SEARCH · Engineering Papers
Results for “Shock Physics”
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.
Magnetic Field Amplification by a Plasma Cavitation Instability in Relativistic Shock Precursors
Abstract Plasma streaming instabilities play an important role in magnetic field amplification and particle acceleration in relativistic shocks and their environments. However, in the far shock precursor region where accelerated particles constitute a highly relativistic and dilute beam, streaming instabilities typically become inefficient and operate at very small scales when compared to the gyroradii of the beam particles. We report on a plasma cavitation instability that is driven by dilute relativistic beams and can increase both the magnetic field strength and coherence scale by orders of magnitude to reach near-equipartition values with the beam energy density. This instability grows after the development of the Weibel instability and is associated with the asymmetric response of background leptons and ions to the beam current. The resulting net inductive electric field drives a strong energy asymmetry between positively and negatively charged beam species. Large-scale particle-in-cell simulations are used to verify analytical predictions for the growth and saturation level of the instability and indicate that it is robust over a wide range of conditions, including those associated with pair-loaded plasmas. These results can have important implications for the magnetization and structure of shocks in gamma-ray bursts, and more generally for magnetic field amplification and asymmetric scattering of relativistic charged particles in plasma astrophysical environments.
Observation of shock-front separation in multi-ion-species collisional plasma shocks
We observe shock-front separation and species-dependent shock widths in multi-ion-species collisional plasma shocks, which are produced by obliquely merging plasma jets of a He/Ar mixture (97% He and 3% Ar by initial number density) on the plasma liner experiment. Visible plasma emission near the He-I 587.6-nm and Ar-II 476.5–514.5-nm lines is simultaneously recorded by splitting a single visible image of the shock into two different fast-framing cameras with different narrow bandpass filters (589 ± 5 nm for observing the He-I line and 500 ± 25 nm for the Ar-II lines). For conditions in these experiments (pre-shock ion and electron densities ≈5×10 14 cm –3 , ion and electron temperatures of ≈2.2 eV, and relative plasma-merging speed of 22 km/s), the observationally inferred magnitude of He/Ar shock-front separation and the shock widths themselves are <1 cm, which correspond to ~50 post-shock thermal ion–ion mean free paths. These experimental length scales are in reasonable qualitative and quantitative agreement with results from 1D multi-fluid simulations using the CHICAGO code. However, we note there are differences between the experimentally inferred and simulation-predicted ionization states and line emission intensities, particularly in the post-shock region. Overall, the experimental and simulation results are consistent with theoretical predictions that the lighter He ions diffuse farther ahead within the overall shock front than the heavier Ar ions.
Magnetohydrodynamic and gasdynamic aspects of solar-wind flow around terrestrial planets: A critical review
Problems in space physics are discussed whose models, in simplified form, reduce to a supersonic flow scheme with a detached shock wave, namely: (1) solar wind interaction with an intrinsic planetary magnetic field; (2) solar wind interaction with the ionized component of the atmosphere of a comet; and (3) solar wind interaction with the ionosphere of a planet which does not possess its own magnetic field. The numerical study of the above problems is performed with the use of magnetogasdynamic equations for an ideal single fluid model. From the physical veiwpoint, the problems are solved in terms of as simple phenomena as possible; the principal objective is to make recently developed methods of numerical analysis of mixed flows applicable to space physics problems.
Robust and Accurate Shock Capturing Method for High-Order Discontinuous Galerkin Methods
A simple yet robust and accurate approach for capturing shock waves using a high-order discontinuous Galerkin (DG) method is presented. The method uses the physical viscous terms of the Navier-Stokes equations as suggested by others; however, the proposed formulation of the numerical viscosity is continuous and compact by construction, and does not require the solution of an auxiliary diffusion equation. This work also presents two analyses that guided the formulation of the numerical viscosity and certain aspects of the DG implementation. A local eigenvalue analysis of the DG discretization applied to a shock containing element is used to evaluate the robustness of several Riemann flux functions, and to evaluate algorithm choices that exist within the underlying DG discretization. A second analysis examines exact solutions to the DG discretization in a shock containing element, and identifies a "model" instability that will inevitably arise when solving the Euler equations using the DG method. This analysis identifies the minimum viscosity required for stability. The shock capturing method is demonstrated for high-speed flow over an inviscid cylinder and for an unsteady disturbance in a hypersonic boundary layer. Numerical tests are presented that evaluate several aspects of the shock detection terms. The sensitivity of the results to model parameters is examined with grid and order refinement studies.
Shock Initiation Microscopy with High Time and Space Resolution
In this paper, we describe studies of shock initiation and shock-to-detonation transitions in energetic materials using a tabletop shock compression microscope with nanosecond time resolution and micrometer spatial resolution. Planar input shocks with durations of 4–20 ns are produced using 0–4.5 km/s laser-launched flyer plates. Emphasis is on measurements of temperature, velocities, pressure, and microstructure using photon Doppler velocimetry (PDV), optical pyrometry and high-speed videography. Techniques are discussed for fabricating disposable shock target arrays of tiny plastic-bonded explosives (PBX), liquid and powder explosives, and single-crystal explosives for high-throughput studies. Optical temperature measurements of shocked triaminotrinitrobenzene (TATB) are discussed. Since TATB is yellow, we developed methods to correct for the blue absorption to obtain more accurate temperatures. Hot spots in shocked polymer-encased HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine) crystals are observed in real-time, showing a hot spot produced in a collapsing void that ignites a deflagration. Despite the small dimensions of our explosive charges (typically 1 mm diameter and 250 μm length), we produced reproducible detonation states in solid and liquid explosives using short-duration shocks near the von Neumann spike (VNS) pressure. We show the VNS pressure is associated with a transition to high-efficiency gas production from the explosives. In studies of NM, prior to detonation, we see reaction originating at hot spots which coalesce to form a superdetonation.
Computational Relativistic Astrophysics Using the Flowfield-Dependent Variation Theory
Theoretical models, observations and measurements have preoccupied astrophysicists for many centuries. Only in recent years, has the theory of relativity as applied to astrophysical flows met the challenges of how the governing equations can be solved numerically with accuracy and efficiency. Even without the effects of relativity, the physics of magnetohydrodynamic flow instability, turbulence, radiation, and enhanced transport in accretion disks has not been completely resolved. Relativistic effects become pronounced in such cases as jet formation from black hole magnetized accretion disks and also in the study of Gamma-Ray bursts (GRB). Thus, our concern in this paper is to reexamine existing numerical simulation tools as to the accuracy and efficiency of computations and introduce a new approach known as the flowfield-dependent variation (FDV) method. The main feature of the FDV method consists of accommodating discontinuities of shock waves and high gradients of flow variables such as occur in turbulence and unstable motions. In this paper, the physics involved in the solution of relativistic hydrodynamics and solution strategies of the FDV theory are elaborated. The general relativistic astrophysical flow and shock solver (GRAFSS) is introduced, and some simple example problems for Computational Relativistic Astrophysics (CRA) are demonstrated.
Analysis of Experimental Shock Propagation through a Diverse Set of Geometric Cavities Embedded in Lab-Scale Polymer Cubes
Experiments were designed and conducted to investigate the impact that geometric cavities have on the transfer of energy from an embedded explosion to the surface of the physical domain. The experimental domains were fabricated as 3-inch polymer cubes, with varying cavity geometries centered in the cubes. The energy transfer, represented as a shock wave, was generated by the detonation of an exploding bridgewire at the center of the cavity. The shock propagation was tracked by schlieren imaging through the optically accessible polymer. The magnitude of energy transferred to the surface was recorded by an array of pressure sensors. A minimum of five experimental runs were conducted for each cavity geometry and statistical results were developed and compared. Results demonstrated the decoupling effect that geometric cavities produce on the energy field at the surface.
Entropy Splitting and Numerical Dissipation
A rigorous stability estimate for arbitrary order of accuracy of spatial central difference schemes for initial-boundary value problems of nonlinear symmetrizable systems of hyperbolic conservation laws was established recently by Olsson and Oliger (1994) and Olsson (1995) and was applied to the two-dimensional compressible Euler equations for a perfect gas by Gerritsen and Olsson (1996) and Gerritsen (1996). The basic building block in developing the stability estimate is a generalized energy approach based on a special splitting of the flux derivative via a convex entropy function and certain homogeneous properties. Due to some of the unique properties of the compressible Euler equations for a perfect gas, the splitting resulted in the sum of a conservative portion and a non-conservative portion of the flux derivative. hereafter referred to as the "Entropy Splitting." There are several potential desirable attributes and side benefits of the entropy splitting for the compressible Euler equations that were not fully explored in Gerritsen and Olsson. The paper has several objectives. The first is to investigate the choice of the arbitrary parameter that determines the amount of splitting and its dependence on the type of physics of current interest to computational fluid dynamics. The second is to investigate in what manner the splitting affects the nonlinear stability of the central schemes for long time integrations of unsteady flows such as in nonlinear aeroacoustics and turbulence dynamics. If numerical dissipation indeed is needed to stabilize the central scheme, can the splitting help minimize the numerical dissipation compared to its un-split cousin? Extensive numerical study on the vortex preservation capability of the splitting in conjunction with central schemes for long time integrations will be presented. The third is to study the effect of the non-conservative proportion of splitting in obtaining the correct shock location for high speed complex shock-turbulence interactions. The fourth is to determine if this method can be extended to other physical equations of state and other evolutionary equation sets. If numerical dissipation is needed, the Yee, Sandham, and Djomehri (1999) numerical dissipation is employed. The Yee et al. schemes fit in the Olsson and Oliger framework.
High-energy-density-physics measurements in implosions using Bayesian inference
Convergent high-energy-density (HED) experimental platforms are used to study matter under some of the most extreme conditions that can be produced on Earth, comparable to the interior of stars. There are many challenges in using these systems for fundamental measurements currently being addressed by new analysis methods, such as the combination of a reduced physics model and Bayesian inference, allowing a self-consistent inference of physical quantities with a robust error analysis. These methods in combination with simple (as compared to inertial confinement fusion implosions) implosion platforms, which can be modified to show sensitivity to different physical mechanisms of interest, are used to study the physical properties of matter under extreme conditions. This work discusses a subset of implosion targets for studying opacity effects, electron-ion equilibration, and thermal conductivity and, as an example, a system consisting of a thick-shelled, gas-filled laser-direct-drive implosion is used to show how a reduced model and Bayesian inference can help inform experimental design decisions such as diagnostic choice. It is shown that for this system that a combination of neutron and x-ray self-emission diagnostics is critical for constraining the details of the thermodynamic states in the system and that the conductivity exponent in a Spitzer like framework can be constrained to the 30% level in deuterium at gigabar conditions. Furthermore, this process can be applied to many HED systems to make underlying model assumptions explicit and facilitate experimental design and analysis.
On the importance of three-dimensional modeling for high-energy-density physics experiments
Laser-driven cylindrical implosion experiments enable direct measurements of hydrodynamic instability growth in convergent geometries, providing a wealth of validation data in the high-energy-density regime. These experiments are designed to be nearly axially invariant, allowing for modeling with complementary two-dimensional slices of the cylinder. Two distinct hydrodynamics codes are employed to model a subset of these experiments, and the results are shown to be in very good agreement with each other and the available experimental data. While this 2D modeling approach adequately captures most of the physics of the implosion and ensuing instability growth, there are crucial aspects from the three-dimensional nature of the experiments that are missed in 2D. The first fully 3D simulations of these experiments are presented, and small but significant differences are found to arise from both the axial and azimuthal non-uniformity in the laser drive. Recent experimental results confirming the drive asymmetry are discussed.
Study of ablation and shock generation across three orders of magnitude of laser intensity with 100 ps laser pulses
The laser ablation and subsequent shock generation in solid targets plays an important role in a variety of research topics from equation of state models for materials to inertial confinement fusion. One of the long-standing issues is the knowledge of ablation depth in the picosecond time regime. Here, we report on a direct technique for determining the ablation depth in aluminum using x-ray diffraction data from Linac Coherent Light Source at the Stanford Linear Accelerator Center. This technique gives a direct measurement of the shock wave propagation in the bulk target, enabling an ability to discern early timescale physics from late timescale effects not available in postmortem analysis. We find that the ablation depths only vary by 0.2 μm across three orders of magnitude of laser intensity, while the pressure increased by a factor of 10 following a square root dependence on laser pulse energy. We further observe that the ablation depth in this intensity range (10 11 –10 13 W/cm 2 in intensity, corresponding to 0.8–80 J/cm 2 in fluence) cannot be modeled by a universal scaling law, given the complexity of the mechanisms governing laser ablation in this intensity regime.
Measurement of Kinetic-Scale Current Filamentation Dynamics and Associated Magnetic Fields in Interpenetrating Plasmas
Here, we present the first local, quantitative measurements of ion current filamentation and magnetic field amplification in interpenetrating plasmas, characterizing the dynamics of the ion Weibel instability. The interaction of a pair of laser-generated, counterpropagating, collisionless, supersonic plasma flows is probed using optical Thomson scattering (TS). Analysis of the TS ion-feature revealed anticorrelated modulations in the density of the two ion streams at the spatial scale of the ion skin depth c/ω pi = 120 μm, and a correlated modulation in the plasma current. The inferred current profile implies a magnetic field amplitude ~30 ± 6 T , corresponding to ~1% of the flow kinetic energy, indicating that magnetic trapping is the dominant saturation mechanism.
Microscopy X-ray imaging enriched with small angle X-ray scattering for few nanometer resolution reveals shock waves and compression in intense short pulse laser irradiation of solids
Understanding how laser pulses compress solids into high-energy-density states requires diagnostics that simultaneously resolve macroscopic geometry and nanometer-scale structure. Here we present a combined X-ray imaging (XRM) and small-angle X-ray scattering (SAXS) approach that bridges this diagnostic gap. Using the Matter in Extreme Conditions end station at LCLS, we irradiated 25 μm copper wires with 45 fs, 0.9 J, 800 nm pulses at 3.5 × 10 19 W/cm 2 while probing with 8.2 keV XFEL pulses. XRM visualizes the evolution of ablation, compression, and inward-propagating fronts with ∼ 200 nm resolution, while SAXS quantifies their nanometer-scale sharpness via the time-resolved evolution of scattering streaks. The joint analysis reveals that an initially smooth compression steepens into a nanometer-sharp shock front after t sh ≈ (18 ± 3) ps, consistent with an analytical steepening model and hydrodynamic simulations. The front reaches a velocity of c sh ≈ 25 k m / s and a lateral width of several tens of microns, demonstrating direct observation of shock formation and decay at solid density for the first time with few-nanometer precision. This integrated XRM–SAXS method establishes a quantitative, multi-scale diagnostic of laser-driven shocks in dense plasmas relevant to inertial confinement fusion, warm dense matter, and planetary physics.
Shock tube instrumentation techniques for study of hypervelocity entry problems.
Shock tube instrumentation to study convective and radiative heat transfer of hot gases simulating hypervelocity reentry
Standard power regulator for the multi-mission modular spacecraft
The Standard Power Regulator Unit (SPRU) which forms the central building block of the Modular Power Subsystem (MPS) for the Multi-mission Modular Spacecraft (MMS) is described. A functional description of the SPRU is presented, detailing key design features, operational characteristics, and internal redundancy as well as giving block and schematic diagrams for all major subassemblies. The description of a qualification test program consisting of ten sequences ranging from a physical examination of the unit to calculating vibration and sine sweep of power on shock tests is presented. The results showed that some significant anomalies were encountered but no catastrophic failures were indicated, although several power module failures were found due to overvoltage malfunctions of the solar array stimulator.
Physics of magnetospheric boundary layers
The central ideas of this grant are that the magnetospheric boundary layers link disparate regions of the magnetosphere together, and the global behavior of the magnetosphere can be understood only by understanding the linking mechanisms. Accordingly the present grant includes simultaneous research on the global, meso-, and micro-scale physics of the magnetosphere and its boundary layers. These boundary layers include the bow shock, magnetosheath, the plasma sheet boundary layer, and the ionosphere. Analytic, numerical and simulation projects have been performed on these subjects, as well as comparison of theoretical results with observational data. Very good progress has been made, with four papers published or in press and two additional papers submitted for publication during the six month period 1 June - 30 November 1993. At least two projects are currently being written up. In addition, members of the group have given papers at scientific meetings. The further structure of this report is as follows: section two contains brief accounts of research completed during the last six months, while section three describes the research projects intended for the grant's final period.
Variability in shell models of GRBs
Many cosmological models of gamma-ray bursts (GRBs) assume that a single relativistic shell carries kinetic energy away from the source and later converts it into gamma rays, perhaps by interactions with the interstellar medium or by internal shocks within the shell. Although such models are able to reproduce general trends in GRB time histories, it is difficult to reproduce the high degree of variability often seen in GRBs. The authors investigate methods of achieving this variability using a simplified external shock model. Since the model emphasizes geometric and statistical considerations, rather than the detailed physics of the shell, it is applicable to any theory that relies on relativistic shells. They find that the variability in GRBs gives strong clues to the efficiency with which the shell converts its kinetic energy into gamma rays.