[Nomis]

Dear Foamers,

today I tested the new buoyancyTurbSource feature using the OF v2106 master branch.
However, using the hotRoom tutorial with the fvOptions file as below, the run crashed.
Without the fvOptions file it worked perfectly.

constant/fvOptions

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

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 ??:?

First of all, everything seems to be fine, but at timestep 58 the solver crashed.
Does anybody have an idea or experience with using buoyancyTurbSource?
Greetings Simon