Investigation of the unsteadiness of a shock-reflection interaction with time-resolved particle image velocimetry

The spatio-temporal dynamics of an impinging shock/boundary layer interaction at Mach 2 and under incipient separation conditions, has been investigated experimentally by means of high-speed particle image velocimetry (PIV). The available PIV acquisition rate of up to 20 kHz permits a time-resolved characterization of the interaction. The dynamics of different flow regions—notably the separation region and the reflected shock—were quantified by means of temporal auto-correlation fields and pseudo-spectral analysis. The PIV data further enable to investigate the relationship between spatially extended flow features, such as shock position and bubble size, as well as the influence of the upstream boundary layer. The results confirm earlier studies that there is an important upstream effect on the present incipient interaction.     

Stability and modal analysis of shock/boundary layer interactions

The dynamics of oblique shock wave/turbulent boundary layer interactions is analyzed by mining a large-eddy simulation (LES) database for various strengths of the incoming shock. The flow dynamics is first analyzed by means of dynamic mode decomposition (DMD), which highlights the simultaneous occurrence of two types of flow modes, namely a low-frequency type associated with breathing motion of the separation bubble, accompanied by flapping motion of the reflected shock, and a high-frequency type associated with the propagation of instability waves past the interaction zone. Global linear stability analysis performed on the mean LES flow fields yields a single unstable zero-frequency mode, plus a variety of marginally stable low-frequency modes whose stability margin decreases with the strength of the interaction. The least stable linear modes are grouped into two classes, one of which bears striking resemblance to the breathing mode recovered from DMD and another class associated with revolving motion within the separation bubble. The results of the modal and linear stability analysis support the notion that low-frequency dynamics is intrinsic to the interaction zone, but some continuous forcing from the upstream boundary layer may be required to keep the system near a limit cycle. This can be modeled as a weakly damped oscillator with forcing, as in the early empirical model by Plotkin (AIAA J 13:1036–1040, 1975).     

Effects of sub-millimeter-bubble injection on transition to turbulence in natural convection boundary layer along a vertical plate in water

Temperature and velocity measurements are performed to clarify the effects of sub-millimeter-bubble injection on the transition to turbulence in the natural convection boundary layer along a vertical plate in water. In particular, we focus on the relationship between the bubble injection position L and the transition to turbulence in the natural convection boundary layer. The bubble injection positions used in our experiments are L = 1.6 and 3.6 mm. Bubble injection at L = 1.6 mm delays the transition to turbulence in the natural convection boundary layer, while that at L = 3.6 mm accelerates the transition to turbulence in the boundary layer. In the case of L = 1.6 mm, the appearance region of the liquid velocity fluctuation in the bubble-induced upward flow in the upstream unheated section is restricted to near the wall, although the peak of the liquid velocity fluctuation is high. In contrast, in the case of L = 3.6 mm, the relatively large liquid velocity fluctuation is distributed widely over the laminar boundary layer width. These results suggest that the effect of the liquid velocity fluctuation on the laminar boundary layer is quite different between L = 1.6 and 3.6 mm. It is therefore expected that the transition to turbulence in the natural convection boundary layer for the case with bubble injection is dependent on the magnitude and appearance region of the liquid velocity fluctuation in the bubble-induced upward flow in the upstream unheated section.     

激波与转捩边界层干扰非定常特性数值分析

激波与边界层干扰的非定常问题是高速飞行器气动设计中基础研究内容之一.以往研究主要针对层流和湍流干扰,在分离激波低频振荡及其内在机理方面存在着上游机制和下游机制两类截然不同的理论解释.分析激波与转捩边界层干扰下非定常运动现象有助于进一步加深理解边界层状态以及分离泡结构对低频振荡特性的影响规律,为揭示其产生机理指出新的方向.采用直接数值模拟方法对来流马赫数2.9,24?压缩拐角内激波与转捩边界层干扰下激波的非定常运动特性进行了数值分析.通过在拐角上游平板特定的流向位置添加吹吸扰动激发流动转捩,使得进入拐角的边界层处于转捩初期阶段.在验证了计算程序可靠性的基础上,详细分析了转捩干扰下激波运动的间歇性和振荡特征,着重研究了分离泡展向三维结构对激波振荡特性的影响规律,最后还初步探索了转捩干扰下激波低频振荡产生的物理机制.研究结果表明:分离激波的非定常运动仍存在强间歇性和低频振荡特征,其时间尺度约为上游无干扰区内脉动信号特征尺度的10倍量级;分离泡展向三维结构不会对分离激波的低频振荡特征产生实质影响.依据瞬态脉动流场的低通滤波结果,转捩干扰下激波低频振荡的诱因来源于拐角干扰区下游,与流场中分离泡的收缩/膨胀运动存在一定的关联.     

Hypersonic interference heating in the vicinity of surface protuberances

The understanding of the behaviour of the flow around surface protuberances in hypersonic vehicles is developed and an engineering approach to predict the location and magnitude of the highest heat transfer rates in their vicinity is presented. To this end, an experimental investigation was performed in a hypersonic facility at freestream Mach numbers of 8.2 and 12.3 and Reynolds numbers ranging from Re _∞/m = 3.35 × 10⁶ to Re _∞/m = 9.35 × 10⁶. The effects of protuberance geometry, boundary layer state, freestream Reynolds number and freestream Mach numbers were assessed based on thin-film heat transfer measurements. Further understanding of the flowfield was obtained through oil-dot visualizations and high-speed schlieren videos. The local interference interaction was shown to be strongly 3-D and to be dominated by the incipient separation angle induced by the protuberance. In interactions in which the incoming boundary layer remains unseparated upstream of the protuberance, the highest heating occurs adjacent to the device. In interactions in which the incoming boundary layer is fully separated ahead of the protuberance, the highest heating generally occurs on the surface just upstream of it except for low-deflection protuberances under low Reynolds freestream flow conditions in which case the heat flux to the side is greater.     

Time-resolved stereo PIV measurements of shock–boundary layer interaction on a supercritical airfoil

Time-resolved stereo particle-image velocimetry (TR-SPIV) and unsteady pressure measurements are used to analyze the unsteady flow over a supercritical DRA-2303 airfoil in transonic flow. The dynamic shock wave–boundary layer interaction is one of the most essential features of this unsteady flow causing a distinct oscillation of the flow field. Results from wind-tunnel experiments with a variation of the freestream Mach number at Reynolds numbers ranging from 2.55 to 2.79 × 10⁶ are analyzed regarding the origin and nature of the unsteady shock–boundary layer interaction. Therefore, the TR-SPIV results are analyzed for three buffet flows. One flow exhibits a sinusoidal streamwise oscillation of the shock wave only due to an acoustic feedback loop formed by the shock wave and the trailing-edge noise. The other two buffet flows have been intentionally influenced by an artificial acoustic source installed downstream of the test section to investigate the behavior of the interaction to upstream-propagating disturbances generated by a defined source of noise. The results show that such upstream-propagating disturbances could be identified to be responsible for the upstream displacement of the shock wave and that the feedback loop is formed by a pulsating separation of the boundary layer dependent on the shock position and the sound pressure level at the shock position. Thereby, the pulsation of the separation could be determined to be a reaction to the shock motion and not vice versa.     

Correlation of the separation region length in shock wave/channel boundary layer interaction

 Results of an experimental investigation of the characteristics of a separation region induced by the interaction of an externally generated oblique shock with the turbulent boundary layer formed in a rectangular half channel are discussed. The experiments were carried out in the supersonic wind tunnel of the Institute of Theoretical and Applied Mechanics SB RAS at a free-stream Mach number M _∞=3.01 over a range of Reynolds numbers Re ₁=(9.7–47.5)×10⁶ m^-1 and at zero incidence and zero yaw of the model. Particular attention is paid to the size of the zone of the upstream propagation of disturbances (upstream influence region) under different experimental conditions: with varied values of the shock wave strength, half channel width, and Reynolds number. It is shown, in particular, that the normalized upstream influence region length as a function of inclination angle of the shock generator in a rectangular half channel is readily approximated by a simple exponential function. In support of the known reference data obtained for supersonic numbers M _∞ and moderate Re in other configurations, it is also shown that the upstream influence region length decreases with increasing Reynolds number. Generalization of experimental data on the length of the upstream influence region formed in similar geometric configurations is possible using an additional reference linear scale which is the distance from the leading edge of the shock generator to the exposed surface. A substantial dependence of the reference dimensions of separation region on the half channel width is also established. Received: 20 January 1997/Accepted: 30 June 1997     

Correlation study in shock wave�Cturbulent boundary layer interaction

Shock wave–turbulent boundary layer interaction is a critical problem in aircraft design. Therefore, a thorough understanding of the processes occurring in such flows is necessary. The most important task is to study the unsteady phenomena, in particular, the low-frequency ones, for this interaction. An experimental study of separated flow has been performed in the zone of interaction of the incident oblique shock wave with a turbulent boundary layer at Mach 2. Two-point correlation data in the separation zone and the upstream flow were obtained and showed that low-frequency oscillations of the reflected shock waves are related to pulsations in the inflow turbulent boundary layer.     

Visualizing supersonic inlet duct unstart using planar laser Rayleigh scattering

Planar laser Rayleigh scattering (PLRS) from condensed CO₂ particles is used to visualize flow structure in a Mach 5 wind tunnel undergoing unstart. Detailed flow features such as laminar/turbulent boundary layers and shockwaves are readily illustrated by the technique. A downstream transverse air jet, inducing flow unchoking downstream of the jet, is injected into the free stream flow of the tunnel, resulting in tunnel unstart. Time sequential PLRS images reveal that the boundary layer growth/separation on a surface with a thick turbulent boundary layer, initiated by the jet injection, propagates upstream and produces an oblique unstart shock. The tunnel unstarts upon the arrival of the shock at the inlet. In contrast, earlier flow separation on the opposite surface, initially supporting a thin laminar boundary layer, is observed when a jet induced bow shock strikes that surface. The resulting disturbance to this boundary layer also propagates upstream and precedes the formation of an unstart shock.     

Shock wave triple-point morphology

The existence and characteristics of shock wave triple points are examined. The analysis utilizes a single flow plane for the three shocks and is local to the triple point. It applies when the flow is unsteady, three-dimensional, and the upstream flow is nonuniform. Under more restrictive conditions, a relation is also derived for the ratio of the curvature of the Mach stem to that of the reflected shock. For given values of the ratio of specific heats, γ, and the upstream Mach number, M ₁, a solution window is established. A parametric set of solutions is generated within the window for γ = 1, 1.4, and 5/3 and for 16 values of M ₁ ranging from solution onset to M ₁ = 6.A solution can be one of three types, these stem from the velocity tangency condition along the slip stream. Topics are addressed such as solution multiplicity, shock wave and slip stream orientation, the nature of the reflected wave (weak, strong, inverted, normal), the nature of the Mach stem (weak, strong, normal), and differences due to changes in γ and M ₁.     

Scanning PIV investigation of the laminar separation bubble on a SD7003 airfoil

A laminar separation bubble occurs on the suction side of the SD7003 airfoil at an angle of attack α =  4–8° and a low Reynolds number less than 100,000, which brings about a significant adverse aerodynamic effect. The spatial and temporal structure of the laminar separation bubble was studied using the scanning PIV method at α =  4° and Re = 60,000 and 20,000. Of particular interest are the dynamic vortex behavior in transition process and the subsequent vortex evolution in the turbulent boundary layer. The flow was continuously sampled in a stack of parallel illuminated planes from two orthogonal views with a frequency of hundreds Hz, and PIV cross-correlation was performed to obtain the 2D velocity field in each plane. Results of both the single-sliced and the volumetric presentations of the laminar separation bubble reveal vortex shedding in transition near the reattachment region at Re = 60,000. In a relatively long distance vortices characterized by paired wall-normal vorticity packets retain their identities in the reattached turbulent boundary layer, though vortices interact through tearing, stretching and tilting. Compared with the restricted LSB at Re = 60,000, the flow at Re = 20,000 presents an earlier separation and a significantly increased reversed flow region followed by “huge” vortical structures.     

High-fidelity simulations of moving and flexible airfoils at low Reynolds numbers

The present paper highlights results derived from the application of a high-fidelity simulation technique to the analysis of low-Reynolds-number transitional flows over moving and flexible canonical configurations motivated by small natural and man-made flyers. This effort addresses three separate fluid dynamic phenomena relevant to small fliers, including: laminar separation and transition over a stationary airfoil, transition effects on the dynamic stall vortex generated by a plunging airfoil, and the effect of flexibility on the flow structure above a membrane airfoil. The specific cases were also selected to permit comparison with available experimental measurements. First, the process of transition on a stationary SD7003 airfoil section over a range of Reynolds numbers and angles of attack is considered. Prior to stall, the flow exhibits a separated shear layer which rolls up into spanwise vortices. These vortices subsequently undergo spanwise instabilities, and ultimately breakdown into fine-scale turbulent structures as the boundary layer reattaches to the airfoil surface. In a time-averaged sense, the flow displays a closed laminar separation bubble which moves upstream and contracts in size with increasing angle of attack for a fixed Reynolds number. For a fixed angle of attack, as the Reynolds number decreases, the laminar separation bubble grows in vertical extent producing a significant increase in drag. For the lowest Reynolds number considered (Re _c  = 10⁴), transition does not occur over the airfoil at moderate angles of attack prior to stall. Next, the impact of a prescribed high-frequency small-amplitude plunging motion on the transitional flow over the SD7003 airfoil is investigated. The motion-induced high angle of attack results in unsteady separation in the leading edge and in the formation of dynamic-stall-like vortices which convect downstream close to the airfoil. At the lowest value of Reynolds number (Re _c  = 10⁴), transition effects are observed to be minor and the dynamic stall vortex system remains fairly coherent. For Re _c  = 4 × 10⁴, the dynamic-stall vortex system is laminar at is inception, however shortly afterwards, it experiences an abrupt breakdown associated with the onset of spanwise instability effects. The computed phased-averaged structures for both values of Reynolds number are found to be in good agreement with the experimental data. Finally, the effect of structural compliance on the unsteady flow past a membrane airfoil is investigated. The membrane deformation results in mean camber and large fluctuations which improve aerodynamic performance. Larger values of lift and a delay in stall are achieved relative to a rigid airfoil configuration. For Re _c = 4.85 × 10⁴, it is shown that correct prediction of the transitional process is critical to capturing the proper membrane structural response.     

Mixed Convection Boundary Layer Flow from a Horizontal Circular Cylinder Embedded in a Porous Medium Filled with a Nanofluid

Steady mixed convection boundary layer flow from an isothermal horizontal circular cylinder embedded in a porous medium filled with a nanofluid has been studied for both cases of a heated and cooled cylinder. The resulting system of nonlinear partial differential equations is solved numerically using an implicit finite-difference scheme. The solutions for the flow and heat transfer characteristics are evaluated numerically for various values of the governing parameters, namely the nanoparticle volume fraction φ and the mixed convection parameter λ. Three different types of nanoparticles are considered, namely Cu, Al₂O₃ and TiO₂. It is found that for each particular nanoparticle, as the nanoparticle volume fraction φ increases, the magnitude of the skin friction coefficient decreases, and this leads to an increase in the value of the mixed convection parameter λ which first produces no separation. On the other hand, it is also found that of all the three types of nanoparticles considered, for any fixed values of φ and λ, the nanoparticle Cu gives the largest values of the skin friction coefficient followed by TiO₂ and Al₂O₃. Finally, it is worth mentioning that heating the cylinder (λ > 0) delays separation of the boundary layer and if the cylinder is hot enough (large values of λ > 0), then it is suppressed completely. On the other hand, cooling the cylinder (λ < 0) brings the boundary layer separation point nearer to the lower stagnation point and for a sufficiently cold cylinder (large values of λ < 0) there will not be a boundary layer on the cylinder.     

Effects of Mach number on turbulent separation behaviours induced by blunt fin

 An experimental study of the interaction between shock wave and turbulent boundary layer induced by blunt fin has been carried out at M _∞=7.8 using oil flow visualization and simultaneous measurements of fluctuating wall pressure and heat transfer. This paper presents the effects of Mach number on turbulent separation behaviours induced by blunt fin. Received: 21 July 1996/Accepted: 4 February 1998     

On the mechanism of unsteady shock oscillation in shock wave/turbulent boundary layer interactions

An experimental investigation into the mechanism of shock wave oscillation in compression ramp-generated shock wave/turbulent boundary layer interactions is presented. Particular emphasis is focused upon documenting the respective roles played by both burst-sweep events in the turbulent boundary layer immediately upstream of the interaction and the downstream separated shear layer upon unsteady shock front motion. Unlike the majority of compression ramp experiments which involve bulk separation and large-scale shock motion, consideration is given here to comparatively “weak” interactions in which the streamwise spatial excursion of the shock front is always less than one boundary layer thickness. In this manner any shock motion due to upstream burst-sweep events should be more apparent in relation to that oscillation associated with the separated region. A discrete Hilbert transform-based conditional sampling technique is used to obtain wall pressure measurements conditioned to burst-sweep events. The conditional sampling technique forms the basis by which the instantaneous shock motion is conditioned to the occurrence of upstream bursting. The relationship between the separation bubble and shock motion is also explored in detail. The results of the experiments indicate that the separation bubble represents a first-order effect on shock oscillation. Although it is demonstrated theoretically that the burst-sweep cycle can also give rise to unsteady shock motion of much lower amplitude, the experiments clearly demonstrate that there is no discernible statistical relationship between burst events and spanwise coherent shock front motion.     

Shock-unsteadiness model applied to oblique shock wave/turbulent boundary-layer interaction

Reynolds-averaged Navier–Stokes prediction of shock wave/turbulent boundary layer interactions can yield significant error in terms of the size of the separation bubble. In many applications, this can alter the shock structure and the resulting surface properties. Shock-unsteadiness modification of Sinha et al. (Physics of Fluids, Vol.15, No.8, 2003) has shown potential in improving separation bubble prediction in compression corner flows. In this article, the modification is applied to oblique shock wave interacting with a turbulent boundary layer. The challenges involved in the implementation of the shock-unsteadiness correction in the presence of multiple shock waves and expansion fans are addressed in detail. The results show that a robust implementation of the model yields appreciable improvement over standard k–ω turbulence model predictions.     

Tomographic and time resolved PIV measurements on a finite cylinder mounted on a flat plate

Tomographic and time resolved PIV measurements were performed to examine the 3D flow topology and the flow dynamic above the upper surface of a low-aspect ratio cylinder at Re ≈ 1 ×  10⁵. This generic experiment is of fundamental interest because it represents flow features which are relevant to many applications such as laminar separation bubbles and turbulent reattachment. At Re  ≈ 1 × 10⁵, laminar separation bubbles arise on the side of the cylinder. Furthermore, on the top of the cylinder a separation with reattachment is of major interest. The tomographic PIV measurement, which allows to determine all three velocity components in a volume instantaneously, was applied to examine the flow topology and interaction between the boundary layer and wake structures on the top of the finite cylinder. In the instantaneous flow fields the tip vortices and the recirculation region becomes visible. However, it is also observed that the flow is quite unsteady due to the large separation occurring on the top of the cylinder. In order to study the temporal behaviour of the separation, time resolved PIV was applied. This technique allows capturing the dynamic processes in detail. The development of vortices in the separated shear layer is observed and in addition regions with different dominant frequencies are identified.     

Development of a turbulent boundary layer after a step from smooth to rough surface

The flow developing downstream of a step change from smooth to rough surface condition is studied in the light of Townsend’s wall similarity hypothesis. Previous studies seem to support the hypothesis for channel and pipe flows, but there are considerable controversies about its application to boundary layers and in particular to surface roughness formed by spanwise bars. It has been suggested that this controversy arises from insufficient separation of scales between the boundary layer thickness and the roughness length scale. An experimental investigation has therefore been undertaken where the flow evolves from a fully developed smooth wall boundary layer at high Reynolds numbers over a step in surface roughness (Re _θ = 13,400 at the step). The flow is mapped through the development of the internal layer until the flow is fully developed over the rough wall. The internal layer is found to grow as δ ∼ X ^0.73, and after about 15 boundary layer thicknesses at the step, the internal layer has reached the outer edge of the incoming layer. At the last rough wall measurement station, the Reynolds number has grown to Re _θ ≈ 32,600 and the ratio of boundary layer to roughness length scales is δ/k ≈ 140. The outer layer differences between the smooth and the rough wall data were found to be sufficiently small to conclude that for this setup the Townsend’s wall similarity hypothesis appears to hold.     

Investigation on onset of shock-induced separation

A great number of experimental data indicating shock wave/boundary layer interactions in internal or external supersonic flows were reviewed to make clear the mechanism of the interaction and to decide the onset of shock-induced separation. The interesting conclusions were obtained for the considerably wide range of flow geometries that the onset of separation is independent of the flow geometries and the boundary layer Reynolds number. It is found that the pressure rise necessary to separate the boundary layer in supersonic external flows could be applied to such internal flows as overexpanded nozzles or diffusers. This is due to the fact that the separation phenomenon caused by shock wave/boundary layer interactions is processed through a supersonic deceleration. The shock-induced separation in almost all of interacting flow fields is governed by the concept of free interaction, and the onset of shock-induced separation is only a function of the Mach number just upstream of shock wave. However, physical scales of the produced separation are not independent of the downstream flow fields.     

Bulk flow parameters in turbulent wall boundary layers near separation

Experiments were conducted in a turbulent boundary layer near separation along a flat plate. The pressure gradient in flow direction was varied such that three significant boundary layer configurations could be maintained. The flow in the test section thus had simultaneously a region of favourable pressure gradient, a region of strong adverse pressure gradient with boundary layer separation and a region of reattached boundary layer. Specially designed fine probes facilitated the measurements of skin friction and velocity distribution very close to the wall. Bulk flow parameters such as skin friction coefficient C _f, Reynold's number Re_δ2 and shape factors H and G, which are significant characteristics of wall boundary layers were evaluated. The dependence of these parameters on the Reynolds number and along the test section was explored and the values were compared with other empirical and analytical formulae known in the literature.