[Zinedine]

Dear foamers,

I have setup a case using OF2.3 using pimpleFoam with
k-epsilon turbulence and this works fine when running using OF2.3.1.

However when attempting to run the same case using OF2.4 - I have the following error:

####################################
Create time

Create mesh for time = 0

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model type RASModel
Selecting RAS turbulence model kEpsilon
bounding k, min: 0 max: 1e-05 average: 1e-05
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
sigmaEps 1.3;
}

No finite volume options present


PIMPLE: Operating solver in PISO mode


Starting time loop

Courant Number mean: 1.21676e-05 max: 0.000519965
deltaT = 1.1999e-05
Time = 1.1999e-05

PIMPLE: iteration 1
smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 4.39177e-09, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 5.70503e-08, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 7.67011e-08, No Iterations 1
GAMG: Solving for p, Initial residual = 1, Final residual = 0.00441416, No Iterations 7
time step continuity errors : sum local = 3.12544e-09, global = -5.85747e-11, cumulative = -5.85747e-11
GAMG: Solving for p, Initial residual = 0.559135, Final residual = 5.71717e-08, No Iterations 14
time step continuity errors : sum local = 8.4225e-13, global = 1.05287e-14, cumulative = -5.85642e-11
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib64/libc.so.6"
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4 void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:perator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#6 Foam::incompressible::RASModels::kEpsilon::correct () at ??:?
#7 ? at ??:?
#8 __libc_start_main in "/lib64/libc.so.6"
#9 ? at ??:?
Floating point exception
##################################################

This is the same case file -I have trying to find out what is the issue
with no success - I am sure this must be something silly.

The error statement in red: Floating point exception
somethig is not well define- that is weird as it is the same case
that runs with OF2.3.1 ....

Your help is really appreciated.

Thanks

Z.

[ppendota]

Hi, if you are using k=0 as a boundary condition, try using k=1e-10 (or a very small value) to avoid division by zero in the solver.

[Zinedine]

Good morning PremchandThanks for your suggestion - just implemented the change and
indeed it seems to function now.
Would you know why this issue doesn't show up on OF2.3.1 but does on OF2.4?

Regards

Zinedine