A stability condition for turbulence model： From EMMS model to EMMS-based turbulence model

The closure problem of turbulence is still a challenging issue in turbulence modeling. In this work, a stability condition is used to close turbulence. Specifically, we regard single-phase flow as a mixture of turbulent and non-turbulent fluids, separating the structure of turbulence. Subsequently, according to the picture of the turbulent eddy cascade, the energy contained in turbulent flow is decomposed into different parts and then quantified. A turbulence stability condition, similar to the principle of the energy-minimization multi-scale （EMMS） model for gas-solid systems, is formulated to close the dynamic constraint equa- tions of turbulence, allowing the inhomogeneous structural parameters of turbulence to be optimized. We name this model as the ＂EMMS-based turbulence model＂, and use it to construct the corresponding turbulent viscosity coefficient. To validate the EMMS-based turbulence model, it is used to simulate two classical benchmark problems, lid-driven cavity flow and turbulent flow with forced convection in an empty room, The numerical results show that the EMMS-hased turbulence model improves the accuracy of turbulence modeling due to it considers the principle of compromise in competition between viscosity and inertia.     

A stability condition for turbulence model: From EMMS model to EMMS-based turbulence model

The closure problem of turbulence is still a challenging issue in turbulence modeling. In this work, a stability condition is used to close turbulence. Specifically, we regard single-phase flow as a mixture of turbulent and non-turbulent fluids, separating the structure of turbulence. Subsequently, according to the picture of the turbulent eddy cascade, the energy contained in turbulent flow is decomposed into different parts and then quantified. A turbulence stability condition, similar to the principle of the energy-minimization multi-scale (EMMS) model for gas–solid systems, is formulated to close the dynamic constraint equations of turbulence, allowing the inhomogeneous structural parameters of turbulence to be optimized. We name this model as the “EMMS-based turbulence model”, and use it to construct the corresponding turbulent viscosity coefficient. To validate the EMMS-based turbulence model, it is used to simulate two classical benchmark problems, lid-driven cavity flow and turbulent flow with forced convection in an empty room. The numerical results show that the EMMS-based turbulence model improves the accuracy of turbulence modeling due to it considers the principle of compromise in competition between viscosity and inertia.     

Non-equilibrium turbulent phenomena in the ?ow over a backward-facing ramp

Non-equilibrium turbulence phenomena have raised great interests in recent years. Significant efforts have been devoted to non-equilibrium turbulence properties in canonical flows, e.g., grid turbulence, turbulent wakes, and homogeneous isotropic turbulence(HIT). The non-equilibrium turbulence in non-canonical flows, however, has rarely been studied due to the complexity of the flows. In the present contribution, a directnumerical simulation(DNS) database of a turbulent flow is analyzed over a backwardfacing ramp, the flow near the boundary is demonstrated, and the non-equilibrium turbulent properties of the flow in the wake of the ramp are presented by using the characteristic parameters such as the dissipation coefficient C and the skewness of longitudinal velocity gradient Sk, but with opposite underlying turbulent energy transfer properties. The equation of Lagrangian velocity gradient correlation is examined, and the results show that non-equilibrium turbulence is the result of phase de-coherence phenomena, which is not taken into account in the modeling of non-equilibrium turbulence. These findings are expected to inspire deeper investigation of different non-equilibrium turbulence phenomena in different flow conditions and the improvement of turbulence modeling.     

磁流体湍流及数值模拟研究综述

自从20世纪90年代初,AJAX计划的提出掀起了磁流体技术的研究热潮.本文总结了国内外关于磁流体湍流研究的成果,包括两方面:磁流体湍流的特点和数值模型.相比于普通湍流,磁流体湍流表现出电磁诱导的各向异性、线性焦耳耗散、湍流重层流化等特性,本文同时对磁流体湍流的结构和脉动抑制现象进行了描述:数值模型研究从雷诺平均方程和大涡模拟两方面作了介绍.当前研究集中在不可压缩磁流体湍流,距离能够应用于航空航天问题还有很长的路要走.     

二维衰减湍流的速度加速度结构函数

湍流场中二阶速度加速度结构函数 (velocity-acceleration structure function, VASF) 被认为与尺度间能量或者拟涡能的传递相关,其正负表明传递的方向. 三维湍流中,能量从大尺度向 小尺度传递,VASF 为负. 在二维湍流中,能量反向传递到大尺度,拟涡能正向传递到小尺度,因此理论上 VASF 无论在反向能量级串区还是在正向拟 涡能级串区均为正. 然而,相对于三维湍流中 VASF 的充分研究,二维湍流中 VASF 的正负性迄今尚无实验或数值模拟数据验证. 本文通过三维二维湍流中普适的公式推导,指出在空间非均匀湍流场中,VASF 除了尺度间传递,还受到非均匀项的影响. 一种常见的空间非均匀湍流场是在实验研究中常用的风洞或水洞中,湍流发生装置 (如栅格) 后的湍流. 该流场中,湍流强度随下游位置增大而逐渐衰减,这种衰减则带来空间上的非均匀性. 本文在基于竖直流动皂膜的二维衰减湍流场中,利用拉格朗日粒子追踪法测得在拟涡能级串区的 VASF,并分析各部分的影响. 结果表明,虽然尺度间传递项为正值,但由于衰减带来的非均匀项为负值,使得 VASF 的值为负,使之失去了表征拟涡能传递方向的意义. 因此,在类似风洞、水洞、水槽等衰减流场中对 VASF 的讨论不应忽略非均匀项. 最后对与 VASF 密切相关的弥散过程进行分析,发现后期弥散过程变慢是由于负的 VASF 导致.      

VELOCITY-ACCELERATION STRUCTURE FUNCTION IN TWO-DIMENSIONAL DECAYING TURBULENCE$^{\bf 1)}$

湍流场中二阶速度加速度结构函数 (velocity-acceleration structure function, VASF) 被认为与尺度间能量或者拟涡能的传递相关,其正负表明传递的方向. 三维湍流中,能量从大尺度向 小尺度传递,VASF 为负. 在二维湍流中,能量反向传递到大尺度,拟涡能正向传递到小尺度,因此理论上 VASF 无论在反向能量级串区还是在正向拟 涡能级串区均为正. 然而,相对于三维湍流中 VASF 的充分研究,二维湍流中 VASF 的正负性迄今尚无实验或数值模拟数据验证. 本文通过三维二维湍流中普适的公式推导,指出在空间非均匀湍流场中,VASF 除了尺度间传递,还受到非均匀项的影响. 一种常见的空间非均匀湍流场是在实验研究中常用的风洞或水洞中,湍流发生装置 (如栅格) 后的湍流. 该流场中,湍流强度随下游位置增大而逐渐衰减,这种衰减则带来空间上的非均匀性. 本文在基于竖直流动皂膜的二维衰减湍流场中,利用拉格朗日粒子追踪法测得在拟涡能级串区的 VASF,并分析各部分的影响. 结果表明,虽然尺度间传递项为正值,但由于衰减带来的非均匀项为负值,使得 VASF 的值为负,使之失去了表征拟涡能传递方向的意义. 因此,在类似风洞、水洞、水槽等衰减流场中对 VASF 的讨论不应忽略非均匀项. 最后对与 VASF 密切相关的弥散过程进行分析,发现后期弥散过程变慢是由于负的 VASF 导致.     

Inlet grid-generated turbulence for large-eddy simulations

A new technique of generating turbulence in large-eddy simulations (LES) has been investigated and results compared with previous studies for validation. The proposed gridInlet technique uses a grid pattern on the inlet boundary patch to produce grid-generated turbulence as used in wind tunnel experiments. This allows the turbulence integral length scale to be controlled by changing the grid size, while the turbulence intensity is controlled by changing the inlet distance. The objective of this paper is to investigate domain and mesh requirements to implement the gridInlet technique. This technique is most suited to studies on the influence of high-intensity isotropic turbulence on objects, particularly if comparisons are to be made to experimental data obtained with grid-generated turbulence.     

Turbulence Intensity and Length Scale Effects on Premixed Turbulent Flame Propagation

The effects of varying turbulence intensity and turbulence length scale on premixed turbulent flame propagation are investigated using Direct Numerical Simulation (DNS). The DNS dataset contains the results of a set of turbulent flame simulations based on separate and systematic changes in either turbulence intensity or turbulence integral length scale while keeping all other parameters constant. All flames considered are in the thin reaction zones regime. Several aspects of flame behaviour are analysed and compared, either by varying the turbulence intensity at constant integral length scale, or by varying the integral length scale at constant turbulence intensity. The turbulent flame speed is found to increase with increasing turbulence intensity and also with increasing integral length scale. Changes in the turbulent flame speed are generally accounted for by changes in the flame surface area, but some deviation is observed at high values of turbulence intensity. The probability density functions (pdfs) of tangential strain rate and mean flame curvature are found to broaden with increasing turbulence intensity and also with decreasing integral length scale. The response of the correlation between tangential strain rate and mean flame curvature is also investigated. The statistics of displacement speed and its components are analysed, and the findings indicate that changes in response to decreasing integral length scale are broadly similar to those observed for increasing turbulence intensity, although there are some interesting differences. These findings serve to improve current understanding of the role of turbulence length scales in flame propagation.

     

Experimental study of 3D turbulence by solving an inverse problem

The problem of the investigation of 3D turbulence from the projected data is studied here. It is shown that statistical information about turbulence can be recovered in such situations through a solution of an inverse problem. Unlike earlier papers, the adopted model of the turbulence envisages the turbulence characteristics to vary inside the 3D volume. Analytical solutions expressing turbulence characteristics (correlation functions, spectra) in terms of experimental data are found for this realistic model.     

On the effects of free-stream turbulence on axisymmetric disc wakes

Wind tunnel experiments have been used to study the effects of free-stream turbulence on the axisymmetric wake behind a disc. The disc and its wake were introduced to various turbulent streams having various levels of turbulence intensity and length scale. It was found that the presence of free-stream turbulence enhances the body??s drag and hence wake momentum deficit, if it is of sufficient strength, changes the far wake??s decay rate and prevents the appearance of self-similarity. The external turbulence causes a significant transformation in the wake??s turbulence structure. This gradually evolves towards the character of the free-stream turbulence itself and thus is characterised by much weaker turbulence (cross-stream) transport processes and a consequent dominance of shear stress production, which acts to maintain the shear stress and mean velocity profiles.     

Evaluation of bubble-induced turbulence using direct numerical simulation

The presented research evaluates the interaction between a single bubble and homogeneous turbulent flow using direct numerical simulation (DNS) approach. The homogeneous single-phase turbulence is numerically generated by passing a uniform flow through grid planes. The turbulence decay rate is compared with experiment-based correlation. The single phase turbulence is then used as an inflow boundary condition for a set of single bubble studies. By estimating the turbulent field around the fully resolved bubble, the effects of bubble deformability, turbulent intensity and relative velocity on the bubble-induced turbulence are investigated. The existence of bubble creates new vortices in the wake region and the enhancement of turbulence is observed in the region behind the bubble. The results show that the magnitude of the turbulence enhancement would increase as the bubble encounters larger liquid turbulent intensity or higher relative velocity. Set of bubble Weber numbers from 0.34 to 3.39 are used to investigate the effect of bubble deformability. The more deformable bubble is the higher the increase in the magnitude of the turbulence enhancement behind the bubble. This research provides systematic insight on the bubble-induced turbulence (BIT) mechanism and is important for multiphase computational fluid dynamics (M-CFD) closure model development.     

The influence of temperature on the measurements of Reynolds stresses in shear free turbulence near a wall

 Temperature changes have a significant influence on the measurements of Reynolds stresses in turbulent boundary layers. As compared to the spanwise velocity fluctuations the streamwise turbulence intensity is especially sensitive to temperature deviations. Although this is a general statement its importance is clearly elucidated in a shear-free turbulence near a solid wall, since the mixing due to turbulence production is minimized in this flow. A consequence of temperature influence on hot-wire measurements is that frictional heating from the wall has produced contradictory results in different experiments on shear-free turbulence. In the current paper, measurements of streamwise and spanwise turbulence intensities have been conducted at different wall temperatures, thereby simulating the contradictory results mentioned above. A simple model has been developed showing that the turbulence intensities are affected by both the rms. value of the temperature fluctuations and the correlation between fluctuating temperature and velocity. These correlations are measured and the developed model is used to explain deviations in earlier measurements on shear-free turbulence. Moreover, the individual magnitudes of the two correlations in the temperature correction are estimated and their individual importance is discussed. Received: 17 February 1997 / Accepted: 18 January 1998     

从光闪烁频谱反演大气湍流谱:原理与数值模拟

建立了一种数值算法,从光波在湍流大气中传播的闪烁时间频谱的高频特征反演大气光学湍流谱.此算法具有两点优越性:（1）反演出的湍流谱直接表征大气折射率的起伏而非大气温度场的起伏;（2）可以获得高空间波数、大部分位于耗散区的湍流谱,而用其它方法不易获得此区域的湍流谱.对Hill湍流谱以及一种复杂的湍流谱进行了大量数值模拟来检验此种算法,结果表明:如果知道大气风速与湍流内尺度就可以相当有效地获得湍流谱.     

Turbulence generation by electric body forces

The effect of an electric body force on vorticity production and turbulence generation in a gas is investigated by examination of the governing electrohydrodynamic equations. The theoretical concepts are illustrated by hot-film anemometer measurements of the electrically induced turbulence in a large scale electrostatic precipitator. The results indicate dramatic increases in turbulence and diffusivity due to the corona discharge and suggest that turbulence control in this application requires a modification in electrode geometry.This paper was presented at the Ninth symposium on turbulence, University of Missouri-Rolla, October 1–3, 1984     

Response of the Sphere Wake to Freestream Fluctuations

Direct numerical simulations have been used to investigate the response of the wake of a sphere to freestream fluctuations. This study has been motivated by the need to understand particle-induced turbulence enhancement in particulate flows. A sequence of simulations of flow past a sphere have been carried out where the frequency and amplitude of the freestream fluctuations and the flow Reynolds number has been varied systematically. It has been suggested that turbulence enhancement is primarily caused by vortex shedding from particles (Gore and Crowe, 1989; Hetsroni, 1989). Our simulations of the forced wake indicate that turbulence enhancement may be attributed to natural vortex shedding only when the freestream fluctuation level is low and the Reynolds number is greater than about 300. In addition to natural vortex shedding, the current simulations also suggest another mechanism for turbulence enhancement. It is found that in the presence of freestream fluctuations, the wake behaves like an oscillator and returns large amounts of kinetic energy to the surrounding fluid at resonance. This mechanism is not associated with natural vortex shedding and is therefore capable of enhancing freestream turbulence even at Reynolds numbers less than 300. Simulations also indicate that when the turbulence intensity of the carrier fluid is high, this resonance mechanism might be solely responsible for turbulence enhancement. Finally, our simulations also suggest a possible explanation for the correlation between turbulence enhancement and the ratio of the particle size to the size of energy containing eddies of turbulence found by Gore and Crowe (1989). Received 5 October 1999 and accepted 14 October 1999     

Constant Reynolds number turbulence downstream of an orificed perforated plate

Quasi-isotropic turbulence was experimentally produced in a wind tunnel via an orificed, perforated plate (OPP) at 10.5 m/s. The OPP consists of a lattice arrangement of 38.1 mm holes occupying 57% of the plate area. The OPP turbulence was found to be homogeneous over the cross section normal to the mean flow with Gaussian-like turbulence fluctuation. The isotropy of the turbulence field as portrayed by the streamwise/lateral turbulence intensity ratio was found to be approximately 1.1. The OPP turbulence is essentially self-preserving wherein the Taylor microscale Reynolds number remains nearly constant and the lateral velocity correlations collapse into a single curve.     

Modern turbulence and new challenges

The paper briefly reviews the progress in turbulence research in the 20th century and a number of issues are addressed based on achievements. The modern theroy of Navier-Stokes equation provides the theoretical basis for the development of turbulence research. The significance and bottle neck of DNS and the physical experiment in exploring turbulent flows are analyzed. The active manipulation of turbulence is directly guided by the knowledge of large-scale coherent structures. The existing problems in the large-eddy simulation are also pointed out. Scalar turbulence, which behaves quite different from fluid turbulence in many aspects, has drawn much attention in recent years. Besides the analysis of the difficulties in turbulence research, a number of examples are also presented to show how to use modern theory, computer and high technology to explore the nature of turbulence. The project supported by the National Natural Science Foundation of China (NSFC) (19572041 and 19732005)     

Measurement of velocity and kinetic energy of turbulence in swirling flows and their numerical prediction

Swirling flows are often employed in gas turbine combustion systems and high intensity industrial furnaces. A detailed analysis of the turbulence in the flow is necessary to achieve optimum combustion conditions. In this paper a method has been described to measure the turbulence levels in three directions using a hot wire anemometer. So far there is no established method available for measurement of turbulence in swirling and recirculating flows. The present method, it is hoped, will bridge the gap. The merit of the present method is the use of a single-wire probe rather than the X-probe. The method has been used for the measurement of turbulence levels in swirling recirculating flows generated by vane swirlers. From the measured turbulence levels, the kinetic energy of turbulence has been calculated and the results are compared with a well-established numerical prediction method. Mean velocity measurements have also been made using a 3-hole Pitot probe. The agreement between the measured and predicted values is quite satisfactory.     

THE EXAMINATION OF TURBULENCE MODELING WITH LES DATABASE

ＴＨＥＥＸＡＭＩＮＡＴＩＯＮＯＦＴＵＲＢＵＬＥＮＣＥＭＯＤＥＬＩＮＧＷＩＴＨＬＥＳＤＡＴＡＢＡＳＥＳｕＭｉｎｇ－ｄｅ（苏铭德）（ＱｉｎｇｈｕａＵｎｉｖｅｒｓｉｔｙ，Ｂｅｉｊｉｎｇ）Ｒ．Ｆｒｉｅｄｒｉｃｈ（ＬｅｈｒｓｉｒｈｌｆｕｒＦｌｕｉｄｓｍｅｃｈａ...     

INTERNAL LEE WAVES AND TURBULENCE MIXING OVER AN ISOLATED SEAMOUNT: RESULTS FROM TURBULENCE ENERGY MODELS

A three-dimensional, primitive equation, baroclinic numerical model incorporating a range of turbulence energy submodels is used to study the generation of internal lee waves over an isolated seamount. Attention is given to the turbulence mixing enhanced by the internal lee waves. The results show that regions of strong turbulence energy appear over the lee side of the seamount associated with the production of the lee waves. The computed vertical eddy viscosity and diffusivity using turbulence models can be as large as 1 m²s⁻¹. A comparison of the magnitude and spatial distribution of the internal lee waves does not reveal any major differences in results computed using different turbulence energy models or mixing determined from a Richards on number formulation. However, the magnitude of the vertical mixing is sensitive to the form of turbulence energy submodel. Also, a study of the relevant importance of the various terms in the turbulence energy equation shows that the term representing the advection of turbulence needs to be retained in order to accurately reproduce the mixing produced by the internal lee waves. Calculations using a range of seamount profiles show that the magnitude of near-bed turbulence is sensitive to the shape of the seamount. The magnitude and spatial distribution of the lee waves and associated flow field are affected by the parametrization of horizontal diffusion, with significant differences between the use of Laplacian and biharmonic forms of horizontal diffusion. The application of biharmonic horizontal diffusion is recommended.