）IRECT NUMERICAL SIMULATION OF DECAYING COMPRESSIBLE ISOTROPIC TURBULENCE

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

湍流能量耗散率标度律实验研究

应用粒子图像测速(PIV)系统对平板湍流边界层内流向和法向的瞬时速度进行了测量。湍流的能量耗散率由轴对称假设得到,同时在研究湍流动能耗散率标度律的过程中采用传统的统计学方法。实验结果显示,对于不同尺度上和不同法向位置湍流耗散率标度律来说,湍流耗散主要发生在小尺度上,也就是说湍动能耗散率标度律在小尺度上具有普适性。另外,根据层次结构理论假设,通过PIV实验数据对最高激发态的标度指数进行了研究,结果发现,最高激发态存在绝对标度指数,并且绝对标度律是由信号中最强耗散涨落的局部结构产生的。     

用子波变换检测壁湍流信号局部标度指数

本文用子波变换检测了刻画壁湍流脉动信号自相似性的局部标度指数，研究了不同尺度的湍流结构的自相似性，发现在湍流边界层猝发过程中，喷射和扫掠发生时刻小尺度脉动速度信号的局部标度指数为负值，说明在大尺度猝发事件发生的时刻小尺度结构具有奇异的自相似性，在猝发过程中其作用不仅仅是对湍能的耗散．     

自由端射流多尺度湍涡结构标度律的实验研究

利用后向接收式激光多普勒测速仪（ＬＤＶ）对自由湍射流进行了测量,我们对采集到的充分发展的湍流信号先求其１－８阶的结构函数及其标度指数,结果验证了Ｋｏｌｍｏｇｏｒｏｖ提出的标度律理论。然后用子波变换将自由湍流脉动速度分解为多尺度湍流结构,研究每一个尺度湍涡速度的结构函数的标度律及其与湍涡速度的自相关函数的关系。     

自由湍射流多尺度湍涡结构标度律的实验研究

利用后向接收式激光多普勒测速仪（ＬＤＶ）对自由湍射流进行了测量,我们对采集到的充分发展的湍流信号先求其１－ ８阶的结构函数及其标度指数,结果验证了Ｋｏｌｍ ｏｇｏｒｏｖ提出的标度律理论． 然后用子波变换将自由湍流脉动速度分解为多尺度湍涡结构, 研究每一个尺度湍涡速度的结构函数的标度律及其与湍涡速度的自相关函数的关系．     

壁湍流边界层奇异标度律的实验研究

采用热线风速仪对平板湍流边界层的流向速度进行测量,用速度结构函数研究不同尺度结构标度律的变化规律,结果显示小尺度区的概率密度曲线尾部明显偏离高斯型,说明高幅值间歇性事件占的份额较多;惯性子区的曲线向高斯型靠近,间歇性事件所占份额减少;大尺度结构的曲线趋于高斯型,间歇性事件所占份额最小。在耗散区、惯性子区和较大的尺度结构区存在大小不同的绝对标度指数,越靠近壁面这些区域的标度指数均越偏离p/3而逐渐变小。绝对标度指数与边界层位置有关,在缓冲层各阶标度指数与线性标度律偏差很大,显示较强的奇异性,当过渡到对数层及外区,标度指数逐渐增大,接近均匀各向同性湍流的状态。缓冲层、对数层及外区具有各异的绝对标度指数增长率,与各层的不同湍流结构特征和运动形式有关。     

金属丝壁面加热湍流标度律的实验研究

通过在固壁表面的平板湍流边界层沿流向平行放置若干通电加热的金属细丝,在平板表面形成沿展向周期性分布的温度场,利用该温度场引起的空气热对流,在湍流边界层近壁区域产生一组沿湍流边界层展向周期分布的流向涡结构。对壁湍流小尺度结构标度律统计特性的研究表明,金属丝加热后形成的规则流向涡结构将壁湍流各种尺度湍涡结构不规则的脉动有序地组织起来,增强了湍流小尺度结构的层次结构相似性,减小了壁湍流中小尺度结构的间歇性和奇异性,抑制了壁湍流中奇异的湍涡结构。     

关于湍流标度律的争鸣

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

壁湍流扩展的自相似标度律的实验研究

对风洞中零压力梯度平板湍流边界层进行了实验研究，用热线风速仪测量了不同法向位置的脉动速度，研究了湍流边界层不同法向位置速度结构函数的扩展的自相似标度律。     

壁面加热边界条件对湍流扩展的自相似性的影响

用热线风速仪测量了风洞中壁面加热和不加热平板湍流边界层不同法向位置的瞬时流向速度的时间序列,研究了壁面加热的边界条件对流向速度增量的p阶结构函数的扩展的自相似性和层次结构模型的影响。实验结果表明,壁面加热的边界条件对流向速度增量的结构函数扩展的自相似标度律的相对标度指数δ（p,3）和层次结构模型中间歇参数β、最奇异指数γ的影响明显存在。     

Statistical study of anisotropic anomalous behavior of a cylinder wake

 An experimental investigation of the statistical properties of the velocity difference in an anisotropic turbulent cylinder wake at moderate Re _λ is conducted by a triple hot-wire anemometer (HWA) probe. The energy spectra show a clear large scale anisotropy due to the presence of the Von Karman vortices. In spite of the low Re _λ and the large scale anisotropy, the extended-self-similarity (ESS) allows location of broad scaling ranges and the calculation of the scaling exponents of the three velocity structure function components in order to examine the intermittency anomalies. More intermittent behavior of the transverse velocity components with respect to the longitudinal one is found both by ESS and by comparison of the longitudinal and transverse velocity difference probability distribution functions (PDF). Received: 12 October 1997 / Accepted: 29 April 1998     

可压缩均匀剪切湍流的二阶矩模拟

基于对可压缩湍流中脉动压力场和脉动速度场特征的理论分析以及ＤＮＳ结果,建立了可均匀剪切湍流中压力－变形率关联的压缩性修正模式,应用这个模式,加上Ｓａｒｋａｒ等建立的脉动体胀率项（ｄｉｌａｔａｔｉｏｎａｌ ｔｅｒｍｓ）的模式,预测可压缩均匀剪切湍流随时间的发展,所得雷诺应力各是性张量的平衡值与Ｂｌａｉｓｄｅｌｌ等的ＤＮＳ数据非常一致。这个模式准确地预测出均匀剪切湍流中压缩性导致的雷诺应力结构的“流向     

Direct simulation of compressible turbulence in a shear flow

The purpose of this study is to investigate compressibility effects on the turbulence in homogeneous shear flow. We find that the growth of the turbulent kinetic energy decreases with increasing Mach number—a phenomenon which is similar to the reduction of turbulent velocity intensities observed in experiments on supersonic free shear layers. An examination of the turbulent energy budget shows that both the compressible dissipation and the pressure-dilatation contribute to the decrease in the growth of kinetic energy. The pressure-dilatation is predominantly negative in homogeneous shear flow, in contrast to its predominantly positive behavior in isotropic turbulence. The different signs of the pressure-dilatation are explained by theoretical consideration of the equations for the pressure variance and density variance. We previously obtained the following results for isotropic turbulence: first, the normalized compressible dissipation is of O(M _t ² ), and, second, there is approximate equipartition between the kinetic and potential energies associated with the fluctuating compressible mode. Both these results have now been substantiated in the case of homogeneous shear. The dilatation field is significantly more skewed and intermittent than the vorticity field. Strong compressions seem to be more likely than strong expansions.Dedicated to Professor J.L. Lumley on the occasion of his 60th birthday.This research was supported by the National Aeronautics and Space Administration under NASA Contract No. NAS1-18605 while the authors were in residence at the Institute for Computer Applications in Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665, U.S.A.     

Analyzing the influence of compressibility on the rapid pressure–strain rate correlation in turbulent shear flows

The influence of compressibility on the rapid pressure–strain rate tensor is investigated using the Green’s function for the wave equation governing pressure fluctuations in compressible homogeneous shear flow. The solution for the Green’s function is obtained as a combination of parabolic cylinder functions; it is oscillatory with monotonically increasing frequency and decreasing amplitude at large times, and anisotropic in wave-vector space. The Green’s function depends explicitly on the turbulent Mach number M _t , given by the root mean square turbulent velocity fluctuations divided by the speed of sound, and the gradient Mach number M _g , which is the mean shear rate times the transverse integral scale of the turbulence divided by the speed of sound. Assuming a form for the temporal decorrelation of velocity fluctuations brought about by the turbulence, the rapid pressure–strain rate tensor is expressed exactly in terms of the energy (or Reynolds stress) spectrum tensor and the time integral of the Green’s function times a decaying exponential. A model for the energy spectrum tensor linear in Reynolds stress anisotropies and in mean shear is assumed for closure. The expression for the rapid pressure–strain correlation is evaluated using parameters applicable to a mixing layer and a boundary layer. It is found that for the same range of M _t there is a large reduction of the pressure–strain correlation in the mixing layer but not in the boundary layer. Implications for compressible turbulence modeling are also explored.      

Hierarchical structures in a turbulent pipe flow

A hierarchical structure (HS) analysis (β-test and γ-test) is applied to a fully developed turbulent pipe flow. Velocity signals are measured at two cross sections in the pipe and at a series of radial locations from the pipe wall. Particular attention is paid to the variation of turbulent statistics at wall units 10<y⁺<3000. It is shown that at all locations the velocity fluctuations satisfy the She–Leveque hierarchical symmetry (Phys. Rev. Lett. 72 (1994) 336). The measured HS parameters, β and γ, are interpreted in terms of the variation of fluid structures. Intense anisotropic fluid structures generated near the wall appear to be more singular than the most intermittent structures in isotropic turbulence and appear to be more outstanding compared to the background fluctuations; this yields a more intermittent velocity signal with smaller γ and β. As turbulence migrates into the logarithmic region, small-scale motions are generated by an energy cascade and large-scale organized structures emerge which are also less singular than the most intermittent structures of isotropic turbulence. At the center, turbulence is nearly isotropic, and β and γ are close to the 1994 She–Leveque predictions. A transition is observed from the logarithmic region to the center in which γ drops and the large-scale organized structures break down. We speculate that it is due to the growing eddy viscosity effects of widely spread turbulent fluctuations in a similar way as in the breakdown of the Taylor vortices in a turbulent Couette–Taylor flow at high Reynolds numbers.     

Direct numerical simulation of compressible turbulent flows

This paper reviews the authors＇ recent studies on compressible turbulence by using direct numerical simulation （DNS）,including DNS of isotropic（decaying） turbulence, turbulent mixing-layer,turbulent boundary-layer and shock/boundary-layer interaction.Turbulence statistics, compressibility effects,turbulent kinetic energy budget and coherent structures are studied based on the DNS data.The mechanism of sound source in turbulent flows is also analyzed. It shows that DNS is a powerful tool for the mechanistic study of compressible turbulence.     

Numerical simulation of expansion/compression effect on shock induced turbulent flow

The interaction of homogeneous and isotropic turbulence with a shock wave is observed by solving the Reynolds-averaged Navier–Stokes equations with the k–? turbulence model. All turbulent fluctuations are measured at the period of expansion in the turbulent field and during compression by the reflected shock on turbulent field, and it is observed that the longitudinal turbulent velocity fluctuation is enhanced more at the period of expansion due to incident shock wave movement far from the turbulent field. The amplification of the turbulent kinetic energy (TKE) level in the shock/turbulence interaction depends on the shock wave strength and the longitudinal velocity difference across the shock wave. On decreasing the longitudinal velocity difference across the shock, the turbulent kinetic energy (TKE) level is less amplified. The TKE level is amplified by the factor of 1.5–1.8 in the shock/turbulence interaction where the dissipation rate of TKE decreases in all cases of shock/turbulence interaction. After the shock/turbulence interaction, the turbulent dissipative-length scale is amplified slightly and the amplification of the length scales decreases when increasing the shock strength. To cite this article: M.A. Jinnah, K. Takayama, C. R. Mecanique 333 (2005).     

Direct Numerical Simulation of a Mach 2 shock interacting with isotropic turbulence

Direct Numerical Simulation (DNS) and linear analysis of a shock interacting with incompressible and compressible isotropic turbulence is conducted. A dependence of amplification ratios on the degree of compressibility of the incoming flow is found. It can be shown that the enhancement of rms values of turbulent quantities across the shock varies according to the ratio of compressible to incompressible kinetic energy (exact definition see eq. 8). Inflow conditions with high values of display reduced amplification ratios of TKE and thermodynamic quantities while vorticity fluctuations are enhanced more strongly. The different behaviour of the turbulent kinetic energy (TKE) is due to the reduced pressure diffusion term in the TKE-equation. Experiments show qualitatively a similar behaviour as the simulation with incompressible inflow conditions, but they could so far not confirm our findings of reduced amplification rates in the compressible case, one of the reasons being the lack of knowledge of all flow parameters upstream of the shock front and the inability to generate isotropic turbulence in real life experiments. For the DNS we use a third order in space shock-capturing scheme based on the ENO algorithm of Harten [10] together with an approximate Riemann solver. This non-TVD scheme turned out to have many advantages over other common Godunov-type high resolution schemes for the specific problem of a shock interacting with turbulent fields.     

Flow-thermodynamics interactions in decaying anisotropic compressible turbulence with imposed temperature fluctuations

The fundamental nature of the non-linear flow-thermodynamics interactions in a compressible turbulent flow with imposed temperature fluctuations is investigated. Direct numerical simulations (DNS) of decaying anisotropic compressible turbulence (turbulent Mach number 0.06–0.6) with imposed temperature fluctuations are performed to examine: (i) interactions between solenoidal and dilatational kinetic energy; (ii) partition between dilatational kinetic energy and thermodynamic potential energy; and (iii) redistribution of solenoidal and dilatational kinetic energy among the various Reynolds stress components. It is found that solenoidal kinetic energy levels and return-to-isotropy are weakly dependent on Mach number but independent of imposed temperature fluctuations in the parameter range studied. The dilatational kinetic energy generated is proportional to the square of the pressure fluctuations associated with the initial solenoidal and temperature fluctuations and thus a strong function of Mach number and heat release intensity. The energy exchange between dilatational kinetic and potential energy is driven by a strong proclivity toward equipartition. Consequently, the dynamics of pressure-dilatation ( ${\overline{pd}}$ ), which is the mechanism of this energy exchange between dilatational and potential energies, is dictated entirely by the requirement to impose energy equipartition. Based on the results, we provide a physical picture of the solenoidal–dilatational–potential energy interactions and the action of pressure-dilatation. The identification of the fundamental precepts underlying the various interactions is of great utility for turbulence closure model development.     

Optimized group velocity control scheme and DNS of decaying compressible turbulence of relative high turbulent Mach number

To overcome the difficulty in the DNS of compressible turbulence at high turbulent Mach number, a new difference scheme called GVC8 is developed. We have succeeded in the direct numerical simulation of decaying compressible turbulence up to turbulent Mach number 0.95. The statistical quantities thus obtained at lower turbulent Mach number agree well with those from previous authors with the same initial conditions, but they are limited to simulate at lower turbulent Mach numbers due to the so‐called start‐up problem. The energy spectrum and coherent structure of compressible turbulent flow are analysed. The scaling law of compressible turbulence is studied. The computed results indicate that the extended self‐similarity holds in decaying compressible turbulence despite the occurrence of shocklets, and compressibility has little effects on relative scaling exponents when turbulent Mach number is not very high. Copyright © 2005 John Wiley & Sons, Ltd.