An explicit algebraic Reynolds stress and heat flux model for incompressible turbulence: Part II Buoyant flow

This paper presents a derivation of an explicit algebraic model for two-dimensional (2-D) buoyant flows. It is an extension of the work reported in Part I (So et al. [27]). The tensor representation method of Jongen and Gatski [14] is extended to derive an explicit algebraic Reynolds stress model (EASM) for 2-D buoyant flow invoking the Boussinesq approximation. The projection methodology is further extended to treat the heat flux transport equation in the derivation of an explicit algebraic heat flux model (EAHFM) for buoyant flow. Again, the weak equilibrium assumption is invoked for the scaled Reynolds stress and scaled heat flux equation. An explicit algebraic model for buoyant flows is then formed with the EASM and EAHFM. From the derived EAHFM, an expression for the thermal diffusivity tensor in buoyant shear flows is deduced. Furthermore, a turbulent Prandtl number (Pr_T) for each of the three heat flux directions is determined. These directional Pr_T are found to be a function of the gradient Richardson number. Alternatively, a scalar Pr_T can be derived; its value is compared with the directional Pr_T. The EASM and EAHFM are used to calculate 2-D homogeneous buoyant shear flows and the results are compared with direct numerical simulation data and other model predictions, where good agreement is obtained. Dedicated to the memory of the late Professor Charles G. Speziale of Boston University     

Accounting for Buoyancy Effects in the Explicit Algebraic Stress Model: Homogeneous Turbulent Shear Flows

Most explicit algebraic stress models are formulated for turbulent shear flows without accounting for external body force effects, such as the buoyant force. These models yield fairly good predictions of the turbulence field generated by mean shear. As for thermal turbulence generated by the buoyant force, the models fail to give satisfactory results. The reason is that the models do not explicitly account for buoyancy effects, which interact with the mean shear to enhance or suppress turbulent mixing. Since applicable, coupled differential equations have been developed describing these thermal turbulent fields, it is possible to develop corresponding explicit algebraic stress models using tensor representation theory. While the procedure to be followed has been employed previously, unique challenges arise in extending the procedure for developing the algebraic representations to turbulent buoyant flows. In this paper the development of an explicit algebraic stress model (EASM) is confined to the homogeneous buoyant shear flow case to illustrate the methodology needed to develop the proper polynomial representations. The derivation is based on the implicit formulation of the Reynolds stress anisotropy at buoyant equilibrium. A five-term representation is found to be necessary to account properly for the effect of the thermal field. Thus derived, external buoyancy effects are represented in the scalar coefficients of the basis tensors, and structural buoyancy effects are accounted for in additional terms in the stress anisotropy tensor. These terms will not vanish even in the absence of mean shear. The performance of the new EASM, together with a two-equation (2-Eq) model, the non-buoyant EASM of Gatski and Speziale (1993) and a full second-order model, is assessed against direct numerical simulations of homogeneous, buoyant shear flows at two different Richardson numbers representing weak and strong buoyancy effects. The calculations show that this five-term representation yields better results than the 2-Eq model and the EASM of Gatski and Speziale where buoyancy effects are not explicitly accounted for. Received 5 March 2001 and accepted 15 January 2002     

Predicting Buoyant Shear Flows Using Anisotropic Dissipation Rate Models

This paper examines the modeling of two-dimensional homogeneous stratified turbulent shear flows using the Reynolds-stress and Reynolds-heat-flux equations. Several closure models have been investigated; the emphasis is placed on assessing the effect of modeling the dissipation rate tensor in the Reynolds-stress equation. Three different approaches are considered; one is an isotropic approach while the other two are anisotropic approaches. The isotropic approach is based on Kolmogorov's hypothesis and a dissipation rate equation modified to account for vortex stretching. One of the anisotropic approaches is based on an algebraic representation of the dissipation rate tensor, while another relies on solving a modeled transport equation for this tensor. In addition, within the former anisotropic approach, two different algebraic representations are examined; one is a function of the Reynolds-stress anisotropy tensor, and the other is a function of the mean velocity gradients. The performance of these closure models is evaluated against experimental and direct numerical simulation data of pure shear flows, pure buoyant flows and buoyant shear flows. Calculations have been carried out over a range of Richardson numbers (Ri) and two different Prandtl numbers (Pr); thus the effect of Pr on the development of counter-gradient heat flux in a stratified shear flow can be assessed. At low Ri, the isotropic model performs well in the predictions of stratified shear flows; however, its performance deteriorates as Ri increases. At high Ri, the transport equation model for the dissipation rate tensor gives the best result. Furthermore, the results also lend credence to the algebraic dissipation rate model based on the Reynolds stress anisotropy tensor. Finally, it is found that Pr has an effect on the development of counter-gradient heat flux. The calculations show that, under the action of shear, counter-gradient heat flux does not occur even at Ri = 1 in an air flow. This revised version was published online in July 2006 with corrections to the Cover Date.     

Applications of URANS on predicting unsteady turbulent separated flows

Accurate prediction of unsteady separated turbulent flows remains one of the toughest tasks and a practi cal challenge for turbulence modeling. In this paper, a 2D flow past a circular cylinder at Reynolds number 3,900 is numerically investigated by using the technique of unsteady RANS （URANS）. Some typical linear and nonlinear eddy viscosity turbulence models （LEVM and NLEVM） and a quadratic explicit algebraic stress model （EASM） are evaluated. Numerical results have shown that a high-performance cubic NLEVM, such as CLS, are superior to the others in simulating turbulent separated flows with unsteady vortex shedding.     

Evaluation a priori d'un modèle non-linéaire de turbulenceA priori evaluation and improvement of a non-linear model for turbulent flows

The aim of this work is a priori evaluation and improvement of a non-linear model for turbulent flows using the results from direct numerical simulation of Navier–Stokes equations. The algebraic explicit non-linear model recently proposed by Rumsey C.L. et al. [1] is studied. The data base used here comes from a direct numerical simulation of a turbulent flow through a square duct. For this flow, this study shows that the hypothesis of equilibrium state for the anisotropic tensor is correct. The analysis is made using the maps of the second and third invariants of the Reynolds stress tensor. The approach used permits to conclude that the model using a wall function improves the numerical prediction of the anisotropy. To cite this article: O. El Yahyaoui et al., C. R. Mecanique 330 (2002) 27–34     

Nonlinear turbulence models for predicting strong curvature effects

Prediction of the characteristics of turbulent flows with strong streamline curvature, such as flows in turbomachines, curved channel flows, flows around airfoils and buildings, is of great importance in engineering applications and poses a very practical challenge for turbulence modeling. In this paper, we analyze qualitatively the curvature effects on the structure of turbulence and conduct numerical simulations of a turbulent Uduct flow with a number of turbulence models in order to assess their overall performance. The models evaluated in this work are some typical linear eddy viscosity turbulence models, nonlinear eddy viscosity turbulence models （NLEVM） （quadratic and cubic）, a quadratic explicit algebraic stress model （EASM） and a Reynolds stress model （RSM） developed based on the second-moment closure. Our numerical results show that a cubic NLEVM that performs considerably well in other benchmark turbulent flows, such as the Craft, Launder and Suga model and the Huang and Ma model, is able to capture the major features of the highly curved turbulent U-duct flow, including the damping of turbulence near the convex wall, the enhancement of turbulence near the concave wall, and the subsequent turbulent flow separation. The predictions of the cubic models are quite close to that of the RSM, in relatively good agreement with the experimental data, which suggests that these models may be employed to simulate the turbulent curved flows in engineering applications.     

Study on anisotropic buoyant turbulence model

Ｂｕｏｙａｎｔｆｌｏｗｉｓｏｎｅｏｆｔｈｅｆｕｎｄａｍｅｎｔａｌｆｌｏｗｓ．Ｔｈｅｄｉｆｆｅｒｅｎｃｅｏｆｄｅｎｓｉｔｙｂｅｔｗｅｅｎｄｉｓｃｈａｒｇｅｄｆｌｕｉｄａｎｄｔｈｅａｍｂｉｅｎｔｆｌｕｉｄｃａｎｃａｕｓｅｂｕｏｙａｎｔｔｕｒｂｕｌｅｎｔｆ...     

Anisotropic Organised Eddy Simulation for the prediction of non-equilibrium turbulent flows around bodies

The unsteady turbulent flow around bodies at high Reynolds number is predicted by an anisotropic eddy-viscosity model in the context of the Organised Eddy Simulation (OES). A tensorial eddy-viscosity concept is developed to reinforce turbulent stress anisotropy, that is a crucial characteristic of non-equilibrium turbulence in the near-region. The theoretical aspects of the modelling are investigated by means of a phase-averaged PIV in the flow around a circular cylinder at Reynolds number 1.4×10⁵. A pronounced stress–strain misalignment is quantified in the near-wake region of the detached flow, that is well captured by a tensorial eddy-viscosity concept. This is achieved by modelling the turbulence stress anisotropy tensor by its projection onto the principal matrices of the strain-rate tensor. Additional transport equations for the projection coefficients are derived from a second-order moment closure scheme. The modification of the turbulence length scale yielded by OES is used in the Detached Eddy Simulation hybrid approach. The detached turbulent flows around a NACA0012 airfoil (2-D) and a circular cylinder (3-D) are studied at Reynolds numbers 10⁵ and 1.4×10⁵, respectively. The results compared to experimental ones emphasise the predictive capabilities of the OES approach concerning the flow physics capture for turbulent unsteady flows around bodies at high Reynolds numbers.     

Numerical modelling of complex turbulent free‐surface flows with the SPH method: an overview

The gridless smoothed particle hydrodynamics (SPH) method is now commonly used in computational fluid dynamics (CFD) and appears to be promising in predicting complex free‐surface flows. However, increasing flow complexity requires appropriate approaches for taking account of turbulent effects, whereas some authors are still working without any turbulence closure in SPH. A review of recently developed turbulence models adapted to the SPH method is presented herein, from the simplistic point of view of a one‐equation model involving mixing length to more sophisticated (and thus realistic) models like explicit algebraic Reynolds stress models (EARSM) or large eddy simulation (LES). Each proposed model is tested and validated on the basis of schematic cases for which laboratory data, theoretical or numerical solutions are available in the general field of turbulent free‐surface incompressible flows (e.g. open‐channel flow and schematic dam break). They give satisfactory results, even though some progress should be made in the future in terms of free‐surface influence and wall conditions. Recommendations are given to SPH users to apply this method to the modelling of complex free‐surface turbulent flows. Copyright © 2006 John Wiley & Sons, Ltd.     

Lower-Dimensional Manifold (Algebraic) Representation of Reynolds Stress Closure Equations

For complex turbulent flows, Reynolds stress closure modeling (RSCM) is the lowest level at which models can be developed with some fidelity to the governing Navier–Stokes equations. Citing computational burden, researchers have long sought to reduce the seven-equation RSCM to the so-called algebraic Reynolds stress model which involves solving only two evolution equations for turbulent kinetic energy and dissipation. In the past, reduction has been accomplished successfully in the weak-equilibrium limit of turbulence. In non-equilibrium turbulence, attempts at reduction have lacked mathematical rigor and have been based on ad hoc hypotheses resulting in less than adequate models.?In this work we undertake a formal (numerical) examination of the dynamical system of equations that constitute the Reynolds stress closure model to investigate the following questions. (i) When does the RSCM equation system formally permit reduced representation? (ii) What is the dimensionality (number of independent variables) of the permitted reduced system? (iii) How can one derive the reduced system (algebraic Reynolds stress model) from the full RSCM equations? Our analysis reveals that a lower-dimensional representation of the RSCM equations is possible not only in the equilibrium limit, but also in the slow-manifold stage of non-equilibrium turbulence. The degree of reduction depends on the type of mean-flow deformation and state of turbulence. We further develop two novel methods for deriving algebraic Reynolds stress models from RSCM equations in non-equilibrium turbulence. The present work is expected to play an important role in bringing much of the sophistication of the RSCM into the realm of two-equation algebraic Reynolds stress models. Another objective of this work is to place the other algebraic stress modeling efforts in the lower-dimensional modeling context. Received 19 November 1999 and accepted 3 August 2000     

Assessing the Role of Curvature Effects in the Representation of Anisotropy Transport in Algebraic Reynolds-Stress Modelling Applied to Separated Flow

A “curvature correction”, in the sense used herein, is model fragment that is intended to account for some proportion of the error that is introduced into explicit algebraic Reynolds-stress models as a consequence of ignoring the convective transport of stress anisotropy within a framework in which the velocity vector and stress tensor are represented in terms of Cartesian components in highly-curved flow conditions. The present paper examines the ability of such corrections to represent this error, based on a-priori investigations of representative 2-D and 3-D massively separated flows for which the ‘exact’ level of the stress transport is known. In essence, the corrections reflect the assumption that the convection of the Reynolds-stress components, or the associated strain-tensor components, expressed in terms of curvature-oriented coordinates, is negligible. The analysis shows, first, that in general recirculating flows, the contribution of anisotropy transport to the stress balance is generally small, so that any form of related correction is of little consequence. Second, the variants of curvature correction examined correlate poorly with the real anisotropy convection. Thus, while these curvature corrections are useful in very particular conditions, such as flow in highly-curved ducts, they are not generality effective – indeed, possibly counterproductive – and cannot be recommended for inclusion in general numerical schemes.     

Experimental assessment of a nonlinear turbulent stress relation in a complex reciprocating engine flow

A nonlinear turbulent stress relationship, based on an explicit algebraic Reynolds stress closure, is compared against experimental data obtained in a swirl-supported, light-duty engine motored at constant speed. The model relationship is applied to measured mean velocity gradients and turbulence scales, and the predictions compared against the measured shear stress and normal stress anisotropy. Significant improvement over the linear stress relationship typically used in two-equation turbulence models is observed. Conditions under which the model predictions are poor are identified and the reasons for the poor performance discussed.     

Objective tensorial representation of the pressure–strain correlations of turbulence

An explicit algebraic model for the fluctuating pressure–strain rate correlations of turbulence is developed by the use of representation theorems for tensor-valued isotropic tensors, and by invoking the principle of objectivity. The resulting model differs from others by the absence of the vorticity tensor from its formulation. The new model is calibrated by reference to data from homogeneous shear flows, and its potential as a practical tool for the analysis of turbulent flows is demonstrated by numerical simulations of a benchmark two-dimensional shear layer.     

Investigation of transitional and turbulent heat and momentum transport in a rotating cavity

The paper gives the results of the DNS/LES which was performed to investigate the transitional and turbulent non-isothermal flows within a rotor/stator cavity. Computations were performed for the cavity of aspect ratio L = 2–35, Rm = 1.8 and for rotational Reynolds numbers up to 290000. The main purpose of the investigations was to analyze the influence of aspect ratio and Reynolds number on the flow structure and heat transfer. The numerical solution is based on a pseudo-spectral Chebyshev–Fourier–Galerkin collocation approximation. The time scheme is semi-implicit second-order accurate, which combines an implicit treatment of the diffusive terms and an explicit Adams–Bashforth extrapolation for the non-linear convective terms. In the paper we analyze distributions of the Reynolds stress tensor components, the turbulent heat flux tensor components, Nusselt number distributions and the turbulent Prandtl number and other structural parameters, which can be useful for modeling purposes. Selected results are compared with the experimental data obtained for single heated rotating disk by Elkins and Eaton (2000).     

Predicting Noninertial Effects with Linear and Nonlinear Eddy-Viscosity, and Algebraic Stress Models

Three types of turbulence models which account for rotational effects in noninertial frames of reference are evaluated for the case of incompressible, fully developed rotating turbulent channel flow. The different types of models are a Coroiolis-modified eddy-viscosity model, a realizable nonlinear eddy-viscosity model, and an algebraic stress model which accounts for dissipation rate anisotropies. A direct numerical simulation of a rotating channel flow is used for the validation of the turbulence models. This simulation differs from previous studies in that significantly higher rotation numbers are investigated. Flows at these higher rotation numbers are characterized by a relaminarization on the cyclonic or suction side of the channel, and a linear velocity profile on the anticyclonic or pressure side of the channel. The predictive performance of the three types of models are examined in detail, and formulation deficiencies are identified which cause poor predictive performance for some of the models. Criteria are identified which allow for accurate prediction of such flows by algebraic stress models and their corresponding Reynolds stress formulations.     

A dynamic global-coefficient mixed subgrid-scale model for large-eddy simulation of turbulent flows

A dynamic global-coefficient mixed subgrid-scale eddy-viscosity model for large-eddy simulation of turbulent flows in complex geometries is developed. In the present model, the subgrid-scale stress is decomposed into the modified Leonard stress, cross stress, and subgrid-scale Reynolds stress. The modified Leonard stress is explicitly computed assuming a scale similarity, while the cross stress and the subgrid-scale Reynolds stress are modeled using the global-coefficient eddy-viscosity model. The model coefficient is determined by a dynamic procedure based on the global-equilibrium between the subgrid-scale dissipation and the viscous dissipation. The new model relieves some of the difficulties associated with an eddy-viscosity closure, such as the nonalignment of the principal axes of the subgrid-scale stress tensor and the strain rate tensor and the anisotropy of turbulent flow fields, while, like other dynamic global-coefficient models, it does not require averaging or clipping of the model coefficient for numerical stabilization. The combination of the global-coefficient eddy-viscosity model and a scale-similarity model is demonstrated to produce improved predictions in a number of turbulent flow simulations.     

Modeling of Three-Dimensional Turbulent Jet and Boundary-Layer Flows

New anisotropic algebraic constitutive relations for the Reynolds stress tensor are formulated. These relations make it possible to model correctly three-dimensional turbulent flows which cannot be described on the basis of traditional modern semi-empirical models of turbulence. Along with the well-known nonlinear Saffman term, these relations contain new nonlinear terms which take the wall effect into account. Several two- and three-dimensional turbulent flows are calculated nu\-merically using the averaged Navier-Stokes equations. The calculation results are compared with known experimental data.