Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Magnetic Fusion Energy”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Cooling flow regime of a plasma thermal quench

A large class of Laboratory, Space, and Astrophysical plasmas is nearly collisionless. When a localized energy or particle sink, for example, in the form of a radiative cooling spot or a black hole, is introduced into such a plasma, it can trigger a plasma thermal collapse, also known as a thermal quench in tokamak fusion. Here we show that the electron thermal conduction in such a nearly collisionless plasma follows the convective energy transport scaling in itself or in its spatial gradient, due to the constraint of ambipolar transport. As a result, a robust cooling flow aggregates mass toward the cooling spot and the thermal collapse of the surrounding plasma takes the form of four propagating fronts that originate from the radiative cooling spot, along the magnetic field line in a magnetized plasma. The slowest one, which is responsible for deep cooling, is a shock front.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Turbulent field fluctuations in gyrokinetic and fluid plasmas

A key uncertainty in the design and development of magnetic confinement fusion energy reactors is predicting edge plasma turbulence. An essential step in overcoming this uncertainty is the validation in accuracy of reduced turbulent transport models. Drift-reduced Braginskii two-fluid theory is one such set of reduced equations that has for decades simulated boundary plasmas in experiment, but significant questions exist regarding its predictive ability. To this end, using a novel physics-informed deep learning framework, we demonstrate the first ever direct quantitative comparisons of turbulent field fluctuations between electrostatic two-fluid theory and electromagnetic gyrokinetic modeling with good overall agreement found in magnetized helical plasmas at low normalized pressure. This framework presents a new technique for the numerical validation and discovery of reduced global plasma turbulence models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Turbulent field fluctuations in gyrokinetic and fluid plasmas

A key uncertainty in the design and development of magnetic confinement fusion energy reactors is predicting edge plasma turbulence. An essential step in overcoming this uncertainty is the validation in accuracy of reduced turbulent transport models. Drift-reduced Braginskii two-fluid theory is one such set of reduced equations that has for decades simulated boundary plasmas in experiment, but significant questions exist regarding its predictive ability. To this end, using a novel physics-informed deep learning framework, we demonstrate the first ever direct quantitative comparisons of turbulent field fluctuations between electrostatic two-fluid theory and electromagnetic gyrokinetic modelling with good overall agreement found in magnetized helical plasmas at low normalized pressure. This framework is readily adaptable to experimental and astrophysical environments, and presents a new technique for the numerical validation and discovery of reduced global plasma turbulence models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Metriplectic foundations of gyrokinetic Vlasov–Maxwell–Landau theory

Here, this Letter reports on a metriplectic formulation of a collisional, nonlinear full-f electromagnetic gyrokinetic theory compliant with energy conservation and monotonic entropy production. In an axisymmetric background magnetic field, the toroidal angular momentum is also conserved. Notably, a new collisional current, contributing to the gyrokinetic Maxwell–Ampère equation and the gyrokinetic charge conservation law, is discovered.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Local conservation of energy in fully implicit PIC algorithms

We consider the issue of strict, fully discrete local energy conservation for a whole class of fully implicit local-charge- and global-energy-conserving particle-in-cell (PIC) algorithms. Earlier studies demonstrated these algorithms feature strict global energy conservation. However, whether a local energy conservation theorem exists (in which the local energy update is governed by a flux balance equation at every mesh cell) for these schemes is unclear. In this study, we show that a local energy conservation theorem indeed exists. We begin our analysis with the 1D electrostatic PIC model without orbit-averaging, and then generalize our conclusions to account for orbit averaging, multiple dimensions, and electromagnetic models (Darwin). In all cases, a temporally, spatially, and particle-discrete local energy conservation theorem is shown to exist, proving that these formulations (as originally proposed in the literature), in addition to being locally charge conserving and globally energy conserving, are strictly locally energy conserving as well. In contrast to earlier proofs of local conservation in the literature, which only considered continuum time, our result is valid for the fully implicit time-discrete version of all models considered, including important features such as orbit averaging. We demonstrate the local-energy-conservation property numerically with a paradigmatic numerical example.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Energy and momentum conservation in the Euler–Poincaré formulation of local Vlasov–Maxwell-type systems

The action principle by Low for the classic Vlasov–Maxwell system contains a mix of Eulerian and Lagrangian variables. This renders the Noether analysis of reparametrization symmetries inconvenient, especially since the well-known energy- and momentum-conservation laws for the system are expressed in terms of Eulerian variables only. While an Euler–Poincaré formulation of Vlasov–Maxwell-type systems, effectively starting with Low's action and using constrained variations for the Eulerian description of particle motion, has been known for a while Cendra et al, it is hard to come by a documented derivation of the related energy- and momentum-conservation laws in the spirit of the Euler–Poincaré machinery. To our knowledge only one such derivation exists in the literature so far, dealing with the so-called guiding-center Vlasov–Darwin system Sugama et al. The present exposition discusses a generic class of local Vlasov–Maxwell-type systems, with a conscious choice of adopting the language of differential geometry to exploit the Euler–Poincaré framework to its full extent. After reviewing the transition from a Lagrangian picture to an Eulerian one, we demonstrate how symmetries generated by isometries in space lead to conservation laws for linear- and angular-momentum density and how symmetry by time translation produces a conservation law for energy density. Here, we also discuss what happens if no symmetries exist. Finally, two explicit examples will be given—the classic Vlasov–Maxwell and the drift-kinetic Vlasov–Maxwell—and the results expressed in the language of regular vector calculus for familiarity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Driven resonant current amplification in self-organized plasma configurations with uniform λ and plasma pressure confinement

Driven resonant current amplification in uniform λ=$^{μoJ}_{B}$, self-organized plasma configurations with plasma pressure confinement is studied. This work extends previous analyses of force-free states (J ∥ B) to equilibria with plasma pressure gradients (∇P = J × B ≠ 0). An analytical treatment shows that in equilibria with uniform λ and a plasma pressure gradient that varies linearly with poloidal magnetic flux, driven resonances are either downshifted or upshifted from their force-free values. Additionally, a modification of off resonance current amplification is observed in pressure confining equilibria. Computational results add support to analytical predictions, and also include driven resonance behavior with equilibrium profiles not considered analytically. Furthermore, these results suggest not only persistence, but also modification of driven resonant current amplification in sustained spheromak configurations with plasma pressure confinement when compared to force-free states. By providing a basis for the existence of driven, high current amplification, pressure confining equilibria, this work provides insight into the accessibility of sustained spheromak states of most interest for potential fusion energy applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Implementation of extreme ultraviolet spectroscopy on a sheared-flow-stabilized Z pinch

A diagnostic for extreme ultraviolet spectroscopy was fielded on the sheared-flow-stabilized (SFS) fusion Z-pinch experiment (FuZE-Q) for the first time. The spectrometer collected time-gated plasma emission spectra in the 5–40 nm wavelength (30–250 eV) range for impurity identification, radiative power studies, and for plasma temperature and density measurements. The unique implementation of the diagnostic included fast (10 ns risetime) pulsed high voltage electronics and a multi-stage differential pumping system that allowed the vacuum-coupled spectrometer to collect three independently timed spectra per FuZE-Q shot while also protecting sensitive internal components. Analysis of line emission identifies oxygen (N-, C-, B-, Be-, Li-, and He-like O), peaking in intensity shortly after maximum current (>500 kA). This work provides a foundation for future high energy spectroscopy experiments on SFS Z-pinch devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Simulations of stand-off runaway electron beam termination by tungsten particulates for tokamak disruption mitigation

Abstract Stand-off runaway electron termination by injected tungsten particulates offers a plausible option in the toolbox of disruption mitigation. Tungsten is an attractive material choice for this application due to large electron stopping power and high melting point. To assess the feasibility of this scheme, we simulate runaway collisions with tungsten particulates using the MCNP program for incident runaway energies ranging from 1 to 10 MeV. We assess runaway termination from energetics and collisional kinematics perspectives. Energetically, the simulations show that 99% of runaway beam energy is removed by tungsten particulates on a timescale of 4–9 µ s. Kinematically, the simulations show that 99% of runaways are terminated by absorption or backscattering on a timescale of 3–4 µ s. By either metric, the runaway beam is effectively terminated before the onset of particulate melting. Furthermore, the simulations show that secondary radiation emission by tungsten particulates does not significantly impact the runaway termination efficacy of this scheme. Secondary radiation is emitted at lower particle energies than the incident runaways and with a broad angular distribution such that the majority of secondary electrons emitted will not experience efficient runaway re-acceleration. Overall, the stand-off runaway termination scheme is a promising concept as a last line of defense against runaway damage in ITER, SPARC, and other future burning-plasma tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Determination of a macro- to micro-scale progression leading to a magnetized plasma disruption

We report here the observations of a plasma jet evolving through a macro- to micro-scale progression sequence. This leads to a fast magnetic reconnection that results in the jet breaking off from its originating electrode and forming a force-free state. A sausage-like pinching occurs first and squeezes an initially fat, short magnetized jet so that it becomes thin. The thin jet then becomes kink unstable. The lengthening of the jet by the kinking thins the jet even more since the kink is an incompressible instability. When the jet radius becomes comparable to the ion-skin depth, Hall and electron inertial physics become important and establish the environment for a fast magnetic reconnection. This fast reconnection occurs, disrupting the jet and establishing a force-free state. X-ray bursts and whistler waves, evidence of a magnetic reconnection, are observed when the plasma jet breaks off from the electrode. This experimentally confirmed sequence of successive thinning from pinching followed by kinking is reproduced in a three-dimensional ideal Magnetohydrodynamic (MHD) numerical simulation. The results of the experiment and the numerical simulation, together demonstrate a viable path from macro-scale MHD physics to micro-scale non-MHD physics where fast reconnection occurs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Approximate symmetries of guiding-centre motion

In a strong, inhomogeneous magnetic field, charged particle dynamics may be studied in the guiding-centre approximation, which is known to be Hamiltonian. When the magnetic field is quasisymmetric, the first-order guiding-centre (FGC) Hamiltonian structure admits a continuous symmetry, and therefore a conserved quantity in addition to the energy. Since the FGC system is only an approximation, it is also interesting to consider approximate symmetries of the guiding-centre Hamiltonian structure. We find that any approximate spatial symmetry coincides with quasisymmetry to leading order. For approximate phase-space symmetries, we derive weaker conditions than quasisymmetry. The latter include 'weak quasisymmetry' as a subcase, recently proposed by Rodríguez et al. Our results, however, show that weak quasisymmetry is necessarily non-spatial at first order. Finally, we demonstrate that if the magnetic field is constrained to satisfy magnetohydrostatic force balance then an approximate symmetry must agree with quasisymmetry to leading order.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

AC compensation of 3D magnetic diagnostic signals in DIII-D and National Spherical Torus Experiment-Upgrade (NSTX-U) for real-time application

A time domain algorithm has been developed to remove the vacuum pickup generated by both coil current (DC) and induced vessel current (AC) in real time from three dimensional (3D) magnetic diagnostic signals in the National Spherical Torus Experiment-Upgrade (NSTX-U) and DIII-D tokamaks. The possibility of detecting 3D plasma perturbations in real time is essential in modern and future tokamaks to avoid and control MHD instabilities. The presence of vacuum field pickup, due to toroidally asymmetric (3D) coils or to misalignment between sensors and axisymmetric (2D) coils, pollutes the measured plasma 3D field, making the detection of the magnetic field produced by the plasma challenging. Although the DC coupling between coils and sensors can be easily calculated and removed, the AC part is more difficult. Here, an algorithm based on a layered low-pass filter approach for the AC compensation and its application for DIII-D and NSTX-U data is presented, showing that this method reduces the vacuum pickup to the noise level. Comparison of plasma response measurements with and without vacuum compensation shows that accurate mode locking detection and plasma response identification require precise AC and DC compensations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Serpent neutronics model of Wendelstein 7-X for 14.1 MeV neutrons

In this work, a Serpent 2 neutronics model of the Wendelstein 7-X (W7-X) stellarator is prepared, and an response function for the Scintillating-Fibre neutron detector (SciFi) is calculated using the model. The neutronics model includes the simplified geometry for the key components of the stellarator itself as well as the torus hall. The objective of the model is to assess the 14.1 MeV neutron flux from deuteron-triton fusions in W7-X, where the neutrons are modelled only until they have slowed down to 1 MeV energy. The key messages of this article are: demonstration of unstructured mesh geometry usage for stellarators, W7-X in particular; technical documentation of the model and first insights in fast neutron behaviour in W7-X, especially related to the SciFi: the model indicates that the superconducting coils are the strongest scatterers and block neutrons from large parts of the plasma. The back-scattering from e.g. massive steel support structures is found to be small. The SciFi will detect neutrons from an extended plasma volume in contrast to having an effective line-of-sight.

14 MeV↗
DOE OSTI2024

Corrigendum: Simulations of stand-off runaway electron beam termination by tungsten particulates for tokamak disruption mitigation (2024 Nucl. Fusion 64 056019)

In the original article, we inadvertently omitted a reference which previously introduced the concept of tungsten injection for 'Disruption Mitigation in Tokamak Reactor via Reducing the Seed Electrons of Avalanche.' In that reference, the authors have proposed injection of a tungsten cylinder of 80 mm by 8 mm dimensions, by rail guns, to absorb the runaway seed population. This can be contrasted with another interesting idea of using small tungsten pellets coated with a low Z material for depleting the runaway seeds as an option for upgrade of the ITER Disruption Mitigation System. Interested readers are referred to those two papers on the tungsten injection for runaway seed removal/reduction, and the related issue of not shortening the current quench excessively. The application focus of was not on reducing the runaway seeds in the initial Ohmic-to-runaway current conversion phase, but on safe termination of a fully formed and likely decaying runaway beam that is about to scrape off against the first wall, which can be accelerated by vertical displacement events. This scheme thus serves as a last line of defense against potential wall damage. The choice here is locally released tungsten particulates, similar to the previous idea of using tungsten particulates as a dust shield for the divertor. The most significant finding of is that strong pitch angle scattering, in addition to energy attenuation and absorption, of high-energy runaways by tungsten particulates, provides a transport mechanism by which runaways damage of the first wall can be mitigated by reducing and spreading the runaway wall load. We regret the omission in the original article, and thank the authors of the cited article for bringing this matter to our attention. We also find the possibility of using the injected tungsten rod from to terminate the runaways, as opposed to the cloud of tungsten particulates in intriguing. In that context, it is of interest to mention another possibility of a retractable tungsten metal arm that can be swung out for deployment in runaway termination, which would remove the need for post-mitigation recovery of tungsten debris in injection schemes.

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

Finite spatial-grid effects in energy-conserving particle-in-cell algorithms

Finite-grid (or aliasing) instabilities are pervasive in particle-in-cell (PIC) plasma simulation algorithms, and force the modeler to resolve the smallest (Debye) length scale in the problem regardless of dynamical relevance. These instabilities originate in the aliasing of interpolation errors between mesh quantities and particles (which live in the space–time continuum). Recently, strictly energy-conserving PIC (EC-PIC) algorithms have been developed that promise enhanced robustness against aliasing instabilities. In this study, we confirm by analysis that EC-PIC is stable against aliasing instabilities for stationary plasmas. For drifting plasmas, we demonstrate by analysis and numerical experiments that, while EC-PIC algorithms are not free from these instabilities in principle, they feature a benign stability threshold for finite-temperature plasmas that make them usable in practice for a large class of problems (featuring ambipolarity and realistic ion-electron mass ratios) without the need to consider the size of the Debye length. We also demonstrate that this threshold is absent for the popular momentum-conserving PIC algorithms, which are therefore unstable for both drifting and stationary plasmas beyond a threshold in cell size compared to Debye length. Finite-grid (or aliasing) instabilities are pervasive in particle-in-cell (PIC) plasma simulation algorithms, and force the modeler to resolve the smallest (Debye) length scale in the problem regardless of dynamical relevance. These instabilities originate in the aliasing of interpolation errors between mesh quantities and particles (which live in the space–time continuum). Recently, strictly energy-conserving PIC (EC-PIC) algorithms have been developed that promise enhanced robustness against aliasing instabilities. In this study, we confirm by analysis that EC-PIC is stable against aliasing instabilities for stationary plasmas. For drifting plasmas, we demonstrate by analysis and numerical experiments that, while EC-PIC algorithms are not free from these instabilities in principle, they feature a benign stability threshold for finite-temperature plasmas that make them usable in practice for a large class of problems (featuring ambipolarity and realistic ion-electron mass ratios) without the need to consider the size of the Debye length. Finally, we also demonstrate that this threshold is absent for the popular momentum-conserving PIC algorithms, which are therefore unstable for both drifting and stationary plasmas beyond a threshold in cell size compared to Debye length.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Further modeling of q 95 windows for the suppression of edge localized modes by resonant magnetic perturbations in the DIII-D tokamak

An improved resonant plasma response model that more accurately captures the physics of the interaction between a tokamak plasma and a resonant magnetic perturbation (RMP) is developed. The model interpolates between the linear and the nonlinear response regimes and takes into account the fact that the slip-frequency is non-zero in the nonlinear regime. The improved model is incorporated into the extended perturbed equilibrium code (EPEC) toroidal asymptotic matching code. The modified EPEC code is used to investigate RMP-induced edge-localized-mode (ELM) suppression in DIII-D H-mode discharge #145380. Somewhat surprisingly, allowing for a finite slip-frequency (i.e., relaxing the so-called no-slip constraint) is found to only slightly facilitate the locking of driven magnetic island chains to the RMP, and, hence, to only slightly facilitate RMP-induced ELM suppression. This is true despite the fact that the nature of non-locked island solutions is radically different when the no-slip constraint is imposed compared to when it is relaxed (in the first case, the widths of the island chains driven at the rational surfaces pulsate, and in the second case, they remain steady). The previously obtained conclusion that the response of a typical H-mode tokamak plasma to an RMP cannot be accurately modeled by linear theory is confirmed. Here, the previously obtained conclusion that the best agreement between theory and observations is achieved by assuming that the natural frequencies of tearing modes, in the absence of the RMP, are determined by the local equilibrium E × B velocity is also confirmed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Large radiation back-flux from Monte Carlo simulations of fusion neutron–material interactions

Abstract Fusion power reactors will generate intense neutron fluxes into plasma-facing and structural materials (SMs). Radiation back-fluxes, generated from neutron–material interactions under these fluxes, can dramatically impact the plasma dynamics, e.g. by seeding runaway electrons during disruptions via Compton scattering of background electrons by wall-emitted gamma radiation. Here, we quantify these back-fluxes, including neutrons, gamma rays, and electrons, using Monte Carlo calculations for a range of SM candidates and first wall (FW) thicknesses. The radiation back-flux magnitudes are remarkably large, with neutron and gamma radiation back-fluxes on the same order of magnitude as the incident fusion neutron flux. Electron back-fluxes are two orders of magnitudes lower, but are emitted at sufficiently high energies to impact the sheath and boundary plasma dynamics. Material configuration plays a key role in determining back-flux magnitudes. The SM chiefly determines the neutron back-flux magnitude, while the FW thickness principally attenuates the gamma ray and electron back-fluxes. In addition to prompt back-fluxes, which are emitted immediately after fusion neutrons impact the surface, significant delayed gamma ray and electron back-fluxes arise from nuclear decay processes in the activated materials. These delayed back-flux magnitudes range from 2% to 7% of the prompt back-fluxes, and remain present during transients when fusion no longer occurs. During disruptions, build-up of delayed gamma radiation back-flux represents potential runaway electron seeding mechanisms, posing additional challenges for disruption mitigation in a power reactor compared with non-nuclear plasma operations. This work highlights the impact of these radiation back-fluxes plasma performance and demonstrates the importance of considering back-flux generation in materials selection for fusion power reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗