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

Possible Bug in pimpleFoam (or createPatch) (or fluent3DMeshToFoam)

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 21, 2013, 09:21
Angry Possible Bug in pimpleFoam (or createPatch) (or fluent3DMeshToFoam)
  #1
Senior Member
 
cfdonline2mohsen's Avatar
 
Mohsen KiaMansouri
Join Date: Jan 2010
Location: CFD Lab
Posts: 118
Rep Power: 7
cfdonline2mohsen is on a distinguished road
Hello to all Foamers

I'm struggling with a problem for about 2 weeks.
I want to simulate a channel flow. I created my mesh in GAMBIT then imported it to OpenFOAM and ran it successfully by using pimpleFoam. after a while, I created the same geometry in GAMBIT and again imported it to OpenFOAM but I couldn't run it this time by pimpleFoam??

I explain my problem in details.
I created two channel geometries with the same geometry and mesh in GAMBIT (TwinLundWithCyclic1.msh and TwinLundTwin.msh)and then import them to OpenFOAM:

1st Step: fluent3DMeshToFoam for My 1st Case (TwinLundWithCyclic2):
Code:
moh3n@moh3n:~$ cd /home/moh3n/Desktop/920729/TwinLundWithCyclic2
moh3n@moh3n:~/Desktop/920729/TwinLundWithCyclic2$ fluent3DMeshToFoam TwinLundWithCyclic1.msh
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.1.1-221db2718bbb
Exec   : fluent3DMeshToFoam TwinLundWithCyclic1.msh
Date   : Oct 21 2013
Time   : 10:48:15
Host   : "moh3n"
PID    : 3029
Case   : /home/moh3n/Desktop/920729/TwinLundWithCyclic2
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

--> FOAM Warning : 
    From function dlOpen(const fileName&, const bool)
    in file POSIX.C at line 1175
    dlopen error : /home/moh3n/OpenFOAM/moh3n-2.1.1/platforms/linux64GccDPOpt/lib/libdynamicSmagorinskyModel.so: undefined symbol: _ZNK4Foam14incompressible9LESModels11GenEddyVisc1BEv
--> FOAM Warning : 
    From function dlLibraryTable::open(const fileName&, const bool)
    in file db/dynamicLibrary/dlLibraryTable/dlLibraryTable.C at line 96
    could not load "libdynamicSmagorinskyModel.so"
Dimension of grid: 3
Number of points: 325710
PointGroup: 1 start: 0 end: 325709.  Reading points...done.
Number of faces: 941323
FaceGroup: 3 start: 0 end: 2989.  Reading mixed faces...done.
FaceGroup: 4 start: 2990 end: 5979.  Reading mixed faces...done.
FaceGroup: 5 start: 5980 end: 6853.  Reading mixed faces...done.
FaceGroup: 6 start: 6854 end: 10717.  Reading mixed faces...done.
FaceGroup: 7 start: 10718 end: 11591.  Reading mixed faces...done.
FaceGroup: 8 start: 11592 end: 15455.  Reading mixed faces...done.
FaceGroup: 9 start: 15456 end: 16690.  Reading mixed faces...done.
FaceGroup: 10 start: 16691 end: 22150.  Reading mixed faces...done.
FaceGroup: 11 start: 22151 end: 23385.  Reading mixed faces...done.
FaceGroup: 12 start: 23386 end: 28845.  Reading mixed faces...done.
FaceGroup: 13 start: 28846 end: 31835.  Reading mixed faces...done.
FaceGroup: 14 start: 31836 end: 34825.  Reading mixed faces...done.
FaceGroup: 16 start: 34826 end: 941322.  Reading mixed faces...done.
Number of cells: 307970
CellGroup: 2 start: 0 end: 307969 type: 1
Zone: 2 name: fluid type: fluid.  Reading zone data...done.
Zone: 3 name: interface1 type: wall.  Reading zone data...done.
Zone: 4 name: inteface0 type: wall.  Reading zone data...done.
Zone: 5 name: back1 type: symmetry.  Reading zone data...done.
Zone: 6 name: back0 type: symmetry.  Reading zone data...done.
Zone: 7 name: front1 type: symmetry.  Reading zone data...done.
Zone: 8 name: front0 type: symmetry.  Reading zone data...done.
Zone: 9 name: plate1 type: wall.  Reading zone data...done.
Zone: 10 name: plate0 type: wall.  Reading zone data...done.
Zone: 11 name: top1 type: symmetry.  Reading zone data...done.
Zone: 12 name: top0 type: symmetry.  Reading zone data...done.
Zone: 13 name: outlet type: pressure-outlet.  Reading zone data...done.
Zone: 14 name: inlet type: velocity-inlet.  Reading zone data...done.
Zone: 16 name: default-interior type: interior.  Reading zone data...done.

FINISHED LEXING

Creating patch 0 for zone: 3 name: interface1 type: wall
Creating patch 1 for zone: 4 name: inteface0 type: wall
Creating patch 2 for zone: 5 name: back1 type: symmetry
Creating patch 3 for zone: 6 name: back0 type: symmetry
Creating patch 4 for zone: 7 name: front1 type: symmetry
Creating patch 5 for zone: 8 name: front0 type: symmetry
Creating patch 6 for zone: 9 name: plate1 type: wall
Creating patch 7 for zone: 10 name: plate0 type: wall
Creating patch 8 for zone: 11 name: top1 type: symmetry
Creating patch 9 for zone: 12 name: top0 type: symmetry
Creating patch 10 for zone: 13 name: outlet type: pressure-outlet
Creating patch 11 for zone: 14 name: inlet type: velocity-inlet
Creating cellZone 0 name: fluid type: fluid
Creating faceZone 0 name: default-interior type: interior
faceZone from Fluent indices: 34826 to: 941322 type: interior
patch 0 from Fluent indices: 0 to: 2989 type: wall
patch 1 from Fluent indices: 2990 to: 5979 type: wall
patch 2 from Fluent indices: 5980 to: 6853 type: symmetry
patch 3 from Fluent indices: 6854 to: 10717 type: symmetry
patch 4 from Fluent indices: 10718 to: 11591 type: symmetry
patch 5 from Fluent indices: 11592 to: 15455 type: symmetry
patch 6 from Fluent indices: 15456 to: 16690 type: wall
patch 7 from Fluent indices: 16691 to: 22150 type: wall
patch 8 from Fluent indices: 22151 to: 23385 type: symmetry
patch 9 from Fluent indices: 23386 to: 28845 type: symmetry
patch 10 from Fluent indices: 28846 to: 31835 type: pressure-outlet
patch 11 from Fluent indices: 31836 to: 34825 type: velocity-inlet

Writing mesh to "/home/moh3n/Desktop/920729/TwinLundWithCyclic2/constant/region0"

End
1st Step: fluent3DMeshToFoam for My 2nd Case(TwinLundWithCyclic3):
Code:
moh3n@moh3n:~/Desktop/920729/TwinLundWithCyclic2$ cd /home/moh3n/Desktop/920729/TwinLundWithCyclic3
moh3n@moh3n:~/Desktop/920729/TwinLundWithCyclic3$ fluent3DMeshToFoam TwinLundTwin.msh
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.1.1-221db2718bbb
Exec   : fluent3DMeshToFoam TwinLundTwin.msh
Date   : Oct 21 2013
Time   : 10:50:35
Host   : "moh3n"
PID    : 3395
Case   : /home/moh3n/Desktop/920729/TwinLundWithCyclic3
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

--> FOAM Warning : 
    From function dlOpen(const fileName&, const bool)
    in file POSIX.C at line 1175
    dlopen error : /home/moh3n/OpenFOAM/moh3n-2.1.1/platforms/linux64GccDPOpt/lib/libdynamicSmagorinskyModel.so: undefined symbol: _ZNK4Foam14incompressible9LESModels11GenEddyVisc1BEv
--> FOAM Warning : 
    From function dlLibraryTable::open(const fileName&, const bool)
    in file db/dynamicLibrary/dlLibraryTable/dlLibraryTable.C at line 96
    could not load "libdynamicSmagorinskyModel.so"
Dimension of grid: 3
Number of points: 325710
PointGroup: 1 start: 0 end: 325709.  Reading points...done.
Number of faces: 941323
FaceGroup: 3 start: 0 end: 2989.  Reading mixed faces...done.
FaceGroup: 4 start: 2990 end: 5979.  Reading mixed faces...done.
FaceGroup: 5 start: 5980 end: 6853.  Reading mixed faces...done.
FaceGroup: 6 start: 6854 end: 10717.  Reading mixed faces...done.
FaceGroup: 7 start: 10718 end: 11591.  Reading mixed faces...done.
FaceGroup: 8 start: 11592 end: 15455.  Reading mixed faces...done.
FaceGroup: 9 start: 15456 end: 16690.  Reading mixed faces...done.
FaceGroup: 10 start: 16691 end: 22150.  Reading mixed faces...done.
FaceGroup: 11 start: 22151 end: 23385.  Reading mixed faces...done.
FaceGroup: 12 start: 23386 end: 28845.  Reading mixed faces...done.
FaceGroup: 13 start: 28846 end: 31835.  Reading mixed faces...done.
FaceGroup: 14 start: 31836 end: 34825.  Reading mixed faces...done.
FaceGroup: 16 start: 34826 end: 941322.  Reading mixed faces...done.
Number of cells: 307970
CellGroup: 2 start: 0 end: 307969 type: 1
Zone: 2 name: fluid type: fluid.  Reading zone data...done.
Zone: 3 name: interface1 type: wall.  Reading zone data...done.
Zone: 4 name: interface0 type: wall.  Reading zone data...done.
Zone: 5 name: black1 type: symmetry.  Reading zone data...done.
Zone: 6 name: back0 type: symmetry.  Reading zone data...done.
Zone: 7 name: front1 type: symmetry.  Reading zone data...done.
Zone: 8 name: front0 type: symmetry.  Reading zone data...done.
Zone: 9 name: plate1 type: wall.  Reading zone data...done.
Zone: 10 name: plate0 type: wall.  Reading zone data...done.
Zone: 11 name: top1 type: symmetry.  Reading zone data...done.
Zone: 12 name: top0 type: symmetry.  Reading zone data...done.
Zone: 13 name: outlet type: pressure-outlet.  Reading zone data...done.
Zone: 14 name: inlet type: velocity-inlet.  Reading zone data...done.
Zone: 16 name: default-interior type: interior.  Reading zone data...done.

FINISHED LEXING

Creating patch 0 for zone: 3 name: interface1 type: wall
Creating patch 1 for zone: 4 name: interface0 type: wall
Creating patch 2 for zone: 5 name: black1 type: symmetry
Creating patch 3 for zone: 6 name: back0 type: symmetry
Creating patch 4 for zone: 7 name: front1 type: symmetry
Creating patch 5 for zone: 8 name: front0 type: symmetry
Creating patch 6 for zone: 9 name: plate1 type: wall
Creating patch 7 for zone: 10 name: plate0 type: wall
Creating patch 8 for zone: 11 name: top1 type: symmetry
Creating patch 9 for zone: 12 name: top0 type: symmetry
Creating patch 10 for zone: 13 name: outlet type: pressure-outlet
Creating patch 11 for zone: 14 name: inlet type: velocity-inlet
Creating cellZone 0 name: fluid type: fluid
Creating faceZone 0 name: default-interior type: interior
faceZone from Fluent indices: 34826 to: 941322 type: interior
patch 0 from Fluent indices: 0 to: 2989 type: wall
patch 1 from Fluent indices: 2990 to: 5979 type: wall
patch 2 from Fluent indices: 5980 to: 6853 type: symmetry
patch 3 from Fluent indices: 6854 to: 10717 type: symmetry
patch 4 from Fluent indices: 10718 to: 11591 type: symmetry
patch 5 from Fluent indices: 11592 to: 15455 type: symmetry
patch 6 from Fluent indices: 15456 to: 16690 type: wall
patch 7 from Fluent indices: 16691 to: 22150 type: wall
patch 8 from Fluent indices: 22151 to: 23385 type: symmetry
patch 9 from Fluent indices: 23386 to: 28845 type: symmetry
patch 10 from Fluent indices: 28846 to: 31835 type: pressure-outlet
patch 11 from Fluent indices: 31836 to: 34825 type: velocity-inlet

Writing mesh to "/home/moh3n/Desktop/920729/TwinLundWithCyclic3/constant/region0"

End
As you can see both geometries are the same (except for the name of the patches inteface0 & black1 that were misspelled and I changed them manually in polyMesh/boundary to interface0 & back1).

Then I did checkMesh for both meshes:

2nd Step: checkMesh -allGeometry for My 1st Case (TwinLundWithCyclic2):
Code:
moh3n@moh3n:~/Desktop/920729/TwinLundWithCyclic2$ checkMesh -allGeometry
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.1.1-221db2718bbb
Exec   : checkMesh -allGeometry
Date   : Oct 21 2013
Time   : 10:56:18
Host   : "moh3n"
PID    : 3420
Case   : /home/moh3n/Desktop/920729/TwinLundWithCyclic2
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

--> FOAM Warning : 
    From function dlOpen(const fileName&, const bool)
    in file POSIX.C at line 1175
    dlopen error : /home/moh3n/OpenFOAM/moh3n-2.1.1/platforms/linux64GccDPOpt/lib/libscaledMappedVelocityFixedValueKataoka.so: undefined symbol: _ZN4Foam12fvPatchFieldINS_6VectorIdEEE31patchMapperConstructorTablePtr_E
--> FOAM Warning : 
    From function dlLibraryTable::open(const fileName&, const bool)
    in file db/dynamicLibrary/dlLibraryTable/dlLibraryTable.C at line 96
    could not load "libscaledMappedVelocityFixedValueKataoka.so"
--> FOAM Warning : 
    From function dlOpen(const fileName&, const bool)
    in file POSIX.C at line 1175
    dlopen error : /home/moh3n/OpenFOAM/moh3n-2.1.1/platforms/linux64GccDPOpt/lib/libdynamicSmagorinskyModel.so: undefined symbol: _ZNK4Foam14incompressible9LESModels11GenEddyVisc1BEv
--> FOAM Warning : 
    From function dlLibraryTable::open(const fileName&, const bool)
    in file db/dynamicLibrary/dlLibraryTable/dlLibraryTable.C at line 96
    could not load "libdynamicSmagorinskyModel.so"
Create polyMesh for time = 0

Time = 0

Mesh stats
    points:           325710
    faces:            941323
    internal faces:   906497
    cells:            307970
    boundary patches: 12
    point zones:      0
    face zones:       1
    cell zones:       1

Overall number of cells of each type:
    hexahedra:     307970
    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: 2
    The mesh has multiple regions which are not connected by any face.
  <<Writing region information to "0/cellToRegion"

Checking patch topology for multiply connected surfaces ...
    Patch               Faces    Points   Surface topology                   Bounding box
    interface1          2990     3102     ok (non-closed singly connected)   (1.981 -2.2728694e-17 -1.296708e-17) (1.981 0.72 0.377)
    interface0          2990     3102     ok (non-closed singly connected)   (1.98 -2.2728694e-17 -1.296708e-17) (1.98 0.72 0.377)
    back1               874      940      ok (non-closed singly connected)   (1.981 0 -1.296708e-17) (2.41 0.72 4.3980468e-17)
    back0               3864     3995     ok (non-closed singly connected)   (1.1541915e-17 0 -2.7755576e-17) (1.98 0.72 4.3980468e-17)
    front1              874      940      ok (non-closed singly connected)   (1.981 0 0.377) (2.41 0.72 0.377)
    front0              3864     3995     ok (non-closed singly connected)   (-1.1541915e-17 0 0.377) (1.98 0.72 0.377)
    plate1              1235     1320     ok (non-closed singly connected)   (1.981 -2.2728694e-17 0) (2.41 2.2728694e-17 0.377)
    plate0              5460     5610     ok (non-closed singly connected)   (-1.1541915e-17 -2.2728694e-17 0) (1.98 2.2728694e-17 0.377)
    top1                1235     1320     ok (non-closed singly connected)   (1.981 0.72 0) (2.41 0.72 0.377)
    top0                5460     5610     ok (non-closed singly connected)   (-1.1541915e-17 0.72 -2.7755576e-17) (1.98 0.72 0.377)
    outlet              2990     3102     ok (non-closed singly connected)   (2.41 -2.2728694e-17 -1.296708e-17) (2.41 0.72 0.377)
    inlet               2990     3102     ok (non-closed singly connected)   (-1.1541915e-17 -2.2728694e-17 -2.7755576e-17) (1.1541915e-17 0.72 0.377)

Checking geometry...
    Overall domain bounding box (-1.1541915e-17 -2.2728694e-17 -2.7755576e-17) (2.41 0.72 0.377)
    Mesh (non-empty, non-wedge) directions (1 1 1)
    Mesh (non-empty) directions (1 1 1)
    Boundary openness (9.8982217e-17 -4.0994532e-15 -2.0078441e-15) OK.
    Max cell openness = 2.0697534e-16 OK.
    Max aspect ratio = 13.698473 OK.
    Minumum face area = 9.9802578e-06. Maximum face area = 0.0012946986.  Face area magnitudes OK.
    Min volume = 2.2534371e-07. Max volume = 7.5092518e-06.  Total volume = 0.65389896.  Cell volumes OK.
    Mesh non-orthogonality Max: 0 average: 0
    Non-orthogonality check OK.
    Face pyramids OK.
    Max skewness = 1.1094164e-08 OK.
    Coupled point location match (average 0) OK.
    Face tets OK.
    Min/max edge length = 0.0017207341 0.054926606 OK.
    All angles in faces OK.
    Face flatness (1 = flat, 0 = butterfly) : average = 1  min = 1
    All face flatness OK.
    Cell determinant (wellposedness) : minimum: 0.027406617 average: 1.9062195
    Cell determinant check OK.
    Concave cell check OK.

Mesh OK.

End
__________________
“If you have an apple and I have an apple and we exchange these apples then you and I will still each have one apple. But if you have an idea and I have an idea and we exchange these ideas, then each of us will have two ideas.”
cfdonline2mohsen is offline   Reply With Quote

Old   October 21, 2013, 09:24
Angry Continued
  #2
Senior Member
 
cfdonline2mohsen's Avatar
 
Mohsen KiaMansouri
Join Date: Jan 2010
Location: CFD Lab
Posts: 118
Rep Power: 7
cfdonline2mohsen is on a distinguished road
2nd Step: checkMesh -allGeometry for My 2nd Case (TwinLundWithCyclic3):
Code:
moh3n@moh3n:~/Desktop/920729/TwinLundWithCyclic3$ checkMesh -allGeometry
 /*---------------------------------------------------------------------------*\
 | =========                 |                                                 |
 | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
 |  \\    /   O peration     | Version:  2.1.1                                 |
 |   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
 |    \\/     M anipulation  |                                                 |
 \*---------------------------------------------------------------------------*/
 Build  : 2.1.1-221db2718bbb
 Exec   : checkMesh -allGeometry
 Date   : Oct 21 2013
 Time   : 10:56:23
 Host   : "moh3n"
 PID    : 3421
 Case   : /home/moh3n/Desktop/920729/TwinLundWithCyclic3
 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
 
 --> FOAM Warning : 
     From function dlOpen(const fileName&, const bool)
     in file POSIX.C at line 1175
     dlopen error :  /home/moh3n/OpenFOAM/moh3n-2.1.1/platforms/linux64GccDPOpt/lib/libscaledMappedVelocityFixedValueKataoka.so:  undefined symbol:  _ZN4Foam12fvPatchFieldINS_6VectorIdEEE31patchMapperConstructorTablePtr_E
 --> FOAM Warning : 
     From function dlLibraryTable::open(const fileName&, const bool)
     in file db/dynamicLibrary/dlLibraryTable/dlLibraryTable.C at line 96
     could not load "libscaledMappedVelocityFixedValueKataoka.so"
 --> FOAM Warning : 
     From function dlOpen(const fileName&, const bool)
     in file POSIX.C at line 1175
     dlopen error :  /home/moh3n/OpenFOAM/moh3n-2.1.1/platforms/linux64GccDPOpt/lib/libdynamicSmagorinskyModel.so:  undefined symbol: _ZNK4Foam14incompressible9LESModels11GenEddyVisc1BEv
 --> FOAM Warning : 
     From function dlLibraryTable::open(const fileName&, const bool)
     in file db/dynamicLibrary/dlLibraryTable/dlLibraryTable.C at line 96
     could not load "libdynamicSmagorinskyModel.so"
 Create polyMesh for time = 0
 
 Time = 0
 
 Mesh stats
     points:           325710
     faces:            941323
     internal faces:   906497
     cells:            307970
     boundary patches: 12
     point zones:      0
     face zones:       1
     cell zones:       1
 
 Overall number of cells of each type:
     hexahedra:     307970
     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: 2
     The mesh has multiple regions which are not connected by any face.
   <<Writing region information to "0/cellToRegion"
 
 Checking patch topology for multiply connected surfaces ...
     Patch               Faces    Points   Surface topology                   Bounding box
     interface1          2990     3102     ok (non-closed singly connected)   (1.981 -2.2728694e-17 -2.7755576e-17) (1.981 0.72 0.377)
     interface0          2990     3102     ok (non-closed singly connected)   (1.98 -2.2728694e-17 -2.7755576e-17) (1.98 0.72 0.377)
     back1               874      940      ok (non-closed singly connected)   (1.981 0 -2.7755576e-17) (2.41 0.72 4.3980468e-17)
     back0               3864     3995     ok (non-closed singly  connected)   (1.1541915e-17 0 -2.7755576e-17) (1.98 0.72 2.7755576e-17)
     front1              874      940      ok (non-closed singly connected)   (1.981 0 0.377) (2.41 0.72 0.377)
     front0              3864     3995     ok (non-closed singly connected)   (-1.1541915e-17 0 0.377) (1.98 0.72 0.377)
     plate1              1235     1320     ok (non-closed singly connected)   (1.981 -2.2728694e-17 0) (2.41 2.2728694e-17 0.377)
     plate0              5460     5610     ok (non-closed singly connected)   (-1.1541915e-17 -2.2728694e-17 2.7755576e-17) (1.98 2.2728694e-17 0.377)
     top1                1235     1320     ok (non-closed singly connected)   (1.981 0.72 -2.7755576e-17) (2.41 0.72 0.377)
     top0                5460     5610     ok (non-closed singly  connected)   (-1.1541915e-17 0.72 -2.7755576e-17) (1.98 0.72 0.377)
     outlet              2990     3102     ok (non-closed singly  connected)   (2.41 -2.2728694e-17 -1.296708e-17) (2.41 0.72 0.377)
     inlet               2990     3102     ok (non-closed singly  connected)   (-1.1541915e-17 -2.2728694e-17 -2.7755576e-17)  (1.1541915e-17 0.72 0.377)
 
 Checking geometry...
     Overall domain bounding box (-1.1541915e-17 -2.2728694e-17 -2.7755576e-17) (2.41 0.72 0.377)
     Mesh (non-empty, non-wedge) directions (1 1 1)
     Mesh (non-empty) directions (1 1 1)
     Boundary openness (6.0085425e-16 -3.2804356e-15 -2.0078441e-15) OK.
     Max cell openness = 2.0697534e-16 OK.
     Max aspect ratio = 13.698473 OK.
     Minumum face area = 9.9802578e-06. Maximum face area = 0.0012946986.  Face area magnitudes OK.
     Min volume = 2.2534371e-07. Max volume = 7.5092518e-06.  Total volume = 0.65389896.  Cell volumes OK.
     Mesh non-orthogonality Max: 0 average: 0
     Non-orthogonality check OK.
     Face pyramids OK.
     Max skewness = 1.1094164e-08 OK.
     Coupled point location match (average 0) OK.
     Face tets OK.
     Min/max edge length = 0.0017207341 0.054926606 OK.
     All angles in faces OK.
     Face flatness (1 = flat, 0 = butterfly) : average = 1  min = 1
     All face flatness OK.
     Cell determinant (wellposedness) : minimum: 0.027406617 average: 1.9062195
     Cell determinant check OK.
     Concave cell check OK.
 
 Mesh OK.
 
 End
As you see there are a little differences between the 2 Meshes That I highlighted them in red for interface1,interface0,back1,plate0,top1 patches and also for Boundary openness. I searched the whole Forums and I found that this e-17 or e-16 are the numerical way of saying 0. e.g. see #2.

My 1st suspicion is using e-17 or e-16 by OpenFOAM instead of 0
Then I used createPatch for creating 2 cyclic patches inside my domain:

3rd Step:
createPatch
for My 1st Case (TwinLundWithCyclic2):
Code:
moh3n@moh3n:~/Desktop/920729/TwinLundWithCyclic2$ createPatch
 /*---------------------------------------------------------------------------*\
 | =========                 |                                                 |
 | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
 |  \\    /   O peration     | Version:  2.1.1                                 |
 |   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
 |    \\/     M anipulation  |                                                 |
 \*---------------------------------------------------------------------------*/
 Build  : 2.1.1-221db2718bbb
 Exec   : createPatch
 Date   : Oct 21 2013
 Time   : 11:10:33
 Host   : "moh3n"
 PID    : 3449
 Case   : /home/moh3n/Desktop/920729/TwinLundWithCyclic2
 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
 
 --> FOAM Warning : 
     From function dlOpen(const fileName&, const bool)
     in file POSIX.C at line 1175
     dlopen error :  /home/moh3n/OpenFOAM/moh3n-2.1.1/platforms/linux64GccDPOpt/lib/libdynamicSmagorinskyModel.so:  undefined symbol: _ZNK4Foam14incompressible9LESModels11GenEddyVisc1BEv
 --> FOAM Warning : 
     From function dlLibraryTable::open(const fileName&, const bool)
     in file db/dynamicLibrary/dlLibraryTable/dlLibraryTable.C at line 96
     could not load "libdynamicSmagorinskyModel.so"
 Reading createPatchDict
 
 Create polyMesh for time = 0
 
 Adding new patch interface0New as patch 12 from 
 {
     type            cyclic;
     neighbourPatch  interface1New;
     matchTolerance  0.0001;
 }
 
 Adding new patch interface1New as patch 13 from 
 {
     type            cyclic;
     neighbourPatch  interface0New;
     matchTolerance  0.0001;
 }
 
 
 Moving faces from patch interface0 to patch 12
 Moving faces from patch interface1 to patch 13
 
 Doing topology modification to order faces.
 
 Dumping interface0New faces to "coupled_interface0New.obj"
 Dumping interface1New faces to "coupled_interface1New.obj"
 Dumping cyclic match as lines between face centres to "coupled_interface0Newinterface1New_match.obj"
 Not synchronising points.
 
 Removing patches with no faces in them.
 
 Removing empty patch interface1 at position 0
 Removing empty patch interface0 at position 1
 Removing patches.
 Dumping interface0New faces to "final_interface0New.obj"
 Dumping interface1New faces to "final_interface1New.obj"
 Dumping cyclic match as lines between face centres to "final_interface0Newinterface1New_match.obj"
 Writing repatched mesh to 1.5e-05
 
 End
3rd Step: createPatch for My 2nd Case (TwinLundWithCyclic3):
Code:
moh3n@moh3n:~/Desktop/920729/TwinLundWithCyclic3$ createPatch
 /*---------------------------------------------------------------------------*\
 | =========                 |                                                 |
 | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
 |  \\    /   O peration     | Version:  2.1.1                                 |
 |   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
 |    \\/     M anipulation  |                                                 |
 \*---------------------------------------------------------------------------*/
 Build  : 2.1.1-221db2718bbb
 Exec   : createPatch
 Date   : Oct 21 2013
 Time   : 11:10:49
 Host   : "moh3n"
 PID    : 3456
 Case   : /home/moh3n/Desktop/920729/TwinLundWithCyclic3
 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
 
 --> FOAM Warning : 
     From function dlOpen(const fileName&, const bool)
     in file POSIX.C at line 1175
     dlopen error :  /home/moh3n/OpenFOAM/moh3n-2.1.1/platforms/linux64GccDPOpt/lib/libdynamicSmagorinskyModel.so:  undefined symbol: _ZNK4Foam14incompressible9LESModels11GenEddyVisc1BEv
 --> FOAM Warning : 
     From function dlLibraryTable::open(const fileName&, const bool)
     in file db/dynamicLibrary/dlLibraryTable/dlLibraryTable.C at line 96
     could not load "libdynamicSmagorinskyModel.so"
 Reading createPatchDict
 
 Create polyMesh for time = 0
 
 Adding new patch interface0New as patch 12 from 
 {
     type            cyclic;
     neighbourPatch  interface1New;
     matchTolerance  0.0001;
 }
 
 Adding new patch interface1New as patch 13 from 
 {
     type            cyclic;
     neighbourPatch  interface0New;
     matchTolerance  0.0001;
 }
 
 
 Moving faces from patch interface0 to patch 12
 Moving faces from patch interface1 to patch 13
 
 Doing topology modification to order faces.
 
 Dumping interface0New faces to "coupled_interface0New.obj"
 Dumping interface1New faces to "coupled_interface1New.obj"
 Dumping cyclic match as lines between face centres to "coupled_interface0Newinterface1New_match.obj"
 Not synchronising points.
 
 Removing patches with no faces in them.
 
 Removing empty patch interface1 at position 0
 Removing empty patch interface0 at position 1
 Removing patches.
 Dumping interface0New faces to "final_interface0New.obj"
 Dumping interface1New faces to "final_interface1New.obj"
 Dumping cyclic match as lines between face centres to "final_interface0Newinterface1New_match.obj"
 Writing repatched mesh to 1.5e-05
 
 End
As you can see createPatch are the same for both meshes.
It seems that createPatch is not a suspect.

Then I replaced my old polyMesh with this new polyMesh that is created in 1.5e-05 and changed the names of the new cyclic patches interface0New & interface1New to interface0&interface1. the I changed the type of the inlet patch to mappedPatch.

__________________
“If you have an apple and I have an apple and we exchange these apples then you and I will still each have one apple. But if you have an idea and I have an idea and we exchange these ideas, then each of us will have two ideas.”
cfdonline2mohsen is offline   Reply With Quote

Old   October 21, 2013, 09:27
Default
  #3
Senior Member
 
cfdonline2mohsen's Avatar
 
Mohsen KiaMansouri
Join Date: Jan 2010
Location: CFD Lab
Posts: 118
Rep Power: 7
cfdonline2mohsen is on a distinguished road
Then I did checkMesh for my Final Meshes before running pimpleFoam:

4th Step: checkMesh -allGeometry for My 1st Case (TwinLundWithCyclic2):
Code:
moh3n@moh3n:~/Desktop/920729/TwinLundWithCyclic2$ checkMesh -allGeometry
  /*---------------------------------------------------------------------------*\
  | =========                 |                                                 |
  | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
  |  \\    /   O peration     | Version:  2.1.1                                 |
  |   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
  |    \\/     M anipulation  |                                                 |
  \*---------------------------------------------------------------------------*/
  Build  : 2.1.1-221db2718bbb
  Exec   : checkMesh -allGeometry
  Date   : Oct 21 2013
  Time   : 11:16:52
  Host   : "moh3n"
  PID    : 3544
  Case   : /home/moh3n/Desktop/920729/TwinLundWithCyclic2
  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
  
  --> FOAM Warning : 
      From function dlOpen(const fileName&, const bool)
      in file POSIX.C at line 1175
      dlopen error :   /home/moh3n/OpenFOAM/moh3n-2.1.1/platforms/linux64GccDPOpt/lib/libscaledMappedVelocityFixedValueKataoka.so:   undefined symbol:   _ZN4Foam12fvPatchFieldINS_6VectorIdEEE31patchMapperConstructorTablePtr_E
  --> FOAM Warning : 
      From function dlLibraryTable::open(const fileName&, const bool)
      in file db/dynamicLibrary/dlLibraryTable/dlLibraryTable.C at line 96
      could not load "libscaledMappedVelocityFixedValueKataoka.so"
  --> FOAM Warning : 
      From function dlOpen(const fileName&, const bool)
      in file POSIX.C at line 1175
      dlopen error :   /home/moh3n/OpenFOAM/moh3n-2.1.1/platforms/linux64GccDPOpt/lib/libdynamicSmagorinskyModel.so:   undefined symbol: _ZNK4Foam14incompressible9LESModels11GenEddyVisc1BEv
  --> FOAM Warning : 
      From function dlLibraryTable::open(const fileName&, const bool)
      in file db/dynamicLibrary/dlLibraryTable/dlLibraryTable.C at line 96
      could not load "libdynamicSmagorinskyModel.so"
  Create polyMesh for time = 0
  
  Time = 0
  
  Mesh stats
      points:           325710
      faces:            941323
      internal faces:   906497
      cells:            307970
      boundary patches: 12
      point zones:      0
      face zones:       1
      cell zones:       1
  
  Overall number of cells of each type:
      hexahedra:     307970
      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                   Bounding box
      back1               874      940      ok (non-closed singly connected)   (1.981 0 -1.296708e-17) (2.41 0.72 4.3980468e-17)
      back0               3864     3995     ok (non-closed singly   connected)   (1.1541915e-17 0 -2.7755576e-17) (1.98 0.72 4.3980468e-17)
      front1              874      940      ok (non-closed singly connected)   (1.981 0 0.377) (2.41 0.72 0.377)
      front0              3864     3995     ok (non-closed singly connected)   (-1.1541915e-17 0 0.377) (1.98 0.72 0.377)
      plate1              1235     1320     ok (non-closed singly  connected)   (1.981 -2.2728694e-17 0) (2.41 2.2728694e-17 0.377)
      plate0              5460     5610     ok (non-closed singly   connected)   (-1.1541915e-17 -2.2728694e-17 0) (1.98 2.2728694e-17   0.377)
      top1                1235     1320     ok (non-closed singly connected)   (1.981 0.72 0) (2.41 0.72 0.377)
      top0                5460     5610     ok (non-closed singly   connected)   (-1.1541915e-17 0.72 -2.7755576e-17) (1.98 0.72 0.377)
      outlet              2990     3102     ok (non-closed singly   connected)   (2.41 -2.2728694e-17 -1.296708e-17) (2.41 0.72 0.377)
      inlet               2990     3102     ok (non-closed singly   connected)   (-1.1541915e-17 -2.2728694e-17 -2.7755576e-17)   (1.1541915e-17 0.72 0.377)
      interface0          2990     3102     ok (non-closed singly   connected)   (1.98 -2.2728694e-17 -1.296708e-17) (1.98 0.72 0.377)
      interface1          2990     3102     ok (non-closed singly   connected)   (1.981 -2.2728694e-17 -1.296708e-17) (1.981 0.72 0.377)
  
  Checking geometry...
      Overall domain bounding box (-1.1541915e-17 -2.2728694e-17 -2.7755576e-17) (2.41 0.72 0.377)
      Mesh (non-empty, non-wedge) directions (1 1 1)
      Mesh (non-empty) directions (1 1 1)
      Boundary openness (-1.0283083e-16 -4.0994532e-15 -2.0078441e-15) OK.
      Max cell openness = 2.0697534e-16 OK.
      Max aspect ratio = 13.698473 OK.
      Minumum face area = 9.9802578e-06. Maximum face area = 0.0012946986.  Face area magnitudes OK.
      Min volume = 2.2534371e-07. Max volume = 7.5092518e-06.  Total volume = 0.65389896.  Cell volumes OK.
      Mesh non-orthogonality Max: 0 average: 0
      Non-orthogonality check OK.
      Face pyramids OK.
      Max skewness = 1.1094164e-08 OK.
      Coupled point location match (average 5.5511152e-16) OK.
      Face tets OK.
      Min/max edge length = 0.0017207341 0.054926606 OK.
      All angles in faces OK.
      Face flatness (1 = flat, 0 = butterfly) : average = 1  min = 1
      All face flatness OK.
      Cell determinant (wellposedness) : minimum: 0.027406617 average: 1.9062195
      Cell determinant check OK.
      Concave cell check OK.
  
  Mesh OK.
  
  End
4th Step: checkMesh -allGeometry for My 2nd Case (TwinLundWithCyclic3):
Code:
moh3n@moh3n:~/Desktop/920729/TwinLundWithCyclic3$ checkMesh -allGeometry
  /*---------------------------------------------------------------------------*\
  | =========                 |                                                 |
  | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
  |  \\    /   O peration     | Version:  2.1.1                                 |
  |   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
  |    \\/     M anipulation  |                                                 |
  \*---------------------------------------------------------------------------*/
  Build  : 2.1.1-221db2718bbb
  Exec   : checkMesh -allGeometry
  Date   : Oct 21 2013
  Time   : 11:17:26
  Host   : "moh3n"
  PID    : 3550
  Case   : /home/moh3n/Desktop/920729/TwinLundWithCyclic3
  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
  
  --> FOAM Warning : 
      From function dlOpen(const fileName&, const bool)
      in file POSIX.C at line 1175
      dlopen error :   /home/moh3n/OpenFOAM/moh3n-2.1.1/platforms/linux64GccDPOpt/lib/libscaledMappedVelocityFixedValueKataoka.so:   undefined symbol:   _ZN4Foam12fvPatchFieldINS_6VectorIdEEE31patchMapperConstructorTablePtr_E
  --> FOAM Warning : 
      From function dlLibraryTable::open(const fileName&, const bool)
      in file db/dynamicLibrary/dlLibraryTable/dlLibraryTable.C at line 96
      could not load "libscaledMappedVelocityFixedValueKataoka.so"
  --> FOAM Warning : 
      From function dlOpen(const fileName&, const bool)
      in file POSIX.C at line 1175
      dlopen error :   /home/moh3n/OpenFOAM/moh3n-2.1.1/platforms/linux64GccDPOpt/lib/libdynamicSmagorinskyModel.so:   undefined symbol: _ZNK4Foam14incompressible9LESModels11GenEddyVisc1BEv
  --> FOAM Warning : 
      From function dlLibraryTable::open(const fileName&, const bool)
      in file db/dynamicLibrary/dlLibraryTable/dlLibraryTable.C at line 96
      could not load "libdynamicSmagorinskyModel.so"
  Create polyMesh for time = 0
  
  Time = 0
  
  Mesh stats
      points:           325710
      faces:            941323
      internal faces:   906497
      cells:            307970
      boundary patches: 12
      point zones:      0
      face zones:       1
      cell zones:       1
  
  Overall number of cells of each type:
      hexahedra:     307970
      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                   Bounding box
      back1               874      940      ok (non-closed singly connected)   (1.981 0 -2.7755576e-17) (2.41 0.72 4.3980468e-17)
      back0               3864     3995     ok (non-closed singly   connected)   (1.1541915e-17 0 -2.7755576e-17) (1.98 0.72 2.7755576e-17)
      front1              874      940      ok (non-closed singly connected)   (1.981 0 0.377) (2.41 0.72 0.377)
      front0              3864     3995     ok (non-closed singly connected)   (-1.1541915e-17 0 0.377) (1.98 0.72 0.377)
      plate1              1235     1320     ok (non-closed singly  connected)   (1.981 -2.2728694e-17 0) (2.41 2.2728694e-17 0.377)
      plate0              5460     5610     ok (non-closed singly   connected)   (-1.1541915e-17 -2.2728694e-17 2.7755576e-17) (1.98   2.2728694e-17 0.377)
      top1                1235     1320     ok (non-closed singly connected)   (1.981 0.72 -2.7755576e-17) (2.41 0.72 0.377)
      top0                5460     5610     ok (non-closed singly   connected)   (-1.1541915e-17 0.72 -2.7755576e-17) (1.98 0.72 0.377)
      outlet              2990     3102     ok (non-closed singly   connected)   (2.41 -2.2728694e-17 -1.296708e-17) (2.41 0.72 0.377)
      inlet               2990     3102     ok (non-closed singly   connected)   (-1.1541915e-17 -2.2728694e-17 -2.7755576e-17)   (1.1541915e-17 0.72 0.377)
      interface0          2990     3102     ok (non-closed singly   connected)   (1.98 -2.2728694e-17 -2.7755576e-17) (1.98 0.72 0.377)
      interface1          2990     3102     ok (non-closed singly   connected)   (1.981 -2.2728694e-17 -2.7755576e-17) (1.981 0.72 0.377)
  
  Checking geometry...
      Overall domain bounding box (-1.1541915e-17 -2.2728694e-17 -2.7755576e-17) (2.41 0.72 0.377)
      Mesh (non-empty, non-wedge) directions (1 1 1)
      Mesh (non-empty) directions (1 1 1)
      Boundary openness (-4.7628464e-16 -3.2804356e-15 -2.0078441e-15) OK.
      Max cell openness = 2.0697534e-16 OK.
      Max aspect ratio = 13.698473 OK.
      Minumum face area = 9.9802578e-06. Maximum face area = 0.0012946986.  Face area magnitudes OK.
      Min volume = 2.2534371e-07. Max volume = 7.5092518e-06.  Total volume = 0.65389896.  Cell volumes OK.
      Mesh non-orthogonality Max: 0 average: 0
      Non-orthogonality check OK.
      Face pyramids OK.
      Max skewness = 1.1094164e-08 OK.
      Coupled point location match (average 5.0394958e-14) OK.
      Face tets OK.
      Min/max edge length = 0.0017207341 0.054926606 OK.
      All angles in faces OK.
      Face flatness (1 = flat, 0 = butterfly) : average = 1  min = 1
      All face flatness OK.
      Cell determinant (wellposedness) : minimum: 0.027406617 average: 1.9062195
      Cell determinant check OK.
      Concave cell check OK.
  
  Mesh OK.
  
  End
As you can see there are a little differences between the 2 Meshes for these this e-17 or e-16 values for interface1,interface0,back1,plate0,top1 patches just like the 2nd step and also for Boundary openness but with different values with respect to the 2nd step.
Also,a new difference has emerged in this step: Coupled point location match.

__________________
“If you have an apple and I have an apple and we exchange these apples then you and I will still each have one apple. But if you have an idea and I have an idea and we exchange these ideas, then each of us will have two ideas.”
cfdonline2mohsen is offline   Reply With Quote

Old   October 21, 2013, 09:28
Default
  #4
Senior Member
 
cfdonline2mohsen's Avatar
 
Mohsen KiaMansouri
Join Date: Jan 2010
Location: CFD Lab
Posts: 118
Rep Power: 7
cfdonline2mohsen is on a distinguished road
and the final step is running pimpleFoam:

5th Step: Running pimpleFoam for My 1st Case (TwinLundWithCyclic2):
Code:
moh3n@moh3n:~/Desktop/920729/TwinLundWithCyclic2$ pimpleFoam
  /*---------------------------------------------------------------------------*\
  | =========                 |                                                 |
  | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
  |  \\    /   O peration     | Version:  2.1.1                                 |
  |   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
  |    \\/     M anipulation  |                                                 |
  \*---------------------------------------------------------------------------*/
  Build  : 2.1.1-221db2718bbb
  Exec   : pimpleFoam
  Date   : Oct 21 2013
  Time   : 11:22:37
  Host   : "moh3n"
  PID    : 3554
  Case   : /home/moh3n/Desktop/920729/TwinLundWithCyclic2
  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
  
  Reading field U
  
  Reading/calculating face flux field phi
  
  Selecting incompressible transport model Newtonian
  Selecting turbulence model type LESModel
  Selecting LES turbulence model dynamicSmagorinsky
  dynamicSmagorinskyCoeffs
  {
      filter          simple;
      ce              1.048;
  }
  
  No field sources present
  
  
  PIMPLE: Operating solver in PISO mode
  
  
  Starting time loop
  
  Courant Number mean: 0.00064134496 max: 0.00066433568
  Time = 1.5e-05
  
  DILUPBiCG:  Solving for Ux, Initial residual = 1, Final residual = 5.4453587e-07, No Iterations 1
  DILUPBiCG:  Solving for Uy, Initial residual = 0.99448474, Final residual = 2.7934192e-08, No Iterations 2
  DILUPBiCG:  Solving for Uz, Initial residual = 0.98434205, Final residual = 5.1759288e-08, No Iterations 2
  DICPCG:  Solving for p, Initial residual = 1, Final residual = 0.041080613, No Iterations 184
  time step continuity errors : sum local = 1.9183968e-07, global = -3.7234788e-10, cumulative = -3.7234788e-10
  DICPCG:  Solving for p, Initial residual = 0.0036086778, Final residual = 9.0427381e-07, No Iterations 255
  time step continuity errors : sum local = 3.3297187e-10, global = -2.0421325e-12, cumulative = -3.7439002e-10
  bounding k, min: -3.4336254e-06 max: 0.93980124 average: 0.001494962
  ExecutionTime = 23.49 s  ClockTime = 24 s
  
  Courant Number mean: 0.00095867724 max: 0.0032413083
  Time = 3e-05
  
  DILUPBiCG:  Solving for Ux, Initial residual = 0.81314678, Final residual = 5.0166127e-06, No Iterations 1
  DILUPBiCG:  Solving for Uy, Initial residual = 0.31087511, Final residual = 2.1087205e-08, No Iterations 1
  DILUPBiCG:  Solving for Uz, Initial residual = 0.32819123, Final residual = 7.1596872e-07, No Iterations 1
  DICPCG:  Solving for p, Initial residual = 0.014583725, Final residual = 0.00072210423, No Iterations 193
  time step continuity errors : sum local = 2.9425828e-07, global = -5.9706547e-09, cumulative = -6.3450447e-09
  DICPCG:  Solving for p, Initial residual = 0.0022751552, Final residual = 9.142224e-07, No Iterations 240
  time step continuity errors : sum local = 3.0511796e-10, global = -8.1926347e-13, cumulative = -6.345864e-09
  bounding k, min: -0.00020715488 max: 0.88381996 average: 0.0024604238
  ExecutionTime = 43.1 s  ClockTime = 43 s
  
  Courant Number mean: 0.0012224251 max: 0.003201858
  Time = 4.5e-05
  
  DILUPBiCG:  Solving for Ux, Initial residual = 0.86187261, Final residual = 6.7931013e-06, No Iterations 1
  DILUPBiCG:  Solving for Uy, Initial residual = 0.45940243, Final residual = 2.8525143e-08, No Iterations 1
  DILUPBiCG:  Solving for Uz, Initial residual = 0.38062296, Final residual = 6.9428457e-07, No Iterations 1
And it works just fine for my 1st case both in Parallel & Serial.

5th Step: Running pimpleFoam for My 2nd Case (TwinLundWithCyclic3):
Code:
moh3n@moh3n:~/Desktop/920729/TwinLundWithCyclic3$ pimpleFoam
  /*---------------------------------------------------------------------------*\
  | =========                 |                                                 |
  | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
  |  \\    /   O peration     | Version:  2.1.1                                 |
  |   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
  |    \\/     M anipulation  |                                                 |
  \*---------------------------------------------------------------------------*/
  Build  : 2.1.1-221db2718bbb
  Exec   : pimpleFoam
  Date   : Oct 21 2013
  Time   : 11:22:54
  Host   : "moh3n"
  PID    : 3555
  Case   : /home/moh3n/Desktop/920729/TwinLundWithCyclic3
  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
  
  Reading field U
  
  Reading/calculating face flux field phi
  
  Selecting incompressible transport model Newtonian
  Selecting turbulence model type LESModel
  Selecting LES turbulence model dynamicSmagorinsky
  dynamicSmagorinskyCoeffs
  {
      filter          simple;
      ce              1.048;
  }
  
  No field sources present
  
  
  PIMPLE: Operating solver in PISO mode
  
  
  Starting time loop
  
  Courant Number mean: 0.00064134496 max: 0.00066433568
  Time = 1.5e-05
  
  #0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
  #1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
  #2   in "/lib/x86_64-linux-gnu/libc.so.6"
  #3  Foam::scaledMappedVelocityFixedValueFvPatchField::updateCoeffs() in  "/home/moh3n/OpenFOAM/moh3n-2.1.1/platforms/linux64GccDPOpt/lib/libscaledMappedVelocityFixedValueKataoka.so"
  #4  Foam::fvMatrix<Foam::Vector<double>  >::fvMatrix(Foam::GeometricField<Foam::Vector<double>,  Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensionSet  const&) in  "/opt/openfoam211/platforms/linux64GccDPOpt/bin/pimpleFoam"
  #5  Foam::fv::gaussLaplacianScheme<Foam::Vector<double>,  double>::fvmLaplacianUncorrected(Foam::GeometricField<double,  Foam::fvsPatchField, Foam::surfaceMesh> const&,  Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField,  Foam::volMesh> const&) in  "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
  #6  Foam::fv::gaussLaplacianScheme<Foam::Vector<double>,  double>::fvmLaplacian(Foam::GeometricField<double,  Foam::fvsPatchField, Foam::surfaceMesh> const&,  Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField,  Foam::volMesh> const&) in  "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
  #7  Foam::fv::laplacianScheme<Foam::Vector<double>,  double>::fvmLaplacian(Foam::GeometricField<double,  Foam::fvPatchField, Foam::volMesh> const&,  Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField,  Foam::volMesh> const&) in  "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
  #8  Foam::tmp<Foam::fvMatrix<Foam::Vector<double> > >  Foam::fvm::laplacian<Foam::Vector<double>,  double>(Foam::GeometricField<double, Foam::fvPatchField,  Foam::volMesh> const&,  Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField,  Foam::volMesh> const&, Foam::word const&) in  "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libincompressibleTurbulenceModel.so"
  #9  Foam::tmp<Foam::fvMatrix<Foam::Vector<double> > >  Foam::fvm::laplacian<Foam::Vector<double>,  double>(Foam::GeometricField<double, Foam::fvPatchField,  Foam::volMesh> const&,  Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField,  Foam::volMesh> const&) in  "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libincompressibleTurbulenceModel.so"
  #10   Foam::incompressible::LESModels::GenEddyVisc::divDevBeff(Foam::GeometricField<Foam::Vector<double>,  Foam::fvPatchField, Foam::volMesh>&) const in  "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libincompressibleLESModels.so"
  #11   Foam::incompressible::LESModel::divDevReff(Foam::GeometricField<Foam::Vector<double>,  Foam::fvPatchField, Foam::volMesh>&) const in  "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libincompressibleLESModels.so"
  #12  
   in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/pimpleFoam"
  #13  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
  #14  
   in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/pimpleFoam"
  Floating point exception
And it gave Floating point exception Error!!!

The 0, system and constant(except polyMesh) folders are the same for both cases. and the geometries are the same for both cases except little differences in checkMesh (using e-17 or e-16 by OpenFOAM instead of 0).
My 1st case runs just fine and my 2nd case gives Floating point exception error!!

I've also uploaded both cases:
http://www.mediafire.com/?m6gd1m06ro88bmj

You can do all the above steps yourself to see that I'm right!
It drove me crazy!
Some body please tell me what's the reason???
__________________
“If you have an apple and I have an apple and we exchange these apples then you and I will still each have one apple. But if you have an idea and I have an idea and we exchange these ideas, then each of us will have two ideas.”

Last edited by cfdonline2mohsen; October 21, 2013 at 14:09.
cfdonline2mohsen 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
pimpleFoam: turbulence->correct(); is not executed when using residualControl hfs OpenFOAM Running, Solving & CFD 3 October 29, 2013 09:35
Differences simpleFoam vs. pimpleFoam / RASModel.H vs turbulenceModel.H uli OpenFOAM Programming & Development 7 January 26, 2013 16:01
Mesh conversion problem (fluent3DMeshToFoam) Aadhavan OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 1 December 12, 2012 11:38
OpenFOAM command from inside MATLAB sega OpenFOAM Post-Processing 18 September 25, 2012 07:35
unv mesh corrupted after createPatch maddalena OpenFOAM Mesh Utilities 1 February 18, 2010 08:43


All times are GMT -4. The time now is 08:52.