
[Sponsors] 
October 26, 2011, 09:49 
Kε Model and SIMPLE?

#1 
New Member
anonymous
Join Date: Aug 2011
Posts: 5
Rep Power: 7 
Hi Guys,
I wanted to get cleared up about how exactly the kε fits into an incompressible solver that doesn't have turbulence handling. I am reading the Versteeg book Introduction to CFD the Finite Volume Method. In the book he describes pretty thoroughly how to implement the SIMPLE method for a laminar solver, but only has a short chapter giving the transport equations for K and ε but he doesn't tell at all how to implement this turbulence modelling. My question is how is kε implemented? Does it remove all swirling vortexes and eddies from the output? Searching online doesn't yield anything about implementation, only a reiteration of what the transport equations mean... Thanks! 

May 25, 2012, 15:33 

#2 
Member
Michael Moor
Join Date: May 2012
Location: Ireland
Posts: 30
Rep Power: 6 
I hear ya! I am currently trying to implement my own 2D cfd code in MATLAB, i'd like to think that I am 80% of the way there, I'm looking for help with it, but not sure where to go!! If anyone is interested my code is very well annotated, I might not be the strongest programmer, but it is fairly explicit!
Did you find out how to do this in the end? I would GREATLY appreciate any help/advice! Best regards, Michael 

May 25, 2012, 15:53 

#3 
New Member
anonymous
Join Date: Aug 2011
Posts: 5
Rep Power: 7 
Hey Michael !
Yes I figured out the situation in the end. The RANS model uses modified versions of the NS equations, and in terms of the SIMPLE algorithm you solve the discretized RANS instead of discretized regular NS. With the RANS you have more variables, and these are simply just scalar variables, so you solve them after you solve momentum, pressure correction ,etc... Note though that in RANS all eddies and these shortscale transient phenomena will be gone. May I ask if you are a student in CFD? A master's student perhaps? Best, Miz 

May 25, 2012, 17:38 

#4 
Member

Why did you choose kepsilon and not other models such as komega or others. What are pros and cons of this model and is it application specific?


May 26, 2012, 07:07 

#5  
Member
Michael Moor
Join Date: May 2012
Location: Ireland
Posts: 30
Rep Power: 6 
Quote:
Those "scalar" variables are in fact the Reynolds stresses, and all zero, one and two equations models use the Boussinesq eddy hypothesis (BEH) to relate the reynolds stresses (basically stresses due to turbulence) to mean flow rates ( since we have 4 equations, continuity and 3 RANS, but now have MEAN unknowns u, v,w, p and 6 independent reynolds stresses (10 unknowns)... therefore, we either need to relate these unknowns back to the mean variables ( which is done through the BEH and a turbulent eddy viscosity mu_t), or solve each REynolds stress with it's own equation to balance the system, which would be the RANS turbulence model, which is more accurate, but computationally expensive. Therefore, all the zero, one and two equation models attempt to do is to describe the turbulent viscosity mu_t... that's about it, where a zero equation e.g. Pandtl mixing length, uses no additional equations, a one equation model, e.g. Spalart Almaras solve one additional equation to describe mu_t etc... nothing too complicated. See Versteeg and Malalasakera, it's pretty good. Be wary though, as the location of your first cell in the domain from a wall is very important, as turbulence is described differently in different portions of a boundary layer, again, see Versteeg, it's pretty clear cut (after a while!) yes I am a masters student, in Aeronautical engineering, it is a taught masters with a thesis as well, and mine is on this development of a code... tough going personally as I had to figure matlab out first! I hope that that clears anything up (if there was anything to be cleared!) Whilst I haven't incorporated turbulence yet, you are right, PATANKAR does say to calculate the turbulence after correcting the velocities, hopefully it will be as easy as that! BR, Michael 

May 26, 2012, 07:23 

#6  
Member
Michael Moor
Join Date: May 2012
Location: Ireland
Posts: 30
Rep Power: 6 
Quote:
I have chosen the Prandtl mixing length, as it is easy to implement, comutationally light, and gives reasonable predictions for velocity gradients in boundary layers. Then I am going to use the kepsilon as it is described (a little) in Versteeg... as well as being generally fairly well rounded ( not as good as komega though). The difference between the two is that there is a term in the epsilon transport equation that cannot be resolved at the wall, and so wall functinos are used to extrapolate the data there, which is a drawback... whereas in the komega model, it solve all the way to the wall, but is very sensitive to the assigned freestream values of omega... this problem however is fixed by the Shear Stress Transport (SST) komega model, which is essentially a kepsilon model in freestream and for the outer region of the boundary layer, and then blends into a komega model at the wall, baically the best of both worlds! Since I am writing the code, I don't want to make life any more difficult than it already is!! As for applicability, the kepsilon is robust, economical, has a fairly wide range of flow, but is not good in: swirling flows, large pressure gradients, separation, jet flows, and with all zero, one, and two equation models, which use the Boussinesq eddy hypothesis, the Reynolds stresses are modelled as being isotropic ( same in all directions), which is not the case in reality. the komega is good external aerodynamic flows, can handle adverse pressure gradients, mixin layers, jet and wake flows, no need for damping functions at the wall, and can be used with or without wall functions! if you need any more help i'd be happy to help. BR, Michael 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
simple model, difficult outlet  Eric  CFX  7  May 23, 2014 08:13 
SIMPLE algorithm and keps turbulence model  zqnwpu  Main CFD Forum  11  May 28, 2012 04:08 
Floworks: simple model of combustion  Cornell  FloEFD, FloWorks & FloTHERM  0  August 18, 2008 19:53 
A simple plastic solid model in FSI  jinwon park  Main CFD Forum  7  April 30, 2008 05:05 
Simple methane kinetics model  Andy  Main CFD Forum  0  December 4, 2006 04:27 