[JeanMichel]

Hello,
I'm trying to simulate the air of the air with a hot plate which is in a box that simulate the atmosphere. But whenever the timestep, openfoam aborts after a few time steps with the explanation that the inital temperature is negative.

I have spent a lot of time but I can't find if it's linked to the mesh, to the boundary conditions or to something else.
You can download my case here: https://epolytechnique-my.sharepoint...m9QIQ?e=rjCm4y

Thanks for your help because I'm starting to desperate !

Courant Number mean: 1.22697e-05 max: 2.55618
deltaT = 1.44231e-05
Time = 0.0318686s

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab:  Solving for h, Initial residual = 0.00344084, Final residual = 1.02413e-09, No Iterations 2
GAMG:  Solving for p_rgh, Initial residual = 0.0375017, Final residual = 5.17781e-08, No Iterations 11
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.19373e-10, global = 1.12428e-10, cumulative = -2.2971e-10
ExecutionTime = 45.1982 s  ClockTime = 47 s

Courant Number mean: 1.652e-05 max: 11.005
deltaT = 3.93176e-06
Time = 0.0318726s

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab:  Solving for h, Initial residual = 0.0108358, Final residual = 2.32118e-09, No Iterations 2


--> FOAM FATAL ERROR:
Negative initial temperature T0: -19.9927

    From function static Foam::scalar Foam::species::thermo<Thermo, Type>::T(const ThermoType&, Foam::scalar, Foam::scalar, Foam::scalar, FType, dFdTType, LimitType, bool) [with ThermoType = Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy>; FType = double (Foam::hConstThermo<Foam::perfectGas<Foam::specie> >::*)(double, double) const; dFdTType = double (Foam::hConstThermo<Foam::perfectGas<Foam::specie> >::*)(double, double) const; LimitType = double (Foam::hConstThermo<Foam::perfectGas<Foam::specie> >::*)(double) const; Thermo = Foam::hConstThermo<Foam::perfectGas<Foam::specie> >; Type = Foam::sensibleEnthalpy; Foam::scalar = double]
    in file /home/ubuntu/OpenFOAM/OpenFOAM-11/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 300.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::error::abort() at ??:?
#2  double Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy>::T<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy>, double (Foam::hConstThermo<Foam::perfectGas<Foam::specie> >::*)(double, double) const, double (Foam::hConstThermo<Foam::perfectGas<Foam::specie> >::*)(double, double) const, double (Foam::hConstThermo<Foam::perfectGas<Foam::specie> >::*)(double) const>(Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> const&, double, double, double, double (Foam::hConstThermo<Foam::perfectGas<Foam::specie> >::*)(double, double) const, double (Foam::hConstThermo<Foam::perfectGas<Foam::specie> >::*)(double, double) const, double (Foam::hConstThermo<Foam::perfectGas<Foam::specie> >::*)(double) const, bool) at ??:?
#3  Foam::heRhoThermo<Foam::rhoThermo::composite, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > >::calculate() at ??:?
#4  virtual thunk to Foam::heRhoThermo<Foam::rhoThermo::composite, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > >::correct() at ??:?
#5  Foam::solvers::fluid::thermophysicalPredictor() at ??:?
#6  ? in "/opt/openfoam11/platforms/linux64GccDPInt32Opt/bin/foamRun"
#7  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#8  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9  ? in "/opt/openfoam11/platforms/linux64GccDPInt32Opt/bin/foamRun"
Aborted