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     

Turbulence Modeling Effects on the Prediction of Equilibrium States of Buoyant Shear Flows

The effects of turbulence modeling on the prediction of equilibrium states of turbulent buoyant shear flows were investigated. The velocity field models used include a two-equation closure, a Reynolds-stress closure assuming two different pressure-strain models and three different dissipation rate tensor models. As for the thermal field closure models, two different pressure-scrambling models and nine different temperature variance dissipation rate ɛ_τ) equations were considered. The emphasis of this paper is focused on the effects of the ɛ_τ-equation, of the dissipation rate models, of the pressure-strain models and of the pressure-scrambling models on the prediction of the approach to equilibrium turbulence. Equilibrium turbulence is defined by the time rate of change of the scaled Reynolds stress anisotropic tensor and heat flux vector becoming zero. These conditions lead to the equilibrium state parameters, given by /ɛ, _τ/ɛ_τ, , Sk/ɛ and G/ɛ, becoming constant. Here, and _τ are the production of turbulent kinetic energy k and temperature variance , respectively, ɛ and ɛ_τ are their respective dissipation rates, R is the mixed time scale ratio, G is the buoyant production of k and S is the mean shear gradient. Calculations show that the ɛ_τ-equation has a significant effect on the prediction of the approach to equilibrium turbulence. For a particular ɛ_τ-equation, all velocity closure models considered give an equilibrium state if anisotropic dissipation is accounted for in one form or another in the dissipation rate tensor or in the ɛ-equation. It is further found that the models considered for the pressure-strain tensor and the pressure-scrambling vector have little or no effect on the prediction of the approach to equilibrium turbulence. Received 21 April 2000 and accepted 21 February 2001     

An explicit algebraic Reynolds stress and heat flux model for incompressible turbulence: Part I Non-isothermal flow

Tensor representation theory is used to derive an explicit algebraic model that consists of an explicit algebraic stress model (EASM) and an explicit algebraic heat flux model (EAHFM) for two-dimensional (2-D) incompressible non-isothermal turbulent flows. The representation methodology used for the heat flux vector is adapted from that used for the polynomial representation of the Reynolds stress anisotropy tensor. Since the methodology is based on the formation of invariants from either vector or tensor basis sets, it is possible to derive explicit polynomial vector expansions for the heat flux vector. The resulting EAHFM is necessarily coupled with the turbulent velocity field through an EASM for the Reynolds stress anisotropy. An EASM has previously been derived by Jongen and Gatski [10]. Therefore, it is used in conjunction with the derived EAHFM to form the explicit algebraic model for incompressible 2-D flows. This explicit algebraic model is analyzed and compared with previous formulations including its ability to approximate the commonly accepted value for the turbulent Prandtl number. The effect of pressure-scrambling vector model calibration on predictive performance is also assessed. Finally, the explicit algebraic model is validated against a 2-D homogeneous shear flow with a variety of thermal gradients. Dedicated to the memory of the late Professor Charles G. Speziale of Boston University     

Correcting Anisotropic Gradient Transport of k

The gradient transport model for k is extended to classes of turbulent flows for which the gradient transport hypothesis is relevant but the anisotropy of the Reynolds stress, to which the eddy diffusivity is proportional, is large and variable. In highly anisotropic turbulence the standard isotropic model used in engineering practice is fundamentally wrong and the conventional anisotropic approximation inadequate. The work is motivated by the important observations that the eddy diffusivity coefficient for a standard gradient transport model for various transported quantities is a factor of 3–10 times larger in highly anisotropic turbulence than that used in standard engineering models. While the conventional anisotropic eddy diffusivity approximation appears adequate for material conserved scalars it is inadequate for k. The problem is solved by addressing the anisotropy of the turbulent transport of k at the level of the underlying third order tensor. It is shown that, unlike the traditional transport models for k, the orientation of the anisotropy with respect to the direction of the gradient plays a crucial role not accounted for in conventional models used in engineering calculations. The new anisotropic eddy diffusivity tensor is quadratic in the anisotropy (the traditional model is linear in the anisotropy). It is shown that the new more rigorous anisotropic eddy diffusivity varies 300% more than the standard model comparing the isotropic limit to the value for the two-dimensional limit. The two-dimensional limit is important in strongly stably stratified flows, in pressure gradient or shock driven flows and in rotating flows. Using the simple shear and the homogeneous non-equilibrium Rayleigh Taylor turbulence the new anisotropic diffusivity tensor is validated in inhomogeneous Rayleigh Taylor turbulence at early and late times.     

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     

The evolution of Hooke's law due to texture development in FCC polycrystals

Initially isotropic aggregates of crystalline grains show a texture-induced anisotropy of both their inelastic and elastic behavior when submitted to large inelastic deformations. The latter, however, is normally neglected, although experiments as well as numerical simulations clearly show a strong alteration of the elastic properties for certain materials. The main purpose of the work is to formulate a phenomenological model for the evolution of the elastic properties of cubic crystal aggregates. The effective elastic properties are determined by orientation averages of the local elasticity tensors. Arithmetic, geometric, and harmonic averages are compared. It can be shown that for cubic crystal aggregates all of these averages depend on the same irreducible fourth-order tensor, which represents the purely anisotropic portion of the effective elasticity tensor. Coupled equations for the flow rule and the evolution of the anisotropic part of the elasticity tensor are formulated. The flow rule is based on an anisotropic norm of the stress deviator defined by means of the elastic anisotropy. In the evolution equation for the anisotropic part of the elasticity tensor the direction of the rate of change depends only on the inelastic rate of deformation. The evolution equation is derived according to the theory of isotropic tensor functions. The transition from an elastically isotropic initial state to a (path-dependent) final anisotropic state is discussed for polycrystalline copper. The predictions of the model are compared with micro–macro simulations based on the Taylor–Lin model and experimental data.     

快速颗粒流碰撞应力的动理模型

颗粒碰撞是快速颗粒流的动量传递,能量传递及耗散的主要机制。因此,碰撞应力是快速颗粒流研究的主要内容。本文基于大小均匀、光滑圆球颗粒的非弹性碰撞模型,类比气体分子动理论中的处理方法得出快速颗粒流积分形式的碰撞应力,能量传递通量及能量耗散率等关系。在颗粒速度分布函数的一级近似式符合Maxwell分布的假设下,对碰撞应力积分得出相应于Boltzmann方程二级近似解的碰撞应力表达式。在简单快速颗粒流动条件下,本文的理论结果与实验资料符合较好。     

Irreversible thermodynamics of liquid crystal interfaces

A macroscopic theory for the dynamics of isothermal compressible interfaces between nematic liquid crystalline polymers and isotropic viscous fluids has been formulated using classical irreversible thermodynamics. The theory is based on the derivation of the interfacial rate of entropy production for ordered interfaces, that takes into account interfacial anisotropic viscous dissipation as well as interfacial anisotropic elastic storage. The symmetry breaking of the interface provides a natural decomposition of the forces and fluxes appearing in the entropy production, and singles out the symmetry properties and tensorial dimensionality of the forces and fluxes. Constitutive equations for the surface extra stress tensor and for surface molecular field are derived, and their use in interfacial balance equations for ordered interfaces is identified. It is found that the surface extra stress tensor is asymmetric, since the anisotropic viscoelasticity of the nematic phase is imprinted onto the surface. Consistency of the proposed surface extra stress tensor with the classical Boussinesq constitutive equation appropriate to Newtonian interfaces is demonstrated. The anisotropic viscoelastic nature of the interface between nematic polymers (NPs) and isotropic viscous fluids is demonstrated by deriving and characterizing the dynamic interfacial tension. The theory provides for the necessary theoretical tools needed to describe the interfacial dynamics of NP interfaces, such as capillary instabilities, Marangoni flows, wetting and spreading phenomena.     

Positivité de la dissipation intrinsèque d'une classe de modèles d'endommagement anisotropes non standardsPositivity of intrinsic dissipation of a class of nonstandard anisotropic damage models

A nonstandard thermodynamics framework ensures the positivity of the dissipation due to degradation mechanisms for damage states represented by a symmetric second order tensor. The proof of the positivity of the intrinsic dissipation is given. An increasing damage, in terms of positive principal values of the damage rate tensor, guaranties this positivity for the considered class of models, extending then to induced anisotropy the isotropic case property of a positive damage rate. This result gives many possibilities of modeling for anisotropic damage. To cite this article: R. Desmorat, C. R. Mecanique 334 (2006).     

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.     

Explicit algebraic stress modelling of homogeneous and inhomogeneous flows

Capability of the explicit algebraic stress models to predict homogeneous and inhomogeneous shear flows are examined. The importance of the explicit solution of the production to dissipation ratio is first highlighted by examining the algebraic stress models performance at purely irrotational strain conditions. Turbulent recirculating flows within sudden expanding pipes are further simulated with explicit algebraic stress model and anisotropic eddy viscosity model. Both models predict better stress–strain interactions, showing reasonable shear layer developments. The anisotropic stress field are also accurately predicted by the models, though the anisotropic eddy viscosity model of Craft et al. returns marginally better results. Copyright © 2005 John Wiley & Sons, Ltd.     

A numerical method for steady and nonisothermal viscoelastic fluid flow for high Deborah and Péclet numbers

A combined finite element/streamline integration method is presented for nonisothermal flows of viscoelastic fluids. The attention is focused on some characteristic problems that arise for numerical simulation of flows with high Deborah and Péclet numbers. The two most important problems to handle are the choice of an outflow boundary condition for not completely developed flow and the treatment of the dissipative term in the temperature equation. The ability of the numerical method to handle high Deborah and Péclet numbers will be demonstrated on a contraction flow of an LDPE melt with isotropic and anisotropic heat conductivity. The influence of anisotropic heat conduction and the difference between the stress work and mechanical dissipation will be discussed for contraction flows. Received:  4 February 1997 Accepted: 23 October 1997     

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.     

Modeling anisotropic Reynolds-stress dissipation in particle- or droplet-laden flows

Particles or droplets dispersed in turbulent flows at sufficiently high volume loadings lead to modifications of turbulence characteristics. More specifically, in a detailed experimental investigation by Poelma and coworkers, where the particle phase moves with a non-zero mean velocity relative to the fluid phase, it was found that anisotropic Reynolds-stress dissipation is induced. Recently, we have proposed a model that can account for this effect in RANS-based and PDF method simulations. In our previous work, however, no simulation results of the RANS model formulation were presented. In the present work, a new compact tensorial RANS formulation is presented and the new formulation is validated against the experimental data of Poelma and coworkers.     

Scale properties of turbulent transport and coherent structure in stably stratified flows

The empirical mode decomposition (EMD) is used to study the scale properties of turbulent transport and coherent structures based on velocity and temperature time series in stably stratified turbulence. The analysis is focused on the scale properties of intermittency and coherent structures in different modes and the contributions of energy-contained coherent structures to turbulent scalar counter-gradient transport (CGT). It is inferred that the velocity intermittency is scattered to more modes with the development of the stratified flow, and the intermittency is enhanced by the vertical stratification, especially in small scales. The anisotropy of the field is presented due to different time scales of coherent structures of streamwise and vertical velocities. There is global counter-gradient heat transport close to the turbulence-generated grid, and there is local counter-gradient heat transport at certain modes in different positions. Coherent structures play a principal role in the turbulent vertical transport of temperature.     

A possible explanation of countergradient fluxes in homogeneous turbulence

This study addresses the phenomenon of persistent countergradient (PCG) fluxes of momentum and heat (density) as observed in homogeneous turbulence forced by shear and stratification. Countergradient fluxes may occur at large scales when stratification is strong. However, they always occur at small scales, independently of stratification. A conceptional model is introduced to explain PCG fluxes at small scales as the result of the collision of large-scale fluid parcels. The large parcels collide under the driving force of inclined vortex structures (in a shear-dominated flow) or of buoyancy (in a strongly stratified shear flow). This collision model also explains the PCG heat flux in an unsheared stratified flow with zero average momentum flux. It is found that the energy of the small-scale PCG motions is provided (i) by quick transport of kinetic energy from the scales of production to relatively slowly dissipating scales if the flow is shear-driven and (ii) by conversion of available potential energy to kinetic energy at small scales when the flow is stratified. The collision mechanism is an inherent property of the turbulence dynamics. Therefore, the PCG fluxes at small scales reflect a universal character of homogeneous turbulence, and are found over a large range of Reynolds numbers. The Prandtl (or Schmidt) number influences the rate of dissipation of temperature (or density) variance but not the dissipation rate of the velocity variance. In stratified flows, therefore, the number directly affects the strength of the PCG heat flux at small scales. It is found, however, that the PCG momentum flux is also altered slightly when the Prandtl number is large enough to sustain small buoyantly moving parcels after collision.     

Some developments in computational modeling of turbulent flows

In this paper, some recent developments of two turbulence closure schemes at ICOMP, NASA Lewis will be discussed. One is the Reynolds-stress algebraic equation model and the other is the Reynolds-stress transport equation model. Various model constraints required by the rapid distortion theory, the invariant theory and the realizability principle, etc. will be described in the model development. The models discussed are for high-turbulent Reynolds number flows, so that the near-wall turbulence and the low-Reynolds-number turbulence are not discussed here.     

Homogenization problems of a hollow cylinder made of elastic materials with discontinuous properties

In this paper, we present the homogenization of an anisotropic hollow layered tube with discontinuous elastic coefficients. We focus on some aspects of technological importance, such as the effective coefficients of anisotropic materials, the behavior of the homogenized displacements and stresses, the discontinuities of in-plane shear, hoop and longitudinal stresses, the homogenization-induced anisotropy in the isotropic case. We conclude that the problem of cylindrically anisotropic tubes under extension, torsion, shearing and pressuring is stable by homogenization and we define the effective tensor of the material elastic coefficients. Some numerical examples confirm the theoretical results.     

Non-symmetric viscoelasticity of anisotropic polymer liquids

The liquid crystalline (LC) polymers are considered as anisotropic viscoelastic liquids with nonsymmetric stresses. A simple constitutive equation for nematic polymers describing the coupled relaxation of symmetric and antisymmetric parts of the stress tensor is formulated. For illustration of non-symmetric anisotropic viscoelasticity, the simplest viscometric flows of polymeric nematics in the magnetic field are considered. The frequency and shear rate dependencies of extended set of Miesowicz viscosities are predicted. Received: 23 March 1999/Accepted: 13 December 1999