CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Programming & Development

implementing k-w to twoPhaseEulerFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 16, 2012, 18:07
Default implementing k-w to twoPhaseEulerFoam
  #1
Ofc
New Member
 
Join Date: Nov 2011
Posts: 9
Rep Power: 5
Ofc is on a distinguished road
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,
Ofc is offline   Reply With Quote

Reply

Thread Tools
Display Modes

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 Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
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


All times are GMT -4. The time now is 03:00.