[chengtun]

I am trying to modify the case of wallBoiling which is OpenFOAM-dev/tutorials/multiphase/reactingTwoPhaseEulerFoam/RAS/wallBoiling
I only changed the fluid properties and geometry parameters and unfortunately got the bug of Negative initial temperature T0 after some iterative steps.
The initial temperature is 303.15K for both gas and liquid. I tried heat flux q in T.liquid for both zero and several thousand. They all ended in bug of Negative initial temperature. More details are as follows.
--> FOAM FATAL ERROR:
Negative initial temperature T0: -35.6173826

From function Foam::scalar Foam::species::thermo<Thermo, Type>::T(Foam::scalar, Foam::scalar, Foam::scalar, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar) const) const [with Thermo = Foam::hRefConstThermo<Foam::rhoConst<Foam::specie> >; Type = Foam::sensibleEnthalpy; Foam::scalar = double; Foam::species::thermo<Thermo, Type> = Foam::species::thermo<Foam::hRefConstThermo<Foam:: rhoConst<Foam::specie> >, Foam::sensibleEnthalpy>]
in file /home/cpx/OpenFOAM/cpx-5.0/OpenFOAM-dev/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 54.

FOAM aborting

#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hRefConstThermo<Foam::rhoConst<Fo am::specie> >, Foam::sensibleEnthalpy> > > >::calculate() at ??:?
#3 Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hRefConstThermo<Foam::rhoConst<Fo am::specie> >, Foam::sensibleEnthalpy> > > >::correct() at ??:?
#4 Foam::phaseSystem::correctThermo() at ??:?
#5 Foam::ThermalPhaseChangePhaseSystem<Foam::HeatAndM assTransferPhaseSystem<Foam::MomentumTransferPhase System<Foam::twoPhaseSystem> > >::correctThermo() at ??:?
#6 ? at ??:?
#7 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#8 ? at ??:?
Aborted (core dumped)

I tried to set temperature constraints in fvOptions for gas and liquid. The following example is for liquid. But it doesn't work.
temperature_constraints
{
type limitTemperature;
selectionMode all;
active true;
min 273;
max 600;
phase liquid;
}
I also tried to delete the judgment for negative temperature in OpenFOAM-dev/src/thermophysicalModels/specie/thermo/thermo/thermoI.H line 52-57
if (T0 < 0)
{
FatalErrorInFunction
<< "Negative initial temperature T0: " << T0
<< abort(FatalError);
}

But it still shows the same error. I am really confused why this bug appear and why the deletion doesn't work. When I open the file /OpenFOAM-dev/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 54, the judgment has already been deleted.
Could anyone give me a hint? Thanks!