CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Native Meshers: blockMesh

blockMesh problem with blocks with vertices

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   March 21, 2016, 12:25
Default blockMesh problem with blocks with vertices
  #1
Member
 
Adrien ORSINI
Join Date: Mar 2016
Location: France
Posts: 85
Rep Power: 2
adrieno is on a distinguished road
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:rintStack(Foam::Ostream&) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 Foam:olyMesh::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:olyMesh::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:olyMesh:olyMesh(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 is offline   Reply With Quote

Old   March 22, 2016, 04:21
Post
  #2
Member
 
Adrien ORSINI
Join Date: Mar 2016
Location: France
Posts: 85
Rep Power: 2
adrieno is on a distinguished road
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).
Attached Images
File Type: jpg DSC_0480.jpg (90.7 KB, 15 views)
Attached Files
File Type: c blockMeshDict.c (2.0 KB, 9 views)
adrieno is offline   Reply With Quote

Old   March 22, 2016, 10:52
Angry
  #3
Member
 
Adrien ORSINI
Join Date: Mar 2016
Location: France
Posts: 85
Rep Power: 2
adrieno is on a distinguished road
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:rintStack(Foam::Ostream&) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 Foam:olyMesh::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:olyMesh::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:olyMesh:olyMesh(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...
Attached Images
File Type: jpg DSC_0483.jpg (67.1 KB, 5 views)
Attached Files
File Type: txt blockMeshDict.txt (1.9 KB, 2 views)
adrieno is offline   Reply With Quote

Old   March 23, 2016, 04:35
Default blockMesh cylinder then square.
  #4
Member
 
Adrien ORSINI
Join Date: Mar 2016
Location: France
Posts: 85
Rep Power: 2
adrieno is on a distinguished road
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
Attached Images
File Type: png 3.png (21.2 KB, 6 views)
File Type: png 4.png (20.9 KB, 6 views)
File Type: png 5.png (22.4 KB, 5 views)
File Type: png 6.png (33.0 KB, 5 views)
Attached Files
File Type: txt blockMeshDict.txt (2.7 KB, 4 views)
adrieno is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Problem with blockMesh and my shape TneurolF OpenFOAM Native Meshers: blockMesh 4 June 25, 2013 13:52
[ICEM] Problem with blocks in 3D Mesh brator ANSYS Meshing & Geometry 3 June 25, 2013 10:38
Blockmesh problem with more than one block sven82 OpenFOAM Pre-Processing 1 June 4, 2013 17:08
Is there any easy way to connect blocks of varying resolutions using blockMesh? murrdpirate OpenFOAM 7 June 4, 2011 09:26
blockMesh problem shangzung OpenFOAM 5 December 9, 2009 09:18


All times are GMT -4. The time now is 14:27.