Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “plasma models”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Transport and confinement physics: Chapter 2 of the special issue: on the path to tokamak burning plasma operation

Progress in physics understanding and theoretical model development of plasma transport and confinement (TC) in the ITPA TC Topical Group since the publication of the ITER Physics Basis (IPB) document (Doyle et al 2007 Nucl. Fusion 47 S18) was summarized focusing on the contributions to ITER and burning plasma prediction and control. This paper provides a general and streamlined overview on the advances that were mainly led by the ITPA TC joint experiments and joint activities for the last 15 years (see JEX/JA table in appendix). This paper starts with the scientific strategy and scope of the ITPA TC Topical group and overall picture of the major progress, followed by the progress of each research field: particle transport, impurity transport, ion and electron thermal turbulent transport, momentum transport, impact of 3D magnetic fields on transport, confinement mode transitions, global confinement, and reduced transport modeling. Cross references with other Topical Groups are given in order to highlight overlapped topics, such as the 3D effect on the plasma transport in the edge and L-H transition physics. The increasing overlap between the topical groups is a reflection of the progress on integrating the known physics into comprehensive models that are better and better able to reproduce the plasma transport. In recent years, such integration has become increasingly prevalent when considering transport from the SOL, through the edge pedestal, and into the plasma core. In the near future, increased collaboration also with the magneto-hydrodynamic and energetic particles community will be important as we approach burning plasma conditions in next-step fusion devices. A summary of remaining challenges and next steps for each research field is given in the Summary section.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Breaking Beams: A Kinetic Treatment of the Filamentation and Forward Scatter of Laser Beams in Plasma (Abbreviated Final Report)

Precise and predictable propagation and deposition of laser energy in high energy density (HED) plasma is threatened by the breakup of the beam into filaments. Such filaments can be deflected and sprayed but also locally intensify the light, increasing the likelihood of other problematic interactions with the plasma. Adequately modeling, anticipating, and avoiding such processes in HED experiments requires a thorough understanding of the thresholds and nonlinear saturation mechanisms, particularly for the weakly collisional, high temperature plasmas achieved at the National Ignition Facility (NIF).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Breaking Beams: A Kinetic Treatment of the Filamentation and Forward Scatter of Laser Beams in Plasma

Precise and predictable propagation and deposition of laser energy in high energy density (HED) plasma is threatened by the breakup of the beam into filaments. Such filaments can be deflected and sprayed but also locally intensify the light, increasing the likelihood of other problematic interactions with the plasma. Adequately modeling, anticipating, and avoiding such processes in HED experiments requires a thorough understanding of the thresholds and nonlinear saturation mechanisms, particularly for the weakly collisional, high temperature plasmas achieved at the National Ignition Facility (NIF).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Integrated modeling of boron powder injection for real-time plasma-facing component conditioning

An integrated modeling framework for investigating the application of solid boron (B) powder injection for real-time surface conditioning of plasma-facing components (PFCs) in tokamak environments is presented. Utilizing the DIII-D impurity powder dropper (IPD) setup, this study simulates B powder injection scenarios ranging from milligrams to tens of milligrams per second, corresponding to boron flux rates of 10 20 –10 21 B/s in standard L-mode conditions. The comprehensive modeling approach combines EMC3-EIRENE for simulating the deuterium plasma background and the Dust Injection Simulator (DIS) for the ablation and transport of the boron powder particles. EMC3 trace impurity fluid modeling results show substantial boron transport to the inboard lower divertor, predominantly influenced by the main ion plasma flow. The dependency on powder particle size (5-250 µm) was found to be insignificant for the scenario considered. The effects of erosion and redeposition were considered to reconcile the discrepancies with experimental observations, which saw substantial deposition on the outer divertor plasma-facing components. For this purpose, the WallDYN3D code was updated to include boron sources within the plasma domain and integrated into the modeling framework. The mixed-material migration modeling shows evolving boron deposition patterns, suggesting the formation of mixed B-C layers or predominantly B coverage depending on the powder mass flow rate. While the modeling outcomes at lower B injection rates tend to align with DIII-D experimental observations, the prediction of near-pure boron layers at higher rates has yet to be experimentally verified in the carbon environment of the DIII-D tokamak. The extensive reach of boron layers found in the modeling suggests the need for modeling that encompasses the entire wall geometry for more accurate experimental correlations. This integrated approach sets a precedent for analyzing and applying real-time in-situ boron coating techniques in advanced tokamak scenarios, potentially extendable to ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Recent progress in L-H transition studies at JET: Tritium, Helium, Hydrogen and Deuterium

In this report we present an overview of results from a series of L-H transition experiments undertaken at JET since the installation of the ITER-like-wall (JET-ILW), with Beryllium wall tiles and a Tungsten divertor. Tritium, Helium and Deuterium plasmas have been investigated. Initial results in Tritium show ohmic L-H transitions at low density and the power threshold for the L-H transition (P LH ) is lower in Tritium plasmas than in Deuterium ones at low densities, while we still lack contrasted data to provide a scaling at high densities. In Helium plasmas there is a notable shift of the density at which the power threshold is minimum (n e,min ) to higher values relative to Deuterium and Hydrogen references. Above n e,min (He) the L-H power threshold at high densities is similar for D and He plasmas. Transport modelling in slab geometry shows that in Helium neoclassical transport competes with interchange-driven transport, unlike in Hydrogen isotopes. Measurements of the radial electric field in Deuterium plasmas show that E r shear is not a good indicator of proximity to the L-H transition. Transport analysis of ion heat flux in Deuterium plasmas show a non-linearity as density is decreased below n e,min . Lastly, a regression of the JET-ILW Deuterium data is compared to the 2008 ITPA scaling law.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

On the non-existence of stepped-pressure equilibria far from symmetry

The Stepped Pressure Equilibrium Code (SPEC) (Hudson et al 2012 Phys. Plasmas 19 112502) has been successful in the construction of equilibria in 3D configurations that contain a mixture of flux surfaces, islands and chaotic magnetic field lines. In this model, the plasma is sliced into sub-volumes separated by ideal interfaces, and in each volume the magnetic field is a Beltrami field. In the cases where the system is far from possessing a continuous symmetry, such as in stellarators, the existence of solutions to a stepped-pressure equilibrium with given constraints, such as a multi-region relaxed MHD minimum energy state, is not guaranteed but is often taken for granted. Using SPEC, we have studied two different scenarios in which a solution fails to exist in a slab with analytic boundary perturbations. Here we found that with a large boundary perturbation, a certain interface becomes fractal, corresponding to the break up of a Kolmogorov–Arnold–Moser (KAM) surface. Moreover, an interface can only support a maximum pressure jump while a solution of the magnetic field consistent with the force balance condition can be found. An interface closer to break-up can support a smaller pressure jump. We discovered that the pressure jump can push the interface closer to being non-smooth through force balance, thus significantly decreasing the maximum pressure it can support. Our work shows that a convergence study must be performed on a SPEC equilibrium with interfaces close to break-up. These results may also provide insights into the choice of interfaces and have applications in finding out the maximum pressure a machine can support.

3D equilibrium↗
DOE OSTI2024

Reconstruction of plasma equilibrium and eddy current induced error fields in the Lithium Tokamak eXperiment-β (Final Report)

This is the final report for the DOE award DE-SC0019239, titled “Reconstruction of plasma equilibrium and eddy current induced error fields in the Lithium Tokamak eXperiment-b” for the period July 1, 2018 – July 2, 2023. This award supported research into computational models for plasma equilibrium and 3D conducting structures in the Lithium Tokamak eXperiment-b (LTX-b) using the PSI-Tri/TokaMaker and PSI-Tet codes. Additional detail on progress in support of these tasks is provided below.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A review on charged-particle transport modeling for laser direct-drive fusion

Inertial confinement fusion (ICF) with the laser-indirect-drive scheme has recently made a tremendous breakthrough recently after decades of intensive research effort. Taking this success to the next step, the ICF community is coming to a general consensus that laser direct-drive (LDD) fusion might be the viable way for enabling inertial fusion energy (IFE) and high-gain targets for other applications. Designing and understanding LDD fusion targets heavily rely on radiation-hydrodynamic code simulations, in which charged-particle transport plays an essential role in modeling laser-target energy coupling and bootstrap heating of fusion-produced α-particles. To better simulate charged-particle transport in LDD targets, over the past four decades the plasma physics community has advanced transport calculations from simple plasma physics models to sophisticated computations based on first-principles methods. In this review, we give an overview of the current status of charged-particle transport modeling for LDD fusion, including what challenges we still face and the possible paths moving forward to advance transport modeling for ICF simulations. We hope this review will provide a summary of exciting challenges to stimulate young minds to enter the field, facilitate further progress in understanding warm-dense matter physics, and ultimately bridge toward the success of reliable LDD fusion designs for IFE and other high-gain ICF applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Thermal conductivity of a laser plasma

Here, we present a model of the electron thermal conductivity of a laser-produced plasma. The model, supported by Vlasov–Fokker–Planck simulations, predicts that laser absorption reduces conductivity by forcing electrons out of a Maxwell–Boltzmann equilibrium, which results in depletion of both low-velocity bulk electrons and high-velocity tail electrons. We show that both the bulk and tail electrons approximately follow super-Gaussian distributions, but with distinct exponents that each depend on the laser intensity and wavelength through the parameter a = Zv$^{2}_{E}$/v$^{2}_{T}$. For a value of a = 0.5, tail depletion reduces the thermal conductivity to half its zero-intensity value. We present our results as simple analytic fits that can be readily implemented in any radiation-hydrodynamics code or used to correct the local limit of non-local conduction models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Electron-impact ionization from excited states method

We present a distorted-wave with exchange (DWE) formulation for computing electron-impact ionization (EII) cross sections within the recently developed excited states method (ESM) for plasmas. The ESM provides a self-consistent quantum-mechanical description of both bound and continuum electronic states in dense plasmas, incorporating finite-temperature and screening effects absent in isolated-atom approaches. Using this framework, we calculate EII cross sections for lithium across a wide range of temperatures and densities. We show that the ESM reproduces the isolated-atom limit at low densities while capturing strong plasma effects, such as pressure ionization and shape resonances, at solid density. Comparisons with isolated-atom DWE and average-atom (AA) calculations reveal that, although AA cross sections can approximate configuration average behavior, the ESM provides more accurate, state-resolved answers. In conclusion, these results demonstrate that the ESM allows one to obtain rate coefficients for collisional-radiative modeling of plasmas out of local thermodynamical equilibrium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A reduced-model (nSOLT) simulation of neutral recycling effects on plasma turbulence in the divertor region of MAST-U

The 2D scrape-off-layer turbulence code (nSOLT), which includes 1D kinetic neutral–plasma interactions, is applied to study effects of neutral recycling on plasma turbulence for parameters illustrative of the MAST-U divertor region. Neutral recycling is modeled by injecting a fraction of the parallel plasma flux to the divertor back into the simulation domain as a source of Franck–Condon-distributed neutrals. Stationary sources, concentrated at the magnetic separatrix, model plasma streaming into the divertor region from the upstream scrape-off-layer and sustain plasma turbulence absent neutral recycling. Starting from one such no-neutrals equilibrium, we initiate recycling in a numerical experiment designed to diagnose and identify the effects of various neutral–plasma interactions on the divertor plasma, divertor turbulence, and plasma exhaust. The onset of recycling triggers an initial burst of enhanced cross field plasma transport that is quelled by ionization cooling and charge–exchange (CX) friction, with growing neutral pressure, leading to a quiescent, turbulence-free state. Diagnosis of this transient burst reveals that (1) the sudden increase in plasma density due to ionization dominates the onset of the burst, (2) electron cooling due to ionization increases collisionality and disconnects blob filaments from the sheath, and (3) CX friction drives tripole polarization of a blob that can dominate the curvature-driven dipole polarization, leading to the stagnation of blob propagation and reduced radial turbulent transport. It is shown that CX friction is negligible compared to sheath physics in determining equilibrium mean flow shearing rates, for parameters considered herein (specifically a short connection length to the divertor target), while it can significantly reduce interchange-instability growth rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Study of Alfvén eigenmode stability in Quasi-Poloidal Stellarator (QPS) plasma using a Landau closure model

The aim of this study is to analyze the linear stability of Alfvén eigenmodes (AE) in the QPS device heated by a tangential neutral beam injector (NBI). The analysis is performed using the gyro-fluid code FAR3d, that solves the reduced MHD equations for the thermal plasma coupled with moments of the kinetic equation for the energetic particles (EP). The AE stability is calculated in several operational regimes of the tangential NBI: EP β between 0.001 and 0.1, EP energy between 12 and 180 keV and different radial locations of the beam. The analysis is performed for vacuum and finite β equilibria as well as QPS configurations with two and three periods. The EP β threshold in the vacuum case is 0.001 and the AE frequency is lower as the energy of the EP population decreases. Toroidal Alfvén eigenmodes with f = 80–120 kHz and elliptical AE between f = 120–350 kHz are triggered between the middle-outer plasma region ($r/a \gt 0.5$). The AE stability improves in the simulations with finite β equilibria and three period configurations with respect to the vacuum case with two periods because the continuum gaps are slender, leading to a higher threshold of the EP β, above 0.03 for the AEs triggered by the helical mode families. Helical effects are not strong enough to destabilize Helical Alfvén eigenmodes, the AEs with the largest growth rates are triggered by the n = 1 and n = 2 toroidal families.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Predicting operational windows of ELMs suppression by resonant magnetic perturbations in the DIII-D and KSTAR tokamaks

A newly developed plasma response model, combining the nonlinear two-fluid MHD code TM1 and toroidal MHD code GPEC run in ideal mode, quantitatively predicts the narrow isolated q 95 windows (Δq 95 ~ 0.1) of ELM suppression by n = 1, 2 and 3 resonant magnetic perturbations (RMPs) in both DIII-D and KSTAR tokamaks across a wide range of plasma parameters. The key physics that unites both experimental observations and our simulations is the close alignment of essential resonant q-surfaces and the location of the top of the pedestal prior to an ELM. This alignment permits an applied RMP to produce field penetration due to the lower E × B rotation at the pedestal top rather than being screened. The model successfully predicts that narrow magnetic islands form when resonant field penetration occurs at the top of pedestal, and these islands are easily screened when q 95 moves off resonance, leading to very narrow windows of ELM suppression (typically Δq 95 ~ 0.1). Furthermore, the observed reduction in the pedestal height is also well captured by the calculated classical collisional transport across the island. We recover observed q 95 , β N and plasma shape dependence of ELM suppression due to the effect of magnetic islands on pedestal transport and Peeling-Ballooning- Mode (PBM) stability. Importantly, experiments do occasionally observe wide windows of ELM suppression (Δq 95 > 0.5). Our model reveals that at low pedestal-top density multiple islands open, leading to wide operational windows of ELM suppression consistent with experiment. The model indicates that wide q95 windows of ELM suppression can be achieved at substantially higher pedestal pressure with less confinement degradation in DIII-D by operating at higher toroidal mode number (n = 4) RMPs. As a result, this can have significant implications for the operation of the ITER ELM control coils for maintaining high confinement together with ELM suppression.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Integrated modelling: Coupling of surface evolution and plasma-impurity transport

During the interaction of the scrape off layer (SOL) plasma with the first wall the evolution of both wall and plasma are tightly coupled: The erosion of the first wall leads to an impurity concentration in the plasma which affects the particle and power balance in the plasma. In turn the impurities, when leaving the plasma via transport, can form deposits and mixed materials, far away from their initial source location. These deposits can be eroded, allowing the impurities to stepwise migrate through the fusion device until they end up at a location where the plasma at the wall is cold enough and no further erosion occurs. To describe these processes an integrated model of surface evolution and plasma transport of impurities is needed. The WallDYN code achieves this required coupling of processes by parameterising the output of surface evolution- and plasma impurity-migration-codes by analytical models. For a given fixed background plasma it evolves the surface composition, derives impurity flux into and from the plasma and can from this derive the impurity densities in the plasma. This paper will show the importance of including this recycling of impurities at the wall in impurity migration modelling: The 13 CH 4 seeding experiment performed in the Wendelstein 7-X Stellerator is modelled using the recent extension of the WallDYN code to 3D plasma and wall geometries. A comparison with post mortem analysis of the 13 C deposition shows both qualitative and quantitative agreement with the WallDYN calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Plasma Emission, Wave Speeds, and Temperatures of Laser-Supported Plasmas in Air

A study of the generation and subsequent expansion of plasmas using a single beam (532 nm) from the second harmonic generator of a Q-switched Nd:YAG laser is described. The plasmas in air are studied using an intensified charge-coupled device camera over periods of 1000 ns following the laser pulses. The plasma experimental results are interpreted in terms of wave speeds and plasma temperatures. Energy balance equations are considered to estimate the plasma temperatures and the laser-supported radiation (LSR) and the laser-supported detonation (LSD) wave speeds. Within 1000 ns observation time, the circular spectral images turn elliptical because of plasma expansion toward the laser beam. The plasma expansion increases with increasing laser irradiance. A laser irradiance range of [Formula: see text] yields LSR wave speeds between 45 and 170 km/s and LSD wave speeds of 90–130 km/s. The measured wave speeds are compared with the LSR and the LSR wave speed models. The plasma temperatures are between 90,000 and 120,000 K at 15 ns after the laser pulse based on measured emission spectra and are between 640,000 and 760,000 K at 3 ns based on the LSR model.

Engineering↗
DOE OSTI2023

Tech-X contributions to the Center for Tokamak Transients (CTTS) Scientific Discovery through Advanced Computing (SciDAC) project (Final Report)

The Center for Tokamak Transients Simulations (CTTS) consists of computational plasma physicists, computer scientists, and applied mathematicians who are together developing and applying simulation software to improve our understanding of and ability to predict and control transient events in tokamak discharges, particularly those that lead to violent disruptions. We have major simulation activities concerning disruptions that are caused by violating ideal MHD stability criteria, vertical displacement events, resistive wall modes, and tearing modes. We are performing comprehensive modeling of disruption mitigation techniques, particularly by shattered pellet injection (SPI). The Center is built around two mature extended MHD codes: NIMROD and M3D-C1. These codes are being extended in a number of ways as required by the demands of disruption modeling, and are being optimized to work well for larger problem size and modern computer architectures as required to model tokamak plasmas with realistic parameters over the required timescales. The codes will be coupled to the FronTier code for state-of-the-art modeling of SPI. Validation with experimental data on DIII-D, NSTX, and JET is being facilitated by rigorous uncertainty quantification (UQ). The high-level objectives of the center are to: (1) Develop better understanding and improved predictive capability of when crossing a linear stability boundary will likely lead to a disruption, (2) Provide a tool for calculating vessel forces, both axisymmetric and non-axisymmetric, for worst case VDE and other disruptions in tokamaks, (3) Develop increased understanding and improved predictive capability for locked mode disruptions and how best to avoid them, (4) Develop, verify, and validate 3D models for disruption mitigation by shattered pellets that can be used to design and optimize a system for future tokamaks, and (5) Increase the efficiency and scalability of the M3D-C1 and NIMROD codes so that they can model disruptions in ITER scale devices using the most powerful DOE HPC hardware, and (6) Demonstrate code-coupling via FronTier coupling with NIMROD and M3D-C1 and by interfacing these codes within a WDM framework.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A Summary of the UCLA HANE-Laser Experiment: 2011-2020

In the early years of this century, there was renewed interest at DTRA in artificial radiation belts, the dynamics of high altitude nuclear explosions that produced them, the development of large-scale kinetic plasma computer models at LLNL and LANL (particularly in the form of “hybrid” {i.e., particle ions, massless fluid electrons} algorithms), and the building of a laser facility at UCLA, under the direction of Prof. Niemann, which was connected to the Large Plasma Device (LAPD), a DOE user facility. These advances sparked the idea for a new laser experimental program to examine early-time HANE issues, reviving the concept from a former program at NRL in the 90’s. The basic motivation for this new effort can be traced back to the first DTRA artificial radiation belt workshop at Stanford in 2009. Subsequent discussions then led to a formal proposal from UCLA that was submitted to DTRA (Grant Jones) in 2010, vigorously reviewed, and finally approved in 2011, with funding begun in 2012. An historical perspective of this development process was presented at a DTRA review last year. In this document we review and summarize the major achievements of the DTRA-sponsored UCLA HANE-laser experiment over the past eight years. Our purpose is to briefly describe the major achievements of this program, as documented in the included extensive list of journal publications [which does not include all the publications nor any of the many invited and contributed presentations of the UCLA group], the role of the national laboratories in this effort, and how this work has impacted (and will continue to improve) our understanding of high altitude nuclear events.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗