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,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
NASA NTRS1976

A numerical solution of the flow field over a transonic airfoil including strong-shock-induced flow separation

An approximate numerical method has been developed for computing the entire transonic flow field over an airfoil, including the boundary layer, the shock-wave/boundary-layer interaction, and the shock-induced-separation. The inviscid flow, the turbulent boundary layer in the weak interaction zone, and the strong interaction elliptic region are solved separately. The inviscid flow is represented by the solution of the transonic small-disturbance equation, while the boundary layer is solved by the generalized Galerkin method, and the strong interaction region (including the wake) is computed from the time-averaged Navier-Stokes equations. A second-order interaction is computed iteratively using the inviscid pressure distribution as a boundary condition for the boundary layer and Navier-Stokes codes, and the viscous displacement thickness to modify the airfoil geometry for the inviscid computation. The computed pressure distributions and shock-wave profiles are compared with experimental data. Agreement is very good except that the pressure recovery on the trailing edge is greater than the experimental data.

Seginer, A.↗
NASA NTRS1985

Supersonic flow with feeding of energy

The present work discusses the results of some experimental studies on the possibility of attenuating shock waves in a supersonic flow. The shock waves were formed by an external source of electrical energy. An electromechanical method is described that permits partial recovery of the expended energy.

Zaremba, W.↗
NASA NTRS2004

High-Order Shock-Capturing Methods for Modeling Dynamics of the Solar Atmosphere

We use one-dimensional high-order central shock capturing numerical methods to study the response of various model solar atmospheres to forcing at the solar surface. The dynamics of the atmosphere is modeled with the Euler equations in a variable-sized flux tube in the presence of gravity. We study dynamics of the atmosphere suggestive of spicule formation and coronal oscillations. These studies are performed on observationally-derived model atmospheres above the quiet sun and above sunspots. To perform these simulations, we provide a new extension of existing second- and third- order shock-capturing methods to irregular grids. We also solve the problem of numerically maintaining initial hydrostatic balance via the introduction of new variables in the model equations and a careful initialization mechanism. We find several striking results: all model atmospheres respond to a single impulsive perturbation with several strong shock waves consistent with the rebound-shock model. These shock waves lift material and the transition region well into the initial corona, and the sensitivity of this lift to the initial impulse depends non-linearly on the details of the atmosphere model. We also reproduce an observed 3-minute coronal oscillation above sunspots compared to 5-minute oscillations above the quiet sun.

Bryson, Steve↗
DOE OSTI2025

Planar shock-induced bubble collapse and jetting in water captured via X-ray phase contrast imaging

Shock wave–bubble interactions in water manifest rich dynamics driven by a combination of strong pressure and density mismatches. They have a wide variety of applications, including the injection of pharmaceuticals, and through scaling, enable the exploration of various aspects of high-energy-density systems such as inertial confinement fusion. In this work, the interaction between a micrometric nitrogen bubble and a planar shock wave, characterized by a Mach number of M = 1.24 and a peak pressure of p max = 0.57 , is experimentally recorded using ultra-high-speed x-ray phase contrast imaging. Highly resolved radiographs provide access to all phase discontinuities along the beam path, offering quantities such as the time-varying bubble size, the speed of a jet produced during the bubble collapse, and the time evolution of the shock wave front, which are critical benchmark data for numerical scheme validation. This study addresses the lack of well-characterized, repeatable, and high spatiotemporal resolution experiments at negative Atwood numbers by providing shock–bubble visualization and corresponding numerical simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1958

Application of the Method of Coordinate Perturbation to Unsteady Duct Flow

The method of coordinate perturbation is applied to the unsteady flow of a compressible fluid in ducts of variable cross section. Solutions, in the form of perturbation series, are obtained for unsteady flows in ducts for which the logarithmic derivative of area variation with respect to the space coordinate is a function of the 'smallness' parameter of the perturbation series. This technique is applied to the problem of the interaction of a disturbance and a shock wave in a diffuser flow. It is found that, for a special choice of the function describing the disturbance, the path of the shock wave can be expressed in closed form to first order. The method is then applied to the determination of the flow field behind a shock wave moving on a prescribed path in the x,t-plane. Perturbation series solutions for quite general paths are developed. The perturbation series solutions are compared with the more exact solutions obtained by the application of the method of characteristics. The approximate solutions are shown to be in reasonably accurate agreement with the solutions obtained by the method of characteristics.

Himmel, Seymour C.↗
NASA NTRS1989

Interplanetary energetic particle observations of the March 1989 events

The IMP-8 spacecraft placed in an elongated orbit of approximately R(sub E) x R(sub E) orbit around the Earth was the only monitor of the energetic particle environment of the near interplanetary space during the period of the solar particle events associated with the Active Region 5395 in March 1989. Measurements of energetic ion and electron intensities were obtained in a series of channels within the energy range: 0.3 to 440 MeV for photons, 0.6 to 52 MeV/nuc for alpha particles, 0.7 to 3.3 MeV/nuc for nuclei with Z greater than or equal to 3, 3 to 9 MeV/nuc with Z greater than or equal to 20, and 0.2 to 2.5 MeV for electrons. The responses of selected energy channels during the period 5 to 23 March 1989 are displayed. It is clearly noted that the most prominent energetic ion intensity enhancements in that time interval were associated with the interplanetary shock wave of March 13 (07:42 UT) as well as that of March 8 (17:56 UT), which have distinct particle acceleration signatures. These shock waves play a major role in determining the near Earth energetic ion intensities during the above period by accelerating and modulating the ambient solar energetic particle population, which was already present in high intensities in the interplanetary medium due to the superposition of a series of solar flare particle events originating in AR 5395. The differential ion intensities at the lowest energy channel of the CPME experiment, which were associated with the March 13 shock wave, reached the highest level in the life of the IMP-8 spacecraft at this energy. At high energies, the shock associated intensity peak was smaller by less than a factor of 3 than the maxima of solar flare particle intensities from some other major flares, in particular from those with sites well connected to the Earth's magnetic flux tubes.

Sarris, E. T.↗
NASA NTRS1983

Passive shock wave/boundary layer control for transonic airfoil drag reduction

An investigation of the passive shock wave/boundary layer interaction control for reducing the drag in 12 percent thick circular arc and 14 percent thick supercritical airfoils was conducted in a 3 in. x 15.4 in. transonic wind tunnel at transonic Mach numbers. A porous surface with a cavity beneath it was positioned on the area of the airfoils, mounted on the test section bottom wall, where the shock wave occurs. The static pressure distributions over the airfoil, the wake impact pressure data for determining the profile drag, and the Schlieren photographs for porous surface airfoils are presented and compared with the results for solid surface airfoils. With the porous surface the normal shock wave for solid surface was changed to a lambda shock wave system, and the wake impact pressure data indicated an appreciable drag reduction for both airfoils with the porous surface at transonic speeds while causing little or no loss of lift. The effect of porosity and cavity size is investigated and off-design performance is discussed.

Bahi, L.↗
NASA NTRS1984

Porosity effect on supercritical airfoil drag reduction by shock wave/boundary layer control

An investigation of the passive shock wave/boundary layer control for reducing the drag of 14 percent-thick supercritical airfoil was conducted in the 3 in. x 15.4 in. RPI Transonic Wind Tunnel at transonic Mach numbers. Various porous surfaces with a cavity beneath it was positioned on the area of the airfoil, mounted on the test section bottom wall, where the shock wave occurs. The static pressure distributions over the airfoil, the wake impact pressure survey for determining the profile drag and the Schlieren photographs for porous surfaces are presented and compared with the results for solid surface airfoil. With a uniform porosity surface the normal shock wave for solid surface was changed to a lambda shock wave, and the wake impact pressure data indicated an appreciable drag reduction at transonic Mach numbers. For a free stream Mach number of 0.81 the profile drag coefficient for the airfoil top surface with uniform porosity was 46 percent lower than for the solid surface airfoil.

Nagamatsu, H. T.↗
NASA NTRS1995

Experimental investigation of the inlet detector configuration variation in the flow field at Mach 1.9

In recent years, active research has been conducted to study the technological feasibility of supersonic laminar flow control on the wing of the High Speed Civil Transport (HSCT). For this study, the F-16XL has been chosen due to its highly swept crank wing planform that closely resembles the HSCT configurations. During flights, it is discovered that the shock wave generated from the aircraft inlet introduces disturbances on the wing where the data acquisition is conducted. The flow field about a supersonic inlet is characterized by a complex three dimensional pattern of shock waves generated by the geometrical configuration of a deflector and a cowl lip. Hence, in this study, experimental method is employed to investigate the effects of the variation of deflector configuration on the flow field, and consequently, the possibility of diverting the incoming shock-disturbances away from the test section. In the present experiments, a model composed of a simple circular tube with a triangular deflector is designed to study the deflector length and the deflector base width variation in the flow field. Experimental results indicate that the lowest external pressure ratio is observed at the junction where the deflector lip and the inlet cowl lip merge. Also, it is noted that the external pressure ratio, the internal pressure ratio, the coefficient of spillage drag, and the shock standoff distance decrease as the deflector length increases. In addition, the Redefined Total Pressure Recovery Ratio (RTPRR) increases with an increase in the deflector length. Results from the study of the effect of the deflector's base width variation on the flow field indicate that the lowest external pressure ratio is observed at the junction between the inlet cowl lip and the deflector lip. As the base width of the deflector increases, the external pressure ratio at 0 rotation increases, whereas the external pressure ratio at 180 rotation decreases. In addition, the internal pressure ratio and the coefficient of spillage drag decrease as the base width of the deflector increases. However, RTPRR and shock standoff distance increase as the base width increases. In conclusion, as deflector dimensions vary, distinctive patterns in the pressure variation around the inlet deflector are observed. With an increase in the deflector length and base width, the magnitude of shock-disturbances are weakened due to a decrease in the external pressure ratio. Also, as the deflector length and base width increase, a smaller bow shock angle is formed. Therefore, the inlet shock wave formation would be significantly altered, and consequently, shock disturbances on the wing test section can be avoided through appropriately designing the deflector.

Hwang, Kyu C.↗
NASA NTRS1985

Two-dimensional kinetics computer program for orbit transfer vehicles

Future Orbit Transfer Vehicles (OTV) presently under consideration require rocket engines delivering a high specific impulse. This high performance can be obtained with large area ratio thrust chambers using oxygen with either hydrogen or hydrocarbon fuels. In the projected nozzles, the combustion products are expanded to low pressure and temperature levels at high Mach number, a domain which has not been experienced with existing rocket engines. In order to assist in the design and evaluation of OTV engines, modifications have recently been incorporated in the Two-Dimensional-Kinetics (TDK) computer program. These modifications include the treatment of shock waves induced by the wall curvature, and shocks which are attached to the wall. The reflection of a shock wave at the flow axis is also treated in an appropriate manner and multiple shock reflections are allowed. A Boundary Layer Module (BLM) has also been incorporated into the TDK computer program in order to provide an automated treatment of the rigorous JANNAF rocket engine performance prediction procedure.

Nickerson, G. R.↗
DOE OSTI2021

Shock generated vorticity in spark discharges

Spark discharges induce a complex flow field consisting of a shock wave at early times (~ 1 μs), a pair of vortex rings, and a hot gas kernel. The vortex rings entrain ambient gas into the hot gas kernel and control its cooling and expansion. In this work, we investigate the shock wave's contribution in producing the vortex ring vorticity. We analyze high-speed (700 kHz) schlieren images of the shock wave for a range of electrical energies to measure the shock properties and estimate shock velocity and curvature. These measurements are combined with a model to calculate the vorticity generated. The measurements show that the highest vorticity is generated near the peak shock curvature location, and the shock curvature and strength increase with electrical energy deposited. A comparison of the vorticity estimated from the model to vorticity measurements from stereoscopic particle image velocimetry shows the results to be statistically equivalent. Furthermore, this suggests that the shock curvature and velocity contribute to the vortex rings induced by spark discharges.

42 ENGINEERING↗
NASA NTRS2015

Analytic Ballistic Performance Model of Whipple Shields

The dual-wall, Whipple shield is the shield of choice for lightweight, long-duration flight. The shield uses an initial sacrificial wall to initiate fragmentation and melt an impacting threat that expands over a void before hitting a subsequent shield wall of a critical component. The key parameters to this type of shield are the rear wall and its mass which stops the debris, as well as the minimum shock wave strength generated by the threat particle impact of the sacrificial wall and the amount of room that is available for expansion. Ensuring the shock wave strength is sufficiently high to achieve large scale fragmentation/melt of the threat particle enables the expansion of the threat and reduces the momentum flux of the debris on the rear wall. Three key factors in the shock wave strength achieved are the thickness of the sacrificial wall relative to the characteristic dimension of the impacting particle, the density and material cohesion contrast of the sacrificial wall relative to the threat particle and the impact speed. The mass of the rear wall and the sacrificial wall are desirable to minimize for launch costs making it important to have an understanding of the effects of density contrast and impact speed. An analytic model is developed here, to describe the influence of these three key factors. In addition this paper develops a description of a fourth key parameter related to fragmentation and its role in establishing the onset of projectile expansion.

Miller, J. E.↗
NASA NTRS1984

Magnetohydrodynamic modelling of interplanetary disturbances between the sun and earth

A time-dependent, nonplanar, two-dimensional magnetohydrodynamic computer model is used to simulate a series, separately examined, of solar flare-generated shock waves and their subsequent disturbances in interplanetary space between the sun and the earth's magnetosphere. The 'canonical' or ansatz series of shock waves include initial velocities near the sun over the range 500 to 3500 km/s. The ambient solar wind, through which they propagate, is taken to be a steady state homogeneous plasma (that is, independent of heliolongitude) with a representative set of plasma and magnetic field parameters. Complete sets of solar wind plasma and magnetic field parameters are presented and discussed. Particular attention is addressed to the MHD model's ability to address fundamental operational questions vis-a-vis the long-range forecasting of geomagnetic disturbances. These questions are: (1) will a disturbance (such as the present canonical series of solar flare shock waves) produce a magnetospheric and ionospheric disturbance, and, if so, (2) when will it start, (3) how severe will it be, and (4) how long will it last? The model's output is used to compute various solar wind indices of current interest as a demonstration of the model's potential for providing 'answers' to these questions.

Dryer, M.↗
NASA NTRS1978

Radiative shock dynamics. II - Hydrogen continua

The interaction between radiation and a shock wave propagating through a stellar atmosphere is investigated. Departures from local thermodynamic equilibrium (LTE) are permitted in the first two levels of a 10-level hydrogen atom; levels 3-10 are in LTE. A piston moving at constant velocity into the bottom of the atmosphere drives a shock wave. This shock produces precursor radiation that diffuses through the gas well ahead of the shock and causes a mild luminosity flash in the emergent Balmer and free-free radiation when it reaches the surface. The precursor wave deposits a large amount of radiative energy in the outer layers of the atmosphere, initiating a radiation-induced pressure wave. The process of energy transfer from the radiation field to the compression wave is similar to the Eddington valve mechanism which drives stellar pulsations. Material is accelerated outward by the radiation-induced wave; eventually it free-falls inward, hits the quasistationary atmosphere, and forms an accretion shock. The piston driven shock is weakened by radiative energy losses. When it reaches the surface, the shock is invisible in the continuum radiation.

Klein, R. I.↗
NASA NTRS2003

Effects of Pre-Existing Target Structure on the Formation of Large Craters

The shapes of large-scale craters and the mechanics responsible for melt generation are influenced by broad and small-scale structures present in a target prior to impact. For example, well-developed systems of fractures often create craters that appear square in outline, good examples being Meteor Crater, AZ and the square craters of 433 Eros. Pre-broken target material also affects melt generation. Kieffer has shown how the shock wave generated in Coconino sandstone at Meteor crater created reverberations which, in combination with the natural target heterogeneity present, created peaks and troughs in pressure and compressed density as individual grains collided to produce a range of shock mineralogies and melts within neighboring samples. In this study, we further explore how pre-existing target structure influences various aspects of the cratering process. We combine experimental and numerical techniques to explore the connection between the scales of the impact generated shock wave and the pre-existing target structure. We focus on the propagation of shock waves in coarse, granular media, emphasizing its consequences on excavation, crater growth, ejecta production, cratering efficiency, melt generation, and crater shape. As a baseline, we present a first series of results for idealized targets where the particles are all identical in size and possess the same shock impedance. We will also present a few results, whereby we increase the complexities of the target properties by varying the grain size, strength, impedance and frictional properties. In addition, we investigate the origin and implications of reverberations that are created by the presence of physical and chemical heterogeneity in a target.

Barnouin-Jha, O. S.↗
NASA NTRS1994

Pressure-sensing performance of upright cylinders in a Mach 10 boundary-layer

An experimental research program to provide basic knowledge of the pressure-sensing performance of upright, flushported cylinders in a hypersonic boundary layer is described. Three upright cylinders of 0.25-, 0.5- and l.0-in. diameters and a conventional rake were placed in the test section sidewall boundary layer of the 31 Inch Mach 10 Wind Tunnel at NASA Langley Research Center, Hampton, Virginia. Boundary-layer pressures from these cylinders were compared to those measured with a conventional rake. A boundary-layer thickness-to-cylinder-diameter ratio of 8 proved sufficient to accurately measure an overall pressure profile and ascertain the boundary-layer thickness. Effects of Reynolds number, flow angularity, and shock wave impingement on pressure measurement were also investigated. Although Reynolds number effects were negligible at the conditions studied, flow angularity above 10 deg significantly affects the measured pressures. Shock wave impingement was used to investigate orifice-to-orifice pressure crosstalk. No crosstalk was measured. The lower pressure measured above the oblique shock wave impingement showed no influence of the higher pressure generated at the lower port locations.

Johnson, Steven↗
NASA NTRS2000

Reconnection Remnants in the Magnetic Cloud of October 18-19, 1995: A Shock, Monochromatic Wave, Heat Flux Dropout and Energetic Ion Beam

Evidence is presented that the WIND spacecraft observed particle and field signatures on October 18-19, 1995 due to reconnection near the footpoints of a magnetic cloud (i.e., between 1 and 5 solar radii). These signatures include: (1) an internal shock traveling approximately along the axis of the magnetic cloud, (2) a simple compression of the magnetic field consistent with the footpoint magnetic fields being thrust outwards at speeds much greater than the solar wind speed, (3) an electron heat flux dropout occurring within minutes of the shock indicating a topological change resulting from disconnection from the solar surface, (4) a very cold 5 keV proton beam and (5) an associated monochromatic wave. We expect that, given observations of enough magnetic clouds, Wind and other spacecraft will see signatures similar to the ones reported here indicating reconnection. However, these observations require the spacecraft to be fortuitously positioned to observe the passing shock and other signatures and will therefore be associated with only a small fraction of magnetic clouds. Consistent with this, a few magnetic clouds observed by Wind have been found to possess internal shock waves.

Collier, Michael R.↗
NASA NTRS1989

Application of a full potential method to AGARD standard airfoils

One of the most important uses of method that calculate unsteady aerodynamic loads is to predict and analyze the aeroelastic responses of flight vehicles. Currently, methods based on transonic small disturbance potential aerodynamics are the primary tools for aeroelastic analysis. Flow solutions obtained using isentropic potential theory can be highly inaccurate and even multivalued, because they do not model the effects of entropy that is produced when shock waves are in the flow field. From the results that are presented, it is concluded that nonisentropic potential methods more accurately model Euler solutions than do isentropic methods. The primary effects of modeling shock generated entropy are: (1) to eliminate mulitple flow solutions when strong shock waves are in the flow field; and (2) to bring the strengths and locations of computed shock waves into better agreement with those calculated using Euler method and those measured during experiments.

Whitlow, Woodrow, Jr.↗