壁湍流相干结构尺度及边界层内的SL标度律

利用高分辨率、高帧率PIV系统对湍流边界层中相干结构的多种空间尺度和边界层内SL 标度律在不同尺度下的具体表达形式进行了实验研究. 实验在两个动量损失厚度雷诺数 (Re_{\theta}=628.5和Re_{\theta}=1032.9)下测量平板湍流边界层中缓冲 层、对数区和外区的二维瞬时速度场. 应用 小波分析以及传统的统计学方法，在垂直于平板和平行于平板的平面内考察平板湍流边界层 中存在的相干结构的流向和展向尺度，并与已知的相干结构尺度实验结果进行了对比分析. 利用在动量损失厚度雷诺数628.5下测得的数据，对多种脉动结构(脉动速度结构等) 的空间关系及其标度律进行了研究. 第2项工作直接利用湍流边界层空间速度分布，对多种 流场尺度结构内部的She-Leveque(简称SL)标度律及自相似律进行了验证. 结果表明，各 单一流场尺度结构内部，流向脉动速度{\pmb u}'、法向脉动速度{\pmb v}'及 脉动涡分量\d {\pmb v}'/\d {\pmb x}的统计结构量均存在明显的标度律，标度 指数的形式与自相似律和SL标度律均非常吻合，只是常数随流场尺度的不同而不同， 且呈现一定的规律性. 但对于结构量的五阶矩随距离l的研究表明，自相似律和SL 标度律成立的范围并不完全一致，同时标度律成立的范围大小与流场尺度有明显关系.     

壁面湍流标度律的PIV实验研究

利用高分辨率、高帧率PIV系统对湍流边界层内的自相似标度律在不同的壁面距离的具体表达形式进行了实验研究。在动量损失雷诺数（Reθ=2167）下测量平板湍流边界层内的二维瞬时速度场。应用小波分析和传统统计学方法,在垂直于平板的平面内,对不同法向位置的自相似的标度律的具体形式及其随尺度的变化情况进行考察,并与已知的She-leveque（简称SL）标度律进行比较分析。结果表明,越接近壁面,结构函数的概率密度曲线尾部越宽,表现出逐渐增强的间歇性;并且随着壁面距离的减小,ESS标度律曲线与SL标度律曲线之间的差别也就越大,而新形式标度律曲线则与SL标度律曲线的差别逐渐减小。在提取的结构尺度一样时,结构函数的标度指数随法向位置的不同变化不大,但是随着提取结构尺度的增加,两种形式的标度律适用的尺度范围也在逐渐增大。     

A RESIDUAL-BASED UNRESOLVED-SCALE FINITE ELEMENT MODELLING FOR IMPLICT LARGE EDDY SIMULATION 1)

对应于湍流的大尺度与小尺度流场信息, 本文在有限元的框架下, 假设Navier-Stokes方程的解的形函数由大尺度和不可解尺度形函数叠加组成, 引入对应的权函数, 将Navier-Stokes方程的有限元变分形式分解为大尺度和不可解尺度系统. 根据不可解尺度系统, 构建基于Navier-Stokes大尺度方程残差的不可解尺度模型, 将其代入Navier-Stokes方程的大尺度系统, 进而数值求解大尺度系统得到Navier-Stokes方程的大尺度解. 该方法无需像传统的大涡模拟方法那样对方程的解进行过滤, 通过对形函数进行尺度分解实现解的尺度分解. 本文使用该方法的自编程序代码开展了槽道湍流的数值模拟. 通过与有限差分大涡模拟、DNS计算结果的比较, 发现在使用较少网格情况下该方法预测的平均流向速度在近壁区与DNS数据吻合, 在黏性外层略偏高; 该方法对雷诺应力预测偏低导致从流向向垂向方向上湍动能输运略偏低. 流向速度等值面图显示该方法有效捕捉到了大尺度旋涡结构; 同时在近壁区可以观察到明显的低速条带结构.     

基于Navier-Stokes方程残差的隐式大涡模拟有限元模型

对应于湍流的大尺度与小尺度流场信息, 本文在有限元的框架下, 假设Navier-Stokes方程的解的形函数由大尺度和不可解尺度形函数叠加组成, 引入对应的权函数, 将Navier-Stokes方程的有限元变分形式分解为大尺度和不可解尺度系统. 根据不可解尺度系统, 构建基于Navier-Stokes大尺度方程残差的不可解尺度模型, 将其代入Navier-Stokes方程的大尺度系统, 进而数值求解大尺度系统得到Navier-Stokes方程的大尺度解. 该方法无需像传统的大涡模拟方法那样对方程的解进行过滤, 通过对形函数进行尺度分解实现解的尺度分解. 本文使用该方法的自编程序代码开展了槽道湍流的数值模拟. 通过与有限差分大涡模拟、DNS计算结果的比较, 发现在使用较少网格情况下该方法预测的平均流向速度在近壁区与DNS数据吻合, 在黏性外层略偏高; 该方法对雷诺应力预测偏低导致从流向向垂向方向上湍动能输运略偏低. 流向速度等值面图显示该方法有效捕捉到了大尺度旋涡结构; 同时在近壁区可以观察到明显的低速条带结构.      

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

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

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

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

关于湍流标度律的争鸣

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

后台阶分离流动中大涡结构演变的数值模拟

本文对后台阶分离流动中涡结构的演变进行了大涡模拟，研究了流场结构的变化规律。详细讨论了随着雷诺数的增加流场结构的典型特征的变化规律，指出流场中的涡结构随着雷诺数的增大变得十分复杂和丰富，回流区的数目、大小及其出现的位置也显著地不同。这些结果与已有的一些实验值和流场显示结果是吻合的。在此基础上，进一步研究了高雷诺数时流场中大尺度涡结构的瞬时发展和演化过程，展示了其中大涡的产生、追随、吸引、合并和破碎等过程。对于高雷诺数情况，对大涡模拟得到的数值结果进行了统计，得到的时均速度分布以及台阶后方的回流区长度与现有的其他实验结果符合得很好。本研究是针对后台阶分离流动深入开展湍流控制以及两相流动研究的基础。     

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

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

大涡模拟的壁模型及其应用

大涡模拟是研究湍流的非定常特性的重要方法.但解析壁面层的大涡模拟所需的计算量与直接数值模拟相当,是大涡模拟在高雷诺数壁湍流数值模拟中所面临的主要困难.解析壁面层所需的网格尺度与壁面黏性长度同量级,是引起壁湍流大涡模拟计算量增加的主要原因.壁模型通过模化近壁流动避免了完全解析壁面层,可以显著地降低壁湍流大涡模拟的计算量,是克服上述困难的有效方法.本文介绍了大涡模拟壁模型的主要类型;详细讨论了常用的壁面应力模型,特别是平衡层模型和双层模型的构建思路和特点;基于近壁流动的特征讨论了应力边界条件的必要性和适用性;指出了壁面应力模型的局限性以及考虑非平衡效应修正的各种方法;讨论了壁面应力模型的研究历史、最新进展和发展趋势,给出了常用的壁面应力模型的分支与发展关系图;并基于Werner-Wengle模型实现了周期山状流的大涡模拟.     

The effect of a systematic change in surface roughness skewness on turbulence and drag

Previous work by the authors (Flack and Schultz, 2010) has identified the root-mean-square roughness height, k_rms, and the skewness, Sk, of the surface elevation distribution as important parameters in scaling the skin-friction drag on rough surfaces. In this study, three surfaces are tested in turbulent boundary layer flow at a friction Reynolds number, Re_τ = 1600–2200. All the surfaces have similar root-mean-square roughness height, while the skewness is varied. Measurements are presented using both two-component LDV and PIV. The results show the anticipated trend of increasing skin-friction drag with increasing skewness. The largest increase in drag occurs going from negative skewness to zero skewness with a more modest increase going from zero to positive skewness. Some differences in the mean velocity and Reynolds stress profiles are observed for the three surfaces. However, these differences are confined to a region close to the rough surface, and the mean velocity and Reynolds stress profiles collapse away from the wall when scaled in outer variables. The turbulence structure as documented through two-point spatial correlations of velocity is also observed to be very similar over the three surfaces. These results support Townsend’s (1976) concept of outer-layer similarity that the wall boundary condition exerts no direct influence on the turbulence structure away from the wall except in setting the velocity and length scales for the outer layer.     

Validation of large eddy simulation method for study of flatness and skewness of decaying compressible magnetohydrodynamic turbulence

In the present work we study potential applicability of large eddy simulation (LES) method for prediction of flatness and skewness of compressible magnetohydrodynamic (MHD) turbulence. The knowledge of these quantities characterizes non-Gaussian properties of turbulence and can be used for verification of hypothesis on Gaussianity for the turbulent flow under consideration. Prediction accuracy of these quantities by means of LES method directly determines efficiency of reconstruction of probability density function (PDF) that depends on used subgrid-scale (SGS) parameterizations. Applicability of LES approach for studying of PDF properties of turbulent compressible magnetic fluid flow is investigated and potential feasibilities of five SGS parameterizations by means of comparison with direct numerical simulation results are explored. The skewness and the flatness of the velocity and the magnetic field components under various hydrodynamic Reynolds numbers, sonic Mach numbers, and magnetic Reynolds numbers are studied. It is shown that various SGS closures demonstrate the best results depending on change of similarity numbers of turbulent MHD flow. The case without any subgrid modeling yields sufficiently good results as well. This indicates that the energy pile-up at the small scales that is characteristic for the model without any subgrid closure, does not significantly influence on determination of PDF. It is shown that, among the subgrid models, the best results for studying of the flatness and the skewness of velocity and magnetic field components are demonstrated by the Smagorinsky model for MHD turbulence and the model based on cross-helicity for MHD case. It is visible from the numerical results that the influence of a choice subgrid parametrization for the flatness and the skewness of velocity is more essential than for the same characteristics of magnetic field.     

基于人工神经网络的亚格子应力建模

亚格子(SGS)应力建模在湍流大涡模拟(LES)中有着极为重要的作用. 传统亚格子应力模型存在相对误差较大、耗散过强等问题. 近年来, 计算机技术的发展使得人工神经网络(ANN)等机器学习方法逐渐成为亚格子应力建模型的新研究范式. 本文着重考虑滤波宽度及雷诺数影响, 在不可压缩槽道湍流中建立了亚格子应力的ANN模型. 该模型以滤波后的直接数值模拟(fDNS)流场物理量及滤波尺度为输入信息, 相应滤波尺度下的亚格子应力为输出量. 通过对不同滤波尺度及不同雷诺数数据的训练, ANN模型能够给出与直接数值模拟(DNS)高度吻合的亚格子应力. 此外, 模型在亚格子耗散等非ANN建模量上也有着优异的预测性能, 与基于DNS获得的对应物理量的相关系数大都在0.9以上, 较梯度模型及Smagorinsky模型有明显提升. 在后验测试中, ANN模型对流向平均速度剖面的预测同样优于梯度模型、Smagorinsky模型及隐式大涡模拟(ILES)等传统LES模型. 在脉动速度均方根预测方面, 除了某些法向位置外ANN模型的性能整体上相对其他3个模型有所提升. 然而, 随着网格尺度的增大ANN模型预测的结果与fDNS结果的偏差逐渐增大. 总之, ANN方法在发展高精度亚格子应力模型上具有很大的潜力.      

External Intermittency Simulation in Turbulent Round Jets

Direct numerical and large eddy simulation (DNS and LES) are applied to study passive scalar mixing and intermittency in turbulent round jets. Both simulation techniques are applied to the case of a low Reynolds number jet with Re = 2,400, whilst LES is also used to predict a high Re = 68,000 flow. Comparison between time-averaged results for the scalar field of the low Re case demonstrate reasonable agreement between the DNS and LES, and with experimental data and the predictions of other authors. Scalar probability density functions (pdfs) for this jet derived from the simulations are also in reasonable accord, although the DNS results demonstrate the more rapid influence of scalar intermittency with radial distance in the jet. This is reflected in derived intermittency profiles, with LES generally giving profiles that are too broad compared to equivalent DNS results, with too low a rate of decay with radial distance. In contrast, good agreement is in general found between LES predictions and experimental data for the mixing field, scalar pdfs and external intermittency in the high Reynolds number jet. Overall, the work described indicates that improved sub-grid scale modelling for use with LES may be beneficial in improving the accuracy of external intermittency predictions by this technique over the wide range of Reynolds numbers of practical interest.     

Direct simulation of turbulent flow and heat transfer in a channel. Part I: Smooth walls

Interest in the use of supercomputers for the direct numerical calculation of turbulence prompts the development of efficient numerical techniques so that calculation at higher Reynolds numbers might be made. This paper presents an efficient pseudo-spectral technique, similar to but different from others that have recently appeared, for the calculation of momentum and heat transfer to a constant-property, turbulent fluid in a two-dimensional channel with walls at different, uniform temperature. The code uses no empiricism, although periodic boundary conditions are used for fluctuating quantities in the streamwise and spanwise directions. Calculations were made for a Prandtl number of 0·72 and Reynolds number based on friction velocity and channel half-height of 180 or 2800 based on channel half-height and average velocity. Calculations of mean velocity profile, turbulence intensities, skewness, flatness, Reynolds stress and eddy diffusivity of heat near a wall compare favourably with experimental results. Representative contour plots of the temperature field near the wall and of the spanwise and streamwise two-point velocity correlations are given. Deficiencies are that the calculation requires many hours on a fast computer with a large high-speed memory and that the grid size in each direction for appropriate resolution is approximately proportional to the square of the Reynolds number and to the Prandtl number raised to some power greater than one.     

Log-layer mismatch and commutation error in hybrid RANS/LES simulation of channel flow

Hybrid approach combining large eddy simulation (LES) with the Reynolds-averaged Navier–Stokes equation (RANS) is expected to accurately simulate wall-bounded turbulent flows at high Reynolds numbers. As an important issue in developing hybrid methods, it is known that the log layers in the RANS and LES regions are not lined up in hybrid RANS/LES simulations of channel flow. Although several methods including additional filtering near the RANS/LES interface have been proposed to eliminate the log-layer mismatch, there is no obvious physical justification for the methods and some ad hoc tuning is necessary. In this work, the commutation error terms in the filtered velocity equations are investigated to justify the method of additional filtering. It is shown that the additional filtering can be considered as a finite difference approximation to extra terms due to the non-commutivity between the hybrid filter and the spatial derivative. Moreover, an expression determining the filter width and its location for the additional filtering is obtained. To validate the expression, a hybrid simulation of channel flow is carried out. The additional filtering with the filter width derived is shown to be effective in eliminating the log-layer mismatch and improving the mean velocity profile.     

Assessment of the modulated gradient model in decaying isotropic turbulence

A recently introduced nonlinear model undergoes evaluations based on two isotropic turbulent cases: a University of Wiscosion-Madison case at a moderate Reynolds number and a Johns Hopkins University case at a high Reynolds number. The model uses an estimation of the subgrid-scale (SGS) kinetic energy to model the magnitude of the SGS stress tensor, and uses the normalized velocity gradient tensor to model the structure of the SGS stress tensor. Testing is performed for the first case through a comparison between direct numerical simulation (DNS) results and large eddy simulation (LES) results regarding resolved kinetic energy and energy spectrum. In the second case, we examine the resolved kinetic energy, the energy spectrum, as well as other key statistics including the probability density functions of velocities and velocity gradients, the skewness factors, and the flatness factors. Simulations using the model are numerically stable, and results are satisfactorily compared with DNS results and consistent with statistical theories of turbulence.     

Large Eddy Simulations: How to evaluate resolution

The present paper gives an analysis of fully developed channel flow at Reynolds number of Re=u_τδ/ν=4000 based on the friction velocity, u_τ, and half the channel height, δ. Since the Reynolds number is high, the LES is coupled to a URANS model near the wall (hybrid LES–RANS) which acts as a wall model. It it found that the energy spectra is not a good measure of LES resolution; neither is the ratio of the resolved turbulent kinetic energy to the total one (i.e. resolved plus modelled turbulent kinetic energy). It is suggested that two-point correlations are the best measures for estimating LES resolution. It is commonly assumed that SGS dissipation takes place at high wavenumbers. Energy spectra of the fluctuating velocity gradients show that this is not true; the major part of the SGS dissipation takes place at low to midrange wavenumbers. Furthermore, the energy spectra of the fluctuating velocity gradients reveals that the accuracy of the predicted velocity gradients at the highest resolved wavenumbers is very poor.     

PROBABILITY DISTRIBUTION FUNCTION OF NEAR-WALL TURBULENT VELOCITY FLUCTUATIONS

IntroductionNear-wall turbulence plays an important role in environmental engineering due to itsgreat impact on the transport of pollutants,nutrients,sediments,etc.Therefore,peoplehave long been studying the wall-bounded turbulence[1,2],in which the proba…     

Statistical Characteristics of an Isothermal, Supersonic Developing Boundary Layer Flow from DNS Data

Direct numerical simulation results for a developing, supersonic boundary layer flow with either an adiabatic wall temperature condition or a cold wall (relative to adiabatic) temperature condition are evaluated to assess the comparative effect on the mean and turbulent fields. Included in the analysis are two-point turbulent spanwise spatial correlations and the corresponding one-dimensional energy spectra distributions as well as higher-order statistics including skewness and flatness factors. In addition to the mean field velocity and temperature behavior, the turbulent Reynolds stress, temperature variance, and heat and mass flux distributions are discussed. An overall focus of the analysis is to both contrast the velocity and thermal field behavior and to provide some additional insight into the dynamic balance of the various statistical correlations and their impact on model development.