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

BlockMesh error in channel flow

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 3, 2008, 16:03
Default Hi I am simulating a flow i
  #1
New Member
 
OpenFOAM Newbie
Join Date: Mar 2009
Posts: 20
Rep Power: 8
lofty is on a distinguished road
Hi

I am simulating a flow in a channel with cyclic BC and four ribs inside which i have specified as walls. I am getting the following error when i execute the blockMesh.

Reading block mesh description dictionary

Creating block mesh

Creating blockCorners

Creating curved edges

Creating blocks

Creating patches

Creating block mesh topology
--> FOAM Warning :
From function cellModel::mag(const labelList&, const pointField&)
in file meshes/meshShapes/cellModel/cellModel.C at line 134
zero or negative pyramid volume: -1.33333 for face 4

Default patch type set to empty


--> FOAM FATAL ERROR : face 1 in patch 0 does not have neighbour cell face: 4(0 18 20 43)#0 Foam::error::printStack(Foam:stream&) in "/home/caelinux/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/caelinux/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/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 "/home/caelinux/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#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&, bool) in "/home/caelinux/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#4 Foam::blockMesh::createTopology(Foam::IOdictionary &) in "/home/caelinux/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/blockMesh "
#5 Foam::blockMesh::blockMesh(Foam::IOdictionary&) in "/home/caelinux/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/blockMesh "
#6 main in "/home/caelinux/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/blockMesh "
#7 __libc_start_main in "/lib/i686/libc.so.6"
#8 Foam::regIOobject::readIfModified() in "/home/caelinux/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/blockMesh "


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

Any idea how to resolve the problem! Thanks in advance
lofty is offline   Reply With Quote

Old   March 3, 2008, 17:00
Default Check your lock definition - o
  #2
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,758
Rep Power: 21
hjasak will become famous soon enough
Check your lock definition - one block is inside-out.

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   March 3, 2008, 18:26
Default Sorry what is lock definition
  #3
New Member
 
OpenFOAM Newbie
Join Date: Mar 2009
Posts: 20
Rep Power: 8
lofty is on a distinguished road
Sorry what is lock definition Jasak. What do you mean by "one block is inside out". All the blocks are inside the channel. Here is my block dic file.

FoamFile
{
version 2.0;
format ascii;

root "/opt/OpenFOAM/caelinux-1.4.1/run";
case "cubes";
instance "constant/polyMesh";
local "";

class dictionary;
object blockMeshDict;
}

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

arguments "/opt/OpenFOAM/caelinux-1.4.1/run/cubes" off;

convertToMeters 0.01;

vertices
(
(0 0 0)
(8 0 0)
(0 2 0)
(8 2 0)
(0 8 0)
(8 8 0)
(0 0 8)
(8 0 8)
(0 2 8)
(8 2 8)
(0 8 8)
(8 8 8)
(0 0 6)
(2 0 8)
(2 0 6)
(0 2 6)
(2 2 8)
(2 2 6)
(0 0 2)
(0 0 4)
(2 0 2)
(2 0 4)
(0 2 2)
(0 2 4)
(2 2 2)
(2 2 4)
(4 0 4)
(4 0 6)
(6 0 4)
(6 0 6)
(4 2 4)
(4 2 6)
(6 2 4)
(6 2 6)
(4 0 0)
(4 0 2)
(6 0 0)
(6 0 2)
(4 2 0)
(4 2 2)
(6 2 0)
(6 2 2)
(2 2 0)
(2 0 0)
(4 0 8)
(4 2 8)
(6 0 8)
(6 2 8)
);

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

edges
(
);

patches
(
wall bottomWall
(
(19 12 14 21)
(0 18 20 43)
(43 13 44 34)
(27 44 46 29)
(35 26 28 37)
(36 46 7 1)
)
wall cube1
(
(12 6 13 14)
(15 3 16 17)
(12 14 17 15)
(6 13 16 3)
(12 6 3 15)
(14 13 16 17)
)
wall cube2
(
(18 19 21 20)
(22 23 25 24)
(18 20 24 22)
(19 21 25 23)
(18 19 23 22)
(20 21 25 24)
)
wall cube3
(
(26 27 29 28)
(30 31 33 32)
(26 27 32 30)
(27 29 33 31)
(26 27 31 30)
(28 29 33 32)
)
wall cube4
(
(34 35 37 36)
(38 39 41 40)
(34 36 40 38)
(35 37 41 39)
(34 35 39 38)
(36 37 41 40)
)
wall topWall
(
(4 10 11 5)
)
symmetryPlane sides1
(
(0 2 3 1)
(6 7 9 8)
)
symmetryPlane sides2
(
(2 4 5 3)
(8 9 11 10)
)
cyclic inout1
(
(1 3 9 7)
(0 6 8 2)
)
cyclic inout2
(
(3 5 11 9)
(2 8 10 4)
)
);

mergePatchPairs
(
);


// ************************************************** *********************** //
lofty is offline   Reply With Quote

Old   March 3, 2008, 18:37
Default Dear OF Newbie, The order
  #4
Senior Member
 
Michael Jaworski
Join Date: Mar 2009
Location: Champaign, IL, USA
Posts: 126
Rep Power: 8
mike_jaworski is on a distinguished road
Dear OF Newbie,
The order in which you specify the vertices of a block matter. If you look at the user-guide section on using blockMesh, it will show you the order you should specify the vertices in and how that order defines local coordinates, etc.

The relevant section of the user guide is here:
http://www.opencfd.co.uk/openfoam/do...#x31-1640006.3

You'll have to examine your block definitions for which specific block is incorrect.

Good Luck,
Mike J.
mike_jaworski is offline   Reply With Quote

Old   March 5, 2008, 09:13
Default Hi Can anyone explain what
  #5
New Member
 
OpenFOAM Newbie
Join Date: Mar 2009
Posts: 20
Rep Power: 8
lofty is on a distinguished road
Hi

Can anyone explain what the following error imply:

FOAM FATAL ERROR : face 0 in patch 11 does not have neighbour cell face: 4(12 14 13 6)#0

I do have a face opposite to face mentioned in error. Why is blockMesh not recognising it? Thank You
lofty 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
BlockMesh error hjasak OpenFOAM Native Meshers: blockMesh 5 August 24, 2012 07:22
maintain a constant flow rate in channel flow ? Lewis Main CFD Forum 2 September 28, 2010 11:35
Error while blockMesh fightingfalcon23 OpenFOAM Native Meshers: blockMesh 0 April 15, 2008 05:41
Blockmesh cavity error message tonitoney OpenFOAM Installation 2 March 17, 2008 12:59
BlockMesh error with growing mesh size kian OpenFOAM Native Meshers: blockMesh 4 September 24, 2007 16:00


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