[PonchO]

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

[abe]

Hi Christoph,

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

ABE

[PonchO]

Originally Posted by abe

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:

//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]

Hi Christoph,

Thank you for sharing.

ABE