
[Sponsors] 
December 15, 2006, 02:46 
Model constants in kepsilon model

#1 
Guest
Posts: n/a

Hi all,
I have a doubt regarding what is the criterion to select/adjust the model constants in standard kepsilon model? Also, for which flow situations/geometries this model won't wor`k? Thanks a lot, Sachin 

December 15, 2006, 11:38 
Re: Model constants in kepsilon model

#2 
Guest
Posts: n/a

Dear Sachin,
Those model constants are not meant to play with.. They have been refined over years based on several flow conditions.. They are not perfect, but they satisfy a lot of consistency checks.. If you change one of them, you must revisit each one of them in order to satisfy a few physical contrains.. You should read a few turbulence modeling papers (usually the older ones) where this consistency checks are better explained. Newer papers assume the reader already knows about. Some books: Mathematical models of turbulence 1972 1972, B. E. Launder and D. B. Spalding. London: Academic Press Inc. (London) Limited. Paper, 169pp., #2.50 Turbulence Models and Their Application in Hydraulics. By W. RODI. International Association for Hydraulic Research, Delft, 1980. Paperback Check early publications by Launder, Rodi, Leschiner, Rotta, Hanjalic, Speziale and if you want to go far in time.. Try the seminal paper by Harlow and Nakayama Harlow, F.H., and P.I. Nakayama, "Turbulence Transport Equations," Phys. Fluids, Vol. 10, pp. 2323Ç2332, 1967. I do not recall if the popular book by Wilcox discuss the origin of the constants.. Turbulence modelling for CFD. By D.C. Wilcox good luck Opaque 

December 17, 2006, 17:56 
Re: Model constants in kepsilon model

#3 
Guest
Posts: n/a

Dear opaque,
Could you briefly describe the consistency you talk about and the physical constrains. I have toyed myself with the coefficients, for instance if I want more turbulence dissipation I would change the Ceps2 to a lower value. In this case what physical constrains would this violate? Best Regards 

December 17, 2006, 23:55 
Re: Model constants in kepsilon model

#4 
Guest
Posts: n/a

Dear Purkmeisster,
It has been quite a few years since I read those papers where the equation constraining Ce2, Ce1, and C_mu was discussed.. Usually in turbulence papers in the late 60's and early 70's.. However, the physical consistency of the model constants must satisfy a few simple turbulent flows. The following is taken from a presentation by Dr. K. Hanjalic (University of Delft) "... The terms have the following physical meanings: P  production of k (of turbulence!) due to meanflow deformation (work of turbulent stresses associated with the mean flow deformation, or transfer of energy from the mean motion to turbulent fluctuations by the action of Reynolds stresses); (note that usually u′iu′j < 0 when ∂Ui/∂xj > 0 so that in most flows P > 0). ε  dissipation of k into heat by viscosity (Note ε > 0 always!) D  diffusive transport of k (viscous, by fluctuating velocity and by fluctuating pressure,respectively) ... a) energy decay in homogeneous turbulence > U1 ∂k/dx1 = âˆ'ε b) rapid distortion: Uj ∂k/dxj = P c) convectiondiffusion equilibrium (pure transport): Uj ∂k/dxj = D d) local turbulence energy equilibrium: 0 = P âˆ' ε e) diffusiondissipation equilibrium: 0 = D âˆ' ε a). In the absence of production (no meanrate of strain, nor body force) in a homogeneous turbulent flows, the turbulence kinetic energy will freely decay and will ultimately be dissipated by viscosity. This is a natural, irreversible process. In the initial period, while the turbulence Reynolds number is still high, the process occurs by energy cascading towards smaller and smaller eddies, independent of viscosity. At a later stage when the Reynolds number becomes sufficiently small, the decay is governed by viscosity. Experiments show (confirmed by DNS) that both, the initial and final periods of decay can be expressed approximately by exponential law, k ∝ tâˆ'n (or ∝ xâˆ'n), with n ≈ 1.2 in the initial period, and n ≈ 2.5 in the final period. This decay law can be used to determine the coefficient associated with the sink term of the scaleproviding equation. b) If a turbulent flow is subjected to a sudden meanrate of strain, the dissipation can be neglected as compared with the imposed turbulence production term. Such is a flow subjected to a sudden contraction, free turbulent flow encountering a solid wall parallel or inclined to the flow direction, or a change in the wall topology, or a boundary layer flow over a stagnant wall encountering transversely moving wall, e.g. turbine rotor. Useful information for modelling and model validation (particularly of secondmoment closures) can be obtained by analyzing the response of turbulence to rapid distortion, since the uncertainties in the model ε equation can be ignored. c) At the free edge of a turbulent flow both production and dissipation can be neglected as compared with the convection and diffusion which are roughly in balance. d) The local energy equilibrium prevails in the inner zone (logarithmic region) of a wall boundary layer at zero or weak pressure gradients. Both the convective and diffusive transport are here small in comparison with the production and dissipation which balance each other locally. This energy balance has a number of additional implications: the shear stress can be regarded as approximately constants, mean velocity obeys the semilogarithmic law and the characteristic turbulence length scale increases linearly with the wall distance. All these features are exploited in tuning the the coefficients in a turbulence model. e) This last case is encountered e.g. in the central zone of a fully developed pipe or channel flow, where the production can be neglected (weak or no mean rate of strain) and convection is zero by definition. Locally dissipated turbulence is supplied by diffusion from other regions in the flow. Although not crucial for most engineering application, satisfying these conditions is often a challenge for a turbulence model. ..." Hope the information is somehow useful .. Regards, Opaque 

December 18, 2006, 07:50 
Re: Model constants in kepsilon model

#5 
Guest
Posts: n/a

Thank you for a very good answer (perhaps good for the wiki?).
May I continue to ask some follow up questions: What I'm thinking of is an applied engineering point of view. That is, in this case, trying to fit the model to a range of experimental data. If a change in a model constant results in a better experimental fit, then what are the pitfalls of the approach, aside from that the constants found are might only be valid over the system investigated (and even only the range investigated)? Best Regards 

December 18, 2006, 09:11 
Re: Model constants in kepsilon model

#6 
Guest
Posts: n/a

What you have to take care about is the interdependance of the constant values one to each other.
Suppose you start with the decay of homogeneus turbulence, then you can fix the value of c_eps2. From the homogeneus shear flow, you can show that c_eps1 depends of c_nu and c_eps2. c_nu is fix from the log law, as for sigma_eps. But sigma_eps also depends on c_nu,c_eps1, c_eps2 and the Karman constant. Finally, sigma_k is given by the value of sigma_eps and c_eps2, so it depends on all the previous test case. So, as an exemple, if you want to change c_eps2, without changing the behaviour of your model in the lag law region, then you also have to change all the other constant to keep the good behavior. I can give you a PDF file  I have written but not finalized and which is in french  which explains all this relationship. ZubenUbi 

December 20, 2006, 03:11 
Re: Model constants in kepsilon model

#7 
Guest
Posts: n/a


December 28, 2006, 10:29 
Re: Model constants in kepsilon model

#8 
Guest
Posts: n/a

Thanks a lot for ur info.!!!


December 28, 2006, 10:45 
Re: Model constants in kepsilon model

#9 
Guest
Posts: n/a

Dear ZubenUbi,
Can I too get the file from you? Also, can I translate it readily into english? My email is: sachintarkunde@yahoo.com 

January 17, 2010, 11:15 
model constants in standard kepsilon model

#10 
New Member
Yogini Patel
Join Date: Jan 2010
Location: Finland
Posts: 15
Rep Power: 9 
Hello ZubenUbi,
I have some confusion regarding to the model constants in standard kepsilon model.So can you give me the pdf file ? My email is yogi.msu@gmail.com Thank you!! Yogini. 

January 18, 2010, 19:32 

#11 
New Member
Join Date: Nov 2009
Posts: 12
Rep Power: 9 
Hi everybody,
would it be possible that anybody could send me the file concerning the model constants, too (sekran@web.de) ? Would be great! Thank you! 

February 18, 2010, 08:39 
model constants in standard kepsilon model

#12 
New Member
Vitor Manuel M. G. da Costa Gomes
Join Date: Feb 2010
Posts: 1
Rep Power: 0 
Dear ZubenUbi
Can I too ask you for the said pdf? thanks in advance! Regards 

March 15, 2011, 15:59 

#13 
New Member
James
Join Date: Mar 2011
Location: GTA Canada
Posts: 1
Rep Power: 0 
Dear ZubenUbi
Sorry to bother you, but can I also have a copy of this PDF. 

June 23, 2011, 19:37 

#14 
Member
RCCS
Join Date: Apr 2011
Posts: 40
Rep Power: 8 
Dear ZubenUbi
Sorry to bother you, but can I also have a copy of this PDF. I have been studying the K Epsilon model constants effect on pipe flow. My email is rsilva@eq.uc.pt Thank you. Cheers. 

March 19, 2012, 08:39 
Pdf copy

#15 
Senior Member
Goutam Saha
Join Date: Dec 2011
Location: UK
Posts: 131
Rep Power: 7 
Can I also have a copy of this PDF (English version).
I have been studying the K Epsilon model constants effect on pipe flow. My Email is g.saha.1@research.gla.ac.uk Thank you. Cheers. 

November 21, 2014, 09:30 
Pdf from Mr ZubenUbi

#16 
New Member
Join Date: Nov 2014
Posts: 1
Rep Power: 0 
Dear all,
Could anyone send me the file of Mr. ZubenUbi?, does not matter if it is in French. Thank you very much in advance (a.gonzalez.sanchez@gmail.com) 

December 4, 2014, 09:13 

#17 
Member
Qureshi M Z I
Join Date: Sep 2013
Posts: 74
Rep Power: 5 
hi
i am working on RANS modeling i have some confusion about constants. could anyone send me the file of Mr. ZubenUbi? english or in french, Regars email. mziqureshi@hotmail.com 

July 18, 2016, 14:35 

#18 
New Member
Ximena Delgadillo
Join Date: Jul 2016
Posts: 2
Rep Power: 0 
Does anyone have the document of Zuben Ubi you could send me ?
A lot of thanks ximena_ax20@hotmail.com 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
about Subgridscale model  impecca  OpenFOAM Running, Solving & CFD  4  December 20, 2013 11:36 
SimpleFoam k and epsilon bounded  nedved  OpenFOAM Running, Solving & CFD  1  November 25, 2008 21:21 
selection of model constants  sachin T  CFDWiki  3  May 22, 2007 04:14 
derivation of epsilon equation of kepsilon model  Ricardo Rezende  Main CFD Forum  5  May 1, 2007 05:38 
changing model constants in ke turbulence model  Sunil  CFX  3  October 3, 2006 12:12 