Search NASA⌕ Search

SEARCH · Search NASA

Results for “fusion energy modeling”

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 OSTI2025

Sensitivities of time-dependent temperature profile predictions for NSTX with the multi-mode model

The Multi-Mode Model (MMM) for turbulent transport was applied to a large set of well-analyzed discharges from the National Spherical Torus Experiment (NSTX) in order to evaluate its sensitivities to a wide range of plasma conditions. MMM calculations were performed for hundreds of milliseconds in each discharge by performing time-dependent predictive simulations with the 1.5D tokamak integrated modeling code TRANSP. A closely related study (Lestz et al 2025 Plasma Phys. Control. Fusion 67 105029) concluded that MMM predicted electron and ion temperature profiles that were in reasonable agreement with NSTX observations, generally outperforming a different reduced transport model, TGLF. This finding motivates the more thorough investigation of the characteristics of the MMM predictions conducted in this work. The simulations with MMM have electron energy transport dominated by electron temperature gradient modes in the examined discharges with relatively low plasma β (ratio of kinetic plasma pressure to magnetic field pressure) and high collisionality, transitioning to a mixture of different modes for higher β and lower collisionality. The thermal ion diffusivity predicted by MMM is much smaller than the neoclassical contribution, in line with previous experimental analysis of NSTX. Nonetheless, the electron and ion temperature profiles are coupled via collisional energy exchange and thus sensitive to which transport channels are predicted. The time-dependent simulations with MMM are robust to the simulation start time, converging to remarkably similar temperature profiles later during the discharge. MMM typically overpredicts confinement relative to NSTX observations, leading to the prediction of overly steep temperature profiles. Plasmas with spatially broader temperature profiles, higher plasma β, and longer energy confinement times tend to be predicted by MMM with better agreement with the experiment. As a result, these findings provide useful context for understanding the regime-dependent tendencies of MMM in anticipation of self-consistent, time-dependent predictive simulations of NSTX-U discharges with these same modeling tools.

MMM↗
DOE OSTI2023

Resonance in radio frequency sheath admittance and enhanced impurity emission near the ion cyclotron frequency

Ion cyclotron resonance heating (ICRH) is of considerable interest among all auxiliary heating techniques, because it transfers power directly to ions, targets the high-density core, and involves the cheapest radio-frequency (RF) components. During ICRH operation, RF sheaths form on the ICRH antenna itself, nearby hardware, and far-field surfaces. These sheaths are associated with large hot-spot formation and impurity emissions. This work presents high-resolution numerical modelling of RF sheaths in nuclear fusion scenarios using hPIC2, a Debye-scale particle-in-cell code. The modelling reveals a new RF sheath phenomenon which occurs when the RF is resonant with a species ion cyclotron frequency or harmonic in an oblique magnetic field. This resonance of the RF sheath causes modifications of the energy-angle distributions of the ions impacting on the walls, with consequent increase in wall impurity emission. Due to the 1/ R scaling of the magnetic field in a tokamak, such cases are possible on divertor or vacuum chamber surfaces, depending on the geometry of the tokamak and the RF heating scenario. A simple physics interpretation using a driven damped harmonic oscillator is proposed, where the ion sheath transit time plays the role of damping the RF sheath admittance. Critically, the resonance leads to increased ion heat flux at the wall as well as increased physical sputtering, despite a lack of increase in RF rectified sheath potential.

ICRH antenna↗
DOE OSTI2025

Q2 Report for FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment

This report describes the work and activities carried out towards the completion of each of the following milestones in FY25 Q2: 1. Demonstrate workflow for generating self-consistent CESOL plasma profiles + first wall and divertor loading prediction and generate the CAT plasma and neutron loading needed for further engineering analysis: $\circ$ Run CESOL with BOUT++/Hermes-3 and immersed boundary condition to directly map to wall: • Run BOUT++/Hermes-3 through the IPS workflow to find radial particle and energy diffusivities to match either the Eich or the physics-based scaling of the SOL heat flux width, and • Expand source of first wall heat flux to include charged particles, neutrals, and radiation from the core+edge. 2. Generate medium fidelity parametrized CAD: $\circ$ Develop the TRACER tool to read an existing CAD, regenerate the geometry based on vertex location and connectivity information, define vertex translation and parameters needed for scaling the CAD, and $\circ$ Utilize the FreeGS code to determine CAT PF coil placement, including minimizing the number of coils, coil current, and electromechanical stresses. 3. Utilize plasma loading for engineering analysis: $\circ$ Couple the plasma loading to input for OpenFOAM and demonstrate initial test of thermal analysis of CAT first wall loading with typical DCLL blanket component cooling boundary conditions. 4. Demonstrate nuclear analysis: $\circ$ Apply initial analysis of tritium transport in DCLL blanket by evaluating spatially resolved tritium generation rates, tritium diffusion and convection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

THz Radiation Generation to Enable Internal Magnetic Field Measurement of Burning Plasmas (Final Scientific Report)

A one-year experimental program was performed using the Multi-TeraWatt (MTW) laser at the Laboratory for Laser Energetics (LLE) to investigate the efficiency and physical process of laser-driven—wire THz generation. This process was evaluated as a potential high-repetition-rate source of powerful THz pulses for use in THz pulsed polarimetry of a moderate-density magnetic-confinement fusion concept developed by TAE. The goal was to develop a source for short pulses in the THz regime (150 μm wavelength) containing 10’s of mJ. The campaign tested the impact of Cu wire diameter (25, 75, 160 μm), length (250 – 750 μm), and laser irradiation orientation on THz production efficiency in comparison to 20-μm thick Cu foils. Of these targets, the 160 µm diameter, 750 µm long wires with sharpened tips and irradiated at the base were found to produce the highest energy and shortest wavelength THz pulses, exceeding 2.5 Joule/sr and meeting the requirements for the polarimetry diagnostic. The brightest THz signal was radiated in the direction of the sharpened wire tip. A combination of two models – a traveling wave antenna combined with coherent transition radiation – describes the emission pattern, pulse energy, and spectrum observed in the experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

DIII-D research to provide solutions for ITER and fusion energy

The DIII-D tokamak has elucidated crucial physics and developed projectable solutions for ITER and fusion power plants in the key areas of core performance, boundary heat and particle transport, and integrated scenario operation, with closing the core-edge integration knowledge gap being the overarching mission. New experimental validation of high-fidelity, multi-channel, non-linear gyrokinetic turbulent transport models for ITER provides strong confidence it will achieve Q ≥ 10 operation. Experiments identify options for easing H-mode access in hydrogen, and give new insight into the isotopic dependence of transport and confinement. Analysis of 2,1 islands in unoptimized low-torque IBS demonstration discharges suggests their onset time occurs randomly in the constant β phase, most often triggered by non-linear 3-wave coupling, thus identifying an NTM seeding mechanism to avoid. Pure deuterium SPI for disruption mitigation is shown to provide favorable slow cooling, but poor core assimilation, suggesting paths for improved SPI on ITER. At the boundary, measured neutral density and ionization source fluxes are strongly poloidally asymmetric, implying a 2D treatment is needed to model pedestal fuelling. Detailed measurements of pedestal and SOL quantities and impurity charge state radiation in detached divertors has validated edge fluid modelling and new self-consistent 'pedestal-to-divertor' integrated modeling that can be used to optimize reactors. New feedback adaptive ELM control minimizes confinement reduction, and RMP ELM suppression with sustained high core performance was obtained for the first time with the outer strike point in a W-coated, compact and unpumped small-angle slot divertor. Advances have been made in integrated operational scenarios for ITER and power plants. Wide pedestal intrinsically ELM-free QH-modes are produced with more reactor-relevant conditions, Low torque IBS with W-equivalent radiators can exhibit predator-prey oscillations in Te and radiation which need control. High-β P scenarios with q min > 2, q 95 –7.9, β N > 4, β T –3.3% and H 98y2 > 1.5 are sustained with high density ($\overline{n}$ = 7E19 m -3 , f G –1) for 6 τ E , improving confidence in steady-state tokamak reactors. Diverted NT plasmas achieve high core performance with a non-ELMing edge, offering a possible highly attractive core-edge integration solution for reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2003

Kinetic Monte Carlo Simulations of Rod Eutectics and the Surface Roughening Transition in Binary Alloys

In this paper we report on two related topics: Kinetic Monte Carlo simulations of the steady state growth of rod eutectics from the melt, and a study of the surface roughness of binary alloys. We have implemented a three dimensional kinetic Monte Carlo (kMC) simulation with diffusion by pair exchange only in the liquid phase. Entropies of fusion are first chosen to fit the surface roughness of the pure materials, and the bond energies are derived from the equilibrium phase diagram, by treating the solid and liquid as regular and ideal solutions respectively. A simple cubic lattice oriented in the {100} direction is used. Growth of the rods is initiated from columns of pure B material embedded in an A matrix, arranged in a close packed array with semi-periodic boundary conditions. The simulation cells typically have dimensions of 50 by 87 by 200 unit cells. Steady state growth is compliant with the Jackson-Hunt model. In the kMC simulations, using the spin-one Ising model, growth of each phase is faceted or nonfaceted phases depending on the entropy of fusion. There have been many studies of the surface roughening transition in single component systems, but none for binary alloy systems. The location of the surface roughening transition for the phases of a eutectic alloy determines whether the eutectic morphology will be regular or irregular. We have conducted a study of surface roughness on the spin-one Ising Model with diffusion using kMC. The surface roughness was found to scale with the melting temperature of the alloy as given by the liquidus line on the equilibrium phase diagram. The density of missing lateral bonds at the surface was used as a measure of surface roughness.

Bentz, Daniel N.↗
DOE OSTI2024

Effects of alpha-ion stopping on ignition and ignition criteria in inertial confinement fusion experiments

With the advent of ignited plasmas at the National Ignition Facility (NIF), alpha physics has become a driving factor in theoretical understanding and experimental behavior. In this communication, we explore aspects of direct alpha-ion heating through comparison of the consequences from the one-fluid and two-fluid models in the hydrodynamic approach. We show that the case with all alpha energy deposited in electrons raises the ignition criteria by ∼4 keV or ∼0.2 g/cm2 in the hotspot relative to the case with all alpha energy deposited in ions. In the case of the recently ignited NIF implosion, 30% of the 3.5 MeV α energy is deposited into the DT fuel ions, for which there is negligible difference between the one-fluid and two-fluid ignition criteria. However, changes in the ion stopping fraction through profile effects and alternate stopping power models could lead to ignition curve shifts of ∼1 keV.

Physics↗
DOE OSTI2020

How data science methods can improve the quality and efficiency of ICF and HEDP research

Data Science methods (many that are Bayesian based) are widely used in the physical sciences to estimate model parameters from experimental data, synthesize heterogeneous data, calibrate models, design experiments, and determine statistical significance of data. These methods provide a wealth of advantages over traditional analysis techniques because: 1) uncertainties are rigorously defined and propagated naturally through complex systems including covariance, 2) prior information is captured within the analysis framework (including rad-MHD and rad-hydro simulations), 3) competing models can be selected and/or ruled out using quantitative criteria, and 4) complex, heterogeneous data can be incorporated simultaneously. While these methods have been widely adopted as the gold standard in fields such as particle physics, astronomy, and biology, they have been slow to catch on in Inertial Confinement Fusion (ICF) and High Energy Density Physics (HEDP) research. Recently, several teams at LLNL, SNL, LANL, and the LLE have been exploring the use of these tools in their research and have found success. Here we propose that a concerted effort to consolidate these independent research efforts by developing and deploying common tools for use across the complex can revolutionize the way we approach data analysis, assimilation of theory and experiment, and decision making. The Bayesian formalism provides a means to accomplish this, but we are lacking certain infrastructure to make it happen on a large scale. Furthermore, once adopted, these techniques can be used to develop standards by which discoveries can be judged, similar to the so-called 5σ rule in high energy particle physics. Such standards may be used in the future to address the issue of unknown reproducibility in ICF and HED experiments caused by low shot rate and high cost per experiment. Our goals as a group are to advance the state of the art in HED measurement science by enabling: 1) better inferences from data with well-defined uncertainties, 2) better use of the data we have and continue to collect, 3) intelligent synthesis of data, 4) evaluation of the statistical significance of our data, and 5) informed decision making regarding the design of new experiments and instruments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Demonstration of low-density, high-performance operation of sustained spheromaks and favorable scalability toward compact, low-cost fusion power plants (Final Scientific/Technical Report)

This project worked to advance the technical viability of a novel method for efficiently sustaining stable, high-performance spheromak plasma configurations to serve as the basis of compact, low-cost fusion power plants. In particular, our group worked to improve the method of Steady Inductive Helicity Injection (SIHI) with Imposed-Dynamo Current Drive (IDCD) for spheromak plasma sustainment. Prior to this project demonstrations of this plasma sustainment technology have achieved plasma performance consistent with entry milestone 3 of the BETHE FOA. Research and development (R&D) activities for this project were focused on increasing plasma performance toward a level consistent with exit milestone 4. To do this the PI and his group worked to increase the performance of sustained spheromaks produced in an existing experimental prototype (HIT-SIU) while improving confidence in projections to and design of future, higher performance devices through three primary R&D activities: 1) Improved control over the density of plasma in the device throughout a discharge to provide a pathway for demonstration of spheromaks Ohmically heating to the Mercier beta limit via: a. Fueling the device directly with plasma through the installation of pre-ionized source on the injectors b. Optimization of electrical current waveforms in the driver circuits to enable low-density plasma formation with a lower fueling rate 2) Computational demonstration of a validated, realistic injector circuit coupled to a dynamic plasma model capable of use as a design tool for SIHI drivers and associated circuits for new experimental design points on the pathway to commercial reactors. The improvements in plasma performance achieved during research activity 1), and the computational projections performed in research activity 2) increased the technological readiness level (TRL) of this fusion energy concept toward a level sufficient to attract early-stage private investment and/or other forms of follow-on investment to pursue required R&D activities required for the eventual fusion power plants based on this novel technical approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A model-independent data assimilation (MIDA) module and its applications in ecology

Abstract. Models are an important tool to predict Earth system dynamics. An accurate prediction of future states of ecosystems depends on not only model structures but also parameterizations. Model parameters can be constrained by data assimilation. However, applications of data assimilation to ecology are restricted by highly technical requirements such as model-dependent coding. To alleviate this technical burden, we developed a model-independent data assimilation (MIDA) module. MIDA works in three steps including data preparation, execution of data assimilation, and visualization. The first step prepares prior ranges of parameter values, a defined number of iterations, and directory paths to access files of observations and models. The execution step calibrates parameter values to best fit the observations and estimates the parameter posterior distributions. The final step automatically visualizes the calibration performance and posterior distributions. MIDA is model independent, and modelers can use MIDA for an accurate and efficient data assimilation in a simple and interactive way without modification of their original models. We applied MIDA to four types of ecological models: the data assimilation linked ecosystem carbon (DALEC) model, a surrogate-based energy exascale earth system model: the land component (ELM), nine phenological models and a stand-alone biome ecological strategy simulator (BiomeE). The applications indicate that MIDA can effectively solve data assimilation problems for different ecological models. Additionally, the easy implementation and model-independent feature of MIDA breaks the technical barrier of applications of data–model fusion in ecology. MIDA facilitates the assimilation of various observations into models for uncertainty reduction in ecological modeling and forecasting.

58 GEOSCIENCES↗
NASA NTRS2021

A model-independent data assimilation (MIDA) module and its applications in ecology

Models are an important tool to predict Earth system dynamics. An accurate prediction of future states of ecosystems depends on not only model structures but also parameterizations. Model parameters can be constrained by data assimilation. However, applications of data assimilation to ecology are restricted by highly technical requirements such as model-dependent coding. To alleviate this technical burden, we developed a model-independent data assimilation (MIDA) module. MIDA works in three steps including data preparation, execution of data assimilation, and visualization. The first step prepares prior ranges of parameter values, a defined number of iterations, and directory paths to access files of observations and models. The execution step calibrates parameter values to best fit the observations and estimates the parameter posterior distributions. The final step automatically visualizes the calibration performance and posterior distributions. MIDA is model independent, and modelers can use MIDA for an accurate and efficient data assimilation in a simple and interactive way without modification of their original models. We applied MIDA to four types of ecological models: the data assimilation linked ecosystem carbon (DALEC) model, a surrogate-based energy exascale earth system model: the land component (ELM), nine phenological models and a stand-alone biome ecological strategy simulator (BiomeE). The applications indicate that MIDA can effectively solve data assimilation problems for different ecological models. Additionally, the easy implementation and model-independent feature of MIDA breaks the technical barrier of applications of data–model fusion in ecology. MIDA facilitates the assimilation of various observations into models for uncertainty reduction in ecological modeling and forecasting.

Earth system dynamics↗
DOE OSTI2024

Modeling the Low-Pressure High-Voltage Branch of the Paschen Curve for Hydrogen and Deuterium

A physical and numerical model of the Townsend discharge in molecular hydrogen and deuterium has been developed to meet the needs of designing a plasma-based switching device for power grid applications. Here, the model allows to predict the low-pressure branch of the Paschen curve for applied voltage in the range of several hundred kilovolts. In the regime of interest, electrons are in a runaway state and ionization by ions and fast neutrals sustains the discharge. It was essential to correctly account for both gas-phase and surface interactions (electron emission and electron backscattering), especially in terms of their dependence on particle energy. The model yields results consistent with prior data obtained for lower voltage. The three-species (electrons, ions, and fast neutrals) model successfully captures the essential physics of the process.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Scalable Implicit Solvers with Dynamic Mesh Adaptation for a Relativistic Drift-Kinetic Fokker–Planck–Boltzmann Model

In this work we consider a relativistic drift-kinetic model for runaway electrons along with a Fokker–Planck operator for small-angle Coulomb collisions, a radiation damping operator, and a secondary knock-on (Boltzmann) collision source. Here, we develop a new scalable fully implicit solver utilizing finite volume and conservative finite difference schemes and dynamic mesh adaptivity. A new data management framework in the PETSc library based on the p4est library is developed to enable simulations with dynamic adaptive mesh refinement (AMR), distributed memory parallelization, and dynamic load balancing of computational work. This framework and the runaway electron solver building on the framework are able to dynamically capture both bulk Maxwellian at the low-energy region and a runaway tail at the high-energy region. To effectively capture features via the AMR algorithm, a new AMR indicator prediction strategy is proposed that is performed alongside the implicit time evolution of the solution. This strategy is complemented by the introduction of computationally cheap feature-based AMR indicators that are analyzed theoretically. Numerical results quantify the advantages of the prediction strategy in better capturing features compared with nonpredictive strategies; and we demonstrate trade-offs regarding computational costs. The robustness with respect to model parameters, algorithmic scalability, and parallel scalability are demonstrated through several benchmark problems including manufactured solutions and solutions of different physics models. We focus on demonstrating the advantages of using implicit time stepping and AMR for runaway electron simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Mixing in ICF implosions on the National Ignition Facility caused by the fill-tube

The micrometer-scale tube that fills capsules with thermonuclear fuel in inertial confinement fusion experiments at the National Ignition Facility is also one of the implosion's main degradation sources. It seeds a perturbation that injects the ablator material into the center, radiating away some of the hot-spot energy. This paper discusses how the perturbation arises in experiments using high-density carbon ablators and how the ablator mix interacts once it enters the hot-spot. Both modeling and experiments show an in-flight areal-density perturbation and localized x-ray emission at stagnation from the fill-tube. Simulations suggest that the fill-tube is degrading an otherwise 1D implosion by ~2×, but when other degradation sources are present, the yield reduction is closer to 20%. Characteristics of the fill-tube assembly, such as the through-hole size and the glue mass, alter the dynamics and magnitude of the degradation. Overall, these aspects point the way toward improvements in the design, some of which (smaller diameter fill-tube) have already shown improvements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Fusion of Experiments and Simulations for Real-Time Identification of Pipeline Defects

In this study, we explored fusion of experiments and simulations for real time identification of pipeline defects across physical and non-physical domains. The challenges associated to data processing were addressed and a combined classification models was presented via CNN models. In addition, regression model based on XGBOOST is built to determine the defect location and defect dimension from data-driven features of guided wave signals captured by SMS fiber optic sensor.

deep learning↗
DOE OSTI2025

Fusion of Experiments and Simulations for Real-Time Identification of Pipeline Defects

In this study, we explored fusion of experiments and simulations for real time identification of pipeline defects across physical and non-physical domains. The challenges associated to data processing were addressed and a combined classification models was presented via CNN models. In addition, regression model based on XGBOOST is built to determine the defect location and defect dimension from data-driven features of guided wave signals captured by SMS fiber optic sensor.

deep learning↗
DOE OSTI2025

Physics-informed transformation toward improving the machine-learned NLTE models of ICF simulations

The integration of machine-learning techniques into inertial confinement fusion (ICF) simulations has emerged as a powerful approach for enhancing computational efficiency. By replacing the costly nonlocal thermodynamic equilibrium (NLTE) model with machine-learning models, significant reductions in calculation time have been achieved. However, determining how to optimize machine-learning-based NLTE models in order to match ICF simulation dynamics remains challenging, underscoring the need for physically relevant error metrics and strategies to enhance model accuracy with respect to these metrics. Thus, we propose novel physics-informed transformations designed to emphasize energy transport, use these transformations to establish new error metrics, and demonstrate that they yield smaller errors within reduced principal-component spaces compared to conventional transformations. Published by the American Physical Society 2025

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗