susp in kEpsilon code in openfoam
 

Hi dears
I'm dealing with k epsilon equations in openfoam to make some desirable change on that. I wonder if somebody could helo to understand this term:


- fvm::SuSp((2.0/3.0)*alpha()*rho()*divU, k_)
in k equation

and also

- fvm::SuSp(((2.0/3.0)*C1_ - C3_)*alpha()*rho()*divU, epsilon_)
in epsilon equation.

so where does 2.0/3.0 comes from? in compaare with the equations n literature.

Regards:)

Hi again !
Problem solved
there is a comperehensive disscutioon about this in here:


https://www.cfd-online.com/Forums/op...arfield-g.html
Regards :)