[quarkz]

Hi,

I need to change the Turbulent kinetic energy and frequency equations in the kOmegaSST formulation. I need to define a slightly different formula for the DkEff(F1) and DomegaEff(F1) and the corresponding laplacian formulation.

I tried to read Kassem's Implement turbulence model on OpenFOAM 1706 guide and being a non-C++ programmer, I got lost half way.

I did create a new model, but I do not know how to change the formulation for the above 2.

Hence, can someone give me a helping hand?

Should I just copy the entire formula to the modified kOmegaSST.C and edit directly from there?

Basically I need to change:

//Turbulent kinetic energy and frequency equations

fvm::laplacian(alpha*rho*DomegaEff(F1), omega_)

to

fvm::laplacian(alpha*sqrt(rho)DomegaEff(F1), sqrt(rho)*omega_)

//Turbulent kinetic energy equation

fvm::laplacian((alpha*rho*DkEff(F1)), k_)

to

fvm::laplacian((alpha*DkEff(F1)), rho*k_)

I have attached the eqns, especially the highlighted laplacian corrections.

How do I represent them in OF?


Thanks!