Problem with dynOneEqEddy
 

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::printStack(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?

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

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

dynOneEqEddy
 

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

Thanks!

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

Thank you Ali, it worked!