CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM (https://www.cfd-online.com/Forums/openfoam/)
-   -   Report "negative initial temperature T0" when using compressibleInterFoam (https://www.cfd-online.com/Forums/openfoam/207550-report-negative-initial-temperature-t0-when-using-compressibleinterfoam.html)

Zhanyuan Wang September 28, 2018 23:22

Report "negative initial temperature T0" when using compressibleInterFoam
 
1 Attachment(s)
Hello,
I m a new learner of OF and I m trying to simulating supercritical (high-temperature and high-pressure) spray characteristics using compressibleInterFoam. When i run the case, it always broke down after only a few steps and report like that:

[0]
[4] [1]
[4]
[4] --> FOAM FATAL ERROR:
[4] Negative initial temperature T0: -1496.0858
[4]
[4] 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::hPolynomialThermo<Foam::icoPolynomial<Foam:: specie> >; Type = Foam::sensibleInternalEnergy; Foam::scalar = double; Foam::species::thermo<Thermo, Type> = Foam::species::thermo<Foam::hPolynomialThermo<Foam ::icoPolynomial<Foam::specie> >, Foam::sensibleInternalEnergy>]
[4] in file /home/wzy/OpenFOAM/OpenFOAM-5.x/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 54.

Because i m a new learner, i just simply copy most part of the tutorial case and modified some details. The tutorial case uses laminar model, const thermo type, hConst transport type and perfectFuild model and i changed them to LES model, polynomial, hPolynomial and icoPolynomial. I dont know if it caused the problem. Did i set wrong BC values, wrong fvSchemes or wrong time steps? Is there any problem with my mesh? Need i modify the source code of this solver?

My case is in the attachment and I will be grateful very very very much for anyone who can offer help.

The King November 14, 2018 03:58

Hi!

did you found a solution? Got the same...

The King

The King November 14, 2018 04:13

Hi,

Found the answwer in on this forum. I changed the pressure from 0 (guage) to 1atm (1e5) and the p_rgh file since the case is compressible and than you need absolute pressures in the calculations.

Simulation running without any problem now.

The King


All times are GMT -4. The time now is 17:19.