[Linea]

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:

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).
--------------------------------------------------------------------------

[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.

[Linea]

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.

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 ;
}