Investigation of laminar and turbulent boundary layers interacting with externally generated shock waves
Laminar and turbulent boundary layer interactions with externally generated shock waves
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Laminar and turbulent boundary layer interactions with externally generated shock waves
An assessment is presented for the experimental data on separated flow in shock wave turbulent boundary layer interactions at hypersonic and supersonic speeds. The data base consists mainly of two dimensional and axisymmetric interactions in compression corners or cylinder-flares, and externally generated oblique shock interactions with boundary layers over flat plates or cylindrical surfaces. The conditions leading to flow separation and the subsequent changes in the flow empirical correlations for incipient separation are reviewed. The effects of the Mach number, Reynolds number, surface cooling and the methods of detecting separation are discussed. The pertinent experimental data for the separated flow characteristics in separated turbulent boundary layer shock interaction are also presented and discussed.
The solar wind disturbances of January 1 to 7, 1978 are reconstructed by a modeling method. First, the interplanetary magnetic field (IMF) background pattern, including a corotating shock, is reproduced using the Stanford source surface map. Then, two solar flares with their onset times on January 1, 0717 UT at S17 deg E10 deg and 2147 UT S17 deg E32 deg, respectively, are selected to generate two interplanetary transient shocks. It is shown that these two shocks interacted with the corotating shock, resulting in a series of interplanetary events observed by four spacecraft, Helios 1 and 2, IMP-8 (Interplanetary Monitoring Platform 8), and Voyager 2. Results show that these three shock waves interact and coalesce in interplanetary space such that Helios 2 and Voyager 2 observed only one shock and Helios 1 and IMP-8 observed two shocks. All shocks observed by the four spacecraft, except the corotating shock at Helios 1, are either a transient shock or a shock which is formed from coalescing of the transient shocks with the corotating shock. The method is useful in reconstructing a very complicated chain of interplanetary events observed by a number of spacecraft.
The solar wind disturbances of January 1 to 7, 1978 are reconstructed by a modeling method. First, the interplanetary magnetic field (IMF) background pattern, including a corotating shock, is reproduced using the Stanford source surface map. Then, two solar flares with their onset times on January 1, 0717 UT at S17 deg E10 deg and 2147 UT S17 deg E32 deg, respectively, are selected to generate two interplanetary transient shocks. It is shown that these two shocks interacted with the corotating shock, resulting in a series of interplanetary events observed by four spacecraft, Helios 1 and 2, IMP-8 (Interplanetary Monitoring Platform 8), and Voyager 2. Results show that these three shock wave interact and coalesce in interplanetary space such that Helios 2 and Voyager 2 observed only one shock and Helios 1 and IMP-8 observed two shocks. All shocks observed by the four spacecraft, except the corotating shock at Helios 1, are either a transient shock or a shock which is formed from coalescing of the transient shocks with the corotating shock. The method is useful in reconstructing a very complicated chain of interplanetary events observed by a number of spacecraft.
A computation procedure is described for predicting the flowfields which develop when successive interactions between oblique shock waves and a turbulent boundary layer occur. Such interactions may occur, for example, in engine inlets for supersonic aircraft. Computations are carried out for axisymmetric internal flows at M 3.82 and 2.82. The effect of boundary layer bleed is considered for the M 2.82 flow. A control volume analysis is used to predict changes in the flow field across the interactions. Two bleed flow models have been considered. A turbulent boundary layer program is used to compute changes in the boundary layer between the interactions. The results given are for flows with two shock wave interactions and for bleed at the second interaction site. In principle the method described may be extended to account for additional interactions. The predicted results are compared with measured results and are shown to be in good agreement when the bleed flow rate is low (on the order of 3% of the boundary layer mass flow), or when there is no bleed. As the bleed flow rate is increased, differences between the predicted and measured results become larger. Shortcomings of the bleed flow models at higher bleed flow rates are discussed.
The swept oblique shock-wave/turbulent-boundary-layer interaction generated by a 20-deg sharp fin at Mach 4 and Reynolds number 21,000 is investigated via a series of computations using both conical and three-dimensional Reynolds-averaged Navier-Stokes equations with turbulence incorporated through the algebraic turbulent eddy viscosity model of Baldwin-Lomax. Results are compared with known experimental data, and it is concluded that the computed three-dimensional flowfield is quasi-conical (in agreement with the experimental data), the computed three-dimensional and conical surface pressure and surface flow direction are in good agreement with the experiment, and the three-dimensional and conical flows significantly underpredict the peak experimental skin friction. It is pointed out that most of the features of the conical flowfield model in the experiment are observed in the conical computation which also describes the complete conical streamline pattern not included in the model of the experiment.
A complex two-dimensional, unsteady, viscous hypersonic shock wave interaction is numerically simulated by a high-resolution, second-order fully implicit shock-capturing scheme. The physical model consists of a nonstationary oblique shock impinging on the bow shock of a blunt body. Studies indicated that the unsteady flow patterns are slightly different from their steady counterparts. However, for the sample cases investigated the peak surface pressures for the unsteady flows seem to occur at very different impingement locations than for the steady flow cases.
A complex two-dimensional, unsteady, viscous hypersonic shock wave interaction is numerically simulated by a high-resolution, second-order fully implicit shock-capturing scheme. The physical model consists of a nonstationary oblique shock impinging on the bow shock of a blunt body. Studies indicate that the unsteady flow patterns are slightly different from their steady counterparts. However, for the sample cases investigated the peak surface pressures for the unsteady flows seem to occur at very different impingement locations than for the steady flow cases.
Abstract In this work, three-dimensional Euler–Lagrange (EL) point-particle simulations of a shock wave interacting with a fixed bed of particles are carried out. The results from the particle-resolved (PR) simulations are used to assess the performance of the point-particle drag models during short time scales. We demonstrate that in a one-way coupled regime, the point-particle simulations recover the dominant gas dynamic features of the flow and are in a good agreement with the exact Riemann solution of a shock traveling through a sudden area contraction. Although the PR simulations are inviscid, we show that a dissipative drag is necessary to predict the mean behavior of the gas. As a model for the inviscid shock-induced (SI) drag two different models are presented in lieu of the quasi-steady drag. Finally, two-way coupled simulations are performed at four different particle volume fractions {0.10, 0.15, 0.20, 0.25} and three different incident shock Mach numbers {1.22, 1.66, 3.0} and compared against the data from PR inviscid simulations. At a lower Mach number (1.22), averaged flow quantities from the two-way coupled simulations agree well with the PR simulations. As the Mach number increases, we observe that the discrepancies between the point-particle and the PR simulations grow. A sensitivity analysis of the drag models involved reveals a strong influence of the inviscid-unsteady force on the gas quantities especially in the case of a strong shock interacting with a dense bed of particles. The use of Mach correlation beyond the subcritical regime coupled with the model for volume fraction correction is identified as a probable cause for the additional drag.
One of the main themes in fluid dynamics at present and in the future is going to be computational fluid dynamics with the primary focus on the determination of drag, flow separation, vortex flows, and unsteady flows. A computation of the flow of a viscous fluid requires an understanding and consideration of the physical aspects of the flow. This is done by identifying the flow regimes and the scales of fluid motion, and the sources of vorticity. Discussions of flow regimes deal with conditions of incompressibility, transitional and turbulent flows, Navier-Stokes and non-Navier-Stokes regimes, shock waves, and strain fields. Discussions of the scales of fluid motion consider transitional and turbulent flows, thin- and slender-shear layers, triple- and four-deck regions, viscous-inviscid interactions, shock waves, strain rates, and temporal scales. In addition, the significance and generation of vorticity are discussed. These physical aspects mainly guide computations of the flow of a viscous fluid.
A test was conducted in the 15 cm x 15 cm supersonic wind tunnel at NASA Glenn Research Center that focused on corner effects of an oblique shock-wave/boundary-layer interaction. In an attempt to control the interaction in the corner region, eight corner fillet configurations were tested. Three parameters were considered for the fillet configurations: the radius, the fillet length, and the taper length from the square corner to the fillet radius. Fillets effectively reduced the boundary-layer thickness in the corner; however, there was an associated penalty in the form of increased boundary-layer thickness at the tunnel centerline. Larger fillet radii caused greater reductions in boundary-layer thickness along the corner bisector. To a lesser, but measureable, extent, shorter fillet lengths resulted in thinner corner boundary layers. Overall, of the configurations tested, the largest radius resulted in the best combination of control in the corner, evidenced by a reduction in boundary-layer thickness, coupled with minimal impacts at the tunnel centerline.
The radio emission observed from SN 1987A can be interpreted as the result of the shock wave interaction with the wind from the blue supergiant progenitor star. There is now evidence at IR and UV wavelengths for the presence of dense circumstellar gas, which may be at a distance of about 10 to the 18th cm from the SN. If this interpretation is correct, the IR dust echo should remain at about constant luminosity for 1 to 2 years. The dense gas is probably the remnant of a dense wind from the progenitor star lost during an earlier red supergiant evolutionary phase. A search for luminous B stars with evidence of dense circumstellar gas at some distance from the star is warranted. The detection of light echoes formed by dust scattering in clouds near SN 1987A is promising. The echoes may appear as luminous overlapping annular regions with a size scale of order 1 arcmin or less.
Nonintrusive 3-D measurements were made of a normal shock wave/turbulent boundary layer interaction. The measurements were made through a quadrant of a square test section of a continuous flow supersonic wind tunnel in which a normal shock wave had been stabilized. Two dimensional measurements were made throughout the interaction region while 3-D measurements were made near the corner in the vicinity of the shock where 3-D flow effects were expected to be significant. Laser Doppler velocimetry, surface static pressure measurements, and flow visualization techniques were used for two freestream nominal Mach number test cases: 1.6 and 1.3. No turbulence information was obtained. The mean velocity measurements were converted to Mach number by recording the wind tunnel total temperature. Some shock oscillation was present during both of the test cases. After startup of the wind tunnel, the shock oscillated with an amplitude of approx. + or - 1 cm, however, after reaching steady condition, the shock oscillation amplitude was greatly reduced, as evidenced by the laser velocimeter results. The Mach 1.3 test case resulted in a nearly uniform flow without secondary shock waves and with no or very isolated corner separation.
An upwind finite element technique that uses cell centered quantities and implicit and/or explicit time marching has been developed for computing hypersonic laminar viscous flows using adaptive unstructured triangular grids. A structured grid of quadrilaterals is laid out near the body surface. For inviscid flows the method is stable at Courant numbers of over 100,000. A first order basic scheme and a higher order flux corrected transport (FCT) scheme have been implemented. This technique has been applied to the problem of predicting type III and IV shock wave interactions on a cylinder, with a view of simulating the pressure and heating rate augmentation caused by an impinging shock on the leading edge of a cowl lip of an engine inlet. The predictions of wall pressure and heating rates compare very well with experimental data. The flow features are very distinctly captured with a sequence of adaptively generated grids. The adaptive mesh generator and the upwind Navier-Stokes solver are combined in a set of programs called LARCNESS, an acronym for Langley Adaptive Remeshing Code and Navier-Stokes Solver.
Large-scale molecular dynamics (MD) simulations were carried out to investigate the shock-induced evolution of microstructure in Fe-based systems comprising single-crystal and layered Cu/Fe alloys with a distribution of interfaces. The shock compression of pure single-crystal Fe oriented along [110] above a threshold pressure results in a BCC (α)→HCP (ε) phase transformation behavior that generates a distribution of ε phase variants in the phase transformed region of the microstructure behind the shock front. The propagation of the release wave through a phase transformed ε phase causes a reverse ε→α phase transformation and renders a distribution of twins for the [110] oriented Fe that serve as void nucleation sites during spall failure. The simulations reveal that the α→ε→α transformation-induced twinning for shock loading along the [110] direction is due to a dominant ε phase variant formed during compression that rotates on the arrival of the release wave followed by a reverse phase transformation to twins in the α phase. The modifications in the evolution of the ε phase variants and twins in Fe behavior are also studied for Cu–Fe layered microstructures due to the shock wave interactions with the Cu/Fe interfaces using a newly constructed Cu–Fe alloy potential. Here, the MD simulations suggest that interfaces affect the observed variants during shock compression and, hence, distributions of twins during shock release that affects the void nucleation stresses in the Fe phase of Cu/Fe microstructures.
Fundamental research for sonic boom reduction is needed to quantify the interaction of shock waves generated from the aircraft wing or tail surfaces with the nozzle exhaust plume. Aft body shock waves that interact with the exhaust plume contribute to the near-field pressure signature of a vehicle. The plume and shock interaction was studied using computational fluid dynamics and compared with experimental data from a coaxial convergent-divergent nozzle flow in an open jet facility. A simple diamond-shaped wedge was used to generate the shock in the outer flow to study its impact on the inner jet flow. Results show that the compression from the wedge deflects the nozzle plume and shocks form on the opposite plume boundary. The sonic boom pressure signature of the nozzle exhaust plume was modified by the presence of the wedge. Both the experimental results and computational predictions show changes in plume deflection.
The three-dimensional oblique shock wave-turbulent boundary layer interaction generated by a sharp fin attached to a flat plate is investigated experimentally and theoretically for Mach 3 and Reynolds number 9 x 10 to the 5th using two different models. Both models employ the three-dimensional compressible Navier-Stokes equations in mass-averaged variables; one model utilizes the algebraic turbulent eddy viscosity model of Baldwin and Lomax (1978), while the other model employs the two-equation turbulence model of Jones and Launder (1972) coupled with the wall function model of Viegas and Rubesin (1985). The computed surface pressure, surface streamlines, pitot pressure, and yaw angle profiles are found to be in good agreement with experimental data. The three-dimensional velocity fields computed by both models are in close agreement, although the eddy viscosity profiles differ significantly within the three-dimensional interaction. This result indicates that the overall structure of this three-dimensional sharp fin interaction is insensitive to the turbulence model.
An experiment, thoroughly documenting the flow field resulting from the interaction of a shock wave with a nonadiabatic hypersonic turbulent boundary layer, is described. Detailed mean flow and surface data are presented for two shock strengths resulting in attached and separated flows, respectively. The surface measurements include continuous pressure, shear and heat-flux distributions upstream, in, and downstream of the interaction regions. At closely spaced intervals along the surface, boundary-layer profiles of static and pitot pressure and total temperature were obtained from which velocity, density and static temperature profiles were derived. The data are presented in both graphical and tabular form. These data are of sufficient detail to validate advanced computer codes and their associated turbulence models.