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

simpleFoam error - "Floating point exception"

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 12, 2014, 11:17
Default simpleFoam error - "Floating point exception"
  #1
New Member
 
James
Join Date: Jun 2012
Posts: 8
Rep Power: 5
mbcx4jc2 is on a distinguished road
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_64-linux-gnu/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_64-linux-gnu/libc.so.6"
#7
in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/simpleFoam"
Floating point exception (core dumped)
Attached Files
File Type: txt simpleFoamLog.txt (5.8 KB, 32 views)
mbcx4jc2 is offline   Reply With Quote

Old   February 12, 2014, 12:01
Default
  #2
Senior Member
 
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,114
Rep Power: 19
alexeym will become famous soon enoughalexeym will become famous soon enough
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.
Attached Files
File Type: gz system.tar.gz (1.1 KB, 27 views)
alexeym is offline   Reply With Quote

Old   February 12, 2014, 13:53
Default
  #3
New Member
 
James
Join Date: Jun 2012
Posts: 8
Rep Power: 5
mbcx4jc2 is on a distinguished road
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?
mbcx4jc2 is offline   Reply With Quote

Old   February 12, 2014, 14:27
Default
  #4
Senior Member
 
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,114
Rep Power: 19
alexeym will become famous soon enoughalexeym will become famous soon enough
Well,

You can try to estimate velocity in the narrow part from cross-section 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.
alexeym is offline   Reply With Quote

Old   March 5, 2014, 22:59
Unhappy error:time step continuty errors increasing and diverged
  #5
New Member
 
卞凤娇
Join Date: Nov 2013
Location: beijing
Posts: 15
Rep Power: 3
jiaojiao is on a distinguished road
Quote:
Originally Posted by mbcx4jc2 View Post
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_64-linux-gnu/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_64-linux-gnu/libc.so.6"
#7
in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/simpleFoam"
Floating point exception (core dumped)
I have a similar problem,I am running the MRFSimplefoam with of2.1 to simulate the 3D wind turbine, but it cause an error as I pasted
Code:
#0  Foam::error::printStack(Foam::Ostream&) in "/home/user5/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64Gcc45DPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/home/user5/OpenFOAM/OpenFOAM-2.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/OpenFOAM-2.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/OpenFOAM-2.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/OpenFOAM-2.1.0/platforms/linux64Gcc45DPOpt/bin/foamToVTK"
#6  main in "/home/user5/OpenFOAM/OpenFOAM-2.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/OpenFOAM-2.1.0/platforms/linux64Gcc45DPOpt/bin/foamToVTK"
浮点数例外
can anybody help me? Where the problems are and how can I fix it? my fvshemes file is
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((1|A(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               ;
}
the log is
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 1e-05
    field U     tolerance 1e-05
    field nuTilda     tolerance 1e-05


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.54463e-06, No Iterations 2
ExecutionTime = 454.86 s  ClockTime = 497 s
when it run to steps21,the errors take placed.
jiaojiao is offline   Reply With Quote

Old   March 6, 2014, 01:05
Default
  #6
New Member
 
Bin Xu
Join Date: Apr 2012
Location: Singapore
Posts: 23
Rep Power: 5
norkistar is on a distinguished road
your first step : time step already diverge, you can have a check about your mesh and your boundary condition settings
norkistar is offline   Reply With Quote

Old   March 7, 2014, 03:55
Default
  #7
New Member
 
卞凤娇
Join Date: Nov 2013
Location: beijing
Posts: 15
Rep Power: 3
jiaojiao is on a distinguished road
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 (non-closed singly connected)
outlet 1841 966 ok (non-closed singly connected)
symmetry 3560 3649 ok (non-closed singly connected)

Checking geometry...
Overall domain bounding box (-12.7106 -6.66098 -12.7139) (12.7102 27.7695 12.7057)
Mesh (non-empty, non-wedge) directions (1 1 1)
Mesh (non-empty) directions (1 1 1)
Boundary openness (-2.3031e-15 -3.09845e-16 5.48772e-16) OK.
Max cell openness = 1.07338e-15 OK.
Max aspect ratio = 30.0527 OK.
Minimum face area = 6.8698e-07. Maximum face area = 0.772395. Face area magnitudes OK.
Min volume = 3.21357e-10. Max volume = 0.154152. Total volume = 17466.7. Cell volumes OK.
Mesh non-orthogonality Max: 70.2416 average: 19.945
*Number of severely non-orthogonal (> 70 degrees) faces: 3.
Non-orthogonality check OK.
<<Writing 3 non-orthogonal 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 non-orthogonal faces or something eles?
jiaojiao is offline   Reply With Quote

Old   March 13, 2014, 21:09
Default
  #8
New Member
 
卞凤娇
Join Date: Nov 2013
Location: beijing
Posts: 15
Rep Power: 3
jiaojiao is on a distinguished road
Quote:
Originally Posted by norkistar View Post
your first step : time step already diverge, you can have a check about your mesh and your boundary condition settings
Hi,I now use Gambit to generate mesh, boundary condition doesnt change and it can run more than 150 steps with the first time steps continuty errors OK. I paste the log file,could you tell me where the problem is?the eorror messege is
Code:
#0  Foam::error::printStack(Foam::Ostream&) in "/home/user5/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64Gcc45DPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/home/user5/OpenFOAM/OpenFOAM-2.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/OpenFOAM-2.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/OpenFOAM-2.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/OpenFOAM-2.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/OpenFOAM-2.1.0/platforms/linux64Gcc45DPOpt/lib/libOpenFOAM.so"
#7  Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/home/user5/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64Gcc45DPOpt/lib/libfiniteVolume.so"
#8  Foam::fvMatrix<double>::solve() in "/home/user5/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64Gcc45DPOpt/bin/MRFSimpleFoam"
#9  main in "/home/user5/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64Gcc45DPOpt/bin/MRFSimpleFoam"
#10  __libc_start_main in "/lib64/libc.so.6"
#11  Foam::UOPstream::write(char) in "/home/user5/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64Gcc45DPOpt/bin/MRFSimpleFoam"
娴偣鏁颁緥澶
jiaojiao is offline   Reply With Quote

Old   March 13, 2014, 21:14
Default
  #9
New Member
 
卞凤娇
Join Date: Nov 2013
Location: beijing
Posts: 15
Rep Power: 3
jiaojiao is on a distinguished road
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.91772e-08, Final residual = 9.91772e-08, 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.47159e-07, 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.85788e-08, Final residual = 9.85788e-08, 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.54768e-07, 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.65561e-08, Final residual = 9.65561e-08, 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.12524e-08, 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.28703e-08, Final residual = 9.28703e-08, 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.31763e-08, 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.84393e-08, Final residual = 8.84393e-08, 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.93949e-08, 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.29655e-08, Final residual = 8.29655e-08, 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.66331e-08, 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
jiaojiao is offline   Reply With Quote

Old   April 21, 2015, 13:35
Default Floating point exception (core dumped)- Immersed boundary- simpleIbFoam
  #10
New Member
 
Anamitra Bagchi
Join Date: Mar 2015
Posts: 3
Rep Power: 2
Anamitra is on a distinguished road
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@anamitra-HP-Pavilion-g6-Notebook-PC:~/IC_Files/Empty_Domain/simpleibfoam_vg$ simpleIbFoam
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 3.1 |
| \\ / A nd | Web: http://www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 3.1-7d8e040bf53d
Exec : simpleIbFoam
Date : Apr 21 2015
Time : 19:30:14
Host : anamitra-HP-Pavilion-g6-Notebook-PC
PID : 8630
CtrlDict : /home/anamitra/foam/foam-extend-3.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@anamitra-HP-Pavilion-g6-Notebook-PC:~/IC_Files/Empty_Domain/simpleibfoam_vg$
Anamitra is offline   Reply With Quote

Old   April 21, 2015, 19:45
Default
  #11
New Member
 
N/A
Join Date: Jul 2010
Posts: 27
Rep Power: 7
kamakura117 is on a distinguished road
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.
kamakura117 is offline   Reply With Quote

Old   June 25, 2015, 13:53
Default
  #12
New Member
 
Anamitra Bagchi
Join Date: Mar 2015
Posts: 3
Rep Power: 2
Anamitra is on a distinguished road
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;
}
}
Anamitra is offline   Reply With Quote

Old   Yesterday, 02:20
Default
  #13
New Member
 
Kittipos
Join Date: Feb 2015
Posts: 7
Rep Power: 2
Kittipos is on a distinguished road
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/i386-linux-gnu/libc.so.6"
#7  
 in "/opt/openfoam221/platforms/linuxGccDPOpt/bin/rhoPorousSimpleFoam"
Floating point exception (core dumped)
Thank you
Attached Files
File Type: txt log.txt (47.7 KB, 0 views)
Kittipos is online now   Reply With Quote

Reply

Tags
error, floating point exception, openfoam, simplefoam

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
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
problem in toposet Ahmed Khattab OpenFOAM Pre-Processing 10 March 24, 2013 16:18
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


All times are GMT -4. The time now is 03:17.