Learning from Each Other: Cross-Cutting Diagnostic Development Activities Between Magnetic and Inertial Confinement Fusion
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
There are programs for high-Z shell experiments at the National Ignition Facility (NIF). For shells made of actinide material, a quantitative fission diagnostic is needed in order to determine how much fission took place and whether the fission was sufficient to produce a non-negligible heat source in the burning capsule. Here, we present a viable coupled experimental and theoretical technique for making quantitative fission measurements possible. The proposed scheme involves using a well-characterized set of depleted uranium foils outside an NIF capsule to verify the conversion of xenon and krypton fission fragments collected at the Radiochemical Analysis of Gaseous Samples (RAGS) facility into total fission yield. We present the calculations needed for this conversion, including the decays in and out of fission fragment chains during the RAGS pump-down of the NIF chamber.
Considerable progress has been made in deuterium-tritium-layered implosion experiments on the OMEGA Laser System, bringing the prospects for thermonuclear ignition in direct-drive configurations with megajoule-class lasers closer to reality. Doing so has required navigating the balance between improved 1D performance and multidimensional stability. Using statistical modeling based on over 350 cryogenic implosions to identify various degradation mechanisms, and combined with multidimensional simulations and experimental techniques such as target offsets to combat residual flows, core conditions have repeatably been achieved that extrapolate to the burning-plasma state when scaled to 2.15 MJ of symmetric laser illumination. Using high implosion velocities (> 450 km/s) and moderately high adiabats (~5), these experiments produced record-high scaled Lawson parameters in direct drive equal to 89 ± 2% of that required for ignition with expected yields of up to 1.5 ± 0.2 MJ. To improve these results still further, focused physics studies are performed to improve physics understanding and identify routes to even greater performance. Recent studies include investigations into the impact of mounting features, laser imprint, reduced fuel temperatures, and greater on-target intensities through subscale experiments. This manuscript gives a summary of the cryogenic direct-drive program on the OMEGA laser, including routes taken to achieve the current best performance, the status of recent focused physics investigations, and future designs—such as target solutions to laser imprint and reducing vapor density to increase convergence—that are expected lead to the demonstration of hydro-equivalent ignition on OMEGA.
The Magnetized Liner Inertial Fusion (MagLIF) experimental platform at Sandia National Laboratories has realized a number of recent improvements in electrical current and laser preheat coupling but fusion yields are below expectations, based on modeling with clean (that is, simulations with no interfacial instabilities and no mix) two-dimensional (2D) magneto-hydrodynamics (MHD) simulations, by a factor of 5 or more. However, magnetized liner implosions are known to produce helical magneto-Rayleigh–Taylor (MRT) modes and complex stagnation structures that cannot be directly modeled in 2D. This paper presents the results of 3D HYDRA MHD simulations, including helical MRT, showing that degradation from these instabilities can readily reproduce experimental yields, but it is difficult to simultaneously match the yield along with inferred fuel temperature, pressure, and burn history (assessed with x-rays), which are larger in the experiments considered here. Additional analysis methods and future experiments are proposed to help address the discrepancies. The 3D simulations also show, without substantive improvements to stability, changing the applied Bz or preheat is unlikely to increase performance of MagLIF to the same degree as clean 2D simulations. Finally, the first 3D HYDRA simulations with the Hall term show the ab initio production of helical MRT and produce a clear change in stagnation morphology compared to pre-seeded simulations without Hall. In spite of the differences, the simulations still produce comparable fusion performance but with longer, lower power, neutron yield history.
An Arrhenius–Wescott–Stewart–Davis (AWSD) reactive flow model for high explosive PBX 9404 is developed. We specifically calibrate an AWSD model for PBX 9404 by fitting equations of state for reactants and detonation products to the results of thermochemical calculations and to experimental data from multiple sources. The calibrated equations of state are then coupled with an Arrhenius rate law based on shock temperature that describes the reaction progress during PBX 9404 detonation. The parameter values in the rate law are calibrated to experimental gas-gun data and diameter effect results. The results of the calibrated AWSD model are in strong agreement with available experimental data for PBX 9404. A similar level of agreement between predicted and experimental results is observed when the calibrated model is validated on data that were not used in the model parameterization procedure. Our results illustrate that the AWSD model is capable of accurately describing the many important properties and observables in the reactive burn of PBX 9404. Because of the historical significance of PBX 9404 in high explosives research and its current use in aging studies, this work provides an important model of a legacy material, which can be used to make comparisons to new high explosive formulations.
Magnetized liner inertial fusion (MagLIF) experiments have demonstrated fusion-relevant ion temperatures up to 3.1 keV and thermonuclear production of up to 1.1 × 1013 deuterium–deuterium neutrons. This performance was enabled through platform development that provided increases in applied magnetic field, coupled preheat energy, and drive current. Advanced coil designs with internal reinforcement enabled an increase from 10 to 20 T. An improved laser pulse shape, beam smoothing, and thinner laser entrance foils increased preheat energy coupling from less than 1 to 2.3 kJ. A redesign of the final transmission line and load region increased peak load current from 16 to 20 MA. The wider range of input parameters was leveraged to study target performance trends with preheat energy, applied magnetic field, and peak load current. Ion temperature and neutron yield generally followed trends in two-dimensional clean Lasnex calculations. Stagnation performance improved with peak load current when other input parameters were also increased such that convergence was maintained. This dataset suggests that reducing convergence to less than 30 would improve predictability of target performance. Lasnex was used to identify a simulation-optimized scaling path, which suggests 10+ kJ of fusion yield is possible on the Z facility with achievable input parameters. This path also indicates >10 MJ could be generated through volume burn on a future facility with a path to high yield (>200 MJ) using cryogenic dense fuel layers. The newly developed MagLIF platform enables exploration of both this simulation optimized scaling path and a recently developed similarity-scaling path.
Magnetized liner inertial fusion (MagLIF) implosions on the Z accelerator require magnetization of the fuel to reduce thermal conduction losses of energy from the laser-preheated deuterium fusion fuel to the surrounding liner material. While external field coils traditionally used to axially magnetize MagLIF targets are limited to 10–20 T (or perhaps ∼ 30 T with technological development), calculations suggest that higher axial magnetic fields would improve thermal insulation of the fuel and improve MagLIF target performance. We present results from three-dimensional radiation-magnetohydrodynamic simulations of MagLIF implosions employing auto-magnetizing helical liners composed of discrete metallic helical conduction paths separated by electrically insulating material. These simulated auto-magnetizing (AutoMag) liners produce axial magnetic fields >30 T inside the fusion fuel prior to implosion. Simulations indicate that higher thermonuclear yields and burn-averaged fuel ion temperatures are attained for implosions using AutoMag liners compared to standard MagLIF implosions. Comparable implosion morphology is evident in synthetic x-ray images of AutoMag–MagLIF and standard MagLIF implosions.
Typically, fluid simulations are used for tokamak divertor design. However, fluid models are only valid if the scrape-off layer (SOL) is highly collisional. This assumption is valid in many present-day experiments but is questionable in the upstream SOL of some high-power scenarios envisioned for burning plasmas and fusion pilot plants. This paper reports on comparisons between fluid and kinetic simulations of the SOL for upstream parameters and geometry representative of the Spherical Tokamak for Energy Production fusion pilot plant. The SOLPS-ITER (fluid) and Gkeyll (gyrokinetic) codes are operated in a two-dimensional axisymmetric mode, which replaces turbulence with ad-hoc diffusivities. In kinetic simulations, we observe that the ions in the upstream SOL experience significant mirror trapping. This substantially increases the upstream temperature and has important implications for impurity dynamics. We show that the mirror force, which is excluded in SOLPS’s fluid equations, enhances the electrostatic potential drop along the field line in the SOL. We also show that the assumption of equal main ion and impurity temperatures, which is made in commonly used fluid codes, is invalid for the regimes explored here. The combination of these effects results in superior confinement of impurities to the divertor region in kinetic simulations, consistent with our earlier predictions [Kotschenreuther et al., in 29th IAEA 29 Fusion Energy Conference (IAEA, London, UK, 2023)]. This effect can be dramatic, reducing the midplane impurity density by orders of magnitude. These results indicate that in lower collisionality SOL’s the tolerable downstream impurity densities may be higher than would be predicted by fluid simulations, allowing for higher radiated power while avoiding unacceptable core contamination. Our results highlight the importance of kinetic simulations for divertor design and optimization for fusion pilot plants.
Using a nondestructive, in situ near-infrared defect imaging technique, we obtained time-resolved images of dark line defect (DLD) growth in the cavities of several high-power 790-nm diode lasers and calculated individual DLD growth rates. A clear longitudinal spatial dependence in the DLD growth rate was observed, in which the DLD growth rate was faster near the anti-reflective (AR) facet. Longitudinal spatial hole burning simulation results indicate that the total photon density increases toward the AR facet, while the carrier density decreases toward the AR facet. Here, these results suggest that photon absorption plays a key role in enhancing DLD growth beyond the contribution of the traditional carrier capture mechanism. Asymmetric DLD growth that was preferentially directed toward the high-reflectivity facet was also observed, providing further support for the role of photon absorption in the growth process. To account for these findings, we propose a photon-accelerated mechanism for the enhancement of DLD growth in which the slower carrier capture process in the traditional recombination-enhanced dislocation glide mechanism is replaced by a faster photon absorption process.
Analytic studies and two-dimensional “clean” radiation-magnetohydrodynamic (rad-MHD) simulations employing dynamical similarity driver-target scaling prescriptions [Ruiz et al., Phys. Plasmas 30, 032708 (2023)] suggest that Magnetized Liner Inertial Fusion (MagLIF) target implosions can scale to > 10 MJ DT fusion yields when peak drive current is increased beyond 60 MA. We present results from three-dimensional (3D) rad-MHD simulations of similarity-scaled MagLIF target implosions at peak drive currents ranging from 15 to 40 MA. Simulations in this study suggest that magneto-Rayleigh–Taylor instability (MRTI) growth and feedthrough to the fuel region are more severe at higher drive current scales, which reduces the fusion yield compared to prior analytic and 2D clean simulation predictions. In contrast to standard MagLIF, simulations of current-scaled MagLIF target implosions driven by a dynamic screw pinch (DSP) demonstrate reduced MRTI feedthrough and greater fuel magnetization, resulting in improved thermonuclear performance and enhanced performance scaling with peak drive current. DSP drive enables additional scaling of the liner mass to increase liner radius but maintain implosion time, resulting in higher implosion velocities at the expense of increased susceptibility to MRTI. We present a current- and mass-scaled simulated DSP-MagLIF target implosion at the ∼ 40 MA peak current level that produces ignition scale performance, demonstrating a burn-averaged Lawson ignition parameter above unity and DT fusion yield above 1 MJ.
Motivated by the need for accurate, timely, and efficient calculations of plasma transport, predictions of plasma turbulence properties made using different TGLF saturation rules are benchmarked against corresponding predictions from linear and nonlinear gyrokinetic CGYRO simulations. This benchmarking is carried out using parameters taken from an inductive burning plasma scenario in a hypothetical compact high-field (R maj = 4 m, B T = 8 T) tokamak, lying in a much different regime of parameter space than either the TGLF calibration regime or current-day experiments. The core turbulent transport in this scenario is predicted to be dominated by ion temperature gradient (ITG) turbulence. In general, the ITG critical gradients predicted by various TGLF saturation rules are quite close to the CGYRO predictions. Both codes predict similar linear ITG growth rates and frequency spectra, as well as their scaling with R/L T i = −Rd ln(T i )/dr. However, TGLF systematically predicts unstable trapped-electron modes (TEMs) above k y ρ s ≃ 0.5 not seen by CGYRO for the same parameters, due to TGLF predicting a lower threshold in R/L T e than CGYRO for TEM onset. It is shown that for this scenario, nonlinear CGYRO simulations predict stiffer ITG turbulence than the TGLF SAT0 and SAT1 saturation rules, with energy fluxes close in magnitude and scaling with R/L T i to what is predicted by the SAT2 saturation rule. Self-consistent core profiles calculated using nonlinear CGYRO flux predictions and the PORTALS transport solver are shown to agree fairly well with corresponding predictions made using the TGLF SAT2 model, including a similar level of density peaking.
Proton-boron 11 (pB11) fusion is relatively safe and clean, but difficult to use for net power production, since bremsstrahlung radiation tends to radiate away power more quickly than it can be generated by fusion power, particularly once poisoning by alpha particles is taken into account. While in magnetic confinement fusion, this problem can be addressed by deconfining the alphas, in inertial confinement fusion (ICF), the alphas that heat the plasma linger for the duration of the reaction. Thus, it becomes essential to trap the bremsstrahlung radiation in the hotspot. Through burn simulations incorporating bremsstrahlung emission and reabsorption, we infer the necessary conditions to capture enough radiation to produce scientific breakeven in a pB11 ICF plasma. We find that breakeven requires a stagnation areal density roughly two orders of magnitude higher than the current state-of-the-art, at pressures three orders of magnitude higher.
The double shell path to inertially confined volumetric thermonuclear burn relies on efficient kinetic energy transfer from an outer ablator shell to an inner shell containing DT fuel at implosion velocities near 250 km/s. In order to experimentally realize adequate compression, target imperfections and assembly features must be tightly controlled. Pulse shape and ablator material are critical design choices that impact all aspects of the implosion. In this article, we evaluate three designs: an Al ablator with single-shock radiation drive, a CH ablator with single-shock drive, and a CH ablator with reduced-adiabat 2-shock drive. All designs use the same W inner (pusher) shell overcoated with a glow discharge polymer tamper and filled with liquid-density DT. One-dimensional HYDRA simulations revealed that the plastic ablator design required the use of the 2-shock drive to match the performance of the Al design. Two-dimensional simulations resolving broadband surface roughness up to modes of several hundred placed only on the ablator showed plastic to greatly outperform Al. The plastic ablator designs did exhibit higher levels of inner shell instability growth when surface roughness was only placed at the pusher/tamper interface, which was mitigated with higher foam cushion densities. A unique and additional benefit of a plastic ablator is also discussed, where heating of the ablator by Au M-band radiation is predicted to rapidly expand the hemi-shell surfaces, closing the equatorial joint gap ahead of the first shock. Taken together, the properties of the plastic ablator for indirect-drive double shell implosions appear promising in removing several long-standing fabrication challenges.
The double shell is a volume-burn inertial confinement fusion concept consisting of two concentric shells: a low-Z outer shell that collides with and transfers momentum to a high-Z inner shell which compresses and heats the thermonuclear fuel. The increased number of capsule interfaces and severe hydrodynamic instability of the high-density pusher during its acceleration phase provide challenges to the success of the double shell. Two-dimensional radiation-hydrodynamics simulations predict the hydrodynamic instability growth on the outer surface of the pusher can be greatly reduced through appropriate timing of two shocks that cross this interface. One of these shocks, unique to multi-shell designs, arises from radiation-driven ablation of the inner shell ahead of the main shock, the second shock of concern. The shock timing is optimized by increasing the thickness of a low-Z tamper layer exterior to the pusher, resulting in only minimal changes to the implosion timing. Reducing the instability growth on the outer surface of the high-Z pusher can dramatically decrease the modulations that feedthrough to the pusher inner surface, improving the efficacy of radiation trapping in the thermonuclear fuel and increasing the predicted thermonuclear yield by ≳20×.
We examine multimode diode laser pump depletion in a graded index multimode (MM) fibre Raman laser generating a high-quality Stokes beam. Spatial ‘hole burning’ in the transmitted pump beam is demonstrated. A balance model is constructed for interaction of the pump and Stokes waves in the MM Raman laser. Analytical formulas are derived and used in calculations demonstrating qualitative agreement with experimental data. The origin of the observed quantitative discrepancies is discussed. (paper)
Background Vegetation fire may change Phosphorus (P) cycling in terrestrial ecosystems through converting biomass into fire residues. Aim The aim of this study was to understand the chemistry and mobility of P in fire residues to help reveal P thermochemistry during biomass burning and post-fire P cycling. Methods A combination of sequential extraction, liquid 31P NMR and P K-edge XANES was used to obtain quantitative P speciation and explain P solubilisation behaviours of charcoal. Key results Despite varying diverse P species existing in raw biomass, only two P structural moieties – orthophosphate and pyrophosphate – were identified in charcoal. However, relative abundance of pyrophosphate differs greatly among charcoal samples from different biomass types, ranging between 0 and 40% of total extractable P. Although P K-edge XANES data indicates abundant soluble phosphate minerals, most P (70–90%) is likely occluded physically in the charcoal. The bicarbonate-extractable P (the Olsen-P) varies significantly and cannot be explained by surface P concentration or elemental stoichiometry alone. Conclusion and implications The results suggest the importance of starting biomass P speciation (i.e. molecular structure and complexation environment) and thermal conditions in controlling P speciation and availability in charcoal. The different P chemistry between charcoal and ash suggests the importance of fire types and severity in disturbing the P cycle.
Soot emissions in combustion are unwanted consequences of burning hydrocarbon fuels. The presence of soot during and following combustion processes is an indication of incomplete combustion and has several negative consequences including the emission of harmful particulates and increased operational costs. Efforts have been made to reduce soot production in combustion engines through utilizing oxygenated biofuels in lieu of traditional nonoxygenated feedstocks. The ongoing Co-Optimization of Fuels and Engines (Co-Optima) initiative from the US Department of Energy (DOE) is focused on accelerating the introduction of affordable, scalable, and sustainable biofuels and high-efficiency, low-emission vehicle engines. The Co-Optima program has identified a handful of biofuel compounds from a list of thousands of potential candidates. In this study, a shock tube was used to evaluate the performance of soot reduction of five high-performance biofuels downselected by the Co-Optima program. Current experiments were performed at test conditions between 1,700 and 2,100 K and 4 and 4.7 atm using shock tube and ultrafast, time-resolve laser absorption diagnostic techniques. The combination of shock heating and nonintrusive laser detection provides a state-of-the-art test platform for high-temperature soot formation under engine conditions. Soot reduction was found in ethanol, cyclopentanone, and methyl acetate; conversely, an α-diisobutylene and methyl furan produced more soot compared to the baseline over longer test times. For each biofuel, several reaction pathways that lead towards soot production were identified. The data collected in these experiments are valuable information for the future of renewable biofuel development and their applicability in engines.
Significance We provide the first assessment of aboveground live tree biomass in a mixed conifer forest over the late Holocene. The biomass record, coupled with local Native oral history and fire scar records, shows that Native burning practices, along with a natural lightning-based fire regime, promoted long-term stability of the forest structure and composition for at least 1 millennium in a California forest. This record demonstrates that climate alone cannot account for observed forest conditions. Instead, forests were also shaped by a regime of frequent fire, including intentional ignitions by Native people. This work suggests a large-scale intervention could be required to achieve the historical conditions that supported forest resiliency and reflected Indigenous influence.