add a source term in incompressible k-epsilon model
I am going to add a new source term (only non-zero at a few cells and not varying with time) in the kinetic energy equation in the kEpsilon code. I have followed a tutorial from the forum to make a new copy of the kEpsilon code and called it mykEpsilon. I also changed simpleFoam to read a new scalar field, mySource. My question is how to pass this scalar field into mykEpsilon so that I can write:
+ fvm::div(phi_, k_)
- fvm::Sp(fvc::div(phi_), k_)
- fvm::laplacian(DkEff(), k_)
- fvm::Sp(epsilon_/k_, k_)
I read somewhere that this can be done by passing the field into the class through a pointer, but not exactly know now. Any hints?
|All times are GMT -4. The time now is 23:31.|