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 55 records · Page 3

Gamma-rays from supernova remnants and the signatures of diffusive shock acceleration

A nonlinear shock acceleration model which generates non-thermal proton distributions and includes a self-consistent determination of shock hydrodynamics, is considered. Gamma ray spectra are obtained for supernova remnants, allowing for the cessation of acceleration to high energies due to the finite ages and the sizes of the remnants. Gamma ray spectral cutoffs can be observed in the TeV range for reasonable remnant parameters and deviations from power law behavior are found at all energies from 1 MeV to cutoff. Correlated observations by the International Gamma Ray Astrophysics Laboratory, Whipple and other instruments may provide stringent constraints to understanding supernova remnants.

Baring, Matthew G.↗

On the hydrodynamic interaction of shock waves with interstellar clouds. 1: Nonradiative shocks in small clouds

The interstellar medium (ISM) is inhomogeneous, with clouds of various temperatures and densities embedded in a tenuous intercloud medium. Shocks propagating through the ISM can ablate or destroy the clouds, at the same time significantly altering the properties of the intercloud medium. This paper presents a comprehensive numerical study of the simplest case of the interaction between a shock wave and a spherical cloud, in which the shock far from the cloud is steady and planar, and in which radiative losses, thermal conduction, magnetic fields, and gravitational forces are all neglected. As a result, the problem is completely specified by two numbers: the Mach number of the shock, M, and the ratio of the density of the cloud to that of the intercloud medium, Chi. For strong shocks we show that the dependence on M scales out, so the primary independent parameter is Chi. Variations from this simple case are also considered: the potential effect of radiative losses is assessed by calculations in which the ratio of specific heats in the cloud is 1.1 instead of 5/3; the effect of the initial shape of the cloud is studied by using a cylindrical cloud instead of a spherical one; and the role of the initial shock is determined by considering the case of a cloud embedded in a wind. Local adaptive mesh refinement techniques with a second-order, two-fluid, two-dimensional Godunov hydrodynamic scheme are used to address these problems, allowing heretofore unobtainable numerical resolution. Convergence studies to be described in a subsequent paper demonstrate that about 100 zones per cloud radius are needed for accurate results; previous calculations have generally used about a third of this number. The results of the calculations are analyzed in terms of global quantities which provide an overall description of te shocked cloud: the size and shape of the cloud, the mean density, the mean pressure, the mean velocity, the velocity dispersion, and the total circulation.

Klein, Richard I.↗

Cosmic-ray hydrodynamics at relativistic shocks

A test particle description of first-order Fermi acceleration of cosmic rays at plane-parallel relativistic shocks is presented in which the pitch angle distribution is obtained explicitly and used to determine the pressure anisotropy or variable Eddington factor as a function of position. The Eddington factor is in turn used to find the hydrodynamical cosmic-ray sound speed, taking the anisotropy into account. For upstream flows speeds above c/sq rt 3, it is shown that the cosmic-ray sound speed exceeds the shock speed in a region just upstream of the shock and decreases montonically with distance further upstream, reaching the value c/sq rt 3 far from the shock. Thus, the shock succeeds in maintaining an anisotropy in the cosmic-ray distribution which is sufficient to keep the flow subsonic with respect to the cosmic-ray sound speed both downstream and in a finite region upstream. The astrophysical implications of these findings are discussed.

Kirk, J. G.↗

Uncertainty Quantification Enabled by Automatic Differentiation for Hydrodynamic Simulation of Shock‐to‐Detonation Transition in High Explosives

Quantifying the effects of uncertainty in a reactive burn model on the run-to-detonation time in high explosives (HEs) provides a robust methodology for assessing the probability of an HE failing the IHE qualification standard. Moreover, uncertainty quantification helps evaluate whether the model calibration accurately represents data outside the calibration set. This study uses a specialized hydrodynamic simulation code for modeling detonation to determine the run-to-detonation time of the HE PBX 9502 for various impact velocities. To quickly approximate uncertainties in the model, a surrogate was constructed using a Taylor series expansion centered at the mean of the input parameters. To obtain the sensitivities required for constructing the Taylor series, HYP-percomplex Automatic Differentiation (HYPAD) was implemented. HYPAD is a methodology for infusing existing codes with automatic differentiation capabilities by augmenting variables with one or more imaginary units to compute step-size independent partial derivatives. These derivatives are accurate to machine precision with respect to the implemented numerical algorithm, meaning their accuracy reflects that of the underlying method (e.g., integration or discretization schemes). Using reduced order modeling techniques, the mean and standard deviation of the run-to-detonation time of a shock within PBX 9502 were computed for a number of initial impact velocities. A weighted least squares regression was then performed to obtain a best fit curve and prediction interval for the computed statistics. Historical data points from explosively driven wedge tests were utilized to validate the prediction interval, ensuring its reliability in predicting future outcomes. With this prediction interval and a known safety constraint curve, the most probable point of failure and the probability of failure for the HE PBX 9502 were determined.

97 MATHEMATICS AND COMPUTING↗

Steady state and dynamical structure of a cosmic-ray-modified termination shock

A hydrodynamic model is developed for the structure of a cosmic-ray-modified termination shock. The model is based on the two-fluid equations of diffuse shock acceleration (Drury and Volk, 1981). Both the steady state structure of the shock and its interaction with outer heliospheric disturbances are considered. Under the assumption that the solar wind is decelerated by diffusing interstellar cosmic rates, it is shown that the natural state of the termination shock is a gradual deceleration and compression, followed by a discontinuous jump to a downstream state which is dominated by the pressure contribution of the cosmic rays. A representative model is calculated for the steady state which incorporates both interstellar cosmic ray mediation and diffusively accelerated anomalous ions through a proposed thermal leakage mechanism. The interaction of large-scale disturbances with the equilibrium termination shock model is shown to result in some unusual downstream structure, including transmitted shocks and cosmic-ray-modified contact discontinuities. The structure observed may be connected to the 2-kHz outer heliospheric radio emission (Cairns et al., 1992a, b). The time-dependent simulations also demonstrate that interaction with solar wind compressible turbulence (e.g., traveling interplanetary shocks, etc.) could induce the termination shock to continually fluctuate between cosmic-ray-dominated and gas-dynamic states. This fluctuation may represent a partial explanation of the galactic cosmic ray modulation effect and illustrates that the Pioneer and Voyager satellites will encounter an evolving shock whose structure and dynamic properties are strongly influence by the mediation of interstellar and anomalous cosmic rays.

Donohue, D. J.↗

Test evaluation of shock buffering concept for hydrodynamic ram induced by yawing projectile impacting a simulated integral fuel tank

A concept for containing the shock inputs due to hydrodynamic ram caused by an impacting projectile within a fuel cell is discussed. This is to provide a buffering layer of foam, flexible, rigid or a combination thereof, which is sealed from the liquid. A program is described in which this buffering concept was tested. The effectiveness of a novel muzzle-mounted, 'tumble', test device is shown.

Zabel, P. H.↗

Steady flow on to a conveyor belt - Causal viscosity and shear shocks

Some hydrodynamical consequences of the adoption of a causal theory of viscosity are explored. Causality is introduced into the theory by letting the coefficient of viscosity go to zero as the flow velocity approaches a designated propagation speed for viscous signals. Consideration is given to a model of viscosity which has a finite propagation speed of shear information, and it is shown that it produces two kinds of shear shock. A 'pure shear shock' corresponds to a transition from a superviscous to a subviscous state with no discontinuity in the velocity. A 'mixed shear shock' has a shear transition occurring at the same location as a normal adiabatic or radiative shock. A generalized version of the Rankine-Hugoniot conditions for mixed shear shocks is derived, and self-consistent numerical solutions to a model 2D problem in which an axisymmetric radially infalling stream encounters a spinning star are presented.

Syer, D.↗

Hydrodynamic simulations of shock-driven chemistry in polyimide

Here, the development for equations of state both for polyimide and its reaction products is presented along with hydrodynamic simulations linking the two EOS through an Arrhenius rate law. The equations of state compare favorably to available data, and the hydrodynamic simulations are able to qualitatively reproduce many of the features seen in both legacy shock data and more modern embedded gauge data. However, quantitative agreement is not achieved, suggesting that new rate laws are needed to fully capture the effect of the reaction over a wide range of conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

X-ray clusters in a cold dark matter + lambda universe: A direct, large-scale, high-resolution, hydrodynamic simulation

A new, three-dimensional, shock-capturing, hydrodynamic code is utilized to determine the distribution of hot gas in a cold dark matter (CDM) + lambda model universe. Periodic boundary conditions are assumed: a box with size 85/h Mpc, having cell size 0.31/h Mpc, is followed in a simulation with 270(exp 3) = 10(exp 7.3) cells. We adopt omega = 0.45, lambda = 0.55, h identically equal to H/100 km/s/Mpc = 0.6, and then, from the cosmic background explorer (COBE) and light element nucleosynthesis, sigma(sub 8) = 0.77, omega(sub b) = 0.043. We identify the X-ray emitting clusters in the simulation box, compute the luminosity function at several wavelength bands, the temperature function and estimated sizes, as well as the evolution of these quantities with redshift. This open model succeeds in matching local observations of clusters in contrast to the standard omega = 1, CDM model, which fails. It predicts an order of magnitude decline in the number density of bright (h nu = 2-10 keV) clusters from z = 0 to z = 2 in contrast to a slight increase in the number density for standard omega = 1, CDM model. This COBE-normalized CDM + lambda model produces approximately the same number of X-ray clusters having L(sub x) greater than 10(exp 43) erg/s as observed. The background radiation field at 1 keV due to clusters is approximately the observed background which, after correction for numerical effects, again indicates that the model is consistent with observations.

Cen, Renyue↗

X-Rays from the Explosion Site: Fifteen Years of Light Curves of SN 1993J

We present a comprehensive analysis of the X-ray light curves of SN 1993J in a nearby galaxy M81. This is the only supernova other than SN 1987A, which is so extensively followed in the X-ray bands. Here we report on SN 1993J observations with the Chandra in the year 2005 and 2008, and Swift observations in 2005, 2006 and 2008. We combined these observations with all available archival data of SN 1993J, which includes ROSAT, ASCA, Chandra, and XMM-Newton, observations from 1993 April to 2006 August. In this paper we report the X-ray light curves of SN 1993J, extending up to fifteen years, in the soft (0.3-2.4 keV), hard (2-8 keV) and combined (0.3-8 keV) bands. The hard and soft-band fluxes decline at different rates initially, but after about 5 years they both undergo a t(sup -1) decline. The soft X-rays, which are initially low, start dominating after a few hundred days. We interpret that most of the emission below 8 keV is coming from the reverse shock which is radiative initially for around first 1000-2000 days and then turn into adiabatic shock. Our hydrodynamic simulation also confirms the reverse shock origin of the observed light curves. We also compare the Ha line luminosity of SN 1993J with its X-ray light curve and note that the Ha line luminosity has a fairly high fraction of the X-ray emission, indicating presence of clumps in the emitting plasma.

Chandra, Poonam↗

Morphology of shocked lateral outflows in colliding hydrodynamic flows

Supersonic interacting flows occurring in phenomena, such as protostellar jets, give rise to strong shocks and have been demonstrated in several laboratory experiments. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in three dimensions. Here, we introduce variations in the flow parameters of density, velocity, and cross-sectional radius of the colliding flows in order to study the propagation and conical shape of the bow shock formed by collisions between two, not necessarily symmetric, hypersonic flows. We find that the motion of the interaction region is driven by imbalances in ram pressure between the two flows, while the conical structure of the bow shock is a result of shocked lateral outflows being deflected from the horizontal when the flows are of differing cross sections.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hydrodynamical constraints on cosmic-ray acceleration in relativistic shocks

A two-fluid hydrodynamical model governing the transport of cosmic rays in a relativistically moving background plasma is developed. The equations are used to discuss the time asymptotic structure of a relativistic, plane-parallel shock wave significantly modified by first-order Fermi acceleration of cosmic rays. The model allows for an anisotropic cosmic-ray pressure tensor with pressure components parallel and perpendicular to the shock normal. It is surmised that the well-known energy threshold problem for diffusive shock acceleration that downstream particles have a large enough velocity to overtake the shock and subsequently scatter in the upstream medium before returning to the downstream medium may be expressed hydrodynamically by the condition V less than V(cr), where V is the fluid velocity relative to the shock and V(cr) is the relativistic cosmic-ray sound speed. Astrophysical implications of the results are briefly discussed.

Webb, G. M.↗

Properties of large scale plasma flow during the early stage of the plasmaspheric refilling

The objective is to better characterize the macroscopic properties of the interhemisphere plasma flow by solving a more complete set of hydrodynamic equations than that solved previously. Specifically, the ion continuity, momentum and energy equations were solved for the plasma flow along the closed magnetic field lines. During the initial stage of the supersonic outflow in the equatorial region, the ions cool substantially. Using the hydrodynamic model for the large-scale plasma flow, the dynamics of shocks was examined which form in the geomagnetic flux tubes during the early stages of refilling. These shocks are more like those forming in neutral gases than the electrostatic shocks driven by microinstabilities involving ion-ion interaction. Therefore, the shocks seen in the hydrodynamic model are termed as hydrodynamic shocks. Such shocks are generally unsteady and therefore the usual shock jump conditions given by Rankine-Hugoniot relations are not strictly applicable to them. The density, flow velocity and temperature structures associated with the shocks are examined for both asymmetrical and symmetrical flows. In the asymmetrical flow the outflow from one of two conjugate ionospheres is dominant. On the other hand, in the symmetrical case outflows from the two ionospheric sources are identical. Both cases are treated by a two-stream model. In the late type of flow, the early-time refilling shows a relaxation type of oscillation, which is driven by the large-scale interactions between the two identical streams. After this early stage, the resulting temperature structure shows some interesting features. In the equatorial region the streams are isothermal, but in the off-equatorial regions the streams have quite different temperatures, and also densities and flow velocities. The dense and slow stream is found to be warmer than the low-density fast stream. In the late stage of refilling, the temperature is found to steadily increase from the conjugate ionospheres towards the equator; the equatorial temperature is found to be as high as about 8000 K compared to the ionospheric temperature of 3600 K.

Singh, Nagendra↗