
[Sponsors] 
December 15, 2009, 21:05 
Floating point exception error

#1 
New Member
Alan Harrland
Join Date: Mar 2009
Posts: 21
Rep Power: 10 
Hi,
I am trying to run a scenario of a three dimensional wedge in supersonic flow, using rhoCentralFoam. The solution will iterate for a while, but then it will crash and I get the following error: #0 Foam::error:rintStack(Foam::Ostream&) in "/home/cdoolan/OpenFOAM/OpenFOAM1.6/lib/linuxGccDPOpt/libOpenFOAM.so" #1 Foam::sigFpe::sigFpeHandler(int) in "/home/cdoolan/OpenFOAM/OpenFOAM1.6/lib/linuxGccDPOpt/libOpenFOAM.so" #2 Uninterpreted: #3 Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&) in "/home/cdoolan/OpenFOAM/OpenFOAM1.6/lib/linuxGccDPOpt/libOpenFOAM.so" #4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&) in "/home/cdoolan/OpenFOAM/OpenFOAM1.6/applications/bin/linuxGccDPOpt/rhoCentralFoam" #5 main in "/home/cdoolan/OpenFOAM/OpenFOAM1.6/applications/bin/linuxGccDPOpt/rhoCentralFoam" #6 __libc_start_main in "/lib/libc.so.6" #7 _start at /usr/src/packages/BUILD/glibc2.9/csu/../sysdeps/i386/elf/start.S:122 Floating point exception I have tried turning off the floating point exception catching, which allows the solver to iterate until the end time, but I cannot view in post processing as the solution space contains NaN values. I am building my mesh in Gmsh, and importing into OpenFOAM (with which I have no problems). My initial condition files are below: p: Code:
/** C++ **\  =========    \\ / F ield  OpenFOAM: The Open Source CFD Toolbox   \\ / O peration  Version: 1.6   \\ / A nd  Web: http://www.OpenFOAM.org   \\/ M anipulation   \**/ FoamFile { version 2.0; format ascii; class volScalarField; object p; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [1 1 2 0 0 0 0]; internalField uniform 1; boundaryField { inlet { type fixedValue; value uniform 1; } outlet { type zeroGradient; } obstacle { type zeroGradient; } defaultFaces { type empty; } } // ************************************************************************* // Code:
/** C++ **\  =========    \\ / F ield  OpenFOAM: The Open Source CFD Toolbox   \\ / O peration  Version: 1.6   \\ / A nd  Web: http://www.OpenFOAM.org   \\/ M anipulation   \**/ FoamFile { version 2.0; format ascii; class volScalarField; object T; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 0 0 1 0 0 0]; internalField uniform 1; boundaryField { inlet { type fixedValue; value uniform 1; } outlet { type zeroGradient; } obstacle { type zeroGradient; } defaultFaces { type empty; } } // ************************************************************************* // Code:
/** C++ **\  =========    \\ / F ield  OpenFOAM: The Open Source CFD Toolbox   \\ / O peration  Version: 1.6   \\ / A nd  Web: http://www.OpenFOAM.org   \\/ M anipulation   \**/ FoamFile { version 2.0; format ascii; class volVectorField; object U; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 1 1 0 0 0 0]; internalField uniform (0 0 0); boundaryField { inlet { type fixedValue; value uniform (5 0 0); } outlet { type zeroGradient; } obstacle { type slip; } defaultFaces { type empty; } } // ************************************************************************* // Code:
/** C++ **\  =========    \\ / F ield  OpenFOAM: The Open Source CFD Toolbox   \\ / O peration  Version: 1.6   \\ / A nd  Web: www.OpenFOAM.org   \\/ M anipulation   \**/ FoamFile { version 2.0; format ascii; class polyBoundaryMesh; location "constant/polyMesh"; object boundary; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 4 ( inlet { type patch; nFaces 156; startFace 14250; } outlet { type patch; nFaces 894; startFace 14406; } obstacle { type patch; nFaces 114; startFace 15300; } defaultFaces { type empty; nFaces 0; startFace 15414; } ) // ************************************************************************* // Thanks, any help would be much appreciated, as I am pulling my hair out here, trying to figure this out. Thanks, Alan 

December 26, 2009, 01:38 

#2 
New Member
Alan Harrland
Join Date: Mar 2009
Posts: 21
Rep Power: 10 
Is there anyone here willing to have a look at this for me? I can't for the life of me figure it out. It would be greatly appreciated as I am pulling my hair out.
Thanks in advance, Alan 

December 27, 2009, 08:10 

#3 
New Member
Kerstin Heinen
Join Date: Mar 2009
Location: Ludwigshafen, Germany
Posts: 27
Rep Power: 10 
Hi Alan,
I would say even as advanced cfd user it is nearly impossible to give the answer to your problem... Just some advice how I would proceed if I were you: Start writing out results very often and look at intermediate results. If crashing is due to bad mesh quality you may see it in the intermediate results , that some single cells get large values which are then later "exploding"  check your direction vector for velocity in the boundary condition, does it really match to the geometry?  how do you initialise the solution? If this is a turbulent flow, don't start from scratch...You can run a "potentialFoam" before or reduce your velocity magnitude, to have less severe starting conditions for the solver...  make a "checkMesh"...with errors in the output there my experience is, that you will not have a stable run in most cases... Good luck for this "hunt". It wouldn't be CFD, if it would never happen that one has exactly the same problems, you have now. Kerstin 

December 27, 2009, 11:32 

#4 
Senior Member

Maybe you could try to modify the pressure initial conditions to
Code:
inlet { type zeroGradient; } outlet { type fixedValue; value uniform 0; } 

December 27, 2009, 11:49 

#5 
Senior Member

I also have a FPE error that keeps coming back and I can't figure out what initiates it...
Code:
Initializing the GGI interpolator between master/shadow patches: insideSlider/outsideSlider Evaluation of GGI weighting factors: Largest slave weighting factor correction : 0.000117224699752 average: 0.000117210464361 Largest master weighting factor correction: 4.4408920985e16 average: 1.09561482693e16 smoothSolver: Solving for Ux, Initial residual = 4.18366053939e05, Final residual = 3.99534809799e11, No Iterations 5 smoothSolver: Solving for Uy, Initial residual = 0.000183494576854, Final residual = 1.58095734727e10, No Iterations 5 PCG: Solving for p, Initial residual = 0.00115761064277, Final residual = 9.69845796082e07, No Iterations 285 PCG: Solving for p, Initial residual = 9.99140357817e05, Final residual = 9.84892398469e07, No Iterations 16 time step continuity errors : sum local = 1.61074737955e11, global = 1.38664307317e13, cumulative = 2.44455566515e08 PCG: Solving for p, Initial residual = 2.7137207064e05, Final residual = 9.65820558314e07, No Iterations 41 PCG: Solving for p, Initial residual = 1.15863553234e05, Final residual = 9.31153715832e07, No Iterations 3 time step continuity errors : sum local = 1.5228622047e11, global = 5.55997394556e14, cumulative = 2.44456122512e08 PBiCG: Solving for omega, Initial residual = 1.98154832644e05, Final residual = 3.46778670816e08, No Iterations 1 PBiCG: Solving for k, Initial residual = 8.33751071488e05, Final residual = 2.53152769991e07, No Iterations 1 ExecutionTime = 16216.54 s ClockTime = 16307 s Courant Number mean: 0.0394310396556 max: 0.499997648141 velocity magnitude: 124.35829345 deltaT = 9.57280153069e06 Time = 0.316468 Initializing the GGI interpolator between master/shadow patches: insideSlider/outsideSlider Evaluation of GGI weighting factors: Largest slave weighting factor correction : 2.36652718848e05 average: 2.36504740638e05 Largest master weighting factor correction: 4.4408920985e16 average: 1.09561482693e16 smoothSolver: Solving for Ux, Initial residual = 4.18646624848e05, Final residual = 4.00745405967e11, No Iterations 5 smoothSolver: Solving for Uy, Initial residual = 0.000183554079661, Final residual = 1.58722402612e10, No Iterations 5 PCG: Solving for p, Initial residual = 0.00113267387872, Final residual = 9.99882920123e07, No Iterations 149 PCG: Solving for p, Initial residual = 9.82117938326e05, Final residual = 9.95912006942e07, No Iterations 15 time step continuity errors : sum local = 1.62857195533e11, global = 1.21882881021e12, cumulative = 2.44443934224e08 PCG: Solving for p, Initial residual = 2.63425356134e05, Final residual = 9.5916336339e07, No Iterations 39 PCG: Solving for p, Initial residual = 1.13546165489e05, Final residual = 9.1377403705e07, No Iterations 3 time step continuity errors : sum local = 1.4942541397e11, global = 1.21024107242e12, cumulative = 2.44431831814e08 PBiCG: Solving for omega, Initial residual = 1.97786652847e05, Final residual = 3.4608657236e08, No Iterations 1 PBiCG: Solving for k, Initial residual = 8.33789411068e05, Final residual = 2.53398294466e07, No Iterations 1 bounding k, min: 0.122378016374 max: 237.949650915 average: 11.0237179132 ExecutionTime = 16216.96 s ClockTime = 16308 s Courant Number mean: 0.0394313290225 max: 0.499998490876 velocity magnitude: 124.35851925 deltaT = 9.57280153069e06 Time = 0.316477 Initializing the GGI interpolator between master/shadow patches: insideSlider/outsideSlider Evaluation of GGI weighting factors: Largest slave weighting factor correction : 0.000160794100478 average: 0.000160779324961 Largest master weighting factor correction: 4.4408920985e16 average: 1.22708860616e16 smoothSolver: Solving for Ux, Initial residual = 4.18648802413e05, Final residual = 4.01555017906e11, No Iterations 5 smoothSolver: Solving for Uy, Initial residual = 0.000183664619068, Final residual = 1.5925205909e10, No Iterations 5 PCG: Solving for p, Initial residual = 0.0012822507928, Final residual = 9.74487717197e07, No Iterations 314 PCG: Solving for p, Initial residual = 0.000183503994712, Final residual = 9.54017668329e07, No Iterations 27 time step continuity errors : sum local = 1.56040920447e11, global = 2.71683168776e14, cumulative = 2.4443156013e08 PCG: Solving for p, Initial residual = 4.54649995932e05, Final residual = 9.7803002942e07, No Iterations 92 PCG: Solving for p, Initial residual = 1.87289344062e05, Final residual = 8.32829759138e07, No Iterations 5 time step continuity errors : sum local = 1.36217381462e11, global = 7.85213463584e13, cumulative = 2.44439412265e08 PBiCG: Solving for omega, Initial residual = 1.97895544736e05, Final residual = 3.43292848071e08, No Iterations 1 PBiCG: Solving for k, Initial residual = 8.33872916715e05, Final residual = 2.5269138733e07, No Iterations 1 ExecutionTime = 16217.56 s ClockTime = 16308 s Courant Number mean: 0.0394316251525 max: 0.499999326467 velocity magnitude: 124.358741937 deltaT = 9.57280153069e06 Time = 0.316487 Initializing the GGI interpolator between master/shadow patches: insideSlider/outsideSlider Evaluation of GGI weighting factors: Largest slave weighting factor correction : 0.000285528627237 average: 0.000285513904142 Largest master weighting factor correction: 4.4408920985e16 average: 9.49532850008e17 smoothSolver: Solving for Ux, Initial residual = 4.18150647223e05, Final residual = 4.02129097578e11, No Iterations 5 smoothSolver: Solving for Uy, Initial residual = 0.000183919043323, Final residual = 1.59663449673e10, No Iterations 5 PCG: Solving for p, Initial residual = 0.00118091660555, Final residual = 9.97058782232e07, No Iterations 280 PCG: Solving for p, Initial residual = 0.000103443967517, Final residual = 8.8633487677e07, No Iterations 12 time step continuity errors : sum local = 1.4493480812e11, global = 2.1001888378e13, cumulative = 2.44437312076e08 PCG: Solving for p, Initial residual = 2.80136398208e05, Final residual = 9.62893451566e07, No Iterations 42 PCG: Solving for p, Initial residual = 1.11499919963e05, Final residual = 9.28004174336e07, No Iterations 3 time step continuity errors : sum local = 1.51747798538e11, global = 4.16825424229e13, cumulative = 2.44433143822e08 PBiCG: Solving for omega, Initial residual = 1.98065657161e05, Final residual = 3.40451890596e08, No Iterations 1 PBiCG: Solving for k, Initial residual = 8.33920881206e05, Final residual = 2.52098556426e07, No Iterations 1 ExecutionTime = 16218.13 s ClockTime = 16309 s Courant Number mean: 0.0394319159774 max: 0.500000147171 velocity magnitude: 124.358960763 deltaT = 9.57170436146e06 Time = 0.316496 Initializing the GGI interpolator between master/shadow patches: insideSlider/outsideSlider Evaluation of GGI weighting factors: Largest slave weighting factor correction : 0.000397873642995 average: 0.000397859002476 Largest master weighting factor correction: 4.4408920985e16 average: 1.18326401309e16 smoothSolver: Solving for Ux, Initial residual = 4.17820269512e05, Final residual = 4.02754195966e11, No Iterations 5 smoothSolver: Solving for Uy, Initial residual = 0.000183916215376, Final residual = 1.60180223308e10, No Iterations 5 PCG: Solving for p, Initial residual = 0.00116209995737, Final residual = 9.91222118143e07, No Iterations 96 PCG: Solving for p, Initial residual = 9.7432771035e05, Final residual = 8.65531722367e07, No Iterations 12 time step continuity errors : sum local = 1.41496821471e11, global = 1.1742061533e12, cumulative = 2.44444885883e08 PCG: Solving for p, Initial residual = 2.64247799783e05, Final residual = 9.26829060343e07, No Iterations 41 PCG: Solving for p, Initial residual = 1.09427368381e05, Final residual = 9.06631309017e07, No Iterations 3 time step continuity errors : sum local = 1.48215256106e11, global = 1.20279142129e12, cumulative = 2.44456913798e08 PBiCG: Solving for omega, Initial residual = 1.98363312799e05, Final residual = 3.4374674519e08, No Iterations 1 PBiCG: Solving for k, Initial residual = 8.33873645469e05, Final residual = 2.51784635071e07, No Iterations 1 ExecutionTime = 16218.51 s ClockTime = 16309 s Courant Number mean: 0.0394276758498 max: 0.499943657681 velocity magnitude: 124.359179205 deltaT = 9.57170436146e06 Time = 0.316506 Initializing the GGI interpolator between master/shadow patches: insideSlider/outsideSlider Evaluation of GGI weighting factors: Largest slave weighting factor correction : 0.000497859420717 average: 0.000497844892901 Largest master weighting factor correction: 4.4408920985e16 average: 8.47275466161e17 smoothSolver: Solving for Ux, Initial residual = 4.17905810822e05, Final residual = 4.03898968375e11, No Iterations 5 smoothSolver: Solving for Uy, Initial residual = 0.000184025898617, Final residual = 1.60739458371e10, No Iterations 5 PCG: Solving for p, Initial residual = 0.00112936194147, Final residual = 9.74360513696e07, No Iterations 274 PCG: Solving for p, Initial residual = 9.76900241161e05, Final residual = 9.7094946801e07, No Iterations 11 time step continuity errors : sum local = 1.58701620076e11, global = 1.80901108399e13, cumulative = 2.44455104787e08 PCG: Solving for p, Initial residual = 2.57632315604e05, Final residual = 9.37704584393e07, No Iterations 41 PCG: Solving for p, Initial residual = 1.09362657999e05, Final residual = 8.88228716267e07, No Iterations 3 time step continuity errors : sum local = 1.45180616517e11, global = 3.24254714623e13, cumulative = 2.44451862239e08 PBiCG: Solving for omega, Initial residual = 1.9899688921e05, Final residual = 3.67212803247e08, No Iterations 1 PBiCG: Solving for k, Initial residual = 8.33898572745e05, Final residual = 2.51738401993e07, No Iterations 1 bounding k, min: 0.802415643715 max: 237.717244127 average: 11.0287760469 ExecutionTime = 16219.02 s ClockTime = 16310 s Courant Number mean: 0.0394279428752 max: 0.499944485086 velocity magnitude: 124.359399357 deltaT = 9.57170436146e06 Time = 0.316516 Initializing the GGI interpolator between master/shadow patches: insideSlider/outsideSlider Evaluation of GGI weighting factors: Largest slave weighting factor correction : 0.000585499543768 average: 0.000585485158807 Largest master weighting factor correction: 4.4408920985e16 average: 9.64141047701e17 smoothSolver: Solving for Ux, Initial residual = 4.17858218617e05, Final residual = 4.05074168974e11, No Iterations 5 smoothSolver: Solving for Uy, Initial residual = 0.000184094269779, Final residual = 1.61390224942e10, No Iterations 5 PCG: Solving for p, Initial residual = 0.00112989517074, Final residual = 9.91499532231e07, No Iterations 152 PCG: Solving for p, Initial residual = 9.69060456898e05, Final residual = 9.08043230719e07, No Iterations 12 time step continuity errors : sum local = 1.48405476369e11, global = 9.5042254985e13, cumulative = 2.44461366465e08 PCG: Solving for p, Initial residual = 2.59766736557e05, Final residual = 9.42002427652e07, No Iterations 41 PCG: Solving for p, Initial residual = 1.1185836139e05, Final residual = 9.04423006686e07, No Iterations 3 time step continuity errors : sum local = 1.47813793588e11, global = 1.04539352081e12, cumulative = 2.444718204e08 PBiCG: Solving for omega, Initial residual = 1.98655390897e05, Final residual = 3.54224756962e08, No Iterations 1 PBiCG: Solving for k, Initial residual = 8.3438685774e05, Final residual = 2.51623930413e07, No Iterations 1 ExecutionTime = 16219.45 s ClockTime = 16310 s Courant Number mean: 0.0394282050904 max: 0.499945308896 velocity magnitude: 124.35961817 deltaT = 9.57170436146e06 Time = 0.316525 Initializing the GGI interpolator between master/shadow patches: insideSlider/outsideSlider Evaluation of GGI weighting factors: Largest slave weighting factor correction : 0.000660805915266 average: 0.000660791703336 Largest master weighting factor correction: 3.33066907388e16 average: 6.42760698467e17 smoothSolver: Solving for Ux, Initial residual = 4.17822084748e05, Final residual = 4.06156755246e11, No Iterations 5 smoothSolver: Solving for Uy, Initial residual = 0.000184161482139, Final residual = 1.619696451e10, No Iterations 5 PCG: Solving for p, Initial residual = 0.00111502624192, Final residual = 9.94145230431e07, No Iterations 280 PCG: Solving for p, Initial residual = 9.81284451381e05, Final residual = 9.83206605956e07, No Iterations 11 time step continuity errors : sum local = 1.6066266089e11, global = 1.6661252468e13, cumulative = 2.44473486525e08 PCG: Solving for p, Initial residual = 2.56735821773e05, Final residual = 9.78307634024e07, No Iterations 43 PCG: Solving for p, Initial residual = 1.10964929717e05, Final residual = 8.87408739619e07, No Iterations 3 time step continuity errors : sum local = 1.45008481624e11, global = 1.02944097058e14, cumulative = 2.44473383581e08 PBiCG: Solving for omega, Initial residual = 1.98451187283e05, Final residual = 2.70582444337e08, No Iterations 1 PBiCG: Solving for k, Initial residual = 8.34347962654e05, Final residual = 2.51644588879e07, No Iterations 1 ExecutionTime = 16219.97 s ClockTime = 16311 s Courant Number mean: 0.0394284702228 max: 0.499946124724 velocity magnitude: 124.359835338 deltaT = 9.57170436146e06 Time = 0.316535 Initializing the GGI interpolator between master/shadow patches: insideSlider/outsideSlider Evaluation of GGI weighting factors: Largest slave weighting factor correction : 0.000723788760952 average: 0.000723774752211 Largest master weighting factor correction: 3.33066907388e16 average: 9.20316454623e17 smoothSolver: Solving for Ux, Initial residual = 4.17873817625e05, Final residual = 4.07204477989e11, No Iterations 5 smoothSolver: Solving for Uy, Initial residual = 0.000184210804886, Final residual = 1.62641723107e10, No Iterations 5 PCG: Solving for p, Initial residual = 0.00110876280886, Final residual = 9.61326185025e07, No Iterations 204 PCG: Solving for p, Initial residual = 9.64025215336e05, Final residual = 9.69402485222e07, No Iterations 11 time step continuity errors : sum local = 1.58389919616e11, global = 1.25756269618e12, cumulative = 2.44485959208e08 PCG: Solving for p, Initial residual = 2.5342476828e05, Final residual = 9.85823081162e07, No Iterations 41 PCG: Solving for p, Initial residual = 1.10396708963e05, Final residual = 9.27978145016e07, No Iterations 3 time step continuity errors : sum local = 1.51621484794e11, global = 1.19091666363e12, cumulative = 2.44497868375e08 [louisdell:03845] *** Process received signal *** [louisdell:03845] Signal: Floating point exception (8) [louisdell:03845] Signal code: (6) [louisdell:03845] Failing at address: 0x3e800000f05 [louisdell:03846] *** Process received signal *** [louisdell:03846] Signal: Floating point exception (8) [louisdell:03846] Signal code: (6) [louisdell:03846] Failing at address: 0x3e800000f06 [louisdell:03845] [ 0] /lib/libc.so.6 [0x7f46d1640530] [louisdell:03845] [ 1] /lib/libc.so.6(gsignal+0x35) [0x7f46d16404b5] [louisdell:03845] [ 2] /lib/libc.so.6 [0x7f46d1640530] [louisdell:03845] [ 3] /home/louis/OpenFOAM/OpenFOAM1.5dev/lib/linux64GccDPOpt/libOpenFOAM.so(_ZNK4Foam5PBiCG5solveERNS_5FieldIdEERKS2_h+0xbd8) [0x7f46d255f418] [louisdell:03845] [ 4] /home/louis/OpenFOAM/OpenFOAM1.5dev/lib/linux64GccDPOpt/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE5solveERNS_7IstreamE+0x13f) [0x7f46d327b71f] [louisdell:03845] [ 5] /home/louis/OpenFOAM/OpenFOAM1.5dev/lib/linux64GccDPOpt/libincompressibleRASModels.so(_ZN4Foam5solveIdEENS_9lduMatrix17solverPerformanceERKNS_3tmpINS_8fvMatrixIT_EEEE+0x4e) [0x7f46d3e7fc0e] [louisdell:03845] [ 6] /home/louis/OpenFOAM/OpenFOAM1.5dev/lib/linux64GccDPOpt/libincompressibleRASModels.so(_ZN4Foam14incompressible9RASModels9kOmegaSST7correctEv+0x123f) [0x7f46d3ebe79f] [louisdell:03845] [ 7] turbDyMFoam [0x41a488] [louisdell:03845] [ 8] /lib/libc.so.6(__libc_start_main+0xfd) [0x7f46d162babd] [louisdell:03845] [ 9] turbDyMFoam [0x415169] [louisdell:03845] *** End of error message *** [louisdell:03846] [ 0] /lib/libc.so.6 [0x7ff0e7f99530] [louisdell:03846] [ 1] /lib/libc.so.6(gsignal+0x35) [0x7ff0e7f994b5] [louisdell:03846] [ 2] /lib/libc.so.6 [0x7ff0e7f99530] [louisdell:03846] [ 3] /home/louis/OpenFOAM/OpenFOAM1.5dev/lib/linux64GccDPOpt/libOpenFOAM.so(_ZNK4Foam5PBiCG5solveERNS_5FieldIdEERKS2_h+0xbd8) [0x7ff0e8eb8418] [louisdell:03846] [ 4] /home/louis/OpenFOAM/OpenFOAM1.5dev/lib/linux64GccDPOpt/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE5solveERNS_7IstreamE+0x13f) [0x7ff0e9bd471f] [louisdell:03846] [ 5] /home/louis/OpenFOAM/OpenFOAM1.5dev/lib/linux64GccDPOpt/libincompressibleRASModels.so(_ZN4Foam5solveIdEENS_9lduMatrix17solverPerformanceERKNS_3tmpINS_8fvMatrixIT_EEEE+0x4e) [0x7ff0ea7d8c0e] [louisdell:03846] [ 6] /home/louis/OpenFOAM/OpenFOAM1.5dev/lib/linux64GccDPOpt/libincompressibleRASModels.so(_ZN4Foam14incompressible9RASModels9kOmegaSST7correctEv+0x123f) [0x7ff0ea81779f] [louisdell:03846] [ 7] turbDyMFoam [0x41a488] [louisdell:03846] [ 8] /lib/libc.so.6(__libc_start_main+0xfd) [0x7ff0e7f84abd] [louisdell:03846] [ 9] turbDyMFoam [0x415169] [louisdell:03846] *** End of error message *** mpirun noticed that job rank 0 with PID 3845 on node louisdell exited on signal 8 (Floating point exception). 1 additional process aborted (not shown) Thanks, Louis PS: I will provide the full logfile if anyone is interested. Last edited by louisgag; December 27, 2009 at 22:14. 

December 28, 2009, 03:54 

#6 
Senior Member
matej forman
Join Date: Mar 2009
Location: Brno, Czech Republic
Posts: 117
Rep Power: 10 
Alan,
I'd support the idea with mesh checking! Although I've never done any supersonic flow computations, I'm not very sure with your boundary cond. First if the obstacle is wall, it should have not a patch, but wall keyword in the boundary file, but I'm not sure if it makes difference. What is more important I think you need some more advanced BC tuned for supersonic pressure then zeroGradient. Looking at the tutorials for your solver I see there are some pressureTransmittive BCs. For example I do not understand your pressure at inlet = 1; Why 1? My guess would be identify the closest tutorial to your case and have a look at the tutorial settings for BC. good luck matej 

December 28, 2009, 19:01 

#7 
New Member
Alan Harrland
Join Date: Mar 2009
Posts: 21
Rep Power: 10 
Thanks for the replies guys. I think the problem may lie in my mesh. There are some highly skew elements along the wedge's surface, and when the shockwave reaches these elements, the simulation crashes.
I don't think there should be a problem with the BC's, as it is based on other supersonic tutorials in the rhoCentralFoam solver, and I have done other simulations without problems. I will try and play around with them though to see if it helps. Thanks again for the responses. Alan 

January 13, 2010, 20:30 

#8  
Member
Mihir
Join Date: Mar 2009
Posts: 40
Rep Power: 10 
Im trying to solve a freejet problem by rhoCentralFoam [OF 1.5]. My domain represents a close chamber in which the jet enters at Supersonic conditions.
My code crashes due to a FPE , I am getting ve pressure & density at the outlet . Quote:
Code:
dimensions [1 1 2 0 0 0 0]; internalField uniform 42000; boundaryField { wall { type slip; } axis { type empty; } outlet { type nonReflective; refValue uniform 6600; gamma 1.4; phi phi; psi psi; valueFraction uniform 0.2; refGradient uniform 0.2; value uniform 6600; } inlet { type totalPressure; U U; phi phi; rho none; psi none; gamma 1.4; p0 uniform 101325; value uniform 101325; } frontAndBackPlanes { type empty; } front { type wedge; } back { type wedge; } } Code:
dimensions [0 1 1 0 0 0 0]; internalField uniform (260.39 0 0); boundaryField { wall { type slip; } axis { type empty; } outlet { type zeroGradient; } inlet { type fixedValue; value uniform (315 0 0); } frontAndBackPlanes { type empty; } front { type wedge; } back { type wedge; } } Code:
thermoType hThermo<pureMixture<sutherlandTransport<specieThermo<hConstThermo<perfectGas>>>>>; mixture perfectGas 1 28.9 1004.3 200.0 1.458e6 110.4; rhoMin rhoMin [1 3 0 0 0] 0.001; Last edited by mihir1310; January 13, 2010 at 21:02. 

February 8, 2012, 04:31 
Floating error

#9 
Member
Join Date: Jan 2012
Posts: 58
Rep Power: 7 
Hello
i am getting floating error:invalid number. I am working on train travelling in tunnel. Can any body help me. My email is sheikhnasir39@gmail.com thanks 

April 6, 2012, 13:33 

#10  
New Member
Join Date: Jan 2010
Posts: 23
Rep Power: 9 
Quote:
Do you have timestep adjustment on? I noticed a similar behavior when using a moving mesh with GGI doing cell layer/addition. Residuals were fine and then suddenly an FPE. I could restart from the latest time step and it would make it past that point, but then at some later time it would happen. I eventually set an maximum timestep that was pretty low, and the problem went away  no more need to baby sit and keep restarting. Don't know if that helps. 

April 6, 2012, 14:02 

#11 
Senior Member

Hi James,
it has been quite a while since I found a way out of this problem, so I don't recall exactly how I did it. Thanks for the tip, I will keep it in mind if I stumble against such a problem again. Louis 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Installation of Netgen in SuSE Linux 92  edvardsenpriv  Open Source Meshers: Gmsh, Netgen, CGNS, ...  23  January 16, 2009 07:12 
Problem with compile the setParabolicInlet  ivanyao  OpenFOAM Running, Solving & CFD  6  September 5, 2008 20:50 
How to get the max value of the whole field  waynezw0618  OpenFOAM Running, Solving & CFD  4  June 17, 2008 05:07 
Gmsh and samplesurface  touf  Open Source Meshers: Gmsh, Netgen, CGNS, ...  2  December 10, 2007 03:27 
Floating point error  Nabil  FLUENT  2  January 19, 2007 06:30 