吹吸扰动对壁湍流相干结构的影响

采用热线风速仪测量受吹吸扰动的壁湍流边界层的流向速度，用傅里叶变换和子波变换研究 吹吸扰动对壁湍流能谱的影响，结果显示施加的低频扰动使边界层内层大尺度结构的能量减 少，小尺度结构的能量有所增强，远离壁面时扰动强度逐步衰减直到在外层中消失；通过VITA 法和子波变换法检测猝发事件，表明该扰动降低了猝发强度，使猝发周期延长，条件平均速 度波形的幅值降低、持续时间变短，说明扰动明显抑制了相干结构的猝发过程. 利用子波变 换可以实现湍谱分析，能有效检测猝发中的湍流结构，是一种客观的分析工具.     

The eduction of structures from flow imagery using wavelets Part I. The mixing layer

 Using two-dimensional Mexican hat wavelets, digitized imagery from the flow visualization pictures of a mixing layer by Brown and Roshko have been analysed. The objective of the present study is to “calibrate” the proposed wavelet technique for educing structures in a flow whose large-scale organization is well established. It is found that, with appropriate thresholds on wavelet-transform coefficients, considerable insight into the structure of the flow at different scales can be obtained. In particular, the analysis reveals that the small-scale structure within the large-scale vortices is approximately homogeneous, and has characteristic scales that do not vary significantly with downstream distance. Received: 5 March 1998 / Accepted: 28 October 1998     

Decay of confined,two-dimensional,spatially periodic arrays of vortices: A numerical investigation

The disarrangement of a perturbed lattice of vortices was studied numerically. The basic state is an exponentially decaying, exact solution of the Navier-Stokes equations. Square arrays of vortices with even numbers of vortex cells along each side were perturbed and their evolution was investigated. Whether the energy in the perturbation grows somewhat before it decays or decays monotonically depends on the initial strength of the vortices of the basic state, the extent of lateral confinement and the structure of the perturbation. The critical condition for temporally local instability, i.e. the critical amplitude of the basic state that must be exceeded to allow energy transfer from the basic state to the perturbation, is discussed. In the strongly confined case of a square lattice of four vortices the appearance of enchancement of global rotation is the result of energy transfer from the basic state to a temporally local unstable mode. Energy is transferred from the basic state to larger-scaled structures (inverse cascade) only if the scales of the larger structures are inherently contained in the initial structure of the perturbation. The initial structure of the double array of vortices is not maintained except for a very special form of perturbation. The facts that large scales decay more slowly than small scales and that, when non-linearities are sufficiently strong, energy is transferred from one scale to another explain the differences in the disarrangement process for different initial strengths of the vortices of the basic state. The stronger vortices, i.e. the vortices perturbed in a manner that increases their strength, tend to dominate the weaker vortices. The pairing and subsequent merging (or capture) of vortices of like sense into larger-scale vortices are described in terms of peaks in the evolution of the square root of the palinstrophy divided by the enstrophy.     

可压缩燃烧反应转捩混合层直接数值模拟

针对三维时间发展可压缩氢/氧非预混燃烧反应平面自由剪切混合层,采用5阶迎风/6阶对称紧致混合差分格式以及3阶显式Runge-Kutta时间推进方法,直接数值模拟了伴随燃烧产物生成和反应能量释放, 流动受扰动激发失稳并转捩的演化过程. 在转捩初期, 获得了${\it\Lambda}$涡、马蹄涡等典型的大尺度拟序结构,观察到了流动失稳后发生双马蹄涡三维对并的现象, 大尺度结构呈较好的对称性.在流动演化后期, 大尺度结构逐次破碎形成小尺度结构, 混合层进入转捩末期,呈明显的不对称性.      

混合层强化混合的数值研究

受 Ｗａｎｇ ＆ Ｆｉｅｄｌｅｒ（１９９７）的实验的启发，采用高阶精度的差分格式，通过数值模拟的方法，研究了二维混合层及限于两平板间的二维混合层（二维受限混合层）入口处加振动对提高混合层混合效率的作用．计算结果表明：对二维混合层，振动的频率越低，在混合层中产生的大尺度涡结构的尺度越大，在频率很低时，涡具有相似性；对限于两平板间的二维混合层，在一定的振动频率下，混合层中产生的涡较大而且破碎得也较好，这将有利于混合．这一结论与 Ｗａｎｇ ＆ Ｆｉｅｄｌｅｒ（１９９７）的实验观测到的结果是一致的．     

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

Numerical study on the turbulent mixing of planar shock-accelerated triangular heavy gases interface

The interaction of a planar shock wave with a triangle-shaped sulfur hexafluoride (\(\mathrm{SF_6}\)) cylinder surrounded by air is numerically studied using a high resolution finite volume method with minimum dispersion and controllable dissipation reconstruction. The vortex dynamics of the Richtmyer–Meshkov instability and the turbulent mixing induced by the Kelvin–Helmholtz instability are discussed. A modified reconstruction model is proposed to predict the circulation for the shock triangular gas–cylinder interaction flow. Several typical stages leading the shock-driven inhomogeneity flow to turbulent mixing transition are demonstrated. Both the decoupled length scales and the broadened inertial range of the turbulent kinetic energy spectrum in late time manifest the turbulent mixing transition for the present case. The analysis of variable-density energy transfer indicates that the flow structures with high wavenumbers inside the Kelvin–Helmholtz vortices can gain energy from the mean flow in total. Consequently, small scale flow structures are generated therein by means of nonlinear interactions. Furthermore, the occasional “pairing” between a vortex and its neighboring vortex will trigger the merging process of vortices and, finally, create a large turbulent mixing zone.     

利用子波分析对平壁湍流猝发现象的研究

利用槽道湍流直接数值模拟的数据库,采用子波分析的方法。对平壁湍流猝发现象的多尺度特性进行了研究,在不同惊讶上对猝发平均周期进行了统计,并利用局部标度指数研究了猝发过程的奇异性。     

混合层中柱状粒子取向分布函数的研究

数值模拟了平面流固混合层中柱状粒子在平面取向的分布函数，提出了描述流场中某点直粒子取向占优程度的参数Dψ的分析，说明在混合层中心区域粒子的取向排列是有序的，在混合层相邻两个大涡之间的区域，是粒子最可几角占优程度最高的地方，而且最可几取向就是沿流场的拉伸方向，说明流场的拉伸作用及方向对粒子的取向排列具有重要影响。与其它模型相比，粒子取向的分布函数模型能定量反映流场大涡对粒子取向控制的强弱程度。     

子波分析辨识壁湍流猝发事件的能量最大准则

用子波分析的方法，对用热膜测速仪得到的平板湍流边界层中流向脉动速度信号，在时域空间和频域空间同时进行时频双局部化分解．用子波系数研究了壁湍流脉动动能随尺度的分布，提出了确定壁湍流猝发事件时间尺度参数的能量最大准则，用子波逆变换得到了猝发事件对应的速度信号波形     

Vortex motion in a rotating barotropic fluid layer

To study vortex motion and the mechanisms of geostrophic adjustment (i.e. the equilibrium between pressure gradient and Coriolis force, which leads to the weakening of inertio-gravity waves) in large scale geophysical flows, we simulate the dynamics of a shallow-water layer in uniform rotation, without any forcing other than the initial injection of energy and potential enstrophy. Such a flow generates inertio-gravity waves which interact with the rotational eddies. We found that both inertio-gravity waves and rotation reduce the non-linear interactions between vortices, namely the condensation of the vorticity field into isolated coherent vortices, corresponding to the inverse rotational energy cascade, and the associated production of vorticity filaments, due to the direct potential enstrophy cascade. Rotation also inhibits the direct inertio-gravitational energy cascade for scales larger than the Rossby deformation radius. Therefore, if inertio-gravity waves are initially excited at large enough scales, they will remain trapped there due to rotation and there will be no geostrophic adjustment. On the contrary, if inertio-gravity waves are only present at scales smaller than the Rossby deformation radius, which are insensitive to the effect of rotation, they will non-linearly interact and cascade towards the dissipative scales, leaving the flow in geostrophic equilibrium.     

On the evolution of large scale structures in three-dimensional mixing layers

In this paper, several mathematical models for the large scale structures in some special kinds of mixing layers, which might be practically useful for enhancing the mixing, are proposed. First, the linear growth rate of the large scale structures in the mixing layers was calculated. Then, using the much improved weakly non-linear theory, combined with the energy method, the non-linear evolution of large scale structures in two special mixing layer configurations is calculated. One of the mixing layers has equal magnitudes of the upstream velocity vectors, while the angles between the velocity vectors and the trailing edge were π/2-φ and π/2+φ, respectively. The other mixing layer was generated by a splitter-plate with a 45-degree-sweep trailing edge. The project supported by the National Natural Science Foundation of China (19642001) and Deutsche Forschungsgemeinschaft (DFG)     

Effects of free‐stream turbulence on large‐scale coherent structures of separated boundary layer transition

It has been well established that large‐scale structures, usually called coherent structures, exist in many transitional and turbulent flows. The topology and range of scales of those large‐scale structures vary from flow to flow such as counter‐rotating vortices in wake flows, streaks and hairpin vortices in turbulent boundary layer. There has been relatively little study of large‐scale structures in separated and reattached transitional flows. Large‐eddy simulation (LES) is employed in the current study to investigate a separated boundary layer transition under 2% free‐stream turbulence on a flat plate with a blunt leading edge. The Reynolds number based on the inlet free stream velocity and the plate thickness is 6500. A dynamic subgrid‐scale model is employed to compute the subgrid‐scale stresses more accurately in the current transitional flow case. Flow visualization has shown that the Kelvin–Helmholtz rolls, which have been so clearly visible under no free‐stream turbulence (NFST) are not as apparent in the present study. The Lambda‐shaped vortical structures which can be clearly seen in the NFST case can hardly be identified in the free‐stream turbulence (FST) case. Generally speaking, the effects of free‐stream turbulence have led to an early breakdown of the boundary layer, and hence increased the randomization in the vortical structures, degraded the spanwise coherence of those large‐scale structures. Copyright © 2005 John Wiley & Sons, Ltd.     

Velocity field after wave breaking

The velocity field in breaking water waves is considered in this paper. A numerical simulation describes in detail the transition from a primary overturning and consequent rebounding jets into a bore front, where the vorticity in the coherent large‐scale eddy structures devolves into turbulence. Spatial changes in the frequency spectra of the kinetic energy and the enstrophy are associated with the production, transport and dissipation of the Reynolds stress and the various wave and turbulent mixing length scales. Mean velocity fields and the wave and kinetic energy in a surf zone are evaluated. Fourier and wavelet spectral analysis is applied to study both the surface elevation and energy changes, and the distinction that must be made between spilling and plunging breakers is clarified in this paper. Copyright © 2002 John Wiley & Sons, Ltd.     

Effect of Streamline Curvature on Intensity of Streamwise Vortices in the Mixing Layer of Supersonic Jets

The structure of supersonic nonisobaric jets with Mach numbers M_a = 1 and 2 is considered experimentally to find the effect of streamline curvature on the evolution of streamwise vortices in the mixing layer. The spatial development of steady streamwise vortices in the mixing layer of supersonic jets is considered. A method for generation of steady streamwise vortices by applying microroughness elements of controlled size onto the inner surface of the nozzle is developed. Radial profiles and azimuthal variations of total pressure are obtained; the mixinglayer thickness and the curvature of streamlines in supersonic jets are determined. A significant effect of microroughness elements of prescribed shape located on the nozzle surface on the behavior of total pressure in the mixing layer of supersonic jets, as compared to natural disturbances, is obtained.     

Coherent structure identification from wavelet analysis of particle image velocimetry data

 A new technique based on wavelet transform is applied to bidimensional velocity fields obtained by particle image velocimetry (PIV) measurements, in order to extract and characterize swirling motion associated with coherent structures. The proposed technique is based on the selectivity property of the wavelet transform and permits the detection of regions of the flow field associated with coherent structures and their spatial localization. Furthermore, being the method based on the analysis of the local energy content at separated scales, it is possible to extract the typical wavenumber associated with structures and therefore the typical length-scale. The procedure is validated by the application to velocity vector fields obtained from PIV measurements in different flow conditions and turbulence levels. Results are compared with those obtained by other more standard procedures, and the advantages and limitations of the proposed method are then discussed. Received: 16 October 2000 / Accepted: 18 June 2001 Published online: 29 November 2001     

Wavelet spatial scaling for educing dynamic structures in turbulent open cavity flows

In this work a methodology was developed for the selection of wavelet spatial scales to educe dynamic structures in turbulent cavity flows. The wavelet transform was applied to both the temporal signal and spatial fields to extract structures from the oscillating shear layer. The dominant frequencies were identified from the temporal transform of the shear layer oscillations, and then the corresponding wavelength was obtained using the relation U_cT=λ at each frequency. The wavelet spatial scaling was examined and a one-to-one relationship was established with respect to the wavelength. At each spatial scale, the transformed images of vorticity, velocity, and pressure fluctuations captured the vortical structures. Using this methodology, the dynamic vortical structures were extracted from the turbulent open cavity flows. Energy analysis was performed to examine the contributions of each structure.     

The structure and dynamics of reacting plane mixing layers

The reacting two-dimensional plane mixing layer has been studied in two configurations: a rearward facing step and a two-stream mixing layer. Observations have been made of the steady state behavior, and the unsteady behavior when the flow was forced by a specific acoustic frequency. The steady behavior of the mean properties of the reacting flows is similar to that of non-reacting free shear flows except for the global effects of thermodynamic property changes. The structure of these flows is qualitatively similar to that of non-reacting flows. Vortices form by the two-dimensional Kelvin-Helmholtz instability and grow by subharmonic combination until the mixing layer interacts with the walls. Entrainment is dominated by the two-dimensional vortex motion. Three-dimensional instabilities give rise to secondary vortices which are coherent over several Kelvin-Helmholtz structures and dominate the fine scale mixing process. The mixing transition corresponds to a loss of coherence in the layer. Unsteady behavior occurs when there are resonant interactions with the Kelvin-Helmholtz instability or the instability associated with the recirculation vortex in the rearward facing step flow. Modeling efforts are reported which show promise of simulating the essential features of plane mixing layers.A version of this paper was presented at the ASME Winter Annual Meeting of 1984 and printed in AMD-Vol. 66     

Active Control of Swirling Coaxial Jet Mixing with Manipulation of Large-Scale Vortical Structures

Large-scale vortical structures and associated mixing in methane/air swirling coaxial jets are actively controlled by manipulating the outer shear layer of the outer swirling coaxial jet with miniature flap actuators. In order to investigate the control mechanisms, stereoscopic particle image verocimetry (stereo-PIV) and plannar laser-induced fluorescence (PLIF) techniques are employed. It is found that intense vortex rings are produced in the outer shear layer in phase with the periodic flap motion regardless of the swirl number examined. The vortical structures in the inner shear layer, however, are strongly dependent on the swirl rate. This is because the central methane jet is accelerated by the negative axial pressure gradient, of which strength is determined by the swirl. As a result, the inner vortex formation is significantly suppressed at a higher swirl rate. On the other hand, at a relatively low swirl rate, the inner vortices are shed continuously and the methane jet is pinched off. This particular mode promotes the mixing of methane and air, so that the flammable mixture can be formed at an earlier stage of the jet flow development. In addition, the evolution of secondary streamwise vortices is prompted by the combination of the periodic vortex ring shedding and the swirl. They also contribute to the mixing enhancement in the downstream region.