CFD Online Discussion Forums - How to avoid negative T0 in rhoCentralFom

Hello,

As far as I see a lot of people have this negative initial T0 error with rhoCentralFoam but I cannot find any clear answer on how to get rid of it.

I am trying to run a compressible case, where I expect a shock wave.

Adjustable time step does not work with sonicFoam.

Limiting Temperature or Pressure does not work with rhoCentralFoam.

A combination of 2 works with rhoPimpleFoam but then I cannot solve for the shock wave.

What would be the way to prevent negative temperatures in rhoCentralFoam ? It is not fvOptions, it cannot be set in fvSolutions PIMPLE like for pimple solvers.

The error message suggests it is something in the thermophysicalProperties file but there is nothing extraordinary in that file.

Code:



--> FOAM FATAL ERROR:

Negative initial temperature T0: -5374.53



    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::hConstThermo<Foam::perfectGas<Foam::specie> >; Type = Foam::sensibleInternalEnergy; Foam::scalar = double; Foam::species::thermo<Thermo, Type> = Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy>]

    in file X:/OpenFOAM-5.x/src/thermophysicalModels/specie/thermo/thermo/thermoI.H at line 54.