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 导致.     

二维槽道湍流数值模拟与白组织特性

从流体力学基本方程出发,讨论了二维槽道湍流的衰减特性,通过对流场施加合适的体积力,采用拟谱方法对二维槽道强制湍流进行了数值模拟.研究了二维槽道衰减湍流的自组织与逆能量级串特性,再现了二维槽道衰减湍流中湍涡的自组织过程,以及不同波数湍流结构所携能量在自组织过程中的变化,并解释了二维槽道湍流平均速度曲线特征以及海洋环流所特有的自然现象.     

关于植被中湍流的研究

本文从实验和模式理论两个方面介绍了近年来在植被湍流研究方面的进展.实验研究表明,植被内湍流是高度间歇性的,大尺度涡在湍流输运中起主导作用.植物枝、干、叶打碎了大尺度涡,产生的迹湍流中的小尺度涡更容易耗散成热。因此,植被内能量不按一般的能量级串理论预示的方式进行,而要发生所谓“短路”现象,使得功率谱曲线的斜率在惯性区比-2/3律更负.现有的高阶封闭模式虽然取得了很大成功,但在模拟湍流强度上高估了其大小,还需根据实验加以改进.      

二维槽道湍流数值模拟与自组织特性

从流体力学基本方程出发,讨论了二维槽道湍流的衰减特性,通过对流场施加合适的体积力,采用拟谱方法对二维槽道强制湍流进行了数值模拟。研究了二维槽道衰减湍流的自组织与逆能量级串特性,再现了二维槽道衰减湍流中湍涡的自组织过程,以及不同波数湍流结构所携能量在自组织过程中的变化,并解释了二维槽道湍流平均速度曲线特征以及海洋环流所特有的自然现象。     

FST方法分析二维Rayleigh-Bénard湍流热对流系统的能量输运

本文应用空间滤波方法：FST(Filter-space technique)方法,研究二维Rayleigh-Bénard(RB)湍流热对流系统中湍动能、热能和拟涡能的能量输运．研究中Rayleigh数(Ra)选取为1x10^8、1x10^9和1x10^10,Prandtl数(Pr)固定为4.38．我们展示了的结果表明,在二维RB系统中,三个Ra数下全场的平均湍动能和平均拟涡能在不同滤波尺度下的能量输运与Kraichnan在1967年预测的二维湍流中的级串理论有所偏差,而中心区域的能量都是向小尺度输运的．结果还揭示了瞬时能量输运的一些局部特性,包括它们在小尺度上不对称的分布．     

提取壁湍流相干结构的数字滤波法

以三丝热线探头测得平板湍流边界层的数据为对象，提出用数字滤波技术将湍流信号分解为接近各向同性的小尺度涡和非各向同性的大尺度涡的方法。并用条件采样技术从大涡信号中提取相干结构。     

壁面加热湍流相干结构的子波分析实验研究

通过在平板湍流边界层沿流向固壁表面平行放置若干条通电加热的金属细丝,在平板表面形成沿展向周期性分布的温度场,利用该温度场引起的空气热对流,在湍流边界层近壁区域产生一组沿湍流边界层展向周期分布的大尺度流向涡结构,改变了平板湍流边界层中不同尺度结构及其能量分布。采用对壁湍流多尺度结构的子波分析表明,在湍流边界层近壁区域产生规则的流向涡结构将壁湍流各种尺度湍涡结构不规则的脉动有序地组织起来,抑制了壁湍流各种尺度湍涡结构脉动,特别抑制了能量最大尺度湍涡结构的脉动,减小由于湍流脉动引起的在湍流边界层法向和展向的动量和能量损耗,从而减小了湍流的阻力。     

湍流边界层外区相干结构的三维波模型

根据流动稳定性理论,提出了一种三维波模型来描述湍流边界层外区大尺度相干结构。计算所得流线图和等涡量线图较罗纪生,周恒（1993）的二维波模型更符合实验结果。说明该三维模型能够较好地反映湍流边界层外区大尺度相干结构的物理特征。     

湍流基础问题研究进展:能量传递,相互作用尺度,各向同性衰减的自保持性

为了更深入地了解湍流的物理过程,本文综述了各向同性湍流的基础问题．在评述了Ｋｏｌｍｏｇｏｒｏｖ能谱及能量级串过程后,深入讨论了Ｋｏｌｍｏｇｏｒｏｖ局部各向同性假设．接着综述了涉及能量传递的以及包括三元组相互作用的各向同性湍流相互作用尺度的详细物理过程．还讨论了惯性区、自相似性以及小尺度对大尺度各向异性的响应和末期衰减过程．之后为了举例说明这些论点,详细讨论了根据各向同性湍流直接模拟及大涡模拟得到的结果（包括对亚格子模型的讨论）.最后,综述了各向同性湍流的自保持性,并展望了今后的研究方向．文末列出了１５５篇参考文献     

展向喷嘴湍流横向射流大涡模拟及其统计分析

采用大涡模拟方法数值模拟了展向椭圆喷嘴的湍流横向射流,对其大尺度结构的时空演化和湍流脉动速度场的时间序列分析、频谱分析、PDF分析以及时、空截面上的统计平均特性进行分析.结果表明,在射流出口附近的下游核心区中速度脉动剧烈,显现出明显的湍流特征.除了三维涡环脱落、扭曲、变形、摆动所对应频率之外,还存在很宽的湍流基频,它与在喷嘴出口附近产生的三维涡环的时空演化过程密切相关.由于展向椭圆喷嘴的湍流横向射流中的三维涡环快速脱落和强相互作用导致射流尾迹中的强湍流脉动,展向椭圆喷嘴湍流横向射流的PDF空间演化特征结构复杂.在射流核心区的湍流偏应力变化平缓,其统计平均值分布接近左右对称.展向椭圆喷嘴的湍流横向射流脉动速度场具有极为复杂的统计行为,与流向椭圆喷嘴相比具有更好的掺混能力.     

Direct simulation of two‐dimensional turbulent flow over a surface‐mounted obstacle

Two‐dimensional turbulent flow over a surface‐mounted obstacle is studied as a numerical experiment that takes place in a wind tunnel. The transient Navier–Stokes equations are solved directly with Galerkin finite elements. The Reynolds number defined with respect to the height of the wind tunnel is 12 518. Instantaneous streamline patterns are shown that give a complete picture of the flow phenomena. Energy and enstrophy spectra yield the dual cascade of two‐dimensional turbulence and the ?1 power law decay of enstrophy. Mean values of velocities and root mean square fluctuations are compared with the available experimental results. Other statistical characteristics of turbulence such as Eulerian autocorrelation coefficients, longitudinal and lateral coefficients are also computed. Finally, oscillation diagrams of computed velocity fluctuations yield the chaotic behaviour of turbulence. Copyright © 2007 John Wiley & Sons, Ltd.     

Systematic tuning of dispersion models for simulation of evaporating sprays

In this study, via an Eulerian–Lagrangian framework, the performance of two recent dispersion models, i.e. a first-order autoregressive process and the PDF model, is compared. The appropriate relations for the turbulence scales and the drift correction term are suggested and the tuned values for the constants of the models are proposed in a systematic approach by starting with the simplest case, i.e. particle-laden stationary isotropic turbulence and adding more complexities in the subsequent cases, including the homogeneous anisotropic shear flow, decaying grid turbulence, and inhomogeneous gas–solid spray. Also, the isotropic relation for the effect of inertia in the Lagrangian turbulence time scale seen by particles is extended to the anisotropic case while it remains consistent in the isotropic limit. Finally, the performance of the tuned models is evaluated for the simulation of an evaporating spray. It is observed that, the tuned constants for the evaporating spray are close to the ones obtained for the homogeneous shear flow.     

Un opérateur de diffusion spatialement sélectif pour le transport d'un traceur passif ou actif par un écoulement de grande échelle

In a fluid flow, fields are measurable up to a cut-off scale at which they are regularized. We show that, for a smooth velocity field, this regularization adds to the advection equation a diffusive term proportional to the strain tensor. We study in two dimensions its effect on the dynamics of velocity and vorticity, and on the conservation of quadratic invariants. Vorticity and energy are still conserved, while enstrophy and tracer variance are dissipated depending on the flow topology. These properties (conservation, dissipation, spatial selectivity) suggest the use of this selective strain–diffusion operator for numerical simulations of inhomogeneous flows in the quasi-two-dimensional approximation.     

Lagrangian Covariance Analysis of β-Plane Turbulence

The effects of Rossby wave–turbulence interactions on particle dispersion are investigated in a Lagrangian analysis of the potential vorticity equation. The analysis produces several exact statistical results for fluid particle dispersion in barotropic turbulence on a β-plane. In the inviscid problem the first integral time scale of the meridional velocity is found to be zero, as might occur in pure wave processes, and the meridional particle dispersion is bounded. The second integral time scale, which determines the magnitude of the bound, is shown to depend explicitly on β, the enstrophy and the energy of the meridional velocity. Expressions relating the autocorrelation of the vorticity to the autocorrelation of the meridional velocity are derived and the Lagrangian integral time scale of the relative vorticity is diagnostically related to the meridional velocity correlation. The applicability of these predictions is verified in a series of numerical simulations. For a range of β values, the meridional extent of quasisteady alternating zonally averaged jets occurring in the numerical solutions scales with a length scale given by the the standard deviation of the meridional particle dispersion. Received 8 March 1999 and accepted 12 December 1999     

On the nature of multitude scalings in decaying isotropic turbulence

We report multitude scaling laws for isotropic fully developed decaying turbulence through group theoretic method employing on the spectral equations both for modelling and without any modelling of nonlinear energy transfer. For modelling, besides the existence of classical power law scalings, an exponential decay of turbulent energy in time is obtained subject to exponentially decaying integral length scale at infinite Reynolds number limit. For the transfer without modelling, at finite Reynolds number, in addition to general power law decay of turbulence intensity with integral length scale growing as a square root of time, an exponential decay of energy in time is explored when integral length scale remains constant. Both the power and exponential decaying laws of energy agree to the theoretical results of George (1992), George and Wang (2009) and experimental results of fractal grid generated turbulence by Hurst and Vassilicos (2007). At infinite Reynolds number limit, a general power law scaling is obtained from which all classical scaling laws are recovered. Further, in this limit, turbulence exhibits a general exponential decaying law of energy with exponential decaying integral length scale depending on two scaling group parameters. The role of symmetry group parameters on turbulence dynamics is discussed in this study.     

粉尘等容燃烧容器内扬尘系统诱导湍流特性的实验研究

采用热线风速仪与系统平均法测定和研究了三种圆柱形封闭容器内的扬尘诱导湍流衰减特性,并与球形爆炸容器内已测定的扬尘诱导湍流衰减特性进行了比较。探讨了扬尘装置、容器体积、形状对扬尘诱导湍流瞬态特性的影响。实验结果显示,在粉尘等容燃烧容器内扬尘诱导湍流强度随时间的衰减特性具有一定相似性,它们均呈负指数关系。     

Couplage d'un modèle stochastique lagrangien sous-maille avec une simulation grandes échelles

A large eddy simulation (LES) is used to compute passive scalar dispersion in a turbulent flow. Instead of resolving the passive scalar transport equation, fluid particles are tracked in a Lagrangian way and a Langevin stochastic modelling is used for the small scale component of the velocity of fluid particles at a subgrid-scale level. The stochastic model is written in terms of subgrid-scale statistics. The specificity of this study resides in the coupling of a LES with a stochastic model using the filtered subgrid-scale statistics in inhomogeneous turbulence. The results are compared with the wind tunnel experiments of Fackrell and Robins [J. Fluid Mech. 117 (1982) 1–26]. To cite this article: I. Vinkovic et al., C. R. Mecanique 333 (2005).     

Effects of the Similarity Model in Finite-Difference LES of Isotropic Turbulence Using a Lagrangian Dynamic Mixed Model

A Lagrangian dynamic formulation of the mixed similarity subgrid (SGS) model for large-eddy simulation (LES) of turbulence is proposed. In this model, averaging is performed over fluid trajectories, which makes the model applicable to complex flows without directions of statistical homogeneity. An alternative version based on a Taylor series expansion (nonlinear mixed model) is also examined. The Lagrangian models are implemented in a finite difference code and tested in forced and decaying isotropic turbulence. As comparison, the dynamic Smagorinsky model and volume-averaged formulations of the mixed models are also tested. Good results are obtained, except in the case of low-resolution LES (32³) of decaying turbulence, where the similarity coefficient becomes negative due to the fact that the test-filter scale exceeds the integral scale of turbulence. At a higher resolution (64³), the dynamic similarity coefficient is positive and good agreement is found between predicted and measured kinetic energy evolution. Compared to the eddy viscosity term, the similarity or the nonlinear terms contribute significantly to both SGS dissipation of kinetic energy and SGS force. In order to dynamically test the accuracy of the modeling, the error incurred in satisfying the Germano identity is evaluated. It is found that the dynamic Smagorinksy model generates a very large error, only 3% lower than the worst-case scenario without model. Addition of the similarity or nonlinear term decreases the error by up to about 50%, confirming that it represents a more realistic parameterization than the Smagorinsky model alone.