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.
A brief overview of the physics thought to determine the properties of the heliospheric termination shock is given and the central role of hot interstellar pick-up ions and cosmic rays is emphasized.
Supernova 1987 A has given us an unprecedented view of the evolution of the explosion debris and its interaction with circumstellar matter. The outer supernova debris, now expanding with velocities approx.8000 km/s, encountered the relatively dense circumstellar ring formed by presupernova mass loss in the early 1990s. The shock interaction is manifested by UV-optical "hotspots", an expanding X-ray ring, an expanding ring of knotty non-thermal radio emission, and a ring of thermal IR emission from silicate dust Recent ultraviolet observations of the emissions from the reverse shock and the ring with the HST/COS reveal new details about the shock interaction. Lyman alpha emission from the reverse shock is much stronger than H alpha and they have different emission morphologies, pointing to different emission mechanisms. The reverse shock was detected for the first time in C IV 1550. The N V to C IV brightness ratio indicates the N/C abundance ratio in the expanding debris is about 100X solar, about 3X N/C in the inner ring.
Supernova 1987 A is the brightest and nearest supernova observed since Kepler's SN1604, and is the only one close enough to resolve and directly observe the temporal growth of the ejecta. Over the past 25 years, intensive observations across the electromagnetic spectrum with observatories on the ground (Australia Telescope Compact Array, Gemini-S, Magellan, VLT) and in space (IUE, KAO, CGRO, Hubble, Chandra, Spitzer, Herschel) have given us an unprecedented view of the evolution of the debris of the supernova and of its shock interaction with circumstellar matter. The outer supernova debris, now expanding with velocities -8000 km/s, encountered the relatively dense circumstellar ring formed by presupernova mass loss starting in 1994. The resulting shock interaction has been manifested by: rapidly brightening UV-optical "hotspots", an expanding X-ray ring. an expanding ring of knotty non-thermal radio emission, and a ring of thermal IR emission from silicate dust. The recent evolution of these emissions reveal new details about the shock interaction, circumstellar material, and the star that exploded. Certain critical problems about SN 1987 A, such as the still undiscovered compact object formed in the explosion and the structure of the central debris, require the capabilities of JWST.
Explore the source record for details and available documents.
The paper analyzes the structure of fast shocks incident upon interstellar gas of ambient density from 10 to the 7th per cu cm, while focusing on the problems of formation and destruction of molecules and infrared emission in the cooling, neutral post shock gas. It is noted that such fast shocks initially dissociate almost all preexisting molecules. Discussion covers the physical processes which determine the post shock structure between 10 to the 4 and 10 to the 2 K. It is shown that the chemistry of important molecular coolants H2, CO, OH, and H2O, as well as HD and CH, is reduced to a relatively small set of gas phase and grain surface reactions. Also, the chemistry follows the slow conversion of atomic hydrogen into H2, which primarily occurs on grain surfaces. The dependence of this H2 formation rate on grain and gas temperatures is examined and the survival of grains behind fast shocks is discussed. Post shock heating and cooling rates are calculated and an appropriate, analytic, universal cooling function is developed for molecules other than hydrogen which includes opacities from both the dust and the lines.
The history and theory of shock acceleration is reviewed, paying particular attention to theories of parallel shocks which include the backreaction of accelerated particles on the shock structure. The work that computer simulations, both plasma and Monte Carlo, are playing in revealing how thermal ions interact with shocks and how particle acceleration appears to be an inevitable and necessary part of the basic plasma physics that governs collisionless shocks is discussed. Some of the outstanding problems that still confront theorists and observers in this field are described.
High resolution diffraction grating spectra of X-ray emission from massive stars obtained with Chandra and XMM-Newton have revolutionized our understanding of their powerful, radiation-driven winds. Emission line shapes and line ratios provide diagnostics on a number of key wind parameters. Modeling of resolved emission line velocity profiles allows us to derive independent constraints on stellar mass-loss rates, leading to downward revisions of a factor of a few from previous measurements. Line ratios in He-like ions strongly constrain the spatial distribution of Xray emitting plasma, confirming the expectations of radiation hydrodynamic simulations that X-ray emission begins moderately close to the stellar surface and extends throughout the wind. Some outstanding questions remain, including the possibility of large optical depths in resonance lines, which is hinted at by differences in line shapes of resonance and intercombination lines from the same ion. Resonance scattering leads to nontrivial radiative transfer effects, and modeling it allows us to place constraints on shock size, density, and velocity structure
Computational fluid dynamic (CFD) analyses of the University of Michigan (UM) Shock/Boundary-Layer Interaction (SBLI) experiments were performed as an extension of the CFD SBLI Workshop held at the 48th AIAA Aerospace Sciences Meeting in 2010. In particular, the UM Mach 2.75 Glass Tunnel with a semi-spanning 7.75deg wedge was analyzed in attempts to explore key physics pertinent to SBLI's, including thermodynamic and viscous boundary conditions as well as turbulence modeling. Most of the analyses were 3D CFD simulations using the OVERFLOW flow solver, with additional quasi-1D simulations performed with an in house MATLAB code interfacing with the NIST REFPROP code to explore perfect verses non-ideal air. A fundamental exploration pertaining to the effects of particle image velocimetry (PIV) on post-processing data is also shown. Results from the CFD simulations showed an improvement in agreement with experimental data with key contributions including adding a laminar zone upstream of the wedge and the necessity of mimicking PIV particle lag for comparisons. Results from the quasi-1D simulation showed that there was little difference between perfect and non-ideal air for the configuration presented.
Particle-based hydrodynamics models offer distinct advantages over Eulerian and Lagrangian hydrocodes in particular shock physics applications. Particle models are designed to avoid the mesh distortion and state variable diffusion problems which can hinder the effective use of Lagrangian and Eulerian codes respectively. However conventional particle-in-cell and smooth particle hydrodynamics methods employ particles which are actually moving interpolation points. A new particle-based modeling methodology, termed Hamiltonian particle hydrodynamics, was developed by Fahrenthold and Koo (1997) to provide an alternative, fully Lagrangian, energy-based approach to shock physics simulations. This alternative formulation avoids the tensile and boundary instabilities associated with standard smooth particle hydrodynamics formulations and the diffusive grid- to-particle mapping schemes characteristic of particle-in-cell methods. In the work described herein, the method of Fahrenthold and Koo has been extended, by coupling the aforementioned hydrodynamic particle model to a hexahedral finite element based description of the continuum dynamics. The resulting continuum model retains all of the features (including general contact-impact effects) of Hamiltonian particle hydrodynamics, while in addition accounting for tensile strength, plasticity, and damage effects important in the simulation of hypervelocity impact on orbital debris shielding. A three dimensional, vectorized, and autotasked implementation of the extended particle method described here has been coded for application to orbital debris shielding design. Source code for the pre-processor (PREP), analysis code (EXOS), post-processor (POST), and rezoner (ZONE), have been delivered separately, along with a User's Guide describing installation and application of the software.
No abstract available
Standard techniques used to model chemically-reacting flows require an artificial viscosity for stability in the presence of strong shocks. The resulting shock is smeared over at least three computational cells, so that the thickness of the shock is dictated by the structure of the overall mesh and not the shock physics. A gas passing through a strong shock is thrown into a nonequilibrium state and subsequently relaxes down over some finite distance to an equilibrium end state. The artificial smearing of the shock envelops this relaxation zone which causes the chemical kinetics of the flow to be altered. A method is presented which can investigate these issues by following the chemical kinetics and flow kinetics of a gas passing through a fully resolved shock wave at hypersonic Mach numbers. A nonequilibrium chemistry model for air is incorporated into a spectral multidomain Navier-Stokes solution method. Since no artificial viscosity is needed for stability of the multidomain technique, the precise effect of this artifice on the chemical kinetics and relevant flow features can be determined.
Standard techniques used to model chemically-reacting flows require an artificial viscosity for stability in the presence of strong shocks. The resulting shock is smeared over at least three computational cells, so that the thickness of the shock is dictated by the structure of the overall mesh and not the shock physics. A gas passing through a strong shock is thrown into a nonequilibrium state and subsequently relaxes down over some finite distance to an equilibrium end state. The artificial smearing of the shock envelops this relaxation zone which causes the chemical kinetics of the flow to be altered. A method is presented which can investigate these issues by following the chemical kinetics and flow kinetics of a gas passing through a fully resolved shock wave at hypersonic Mach numbers. A nonequilibrium chemistry model for air is incorporated into a spectral multidomain Navier-Stokes solution method. Since no artificial viscosity is needed for stability of the multidomain technique, the precise effect of this artifice on the chemical kinetics and relevant flow features can be determined.
Although collisionless shocks are ubiquitous in astrophysics, certain key aspects of them are not well understood. In particular, the process known as collisionless electron heating, whereby electrons are rapidly energized at the shock front, is one of the main open issues in shock physics. Here, we present the first clear evidence for efficient collisionless electron heating at the reverse shock of Tycho's supernova remnant (SNR), revealed by Fe K diagnostics using high-quality X-ray data obtained by the Suzaku satellite. We detect K beta (3p yields 1s) fluorescence emission from low-ionization Fe ejecta excited by energetic thermal electrons at the reverse shock front, which peaks at a smaller radius than Fe K alpha (2p yields 1s) emission dominated by a relatively highly ionized component. Comparisons with our hydrodynamical simulations imply instantaneous electron heating to a temperature 1000 times higher than expected from Coulomb collisions alone. The unique environment of the reverse shock, which is propagating with a high Mach number into rarefied ejecta with a low magnetic field strength, puts strong constraints on the physical mechanism responsible for this heating and favors a cross-shock potential created by charge deflection at the shock front. Our sensitive observation also reveals that the reverse shock radius of this SNR is about 10% smaller than the previous measurement using the Fe K alpha morphology from the Chandra observations. Since strong Fe K beta fluorescence is expected only from low-ionization plasma where Fe ions still have many 3p electrons, this feature is key to diagnosing the plasma state and distribution of the immediate postshock ejecta in a young SNR.
Observations of the plasma waves associated with collisionless shocks are reviewed, and the understanding of their generation mechanisms and their importance to shock physics are summarized. The emphasis is on waves generated directly at the shock, especially ion acoustic and lower-hybrid-like modes. The observations are discussed in the context of shock structure, with attention given to the distinctions between waves generated in the shock foot and ramp. The behavior of resistive, dispersive, and supercritical quasi-perpendicular shocks is contrasted. Evidence for the operation of various generation mechanisms, including interactions with cross-field currents, gyrating reflected ions, and field-aligned electron beams, are summarized. The various forms of plasma heating which are actually observed are outlined, and the role of the various wave modes in this heating is discussed. Conclusion, it is argued that, while plasma wave turbulence may play a vital role in plasma heating for some special shocks, it is of second-order importance in most cases.
Fourth order difference approximations of initial-boundary value problems for hyperbolic partial differential equations are considered. We use the method of lines approach with both explicit and compact implicit difference operators in space. The explicit operator satisfies an energy estimate leading to strict stability. For the implicit operator we develop boundary conditions and give a complete proof of strong stability using the Laplace transform technique. We also present numerical experiments for the linear advection equation and Burgers' equation with discontinuities in the solution or in its derivative. The first equation is used for modeling contact discontinuities in fluid dynamics, the second one for modeling shocks and rarefaction waves. The time discretization is done with a third order Runge-Kutta TVD method. For solutions with discontinuities in the solution itself we add a filter based on second order viscosity. In case of the non-linear Burger's equation we use a flux splitting technique that results in an energy estimate for certain different approximations, in which case also an entropy condition is fulfilled. In particular we shall demonstrate that the unsplit conservative form produces a non-physical shock instead of the physically correct rarefaction wave. In the numerical experiments we compare our fourth order methods with a standard second order one and with a third order TVD-method. The results show that the fourth order methods are the only ones that give good results for all the considered test problems.
This paper reports on a Workshop on Supernova Remnants and the Physics of Strong Shock Waves hosted by North Carolina State University at Raleigh, North Carolina, September 16-18, 1993. The workshop brought together observers, shock theorists, cosmic-ray specialists, and simulators to address the role supernova remnants can play in furthering our understanding of the complex plasma physics associated with collisionless shocks and particle acceleration. Over fifty scientists presented papers on various aspects of supernova remnants. In lieu of a proceedings volume, we present here a synopsis of the workshop, in the form of brief summaries of each workshop session.
One of the outstanding problems of collisionless shock physics is the microphysical description of electron thermalization and associated wave processes. Observations at the Earth's bow shock have shown various features in the electron velocity distribution such as a loss-cone, temperature anisotropies, and beams before a thermalized, flat-topped distribution is reached downstream. It is generally believed that in addition to the macroscopic fields, which determine the overall mapping of the distribution function, electromagnetic whistler waves and lower hybrid waves are responsible for the pitch-angle scattering and thermalization of the electrons. However, a genuine physical explanation requires analysis of the evolution of the electron distribution in the self-consistent wave fields. Particle simulations are ideally suited to do this, and include processes in which both electrons and ions are involved, such as cross-field streaming instabilities. We show results from l-D and 2-D simulations using an implicit code, which allows time steps large compared to the inverse electron plasma frequency. It is shown that in order to capture the correct propagation characteristics of the waves, a large mass ratio between the protons and electrons is generally necessary. Detailed diagnostics of the waves as well as the evolution of the electron velocity distribution will be presented and will be discussed in the frame work of satellite observations.
Collisionless shocks in astrophysical plasmas are important thermalizers, converting some of the incident flow energy into thermal energy, and non-thermalizers, partitioning that energy in unequal ways to different particle species, sub-populations thereof, and field components. This partition problem, or equivalently the shock equation of state, lies at the heart of shock physics. Here we employ systematically a framework to capture all the incident and downstream energy fluxes at two example traversals of the Earth’s bow shock by the Magnetospheric Multiscale Mission. Here and traditionally such data has to be augmented by information from other spacecraft, e.g., to provide more accurate measurements of the cold solar wind beam. With some care and fortuitous choices, the energy fluxes are constant, including instantaneous measurements through the shock layer. The dominant incident proton ram energy is converted primarily into downstream proton enthalpy flux, the majority of which is actually carried by a small fraction of suprathermal protons. Fluctuations include both real and instrumental effects. Separating these, resolving the solar wind beam, and other considerations point the way to a dedicated mission to solve this energy partition problem across a full range of plasma and shock conditions.