湍流场中被动标量的细微结构

本文采用直接数值模拟方法,在具有平均标量梯度的各向同性湍流中,研究被动标量的小尺度结构特性及其与湍流场中应变与涡量的关系.对欧拉统计量及拉格朗日统计量的统计表明:标量耗散的形成主要是由于标量梯度同流场的应变张量压缩主轴耦合的结果,而涡量对标量梯度的形成只有较弱的影响,然而它可以间接影响大强度标量耗散的产生.强标量耗散的细微片状结构的形成时间尺度大约为10倍Kolmogrov时间尺度;在形成强标量梯度的细微片状结构过程中,应变强度随标量梯度同步增大,而涡量则先减小后增大,并在5倍Kolmogorov时间尺度时达到最大.     

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

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

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

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

旋转槽道湍流标量场输运特性数值模拟研究

展向旋转槽道湍流中的标量场输运过程是与众多工程流动直接相关的模型问题。基于直接数值模拟工作对该问题开展系统性的研究。由壁面摩擦速度定义的流动Reynolds数固定在180,重点考察Schmidt数和旋转数的影响。结果表明,较弱旋转强度即可对主导流动结构形态产生明显的影响：此时槽道不稳定侧产生流向大尺度结构,由此导致标量场出现条带状结构。强旋转时不稳定侧出现被湍流充分混合的区域,而在稳定侧流动层流化并出现近似传导区。平均标量剖面在湍流区和层流区呈现斜率不同的线性分布。Schmidt数小于1时,湍流区标量场脉动和湍流输运随旋转数出现非单调变化,而Schmidt数大于等于1时两者都随旋转数单调下降。由此导致总标量传输率在Schmidt数小于1时随旋转数先上升后下降,而当Schmidt数大于1时单调下降且在弱旋转区域下降趋势最快。     

雷诺数对圆柱尾流中被动标量场的影响

本文通过实验研究雷诺数对加热圆柱尾流中温度场的影响.实验中雷诺数Re(≡U∞d/v,其中U∞为来流速度、d为圆柱直径、v为流体黏度)的取值范围为1200-8600.实验中温度是由直径为0.63 μm的冷线探针测量的.实验结果表明,一般而言,雷诺数对整个尾流的标量混合特性有着显著的影响.随着雷诺数的增加,平均标量场向外的扩散速度加快、标量脉动强度增加了但衰减也加快.本文还发现:尾流中似乎存在两个区域,一个位于卡门涡街下游靠后,另一个就是传统的远场自相似区;在这两个区域,某些描述标量和动量的相似关系式近似成立.     

惯性颗粒所见被动标量统计特性的DNS与模型研究

本文采用DNS方法,对惯性颗粒所见各向同性湍流中具有平均标量梯度的被动标量场统计特性进行了研究。结果表明：惯性对颗粒温度脉动强度,两相温度关联,自相关特性以及颗粒热流与两相交叉热流的统计特性具有明显的影响。在PDF方法的框架下,系统地推导了非等温气固两相流的PDF方程,且基于朗之万随机体系对方程进行了封闭,并利用前面的...     

超音速平板混合层被动标量的数值模拟

采用空间五阶精度紧致-WENO混合格式,数值模拟了Mc=0.2和1.0的二维空间发展超音速平板混合层被动标量混合过程.并与Clemens的实验结果进行比较,结果表明,被动标量时均值与实验结果符合很好,而被动标量脉动均方根值与实验值相比较小.还研究了可压缩性对标量混合过程的影响,随着对流马赫数的增加,被动标量时均值剖面变化不大,而脉动均方根峰值增大.     

三维均匀各向同性两相湍流的直接模拟

本文对三维气粒两相均匀各向同性湍流进行了直接数值模拟。气相控制方程组采用分布投影方法进行求解,微分方程采用六紧致阶差分格式和快速Fourier变换结合求解;计算颗粒场时,采用Lagrangian方法。由该方法得到的能谱和各统计量与由谱方法得到的对应值进行了比较,吻合十分理想,对不同Stocks数颗粒在流场内的瞬态分布也进行了初步模拟,并观察到局部富集现象,证明该方法是进行两相湍流直接数值模拟行之有效的方法。     

圆柱尾流温度标量场的小波分析*

通过使用一排16根冷线探头排在多个空间点同时测量微加热圆柱的尾流温度场,用小波分析技术对瞬时温度场的时间序列信号进行多尺度分析,目的是研究不同尺度脉动温度对总体温度场的贡献.直径为d=12.7 mm的圆柱产生了被测的尾流,对应的雷诺数为5500,测量区域位于下游距离为2d和20d之间.基于小波多尺度分辨技术,尾流温度场被分解为不同温度脉动特征尺度的小波分量.通过分析这些小波分量的瞬时温度等值线图,能够直接观测到不同特征尺度的涡结构运动特征和湍流间歇过程.特别地,我们在近场区从原始信号分解获得的高频区域中发现了K-H涡的存在.不同尺度的温度方差沿流向的变化表明,在下游距离为x=3d和20d之间,中等尺度的结构比大尺度和小尺度结构对总的温度均方根的贡献更大.不同尺度的自相关函数表明,大尺度和中等尺度的结构显示出较大的相关性,而高频的小波分量则更快地失去了原有的拟序性.     

Scalar Statistics along Inertial Particle Trajectory in Isotropic Turbulence

The statistics of a passive scalar along inertial particle trajectory in homogeneous isotropic turbulence with a mean scalar gradient is investigated by using direct numerical simulation. We are interested in the influence of particle inertia on such statistics, which is crucial for further understanding and development of models in non-isothermal gas-particle flows. The results show that the scalar variance along particle trajectory decreases with the increasing particle inertia firstly; when the particle＇s Stokes number St is less than 1.0, it reaches the minimal value when St is around 1.0, then it increases if St increases further. However, the scalar dissipation rate along the particle trajectory shows completely contrasting behavior in comparison with the scalar variance. The mechanical-to-thermal time scale ratios averaged along particle, （r）p, are approximately two times smaller than that computed in the Eulerian frame r, and stay at nearly 1.77 with a weak dependence on particle inertia. In addition, the correlations between scalar dissipation and flow structure characteristics along particle trajectories, such as strain and vorticity, are also computed, and they reach their maximum and minimum, 0.31 and 0.25, respectively, when St is around 1.0.     

平板通道内充分发展湍流直接数值模拟

采用有限差分法对平板通道内充分发展湍流进行了直接模拟(DNS),Reτ=150,网格数64×64×64。得到了湍流脉动速度沿y 分布的等值线图,分析了近壁面脉动特性研究的重要性,为理解湍流物理机制提供数据。获得了雷诺应力的变化、雷诺应力输运方程中各项的平衡关系以及湍动能方程中各项的变化曲线,为进一步评价和改进湍流模式提供依据。     

Scalar dispersion by coherent structures in uniformly sheared flow generated in a water tunnel

In order to investigate the role of coherent structures as mechanisms of scalar dispersion, we studied measurements of a passive scalar plume released in a uniformly sheared turbulent flow generated in a water tunnel. The flow had homogeneous turbulence properties in the measurement domain and contained hairpin vortices similar to those in boundary layers, and so was an ideal test bed to study the effects of coherent structures on turbulent dispersion, free from the effects of inhomogeneities or boundaries. Measurements of the velocity and concentration fields were acquired simultaneously using stereo particle image velocimetry and planar laser-induced fluorescence. We found that dye was preferentially located far away from vortices and was less likely to appear in close proximity to vortices, which is attributed to the high dissipation at the periphery of the vortices. However, we also found that dye was not directly correlated with the uniform momentum zones in the flow, suggesting a more complex relationship exists between these zones, the locations of vortices, and dye transport. Considering scalar flux events rather than simply the presence of dye as our condition of interest, a conditional eddy analysis demonstrated that hairpin vortices are responsible for the large scalar flux events as well as the large Reynolds stress events in the flow. The fact that the Reynolds stress was correlated with the scalar flux further confirmed that coherent structures are dominant mechanisms for scalar transport. Furthermore, we found that the scalar flux vector was preferentially inclined by 155° and ?25° with respect to the streamwise direction, and was thus approximately orthogonal to the planes of the legs of the most common upright and inverted hairpin structures in the flow. These findings demonstrate that coherent structures play an important and intricate role in turbulent diffusion.     

自定义标量法模拟丙烷燃烧过程

采用自定义标量法模拟丙烷扩散燃烧,该方法通过把反应组分定义为Fluent程序的自定义标量、化学反应速率作为源项求解质量、动量、能量和组分守恒方程,并用化学反应引起的能量变化修正能量方程.考虑了详细的化学反应机理,整个燃烧反应机理包括27种化学物质(不含N2)和83个基元反应.合理地模拟出了丙烷的燃烧过程,并将火焰的长度、温度、丙烷、氧气以及中间产物的分布与实验数据进行比较.     

Decaying turbulence in the presence of a shearless uniform kinetic energy gradient

The decay of turbulent kinetic energy in nearly isotropic grid turbulence has been studied extensively as a fundamental point of reference for turbulence theories and numerical simulations. Most studies have focused on nearly homogeneous turbulence characterised by power-law decay. Other studies have focused on so-called shearless mixing layers, in which two regions with the same mean velocity but distinctly different kinetic energy levels slowly diffuse into each other downstream thus providing information about spatial transport of turbulence. Here, we introduce and study another type of shearless turbulent flow. It has initially a nearly uniform spatial gradient of kinetic energy of the form k ～ β(y ? y₀), where y is the spanwise position. In the experiments, this gradient is generated with the use of an active grid and screens mounted upstream of the wind-tunnel’s test section, iteratively designed to produce a uniform gradient of turbulent kinetic energy without mean velocity shear. Data are acquired using X-wire thermal anemometry at different spanwise and downstream locations. Profile measurements are used to quantify the constancy of the mean velocity and the linearity of the initial profile of kinetic energy. Measurements show that at all spanwise locations, the decay in the streamwise direction follows a power-law but with exponents n(y) that depend upon the spanwise location. The results are consistent with a decay of the form k/?u?² = β(x/x_ref)^?n(y)(y ? y₀)/M. Results for the development of integral length scale, and for velocity skewness and flatness factors are also presented. Significant deviations from Gaussianity are observed especially for the spanwise velocity component in the lower kinetic energy region. Future experiments will be needed including measurements of the dissipation rate ? at sufficient accuracy, in order to unambiguously partition the energy decay into dissipation and spatial diffusion.     

Reassessment of the classical closures for scalar turbulence

In deducing the consequences of the Direct Interaction Approximation, Kraichnan was sometimes led to consider the properties of special classes of nonlinear interactions in degenerate triads in which one wavevector is very small. Such interactions can be described by simplified models closely related to elementary closures for homogeneous isotropic turbulence such as the Heisenberg and Leith models. These connections can be exploited to derive considerably improved versions of the Heisenberg and Leith models that are only slightly more complicated analytically. This paper applies this approach to derive some new simplified closure models for passive scalar advection and investigates the consistency of these models with fundamental properties of scalar turbulence. Whereas some properties, such as the existence of the Kolmogorov–Obukhov range and the existence of thermal equilibrium ensembles, follow the velocity case closely, phenomena special to the scalar case arise when the diffusive and viscous effects become important at different scales of motion. These include the Batchelor and Batchelor–Howells–Townsend ranges pertaining, respectively, to high and low molecular Schmidt number. We also consider the spectrum in the diffusive range that follows the Batchelor range. We conclude that improved elementary models can be made consistent with many nontrivial properties of scalar turbulence, but that such models have unavoidable limitations.     

A new assessment of the second-order moment of Lagrangian velocity increments in turbulence

The behaviour of the second-order Lagrangian structure functions on state-of-the-art numerical data both in two and three dimensions is studied. On the basis of a phenomenological connection between Eulerian space-fluctuations and the Lagrangian time-fluctuations, it is possible to rephrase the Kolmogorov 4/5-law into a relation predicting the linear (in time) scaling for the second-order Lagrangian structure function. When such a function is directly observed on current experimental or numerical data, it does not clearly display a scaling regime. A parameterisation of the Lagrangian structure functions based on Batchelor model is introduced and tested on data for 3d turbulence, and for 2d turbulence in the inverse cascade regime. Such parameterisation supports the idea, previously suggested, that both Eulerian and Lagrangian data are consistent with a linear scaling plus finite-Reynolds number effects affecting the small- and large timescales. When large-time saturation effects are properly accounted for, compensated plots show a detectable plateau already at the available Reynolds number. Furthermore, this parameterisation allows us to make quantitative predictions on the Reynolds number value for which Lagrangian structure functions are expected to display a scaling region. Finally, we show that this is also sufficient to predict the anomalous dependency of the normalised root mean squared acceleration as a function of the Reynolds number, without fitting parameters.