|
[Sponsors] |
Floating point exception (core dumped) for GAMG solver |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
March 24, 2015, 08:07 |
Floating point exception (core dumped) for GAMG solver
|
#1 |
New Member
Yolanda Hou
Join Date: Mar 2015
Location: Southampton
Posts: 2
Rep Power: 0 |
Hi, I am trying to simulate the flow field around a circular cylinder using Detached-Eddy Simulation and would like to run a RANS simulation with realizable k-e model as its initial conditions. However the following floating point exception at the very beginning of the simulation:
/*---------------------------------------------------------------------------*\ | ========= | | | \\ / 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 : simpleFoam Date : Mar 24 2015 Time : 10:52:26 Host : "cyan01" PID : 16434 Case : /scratch/yh25g11/yaw0DES-rke 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 RAS turbulence model SpalartAllmaras SpalartAllmarasCoeffs { sigmaNut 0.66666; kappa 0.41; Cb1 0.1355; Cb2 0.622; Cw2 0.3; Cw3 2; Cv1 7.1; Cv2 5; } Employing Ashford correction No finite volume options present SIMPLE: no convergence criteria found. Calculations will run for 10000 steps. Starting time loop forceCoeffs forcesCylinder: Not including porosity effects Time = 1 DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 0.04933216, No Iterations 2 DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 0.05075823, No Iterations 2 DILUPBiCG: Solving for Uz, Initial residual = 0.9999973, Final residual = 0.02594695, No Iterations 1 #0 Foam::error:rintStack(Foam::Ostream&) in "/home/yh25g11/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #1 Foam::sigFpe::sigHandler(int) in "/home/yh25g11/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #2 at sigaction.c:0 #3 Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const in "/home/yh25g11/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #4 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const in "/home/yh25g11/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #5 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/yh25g11/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #6 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) in "/home/yh25g11/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so" #7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/home/yh25g11/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/bin/simpleFoam" #8 Foam::fvMatrix<double>::solve() in "/home/yh25g11/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/bin/simpleFoam" #9 in "/home/yh25g11/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/bin/simpleFoam" #10 __libc_start_main in "/lib64/libc.so.6" #11 in "/home/yh25g11/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/bin/simpleFoam" Floating point exception (core dumped) And the boundary conditions are set as attached: |
|
March 24, 2015, 08:08 |
|
#2 |
New Member
Yolanda Hou
Join Date: Mar 2015
Location: Southampton
Posts: 2
Rep Power: 0 |
And the schemes are set as follows:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.1.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) cellLimited Gauss linear 1; grad(nuTilda) Gauss linear; } divSchemes { default none; div(phi,U) bounded Gauss linearUpwindV grad(U); div(phi,k) bounded Gauss upwind; div(phi,epsilon) bounded Gauss upwind; //div(phi,nuTilda) bounded Gauss upwind; div((nuEff*dev(T(grad(U))))) Gauss linear; } laplacianSchemes { default Gauss linear corrected; laplacian(nuEff,U) Gauss linear corrected; laplacian((1|A(U)),p) Gauss linear corrected; laplacian(DnuTildaEff,nuTilda) Gauss linear corrected; } interpolationSchemes { default linear; interpolate(HbyA) linear; } snGradSchemes { default corrected; } fluxRequired { default no; p; } // ************************************************** *********************** // |
|
March 24, 2015, 20:28 |
core dumped
|
#3 | |
New Member
QP Dong
Join Date: Feb 2015
Posts: 7
Rep Power: 11 |
Quote:
Have you fixed this problem? How? Thanks |
||
Tags |
simplefoam gamg |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
2D NACA0012 O-Mesh GAMG Solver Floating Point Exception | mdeaves | OpenFOAM Running, Solving & CFD | 4 | September 23, 2019 03:51 |
simpleFoam: floating point exception | student666 | OpenFOAM Running, Solving & CFD | 9 | March 13, 2016 19:23 |
[snappyHexMesh] snappyHexMesh and cyclic boundaries | Ruli | OpenFOAM Meshing & Mesh Conversion | 2 | December 9, 2013 07:51 |
ERROR : Floating point exception (core dumped) | sam.ho | OpenFOAM Running, Solving & CFD | 4 | December 3, 2013 06:25 |
floating point exception !! | SamerAli | OpenFOAM Programming & Development | 1 | May 18, 2013 04:43 |