buoyantPisoFam error
 

Hello Foamers,

I try to run a LES case with buoyantPisoFoam but whatever I do it always returns me this error :

************************************************** ************************************************** *
Create time

Create mesh for time = 0


Reading g
Reading thermophysical properties

Selecting thermodynamics package hRhoThermo<pureMixture<constTransport<specieThermo <hConstThermo<perfectGas>>>>>
Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model type LESModel
Selecting LES turbulence model oneEqEddy
oneEqEddyCoeffs
{
ce 1.048;
Prt 1;
ck 0.094;
}

Creating field DpDt

#0 Foam::error:: printStack(Foam::Ostream&) in "/home/norbertb/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/norbertb/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted:
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/home/norbertb/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#4 void Foam::divide<Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/home/norbertb/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxGccDPOpt/buoyantPisoFoam"
#5 Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > Foam:: operator/<Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > const&) in "/home/norbertb/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxGccDPOpt/buoyantPisoFoam"
#6 main in "/home/norbertb/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxGccDPOpt/buoyantPisoFoam"
#7 __libc_start_main in "/lib/i686/cmov/libc.so.6"
#8 _start at /usr/src/packages/BUILD/glibc-2.9/csu/../sysdeps/i386/elf/start.S:122
Floating point exception
************************************************** ************************************************** *

Since there is no issue with RANS tutos, I suppose it is caused by LES, but I do not know where exactly.
Any help would be greatly appreciated :)

It seems that error happens at this level, in the createFields.h file :

Info<< "Creating field DpDt\n" << endl;
volScalarField DpDt
(
"DpDt",
fvc:: DpDt(surfaceScalarField("phiU", phi/fvc::interpolate(rho)), p)
);

But I still don't understand how to avoid it ...
Does anyone got the same issue / an idea ?