[shv.al]

Hi all, I am trying to simulate a dam break with sloping sea wall using the compressibleMultiphaseInterFoam solver. After running blockMesh and setFields, I get this error when compressibleMultiphaseInterFoam is running:

PIMPLE: iteration 1
MULES: Solving for alpha.water
water volume fraction, min, max = 0.151061 -2.18497e-34 1
MULES: Solving for alpha.air
air volume fraction, min, max = 0.848939 0 1
Phase-sum volume fraction, min, max = 1 1 1
MULES: Solving for alpha.water
water volume fraction, min, max = 0.151061 -5.92092e-30 1
MULES: Solving for alpha.air
air volume fraction, min, max = 0.848939 -4.04094e-43 1
Phase-sum volume fraction, min, max = 1 1 1
MULES: Solving for alpha.water
water volume fraction, min, max = 0.151061 -4.93459e-29 1
MULES: Solving for alpha.air
air volume fraction, min, max = 0.848939 -8.0895e-43 1
Phase-sum volume fraction, min, max = 1 1 1
MULES: Solving for alpha.water
water volume fraction, min, max = 0.151061 -1.69187e-28 1
MULES: Solving for alpha.air
air volume fraction, min, max = 0.848939 -1.21457e-42 1
Phase-sum volume fraction, min, max = 1 1 1
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for T, Initial residual = 8.19455e-05, Final residual = 6.16216e-09, No Iterations 1


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

From 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:erfectFluid<Foam::specie > >; Type = Foam::sensibleInternalEnergy; Foam::scalar = double; Foam::species::thermo<Thermo, Type> = Foam::species::thermo<Foam::hConstThermo<Foam:er fectFluid<Foam::specie> >, Foam::sensibleInternalEnergy>]
in file /home/pawan/OpenFOAM/OpenFOAM-v2006/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 56.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::error::exitOrAbort(int, bool) at ??:?
#2 Foam::heRhoThermo<Foam::rhoThermo, Foam:ureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hConstThermo<Foam:erfectFluid<F oam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, bool) at ??:?
#3 Foam::heRhoThermo<Foam::rhoThermo, Foam:ureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hConstThermo<Foam:erfectFluid<F oam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() at ??:?
#4 Foam:haseModel::correct() at ??:?
#5 Foam::multiphaseMixtureThermo::correct() at ??:?
#6 ? at ??:?
#7 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#8 ? at ??:?
Aborted (core dumped)

I've checked the boundary conditions and they all seem fine. Could use some help, thank you.

https://drive.google.com/file/d/15Im...ew?usp=sharing

Here's a link to download my case for more information.