[CeliaSanchez]

I'm trying to simulate the heat transfer using chtMultiRegionTwoPhaseEulerFoam solver of a chamber half filled with liquid. For this I have adapt the tutorial heatTransfer/chtMultiRegionTwoPhaseEulerFoam/solidQuenching2D.
At the second time step I get the following error:

[stack trace]
=============
#1 Foam::sigFpe::sigHandler(int) in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#2 ? in /lib/x86_64-linux-gnu/libc.so.6
#3 ? in /lib/x86_64-linux-gnu/libm.so.6
#4 Foam:: pow025(Foam::Field<double>&, Foam::UList<double> const&) in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#5 Foam:: pow025(Foam::tmp<Foam::Field<double> > const&) in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#6 Foam::wallBoilingModels::CHFModels::HuaXu::CHFSubC ool(Foam:: phaseModel const&, Foam:: phaseModel const&, int, Foam::Field<double> const&, Foam::Field<double> const&, Foam::Field<double> const&) const in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libreactingMultiphaseSystem.so
#7 Foam::compressible::alphatWallBoilingWallFunctionF vPatchScalarField::updateCoeffs() in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libreactingMultiphaseSystem.so
#8 ? in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libreactingMultiphaseSystem.so
#9 ? in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libtwoPhaseReactingTurbulenceModels.so
#10 ? in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libtwoPhaseReactingTurbulenceModels.so
#11 Foam:: phaseSystem::correctTurbulence() in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libreactingMultiphaseSystem.so
#12 ? in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/bin/chtMultiRegionTwoPhaseEulerFoam
#13 ? in /lib/x86_64-linux-gnu/libc.so.6
#14 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#15 ? in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/bin/chtMultiRegionTwoPhaseEulerFoam
=============
Floating point exception (core dumped)

The boundary conditions that I'm using for alphat are attached.

I would be very grateful if someone who has successfully adapted the tutorial for a fluid region with liquid and gas can help me.