CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   dynOneEqEddy for Compressible Flow (http://www.cfd-online.com/Forums/openfoam/81988-dynoneeqeddy-compressible-flow.html)

deji November 12, 2010 09:38

dynOneEqEddy for Compressible Flow
 
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::printStack(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::operator/<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 November 14, 2010 16:21

What is your initial condition for k and other fields?

Fransje November 17, 2010 06:57

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:
Code:

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. :rolleyes:

Kind regards,

Francois.


All times are GMT -4. The time now is 12:41.