CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Floating point exception (core dumped) for GAMG solver

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 24, 2015, 08:07
Default Floating point exception (core dumped) for GAMG solver
  #1
New Member
 
Yolanda Hou
Join Date: Mar 2015
Location: Southampton
Posts: 2
Rep Power: 0
yuhou1989 is on a distinguished road
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:
Attached Files
File Type: txt U.txt (1.4 KB, 4 views)
File Type: txt p.txt (1.3 KB, 4 views)
File Type: txt nut.txt (1.6 KB, 2 views)
File Type: txt k.txt (1.4 KB, 5 views)
File Type: txt epsilon.txt (1.4 KB, 3 views)
yuhou1989 is offline   Reply With Quote

Old   March 24, 2015, 08:08
Default
  #2
New Member
 
Yolanda Hou
Join Date: Mar 2015
Location: Southampton
Posts: 2
Rep Power: 0
yuhou1989 is on a distinguished road
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;
}

// ************************************************** *********************** //
yuhou1989 is offline   Reply With Quote

Old   March 24, 2015, 20:28
Default core dumped
  #3
New Member
 
QP Dong
Join Date: Feb 2015
Posts: 7
Rep Power: 9
high_er is on a distinguished road
Quote:
Originally Posted by yuhou1989 View Post
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;
}

// ************************************************** *********************** //
Hi
Have you fixed this problem? How?
Thanks
high_er is offline   Reply With Quote

Reply

Tags
simplefoam gamg

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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 Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


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


All times are GMT -4. The time now is 14:26.