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

printStack(Foam::Ostream& - Debug

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 15, 2020, 10:18
Cool printStack(Foam::Ostream& - Debug
  #1
New Member
 
Marlen
Join Date: Dec 2019
Posts: 4
Rep Power: 6
Linea is on a distinguished road
Hello,


I am using buoyantBoussinesqPimple trying to simulate natural ventilation in a generic room. I have been trying all different kind of settings but my simulation keeps on crashing with the following Error (see CODE). The iteration just stops and the logfile shows following:

Code:
PIMPLE: iteration 1
DILUPBiCGStab:  Solving for T, Initial residual = 0.9991967, Final residual = 1.533909e-06, No Iterations 5
DICPCG:  Solving for p_rgh, Initial residual = 0.8879945, Final residual = 4.709069e+12, No Iterations 1001
time step continuity errors : sum local = 1.455111e+36, global = -3.750831e+21, cumulative = -3.750831e+21
DICPCG:  Solving for p_rgh, Initial residual = 2.059018e-12, Final residual = 2.059018e-12, No Iterations 0
time step continuity errors : sum local = 6.534109e+68, global = 3.122948e+53, cumulative = 3.122948e+53
DILUPBiCGStab:  Solving for epsilon, Initial residual = 0.0640009, Final residual = 2.124775e-10, No Iterations 1
bounding epsilon, min: -5.432239e+53 max: 1.125331e+116 average: 2.376524e+110
DILUPBiCGStab:  Solving for k, Initial residual = 8.011081e-38, Final residual = 8.011081e-38, No Iterations 0
ExecutionTime = 2110.64 s  ClockTime = 2117 s

Courant Number mean: 3.267282e+68 max: 1.096909e+75
Time = 608.75

PIMPLE: iteration 1
DILUPBiCGStab:  Solving for T, Initial residual = 8.475585e-06, Final residual = 8.475585e-06, No Iterations 0
DICPCG:  Solving for p_rgh, Initial residual = 1, Final residual = 0.009356768, No Iterations 11
time step continuity errors : sum local = 1.135053e+124, global = 2.231873e+123, cumulative = 2.231873e+123
DICPCG:  Solving for p_rgh, Initial residual = 0.9999535, Final residual = 0.006464629, No Iterations 1
time step continuity errors : sum local = 1.192482e+140, global = -1.145045e+127, cumulative = -1.144822e+127
--------------------------------------------------------------------------
A process has executed an operation involving a call to the
"fork()" system call to create a child process.  Open MPI is currently
operating in a condition that could result in memory corruption or
other system errors; your job may hang, crash, or produce silent
data corruption.  The use of fork() (or system() or other calls that
create child processes) is strongly discouraged.

The process that invoked fork was:

  Local host:          [[863,1],1] (PID 7251)

If you are *absolutely sure* that your application will successfully
and correctly survive a call to fork(), you may disable this warning
by setting the mpi_warn_on_fork MCA parameter to 0.
--------------------------------------------------------------------------
--------------------------------------------------------------------------
mpirun noticed that process rank 3 with PID 0 on node n-62-29-24 exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------


Code:
[1] #0  [3] Foam::error::printStack(Foam::Ostream&)#0  Foam::error::printStack(Foam::Ostream&) at ??:?
[1] #1  Foam::sigFpe::sigHandler(int) at ??:?
[3] #1  Foam::sigFpe::sigHandler(int) at ??:?
[1] #2  ? at ??:?
[3] #2  ? in "/lib64/libc.so.6"
[1] #3  Foam::multiply(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/lib64/libc.so.6"
[3] #3  Foam::multiply(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
[1] #4  Foam::tmp<Foam::DimensionedField<double, Foam::volMesh> > Foam::operator*<Foam::volMesh>(Foam::tmp<Foam::DimensionedField<double, Foam::volMesh> > const&, Foam::DimensionedField<double, Foam::volMesh> const&) at ??:?
[3] #4  Foam::tmp<Foam::DimensionedField<double, Foam::volMesh> > Foam::operator*<Foam::volMesh>(Foam::tmp<Foam::DimensionedField<double, Foam::volMesh> > const&, Foam::DimensionedField<double, Foam::volMesh> const&) at ??:?
[1] #5  Foam::RASModels::kEpsilon<Foam::IncompressibleTurbulenceModel<Foam::transportModel> >::correct() at ??:?
[3] #5  Foam::RASModels::kEpsilon<Foam::IncompressibleTurbulenceModel<Foam::transportModel> >::correct() at ??:?
[1] #6   at ??:?
[3] #6  ?? at ??:?
[1] #7  __libc_start_main at ??:?
[3] #7  __libc_start_main in "/lib64/libc.so.6"
[1] #8   in "/lib64/libc.so.6"
[3] #8  ?? at ??:?
[n-62-29-24:07251] *** Process received signal ***
[n-62-29-24:07251] Signal: Floating point exception (8)
[n-62-29-24:07251] Signal code:  (-6)
[n-62-29-24:07251] Failing at address: 0x3863a00001c53
[n-62-29-24:07251] [ 0] /lib64/libc.so.6(+0x363f0)[0x7fe7d4ad63f0]
[n-62-29-24:07251] [ 1] /lib64/libc.so.6(gsignal+0x37)[0x7fe7d4ad6377]
[n-62-29-24:07251] [ 2] /lib64/libc.so.6(+0x363f0)[0x7fe7d4ad63f0]
[n-62-29-24:07251] [ 3] /appl/OpenFOAM/5.0/XeonE5-2660v3/gcc-6.3.0//OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam8multiplyERNS_5FieldIdEERKNS_5UListIdEES6_+0x1c5)[0x7fe7d5fbbcc5]
[n-62-29-24:07251] [ 4] /appl/OpenFOAM/5.0/XeonE5-2660v3/gcc-6.3.0//OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(_ZN4FoammlINS_7volMeshEEENS_3tmpINS_16DimensionedFieldIdT_EEEERKS6_RKS5_+0x16a)[0x7fe7d9badefa]
[n-62-29-24:07251] [ 5] /appl/OpenFOAM/5.0/XeonE5-2660v3/gcc-6.3.0//OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(_ZN4Foam9RASModels8kEpsilonINS_29IncompressibleTurbulenceModelINS_14transportModelEEEE7correctEv+0x437)[0x7fe7d9c1cef7]
[n-62-29-24:07251] [ 6] buoyantBoussinesqPimpleFoam[0x433ce3]
[n-62-29-24:07251] [ 7] /lib64/libc.so.6(__libc_start_main+0xf5)[0x7fe7d4ac2545]
[n-62-29-24:07251] [ 8] buoyantBoussinesqPimpleFoam[0x4381b1]
[n-62-29-24:07251] *** End of error message ***
 at ??:?
[n-62-29-24:07253] *** Process received signal ***
[n-62-29-24:07253] Signal: Floating point exception (8)
[n-62-29-24:07253] Signal code:  (-6)
[n-62-29-24:07253] Failing at address: 0x3863a00001c55
[n-62-29-24:07253] [ 0] /lib64/libc.so.6(+0x363f0)[0x7eff6144c3f0]
[n-62-29-24:07253] [ 1] /lib64/libc.so.6(gsignal+0x37)[0x7eff6144c377]
[n-62-29-24:07253] [ 2] /lib64/libc.so.6(+0x363f0)[0x7eff6144c3f0]
[n-62-29-24:07253] [ 3] /appl/OpenFOAM/5.0/XeonE5-2660v3/gcc-6.3.0//OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam8multiplyERNS_5FieldIdEERKNS_5UListIdEES6_+0x1c5)[0x7eff62931cc5]
[n-62-29-24:07253] [ 4] /appl/OpenFOAM/5.0/XeonE5-2660v3/gcc-6.3.0//OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(_ZN4FoammlINS_7volMeshEEENS_3tmpINS_16DimensionedFieldIdT_EEEERKS6_RKS5_+0x16a)[0x7eff66523efa]
[n-62-29-24:07253] [ 5] /appl/OpenFOAM/5.0/XeonE5-2660v3/gcc-6.3.0//OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(_ZN4Foam9RASModels8kEpsilonINS_29IncompressibleTurbulenceModelINS_14transportModelEEEE7correctEv+0x437)[0x7eff66592ef7]
[n-62-29-24:07253] [ 6] buoyantBoussinesqPimpleFoam[0x433ce3]
[n-62-29-24:07253] [ 7] /lib64/libc.so.6(__libc_start_main+0xf5)[0x7eff61438545]
[n-62-29-24:07253] [ 8] buoyantBoussinesqPimpleFoam[0x4381b1]
[n-62-29-24:07253] *** End of error message ***
[n-62-29-24:07245] 1 more process has sent help message help-opal-runtime.txt / opal_init:warn-fork
[n-62-29-24:07245] Set MCA parameter "orte_base_help_aggregate" to 0 to see all help / error messages

I've attached 0/U and 0/T, as well as the controldDict. I haven't made any other changes than velocity and temperature.

Let me know if more information are needed. I really tried to google a lot but nothing seemed to work.
Attached Files
File Type: txt U.txt (1.6 KB, 1 views)
File Type: txt T.txt (1.2 KB, 0 views)
File Type: txt controlDict.txt (1.0 KB, 0 views)
Linea is offline   Reply With Quote

Old   March 18, 2020, 08:59
Default Update: High Courant Number & 1001 Iterations
  #2
New Member
 
Marlen
Join Date: Dec 2019
Posts: 4
Rep Power: 6
Linea is on a distinguished road
Just a little update:

I have realised that my courant number gets really high just before the simulation crashes and that in the iteration for p_rgh it does not converge before reaching the limit of 1000 iterations.

Does anyone know how to fix it or what would need to be changed in order to stop it from crashing?

I have tried to increase the timesteps, but it does not seem to help.

I read I might need to change my schemes but I am unsure how. I posted my fvSchemes file.


Code:
FoamFile
{
version     2.0;
format      ascii;
class       dictionary;
object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
ddtSchemes { default Euler; } 
gradSchemes { default Gauss linear; } 
divSchemes
{
    default Gauss upwind;
    div(phi,U) Gauss upwind;
    div(phi,k) Gauss upwind;
    div(phi,epsilon) Gauss upwind;
    div(phi,omega) Gauss upwind;
    div((nuEff*dev2(T(grad(U))))) Gauss linear;
    div((nuEff*dev2(grad(U).T()))) Gauss linear;
    div(phi,h) Gauss upwind;
    div(R) Gauss linear;
    div((muEff*dev2(grad(U).T()))) Gauss linear;
    div((muEff*dev2(T(grad(U))))) Gauss linear;
    div(phi,R) Gauss limitedLinear 1;
    div(phi,nuTilda) Gauss limitedLinear 1;
    div(phiU,p) Gauss upwind;
    div((nuEff*dev(grad(U).T()))) Gauss linear;
    div((nuEff*dev(T(grad(U))))) Gauss linear;
    div(phi,T) Gauss upwind;
} 
laplacianSchemes
{
    default Gauss linear corrected;
    laplacian(1,p) Gauss linear corrected;
} 
interpolationSchemes { default linear; } 
snGradSchemes { default uncorrected; } 
fluxRequired
{
    default no;
    p_rgh ;
}
Linea is offline   Reply With Quote

Reply

Tags
buoyantboussinesqpimple, debug, fork, sigfpe, ventilation

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
[General] Extracting ParaView Data into Python Arrays Jeffzda ParaView 30 November 6, 2023 21:00
Finding the corresponding debug switch Tobi OpenFOAM Programming & Development 8 October 14, 2020 09:23
create the file *.foam phongstar OpenFOAM 12 October 14, 2018 18:06
Recursive output of Foam::word in gdb print command chengdi OpenFOAM Programming & Development 1 July 19, 2017 20:19
Eclipse - case debug error Bufacchi OpenFOAM 1 February 7, 2012 15:15


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