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