|
[Sponsors] |
Possible Bug in pimpleFoam (or createPatch) (or fluent3DMeshToFoam) |
![]() |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
![]() |
![]() |
#1 |
Senior Member
Mohsen KiaMansouri
Join Date: Jan 2010
Location: CFD Lab
Posts: 118
Rep Power: 17 ![]() |
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 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 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.” ![]() |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Senior Member
Mohsen KiaMansouri
Join Date: Jan 2010
Location: CFD Lab
Posts: 118
Rep Power: 17 ![]() |
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 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 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 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.” ![]() |
|
![]() |
![]() |
![]() |
![]() |
#3 |
Senior Member
Mohsen KiaMansouri
Join Date: Jan 2010
Location: CFD Lab
Posts: 118
Rep Power: 17 ![]() |
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 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 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.” ![]() |
|
![]() |
![]() |
![]() |
![]() |
#4 |
Senior Member
Mohsen KiaMansouri
Join Date: Jan 2010
Location: CFD Lab
Posts: 118
Rep Power: 17 ![]() |
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 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 ![]() ![]() ![]() 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 15:09. |
|
![]() |
![]() |
![]() |
Thread Tools | Search this Thread |
Display Modes | |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
[Commercial meshers] Mesh conversion problem (fluent3DMeshToFoam) | Aadhavan | OpenFOAM Meshing & Mesh Conversion | 2 | March 8, 2018 02:47 |
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 |
OpenFOAM command from inside MATLAB | sega | OpenFOAM Post-Processing | 18 | September 25, 2012 08:35 |
[Salome] unv mesh corrupted after createPatch | maddalena | OpenFOAM Meshing & Mesh Conversion | 1 | February 18, 2010 08:43 |