Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock hydrodynamics”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Accretion Disks Around Binary Black Holes of Unequal Mass: GRMHD Simulations Near Decoupling

We report on simulations in general relativity of magnetized disks onto black hole binaries. We vary the binary mass ratio from 1:1 to 1:10 and evolve the systems when they orbit near the binary disk decoupling radius. We compare (surface) density profiles, accretion rates (relative to a single, non-spinning black hole), variability, effective alpha-stress levels and luminosities as functions of the mass ratio. We treat the disks in two limiting regimes: rapid radiative cooling and no radiative cooling. The magnetic field lines clearly reveal jets emerging from both black hole horizons and merging into one common jet at large distances. The magnetic fields give rise to much stronger shock heating than the pure hydrodynamic flows, completely alter the disk structure, and boost accretion rates and luminosities. Accretion streams near the horizons are among the densest structures; in fact, the 1:10 no-cooling evolution results in a refilling of the cavity. The typical effective temperature in the bulk of the disk is approx. 10(exp5) (M / 10(exp 8)M solar mass (exp -1/4(L/L(sub edd) (exp 1/4K) yielding characteristic thermal frequencies approx. 10 (exp 15) (M /10(exp 8)M solar mass) (exp -1/4(L/L (sub edd) (1+z) (exp -1)Hz. These systems are thus promising targets for many extragalactic optical surveys, such as LSST, WFIRST, and PanSTARRS.

Supermassive↗
NASA NTRS1975

Hydrodynamic effects in the atmosphere of variable stars

Numerical models of variable stars are established, using a nonlinear radiative transfer coupled hydrodynamics code. The variable Eddington method of radiative transfer is used. Comparisons are for models of W Virginis, beta Doradus, and eta Aquilae. From these models it appears that shocks are formed in the atmospheres of classical Cepheids as well as W Virginis stars. In classical Cepheids, with periods from 7 to 10 days, the bumps occurring in the light and velocity curves appear as the result of a compression wave that reflects from the star's center. At the head of the outward going compression wave, shocks form in the atmosphere. Comparisons between the hydrodynamic motions in W Virginis and classical Cepheids are made. The strong shocks in W Virginis do not penetrate into the interior as do the compression waves formed in classical Cepheids. The shocks formed in W Virginis stars cause emission lines, while in classical Cepheids the shocks are weaker.

Davis, C. G., Jr.↗
NASA NTRS1983

Monte Carlo simulation of steady state shock structure including cosmic ray mediation and particle escape

Both hydrodynamic calculations (Drury and Volk, 1981, and Axford et al., 1982) and kinetic simulations imply the existence of thermal subshocks in high-Mach-number cosmic-ray-mediated shocks. The injection efficiency of particles from the thermal background into the diffusive shock-acceleration process is determined in part by the sharpness and compression ratio of these subshocks. Results are reported for a Monte Carlo simulation that includes both the back reaction of accelerated particles on the inflowing plasma, producing a smoothing of the shock transition, and the free escape of particles allowing arbitrarily large overall compression ratios in high-Mach-number steady-state shocks. Energy spectra and estimates of the proportion of thermal ions accelerated to high energy are obtained.

Ellison, D. C.↗
NASA NTRS1998

High-Velocity Ly(Alpha) Emission from SMR 1987A

The high-velocity Ly(Alpha) emission from SN 1987A observed with the Space Telescope Imaging Spectrograph (STIS) evidently comes from a reverse shock formed where the outer envelope of SN 1987A strikes ionized gas inside the inner circumstellar ring. The observations can be explained by a simple kinematic model, in which the Ly(Alpha) emission comes from hydrogen atoms with radial velocity approximately 15,000 km s(exp -1) crossing a reverse shock in the shape of a slightly prolate ellipsoid with equatorial radius 4.8 x 10(exp 17) cm or approximately 80% of the distance to the inner surface of the inner ring. N v double Lambda 1239, 1243 emission, if present, has a net luminosity approximately less than 30% times that of the Ly(Alpha) emission. Future STIS observations should enable us to predict the time of impact with the inner ring and to determine unambiguously whether or not N v emission is present. These observations will offer a unique opportunity to probe the structure of SN 1987A's circumstellar environment and the hydrodynamics and kinetics of very fast shocks.

Michael, Eli↗
NASA NTRS1978

Shock tube and shock wave research; Proceedings of the Eleventh International Symposium, University of Washington, Seattle, Wash., July 11-14, 1977

Papers are presented on the applications of shock-wave technology to the study of hydrodynamics, the use of the pressure-wave machine for charging diesel engines, and measurements of the heat-transfer rate in gas-turbine components. Consideration is given to shock propagation along 90-degree bends, the explosive dissemination of liquids, and rotational and vibrational relaxation behind weak shock waves in water vapor. Shock phenomena associated with expansion flows are described and stratospheric-related research using the shock tube is outlined. Attention is given to shock-wave ignition of magnesium powders, Mach reflection and boundary layers, and transition in the shock-induced unsteady boundary layer on a flat plate. Shock-tube measurements of induction and post-induction rates for low-Btu gas mixtures are presented and shock-initiated ignition in COS-N2O-Ar mixtures is described. Cluster growth rates in supersaturated lead vapor are presented and a study of laser-induced plasma motion in a solenoidal magnetic field is reviewed.

Ahlborn, B.↗
NASA NTRS1992

Shock breakout in SN 1987A

Detailed models for the UV burst from SN 1987A are calculated using a two-temperature radiation hydrodynamics code. The sequence of events which occur during shock breakout is described, including the formation of a radiative precursor, the disappearance of the shock, and the growth of a thin, dense shell above the photosphere which gives rise to a new viscous, isothermal shock. Continuum spectra are calculated under the assumption that the emergent radiation is a dilute blackbody with a color temperature equal to the radiation temperature at the thermalization depth. This color temperature is 2-3 times higher than the effective temperature. Light-travel-time corrections have been included in the light curves and the spectra. Quantities most relevant to accurate modeling of the fluorescing ring around SN 1987A are highlighted.

Ensman, Lisa↗
NASA NTRS1984

On the theory of cosmic-ray-mediated shocks with variable compression ratio

Cosmic-ray-mediated shocks may accelerate enough cosmic rays to high enough energies that they escape the shock, carrying an appreciable amount of energy before being convected to downstream infinity. Under such conditions, it is noted, the overall compression ratio cannot be determined from the conservation equations as in conventional hydrodynamic treatments, and the standard equations for shock acceleration admit arbitrarily high compression ratios. A procedure is outlined for obtaining the structure of high Mach number, cosmic-ray-mediated shocks, including their overall compresion ratio, around a low Mach number viscous subshock. Analytic solutions are obtained by quardrature for an energy-dependent diffusion coefficient in the limit of extreme sensitivity to energy, which, unlike previous solutions, include the finite thermal pressure of the preshock gas.

Eichler, D.↗
NASA NTRS1995

Motion of the termination shock in response to an 11 year variation in the solar wind

A two-dimensional hydrodynamic numerical model was used to study the motion of the termination shock in response to an 11 year variation in the solar wind ram pression. For a total variation in the ram pressue by a factor of 2, a termination shock at 89 AU was found to move inward and outward about +/- 8% of its distance with a typical velocity of 12 km/sec. This movement may be understood in terms of the various time scales associated with the response of the termination shock and heliopause to variations in the solar wind ram pressure.

Karmesin, S. Roy↗
NASA NTRS2010

Observing Supernova 1987A with the Refurbished Hubble Space Telescope

The young remnant of supernova 1987A (SN 1987A) offers an unprecedented glimpse into the hydrodynamics and kinetics of fast astrophysical shocks. We have been monitoring SN 1987A with the Hubble Space Telescope (HST) since it was launched. The recent repair of the Space Telescope Imaging Spectrograph (STIS) allows us to compare observations in 2004, just before its demise, with those in 2010, shortly after its resuscitation by NASA astronauts. We find that the Ly-alpha and H-alpha lines from shock emission continue to brighten, while their maximum velocities continue to decrease. We report evidence for nearly coherent, resonant scattering of Lya photons (to blueshifts approximately -12,000 km /s) from hotspots on the equatorial ring. We also report emission to the red of Ly-alpha that we attribute to N v lambda lambda 1239,1243 Angstrom line emission. These lines are detectable because, unlike hydrogen atoms, N4+ ions emit hundreds of photons before they are ionized. The profiles of the N v lines differ markedly from that of H-alpha. We attribute this to scattering of N4+ ions by magnetic fields in the ionized plasma. Thus, N v emission provides a unique probe of the isotropization zone of the collisionless shock. Observations with the recently installed Cosmic Origins Spectrograph (COS) will enable us to observe the N v lambda lambda 1239,1243 Angstrom line profiles with much higher signal-to-noise ratios than possible with STIS and may reveal lines of other highly ionized species (such as C IVlambda lambda 1548,1551 Angstrom) that will test our explanation for the N v emission

France, Kevin↗
NASA NTRS2004

High-Order Simulation of Non-Linear Oscillations and Shocks in the Solar Atmosphere

The solar atmosphere presents a rich source of highly non-linear magneto-hydrodynamic phenomena: strong gradients and forcing terms result in both large shocks and oscillations. The additional requirements of energy balance and initialization in hydrostatic equilibrium compound the challenge of this problem. A wealth of observational data allows us to check the results of our simulations. The problem of simulating the solar atmosphere provides, in addition to an interesting system in its own right, a challenging testbed for high-order shock-capturing methods. We discuss the challenge of simulating solar atmospheric phenomena, concentrating on various high-order central methods ranging from second to fourth order. Our method is based on the central-upwind scheme of Kurganov, Noelle and Petrova, which we extend to high order via various interpolants. We investigate various initial data for our simulations, corresponding to observed conditions in different regions of the solar surface: the normal quiet sun and sunspots. When non-oscillatory using second- and third-order methods, we are able to reproduce non-trivial observational results. In particular we find a correlation between initial data and both the shock speeds and particle oscillation spectra that match observations in the corresponding regions. When using fourth-order WENO interpolants, we find that while the individual shock profiles at any given time appear non-oscillatory, spurious oscillations appear in the fields after long time integrations. The issue of initialization in hydrostatic equilibrium raises difficult issues. Careful treatment of the gravitational source term can reduce violations of hydrostatic equilibrium, but difficulties remain, primarily due to discontinuities in the piecewise-polynomial reconstructions. In the case of the Euler equations in gravity, the use of high-order methods reduces the violation of hydrostatic equilibrium to a sufficiently low level for the above described results to be obtained. Including various energy flux terms in the equations, however, leads to severe instability when hydrostatic balance is violated. We discuss various strategies to address this problem, including adaptive grid methods.

Bryson, S.↗
NASA NTRS1990

Radiative shocks and nonequilibrium chemistry in the early universe - Galaxy and primordial star formation

The nonequilibrium radiative cooling, recombination, and molecule formation behind steady-state shock waves in primordial composition gas is studied. The hydrodynamical conservation equations are solved as well as equations for nonequilibrium ionization, recombination, and molecule formation and the equation of radiative transfer. The shocked gas is found to cool faster than it can recombine; as a result it is able to form an H2 concentration as high as 0.001 or higher via the formation of H(-) and H2(+) intermediaries due to the enhanced nonequilibrium ionization at 10,000 K. With such an H2 concentration, the gas cools by rotational-vibrational line excitation of H2 molecules to well below the canonical final temperature of 10,000 K for a molecule-free gas without metals. It is shown that, as the level of external ionizing and dissociating radiation flux is increased, the formation of and cooling by H2 molecules can be inhibited and delayed.

Shapiro, Paul R.↗
NASA NTRS1995

Collapse and fragmentation of molecular cloud cores. 2: Collapse induced by stellar shock waves

The standard scenario for low-mass star formation involves 'inside-out' collapse of a dense molecular cloud core following loss of magnetic field support through ambipolar diffusion. However, isotopic anomalies in presolar grains and meteoritical inclusions imply that the collapse of the presolar cloud may have been triggered by a stellar shock wave. This paper explores 'outside-in' collapse, that is, protostellar collapse initiated directly by the compression of quiescent dense cloud cores impacted by relatively slow stellar shock waves. A second-order accurate, gravitational hydrodynamics code has been used to study both the spherically symmetrical and three-dimensional evolution of initially centrally condensed, isothermal, self-gravitating, solar-mass cloud cores that are struck by stellar shock waves with velocities up to 25 km/s and postshock temperatures of 10 to 10,000 K. The models show that such mild shock waves do not completely shred and destroy the cloud, and that the dynamical ram pressure can compress the cloud to the verge of self-gravitational collapse. However, compression caused by a high postshock temperature is a considerably more effective means of inducing collapse. Shock-induced collapse produces high initial mass accretion rates (greater than 10(exp -4) solar mass/yr in a solar-mass cloud) that decline rapidly to much lower values, depending on the presence (approximately 10(exp -6) solar mass/yr) or absence (approximately 10(exp -8) to 10(exp -7) solar mass/yr) of an infinite reservoir of mass. Stellar mass accretion rates approximately 10(exp -7) solar mass/yr have been previously inferred from the luminosities of T Tauri stars; balanced mass accretion (stellar rate = envelope rate) at approximately 10(exp -7) solar mass/yr could then be possible if accretion occurs from a finite mass reservoir. Fluid tracers are used to determine what fraction of the stellar shock material is incorporated into the resulting protostellar object and disk; roughly half the impinging material is injected into the collapsing cloud core when there is a high postshock temperature. The models are consistent with a scenario where an AGB star wind triggered the collapse of the presolar cloud while injecting about 0.01 solar mass of matter derived from the AGB star envelope, as has been separately inferred on the basis of nucleosynthesis calculations.

Boss, Alan P.↗
NASA NTRS1989

Shock-generated X-ray emission in radiatively driven winds - A model for Tau Scorpii

A one-dimensional radiation-hydrodynamics code is used to numerically investigate the structure and evolution of shocks in the winds of hot stars. Results are presented for the specific case of Tau Sco, a well-studied main-sequence B star for which there are X-ray data from the Einstien satellite's Solid State Spectrometer. A phenomenological radiative acceleration term and a mass-loss rate consistent with UV observations, are used to determine the time dependence of the temperatures within and X-ray emission from an isolated shock region. The driving acceleration leads to the formation of a two-component shock zone with 'forward' and 'reverse' shocks, each with their own characteristic temperature. A denser cold region forms between the two shocks, which could potentially account for the presence of narrow absorption features that are observed in the UV P Cygni profiles of many hot stars. The X-ray emission spectra from the shocks in the calculations are in good general agreement with two-temperature model fits to Einstein X-ray observations.

Macfarlane, Joseph J.↗
NASA NTRS1992

Radiative shocks and hydrogen molecules in pregalactic gas - The effects of postshock radiation

Previous calculations of radiative shocks in pregalactic gas are generalized and substantially improved. The hydrodynamical conservations are solved, along with the rate equations for nonequilibrium ionization, recombination, and molecule formation, and the equation of radiative transfer, for steady-state shocks in a gas of primordial compositions. This calculation self-consistently includes the effects of the diffuse postshock emission as well as a possible external radiation flux on the postshock flow and the preshock ionization levels. A previous result is confirmed that the shocked gas cools faster than it can recombine and, as a result, is able to form an H2 concentration as high as 0.003 via the formation of H(-) and H2(+) intermediaries due to the enhanced nonequilibrium ionization fraction at 10,000 K. It is shown that, as the level of the external radiation flux is increased, the formation of and cooling by H2 molecules can be inhibited and delayed. The detailed radiation spectra and ionizing photon number fluxes emergent from these shocks are presented.

Kang, Hyesung↗
NASA NTRS1997

Coupling of Sph and Finite Element Codes for Multi-Layer Orbital Debris Shield Design

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.

Fahrenthold, Eric P.↗
NASA NTRS1975

Shock waves in a beta Doradus model

The accumulation of observational evidence and interpretation supports the conclusion that beta Doradus has running shock waves in its atmosphere. It is the work by Hutchinson with its predictions for the phases of the shock waves that prompted this author to construct a numerical hydrodynamical model for the atmosphere of beta Doradus. A short description of the procedure of the calculation follows. The results for the beta Doradus model contain a comparison with the observations. It is of particular interest to verify that the atmosphere contains multiple shock waves as suggested by Hutchinson for the ultraviolet observations.

Hill, S. J.↗
NASA NTRS1985

The cosmic-ray shock structure problem for relativistic shocks

The time asymptotic behaviour of a relativistic (parallel) shock wave significantly modified by the diffusive acceleration of cosmic-rays is investigated by means of relativistic hydrodynamical equations for both the cosmic-rays and thermal gas. The form of the shock structure equation and the dispersion relation for both long and short wavelength waves in the system are obtained. The dependence of the shock acceleration efficiency on the upstream fluid spped, long wavelength Mach number and the ratio N = P sub co/cP sub co+P sub go)(Psub co and P sub go are the upstream cosmic-ray and thermal gas pressures respectively) are studied.

Webb, G. M.↗