圆柱绕流尾迹中相干结构对湍流特性的影响

实验研究了圆柱尾迹中相干结构对湍流平均量的影响.用一个X热线在距离圆柱体y/d=4测量参考信号,用X热线和冷线相结合的三线探头测量从x/d=10～40的圆柱绕流尾迹中的速度和温度脉动,用条件平均的方法研充圆柱绕流近尾迹中的相干结构,并对相干结构对动量和热量的湍流输运过程的影响进行初步分析.结果表明:相干结构对尾迹中速度的横向脉动影响最大;相干结构对湍流参数的影响随x/d的变化明显.该工作还对相干结构引起的湍动能产生率变成负和逆温度梯度输运现象做了定性的解释.     

定常流作用于波浪形边界附近的圆柱上的水动力

本文介绍在水洞中做的作用于波浪形边界附近的水动力特性的实验研究。在基于圆柱模型直径的 R_c=10~4～1.9×10~4 范围内,测量了圆柱在波谷、波峰和不同距离上的阻力、升力脉动变化的频率。流谱显示实验揭示了尾流结构随距离-直径比 G/D 的变化情况及圆柱与边界相互作用的机制。     

隔离板对流向振荡圆柱尾流旋涡脱落的抑制

用隔离板对直径为D, 沿流向振荡的圆柱后涡脱落进行抑制. 隔离板放于圆柱尾流中心线上,控制参数包括隔离板长度L/D以及隔离板前缘到柱体振荡中心的距离G/D. 实验的雷诺数范围Re=V_∞D/v=1.01×10⁴～1.69×10⁴,柱体折减振频范围f_eD/V_∞=0～0.03, 柱体振幅固定为A/D=0.2. 风洞烟线显示和热线测量结果表明:当 G/D位于一个有效区域内时,可有效抑制振荡柱体尾流的旋涡脱落. 该有效区的大小随着隔离板板长的增大而增大, 随着Re数和圆柱振荡频率的增大而减小.     

狭缝射流撞击圆柱表面的湍流特性实验研究

实验研究了狭缝射流撞击圆柱表面后壁面射流区的平均流动和湍流特 性. 考察了雷诺数 Re (6000-20000), 喷口到受撞表面距 离 Y/W (5-13), 喷口宽度 W (6.25mm, 9.38mm), 受撞表 面曲率(半圆柱体直径 D = 150mm)对流动和湍流结构的影响. 通过分析 X 热线 在壁面射流区的测量结果发现，在近壁区域，表面曲率、 Re_{w} , Y/W 和 S/W 等 参数对 \sqrt {\overline{u^2}} / U_m 的影响比对 \sqrt {\overline{v^2}} / U_m 强，并且切 应力 \overline {uv} /U_m^2 对表面曲率变化最敏感. 当喷口与受撞击表面之间的距 离 Y/W 在一定范围内增加时， 沿圆柱表面流动的流向和横向的湍流强度增强. 用平板射流和圆柱体表面壁面射流的数据进行比较，从而得到表面曲率对壁面射流特 性的影响. 结果表明，曲率对壁面射流的影响较强， 并随着 S/W 的增大而增强. 随着雷诺数的增大，壁面曲率的影响也有强化的趋势.     

圆柱尾迹影响旁路转捩末期发卡涡涡包的研究

基于TR-PIV技术, 通过侧视和俯视两种情况对圆柱尾迹影响下旁路转 捩末期发卡涡涡包的结构及特征尺寸进行了实验研究. 结合二维空间子波变换和\lambda _{ci}准则, 运用线性随机估计方法对速度信号进行条件平均. 在侧视情况下, 条件平均结 果显示, 在边界层中一系列发卡涡涡头与壁面构成17^{\circ}的倾角, 并且被尾迹涡所占据的低 速区域出现在涡包上方的主流区中. 在俯视的结果中, 沿流向方向拉伸(流向尺度 3\delta, 展向尺度0.55δ)的低速条带结构出现在法向高 度为y/δ =0.2的流向-展向平面中, 并且在该低速条带的两侧对称地出 现了沿流向分布的反向旋转的涡结构. 可以得出: 在圆柱尾迹影响下旁路转捩的末期, 由于 尾迹涡诱导作用的影响, 发卡涡涡包在形态上显示出了更大尺度的特征.     

结构振动对湍流近尾迹的影响

研究了圆柱绕流中流体与结构的相互作用,侧重结构振动对湍流尾迹的影响,用激光测振仪测量圆柱在升力方向的位移;用热线和LDA（二维）测量湍流的近尾迹,通过变化自由流的速度和圆柱体直径（特征尺寸）来变化雷诺数,用两个振动特性不同的（一个相对刚性,一个相对弹）圆柱来产生尾迹,研究固体结构振动对湍流近尾迹的平均速度场和湍流场的影响,结果表明,结构自由振动对湍流近尾迹场影响明显,该影响随雷诺数的变化不明显。     

流向磁场作用下圆柱绕流的直接数值模拟

绕流是托卡马克装置中液态包层内常见的流动形态,对流场与热量分布有着重要的影响.本文通过直接数值模拟(DNS),研究了不同磁场强度下$Re=3900$的圆柱绕流,分析了磁场强度对于湍流尾迹的影响.无磁场情况下,直接数值模拟的结果与前人的实验及模拟结果吻合很好.圆柱下游的尾迹中,随着流向距离的增大, 流向速度剖面逐渐从U型进化呈V型, 并慢慢趋于平缓,这表明尾迹中的流动结构受圆柱影响逐渐减小.圆柱后方两侧的剪切层中,由于Kelvin-Helmholtz不稳定性的影响,可以清晰地看到小尺度剪切层涡的脱落.通过对无磁场的计算结果施加流向磁场,本文计算了哈特曼数($Ha$)分别为20, 40和80的工况,以研究磁场效应对于湍流的影响.结果表明磁场较弱时,流动依然呈三维湍流状态.随着磁场增强, 近圆柱尾流区受磁场抑制明显,回流区被拉长,剪切层失稳位置向下游转移.圆柱后方的涡结构由于受到竖直方向洛伦兹力的挤压作用,随着哈特曼数的增加尾迹区域逐渐变窄.相比于无磁场情况的涡结构,由于磁场的耗散作用,相应的涡结构尺度变小.该研究不仅扩展了现有磁场下湍流运动的参数范围,对于液态包层的设计及安全运行同样具有重要的理论指导意义和工程应用价值.      

DIRECT NUMERICAL SIMULATIONS ON THE TURBULENT FLOW PAST A CONFINED CIRCULAR CYLINDER WITH THE INFLUENCE OF THE STREAMWISE MAGNETIC FIELDS1)

绕流是托卡马克装置中液态包层内常见的流动形态,对流场与热量分布有着重要的影响.本文通过直接数值模拟(DNS),研究了不同磁场强度下$Re=3900$的圆柱绕流,分析了磁场强度对于湍流尾迹的影响.无磁场情况下,直接数值模拟的结果与前人的实验及模拟结果吻合很好.圆柱下游的尾迹中,随着流向距离的增大, 流向速度剖面逐渐从U型进化呈V型, 并慢慢趋于平缓,这表明尾迹中的流动结构受圆柱影响逐渐减小.圆柱后方两侧的剪切层中,由于Kelvin-Helmholtz不稳定性的影响,可以清晰地看到小尺度剪切层涡的脱落.通过对无磁场的计算结果施加流向磁场,本文计算了哈特曼数($Ha$)分别为20, 40和80的工况,以研究磁场效应对于湍流的影响.结果表明磁场较弱时,流动依然呈三维湍流状态.随着磁场增强, 近圆柱尾流区受磁场抑制明显,回流区被拉长,剪切层失稳位置向下游转移.圆柱后方的涡结构由于受到竖直方向洛伦兹力的挤压作用,随着哈特曼数的增加尾迹区域逐渐变窄.相比于无磁场情况的涡结构,由于磁场的耗散作用,相应的涡结构尺度变小.该研究不仅扩展了现有磁场下湍流运动的参数范围,对于液态包层的设计及安全运行同样具有重要的理论指导意义和工程应用价值.     

基于LBM-LES的三维串列多圆柱绕流尾涡结构研究

为了研究三维串联多圆柱绕流在不同圆柱间距比L/D(L为圆柱的间距,D为圆柱直径)情况下的尾涡结构,基于格子玻尔兹曼(LBM)方法,采用大涡模拟(LES)方法计算了高雷诺数下三维单圆柱及三维串联多圆柱在不同L/D情况下的尾涡流场.数值结果表明:LBM-LES方法可以较好地模拟出三维单圆柱以及串联多圆柱的尾涡结构;高雷诺数...     

高速碰撞下圆柱壳自由梁的穿孔特性

利用二级轻气炮加载,进行了球状2A12铝弹丸垂直撞击圆柱壳自由梁实验。并进行了弹丸速度、圆柱壳直径和壁厚等因素对穿孔直径影响的数值模拟,数值模拟结果和实验结果基本吻合。通过量纲分析和数值模拟结合,推导了穿孔直径与相关影响参数的经验关系式。研究结果表明:当圆柱壳直径和厚度不变时,高速撞击产生的穿孔直径在径向和轴向都随着弹丸速度增大而增大;当弹丸速度和圆柱壳厚度不变时,高速撞击产生的穿孔直径随着圆柱壳自由梁直径的增大而减小;当弹丸速度和圆柱壳直径不变时,穿孔直径随着圆柱壳厚度的增大而减小。     

DNS of a Laminar Separation Bubble in the Presence of Oscillating External Flow

A three-dimensional Direct Numerical Simulation (DNS) of a laminar separation bubble in the presence of oscillating flow is performed. The oscillating flow induces a streamwise pressure gradient varying in time. The special shape of the upper boundary of the computational domain, together with the oscillating pressure gradient causes the boundary layer flow to alternately separate and re-attach. When the inflow decelerates, the shear layer starts to separate and rolls up. Simultaneously the flow becomes 3D. After a transient period, the phase-averaged reverse flow inside the separation bubble reaches speeds ranging from 20 up to 150% of the free-stream velocity. During these phases, the flow is absolutely unstable and self-sustained turbulence can exist. When the inflow starts to accelerate, a spanwise roll of turbulent flow is shed from the shear layer. Shortly after this, the remainder of the separation bubble moves downstream and rejoins with the shed turbulent roll. During the flow-acceleration phase, a patch of laminar boundary layer flow is obtained. Along the flat plate, a series of turbulent patches of flow travelling downstream, separated by laminar flow can be observed, reminiscent of boundary layer flow in a turbine cascade with periodically appearing free-stream disturbances.     

Influence of bubble size on micro-bubble drag reduction

Micro-bubble drag reduction experiments were conducted in a turbulent water channel flow. Compressed nitrogen was used to force flow through a slot injector located in the plate beneath the boundary layer of the tunnel test section. Gas and bubbly mixtures were injected into a turbulent boundary layer (TBL), and the resulting friction drag was measured downstream of the injector. Injection into tap water, a surfactant solution (Triton X-100, 20 ppm), and a salt-water solution (35 ppt) yielded bubbles of average diameter 476, 322 and 254 μm, respectively. In addition, lipid stabilized gas bubbles (44 μm) were injected into the boundary layer. Thus, bubbles with d ⁺ values of 200 to 18 were injected. The results indicate that the measured drag reduction by micro-bubbles in a TBL is related strongly to the injected gas volumetric flow rate and the static pressure in the boundary layer, but is essentially independent of the size of the micro-bubbles over the size range tested.     

用三分量LDV测量壁湍流边界层

本文用三分量激光多普勒测速仪对壁湍流边界层进行了测量，经过数据处理后，给出了各湍流参数的分布曲线，并与Ｄｊｅｎｉｄｉ等的测量结果进行了比较，吻合得很好。结果表明，用３Ｄ－ＬＤＶ研究湍流是可行的     

Large-eddy simulation study of upstream boundary conditions influence upon a backward-facing step flowÉtude par simulation des grandes échelles de l'influence des conditions aux limites amont sur l'écoulement derrière une marche descendante

We use Large Eddy Simulation to investigate the influence of upstream boundary conditions on the development of a backward facing step flow. The first inlet condition consists of a mean turbulent boundary layer velocity profile perturbed by a white noise. The second relies upon a precursor calculation where the development of a quasi-temporal turbulent boundary layer is simulated. In this case, the quasi-longitudinal vortices in the upstream turbulent boundary-layer trigger the destabilization of the shear layer just behind the step, resulting in a shortening of the recirculation length and an increase of the characteristic frequency associated to the Kelvin–Helmholtz vortices. The mean flow and the characteristic frequencies of pressure fluctuations are strongly dependent of the upstream flow. It demonstrates the importance of realistic boundary conditions for the simulation of complex 3D flows or for flow control simulations. To cite this article: J.-L. Aider, A. Danet, C. R. Mecanique 334 (2006).     

Vortex reynolds number in turbulent boundary layers

The effects of vortex Reynolds number on the statistics of turbulence in a turbulent boundary layer have been investigated. Vortex Reynolds number is defined as the ratio of circulation around the vortex structure to the fluid viscosity. The vortex structure of the outer region was modeled and a full numerical simulation was then conducted using a high-order spectral method. A unit domain of the outer region of a turbulent boundary layer was assumed to be composed of essentially three elements: a wall, a Blasius mean shear, and an elliptic vortex inclined at 45° to the flow direction. The laminar base-flow Reynolds number is roughly in the same range as that of a turbulent boundary layer based on eddy viscosity, and the vortex-core diameter based on the boundary-layer thickness is nearly the same as the maximum mixing length in a turbulent boundary layer. The computational box size, namely, 500, 150, and 250 wall units in the streamwise, surface-normal, and spanwise directions, respectively, is approximately the same as the measured quasi-periodic spacings of the near-wall turbulence-producing events in a turbulent boundary layer. The effects of vortex Reynolds number and the signs of the circulation on the moments of turbulence were examined. The signs mimic the ejection and sweep types of organized motions of a turbulent boundary layer. A vortex Reynolds number of 200 describes the turbulence moments in the outer layer reasonably well.     

Response of a skewed turbulent boundary layer to favourable pressure gradient

Experimental results are reported for the response to a favourable pressure gradient of an initially turbulent boundary layer (Re _θ?≈?1600) developing on a flat plate with its leading edge skewed at 60° to the approach flow. The pressure gradient orthogonal to the leading edge is nominally the same as that which was shown by Escudier et?al. [(1998) Exp Fluids 25: 491–502] to cause extreme thinning of a two-dimensional (2D) (i.e. unskewed) turbulent boundary layer and the intermittency in the immediate vicinity of the surface to fall to zero, i.e. an apparent laminarisation of the boundary layer. In the case of the skewed boundary layer, the responses of the turbulence and mean-flow structures are qualitatively similar to those for the 2D situation. However, the streamwise pressure gradient is much weaker than for the 2D experiment and the extent of the changes it produces is much reduced. Even so, the changes are considerably greater than would be expected from the magnitude of the streamwise pressure gradient.     

DDT in a smooth tube filled with a hydrogen–oxygen mixture

Results of experimental study on DDT in a smooth tube filled with sensitive mixtures having detonation cell size from 1 to 3 orders of magnitude smaller than the tube diameter are presented. Stoichiometric hydrogen–oxygen mixtures were used in the tests with initial pressure ranging from 0.2 to 8 bar. A dependence of the run-up distance to DDT on the initial pressure is studied. This dependence is found to be close to the inverse proportionality. It is suggested that the flow ahead of the flame results in formation of the turbulent boundary layer. This boundary layer controls the scale of turbulent motions in the flow. A simple model to estimate the maximum scale of the turbulent pulsations (boundary layer thickness) at flame positions along the tube is presented. The largest scale of the turbulent motions at the location of the onset of detonation is shown to be 1 order of magnitude greater than the detonation cell widths, λ, in all the tests. It is suggested that the onset of detonation is triggered during flame acceleration as soon as the maximum scale of the turbulent pulsations increases up to about 10 λ. The model to estimate the maximum size of turbulent motions, δ, and the correlation δ≈ 10λ, give a basis for estimations of the run-up distances to DDT in tubes with internal diameter D > 20λ. PACS 47.40.-x; 47.27.Nz This paper was based on work that was presented at the 19th Inter-national Colloquium on the Dynamics of Explosions and Reactive Systems, Hakone, Japan, July 27 - August 1, 2003     

The turbulent boundary layer on the moving surface of a cylindrical body

A study is made of the problem of a two-dimensional turbulent boundary layer on the moving surface of a cylindrical body (a Rankine oval with a relative elongation of four) moving at constant velocity in an incompressible fluid. For the numerical simulation of the turbulent flow of the fluid, the boundary layer is divided into exterior and interior regions in accordance with a two-layer model, using different expressions for the coefficients of turbulent transfer for each region. A study was nade of the development of the boundary layer on the body at different speeds of the body surface and different Reynolds numbers. The following integral characteristics were found by numerical calculation: the work of friction as the body is displaced; the work expended on the movement of its surface; and, for a flow regime with separation, the work of the pressure force. In this case the following model of separation flow is assumed: beyond the singular point in the solution of the boundary layer equations that indicates the appearance of a region of reverse flow, the pressure and friction stress on the wall are constant and are determined by their values at the singular point.Translated from Izvestiya Akademii Nauk SSSH, Mekhanika Zhidkosti i Gaza, No. 5, pp. 61–67, September–October, 1984.Finally, the author would like to thank G. G. Chernyi and Yu. D. Shevelev for useful discussions and for their interest in this work.     

Numerical prediction of turbulent boundary layer noise from a sharp-edged flat plate

An efficient hybrid uncorrelated wall plane waves–boundary element method (UWPW-BEM) technique is proposed to predict the flow-induced noise from a structure in low Mach number turbulent flow. Reynolds-averaged Navier-Stokes equations are used to estimate the turbulent boundary layer parameters such as convective velocity, boundary layer thickness, and wall shear stress over the surface of the structure. The spectrum of the wall pressure fluctuations is evaluated from the turbulent boundary layer parameters and by using semi-empirical models from literature. The wall pressure field underneath the turbulent boundary layer is synthesized by realizations of uncorrelated wall plane waves (UWPW). An acoustic BEM solver is then employed to compute the acoustic pressure scattered by the structure from the synthesized wall pressure field. Finally, the acoustic response of the structure in turbulent flow is obtained as an ensemble average of the acoustic pressures due to all realizations of uncorrelated plane waves. To demonstrate the hybrid UWPW-BEM approach, the self-noise generated by a flat plate in turbulent flow with Reynolds number based on chord Re_c = 4.9 × 10⁵ is predicted. The results are compared with those obtained from a large eddy simulation (LES)-BEM technique as well as with experimental data from literature.