cramsdead |
November 10, 2014 05:53 |
error messages about blockmeshing and solving
1 Attachment(s)
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 ! :confused:
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 ! :confused:
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 :(
|