implementing k-w to twoPhaseEulerFoam
 

Hello Foamers,

right now I try to implement the k-w-model into the twoPhaseEulerFoam solver by replacing the kEpsilon.H by my own one. Some questions came up, I cannot solve by myself:

1. In the solver the continuous phase turbulent viscosity is given by
nutb = sqrt(k)/epsilon instead of sqr(k)/epsilon. Even though the units should not fit, they do. If I replace them by nutb = k/omega, the solver answers after compiling & running a testcase:

Calculating field nuEffa

LHS and RHS of + have different dimensions
dimensions : [0 0 1 0 0 0 0] + [0 2 -1 0 0 0 0]

(nuEffa = sqrt(Ct)*nutb + nub)

What did I do wrong?

2. In the original k-w model for one phase I can find the term
fvm::Sp(fvc::div(phi_),omega_)
This part I cannot identifiy in the equation for the k-w model. To me it seems this is an additional source? Where does it come from? Do I have to include it to my code?

3. In the solver there is a calculation for the tensor G, which is done
if (yPlus > 11.6)
Can I replace this calculation with the calculation of the one-phase solver and keep the if? Or do I have to delete it?


Thanks for your time & help,