[quarkz]

Hi,

I need to modify the turbulent kinetic energy eqn. The new k formulation involves the speed of sound (c). I can compute c in the eqn using the formula c = sqrt(gamma*R*T), where gamma, R and T are specific heat ratios, gas constant and temperature. I also realise that there is a computation of c in rhoCentralFoam.C:

volScalarField c("c", sqrt(thermo.Cp()/thermo.Cv()*rPsi));

So my question is how do I call the value of c from rhoCentralFoam?

Or if using the 1st mtd, how do I call gamma, R and T?

What are the files which I must modify?

Being a newbie, I tried to add:

#include "psiThermo.H"

to my modified kOmegaSST.C (obtained using Prof Nilsson's tutorial)

and also add directories to my Make/options:

-I$(LIB_SRC)/thermophysicalModels/basic/lnInclude \
    -I$(LIB_SRC)/thermophysicalModels/radiation/lnInclude \
    -I$(LIB_SRC)/thermophysicalModels/specie/lnInclude \
    -I$(LIB_SRC)/thermophysicalModels/solidThermo/lnInclude \
    -I$(LIB_SRC)/thermophysicalModels/solidSpecie/lnInclude \

Although it compiled successfully, once I add things like thermo.Cp into the code, it complains:

error: ‘thermo’ was not declared in this scope

I also tried to add the line:

volScalarField c("c", sqrt(thermo.Cp()/thermo.Cv()*rPsi));

directly into my code but got the same error.

Any idea? Of cos, if possible, the 1st mtd seems much easier but I also have no idea how it can be implemented.

[HPE]

what did you mean by calling c? you compute c inside rhoFoam and you it thereat?

[quarkz]

Sorry I mean I need to obtain the value of c for calculation of turbulent kinetic energy eqn. Since I'll be running rhoCentralFoam with kOmegaSST, is it possible to get the current value of c while running rhoCentralFoam and feed it into the turbulent kinetic energy eqn?

Else, I would need to calculate c using the formula in the source file of the turbulent kinetic energy eqn:

thermo.Cp()/thermo.Cv()*rPsi