CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Mesh Utilities (http://www.cfd-online.com/Forums/openfoam-meshing-utilities/)
-   -   Doesnbt have neighbor cells (http://www.cfd-online.com/Forums/openfoam-meshing-utilities/61545-doesnbt-have-neighbor-cells.html)

ivanyao July 15, 2008 23:19

hi when i simulation a cube
 
hi
when i simulation a cube 6*6*6,i define the bottom (0 1 2 3),when i excute the blockmesh and it show bottom (0 1 2 3)doesn't have neighbor cells.my blockmesh.file follow:
arguments "/home/ivan/OpenFOAM/ivan-1.4/run/tutorials/rhoTurbFoam/building" off;

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)
);

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
(
);


// ************************************************** *********************** //
could anyone help me?

ngj July 16, 2008 06:49

Could you please post the erro
 
Could you please post the error message from blockMesh?

Best regards

Niels

ivanyao July 16, 2008 21:14

hi sorry about that the err
 
hi
sorry about that
the error message is:
FOAM FATAL ERROR : face 0 in patch 8 does not have neighbour cell face: 4(0 1 2 3)#0 Foam::error::printStack(Foam:http://www.cfd-online.com/OpenFOAM_D...part/proud.gifstream&)
#1 Foam::error::abort()
#2 Foam::polyMesh::facePatchFaceCells(Foam::List<foam ::face> const&, Foam::List<foam::list<int> > const&, Foam::List<foam::list<foam::face> > const&, int) const
#3 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&)
#4 Foam::blockMesh::createTopology(Foam::IOdictionary &)
#5 Foam::blockMesh::blockMesh(Foam::IOdictionary&)
#6 main
#7 __libc_start_main
#8 __gxx_personality_v0 at /usr/src/packages/BUILD/glibc-2.3/csu/../sysdeps/i386/elf/start.S:122


From function polyMesh::facePatchFaceCells(const faceList& patchFaces,const labelListList& pointCells,const faceListList& cellsFaceShapes,const label patchID)
in file meshes/polyMesh/polyMeshFromShapeMesh.C at line 127.

FOAM aborting

ngj July 17, 2008 04:17

Hi Weihong What you are tol
 
Hi Weihong

What you are told, is that the first face (face no. 0) in patch 8 is wrong. It could either be wrong orientation, but in your case, I cannot see any block which include all the points (0 1 2 3).

The specific patch which cause problems is the wall patch in the very bottom after the empty patches.

/ Niels


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