
[Sponsors] 
February 12, 2014, 11:17 
simpleFoam error  "Floating point exception"

#1 
New Member
James
Join Date: Jun 2012
Posts: 8
Rep Power: 6 
Hi everyone,
I'm running into an error running simpleFoam (this problem also happens with icoFoam). Casefile can be found at the following address: https://drive.google.com/file/d/0B5R...it?usp=sharing I'm getting the error pasted below. The log file also reveals that the solution is diverging which I suspect why this is happening. checkMesh looks good, boundary conditions appear sensible. Any ideas would be much appreciated! #0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #2 in "/lib/x86_64linuxgnu/libc.so.6" #3 Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #4 at simpleFoam.C:0 #5 in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/simpleFoam" #6 __libc_start_main in "/lib/x86_64linuxgnu/libc.so.6" #7 in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/simpleFoam" Floating point exception (core dumped) 

February 12, 2014, 12:01 

#2 
Senior Member
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,438
Rep Power: 25 
Hi,
Obviously you've got diverging solution. Look at all those 1001 in log file, it is default maximum number of iteration of the matrix solver after which it will give up even it doesn't get the solution within given tolerances. I'm not quite get from the case files what you are trying to do but I've attached archive with the settings that can bear your setup. But the velocities in the narrow part of the mesh are around 4 km/s so I guess it doesn't make much sense. I've attached only archive of the system folder as I've made no changes in other folders. 

February 12, 2014, 13:53 

#3 
New Member
James
Join Date: Jun 2012
Posts: 8
Rep Power: 6 
Thanks alexeym!
It's running for me up to timestep 19.8 so obviously this is a big improvement. Basically, I'm just trying to measure permeability across this mesh by sampling pressure difference and flux through the two planes that are my inlets and outlets. Perhaps reducing the velocity bc at the inlet will reduce the absolute velocity in the narrow section but I imagine relatively speaking, the velocity would still be this high. Would reducing the initial velocity amount help the solution converge? 

February 12, 2014, 14:27 

#4 
Senior Member
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,438
Rep Power: 25 
Well,
You can try to estimate velocity in the narrow part from crosssection areas ratio. As you are trying to run simulation in laminar case, you should also estimate Reynolds number to see if it is really so. And answering your question: yes, it will help. 

March 5, 2014, 22:59 
error:time step continuty errors increasing and diverged

#5  
Member
Fengjiao Bian
Join Date: Nov 2013
Location: beijing
Posts: 30
Rep Power: 5 
Quote:
Code:
#0 Foam::error::printStack(Foam::Ostream&) in "/home/user5/OpenFOAM/OpenFOAM2.1.0/platforms/linux64Gcc45DPOpt/lib/libOpenFOAM.so" #1 Foam::sigFpe::sigHandler(int) in "/home/user5/OpenFOAM/OpenFOAM2.1.0/platforms/linux64Gcc45DPOpt/lib/libOpenFOAM.so" #2 __restore_rt at sigaction.c:0 #3 Foam::writeFuns::insert(double, Foam::DynamicList<float, 0u, 2u, 1u>&) in "/home/user5/OpenFOAM/OpenFOAM2.1.0/platforms/linux64Gcc45DPOpt/bin/foamToVTK" #4 void Foam::writeFuns::write<double>(std::ostream&, bool, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::pointPatchField, Foam::pointMesh> const&, Foam::vtkMesh const&) in "/home/user5/OpenFOAM/OpenFOAM2.1.0/platforms/linux64Gcc45DPOpt/bin/foamToVTK" #5 void Foam::internalWriter::write<double>(Foam::volPointInterpolation const&, Foam::PtrList<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) in "/home/user5/OpenFOAM/OpenFOAM2.1.0/platforms/linux64Gcc45DPOpt/bin/foamToVTK" #6 main in "/home/user5/OpenFOAM/OpenFOAM2.1.0/platforms/linux64Gcc45DPOpt/bin/foamToVTK" #7 __libc_start_main in "/lib64/libc.so.6" #8 Foam::regIOobject::writeObject(Foam::IOstream::streamFormat, Foam::IOstream::versionNumber, Foam::IOstream::compressionType) const in "/home/user5/OpenFOAM/OpenFOAM2.1.0/platforms/linux64Gcc45DPOpt/bin/foamToVTK" 浮点数例外 Code:
ddtSchemes { default steadyState; } gradSchemes { default cellLimited Gauss linear 1; grad(p) Gauss linear; grad(U) Gauss linear; } divSchemes { default Gauss upwind; div(phi,U) Gauss limitedLinearV 1;//Gauss upwind; div(phi,k) Gauss upwind; div(phi,epsilon) Gauss upwind;//linear; div(phi,nuTilda) Gauss linearUpwind grad(nuTilda); div((nuEff*dev(T(grad(U))))) Gauss linear; } laplacianSchemes { default none; laplacian(nuEff,U) Gauss linear corrected; laplacian((1A(U)),p) Gauss linear corrected; laplacian(DkEff,k) Gauss linear corrected; laplacian(DnuTildaEff,nuTilda) Gauss linear corrected; laplacian(1,p) Gauss linear corrected; laplacian(DepsilonEff,epsilon) Gauss linear corrected; } interpolationSchemes { default linear; interpolate(U) linear; } snGradSchemes { default corrected; } fluxRequired { default no; p ; } Code:
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 RAS turbulence model kEpsilon kEpsilonCoeffs { Cmu 0.09; C1 1.44; C2 1.92; sigmaEps 1.3; } No field sources present SIMPLE: convergence criteria field p tolerance 1e05 field U tolerance 1e05 field nuTilda tolerance 1e05 Starting time loop Time = 1 smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.0116175, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.0478811, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 0.0116854, No Iterations 2 GAMG: Solving for p, Initial residual = 1, Final residual = 0.0747252, No Iterations 1 time step continuity errors : sum local = 1222.25, global = 625.783, cumulative = 625.783 smoothSolver: Solving for epsilon, Initial residual = 0.999886, Final residual = 0.00670035, No Iterations 4 smoothSolver: Solving for k, Initial residual = 1, Final residual = 0.00694971, No Iterations 2 ExecutionTime = 39.72 s ClockTime = 40 s Time = 2 smoothSolver: Solving for Ux, Initial residual = 0.255939, Final residual = 0.00653363, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.207718, Final residual = 0.00447832, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.253055, Final residual = 0.00686111, No Iterations 2 GAMG: Solving for p, Initial residual = 1, Final residual = 0.16482, No Iterations 1 time step continuity errors : sum local = 1507.88, global = 682.725, cumulative = 1308.51 smoothSolver: Solving for epsilon, Initial residual = 0.36138, Final residual = 0.00191422, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.577132, Final residual = 0.00949441, No Iterations 2 ExecutionTime = 60.56 s ClockTime = 66 s Time = 3 smoothSolver: Solving for Ux, Initial residual = 0.199729, Final residual = 0.00633206, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.27506, Final residual = 0.00436556, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.200162, Final residual = 0.00666538, No Iterations 2 GAMG: Solving for p, Initial residual = 0.486066, Final residual = 0.0946856, No Iterations 1 time step continuity errors : sum local = 824.775, global = 340.359, cumulative = 1648.87 smoothSolver: Solving for epsilon, Initial residual = 0.298414, Final residual = 0.00162272, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.23253, Final residual = 0.00363581, No Iterations 2 ExecutionTime = 80.71 s ClockTime = 86 s Time = 4 smoothSolver: Solving for Ux, Initial residual = 0.206463, Final residual = 0.00682647, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.447125, Final residual = 0.00887106, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.206458, Final residual = 0.00696421, No Iterations 2 GAMG: Solving for p, Initial residual = 0.443974, Final residual = 0.0504603, No Iterations 1 time step continuity errors : sum local = 661.598, global = 177.288, cumulative = 1471.58 smoothSolver: Solving for epsilon, Initial residual = 0.432129, Final residual = 0.00656384, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.43187, Final residual = 0.00188758, No Iterations 2 ExecutionTime = 101.48 s ClockTime = 112 s Time = 5 smoothSolver: Solving for Ux, Initial residual = 0.180274, Final residual = 0.00628132, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.279662, Final residual = 0.00753944, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.181797, Final residual = 0.00637471, No Iterations 2 GAMG: Solving for p, Initial residual = 0.615593, Final residual = 0.0888842, No Iterations 1 time step continuity errors : sum local = 1278.21, global = 675.965, cumulative = 795.615 smoothSolver: Solving for epsilon, Initial residual = 0.310979, Final residual = 0.00371577, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.273199, Final residual = 0.00232058, No Iterations 2 ExecutionTime = 123.77 s ClockTime = 134 s Time = 6 smoothSolver: Solving for Ux, Initial residual = 0.134419, Final residual = 0.00424664, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.13053, Final residual = 0.00380063, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.135395, Final residual = 0.00435766, No Iterations 2 GAMG: Solving for p, Initial residual = 0.784403, Final residual = 0.141384, No Iterations 1 time step continuity errors : sum local = 1529.09, global = 741.707, cumulative = 53.9077 smoothSolver: Solving for epsilon, Initial residual = 0.231691, Final residual = 0.00282583, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.166224, Final residual = 0.00223461, No Iterations 2 ExecutionTime = 145.56 s ClockTime = 161 s Time = 7 smoothSolver: Solving for Ux, Initial residual = 0.136914, Final residual = 0.00399246, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.105305, Final residual = 0.00289856, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.13624, Final residual = 0.00405558, No Iterations 2 GAMG: Solving for p, Initial residual = 0.721444, Final residual = 0.143873, No Iterations 1 time step continuity errors : sum local = 1409.36, global = 557.136, cumulative = 503.228 smoothSolver: Solving for epsilon, Initial residual = 0.209118, Final residual = 0.0031759, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.156912, Final residual = 0.00272508, No Iterations 2 ExecutionTime = 167.51 s ClockTime = 183 s Time = 8 smoothSolver: Solving for Ux, Initial residual = 0.142662, Final residual = 0.00416516, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.106274, Final residual = 0.0027194, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.142317, Final residual = 0.00423871, No Iterations 2 GAMG: Solving for p, Initial residual = 0.50711, Final residual = 0.0923087, No Iterations 1 time step continuity errors : sum local = 1146.77, global = 448.292, cumulative = 951.52 smoothSolver: Solving for epsilon, Initial residual = 0.225749, Final residual = 0.00427787, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.174881, Final residual = 0.00347794, No Iterations 2 ExecutionTime = 189.73 s ClockTime = 210 s Time = 9 smoothSolver: Solving for Ux, Initial residual = 0.132193, Final residual = 0.00415394, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.0937998, Final residual = 0.00250987, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.133503, Final residual = 0.00419109, No Iterations 2 GAMG: Solving for p, Initial residual = 0.41403, Final residual = 0.0826546, No Iterations 1 time step continuity errors : sum local = 1231.54, global = 483.631, cumulative = 1435.15 smoothSolver: Solving for epsilon, Initial residual = 0.238165, Final residual = 0.00537919, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.18193, Final residual = 0.00424883, No Iterations 2 ExecutionTime = 211.34 s ClockTime = 232 s Time = 10 smoothSolver: Solving for Ux, Initial residual = 0.114206, Final residual = 0.00371409, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.0774936, Final residual = 0.00227964, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.116838, Final residual = 0.00371665, No Iterations 2 GAMG: Solving for p, Initial residual = 0.390742, Final residual = 0.0831381, No Iterations 1 time step continuity errors : sum local = 1287.78, global = 275.8, cumulative = 1710.95 smoothSolver: Solving for epsilon, Initial residual = 0.203471, Final residual = 0.00501127, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.169364, Final residual = 0.00457364, No Iterations 2 ExecutionTime = 233.55 s ClockTime = 259 s Time = 11 smoothSolver: Solving for Ux, Initial residual = 0.103788, Final residual = 0.00356374, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.0637353, Final residual = 0.00211766, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.104148, Final residual = 0.0034536, No Iterations 2 GAMG: Solving for p, Initial residual = 0.380292, Final residual = 0.0869304, No Iterations 1 time step continuity errors : sum local = 1302.78, global = 205.145, cumulative = 1916.1 smoothSolver: Solving for epsilon, Initial residual = 0.160339, Final residual = 0.00424108, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.14831, Final residual = 0.00447062, No Iterations 2 ExecutionTime = 255.03 s ClockTime = 281 s Time = 12 smoothSolver: Solving for Ux, Initial residual = 0.0950877, Final residual = 0.0033789, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.0547081, Final residual = 0.0018934, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.0942434, Final residual = 0.00320187, No Iterations 2 GAMG: Solving for p, Initial residual = 0.36055, Final residual = 0.0855903, No Iterations 1 time step continuity errors : sum local = 1272.49, global = 177.985, cumulative = 2094.08 smoothSolver: Solving for epsilon, Initial residual = 0.141591, Final residual = 0.0039718, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.139064, Final residual = 0.00437892, No Iterations 2 ExecutionTime = 276.63 s ClockTime = 307 s Time = 13 smoothSolver: Solving for Ux, Initial residual = 0.08845, Final residual = 0.00326364, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.0481885, Final residual = 0.00174469, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.0879515, Final residual = 0.00315651, No Iterations 2 GAMG: Solving for p, Initial residual = 0.346286, Final residual = 0.081073, No Iterations 1 time step continuity errors : sum local = 1200.51, global = 165.516, cumulative = 2259.6 smoothSolver: Solving for epsilon, Initial residual = 0.13731, Final residual = 0.00392148, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.13561, Final residual = 0.0043371, No Iterations 2 ExecutionTime = 296.61 s ClockTime = 327 s Time = 14 smoothSolver: Solving for Ux, Initial residual = 0.0826409, Final residual = 0.00302597, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.0436394, Final residual = 0.0016168, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.0823954, Final residual = 0.00302672, No Iterations 2 GAMG: Solving for p, Initial residual = 0.334584, Final residual = 0.0735457, No Iterations 1 time step continuity errors : sum local = 1086.94, global = 161.867, cumulative = 2421.46 smoothSolver: Solving for epsilon, Initial residual = 0.13495, Final residual = 0.00391054, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.134138, Final residual = 0.00437452, No Iterations 2 ExecutionTime = 316.77 s ClockTime = 354 s Time = 15 smoothSolver: Solving for Ux, Initial residual = 0.0779784, Final residual = 0.00287186, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.0403581, Final residual = 0.00151513, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.0772994, Final residual = 0.00287905, No Iterations 2 GAMG: Solving for p, Initial residual = 0.329316, Final residual = 0.0656906, No Iterations 1 time step continuity errors : sum local = 972.644, global = 160.998, cumulative = 2582.46 smoothSolver: Solving for epsilon, Initial residual = 0.134529, Final residual = 0.0038749, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.139388, Final residual = 0.00462282, No Iterations 2 ExecutionTime = 338.2 s ClockTime = 375 s Time = 16 smoothSolver: Solving for Ux, Initial residual = 0.075164, Final residual = 0.00285694, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.038141, Final residual = 0.00143983, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.0761046, Final residual = 0.00311105, No Iterations 2 GAMG: Solving for p, Initial residual = 0.325418, Final residual = 0.0581975, No Iterations 1 time step continuity errors : sum local = 870.478, global = 83.0109, cumulative = 2665.47 smoothSolver: Solving for epsilon, Initial residual = 0.439395, Final residual = 0.0130362, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.552394, Final residual = 0.0347422, No Iterations 2 ExecutionTime = 359.51 s ClockTime = 402 s Time = 17 smoothSolver: Solving for Ux, Initial residual = 0.464726, Final residual = 0.00644288, No Iterations 4 smoothSolver: Solving for Uy, Initial residual = 0.129872, Final residual = 0.00153488, No Iterations 4 smoothSolver: Solving for Uz, Initial residual = 0.773062, Final residual = 0.0125177, No Iterations 4 GAMG: Solving for p, Initial residual = 0.720437, Final residual = 0.12145, No Iterations 1 time step continuity errors : sum local = 4831.8, global = 1203.49, cumulative = 3868.96 smoothSolver: Solving for epsilon, Initial residual = 0.99852, Final residual = 0.0378969, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.99975, Final residual = 0.000130416, No Iterations 4 ExecutionTime = 381.33 s ClockTime = 424 s Time = 18 smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.00994666, No Iterations 4 smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.0101358, No Iterations 4 smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 0.0100816, No Iterations 4 GAMG: Solving for p, Initial residual = 1, Final residual = 0.28906, No Iterations 2 time step continuity errors : sum local = 1.85649e+11, global = 7.84046e+10, cumulative = 7.84047e+10 smoothSolver: Solving for epsilon, Initial residual = 1, Final residual = 0.0601755, No Iterations 2 smoothSolver: Solving for k, Initial residual = 1, Final residual = 0.0695155, No Iterations 2 ExecutionTime = 403.73 s ClockTime = 446 s Time = 19 smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.00701573, No Iterations 4 smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.00701791, No Iterations 4 smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 0.0070213, No Iterations 4 GAMG: Solving for p, Initial residual = 1, Final residual = 0.203189, No Iterations 2 time step continuity errors : sum local = 2.55152e+32, global = 1.23564e+32, cumulative = 1.23564e+32 smoothSolver: Solving for epsilon, Initial residual = 1, Final residual = 0.0624321, No Iterations 2 smoothSolver: Solving for k, Initial residual = 1, Final residual = 0.0115553, No Iterations 2 ExecutionTime = 426 s ClockTime = 468 s Time = 20 smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.00687755, No Iterations 4 smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.00687755, No Iterations 4 smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 0.00687755, No Iterations 4 GAMG: Solving for p, Initial residual = 1, Final residual = 0.246173, No Iterations 7 time step continuity errors : sum local = 3.63847e+82, global = 9.72572e+81, cumulative = 9.72572e+81 smoothSolver: Solving for epsilon, Initial residual = 1, Final residual = 0.00747124, No Iterations 8 smoothSolver: Solving for k, Initial residual = 1, Final residual = 5.54463e06, No Iterations 2 ExecutionTime = 454.86 s ClockTime = 497 s 

March 6, 2014, 01:05 

#6 
New Member
Bin Xu
Join Date: Apr 2012
Location: Singapore
Posts: 23
Rep Power: 6 
your first step : time step already diverge, you can have a check about your mesh and your boundary condition settings


March 7, 2014, 03:55 

#7 
Member
Fengjiao Bian
Join Date: Nov 2013
Location: beijing
Posts: 30
Rep Power: 5 
i run checkMesh and it shows:
Create polyMesh for time = 0 Time = 0 Mesh stats points: 644339 faces: 6657993 internal faces: 6626605 cells: 3261204 faces per cell: 4.07353 boundary patches: 4 point zones: 0 face zones: 1 cell zones: 2 Overall number of cells of each type: hexahedra: 100740 prisms: 3984 wedges: 0 pyramids: 34318 tet wedges: 0 tetrahedra: 3122162 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 blade 24146 21942 ok (closed singly connected) inlet 1841 966 ok (nonclosed singly connected) outlet 1841 966 ok (nonclosed singly connected) symmetry 3560 3649 ok (nonclosed singly connected) Checking geometry... Overall domain bounding box (12.7106 6.66098 12.7139) (12.7102 27.7695 12.7057) Mesh (nonempty, nonwedge) directions (1 1 1) Mesh (nonempty) directions (1 1 1) Boundary openness (2.3031e15 3.09845e16 5.48772e16) OK. Max cell openness = 1.07338e15 OK. Max aspect ratio = 30.0527 OK. Minimum face area = 6.8698e07. Maximum face area = 0.772395. Face area magnitudes OK. Min volume = 3.21357e10. Max volume = 0.154152. Total volume = 17466.7. Cell volumes OK. Mesh nonorthogonality Max: 70.2416 average: 19.945 *Number of severely nonorthogonal (> 70 degrees) faces: 3. Nonorthogonality check OK. <<Writing 3 nonorthogonal faces to set nonOrthoFaces Face pyramids OK. Max skewness = 3.90819 OK. Coupled point location match (average 0) OK. Mesh OK. how ablout the mesh quality? if bad,does the problem lie in the volume range or the 3 nonorthogonal faces or something eles? 

March 13, 2014, 21:09 

#8  
Member
Fengjiao Bian
Join Date: Nov 2013
Location: beijing
Posts: 30
Rep Power: 5 
Quote:
Code:
#0 Foam::error::printStack(Foam::Ostream&) in "/home/user5/OpenFOAM/OpenFOAM2.1.0/platforms/linux64Gcc45DPOpt/lib/libOpenFOAM.so" #1 Foam::sigFpe::sigHandler(int) in "/home/user5/OpenFOAM/OpenFOAM2.1.0/platforms/linux64Gcc45DPOpt/lib/libOpenFOAM.so" #2 __restore_rt at sigaction.c:0 #3 Foam::GAMGSolver::scalingFactor(Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double> const&, Foam::Field<double> const&) const in "/home/user5/OpenFOAM/OpenFOAM2.1.0/platforms/linux64Gcc45DPOpt/lib/libOpenFOAM.so" #4 Foam::GAMGSolver::scalingFactor(Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double>&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const in "/home/user5/OpenFOAM/OpenFOAM2.1.0/platforms/linux64Gcc45DPOpt/lib/libOpenFOAM.so" #5 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const in "/home/user5/OpenFOAM/OpenFOAM2.1.0/platforms/linux64Gcc45DPOpt/lib/libOpenFOAM.so" #6 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/user5/OpenFOAM/OpenFOAM2.1.0/platforms/linux64Gcc45DPOpt/lib/libOpenFOAM.so" #7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/home/user5/OpenFOAM/OpenFOAM2.1.0/platforms/linux64Gcc45DPOpt/lib/libfiniteVolume.so" #8 Foam::fvMatrix<double>::solve() in "/home/user5/OpenFOAM/OpenFOAM2.1.0/platforms/linux64Gcc45DPOpt/bin/MRFSimpleFoam" #9 main in "/home/user5/OpenFOAM/OpenFOAM2.1.0/platforms/linux64Gcc45DPOpt/bin/MRFSimpleFoam" #10 __libc_start_main in "/lib64/libc.so.6" #11 Foam::UOPstream::write(char) in "/home/user5/OpenFOAM/OpenFOAM2.1.0/platforms/linux64Gcc45DPOpt/bin/MRFSimpleFoam" 娴偣鏁颁緥澶 

March 13, 2014, 21:14 

#9 
Member
Fengjiao Bian
Join Date: Nov 2013
Location: beijing
Posts: 30
Rep Power: 5 
It seems the log file didnt upload, I am so confused. I can only cut the part where error happens.
Time = 75 smoothSolver: Solving for Ux, Initial residual = 0.0602954, Final residual = 0.00373723, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.0653118, Final residual = 0.00486974, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.0597028, Final residual = 0.00373956, No Iterations 2 GAMG: Solving for p, Initial residual = 9.91772e08, Final residual = 9.91772e08, No Iterations 0 time step continuity errors : sum local = 0.508126, global = 0.00403281, cumulative = 0.0725249 smoothSolver: Solving for epsilon, Initial residual = 0.111438, Final residual = 2.47159e07, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.0644941, Final residual = 0.00196909, No Iterations 2 ExecutionTime = 1101.44 s ClockTime = 1104 s Time = 76 smoothSolver: Solving for Ux, Initial residual = 0.0594781, Final residual = 0.00349262, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.060454, Final residual = 0.00433672, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.0589499, Final residual = 0.00349089, No Iterations 2 GAMG: Solving for p, Initial residual = 9.85788e08, Final residual = 9.85788e08, No Iterations 0 time step continuity errors : sum local = 0.536801, global = 0.00461372, cumulative = 0.0679112 smoothSolver: Solving for epsilon, Initial residual = 0.119377, Final residual = 1.54768e07, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.0767105, Final residual = 0.00221075, No Iterations 2 ExecutionTime = 1115.63 s ClockTime = 1118 s Time = 77 smoothSolver: Solving for Ux, Initial residual = 0.0592163, Final residual = 0.00329118, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.0563358, Final residual = 0.00385813, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.0588684, Final residual = 0.00329036, No Iterations 2 GAMG: Solving for p, Initial residual = 9.65561e08, Final residual = 9.65561e08, No Iterations 0 time step continuity errors : sum local = 0.565758, global = 0.00525311, cumulative = 0.0626581 smoothSolver: Solving for epsilon, Initial residual = 0.124077, Final residual = 9.12524e08, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.0872853, Final residual = 0.00261114, No Iterations 2 ExecutionTime = 1129.82 s ClockTime = 1133 s Time = 78 smoothSolver: Solving for Ux, Initial residual = 0.0601651, Final residual = 0.00311845, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.0534155, Final residual = 0.00342506, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.0595368, Final residual = 0.0031283, No Iterations 2 GAMG: Solving for p, Initial residual = 9.28703e08, Final residual = 9.28703e08, No Iterations 0 time step continuity errors : sum local = 0.594105, global = 0.00596336, cumulative = 0.0566947 smoothSolver: Solving for epsilon, Initial residual = 0.130123, Final residual = 5.31763e08, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.0996443, Final residual = 0.00317308, No Iterations 2 bounding k, min: 27.2023 max: 6.08394e+06 average: 47.9263 ExecutionTime = 1144.26 s ClockTime = 1147 s Time = 79 smoothSolver: Solving for Ux, Initial residual = 0.0614544, Final residual = 0.00297322, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.0514696, Final residual = 0.00306002, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.0607817, Final residual = 0.00298804, No Iterations 2 GAMG: Solving for p, Initial residual = 8.84393e08, Final residual = 8.84393e08, No Iterations 0 time step continuity errors : sum local = 0.620817, global = 0.00674595, cumulative = 0.0499488 smoothSolver: Solving for epsilon, Initial residual = 0.13737, Final residual = 2.93949e08, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.113938, Final residual = 0.00343359, No Iterations 2 bounding k, min: 216.755 max: 8.91232e+06 average: 61.9746 ExecutionTime = 1158.68 s ClockTime = 1161 s Time = 80 smoothSolver: Solving for Ux, Initial residual = 0.063004, Final residual = 0.00283814, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.0505971, Final residual = 0.00273304, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.0624865, Final residual = 0.0028546, No Iterations 2 GAMG: Solving for p, Initial residual = 8.29655e08, Final residual = 8.29655e08, No Iterations 0 time step continuity errors : sum local = 0.645584, global = 0.00760253, cumulative = 0.0423462 smoothSolver: Solving for epsilon, Initial residual = 0.143916, Final residual = 1.66331e08, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.127062, Final residual = 0.00356757, No Iterations 2 bounding k, min: 19.5544 max: 1.31047e+07 average: 81.9192 ExecutionTime = 1173.05 s ClockTime = 1176 s 

April 21, 2015, 13:35 
Floating point exception (core dumped) Immersed boundary simpleIbFoam

#10 
New Member
Anamitra Bagchi
Join Date: Mar 2015
Posts: 4
Rep Power: 3 
Hi I am using simpleIbFoam, each time I run it shows this error. The checkMesh looks ok. Can someone help with this issue ? Thanks in advance !!
anamitra@anamitraHPPaviliong6NotebookPC:~/IC_Files/Empty_Domain/simpleibfoam_vg$ simpleIbFoam /**\  =========    \\ / F ield  foamextend: Open Source CFD   \\ / O peration  Version: 3.1   \\ / A nd  Web: http://www.extendproject.de   \\/ M anipulation   \**/ Build : 3.17d8e040bf53d Exec : simpleIbFoam Date : Apr 21 2015 Time : 19:30:14 Host : anamitraHPPaviliong6NotebookPC PID : 8630 CtrlDict : /home/anamitra/foam/foamextend3.1/etc/controlDict Case : /home/anamitra/IC_Files/Empty_Domain/simpleibfoam_vg nProcs : 1 SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 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 RAS turbulence model kOmegaSST External flow Number of IB cells: 76 kOmegaSSTCoeffs { alphaK1 0.85034; alphaK2 1; alphaOmega1 0.5; alphaOmega2 0.85616; gamma1 0.5532; gamma2 0.4403; beta1 0.075; beta2 0.0828; betaStar 0.09; a1 0.31; c1 10; } Create immersed boundary cell mask Create immersed boundary face mask Found immersed boundary patch 0 named ibsimplevg Starting time loop Time = 1 DILUPBiCG: Solving for Ux, Initial residual = 0.357069, Final residual = 0.000767327, No Iterations 2 DILUPBiCG: Solving for Uy, Initial residual = 0.999927, Final residual = 0.00105008, No Iterations 2 DILUPBiCG: Solving for Uz, Initial residual = 0.285562, Final residual = 0.000466337, No Iterations 2 DICPCG: Solving for p, Initial residual = 0.984177, Final residual = 0.0909034, No Iterations 85 time step continuity errors : sum local = 0.175394, global = 0.00155391, cumulative = 0.00155391 Floating point exception (core dumped) anamitra@anamitraHPPaviliong6NotebookPC:~/IC_Files/Empty_Domain/simpleibfoam_vg$ 

April 21, 2015, 19:45 

#11 
New Member
N/A
Join Date: Jul 2010
Posts: 29
Rep Power: 8 
Anamitra:
It looks like you need to double check your BCs and ICs. Also, it's taking a lot of iterations to converge your pressure, you might consider a different numerical method if your BC/ICs are okay. 

June 25, 2015, 13:53 

#12 
New Member
Anamitra Bagchi
Join Date: Mar 2015
Posts: 4
Rep Power: 3 
Hi,
Thanks for the reply!! Could you suggest what changes are possible in the boundary conditions, I am still looking for it. I am not very sure what factors can cause the floating point exception. I am solving a case of external flow around a cuboid (with Immersed boundary) placed around the centre of the flow domain. Pressure looks like this: internalField nonuniform List<scalar> 256000 (long list found from mapping fields) boundaryField { ibsimplevg { type immersedBoundary; refValue uniform 0; refGradient uniform 0; fixesValue no; setDeadCellValue yes; deadCellValue 0; value nonuniform 0(); } inlet { type zeroGradient; } outlet { type zeroGradient; } sides { type symmetryPlane; } wall_plate { type zeroGradient; } top { type zeroGradient; } } 

August 4, 2015, 02:20 

#13 
New Member
Kittipos
Join Date: Feb 2015
Posts: 7
Rep Power: 3 
Hi Foamer
I have an error with running rhoPorousSimpleFoam in my simple model. Code:
#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam221/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam221/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #2 Uninterpreted: #3 Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() in "/opt/openfoam221/platforms/linuxGccDPOpt/lib/libfluidThermophysicalModels.so" #4 Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() in "/opt/openfoam221/platforms/linuxGccDPOpt/lib/libfluidThermophysicalModels.so" #5 in "/opt/openfoam221/platforms/linuxGccDPOpt/bin/rhoPorousSimpleFoam" #6 __libc_start_main in "/lib/i386linuxgnu/libc.so.6" #7 in "/opt/openfoam221/platforms/linuxGccDPOpt/bin/rhoPorousSimpleFoam" Floating point exception (core dumped) 

Tags 
error, floating point exception, openfoam, simplefoam 
Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
problem in toposet  Ahmed Khattab  OpenFOAM PreProcessing  11  October 29, 2015 09:26 
ActuatorDiskExplicitForce in OF2.1. Help  be_inspired  OpenFOAM Programming & Development  8  July 3, 2014 10:54 
3Dmixer on simpleFoam and MRFInterFoam.  sharonyue  OpenFOAM Running, Solving & CFD  1  May 12, 2013 01:04 
error EOF in blockMesh  Ahmed Khattab  OpenFOAM Meshing & Mesh Conversion  7  May 17, 2012 00:37 
CFX Post: Problems with moving point cloud for changing time steps  spatialtime  CFX  0  December 7, 2009 05:56 