CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Meshing & Mesh Conversion (https://www.cfd-online.com/Forums/openfoam-meshing/)
-   -   [blockMesh] blockMesh problem with blocks with vertices (https://www.cfd-online.com/Forums/openfoam-meshing/168408-blockmesh-problem-blocks-vertices.html)

adrieno March 21, 2016 12:25

blockMesh problem with blocks with vertices
 
hello all,
I'm trying a mesh with block with only 6 vertices, but I allways have this error message

--> FOAM FATAL ERROR:
face 0 in patch 1 does not have neighbour cell face: 4(1 3 6 1)

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

#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 Foam::polyMesh::facePatchFaceCells(Foam::List<Foam ::face> const&, Foam::List<Foam::List<int> > const&, Foam::List<Foam::List<Foam::face> > const&, int) const in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#3 Foam::polyMesh::setTopology(Foam::List<Foam::cellS hape> const&, Foam::List<Foam::List<Foam::face> > const&, Foam::List<Foam::word> const&, Foam::List<int>&, Foam::List<int>&, int&, int&, Foam::List<Foam::cell>&) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4 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::PtrList<Foam::dictionary> const&, Foam::word const&, Foam::word const&, bool) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5 Foam::blockMesh::createTopology(Foam::IOdictionary const&, Foam::word const&) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libblockMesh.so"
#6 Foam::blockMesh::blockMesh(Foam::IOdictionary const&, Foam::word const&) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libblockMesh.so"
#7
in "/opt/openfoam222/platforms/linux64GccDPOpt/bin/blockMesh"
#8 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9
in "/opt/openfoam222/platforms/linux64GccDPOpt/bin/blockMesh"
Abandon (core dumped)

Any Idea ?

Thanks




adrieno March 22, 2016 04:21

2 Attachment(s)
Here is a photo of the geometry expected with the vertices and here is also the blockMeshDict (I juste add ".c" to upload the file, otherwise it didn't worked).

adrieno March 22, 2016 10:52

2 Attachment(s)
Hi all,

I've simplified my geometry as possible in order to find my error.

The thing is, with only 2 blocks, I still have this fatal error.


Creating block mesh topology
--> FOAM Warning :
From function cellModel::mag(const labelList&, const pointField&)
in file meshes/meshShapes/cellModel/cellModel.C at line 128
zero or negative pyramid volume: -2.04827e+08 for face 1
--> FOAM Warning :
From function cellModel::mag(const labelList&, const pointField&)
in file meshes/meshShapes/cellModel/cellModel.C at line 128
zero or negative pyramid volume: -2.04827e+08 for face 2
--> FOAM Warning :
From function cellModel::mag(const labelList&, const pointField&)
in file meshes/meshShapes/cellModel/cellModel.C at line 128
zero or negative pyramid volume: -1.69684e+08 for face 0

Check topology

Basic statistics
Number of internal faces : 1
Number of boundary faces : 10
Number of defined boundary faces : 10
Number of undefined boundary faces : 0
Checking patch -> block consistency

Creating block offsets
Creating merge list .

Creating polyMesh from blockMesh
Creating patches
Creating cells
Creating points with scale 0.1


--> FOAM FATAL ERROR:
face 100 in patch 0 does not have neighbour cell face: 4(0 111 122 1)

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

#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 Foam::polyMesh::facePatchFaceCells(Foam::List<Foam ::face> const&, Foam::List<Foam::List<int> > const&, Foam::List<Foam::List<Foam::face> > const&, int) const in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#3 Foam::polyMesh::setTopology(Foam::List<Foam::cellS hape> const&, Foam::List<Foam::List<Foam::face> > const&, Foam::List<Foam::word> const&, Foam::List<int>&, Foam::List<int>&, int&, int&, Foam::List<Foam::cell>&) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4 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::PtrList<Foam::dictionary> const&, Foam::word const&, Foam::word const&, bool) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5
in "/opt/openfoam222/platforms/linux64GccDPOpt/bin/blockMesh"
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7
in "/opt/openfoam222/platforms/linux64GccDPOpt/bin/blockMesh"
Abandon (core dumped)



I've joined a picture of my geometry expected (a draw of the face up z=0 and down z=3976) and my blockMeshDict file.

(In fact I would like at the end to mesh a cylinder, but as it doesn't work I go step by step...)

I've allready done a simple mesh with two blocks, so I don't understand why this one doesn't work at all...

adrieno March 23, 2016 04:35

blockMesh cylinder then square.
 
5 Attachment(s)
Hi all,

I've fixed my problems.
In fact I don't really fixed it but I've started from the very beginning... one block at a time as show on the pictures.

Because that's a quit standard case I think, I add the blockMeshDict if someone want's to work on it, adapting his own dimension..

By the way if someone has any other idea of blocking for this geometry, let me know. I'm really new, so there must be many ways to do it better.

Have a nice project !
Adrien


All times are GMT -4. The time now is 06:09.