scispace - formally typeset
Search or ask a question
Journal ArticleDOI

Development of turbulence models for shear flows by a double expansion technique

01 Jul 1992-Physics of Fluids (PHYSICS OF FLUIDS)-Vol. 4, Iss: 7, pp 1510-1520
TL;DR: In this article, a two-equation model and Reynolds stress transport model are developed for turbulent shear flows and tested for homogeneous shear flow and flow over a backward facing step.
Abstract: Turbulence models are developed by supplementing the renormalization group (RNG) approach of Yakhot and Orszag [J. Sci. Comput. 1, 3 (1986)] with scale expansions for the Reynolds stress and production of dissipation terms. The additional expansion parameter (η≡SK/■) is the ratio of the turbulent to mean strain time scale. While low‐order expansions appear to provide an adequate description for the Reynolds stress, no finite truncation of the expansion for the production of dissipation term in powers of η suffices−terms of all orders must be retained. Based on these ideas, a new two‐equation model and Reynolds stress transport model are developed for turbulent shear flows. The models are tested for homogeneous shear flow and flow over a backward facing step. Comparisons between the model predictions and experimental data are excellent.

Summary (1 min read)

Jump to:  and [According]

According

  • It was not an accident that by the early 20th century the ideas of kinetic theory were adapted to describeturbulence with turbulent eddiesas the molecules or building blocks of this "gas.".
  • The assumption (13) was not rigorously justified but it has led to simple equations useful for practical calculations.
  • On the other hand, in the limit of strong strains where r/--, oo: EQUATION /_ _0K which becomes singular.

Did you find this useful? Give us your feedback

Content maybe subject to copyright    Report

NASA Contractor Report 187611
ICASE Report No. 91-65
ICASE
DEVELOPMENT OF TURBULENCE MODELS FOR SHEAR
FLOWS BY A DOUBLE EXPANSION TECHNIQUE
V. Yakhot
S. Thangam
T. B. Gatski
S. A. Orszag
C. G. Speziale
Contract No. NAS1-18605
July 1991
Institute for Computer Applications in Science and Engineering
NASA Langley Research Center
Hampton, Virginia 23665-5225
Operated by the Universities Space Research Association
National Aeronaulic._ and
Space Administration
Langley Re_eorch Center
Hampton, Virginia 23665-5225
'4O
4" W_
O _ _O
I ,--t_
Z Z30


DEVELOPMENT OF TURBULENCE MODELS FOR
SHEAR FLOWS BY A DOUBLE EXPANSION TECHNIQUE
V. Yakhot*, S. Thangam t_, T. B. Gatski§,
S. A. Orszag* and C. G. Speziale t
*Applied &: Computational Mathematics, Princeton University, Princeton, NJ 08544
tICASE, NASA Langley Research Center, Hampton, VA 23665
§NASA Langley Research Center, Hampton, VA 23665
ABSTRACT
Turbulence models are developed by supplementing the renormalization group (RNG)
approach of Yakhot L: Orszag with scale expansions for the Reynolds stress and production
of dissipation terms. The additional expansion parameter (r/ - Sl-(/-g) is the ratio of the
turbulent to mean strain time scale. While low-order expansions appear to provide an
adequate description for the Reynolds stress, no finite truncation of the expansion for the
production of dissipation term in powers of 7/suffices - terms of all orders must be retained.
Based on these ideas, a new two-equation model and Reynolds stress transport model are
developed for turbulent shear flows. The models are tested for homogeneous shear flow
and flow over a backward facing step. Comparisons between the model predictions and
experimental data are excellent.
tPermanent Address: Stevens Institute of Technology, Hoboken, NJ 07030.
tThis research was supported by the National Aeronautics and Space Administration under NASA
Contract No. NAS1-18605 while the authors were in residence at the Institute for Computer Applications
in Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665.


I. Introduction
According to the Kolmogorov theory 1 of turbulence, the dynamics of velocity fluctuations
v(k, w) at the scale l = 2_r/k depend only on the mean rate of dissipation of kinetic energy
g and the length scale I. Neither the integral (L) nor dissipation (ld = 27r/kd) scales enter
the dynamics of the inertial range of scales where ld _<: l <4 L. Based on this theory, the
velocity of eddies of size I scales as gl/31a/3 so that the characteristic (turnover) time of these
eddies is
l ,.., -g_1131213"
T(1) _ v(l----)- (1)
The velocity correlation function can then be represented in the scaling form
<vi(k,w)vi(k,,j)>_k_13/3¢( w )
__113k_t3 - C(k,w) (2)
where the scaling function ¢(x) is to be determined from other considerations. The energy
spectrum E(D) is of the form
E(k) _ k 2 f C(k,_o)aw = CK-g2/ak -5/3 (3)
where C_- is the so called Kolmogorov constant. These results also yield the effective (tur-
bulent) viscosity UT at the scale l = 2rr/k:
I]T(lg ) "" VII "" _1/3]¢ -4/3 . (4)
This effective viscosity plays a profound role in turbulence modelling. For example,
it is quite reasonable to assume that the large-scale properties of the flow are governed
by effective equations of motion similar to the Navier-Stokes equations with the effects of
strong interactions between the velocity fluctuations taken into account through an effective
viscosity (see Yakhot and Orszag 2)
UT _ -_1/3L4/3. (5)
However, the strict notion of eddy viscosity requires the existence of a small parameter
I/L << 1 which is absent in turbulent flow. Still, the eddy viscosity concept proves to be
extremely useful, working much better than expected. This situation is not unique: in a
fluid close to thermal equilibrium, the molecular viscosity representation is very accurate
even when the mean free path A is not that small (A/L __ 1).
Many years before Kolmogorov's work, Osborne Reynolds realized that turbulent flow was
a random system which had to be treated using statistical methods. The work of Reynolds
was similar in concept to the work of Boltzmann, Gibbs and others who formulated the

Citations
More filters
Journal ArticleDOI
TL;DR: The implementation of various types of turbulence modeling in a FOAM computational-fluid-dynamics code is discussed, and calculations performed on a standard test case, that of flow around a square prism, are presented.
Abstract: In this article the principles of the field operation and manipulation (FOAM) C++ class library for continuum mechanics are outlined. Our intention is to make it as easy as possible to develop reliable and efficient computational continuum-mechanics codes: this is achieved by making the top-level syntax of the code as close as possible to conventional mathematical notation for tensors and partial differential equations. Object-orientation techniques enable the creation of data types that closely mimic those of continuum mechanics, and the operator overloading possible in C++ allows normal mathematical symbols to be used for the basic operations. As an example, the implementation of various types of turbulence modeling in a FOAM computational-fluid-dynamics code is discussed, and calculations performed on a standard test case, that of flow around a square prism, are presented. To demonstrate the flexibility of the FOAM library, codes for solving structures and magnetohydrodynamics are also presented with appropriate test case results given. © 1998 American Institute of Physics.

3,987 citations

Journal ArticleDOI
TL;DR: In this paper, the authors present guidelines for using computational fluid dynamics (CFD) techniques for predicting pedestrian wind environment around buildings in the design stage, based on cross-comparison between CFD predictions, wind tunnel test results and field measurements.

1,619 citations

Dissertation
01 Jan 1996
TL;DR: An automatic error-controlled adaptive mesh refinement algorithm is set up in order to automatically produce a solution of pre-determined accuracy, based on a new stabilised and bounded second-order differencing scheme proposed.
Abstract: The accuracy of numerical simulation algorithms is one of main concerns in modern Computational Fluid Dynamics. Development of new and more accurate mathematical models requires an insight into the problem of numerical errors. In order to construct an estimate of the solution error in Finite Volume calculations, it is first necessary to examine its sources. Discretisation errors can be divided into two groups: errors caused by the discretisation of the solution domain and equation discretisation errors. The first group includes insufficient mesh resolution, mesh skewness and non-orthogonality. In the case of the second order Finite Volume method, equation discretisation errors are represented through numerical diffusion. Numerical diffusion coefficients from the discretisation of the convection term and the temporal derivative are derived. In an attempt to reduce numerical diffusion from the convection term, a new stabilised and bounded second-order differencing scheme is proposed. Three new methods of error estimation are presented. The Direct Taylor Series Error estimate is based on the Taylor series truncation error analysis. It is set up to enable single-mesh single-run error estimation. The Moment Error estimate derives the solution error from the cell imbalance in higher moments of the solution. A suitable normalisation is used to estimate the error magnitude. The Residual Error estimate is based on the local inconsistency between face interpolation and volume integration. Extensions of the method to transient flows and the Local Residual Problem error estimate are also given. Finally, an automatic error-controlled adaptive mesh refinement algorithm is set up in order to automatically produce a solution of pre-determined accuracy. It uses mesh refinement and unrefinement to control the local error magnitude. The method is tested on several characteristic flow situations, ranging from incompressible to supersonic flows, for both steady-state and transient problems.

1,418 citations

Journal ArticleDOI
TL;DR: In this article, the RNG κ-e turbulence model derived by Yakhot and Orszag (1986) based on the Renormalization Group theory has been modified and applied to variable-density engine flows in the present study.
Abstract: The RNG κ-e turbulence model derived by Yakhot and Orszag (1986) based on the Renormalization Group theory has been modified and applied to variable-density engine flows in the present study. The original RNG-based turbulence transport approximations were developed formally for an incompressible flow. In order to account for flow compressibility the RNG e-equation is modified and closed through an isotropic rapid distortion analysis. Computations were made of engine compressing/expanding flows and the results were compared with available experimental observations in a production diesel engine geometry. The modified RNG κ-e model was also applied to diesel spray combustion computations. It is shown that the use of the RNG model is warranted for spray combustion modeling since the ratio of the turbulent to mean-strain time scales is appreciable due to spray-generated mean flow gradients, and the model introduces a term to account for these effects. Large scale flow structures are predicted which ar...

1,200 citations

Journal ArticleDOI
TL;DR: In this paper, the main issues and related closures of turbulent combustion modeling are reviewed and a review of the models for non-premixed turbulent flames is given, along with examples of numerical models for mean burning rates for premixed turbulent combustion.

1,069 citations

References
More filters
Journal ArticleDOI
TL;DR: In this article, an extended period numerical integration of a baroclinic primitive equation model has been made for the simulation and the study of the dynamics of the atmosphere's general circulation, and the solution corresponding to external gravitational propagation is filtered by requiring the vertically integrated divergence to vanish identically.
Abstract: An extended period numerical integration of a baroclinic primitive equation model has been made for the simulation and the study of the dynamics of the atmosphere's general circulation. The solution corresponding to external gravitational propagation is filtered by requiring the vertically integrated divergence to vanish identically. The vertical structure permits as dependent variables the horizontal wind at two internal levels and a single temperature, with the static stability entering as a parameter. The incoming radiation is a function of latitude only corresponding to the annual mean, and the outgoing radiation is taken to be a function of the local temperature. With the requirement for thermal equilibrium, the domain mean temperature is specified as a parameter. The role of condensation is taken into account only as it effectively reduces the static stability. All other external sources and sinks of heat are assumed to balance each other locally, and are thus omitted. The kinematics are th...

12,952 citations

Journal ArticleDOI
TL;DR: In this paper, the authors present a review of the applicability and applicability of numerical predictions of turbulent flow, and advocate that computational economy, range of applicability, and physical realism are best served by turbulence models in which the magnitudes of two turbulence quantities, the turbulence kinetic energy k and its dissipation rate ϵ, are calculated from transport equations solved simultaneously with those governing the mean flow behaviour.

11,866 citations

Book
01 Jan 1972
TL;DR: In this paper, the authors present a reference record created on 2005-11-18, modified on 2016-08-08 and used for the analysis of turbulence and transport in the context of energie.
Abstract: Keywords: turbulence ; transport ; contraintes ; transport ; couche : limite ; ecoulement ; tourbillon ; energie Reference Record created on 2005-11-18, modified on 2016-08-08

8,276 citations

Journal ArticleDOI
TL;DR: In this article, the authors developed a model of turbulence in which the Reynolds stresses are determined from the solution of transport equations for these variables and for the turbulence energy dissipation rate E. Particular attention is given to the approximation of the pressure-strain correlations; the forms adopted appear to give reasonably satisfactory partitioning of the stresses both near walls and in free shear flows.
Abstract: The paper develops proposals for a model of turbulence in which the Reynolds stresses are determined from the solution of transport equations for these variables and for the turbulence energy dissipation rate E. Particular attention is given to the approximation of the pressure-strain correlations; the forms adopted appear to give reasonably satisfactory partitioning of the stresses both near walls and in free shear flows. Numerical solutions of the model equations are presented for a selection of strained homogeneous shear flows and for two-dimensional inhomogeneous shear flows including the jet, the wake, the mixing layer and plane channel flow. In addition, it is shown that the closure does predict a very strong influence of secondary strain terms for flow over curved surfaces.

3,855 citations

Journal ArticleDOI
TL;DR: In this article, a dynamic renormalization group (RNG) method for hydrodynamic turbulence was developed, which uses dynamic scaling and invariance together with iterated perturbation methods, allowing us to evaluate transport coefficients and transport equations for the large scale (slow) modes.
Abstract: We develop the dynamic renormalization group (RNG) method for hydrodynamic turbulence. This procedure, which uses dynamic scaling and invariance together with iterated perturbation methods, allows us to evaluate transport coefficients and transport equations for the large-scale (slow) modes. The RNG theory, which does not include any experimentally adjustable parameters, gives the following numerical values for important constants of turbulent flows: Kolmogorov constant for the inertial-range spectrumCK=1.617; turbulent Prandtl number for high-Reynolds-number heat transferPt=0.7179; Batchelor constantBa=1.161; and skewness factor¯S3=0.4878. A differentialK-\(\bar \varepsilon \) model is derived, which, in the high-Reynolds-number regions of the flow, gives the algebraic relationv=0.0837 K2/\(\bar \varepsilon \), decay of isotropic turbulence asK=O(t−1.3307), and the von Karman constantκ=0.372. A differential transport model, based on differential relations betweenK,\(\bar \varepsilon \), andν, is derived that is not divergent whenK→ 0 and\(\bar \varepsilon \) is finite. This latter model is particularly useful near walls.

3,342 citations