Engineering Papers⌕ Search

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.

At least 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Refractive index of lithium fluoride at high dynamic stresses

Alkali halides are prototypical ionic solids whose refractive index is a fundamental property related to their lattice structure and ionic polarizability. In particular, lithium fluoride (LiF) has the largest band gap of any known transparent material and maintains transparency at high pressures, making it well suited for use as an optical window in dynamic compression experiments. While empirical fits to the measured density dependence of the refractive index have been provided, a model that is valid over the entire pressure range of experiments and based on a polarization-based theoretical description is lacking. We present a refractive index model for dynamically compressed LiF based on the Lorentz-Lorenz equation, where the molecular polarizability is determined using a single oscillator model with a strain polarizability parameter Λ=0.73. Here, we show that our modeling approach provides an excellent match to the LiF refractive index data for both shock and ramp compression experiments to 900 GPa (density compression greater than threefold). Additionally, updated fits are provided to determine the refractive index correction for LiF windows used in dynamic compression experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Shock state distributions in porous tantalum and characterization with multipoint velocimetry

Heterogenous materials under shock compression can be expected to reach different shock states throughout the material according to local differences in microstructure and the history of wave propagation. Here, in this work, a compact, multiple-beam focusing optic assembly is used with high-speed velocimetry to interrogate the shock response of porous tantalum films prepared through thermal-spray deposition. The distribution of particle velocities across a shocked interface is compared to results obtained using a set of defocused interferometric beams that sampled the shock response over larger areas. The two methods produced velocity distributions along the shock plateau with the same mean, while a larger variance was measured with narrower beams. The finding was replicated using three-dimensional, mesoscopically resolved hydrodynamics simulations of solid tantalum with a pore structure mimicking statistical attributes of the material and accounting for radial divergence of the beams, with agreement across several impact velocities. Accounting for pore morphology in the simulations was found to be necessary for replicating the rise time of the shock plateau. The validated simulations were then used to show that while the average velocity along the shock plateau could be determined accurately with only a few interferometric beams, accurately determining the width of the velocity distribution, which here was approximately Gaussian, required a beam dimension much smaller than the spatial correlation lengthscale of the velocity field, here by a factor of ∼30×, with implications for the study of other porous materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Numerical simulations of the Los Alamos gapstick experiment

The gapstick is a high explosive (HE) sensitivity test recently developed at Los Alamos. The experiment design was motivated by the traditional gap test and consists of a series of HE and inert pellets in a rate stick configuration. The inert pellets are made increasingly longer so that eventually the attenuated shock is unable to initiate the next HE pellet. While the gapstick is a mechanically simple experiment it poses several challenges for numerical simulation. In particular, accurately modeling the HE initiation and detonation phenomena requires a reactive burn model with sufficient mesh resolution to capture the reaction scales. In this work, the Scaled Uniform Reactive Front (SURF) and Arrhenius Wescott-Stewart-Davis (AWSD) burn models are used for simulations of a PBX 9501 (95 wt% HMX, 5% binder) gapstick. Furthermore, the material model for the inert pellet material, Polyvinylidene Fluoride (PVDF), has a direct influence on the shock propagation and a new equation of state (EOS) calibration for PVDF is developed using available Hugoniot data. Despite some challenges, the simulations are able to reasonably predict shock transit velocities and detonation failure in the gapstick.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Anharmonicity and Phase Diagram of Magnesium Oxide in the Megabar Regime

With density functional molecular dynamics simulations, we computed the phase diagram of MgO from 50 to 2000 GPa up to 20000 K. Via thermodynamic integration (TDI), we derive the Gibbs free energies of the B1, B2, and liquid phases and determine their phase boundaries. With TDI and a pseudo-quasi-harmonic approach, we show that anharmonic effects are important and stabilize the B1 phase in particular. As a result, the B1-B2 transition boundary in the pressure-temperature plane exhibits a steep slope. We predict the B1-B2-liquid triple point to occur at approximately T= 10000 K and P = 370 GPa, which is higher in pressure than was inferred with quasi-harmonic methods alone. We predict the principal shock Hugoniot curve to enter the B2 phase stability domain but only over a very small range of parameters. Here, this may render it difficult to observe this phase with shock experiments because of kinetic effects.

3-dimensional systems↗
DOE OSTI2023

New Sesame Equations of State for Polypropylene, Polyvinylchloride (PVC), and Nylon 66

I built new Sesame equations of state (EOS) for polypropylene (PP), polyvinylchloride (PVC), and nylon 66. PP has the same stoichiometry as polyethylene and TPX but with a density intermediate between the two. Nylon is a polyamide manufactured chiefly in the form of fibers. PP and PVC are the second- and third-most widely produced polymers in the world. All three are thermoplastics. Their monomer units are depicted in Figure 1. I discuss the theory on which the EOS are based in the following section, then compare results to thermal and shock compression data in Section III.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Using neutrons and x rays to measure plasma conditions in a solid sphere of deuterated polyethylene compressed to densities of 35 g/cc at temperatures of 2 keV and pressures of 40 Gbar

This paper describes an experiment that shock compresses the center of a solid deuterated polyethylene sphere, CD 2 , to densities of 35 g/cc and temperatures of 2 keV with corresponding pressure of 40 Gbar. The design employs a strong spherically converging shock launched through a solid ball of material using a Hohlraum radiation drive. As the shock coalesces at the center it produces a hot spot that we characterize by measuring the x-ray self-emission and 2.45 MeV neutrons emitted. Two-dimensional images and time-resolved measurements of the x rays emitted determine the size and time duration of the hot spot, leading to an estimated 2k eV electron temperature. The neutron time of flight spectrometer measures an average ion temperature of 1.06 ± 0.15 keV and neutron yield of 7.0 (±0.5) × 10 9 DD neutrons. Our new distribution function tool enables us to create a forward model of the experimental data based on 1D radiation-hydrodynamic simulations, leading to a better understanding of the plasma conditions that produce the measured neutrons and x rays. Furthermore, our simulations indicate that the x rays are produced in a short-lived hot-dense core over tens of picoseconds, whereas the neutron emission continues for about 200 ps, as the hot core starts to expand, thereby leading to a lower mean temperature of the plasma during neutron production. This finding is in agreement with the experimental data, and we therefore conclude that the forward-modeling is a useful tool for inferring the conditions of the hot spot in a laser-driven implosion during burn.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

RELAP-7 Application and Enhancement for FLEX Strategies and ATF Behavior under Extended Loss of AC Power Conditions

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.

42 ENGINEERING↗
DOE OSTI2024

An exploration of anomalous electrical noise in shocked cyclotrimethylenetrinitramine (RDX)-based explosives

Gas gun shock experiments on cyclotrimethylenetrinitramine (RDX)-based explosive compositions that employ embedded gauge particle velocity tracers have noted a significant amount of electrical noise when compared to other explosive formulations. This paper reexamines previously published embedded gauge data on Cyclotols (60–80 wt. % RDX) to quantify the electromagnetic behavior of these materials. The primary observation is a fourfold increase in the electrical noise when Cyclotols are shocked above 4.22 ± 0.08 GPa. Electromagnetic gauge noise is also observed within particle velocity traces in reactive growth and off-Hugoniot shocks, although at higher pressures than the direct shock case, suggesting a temperature- or kinetically dependent transition. In all cases, the electrical noise disappears upon detonation. By comparing with the static high-pressure phase diagram of RDX, we interpret this change in electromagnetic behavior to be a change in the RDX crystal structure to a piezoelectric phase, although it is uncertain whether the γ or ε phase is responsible for the observed behavior.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

A detailed chemical kinetic modeling and experimental investigation of the low- and high-temperature chemistry of n-butylcyclohexane

Chemical kinetic models of gasoline, jet, and diesel fuels and their mixtures with bioderived fuels are needed to assess fuel property effects on efficiency, emissions, and other performance metrics in internal combustion and gas turbine engines. As these real fuels have too many fuel components to be included in a chemical kinetic model, surrogate fuels containing fewer components are used to represent them. These surrogate fuels mimic the chemical classes or molecular structures contained in the real fuel. One of the important chemical classes in gasoline, jet, and diesel fuels comprises cyclohexanes. Cyclohexanes comprise about 30% or more by weight in diesel fuel. Also, Mueller et al proposed n-butylcyclohexane (nBCH) as a component in a nine-component surrogate palette to represent the ignition properties, distillation curve, density, and molecular structures of a diesel certification fuel. Herein, experimental measurements of the ignition delay times (IDTs) of nBCH in a shock tube and in a rapid compression machine are reported over a wide range of temperature, pressure, and equivalence ratio important for enabling the validation of a chemical kinetic model for nBCH for combustion in diesel engines. The range of conditions are temperatures of 630–1420 K, pressures of 10, 30, and 50 bar, and equivalence ratios of 0.3, 0.5, 1.0, and 2.0 in ‘air’. A detailed chemical model is developed for nBCH to simulate its ignition at both low and high-temperature conditions and at relevant elevated pressures. The experimentally measured IDTs are used to improve and validate the chemical kinetic model.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Synthesis and investigation into explosive sensitivity for a series of new picramide explosives

Tailoring the molecular properties that govern energetic material sensitivity is essential to improve safety and help develop new energetic materials. Despite this need, understanding the complex chemistry and physics of explosive initiation and propagation is still a challenge. Recent work by our group has reinforced the view that explosive sensitivity under sub-shock conditions is connected to the strength of the weakest covalent bond in the molecule, that is, its “trigger linkage.” These correlations have been observed with different classes of energetic molecules and indicate that “trigger linkage” bond breaking, and heat of explosion are good indicators for the sensitivity trends. Herein we report the synthesis of aliphatic energetic materials with ethane, propane and neopentane backbones. Experimental and computational studies show that the trigger linkage model, based on results from quantum molecular dynamics simulations, correctly predicts trends observed in the impact sensitivity of the molecules. However, while the model predicts the impact sensitivities of the ethane series, the neopentane series has higher impact sensitivities than predicted, which is presumably influenced by crystal packing effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

A theoretical approach for transient shock strengthening in high-energy-density laser compression experiments

In high-energy-density shock compression experiments, the desired state of compression is typically achieved by passing shock waves through a sequence of different materials. Here, a theoretical approach for transiently strengthening such shocks passing through interfaces of experimentally relevant impedance ratios is examined. A semi-analytical method based on characteristics analysis is used to solve the problem of a shock passing from one material to another through an intermediate region of non-uniform impedance; by appropriately designing this region, a greater shock strength can be achieved in the second material for a finite duration than in the absence of this region. When a shock passes into a material of higher impedance, the shock strengthening increases to a maximum before decreasing to an asymptotic value as the strength of the incident shock is increased. For shocks passing into materials of decreasing impedance, the shock strengthening increases monotonically as the strength of the incident shock is increased when the impedance ratio is above a critical threshold. Incorporating multiple intermediate materials can further increase the strength of the transmitted shock, with an exponential discretization of intermediate material impedances being the most effective distribution for strengthening strong shocks. The results suggest that up to 25% and 9% increases in pressure behind the leading shock can be achieved for materials of increasing and decreasing impedance, respectively. The technique is applied to the design of laser-driven dynamic compression experiments, and the results of the analysis are verified via comparison to simulations performed with the HYADES hydrodynamics code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Shock wave formation in radiative plasmas

The temporal evolution of weak shocks in radiative media is theoretically investigated in this work. The structure of radiative shocks has traditionally been studied in a stationary framework. Their systematic classification is complex because layers of optically thick and thin regions alternate to form a radiatively-driven precursor and a temperature-relaxation layer, between which the hydrodynamic shock is embedded. In this work, we analyze the formation of weak shocks when two radiative plasmas with different pressures are put in contact. Applying a reductive perturbative method yields a Burgers-type equation that governs the temporal evolution of the perturbed variables including the radiation field. The conditions upon which optically thick and thin solutions exist have been derived and expressed as a function of the shock strength and Boltzmann number. Below a certain Boltzmann number threshold, weak shocks always become optically thick asymptotically in time, while thin solutions appear as transitory structures. The existence of an optically thin regime is related to the presence of an overdense layer in the compressed material. Scaling laws for the characteristic formation time and shock width are provided for each regime. The theoretical analysis is supported by FLASH simulations, and a comprehensive testcase has been designed to benchmark radiative hydrodynamic codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Design of a High Toughness Epoxy for Superconducting Magnets and Its Key Properties

Nb 3 Sn accelerator magnets are poised to enable the luminosity upgrade of the Large Hadron Collider (LHC) at CERN, improving its potential for exploring physics beyond the standard model of particle physics. The prototype Nb 3 Sn magnets consistently need 10–25 quenches to achieve their best performance. A hypothesis is that the long training of these magnets may at least be partially induced by epoxy cracking and bonding failures. In 2018, we showed that several existing epoxies have a higher toughness and less tendency to crack at low temperatures than CTD-101K, the epoxy resin with which almost all Nb 3 Sn accelerator magnets have been impregnated. Here we explore a new high toughness formulation for reducing quench training of Nb 3 Sn accelerator magnets, through combining two amine curing agents to achieve a good glass transition temperature (T g ), adding a viscosity reducer to achieve low viscosity, and experimenting a coupling agent for improved bonding strengths. In this work, we report results of comprehensive materials tests, including thermal shock, tensile, compressive, shear tests, viscosity and T g tests.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Density variance dynamics in disparate shock tubes

This report discusses two experiments which investigate thin layer, heavy curtain fragmentation from the perspective of a Reynolds-Averaged mix model, in drastically disparate experimental regimes. The first, the centimeter/millisecond-scale “Horizontal Shock Tube” (HST) is a compressed-gas piston-driven shock tube experiment. The second, the “Multishock thin layer” (Mshock) experiment performed at the National Ignition Facility, is a micrometer/nanosecond-scale laser-driven shock tube experiment. Both are situations in which a heavy plane layer (a ‘curtain’) is initially suspended in a lighter medium. After being shocked from at least one side, the layer translates while its interfaces evolve due to the excitation of the Richtmyer-Meshkov instability at its surfaces. The evolution of density variance, which initially exists only on the surface of the layer, as it comes to encompass the whole layer interior is used as a description of layer fragmentation and dissolution. These experiments have each been simulated in the Los Alamos National Laboratory multi-physics code xRAGE, which includes fundamental hydrodynamics, extended plasma physics and radiation effects which are important to drive the high-energy density experiment, and the Besnard-Harlow-Rauenzahn (BHR) turbulence model. In each, the principal diagnostic for comparison is an experimental metric for the density (co)variance, b, which tracks the moments of the density field at the curtain interfaces and body. Due to experimental constraints in different regimes (i.e. optical diagnostics can be deployed on conventional shock tubes, while the plasma shock tubes must be imaged by x-rays; interfaces can be imposed to specification on laser-driven experiments, which are stored in the solid phase, while conventional experiments have imperfect control of the flow fields which separate the layer, etc.) the experiments are not designed to be perfect scaled cognates of one another. However, despite the separation of six orders of magnitude of scaling in time, and four in space, we are able to demonstrate that the same turbulence model, operating in the same fashion in the same computer code, is able to reproduce results in each experiment, by tracking evolution due to common relevant physics. Additionally, we will present preliminary work toward density variance comparisons in a single-interface Richtmyer-Meshkov configuration, the conventional fluid “Vertical Shock Tube” (VST) experiment, and the Modal Initial Conditions (ModCons) campaign fielded at the OMEGA-EP laser facility.

42 ENGINEERING↗
DOE OSTI2025

Development, Verification, and Validation of an OpenFOAM-Based Solver for Modeling Inertial Fusion Energy Chambers

Our work seeks to introduce a computational tool tailored to the physics of inertial fusion energy chambers, in particular, those concepts based on thick liquid walls. In this approach, the structural materials are protected by several neutron mean-free-paths of renewable liquid and thus will be able to survive much longer than un-shielded walls, with virtually all structures lasting for the life of the plant and enabling the use of commercially available and qualified materials. The OpenFOAM-based solver named rhoCentralFoam has been used as a starting point. rhoCentralFoam belongs to the standard OpenFOAM solver toolset. It is a high-speed, explicit compressible flow solver with shock-capturing capability. While the main features have been retained, the solver had to be restructured to make use of tabular data for equations of states, a necessary addition to model the complex thermo-physical properties of ionized gasses. This entailed the need to change the independent state variables used by the solver, resulting in a new thermodynamic library and slightly different solution algorithm. Moreover, a radiation heat transfer model based on the P-1 approximation was added to the solver. The solver is verified against an analytical solution from the Sedov-Taylor-Neumann test problem to showcase the ability of the hydrodynamic solvers to handle strong shocks, whereas the P-1 model was verified using a simple one-dimensional problem with an analytical solution. Additionally, a validation case involving shock-wave propagation through jet array is presented, and the results are compared with experimental data from the open literature. Lastly, in order to showcase the utility of the solver for practical cases, we applied the refined solver to two representative scenarios: gas venting within the HYLIFE-II chamber and the compression of the gas following the partial ablation of the liquid wall.

Chamber dynamics↗
DOE OSTI2022

Reactions of hydrazine with the amidogen radical and atomic hydrogen

The rate coefficient k 1 for NH 2 + N 2 H 4 was measured to be (5.4 ± 0.4) × 10 –14 cm 3 molecule –1 s –1 at 296 K. NH 2 was generated by pulsed laser photolysis of NH 3 at 193 nm, and monitored as a function of time by pulsed laser-induced fluorescence excited at 570.3 nm under pseudo-first order conditions in the presence of excess N 2 H 4 in an Ar bath gas. This reaction was also investigated computationally, with geometries and scaled frequencies obtained with M06-2X/6-311+G(2df,2p) theory, and single-point energies from CCSD(T)-F12b/cc-pVTZ-F12 theory, plus a term to correct approximately for electron correlation through CCSDT(Q). Three connected transition states are involved and rate constants were obtained via Multistructural Improved Canonical Variational Transition State Theory with Small Curvature Tunneling. Combination of experiment and theory leads to a recommended rate coefficient for hydrogen abstraction of k 1 = 6.3 × 10 –23 T 3.44 exp(+289 K/T) cm 3 molecule –1 s –1 . The minor channel for H + N 2 H 4 forming NH 2 + NH 3 was characterized computationally as well, to yield 5.0 × 10 –19 T 2.07 exp(-4032 K/T) cm 3 molecule –1 s –1 . Finally, these results are compared to several discordant prior estimates, and are employed in an overall mechanism to compare with measurements of half-lives of hydrazine in a shock tube.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Damage to Polystyrene Polymer Film by Shock Wave Induced Bubble Collapse

Metallic surfaces that are in contact with solutions are commonly used in numerous applications where these surfaces can be damaged by shock wave induced bubble collapse. Use of polymer films that coat such surfaces to prevent them from damage requires a better understanding of how much harm collapsing bubbles produce in the films. In this study, we report the results from coarse-grained molecular dynamics simulations to study the damage to polystyrene (PS) films coating a hard surface. The damage was caused by a collapsing nanobubble located in the proximity of the film and interacting with an impinging shock wave. This collapse produces a high-speed water jet that impacts the PS film with a greater force than the shock front and creates cavities/pits in the PS film. We observed that polymer molecules located in the jet vicinity undergo conformational extension in the direction perpendicular to the jet motion, while chain molecules in the rest of the film undergo compression. Furthermore, we also observed that damage to the film is sensitive to the strength of the shock wave.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Scaling laser preheat for MagLIF with the Z-Beamlet laser

Optimizing the performance of the Magnetized Liner Inertial Fusion (MagLIF) platform on the Z pulsed power facility requires coupling greater than 2 kJ of preheat energy to an underdense fuel in the presence of an applied axial magnetic field ranging from 10 to 30 T. Achieving the suggested optimal preheat energies has not been experimentally achieved so far. In this work, we explore the preheat design space for cryogenically cooled MagLIF targets, which represent a viable candidate for increasing preheat energies. Using 2D and 3D HYDRA MHD simulations, we first discuss the various physical effects that occur during laser preheat, such as laser energy deposition, self-focusing, and filamentation. After identifying the changes that different phase plates, gas-fill densities, and magnetic fields bring to the aforementioned physical effects, we, then, consider higher laser energies that are achievable with modest upgrades to the Z Beamlet laser. Lastly, with a 6.0-kJ upgraded laser, 3D calculations suggest that it is possible to deliver 4.25 kJ into the MagLIF fuel, resulting in an expected deuterium neutron yield of Y DD ≃ 1.5 × 10 14 , or roughly 50 kJ of DT equivalent yield, at 20-MA current drive. This represents a 10-fold increase in the currently achieved yields for MagLIF.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗