[gu1]

Hi,

I'm trying to use the kEpsilon model and it's giving me this error. Interestingly, for the mixtureKEpsilon model, the simulation runs fine.

Courant Number mean: 0.0375 max: 0.0375
Time = 0.005

PIMPLE: Iteration 1
MULES: Solving for alpha.air
MULES: Solving for alpha.water
air fraction, min, max = 0.293333 0 1
water fraction, min, max = 0.706667 0 1
Phase-sum volume fraction, min, max = 1 1 1
MULES: Solving for alpha.air
MULES: Solving for alpha.water
air fraction, min, max = 0.293333 0 1
water fraction, min, max = 0.706667 0 1
Phase-sum volume fraction, min, max = 1 1 1
Constructing momentum equations
smoothSolver:  Solving for e.air, Initial residual = 1, Final residual = 6.71348e-09, No Iterations 2
smoothSolver:  Solving for e.water, Initial residual = 1, Final residual = 3.81242e-09, No Iterations 4
air min/max T 300 - 350
water min/max T 300 - 350.04
GAMG:  Solving for p_rgh, Initial residual = 0.000188923, Final residual = 7.61464e-09, No Iterations 20
PIMPLE: Iteration 2
MULES: Solving for alpha.air
MULES: Solving for alpha.water
air fraction, min, max = 0.293327 0 1
water fraction, min, max = 0.706673 -1.12253e-16 1
Phase-sum volume fraction, min, max = 1 1 1
MULES: Solving for alpha.air
MULES: Solving for alpha.water
air fraction, min, max = 0.293329 0 1
water fraction, min, max = 0.706671 -8.69992e-18 1
Phase-sum volume fraction, min, max = 1 1 1
Constructing momentum equations
smoothSolver:  Solving for e.air, Initial residual = 0.0385992, Final residual = 3.57737e-09, No Iterations 3
smoothSolver:  Solving for e.water, Initial residual = 8.92925e-05, Final residual = 2.18232e-10, No Iterations 1
air min/max T 299.571 - 350
water min/max T 301.108 - 350
GAMG:  Solving for p_rgh, Initial residual = 4.16813e-05, Final residual = 1.06388e-09, No Iterations 13
PIMPLE: Iteration 3
MULES: Solving for alpha.air
MULES: Solving for alpha.water
air fraction, min, max = 0.293328 0 1
water fraction, min, max = 0.706672 0 1
Phase-sum volume fraction, min, max = 1 1 1
MULES: Solving for alpha.air
MULES: Solving for alpha.water
air fraction, min, max = 0.293331 0 1
water fraction, min, max = 0.706669 0 1
Phase-sum volume fraction, min, max = 1 1 1
Constructing momentum equations
smoothSolver:  Solving for e.air, Initial residual = 0.0273701, Final residual = 4.29279e-10, No Iterations 3
smoothSolver:  Solving for e.water, Initial residual = 4.02651e-06, Final residual = 8.14158e-12, No Iterations 1
air min/max T 299.65 - 350
water min/max T 299.571 - 350
GAMG:  Solving for p_rgh, Initial residual = 2.15398e-05, Final residual = 8.44514e-09, No Iterations 11
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::symGaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:?
#4  Foam::symGaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const at ??:?
#5  Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
#7  Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam"
#8  Foam::fvMatrix<double>::solve() in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam"
#9  Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<double> > const&) at ??:?
#10  Foam::RASModels::kEpsilon<Foam::PhaseCompressibleMomentumTransportModel<Foam::dynamicTransportModel> >::correct() at ??:?
#11  Foam::phaseSystem::correctTurbulence() at ??:?
#12  ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam"
#13  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#14  ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam"
Floating point exception (core dumped)

Anyone can help me?