# implementing k-w to twoPhaseEulerFoam

 Register Blogs Members List Search Today's Posts Mark Forums Read

 January 16, 2012, 18:07 implementing k-w to twoPhaseEulerFoam #1 New Member   Join Date: Nov 2011 Posts: 9 Rep Power: 7 Sponsored Links 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,

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post raagh77 OpenFOAM Running, Solving & CFD 43 November 11, 2013 05:07 cheng1988sjtu OpenFOAM 2 June 24, 2011 10:48 freemankofi OpenFOAM 0 May 23, 2011 16:24 karthik1414 OpenFOAM 0 April 12, 2011 09:57 hemph OpenFOAM 3 December 5, 2009 05:19