Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock jumps”

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 181 records · Page 10

Ion and electron heating at the low-Mach-number, quasi-parallel bow shock

The study examines the magnetic structure and energy dissipation at low-Mach-number, quasi-parallel collisionless shocks on the basis of observations obtained by the ISEE 1 and 2 spacecraft at seven crossings of the earth's dayside bow shock. All the shocks exhibit a fairly short-scale-length ramp where the principal jumps in the electron temperature, ion temperature, and magnetic field strength take place. Large-amplitude LF transverse waves in the magnetic field are present in and downstream from the ramp, with more modest waves upstream. The amplitude of the downstream waves is typically larger than can be accounted for by shock compression of the upstream waves. The electron heating represents about 6 percent of the dissipated bulk flow energy, consistent with observations in other parameter regimes, and suggesting that the dominant electron heating process is the same. It is suggested that the low-Mach-number shocks are reforming but do not exhibit downstream variability because of the relatively low levels of reflected ions.

Thomsen, M. F.↗

Development of a computer code for calculating the steady super/hypersonic inviscid flow around real configurations. Volume 1: Computational technique

A numerical procedure has been developed to compute the inviscid super/hypersonic flow field about complex vehicle geometries accurately and efficiently. A second order accurate finite difference scheme is used to integrate the three dimensional Euler equations in regions of continuous flow, while all shock waves are computed as discontinuities via the Rankine Hugoniot jump conditions. Conformal mappings are used to develop a computational grid. The effects of blunt nose entropy layers are computed in detail. Real gas effects for equilibrium air are included using curve fits of Mollier charts. Typical calculated results for shuttle orbiter, hypersonic transport, and supersonic aircraft configurations are included to demonstrate the usefulness of this tool.

Marconi, F.↗

Pioneer Venus plasma observations of the solar wind-Venus interaction

Collection of data from the Ames plasma analyzer on the Pioneer Venus orbiter has permitted long-term measurements of the interaction of the solar wind with Venus. The paper presents a mapping of the ionosheath flow field, plasma measurements in the distant ionosheath and near the distant plasma cavity, and a summary of observations of jumps in the solar wind proton parameters across Venus' bow shock. Also, the apparent detection of ionospheric O(+) accelerated up to solar wind speeds downstream in Venus' ionosheath is discussed.

Mihalov, J. D.↗

Unusual Bow Shock and Magnetopause Locations During May 11-12, 1999

On May 11 and 12, 1999, when the solar wind proton density dropped to only about 3% of its nominal value, WIND and IMP 8 (and possibly Lunar Prospector) encountered the Earth's bow shock at unusually large distances. These observations, at 57-63 Re, correspond to the most distant crossings of the bow shock ever recorded. This paper will present observations from Interball-1 and Geotail, showing the bow shock near its nominal position, and from IMP 8 and WIND (and perhaps Lunar Prospector) crossing the bow shock at increasing radial distances. The bow shock is extremely weak by the time of its observation by WIND. However, a clear jump in both the magnetic field strength and solar wind bulk velocity is quite apparent, although with no significant changes in solar wind proton density. The nature of this dissipating shock will be discussed. Also, interball-1 has observed the magnetopause at an unusually large standoff distance. The observed motion of the magnetopause and bow shock will be compared with model predictions. This set of observations provide an excellent test case for the behavior of gas dynamic and MHD models under extreme conditions.

Szabo, Adam↗

Gas Density Discontinuities in Merging Clusters

In a number of galaxy clusters, Chandra has discovered a new phenomenon, the sharp gas density edges. Depending on the sign of the temperature jump across the edge, these features may either be bow shocks or cold fronts. While the merger origin of the bow shocks is unambiguous, what causes cold fronts is not entirely clear, as they are observed both in mergers and in relaxed clusters. The purpose of the study of A3376, an interesting cluster with density edges, is to understand the origin of cold fronts and to look for possible shocks. This cluster also has a radio halo, and comparison of an X-ray temperature map and radio image may shed light on the nature of the cluster radio halos, which is our secondary goal. The XMM data for A3376 have been obtained, processed and cleaned. Most of the preliminary analysis, including deriving a temperature map, is complete. We also have Chandra data for this cluster, and intend to compare and combine the two datasets. We are well underway toward determining the dynamical state of A3376 and the nature of its gas density discontinuities.

Mushotzky, Richard↗

Gas Density Discontinuities in Merging Clusters

Chandra has discovered a new phenomenon in galaxy clusters, the sharp gas density edges. Depending on the sign of the temperature jump across the edge, these features may either be bow shocks or cold fronts. While bow shocks obviously are driven by merging sub-clusters, what causes cold fronts is not entirely clear, as they are observed both in mergers and in relaxed clusters. The purpose of the XMM study of A3376, an interesting cluster with density edges, is to understand the origin of cold fronts and to look for possible shocks. The XMM data for A3376 have been mostly analyzed (the X-ray edge turned out to be a cold front). Preliminary results have been shown at a conference and a paper is in preparation. We also have Chandra data for this cluster, and are comparing and combining the two datasets. In the course of analyzing the X-ray data for this cluster as well as several others, it has become apparent that we need the help of hydrodynamic simulations to study the precise mechanism by which cold fronts are formed, the main goal of the present project. A postdoc (Yago Ascasibar) is currently running SPH simulations of an idealized sub- cluster merger. These advanced simulations are nearing completion and two papers with their results are in preparation.

Mushotzky, Richard↗

Optimizing ablator thickness for laser shock experiments

In laser shock experiments, a well-defined, flat-top shock wave at the ablator/sample interface is important for accurately probing material response under uniaxial strain compression. However, the relationship between the ablator thickness and the resulting shock wave characteristics is insufficiently understood, limiting the ability to design optimal experiments. To address this need, we conducted a systematic experimental study using a 100 J laser to shock-compress polyimide ablators to peak stresses ranging from 20.4 to 111.6 GPa. Laser interferometry diagnostics measured the transmitted wave profiles at the ablator/sample interface, consistently showing a single jump followed by a constant peak state before the arrival of release waves. Here, by analyzing shock transit time, flat-top duration, and stress, our results establish a framework for selecting ablator thickness to maximize the flat-top duration, improving the precision and reproducibility of laser shock experiments.

design of experiments↗

Development of a computer code for calculating the steady super/hypersonic inviscid flow around real configurations. Volume 2: Code description

A numerical procedure was developed to compute the inviscid super/hypersonic flow field about complex vehicle geometries accurately and efficiently. A second-order accurate finite difference scheme is used to integrate the three-dimensional Euler equations in regions of continuous flow, while all shock waves are computed as discontinuities via the Rankine-Hugoniot jump conditions. Conformal mappings are used to develop a computational grid. The effects of blunt nose entropy layers are computed in detail. Real gas effects for equilibrium air are included using curve fits of Mollier charts. Typical calculated results for shuttle orbiter, hypersonic transport, and supersonic aircraft configurations are included to demonstrate the usefulness of this tool.

Marconi, F.↗

An improved supersonic, three-dimensional, external, inviscid flow field code

A numerical procedure was developed to compute the inviscid super/hypersonic flow fields about complex vehicle geometries accurately and efficiently. A second-order accurate finite difference scheme is used to integrate the three-dimensional Euler equations in regions of continuous flow, while all shock waves are computed as discontinuities via the Rankine-Hugoniot jump conditions. Conformal mappings are used to develop a computational grid. The effects for equilibrium air are included using curve fits of Mollier charts. This report deals only with modifications to these procedures in four specific areas: inlet mass ingestion, subsonic axial Mach number, improved conformal mappings, and vehicles flying at yaw. In each area both the modifications to the computational procedures and computer code are discussed.

Marconi, F.↗

Direct simulation of transitional flow for hypersonic reentry conditions

This paper presents results of flowfield calculations for typical hypersonic reentry conditions encountered by the nose region of the Space Shuttle Orbiter. Most of the transitional flow regime is covered by the altitude range of 150 to 92 km. Calculations were made with the Direct Simulation Monte Carlo (DSMC) method that accounts for translational, rotational, vibrational, and chemical nonequilibrium effects. Comparison of the DSMC heating results with both Shuttle flight data and continuum predictions showed good agreement at the lowest altitude considered. However, as the altitude increased, the continuum predictions, which did not include slip effects, departed rapidly from the DSMC results by overpredicting both heating and drag. The results demonstrate the effects of rarefaction on the shock and the shock layer, along with the extent of the slip and temperature jump at the surface. Also, the sensitivity of the flow structure to the gas-surface interaction model, thermal accommodation, and surface catalysis are studied.

Moss, J. N.↗

A Spectacular Radio Flare from XRF 050416a at 40 Days and Implications for the Nature of X-Ray Flashes

We present detailed optical, near-infrared, and radio observations of the X-ray flash 050416a obtained with Palomar and Siding Springs Observatories as well as HST and the VLA, placing this event among the best-studied X-ray flashes to date. In addition, we present an optical spectrum from Keck LRIS from which we measure the redshift of the burst, Z=0.6528. At this redshift the isotropic-equivalent prompt energy release was about 10(exp 51) erg, and using a standard afterglow synchrotron model we find that the blastwave kinetic energy is a factor of 10 larger, E-K,iso approximately equals 10 (exp 52) erg. The lack of an observed jet break to t - 20 days indicates that the opening angle is larger than 7 deg and the total beaming-corrected relativistic energy is larger than 10 exp (50) erg. We further show that the burst produced a strong radio flare at t is similar to 40 days accompanied by an observed flattening in the X-ray band which we attribute to an abrupt circumburst density jump or an episode of energy injection (either from a refreshed shock or off-axis ejecta). Late-time observations with HST show evidence for an associated supernova with peak optical luminosity roughly comparable to that of SN 1998bw. Next, we show that the host galaxy of XRF 050416a is actively forming stars at a rate of at least 2 M-solar per year with a luminosity of L-B is similar to 0.5L* and metallicity of Z is similar to 0.2-0.8 Z-solar. Finally, we discuss the nature of XRF 050416a in the context of short-hard gamma-ray bursts and under the framework of off-axis and dirty fireball models for X-ray flashes.

Soderberg, A. M.↗

Narrow Radiative Recombination Continua: A Signature of Ions Crossing the Contact Discontinuity of Astrophysical Shocks

X-rays from planetary nebulae (PNs) are believed to originate from a shock driven into the fast stellar wind (v ~ 1000 kilometers per second) as it collides with an earlier circumstellar slow wind (v ~ 10 kilometers per second). In theory, the shocked fast wind (hot hubble) and the ambient cold nebula can remain separated by magnetic fields along a surface referred to as the contact discontinuity (CD) that inhibits diffusion and heat conduction. The CD region is extremely difficult to probe directly owing to its small size and faint emission. This has largely left the study of CDs, stellar-shocks, and the associated micro-physics in the realm of theory. This paper presents spectroscopic evidence for ions from the hot bubble (kT approximately equal to 100 eV) crossing the CD and penetrating the cold nebular gas (kT approximately equal to 1 eV). Specifically, a narrow radiative recombination continuum (RRC) emission feature is identified in the high resolution X-ray spectrum of the PN BD+30degree3639 indicating bare C VII ions are recombining with cool electrons at kT(sub e) = 1.7 plus or minus 1.3 eV. An upper limit to the flux of the narrow RRC of H-like C VI is obtained as well. The RRCs are interpreted as due to C ions from the hot bubble of BD+30degree3639 crossing the CD into the cold nebula, where they ultimately recombine with its cool electrons. The RRC flux ratio of C VII to C VI constrains the temperature jump across the CD to deltakT greater than 80 eV, providing for the first time direct evidence for the stark temperature disparity between the two sides of an astrophysical CD, and constraining the role of magnetic fields and heat conduction accordingly. Two colliding-wind binaries are noted to have similar RRCs suggesting a temperature jump and CD crossing by ions may be common feature of stellar wind shocks.

Behar, Ehud↗

Evolution of the solar wind structure in the outer heliosphere

Shocks and interaction regions play very important roles in the evolution of large-scale solar wind structure in the outer heliosphere. This study is based on (1) plasma and magnetic field data observed from Voyager and Pioneer spacecraft, and (2) a quantitative magnetohydrodynamic simulation model. Interaction regions bounded by a forward and a reverse shock begin to form near 1 AU at the leading edges of a large-scale stream. The total pressure in the region is greater than the ambient pressure by a factor of ten or more. Large jumps in pressure remain as a prominant feature of the interplanetary structure even as the jumps in flow speed become less visible in the outer heliosphere. The propagation of the forward and reverse shocks widens the dimension of an interaction region. As a result, two interaction regions belonging to neighboring streams coalesce to form a merged interaction region (MIR). Collision and merging of shocks take place during the coalescence process. Two MIRs can themselves merge again at greater heliocentric distances. Simulation results agree well with spacecraft observations, and they explain major restructuring of the solar wind in the outer heliosphere.

Whang, Y. C.↗

Biermann-Battery-Driven Magnetized Collisionless Shock Precursors in Laser-Produced Plasmas

This Letter reports the first complete observation of magnetized collisionless shock precursors formed through the compression of Biermann-battery magnetic fields in laser produced plasmas. At OMEGA, lasers produce a supersonic CH plasma flow which is magnetized with Biermann-battery magnetic fields. The plasma flow collides with an unmagnetized hydrogen gas jet plasma to create a magnetized shock precursor. The situation where the flowing plasma carries the magnetic field is similar to the Venusian bow shock. Imaging 2⁢𝜔 Thomson scattering confirms that the interaction is collisionless and shows density and temperature jumps. Proton radiographs have regions of strong deflections and FLASH magnetohydrodynamic (MHD) simulations show the presence of Biermann fields in the Thomson scattering region. Electrons are accelerated to energies of up to 100 keV in a power-law spectrum. Here, OSIRIS particle-in-cell (PIC) simulations, initialized with measured parameters, show the formation of a magnetized shock precursor and corroborate the experimental observables.

High-energy-density plasmas↗

Large-scale fluctuations in the solar wind at 1 AU - 1978-1982

The large-scale fluctuations in the strength of the magnetic field B and the bulk speed V at 1 AU were investigated on the basis of hourly averages of IMP 8 and ISEE 3 data compiled by Couzens and King (1986). The power spectra of the intermediate frequency fluctuations in B at 1 AU had the form of a power law with exponent 1.92 + or - 0.06 during the 1978-1982 interval. The contributions of jumps to the spectra of the large-scale fluctuations in B and V were studied for the interval from days 226 to 365, 1978, which is representative of other such intervals. Of the jumps in B, 22 percent have the signature of a forward or reverse shock.

Burlaga, L. F.↗

A High-Order Eulerian–Lagrangian Runge–Kutta Finite Volume (EL–RK–FV) Method for Scalar Nonlinear Conservation Laws

Abstract We present a class of high-order Eulerian–Lagrangian Runge–Kutta finite volume methods that can numerically solve Burgers’ equation with shock formations, which could be extended to general scalar conservation laws. Eulerian–Lagrangian (EL) and semi-Lagrangian (SL) methods have recently seen increased development and have become a staple for allowing large time-stepping sizes. Yet, maintaining relatively large time-stepping sizes post shock formation remains quite challenging. Our proposed scheme integrates the partial differential equation on a space-time region partitioned by linear approximations to the characteristics determined by the Rankine–Hugoniot jump condition. We trace the characteristics forward in time and present a merging procedure for the mesh cells to handle intersecting characteristics due to shocks. Following this partitioning, we write the equation in a time-differential form and evolve with Runge–Kutta methods in a method-of-lines fashion. High-resolution methods such as ENO and WENO-AO schemes are used for spatial reconstruction. Extension to higher dimensions is done via dimensional splitting. Numerical experiments demonstrate our scheme’s high-order accuracy and ability to sharply capture post-shock solutions with large time-stepping sizes.

Chen, Jiajie↗

Symmetry Determining Equations of the Rankine-Hugoniot Equations for Variable Velocity Shock Waves

The “constant velocity piston” problem (Fig. 1), also known as the “piston problem,” is a standard model for a one dimensional, in our case linear, symmetric shock wave moving through an inviscid, perfect gas. The model can be divided into two regions - a perturbed section on the left and an unperturbed section on the right - by a moving shock wave moving left to right. Both the perturbed and unperturbed sections, i.e. the shocked and unshocked regions, respectively, obey the Eulerian conservation equations; however, at the exact location of the shock, there is a mathematical discontinuity not satisfied by the Euler equations. To ensure continuity and conservation of certain quantities when crossing between the unshocked and shocked regions, we evoke a series of equations derived from the Eulerian conservation equations, called the Rankine-Hugoniot equations, or “jump” equations as it is often referred to in the literature on the topic. The classical constant-velocity piston problem assumes the piston features a constant driving velocity (among many other willing suspensions of belief required in the pursuit of a first principles equation model); consequent to this assumption is a constant-velocity shock and a constant-velocity shocked flow state. However, using Lie Group Theory (LGT), also known as symmetry analysis, we can attempt to reinterpret the model with a shock wave of variable velocity in time and space. An extension of the model in this way opens up the possibility for obtaining new analytical solutions to the piston problem for certain shock velocity models. In this report, we use LGT to derive the symmetry determining equations (SDEs), whose solutions are Lie groups, which permit analytical solutions. In the future, we can then use the SDEs to define constraint equations on the shock velocity model and what the successive solutions to the Euler equations might be based off such constraints. This report is structured as follows: Section 2 provides a brief derivation of the Rankine-Hugoniot (“jump”) equations; Section 3 gives an overview of Lie group theory; Section 4 derives the SDEs of the jump equations; Section 5 derives the Euler conservation equations for fluids; and Section 6 presents concluding remarks and opportunities for future studies.

42 ENGINEERING↗

Shock system of February 2, 1969

The shock system observed in the solar wind by Pioneer 9 and Ogo 5 on Feb. 2, 1969, consisted of the following major discontinuities: a forward slow shock; a forward fast shock; a tangential discontinuity at which the density dropped sharply and the flow direction changed; a tangential discontinuity at which the magnetic field strength jumped to an unusually high value; two closely spaced tangential discontinuities that bracketed a region of even greater field strength and that fronted a region of very cool, very dense, helium-enriched plasma; a reverse fast shock of low Mach number; and a second reverse fast shock of very low Mach number. The event had aspects of both corotating and flare-induced shock systems; it is suggested that the source of the disturbances was a flare occurring at or near an M region.

Unti, T.↗