[shadabdyn]

Hello everyone,

I am using my self-made code to run the simulation in annulus geometry. My simulation was core dumping at this time even though my courant number is still less than 1. please help me to come out of this problem.The output of terminal is given below.

Time = 240.999999999998

Courant Number mean: 0.0249712 max: 0.0875329
DICPCG:  Solving for Ux, Initial residual = 0.105497, Final residual = 0.00884902, No Iterations 1
DICPCG:  Solving for Uy, Initial residual = 0.469317, Final residual = 0.0347587, No Iterations 1
DICPCG:  Solving for p, Initial residual = 0.0102862, Final residual = 0.000474484, No Iterations 8
time step continuity errors : sum local = 2.46663e-06, global = -1.54741e-20, cumulative = 6.03058e-20
DICPCG:  Solving for Ux, Initial residual = 0.0414099, Final residual = 0.000467468, No Iterations 2
DICPCG:  Solving for Uy, Initial residual = 0.44559, Final residual = 0.0268967, No Iterations 1
DICPCG:  Solving for p, Initial residual = 0.00748688, Final residual = 0.000345887, No Iterations 142
time step continuity errors : sum local = 1.3637e-06, global = -1.29084e-21, cumulative = 5.90149e-20
DILUPBiCG:  Solving for c, Initial residual = 0.000474526, Final residual = 2.22852e-09, No Iterations 1
ExecutionTime = 11.78 s


Time = 241.099999999998

Courant Number mean: 0.0249315 max: 0.0873926
DICPCG:  Solving for Ux, Initial residual = 0.234292, Final residual = 0.0211142, No Iterations 1
DICPCG:  Solving for Uy, Initial residual = 0.520758, Final residual = 0.00558537, No Iterations 2
DICPCG:  Solving for p, Initial residual = 0.0389959, Final residual = 0.00179561, No Iterations 35
time step continuity errors : sum local = 7.66976e-06, global = -7.74512e-21, cumulative = 5.12698e-20
DICPCG:  Solving for Ux, Initial residual = 0.112113, Final residual = 0.00129323, No Iterations 2
DICPCG:  Solving for Uy, Initial residual = 0.161259, Final residual = 0.0148576, No Iterations 1
DICPCG:  Solving for p, Initial residual = 0.012724, Final residual = 0.000543076, No Iterations 37
time step continuity errors : sum local = 2.841e-06, global = -8.35252e-22, cumulative = 5.04346e-20
DILUPBiCG:  Solving for c, Initial residual = 0.000661722, Final residual = 2.16062e-09, No Iterations 1
ExecutionTime = 11.81 s


Time = 241.199999999998

Courant Number mean: 0.0249643 max: 0.0874092
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2   in "/lib/x86_64-linux-gnu/libc.so.6"
#3   in "/lib/x86_64-linux-gnu/libm.so.6"
#4  pow in "/lib/x86_64-linux-gnu/libm.so.6"
#5  Foam::pow(Foam::Field<double>&, Foam::UList<double> const&, double const&) at ??:?
#6  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::pow<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::dimensioned<double> const&) at ??:?
#7  
 at ??:?
#8  
 at ??:?
#9  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10  
 at ??:?
Floating point exception (core dumped)
shadab@shadab-HP-ProDesk-600-G1-TWR:~/OpenFOAM/OpenFOAM-2.3.0/Project/cavity_settling2$

Thanks in advace.

[KUMAR7491]

I have added bottom friction as a source term in shallowwaterfoam.
bottom stress= Sp(2*mag(U), U)

After compiling and running i get error message as below.

Creating field hTotal for post processing

Reading/calculating face flux field phi

Creating Coriolis Force

PIMPLE: no residual control data found. Calculations will employ 3 corrector loops


Starting time loop


Time = 0.1

Courant number mean: 0.200006 max: 0.205556
Gravity wave Courant number max: 0.626418
PIMPLE: iteration 1


--> FOAM FATAL ERROR:
incompatible fields for operation
[hU] + [U]

From function checkMethod(const fvMatrix<Type>&, const fvMatrix<Type>&)
in file /opt/openfoam211/src/finiteVolume/lnInclude/fvMatrix.C at line 1303.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2
in "/home/caelinux/OpenFOAM/caelinux-2.1.1/platforms/linux64GccDPOpt/bin/my_shallowWaterFoam"
#3
in "/home/caelinux/OpenFOAM/caelinux-2.1.1/platforms/linux64GccDPOpt/bin/my_shallowWaterFoam"
#4
in "/home/caelinux/OpenFOAM/caelinux-2.1.1/platforms/linux64GccDPOpt/bin/my_shallowWaterFoam"
#5 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#6
in "/home/caelinux/OpenFOAM/caelinux-2.1.1/platforms/linux64GccDPOpt/bin/my_shallowWaterFoam"
Aborted (core dumped)
caelinux@caelinux:~/OpenFOAM/caelinux-2.1.1/squareBump$

[esujby]

Hello Shabad and Kumar,

Have you checked your mesh to ensure its ok to run the simulation - you can use the following to do just that:

https://openfoamwiki.net/index.php/CheckMesh

it is a print stack error which means maybe one of your fields is been divided by 0 or some other unrecognized character.

lastly, i will advice you both follow the method outlined in the first post of the link bellow to effectively report your problems and get better help. it helps the readers understand what exactly you are trying to achieve. i know its a pain in the *** but it would help you allot.

How to give enough info to get help

[andreamc]

Hello
Sorry to activate this old post. But I'm having a similar problem.
I am running a case with rhoSimpleFoam and SA turbulence model on OF v. 3.1. I've ran the same case on three coarse grids; but whem I'm gonna run the exact same case on a refined mesh it crashes and I get an error similar to the posted above:

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::hePsiThermo<Foam:siThermo, Foam:ureMixture<Foam::sutherlandTransport<Foam:: species::thermo<Foam::hConstThermo<Foam:erfectGa s<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() at ??:?
#4 Foam::hePsiThermo<Foam:siThermo, Foam:ureMixture<Foam::sutherlandTransport<Foam:: species::thermo<Foam::hConstThermo<Foam:erfectGa s<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() at ??:?
#5 ? at ??:?
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7 ? at ??:?

I've tried to:
-change boundary conditions of pressure and nut, but this has not have any effect on the simulation
-change schemes but this doesn't show any effect on the simulation.
-I map the last time step field from the coarse mesh simulation to the refined mesh but the error keeps happening.
-I've change thermophysicalProperties file, the transport method but then I got a different error of MAX NUMBER OF ITERATIONS EXCEED

I'm totally clueless about what is happening. And I do not know how can I solve this. I would be very thankful if anyone can give me a hint.

Regards
Andrea

[iamranjitsingh]

Hello esujby

thanks for reply,, i got my solution after that problem,,,,,,

[haydii]

Hi
I want to simulate 3d flow over weir but when i try to run it with interFoam i see "Floating point exception (core dumped)" error, what shuold i do to fixed it?

[esujby]

hmmmmm!

this could be caused by allot of things however, it seems like something is been divided by Zero.

if i were you, i will install the debug version. if nothing helpfull comes up, i will get a debug version.

have you been able to runt he case before? did you modify something...?

[haydii]

Originally Posted by esujby

hmmmmm!

this could be caused by allot of things however, it seems like something is been divided by Zero.

if i were you, i will install the debug version. if nothing helpfull comes up, i will get a debug version.

have you been able to runt he case before? did you modify something...?

thank you for your attention
yes i modified it for a 3d simulation with totally different boundary and initial condition,please give me some information about the debug version.

[sheaker]

Hello.
http://www.tfd.chalmers.se/~hani/kur.../debugging.pdf
Is this what Your are looking for?

I think your error is somewhere above statement of "Floating point exception (core dumped)". Look through Your log file for other errors.

[haydii]

Originally Posted by sheaker

Hello.
http://www.tfd.chalmers.se/~hani/kur.../debugging.pdf
Is this what Your are looking for?

I think your error is somewhere above statement of "Floating point exception (core dumped)". Look through Your log file for other errors.

Hi
thanks for reply
but I cant download the file

[esujby]

Hello Haydii,

to get the best help i recommend the following advice.

Originally Posted by esujby

Hello Rsingh,

details, details, details! the most important in order to get the most effective help. what are you trying to achieve?

in order to get the best help possible, please go through this thread.

http://www.cfd-online.com/Forums/ope...-get-help.html

kind regards

[haydii]

Originally Posted by esujby

Hello Haydii,

to get the best help i recommend the following advice.

Hi esujby
thanks so much for replay,i'll try to fixed it

Originally Posted by iamranjitsingh

Hello esujby

thanks for reply,, i got my solution after that problem,,,,,,

Hi iamranjitsingh
I have same problem ,could you help me about it?

[Alisa_W]

Originally Posted by ebah6

Hello all,

I am having a similar error in using pimpleDyMFoam.
Below is the error output:

------------------------------------
Courant Number mean: 0.00997899 max: 0.807965
deltaT = 1.32295e-104
--> FOAM Warning :
From function Time:perator++()
in file db/Time/Time.C at line 982
Increased the timePrecision from 267 to 268 to distinguish between timeNames at time 1.97982e-05
Time = 1.979823486337861903608045799352055382769322022795 67718505859375e-05

solidBodyMotionFunctions::rotatingMotion::transfor mation(): Time = 1.97982e-05 transformation: ((0 0 0) (1 (0 0 0.000103663)))
AMI: Creating addressing and weights between 16 source faces and 16 target faces
AMI: Patch source weights min/max/average = 1, 1.0007, 1.00035
AMI: Patch target weights min/max/average = 0.986951, 0.987248, 0.987099
smoothSolver: Solving for Ux, Initial residual = 0.140328, Final residual = 5.75579e-08, No Iterations 3
smoothSolver: Solving for Uy, Initial residual = 0.140962, Final residual = 5.70666e-08, No Iterations 3
GAMG: Solving for p, Initial residual = 0.814891, Final residual = 0.00591061, No Iterations 3
time step continuity errors : sum local = 0.00031744, global = 5.68686e-06, cumulative = 0.00082973
#0 Foam::error:rintStack(Foam::Ostream&) in "/home/alpha/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/home/alpha/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::GAMGSolver::scalingFactor(Foam::Field<double >&, Foam::Field<double> const&, Foam::Field<double> const&, Foam::Field<double> const&) const in "/home/alpha/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/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/alpha/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5 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::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const in "/home/alpha/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#6 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/alpha/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/home/alpha/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#8
in "/home/alpha/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/bin/pimpleDyMFoam"
#9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10
in "/home/alpha/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/bin/pimpleDyMFoam"
--------------------------

Can you tell how you went along in solving the issue?

Thank you for your help.

Have you solved your problem?I just met with the same questions as yours, After solving for U,the error appeared...Can you told me what may cause this problems?

[esujby]

hmm interesting! it could be allot of things but so long as your mesh is spot on, boundaries well defined and case set up well, then i think it could be linked to your tolerances in fvsolutions. so for you based on the attachement, i would say U is not getting solved correctly......tolerances should be smaller than 1e-07 i reckon.

[Alisa_W]

Originally Posted by esujby

hmm interesting! it could be allot of things but so long as your mesh is spot on, boundaries well defined and case set up well, then i think it could be linked to your tolerances in fvsolutions. so for you based on the attachement, i would say U is not getting solved correctly......tolerances should be smaller than 1e-07 i reckon.

Thank you for your nice advice~I have set the tolerances to around e-11 but the error is still the same. Are there any other possibilities?

[esujby]

heyya sorry for the delayed reply, been busy with my PhD thesis. well, from what i see it could be ur relaxation factors 0.2 seems quite low. firstly have you checked ur mesh to make sure its en pointe using "checkMesh" on terminal

[haydii]

Originally Posted by esujby

heyya sorry for the delayed reply, been busy with my PhD thesis. well, from what i see it could be ur relaxation factors 0.2 seems quite low. firstly have you checked ur mesh to make sure its en pointe using "checkMesh" on terminal

your welcome dear esujby, that problem has been solved,anyway so pleased for your response, thank you

[gouravjee]

Originally Posted by wyldckat

Greetings Almir,

At the risk of sending you off in the wrong direction, you can try this answer: My program stops with an output that starts with #0 Foam::error:: PrintStack(Foam::Ostream&)

But in an attempt to send you in the right direction:

1. You should pay closer attention to the output. For example:

   Build  : 1.7.x-3776603e4c6c
   Exec   : buoyantSimpleFoam
   Date   : Jun 15 2011
   Time   : 12:26:48
   Host   : ubuntu
   PID    : 5430
   Case   : /home/almir/OpenFOAM/zylinder
   nProcs : 1
    SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

   The line in bold tells you what SigFpe is: Floating point exception @wikipedia
2. The second line in the print stack says this:

   #1  Foam::sigFpe::sigFpeHandler(int) in "/opt/openfoam171/lib/linux64GccDPOpt/libOpenFOAM.so"

   Which means that some bad math went about doing something wrong... in other words, division by infinite or by zero or something like that.
3. Examining the iteration outputs, you will see the following breadcrumbs about the impeding doom that looms in the solver's horizon:

   Time = 1
   
   DILUPBiCG:  Solving for Ux, Initial residual = 1, Final residual = 0.00987294, No Iterations 1
   DILUPBiCG:  Solving for Uy, Initial residual = 1, Final residual = 0.0157, No Iterations 1
   DILUPBiCG:  Solving for Uz, Initial residual = 1, Final residual = 0.00987846, No Iterations 1
   DILUPBiCG:  Solving for h, Initial residual = 1, Final residual = 0.0105804, No Iterations 1
   GAMG:  Solving for p_rgh, Initial residual = 0.899378, Final residual = 0.00305647, No Iterations 4
   time step continuity errors : sum local = 18.9686, global = -1.40909e-15, cumulative = -1.40909e-15
   rho max/min : 1.22108 1.13449
   DILUPBiCG:  Solving for omega, Initial residual = 0.999913, Final residual = 0.0105502, No Iterations 2
   bounding omega, min: -902.694 max: 24331.5 average: 741.476
   DILUPBiCG:  Solving for k, Initial residual = 1, Final residual = 0.0973743, No Iterations 1
   bounding k, min: -0.000335494 max: 0.0031317 average: 0.000848312
   ExecutionTime = 0.11 s  ClockTime = 0 s
   
   Time = 2
   
   DILUPBiCG:  Solving for Ux, Initial residual = 0.11425, Final residual = 0.00308422, No Iterations 1
   DILUPBiCG:  Solving for Uy, Initial residual = 0.0715711, Final residual = 2.95363e-05, No Iterations 2
   DILUPBiCG:  Solving for Uz, Initial residual = 0.109402, Final residual = 0.00196119, No Iterations 1
   DILUPBiCG:  Solving for h, Initial residual = 0.177633, Final residual = 0.00377685, No Iterations 1
   GAMG:  Solving for p_rgh, Initial residual = 0.997056, Final residual = 0.00735437, No Iterations 4
   time step continuity errors : sum local = 11.1294, global = 4.01181e-15, cumulative = 2.60272e-15
   rho max/min : 309747 -321318
   DILUPBiCG:  Solving for omega, Initial residual = 0.594095, Final residual = 0.0340323, No Iterations 1
   bounding omega, min: -7.04999e+17 max: 1.99134e+09 average: -1.78979e+14
   DILUPBiCG:  Solving for k, Initial residual = 0.999984, Final residual = 0.0558935, No Iterations 2
   bounding k, min: -7.30668e+08 max: 1.51545e+08 average: -1.226e+06
   ExecutionTime = 0.15 s  ClockTime = 0 s

   As I've noted in bold+underline: rho, omega, k and continuity errors indicate that something very not physical is happening!! Expressions like über compression and super turbulence come to mind!

So, to sum up: you are giving very bad boundary conditions to your case!

Best regards,
Bruno

hello,
i am trying to convert a mesh from fluent to openfoam and i have followed a procedure as shown below ::
dos2unix mesh.msh (i am using fluent in windows)
fluentMeshToFoam mesh.msh but ige the following error but i get the following error.

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

number of faces: 3486
Number of points: 614
Reading uniform faces
Reading uniform faces
Reading points


FINISHED LEXING


#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigSegv::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3  ? at ??:?
#4  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#5  ? at ??:?
Segmentation fault (core dumped)

can you tell me where might be the problem??

[einstein_zee]

Originally Posted by gouravjee

hello,
i am trying to convert a mesh from fluent to openfoam and i have followed a procedure as shown below ::
dos2unix mesh.msh (i am using fluent in windows)
fluentMeshToFoam mesh.msh but ige the following error but i get the following error.

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

number of faces: 3486
Number of points: 614
Reading uniform faces
Reading uniform faces
Reading points


FINISHED LEXING


#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigSegv::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3  ? at ??:?
#4  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#5  ? at ??:?
Segmentation fault (core dumped)

can you tell me where might be the problem??

Hey there,

Just a quick hint before somebody answers your question. It is a good practice to use fluent3DMeshToFoam instead of the version you are using currently. You also might find more information here http://openfoamwiki.net/index.php/Fluent3DMeshToFoam

[Lighto]

Hi everyone,

I have a similar problem in OpenFOAM while running a modified pimplefoam with lagrangian particle tracking. The solver is working in PISO mode, and after a few time steps it stops with FPE error.

I have time-varying boundary conditions (I specified them with a txt file) for velocity at the inlet and noSlip conditions at walls and zeroGradient at outlet, while for the pressure I imposed a zeroGradient conditions in all the patches and uniform 0 at the outlet.

I am also running in parallel on a machine with 16 cores.

The error log is the following:

Foam::error::printStack(Foam::Ostream&) at ??:?
[13] #1  Foam::sigFpe::sigHandler(int) at ??:?           
[13] #2  ? at sigaction.c:?m     
[13] #3  Foam::symGaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, long) at ??:?
[13] #4  Foam::symGaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, long) const at ??:?
[13] #5  Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[13] #6  Foam::fvMatrix<Foam::Vector<double> >::solveSegregated(Foam::dictionary const&) at ??:?
[13] #7  Foam::fvMatrix<Foam::Vector<double> >::solve(Foam::dictionary const&) at ??:?
[13] #8  ? at ??:?
[13] #9  ? at ??:?
[13] #10  __libc_start_main at ??:?
[13] #11  ? at ??:?
[nodo-b06:21640] *** Process received signal ***
[nodo-b06:21640] Signal: Floating point exception (8)
[nodo-b06:21640] Signal code:  (-6)

I think that the problem is the courant number, that at a certain time goes above 1... I have already refined the time step and the problem just shifted from 0.036 to 0.040 (of the "real" time of the phenomenon that I have to simulate). It took me a week of simulation, so I can't refine that much...

I attach here two file showing the Co values with respect to the inlet velocity values.