|January 16, 2012, 18:07||
implementing k-w to twoPhaseEulerFoam
Join Date: Nov 2011
Posts: 9Rep Power: 7
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,
|Thread||Thread Starter||Forum||Replies||Last Post|
|twoPhaseEulerFoam and LES||raagh77||OpenFOAM Running, Solving & CFD||43||November 11, 2013 05:07|
|Something wrong in UEqns.H within twoPhaseEulerFoam||cheng1988sjtu||OpenFOAM||2||June 24, 2011 10:48|
|twoPhaseEulerFoam||freemankofi||OpenFOAM||0||May 23, 2011 16:24|
|stratified horizontal two phase flow usinfg twoPhaseEulerFoam||karthik1414||OpenFOAM||0||April 12, 2011 09:57|
|TwoPhaseEulerFoam bed tutorial case stable in 1.5, crashes in 1.6||hemph||OpenFOAM||3||December 5, 2009 05:19|