[deji]

Good morning and good day FOAMERS. I am trying to utilize the compressible flow version of the LES dynamic K equation model and I am getting some error messages. I was utilizing the Smagorinsky LES model before and in need to test the dynOneEqEddy LES model. Does anyone have a sample case just so I can look at the necessary coefficients needed or can anyone give me any feedback as to how to utilize this model. Thanks. Here is the error message. Any help??



Selecting LES turbulence model dynOneEqEddy
#0 Foam::error:rintStack(Foam::Ostream&) in "/home/ayodeji/OpenFOAM/OpenFOAM-1.6.x/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/ayodeji/OpenFOAM/OpenFOAM-1.6.x/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted:
#3 Foam::dimensioned<double> Foam:perator/<double>(Foam::dimensioned<double> const&, Foam::dimensioned<double> const&) in "/home/ayodeji/OpenFOAM/ayodeji-1.6.x/applications/bin/linuxGccDPOpt/fireFoam"
#4 Foam::compressible::LESModels::dynOneEqEddy::ck_(F oam::GeometricField<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh> const&) const in "/home/ayodeji/OpenFOAM/OpenFOAM-1.6.x/lib/linuxGccDPOpt/libcompressibleLESModels.so"
#5 Foam::compressible::LESModels::dynOneEqEddy::updat eSubGridScaleFields(Foam::GeometricField<Foam::Sym mTensor<double>, Foam::fvPatchField, Foam::volMesh> const&) in "/home/ayodeji/OpenFOAM/OpenFOAM-1.6.x/lib/linuxGccDPOpt/libcompressibleLESModels.so"
#6 Foam::compressible::LESModels::dynOneEqEddy::dynOn eEqEddy(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::basicThermo const&) in "/home/ayodeji/OpenFOAM/OpenFOAM-1.6.x/lib/linuxGccDPOpt/libcompressibleLESModels.so"
#7 Foam::compressible::LESModel::adddictionaryConstru ctorToTable<Foam::compressible::LESModels::dynOneE qEddy>::New(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::basicThermo const&) in "/home/ayodeji/OpenFOAM/OpenFOAM-1.6.x/lib/linuxGccDPOpt/libcompressibleLESModels.so"
#8 Foam::compressible::LESModel::New(Foam::GeometricF ield<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::basicThermo const&) in "/home/ayodeji/OpenFOAM/OpenFOAM-1.6.x/lib/linuxGccDPOpt/libcompressibleLESModels.so"
#9 main in "/home/ayodeji/OpenFOAM/ayodeji-1.6.x/applications/bin/linuxGccDPOpt/fireFoam"
#10 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#11 _start at /build/buildd/glibc-2.9/csu/../sysdeps/i386/elf/start.S:122
Floating point exception

[jason_wy]

What is your initial condition for k and other fields?

[Fransje]

Dear Deji,

This doesn't necessarily look like a wrong specification of the turbulence model, but here is a definition of a ./constant/LESProperties file which should work:

LESModel        dynOneEqEddy;

turbulence      on;

printCoeffs     on;

delta           cubeRootVol;

dynOneEqEddyCoeffs
{
    cubeRootVolCoeffs
    {
        deltaCoeff    1;
    }
    filter      simple;
}

Of course, you can change the delta and the filter to suit your needs.
As it is a "dynamic" LES model, you have no other coefficients to specify since the coefficients used are determined "on the fly" by the turbulent model.

You should be aware of the fact that, to my knowledge, none of the OpenFOAM LES models are truly dynamic in a local sense, not even the locDynOneEqEddy, but are merely domain averages computed "dynamically" at every time step.

Kind regards,

Francois.