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

defining constant transport property inside a turbulence model

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

Like Tree1Likes
  • 1 Post By PonchO

Reply
 
LinkBack Thread Tools Display Modes
Old   May 2, 2013, 05:59
Default defining constant transport property inside a turbulence model
  #1
New Member
 
Christoph Kratzsch
Join Date: Nov 2011
Location: Freiberg
Posts: 20
Rep Power: 5
PonchO is on a distinguished road
Hello Foamers,

i want to implement a turbulence model based on a k-epsilon-model for turbulent magnetohydrodynamics-flow.
Therfore i have to add additional dissipation terms in the epsilon-equation and in the k-equation.
These terms containing the density "rho" and the electrical conductivity "sigma".

Because that properties are constants in my incompressible case, i can't introduce them like it's done in the compressible turbulence-model.

So how can i introduce or address them inside my turbulence-model-files?

Best Regards,

Christoph
PonchO is online now   Reply With Quote

Old   June 19, 2013, 05:33
Default
  #2
abe
Member
 
ABE
Join Date: Jul 2012
Posts: 37
Rep Power: 4
abe is on a distinguished road
Hi Christoph,

Have you find a way to call rho inside turbulence model for incompressible case? if yes, please share.

ABE
abe is offline   Reply With Quote

Old   June 21, 2013, 07:48
Default
  #3
New Member
 
Christoph Kratzsch
Join Date: Nov 2011
Location: Freiberg
Posts: 20
Rep Power: 5
PonchO is on a distinguished road
Quote:
Originally Posted by abe View Post
Hi Christoph,

Have you find a way to call rho inside turbulence model for incompressible case? if yes, please share.

ABE
Hi abe,

after some trial and error i've found a solution.

In my .C-File i read the properties as follows:

Code:
//that part you should find in your turbulence-model-code
void kEpsilonZeta::correct()
{
    RASModel::correct();

    if (!turbulence_)
    {
        return;
    }

//Here the reading of the properties starts
    Info<< "\nReading transportProperties\n" << endl;
    IOdictionary transportProperties
    (
        IOobject
        (
            "transportProperties",
            runTime_.constant(),
            mesh_,
            IOobject::MUST_READ,
            IOobject::NO_WRITE,
            false
        )
    );


    dimensionedScalar nu
    (
        transportProperties.lookup("nu")
    );

    
    dimensionedScalar rho
    (
        transportProperties.lookup("rho")
    );

    // electromagnetic properties
    dimensionedScalar mu
    (
        transportProperties.lookup("mu")
    );

    dimensionedScalar sigma
    (
        transportProperties.lookup("sigma")
    );

    dimensionedScalar Cm
    (
        transportProperties.lookup("Cm")
    );

//Now the Solution of the turbulence transport equations follows
Hope that helps you too.

Best regards,

Christoph
abe likes this.
PonchO is online now   Reply With Quote

Old   June 21, 2013, 10:00
Default
  #4
abe
Member
 
ABE
Join Date: Jul 2012
Posts: 37
Rep Power: 4
abe is on a distinguished road
Hi Christoph,

Thank you for sharing.

ABE
abe 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
Modified k-e turbulence model UDF - help! Travis FLUENT 5 January 20, 2014 08:34
Water subcooled boiling Attesz CFX 7 January 5, 2013 04:32
K - epsilon VS SST turbulence model Maicol Main CFD Forum 0 November 30, 2012 17:25
build your own turbulence model with buoyancy Thomas Baumann OpenFOAM 11 November 23, 2009 09:53
Two-Phase Buoyant Flow Issue Miguel Baritto CFX 4 August 31, 2006 13:02


All times are GMT -4. The time now is 13:03.