CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   Possible Bug in pimpleFoam (or createPatch) (or fluent3DMeshToFoam) (http://www.cfd-online.com/Forums/openfoam/125200-possible-bug-pimplefoam-createpatch-fluent3dmeshtofoam.html)

cfdonline2mohsen October 21, 2013 09:21

Possible Bug in pimpleFoam (or createPatch) (or fluent3DMeshToFoam)
 
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??:confused::confused:

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


cfdonline2mohsen October 21, 2013 09:24

Continued
 
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.


cfdonline2mohsen October 21, 2013 09:27

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.


cfdonline2mohsen October 21, 2013 09:28

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!!!:mad::mad::mad:

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???


All times are GMT -4. The time now is 20:18.