Skip to content

bouyancyTurbSource with hotRoom tutorial fails

Summary

bouyancyTurbSource with hotRoom tutorial fails

Steps to reproduce

Using the hotRoom tutorial with the fvOptions file as below, the run crashed. Without the fvOptions file it worked perfectly.

constant/fvOptions Code:

buoyancyTurbSource1
{
    type            buoyancyTurbSource;
    beta            3.3e-03; // Thermal expansion coefficient for incompressible cases
    rho             rho;
    alphat          alphat;
    T               T;
    selectionMode    all;
}

Example case

hotRoom tutorial and fvOptions as described above.

What is the current bug behaviour?

First of all, everything seems to be fine, but at timestep 58 the solver crashed.

What is the expected correct behavior?

A simulation running ;-)

Relevant logs and/or images

LOG:
...
Selecting finite volume options type buoyancyTurbSource
Source: buoyancyTurbSource1
- selecting all cells
- selected 4000 cell(s) with volume 500
Applying buoyancyTurbSource to: epsilon and k
Courant Number mean: 0 max: 0
....

Courant Number mean: 1.22152e+25 max: 2.19919e+26
Time = 56

PIMPLE: iteration 1
DILUPBiCGStab: Solving for T, Initial residual = 1, Final residual = 1.34532e-07, No Iterations 2
DICPCG: Solving for p_rgh, Initial residual = 0.999025, Final residual = 0.00123228, No Iterations 1
time step continuity errors : sum local = 4.44339e+23, global = 6.50956e+06, cumulative = 6.89475e+09
DICPCG: Solving for p_rgh, Initial residual = 6.37374e-11, Final residual = 6.37374e-11, No Iterations 0
time step continuity errors : sum local = 3.33592e+24, global = 3.54335e+07, cumulative = 6.93018e+09
DILUPBiCGStab: Solving for epsilon, Initial residual = 0.0435253, Final residual = 7.14145e-07, No Iterations 58
bounding epsilon, min: -1.09797e+67 max: 5.40352e+68 average: 1.83296e+66
DILUPBiCGStab: Solving for k, Initial residual = 0.843783, Final residual = 4.15889e-07, No Iterations 35
bounding k, min: -2.3265e+59 max: 1.81813e+59 average: 1.87383e+58
ExecutionTime = 0.42 s ClockTime = 0 s

Courant Number mean: 1.28264e+25 max: 1.89296e+26
Time = 58

PIMPLE: iteration 1
DILUPBiCGStab: Solving for T, Initial residual = 1, Final residual = 3.74953e-07, No Iterations 2
DICPCG: Solving for p_rgh, Initial residual = 0.0478055, Final residual = 0.000377563, No Iterations 6
time step continuity errors : sum local = 8.32321e+30, global = -1.83838e+14, cumulative = -1.83831e+14
DICPCG: Solving for p_rgh, Initial residual = 9.1225e-72, Final residual = 9.1225e-72, No Iterations 0
time step continuity errors : sum local = 6.17026e+30, global = -3.54588e+14, cumulative = -5.3842e+14
DILUPBiCGStab: Solving for epsilon, Initial residual = 1.70957e-08, Final residual = 1.70957e-08, No Iterations 0
bounding epsilon, min: 5.06942e-23 max: 1.2042e+89 average: 1.62094e+88
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in /lib64/libpthread.so.0
#3 Foam::scalarProduct<double, double>::type Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4 Foam::PBiCGStab::scalarSolve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#5 Foam::PBiCGStab::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
#7 Foam::fvMatrix<double>::solveSegregatedOrCoupled(F oam::dictionary const&) at ??:?
#8 Foam::fvMesh::solve(Foam::fvMatrix<double>&, Foam::dictionary const&) const at ??:?
#9 Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<doubl e> > const&) at ??:?
#10 Foam::RASModels::kEpsilon<Foam::IncompressibleTurb ulenceModel<Foam::transportModel> >::correct() at ??:?
#11 ? at ??:?
#12 __libc_start_main in /lib64/libc.so.6
#13 ? at ??:?

Environment information

  • OpenFOAM version : v2106, branch=master, build=f815a12b-20210902
  • Operating system : CentOS Linux Version 8
  • Hardware info : Cluster
  • Compiler : gcc

Possible fixes

Edited by Simon Holz