CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Foam::error::PrintStack (https://www.cfd-online.com/Forums/openfoam-solving/89644-foam-error-printstack.html)

shadabdyn January 18, 2017 09:34

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.


Code:

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 January 20, 2017 01:41

error after adding source term
 
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::printStack(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 January 20, 2017 07:28

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.

https://www.cfd-online.com/Forums/op...-get-help.html

andreamc April 8, 2017 14:56

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::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam:: species::thermo<Foam::hConstThermo<Foam::perfectGa s<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() at ??:?
#4 Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam:: species::thermo<Foam::hConstThermo<Foam::perfectGa 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 April 18, 2017 00:50

Hello esujby

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

haydii April 26, 2017 13:21

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 April 26, 2017 15:29

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 April 28, 2017 04:46

dear esujby
 
Quote:

Originally Posted by esujby (Post 646733)
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 April 28, 2017 05:24

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 April 28, 2017 11:53

Quote:

Originally Posted by sheaker (Post 646931)
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 April 29, 2017 05:52

Hello Haydii,

to get the best help i recommend the following advice.

Quote:

Originally Posted by esujby (Post 626624)
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 April 29, 2017 08:17

Quote:

Originally Posted by esujby (Post 647024)
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

Quote:

Originally Posted by iamranjitsingh (Post 645257)
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 July 6, 2017 10:08

1 Attachment(s)
Quote:

Originally Posted by ebah6 (Post 345030)
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::operator++()
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::printStack(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 July 6, 2017 11:40

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 July 7, 2017 07:31

error still
 
5 Attachment(s)
Quote:

Originally Posted by esujby (Post 656065)
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 July 21, 2017 10:03

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 July 22, 2017 08:26

say thanks
 
Quote:

Originally Posted by esujby (Post 657942)
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 March 13, 2018 08:55

problem in fluentMeshToFoam mesh conversion
 
Quote:

Originally Posted by wyldckat (Post 312584)
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:
    Code:

    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:
    Code:

    #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:
    Code:

    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.

Code:


// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
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 March 13, 2018 09:50

Quote:

Originally Posted by gouravjee (Post 685012)
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.

Code:


// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
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 August 7, 2018 05:32

FPE error using pimpleFoam
 
2 Attachment(s)
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:
Code:

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.


All times are GMT -4. The time now is 23:30.