CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Meshing & Mesh Conversion

Specifying boundary faces failes in blockMesh

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   May 10, 2010, 03:56
Default Specifying boundary faces failes in blockMesh
  #1
Member
 
P.A.
Join Date: Mar 2009
Location: Germany
Posts: 47
Rep Power: 8
blaise is on a distinguished road
Hi Foamers,

I have a problem with the grid generation using blockMesh. With the setup included below I get the following output and error message from blockMesh:

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time


Creating block mesh from
"/daten/ShipLES/geom_bodies/quader/LES/calc/snappy_LES/constant/polyMesh/blockMeshDict"


Creating blockCorners

Creating curved edges

Creating blocks

Creating patches

Creating block mesh topology

Default patch type set to empty


Trying to specify a boundary face 4(2 3 7 6) on the face on cell 0 which is either an internal face or already belongs to some other patch. This is face 0 of patch 3 named Q_TOP.#0 Foam::error:rintStack(Foam::Ostream&) in "/usr/local/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/usr/local/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 Foam:olyMesh:olyMesh(Foam::IOobject const&, Foam::Xfer<Foam::Field<Foam::Vector<double> > > const&, Foam::List<Foam::cellShape> const&, Foam::List<Foam::List<Foam::face> > const&, Foam::List<Foam::word> const&, Foam::List<Foam::word> const&, Foam::word const&, Foam::word const&, Foam::List<Foam::word> const&, bool) in "/usr/local/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libOpenFOAM.so"
#3 Foam::blockMesh::createTopology(Foam::IOdictionary &) in "/usr/local/OpenFOAM/OpenFOAM-1.6.x/applications/bin/linux64GccDPOpt/blockMesh"
#4 Foam::blockMesh::blockMesh(Foam::IOdictionary&) in "/usr/local/OpenFOAM/OpenFOAM-1.6.x/applications/bin/linux64GccDPOpt/blockMesh"
#5 main in "/usr/local/OpenFOAM/OpenFOAM-1.6.x/applications/bin/linux64GccDPOpt/blockMesh"
#6 __libc_start_main in "/lib64/libc.so.6"
#7 Foam::regIOobject::writeObject(Foam::IOstream::str eamFormat, Foam::IOstream::versionNumber, Foam::IOstream::compressionType) const in "/usr/local/OpenFOAM/OpenFOAM-1.6.x/applications/bin/linux64GccDPOpt/blockMesh"


From function polyMesh:olyMesh
(
const IOobject&,
const Xfer<pointField>&,
const cellShapeList& cellsAsShapes,
const faceListList& boundaryFaces,
const wordList& boundaryPatchTypes,
const wordList& boundaryPatchNames,
const word& defaultBoundaryPatchType
)
in file meshes/polyMesh/polyMeshFromShapeMesh.C at line 484.

FOAM aborting

Aborted

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //


What am I doing wrong here? Any hint is appreciated!

When I look at the blocking with simpleViewer, everything seems to be well arranged.
The blockMeshDict is as follows (it is a cuboid placed in the middle of a wide cylinder; I hope everything is commented sufficiently):

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: http://www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
// Cuboid in Cylinder

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

convertToMeters 0.001;


// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

vertices
(
// Base Cuboid
(-70.8 -70.8 0) // 0
(70.8 -70.8 0) // 1
(70.8 70.8 0) // 2
(-70.8 70.8 0) // 3

// Lid Cuboid
(-70.8 -70.8 283) // 4
(70.8 -70.8 283) // 5
(70.8 70.8 283) // 6
(-70.8 70.8 283) // 7

// Empty Cuboid
(-70.8 -70.8 850) // 8
(70.8 -70.8 850) // 9
(70.8 70.8 850) // 10
(-70.8 70.8 850) // 11

// Cylinder base
(1414.2 1414.2 0) // 12
(0 2000 0) // 13
(-1414.2 1414.2 0) // 14
(-2000 0 0) // 15
(-1414.2 -1414.2 0) // 16
(0 -2000 0) // 17
(1414.2 -1414.2 0) // 18
(2000 0 0) // 19

// Middle plane cylinder
(1414.2 1414.2 283) // 20
(0 2000 283) // 21
(-1414.2 1414.2 283) // 22
(-2000 0 283) // 23
(-1414.2 -1414.2 283) // 24
(0 -2000 283) // 25
(1414.2 -1414.2 283) // 26
(2000 0 283) // 27

// Cylinder Top
(1414.2 1414.2 850) // 28
(0 2000 850) // 29
(-1414.2 1414.2 850) // 30
(-2000 0 850) // 31
(-1414.2 -1414.2 850) // 32
(0 -2000 850) // 33
(1414.2 -1414.2 850) // 34
(2000 0 850) // 35
);

blocks
(
// Cuboid
hex (0 1 2 3 4 5 6 7)
quader
(5 5 5)
simpleGrading (1 1 1)

// Empty Cuboid
hex (4 5 6 7 8 9 10 11)
(5 5 5)
simpleGrading (1 1 1)

// North-Quadrant on Cylinder Basis
hex (3 2 12 14 7 6 20 22)
(5 5 5)
simpleGrading (1 1 1)

// West-Quadrant on Cylinder Basis
hex (0 3 14 16 4 7 22 24)
(5 5 5)
simpleGrading (1 1 1)

// South-Quadrant on Cylinder Basis
hex (1 0 16 18 5 4 24 26)
(5 5 5)
simpleGrading (1 1 1)

// East-Quadrant on Cylinder Basis
hex (2 1 18 12 6 5 26 20)
(5 5 5)
simpleGrading (1 1 1)

// North-Quadrant on Cylinder Top
hex (7 6 20 22 11 10 28 30)
(5 5 5)
simpleGrading (1 1 1)

// West-Quadrant on Cylinder Top
hex (4 7 22 24 8 11 30 32)
(5 5 5)
simpleGrading (1 1 1)

// South-Quadrant on Cylinder Top
hex (5 4 24 26 9 8 32 34)
(5 5 5)
simpleGrading (1 1 1)

// East-Quadrant on Cylinder Top
hex (6 5 26 20 10 9 34 28)
(5 5 5)
simpleGrading (1 1 1)
);

edges
(
// Arc Basis Cylinder
arc 12 14 (0 2000 0) // arc 0
arc 14 16 (-2000 0 0) // arc 1
arc 16 18 (0 -2000 0) // arc 3
arc 18 12 (2000 0 0) // arc 4

// Arc middle plane
arc 20 22 (0 2000 283) // arc 5
arc 22 24 (-2000 0 283) // arc 6
arc 24 26 (0 -2000 283) // arc 7
arc 26 20 (2000 0 283) // arc 8

// Arc Top Cylinder
arc 28 30 (0 2000 850) // arc 9
arc 30 32 (-2000 0 850) // arc 10
arc 32 34 (0 -2000 850) // arc 11
arc 34 28 (2000 0 850) // arc 12
);

patches
(
// Cylinder shell
patch ZYL_AUSSEN
(
// Cylinder shell, on cylinder Basis
(12 14 22 20) // North
(14 16 24 22) // West
(16 18 26 24) // South
(18 12 20 26) // East

// Cylinder shell, on Cylinder Top
(20 22 30 28) // North
(22 24 32 30) // West
(24 26 34 32) // South
(26 20 28 34) // East
)

// Cylinder Top
patch ZYL_TOP
(
(28 30 11 10) // North
(30 32 8 11) // West
(32 34 9 8) // South
(34 28 10 9) // East
(8 9 10 11) // Center Cylinder Top
)

// Cylinder Basis
patch ZYL_BASE
(
(2 3 14 12) // North
(3 0 16 14) // West
(0 1 18 16) // South
(1 2 12 18) // East
(3 2 1 0) // Center Cylinder Basis
)

wall Q_TOP
(
(2 3 7 6)
)

wall Q_FRONT
(
(1 2 6 5)
)

wall Q_BOTTOM
(
(0 1 5 4)
)

wall Q_BACK
(
(0 4 7 3)
)

wall Q_LID
(
(4 5 6 7)
)

);

// ************************************************** *********************** //

Cheers,

Pascal.
blaise is offline   Reply With Quote

Reply

Tags
blockmeshdict

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Implementation of boundary conditions for FVM Tom Main CFD Forum 7 August 26, 2014 05:58
RPM in Wind Turbine Pankaj CFX 9 November 23, 2009 05:05
Problem in IMPORT of ICEM input file in FLUENT csvirume FLUENT 2 September 9, 2009 01:08
Unaligned accesses on IA64 andre OpenFOAM 5 June 23, 2008 10:37
Boundary conditions? Tom Main CFD Forum 0 November 5, 2002 02:54


All times are GMT -4. The time now is 03:42.