CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM (https://www.cfd-online.com/Forums/openfoam/)
-   -   Problem related to topology (https://www.cfd-online.com/Forums/openfoam/101324-problem-related-topology.html)

Ank May 1, 2012 06:31
Problem related to topology
 
Hey
I am trying to mesh a rectangular tank with 5 pipes, but I am having the following problem with my blockMesh :


/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 1.6-f802ff2d6c5a
Exec : blockMesh
Date : May 01 2012
Time : 15:27:17
Host : karya101.hbni.ac.in
PID : 14790
Case : /home/usr0201/ankur_a/hotRoom
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time


Creating block mesh from
"/home/usr0201/ankur_a/hotRoom/constant/polyMesh/blockMeshDict"


Creating blockCorners

Creating curved edges

Creating blocks
--> FOAM Warning :
From function bool Foam::blockMesh::blockLabelsOK(const label blockLabel, co nst pointField& points, const cellShape& blockShape)
in file createTopology.C at line 67
block 31 point label 112 larger than 111 the largest defined point label

Creating patches


Cannot create mesh due to errors in topology, exiting !


From function blockMesh::createTopology(IOdictionary&)
in file createTopology.C at line 376.

FOAM exiting



I am new to the openfoam so could not figure out the problem, can anyone please help me..

winden May 1, 2012 08:02
Hi.

That means that you have defined 112 points in your vertex list (so the last one has the label 111) but in block 31 you have tried to access point number 113 (with label 112.) Remember that the first point/vertex is number 0 not number 1.

//Björn

Ank May 2, 2012 04:02
hey thanks alot..it worked but now I am having a new error.. I just could not debug it..if somone know the solution please let me know..






/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 1.6-f802ff2d6c5a
Exec : blockMesh
Date : May 02 2012
Time : 13:13:20
Host : karya101.hbni.ac.in
PID : 28798
Case : /home/usr0201/ankur_a/hotRoom
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time


Creating block mesh from
"/home/usr0201/ankur_a/hotRoom/constant/polyMesh/blockMeshDict"


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 129
zero or negative pyramid volume: -1.17853e-06 for face 1
--> FOAM Warning :
From function cellModel::mag(const labelList&, const pointField&)
in file meshes/meshShapes/cellModel/cellModel.C at line 129
zero or negative pyramid volume: -7.07108e-06 for face 2
--> FOAM Warning :
From function cellModel::mag(const labelList&, const pointField&)
in file meshes/meshShapes/cellModel/cellModel.C at line 129
zero or negative pyramid volume: -9.46363e-06 for face 0
--> FOAM Warning :
From function cellModel::mag(const labelList&, const pointField&)
in file meshes/meshShapes/cellModel/cellModel.C at line 129
zero or negative pyramid volume: -0.000270357 for face 1
--> FOAM Warning :
From function cellModel::mag(const labelList&, const pointField&)
in file meshes/meshShapes/cellModel/cellModel.C at line 129
zero or negative pyramid volume: -0.000119752 for face 2
--> FOAM Warning :
From function cellModel::mag(const labelList&, const pointField&)
in file meshes/meshShapes/cellModel/cellModel.C at line 129
zero or negative pyramid volume: -0.000160068 for face 3
--> FOAM Warning :
From function cellModel::mag(const labelList&, const pointField&)
in file meshes/meshShapes/cellModel/cellModel.C at line 129
zero or negative pyramid volume: -0.00013991 for face 4
--> FOAM Warning :
From function cellModel::mag(const labelList&, const pointField&)
in file meshes/meshShapes/cellModel/cellModel.C at line 129
zero or negative pyramid volume: -0.00013991 for face 5
--> FOAM Warning :
From function blockMesh::createTopology(IOdictionary&)
in file createTopology.C at line 398
negative volume block : 25, probably defined inside-out
--> FOAM Warning :
From function cellModel::mag(const labelList&, const pointField&)
in file meshes/meshShapes/cellModel/cellModel.C at line 129
zero or negative pyramid volume: -2.24473e-05 for face 0
--> FOAM Warning :
From function cellModel::mag(const labelList&, const pointField&)
in file meshes/meshShapes/cellModel/cellModel.C at line 129
zero or negative pyramid volume: -2.33901e-05 for face 1
--> FOAM Warning :
From function cellModel::mag(const labelList&, const pointField&)
in file meshes/meshShapes/cellModel/cellModel.C at line 129
zero or negative pyramid volume: -4.72524e-05 for face 3
--> FOAM Warning :
From function cellModel::mag(const labelList&, const pointField&)
in file meshes/meshShapes/cellModel/cellModel.C at line 129
zero or negative pyramid volume: -2.29187e-05 for face 4
--> FOAM Warning :
From function cellModel::mag(const labelList&, const pointField&)
in file meshes/meshShapes/cellModel/cellModel.C at line 129
zero or negative pyramid volume: -2.29187e-05 for face 5
--> FOAM Warning :
From function blockMesh::createTopology(IOdictionary&)
in file createTopology.C at line 398
negative volume block : 26, probably defined inside-out

Default patch type set to empty


face 10 in patch 0 does not have neighbour cell face: 4(30 31 87 86)#0 Foam::error::printStack(Foam::Ostream&) in "/home/usr0201/ankur_a/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/usr0201/ankur_a/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 Foam::Ostream& Foam::operator<< <Foam::error>(Foam::Ostream&, Foam::errorManip<Foam::error>) in "/home/usr0201/ankur_a/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/blockMesh"
#3 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/usr0201/ankur_a/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/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::List<Foam::word> const&, Foam::word const&, Foam::word const&, Foam::List<Foam::word> const&, bool) in "/home/usr0201/ankur_a/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#5 Foam::blockMesh::createTopology(Foam::IOdictionary &) in "/home/usr0201/ankur_a/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/blockMesh"
#6 Foam::blockMesh::blockMesh(Foam::IOdictionary&) in "/home/usr0201/ankur_a/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/blockMesh"
#7 main in "/home/usr0201/ankur_a/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/blockMesh"
#8 __libc_start_main in "/lib64/libc.so.6"
#9 _start at /usr/src/packages/BUILD/glibc-2.9/csu/../sysdeps/x86_64/elf/start.S:116


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

FOAM aborting

Aborted

winden May 2, 2012 05:08
That probably means that you have defined a block inside out, have you read the description of how to arrange the vertices within a block? It is very important that the vertices come in the right order when defining the blocks and the patches.

http://www.openfoam.org/docs/user/blockMesh.php

//Björn


All times are GMT -4. The time now is 11:51.