Add variable coefficient in turbulence model
 

Hi everyone,


I need to add a variable parameter in front of the production term of the omega equation in the k-omega turbulence model.
It should be like


D omega / Dt = beta(x,y) * P(k,omega,U) - D(k,omega;U) + T(k,omega,U)


I'm new to OpenFoam so i thought to add a volScalarField betaNew






kOmega.H as


volScalarField betaNew = 1;






kOmega.C as


betaNew
(
IOobject
(
"betaNew",
this->runTime.timeName(),
this->mesh_,
IOobject::NO_READ,
IOobject::AUTO_WRITE
)
this->mesh_
)





What I'm trying to do is:
- initialize betaNew as unit over the domain
- print its distribution on a file
- modify its distribution somehow
- run the case providing a different distribution



Thank you all

Hi Riccardo!

Did you find a solution about what you trying to do?

I'm working on something similar, but with the k-epsilon model, and I'm still working on it.

I found something on GitHub about adding a source term in a turbulence model, but the source term is a scalar not function of 'x', 'y' and 'z'. It is relatively simple to understand, but because our variable is function of 'x', 'y' (and possibly 'z'), it's less complicated compare of what we're trying to attempt.

In case you're interested, here is the link:

https://github.com/AndreaDesan/LaunderSharmaKEYap



https://www.cfd-online.com/Forums/op...-equation.html

https://www.cfd-online.com/Forums/op...l-2d-mesh.html