关于湍流标度律的争鸣

研究湍流结构函数的标度律。实验或数值模拟得到的湍流结构函数的标度指数是奇异的。很多学者认为:这一实验事实否定Kolmogorov1941年(K41)提出的正常标度律,各向同性湍流惯性区的标度律是奇异的。近年来作者发表一系列文章,提出不同的观点:由于有限雷诺数效应,有限雷诺数湍流的标度指数不等于真正的惯性区标度指数,湍流结构函数的标度指数的实验数据并不否定K41正常标度律,各向同性湍流惯性区的标度律可能是正常的。惯性区奇异标度律和正常标度律对应的湍流物理本质是完全不同的,因而研究解决这个争论具有重要的意义。      

Jet Noise: Since 1952

Jet noise research was initiated by Sir James Lighthill in 1952. Since that time, the development of jet noise theory has followed a very tortuous path. This is, perhaps, not surprising for the understanding of jet noise is inherently tied to the understanding of turbulence in jet flows. Even now, our understanding of turbulence is still tenuous. In the fifties, turbulence was regarded as consisting of a random assortment of small eddies. As a result, the primary focus of jet noise research was to quantify the noise from fine-scale turbulence. This line of work persisted into the eighties. The discovery of large turbulence structures in free shear flows in the early seventies led some investigators to begin questioning the validity of the then established theories. Some went further to suggest that, for high-speed jets, it was the large turbulence structures/instability waves of the flow that were responsible for the dominant part of jet mixing noise. Development of a quantitative theory of noise from large turbulence structures/instability waves took place during the next 15 years. Precision instrumentation and facilities for jet noise measurements became available in the mid-eighties. This permitted a large bank of high-quality narrow band jet noise data to be gathered over the subsequent years. Recent analysis of these data has provided irrefutable evidence that jet noise, in fact, is made up of two basic components; one from the large turbulence structures/instability waves, the other from the fine-scale turbulence. This is true even for subsonic jets. In this paper, some of the crucial research results of the past 44 years, that form the basis of our present understanding of jet noise generation and propagation, are discussed. Received 14 November 1996 and accepted 14 March 1997     

旋涡、湍流与自由表面的相互作用

旋涡及自由表面湍流与自由表面的非定常、非线性相互作用是当前流体力学中一个十分活跃的前沿领域的研究课题,它具有深刻的理论意义和重要的实用背景．旋涡、湍流与自由表面的相互作用是一个非常复杂的非线性、非定常过程,涉及到涡－涡、波－涡相互作用,旋涡在自由表面处的断裂、重联、合并及自由表面湍流中的准拟序结构的形成、发展等复杂动力学过程．关于这一领域的研究已经取得了一些重要成果．对旋涡与自由表面的相互作用的基本过程有了较为清楚的认识;但对湍流与自由表面的相互作用中的许多基本物理过程还知之甚少．本文综述此领域的实验、理论和数值模拟研究的最新进展和主要结果,讨论若干有待解决的问题,并指出进一步的研究方向      

Thirty years of turbulence study in China

This paper provides a comprehensive introduction to the major progress in transition and turbulence studies in China over the past thirty years. Despite the author's best efforts, there are some unavoidable omissions. We expect that with the continued deepening of turbulence research, Chinese scholars will make increasing contributions in this field, which will certainly have strong influence on the establishment of a China School of Turbulence.     

Using a Fuzzy Piecewise Regression Analysis to Predict the Nonlinear Time-Series of Turbulent Flows with Automatic Change-Point Detection

Research has already shown that turbulent flow consists of some coherent time- and space-organized vortical structures. Some dynamic systems and experimental models are employed to understand the turbulent generation mechanism. However, these approaches still cannot provide a good nonlinear analysis of turbulent time-series. In the real turbulent flow, very complicated nonlinear behaviors, which are affected by many vague factors are present. Based on the nonlinear behavior and the results of from this traditional research, we introduce multivariate statistical analysis of an experimental study to explain practical phenomenon. In this paper, a new approach of fuzzy piecewise regression analysis with automatic change-point detection is proposed to predict the nonlinear time-series of turbulent flows. In order to show the practicality and usefulness of this model, we present an example of predicting the near-wall turbulence time-series as a verifiable model. The results of practical applications show that the proposed method is appropriate and appears to be useful in nonlinear analysis and in fuzzy environments to predict the turbulence time-series. This revised version was published online in July 2006 with corrections to the Cover Date.     

高雷诺数壁湍流的研究进展及挑战

高雷诺数壁湍流(high Reynolds number wall-bounded turbulence,HRNWT)是目前湍流科学研究的一个热点也是一个难点,对其现象、规律及机制的认知不足,理论体系远未建立而且研究手段受到各种限制.本文基于对HRNWT主要研究手段的介绍,针对HRNWT中的湍流统计量、超大尺度结构(very large scale motions,VLSMs)的尺度和形态以及起源和影响及其与颗粒的相互作用,总结了HRNWT的研究现状和最新进展,特别梳理了近年来本文作者团队在HRNWT特别是高雷诺数颗粒两相壁湍流方面的研究成果,并对HRNWT的进一步研究给出了建议及展望.      

Experiments in Turbulent Pipe Flow with Polymer Additives at Maximum Drag Reduction

In this paper we report on (two-component) LDV experiments in a fully developed turbulent pipe flow with a drag-reducing polymer (partially hydrolyzed polyacrylamide) dissolved in water. The Reynolds number based on the mean velocity, the pipe diameter and the local viscosity at the wall is approximately 10000. We have used polymer solutions with three different concentrations which have been chosen such that maximum drag reduction occurs. The amount of drag reduction found is 60–70%. Our experimental results are compared with results obtained with water and with a very dilute solution which exhibits only a small amount of drag reduction. We have focused on the observation of turbulence statistics (mean velocities and turbulence intensities) and on the various contributions to the total shear stress. The latter consists of a turbulent, a solvent (viscous) and a polymeric part. The polymers are found to contribute significantly to the total stress. With respect to the mean velocity profile we find a thickening of the buffer layer and an increase in the slope of the logarithmic profile. With respect to the turbulence statistics we find for the streamwise velocity fluctuations an increase of the root mean square at low polymer concentration but a return to values comparable to those for water at higher concentrations. The root mean square of the normal velocity fluctuations shows a strong decrease. Also the Reynolds (turbulent) shear stress and the correlation coefficient between the stream wise and the normal components are drastically reduced over the entire pipe diameter. In all cases the Reynolds stress stays definitely non-zero at maximum drag reduction. The consequence of the drop of the Reynolds stress is a large polymer stress, which can be 60% of the total stress. The kinetic-energy balance of the mean flow shows a large transfer of energy directly to the polymers instead of the route by turbulence. The kinetic energy of the turbulence suggests a possibly negative polymeric dissipation of turbulent energy. This revised version was published online in July 2006 with corrections to the Cover Date.     

可压缩各向同性衰减湍流直接数值模拟研究

采用五阶有限差分WENO格式直接模拟了高初始湍流Mach数的可压缩均匀各向同性湍流,主要分析了湍流的统计特性 和压缩性的影响,包括能谱特征、激波串、耗散率、标度律等. 研究表明,湍动能主要来自于速度场螺旋分量的贡献;各向同性湍流的小尺度脉动对压缩性更为敏感,并且压缩性的增强加快了湍流大 尺度脉动向小尺度脉动的湍动能输运;随着湍流Mach数的升高,胀量(压缩)耗散率所占比率也显著增长. 标度律分析表明,强可压缩湍流的横向速度结构函数仍然具有扩展自相似性;当阶数较高(p ≥ 5)时,纵向速度结构函数的扩展自相似性则不再成立. 对于压缩性较弱的湍流,与不可压缩湍流一致,横向湍流脉动的间歇性要强于纵向湍流脉动;而对于强可压缩湍流,纵向湍流脉动的 间歇性要强于横向湍流脉动.     

基于湍流类比的金属橡胶吸声特性定量分析

金属橡胶材料从表到里都具有大量互相贯通且混乱的孔隙, 经分析认为, 这种混乱性及一定的周期性与湍流中的不规则性和准周期性极其相似. 由于金属橡胶材料内部芯材结构的不规则性, 即便此时的雷诺数很小, 波在金属橡胶中的传播仍是以湍流流动为主. 因此, 引入Kolmogorov的关于湍流的局部各向同性概念, 同时对Kolmogorov关于湍流局部各向同性的两个假设进行类比. 从湍流物理模型出发, 借鉴了湍流的统计处理方法, 对金属橡胶材料的吸声特性进行定量分析, 得到金属橡胶材料的能量耗散率与其结构参数之间的表达式. 研究结果表明, 湍流统计方法的引入, 为基于金属橡胶材料的减振器、阻尼器、消声器等的优化设计提供了可靠的理论依据, 也为超轻金属多孔材料的工程应用提供了一种有效的定量分析方法.      

Behavior of homogeneous turbulence in channels with variable cross-sectional area

On the basis of the equation describing the behavior of the spectral tensor of the energy of the pulsation velocity for one-dimensional flow in a channel with variable cross-sectional area, we obtained a system of equations for the meansquare components of the pulsation velocity vector and integral scales of the turbulence length in various directions, enabling us to use these parameters in the initial section of the channel for determining their behavior along the channel. We made use of some ideas of A. N. Kolmogorov and J. Rotta concerning the possibility of describing viscous and nonlinear terms in the equations for the components of the tensor of Reynolds stresses in terms of the energy of pulsation motion and the integral scale of the turbulence length. The resulting system, in the special cases of very low intensity of turbulence, leads to the results of the linear theory; for constant cross-sectional area (or for a very high intensity of turbulence, when the damping affects the turbulence much more strongly than does the deformation effect) it describes the known empirical laws of the degeneration of turbulence beyond grids. We made a comparison with the data available in the literature on the behavior of the characteristics of turbulence in channels with variable cross-sectional area.     

弹性湍流及其对微混合效率影响研究的若干进展

介绍和评述弹性湍流的产生及其对于微混合效率的影响等问题上的若干研究进展. 弹性 湍流和惯性湍流具有类似的流场特征,但引发机理有所不同. 惯性湍流产生的原因 是惯性引起的Reynolds应力,而弹性湍流则是由弹性应力所引起的. 鉴于在微流动 中,惯性力可忽略不计,因此牛顿流体的混合变得十分困难. 此时可在流体中加入 微量高分子聚合物以生成黏弹性流体. 由黏弹性流体所引发的弹性湍流在提高流体 微混合的效率上可发挥重要作用.     

A parametric study of adverse pressure gradient turbulent boundary layers

There are many open questions regarding the behaviour of turbulent boundary layers subjected to pressure gradients and this is confounded by the large parameter space that may affect these flows. While there have been many valuable investigations conducted within this parameter space, there are still insufficient data to attempt to reduce this parameter space. Here, we consider a parametric study of adverse pressure gradient turbulent boundary layers where we restrict our attention to the pressure gradient parameter, β, the Reynolds number and the acceleration parameter, K. The statistics analyzed are limited to the streamwise fluctuating velocity. The data show that the mean velocity profile in strong pressure gradient boundary layers does not conform to the classical logarithmic law. Moreover, there appears to be no measurable logarithmic region in these cases. It is also found that the large-scale motions scaling with outer variables are energised by the pressure gradient. These increasingly strong large-scale motions are found to be the dominant contributor to the increase in turbulence intensity (scaled with friction velocity) with increasing pressure gradient across the boundary layer.     

Turbulent thermal diffusion in a multi-fan turbulence generator with imposed mean temperature gradient

We studied experimentally the effect of turbulent thermal diffusion in a multi-fan turbulence generator which produces a nearly homogeneous and isotropic flow with a small mean velocity. Using particle image velocimetry and image processing techniques, we showed that in a turbulent flow with an imposed mean vertical temperature gradient (stably stratified flow) particles accumulate in the regions with the mean temperature minimum. These experiments detected the effect of turbulent thermal diffusion in a multi-fan turbulence generator for relatively high Reynolds numbers. The experimental results are in compliance with the results of the previous experimental studies of turbulent thermal diffusion in oscillating grid turbulence (Buchholz et al. 2004; Eidelman et al. 2004). We demonstrated that the turbulent thermal diffusion is an universal phenomenon. It occurs independently of the method of turbulence generation, and the qualitative behavior of particle spatial distribution in these very different turbulent flows is similar. Competition between turbulent fluxes caused by turbulent thermal diffusion and turbulent diffusion determines the formation of particle inhomogeneities.     

An investigation of forced structures in turbulent jet flows

In many cases, turbulence is superimposed on an unsteady organized motion of a mean flow. In the past, these turbulent flows have been studied by time or ensemble averaging methods and some decomposition techniques such as proper orthogonal decomposition (POD). In this study, a new decomposition technique called the turbulence filter will be used to decompose the forced turbulent jet flows. By using the turbulence filtering technique, the fluctuating (turbulent) part and the more organized (forced) part of the velocity field are analyzed. Within this context several experiments on organized turbulent jet flow have been carried out. In the experiments, variable frequency and amplitude oscillation are imposed on a 1D jet. An elliptical plate was used in order to obtain sinusoidal forcing. The axial distance, Reynolds number and the forcing frequency of the signal were varied. The multiple hot wires (six probes) were used to investigate the evolution of the signal along the radial distance. The obtained results of the turbulence filter are compared with those of phase-averaging and POD techniques. The eigenmodes of the data are also evaluated by using the POD method. Received: 31 July 1998/Accepted: 19 January 2000     

The Relevance of a Back-Scatter Model for Depth-Averaged Flow Simulation

This study demonstrates the importance of a sophisticated sub-grid model when performing a depth-averaged unsteady RANS simulation of a shallow flow. The reduction of resolution and the spatial dimensions exclude important physical processes as present in three-dimensional turbulence. Especially the effect of the bottom turbulence on the formation of horizontal eddies appears of key importance. A method is proposed to incorporate these effects by means of a kinematic simulation that mimics the residual turbulent fluctuations in a straight channel flow after depth-averaging. This method is developed in the context of the evolution of large eddies in a shallow mixing layer. A comparison with experiments shows that the proposed method works satisfactory. Naturally, it does not fully account for the omission of all 3D-effects.     

The artificial buffer layer and the effects of forcing in hybrid LES/RANS

The present study is focused on large eddy simulations (LES) that use a statistical (RANS) turbulence model near solid walls, and on the artificial buffer layer that is formed at the interface between these two modeling regions. Additional forcing is used to trigger resolved motions in the LES region more quickly, and leads to improved results in several ways. The study investigates the artificial buffer layer and how it changes with the use of forcing in an in-depth manner, with the purpose of increased understanding of the increasingly popular hybrid LES/RANS group of methods.

The artificial buffer layer is shown to extend from below the modeling interface to well above it, in fact up to 20% of the boundary layer thickness for the cases studied here. The artificial buffer layer is found to be similar to the true buffer layer in many aspects, including a high correlation between the streamwise and wall normal velocity components in the ‘superstreaks’. This indicates that while the superstreaks are highly anisotropic and have unphysical length scales, they still contribute to the resolved shear stress. The forcing does not remove the artificial buffer layer, but it does reduce its extent and increases the resolved shear stress. This increase is mainly associated with increased fluctuations of the wall normal velocity.

A simple, low-dimensional forcing model is proposed and tested, with favorable results. The model is simple to implement and easily generalized to more complex geometries.     

Flow through a diverging channel: instability and bifurcation

The flow of a viscous fluid along a diverging channel is important practically. The mechanics of the flow also poses many theoretical problems. Some of these problems were considered before the foundation of the Japan Society of Fluid Mechanics, so it is perhaps surprising that their solutions are still so poorly understood. Nonetheless, a revival of interest in the last 20 years has led to progress in our understanding of the early stages of transition to turbulence in a two-dimensional channel. The present paper is a review of this progress in the theory of instabilities and bifurcations as the flux driven steadily along the channel is increased.     

Dissipation rate correction methods

Dissipation rates of the turbulent kinetic energy and of the scalar variance are underestimated when the measurement resolution of the small scales of a turbulent flow field are insufficient. Results are presented of experiments conducted in a salt-stratified water tunnel (Schmidt number ∼700). Dissipation rates are determined to be underestimated, and thus correction techniques based on velocity structure functions and mixed-moment functions are proposed. Dissipation rates in laboratory experiments of shear-free, grid-generated turbulence are determined from balance calculations of the kinetic energy and scalar variance evolution equations. Comparisons between the structure function and balance estimates of dissipation show that the corrections are O(1) for the kinetic energy dissipation rate, and are O(100) for the scalar variance dissipation rate. This difference is due to the lack of resolution down to the Batchelor scales that is required for a high Schmidt number flow. Simple correction functions based on microscale Reynolds numbers are developed for both turbulent kinetic energy and scalar variance dissipation rates. Application of the technique to the results of laboratory experiments of density stratified turbulence, sheared turbulence, and sheared density stratified turbulence yields successful corrections. It is also demonstrated that the Karman–Howarth equality (and the analogous Yaglom equation) that relates second and third-order structure functions to dissipation rates is valid for both unstrained (decaying grid-generated turbulence) and density stratified and sheared turbulence at least up to the magnitudes of strains of the current experiments Nt∼10, St∼10, respectively. This is helpful for it allows the use of these equations in the analysis of turbulence even when the large scale background profiles of velocity and scalar are unknown.     

Detached-Eddy Simulations Past a Circular Cylinder

The flow is calculated with laminar separation (LS) at Reynolds numbers 50,000 and 140,000, and with turbulent separation (TS) at140,000 and 3 × 10₆. The TS cases are effectively tripped, but compared with untripped experiments at very high Reynolds numbers. The finest grid has about 18,000 points in each of 56 grid planes spanwise; the resolution is far removed from Direct Numerical Simulations, and the turbulence model controls the separation if turbulent. The agreement is quite good for drag, shedding frequency, pressure, and skin friction. However the comparison is obscured by large modulations of the vortex shedding and drag which are very similar to those seen in experiments but also, curiously, durably different between cases especially of the LS type. The longest simulations reach only about 50 shedding cycles. Disagreement with experimental Reynolds stresses reaches about 30%, and the length of the recirculation bubble is about double that measured. The discrepancies are discussed, as are the effects of grid refinement, Reynolds number, and a turbulence-model curvature correction. The finest grid does not give the very best agreement with experiment. The results add to the validation base of the Detached-Eddy Simulation (DES) technique for smooth-surface separation. Unsteady Reynolds-averaged simulations are much less accurate than DES for LS cases, but very close for TS cases. Cases with a more intricate relationship between transition and separation are left for future study. This revised version was published online in July 2006 with corrections to the Cover Date.     

Experimental study on oscillating grid turbulence and free surface fluctuation

This paper analyses the interaction between the turbulence and free surface. The phenomenon takes place in many natural flows and industrial processes. In the present experiments, turbulence is generated by a vertically oscillating grid moving beneath the free surface. Fluid velocity has been measured through a hot-film anemometer, and the free surface elevation has been measured by an ultrasonic sensor. Integral length scales and several turbulence estimators have been computed. In order to detect the generation of turbulence near the free surface, the correlation between free surface elevation and the underneath flow velocity has been studied, as well as the time lag between turbulence and free surface. The free surface dynamics has been characterized by a velocity scale and a length scale. The kinetic energy associated with the free surface fluctuations increases with the Reynolds number at a rate depending on the frequency of the grid movement. For Reynolds number larger than ≈1000, however, the relationships collapse to a single curve characterized by a lower rate. The present experiments do not achieve the inertial sub-range in the vertical velocity fluctuations, and the estimated spectrum decays with an exponent smaller than ?3, which is the typical value for the two-dimensional turbulence in the inertial sub-range. The macro length scale, estimated by using the Taylor’s frozen turbulence hypothesis, experiences a decay away from the grid, which follows reasonably well the profile of Thompson and Turner (J Fluid Mechanics 67: 349–368, 1975). The micro length scale reduces immediately beneath the free surface, which can be interpreted by the increase of dissipation rate in the subsurface layer. The classification diagram by Brocchini and Peregrine (J Fluid Mech 449: 225–254, 2001) indicates that most tests fall in the weak turbulence domain, but some tests fall in the wavy domain. The vertical velocity fluctuations and the free surface level show a significant correlation with a negative phase lag, that is, the free surface fluctuations are ahead of the vertical velocity fluctuations.