[Zhanyuan Wang]

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]

Hi!

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

The King

[The King]

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