[lillberg]

Hi forum and Henry!

I've run into this problem when running coodles in parallel using the Spalart-Allmaras DES model.

[3] #0 Foam::error::printStack(Foam:stream&) in "/local/users/exterli/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpe
nFOAM.so"
[3] #1 Foam::sigFpe::sigFpeHandler(int) in "/local/users/exterli/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.s
o"
[3] #2 ?? in "/lib64/tls/libc.so.6"
[3] #3 ?? in "/lib64/tls/libm.so.6"
[3] #4 ?? in "/lib64/tls/libm.so.6"
[3] #5 pow in "/lib64/tls/libm.so.6"
[3] #6 Foam::pow(Foam::Field<double>&, Foam::UList<double> const&, double const&) in "/local/users/exterli/OpenFOAM/OpenFO
AM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
[3] #7 void Foam::pow<foam::fvpatchfield,>(Foam::GeometricFiel d<double,>&,
Foam::GeometricField<double,> const&, Foam::dimensioned<double> const&) in "/local/users
/exterli/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libcompressibleLESmodels.so "
[3] #8 Foam::tmp<foam::geometricfield<double,> > Foam::pow<Foam::fvPatchField, Foam::vol
Mesh>(Foam::GeometricField<double,> const&, Foam::dimensioned<double> const&) in "/local/
users/exterli/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libcompressibleLESmode ls.so"
[3] #9 Foam::tmp<foam::geometricfield<double,> > Foam::pow<Foam::fvPatchField, Foam::vol
Mesh>(Foam::GeometricField<double,> const&, double const&) in "/local/users/exterli/OpenF
OAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libcompressibleLESmodels.so"
[3] #10 Foam::compressible::LESmodels::SpalartAllmaras::fw (Foam::GeometricField<double,>
const&) const in "/local/users/exterli/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libcompressibl eLESmodels.so"
[3] #11 Foam::compressible::LESmodels::SpalartAllmaras::co rrect(Foam::tmp<foam::geometri cfield<foam::tensor<double>, Foam:
:fvPatchField, Foam::volMesh> > const&) in "/local/users/exterli/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libcompressibl
eLESmodels.so"
[3] #12 Foam::compressible::LESmodel::correct() in "/local/users/exterli/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libco
mpressibleLESmodels.so"
[3] #13 main in "/local/users/exterli/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/c oodles"
[3] #14 __libc_start_main in "/lib64/tls/libc.so.6"
[3] #15 Foam::lduMatrix::lduMatrix(Foam::lduMatrix&, bool) in "/local/users/exterli/OpenFOAM/OpenFOAM-1.4.1/applications/b
in/linux64GccDPOpt/coodles"
srun: error: n107: task3: Floating point exception (core dumped)

Running the same mesh with the same code using any other LES model works just fine. Tried setting floatTransfer to 0 but it didn't change anything. Printed the max/min values of what goes into fw but couldn't see anything unusual. The simulation runs for some 10.000 timesteps and the solution look Ok.

The cluster has an older kernel and glibc and the problem does not occur on my laptop which is running openSuSE 10.3. Hence, this is likely a system specific problem. Anyhow, I would like to hear if someone has any idea how to overcome this. (except upgrading the cluster nodes, 400 cpu's and I'm not admin :-( )

Best Regards

//Eric

[henry]

Hi Eric,

Try with this SpalartAllmaras.C

Henry

[lillberg]

Works like a charm, thanks a lot Henry!

Just a question: Why bounding r by 10?

//Eric

[henry]

r must be clipped otherwise the (r^6)^6 causes a floating point exception. Clipping r with 10 ensures that (r^6)^6 does not cause over-flow in either single or double precision.

Henry

[lillberg]

Thanks again Henry!

//Eric