Numerical modeling of the receptivity of a supersonic boundary layer to entropy disturbances

Numerical modeling of the receptivity of a two-dimensional flat-plate boundary layer to entropy disturbances is carried out at the freestream Mach number M_∞ = 6. Low-intensity perturbations considered are in the form of temperature spots of various shapes and with different initial positions downstream of the shock. They are shown to be able to generate unstable disturbances in the boundary layer. This receptivity mechanism is relatively weak as compared with the receptivity to acoustic waves. When the entropy perturbations are introduced upstream of the bow shock, they first pass across the shock. Downstream of the shock this interaction generates acoustic waves which, in turn, penetrate into the boundary layer thus exciting unstable disturbances of a considerably greater amplitude than the temperature spots. Thus, the bow shock can change the receptivity mechanism.     

Particle image velocimetry measurements of a shock wave/turbulent boundary layer interaction

Particle image velocimetry is used to investigate the interaction between an incident shock wave and a turbulent boundary layer at Mach 2.1. A particle response assessment establishes the fidelity of the tracer particles. The undisturbed boundary layer is characterized in detail. The mean velocity field of the interaction shows the incident and reflected shock wave pattern, as well as the boundary layer distortion. Significant reversed flow is measured instantaneously, although, on average no reversed flow is observed. The interaction instantaneously exhibits a multi-layered structure, namely, a high-velocity outer region and a low-velocity inner region. Flow turbulence shows the highest intensity in the region beneath the impingement of the incident shock wave. The turbulent fluctuations are found to be highly anisotropic, with the streamwise component dominating. A distinct streamwise-oriented region of relatively large kinematic Reynolds shear stress magnitude appears within the lower half of the redeveloping boundary layer. Boundary layer recovery towards initial equilibrium conditions appears to be a gradual process.     

Interaction of a turbulent boundary layer with a shock wave at transonic flow velocities

An asymptotic theory of the interaction of a turbulent boundary layer on a plate with a normal shock wave of low intensity has been constructed in various studies [1–4] under the assumption that the averaged velocity of the particles in the boundary layer in front of the interaction region satisfies a logarithmic law. In the present paper a different approach to this problem is proposed based on a power law of the velocity in the undisturbed boundary layer. The obtained results give different estimates for not only the sizes of the characteristic flow regions in the interaction region but also for the shock intensity leading to boundary layer separation.     

Shock boundary layer interaction under the influence of a normal suction slot

The phenomenon of shock boundary layer interaction of a shock train under the influence of a normal suction slot is studied. In previous work, it was found that a normal, circumferential suction slot is sufficient to stabilize the primary shock of a shock train in as much as that the back pressure of the shock train can be increased until the shock train gradually changes into a single normal shock. Based on the experimental and numerical results, a flow model was derived which explains the transition of a shock train into a single shock under the influence of boundary layer suction. In this work, the normal shock boundary layer interaction model is validated against flow cases with different upstream Mach and Reynolds numbers. For that purpose three different nozzle flows are investigated at various total pressure levels. In a second step, the flow model is extended to the oblique shock case, correlating the suction mass flow with the total pressure distribution of the incoming boundary layer and the static pressure downstream of the oblique shock. Finally, the influence of the suction cavity pressure onto the shock boundary layer interaction is considered.     

Characteristics of unsteady type IV shock/shock interaction

Characteristics of the unsteady type IV shock/shock interaction of hypersonic blunt body flows are investigated by solving the Navier–Stokes equations with high-order numerical methods. The intrinsic relations of flow structures to shear, compression, and heating processes are studied and the physical mechanisms of the unsteady flow evolution are revealed. It is found that the instantaneous surface-heating peak is caused by the fluid in the “hot spot” generated by an oscillating and deforming jet bow shock (JBS) just ahead of the body surface. The features of local shock/boundary layer interaction and vortex/boundary layer interaction are clarified. Based on the analysis of flow evolution, it is identified that the upstream-propagating compression waves are associated with the interaction of the JBS and the shear layers formed by a supersonic impinging jet, and then the interaction of the freestream bow shocks and the compression waves results in entropy and vortical waves propagating to the body surface. Further, the feedback mechanism of the inherent unsteadiness of the flow field is revealed to be related to the impinging jet. A feedback model is proposed to reliably predict the dominant frequency of flow evolution. The results obtained in this study provide physical insight into the understanding of the mechanisms relevant to this complex flow.     

A calculating method of shock wave oscillating frequency due to turbulent shear layer fluctuations in supersonic flow

One of the more severe fluctuating pressure environments encountered in supersonic orhypersonic flows is the shock wave oscillation driven by interaction of a shock wave withboundary layer.The high intensity oscillating shock wave may induce structure resonanceof a high speed vehicle.The research for the shock oscillation used to adopt empirical orsemiempirical methods because the phenomenon is very complex.In this paper atheoretical solution on shock oscillating frequency due to turbulent shear layer fluctuationshas been obtained from basic conservation equations.Moreover,we have attained theregularity of the frequency of oscillating shock varying with incoming flow Mach numbersM_∞and turning angleθ.The calculating results indicate excellent agreement withmeasurements.This paper has supplied a valuable analytical method to study aeroelasticproblems produced by shock wave oscillation.     

Aerodynamic heating in the region of shock and turbulent boundary layer interaction induced by a cylinder

Detailed distributions of heat flux in the region of shock wave and turbulent boundary layer interaction induced by a cylinder were measured in the shock tunnel. Oil flow patterns and Schlieren photographs were taken. Empirical relations were given for determining separation shock angle, peaks of heat flux and their locations on both cylinder leading edge and flat plate surface, and other characteristic parameters of the interaction region.     

Interaction between the boundary layer and a supersonic flow when part of the wall is rapidly heated

There have been many theoretical studies of aspects of the unsteady interaction of an exterior inviscid flow with a boundary layer [1–9]. The mathematical flow models obtained in these studies by the method of matched asymptotic expansions describe a wide range of phenomena observed experimentally. These include boundary layer separation near the hinge of a flap, the flow in the neighborhood of the trailing edge of an oscillating airfoil [1–2], and the development and propagation of perturbations in a boundary layer excited by an oscillating wall or some other way [3–5]. The present paper studies the interaction of an unsteady boundary layer with a supersonic flow when a small part of the surface of a body in the flow is rapidly heated.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 66–70, January–February, 1984.     

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.     

Numerical Modeling of the Development of a Streaky Structure in a Boundary Layer at a High Freestream Turbulence Level

The disturbances generated by external turbulence in the boundary layer on a flat plate set suddenly in motion are determined. A turbulent flow calculated by direct numerical simulation is taken as the initial conditions. The solution obtained simulates the initial stage of laminar-turbulent transition in the flat-plate boundary layer at a high turbulence level in the oncoming flow. The solution makes it possible to estimate the effects of different factors, such as nonstationarity, nonlinearity, and the parameters of the freestream velocity fluctuation spectrum, on disturbance enhancement in the boundary layer.     

Study of the shock motion in a hypersonic shock system/turbulent boundary layer interaction

Wall pressure fluctuations and surface heat transfer signals have been measured in the hypersonic turbulent boundary layer over a number of compression-corner models. The distributions of the separation shock oscillation frequencies and periods have been calculated using a conditional sampling algorithm. In all cases the oscillation frequency distributions are of broad band, but the most probable frequencies are low. The VITA method is used for deducing large scale disturbances at the wall in the incoming boundary layer and the separated flow region. The results at present showed the existence of coherent structures in the two regions. The zero-cross frequencies of the large scale structures in the two regions are of the same order as that of the separation shock oscillation. The average amplitude of the large scale structures in the separated region is much higher than that in the incoming boundary layer. The length scale of the separation shock motion region is found to increase with the disturbance strength. The results show that the shock oscillation is of inherent nature in the shock wave/turbulent boundary layer interaction with separation. The shock oscillation is considered to be the consequence of the coherent structures in the separated region.This work was supported by the Chinese National Science Foundation. Thanks for Prof. Z. B. Lin and Miss X. Y. Feng for their helps. The authors wish to express thanks to Professor W. Merzkirch who has helped us to check the paper again and again.     

Supersonic flow past a wedge on a flat plate. Laminar boundary layer separation and reattachment

Supersonic laminar flow past a two-dimensional “flat-plate/wedge“ configuration is numerically investigated. The pressures at the boundary layer separation and reattachment points are calculated over wide Mach and Reynolds number ranges. The minimum angles of the wedge surface inclination at which a return flow occurs are determined. The results are presented in the form of generalized Mach-number-dependences of the theoretical pressure on the wedge surface initiating boundary layer separation and the pressure at the boundary layer reattachment point.     

The effects of wall roughness on the separated flow over a smoothly contoured ramp

Experimental measurements address the effects on a turbulent boundary layer of wall roughness on a flat plate and a ramp that produces a separation bubble over the ramp trailing edge. A fully rough flow condition is achieved on the upstream flat plate. The main effect of the wall roughness on the outer layer turbulence on a flat plate is to change the friction velocity. The separation region is substantially larger for the rough-wall case. The rough-wall boundary layer turbulence is less sensitive to the onset of an adverse pressure gradient over the ramp, producing substantially smaller Reynolds stress peaks in upstream flat-plate, wall-unit coordinates.     

Reynolds number effects on a turbulent boundary layer with separation,reattachment, and recovery

The present paper addresses experimental studies of Reynolds number effects on a turbulent boundary layer with separation, reattachment, and recovery. A momentum thickness Reynolds number varies from 1,100 to 20,100 with a wind tunnel enclosed in a pressure vessel by varying the air density and wind tunnel speed. A custom-built, high-resolution laser Doppler anemometer provides fully resolved turbulence measurements over the full Reynolds number range. The experiments show that the mean flow is at most a very weak function of Reynolds number while turbulence quantities strongly depend on Reynolds number. Roller vortices are generated in the separated shear layer caused by the Kelvin–Helmholtz instability. Empirical Reynolds number scalings for the mean velocity and Reynolds stresses are proposed for the upstream boundary layer, the separated region, and the recovery region. The inflectional instability plays a critical role in the scaling in the separated region. The near-wall flow recovers quickly downstream of reattachment even if the outer layer is far from an equilibrium state. As a result, a stress equilibrium layer where a flat-plate boundary layer scaling is valid develops in the recovery region and grows outward moving downstream.     

Experimental investigation of three-dimensional supersonic flow of a gas in the interference region of intersecting surfaces

A contemporary high-speed aircraft represents a complex three-dimensional configuration, where supersonic gas flow is accompanied by numerous local flow interaction zones, in particular, near the intersection of different surfaces. Such a flow is characterized by three-dimensional systems of shock and expansion waves, and close to the surfaces one finds interaction of boundary layers and, above all, interaction of shock waves with the boundary layer. In general, the angular configurations are formed by intersection or contact of nonplanar surfaces with swept-back or blunted leading edges. This makes it practically impossible to obtain a rigorous theoretiical solution to the problem of gas flow over these surfaces, and presents considerable difficulty in an experimental investigation. It is therefore of interest to study the physical features of gas flow in corner configurations of very simple form [1–3]. The present paper examines the results of an experimental investigation of typical features of symmetric and asymmetric interaction of compressive, expansive, and mixed flows in the interference region of planar surfaces intersecting at an angle of less than 180?.     

Global stability analysis of axisymmetric boundary layer over a circular cylinder

This paper presents a linear global stability analysis of the incompressible axisymmetric boundary layer on a circular cylinder. The base flow is parallel to the axis of the cylinder at inflow boundary. The pressure gradient is zero in the streamwise direction. The base flow velocity profile is fully non-parallel and non-similar in nature. The boundary layer grows continuously in the spatial directions. Linearized Navier–Stokes (LNS) equations are derived for the disturbance flow quantities in the cylindrical polar coordinates. The LNS equations along with homogeneous boundary conditions forms a generalized eigenvalues problem. Since the base flow is axisymmetric, the disturbances are periodic in azimuthal direction. Chebyshev spectral collocation method and Arnoldi’s iterative algorithm is used for the solution of the general eigenvalues problem. The global temporal modes are computed for the range of Reynolds numbers and different azimuthal wave numbers. The largest imaginary part of the computed eigenmodes is negative, and hence, the flow is temporally stable. The spatial structure of the eigenmodes shows that the disturbance amplitudes grow in size and magnitude while they are moving towards downstream. The global modes of axisymmetric boundary layer are more stable than that of 2D flat-plate boundary layer at low Reynolds number. However, at higher Reynolds number they approach 2D flat-plate boundary layer. Thus, the damping effect of transverse curvature is significant at low Reynolds number. The wave-like nature of the disturbance amplitudes is found in the streamwise direction for the least stable eigenmodes.     

The dynamics of coherent structures in a flat-plate boundary layer

From the data of a direct numerical simulation the three-dimensional coherent structures of an incompressible, spatially evolving flat-plate boundary layer have been calculated using the POD method. By Galerkin projection of the Navier-Stokes equations onto the corresponding system of eigenfunctions then a low-dimensional model of the flow in the form of a system of ODE's has been derived. In a region of the boundary layer just beyond the spike stages of transition this system displays deterministic chaos that has been quantified by determining its Lyapunov exponents.     

Numerical prediction of locally forced turbulent boundary layer

An unsteady numerical simulation was performed to analyze flow structure behind a local suction/blowing in a flat-plate turbulent boundary layer. The local forcing was given to the boundary layer flow by means of a sinusoidally oscillating jet. A version of the unsteady k––f_μ model [Fluid Dyn. Res. 26 (6) (2000) 421] was employed. The Reynolds number based on the momentum thickness was about Re_θ=1700. The forcing frequency was varied in the range 0.011f⁺0.044 with a fixed forcing amplitude A_o=0.4. The predicted results were compared and validated with the experimental data. It was shown that the unsteady locally forced boundary layer flow is predicted well by the k––f_μ model. The time-dependent numerical flow visualizations were demonstrated during one period of the local forcing. The effect of the pitch angle of local forcing on the reduction of skin friction was examined.     

Effects of micro-ramps on a shock wave/turbulent boundary layer interaction

Stereoscopic particle image velocimetry is used to investigate the effects of micro-ramp sub-boundary layer vortex generators, on an incident shock wave/boundary layer interaction at Mach 1.84. Single- and double-row arrangements of micro-ramps are considered. The micro-ramps have a height of 20% of the unperturbed boundary layer thickness and the measurement planes are located 0.1 and 0.6 boundary layer thicknesses from the wall. The micro-ramps generate packets of individual vortex pairs downstream of their vertices, which produce counter-rotating longitudinal streamwise vortex pairs in a time-averaged view. These structures induce a pronounced spanwise variation of the flow properties, namely the mixing across the boundary layer interface. The probability of reversed-flow occurrence is decreased by 20 and 30% for the single- and double-row configurations, respectively. Both configurations of micro-ramps stabilize the shock motion by reducing the length of its motion by about 20% in the lower measurement plane. The results are summarized by a conceptual model describing the boundary layer’s and interaction’s flow pattern under the effect of the micro-ramps.     

激波风洞内超燃冲压发动机三面压缩进气道流场实验观测

主要进行了超燃冲压发动机三面压缩进气道的实验观测。利用来流马赫数4.5的直通式激波风洞,考察了三组具有不同压缩角度的进气道模型内部的流场情况。实验观测手段为油流法、丝线法和高速纹影,同时,辅以数值模拟以有助于流场细节分析。纹影照片展示了进气道内部以激波边界层相互作用为主要影响因素的流场复杂结构,数值模拟也显示了相近的结果。油流技术与丝线法显示了近壁面处的流动图像,照片中可见激波、分离线、再附线等分界线位置。根据实验结果,可以推测唇口激波与进气道内边界层的相互作用及其引起的壁面分离是影响进气道内流动的主要因素。同时,尝试了利用抽吸方法减弱激波与边界层相互作用诱发的壁面流动分离,并取得一定结果。