Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock wave”

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 1,261 records · Page 70

Evidence of Reynolds number sensitivity in supersonic turbulent shocklets

A local shock wave can destroy the evolution of supersonic combustion in a scramjet. The Reynolds number correlation of the systems of internal shock waves produced in such turbulent free shear layers are presently noted to imply that a direct connection with turbulence intensity exists. By way of this connection, as well as the adjustability of the Reynolds number, opportunities emerge for the control of turbulent shocklets; this will in turn allow control of their influence on supersonic mixing.

Johnson, J. A., III↗

On the theory of the propagation of detonation in gaseous systems

The existing theory of detonation is critically examined. It is shown that the considerations with which the steady value of the velocity of detonation is chosen are not convincing. In connection with the problem of the process of the chemical reaction in a detonation wave, the objections raised against the conceptions of Le Chatelier and Vieille of the 19th century with regard to the ignition of the gas by the shock wave are refuted. On the basis of this concept, it is possible to give a rigorous foundation for the existing method of computing the detonation velocity. The distributions of the temperature, the pressure, and the velocity in the detonation wave front as the chemical reaction proceeds, are considered. On the assumption of the absence of losses, the pure compression of the gas in the shock wave at the start of the chemical reaction develops a temperature that is near the temperature of combustion of the given mixture at constant pressure.

THEORIES - DETONATION WAVES↗

Explosions of infalling comets in Jupiter's atmosphere

In view of the expected collision of comet Shoemaker-Levy 9 (1993e) with Jupiter in 1994 July, we calculate basic properties of the initial interaction for a simplified Jovian atmosphere. The comet is expected to impact Jupiter at 60 km/sec and at an angle of 45 deg to the zenith. The shock wave generated by the bolide should be optically thick once it has penetrated to an atmospheric density approximately 10(exp -6) gr/cm(exp 3), and we calculate the post-shock conditions assuming local thermodynamic equilibrium (LTE) for shock velocities v(sub sh) in the range 10 to 60 km/sec and preshock densities rho(sub a) = 10(exp -6) to 10(exp -2) gr/cm(exp 3). Our shock calculations include molecular hydrogen, atomic hydrogen, ionized hydrogen, neutral helium, and singly ionized helium. Even at the highest shock velocity, the gas is only partially ionized and the postshock temperature rises with preshock density in order to maintain the ionization. The value of the effective shock adiabatic index gamma(sub sh) varies from 1.17 (at low v(sub sh) and rho(sub a) to 1.40 (at high v(sub sh) and rho(sub a). The ablation rate is limited by the radiative flux that reaches the bolide surface. We argue that the ablated gas does not efficiently transfer its kinetic energy to the atmosphere, and it ultimately slows in a similar fashion to the comet material. As the bolide initially falls through the atmosphere, the character of the shock emission changes. At rho(sub a) approximately 10(exp -8) gr/cm (exp 3), the gas is optically thin and we expect line emission; in the optical spectrum, Balmer emission is expected from the shocked atmosphere and low-ionization metal lines from ablated cometary material. At rho(sub a) approximately 10(exp -6) gr/m(exp 3), the shocked gas is optically thick and the shock front near the bolide produces a blackbody spectrum. The temperature is favorable for ultraviolet (1000 to 3000 A) emission and the luminosity may be approximately 5 x 10(exp 23) ergs/sec for approximately 0.6 sec for a bolide 1 km in radius. At rho(sub a) approximately 10(exp -4) gr/cm(exp 3), the bolide has passed below the ultraviolet photosphere. The shock front emits considerable ionizing radiation, but it is absorbed in a narrow preshock region. The bolometric correction for the optical luminosity is large and we expect a 3000 to 8000 A luminosity of approximately 3 x 10(exp 23) ergs/sec for approximately 1 sec. The optical emission is strongly peaked in the vicinity of the bolide. The bolide does have a somewhat less luminous, optically thick trail extending greater than or equal to 10 km, but the radiation is characterized by a temperature of 4000 to 5000 K. From the fragmentation model of Chyba, Thomas, & Zahnle (1993), the bolide deposits most of its kinetic energy at rho(sub a) approximately 10(exp -3) gr/cm(exp 3) and this is the effective explosion site. The shock wave from such an explosion can move up about one density scale height. We examine the breakout of the shock front from the Jovian atmosphere and find that the shock acceleration in the decreasing density region is slow, so that the energy flux in the shock front is small. Higher velocities might be generated by shock acceleration along the channel left by the bolide if the shock motion can occur before the channel closes off as a result of radiative cooling. Hot gas created by the explosion ultimately rises due to buoyancy on a timescale of a minute. The luminosity is highest when the bubble first rises into the optically thin part of the atmosphere and may be approximately 1 x 10(exp 25) ergs/sec in the near-infrared. Roughly 1% of the initial bolide energy may be radiated in this way; the rest of the energy is lost to sound waves from the initial explosion and to work done by the bubble on the surrounding atmosphere.

Chevalier, Roger A.↗

Calculation of Velocities from Explosive Shot Test Fast-Frame Imagery

A cylinder of the experimental HE is detonated inside of an aquarium and sequential fast-frame images are taken of the resulting shock wave. We wish to use the sequence of images to extract quantitative data regarding the detonation velocity and velocity of the shock wave which can then be used to tune the parameters of the equation of state for the experimental HE material.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Computation of sharp-fin-induced shockwave/turbulent boundary layer interactions

Solutions of the Reynolds-averaged Navier-Stokes equations are presented and are compared with a family of experimental results for the three-dimensional interaction of a sharp-fin-induced shock wave with a turbulent boundary layer. The solutions predict most of the essential features of the flow fields for various shock-wave strengths. However, some features of the measured flow fields, such as secondary separation and size of the largest separated zones were not accurately computed. The computed flow fields, aided by particle tracing techniques, display a prominent vortical structure which can be correlated with the observed surface phenomena.

Horstman, C. C.↗

Calculation of Unsteady Transonic Flow over an Airfoil

An implicit finite-difference solver for either the Euler equations or the "thin-layer" Navier-Stokes equations was used to calculate a transonic flow over the NACA 64A010 airfoil pitching about its one-quarter chord. An unsteady automatic grid-generation procedure that will improve significantly the computational efficiency of various unsteady flow problems is described. The calculated results for both inviscid and viscous flows at Much number 0.8 over the airfoil oscillating with reduced frequency referenced to one-half chord, 0.2, are compared with experimental data measured in the Ames 11 x 11-ft Transonic Wind Tunnel. Nonlinear, unsteady effects of the flow on the surface pressure variations, shock-wave excursions, and overall airloads are examined. Good agreements between the results of computations and experiments were obtained. In the shock-wave region, however, the results of the viscous-flow computations showed closer agreement with the experimental data.

Chyu, W. J.↗

Calculation of Unsteady Transonic Flow Over an Airfoil

An implicit finite-difference solver for either the Euler equations or the "thin-layer" Navier-Stokes equations was used to calculate a transonic flow over the NACA 64A010 airfoil pitching about its one-quarter chord. An unsteady automatic grid-generation procedure that will improve significantly the computational efficiency of various unsteady flow problems is described. The calculated results for both inviscid and viscous flows at Mach number 0.8 over the airfoil oscillating with reduced frequency referenced to one-half chord, 0.2, are compared with experimental data measured in the Ames 11 x 11 ft Transonic Wind Tunnel. Nonlinear, unsteady effects of the flow on the surface pressure variations, shock-wave excursions, and overall airloads are examined. Good agreements between the results of computations and experiments were obtained. In the shock-wave region, however, the results of the viscous-flow computations showed closer agreement with the experimental data.

Chyu, W. J.↗

Energization in regions of CIRs unconnected to shocks are probably not the result of cross-field transport

Corotating energetic ion populations (CEIPs) associated with the forward and reverse shocks of corotating interaction regions (CIRs) are observed in CIRs at places where models say are magnetically unconnected to either shock. Such disconnections between CEIPs and shocks are common and have been documented with data from Pioneers 10 and 11 and confirmed with data from Ulysses. They pose a problem for models that account for these CEIPs in terms of ion energization at the shocks followed by ion propagation along field lines. Two possible resolutions to this problem have been suggested: diffusion of the ions across field lines and extension of the ion energization process to regions beyond the shock waves. Here we quantitatively examine the first of these possibilities. We give the Green's function solution to the convection-diffusion equation applied to idealized CIR geometry, with a source at the reverse shock -- the main producer of CEIPs. Two kinds of diffusion are considered: resonant diffusion and stochastic field line diffusion. We find that for resonant diffusion the computed ratio is many orders of magnitudes below the observed ratio. For stochastic field line diffusion, the computed ratio approximately equals the observed ratio if a diffusion coefficient appropriate to the free solar wind is used. It is several orders of magnitude below the observed ratio, however, if a diffusion coefficient appropriate to CIRs is used. We conclude that cross-field diffusion probably does not account for the presence of energetic ions in regions of CIRs that are magnetically unconnected to its shock waves. We suggest that the alternative possibility -- the energetic ions in regions magnetically unconnected to shocks result from an acceleration process that is independent of shocks -- be pursued to the point where quantitative tests can be performed.

Intriligator, Devrie S.↗

Optical Response of Soda-Lime Glass Shocked to 14 GPa

To investigate the potential use of soda-lime glass (SLG) as an optical window in shock wave experiments, laser interferometry (VISAR) was used to determine the window correction and the refractive index in SLG samples subjected to plane shock wave compression. Symmetric plate impact experiments were conducted at peak stresses varying between 4 and 14 GPa and laser interferometry ($λ_0$ = 532$\textit {nm}$ wavelength) measurements were used to monitor the particle velocity histories at the impact surface. The higher stress experiments displayed a constant state prior to stress unloading. In contrast, the lower stress experiments (4.3 and 7.8 GPa) displayed an initial constant state followed by weakly time-dependent response prior to stress unloading. Because the actual particle velocities were independently known from the projectile velocity measurements, the window correction factor and the refractive index changes were determined over the stress range examined. Above 7.5 GPa, our window correction and refractive index results show significant differences with previously published data for the same material. Potential reasons for these differences are indicated.

36 MATERIALS SCIENCE↗

Synthesis of diamonds

Shock wave transmitted through a solid exponential horn generates heat and pressure to convert part of a charge of graphite to diamonds. The shock wave is generated in the apparatus by a complex of magnetic fields and eddy currents.

Estes, M. F.↗

Optimizing infrasound observations for sample return capsule re-entry: Insights from OSIRIS-REx and Hayabusa2

The atmospheric entry of meteoroids presents a rare and unpredictable phenomenon, posing challenges for systematic observation and detailed characterization. Such events are nonetheless critical for advancing understanding of acoustic wave propagation, atmospheric structure, and entry dynamics. In contrast, sample return capsules (SRCs) from space missions follow well constrained re-entry trajectories, enabling planned observations of shock wave generation and propagation under controlled conditions. This study compares two SRC atmospheric entries, Hayabusa2 in 2020 and OSIRIS-REx in 2023, to assess how different infrasound array configurations influence shock wave detection and trajectory validation. Hayabusa2’s re-entry was monitored using a distributed network of 28 portable infrasound sensors across seven arrays in Woomera, Australia, permitting three-dimensional reconstruction of the trajectory and analysis of wave characteristics. For OSIRIS-REx, a compact four-sensor array deployed near Eureka Airport provided trajectory confirmation through arrival-time differences and back-azimuth estimates. Spectral and waveform analyses revealed differences in signal properties associated with variations in entry angle and velocity. The results illustrate both the strengths and the limitations of the deployed array configurations. The dense and distributed arrays during Hayabusa2’s re-entry enabled detailed trajectory reconstruction, whereas the compact array at Eureka primarily provided confirmation of signal coherence and back-azimuth consistency. These case studies highlight that even relatively small arrays, if located close to the predicted ground track, can still capture useful information on arrival direction and timing. Such insights provide practical guidance for planning future observational campaigns of SRC returns and other controlled atmospheric entries, and may also inform approaches to opportunistic observations of meteoroid events and other atmospheric acoustic phenomena.

celestial mechanics: orbit determination↗

Wind-tunnel investigation of surface-pressure fluctuations associated with aircraft buffet

Fluctuating pressures and forces that cause aircraft buffeting have been measured on a semispan rigid-wing model of a typical variable-sweep fighter-type aircraft at transonic speeds. The rms spectral and spatial correlation characteristics of wing fluctuating pressures, fluctuating pressure summations, and structural responses are presented and discussed for a Mach number of 0.85, wing sweep angles of 26 and 72 deg, and angles of attack up to 12 deg. The fluctuating pressure characteristics beneath wing shock waves and leading-edge vortices and in regions of attached and separated flows are presented. Results indicate that: (1) the mean and fluctuating static pressure characteristics are related; (2) a circulation oscillation exists for attached flow conditions below buffet onset; and (3) a significant coupling exists between the wing shock-wave oscillation and the wing first torsional mode when shock-induced separation is present.

Riddle, D. W.↗

The Effect of the Ambient Solar Wind Medium on a CME-driven Shock and the Associated Gradual Solar Energetic Particle Event

We present simulation results of a gradual solar energetic particle (SEP) event detected on 2021 October 9 by multiple spacecraft, including BepiColombo (Bepi) and near-Earth spacecraft such as the Advanced Composition Explorer (ACE). A peculiarity of this event is that the presence of a high-speed stream (HSS) affected the low-energy ion component (≲5 MeV) of the gradual SEP event at both Bepi and ACE, despite the HSS having only a modest solar wind speed increase. Using the EUHFORIA (European Heliospheric FORecasting Information Asset) magnetohydrodynamic model, we replicate the solar wind during the event and the coronal mass ejection (CME) that generated it. We then combine these results with the energetic particle transport model PARADISE (PArticle Radiation Asset Directed at Interplanetary Space Exploration). We find that the structure of the CME-driven shock was affected by the nonuniform solar wind, especially near the HSS, resulting in a shock wave front with strong variations in its properties such as its compression ratio and obliquity. By scaling the emission of energetic particles from the shock to the solar wind compression at the shock, an excellent match between the PARADISE simulation and in situ measurements of ≲5 MeV ions is obtained. Our modeling shows that the intricate intensity variations observed at both ACE and Bepi were influenced by the nonuniform emission of energetic particles from the deformed shock wave and demonstrates the influence of even modest background solar wind structures on the development of SEP events.

Solar energetic particles↗

Planetesimal Impact Vapor Plumes and Nebular Shocks Form Chondritic Mixtures

The origin of chondrules and the chondritic sedimentary rocks that dominate the meteoritic record is a long-standing problem in planetary science. Here, we develop a physical model for the formation of chondritic mixtures as an outcome of vaporizing collisions between planetesimals that were dynamically excited by the growth and migration of planets. We present calculations of nebular shock waves generated by impact vapor plumes and focus on aspects of the plume interaction with the nebular gas and dust that have been neglected in previous studies of impact ejecta. We find that when water dominates the vapor, the plumes are relatively cool. However, the plume expansion is supersonic and can drive strong shock waves in the dusty nebular gas. Portions of these nebular shock fronts initially melt nebular dust, forming chondrules that are coupled to the moving front. As the shock front expands and cools, the chondrules solidify, while the shock front entrains additional dust. Eventually, the plume expansion stalls and then hydrodynamically collapses, turbulently mixing variably processed dust and size-sorted chondrules. For probable impact parameters and nebular conditions during giant planet growth and migration, the impact-generated mixtures have characteristics that span the range observed in chondritic meteorites, providing an environment for rapid formation of chondritic assemblages after chondrule formation. Our impact vapor and nebular shocks model, named the IVANS model, links chondrule formation to the overall context of planet formation and provides a framework for interpreting the detailed chronological and geochemical record contained in chondritic meteorites.

79 ASTRONOMY AND ASTROPHYSICS↗

A study to determine the feasibility of a low sonic boom supersonic transport

A study was made to determine the feasibility of supersonic transport configurations designed to produce a goal sonic boom signature with low overpressure. The results indicate that, in principle, such a concept represents a potentially realistic design approach assuming technology of the 1985 time period. Two sonic boom goals were selected which included: (1) A high speed design that would produce shock waves no stronger than 48 Newtons per square meter (1.0 psf); and an intermediate Mach number (mid-Mach) design that would produce shock waves no stronger than 24 Newtons per square meter. The high speed airplane design was a Mach 2.7 blended arrow wing configuration which was capable of carrying 183 passengers a distance of 7000 km (3780 nmi) while meeting the signature goal. The mid-Mach airplane designed was a Mach 1.5 low arrow wing configuration with a horizontal tail which could carry 180 passengers a distance of 5960 km (3220 nmi).

Kane, E. J.↗

Magma ocean formation due to giant impacts

The thermal effects of giant impacts are studied by estimating the melt volume generated by the initial shock wave and corresponding magma ocean depths. Additionally, the effects of the planet's initial temperature on the generated melt volume are examined. The shock pressure required to completely melt the material is determined using the Hugoniot curve plotted in pressure-entropy space. Once the melting pressure is known, an impact melting model is used to estimate the radial distance melting occurred from the impact site. The melt region's geometry then determines the associated melt volume. The model is also used to estimate the partial melt volume. Magma ocean depths resulting from both excavated and retained melt are calculated, and the melt fraction not excavated during the formation of the crater is estimated. The fraction of a planet melted by the initial shock wave is also estimated using the model.

Tonks, W. B.↗

Wall pressure fluctuations in the reattachment region of a supersonic free shear layer

The primary aim of this research program was to investigate the mechanisms which cause the unsteady wall-pressure fluctuations in shock wave turbulent shear layer interactions. The secondary aim was to find means to reduce the magnitude of the fluctuating pressure loads by controlling the unsteady shock motion. The particular flow under study is the unsteady shock wave interaction formed in the reattachment zone of a separated supersonic flow. Similar flows are encountered in many practical situations, and they are associated with high levels of fluctuating wall pressure. The free shear layer is formed by the flow over a backward facing step, using an existing model, with the base pressure on the step adjusted so that there is no pressure discontinuity at the lip. The shear layer therefore develops in a zero pressure gradient. The primary advantage of this flow configuration is that the reattachment process can be studied in the absence of a separation shock. The mean flow data, and some preliminary hot-wire measurements of the mass-flux fluctuations were made by Baca and Settles, Baca, Williams and Bogdonoff, who showed that the shear layer became self-similar at about 17 delta(sub 0) downstream of the lip, and that it grew at a rate typical of the observed Mach number difference (about 1/3rd the incompressible growth rate). The turbulence measurements were later extended by Hayakawa, Smits and Bogdonoff under NASA Headquarters support.

Smits, Alexander J.↗

The dynamics and control of fluctuating pressure loads in the reattachment region of a supersonic free shear layer

The primary aim is to investigate the mechanisms which cause the unsteady wall-pressure fluctuations in shock wave turbulent shear layer interactions. The secondary aim is to find means to reduce the magnitude of the fluctuating pressure loads by controlling the unsteady shock motion. The particular flow proposed for study is the unsteady shock wave interaction formed in the reattachment zone of a separated supersonic flow. Similar flows are encountered in many practical situations, and they are associated with high levels of fluctuating wall pressure. Wall pressure fluctuations were measured in the reattachment region of the supersonic free shear layer. The free shear layer was formed by the separation of a Mach 2.9 turbulent boundary layer from a backward facing step. Reattachment occurred on a 20 deg ramp. By adjusting the position of the ramp, the base pressure was set equal to the freestream pressure, and the free shear layer formed in the absence of a separation shock. An array of flush-mounted, miniature, high-frequency pressure transducers was used to make multichannel measurements of the fluctuating wall pressure in the vicinity of the reattachment region. Contrary to previous observations of this flow, the reattachment region was found to be highly unsteady, and the pressure fluctuations were found to be significant. The overall behavior of the wall pressure loading is similar in scale and magnitude to the unsteadiness of the wall pressure field in compression ramp flows at the same Mach number. Rayleigh scattering was used to visualize the instantaneous shock structure in the streamwise and spanwise direction. Spanwise wrinkles on the order of half the boundary layer thickness were observed.

Smits, A. J.↗