# how to decide boundary condition of turbulence energy dissipation rate?

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

 October 27, 1999, 10:49 how to decide boundary condition of turbulence energy dissipation rate? #1 Vincent B. Guest   Posts: n/a During computation, I found that the equation of dissipation rate very difficult to converge. Can anybody tell me how to difinite the boundary condition of this eqation to get good convergence?

 October 27, 1999, 11:33 Re: how to decide boundary condition of turbulence energy dissipation rate? #2 John C. Chien Guest   Posts: n/a (1). I am not writing the program right now. So, I don't have the exact definition for you. (2). There are two kinds of definition for the turbulence kinetic energy dissipation function (epsilon). They require different wall boundary conditions for the low Reynolds number model. You have to dig out some papers on low Reynolds number modeling. The review paper by W. Rodi in AIAA Journal is a good place to start. (3). Then, there is also high Reynolds number model with wall function treatment. This is based on the equilibrium boundary layer assumptions. So, there is a simple relationship between the wall skin friction and the turbulence kinetic energy dissipation function (epsilon). There are also several forms available. For this, you also need to check out some papers on the two-equation k-epsilon models. (4). The convergence issue depends on several factors,: the governing equations, the numerical formulations (FD,FV,FE,...), the solution procedures, the boundary conditions, and the computational mesh used. (5). So, the problem is very specific. And there is no general solution to it. It has to be studied based on above factors.

 October 27, 1999, 13:19 Re: how to decide boundary condition of turbulence energy dissipation rate? #3 Stacey Guest   Posts: n/a Vincent, I use the following for k-e model and external flows. I set k by k_inf = 0.5(%turb * U_inf)**2 where %turb is the percent f.s. turbulence. A value of 1% (0.01) is usually good. I then set the dissipation by assuming the turbulent viscosity is larger (say 4x) than the laminar viscosity yielding e = Cmu * rho_inf * k_inf**2 / (4* mu_laminar). A typical value of Cmu is 0.09 and mu_laminar can be had from Sutherland's Law. For internal flows you can use a length scale to get e. Something similar to e = Cmu**0.75 * k_inf**1.5/(0.03*L) where L is a mixing length for the problem. Stacey

 October 27, 1999, 20:50 Re: how to decide boundary condition of turbulence energy dissipation rate? #4 mahesh prakash Guest   Posts: n/a Dear Vincent, You will have to explain your turbulence model and the boundary condition used for epsilon so that we can give some useful advice. For example, if you are solving for the pseudo energy dissipation (see the low reynolds number model of Launder and Spalding), its value is set to zero at the wall. This specification makes the problem stiff and there is bound to be problems with convergence. One can get over such problems by making the grid very fine close to the wall. Using a high Reynolds number model should not give you such problems although in some cases grid refinement blows up the value of epsilon close to the wall (convergent solutions may still be obtained). Cheers

 October 28, 1999, 01:01 Re: how to decide boundary condition of turbulence energy dissipation rate? #5 Vincent B. Guest   Posts: n/a Thank all of you for your valuable advice. I did not introduce the detailed computaion case. I mean, I am calculating flowfield inside an experiment funace, using FD and SIMPLE algorithm, and employing full reynolds stress model as the turbulennce model. In fact, the boundary condition that I am using is the latter one for internal flows. But I got very large residual source for the e equation at nodes near wall, and the total residual source for e equation kept at a fixed high value while other equations converged nicely. SO, what do you think is the most applicable reason?

 October 28, 1999, 04:07 Re: how to decide boundary condition of turbulence energy dissipation rate? #6 mahesh prakash Guest   Posts: n/a Dear Vincent, If your other residuals seem to have converged then most probably your final solution might be correct. Since epsilon values tend to blow up near the wall the epsilon residuals may show up as a very high value in some cases. You can check whether your final solution is correct by carrying out a grid refinement study. Cheers

 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 peob OpenFOAM Running, Solving & CFD 2 August 14, 2014 09:07 therockyy FLOW-3D 0 May 23, 2011 14:19 blackbody Main CFD Forum 1 May 12, 2010 18:54 Sima Phoenics 1 December 1, 2007 19:55 Zhong Lei Main CFD Forum 23 May 14, 1999 13:22

All times are GMT -4. The time now is 05:21.