均匀各向同性湍流平衡区的标量谱

本文把von Kármán关于波数间相互作用的假设推广到标量的情况,得到标量谱传递函数的模化式,从而封闭了均匀各向同性湍流平衡区的谱方程。考虑到在平衡区的前部受到大尺度运动的影响,计入这一剪切效应可以合理地假设一维纵向动量谱的对数斜率达到—1时进入平衡区,由此可以根据实验数据及各向同性湍流一维谱和三维谱相互关系决定边界条件。通过数值方法解出的平衡区谱函数和实验结果相当一致。     

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

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

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

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

湍流数值计算的最近发展

湍流是有大量自由度的非线性力学系统,看来是流体的复杂的宏观不规则运动。因此,就是对于最简单的理想湍流——均匀各向同性湍流进行严格的理论分析,也是很困难的。湍流数值计算长期停留在半经验理论阶段。最近十几年来,高速电子计算机的应用,使湍流数值计算发生了巨大变化,根本改变了湍流问题的可解性。现在求出误差在百分之几以内的湍流精确数值结果已变成了现实。相应于计算机的发      

剪切湍流大尺度相干结构的模式研究

发展了一种计算剪切湍流大尺度相干结构的新模式．该模式的基础是认为大尺度相干结构为湍流场中流体脉动能量增长最快的那部分，且包含大部分的湍流脉动能量．在此基础上。通过对湍流相干能量方程的推演。建立了描述大尺度相干结构的特征控制方程，并应用Ｃｈｅｂｙｓｈｅｖ多项式方法求得湍流相干能量的最大增长率在波数空间的分布，从而获得对应的大尺度相干结构．应用该模式研究了槽流和一自然对流中的大尺度相干结构，得到的近壁区流动结构与实验现象十分接近．     

平板湍流边界层的声辐射

基于平板湍流边界层的壁压起伏波数—频率谱 ,给出了一种湍流边界层声辐射的估算方法 ,并对光滑平板湍流边界层和平板表面粗糙度引起的湍流边界层声辐射进行了分析。结果表明 :湍流边界层声辐射是一种四极子声辐射 ,且其辐射声能集中于平板表面粗糙度引起的湍流边界层声辐射 ;光滑平板湍流边界层的声辐射也不可忽略。     

均匀各向同性湍流的涡旋结构统计理论的一级近似数值解

本文数值计算了均匀各向同性湍流内涡旋结构统计理论的一级近似解。一级近似的结果比零级近似有改进,计算求得的二元纵向速度的相关系数和横向相关系数,都与实验值更吻合。     

8阶群速度控制格式及其应用

构造了8阶精度的群速度控制型差分格式.根据激波捕捉能力及对小尺度分辨能力对格式系数进行了优化.采用该格式对可压缩均匀各向同性湍流进行了直接数值模拟,所计算的最大湍流马赫数达到0.95.相同计算条件下的结果与他人结果吻合较好,说明了格式在捕捉激波的同时对湍流小尺度有较好的分辨能力.     

一种改进的均匀各向同性湍流初始化方法

均匀各向同性湍流是一种最简单的湍流理想状态,也是湍流基础理论研究的最重要对象之一.为了用数值方法产生均匀各向同性湍流场,一般采用Rogallo提出的方法在谱空间生成初始场,然后再转换到物理空间.研究表明,由该方法生成的初始湍流场在3个棱向上呈各向异性,在结构函数和速度概率密度分布上均有体现.尽管在初始场样本很多时,这种各向异性可以在平均意义上消除,但作为数值模拟采用的单个流场则波动较大,不利于在实际计算中作为单个初始场生成各向同性湍流.在此基础上提出一种改进的Rogallo初始化方法,称为模量平均法,将Rogallo方法在3个轴向分别进行,并进行模量平均,最后采用能谱进行模量控制.这种方法可以一方面保持初始场能谱,另一方面减小单个流场的各向异性波动,以产生各向同性程度更佳的单个初始场.在统计意义上,新方法可以分别将结构函数和速度概率密度的相对标准差减小约10%.     

上海陆域古河道分布及对工程建设影响研究

采用直接数值模拟方法, 对槽道湍流中展向振动流向传播的波动壁面的流动 控制和减阻问题进行了研究, 讨论了流向参数k_{x}对Stokes层、湍流拟序结构、湍流猝 发事件以及壁面阻力的影响, 并对此类波动壁面的湍流控制和减阻机理进行了讨论. 结果表 明, 当此类波动壁面被用来调制近壁流动时, 仅低频波对湍流流场具有显著影响, 可导致湍 流猝发事件的频率和强度的显著变化; 波数k_{x}的增大对于湍流猝发事件的频率和强度增 减的影响并不同步, 存在一个最优的波数k_{x}, 在其调制下, 固有流场对诱导流场的影响 最弱, 而诱导流场对固有流场的影响显著, 减阻效果最好.     

Modulation of homogeneous turbulence seeded with finite size bubbles or particles

The dynamics of homogeneous, isotropic turbulence seeded with finite sized particles or bubbles is investigated in a series of numerical simulations, using the force-coupling method for the particle phase and low wavenumber forcing of the flow to sustain the turbulence. Results are given on the modulation of the turbulence due to massless bubbles, neutrally buoyant particles and inertial particles of specific density 1.4 at volumetric concentrations of 6%. Buoyancy forces due to gravity are excluded to emphasize finite size and inertial effects for the bubbles or particles and their interactions with the turbulence. Besides observing the classical entrapment of bubbles and the expulsion of inertial particles by vortex structures, we analyze the Lagrangian statistics for the velocity and acceleration of the dispersed phase. The turbulent fluctuations are damped at mid-range wavenumbers by the bubbles or particles while the small-scale kinetic energy is significantly enhanced. Unexpectedly, the modulation of turbulence depends only slightly on the dispersion characteristics (bubble entrapment in vortices or inertial sweeping of the solid particles) but is closely related to the stresslet component (finite size effect) of the flow disturbances. The pivoting wavenumber characterizing the transition from damped to enhanced energy content is shown to vary with the size of the bubbles or particles. The spectrum for the energy transfer by the particle phase is examined and the possibility of representing this, at large scales, through an additional effective viscosity is discussed.     

On the evolution of decaying isotropic turbulence

A direct numerical simulation is performed on 256³ grids for decaying isotropic turbulence. The total kinematic energy, Taylor micro-scale, Taylor micro-scale Reynolds number and the velocity derivative skewness are calculated. The snapshots of energy spectra and energy transfer spectra are plotted. These measurements verify the DIA predictions: decaying isotropic turbulence has the energy propagation and occupies the final decay periods. The skewness remains to some level with small variation even in the final decay period.     

Band-pass filtered velocity statistics in decaying turbulent box

For homogeneous isotropic turbulence study, the acquisition of band-pass filtered velocity increments (FVI) in a non-forced turbulent box is still a challenge both experimentally and numerically. Turbulence and associated physical processes, at a given instant, are permanently contaminated by a forcing process which can seldom be universal. The situation tends to be the origin of intermittency and the non-Gaussian probability density distribution for acceleration and velocity gradients. To reveal implied mechanism, grid turbulence is adapted to observe non-perturbed homogeneous isotropic turbulence. The velocity increments (VI) can be obtained following Comte-Bellot and Corrsin (GCBC) by means of two point-two time shifted velocity measurements. It is difficult to obtain decaying turbulence (DT) at large turbulent Reynolds number without pollution coming from walls. Nevertheless it is also significant to investigate DT in low Reynolds number regimes to determine non-polluted tendencies. The similarity of DT between particle image velocimetry (PIV) and hot wire anemometry measurements by GCBC are presented. Here we focus our tendency on VI and FVI probability density function (PDF) shapes in this letter. In conclusion, the tendency to Gaussian shape in inertial zone wavenumbers, demonstrates that there will be no intermittency if turbulent cascade is not perturbed.     

数值耗散对压气机流动分离涡模拟的影响研究

分离涡模拟DES是压气机流动模拟中常用的高保真湍流模式。为了使DES准确解析湍流,数值耗散必须限制在合理范围内。然而,当前的压气机流动DES类研究工作中仍然普遍采用高耗散的迎风格式。本文首先基于DES类方法计算的各向同性衰减湍流结果,定量比较了多种不同数值格式的耗散,证实了高耗散迎风格式严重低估中高波数湍流能量。高阶重构格式可以一定程度上改善该问题,但能量耗散仍然过高。本文在高阶重构的基础上,进一步引入自适应耗散函数修改Riemann求解器,构造了自适应耗散格式。该格式在全波数范围都能准确地预测湍流能谱。将该格式配合DES类方法模拟跨声速离心压气机流动,其预测的压比相比于三阶迎风格式,更加接近实验结果。此外,自适应耗散格式显著提高了中小尺度流动结构的分辨率。分析表明,在使用DES类方法模拟压气机流动时,有必要采用数值耗散较低的离散格式,以准确预测压气机总体性能和流动结构。本文构造的自适应耗散格式是一种良好选择。     

Experimental estimation of fluctuating velocity and scalar gradients in turbulence

The effect of numerical differentiation is investigated in the context of evaluating fluctuating velocity and scalar quantities in turbulent flows. In particular, 2-point forward-difference and 3-, 5-, 7-, and 9-point centred-difference schemes are investigated. The spectral technique introduced by Wyngaard (in J Sci Instr 1(2):1105–1108, 1968) for homogeneous turbulence is used to quantify the effects of the schemes. Numerical differentiation is shown to attenuate gradient spectra over a range of wavenumbers. The spectral attenuation, which varies with the order of the scheme, results in a reduction in the measured mean-squared gradients. High-order schemes (e.g. 7- or 9-point) are shown to significantly decrease the attenuation at all wavenumbers and as a result produce more accurate gradients. Hot-wire measurements and direct numerical simulations of decaying homogeneous, isotropic turbulence are found to be in good agreement with the predictions of the analysis, which suggests that high-order schemes can be used to improve empirical gradient estimates. The shape of the probability density functions is also found to be sensitive to the choice of numerical differentiation scheme. The effect of numerical differentiation is also discussed with respect to particle image velocimetry (PIV) measurements of a nominally two-dimensional planar mixing layer. It is found that the relatively low signal-to-noise ratio inherent in typical PIV measurements necessitates the use of low-order schemes to prevent excessive noise amplification, which increases with the order of the scheme. The results of the present work demonstrate that high-order numerical differentiation schemes can be employed to more accurately resolve gradients measured at a given resolution provided the measurements have an adequate signal-to-noise ratio.     

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

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

The reduction of spatial aliasing by long hot-wire anemometer probes

 Experiment and numerical analysis are presented to demonstrate that a hot-wire anemometer probe reduces spatial aliasing of turbulent velocity fluctuations because of the filtering property of the probe sensing element. The experiment focuses on the one-dimensional turbulent velocity spectrum and utilizes a long sensing length hot-wire probe to exaggerate the effect of the sensing element on the turbulent field. The numerical analysis utilizes a model of the hot-wire probe from Wyngaard (1968) along with isotropic turbulence relations to obtain an equation for the hot-wire response in a turbulent velocity field. The model can be used to determine the effect of hot-wire length on the one and three-dimensional turbulent spectra. The experimental study demonstrates that the finite length, hot-wire probe filters out energy in the high wave number region of the one-dimensional spectrum thereby verifying its ability to reduce spatial aliasing. Interestingly, the study also shows that energy in the low wave numbers of the one-dimensional spectrum is attenuated. The numerical study of the hot-wire probe demonstrates that this low wave-number attenuation is purely an artifact of the one-dimensional spectrum and not an effect of the hot-wire probe. Received: 20 May 996 / Accepted: 14 November 1996     

Performance of a probe for measuring turbulent energy and temperature dissipation rates

A probe consisting of four X-wires (a total of eight hot wires operated in constant temperature mode) and two pairs of parallel cold wires (operated in constant current mode) is in principle capable of providing simultaneous data for the instantaneous energy and temperature dissipation rates. To evaluate the performance of this probe, measurements have been made in decaying turbulence downstream of a grid/screen combination. The directly measured mean values of the energy and temperature dissipation rates are compared with those obtained from the streamwise decay rates of the mean turbulent energy and temperature variance. The probe also yields all three fluctuating vorticity components; after applying spatial resolution corrections, their spectra are in close agreement with isotropic calculations over nearly all wavenumbers. Both the vorticity variance and the mean energy dissipation rate exhibit the same power-law decay rate. Transport equations for the mean energy and temperature dissipation rates are satisfied to within ᆞ%.     

On the correlation between enstrophy and energy dissipation rate in a turbulent wake

All three components of the vorticity fluctuation have been measured simultaneously in a turbulent wake using a new eight-sensor vorticity probe. The vorticity fluctuation spectra agree reasonably well with those from a direct numerical simulation of a turbulent channel flow at high wavenumbers. The similarity between the instantaneous energy dissipation rate ε and the instantaneous enstrophy ω² is examined using spectra and probability density functions. The correlation between ω² and ε is evaluated in some detail. The homogeneous value of ε is strongly correlated with ω². The full value of ε and, more especially its isotropic value, are less well correlated with the enstrophy. Conditional averaging indicates that high enstrophy regions are associated with high energy dissipation rate regions.     

Error quantification of 3D homogeneous and isotropic turbulence measurements using 2D PIV

Particle image velocimetry (PIV) has become a popular non-intrusive tool for measuring various types of flows. However, when measuring three dimensional flows with 2D PIV, there is inherent measurement error due to out-of-plane motion. Errors in the measured velocity field propagate to turbulence statistics. Since this can distort the overall flow characteristics, it is important to understand the effect of this out-of-plane error. In this study, the effect of out-of-plane motion on turbulence statistics is quantified. Using forced isotropic turbulence direct numerical simulation (DNS) flow field data provided by the Johns Hopkins turbulence database (JHTDB), synthetic image tests are performed. Turbulence statistics such as turbulence kinetic energy, dissipation rate, Taylor microscale, Kolmogorov scale, and velocity correlations are calculated. Various test cases were simulated while controlling three main parameters which affect the out-of-plane motion: PIV interrogation window size, camera inter-frame time, and laser sheet thickness. The amount of out-of-plane motion was first quantified, and then the error variation according to these parameters was examined. This information can be useful when examining fully three dimensional flows such as homogeneous and isotropic turbulence via 2D PIV.