[hannes]

Description:
While starting the Xoodles/pitzDaily3D tutorial, Xoodles fails and reports a floating point error.
The reason seems to be in XiFoam/bEqn.H in line 79. There, muu is divided by (rhou*epsilon). It seems, that epsilon may be zero in some cases.
Change the statement to:

volScalarField tauEta = sqrt(thermo->muu()/(
rhou*epsilon + dimensionedScalar("1e-6", rhou.dimensions()*epsilon.dimensions(), 1e-6)
));


Solver/Application:
Xoodles, XiFoam, engineFoam

Source file:
$FOAM_APP/solvers/combustion/XiFoam/bEqn.H

Testcase:
$FOAM_TUTORIALS/Xoodles/pitzDaily3D

Platform:
x86_64

Version:
1.4 (also 1.3)

Notes:

[henry]

Thanks for the bug report, the problem is that in the initial k field the internal and wall boundary values are set to zero from which a zero value of epsilon is calculated. The best solution to this problem is to initialise the k field to a small but non-zero value and also select a non-zero value for the walls.

Henry

[gschaider]

Hi!

It seems that this problem is fixed in 1.5 for the Xoodles/pitzDaily (default value set to 1e-8) but NOT in the Xoodles/pitzDaily3D-tutorial-case (values are still 0)

Bernhard

[henry]

Thanks for the report, I will fix it in pitzDaily3D also.

H