|
[Sponsors] |
Possible Bug in pimpleFoam (or createPatch) (or fluent3DMeshToFoam) |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
October 21, 2013, 09:21 |
Possible Bug in pimpleFoam (or createPatch) (or fluent3DMeshToFoam)
|
#1 |
Senior Member
Mohsen KiaMansouri
Join Date: Jan 2010
Location: CFD Lab
Posts: 118
Rep Power: 16 |
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.” |
|
October 21, 2013, 09:24 |
Continued
|
#2 |
Senior Member
Mohsen KiaMansouri
Join Date: Jan 2010
Location: CFD Lab
Posts: 118
Rep Power: 16 |
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.” |
|
October 21, 2013, 09:27 |
|
#3 |
Senior Member
Mohsen KiaMansouri
Join Date: Jan 2010
Location: CFD Lab
Posts: 118
Rep Power: 16 |
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.” |
|
October 21, 2013, 09:28 |
|
#4 |
Senior Member
Mohsen KiaMansouri
Join Date: Jan 2010
Location: CFD Lab
Posts: 118
Rep Power: 16 |
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 14:09. |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[Commercial meshers] Mesh conversion problem (fluent3DMeshToFoam) | Aadhavan | OpenFOAM Meshing & Mesh Conversion | 2 | March 8, 2018 01:47 |
pimpleFoam: turbulence->correct(); is not executed when using residualControl | hfs | OpenFOAM Running, Solving & CFD | 3 | October 29, 2013 08:35 |
Differences simpleFoam vs. pimpleFoam / RASModel.H vs turbulenceModel.H | uli | OpenFOAM Programming & Development | 7 | January 26, 2013 15:01 |
OpenFOAM command from inside MATLAB | sega | OpenFOAM Post-Processing | 18 | September 25, 2012 07:35 |
[Salome] unv mesh corrupted after createPatch | maddalena | OpenFOAM Meshing & Mesh Conversion | 1 | February 18, 2010 07:43 |