Search NASA⌕ Search

SEARCH · Search NASA

Results for “Shock Interactions”

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.

At least 109 records · Page 6

Interference heating from interactions of shock waves with turbulent boundary layers at Mach 6

An experimental investigation of interference heating resulting from interactions of shock waves and turbulent boundary layers was conducted. Pressure and heat-transfer distributions were measured on a flat plate in the free stream and on the wall of the test section of the Langley Mach 6 high Reynolds number tunnel for Reynolds numbers ranging from 2 million to 400 million. Various incident shock strengths were obtained by varying a wedge-shock generator angle (from 10 deg to 15 deg) and by placing a spherical-shock generator at different vertical positions above the instrumented flat plate and tunnel wall. The largest heating-rate amplification factors obtained for completely turbulent boundary layers were 22.1 for the flat plate and 11.6 for the tunnel wall experiments. Maximum heating correlated with peak pressures using a power law with a 0.85 exponent. Measured pressure distributions were compared with those calculated using turbulent free-interaction pressure rise theories, and separation lengths were compared with values calculated by using different methods.

Johnson, C. B.↗

Shock interference prediction using direct simulation Monte Carlo

The shock interaction produced when an incident shock impinges on an inlet cowl lip is investigated for a cowl lip radius of 0.1 in. and flight conditions of Mach 15 and 35-km altitude. This problem is of interest in the design of hypersonic flight vehicles because the interference heating at the cowl lip surface from such a shock interaction is expected to be very high and the design of adequate thermal protection in this region will be challenging. Noncontinuum effects may be significant at this combination of altitude and small dimensions. If this is so, a computational solution based on the Navier-Stokes equations could significantly overpredict the surface heating on the cowl lip. Therefore, a kinetic theory solution technique, direct simulation Monte Carlo (DSMC), is used to predict the flowfield and surface heating. The details of modeling this problem and some preliminary results are presented. The same qualitative flow patterns are seen with DSMC and a continuum solution of the shock interference case, but the peak heating predicted with DSMC is somewhat lower.

Carlson, Ann B.↗

Numerical Study of Boundary Layer Interaction with Shocks: Method Improvement and Test Computation

The objective is the development of a high-order and high-resolution method for the direct numerical simulation of shock turbulent-boundary-layer interaction. Details concerning the spatial discretization of the convective terms can be found in Adams and Shariff (1995). The computer code based on this method as introduced in Adams (1994) was formulated in Cartesian coordinates and thus has been limited to simple rectangular domains. For more general two-dimensional geometries, as a compression corner, an extension to generalized coordinates is necessary. To keep the requirements or limitations for grid generation low, the extended formulation should allow for non-orthogonal grids. Still, for simplicity and cost efficiency, periodicity can be assumed in one cross-flow direction. For easy vectorization, the compact-ENO coupling algorithm as used in Adams (1994) treated whole planes normal to the derivative direction with the ENO scheme whenever at least one point of this plane satisfied the detection criterion. This is apparently too restrictive for more general geometries and more complex shock patterns. Here we introduce a localized compact-ENO coupling algorithm, which is efficient as long as the overall number of grid points treated by the ENO scheme is small compared to the total number of grid points. Validation and test computations with the final code are performed to assess the efficiency and suitability of the computer code for the problems of interest. We define a set of parameters where a direct numerical simulation of a turbulent boundary layer along a compression corner with reasonably fine resolution is affordable.

Adams, N. A.↗

Numerical study of boundary layer interaction with shocks: Method and code validation

A major problem in modeling of turbulent supersonic flows is the correct assessment of viscous-inviscid interaction problems. Of particular interest is the interaction of boundary layers with shocks. Present turbulence models give in most cases unsatisfactory results in the region of rapid distortion and in the separation region (if one is present) in particular with regard to mean flow profiles and turbulence quantities. The objective of the present work is the direct numerical simulation of shock boundary layer interaction. This report summarizes the first phase during which a numerical method suitable for this problem has been developed and a computer code has been written and tested.

Adams, Nikolaus A.↗

Self-sustained oscillations of a shock wave interacting with a boundary layer on a supercritical airfoil

A theory is proposed of the self-sustaining oscillations of a weak shock on an airfoil in steady, transonic flow. The interaction of the shock with the boundary layer on the airfoil produces displacement thickness fluctuations which convect downstream and generate sound by interaction with the trailing edge. A feedback loop is established when this sound impinges on the shock wave, resulting in the production of further fluctuations in the displacement thickness. The details are worked out for an idealized mean boundary layer velocity profile, but strong support for the basic hypotheses of the theory is provided by a comparison with recent experiments involving the generation of acoustic "tone bursts' by a supercritical airfoil section.

Ventres, C. S.↗

Self-sustained oscillations of a shock wave interacting with a boundary layer on a supercritical airfoil

A theory is proposed of the self-sustaining oscillations of a weak shock on an airfoi in steady, transonic flow. The interaction of the shock with the boundary layer on the airfoil produces displacement thickness fluctuations which convect downstream and generate sound by interaction with the trailing edge. A feedback loop is established when this sound impinges on the shock wave, resulting in the production of further fluctuations in the displacement thickness. The details are worked out for an idealized mean boundary layer velocity profile, but strong support for the basic hypotheses of the theory is provided by a comparison with recent experiments involving the generation of acoustic 'tone bursts' by a supercritical airfoil section.

Ventres, C. S.↗

Statistical Properties of Line Centroid Velocity Increments in the rho Ophiuchi Cloud

We present a comparison of histograms of CO (2-1) line centroid velocity increments in the rho Ophiuchi molecular cloud with those computed for spectra synthesized from a three-dimensional, compressible, but non-starforming and non-gravitating hydrodynamic simulation. Histograms of centroid velocity increments in the rho Ophiuchi cloud show clearly non-Gaussian wings, similar to those found in histograms of velocity increments and derivatives in experimental studies of laboratory and atmospheric flows, as well as numerical simulations of turbulence. The magnitude of these wings increases monotonically with decreasing separation, down to the angular resolution of the data. This behavior is consistent with that found in the phase of the simulation which has most of the properties of incompressible turbulence. The time evolution of the magnitude of the non-Gaussian wings in the histograms of centroid velocity increments in the simulation is consistent with the evolution of the vorticity in the flow. However, we cannot exclude the possibility that the wings are associated with the shock interaction regions. Moreover, in an active starforming region like the rho Ophiuchi cloud, the effects of shocks may be more important than in the simulation. However, being able to identify shock interaction regions in the interstellar medium is also important, since numerical simulations show that vorticity is generated in shock interactions.

Lis, D. C.↗

Tetrahedral-Mesh Simulations of Shock-Turbulence Interaction

Despite decades of development of unstructured mesh methods, direct numerical simulations (DNS) of turbulent flows are still predominantly performed on structured or unstructured hexahedral meshes with high-order finite-difference methods, weighted essentially nonoscillatory (WENO) schemes, or hybrid schemes formed by their combinations. Tetrahedral meshes offer easy mesh generation and adaptation around complex geometries and the potential of an orientation-free grid that would benefit the isotropic nature of small-scale dissipation, as well as the solution accuracy of intermediate scales. To advance the state of the art of unstructured-mesh simulation capabilities for shock/turbulence interaction, DNS using pure tetrahedral meshes are carried out with the space-time conservation element, solution element (CESE) method in this research. By its design, the CESE method is constructed based on a non-dissipative scheme and is a genuinely multidimensional numerical framework that is free from the use of an approximate Riemann-solver. The numerical framework also provides the ability to add numerical dissipation (the nondissipative scheme acts as the reference state like that of the reversible state in thermodynamics) when needed (with justification from mathematics/physics). The above-mentioned features along with the CESE method's consistent shock-capturing approach and strong enforcement of flux conservation in spacetime offers a novel method to accurately simulate turbulent flows and their interaction with shocks using tetrahedral meshes. Two canonical problems, namely, isotropic turbulence interaction with a normal shock and a Mach 2.9 turbulent boundary layer flow over a 24deg compression corner are investigated in this study. Computational results show reasonably good agreement with experimental data and results from structured-mesh, high-order simulations available in the literature. Successful validation of these canonical problems demonstrated here paves the way for future high-fidelity supersonic flow simulations involving complex-geometries.

Venkatachari, Balaji Shankar↗

Internal corner flow fields

A computational procedure has been developed to predict the inviscid super/hypersonic flow field of conical internal corners. The prediction of internal corner flow fields can be important in the design of supersonic 'box' type inlets. The computational procedure utilizes a second order finite difference marching technique to asymptote to the conical corner flow solution of Euler's equations. These flow fields are dominated by complex shock interactions. All discontinuities, shocks and slip surfaces are fitted with the appropriate jump conditions. The 'triple' points (the interaction of two shocks and a slip surface) are also computed exactly. Computed results are compared with experimental data and the computational results of other investigators. In addition, the sensitivity of these flow fields to a number of geometric parameters is studied, and the impact of these flows on inlet performance is assessed.

Marconi, F.↗

Direct numerical simulation and analysis of shock turbulence interaction

Two kinds of linear analysis, rapid distortion theory (RDT) and linear interaction analysis (LIA), were used to investigate the effects of a shock wave on turbulence. Direct numerical simulations of two-dimensional isotropic turbulence interaction with a normal shock were also performed. The results from RDT and LIA are in good agreement for weak shock waves, where the effects of shock front curvature and shock front unsteadiness are not significant in producing vorticity. The linear analyses predict wavenumber-dependent amplification of the upstream one-dimensional energy spectrum, leading to turbulence scale length scale decrease through the interaction. Instantaneous vorticity fields show that vortical structures are enhanced while they are compressed in the shock normal direction. Entrophy amplfication through the shock wave compares favorably with the results of linear analyses.

Lee, Sangsan↗

Viscous computation of a space shuttle flow field

A procedure is presented, as well as some results, to calculate the flow over the winged orbiter. This necessitates the use of two computer codes. A parabolized marching Navier-Stokes code is used to obtain the solution up to the bow shock-wing shock interaction region and for the region after the interaction. An unsteady Navier-Stokes code is to be used in the region of the shock interaction. Only resuls for the marching code are presented. For the flow conditions calculated, M infinity = 7.9, alpha = 25 deg, T(wall) = 540 R, Re(L) = 60728 per inch, laminar or turbulent, the PNS code was marched up to an X/L = 0.7 which is where the bow shock-wing shock interaction region occurs.

Chaussee, D. S.↗

Interaction of global merged interaction region shock with the heliopause and its relation to the 2- and 3-kHz radio emissions

We use the Voyager 2 plasma and magnetic field data together with a one-fluid magnetohydrodynamics model to study the interactions of the 1991 global merge interaction region (GMIR) shock with the termination shock and the heliopause. The 1991 GMIR is an extraordinarily large global solar wind structure in radial, longitudinal, and latitudinal extents. It has a strong shock at the leading edge. After its penetration through the termination shock, the GMIR shock first propagates through the subsonic solar wind, then interacts with the heliopause. The interaction produces a transmitted shock propagating outward in the interstellar medium, and a reflected shock propagating backward in the subsonic solar wind. We identify the reflected shock and the transmitted shock as the possible source of the radio noise detected at Voyagers (Gurnett, et al., 1993). The plasma frequency behind the reflected and the transmitted shock can be responsible for the 2- and 3-kHz radio emissions, respectively. The two bands of radio noise are emitted from sources on both sides of the heliopause starting at about the same time. If the emissions picked up by Voyager are due to f(sub P) radiation, then the heliopause is located at R is approximately = 130 AU. If the emissions are due to 2f(sub P) radiation, the R is approximately 150 AU. Because the relative speed of the interstellar plasma with respect to the sun appears to be sub-Alfvenic, it is very unlikely there is a fast-mode bow shock of the heliosphere.

Whang, Y. C.↗

Experimental Space Shuttle Orbiter Studies to Acquire Data for Code and Flight Heating Model Validation

In an experimental study to obtain detailed heating data over the Space Shuttle Orbiter, CUBRC has completed an extensive matrix of experiments using three distinct models and two unique hypervelocity wind tunnel facilities. This detailed data will be employed to assess heating augmentation due to boundary layer transition on the Orbiter wing leading edge and wind side acreage with comparisons to computational methods and flight data obtained during the Orbiter Entry Boundary Layer Flight Experiment and HYTHIRM during STS-119 reentry. These comparisons will facilitate critical updates to be made to the engineering tools employed to make assessments about natural and tripped boundary layer transition during Orbiter reentry. To achieve the goals of this study data was obtained over a range of Mach numbers from 10 to 18, with flight scaled Reynolds numbers and model attitudes representing key points on the Orbiter reentry trajectory. The first of these studies were performed as an integral part of Return to Flight activities following the accident that occurred during the reentry of the Space Shuttle Columbia (STS-107) in February of 2003. This accident was caused by debris, which originated from the foam covering the external tank bipod fitting ramps, striking and damaging critical wing leading edge heating tiles that reside in the Orbiter bow shock/wing interaction region. During investigation of the accident aeroheating team members discovered that only a limited amount of experimental wing leading edge data existed in this critical peak heating area and a need arose to acquire a detailed dataset of heating in this region. This new dataset was acquired in three phases consisting of a risk mitigation phase employing a 1.8% scale Orbiter model with special temperature sensitive paint covering the wing leading edge, a 0.9% scale Orbiter model with high resolution thin-film instrumentation in the span direction, and the primary 1.8% scale Orbiter model with detailed thin-film resolution in both the span and chord direction in the area of peak heating. Additional objectives of this first study included: obtaining natural or tripped turbulent wing leading edge heating levels, assessing the effectiveness of protuberances and cavities placed at specified locations on the orbiter over a range of Mach numbers and Reynolds numbers to evaluate and compare to existing engineering and computational tools, obtaining cavity floor heating to aid in the verification of cavity heating correlations, acquiring control surface deflection heating data on both the main body flap and elevons, and obtain high speed schlieren videos of the interaction of the orbiter nose bow shock with the wing leading edge. To support these objectives, the stainless steel 1.8% scale orbiter model in addition to the sensors on the wing leading edge was instrumented down the windward centerline, over the wing acreage on the port side, and painted with temperature sensitive paint on the starboard side wing acreage. In all, the stainless steel 1.8% scale Orbiter model was instrumented with over three-hundred highly sensitive thin-film heating sensors, two-hundred of which were located in the wing leading edge shock interaction region. Further experimental studies will also be performed following the successful acquisition of flight data during the Orbiter Entry Boundary Layer Flight Experiment and HYTHIRM on STS-119 at specific data points simulating flight conditions and geometries. Additional instrumentation and a protuberance matching the layout present during the STS-119 boundary layer transition flight experiment were added with testing performed at Mach number and Reynolds number conditions simulating conditions experienced in flight. In addition to the experimental studies, CUBRC also performed a large amount of CFD analysis to confirm and validate not only the tunnel freestream conditions, but also 3D flows over the orbiter acreage, wing leading edge, and controlurfaces to assess data quality, shock interaction locations, and control surface separation regions. This analysis is a standard part of any experimental program at CUBRC, and this information was of key importance for post-test data quality analysis and understanding particular phenomena seen in the data. All work during this effort was sponsored and paid for by the NASA Space Shuttle Program Office at the Johnson Space Center in Houston, Texas.

Wadhams, T. P.↗

Hybrid simulation of comet Shoemaker-Levy 9 interaction with Jovian bow shock

The interaction of the solar wind with comet Shoemaker-Levy 9 leading to the formation of the cometary magnetosphere and its interaction with the Jovian bow shock is simulated using a one dimensional hybrid code. The mass loading of the solar wind by the cometary ions leads to the formation of a bow shock behind which the plasma density is 2-3/cu cm and the electron temperature is 4 eV. The interaction of this system with the Jovian bow shock yields local enhancements of the magnetic field and the plasma density by factors of 4-5 and the electron temperature by 2-3.

Lipatov, A. S.↗

Spectra of noise and amplified turbulence emanating from shock-turbulence interaction

This work is a small extension of NACA studies of the early fifties that predicted amplification of turbulence on passing through a shock wave (observed for turbulent boundary layers), as well as the generation of intense noise (observed for supersonic jets). The first solved the basic gasdynamics problem of the interaction of an infinite planar shock with a single three-dimensional spectrum component of turbulence (an oblique sinusoidal shear wave). The second developed the comprehensive 3D spectrum analysis necessary to generalize the scenario to the interaction of a shock wave with convected homogeneous turbulence. Numerical calculations were carried out to yield curves (vs. Mach number) of rms sound pressure, temperature fluctuation, and two components of turbulent velocity downstream of the shock, for two cases of preshock turbulence. The present numerical study reproduces these for one case and provides in addition their one-dimensional power spectra (vs. wavenumber of frequency). Ratios of the several postshock spectra to the longitudinal preshock turbulence spectrum (4D) have been computed for a wide range of Mach numbers; curves vs. wavenumber are presented for two scenarios of preshock turbulence: isotropy and axisymmetry, both based on the von Karman 3D spectrum.

Ribner, Herbert S.↗

Motion of shocks through interplanetary streams

A model for the motion of flare-generated shocks through interplanetary streams is presented which illustrates the effects of a stream-shock interaction on the shock strength and geometry. It is a gasdynamic calculation based on Whitham's (1958) method and on an empirical approximation for the relevant characteristics of streams. The results show that the Mach number of a shock can decrease appreciably to near unity in the interaction region ahead of streams and that the interaction of a spherically symmetric shock with a spiral-shaped corotating stream can cause significant distortions of the initial shock-front geometry. The geometry of a shock discussed by Lepping and Chao (1972) is qualitatively explained by this model.

Burlaga, L. F.↗