CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Pre-Processing (https://www.cfd-online.com/Forums/openfoam-pre-processing/)
-   -   How to handle it (https://www.cfd-online.com/Forums/openfoam-pre-processing/62045-how-handle.html)

ivanyao August 20, 2008 02:51

hi, i am doing a simulation a
 
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:http://www.cfd-online.com/OpenFOAM_D...part/proud.gifstream&)
#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 August 20, 2008 04:03

Run blockMesh with -blockTopol
 
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 August 20, 2008 05:39

I have written an application
 
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 August 20, 2008 22:36

hi Mattijs, it can't work.whe
 
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 August 21, 2008 16:21

I just went through your block
 
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 August 22, 2008 21:35

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

ivanyao August 24, 2008 21:27

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


All times are GMT -4. The time now is 16:25.