Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Plasma confinement”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

The separatrix electron density in JET, ASDEX upgrade and alcator C-Mod H-mode plasmas: A common evaluation procedure and correlation with engineering parameters

The separatrix electron density is an important parameter for core-edge scenario integration in tokamak devices, as it influences plasma confinement, divertor detachment and disruption avoidance. This quantity has been measured in H-mode discharges on JET, ASDEX Upgrade and Alcator C-Mod by applying the same fitting function to Thomson scattering measurements, and by employing the same analysis technique based on scrape-off layer power balance. To estimate the power crossing the separatrix, the inter-ELM time derivative of the plasma energy d W /d t has been experimentally evaluated and found to be approximately a constant fraction of the absorbed heating power. Correlations between n e,sep and engineering parameters have been investigated, revealing that n e,sep scales with the divertor neutral pressure p 0,div in a similar manner across all devices. Additionally, when n e,sep is normalized to the obtained n e,sep dependency, no clear correlation with the plasma current is found. These observations are in agreement with the 2-point model, which suggests that the upstream separatrix density is mainly set by the recycling at the divertor target.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A framework for discrete optimization of stellarator coils

Designing magnets for three-dimensional plasma confinement is a key task for advancing the stellarator as a fusion reactor concept. Stellarator magnets must produce an accurate field while leaving adequate room for other components and being reasonably simple to construct and assemble. In this paper, a framework for coil design and optimization is introduced that enables the attainment of sparse magnet solutions with arbitrary restrictions on where coils may be located. The solution space is formulated as a 'wireframe' consisting of a mesh of interconnected wire segments enclosing the plasma. Two methods are developed for optimizing the current distribution on a wireframe: Regularized Constrained Least Squares, which uses a linear least-squares approach to optimize the currents in each segment, and Greedy Stellarator Coil Optimization, a fully discrete procedure in which loops of current are added to the mesh one by one to achieve the desired magnetic field on the plasma boundary. Examples are presented of solutions obtainable with each method, some of which achieve high field accuracy while obeying spatial constraints that permit easy assembly.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2023

Self-Similar Outflows at the Source of the Fast Solar Wind: A Smoking Gun of Multiscale Impulsive Reconnection?

We present results of a quantitative analysis of structured plasma outflows above a polar coronal hole observed by the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO) spacecraft. In a 6 hr interval of continuous high-cadence SDO/AIA images, we identified more than 2300 episodes of small-scale plasma flows in the polar corona. The mean upward flow speed measured by the surfing transform technique is estimated to be 122 ± 34 km s−1, which is comparable to the local sound speed. The typical recurrence period of the flow episodes is 10–30 minutes, and the mean duration and transverse size of each episode are about 3–5 minutes and 3–4 Mm, respectively. The largest identifiable episodes last for tens of minutes and reach widths up to 40 Mm. For the first time, we demonstrate that the polar coronal-hole outflows obey a family of power-law probability distributions characteristic of impulsive interchange magnetic reconnection. Turbulent photospheric driving may play a crucial role in releasing magnetically confined plasma onto open field. The estimated occurrence rate of the detected self-similar coronal outflows is sufficient for them to make a dominant contribution to the fast-wind mass and energy fluxes and to account for the wind's small-scale structure.

jets↗
DOE OSTI2024

Survey of tungsten gross erosion from main plasma facing components in WEST during a L-mode high fluence campaign

An initial high fluence campaign was performed in WEST, in 2023, on the newly installed actively cooled tungsten divertor composed of ITER-grade monoblocks. The campaign consisted in the repetition of a 60 s long Deuterium L-mode pulse in attached divertor conditions, cumulating over 10000s of plasma exposure. A maximum deuterium fluence of approximately 5⋅1⁢026 m−2 was reached in the outer strike point region, representative of a few high performance ITER pulses. Gross tungsten erosion inferred from visible spectroscopy shows that the most eroded plasma facing component is the inner divertor target with rates ten times larger than on the outer divertor target. The outer midplane tungsten bumpers, located a few centimeters from the plasma, show gross erosion rates two times lower than at the outer divertor. We conclude that the outer midplane bumpers have a negligible contribution to the long range tungsten migration and deposition onto the lower divertor. The cumulated gross erosion rate on the inner divertor translates in an effective gross erosion thickness of about 20μ⁢m, while it is about 2μ⁢m for the outer divertor. Strikingly, these orderings coincide with the thickness of deposits found locally on the divertor: the exposed surfaces of high field side monoblocks are covered with several tens of μ⁢m tungsten deposits, while on the lower field side, few μ⁢m thin tungsten deposits are only found on the magnetically shadowed parts of monoblocks. The strong impact of those deposits on WEST operation, namely perturbation of surface temperature measurement with infra-red thermography, and the emission of flakes causing radiative perturbation of the confined plasma, calls for anticipating similar issues in ITER. In particular, the start of research operation shall consider the definition of a divertor erosion budget in order to anticipate the formation of deleterious deposits.

Fedorczak, N.↗
DOE OSTI2024

Characterization and controllability of radiated power via extrinsic impurity seeding in strongly negative triangularity plasmas in DIII-D

Experiments with extrinsic impurity seeding in strongly negative triangularity shapes in DIII-D achieved radiated power fractions (relative to input power) of up to ≈85% total radiation and ≈55% core radiation in steady-operating conditions. The relationship between core and total radiation was sensitive to impurity species and input power. Attempts to reach higher radiation levels via higher impurity flows resulted in radiative collapse disruptions. Nitrogen, neon, argon, and krypton were tested. Injection was by gas puffing, usually controlled by feeding back real-time estimates for core or total radiating fractions ($P_\textrm{rad} / P_\textrm{input}$). Argon and krypton were controlled by feeding back total radiating fraction, whereas neon was controlled by feeding back the core radiating fraction, due to the lower efficiency of neon as a divertor radiator. Nitrogen flows were pre-programmed. Reasonable total $P_\textrm{rad}$ control target following was achieved with argon or krypton. Control was more challenging with the neon/core radiation configuration, which was more prone to slow response and overshooting of the control target. Poor particle removal contributed to the control challenge: neon particle inventory within the last closed flux surface was roughly constant for up to 1 s (the longest duration tested) after neon injection was halted. However, separate experiments with laser-blow off of non-recycling impurities measured a short impurity confinement time, on the scale of the energy confinement time of ∼100 ms. Modeling with the Aurora impurity transport simulation matched experimental neon density profiles with full recycling (R = 1.0) and weak pumping but predicted rapid decreases in neon inventory if pumping were increased or recycling decreased. This indicates that changes outside the confined plasma (adding a well-placed pump) would improve controllability for all highly recycling species.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Progress in Extending the High-Temperature Plasma Regime in the LHD

In this paper, we show the expansion of the high-temperature regime and related physical phenomena in the Large Helical Device (LHD). At the LHD, we have developed a high-temperature operational regime for both ion (T i ) and electron (T e ) temperatures through understanding the physical characteristics of high-temperature plasma, the enhancement of plasma heating, and the development of technical methods such as wall recycling control and impurity control. In particular, the formation of the ion internal transport barrier (i-ITB) and electron internal transport barrier (e-ITB) in the plasma core region is the key to improving plasma performance. In addition, the deuterium experiment that began in 2017 showed that, due to the isotope effect, the thermal transport of both ions and electrons is suppressed in deuterium plasma compared to light hydrogen plasma, and that a high central temperature can be achieved efficiently. Furthermore, the combination of neutral beam injection (NBI) and electron cyclotron heating (ECH) has extended the operating conditions that simultaneously maintain high T i and T e , and comprehensive research on plasma confinement characteristics with an eye to future fusion reactor conditions has progressed.

Large Helical Device↗
DOE OSTI2026

Simultaneous ELM suppression and divertor detachment via synergistic boron powder and neon injection in EAST

A novel approach for simultaneous power exhaust and edge-localized mode (ELM) control is presented in the Experimental Advanced Superconducting Tokamak discharges, which utilize an ITER-like tungsten divertor. Real-time injection of boron (B) powder and neon (Ne) gas overcomes their limitations encountered when used separately. Pure Ne seeding leads to a narrow operational window constrained by core impurity accumulation and H-mode to L-mode back transitions, while pure solid B injection (SBI) is insufficient for effective divertor cooling. In comparison, their combined use achieves a stable, stationary, ELM-suppressed H-mode with adequate power exhaust. This synergistic scenario features partial energy detachment at the outer divertor while maintaining good plasma confinement (H 98 ∼ 1) with minimal degradation. Two key features of this scenario are: (1) the SBI triggers a persistent Edge Harmonic Mode (EHM), which provides a crucial continuous particle transport channel, preventing Ne and tungsten/molybdenum accumulation without flushing out by ELM, and (2) the B + Ne mixture allows for active optimization of the radiated power profile. Core radiation can be reduced by substituting a portion of the Ne with B, leveraging their complementary non-coronal equilibrium radiation efficiencies. This combined B + Ne injection scheme presents a promising pathway toward integrated core-edge scenarios, offering the potential to minimize total impurity throughput while leveraging an actuator (powder injection) already being considered for ITER.

Boron↗
DOE OSTI2024

Compatibility of divertor detachment and ELM suppression in DIII-D high- β p plasmas with ITER-similar shape

Abstract Integration of transient and steady-state divertor heat fluxes control with a high-performance core is necessary for future fusion reactors. In recent DIII-D high- β p experiments, divertor detachment and simultaneous edge localized mode (ELM) suppression are demonstrated while the plasma confinement quality is maintained high in ITER-similar shape. By optimizing the neon injection in high- β p scenario with ITER-similar shape, deep detachment and ELM suppression are achieved with a high-performance core ( β N ∼ 2.8, β p ∼ 2.3) at q 95 ∼ 7.5. Partial divertor detachment and suppression of large ELMs are achieved at q 95 ∼ 6. The stability analyses suggest that with low neon injection, the density pedestal becomes higher and steeper and the T i profile also increases, therefore the increased edge pressure and higher current density destabilize the Peeling-Ballooning mode (PBM), which would lead to a large ELM collapse. With strong neon gas puffing, the significantly reduced pedestal pressure and current density, due to the degraded T e pedestal, lead to the stabilization of PBM and ELMs are suppressed. For both cases, the coupling between the large radius internal transport barrier (ITB) and edge pedestal is the key reason for maintaining high global performance. The formation of large radius ITB compensates for pedestal degradation. Such results could provide an attractive scenario to well control the transient and steady-state heat flux onto the divertor plates while maintaining good plasma performance, which is an important step toward the steady-state operation of future fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Design of radial interferometer–polarimeter for internal magnetic and density fluctuation measurements at multiple space–time scales in the National Spherical Torus Experiment-Upgrade (NSTX-U)

A Faraday-effect radial interferometer–polarimeter is designed for the National Spherical Torus Experiment-Upgrade (NSTX-U) to measure multiscale magnetic and density fluctuations critical to understanding fusion plasma confinement and stability, including those originating from magnetohydrodynamic instabilities, energetic particle-driven modes, and turbulence. The diagnostic will utilize the three-wave technique with 5 MHz bandwidth to simultaneously measure line-integrated magnetic and density fluctuations up to the ion-cyclotron frequency. Probe beams will be launched radially from the low-field side at the NSTX-U midplane, where the measured Faraday fluctuations mainly correspond to radial magnetic fluctuations that directly link to magnetic transport. A correlation technique will be employed to reduce the measurement noise to below 0.01° enabling detection of small amplitude fluctuations. Two toroidally displaced chords with 7° separation will be installed to measure toroidal mode numbers up to n = 25 for mode identification. Finally, solid-state microwave sources operating at 321 μm (935 GHz) will be used to minimize the impact of the Cotton–Mouton effect.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Project Title: Demonstration High Temperature Superconducting NonPlanar Stellarator Magnet with Advanced Manufactured Assemblies

This is the final report for the project “Demonstration High Temperature Superconducting Non- Planar Stellarator Magnet with Advanced Manufactured Assemblies”, funded by DOE, and performed by Type One Energy from September, 2020, to March 2024 involving the Fusion Technology Institute at the University of Wisconsin–Madison, the Plasma Science and Fusion Center (PSFC) at the Massachusetts Institute of Technology (MIT) and Commonwealth Fusion Systems (CFS) to design and fabricate the first non-planar HTS (REBCO) coil for a high-field stellarator based on the SPARC tokamak’s VIPER cable concept. Stellarators at high fields make high-temperature superconducting magnets necessary for a compact fusion device. But the asymmetric and non-planar nature of its components, especially the magnets make it difficult for scalable producibility. To address these challenges, two promising technologies have emerged: advanced manufacturing (AM) for the supporting plates for forming the magnets, and high-temperature superconducting (HTS) cables inside the plates. AM has advanced enough to produce stellarator components with the necessary geometric complexity, size, and the precision, leading to potentially significant reductions in production time, cost, and waste. The cost of HTS tape has decreased dramatically, and progress in HTS planar magnet development has reached a point where it can be proposed for application to complex 3D non-planar magnets. The main objective of this project is to develop, demonstrate and pre-commercialize a novel, non-planar HTS coil shape that remains superconducting to achieve production scalable reductions in time and cost and performance. The proposed technology is based on the novel concept of a precision sub-scale HTS nonplanar coil assembly. This project focuses on the design, fabrication, material optimization of cable design, and validation and demonstration of the high current carrying capability of superconducting magnets and their support in a complex 3D shape needed for application to stellarator magnetic plasma confinement. The specific objectives of this research program include: (1) The successful application of metal AM to build a precision sub-scale HTS nonplanar coil, (2) An HTS cable and cross-section design that can conform to the required nonplanar coil shape (bend radii as tight as 10-cm) and remains superconducting at an engineering current density of 1.35 kA/cm 2 at 77 K and 1 tesla at the conductor (5 kA in the cable). To achieve the above challenging goals, we have formed a multidisciplinary research team consisting of members from Type One Energy and UW-Madison, MIT PSFC and CFS with complementary skills and strong facilities. The team worked collaboratively on fundamental and applied research on the following three major technical areas: (1) Design, fabrication, and optimization of non-planar HTS Cable The ultimate goal of the project is to determine if commercial REBCO tapes and additive manufacturing can be used to fabricate high field (≥ 10T) non-planar coils with tight bending radii (≃ 100mm) and with a degradation of the critical current (Ic) smaller than 20% with respect to the expected performance. We started with shorter length cable to evaluate the scalability of the production process and eventually reached multiple turns for higher magnetic fields. Our findings suggest that a stellarator coil system of a relevant size, characterized by its asymmetric and non-planar components, can be fabricated using a formed cable in plate method. This system can be simulated using a large-scale modeling approach. The use of hybrid modeling 3 techniques will be pivotal in reducing the complexity of the model and in assessing expected performance in designs. (2) Modeling and simulation of the non-planar HTS Cable Multiphysics simulations are performed using the commercial software and are carried out in self-field conditions, involving 2D and 3D models and twisted around one slot of twist-pitched VIPER cable. Multiphysics simulations are mainly focused on the coil for the critical current evaluation, the magnetic field map, self-Lorentz forces and mechanical, and magnetothermal behavior and the quench dynamics. The detailed model and prediction of the superconducting performance of a stellarator-relevant demonstration cable from numerical simulations supports the results from the actual testing backing the results. A detailed description and results are provided in the later sections. (3) Design, fabrication, and optimization of support for the non-planar HTS Cable The team developed an additive manufactured (AM) coil positioning plate that formed into the required non-planar geometry (with bend radii as tight as 10-cm) and to acceptable tolerances required for a stellarator magnet: (+0.25-mm from ideal on dimensions of coil positioning plates and up to +1-mm from ideal for position of wound coil). The plate materials is also included in this selection process from fabrication and 3D printing perspective and commensurate with eventual application to a fusion reactor. From the cost effectiveness point of view, the HTS coil and plate has the potential to cost less than that made in conventional methods with less waste (<75% waste) reducing time (<50%) and cost (<50%), especially as the AM field matures. The application of advanced manufacturing in the construction of the support plates will also lead to cost reduction, as the cables can be easily replaced, thereby making the assembly modular. With such high primary cost and time savings, high current densities and magnetic field, the funded R&D work has validated the designs, proven the feasibility, and characterized the performance of the HTS coil and plate assembly, paving the way for a relevant-size stellarator coil system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

A Full-Induction Magnetohydrodynamics Solver for Liquid Metal Fusion Blankets in Vertex-CFD

Multiphysics modeling of liquid metal fusion blankets, which produce tritium and convert energy of neutrons created via fusion reactions into heat, is crucial for predicting performance, ensuring structural integrity, and optimizing energy production. While traditional blanket modeling of liquid metal flows during normal steady operating conditions commonly employs the inductionless approximation of the magnetohydrodynamics (MHD) equations, transient scenarios, when the plasma-confining magnetic field varies on millisecond time scales, require a full-induction MHD approach that dynamically evolves the magnetic field via the time-dependent induction equation. This paper presents the formulation, implementation, and initial verification of a full-induction MHD solver integrated within the open-source Vertex-CFD framework, which aims to achieve tight multiphysics coupling, a flexible software design enabling easy extension and addition of physics models, and performance portability across computing platforms. The solver utilizes finite element spatial discretization, implicit Runge–Kutta time integration, and an inexact Newton method to solve the resulting discrete nonlinear system, leveraging Trilinos packages for efficient computation. Verification against selected benchmark problems demonstrates accuracy and robustness of the solver. Furthermore, when the solver is applied to an idealized blanket model in 2.5D and full 3D, results obtained with Vertex-CFD are in good agreement with recently published quasi-2D simulations. These findings establish a computational foundation for future simulations of transient MHD phenomena in liquid metal blankets with Vertex-CFD, and open avenues for future extensions and performance optimizations.

Endeve, Eirik [ORNL] (ORCID:0000000312519507)↗
DOE OSTI2024

Stellarator optimization with constraints

In this work we consider the problem of optimizing a stellarator subject to hard constraints on the design variables and physics properties of the equilibrium. We survey current numerical methods for handling these constraints, and summarize a number of methods from the wider optimization community that have not been used extensively for stellarator optimization thus far. We demonstrate the utility of new methods of constrained optimization by optimizing a quasi-axisymmetric stellarator for favourable physics properties while preventing strong shaping of the plasma boundary, which can be difficult to create with external current sources.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Exploration of the parameter space of quasisymmetric stellarator vacuum fields through adjoint optimisation

Optimising stellarators for quasisymmetry leads to strongly reduced collisional transport and energetic particle losses compared with unoptimised configurations. Although stellarators with precise quasisymmetry have been obtained in the past, it remains unclear how broad the parameter space is where good quasisymmetry may be achieved. We study the range of aspect ratios and rotational transform values for which stellarators with excellent quasisymmetry on the boundary can be obtained. A large number of Fourier harmonics is included in the boundary representation, which is made computationally tractable by the use of adjoint methods to enable fast gradient-based optimisation and by the direct optimisation of vacuum magnetic fields, which converge more robustly compared with solutions from magnetohydrostatics. Several novel configurations are presented, including stellarators with record levels of quasisymmetry on a surface, three field period quasiaxisymmetric stellarators with substantial magnetic shear, and compact quasisymmetric stellarators at low aspect ratios similar to tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A quasi-linear model of electromagnetic turbulent transport and its application to flux-driven transport predictions for STEP

A quasi-linear reduced transport model is developed from a database of high-β electromagnetic nonlinear gyrokinetic simulations performed with spherical tokamak for energy production (STEP) relevant parameters. The quasi-linear model is fully electromagnetic and accounts for the effect of equilibrium flow shear using a novel approach. Its flux predictions are shown to agree quantitatively with predictions from local nonlinear gyrokinetic simulations across a broad range of STEP-relevant local equilibria. This reduced transport model is implemented in the T3D transport solver that is used to perform the first flux-driven simulations for STEP to account for transport from hybrid kinetic ballooning mode turbulence, which dominates over a wide region of the core plasma. Nonlinear gyrokinetic simulations of the final transport steady state from T3D return turbulent fluxes that are consistent with the reduced model, indicating that the quasi-linear model may also be appropriate for describing the transport steady state. Within the assumption considered here, our simulations support the existence of a transport steady state in STEP with a fusion power comparable to that in the burning flat top of the conceptual design, but do not demonstrate how this state can be accessed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Constructing nested coordinates inside strongly shaped toroids using an action principle

A new approach for constructing polar-like boundary-conforming coordinates inside a toroid with strongly shaped cross-sections is presented. A coordinate mapping is obtained through a variational approach, which involves identifying extremal points of a proposed action in the mapping space from [0,2π] 2 ×[0,1] to a toroidal domain in $\mathbb{R}$ 3 . This approach employs an action built on the squared Jacobian and radial length. Extensive testing is conducted on general toroidal boundaries using a global Fourier–Zernike basis via action minimisation. The results demonstrate successful coordinate construction capable of accurately describing strongly shaped toroidal domains. The coordinate construction is successfully applied to the computation of three-dimensional magnetohydrodynamic equilibria in the GVEC code where the use of traditional coordinate construction by interpolation from the boundary failed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The shear Alfvén continuum of quasisymmetric stellarators

The shear Alfvén wave (SAW) continuum plays a critical role in the stability of energetic particle-driven Alfvén eigenmodes. We develop a theoretical framework to analyze the SAW continuum in three-dimensional (3-D) quasisymmetric magnetic fields, focusing on its implications for stellarator design. By employing a near-axis model and degenerate perturbation theory, the continuum equation is solved, highlighting unique features in 3-D configurations, such as the interactions between spectral gaps. Numerical examples validate the theory, demonstrating the impact of flux-surface shaping and quasisymmetric field properties on continuum structure. The results provide insights into optimizing stellarator configurations to minimize resonance-driven losses of energetic particles. This work establishes a basis for incorporating Alfvénic stability considerations into the stellarator design process, demonstrated through optimization of a quasihelical configuration to avoid high-frequency spectral gaps.

fusion plasma↗
DOE OSTI2025

How does ion temperature gradient turbulence depend on magnetic geometry? Insights from data and machine learning

Magnetic geometry has a significant effect on the level of turbulent transport in fusion plasmas. Here, we model and analyse this dependence using multiple machine learning methods and a dataset of >200 000 nonlinear gyrokinetic simulations of ion-temperature-gradient turbulence in diverse non-axisymmetric geometries. The dataset is generated using a large collection of both optimised and randomly generated stellarator equilibria. At fixed gradients and other input parameters, the turbulent heat flux varies between geometries by several orders of magnitude. Trends are apparent among the configurations with particularly high or particularly low heat flux. Regression and classification techniques from machine learning are then applied to extract patterns in the dataset. Due to a symmetry of the gyrokinetic equation, the heat flux and regressions thereof should be invariant to translations of the raw features in the parallel coordinate, similar to translation invariance in computer vision applications. Multiple regression models including convolutional neural networks (CNNs) and decision trees can achieve reasonable predictive power for the heat flux in held-out test configurations, with highest accuracy for the CNNs. Using Spearman correlation, sequential feature selection and Shapley values to measure feature importance, it is consistently found that the most important geometric lever on the heat flux is the flux surface compression in regions of bad curvature. The second most important geometric feature relates to the magnitude of geodesic curvature. These two features align remarkably with surrogates that have been proposed based on theory, while the methods here allow a natural extension to more features for increased accuracy. The dataset, released with this publication, may also be used to test other proposed surrogates, and we find that many previously published proxies do correlate well with both the heat flux and stability boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Bounce-averaged theory in arbitrary multi-well plasmas: solution domains and the graph structure of their connections

Bounce-averaged theories provide a framework for simulating relatively slow processes, such as collisional transport and quasilinear diffusion, by averaging these processes over the fast periodic motions of a particle on a closed orbit. This procedure dramatically increases the characteristic time scale and reduces the dimensionality of the modelled system. The natural coordinates for such calculations are the constants of motion (COM) of the fast particle motion, which by definition do not change during an orbit. However, for sufficiently complicated fields – particularly in the presence of local maxima of the electric potential and magnetic field – the COM are not sufficient to specify the particle trajectory. In such cases, multiple domains in COM space must be used to solve the problem, with boundary conditions enforced between the domains to ensure continuity and particle conservation. Previously, these domains have been imposed by hand, or by recognising local maxima in the fields, limiting the flexibility of bounce-averaged simulations. Here, we present a general set of conditions for identifying consistent domains and the boundary condition connections between the domains, allowing the application of bounce-averaged theories in arbitrarily complicated and dynamically evolving electromagnetic field geometries. We also show how the connections between the domains can be represented by a directed graph, which can help to succinctly represent the trajectory bifurcation structure.

fusion plasma↗