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

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

不可压N-S方程高效算法及二维槽道湍流分析

构造了基于非等距网格的迎风紧致格式,并将其与三阶精度的Adams半隐方法相结合,构造了求解不可压N－S方程高效算法。该算法利用基于交错网格的离散形式的压力Poisson方程求解压力项,解决了边界处的残余散度问题;同时还利用快速Fourier变换将方程的隐式部分解耦,离散后的代数方程组利用追赶法求解,大大减少了计算量。通过对二维槽道流动的数值模拟,证实了所构造的数值方法具有精度高,稳定性好,能抑制混淆误差等优点,同时具有很高的计算效率,是进行壁湍流直接数值模拟的有效方法。在数值模拟的基础上对二维槽道流动进行了分析,得到了Reynolds数从6000到15000的二维流动饱和态解（所谓“二维槽道湍流”）;定性及定量结果均与他人的数值计算结果吻合十分理想。对流场进行了分析,指出了“二维湍流”与三维湍流统计特性的区别。     

壁湍流相干结构尺度及边界层内的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 标度律成立的范围并不完全一致，同时标度律成立的范围大小与流场尺度有明显关系.     

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

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

关于湍流标度律的争鸣

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

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

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

深槽表面气体静压润滑低Reynolds数层流节流效应实验研究

气体介质在润滑间隙流动过程中,沿气体流动方向开设的表面矩形深槽结构内部产生旋涡阻碍气体流动,形成节流效应.通过矩形深沟槽表面静压润滑实验,开展深槽表面和光滑表面的流量特性对比测试,研究了低Reynolds数层流状态下深槽表面间隙的气体节流效应.结果表明:深槽表面产生显著的节流效应,并且随Reynolds数增加而增强;在低Reynolds数层流状态下,深槽表面可以产生节流效应,但是节流效应强度不稳定,当处于完全湍流状态时,节流效应维持定值;间隙尺寸、槽深等结构参数对节流效应影响明显.     

可压缩各向同性衰减湍流直接数值模拟研究

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

）IRECT NUMERICAL SIMULATION OF DECAYING COMPRESSIBLE ISOTROPIC TURBULENCE

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

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

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

Symmetries, Invariance and Scaling-Laws in Inhomogeneous Turbulent Shear Flows

An approach to derive turbulent scaling laws based on symmetry analysis is presented. It unifies a large set of scaling laws for the mean velocity of stationary parallel turbulent shear flows. The approach is derived from the Reynolds averaged Navier–Stokes equations, the fluctuation equations, and the velocity product equations, which are the dyad product of the velocity fluctuations with the equations for the velocity fluctuations. For the plane case the results include the logarithmic law of the wall, an algebraic law, the viscous sublayer, the linear region in the centre of a Couette flow and in the centre of a rotating channel flow, and a new exponential mean velocity profile that is found in the mid-wake region of high Reynolds number flat-plate boundary layers. The algebraic scaling law is confirmed in both the centre and the near wall regions in both experimental and DNS data of turbulent channel flows. For a non-rotating and a moderately rotating pipe about its axis an algebraic law was found for the axial and the azimuthal velocity near the pipe-axis with both laws having equal scaling exponents. In case of a rapidly rotating pipe, a new logarithmic scaling law for the axial velocity is developed. The key elements of the entire analysis are two scaling symmetries and Galilean invariance. Combining the scaling symmetries leads to the variety of different scaling laws. Galilean invariance is crucial for all of them. It has been demonstrated that two-equation models such as the k– model are not consistent with most of the new turbulent scaling laws.     

An Attempt to Realize Experimental Isotropic Turbulence at Low Reynolds Number

This paper presents an attempt to realize experimental isotropicturbulence at low Reynolds number. For this aim an experimentalapparatus, a turbulence chamber “Box”, was designed and built togenerate a turbulent flow field in the center of the chamber. Theturbulent airflow field was generated by eight electrical fans placedsymmetrically at the eight internal corners of the externally cubicchamber. The turbulence intensity was controlled by the fans speed.Laser Doppler velocimeter (LDV) in single and two-point velocitymeasurements was used to fully characterize the turbulent field insidethe chamber. The main results indicate that the turbulence ishomogeneous and isotropic with a quasi-zero mean velocity within aspherical region of 20 mm radius from the center of the chamber. Themeasured turbulent integral length scale was found to be constant andindependent of the turbulence intensity (or fans speed). Furthermore, anoticeable spectral inertial subrange as prescribed by the Kolmogorovtheory has not been observed at the range of Reynolds number exploredhere, where Re_λ < 100. But rather a scaling region characterized by anexponent that is lower than the Kolmogorov value, ?5/3, has beenidentified. Moreover, the value of this exponent showed no definedtrend, while the width of the inertial scaling region expands as themicroscale Reynolds number increases.     

Multipoint correlations of concentration fluctuations in a turbulent passive scalar field

This study experimentally analyzes a turbulent passive scalar field using two-point and three-point correlation functions of the concentration fluctuations. The scalar field was created by the iso-kinetic release of a high Schmidt number dye into a fully developed turbulent boundary layer of an open channel flow. Concentration data at spatially separated measurement points were collected via the planar laser induced fluorescence (PLIF) technique. The current study complements previous research efforts by examining three-point correlations for several configurations in a turbulent shear flow. In the case of two streamwise-aligned points combined with one transversely separated point, contours of the three-point correlation function exhibit the symmetric properties reported in an earlier study of non-shear flow. In a second set of three-point configurations consisting of isosceles and collinear geometries, the influence of the orientation angle between the three-point configuration and the mean concentration gradient varies depending on the specified three-point geometry. The results also suggest that the scaling exponent in the inertial-convective regime is dependent on the injection length scale with weak dependence on the Reynolds number.     

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.     

Turbulent separated reattached flow in a two-dimensional curved-wall diffuser

A turbulent separation-reattachment flow in a two-dimensional asymmetrical curved-wall diffuser is studied by a two-dimensional laser doppler velocimeter. The turbulent boundary layer separates on the lower curved wall under strong pressure gradient and then reattaches on a parallel channel. At the inlet of the diffuser, Reynolds number based on the diffuser height is 1.2×10⁵ and the velocity is 25.2m/s. The results of experiments are presented and analyzed in new defined streamline-aligned coordinates. The experiment shows that after Transitory Detachment Reynolds shear stress is negative in the near-wall backflow region. Their characteristics are approximately the same as in simple turbulent shear layers near the maximum Reynolds shear stress. A scale is formed using the maximum Reynolds shear stresses. It is found that a Reynolds shear stress similarity exists from separation to reattachment and the Schofield-Perry velocity law exists in the forward shear flow. Both profiles are used in the experimental work that leads to the design of a new eddy-viscosity model. The length scale is taken from that developed by Schofield and Perry. The composite velocity scale is formed by the maximum Reynolds shear stress and the Schofield-Perry velocity scale as well as the edge velocity of the boundary layer. The results of these experiments are presented in this paper.     

Large-scale computer simulation of fully developed turbulent channel flow with heat transfer

Recently, with the advent of supercomputers, there has been considerable interest in the use of direct numerical simulation to obtain information about turbulent shear flow at low Reynolds number. This paper presents a pseudospectral technique to solve the full three-dimensional time-dependent Navier-Stokes and advection-diffusion equations without the use of subgrid-scale modelling. The technique has not been previously used for fully developed turbulent channel flow simulation and is based on methods applied in other contexts. The emphasis of this paper is to provide a reasonably detailed account of how the simulation is done rather than to present new calculations of turbulence. The details of an algorithm for turbulent channel flow simulation and the grid and time step sizes needed to integrate through transient behaviour to steady state turbulence have not been published before and are presented here. Results from a Cray-2 simulation of fully developed turbulent flow in a channel with heat transfer are presented along with a critical comparison between experiment and computation. The first- and second-order moments agree well with experimental measurements; the agreement is poor for higher-order moments such as the skewness and flatness near the walls of the channel. Detailed information given about the effects of spatial grid resolution on a computed results is important for estimating the size of the computation required to study various aspects of a turbulent flow.     

Low Reynolds number effects in open-channel turbulent boundary layers

The present study reports measurements of a turbulent boundary layer in an open-channel flow using fiber-optic laser Doppler anemometry. The Reynolds numbers based on momentum thickness and depth of flow are in the range 750≤Re _θ ≤2,400 and 15,300≤Re _h ≤54,200, respectively. It is shown that an accurate estimate of the wall shear stress can be made by fitting a fifth-order polynomial to the near-wall data. The effect of Reynolds number on the mean turbulence intensity and triple correlation is examined using both conventional scaling laws and the recent scaling laws proposed by George and Castillo. The present results show that different scaling laws lead to different conclusions on low Reynolds number effects.     

Reynolds number effects on a turbulent boundary layer with separation,reattachment, and recovery

The present paper addresses experimental studies of Reynolds number effects on a turbulent boundary layer with separation, reattachment, and recovery. A momentum thickness Reynolds number varies from 1,100 to 20,100 with a wind tunnel enclosed in a pressure vessel by varying the air density and wind tunnel speed. A custom-built, high-resolution laser Doppler anemometer provides fully resolved turbulence measurements over the full Reynolds number range. The experiments show that the mean flow is at most a very weak function of Reynolds number while turbulence quantities strongly depend on Reynolds number. Roller vortices are generated in the separated shear layer caused by the Kelvin–Helmholtz instability. Empirical Reynolds number scalings for the mean velocity and Reynolds stresses are proposed for the upstream boundary layer, the separated region, and the recovery region. The inflectional instability plays a critical role in the scaling in the separated region. The near-wall flow recovers quickly downstream of reattachment even if the outer layer is far from an equilibrium state. As a result, a stress equilibrium layer where a flat-plate boundary layer scaling is valid develops in the recovery region and grows outward moving downstream.     

Development of a turbulent near-wall temperature model and its application to channel flow

A numerical study of fluid flow and heat transfer in a two-dimensional channel under fully developed turbulent conditions is reported. A computer program which is capable of treating both forced and natural convection problems under turbulent conditions has been developed. The code uses the high-Reynolds-number form of the two equation turbulent model(k-?) in which a turbulent kinetic energy near-wall model is incorporated in order to accurately represent the behavior of the flow near the wall, particularly in the viscous sublayer where the turbulent Reynolds number is small. A near-wall temperature model has been developed and incorporated into the energy equation to allow accurate prediction of the temperature distribution near the wall and, therefore, accurate calculation of heat transfer coefficients. The sensitivity of the prediction of flow and heat transfer to variations in the coefficients used in the turbulence model is investigated. The predictions of the model are compared to available experimental and theoretical results; good agreement is obtained. The inclusion of the near-wall temperature model has further improved the predictions of the temperature profile and heat transfer coefficient. The results indicate that the turbulent kinetic energy Prandtl number should be a function of Reynolds number.