标量脉动“能量”传递函数与—5/3次方律

本文在Heisenberg关于速度脉动能量传递函数假设的基础上,提出标量脉动“能量”传递函数的相应形式.讨论了具有常值横向平均速度梯度与标量梯度的均匀湍流流场中关于大波数情形的标量脉动“能”谱性质。     

惯性颗粒所见被动标量统计特性的直接数值模拟研究

本文采用直接数值模拟方法,对惯性颗粒所见均匀各向同性湍流中具有平均标量梯度的被动标量场统计特性进行了研究。模拟结果表明:与颗粒所见流体速度的自相关特性不同,颗粒所见标量的自相关特性随颗粒惯性的增加而单调减少;颗粒所见标量脉动能随颗粒惯性的增加先减少再增大,在St≈1.0的临界颗粒附近达到最小值,而颗粒所见标量耗散率随颗粒惯性的变化行为与颗粒所见标量脉动能的变化行为相反;数值模拟的结果进一步揭示,在St≈1.0时颗粒所见流体标量脉动能和耗散率的极值是因为St≈1.0的临界颗粒聚集于低涡量、高应变区域和标量场在高应变区域形成强耗散的片状结构所致。     

颗粒惯性对标量混合的影响

本文采用直接数值模拟的方法,研究了颗粒惯性对均匀各向同性湍流场中被动标量混合特性的影响.结果表明:与无颗粒算例相比,Tp/Tk<1的小颗粒,增大了标量脉动在高波数上的能量,增强了标量耗散率,标量耗散率最大高出了25%;Tp/Tk>1的大颗粒,在初期显著削弱了标量脉动强度,同时也加快了标量耗散率的衰减,标量脉动强度降低了约15%,标量耗散减少了约5%;而Tp/Tk=1的临界颗粒,对标量场混合的影响较弱,标量脉动强度和标量耗散偏差在5%以内;此外,随着颗粒惯性的增大,标量耗散率的概率密度函数向左平移,但分布特征的改变并不显著.     

氢气扩散火焰中辐射源项湍流脉动特征的PDF模拟

采用κ-ε湍流模型、标量联合的概率密度函数(PDF)输运方程和层流火焰面模型相结合,模拟氢气自由扩散火焰中辐射源项湍流脉动特征.给出了主燃区内辐射源项湍流脉动的频率图.辐射源项的样本点分布集中,大约95%以上的样本落在其系综的±3倍方差以内,频谱图为单峰.     

湍流热对流DNS多分辨率并行直接求解方法

依据温度标量场与动量计算的空间和时间计算分辨率不同的特点,采用两套网格,建立多分辨率双网格并行直接求解方法,用以解决极高Ra数湍流热对流DNS模拟巨大计算工作量的难题.在两套网格的数据交换上,根据每个细网格都满足连续方程,设计了速度的守恒平移插值方法.二维极高Ra数湍流热对流的计算结果表明,采用多分辨率双网格并行直接求解方法的DNS计算,可以使计算工作量降低近一个量级.瞬时温度场显示,双网格方法的计算结果可以很好地描述极高Ra数下快速运动的小尺寸漩涡团状羽流,得到的结果与原网格一致,不同方法计算得到的传热Nu数误差不超过1%.     

沟槽壁面减阻机理实验研究

利用IFA300型热线风速仪,测量了光滑壁面和沟槽减阻壁面湍流边界层内的瞬时速度,利用自行设计的阻力天平仪测量了壁面摩擦力。得到了边界层无量纲速度分布和平均湍动能分布。对测得的脉动速度信号,利用离散正交小波变换按时间和尺度分解,得到各尺度分量的湍动能,并且发现其分布在湍流惯性区具有极大值。分析表明,当沟槽有减阻效果时,边界层内的平均湍动能减小,湍流惯性区各分量的湍动能极大值亦减小。     

QUICK格式在湍流旋流流动数值模拟中的应用

为研究不同精度的离散格式对湍流旋流流动数值模拟结果的影响,同时应用QUICK格式和混合格式对同轴射流旋流燃烧室内的湍流流动进行了数值模拟.在采用k-ε湍流模型的条件下,QUICK格式计算得到的燃烧室内气体轴向与切向速度及轴向脉动速度均方根值分布与实验数据符合较好,而混合格式给出的数值模拟结果则与实验有一定的偏差.     

尾流的双尺度简化模型及其大涡模拟建模

叶轮机中存在着复杂的湍流流动结构,如何准确地模拟湍流流动一直是工程界关注的核心问题。由于湍流在大部分区域都处于强非均衡态,根据均衡能量传输假设建立起来的RANS湍流模型和大涡模拟亚格子模型都不再适用。非均衡湍流的一个经典的模型流动是尾流,但通常意义上的尾流往往是无数不同尺度相互作用的结果,难以进行简化分析。本文构造了一个极度简化的双尺度尾流模型进行了直接数值模拟(DNS),结果显示在转捩前期就存在反向的非均衡能量传输。分别采用各种传统的RANS湍流模型和大涡模拟亚格子模型进行了数值模拟,发现各种模型都无法准确预测尾流转捩。在考虑非均衡能量传输机制后,通过流向速度梯度扭率构造的理性亚格子应力模型可以使大涡模拟模拟的转捩位置相对于标准Smagorinsky模型有所改善;不采用涡粘假设的改进速度增量(IVI)模型可以得到更加接近DNS的结果。这些结果为湍流模型进一步的改进提供了思路。     

大气湍流中光束的高阶强度矩

研究了光束通过大气湍流传输的高阶强度矩, 提出了大气湍流中光束高阶强度矩的推导方法, 并推导出了一至四阶光束强度矩传输的解析表达式. 所得结果具有一般性,任意某一光束在自由空间和大气湍流中传输的高阶强度矩均可作为本文结果的特例. 另一方面, 以高斯光束为例, 研究了其K参数在湍流大气中的传输规律. 研究表明,高斯光束在大气湍流中其K参数并不是一个传输不变量,它与传输距离、束腰半径、湍流内外尺度以及湍流强度均有关.这个结论与采用Rytov相位结构函数二次近似或强湍流近似下的结论不同,本文给出了合理解释. 关键词： 高阶强度矩 大气传输 大气湍流 K参数')" href="#">K参数     

声波扰动对大气湍流内外尺度与折射率功率谱函数的影响分析

声波是一种机械波,作为能量的载体,在大气湍流环境下传输会“扰动”湍流耗散率的变化,从而会影响湍流物理结构演化.本文基于声波能量和湍流能量平衡方程,结合湍流内外尺度和大气折射率功率谱函数,研究了在不同声波扰动下大气湍流的内外尺度和折射率功率谱函数的变化特征.结果表明:不同声波的传播会使得湍流的内外尺度发生变化,声源功率越大,对湍流尺度的影响越大,然而声源频率越大,对湍流尺度的影响并不是特别明显;不同声波的传播会使大气折射率功率谱函数发生改变,在惯性区内,考虑到声波对湍流内外尺度的影响,不同声源对大气折射率功率谱的影响程度不同,在耗散区内,大气折射率功率谱都出现随声波传输距离波动的情况.本文探索声波扰动对大气湍流折射率功率谱函数特征参数的变化规律,为激光在声波扰动大气湍流中传输特性以及声光耦合研究提供理论依据.     

Generalised scale-by-scale energy-budget equations and large-eddy simulations of anisotropic scalar turbulence at various Schmidt numbers

A necessary condition for the accurate prediction of turbulent flows using large-eddy simulation (LES) is the correct representation of energy transfer between the different scales of turbulence in the LES. For scalar turbulence, transfer of energy between turbulent length scales is described by a transport equation for the second moment of the scalar increment. For homogeneous isotropic turbulence, the underlying equation is the well-known Yaglom equation. In the present work, we study the turbulent mixing of a passive scalar with an imposed mean gradient by homogeneous isotropic turbulence. Both direct numerical simulations (DNS) and LES are performed for this configuration at various Schmidt numbers, ranging from 0.11 to 5.56. As the assumptions made in the derivation of the Yaglom equation are violated for the case considered here, a generalised Yaglom equation accounting for anisotropic effects, induced by the mean gradient, is derived in this work. This equation can be interpreted as a scale-by-scale energy-budget equation, as it relates at a certain scale r terms representing the production, turbulent transport, diffusive transport and dissipation of scalar energy. The equation is evaluated for the conducted DNS, followed by a discussion of physical effects present at different scales for various Schmidt numbers. For an analysis of the energy transfer in LES, a generalised Yaglom equation for the second moment of the filtered scalar increment is derived. In this equation, new terms appear due to the interaction between resolved and unresolved scales. In an a-priori test, this filtered energy-budget equation is evaluated by means of explicitly filtered DNS data. In addition, LES calculations of the same configuration are performed, and the energy budget as well as the different terms are thereby analysed in an a-posteriori test. It is shown that LES using an eddy viscosity model is able to fulfil the generalised filtered Yaglom equation for the present configuration. Further, the dependence of the terms appearing in the filtered energy-budget equation on varying Schmidt numbers is discussed.     

Conserved scalar mixing in a confined-opposed-jets flow

The focus of this paper is on the mixing of a conserved passive scalar for Sc = 1 (Sc is the Schmidt number) in axisymmetric turbulence for which the initial injections of turbulent kinetic energy and scalar variance are similar. Two confined-opposed-jets (COJ) are experimentally studied through simultaneous PIV (particle image velocimetry) and PLIF (planar laser induced fluorescence) measurements, for different flow regimes. One-point transport equation for the scalar variance is assessed through experimental data, along the common axis of the two opposed jets, and different physical phenomena are revealed (production, diffusion, dissipation). The production of scalar variance is equilibrated by the diffusion term (～75%) and the mean dissipation of the scalar variance (～25%). To further assess the scalar behaviour at each scale in this anisotropic, but axisymmetric, flow, a scale-by-scale scalar variance budget equation is derived for axisymmetric turbulence. This equation reduces to Yaglom's 4/3 law, under additional restrictions. The equation is assessed through experimental data, in the impingement region between the two COJ. In particular, the anisotropic energy transfer along different directions is quantified. It is shown that for scales smaller than the size of the central region, Δ, the cascade of the scalar variance is completely inhibited, independently of the particular direction. For scales larger than Δ, the apparent aspect of the energy transfer is that of an inverse cascade, with positive values of the scalar variance transfer. Nonetheless, inhomogeneity of the flow and mixing at those scales is directly responsible for these positive values.     

The energy and dissipation of turbulent fluctuations of wind velocity and temperature in the boundary layer of the atmosphere

The relationships between the energy of small-scale turbulence and its dissipation rate are studied based on the data of long-term high-frequency measurements of temperature and wind velocity fluctuations in urban area. It is shown that the energy of wind velocity turbulent fluctuations is linearly related to the dissipation rate ɛ. The proportionality coefficient between turbulent kinetic energy (TKE) and ɛ is dimensional and does not depend on the stratification of the atmosphere, the Richardson number, or the Monin-Obukhov scale. Measurements in different seasons show that this coefficient can be related to the mean velocity of adiabatic motions (sound speed or air temperature), which enables one to select a more universal constant, γ. A linear relationship between the temperature fluctuations variance (the characteristic of the inner energy of turbulence) and their dissipation rate is also shown. The revealed proportionality is confirmed by measurements in urban and forest conditions, as well as in the surface layer over a flat desert terrain.     

Etude numérique des jets turbulents à température élevée

Numerical prediction of the structure of high temperature axisymmetric turbulent jets. Turbulent axisymmetric jets at high temperature are studied numerically by using first and second order turbulence models. Regarding the temperature fields, on which we concentrate in this work, predictions with both types of models do not show large differences. In general, predictions agree well with the measurements; the existing differences are usually favorable for the second order model. The effect of solving a transport equation for the scalar dissipation rate on the prediction of the mechanical to scalar time scale ratio and on the prediction of the scalar fluctuations is studied. The influence of varying the density ratio on parameters such as the axial decay rates of the temperature and velocity and the turbulence intensity are studied. Two definitions of the mixing efficiency are introduced. According to both definitions, the mixing efficiency decreases with increasing effects of buoyancy.     

Direct numerical simulation of passive scalar in decaying compressible turbulence

1IntroductionDirectnumericalsimulation(DNS)becomesanimportanttoolinrecentresearchofturbulence[1].DNSofcompressibleturbulenceismoredifficultthanthatoftheincompressibleturbulence.WhentheturbulentMachnumberisgreaterthan0.3theshockletsmayappearinthecompressibleturbulentflowfields.Thereasonandmechanismofshockletsexistencearenotclearyet.TheturbulentMachnumberinDNScannotbeveryhighwiththepresentexistingnumericalmethodsandcomputerresource.Fortheproblemofcompressibleisotropicturbulencewiththeinitia…     

Restricted scaling range models for turbulent velocity and scalar energy transfers in decaying turbulence

The effect of finite Reynolds numbers and/or internal intermittency on the total kinetic energy and scalar energy transfers is examined in detail. For this purpose, two distinct models for velocity and scalar energy transfer are proposed in the specific context of freely decaying isotropic turbulence. The first one extends the already existing dynamical models (hereafter DYM, i.e. based on transport equations originated in Navier–Stokes and advection-diffusion transport equations). The second one relies on the characteristic time of the strain at a specific scale (hereafter SBM). Both models account for the Reynolds number dependence of the scaling exponent of the second-order structure functions, over a range of scales where such exponents may be defined, i.e. a restricted scaling range (RSR). Therefore, the models developed aim at reproducing the energy transfer over the RSR. The predicted energy transfer is very sensible to variations of the scaling exponent, especially at low Reynolds numbers. The approach towards the asymptotic 4/3 law is closely reproduced by the two models. The dynamical model reproduces the experimental results accurately especially in the vicinity of the Taylor microscale, while the SBM agrees almost perfectly with measurements at nearly all scales.