Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamic codes”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Impact of flow-induced beam deflection on beam propagation in ignition scale hohlraums

Experiments examining the amount of specular reflection (or “glint”) within hohlraums containing different gas fill densities have recently been performed. Simulations of these experiments are presented that show using a single flux limiter cannot explain the decrease in glinted power with increasing gas fill density. The hypothesis that flow-induced beam deflection alters laser absorption is presented. A model is proposed that can be implemented into a ray tracing description of the laser commonly used in radiation hydrodynamic codes. It is shown that simulations using this model capture the trend with gas fill density improving agreement with measurements. This formulation is then applied to an ensemble of laser-driven inertial confinement fusion experiments performed at the National Ignition Facility. The proposed model shows little impact on the total x-ray drive on the capsule but a large impact on the resulting implosion symmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Time-dependent density-functional theory study on nonlocal electron stopping for inertial confinement fusion

Understanding laser–target coupling is of the utmost importance for achieving high performance in laser-direct-drive (LDD) inertial confinement fusion (ICF) experiments. Thus, accurate modeling of electron transport and deposition through ICF-relevant materials and conditions is necessary to quantify the total thermal conduction and ablation. The stopping range is a key transport quantity used in thermal conduction models; in this work, we review the overall role that the electron mean free path (MFP) plays in thermal conduction and hydrodynamic simulations. The currently used modified Lee–More model employs various physics approximations. We discuss a recent model that uses time-dependent density functional theory (TD-DFT) to eliminate these approximations in both the calculation of the electron stopping power and corresponding MFP in conduction zone polystyrene (CH) plasma. In general, the TD-DFT calculations showed a larger MFP (lower stopping power) than the standard modified Lee–More model. Using the TD-DFT results, an analytical model for the electron deposition range, λTD−DFT(ρ,T,K), was devised for CH plasmas between ρ=[0.05−1.05] g/cm3, kBT=[100−1000] eV. We implemented this model into LILAC, for simulations of a National Ignition Facility-scale LDD implosion and compared key physics quantities to ones obtained by simulations using the standard model. The implications of the obtained results and the path moving forward to calculate this same quantity in conduction-zone deuterium–tritium plasmas are further discussed, to hopefully close the understanding gap for laser target coupling in LDD-ICF simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Numerical simulations of laser-driven experiments of ion acceleration in stochastic magnetic fields

We present numerical simulations used to interpret laser-driven plasma experiments at the GSI Helmholtz Centre for Heavy Ion Research. The mechanisms by which non-thermal particles are accelerated in astrophysical environments, e.g., the solar wind, supernova remnants, and gamma ray bursts, is a topic of intense study. When shocks are present, the primary acceleration mechanism is believed to be first-order Fermi, which accelerates particles as they cross a shock. Second-order Fermi acceleration can also contribute, utilizing magnetic mirrors for particle energization. Despite this mechanism being less efficient, the ubiquity of magnetized turbulence in the universe necessitates its consideration. Another acceleration mechanism is the lower-hybrid drift instability, arising from gradients of both density and magnetic field, which produce lower-hybrid waves with an electric field that energizes particles as they cross these waves. With the combination of high-powered laser systems and particle accelerators, it is possible to study the mechanisms behind cosmic-ray acceleration in the laboratory. In this work, we combine experimental results and high-fidelity three-dimensional simulations to estimate the efficiency of ion acceleration in a weakly magnetized interaction region. We validate the FLASH magneto-hydrodynamic code with experimental results and use OSIRIS particle-in-cell code to verify the initial formation of the interaction region, showing good agreement between codes and experimental results. We find that the plasma conditions in the experiment are conducive to the lower-hybrid drift instability, yielding an increase in energy ΔE of ~ 264 keV for 242 MeV calcium ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Shock-driven three-fluid mixing with various chevron interface configurations

When a shock wave crosses a density interface, the Richtmyer–Meshkov instability causes perturbations to grow. Richtmyer–Meshkov instabilities arise from the deposition of vorticity from the misaligned density and pressure gradients at the shock front. In many engineering applications, microscopic surface roughness will grow into multi-mode perturbations, inducing mixing between the fluid on either side of an initial interface. Applications often have multiple interfaces, some of which are close enough to interact in the later stages of instability growth. In this study, we numerically investigate the mixing of a three-layer system with periodic zigzag (or chevron) interfaces, calculating the dependence of the width and mass of mixed material on properties such as the shock timing, chevron amplitude, multi-mode perturbation spectrum, density ratio, and shock mach number. The multi-mode case is also compared with a single-mode perturbation. The Flash hydrodynamic code is used to solve the Euler equations in three dimensions with adaptive grid refinement. Key results include a significant increase in mixed mass when changing from a single-mode to a multi-mode perturbation on one of the interfaces. The mixed width is mainly sensitive to the density ratio and chevron amplitude, whereas the mixed mass also depends on the multi-mode spectrum. In conclusion, steeper initial perturbation spectra have lower mixed mass at early times but a greater mixed mass after the reflected shock transits back across the layer.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Shock tube simulations for the three-layer Richtmyer–Meshkov instability with single-mode and multimode perturbations

While the canonical two-component, single-mode Richtmyer–Meshkov instability (RMI) has been extensively studied, relatively less work has focused on the effects of an additional intermediate-density middle layer. This work investigates such three-material RMI configurations at two Atwood number scenarios using the ARES hydrodynamics code. After validation against previous experimental and computational studies, setups corresponding to recent three-layer shock tube experiments are simulated. Cases with both single-mode and multimode perturbations are studied to quantify mixing across the interface between the materials with highest and intermediate density. In particular, this work is able to comprehensibly examine differences between two- and three-dimensional setups for the single-mode and multimode problems. Observations from previous two-layer investigations still apply in the three-layer setup, but over the time horizons considered, there appears to be insufficient nonlinear mode coupling to create significant differences between two- and three-dimensional simulations following the first passage of a shock. Finally, additional reshock simulations have additional nonlinear growth that does result in expected differences between two- and three-dimensional cases in this three-layer setup, but significant differences do not manifest during the time horizon studied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

A future of inertial confinement fusion without laser-plasma instabilities

From the beginning of inertial confinement fusion (ICF) research, laser-plasma instabilities excited by narrowband lasers have limited the hydrodynamic design space of all laser-based approaches to inertial fusion energy (IFE). With advances in broadband laser technologies, the next generation of ICF drivers will likely have large bandwidth that is engineered to mitigate laser-plasma instabilities, thereby expanding the hydrodynamic design space to include both robust high yields (>200−MJ) with large-energy laser systems (>4 MJ) and high gains (>10) with moderate-energy laser facilities (<2 MJ). State-of-the-art simulations indicate that laser bandwidths of a few percent are required to mitigate instabilities for IFE-relevant conditions. To test these models and demonstrate that high-bandwidth lasers mitigate laser-plasma instabilities, the Fourth-generation Laser for Ultrabroadband eXperiments will be used with the OMEGA Laser System. The goal is to provide the community with the confidence to invest in a multiple-beam high-bandwidth laser facility that will demonstrate the necessary ablation pressures for robust direct-drive ignition without detrimental levels of hot electrons that degrade the implosion performance. It is important to recognize that the highest performing ICF implosions will likely never be completely LPI free because of the significant advantages to maximizing the laser intensity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Detailed simulations of the first deuterium-tritium-filled double shell implosions on the National Ignition Facility

The first indirectly driven, liquid DT-filled double shell inertial confinement fusion (ICF) implosions have recently been successfully performed on the National Ignition Facility (NIF). Double shells are a class of alternative designs that use a low-Z outer shell to compress a foam cushion that accelerates a high-Z inner shell to efficiently compress a liquid DT core. Double shells are challenging to fabricate, field, and model. Important engineering features enabling double shell fabrication include a fill-tube penetrating all shells and a carefully designed and very narrow (few μm) step-joint in the ablator. Due to the higher density materials involved, high Atwood number instabilities are also important at many material interfaces. In this paper, numerical simulations of double shell implosions using the Los Alamos National Laboratory multi-physics radiation-hydrodynamics code xRAGE will be discussed. An extensive effort has been under way for several years to develop the code capabilities for ICF simulations in a common modeling framework to allow ease of simulation setup and standardization of the computational methodology. This paper will present a wide range of simulation results capturing, quantifying, and comparing the impact of all these degradation mechanisms on implosion performance. Brief comparisons with recent experimental results and suggestions for future improvements will also be discussed. Our results suggest that capsule surface roughness and the step-joint gap have the largest impact on implosion performance. Initial experimental data may suggest that the sensitivity to the step-joint gap could provide the dominant explanation for DT-filled double shell experiments that have been fielded on NIF thus far.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Hot electron preheat effects in two-photon polymerization plastic lattice laser-driven shock-tube experiments on the OMEGA laser

We present simulations of Omega shock tube experiments designed to investigate hot electron preheat effects in 3D-printed, two-photon polymerization (2PP) plastic lattices. Preheat is inferred in the experiments from the expansion of a plastic witness disk embedded in the lattice. Using the Eulerian radiation-hydrodynamics code xRAGE, we model shock propagation and preheat from both radiative and hot electron energy sources to evaluate their relative impact. To simulate the transport of laser-generated hot electrons, the nonlocal electron heat transport model proposed by Schurtz, Nicolaï, and Busquet (SNB) is extended with a hot electron source term and an energy cascade algorithm. We explore how variations in ablator, lattice geometry, and laser drive affect the shock velocity and witness disk expansion. Simulations show that the inclusion of a 5 μm gold layer reduces shock pressure by 60% and shock speeds by 30%–40% but does not significantly reduce the hot electron preheat, and that different lattice geometries lead to enhanced shock velocities—up to 40% faster than in homogeneous foams. However, radiative and conductive preheat from classical mechanisms alone fail to match experiment. By including a hot electron source term, we reproduce experimental observables such as disk expansion rates and spatial radiographic features. We find that a hot electron population corresponding to 4%–8% of the incident laser energy with T hot = 50 keV produces expansion which agrees with the experimental data, suggesting hot electron preheat is the most plausible explanation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Simulation studies of the interaction of laser radiation with additively manufactured foams

The interaction of laser radiation with foams of various porosities and low densities has been the subject of several numerical and experimental studies. In all cases, the modeling of low-Z under-dense foams as uniform gases of equivalent average density using standard radiation-hydrodynamics codes has resulted in heat-front velocities that are considerably faster than those observed experimentally. It has been theoretically conjectured that this difference may be attributed to the breakdown of the foam’s morphology, leading to a dynamics of filament expansion where the ion and electron energy partitions are significantly different from those calculated using the uniform gas model. We found that 3D computer simulations employing a disconnected representation of the foam’s microstructure which allowed for the dynamics of foam element heating, expansion, and stagnation largely supported the theoretical picture. Simulations using this model for laser experiments on under-dense 2 mg cc -1 SiO 2 aerogel foams reproduced the experimental data fairly well. We used the validated model in simulations of low-density structured foam-like materials (produced via additive manufacturing) with a variety of morphologies. We found that the log-pile configurations were consistent with the analytical propagation model of Gus’kov et al. Further validation of the model was obtained by simulating experiments performed at the Jupiter Laser Facility using the log-pile and octet-truss foam morphologies. Simulations of the foam-laser interaction using a wave propagation code showed that the microstructure was able to enhance stimulated Brillouin scattering (SBS) by concentrating the light energy into density holes. In turn, this promotes laser filamentation, reducing SBS and bringing the predicted values closer to the experimental data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The collapse and three-dimensional explosion of three-dimensional massive-star supernova progenitor models

ABSTRACT The explosion outcome and diagnostics of core-collapse supernovae depend sensitively on the nature of the stellar progenitor, but most studies to date have focused exclusively on one-dimensional, spherically symmetric massive star progenitors. We present some of the first core-collapse supernovae simulations of three-dimensional massive star supernovae progenitors, a 12.5- and a 15-M⊙ model, evolved in three dimensions from collapse to bounce through explosion with the radiation-hydrodynamic code fornax. We compare the results using those starting from three-dimensional progenitors to three-dimensional simulations of spherically symmetric, one-dimensional progenitors of the same mass. We find that the models evolved in three dimensions during the final stages of massive star evolution are more prone to explosion. The turbulence arising in these multidimensional initial models serves as seed turbulence that promotes shock revival. Detection of gravitational waves and neutrinos signals could reveal signatures of pre-bounce turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

flash -light on the ring : hydrodynamic simulations of expanding supernova shells near supermassive black holes

ABSTRACT The way supermassive black holes (SMBHs) in Galactic Centres (GCs) accumulate their mass is not completely determined. At large scales, it is governed by galactic encounters, mass inflows connected to spirals arms and bars, or due to expanding shells from supernova (SN) explosions in the central parts of galaxies. The investigation of the latter process requires an extensive set of gas dynamical simulations to explore the multidimensional parameter space needed to frame the phenomenon. The aims of this paper are to extend our investigation of the importance of SNe for inducing accretion on to an SMBH and carry out a comparison between the fully hydrodynamic code flash and the much less computationally intensive code ring, which uses the thin shell approximation. We simulate 3D expanding shells in a gravitational potential similar to that of the GC with a variety of homogeneous and turbulent environments. In homogeneous media, we find convincing agreement between flash and ring in the shapes of shells and their equivalent radii throughout their whole evolution until they become subsonic. In highly inhomogeneous, turbulent media, there is also a good agreement of shapes and sizes of shells, and of the times of their first contact with the central 1-pc sphere, where we assume that they join the accretion flow. The comparison supports the proposition that an SN occurring at a galactocentric distance of 5 pc typically drives 1–3 M⊙ into the central 1 pc around the GC.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

r -process nucleosynthesis and kilonovae from hypermassive neutron star post-merger remnants

ABSTRACT We investigate r-process nucleosynthesis and kilonova emission resulting from binary neutron star (BNS) mergers based on a three-dimensional (3D) general-relativistic magnetohydrodynamic (GRMHD) simulation of a hypermassive neutron star (HMNS) remnant. The simulation includes a microphysical finite-temperature equation of state (EOS) and neutrino emission and absorption effects via a leakage scheme. We track the thermodynamic properties of the ejecta using Lagrangian tracer particles and determine its composition using the nuclear reaction network SkyNet. We investigate the impact of neutrinos on the nucleosynthetic yields by varying the neutrino luminosities during post-processing. The ejecta show a broad distribution with respect to their electron fraction Ye, peaking between ∼0.25–0.4 depending on the neutrino luminosity employed. We find that the resulting r-process abundance patterns differ from solar, with no significant production of material beyond the second r-process peak when using luminosities recorded by the tracer particles. We also map the HMNS outflows to the radiation hydrodynamics code SNEC and predict the evolution of the bolometric luminosity as well as broadband light curves of the kilonova. The bolometric light curve peaks on the timescale of a day and the brightest emission is seen in the infrared bands. This is the first direct calculation of the r-process yields and kilonova signal expected from HMNS winds based on 3D GRMHD simulations. For longer-lived remnants, these winds may be the dominant ejecta component producing the kilonova emission.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Convergence of small scale Ly α structure at high- z under different reionization scenarios

ABSTRACT The Ly α forest (LAF) at z > 5 probes the thermal and reionization history of the intergalactic medium (IGM) and the nature of dark matter, but its interpretation requires comparison to cosmological hydrodynamical simulations. At high-z, convergence of these simulations is more exacting since transmission is dominated by underdense voids that are challenging to resolve. With evidence mounting for a late end to reionization, small structures down to the sub-kpc level may survive to later times than conventionally thought due to the reduced time for pressure smoothing to impact the gas, further tightening simulation resolution requirements. We perform a suite of simulations using the Eulerian cosmological hydrodynamics code Nyx, spanning domain sizes of 1.25 − 10 h−1 Mpc and 5 − 80 h−1 kpc cells, and explore the interaction of these variables with the timing of reionization on the properties of the matter distribution and the simulated LAF at z = 5.5. In observable Ly α power, convergence within 10 per cent is achieved for k < 0.1 s km–1, but larger k shows deviation of up to 20 per cent. While a later reionization retains more small structure in the density field, because of the greater thermal broadening there is little difference in the convergence of LAF power between early (z = 9) and later (z = 6) reionizations. We conclude that at z ∼ 5.5, resolutions of 10 kpc are necessary for convergence of LAF power at k < 0.1 s km–1, while higher-k modes require higher resolution, and that the timing of reionization does not significantly impact convergence given realistic photoheating.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Development of an ab initio learned model of electron deposition range in deuterium-tritium plasmas through time-dependent density functional theory calculations and machine learning

Accurate hydrodynamic modeling for laser-direct-drive (LDD) inertial-confinement-fusion (ICF) relies on precise calculations of the electron thermal conduction in all target materials. The nonlocal stopping range of electrons in ICF plasmas directly influences thermal conduction; yet, no first principles model exists for the electron mean free path in the conduction-zone regime. This work utilized time-dependent stochastic density-functional theory (TD-sDFT) to calculate the electron stopping power in deuterium-tritium (DT) plasmas at (ρ, T) conditions relevant to the conduction zone and the compressed shell in ICF. Using a combination of our TD-sDFT data and already established analytical models, we developed and trained an artificial neural network to create a global model for the nonlocal electron deposition range, λ E . We compared our machine-learning (ML) based model for λ E to the currently-used modified-Lee-More model in LDD radiation-hydrodynamic codes, such as lilac, and saw an overall decrease in the deposition range. To understand the effects of λ E on LDD ICF implosion dynamics, we implemented the ML-based model into lilac; specifically, we looked at designs consistent with a current experiment on the OMEGA laser and for a newly designed LDD-ICF target for the future OMEGA-Next facility. In both cases, we saw an overall drop in predicted ablation pressure, peak areal density, and neutron yield due to the reduced thermal conduction (smaller λ E ) in DT plasmas. Comparisons with the experiment on OMEGA are also made.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Gamma-ray signatures of 𝑟-process radioactivity from the collapse of magnetized white dwarfs

We predict the gamma-ray line emission from 𝑟-process nuclei synthesized in the ejecta of the accretion-induced collapse (AIC) of a magnetized, rapidly rotating white dwarf. Using ejecta from a two-dimensional general-relativistic neutrino-magnetohydrodynamic simulation, further evolved with a radiation-hydrodynamics code coupled to an in situ nuclear reaction network, we construct angle-dependent gamma-ray spectra in the 0.01–10 MeV band via composition-dependent ray tracing through the ejecta. The emission between ∼1 and 10 d is dominated by 132 I (𝑡 1/2 = 2.3 h), continuously replenished by the decay of its parent 132 Te (𝑡 1/2 = 3.2 d), with additional contributions from 131 I, 133 Xe, and 132 Te . At 𝑡 ≳ 20 d, 56 Co (from 56 Ni decay) becomes the primary emitter. The simultaneous presence of 𝑟 process and iron-peak gamma-ray lines is distinctive of AIC ejecta and absent in binary neutron star mergers, where iron-peak nuclei are generally not synthesized. Comparing with the 3⁢𝜎 continuum sensitivities of planned MeV gamma-ray telescopes (COSI, AMEGO-X, e-ASTROGAM, GRAMS, GammaTPC), we find the brightest 𝑟-process lines detectable to ∼10 Mpc by GammaTPC and GRAMS, with the signal approaching their sensitivity threshold at 30 Mpc. As a result, the 𝑟-process spectral features survive time integration over ∼30 d exposures, demonstrating robustness against the long observation times required by gamma-ray detectors.

Nuclear reactions↗
DOE OSTI2022

Model study of the energy dependence of the correlation between anisotropic flow and the mean transverse momentum in Au + Au collisions

A hybrid model that employs the hadron-string transport model UrQMD and the (3+1)-dimensional relativistic viscous hydrodynamic code vHLLE, is used to investigate the beam energy dependence of the correlation coefficient ρ(v$^{2}_{2}$,[p T ]) between the average transverse momentum [p T ] of hadrons emitted in an event and the square of the anisotropic flow coefficient v$^{2}_{2}$. For Au+Au collisions, the model predicts characteristic patterns for the energy and event-shape dependence of the variances for [p T ] and v$^{2}_{n}$ [Var([p T ] and Var(v$^{2}_{2}$)], and the covariance of v$^{2}_{n}$ and [p T ][cov(v$^{2}_{2}$,[p T ])], consistent with the attenuation effects of the specific shear viscosity η/s. Here in contrast, ρ(v$^{2}_{2}$,[p T ]) is predicted to be insensitive to the beam energy but sensitive to the initial-state geometry of the collisions. These observations suggest that a precise set of measurements for Var([p T ]), Var(v$^{2}_{2}$), cov(v$^{2}_{2}$,[p T ]) and ρ(v$^{2}_{2}$,[p T ]) as a function of beam energy and event shape, could serve to constrain better the eccentricities and their fluctuations over the entire span of beam energies, and consequently aid precision extraction of the temperature and baryon chemical-potential dependence of η/s from the wealth of Au+Au data obtained in the RHIC beam energy scan.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Gravitational-wave signal of a core-collapse supernova explosion of a 15M ⊙ star

In this paper we report on the gravitational-wave signal computed in the context of a three-dimensional simulation of a core-collapse supernova explosion of a 15 M ⊙ star. The simulation was performed with our neutrino hydrodynamics code chimera . We detail the gravitational wave strains as a function of time, for both polarizations, and discuss their physical origins. We also present the corresponding spectral signatures. Gravitational wave emission in our model has two key features: low-frequency emission (less than 200 Hz) emanates from the gain layer as a result of neutrino-driven convection and the standing accretion shock instability (SASI), and high-frequency emission (greater than 600 Hz) emanates from the proto–neutron star due to Ledoux convection within it. The high-frequency emission dominates the gravitational wave emission in our model and emanates largely from the convective layer itself, not from the convectively stable layer above it, due to convective overshoot. Moreover, the low-frequency emission emanates from the gain layer itself, not from the proto–neutron star, due to accretion onto it. We provide evidence of the SASI in our model and demonstrate that the peak of our low-frequency gravitational wave emission spectrum corresponds to it. Given its origin in the gain layer, we classify the SASI emission in our model as p-mode emission and assign a purely acoustic origin, not a vortical–acoustic origin, to it. We compare the results of our three-dimensional model analysis with those obtained from the model’s two-dimensional counterpart and find a significant reduction in the strain amplitudes in the former case, as well as significant reductions in all related quantities. Our dominant proto–neutron star gravitational wave emission is not well characterized by emission from surface g modes, complicating the relationship between peak frequencies observed and the mass and radius of the proto–neutron star expressed by analytic estimates under the assumption of surface g-mode emission. We present our frequency normalized characteristic strain along with the sensitivity curves of current- and next-generation gravitational wave detectors. This simple analysis indicates that the spectrum of gravitational wave emission between approximately 20 Hz and approximately 1 kHz, stemming from neutrino-driven convection, the SASI, accretion onto the proto–neutron star, and proto–neutron star convection, will be accessible for a Galactic event.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Direct-drive double-shell implosion: A platform for burning-plasma physics studies

Double-shell ignition designs have been investigated with the indirect-drive inertial confinement fusion (ICF) scheme in both simulations and experiments in which the inner-shell kinetic energy was limited to ~10 to 15 kJ, even driven by megajoule-class lasers such as the National Ignition Facility. Since direct-drive ICF can couple more energy to the imploding shells, we have performed a detailed study on direct-drive double-shell (D 3 S) implosions with state-of-the-art physics models implemented in radiation-hydrodynamic codes (LILAC and DRACO), including nonlocal thermal transport, cross-beam energy transfer (CBET), and first-principles-based material properties. To mitigate classical unstable interfaces, we have proposed the use of a tungsten/beryllium–mixed inner shell with gradient-density layers that can be made by magnetron sputtering. In our D 3 S designs, a 70-μm-thick beryllium outer shell is driven symmetrically by a high-adiabat (α ≥ 10), 1.9-MJ laser pulse to a peak velocity of ~240 km/s. Upon spherical impact, the outer shell transfers ~30 to 40 kJ of kinetic energy to the inner shell filled with deuterium–tritium gas or liquid, giving neutron-yield energies of ~6 MJ in 1-D simulations. Two-dimensional high-mode DRACO simulations indicated that such high-adiabat D 3 S implosions are not susceptible to laser imprint, but the long-wavelength perturbations from the laser port configuration along with CBET can be detrimental to the target performance. Yet, neutron yields of ~0.3- to 1.0-MJ energies can still be obtained from our high-mode DRACO simulations. The robust α-particle bootstrap is readily reached, which could provide a viable platform for burning-plasma physics studies. Once CBET mitigation and/or more laser energy becomes available, we believe that breakeven or moderate energy gain might be feasible with the proposed D 3 S scheme.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗