tutorial case for simpleReactingParcelFoam crashes

Summary

Running the tutorials/lagrangian/simpleReactingParcelFoam/verticalChannel case out-of-the-box, the simulation crashes after 41 iterations, complaining of a negative temperature:

--> FOAM FATAL ERROR: (openfoam-2112) Negative initial temperature T0: -3.334796867

Steps to reproduce

Run the tutorial case.

Environment information

OpenFOAM version : v2112 Operating system : ubuntu Hardware info : n/a Compiler : n/a

Possible fixes

A simple fix is to include a temperature limiter, using the following constant/fvOptions file (the header has been omitted for the sake of brevity):

limitT { type limitTemperature; min 101; max 600; selectionMode all; }

This helps the solution get over the initial low-T spike, and once past this the temperature limiter has no impact (i.e. it is not constraining the final solution).