implementing k-w to twoPhaseEulerFoam
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
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,
|All times are GMT -4. The time now is 11:16.|