[kkpal]

hi everyone!
I've came across a very strange error in my LES simulation with dynOneEqEddy model. When my velocity at the initial condition is set to 1m/s, it runs smoothly. However, when velocity is set to be 1.1m/s, or 0.9m/s, or even 0.999 and 1.001, the simulation breaks with the following error:

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 LESModel
Selecting LES turbulence model dynOneEqEddy
Selecting LES delta type cubeRootVol
bounding k, min: 0 max: 0.03 average: 0.03
#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#2 Uninterpreted:
#3 Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libincompressibleRASModels.so"
#5 Foam::incompressible::LESModels::dynOneEqEddy::ck( Foam::GeometricField<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libincompressibleLESModels.so"
#6 Foam::incompressible::LESModels::dynOneEqEddy::upd ateSubGridScaleFields(Foam::GeometricField<Foam::S ymmTensor<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libincompressibleLESModels.so"
#7 Foam::incompressible::LESModels::dynOneEqEddy::dyn OneEqEddy(Foam::GeometricField<Foam::Vector<double >, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&, Foam::word const&) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libincompressibleLESModels.so"
#8 Foam::incompressible::LESModel::adddictionaryConst ructorToTable<Foam::incompressible::LESModels::dyn OneEqEddy>::New(Foam::GeometricField<Foam::Vector< double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libincompressibleLESModels.so"
#9 Foam::incompressible::LESModel::New(Foam::Geometri cField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libincompressibleLESModels.so"
#10 Foam::incompressible::turbulenceModel::addturbulen ceModelConstructorToTable<Foam::incompressible::LE SModel>::NewturbulenceModel(Foam::GeometricField<F oam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libincompressibleLESModels.so"
#11 Foam::incompressible::turbulenceModel::New(Foam::G eometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libincompressibleTurbulenceModel.so"
#12
in "/opt/openfoam222/platforms/linuxGccDPOpt/bin/pimpleFoam"
#13 __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
#14
in "/opt/openfoam222/platforms/linuxGccDPOpt/bin/pimpleFoam"
Floating point exception (core dumped)

I tried on several different computers, without any luck.
From the message it seems that the error comes from the implementation of the dynOneEqEddy. Can someone explain this?

[ashim]

Hi kai Zhang,

Have you figured out the problem? I am facing same problem as you. I need to use 72 m/s. when I am using 1 it's ok. but at 72 it's fail. Could you share your experience.

Regards,

Ali

[kkpal]

sorry to response late. I haven't solved it yet. I changed to use another similar model "homogeneousdynoneeq" and it worked out fine.

[cabrera]

Hey! I have exactly the same problem, anyone got to solve it? I am using OF230.

Thanks!

[ashim]

Hey Cabrera,

I have used " internalField uniform (0 0 0);" in my U file, which solved my problem. I hope it also will work for you.

cheers,

Ali

[cabrera]

Thank you Ali, it worked!