Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FUSION”

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 1,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Neutron time of flight (nToF) detectors for inertial fusion experiments

Neutrons generated in Inertial Confinement Fusion (ICF) experiments provide valuable information to interpret the conditions reached in the plasma. The neutron time-of-flight (nToF) technique is well suited for measuring the neutron energy spectrum due to the short time (100 ps) over which neutrons are typically emitted in ICF experiments. By locating detectors 10s of meters from the source, the neutron energy spectrum can be measured to high precision. Here, we present a contextual review of the current state of the art in nToF detectors at ICF facilities in the United States, outlining the physics that can be measured, the detector technologies currently deployed and analysis techniques used.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Required toroidal confinement for fusion and omnigeneity

Deuterium–tritium (DT) burning requires a long energy confinement times compared to collision times, so the particle distribution functions must approximate local Maxwellians. Non-equilibrium thermodynamics is applicable, which gives relations among transport, entropy production, the collision frequency, and the deviation from a Maxwellian. The distribution functions are given by the Fokker–Planck equation, which is an advection–diffusion equation. A large hyperbolic operator, the Vlasov operator with the particle trajectories as its characteristics, equals a small diffusive operator, the collision operator. The collisionless particle trajectories would be chaotic in stellarators without careful optimization. This would lead to rapid entropy production and transport—far beyond what is consistent with a self-sustaining DT burn. Omnigeneity is the weakest general condition that is consistent with a sufficiently small entropy production associated with the thermal particle trajectories. Omnigeneity requires that the contours of constant magnetic field strength be unbounded in at least one of the two angular coordinates in magnetic surfaces and that there be a symmetry in the field-strength wells along the field lines. Even in omnigenous plasmas, fluctuations due to microturbulence can produce chaotic particle trajectories and the gyro-Bohm transport is seen in many stellarator and tokamak experiments. The higher the plasma temperature above 10 keV, the smaller the transport must be compared to gyro-Bohm for a self-sustaining DT burn. The hot alphas of DT fusion heat the electrons. When the ion–electron equilibration time is long compared to the ion energy confinement time, a self-sustaining DT burn is not possible, which sets a limit on the electron temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Mitigation of deceleration-phase Rayleigh–Taylor instability growth in inertial confinement fusion implosions

Rayleigh–Taylor growth during shell deceleration is one of the main limiting factors for target performance in inertial confinement fusion implosions. Here, using analytical scaling laws and hydrodynamic simulations, we show that such amplification can be mitigated by reducing the initial mass density in the central target region. The perturbation growth reduction is caused by a smaller hot-spot convergence ratio during deceleration, increased density scale length, and enhanced ablation stabilization. The required central density reduction can be achieved using the dynamic shell formation concept.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Modeling stimulated Brillouin backscatter from outer-cone quads across multiple inertial confinement fusion hohlraum designs

Stimulated Brillouin scattering (SBS) is a potential risk for laser damage in the experiments carried out at the National Ignition Facility (NIF), and by altering the energy deposition pattern in hohlraums, it affects the symmetry of indirect-drive inertial confinement fusion implosions. We have surveyed backward SBS on outer-cone quads across NIF integrated hohlraums of various platforms numerically, using three-dimensional (3D) simulations with the backscatter code pF3D [Berger et al. Phys. Plasmas 5, 4337 (1998)] and ray-based gain calculations. Simulated reflected powers and energies, as well as the spectrum of reflected light all compare favorably with measurements. Ray-based calculations of exponential SBS amplification (“gain”), which assume a strongly damped plasma wave and steady-state response, are performed using a novel method that includes the 3D speckled field of the laser that drives SBS. This approach is useful for understanding qualitative differences between hohlraum designs and identifying regions susceptible to SBS within hohlraums. Quantitatively, gains are not found to correlate with SBS reflectivities in 3D, necessitating fully wave-based calculations that naturally include diffraction and various temporal dependencies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A multifidelity Bayesian optimization method for inertial confinement fusion design

Due to their cost, experiments for inertial confinement fusion (ICF) heavily rely on numerical simulations to guide design. As simulation technology progresses, so too can the fidelity of models used to plan for new experiments. However, these high-fidelity models are by themselves insufficient for optimal experimental design, because their computational cost remains too high to efficiently and effectively explore the numerous parameters required to describe a typical experiment. Therefore, traditionally, ICF design has relied on low-fidelity modeling to initially identify potentially interesting design regions, which are then subsequently explored via selected high-fidelity modeling. In this paper, we demonstrate that this two-step approach can be insufficient: even for simple design problems, a two-step optimization strategy can lead high-fidelity searching toward incorrect regions and consequently waste computational resources on parameter regimes far away from the true optimal solution. We reveal that a primary cause of this behavior in ICF design problems is the presence of low-fidelity optima in different regions of the parameter space far away from high-fidelity optima. To address this issue, we propose an iterative multifidelity Bayesian optimization method based on Gaussian Process Regression that leverages both low- and high-fidelity models simultaneously. We demonstrate, using both two- and eight-dimensional ICF test problems, that our algorithm can effectively utilize both low-fidelity and high-fidelity models to refine the designs. This approach proves to be more efficient than relying solely on high-fidelity modeling for optimization.

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 OSTI2025

Burn propagation in magnetized high-yield inertial fusion

Recent experiments at the National Ignition Facility (NIF) have demonstrated ignition for the first time in an inertial confinement fusion (ICF) experiment, a major milestone allowing the possibility of high energy gain through burn propagation. Use of external magnetic fields, applied primarily to reduce thermal losses, could increase hotspot temperature and ease requirements for ignition, opening up the capsule design space for high energy gain. However, this same restriction of thermal transport has the potential to inhibit burn propagation, which is vital in the attainment of high gain. In this work, radiation-magnetohydrodynamics (MHD) simulations carried out using the code Chimera are used to investigate the effect of a pre-imposed magnetic field on ignition and burn propagation. This paper studies the propagation of burn using both an idealized planar model and in fully integrated 2D MHD simulations of an igniting NIF capsule. A study of magnetized burn propagation in the idealized planar model identifies three regimes of magnetized burn propagation: (1) thermal conduction driven; (2) alpha transport driven; and (3) fully suppressed burn. Simulations of NIF shot N210808 with an applied 40 T axial field show clear indication of burn suppression perpendicular to field lines, with rapid burn observed along field lines. Implosion shape is altered by the field, and anisotropic conduction causes significant modification to the rate of ablation during stagnation. These results highlight the fundamental changes to implosion dynamics in high-yield magnetized ICF and motivate further study to better optimize future magnetized target designs for high gain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

An experimentally informed design process for future inertial confinement fusion facilities

The achievement of ignition in the laboratory has renewed interest in defining the requirements for a future high-gain inertial confinement fusion (ICF) facility. Our best chance of predicting future ICF performance is with 3-D radiation hydrodynamic simulations that have been benchmarked against experimental data, but their high computational cost is prohibitive for use in practical design studies. We introduce a hierarchical approach where 3-D simulations are tuned to match experimental measurements and used to train 3-D degradation models in 1-D simulations allowing for accurate predictions over the entire OMEGA direct-drive database. A genetic algorithm was used in combination with the trained 1-D simulations to search for optimal direct-drive implosion designs at driver energies ranging from 20 kJ to 10 MJ. As the fidelity of 3-D codes improves, this approach will provide a viable experimentally informed tool for defining the next ICF facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Modeling of Transport Processes in Liquid-Metal Fusion Blankets: Past, Present, and Future

The successful development of robust breeding blanket systems will strongly rely on computational tools for predicting the complex behavior of the electrically conducting liquid-metal (LM) breeder flowing in the complex-shaped blanket ducts in the presence of a strong plasma-confining magnetic field, volumetric heating, and tritium generation. Associated transport processes involve magnetohydrodynamic (MHD) flows, heat transfer, corrosion, and tritium transport. This paper is an overview of past and present efforts in the development, application, and verification and validation (V&V) of such computational tools. As a result of the ongoing campaign on V&V of computer codes for LM blankets, the international fusion community has identified several candidates that promise to become real blanket design and analysis tools in the near future. Among them are HIMAG, MHD-UCAS, COMSOL Multiphysics, ANSYS FLUENT, ANSYS CFX, and OpenFOAM. The progress, over the last decade, in the application of such codes in blanket studies is tremendous. This is illustrated with two examples for a dual-coolant lead-lithium (DCLL) blanket: (1) integrated computer modeling for the recently designed DCLL blanket in the United States and (2) application of the code MHD-UCAS to the analysis of PbLi flows and heat transfer in a generic DCLL blanket prototype at high Hartmann (Ha ~ 10 4 ) and Grashof numbers (Gr ~ 10 12 ). Here, this paper also presents an approach to the development of a new integrated computational tool called the virtual dual-coolant lead-lithium (VDCLL) blanket, which elaborates the existing U.S. MHD code HIMAG.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Ultrathick Boron Carbide Coatings for Nuclear Fusion Targets

Boron carbide is an attractive ablator for next-generation inertial confinement fusion (ICF) targets. Here we describe several aspects of our ongoing systematic studies of the deposition and processing of B 4 C coatings for ICF targets. We show that residual compressive stress in films can be reduced and the deposition rate increased by N-doping. Further, we also demonstrate successful Si substrate etching and surface polishing and discuss remaining challenges and offer potential solutions to the buildup of particulates in the deposition chamber during prolonged coating runs, control of nodular growth defects, and lateral nonuniformity of film properties for deposition conditions with relatively low target-to-substrate distances.

36 MATERIALS SCIENCE↗
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 OSTI2026

Fully Additively Manufactured Wetted Foam Capsules for Inertial Confinement Fusion

In the pursuit of realizing reliable clean energy generation via inertial confinement fusion (ICF), wetted foam (WF) capsule targets have long been coveted due to their potential to simplify the target fielding process and suppress hydrodynamic instabilities and material mixing that limit achievable energy output, yet producing and deploying a WF target has proven challenging. Here, in this work, we demonstrate the design, fabrication, metrology, and testing of fully additively manufactured (AM) foam-lined capsules using two-photon polymerization (2PP) for ICF. We successfully fielded an AM polymeric foam capsule with a 3-mm outer diameter, a nominally 15-µm-thick solid outer layer, a 120-µm-thick inner foam layer, and a 250-µm outer diameter copper fill tube on the National Ignition Facility for a polar direct-drive shot, and we showcase deuterium wetting of the capsule foam layer inside an ignition target proofing station. Our exploration showed that 2PP can produce fieldable targets with complex geometries and potentially shorten the design iteration turnaround time and the overall target fabrication time.

Target↗
DOE OSTI2025

Inertial Confinement Fusion Design Search Using Bayesian Optimization

Inertial confinement fusion (ICF) experiments rely on complex multi-physics simulation codes such as the Lawrence Livermore National Laboratory-developed HYDRA to guide design work. However, these simulations have several dozen tunable parameters and can be computationally expensive. This makes searching the parameter space challenging and time-consuming. Recently developed automated tools utilize Bayesian optimization to search these high-dimensional parameter spaces for optimal designs. The optimization tools run 2D integrated simulations in HYDRA to converge on a design that produces specified scalar or vector outputs. In this paper, we apply the Bayesian optimization tools to two common tuning scenarios. First, we tune simulation inputs to match measurements of a well-characterized experiment at the National Ignition Facility. This type of tuning is commonly performed to compensate for the use of simplified simulation settings (e.g. reduced resolution) or to account for missing physics in the simulations. Second, we search for an ICF simulation design that has a particular radiation drive profile. These optimizations replicate the kinds of tuning researchers routinely perform, but do so with significantly reduced manual effort. This approach demonstrates a powerful and efficient pathway toward autonomous, high-fidelity design optimization for future ICF experiments.

Bayesian optimization↗
DOE OSTI2019

Efficient high-order singular quadrature schemes in magnetic fusion

Several problems in magnetically confined fusion, such as the computation of exterior vacuum fields or the decomposition of the total magnetic field into separate contributions from the plasma and the external sources, are best formulated in terms of integral equation expressions. Based on Biot–Savart-like formulae, these integrals contain singular integrands. The regularization method commonly used to address the computation of various singular surface integrals along general toroidal surfaces is low-order accurate, and therefore requires a dense computational mesh in order to obtain sufficient accuracy. In this work, we present a fast, high-order quadrature scheme for the efficient computation of these integrals. Several numerical examples are provided demonstrating the computational efficiency and the high-order accurate convergence. A corresponding code for use in the community has been publicly released.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Review: Turbulence dynamics during the pedestal evolution between edge localized modes in magnetic fusion devices

Fusion performance has been shown to be related to the H-mode pedestal structure. The pedestal is associated with steep gradients that are the source of free energy for microinstabilities. A variety of instabilities have been shown to co-exist in the pedestal. This paper reviews the experimentally observed signatures of instabilities during the pedestal parameters' evolutions, with a focus on turbulence observations, made in between edge localized modes (ELMs), in multiple tokamaks. ELMs are cyclic events associated with bursty relaxations of the pedestal. The multiple machine results point to very similar pedestal localized modes suggesting the existence of a unifying mechanism governing the modes' onset and dynamics. Modeling efforts to identify the generation mechanism of these instabilities are introduced. Several remaining challenges include the study of the instability saturation mechanisms and impact on the pedestal structure formation. Finally, potential future research avenues will require multiscale-non-linear-gyrokinetic analyses to study the non-linear interaction between these instabilities, and which pedestal conditions facilitate such non-linear interactions, as well as the relationship with the formation of the pedestal through transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Normalizing flows for likelihood-free inference with fusion simulations

Fluid-based scrape-off layer transport codes, such as UEDGE, are heavily utilized in tokamak analysis and design, but typically require user-specified anomalous transport coefficients to match experiments. Determining the uniqueness of these parameters and the uncertainties in them to match experiments can provide valuable insights to fusion scientists. Here, we leverage recent work in the area of likelihood-free inference (‘simulation-based inference’) to train a neural network, which enables accurate statistical inference of the anomalous transport coefficients given experimental plasma profile input. UEDGE is treated as a black-box simulator and runs multiple times with anomalous transport coefficients sampled from priors, and the neural network is trained on these simulations to emulate the posterior. The neural network is trained as a normalizing flow model for density estimation, allowing it to accurately represent complicated, high-dimensional distribution functions. With a fixed simulation budget, we compare a single-round procedure to a multi-round approach that guides the training simulations toward a specific target observation. Finally, we discuss the future possibilities for use of amortized models, which train on a wide range of simulations and enable fast statistical inference for results during experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

3D simulations of inertial confinement fusion implosions part 2: systematic flow anomalies and impact of low modes on performances in OMEGA experiments

We present the first 3D radiation-hydrodynamic simulations of directly driven inertial confinement fusion implosions with an inline package for polarized crossed-beam energy transfer, which were used to assess the impact of the current distributed polarization rotators (DPRs) on OMEGA as well as other known sources of asymmetry. Applied to OMEGA implosions, the simulations predict bang times with no need for ad hoc multipliers, as well as yields—if you separately account for the impacts of imprint and fuel age. The magnitude of the flow is well reproduced when the low mode sources are large, whereas the modeling of the stalk is thought to be required to match the flow magnitude in the remaining cases. For the cases explored in more detail, polarized cross-beam energy transfer (CBET)—the only known systematic drive asymmetry, brought the results closest to the measured flow vectors. The remaining discrepancies are shown to be stemming from the limited knowledge of the laser pointing modes. For typical current levels of beam mispointing, power imbalance, target offset, and asymmetry caused by polarized CBET, low modes degrade the yield by more than 40%. The current strategy of attempting to compensate the mode-1 asymmetry with a preimposed target offset recovers only about one-third of the losses caused by the low modes due to the dynamic nature of the multiple asymmetries and the presence of low modes other than l = 1. Furthermore, addressing the root causes of the drive asymmetries is apt to be more beneficial. To that end, one possible solution to the specific issue of polarized CBET (10 µm DPRs) is shown to work well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗