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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS1981

Heating of the solar chromosphere by ionization pumping

A new theory is proposed to explain the heating of the solar chromosphere, and possibly the corona, by the dissipation of hydrodynamic compression waves. The basis of the dissipative mechanism, here referred to as ionization pumping, is hysteresis caused by irreversible relaxation of the chromospheric medium to ionization equilibrium following pressure perturbations. In the middle chromosphere, where hydrogen is partially ionized, it is shown that ionization pumping will cause strong dissipation of waves whose periods are 200s or less. This could cause heating of the chromosphere sufficient to compensate for the radiative losses. The mechanism retains a high efficiency for waves of arbitrarily small amplitude and, thus, can be more efficient than shock dissipation for small perturbations in pressure. The formation of shocks therefore is not required for the dissipation of waves whose periods are several minutes or less.

Lindsey, C. A.↗
NASA NTRS1999

High-Energy Space Propulsion Based on Magnetized Target Fusion

Magnetized target fusion is an approach in which a magnetized target plasma is compressed inertially by an imploding material wall. A high energy plasma liner may be used to produce the required implosion. The plasma liner is formed by the merging of a number of high momentum plasma jets converging towards the center of a sphere where two compact toroids have been introduced. Preliminary 3-D hydrodynamics modeling results using the SPHINX code of Los Alamos National Laboratory have been very encouraging and confirm earlier theoretical expectations. The concept appears ready for experimental exploration and plans for doing so are being pursued. In this talk, we explore conceptually how this innovative fusion approach could be packaged for space propulsion for interplanetary travel. We discuss the generally generic components of a baseline propulsion concept including the fusion engine, high velocity plasma accelerators, generators of compact toroids using conical theta pinches, magnetic nozzle, neutron absorption blanket, tritium reprocessing system, shock absorber, magnetohydrodynamic generator, capacitor pulsed power system, thermal management system, and micrometeorite shields.

Thio, Y. C. F.↗
DOE OSTI2022

Particle-in-cell modeling of plasma jet merging in the large-Hall-parameter regime

The merging process of magnetized plasma jets with parameters relevant to the plasma-jet-driven magneto-inertial fusion (PJMIF) design and the plasma liner experiment (PLX) is modeled by fully kinetic particle-in-cell (PIC) simulations in one and two spatial dimensions. Here, the modified two-stream instability is identified to be the main mechanism responsible for stopping the plasma jets and preventing species interpenetration. The electron and ion Hall parameters of the merged plasma are greater than unity, and the plasma β is close to unity, which is the desired characteristic of planned experiments at PLX. Our 2D PIC simulations validate the results of the radiation magneto-hydrodynamics code FLASH, which will be the primary tool for modeling various stages of future PJMIF experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1990

Regularized Chapman-Enskog expansion for scalar conservation laws

Rosenau has recently proposed a regularized version of the Chapman-Enskog expansion of hydrodynamics. This regularized expansion resembles the usual Navier-Stokes viscosity terms at law wave-numbers, but unlike the latter, it has the advantage of being a bounded macroscopic approximation to the linearized collision operator. The behavior of Rosenau regularization of the Chapman-Enskog expansion (RCE) is studied in the context of scalar conservation laws. It is shown that thie RCE model retains the essential properties of the usual viscosity approximation, e.g., existence of traveling waves, monotonicity, upper-Lipschitz continuity..., and at the same time, it sharpens the standard viscous shock layers. It is proved that the regularized RCE approximation converges to the underlying inviscid entropy solution as its mean-free-path epsilon approaches 0, and the convergence rate is estimated.

Schochet, Steven↗
NASA NTRS1994

An Approximate Solution for the Circumsolar Flow Field of a Sun Moving Through the Local Interstellar Medium (LISM)

The solar system is in relative motion with respect to the ambient interstellar medium. The supersonic solar wind is expected to pass through the termination shock, thus the solar wind plasma eventually has to enter into an asymptotic outflow geometry appropriately adopted to this counterflow situation. Many attempts have been done to simulate the interaction between the solar wind and the LISM numerically. In this paper we generalize a Parker type analytical solution of the counterflow. The idea is to introduce a special kind of compressibility of the solar wind flow. With the assumption that only a transversal component of the density gradient normal to the flow lines exists we are able to calculate a full set of hydrodynamical quantities describing the circumsolar flow field of a Sun moving through the LISM. The equations governing the velocity and density fields lead to analytical solutions which can be taken as good approximations to the more general case of compressible plasma flows.

Ratkiewicz, Romana E.↗
DOE OSTI2026

Modeling thermal radiation waves in silica plasmas for the Mooncat NIF experiment

The Mooncat experiment on the National Ignition Facility uses a laser-driven hohlraum to create a thermal radiation wave in a titanium-doped silica plasma. The titanium dopant enables absorption spectroscopy measurements to infer the temperature of the wave as it propagates. This measurement can be used to constrain multi-physics simulation codes to better understand when simulations do not match an experiment. In this paper, we present radiation-hydrodynamics simulations of the thermal radiation wave in the first full-platform shots of the Mooncat experiment. We examine the important parameters of the simulation, focusing on the radiation temperature source, the material model of the silica plasma as it pertains to radiation transport, and lateral leakage through a beryllium tube enclosing the silica. We compare different simulation modeling strategies to an analytic model of diffusive radiation transport and find that the simulation agrees with the analytic model when it is sufficiently simplified. These simulations show how radiation energy couples to matter to develop a shock wave in a radiative heat wave, an important topic in astrophysics and nuclear fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1982

The dynamics of the Venus ionosphere. II - The effects of the time scale of the solar wind dynamic pressure variations

The effects on the upper dayside Venus ionosphere of a slow increase in solar wind dynamic pressure are simulated numerically with a one-dimensional (spherically symmetric) Lagrangian hydrodynamical code. The simulation is started with an extended ionosphere in pressure equilibrium with the solar wind at the ionopause. The pressure at the ionopause is gradually increased to five times the initial pressure with rise times of 5, 15, and 30 min. It is found that, for rise times greater than about 10 min, the compression of the ionopause is nearly adiabatic, with the ionopause moving downward at velocities of approximately 1-2 km/sec until it reaches a maximally compressed state, at which time the motion reverses. For short rise times the compression produces a shock wave similar to that occurring in the case of a sudden increase in pressure. The global implications of these processes are discussed within the context of Pioneer Venus observations and future theoretical work on this problem is outlined.

Stein, R. F.↗
NASA NTRS1986

The collision of high-velocity clouds with a galactic disk

Two-dimensional hydrodynamic simulations for the interaction of high-velocity clouds with a galactic disk are presented. The impinging clouds are assumed to be spherical and the target disk is represented by a constant density slab, n(g) = 1/cu cm, with a total width W(g) = 200 pc. The numerical experiments cover a wide range of cloud densities, between 0.1 and 100/cu cm, and velocities between 100 and 300 km/s. At a time approximately 10 to the 7th yr after impact, two types of final configurations are found. In the first case, the infalling cloud is completely shocked in a time short compared with the crossing time of the disk. Then, the generated cavity has time to grow sideways and large scale structures with a round shape, and in some cases nearly spherical, are produced. In the second case, which occurs for high density clouds, the cloud is shocked on a time scale longer than or comparable to the crossing time. The resultant cylindrical holes drilled across the entire disk have the dimensions of the impinging cloud. Cloud-galaxy interactions are compared with other energy sources and the morphologies of the resultant structures are suggested to resemble the large scale structures observed in H I.

Tenorio-Tagle, G.↗
NASA NTRS1994

Two-dimensional hydrodynamical simulations of wind-compressed disks around rapidly rotating B stars

We use a two-dimensional piecewise parabolic method (PPM) code to simulate numerically the hydrodynamics of a radiation-driven stellar wind from a rapidly rotating Be star. The results generally confirm predictions of the semianalytic 'wind-compressed disk' model recently proposed by Bjorkman and Cassinelli to explain the circumstellar disks inferred observationally to exist around such rapidly rotating stars. However, this numerical simulation is able to incorporate several important effects not accounted for in the simple model, including a dynamical treatment of the outward radiative driving and gas pressure, as well as a rotationally distorted, oblate stellar surface. This enables us to model quantitatively the compressed wind and shock that forms the equatorial disk. The simulation results thus do differ in several important details from the simple method, showing, for example, an inner disk inflow not possible in the heuristic approach of assuming a fixed outward velocity law. There is also no evidence for the predicted detachment of the disk that arises in the fixed outflow picture. The peak equatorward velocity in the dynamical models is furthermore about a factor of 2 smaller than the lytically predicted value of approximately 50% of the stellar equatorial rotation speed. As a result, the dynamical disks are somewhat weaker than predicted, with a wider opening angle, lower disk/pole density ratio, and smaller shock velocity jump. The principal cause of these latter differences appears to be an artificially strong equatorward drift of the subsonic outflow in the original analytic model. Much better agreement with the dynamical results can be obtained, however, from a slightly modified, analytic wind-compression model with a more detailed specification of the fixed wind outflow and a lower boundary set to the sonic radius along a rotationally oblate stellar surface. Hence, despite these detailed differences, the general predicted effect of disk formation by wind compression toward the equator is substantially confirmed.

Owocki, Stanley P.↗
DOE OSTI2021

Determining the Structure of Rotating Massive Stellar Cores with Gravitational Waves

The gravitational wave (GW) signal resulting from stellar core collapse encodes a wealth of information about the physical parameters of the progenitor star and the resulting core-collapse supernova (CCSN). We present a novel approach to constrain CCSN progenitor properties at collapse using two of the most detectable parts of the GW signal: the core-bounce signal and evolution of the dominant frequency mode from the protoneutron star. We focus on the period after core bounce but before explosion and investigate the predictive power of GWs from rotating CCSNe to constrain properties of the progenitor star. We analyze 34 2D and four 3D neutrino-radiation-hydrodynamic simulations of stellar core collapse in progenitors of varied initial mass and rotation rate. Extending previous work, we verify the compactness of the progenitor at collapse to correlate with the early ramp-up slope, and in rotating cases, also with the core angular momentum. Combining this information with the bounce signal, we present a new analysis method to constrain the pre-collapse core compactness of the progenitor. Because these GW features occur less than a second after core bounce, this analysis could allow astronomers to predict electromagnetic properties of a resulting CCSN even before shock breakout.

Astronomy & Astrophysics↗
NASA NTRS1990

Thermal instabilities in protogalactic clouds

The means by which a protogalaxy can fragment to form the first generation of stars and globular clusters remains an important problem in astrophysics. Gravitational instabilities grow on timescales too long to drive fragmentation before the background density grows by many orders of magnitude (see Murray and Lin 1989a, and references therein). Thermal instability provides a much more likely mechanism. After its initial collapse, a protogalactic cloud is expected to be shock heated to its virial temperature approx. 10(exp 6) K. Cooling by H and He+ below 10(exp 6) K has a negative slope, so that the cloud is subject to strong thermal instabilities. Density enhancements may then grow rapidly, fragmenting the protogalaxy as it cools to lower temperatures. The role of dynamical effects upon the growth of perturbations is considered here. The method used is similar to that used in Murray and Lin (1989a; see also the Erratum to appear September 15), which examined the growth of thermal instabilities with a one-dimensional Lagrangian hydrodynamics code, written for spherical symmetry. Perturbed regions therefore take the form of shells. The dynamical variables are integrated explicitly, while the temperature, ionization fraction, and molecular fraction are integrated implicitly, and account is taken for non-equilibrium values of these quantities.

Murray, S. D.↗
NASA NTRS1993

Precision linear shaped charge severance of graphite-epoxy materials

This paper presents Precision Linear Shaped Charge (PLSC) components designed to sever a variety of target materials. Recent data for the severance of graphite-epoxy panels or targets with PLSC's are presented. A brief history of the requirement to originate the development of PLSC's for weapon components at Sandia National Laboratories is presented. The Department of Energy's (DOE) nuclear weapon systems have continually decreased in size. Today's relatively small weapons require the design of much more efficient, lighter, and smaller explosive components because fragments, air shocks, and pyro-shocks associated with the function of these components can damage electrical and other sensitive components located nearby. The DOE requirements for PLSC's are listed. Therefore, linear shaped charge (LSC) components for weapon systems can no longer be empirically or experimentally designed for a given application. Many of today's designs require severing concentric cylinders, for example, where the LSC jet is designed to sever only one of the two cylinders as was the case for the B90/Nuclear Depth Strike Bomb. Therefore, code modeling and simulation technology must be utilized to obtain a better understanding of the LSC jet hydrodynamic penetration, fracture, shear, and spall mechanisms associated with the severance of metallic as well as composite targets.

Vigil, Manuel G.↗
NASA NTRS1995

The formation of protostellar disks. 2: Disks around intermediate-mass stars

Hydrodynamical calculations of the evolution of a collapsing, rotating axisymmetric 10 solar masses molecular clump, including the effects of radiative acceleration but without magnetic fields, are represented. The initial cloud is assumed to be uniformly rotating, centrally condensed sphere with rho is proportional to r(exp -2). Several cases are considered, in which both the overall clump size and the total amount of angular momentum are varied. The calculations show how a warm, quasi-hydrostatic disk surrounding a central unresolved core of only a few solar masses forms and grows in size and mass. The disk is encased in two distinct accretion shock fronts, both of which are several scale heights above the equatorial plane. At the end of the calculation of our standard case, the central unresolved region is found to have a mass of 2.7 solar masses and a ratio of rotational to gravitational energy of approximately 0.45, sufficiently large to be unstable to nonaxisymmetric perturbations. In addition, the inner portions of the disk containing most of the mass are unstable according to the local Toomre criterion, implying that also in this region nonaxisymmetric perturbations will lead to rapid evolution. Under the assumption that gravitational torques would transport angular momentum out of this region, a central core of less than or approximately 8 solar masses with a stable disk of greater than or approximately = 2 solar masses should result. Frequency-dependent radiative transfer calculations of the standard case at selected ages show how the continuum spectrum of the structure depends on the disk's orientation and age and how the observed isophotal contours vary with wavelength. Because of the strong dependence on viewing angle, continuum spectra alone should not be used to estimate the evolutionary stage of development of these objects. Comparable results were obtained for the other cases considered.

Yorke, Harold W.↗
NASA NTRS1990

Hydrodynamic simulations of stellar wind disruption by a compact X-ray source

This paper presents two-dimensional numerical simulations of the gas flow in the orbital plane of a massive X-ray binary system, in which the mass accretion is fueled by a radiation-driven wind from an early-type companion star. These simulations are used to examine the role of the compact object (either a neutron star or a black hole) in disturbing the radiatively accelerating wind of the OB companion, with an emphasis on understanding the origin of the observed soft X-ray photoelectric absorption seen at late orbital phases in these systems. On the basis of these simulations, it is suggested that the phase-dependent photoelectric absorption seen in several of these systems can be explained by dense filaments of compressend gas formed in the nonsteady accreation bow shock and wake of the compact object.

Blondin, John M.↗
DOE OSTI2022

Experimental evidence of early-time saturation of the ion-Weibel instability in counterstreaming plasmas of CH, Al, and Cu

The collisionless ion-Weibel instability is a leading candidate mechanism for the formation of collisionless shocks in many astrophysical systems, where the typical distance between particle collisions is much larger than the system size. Multiple laboratory experiments aimed at studying this process utilize laser-driven (I≳10 15 W/cm 2 ), counterstreaming plasma flows (V≲ 2000 km/s) to create conditions unstable to Weibel-filamentation and growth. This technique intrinsically produces temporally varying plasma conditions at the midplane of the interaction where Weibel-driven B fields are generated and studied. Experiments discussed herein demonstrate robust formation of Weibel-driven B fields under multiple plasma conditions using CH, Al, and Cu plasmas. Linear theory based on benchmarked radiation-hydrodynamic FLASH calculations is compared with Fourier analyses of proton images taken ~5–6 linear growth times into the evolution. Our new analyses presented here indicate that the low-density, high-velocity plasma-conditions present during the first linear-growth time (~300–500 ps) sets the spectral characteristics of Weibel filaments during the entire evolution. It is shown that the dominant wavelength (~300 μm) at saturation persists well into the nonlinear phase, consistent with theory under these experimental conditions. However, estimates of B-field strength, while difficult to determine accurately due to the path-integrated nature of proton imaging, are shown to be in the ~10–30 T range, an order of magnitude above the expected saturation limit in homogenous plamas but consistent with enhanced B fields in the midplane due to temporally varying plasma conditions in experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1972

A numerical hydrodynamic study of coalescence in head-on collisions of identical stars.

The specific results reported refer to head-on collisions between identical polytropes of index 3 having solar mass and radius. If the polytropes were initially at rest at infinity, then about 5% of the combined mass is lost by ejection following collision. The volatilized mass fraction rises to about 18% for an initial relative collision velocity of 1000 km/sec at infinite separation, and to about 60% for the 2000 km/sec case. Since the initial kinetic and gravitational energies balance for a relative velocity of 1512 km/sec at infinity, it may be seen that net coalescence persists to velocities somewhat in excess of this figure. Mass ejection takes place in two ways simultaneously: (1) by a rapid sideward expulsion of fluid in a massive lateral sheet normal to the collision axis, and (2) as a result of two recoil shocks which lead momentum flows backward along this axis. The lateral effect has similarities to the expansion of gas into a vacuum i.e., shocks are not involved. However, the ejection of material from the rear colliding hemisphere due to the recoil shocks predominates at low collision velocities. As the velocity increases, both effects strengthen, but the lateral expulsion intensifies more rapidly than the recoil shocks.

Seidl, F. G. P.↗
DOE OSTI2020

Experimental characterization of a section of a spherically imploding plasma liner formed by merging hypersonic plasma jets

In this work, we report experimental results on merging of hypersonic plasma jets, which is the fundamental building block for forming spherically imploding plasma liners as a potential standoff compression driver for magneto-inertial fusion. Jets are formed and launched by contoured-gap coaxial plasma guns mounted at the six vertices and the center of a hexagon covering approximately one-tenth of the surface area of a 9-ft-diameter spherical chamber. First, from experiments with two and three merging jets of four different species (N, Ar, Kr, and Xe), we show that (1) density spatial non-uniformities can be large (with electron-density jumps ranging from 2.9 for N to 6.6 for Xe) when shocks form upon jet merging, but smaller (density jumps <2) when shocks do not form; (2) jet impurities (20% Ti in these experiments) can increase the level of density spatial non-uniformity by increasing the collisionality of jet merging, leading to shock formation rather than potentially more desirable shockless jet merging; and (3) the liner Mach number can remain high (≳10), as required for plasma liners to be an effective compression driver. Secondly, from experiments with six and seven merging jets using Ar, we present results with improved jet-to-jet mass balance of <2% across jets, including (1) evidence of substantially increased balance in the jet merging and symmetry of the liner structure and (2) potentially favorable changes in the jet-merging morphology with the addition of the seventh jet. Ultimately, for both experiments, we present comparisons between experimental and synthetic data from three-dimensional hydrodynamic codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Collision of localized shocks in AdS 5 as a series expansion in transverse gradients

We introduce a computational framework to more efficiently calculate the collision of localized shocks in five dimensional asymptotically Anti-de Sitter space. We expand the Einstein equations in transverse gradients and find that our numerical results agree well with exact solutions already at first order in the expansion. Moreover, the Einstein equations at first order in transverse gradients can be decoupled into two sets of differential equations. Here, the bulk fields of one of these sets has only a negligible contribution to boundary observables, such that the computation on each time slice can be simplified to the solution of several planar shockwave equations plus four further differential equations for each transverse plane ‘pixel’. At the cost of errors of ≲ 10% at the hydrodynamization time and for low to mid rapidities, useful numerical solutions can be sped up by roughly one order of magnitude.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗