error messages about blockmeshing and solving

cramsdead · November 10, 2014, 06:53

Good morning everybody !

I´m new on OpenFOAM, and I´m doing my master Thesis on this software. My Topic is about create an RTM simulation and validate the best solver I can use with.

I created a tank (at the bottom), a tube and plate (at the top) which is my composite.

http://i.imgur.com/Z2wUMpW.png => global view
http://i.imgur.com/hPLnPAU.png => top view (Composite plate)
http://i.imgur.com/T5ID7Mx.png => top view zoom (Composite plate)
http://i.imgur.com/SyjQN9l.png => bottom view (tank)

and here my blockMeshDict :

Code:

convertToMeters 0.001;
vertices
(
//---------------composite plate---------------//
//top composite
(0 0 136.5) //0
(150 0 136.5) //1
(300 0 136.5) //2
(300 150 136.5) //3
(300 300 136.5) //4
(150 300 136.5) //5
(0 300 136.5) //6
(0 150 136.5) //7
//top big circle composite
(44.641057 44.641057 136.5) //8
(255.358943 44.641057 136.5) //9
(255.358943 255.358943 136.5) //10
(44.641057 255.358943 136.5) //11
//top little composite hole
 
(145.58058125 145.58058125 136.5) //12
(154.41941875 145.58058125 136.5) //13
(154.41941875 154.41941875 136.5) //14
(145.58058125 154.41941875 136.5) //15
//top little composite square
(147.90290625 147.90290625 136.5) //16
(152.20970938 147.90290625 136.5) //17
(152.20970938 152.20970938 136.5) //18
(147.90290625 152.20970938 136.5) //19
//bottom composite
(0 0 136) //20
(150 0 136) //21
(300 0 136) //22
(300 150 136) //23
(300 300 136) //24
(150 300 136) //25
(0 300 136) //26
(0 150 136) //27
//bottom big circle composite
(44.641057 44.641057 136) //28
(255.358943 44.641057 136) //29
(255.358943 255.358943 136) //30
(44.641057 255.358943 136) //31
//bottom little composite hole (top tube)
 
(145.58058125 145.58058125 136) //32
(154.41941875 145.58058125 136) //33
(154.41941875 154.41941875 136) //34
(145.58058125 154.41941875 136) //35
//bottom little composite square (top tube)
(147.90290625 147.90290625 136) //36
(152.20970938 147.90290625 136) //37
(152.20970938 152.20970938 136) //38
(147.90290625 152.20970938 136) //39
//---------------tank cube---------------//
//top little circle hole (bottom tube)
 
(145.58058125 145.58058125 36) //40
(154.41941875 145.58058125 36) //41
(154.41941875 154.41941875 36) //42
(145.58058125 154.41941875 36) //43
//top little tank square (bottom tube)
(147.90290625 147.90290625 36) //44
(152.20970938 147.90290625 36) //45
(152.20970938 152.20970938 36) //46
(147.90290625 152.20970938 36) //47
//bottom little circle hole
 
(145.58058125 145.58058125 0) //48
(154.41941875 145.58058125 0) //49
(154.41941875 154.41941875 0) //50
(145.58058125 154.41941875 0) //51
//bottom little tank square
(147.90290625 147.90290625 0) //52
(152.20970938 147.90290625 0) //53
(152.20970938 152.20970938 0) //54
(147.90290625 152.20970938 0) //55
//bottom tank cube
(132 132 0) //56
(168 132 0) //57
(168 168 0) //58
(132 168 0) //59
//top tank cube
(132 132 36) //60
(168 132 36) //61
(168 168 36) //62
(132 168 36) //63
　
　
);
blocks
(
//---------------composite plate---------------//
//composite hole circle
hex (11 8 12 15 31 28 32 35) composite (620 5) simpleGrading (1 1 1) //0
hex (8 9 13 12 28 29 33 32) composite (6 20 5) simpleGrading (1 1 1) //1
hex (9 10 14 13 29 30 34 33) composite (6 20 5) simpleGrading (1 1 1) //2
hex (10 11 15 14 30 31 35 34) composite (6 20 5) simpleGrading (1 1 1) //3
hex (12 16 19 15 32 36 39 35) composite (6 6 5) simpleGrading (1 1 1) //4
hex (12 13 17 16 32 33 37 36) composite (6 6 5) simpleGrading (1 1 1) //5
hex (13 14 18 17 33 34 38 37) composite (6 6 5) simpleGrading (1 1 1) //6 
hex (19 18 14 15 39 38 34 35) composite (6 6 5) simpleGrading (1 1 1) //7
//composite cube center
hex (16 17 18 19 36 37 38 39) composite (6 6 5) simpleGrading (1 1 1) //8
//composite sides
hex (6 0 8 11 26 20 28 31) composite (62 5) simpleGrading (1 1 1) //9
hex (0 2 9 8 20 22 29 28) composite (6 2 5) simpleGrading (1 1 1) //10
hex (2 4 10 9 22 24 30 29) composite (6 2 5) simpleGrading (1 1 1) //11
hex (4 6 11 10 24 26 31 30) composite (6 2 5) simpleGrading (1 1 1) //12
　
//---------------tube---------------//
//tube 
hex (35 32 36 39 43 40 44 47) tube (6 6 50) simpleGrading (1 1 1) //13
hex (32 33 37 36 40 41 45 44) tube (6 6 50) simpleGrading (1 1 1) //14
hex (33 34 38 37 41 42 46 45) tube (6 6 50) simpleGrading (1 1 1) //15
hex (34 35 39 38 42 43 47 46) tube (6 6 50) simpleGrading (1 1 1) //16
//tube square center
hex (36 37 38 39 44 45 46 47) tube (6 6 50) simpleGrading (1 1 1) //17
　
//---------------tank---------------//
//tank hole circle
hex (43 40 44 47 51 48 52 55) tank (6 6 10) simpleGrading (1 1 1) //18
hex (40 41 45 44 48 49 53 52) tank (6 6 10) simpleGrading (1 1 1) //19
hex (41 42 46 45 49 50 54 53) tank (6 6 10) simpleGrading (1 1 1) //20
hex (42 43 47 46 50 51 55 54) tank (6 6 10) simpleGrading (1 1 1) //21
//tank square center
hex (44 45 46 47 52 53 54 55) tank (6 6 10) simpleGrading (1 1 1) //22
//tank sides
hex (63 60 40 43 59 5648 51) tank (6 6 10) simpleGrading (1 1 1) //23
hex (60 61 41 40 56 57 49 48) tank (6 6 10) simpleGrading (1 1 1) //24
hex (61 62 42 41 57 58 50 49) tank (6 6 10) simpleGrading (1 1 1) //25
hex (62 63 43 42 58 59 51 50) tank (6 6 10) simpleGrading (1 1 1) //26
 
);
edges
(
　
//big circle top composite
arc 11 8 (1 150 136.5)
arc 8 9 (150 1 136.5)
arc 9 10 (299 150 136.5)
arc 10 11 (150 299 136.5)
//big circle bottom composite
arc 31 28 (1 150 136)
arc 28 29 (150 1 136)
arc 29 30 (299 150 136)
arc 30 31 (150 299 136)
//little circle top composite
arc 15 12 (143.75 150 136.5)
arc 12 13 (150 143.75 136.5)
arc 13 14 (156.25 150 136.5)
arc 14 15 (150 156.25 136.5)
//little circle bottom composite-top cylinder
arc 35 32 (143.75 150 136)
arc 32 33 (150 143.75 136)
arc 33 34 (156.25 150 136)
arc 34 35 (150 156.25 136)
//little circle bottom cylinder-top tank
arc 43 40 (143.75 150 36)
arc 40 41 (150 143.75 36)
arc 41 42 (156.25 150 36)
arc 42 43 (150 156.25 36)
//little circle bottom tank
arc 51 48 (143.75 150 0)
arc 48 49 (150 143.75 0)
arc 49 50 (156.25 150 0)
arc 50 51 (150 156.25 0)
　
);
boundary
(
　
compositeWall
{
type patch;
faces
(
//top
(8 12 15 11)
(8 9 13 12)
(9 10 14 13)
(10 11 15 14)
(0 8 11 6)
(0 2 9 8)
(2 4 10 9)
(4 6 11 10)
(12 13 17 16)
(13 14 18 17)
(14 15 19 18)
(15 12 16 19)
(16 17 18 19)
//bottom
(28 32 35 31)
(28 29 33 32)
(29 30 34 33)
(30 31 35 34)
(20 28 31 26)
(20 22 29 28)
(22 24 30 29)
(24 26 31 30)
//(32 33 37 36)
//(33 34 38 37)
//(34 35 39 38)
//(35 32 36 39)
//(36 37 38 39)
 
);
}
outletWall //side composite wall
{
type patch;
faces
(
(0 2 22 20)
(2 4 24 22)
(4 6 26 24)
(6 0 20 26)
　
);
}
cylinderWall
{
type patch;
faces
(
(35 32 40 43)
(32 33 41 40)
(33 34 42 41)
(34 35 43 42)
);
}
tankWall
{
type patch;
faces
(
 
//top tank wall
(60 61 41 40)
(61 62 42 41)
(62 63 43 42)
(63 60 40 43)
//side tank wall
(56 57 61 60)
(57 58 62 61)
(58 59 63 62)
(59 56 60 63)
 
);
}
inletWall //bottom tank wall
{
type patch;
faces
(
(56 57 49 48)
(57 58 50 49)
(58 59 51 50)
(59 56 48 51)
 
(48 49 53 52)
(49 50 54 53)
(50 51 55 54)
(51 48 52 55)
 
(52 53 54 55)
);
}
);
mergePatchPairs
(
);
// ************************************************************************* //

My Problems are :

during the blockMeshing, I have a lot of Foam Warning Messages :

Code:

Create time
Creating block mesh from
    "/home/lodato/OpenFOAM/lodato-2.2.x/run/RTMSimulationMesh1/constant/polyMesh/blockMeshDict"
Creating curved edges
Creating topology blocks
Creating topology patches
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: -923.415 for face 0
--> 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: -923.415 for face 1
 
 
[...]
 
 
--> FOAM Warning :
    From function blockMesh::createTopology(IOdictionary&)
    in file blockMesh/blockMeshTopology.C at line 255
    negative volume block : 26, probably defined inside-out
Check topology
        Basic statistics
                Number of internal faces : 58
                Number of boundary faces : 46
                Number of defined boundary faces : 46
                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.001
Adding cell zones
    0   composite
    1   tube
    2   tank
Writing cell zones as cellSets
Writing polyMesh
----------------
Mesh Information
----------------
  boundingBox: (0 0 0) (0.3 0.3 0.1365)
  nPoints: 17490
  nCells: 15780
  nFaces: 48972
  nInternalFaces: 45708
----------------
Patches
----------------
  patch 0 (start: 45708 size: 1236) name: compositeWall
  patch 1 (start: 46944 size: 120) name: outletWall
  patch 2 (start: 47064 size: 1200) name: cylinderWall
  patch 3 (start: 48264 size: 384) name: tankWall
  patch 4 (start: 48648 size: 324) name: inletWall
End

and here the checkMesh

Code:

Create time
Create polyMesh for time = 0
Time = 0
Mesh stats
    points:           17490
    faces:            48972
    internal faces:   45708
    cells:            15780
    faces per cell:   6
    boundary patches: 5
    point zones:      0
    face zones:       0
    cell zones:       3
Overall number of cells of each type:
    hexahedra:     15780
    prisms:        0
    wedges:        0
    pyramids:      0
    tet wedges:    0
    tetrahedra:    0
    polyhedra:     0
Checking topology...
    Boundary definition OK.
    Cell to face addressing OK.
    Point usage OK.
    Upper triangular ordering OK.
    Face vertices OK.
    Number of regions: 1 (OK).
Checking patch topology for multiply connected surfaces...
    Patch               Faces    Points   Surface topology
    compositeWall       1236     1273     ok (non-closed singly connected)
    outletWall          120      144      ok (non-closed singly connected)
    cylinderWall        1200     1224     ok (non-closed singly connected)
    tankWall            384      408      ok (non-closed singly connected)
    inletWall           324      337      ok (non-closed singly connected)
Checking geometry...
    Overall domain bounding box (0 0 0) (0.3 0.3 0.1365)
    Mesh (non-empty, non-wedge) directions (1 1 1)
    Mesh (non-empty) directions (1 1 1)
    Boundary openness (-2.46988e-19 1.37608e-18 -6.09032e-16) OK.
 ***High aspect ratio cells found, Max aspect ratio: 1.89983e+196, number of cells 15780
  <<Writing 15780 cells with high aspect ratio to set highAspectRatioCells
    Minimum face area = 5e-08. Maximum face area = 0.000896005.  Face area magnitudes OK.
    Min volume = 2e-300. Max volume = 2e-300.  Total volume = 3.156e-296.  Cell volumes OK.
#0  Foam::error::printStack(Foam::Ostream&) in "/usr/local/IFB/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/usr/local/IFB/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2  __restore_rt at sigaction.c:0
#3  __ieee754_acos at interp.c:0
#4  acos in "/lib64/libm.so.6"
#5  Foam::polyMesh::checkFaceOrthogonality(Foam::Field<Foam::Vector<double> > const&, Foam::Field<Foam::Vector<double> > const&, bool, bool, Foam::HashSet<int, Foam::Hash<int> >*) const in "/usr/local/IFB/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#6  Foam::checkGeometry(Foam::polyMesh const&, bool) in "/usr/local/IFB/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/bin/checkMesh"
#7  main in "/usr/local/IFB/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/bin/checkMesh"
#8  __libc_start_main in "/lib64/libc.so.6"
#9  Foam::UOPstream::write(char) in "/usr/local/IFB/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/bin/checkMesh"
Floating point exception

Then, when I want to compute it by running porousInterFoam, I have this error message :

Code:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.2.x                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.2.x-1ef4d132e582
Exec   : porousInterFoam
Date   : Nov 10 2014
Time   : 10:50:53
Host   : "ideefix.ifb.loc"
PID    : 22284
Case   : /home/lodato/OpenFOAM/lodato-2.2.x/run/RTMSimulationMesh1
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create mesh for time = 0
Reading field p_rgh
Reading field U
Reading/calculating face flux field phi
Reading transportProperties
Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
#0  Foam::error::printStack(Foam::Ostream&) in "/usr/local/IFB/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/usr/local/IFB/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2  __restore_rt at sigaction.c:0
#3  void Foam::mag<Foam::Vector<double> >(Foam::Field<double>&, Foam::UList<Foam::Vector<double> > const&) in "/usr/local/IFB/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libinterfaceProperties.so"
#4  void Foam::mag<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::GeometricField<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/usr/local/IFB/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libinterfaceProperties.so"
#5  Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > Foam::mag<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/usr/local/IFB/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libinterfaceProperties.so"
#6  Foam::interfaceProperties::calculateK() in "/usr/local/IFB/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libinterfaceProperties.so"
#7  Foam::interfaceProperties::interfaceProperties(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::IOdictionary const&) in "/usr/local/IFB/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libinterfaceProperties.so"
#8  main in "/usr/local/IFB/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/bin/porousInterFoam"
#9  __libc_start_main in "/lib64/libc.so.6"
#10  __gxx_personality_v0 in "/usr/local/IFB/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/bin/porousInterFoam"
Floating point exception

I´m pretty sure this error has his root into the blockMeshDict but I don´t know where !

I´m on it know for more than 3 weeks and I´m totally blocked

I used to go on this usefull Forum to fix my problems, but this time I need to ask your help !

This is my first post ever here, I hope it is clear ! I upload the whole file here.

Thank you a lot for help, I begin to be desperate about my job

alexeym · November 10, 2014, 07:29

Hi,

there's explanation in blockMesh output:

Code:

negative volume block : 26, probably defined inside-out

i.e. you've messed up node numbering in block description. Check if order of nodes in block description is in accordance with http://openfoam.org/docs/user/mesh-d...hp#x23-1350072.

cramsdead · November 11, 2014, 11:53

Hi,

Thanks for your help, everything was perfectly good in my blockMeshDict even after 3 checks. I´m may be new but not noob.

So I retried from the begining step by step and I fixed my Problem !

It seems the origin point (0 0 0) has to be a point of your block or it occurs some meshing Errors.

So I just translated my points and here it is !
no more meshing problem and the solver compute it.

November 10, 2014, 07:29		#2
alexeym Senior Member Alexey Matveichev Join Date: Aug 2011 Location: Nancy, France Posts: 1,938 Rep Power: 39	Hi, there's explanation in blockMesh output: Code: negative volume block : 26, probably defined inside-out i.e. you've messed up node numbering in block description. Check if order of nodes in block description is in accordance with http://openfoam.org/docs/user/mesh-d...hp#x23-1350072.

November 11, 2014, 11:53		#3
cramsdead New Member otadol Join Date: Nov 2014 Posts: 14 Rep Power: 12	Hi, Thanks for your help, everything was perfectly good in my blockMeshDict even after 3 checks. I´m may be new but not noob. So I retried from the begining step by step and I fixed my Problem ! It seems the origin point (0 0 0) has to be a point of your block or it occurs some meshing Errors. So I just translated my points and here it is ! no more meshing problem and the solver compute it. Last edited by cramsdead; November 12, 2014 at 04:26.