[SHY22]

Hi I encounter a problem about openfoam7 reactingTwoPhaseEulerFoam.In my simulation I want to simulate a Bartolomej experiment of suncooled flow boiling at 45bar using water ,but when i turn the thermophsyical.liquid from const to polynomial it show the erro:
4] #1 Foam::sigFpe::sigHandler(int) at ??:?
[4] #2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
[4] #3 ? in "/lib/x86_64-linux-gnu/libm.so.6"
[4] #4 powf64 in "/lib/x86_64-linux-gnu/libm.so.6"
[4] #5 Foam:ow(Foam::Field<double>&, Foam::UList<double> const&, double const&) at ??:?
[4] #6 Foam:ow(Foam::UList<double> const&, double const&) at ??:?
[4] #7 Foam::compressible::alphatPhaseChangeJayatillekeWa llFunctionFvPatchScalarField::Psmooth(Foam::Field< double> const&) const at ??:?
[4] #8 Foam::compressible::alphatWallBoilingWallFunctionF vPatchScalarField::updateCoeffs() at ??:?
[4] #9 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::Boundary::evaluate() in "/home/dyfluid/OpenFOAM/dyfluid-7/platforms/linux64GccDPInt32Opt/bin/reactingTwoPhaseEulerFoam_Modif1"
[4] #10 Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foa m::PhaseCompressibleTurbulenceModel<Foam:haseMod el> > >::correctNut() at ??:?
[4] #11 Foam:haseSystem::correctEnergyTransport() at ??:?
[4] #12 Foam::TwoResistanceHeatTransferPhaseSystem<Foam::M omentumTransferPhaseSystem<Foam::twoPhaseSystem> >::correctEnergyTransport() at ??:?
[4] #13 ? in "/home/dyfluid/OpenFOAM/dyfluid-7/platforms/linux64GccDPInt32Opt/bin/reactingTwoPhaseEulerFoam_Modif1"
[4] #14 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[4] #15 ? in "/home/dyfluid/OpenFOAM/dyfluid-7/platforms/linux64GccDPInt32Opt/bin/reactingTwoPhaseEulerFoam_Modif1"

my thermophsyical.liquid is :

const:
thermoType
{
type heRhoThermo;
mixture pureMixture;
transport const;
thermo hRefConst;
equationOfState rhoConst;
specie specie;
energy sensibleEnthalpy;
}

dpdt no;

mixture
{
specie
{
molWeight 18.02;
}
equationOfState
{
rho 787.9;
}
thermodynamics
{
Hf 0;
Cp 4922;
Tref 530.6;
Href 0;
}
transport
{
mu 1.07e-4;
Pr 0.8680132;
}
}

polynimial:
thermoType
{
type heRhoThermo;
mixture pureMixture;
transport polynomial;
thermo hPolynomial;
equationOfState icoPolynomial;
specie specie;
energy sensibleEnthalpy;
}

dpdt no;//no;



mixture
{
specie
{
nMoles 1;
molWeight 18.0153;
}
equationOfState // rho = rho(T)
{
rhoCoeffs<8> ( 539.1576625 2.541187391 -0.003906805 0 0 0 0 0 ); // rho = a + bT
}
thermodynamics
{
Sf 9.103679;//2.8615;//Standard entropy
Hf 0;//Heat of formation
Pref 0.001e6;
Tref 300;//530.6;//
Href 2551.013479;//1122200;//
CpCoeffs<8> ( 16967.87383 -57.32039084 0.065357618 0 0 0 0 0 ); // Specific heat at constant pressure polynomial coeffs
}
transport
{
muCoeffs<8> ( 0.000951045 -2.75795E-06 2.18533E-09 0 0 0 0 0 ); // mu = a + bT
kappaCoeffs<8> ( -0.081242398 0.003807896 -4.72596E-06 0 0 0 0 0); // kappa = a + bT
PrCoeffs<8> ( 7.909452932 -0.026764627 2.53238E-05 0 0 0 0 0); //Pr 0.8304;
}
}

I cannot find any solution about the erro?Could anyone give me some advice ,please?
Best whishes.