STRUCTURE OF A PLANE SHOCK LAYER
Discussion of plane shock wave structure computation of shock profiles and application of navier-stokes solution
SEARCH · Engineering Papers
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.
Discussion of plane shock wave structure computation of shock profiles and application of navier-stokes solution
Experiments on shock wave structure were conducted in the 12 inch shock tube facility, and it is shown that this work was motivated by the need for improved kinetics for the modeling of NO production in space shuttle flow fields. The experiments on shock structure involved simultaneous observations of pressure and temperature on the wall of the shock tube during reflection of normal shock waves in inert gases.
A two-dimensional symmetric wedge configuration representative of a single high-speed intake in steady flow was investigated. The analysis presented here is intended as an engineering approach for estimating certain features of the internal shock system. The primary interest here is the prediction of the size and location of the almost-normal shock wave that develops when the leading-edge shocks intersect at angles above a certain critical value that is less than the wedge detachment angle. The almost-normal shock wave is frequently referred to as the 'Mach stem', Parametric studies enabled the sensitivity of the Mach stem height to various flowfield parameters to be examined, thus indicating how accurately these parameters must be measured in a given experiment. Results of these predictions were compared with those of a steady-flow experiment performed at nominal freestream Mach numbers from 2.8 to 5. The predicted stem heights were consistently lower than the mean experimental values, attributable both to experimental uncertainties and to certain simplifying assumptions used in the analysis. Modification of these assumptions to better represent the test environment improved the analytical results.
The flow behind the shock wave formed around objects which fly at hypervelocity behaves differently from that of a perfect gas. Molecules become vibrationally excited, dissociated, and ionized. The hot gas may emit or absorb radiation. When the atoms produced by dissociation reach the wall surface, chemical reactions, including recombination, may occur. The thermochemical phenomena of vibration, dissociation, ionization, surface chemical reaction, and radiation are referred to commonly as high-temperature real-gas phenomena. The phenomena cause changes in the dynamic behavior of the flow and the surface pressure and heat transfer distribution around the object. The character of a real gas is described by the internal degrees of freedom and state of constituent molecules; nitrogen and oxygen for air. The internal energy states, rotation, vibration and electronic, of the molecules are excited and, in the limit, the molecular bonds are exceeded and the gas dissociated into atomic and, possibly, ionic constituents. The process of energy transfer causing excitation, dissociation and recombination is a rate process controlled by particle collisions. Binary, two-body, collisions are sufficient to cause internal excitation, dissociation and ionization while three-body collisions are required to recombine the particles into molecular constituents. If the rates of energy transfer are fast with respect to the local fluid dynamic time scale the gas is in, or nearly in, equilibrium. If the energy transfer rates are very slow the gas can be described as frozen. In all other instances, wherein any of the energy exchange rates are comparable to the local fluid time scale, the gas will be thermally or chemically reacting and out of equilibrium. Real gas thermochemical nonequilibrium processes are important in the determination of aerodynamic heating; both convective (including wall catalytic effects) and radiative heating. To illustrate this we consider the hypervelocity flow over a bluff body typical of an atmospheric entry vehicle or an aerospace transfer vehicle (ASTV). The qualitative aspects of a hypersonic flow field over a bluff body are discussed in two parts, forebody and afterbody, with attention to which particular physical effects must be included in an analysis. This will indicate what type of numerical modeling will be adequate in each region of the flow. A bluff forebody flow field is dominated by the presence of the strong bow shock wave and the consequent heating, and chemical reaction of the gas. At high altitude hypersonic flight conditions the thermal excitation and chemical reaction of the gas occur slowly enough that a significant portion of the flow field is in a state of thermochemical nonequilibrium. A second important effect is the presence of the thick boundary layer along the forebody surface. In this region there are large thermal and chemical species gradients due to the interaction of the gas with the wall. Also at high altitudes the shock wave and the boundary layer may become so thick that they merge; in this case the entire shock layer is dominated by viscous effects.
Dynamic compression of composite materials is of scientific interest because the mechanical mismatch between internal phases challenges continuum theories. Typical assumptions about steady wave propagation and quasi-instantaneous state changes require reexamination along with the need for time-dependent models. To that end, data and models are presented here for the shock compression of an idealized particulate composite. To serve as a generic representative of this material class, a polymer matrix was filled with tungsten particles, ranging from 1 to 50 vol. %. This creates a simple microstructure containing randomly scattered particles with an extreme impedance mismatch to the binding matrix. These materials were parallel plate impact loaded by Al flyers traveling at 1.8–5.0 km/s. Velocimetry provided records of the equilibrium state and the compression wave structure for each case with trends quantified by an empirical fit. The same quantities were also studied as a function of the wave's propagation distance. A homogenized viscoelastic model then made it possible to progress from cataloging trends to making predictions. Starting from a Mie–Gruneisen equation of state, additional time varying terms were added to capture the transient response. After calibration, accurate predictions of the steady wave structure were possible.
The unsteady laminar boundary layer induced by the flow-initiating shock wave passing over a flat plate mounted in a shock tube was theoretically and experimentally studied in terms of heat transfer rates to the plate for shock speeds ranging from 1.695 to 7.34 km/sec. The theory presented by Cook and Chapman for the shock-induced unsteady boundary layer on a plate is reviewed with emphasis on unsteady heat transfer. A method of measuring time-dependent heat-transfer rates using thin-film heat-flux gages and an associated data reduction technique are outlined in detail. Particular consideration is given to heat-flux measurement in short-duration ionized shocktube flows. Experimental unsteady plate heat transfer rates obtained in both air and nitrogen using thin-film heat-flux gages generally agree well with theoretical predictions. The experimental results indicate that the theory continues to predict the unsteady boundary layer behavior after the shock wave leaves the trailing edge of the plate even though the theory is strictly applicable only for the time interval in which the shock remains on the plate.
The rms field strengths of electrostatic and electromagnetic turbulence in the earth's bow shock, measured in the frequency range 20 Hz to 200 kHz with the Imp 6 satellite, are examined. The largest strengths of electrostatic turbulence occur when the upstream electron-to-proton temperature ratio is large and the proton temperature is small. No substantial correlation is found among the rms field strengths of electrostatic turbulence, the Alfven Mach number, the ratio of particle to magnetic-field pressure, and the shock normal angle. These results indicate that the strength of electrostatic turbulence in the bow shock is determined by the kinetic properties of the solar-wind plasma. The largest strengths for electromagnetic turbulence occur when the upstream particle density is large and when the shock normal angle is closer to 90 deg, supporting a previous conclusion that whistler waves comprise the electromagnetic turbulence in the bow shock.
Chordwise pressure distributions were obtained in-flight on the upper and lower surfaces of the F-16XL ship 2 aircraft wing between Mach 1.4 and Mach 2.0. This experiment was conducted to determine the location of shock waves which could compromise or invalidate a follow-on test of a large chord laminar flow control suction panel. On the upper surface, the canopy closure shock crossed an area which would be covered by a proposed laminar flow suction panel. At the laminar flow experiment design Mach number of 1.9, 91 percent of the suction panel area would be forward of the shock. At Mach 1.4, that value reduces to 65 percent. On the lower surface, a shock from the inlet diverter would impinge on the proposed suction panel leading edge. A chordwise plate mounted vertically to deflect shock waves, called a shock fence, was installed between the inlet diverter and the leading edge. This plate was effective in reducing the pressure gradients caused by the inlet shock system.
A multi-point focused laser differential interferometer (FLDI) has been developed to measure density fluctuations at 16 points along a line. A pair of cylindrical lenses on the transmitter side of a conventional single-point FLDI instrument form two closely spaced (≤200 μm), orthogonally polarized, parallel laser lines at the instrument’s focus. On the receiver side of the instrument, the interference of the beams on a 16-element photodiode array results in a single line of measurements. The further addition of a Nomarski prism creates two separate measurement lines, and the addition of a second photodiode array to the instrument enables simultaneous measurements of density fluctuations along the two lines separated by several millimeters. These two lines of measurement can be conveniently oriented at any azimuthal angle relative to the instrument’s optical axis on the measurement plane, coinciding with the instrument’s focus. Two experiments were performed to demonstrate the capabilities of the instrument. In the first experiment, a laser-induced breakdown spark generated a traveling spherical shock wave, and measurements of the resulting density disturbance and wave velocity were obtained. These results were compared to high-speed schlieren images of the shock wave acquired at 400 kHz. In the second experiment, the multi-point FLDI instrument was used to measure density disturbances in the boundary layer of a flat plate in a Mach 6 freestream flow. The measurements were made along two lines, both approximately 6 mm in length, extending from the surface of the plate through the boundary layer. High-speed schlieren images were acquired at 100 kHz during separate wind tunnel runs at matching unit Reynolds numbers to visualize the unsteady boundary layer flow and compare to the FLDI measurements.
The present examination of the response of anisotropic turbulent flows to a shock wave uses incompressible, homogeneous rapid distortion theory and idealizes the shock-wave as 1D compression. Attention is given to the shock's effect on axisymmetric flow and both the normal and oblique cases for shear flow. The oblique angle between the directions of compression and shear significantly affects turbulence evolution. Compression amplifies all components of turbulence intensity.
To accurately model inviscid flow with shock waves using staggered-grid Lagrangian hydrodynamics, artificial viscosity is introduced to convert kinetic energy into internal energy, thereby providing a mechanism to generate the required entropy increase across shocks. In this paper, we propose a new method for constructing an adaptive, artificial viscosity in the context of one-dimensional, staggered-grid Lagrangian hydrodynamics. Our adaptive, artificial viscosity is defined in terms of two parameters that depend on density locally in a neighborhood around a shock, and hence, vary cell-by-cell. Our methodology is based on building a reference set of pre-computed optimal, globally constant artificial viscosity coefficients for a family of isolated shock test problems. For arbitrary flows, the evaluation of the unknown coefficients is automated by first estimating shock intensity locally, and second computing the corresponding adaptive parameter value via interpolation over the data from the pre-computed reference set of optimal, globally constant coefficient values. To illustrate the performance of our new approach, we compare our results against two existing methods. The first method is a limiter-based approach, which relies on estimating velocity gradients of the flow, and the second method is utilized in several commercial codes. Finally, we demonstrate that our new adaptive methodology produces more accurate results for a variety of tests with propagating shock waves, as well as for the aforementioned family of isolated shock problems.
Here solid-solid phase transitions are investigated in Zr thin films using shock compression induced by a short laser pulse (<1 ns). Shock wave profiles are measured at free surfaces for films of different thicknesses (a few micrometers) using chirped-pulse line velocimetry with 10 ps time resolution. Experiments are performed at pressures up to ~50 GPa, which is sufficient to reach the ω and β phases under equilibrium conditions. The shock wave structures are analyzed using a general Lagrangian analysis method, which allows for the calculation of stress–strain paths and assessments of phase transition behavior. In agreement with recent short laser pulse experiments using ultra-fast x-ray diffraction, we do not find any clear evidence of the α–ω transition, though this would be expected from the phase diagram. Instead, we infer a direct transformation to a metastable β-phase at lower shock pressures (<20 GPa) and equilibrium β at higher pressures. Through the velocimetry analysis, we find α–β transformation onset times of less than ~100 ps and completion times of less than ~200 ps.
Relaxation solutions to classical three-dimensional small-disturbance (CSD) theory for transonic flow about lifting swept wings are reported. For such wings, the CSD theory was found to be a poor approximation to the full potential equation in regions of the flow field that are essentially two-dimensional in a plane normal to the sweep direction. The effect of this deficiency on the capture of embedded shock waves in terms of (1) the conditions under which shock waves can exist and (2) the relations they must satisfy when they do exist is emphasized. A modified small-disturbance (MSD) equation, derived by retaining two previously neglected terms, was proposed and shown to be a consistent approximation to the full potential equation over a wider range of sweep angles. The effect of these extra terms is demonstrated by comparing CSD, MSD, and experimental wing surface pressures.
Numerical solutions of the time-averaged Navier-Stokes equations employing a simple eddy-viscosity model have been obtained for three dimensional turbulent flow fields at supersonic speeds. The computer results are compared with a series of experimental test flows describing the interaction of a swept shock wave with a turbulent boundary layer for various shock-wave strengths. Very good agreement is obtained between the computed and experimental surface and flow-field results. The computed flow fields are examined in detail to investigate the physics of this type of flow field. Questions concerning the existence of a vortex and the relationship between converging surface oil streaks and the resulting flow field are addressed.
A change in flow pattern that was observed as the free-stream Mach number was increased in the vicinity of 0.8 was described in NACA Technical Note 1211 by Lindsey, Daley, and Humphreys. The flow on the upper surface behind the leading edge of an airfoil at an angle of attack changed abruptly from detached flow with an extensive region of separation to attached supersonic flow terminated by a shock wave. In the present paper, the consequences of shock-wave - boundary layer interaction are proposed as a factor that may be important in determining the conditions under which the change in flow pattern occurs. Some experimental evidence in support of the importance of this factor is presented.
A brief review is given of the unsteadiness of shock wave/turbulent boundary layer interaction. The focus is on interactions generated by swept and unswept compression ramps, by flares, steps and incident shock waves, by cylinders and blunt fins, and by glancing shock waves. The effects of Mach number, Reynolds number, and separated flow scale are discussed as are the physical causes of the unsteadiness. The implications that the unsteadiness has for interpreting time-average surface and flowfield data, and for comparisons of such experimental data with computation, is also briefly discussed. Finally, some suggestions for future work are given. It is clear that there are large gaps in the data base and that many aspects of such phenomena are poorly understood. Much work remains to be done.
Active separation control, using periodic excitation, was studied experimentally at high Reynolds numbers. The effects of compressibility, mild sweep, location o excitation slot and steady momentum transfer on the efficacy of the method were identified. Tests conducted at chord Reynolds numbers as high as 40 x 10(exp 6) demonstrated that active control using oscillatory flow excitation can effectively delay flow separation from and reattach separated flow to aerodynamic surfaces at flight conditions. The effective frequencies generate one to four vortices over the controlled region at all times, regardless of the Reynolds number. The vortices are initially amplified by the separated shear-layer, and after initiating reattachment, the strength of the vortices decay as they are convected downstream. Large amplitude, low frequency vortices break down to smaller ones upon introduction at the excitation slot. The effects of steady mass transfer were compared to those of periodic excitation. It was found that steady blowing is significantly inferior to periodic excitation in terms o performance benefits and that the response to steady blowing is abrupt, and therefore undesirable from a control point of view. Steady suction and periodic excitation are comparable in effectiveness and both exhibit a gradual response to changes in the magnitude of the control input. The combination of weak steady suction and periodic excitation is extremely effective while the addition of steady blowing could be detrimental. Compressibility effects are weak as long as separation is not caused by a shock-wave/boundary-layer interaction The undesirable effects of the shock-induced separation could be alleviated by the introduction of periodic excitation upstream of the shock wave, inside the region of supersonic flow. The effects of mild sweep were also studied and periodic excitation was found to be very effective in reattaching three-dimensional separated flow. Scaling laws that correlate 2D and 3D controlled flows were tested and verified. Several performance benefits could be gained by applying the method to existing configurations, but it is expected that the full potential of the method can only be realized through the design of new configurations. A comprehensive, fully turbulent, database was generated in order to guide the development, and enable validation, of candidate unsteady CFD design tools.
Active separation control, using periodic excitation, was studied experimentally at high Reynolds numbers. The effects of compressibility, mild sweep, location of excitation slot and steady momentum transfer on the efficacy of the method were identified. Tests conducted at chord Reynolds numbers as high as 40 x 10(exp 6) demonstrated that active control using oscillatory flow excitation can effectively delay flow separation from, and reattach separated flow to aerodynamic surfaces at flight conditions. The effective frequencies generate one to four vortices over the controlled region at all times, regardless of the Reynolds number. The vortices are initially amplified by the separated shear-layer, and after initiating reattachment, the strength of the vortices decay as they are convected downstream. Large amplitude, low frequency vortices break down to smaller ones upon introduction at the excitation slot. The effects of steady mass transfer were compared to those of periodic excitation. It was found that steady blowing is significantly inferior to periodic excitation in terms of performance benefits and that the response to steady blowing is abrupt, and therefore undesirable from a control point of view. Steady suction and periodic excitation are comparable in effectiveness and both exhibit a gradual response to changes in the magnitude of the control input. The combination of weak steady suction and periodic excitation is extremely effective while the addition of steady blowing could be detrimental. Compressibility effects are weak as long as separation is not caused by a shock-wave/boundary-layer interaction. The undesirable effects of the shock-induced separation could be alleviated by the introduction of periodic excitation upstream of the shock wave, inside the region of supersonic flow. The effects of mild sweep were also studied and periodic excitation was found to be very effective in reattaching three-dimensional separated flow. Scaling laws that correlate 2D and 3D controlled flows were tested and verified. Several performance benefits could be gained by applying the method to existing configurations, but it is expected that the full potential of the method can only be realized through the design of new configurations. A comprehensive, fully turbulent, database was generated in order to guide the development, and enable validation, of candidate unsteady CFD design tools.