CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Meshing & Mesh Conversion (http://www.cfd-online.com/Forums/openfoam-meshing/)
-   -   Specifying boundary faces failes in blockMesh (http://www.cfd-online.com/Forums/openfoam-meshing/75930-specifying-boundary-faces-failes-blockmesh.html)

blaise May 10, 2010 03:56

Specifying boundary faces failes in blockMesh
 
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::printStack(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::polyMesh::polyMesh(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::polyMesh
(
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.


All times are GMT -4. The time now is 18:47.