[ivanyao]

hi,
i am doing a simulation about wind pressure on a 6m cube,i put the 6m*6m*6m cube into 192m*120m*60m box.when i extrudeMesh,it show follow message;
FOAM FATAL ERROR : Trying to specify a boundary face 4(0 3 7 4) on the face on cell 14 which is either an internal face or already belongs to some other patch. This is face 0 of patch 2 named front.#0 Foam::error::printStack(Foam:stream&)
#1 Foam::error::abort()
#2 Foam::polyMesh::polyMesh(Foam::IOobject const&, 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::List<foam::word> const&)
#3 Foam::blockMesh::createTopology(Foam::IOdictionary &)
#4 Foam::blockMesh::blockMesh(Foam::IOdictionary&)
#5 main
#6 __libc_start_main
#7 __gxx_personality_v0 at /usr/src/packages/BUILD/glibc-2.3/csu/../sysdeps/i386/elf/start.S:122


From function polyMesh::polyMesh
(
const IOobject& io,
const pointField& points,
const cellShapeList& cellsAsShapes,
const faceListList& boundaryFaces,
const wordList& boundaryPatchTypes,
const wordList& boundaryPatchNames,
const word& defaultBoundaryPatchType
)
in file meshes/polyMesh/polyMeshFromShapeMesh.C at line 482.
FOAM aborting

my blockMesh file are:
convertToMeters 1;

vertices
(
(-3 -3 0)
(3 -3 0)
(3 3 0)
(-3 3 0)
(-3 -3 6)
(3 -3 6)
(3 3 6)
(-3 3 6)
(-39 -3 0)
(-39 3 0)
(-39 -3 6)
(-39 3 6)
(153 -3 0)
(153 3 0)
(153 -3 6)
(153 3 6)
(-3 -60 0)
(3 -60 0)
(-3 -60 6)
(3 -60 6)
(-3 60 0)
(3 60 0)
(-3 60 6)
(3 60 6)
(-39 -60 0)
(153 -60 0)
(153 60 0)
(-39 60 0)
(-39 -60 6)
(153 -60 6)
(153 60 6)
(-39 60 6)
(-39 -60 60)
(153 -60 60)
(153 60 60)
(-39 60 60)
(-39 -3 60)
(-39 3 60)
(-3 -60 60)
(3 -60 60)
(153 -3 60)
(153 3 60)
(-3 60 60)
(3 60 60)
(-3 -3 60)
(3 -3 60)
(3 3 60)
(-3 3 60)
);

blocks
(
hex (24 16 0 8 28 18 4 10) (36 57 6) simpleGrading (1 1 1)
hex (28 18 4 10 32 38 44 36) (36 57 54) simpleGrading (1 1 1)
hex (16 17 1 0 18 19 5 4) (6 57 6) simpleGrading (1 1 1)
hex (18 19 5 4 38 39 45 44) (6 57 54) simpleGrading (1 1 1)
hex (17 25 12 1 19 29 14 5) (150 57 6) simpleGrading (1 1 1)
hex (19 29 14 5 39 33 40 45) (150 57 54) simpleGrading (1 1 1)
hex (1 12 13 2 5 14 15 6) (150 6 6) simpleGrading (1 1 1)
hex (5 14 15 6 45 40 41 46) (150 6 54) simpleGrading (1 1 1)
hex (2 13 26 21 6 15 30 23) (150 57 6) simpleGrading (1 1 1)
hex (6 15 30 23 46 41 34 43) (150 57 54) simpleGrading (1 1 1)
hex (3 2 21 20 7 6 23 22) (6 57 6) simpleGrading (1 1 1)
hex (7 6 23 22 47 46 43 42) (6 57 54) simpleGrading (1 1 1)
hex (9 3 20 27 11 7 22 31) (36 57 6) simpleGrading (1 1 1)
hex (11 7 22 31 37 47 42 35) (36 57 54) simpleGrading (1 1 1)
hex (8 0 3 9 10 4 7 11) (36 6 6) simpleGrading (1 1 1)
hex (10 4 7 11 36 44 47 37) (36 6 54) simpleGrading (1 1 1)
hex (4 5 6 7 44 45 46 47) (6 6 54) simpleGrading (1 1 1)
hex (0 1 2 3 4 5 6 7) (6 6 6) edgeGrading (1 1 1 1 1 1 1 1 1 1 1 1)
);

edges
(
);

patches
(
patch inlet
(
(27 9 11 31)
(37 11 31 35)
(9 8 10 11)
(11 10 36 37)
(8 24 28 10)
(10 28 32 36)
)
patch outlet
(
(25 12 14 29)
(29 14 40 33)
(12 13 15 14)
(14 15 41 40)
(13 26 30 15)
(15 30 34 41)
)
wall front
(
(0 3 7 4)
)
wall <back>
(
(1 2 6 5)
)
wall <left>
(
(0 1 5 4)
)
wall <right>
(
(2 3 7 6)
)
wall <top>
(
(4 5 6 7)
)
empty <wall>
(
(24 16 18 28)
(28 18 38 32)
(16 17 19 18)
(18 19 39 38)
(17 25 29 19)
(19 29 33 39)
(26 21 23 30)
(30 23 43 34)
(21 20 22 23)
(23 22 42 43)
(20 27 31 22)
(22 31 35 42)
(32 38 44 36)
(44 36 37 47)
(47 37 35 42)
(38 39 45 44)
(44 45 46 47)
(47 46 43 42)
(39 33 40 45)
(40 45 46 41)
(41 46 43 34)
(16 24 8 0)
(0 8 9 3)
(3 9 27 20)
(17 16 0 1)
(2 3 20 21)
(25 17 1 12)
(12 1 2 13)
(13 2 21 26)
)
wall <bottom>
(
(0 1 2 3)
)
);

mergePatchPairs
(
);
would anyone help me?

[mattijs]

Run blockMesh with -blockTopology option. This will write a .obj file which can be visualised with e.g. javaview or (with objToVTK) with paraview. Both have the option of displaying vertex numbers (with selecting and selection inspector in paraview). Very useful to debug blockMeshes.

[grtabor]

I have written an application called simpleViewer based on blockMesh, which visualises the block structure of a blockMeshDict file for you. Uses FOX/OpenGL. It works for OF1.4.1; I'm guessing making it work for 1.5 would not be too hard. Email me at g.r.tabor@ex.ac.uk if you would like to try with this.

Gavin

[ivanyao]

hi Mattijs,
it can't work.when i type blockMesh <root> <case> -blockTopology.it show:
FOAM FATAL ERROR : Trying to specify a boundary face 4(0 3 7 4) on the face on cell 14 which is either an internal face or already belongs to some other patch. This is face 0 of patch 2 named front.#0 Foam::error::printStack(Foam:stream&)
#1 Foam::error::abort()
#2 Foam::polyMesh::polyMesh(Foam::IOobject const&, 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::List<foam::word> const&)
#3 Foam::blockMesh::createTopology(Foam::IOdictionary &)
#4 Foam::blockMesh::blockMesh(Foam::IOdictionary&)
#5 main
#6 __libc_start_main
#7 __gxx_personality_v0 at /usr/src/packages/BUILD/glibc-2.3/csu/../sysdeps/i386/elf/start.S:122


From function polyMesh::polyMesh
(
const IOobject& io,
const pointField& points,
const cellShapeList& cellsAsShapes,
const faceListList& boundaryFaces,
const wordList& boundaryPatchTypes,
const wordList& boundaryPatchNames,
const word& defaultBoundaryPatchType
)
in file meshes/polyMesh/polyMeshFromShapeMesh.C at line 482.
FOAM aborting
so how do i fix it?

[mattijs]

I just went through your blockMeshDict. Set all block sizes to (1 1 1) so there are no matching errors. Then commented out all patch specifications (front, back, left, right, top) that give errors. Now with forementioned -blockTopology start putting patches and block sizes back in.

[ivanyao]

hi,
thank you!but now what can I do?

[ivanyao]

hi Mattijs,
i am poor in English.so i don't know what to do.
would you tell me the step in detail.please...