# floating point exception

 Register Blogs Members List Search Today's Posts Mark Forums Read

 March 7, 2014, 20:31 floating point exception #1 New Member   J. Brindley Join Date: Mar 2014 Posts: 11 Rep Power: 4 Not sure what's going on. Can't get simpleFoam to run Code: ```Time = 1 smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.081410493, No Iterations 4 smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.061596729, No Iterations 3 smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 0.091684187, No Iterations 3 GAMG: Solving for p, Initial residual = 1, Final residual = 0.0059667553, No Iterations 5 GAMG: Solving for p, Initial residual = 0.40650934, Final residual = 0.0012374181, No Iterations 3 GAMG: Solving for p, Initial residual = 0.036894522, Final residual = 0.00027957842, No Iterations 3 GAMG: Solving for p, Initial residual = 0.0062564286, Final residual = 4.8770264e-05, No Iterations 3 GAMG: Solving for p, Initial residual = 0.0013606213, Final residual = 1.2695641e-05, No Iterations 3 GAMG: Solving for p, Initial residual = 0.00035198339, Final residual = 3.117196e-06, No Iterations 3 GAMG: Solving for p, Initial residual = 0.00010876533, Final residual = 1.0124916e-06, No Iterations 3 GAMG: Solving for p, Initial residual = 3.7107364e-05, Final residual = 3.2700493e-07, No Iterations 3 GAMG: Solving for p, Initial residual = 1.4011258e-05, Final residual = 1.2603484e-07, No Iterations 3 GAMG: Solving for p, Initial residual = 5.5796704e-06, Final residual = 3.7850784e-08, No Iterations 6 GAMG: Solving for p, Initial residual = 2.3340101e-06, Final residual = 6.442118e+119, No Iterations 1000 time step continuity errors : sum local = 1.2632859e+123, global = 4.4873194e+101, cumulative = 4.4873194e+101 #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 in "/lib/x86_64-linux-gnu/libc.so.6" #3 Foam::GaussSeidelSmoother::smooth(Foam::word const&, Foam::Field&, Foam::lduMatrix const&, Foam::Field const&, Foam::FieldField const&, Foam::UPtrList const&, unsigned char, int) at ??:? #4 Foam::GaussSeidelSmoother::smooth(Foam::Field&, Foam::Field const&, unsigned char, int) const at ??:? #5 Foam::smoothSolver::solve(Foam::Field&, Foam::Field const&, unsigned char) const at ??:? #6 Foam::fvMatrix::solveSegregated(Foam::dictionary const&) at ??:? #7 Foam::fvMatrix::solve(Foam::dictionary const&) at ??:? #8 Foam::SolverPerformance Foam::solve(Foam::tmp > const&) at ??:? #9 Foam::incompressible::RASModels::kOmegaSST::correct() at ??:? #10 at ??:? #11 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #12 at ??:? Floating point exception (core dumped)```

 March 7, 2014, 21:06 #2 New Member   J. Brindley Join Date: Mar 2014 Posts: 11 Rep Power: 4 [ Code: ```/*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.3.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : 2.3.0-f5222ca19ce6 Exec : checkMesh Date : Mar 07 2014 Time : 20:04:09 Host : "AVA-411296" PID : 26818 Case : /home/kathyb/OpenFOAM/kathyb-2.3.0/run/SplitTipTesting/AOA/simplesolvefixed(donotuse) 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 polyMesh for time = 0 Time = 0 Mesh stats points: 30927 faces: 72559 internal faces: 60772 cells: 21783 faces per cell: 6.1208741 boundary patches: 5 point zones: 0 face zones: 0 cell zones: 0 Overall number of cells of each type: hexahedra: 15749 prisms: 1429 wedges: 0 pyramids: 0 tet wedges: 128 tetrahedra: 2 polyhedra: 4475 Breakdown of polyhedra by number of faces: faces number of cells 4 1186 5 814 6 740 7 218 8 90 9 678 10 86 11 47 12 401 13 3 14 6 15 150 18 56 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 inlet 0 0 ok (empty) outlet 0 0 ok (empty) walls 40 67 ok (non-closed singly connected) frontBack 0 0 ok (empty) Foil 11747 13554 ok (non-closed singly connected) Checking geometry... Overall domain bounding box (-0.020220271 -0.010812593 -0.00023150215) (0.020260984 0.010818352 0.15183358) Mesh (non-empty, non-wedge) directions (1 1 1) Mesh (non-empty) directions (1 1 1) Boundary openness (-5.8212676e-17 4.5742309e-17 1.775036e-17) OK. Max cell openness = 2.9738872e-16 OK. Max aspect ratio = 21.38308 OK. Minimum face area = 6.6009305e-09. Maximum face area = 1.3716404e-05. Face area magnitudes OK. Min volume = 1.1250083e-12. Max volume = 2.2923998e-08. Total volume = 2.8560307e-05. Cell volumes OK. Mesh non-orthogonality Max: 62.160689 average: 13.151348 Non-orthogonality check OK. Face pyramids OK. Max skewness = 2.6576963 OK. Coupled point location match (average 0) OK. Mesh OK. End``` Code: ```/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.2.1 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; object fvSolution; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // solvers { p { solver GAMG; tolerance 0.02; relTol 0; smoother GaussSeidel; nPreSweeps 0; nPostSweeps 2; cacheAgglomeration on; agglomerator faceAreaPair; nCellsInCoarsestLevel 10; mergeLevels 1; } U { solver PBiCG; preconditioner DILU; tolerance 1e-6; relTol 0.1; } k { solver PBiCG; preconditioner DILU; tolerance 1e-6; relTol 0.1; } omega { solver PBiCG; preconditioner DILU; tolerance 1e-6; relTol 0.1; } cellDisplacement { solver GAMG; tolerance 1e-5; relTol 0; smoother GaussSeridel; cacheAgglomeration true; nCellsInCoursestLevel 10; agglomerator faceAreaPair; mergeLevels 1; } } SIMPLE { nNonOrthogonalCorrectors 0; pReffCell 0; pRefValue 0; residualControl { p 1e-4; U 1e-4; } } potentialFlow { nNonOrthogonalCorrectors 10; } relaxationFactors { fields { p 0.3; } equations { U 0.7; k 0.7; omega 0.7; } } cache { grad(U); } // ************************************************************************* //``` Code: ```/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.2.1 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; object fvSchemes; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // ddtSchemes { default steadyState; } gradSchemes { default Gauss linear; grad(p) Gauss linear; grad(U) Gauss linear; grad(k) cellLimited Gauss linear 1; grad(omega) cellLimited Gauss linear 1; } divSchemes { default none; div(phi,U) bounded Gauss upwind;; div(phi,k) bounded Gauss upwind; div(phi,omega) bounded Gauss upwind; div((nuEff*dev(T(grad(U))))) Gauss linear; } laplacianSchemes { default Gauss linear limited corrected 0.5; laplacian(DomegaEff,omega) Gauss linear limited 0.33; } interpolationSchemes { default linear; } snGradSchemes { default limited 0.33; } fluxRequired { default no; p; } // ************************************************************************* //```

 March 8, 2014, 17:56 #3 Senior Member   Lieven Join Date: Dec 2011 Location: Leuven, Belgium Posts: 297 Rep Power: 15 Hi olddirtybrindly, Do you get the same crash when you run the case with pisofoam? Cheers, L

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post ebah6 OpenFOAM Running, Solving & CFD 8 May 13, 2013 09:01 Janshi STAR-CCM+ 4 March 14, 2012 11:21 sudhasran OpenFOAM Running, Solving & CFD 3 March 12, 2012 17:23 jochemvandenbosch OpenFOAM Running, Solving & CFD 4 February 16, 2012 04:24 Robert@cfd ANSYS Meshing & Geometry 20 November 11, 2011 05:59

All times are GMT -4. The time now is 02:36.