超声速湍流边界层湍流统计与湍流结构分析

可压缩边界层转捩问题与湍流问题一直是制约高超声速飞行器发展的关键基础问题,也是近年来流体力学领域研究的热点问题.采用直接数值模拟方法,获得了空间发展的Ma=2.25超声速湍流边界层流场,通过对湍流边界层的发展状态进行评估,得出有效的Reynolds数Re_θ范围约为2 600~4 600.对壁面摩阻系数开展了分解,获得了各分量的占比,对充分发展的湍流边界层进行1阶和高阶统计分析,包括形状因子、壁面律、平坦因子与偏斜因子、Reynolds应力、脉动涡量等,得到了剪切Reynolds数与动量Reynolds数之间的关系式,分析了湍流边界层壁面律的分层特性,发现湍流的间歇特性主要分布在y+ < 30的区域,并且法向速度脉动的间歇性远高于另外两者,3个方向上的Reynolds应力分布和涡量分布都存在较大差异.通过两点相关性分析和Lagrange涡结构,对近壁区湍流结构进行了分析,包括流向平面和展向平面,发现流向脉动速度的相关区域流向尺度较长,呈现狭长的特性,并且流向平面的相关系数与壁面存在一定的夹角;而在边界层外层,流向速度脉动相关区域的流向尺度变短而展向尺度增加,呈现宽胖型.研究结果进一步加深了对超声速湍流边界层的认识,为下一步湍流边界层的流动控制奠定了基础.      

Incomplete similarity of a plane turbulent wall jet on smooth and transitionally rough surfaces

This study assesses the hypothesis of incomplete similarity for a plane turbulent wall jet on smooth and transitionally rough surfaces. Typically, a wall jet is considered to consist of two regions: an inner layer and an outer layer. The degree to which these two regions reach equilibrium with each other and interact to produce the property of self-similarity remains an open question. In this study, the analysis of the outer and inner regions indicates that each region is characterised by a half-width which exhibits its own distinct dependence on the streamwise distance x from the slot, and a single self-similar structure for both regions does not exist. More specifically, the inner and outer layers of the wall jet exhibit different scaling laws, which results in two self-similar mean velocity profiles, both of which retain a dependence on the slot height H. As such, incomplete similarity of the wall jet on smooth and transitionally rough surfaces is confirmed by this study. In addition, comparison of the experimental results for the transitionally rough surface with the smooth wall case indicates that the surface roughness modifies the development of the mean velocity profile in both the inner and outer regions, although the effect on the outer region is relatively small and close to the experimental uncertainty.     

Investigation of scaling laws in a turbulent boundary layer flow with adverse pressure gradient using PIV

We present an experimental investigation and data analysis of a turbulent boundary layer flow at a significant adverse pressure gradient at Reynolds number up to Re_θ = 10, 000. We combine large-scale particle image velocimetry (PIV) with microscopic PIV for measuring the near wall region including the viscous sublayer. We investigate scaling laws for the mean velocity and for the total shear stress in the inner part of the boundary layer. In the inner part the mean velocity can be fitted by a log-law. In the outer part of the inner layer the log-law ceases to be valid. Instead, a modified log-law provides a good fit, which is given in terms of the pressure gradient parameter and a parameter for the mean inertial effects. Finally we describe and assess a simple quantitative model for the total shear stress distribution which is local in wall-normal direction without streamwise history effects.     

出口雷诺数对平面射流自保持性的影响

通过实验研究出口雷诺数对平面湍流射流自保持性的影响. 测量的射流来自相同的喷嘴但不同的雷诺数Re（≡U_jh/ν,其中U_j是出口平均速度、h是窄缝出口的厚度和ν是黏性系数）,其变化范围是Re=4582—57735.所得的数据包括沿轴线的平均速度、湍流强度、积分尺度、高阶矩和能谱. 实验发现,随着Re的增大,平面射流发展减慢,平均速度和湍流强 关键词： 平面射流 雷诺数 自保持性     

LES of spatially developing turbulent boundary layer over a concave surface

We revisit the problem of a spatially developing turbulent boundary layer over a concave surface. Unlike previous investigations, we simulate the combined effects of streamline curvature as well as curvature-induced pressure gradients on the turbulence. Our focus is on investigating the response of the turbulent boundary layer to the sudden onset of curvature and the destabilising influence of concave surface in the presence of pressure gradients. This is of interest for evaluating the turbulence closure models. At the beginning of the curve, the momentum thickness Reynolds number is 1520 and the ratio of the boundary layer thickness to the radius of curvature is δ₀/R = 0.055. The radial profiles of the mean velocity and turbulence statistics at different locations along the concave surface are presented. Our recently proposed curvature-corrected Reynolds Averaged Navier-Stokes (RANS) model is assessed in an a posteriori sense and the improvements obtained over the base model are reported. From the large Eddy simulation (LES) results, it was found that the maximum influence of concave curvature is on the wall-normal component of the Reynolds stress. The budgets of wall-normal Reynolds stress also confirmed this observation. At the onset of curvature, the effect of adverse pressure gradient is found to be predominant. This decreases the skin friction levels below that in the flat section.     

Large eddy simulation of turbulent horizontal buoyant jets

Large eddy simulations (LESs) of turbulent horizontal buoyant jets are carried out using a high-order numerical method and Sigma subgrid-scale (SGS) eddy-viscosity model, for a number of different Reynolds (Re) and Richardson (Ri) numbers. Simulations at previous experimental flow conditions (Re = 3200, 24, 000 and Ri = 0, 0.01) are carried out first, and the results are found to be qualitatively and quantitatively similar to the experimental results, thus validating the numerical methodology. The effect of varying Ri (values 2×10^?4, 0.001, 0.005, and 0.01) and Re (3200 and 24, 000) is studied next. The presence of stable stratification on one side and unstable stratification on the other side of the jet centreline leads to an asymmetric development of horizontal buoyant jets. It is found that this asymmetry, the total radial spread and the vertical deflection are significantly affected by Ri, while Re affects only the radial asymmetry. The need for developing improved integral models, accounting for this asymmetry, is pointed out. Turbulent production and dissipation rates are investigated, and are found to be symmetric in the horizontal plane, but asymmetric in the mid-vertical plane. A previously proposed model, for correlation between the vertical component of the fluctuating scalar flux vector and the vertical cross-correlation component of the Reynolds tensor, is modified based on the current LES results. Instantaneous scalar and velocity fields are analysed to reveal the structure of horizontal buoyant jets. Similar to the developed turbulent jet, the flow close to the nozzle too is found to be markedly different in the stable and unstable stratification regions. Persistent coherent vortex rings are found in the stable stratification region, while intermittent breakdown of vortex rings into small-scale structures is observed in the unstable stratification region. Similarities and differences between the flow structures in the horizontal buoyant jet configuration and those in the jet in crossflow configuration are discussed. Finally, a dynamic mode decomposition analysis is carried out, which indicates that the flow in the unstable stratification region is more energetic and prone to instabilities, as compared to the flow in the stable stratification region.     

Surface roughness effects on the turbulence statistics in a low Reynolds number channel flow

A high-resolution particle image velocimetry was used to characterize a low Reynolds number turbulent flow in a channel. Experiments were conducted over a sand grain-coated surface of large relative roughness, and the results were compared with measurements over a smooth surface. The roughness perturbation significantly modified the outer layer. Even though the streamwise Reynolds stress shows less sensitivity in the outer layer to the boundary condition, significant enhancements were observed in the wall-normal Reynolds stress and the Reynolds shear stress. These modifications were considered as footprints of the larger-scale eddies transporting intense wall-normal motions away from the rough wall. A quadrant decomposition shows that strong and more frequent ejections are responsible for the larger values of the mean Reynolds shear stress over the rough wall. The results also indicate that spanwise vortex cores with mean vorticity of the same sign as the mean shear are the dominant smaller-scale vortical structures over the smooth and rough walls. A linear stochastic estimation-based analysis shows that the average larger-scale structure associated with these vortices is a shear layer that strongly connects the outer layer flow to the near-wall flow. A proper orthogonal decomposition of the flow suggests that the large-scale eddy is more energetic for the rough wall, and contributes more significantly to the resolved turbulent kinetic energy and the Reynolds shear stress than the smooth wall.     

Turbulent plane wall jets over smooth and rough surfaces

Large-eddy simulations were carried out to study the effects of surface roughness on a plane wall-jet using the Lagrangian dynamic eddy-viscosity subgrid-scale model, at Re = 7500 (based on the jet bulk velocity and height). Results over both smooth and rough surfaces were validated by experimental data at the same Reynolds number. As the jet is injected into the still environment, large-scale rollers are generated in the shear layer between the high-momentum fluid of the jet and the surrounding and are convected downstream with the flow. To understand the extent to which the outer-layer structures modify the flow in the inner layer and the extent to which the effect of roughness spreads away from the wall, both instantaneous and mean flow fields were investigated. The results revealed that, for the Reynolds number and roughness height considered in this study, the effect of roughness is mostly confined to the near-wall region of the wall jet. There is no structural difference between the outer layer of the wall jet over the smooth and rough surfaces. Roughness does not affect the size of the outer-layer structures or the scaling of the profiles of Reynolds stresses in the outer layer. However, in the inner layer, roughness redistributes stresses from streamwise to wall-normal and spanwise directions toward isotropy. Contours of joint probability-density function of the streamwise and wall-normal velocity fluctuations at the bottom of the logarithmic region match those of the turbulent boundary layer at the same height; while the traces of the outer-layer structure were detected at the top of the logarithmic region, indicating that they do not affect the flow very close to the wall, but still modify a major portion of the inner layer. This modification must be taken into consideration when the inner layer of a wall jet is compared with the conventional turbulent boundary layer.     

Large-scale flow structures in a planar offset attaching jet with a co-flowing wall jet

Conditional averaging techniques were used to examine the periodic motions that were observed in flows consisting of an offset planar jet and a co-flowing planar wall jet. The offset jet is one jet height (H_j) away from the wall and has a Reynolds number of approximately 40, 000, based on H_j and flow-rate averaged velocity U_o; for the co-flowing jets, different heights (i.e., 0.18H_j and 0.5H_j) and velocities (i.e., 0.56U_o and 0.36U_o) were considered. The flows had periodic motions with frequencies fH_j/U_o = 0.28 and 0.49 or fH_c/U_o = 0.23 and 0.25, where H_c is the distance between the jets. The periodic motions were present in both the inner shear layer of the offset jet above the re-circulation region and the outer shear layer of the wall jet below the re-circulation region. The motions from the inner shear layer of the offset jet persisted in the shear layer that formed downstream of the re-circulating region. There were periodic motions in the outer shear layer of the offset jet particularly in the flow with the smaller wall jet. The relative contribution of the motions to the total fluctuations increased as the flow evolved downstream reaching a maximum approximately 4H_c downstream of the flow exit. The relative contribution of the periodic motions to the turbulent fluctuations was similar in the two flows but the periodic motions had a much larger impact on the near-wall velocity and pressure fluctuations in the flow with the smaller wall jet due to the trajectory of the periodic structures, the distance of the periodic structures to the wall and the size of these structures.     

Large-eddy simulation of impinging jets on smooth and rough surfaces

We performed large-eddy simulations (LES) of forced impinging jets over smooth and rough surfaces, containing large-scale, azimuthal vortices generated by the enhanced primary instability in the jet shear layer. The interaction between these vortices and the turbulence in the wall jet that is formed downstream of the impingement region determines their rate of decay. To explore the surface-roughness effects on the evolution of the vortices, sand-grain-like surfaces are generated using uniformly distributed but randomly oriented ellipsoids. The flow is compared to our previous LES of jets impinging on a smooth surface. In spite of the severe modification caused by the roughness on the near-wall flow, the vortex development is not significantly altered. Slightly faster decay of the primary vortices is observed in the rough-wall case compared to the smooth-wall one; the secondary vortex that detaches from the wall and is lifted up has larger vorticity. The highly disturbed near-wall flow is advected outward and affects the evolution of the primary vortex for a longer period during the vortex interaction. The robust turbulent generation mechanism in the outer shear layer, however, mitigates the changes in vortex behaviour. The momentum deficit and the enhancement of turbulence due to the surface roughness play a key role during this process.     

Level-crossing statistics and production in low Reynolds number wall turbulence

Statistical properties at velocity level crossings are analysed in a low Reynolds number fully developed turbulent channel flow. Emphasis is placed on local production statistics conditioned by fixed amplitudes of the streamwise u and wall-normal v velocity components. Direct numerical simulations performed in large computational domains are used for this purpose. The Reynolds number based on the channel half width and the shear velocity varies from 180 to 1100. Particular attention is paid to correctly determine the conditional quantities at level crossings of u and v velocity fluctuations, to prevent biasing effects. Level crossings along the longitudinal x and spanwise z directions in homogeneous planes are introduced together with different characterisations, such as directional and contour crossings. The frequency of events detected at fixed thresholds follows a Gaussian model acceptably well. There is more level-crossing activity in the spanwise direction than in the streamwise. The mean conditional production distributions are dissymmetrical with respect to the level-crossing thresholds in the low buffer and viscous sublayers. These statistics differ significantly from a statistical model that assumes joint normality between u and v. There is a clear but relatively mild Reynolds number dependence of the conditional expected means of u and v that do not scale with inner, outer or mixed variables. The Reynolds number sensitivity of these statistics increases towards the edge of the wall layer.     

Effects of upstream roughness and Reynolds number on separated and reattached turbulent flow

This paper presents experimental investigation of upstream roughness and Reynolds number effects on the recirculation region over a smooth forward facing step. The upstream rough wall was produced from 1.5 mm sand grains and the Reynolds number based on step height, Re_h, was varied from 2040 to 9130 for both the upstream smooth and rough walls. For the smooth wall, the reattachment length increased monotonically with Re_h to an asymptotic value of 2.2 step heights for Re_h ≥ 6380. Upstream roughness reduced the reattachment length by 44% because of larger momentum deficit and higher turbulence level in the rough wall boundary layer. The mean velocities and Reynolds stresses were also reduced by roughness. The Reynolds shear stress and production of turbulent kinetic energy showed high negative values at the leading edge of the step indicating counter-gradient diffusion. The implications of these results for standard eddy viscosity models are discussed.     

Numerical investigation of the interaction of the turbulent dual-jet and acoustic propagation

In order to study the interaction between two independent jets, a three-dimensional(3D) transient mathematical model is developed to investigate the flow field and acoustic properties of the two-stream jets. The results are compared with those of the single-stream jet at Mach number 0.9 and Reynolds number 3600. The large eddy simulation(LES) with dynamic Smagorinsky sub-grid scale(SGS) approach is used to simulate the turbulent jet flow structure. The acoustic field is evaluated by the Ffowcs Williams–Hawkings(FW-H) integral equation. Considering the compressibility of high-speed gas jets, the density-based explicit formulation is adopted to solve the governing equations. Meanwhile, the viscosity is approximated by using the Sutherland kinetic theory. The predicted flow characteristics as well as the acoustic properties show that they are in good agreement with the existing experimental and numerical results under the same flow conditions available in the literature. The results indicate that the merging phenomenon of the dual-jet is triggered by the deflection mechanism of the Coanda effect, which sequentially introduces additional complexity and instability of flow structure. One of the main factors affecting the dual-jet merging is the aperture ratio, which has a direct influence on the potential core and surrounding flow fluctuation. The analysis on the noise pollution reveals that the potential core plays a fundamental role in noise emission while the additional mixing noise makes less contribution than the single jet noise. The overall sound pressure level(OASPL) profiles have a directive property, suggesting an approximate 25° deflection from the streamwise direction, however, shifting toward lateral direction of about 10° to 15° in the dual-jet. The conclusion obtained in this study can provide valuable data to guide the development of manufacturing-green technology in the multi-jet applications.     

Recycling inflow method for simulations of spatially evolving turbulent boundary layers over rough surfaces

The technique by Lund et al. to generate turbulent inflow for simulations of developing boundary layers over smooth flat plates is extended to the case of surfaces with roughness elements. In the Lund et al. method, turbulent velocities on a sampling plane are rescaled and recycled back to the inlet as inflow boundary condition. To rescale mean and fluctuating velocities, appropriate length scales need be identified and for smooth surfaces, the viscous scale l_ν = ν/u_τ (where ν is the kinematic viscosity and u_τ is the friction velocity) is employed for the inner layer. Different from smooth surfaces, in rough wall boundary layers the length scale of the inner layer, i.e. the roughness sub-layer scale l_d, must be determined by the geometric details of the surface roughness elements and the flow around them. In the proposed approach, it is determined by diagnosing dispersive stresses that quantify the spatial inhomogeneity caused by the roughness elements in the flow. The scale l_d is used for rescaling in the inner layer, and the boundary layer thickness δ is used in the outer region. Both parts are then combined for recycling using a blending function. Unlike the blending function proposed by Lund et al. which transitions from the inner layer to the outer layer at approximately 0.2δ, here the location of blending is shifted upwards to enable simulations of very rough surfaces in which the roughness length may exceed the height of 0.2δ assumed in the traditional method. The extended rescaling–recycling method is tested in large eddy simulation of flow over surfaces with various types of roughness element shapes.     

Research on the jet characteristics of the deflector–jet mechanism of the servo valve

In view of the complicated structure of the deflector-jet mechanism,a mathematical model based on the turbulent jet flow theory in the deflector-jet amplifier is proposed.Considering the energy transformation and momentum variation,an equation of the flow velocity distribution at the key fluid region is established to describe the morphological changes of the fluid when it passes through the deflector and jets into the receiver.Moreover,the process is segmented into four stages.According to the research results,the oil enters the deflector and impinges with the side wall.Then one part of the oil's flow velocity decreases and a high pressure zone is formed by the oil accumulation,the other part of the oil reverses out of the deflector along the side wall.Prior to entering the receiver,the flow is a kind of plane impinging jet.Virtually,the working pressure of the receiver is generated by the impact force,while the high speed fluid flows out of the receiver and forms a violent vortex,which generates negative pressure and causes the oil to be gasified.Compared with the numerical simulation results,the turbulent jet model that can effectively describe the characteristics of the deflector-jet mechanism is accurate.In addition,the calculation results of the prestage pressure characteristic have been verified by experiments.     

同轴流动聚焦中射流不稳定性的理论研究

壳核结构的微胶囊在医学药学材料食品农业等领域具有广泛的应用前景, 其制备方法一直是相关领域关注的焦点.同轴流动聚焦(co-flow focusing)是一种新型制备技术, 利用复合射流的破碎制备微胶囊具有包裹率高过程量化可控参数域广产率高等诸多优势.在实验中, 复合射流的破碎受到多个过程参数的影响, 并涉及了多层界面的耦合效应.利用简化的物理模型, 在时间和时空域中分析了三相水-油-水复合射流不稳定性的发展和演化.在黏性流体线性稳定性理论中, 同轴射流和驱动液体的基本速度型分别基于管流和误差函数构造, 并通过数值方法求解满足相应边界条件下的线化小扰动控制方程.结果表明:增加内外层界面的界面张力均有利于射流的破碎; 流体的黏性对同轴射流的稳定性均有着促进作用; 越大的黏性越小的内界面张力对应着越大的射流破碎波长; 内外界面的耦合作用以及复合液滴的包裹情况均与内外射流的半径比息息相关; 绝对-对流不稳定性转换的临界Weber数随Reynolds数内层界面张力的增大而增大, 随内层和驱动流体的黏性增大而减小.这些结果将有助于提高液体驱动下同轴流动聚焦技术的过程控制, 为实际应用提供理论指导.      

Flow Characteristics of Flapping Motion of a Plane Water Jet Impinging onto Free Surface

A submerged turbulent plane jet in shallow water impinging vertically onto the free surface will produce a large-scale flapping motion when the jet exit velocity is larger than a critical one. The flapping phenomenon is verified in this paper through a large eddy simulation where the free surface is modeled by volume of fluid approach. The quantitative results for flapping jet are found to be in good agreement with available experimental data in terms of mean velocity, flapping-induced velocity and turbulence intensity. Results show that the flapping motion is a new flow pattern with characteristic flapping frequency for submerged turbulent plane jets, the mean centerline velocity decay is considerably faster than that of the stable impinging jet without flapping motion, and the flapping-induced velocities are as important as the turbulent fluctuations.     

Numerical analysis of shock wave and supersonic turbulent boundary interaction between adiabatic and cold walls

Direct numerical simulations of shock wave and supersonic turbulent boundary layer interaction in a 24° compression ramp with adiabatic and cold-wall temperatures are conducted. The wall temperature effects on turbulence structures and shock motions are investigated. The results are validated against previous experimental and numerical data. The effects of wall cooling on boundary layer characteristics are analysed. Statistical data show that wall cooling has a significant effect on the logarithmic region of mean velocity profile downstream the interaction region. Moreover, the influence of wall temperature on Reynolds stress anisotropy is mainly limited in the near-wall region and has little change on the outer layer. As the wall temperature decreases, the streamwise coherency of streaks increases. Based on the analysis of instantaneous Lamb vector divergence, the momentum transport between small-scale vortices on cold-wall condition is significantly enhanced. In addition, spectral analysis of wall pressure signals indicates that the location of peak of low-frequency energy shifts toward higher frequencies in cold case. Furthermore, the dynamic mode decomposition results reveal two characteristic modes, namely a low-frequency mode exhibiting the breathing motion of separation bubble and a high-frequency mode associated with the propagation of instability waves above separation bubble. The shape of dynamic modes is not sensitive to wall temperature.     

A mixed-time-scale low-Reynolds-number one-equation turbulence model

In this paper, a new low-Reynolds-number (LRN) one-equation turbulence model for eddy viscosity is proposed. A mixed time scale, representing a combination of three time scales: two time scales made of strain-rate parameter S and vorticity parameter Ω and the turbulent time scale k/?, is introduced into this model. The proposed model is derived from an LRN k?? two-equation model where the mixed time scale has been proved to be very effective for predicting local flows over complex terrains. In the transport equation of the model, the mixed time scale is included in the production and the dissipation terms. The new model is evaluated in channel flows at various Reynolds numbers, boundary layer flows with or without pressure gradient and backward-facing step flows with different expansion ratios and Reynolds numbers. Then the grid convergence of the model is investigated. Finally, the model performance for different values of the weighting constant C_s in the mixed time scale is assessed. The results show that the proposed model reproduces the correct wall-limiting behaviour of turbulent quantities and performs well in the near-wall region of turbulent flows. The model could be expected to be adopted in hybrid Reynolds averaged Navier–Stokes/large eddy simulation methodology for complex wall-bounded flows at high Reynolds numbers.     

Large-eddy simulation of impinging jets with embedded azimuthal vortices

We have carried out large-eddy simulations of an impinging jet with embedded azimuthal vortices, a model of the wake of a helicopter hovering in ground effect. The azimuthal vortices are generated by sinusoidal forcing of the velocity at the jet exit. They strengthen while they are advected towards the ground; when they are close to the solid surface, a layer of opposite-sign vorticity is formed at the wall, and lifted up to form a secondary vortex that interacts with the primary one. Regions of reversed flow are caused by the strong, localised, adverse pressure gradient. After this interaction, the primary vortices begin to decay, mostly due to the Reynolds shear stresses, which contribute to the turbulent diffusion of vorticity term in the budget of the phase-averaged azimuthal vorticity. This mechanism is extremely robust, and plays the most important role in the vortex decay even if no turbulence is initially present in the jet, or if the no-slip condition is removed. A three-dimensional instability also plays a role: removing it leads to slower decay. Our results also point out some challenges for turbulence models for the unsteady Reynolds-averaged Navier–Stokes equations.