A unified theory of turbulent flow

Summary A critical analysis is made of the assumptions underlying Reynolds' equations for turbulent flow. It is shown that there are regions in a flow field where these assumptions break down, and it is therefore necessary to separate the flow into a turbulent core and a laminar sublayer. The importance of the boundary conditions to be imposed on the mean velocities and Reynolds stresses at the junction is emphasized as this is the way in which the effect of surface roughness enters the theory. A set of equations for calculating turbulent flows is proposed. The distinctive feature is that the turbulent stresses are represented as the difference between viscous terms with a large eddy viscosity and terms satisfying auxiliary differential equations proposed by Broszko. These terms may be associated with the free and wall turbulence respectively. The theory enables the idea of a large eddy viscosity to be applied even where the velocity gradient is large. The results obtained for specific configurations, which will be reported in detail in future papers, are previewed.     

Compound wall treatment with low‐Re turbulence model

A generalized treatment for the wall boundary conditions relating to turbulent flows is developed that blends the integration to a solid wall with wall functions. The blending function ensures a smooth transition between the viscous and turbulent regions. An improved low Reynolds number k?ε model is coupled with the proposed compound wall treatment to determine the turbulence field. The eddy viscosity formulation maintains the positivity of normal Reynolds stresses and Schwarz' inequality for turbulent shear stresses. The model coefficients/functions preserve the anisotropic characteristics of turbulence. Computations with fine and coarse meshes of a few flow cases yield appreciably good agreement with the direct numerical simulation and experimental data. The method is recommended for computing the complex flows where computational grids cannot satisfy a priori the prerequisites of viscous/turbulence regions. Copyright © 2011 John Wiley & Sons, Ltd.     

Direct numerical simulation of a freely decaying turbulent interfacial flow

Whereas Large Eddy Simulation (LES) of single-phase flows is already widely used in the CFD world, even for industrial applications, LES of two-phase interfacial flows, i.e. two-phase flows where an interface separates liquid and gas phases, still remains a challenging task. The main issue is the development of subgrid scale models well suited for two-phase interfacial flows. The aim of this work is to generate a detailed data base from direct numerical simulation (DNS) of two-phase interfacial flows in order to clearly understand interactions between small turbulent scales and the interface separating the two phases. This work is a first contribution in the study of the interface/turbulence interaction in the configuration where the interface is widely deformed and where both phases are resolved by DNS. To do this, the interaction between an initially plane interface and a freely decaying homogeneous isotropic turbulence (HIT) is studied. The densities and viscosities are the same for both phases in order to focus on the effect of the surface tension coefficient. Comparisons with existing theories built on wall-bounded or free-surface turbulence are carried out. To understand energy transfers between the interfacial energy and the turbulent one, PDFs of the droplet sizes distribution are calculated. An energy budget is carried out and turbulent statistics are performed including the distance to the interface as a parameter. A spectral analysis is achieved to highlight the energy transfer between turbulent scales of different sizes. The originality of this work is the study of the interface/turbulence interactions in the case of a widely deformed interface evolving in a turbulent flow.     

A qualitative assessment of the role of a viscosity depending on the third invariant of the rate-of-deformation tensor upon turbulent non-Newtonian flow

The numerical simulation of some non-Newtonian effects in wall and wall-free turbulent flows, such as drag reduction in pipe flows or the decrease in transverse normal Reynolds stresses, has been attempted in the past with a limited degree of success on the basis of modified wall functions applied to traditional turbulence models (k–ε), rather than through more realistic rheological constitutive equations. In this work, it is qualitatively shown that if the viscosity function of a generalised Newtonian fluid is assumed to depend on the third invariant of the rate of deformation tensor, there is an increase of the viscous diffusion terms, but especially, of the dissipation of turbulence kinetic energy by a factor equal to the Trouton ratio of the fluid, divided by the Trouton ratio of the solvent, thus indicating a possible way to improve rheological–turbulence modelling.     

A RESIDUAL-BASED UNRESOLVED-SCALE FINITE ELEMENT MODELLING FOR IMPLICT LARGE EDDY SIMULATION 1)

对应于湍流的大尺度与小尺度流场信息, 本文在有限元的框架下, 假设Navier-Stokes方程的解的形函数由大尺度和不可解尺度形函数叠加组成, 引入对应的权函数, 将Navier-Stokes方程的有限元变分形式分解为大尺度和不可解尺度系统. 根据不可解尺度系统, 构建基于Navier-Stokes大尺度方程残差的不可解尺度模型, 将其代入Navier-Stokes方程的大尺度系统, 进而数值求解大尺度系统得到Navier-Stokes方程的大尺度解. 该方法无需像传统的大涡模拟方法那样对方程的解进行过滤, 通过对形函数进行尺度分解实现解的尺度分解. 本文使用该方法的自编程序代码开展了槽道湍流的数值模拟. 通过与有限差分大涡模拟、DNS计算结果的比较, 发现在使用较少网格情况下该方法预测的平均流向速度在近壁区与DNS数据吻合, 在黏性外层略偏高; 该方法对雷诺应力预测偏低导致从流向向垂向方向上湍动能输运略偏低. 流向速度等值面图显示该方法有效捕捉到了大尺度旋涡结构; 同时在近壁区可以观察到明显的低速条带结构.     

基于Navier-Stokes方程残差的隐式大涡模拟有限元模型

对应于湍流的大尺度与小尺度流场信息, 本文在有限元的框架下, 假设Navier-Stokes方程的解的形函数由大尺度和不可解尺度形函数叠加组成, 引入对应的权函数, 将Navier-Stokes方程的有限元变分形式分解为大尺度和不可解尺度系统. 根据不可解尺度系统, 构建基于Navier-Stokes大尺度方程残差的不可解尺度模型, 将其代入Navier-Stokes方程的大尺度系统, 进而数值求解大尺度系统得到Navier-Stokes方程的大尺度解. 该方法无需像传统的大涡模拟方法那样对方程的解进行过滤, 通过对形函数进行尺度分解实现解的尺度分解. 本文使用该方法的自编程序代码开展了槽道湍流的数值模拟. 通过与有限差分大涡模拟、DNS计算结果的比较, 发现在使用较少网格情况下该方法预测的平均流向速度在近壁区与DNS数据吻合, 在黏性外层略偏高; 该方法对雷诺应力预测偏低导致从流向向垂向方向上湍动能输运略偏低. 流向速度等值面图显示该方法有效捕捉到了大尺度旋涡结构; 同时在近壁区可以观察到明显的低速条带结构.      

粘流-无粘干扰流动(IF)理论

对不可压缩层流二维干扰流动,本文提出一个干扰流动(IF)理论。IF理论要点为:1)干扰流动沿主流的法向被分为三层即粘性层、干扰层和无粘层,引进了法向动量交换为主导过程的干扰层概念。2)利用力学守恒律、三层匹配关系及文中引进的干扰模型,把三层的空间尺度及惯性-粘性诸力的数置级表示为单参数m的函数,m<1/2·3)导出描述各层流动的控制方程、导出描述全城流动的控制方程为简化Navie-Stokes(SNS)方程。IF理论适用于不存在分离的附着干扰流动以及存在分离的大范围干扰流动,经典边界层(CBL)理论和流动分离局部区域Triple-Deck(TD)理论分别是本文理论在参数m=O和1/4时的两个特例,本文理论容易推广到可压缩、三维及湍流流动。     

Transonic viscous-inviscid interaction by a finite element method

A method is outlined for solving two-dimensional transonic viscous flow problems, in which the velocity vector is split into the gradient of a potential and a rotational component. The approach takes advantage of the fact that for high-Reynolds-number flows the viscous terms of the Navier-Stokes equations are important only in a thin shear layer and therefore solution of the full equations may not be needed everywhere. Most of the flow can be considered inviscid and, neglecting the entropy and vorticity effects, a potential model is a good approximation in the flow core. The rotational part of the flow can then be calculated by solution of the potential, streamfunction and vorticity transport equations. Implementation of the no-slip and no-penetration boundary conditions at the walls provides a simple mechanism for the interaction between the viscous and inviscid solutions and no extra coupling procedures are needed. Results are presented for turbulent transonic internal choked flows.     

Toward a Turbulence Constitutive Relation for Geophysical Flows

Rapidly rotating turbulent flows are frequently in approximate geostrophic balance. Single-point turbulence closures, in general, are not consistent with a geostrophic balance. This article addresses and resolves the possibility of a constitutive relation for single-point second-order closures for classes of rotating and stratified flows relevant to geophysics. Physical situations in which a geostrophic balance is attained are described. Closely related issues of frame-indifference, horizontal divergence, and the Taylor–Proudman theorem are discussed. It is shown that, in the absence of vortex stretching along the axis of rotation, turbulence is frame-indifferent. Unfortunately, no turbulence closures are consistent with this frame-indifference that is frequently an important feature of rotating or quasi-geostrophic flows. A derivation and discussion of the geostrophic constraint which ensures that the modeled second-moment equations are frame-invariant, in the appropriate limit, is given. It is shown that rotating, stratified, and shallow water flows are situations in which such a constitutive relation procedure is useful. A nonlinear nonconstant coefficient representation for the rapid-pressure strain covariance appearing in the Reynolds stress and heat flux equations, consistent with the geostrophic balance, is described. The rapid-pressure strain closure features coefficients that are not constants determined by numerical optimization but are functions of the state of turbulence as parametrized by the Reynolds stresses and the turbulent heat fluxes as is required by tensor representation theory. These issues are relevant to baroclinic and barotropic atmospheric and oceanic flows. The planetary boundary layers in which there is a transition, with height or depth, from a thermally or shear driven turbulence to a geostrophic turbulence is a classic geophysical example to which the considerations in this article are relevant. Received 14 October 1996 and accepted 9 June 1997     

Calculation of Reynolds stresses in turbulent supersonic flows

The baseline numerical procedure of interest in this study combines flux vector splitting, flux difference splitting and an explicit treatment of the diffusion terms of the flow equations. The viscous terms are treated explicitly to preserve the wave propagation properties of the Euler fluxes and permit splitting. The experience with this scheme has been limited to laminar or, at best, ‘eddy viscosity’ flows. In this paper the applicability of the scheme is extended to include the calculation of turbulent Reynolds stresses in supersonic flows. The schemes and our implementation are discussed. Both laminar and turbulence subsets of the Reynolds/Favre-averaged equations are tested, with a discussion of relative performance. The test problem for turbulence consists of a zero-pressure-gradient supersonic boundary layer as well as a supersonic boundary layer experiencing the combined effects of adverse pressure gradient, bulk compression and a concave streamline curvature. Excellent agreement with experimental measurements is observed for most of the quantities compared, which suggests that the numerical procedures presented in this paper are potentially very useful.     

气体动理学格式与多尺度流动模拟

介绍了气体动理学格式(GKS)的基本构造原理及其在两种典型多尺度流动模拟中的应用。GKS利用介观BGK方程的跨尺度演化解来构造网格界面上的数值通量,从而发展出能随计算网格尺度变化自动切换物理模型的多尺度方法。对湍流这种宏观多尺度流动,发展了高精度GKS方法并成功用于低雷诺数湍流的直接数值模拟;为实现对高雷诺数湍流的高效精细模拟,基于拓展BGK方程和已有的RANS,LES模型建立了新型多尺度模拟框架。对跨流域稀薄流动,发展了适合大规模并行的三维统一气体动理学格式(UGKS),并建立了适合轴对称稀薄流动的UGKS。研究表明,GKS在多尺度流动高效模拟中的优异性能,具有很好的发展前景。     

A segregated explicit algebraic structure-based model for wall-bounded turbulent flows

The Algebraic Structure-based Model (ASBM) is an engineering model of turbulence that incorporates information about the morphology of turbulent structures, and thus characterizes additional turbulence information beyond what is provided solely by the Reynolds stresses. In this work, a new formulation of the ASBM is developed that exhibits three defining characteristics. The first is a near-wall correction segregated from the homogeneous component of the model, which allows for a new paradigm for model development and comparison. The second characteristic is a set of fully-explicit equations that replace the original implicit formulation and make evident the highly nonlinear nature of the model. Finally, the last defining characteristic is based on the coupling of the ASBM with turbulent transport equations that provide the required turbulent scales. Using arguments of consistency, stability, and accuracy we arrive at a model that uses transport equations for k and ϵ. The resulting model is simpler to implement, more robust, and amenable to further development. The analysis of the new ASBM formulation is carried out by computing fully-developed channel, pipe and squared duct flows.     

Control of the aerodynamic characteristics of a profile oscillating with respect to the incidence angle in subsonic viscous flow

The results of a numerical simulation of the unsteady subsonic viscous gas flow around a two-dimensional profile oscillating with respect to the incidence angle are presented and the possibility of controlling the nonstationary aerodynamic characteristics is considered. The hysteresis phenomena typical of oscillatory profile motions are investigated, the dependence of the lift force and drag is found for various laws of periodic variation of the incidence angle with time, and the effect of the frequency and amplitude of the angular profile oscillations on the shape of the hysteresis curves is studied. The calculations were based on the numerical solution of the nonstationary Navier-Stokes equations averaged in the Reynolds sense (Reynolds equations) which were closed using the k-ω turbulence model with modeling of the laminar/turbulent transition.     

Resolution-optimised nonlinear scheme for secondary derivatives

A 5-point-stencil optimised nonlinear scheme with spectral-like resolution within the whole wave number range for secondary derivatives is devised. The proposed scheme can compensate for the dissipation deficiency of traditional linear schemes and suppress the spurious energy accumulation that occurs at high wave numbers, both of which are frequently encountered in large eddy simulation. The new scheme is composed of a linear fourth-order central scheme term and an artificial viscosity term. These two terms are connected by a nonlinear weight. The proposed nonlinear weight is designed based on Fourier analysis, rather than Taylor analysis, to guarantee a spectral-like resolution. Moreover, the accuracy is not affected by the optimisation, and the new scheme reaches fourth-order accuracy. The new scheme is tested numerically using the one-dimensional diffusion problem, one-dimensional steady viscous Burger’s shock, two-dimensional vortex decaying, three-dimensional isotropic decaying turbulence and fully developed turbulent channel flow. All the tests confirm that the new scheme has spectral-like resolution and can improve the accuracy of the energy spectrum, dissipation rate and high-order statistics of turbulent flows.     

Subgrid-Scale Models for Compressible Large-Eddy Simulations

An a priori study of subgrid-scale (SGS) models for the unclosed terms in the energy equation is carried out using the flow field obtained from the direct simulation of homogeneous isotropic turbulence. Scale-similar models involve multiple filtering operations to identify the smallest resolved scales that have been shown to be the most active in the interaction with the unresolved SGSs. In the present study these models are found to give more accurate prediction of the SGS stresses and heat fluxes than eddy-viscosity and eddy-diffusivity models, as well as improved predictions of the SGS turbulent diffusion, SGS viscous dissipation, and SGS viscous diffusion.     

Effects of free‐stream turbulence on large‐scale coherent structures of separated boundary layer transition

It has been well established that large‐scale structures, usually called coherent structures, exist in many transitional and turbulent flows. The topology and range of scales of those large‐scale structures vary from flow to flow such as counter‐rotating vortices in wake flows, streaks and hairpin vortices in turbulent boundary layer. There has been relatively little study of large‐scale structures in separated and reattached transitional flows. Large‐eddy simulation (LES) is employed in the current study to investigate a separated boundary layer transition under 2% free‐stream turbulence on a flat plate with a blunt leading edge. The Reynolds number based on the inlet free stream velocity and the plate thickness is 6500. A dynamic subgrid‐scale model is employed to compute the subgrid‐scale stresses more accurately in the current transitional flow case. Flow visualization has shown that the Kelvin–Helmholtz rolls, which have been so clearly visible under no free‐stream turbulence (NFST) are not as apparent in the present study. The Lambda‐shaped vortical structures which can be clearly seen in the NFST case can hardly be identified in the free‐stream turbulence (FST) case. Generally speaking, the effects of free‐stream turbulence have led to an early breakdown of the boundary layer, and hence increased the randomization in the vortical structures, degraded the spanwise coherence of those large‐scale structures. Copyright © 2005 John Wiley & Sons, Ltd.     

Study on anisotropic buoyant turbulence model

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

A K-ε TWO-EQUATION TURBULENCE MODEL FOR THE SOLID-LIQUID TWO-PHASE FLOWS

A two-equation turbulence model has been dereloped for predicting two-phase flow the two equations describe the conserration of turbulence kinetic energy and dissipation rate of that energy for the incompressible carrier fluid in a two-phase flow The continuity, the momentum, K and ε equations are modeled. In this model,the solid-liquid slip veloeites, the particle-particte interactions and the interactions between two phases are considered,The sandy water pipe turbulent flows are sueeessfuly predicted by this turbulince model.