CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

dynOneEqEddy for Compressible Flow

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   November 12, 2010, 09:38
Default dynOneEqEddy for Compressible Flow
  #1
Senior Member
 
n/a
Join Date: Sep 2009
Posts: 198
Rep Power: 7
deji is on a distinguished road
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
deji is offline   Reply With Quote

Old   November 14, 2010, 16:21
Default
  #2
New Member
 
Yi Wang
Join Date: Mar 2009
Posts: 17
Rep Power: 8
jason_wy is on a distinguished road
What is your initial condition for k and other fields?
jason_wy is offline   Reply With Quote

Old   November 17, 2010, 06:57
Default
  #3
Senior Member
 
Francois
Join Date: Jun 2010
Posts: 107
Rep Power: 7
Fransje is on a distinguished road
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.

Kind regards,

Francois.
Fransje is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Flow meter Design CD adapco Group Marketing CD-adapco 3 June 21, 2011 08:33
mass flow in is not equal to mass flow out saii CFX 2 September 18, 2009 08:07
What is the difference between liquid reactive flow and gas reactive flow? James Main CFD Forum 6 May 15, 2009 12:14
potential flow vs. Euler flow curious ... Main CFD Forum 23 July 21, 2006 07:40
Plug Flow Franck Main CFD Forum 3 September 4, 2003 05:57


All times are GMT -4. The time now is 22:59.