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

Foam::error::PrintStack

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

Like Tree60Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 7, 2018, 11:18
Smile
  #81
Member
 
Join Date: Oct 2017
Posts: 52
Rep Power: 4
gouravjee is on a distinguished road
Quote:
Originally Posted by Lighto View Post
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.
Thanks for helping me
gouravjee is offline   Reply With Quote

Old   August 7, 2018, 11:21
Default
  #82
Member
 
Join Date: Oct 2017
Posts: 52
Rep Power: 4
gouravjee is on a distinguished road
Quote:
Originally Posted by Lighto View Post
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.
Refer to this link for identifying error:
https://en.wikipedia.org/wiki/Signal_(IPC)#SIGFPE
gouravjee is offline   Reply With Quote

Old   September 2, 2018, 12:03
Default Foam::error::PrintStack
  #83
New Member
 
ceren's Avatar
 
ceren cilavdaroğlu
Join Date: Aug 2018
Location: Turkey
Posts: 6
Rep Power: 3
ceren is on a distinguished road
Hello all,

I am having a similar error in using weirOwerflow
Below is the log folder:


Code:
Courant Number mean: 5.66695e+09 max: 3.11881e+15
Interface Courant Number mean: 0.0013693 max: 334.573
deltaT = 3.76816e-51
Time = 0.615178

PIMPLE: iteration 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 1  Min(alpha.water) = 0.656874  Max(alpha.water) = 1.29113
MULES: Solving for alpha.water
Phase-1 volume fraction = 1  Min(alpha.water) = 0.656874  Max(alpha.water) = 1.36343
DICPCG:  Solving for p_rgh, Initial residual = 0.991688, Final residual = 0.0386321, No Iterations 2
time step continuity errors : sum local = 3.17812e+17, global = 1.31038e+17, cumulative = 1.31038e+17
DICPCG:  Solving for p_rgh, Initial residual = 9.85132e-44, Final residual = 9.85132e-44, No Iterations 0
time step continuity errors : sum local = 3.17915e+17, global = 1.31083e+17, cumulative = 2.62121e+17
DICPCG:  Solving for p_rgh, Initial residual = 9.8524e-44, Final residual = 9.8524e-44, No Iterations 0
time step continuity errors : sum local = 3.1795e+17, global = 1.31132e+17, cumulative = 3.93252e+17
--------------------------------------------------------------------------
mpirun noticed that process rank 0 with PID 10107 on node ceren-CASPER-BILGISAYAR-SISTEMLERI-A-S exited on signal 8 (Floating point exception).

i am the beginner. please help me..


Thanks in advance for your reply
Ceren
ceren is offline   Reply With Quote

Old   September 7, 2018, 11:16
Default
  #84
Senior Member
 
Deep
Join Date: Oct 2017
Posts: 180
Rep Power: 4
deepbandivadekar is on a distinguished road
Quote:
Originally Posted by ceren View Post
Hello all,

I am having a similar error in using weirOwerflow
Below is the log folder:


Code:
Courant Number mean: 5.66695e+09 max: 3.11881e+15
Interface Courant Number mean: 0.0013693 max: 334.573
deltaT = 3.76816e-51
Time = 0.615178

PIMPLE: iteration 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 1  Min(alpha.water) = 0.656874  Max(alpha.water) = 1.29113
MULES: Solving for alpha.water
Phase-1 volume fraction = 1  Min(alpha.water) = 0.656874  Max(alpha.water) = 1.36343
DICPCG:  Solving for p_rgh, Initial residual = 0.991688, Final residual = 0.0386321, No Iterations 2
time step continuity errors : sum local = 3.17812e+17, global = 1.31038e+17, cumulative = 1.31038e+17
DICPCG:  Solving for p_rgh, Initial residual = 9.85132e-44, Final residual = 9.85132e-44, No Iterations 0
time step continuity errors : sum local = 3.17915e+17, global = 1.31083e+17, cumulative = 2.62121e+17
DICPCG:  Solving for p_rgh, Initial residual = 9.8524e-44, Final residual = 9.8524e-44, No Iterations 0
time step continuity errors : sum local = 3.1795e+17, global = 1.31132e+17, cumulative = 3.93252e+17
--------------------------------------------------------------------------
mpirun noticed that process rank 0 with PID 10107 on node ceren-CASPER-BILGISAYAR-SISTEMLERI-A-S exited on signal 8 (Floating point exception).
i am the beginner. please help me..


Thanks in advance for your reply
Ceren

Your Courant Numbers (entries in bold) and timestep continuity errors are huge. This suggests:
  1. Your time step for the cell sizes are not compatible. Try reducing the timestep.
  2. There's something that's causing imbalance in the system and increasing the continuity error. Take hard look at the boundary conditions again. See if all patches are correctly defined.
deepbandivadekar is offline   Reply With Quote

Old   September 18, 2019, 06:53
Default buoyantSimpleFoam crash - turbulence related
  #85
Member
 
Gareth
Join Date: Jun 2010
Posts: 48
Rep Power: 11
bullmut is on a distinguished road
Hi All


this is a great thread, old though it may be, and very informative.
I will start with my problem. I am trying to run buoyantSimpleFoam on a engine test cell. I am not modelling flow in the engine. Instead I have no mesh where the engine sis. Instead i have inlet and outlet boundaries for the engine inlet and outlet. I am trying to see re-circulation within the test cell (room).

I started out the simulation with just simpleFoam, to see if a cold flow run would work,. yeah success.
Then i added turbulence to the mix and also success.
The final step was adding heat, so i switched to buoyantSimplefoma as the temperature difference between the room and engine exhaust was greater than 10 degrees. Below you can see the last section of my log file. I am getting a floating point exception. So there are some odd maths going on.

Code:
Time = 5

DILUPBiCGStab:  Solving for Ux, Initial residual = 0.181033, Final residual = 3.19602e-05, No Iterations 1
DILUPBiCGStab:  Solving for Uy, Initial residual = 0.236924, Final residual = 6.02134e-05, No Iterations 1
DILUPBiCGStab:  Solving for Uz, Initial residual = 0.210526, Final residual = 3.29874e-05, No Iterations 1
DILUPBiCGStab:  Solving for h, Initial residual = 0.0253854, Final residual = 6.01331e-05, No Iterations 2
#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::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 at ??:?
#4  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::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
#5  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
#7  Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#8  Foam::fvMatrix<double>::solve() at ??:?
#9  ? at ??:?
#10  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11  ? at ??:?
Floating point exception (core dumped)
My first stop was checkMesh which said i had a good mesh (created with snappy).

After reading through this thread i followed the crumbs and saw my epsilon was going negative. And this knocked on to how my p_rgh was being calculated. So i now assume i have the wrong boundary condition epsilon and k. I have attached the patch summary and the log file of the simulation run. If anyone could please advise on what i am supposed to use for boundary conditions i would be very grateful.
For reference my patches are named:
Quote:
Box_Wall, Box_rig - walls
Box_inlet - air inlet from atmosphere
Box_outlet - air outlet, suction based with volume flow rate
Box_port - inlet/outlet not sure - exposed to atmosphere
Box_t_in - engine inlet, modeled as suction outlet
box_t_out - engine outlet,
Box_exhaust - outlet exposed to atmosphere,

Look fwd to your help
Attached Files
File Type: txt patch.txt (3.1 KB, 0 views)
File Type: txt log.txt (10.5 KB, 2 views)
bullmut is offline   Reply With Quote

Old   September 27, 2019, 14:11
Default
  #86
New Member
 
Join Date: Sep 2019
Posts: 5
Rep Power: 2
ASurendran is on a distinguished road
Hi Gareth,


I am also experiencing a similar error when I try to run buoyantSimpleFoam. I am very new to OpenFOAM. My problem is that of a heat exchanger (cold cylinder) in hot flow. For the moment, I am only looking at a 2D problem. In order to familiarise myself with OpenFOAM, I first tried to work on the cold flow case using simpleFoam and it works fine. The moment I add heat transfer, it goes bonkers and I am not able to understand the error.


From what Bruno mentioned in the beginning of the thread, all I could gather was that somewhere in my code there is some nasty math happening. But I am not able to figure out where. I am also not very confident with my BCs. I thought flow over a cylinder with heat transfer is sort of a standard problem that many people must have worked on. Unfortunately, the countless hours spend on the internet trying to figure out the BCs didn't help me much. Eventually, I decided to follow the example of CicuitBoardCooling in the tutorials. But, still going nowhere.


If somebody wants to see how my domain looks like, it is nothing but a half cylinder with symmetry conditions on the top and bottom wall. There is a velocity inlet and the pressure is assumed to be constant. The flow is hot and the tube is cold with a fixed temperature prescribed to it. I would have included the schematic if I knew how to include as attachments. Hope somebody can help me with that as well.


If somebody is willing to have a look at my codes, do let me know, I can post them here.


Cheers!
ASurendran is offline   Reply With Quote

Old   September 27, 2019, 14:29
Default
  #87
New Member
 
Rsingh
Join Date: Mar 2016
Posts: 3
Rep Power: 6
iamranjitsingh is on a distinguished road
Quote:
Originally Posted by ASurendran View Post
Hi Gareth,


I am also experiencing a similar error when I try to run buoyantSimpleFoam. I am very new to OpenFOAM. My problem is that of a heat exchanger (cold cylinder) in hot flow. For the moment, I am only looking at a 2D problem. In order to familiarise myself with OpenFOAM, I first tried to work on the cold flow case using simpleFoam and it works fine. The moment I add heat transfer, it goes bonkers and I am not able to understand the error.


From what Bruno mentioned in the beginning of the thread, all I could gather was that somewhere in my code there is some nasty math happening. But I am not able to figure out where. I am also not very confident with my BCs. I thought flow over a cylinder with heat transfer is sort of a standard problem that many people must have worked on. Unfortunately, the countless hours spend on the internet trying to figure out the BCs didn't help me much. Eventually, I decided to follow the example of CicuitBoardCooling in the tutorials. But, still going nowhere.


If somebody wants to see how my domain looks like, it is nothing but a half cylinder with symmetry conditions on the top and bottom wall. There is a velocity inlet and the pressure is assumed to be constant. The flow is hot and the tube is cold with a fixed temperature prescribed to it. I would have included the schematic if I knew how to include as attachments. Hope somebody can help me with that as well.


If somebody is willing to have a look at my codes, do let me know, I can post them here.


Cheers!
Dear ASurendran, Just check out the boundary conditions for pressure and velocity you have given. The solver you had mentioned here is default available in OpenFOAM, so there is no doubt regarding accuracy. the only thing here can be modified is pressure velocity, and temperature BC. Therefore have a look on in-built bouyantSimpleFOAM tutorial.
iamranjitsingh is offline   Reply With Quote

Old   September 27, 2019, 14:44
Default
  #88
New Member
 
Join Date: Sep 2019
Posts: 5
Rep Power: 2
ASurendran is on a distinguished road
Hi Ranjit,


Thank you for the reply. The closest tutorial that I can compare with my case is that of CircuitCooling. Shown below are my P, T and U codes. I am really doubtful of how to define p_rgh though, as this concept is a bit tricky for me to grasp. Is p_rgh really pressure? what value should we give? etc.


Code for P:
Quote:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [1 -1 -2 0 0 0 0]; // P
//dimensions [0 2 -2 0 0 0 0]; // P/rho

internalField uniform 1e5;//

boundaryField
{
inlet
{
type zeroGradient;
//value uniform 1e5;
}
outlet
{
type fixedValue; //
value uniform 1e5;// $internalField; //
}
top
{
type symmetry;
}
bottom
{
type symmetry;
}
cylinder
{
type zeroGradient;
}
defaultFaces
{
type empty;
}
}

// ************************************************** *********************** //

Code for T:


Quote:
/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 0 0 1 0 0 0];

internalField uniform 1500;

boundaryField
{
top
{
type symmetry;
}

bottom
{
type symmetry;
}

inlet
{
type fixedValue;
value uniform 1500;
}

outlet
{
//type zeroGradient;
type inletOutlet;
inletValue uniform 1500;
value uniform 1500;
}

cylinder
{
type fixedValue;
value uniform 300;
}

defaultFaces
{
type empty;
}


}

// ************************************************** *********************** //

Code for U:


Quote:
/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
location "0";
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 1 -1 0 0 0 0];

internalField uniform (21.0 0 0); // 0.627

boundaryField
{
inlet
{
type fixedValue;
value uniform (21.0 0 0);
}
outlet
{
type inletOutlet;
inletValue uniform (0 0 0);
value uniform (21.0 0 0);
}
top
{
type symmetry;
}
bottom
{
type symmetry;
}
cylinder
{
type noSlip;
//type fixedValue;
//value uniform (0 0 0);
}
defaultFaces
{
type empty;
}
}


// ************************************************** *********************** //

And finally Code for p_rgh:


Quote:
/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p_rgh;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [1 -1 -2 0 0 0 0]; // P
//dimensions [0 2 -2 0 0 0 0]; // P/rho

internalField uniform 1e5;

boundaryField
{
inlet
{
type zeroGradient;
//type fixedFluxPressure;
//gradient uniform 0;
//value uniform 1e5;
}
outlet
{
type fixedValue;
value uniform 1e5;
}
top
{
type symmetry;
}
bottom
{
type symmetry;
}
cylinder
{
type zeroGradient;
//type fixedFluxPressure;
//gradient uniform 0;
//value uniform 1e5;
}
defaultFaces
{
type empty;
}
}

// ************************************************** *********************** //

My simulation crashes after the 9th step. Don't know why.
ASurendran is offline   Reply With Quote

Old   September 30, 2019, 03:57
Default
  #89
Member
 
Gareth
Join Date: Jun 2010
Posts: 48
Rep Power: 11
bullmut is on a distinguished road
Quote:
Originally Posted by iamranjitsingh View Post
Dear ASurendran, Just check out the boundary conditions for pressure and velocity you have given. The solver you had mentioned here is default available in OpenFOAM, so there is no doubt regarding accuracy. the only thing here can be modified is pressure velocity, and temperature BC. Therefore have a look on in-built bouyantSimpleFOAM tutorial.

Hi ASurendran
In my case there were some poor cells causing me trouble. while checkMesh gave me the all clear, a personal inspection of the mesh showed me what to fix.
But i think the boundary conditions are usually a culprit. If you want to eliminate meshing as a possible source try a very simple setup with well defines cells...



regards
bullmut is offline   Reply With Quote

Old   September 30, 2019, 04:55
Default
  #90
New Member
 
Join Date: Sep 2019
Posts: 5
Rep Power: 2
ASurendran is on a distinguished road
Hi Gareth,


Thank you for your suggestion. I will check the mesh and BCs as suggested.


Kind regards
ASurendran is offline   Reply With Quote

Old   October 3, 2019, 07:01
Default
  #91
New Member
 
Join Date: Sep 2019
Posts: 5
Rep Power: 2
ASurendran is on a distinguished road
Hi All,


I forgot to mention that I work on High pressure flows with high temperature differences i.e., my domain pressure is of the order of 10bars and the difference in temperature between the flow and the heat exchanger is around 1100K. Has anyone encountered such situations?


To Gareth - I tried other possible BCs. Sadly, none of them worked.
ASurendran is offline   Reply With Quote

Reply

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
FoamerrorprintStack mayank OpenFOAM Running, Solving & CFD 38 November 25, 2011 22:58


All times are GMT -4. The time now is 08:34.