SEARCH · Search NASA

Results for “Fusion technology”

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.

SourceResearch connectionFrom yearThrough year

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

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Finite orbit width effects in large aspect ratio stellarators

New orbit-averaged equations for low collisionality neoclassical fluxes in large aspect ratio stellarators with mirror ratios close to unity are derived. The equations retain finite orbit width effects by employing the second adiabatic invariant $J$ as a velocity-space coordinate and they have been implemented in the orbit-averaged neoclassical code KNOSOS (Velasco et al. , J. Comput. Phys. , vol. 418, 2020, 109512; Velasco et al. , Nucl. Fusion , vol. 61, 2021, 116013). The equations are used to study the $1/\nu$ regime and the lower collisionality regimes. For generic large aspect ratio stellarators with mirror ratios close to unity, as the collision frequency decreases, the $1/\nu$ regime transitions directly into the $\nu$ regime, without passing through a $\sqrt {\nu }$ regime. An explicit formula for the neoclassical fluxes in the $\nu$ regime is obtained. The formula includes the effect of particles that transition between different types of wells. While these transitions produce stochastic scattering independent of the value of the collision frequency in velocity space, the diffusion in real space remains proportional to the collision frequency. The $\sqrt {\nu }$ regime is only recovered in large aspect ratio stellarators close to omnigeneity: large aspect ratio stellarators with large mirror ratios and optimized large aspect ratio stellarators with mirror ratios close to unity. Neoclassical transport in large aspect ratio stellarators with large mirror ratios can be calculated with the orbit-averaged equations derived by Calvo et al. ( Plasma Phys. Control. Fusion , vol. 59, 2017, 055014). In these stellarators, the $\sqrt {\nu }$ regime exists in the collisionality interval $(a/R) \ln (R/a) \ll \nu _{ii} R a/\rho _i v_{ti} \ll R/a$ . In optimized large aspect ratio stellarators with mirror ratios close to unity, the $\sqrt {\nu }$ regime occurs in an interval of collisionality that depends on the deviation from omnigeneity $\delta$ : $\delta ^{2} |\ln \delta | \ll \nu _{ii} R a/\rho _i v_{ti} \ll 1$ . Here, $\nu _{ii}$ is the ion–ion collision frequency, $\rho _i$ and $v_{ti}$ are the ion gyroradius and thermal speed, and $a$ and $R$ are the minor and major radii.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Modelling the nonlinear plasma response to externally applied three-dimensional fields with the Stepped Pressure Equilibrium Code

Small-amplitude, symmetry-breaking magnetic field perturbations, including resonant magnetic perturbations (RMPs) and error fields, can profoundly impact plasma properties in both tokamaks and stellarators. In this work, we perform the first comparison between the Stepped Pressure Equilibrium Code (SPEC) (a comparatively fast and efficient equilibrium code based on energy-minimisation principles) and M3D-C 1 (a high-fidelity albeit computationally expensive initial-value extended-magnetohydro-dynamic (MHD) code) to assess the conditions under which SPEC can be used to model the nonlinear, non-ideal plasma response to an externally applied (m = 2,n = 1) RMP field in an experimentally relevant geometry. We find that SPEC is able to capture the plasma response in the weakly nonlinear regime – meaning perturbation amplitudes below the threshold for break up of the separatrix and onset of secondary magnetic island formation – when around half of the total toroidal flux is enclosed in the volume containing the q = 2 resonant surface. Further, the observed dependence of SPEC solutions on input parameters, including toroidal flux and the number of volumes into which the plasma is partitioned, indicates that additional exploration of the underlying Multi-Region Relaxed MHD physics model is needed to constrain the choice of parameters. Nonetheless, this work suggests promising applications of SPEC to optimisation and fusion plasma design.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Thermal ion kinetic effects and Landau damping in fishbone modes

We report the kinetic–magnetohydrodynamic (MHD) hybrid simulation approach for macroscopic instabilities in plasmas can be extended to include the kinetic effects of both thermal ions and energetic ions. The new coupling scheme includes synchronization of the density and parallel velocity between thermal ions and MHD, in addition to pressure coupling, to ensure the quasineutrality condition and avoid numerical errors. The new approach has been implemented in the kinetic-MHD code M3D-C1-K, and was used to study the thermal ion kinetic effects and Landau damping in fishbone modes in both DIII-D and NSTX. It is found that the thermal ion kinetic effects can cause an increase of the frequencies of the non-resonant n = 1 fishbone modes driven by energetic particles for q min >1, and Landau damping can provide additional stabilization effects. A nonlinear simulation for n = 1 fishbone mode in NSTX is also performed, and the perturbation on magnetic flux surfaces and the transport of energetic particles are calculated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Database-wide hazard modelling of the onset of DIII-D tearing modes with field features

The rate of onset (hazard) of tearing modes is modelled probabilistically using statistical learning algorithms. Axisymmetric energy-density equilibrium fields are taken as raw high-dimensional input features which are reduced with principal component analysis. Signal processing of non-axisymmetric magnetics fluctuation array data provides the target information from which to learn. Model selection, visualization and calibration assessment procedures are detailed. Here, the analysis is deployed at large scale across the DIII-D tokamak database. Standard model selection criteria suggest that the energy-density post-processed feature is a better choice for modelling the onset rate compared to the non-processed equilibrium reconstruction solution. Two example applications of the learned rate function are demonstrated: (i) proximity-to-onset discharge monitoring and (ii) database analysis showing an (expected) observational global trend that the general hazard increases as a plasma performance metric increases. An important connection between the hazard function and its use as a conditional probability generator is reviewed in the Appendix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Momentum transfer and plasma rotation caused by destabilized eigenmodes in tokamaks

In this work, the influence of magnetohydrodynamic eigenmodes destabilized by energetic ions on the momentum of these ions and concomitant sheared plasma rotation are studied. Two mechanisms affecting rotation are revealed: (i) spatial channelling (SC) – radially separated emission and absorption of the momentum; (ii) mode induced redistribution (MIR) across the magnetic field of the momentum of energetic ions. Forces arising during SC and MIR produce both toroidal and poloidal rotations. In addition, the momentum emission during SC leads to a radial flux of fast ions and generation of a radial electric field. Using the developed theory, estimates were made for the ITER (International Thermonuclear Experimental Reactor) 15 MA baseline scenario. They show that a global toroidicity-induced Alfvén eigenmode destabilized by alpha particles and neutral beam injection can result in significant radial electric field and forces applied to plasma. However, available data are not sufficient for a reliable prediction of the effects of SC and MIR in ITER. In general, one can expect that sheared rotation arising after destabilization of Alfvén modes and fast magnetoacoustic modes by energetic ions will tend to suppress the turbulence and improve plasma performance. The importance of plasma rotation is supported, in particular, by the fact that during the JET DTE1 experimental campaign the best parameters were achieved in a deuterium–tritium discharge where the rotation frequency was largest.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

A general infrastructure for data-driven control design and implementation in tokamaks

A general infrastructure for tokamak controllers based on data-driven neural net models is presented. The paradigm allows for more flexible choices of both the underlying model and the desired controlled variables and targets. The system is implemented and tested on the DIII-D tokamak, enacting simultaneous pressure and temperature control via a finite-set model-predictive controller. Traditional control methods such as proportional–integral–derivative (PID) have proven effective for decoupled control tasks, but scale poorly when trying to achieve more complicated goals such as full state control. This is exactly where model-based controllers succeed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Merging of the superbanana plateau and $\sqrt{v}$ transport regimes in nearly quasisymmetric stellarators

We report alpha particle confinement is one of the most demanding issues for stellarators. It now seems clear that it is possible to design optimized stellarators that confine the background plasma at near tokamak radial transport levels. Moreover, adequate collisionless alpha particle confinement is possible in the core of a highly optimized stellarator. Here, the collisional confinement of barely trapped alphas in an optimized stellarator is considered by accounting for the resonance due to the reversal in direction of the drift within a flux surface and investigating the sensitive role of magnetic shear in keeping this resonance close to the passing boundary in some nearly quasisymmetric stellarator configurations. The treatment relies on a narrow collisional boundary layer formulation that combines the responses of both these resonant pitch angle alphas and the remaining barely trapped alphas. A novel merged regime treatment leads to explicit expressions for the energy diffusivity for both superbanana plateau (or resonant plateau) and $\sqrt{v}$ transport in the large aspect ratio limit for a slowing down tail alpha distribution function, where v is the effective pitch angle scattering collision frequency of the trapped alphas off the background ions. Depending on the details of the optimization scheme and the sign of the magnetic shear, modest magnetic shear can be used to reduce superbanana (or resonant) plateau transport to below the $\sqrt{v}$ transport level. In addition, a quasilinear equation retaining spatial diffusion is derived for a general alpha distribution function that allows the radial alpha transport to modify the distribution so it is no longer isotropic in velocity space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Microstability of β ~ 1 tokamak equilibria

High-power-density tokamaks offer a potential solution to design cost-effective fusion devices. One way to achieve high power density is to operate at a high β value (the ratio of thermal to magnetic pressure), i.e. β~1 . However, a β~1 state may be unstable to various pressure- and current-driven instabilities or have unfavourable microstability properties. To explore these possibilities, we generate β~1 equilibria and investigate their stability. First, we demonstrate the generation of high-β equilibria with the computer code VMEC. We then analyse these equilibria to determine their stability against the infinite-n ideal-ballooning mode. We follow that by engaging in a detailed microstability study using the GS2 code, beginning with assessments of electrostatic ion-temperature-gradient and trapped election mode instabilities. We observe interesting behaviour for the high-β equilibria – stabilization of these modes through two distinct mechanisms – large negative local shear and reversal of electron precession drift. Finally, we perform electromagnetic gyrokinetic simulations and observe enhanced stability in the outer core of high-β equilibria and absence of kinetic ballooning modes in the negative-triangularity, high-β equilibria. The enhanced outer-core stability of high-β equilibria is different from their lower-β counterparts and offers an alternative, potentially favourable regime of tokamak operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Electron kinetics in a high- Z plasmoid

The problem of the electron dynamics on a closed magnetic field line passing through a high-Z plasmoid is considered. The electron kinetic equation is integrated over bounce motion and pitch angle, reducing the independent variables to a single adiabatic invariant plus time. Integration of the full Landau self-collision operator is carried out exactly, resulting in a nonlinear integro-differential operator in the new invariant. Conservation laws and the H theorem of the integrated self-collision operator are proven. Numerical solutions of the integrated kinetic equation are obtained with a self-consistent quasineutral electric potential, given the initial condition of a cold plasmoid immersed in a hot ambient plasma. The fact that cold electrons are deeply trapped in a potential with a parabolic peak leads to exactly 3/4 the usual rate of collisional heating by the ambient plasma, independent of any other parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

A fast, accurate and easy to implement Kapur–Rokhlin quadrature scheme for singular integrals in axisymmetric geometries

Many applications in magnetic confinement fusion require the efficient calculation of surface integrals with singular integrands. The singularity subtraction approaches typically used to handle such singularities are complicated to implement and low-order accurate. In contrast, we demonstrate that the Kapur–Rokhlin quadrature scheme is well-suited for the logarithmically singular integrals encountered for a toroidally axisymmetric confinement system, is easy to implement and is high-order accurate. As an illustration, we show how to apply this quadrature scheme for the efficient and accurate calculation of the normal component of the magnetic field due to the plasma current on the plasma boundary, via the virtual-casing principle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

The DESC stellarator code suite Part 3: Quasi-symmetry optimization

The DESC stellarator optimization code takes advantage of advanced numerical methods to search the full parameter space much faster than conventional tools. Only a single equilibrium solution is needed at each optimization step thanks to automatic differentiation, which efficiently provides exact derivative information. A Gauss–Newton trust-region optimization method uses second-order derivative information to take large steps in parameter space and converges rapidly. With just-in-time compilation and GPU portability, high-dimensional stellarator optimization runs take orders of magnitude less computation time with DESC compared to other approaches. This paper presents the theory of the DESC fixed-boundary local optimization algorithm along with demonstrations of how to easily implement it in the code. Example quasi-symmetry optimizations are shown and compared to results from conventional tools. Three different forms of quasi-symmetry objectives are available in DESC, and their relative advantages are discussed in detail. In the examples presented, the triple product formulation yields the best optimization results in terms of minimized computation time and particle transport. This paper concludes with an explanation of how the modular code suite can be extended to accommodate other types of optimization problems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

The DESC stellarator code suite. Part 1. Quick and accurate equilibria computations

Three-dimensional equilibrium codes are vital for stellarator design and operation, and high-accuracy equilibria are also necessary for stability studies. This paper details comparisons of two three-dimensional equilibrium codes: VMEC, which uses a steepest-descent algorithm to reach a minimum-energy plasma state, and DESC, which minimizes the magnetohydrodynamic (MHD) force error in real space directly. Accuracy as measured by satisfaction of MHD force balance is presented for each code, along with the computation time. It is shown that DESC is able to achieve more accurate solutions, especially near axis. The importance of higher-accuracy equilibria is shown in DESC's better agreement of stability metrics with asymptotic formulae. DESC's global Fourier–Zernike basis also yields solutions with analytic derivatives explicitly everywhere in the plasma volume, provides improved accuracy in the radial direction versus conventional finite differences and allows for exponential convergence. Further, DESC can compute a solution with the same accuracy as VMEC in order-of-magnitude less time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Benchmarking magnetised three-wave coupling for laser backscattering: analytic solutions and kinetic simulations

Understanding magnetised laser–plasma interactions is important for controlling magneto-inertial fusion experiments and developing magnetically assisted radiation and particle sources. For nanosecond pulses at non-relativistic intensities, interactions are dominated by coherent three-wave interactions, whose nonlinear coupling coefficients became known only recently when waves propagate at oblique angles with the magnetic field. In this paper, backscattering coupling coefficients predicted by warm-fluid theory are benchmarked using particle-in-cell simulations in one spatial dimension, and excellent agreements are found for a wide range of plasma temperatures, magnetic field strengths and laser propagation angles, when the interactions are mediated by electron-dominant hybrid waves. Systematic comparisons between theory and simulations are made possible by a rigorous protocol. On the theory side, the initial boundary value problem of linearised three-wave equations is solved, and the transient-time solutions allow the effects of growth and damping to be distinguished. On the simulation side, parameters are carefully chosen and calibration runs are performed to ensure that comparisons are well controlled. Fitting simulation data to analytical solutions yields numerical growth rates that match theory predictions within error bars. Although warm-fluid theory is found to be valid for a wide parameter range, genuine kinetic effects have also been observed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Neoclassical transport in strong gradient regions of large aspect ratio tokamaks

We present a new neoclassical transport model for large aspect ratio tokamaks where the gradient scale lengths are of the size of the ion poloidal gyroradius. Previous work on neoclassical transport across transport barriers assumed large density and potential gradients but a small temperature gradient, or neglected the gradient of the mean parallel flow. Using large aspect ratio and low collisionality expansions, we relax these restrictive assumptions. We define a new set of variables based on conserved quantities, which simplifies the drift kinetic equation whilst keeping strong gradients, and derive equations describing the transport of particles, parallel momentum and energy by ions in the banana regime. The poloidally varying parts of density and electric potential are included. Studying contributions from both passing and trapped particles, we show that the resulting transport is dominated by trapped particles. We find that a non-zero neoclassical particle flux requires parallel momentum input which could be provided through interaction with turbulence or impurities. We derive upper and lower bounds for the energy flux across a transport barrier in both temperature and density and present example profiles and fluxes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Design and construction of the MUSE permanent magnet stellarator

This paper documents the design and construction of MUSE, the world's first permanent magnet (PM) stellarator and the first quasi-axisymmetric experiment. The purpose of MUSE is to develop and assess a new way of building optimised stellarators that uses simple planar coils PMs. Our PM optimisation algorithm consists of initialising a geometry to pack dipoles densely, running the FAMUS code to minimise surface field error subject to PM constraints and applying discrete jumps to reach a physically realisable solution. FAMUS treats the PM system as a set of ideal point dipoles. From there we construct finite-volume magnet towers to be housed in 3D-printed PM holders. We describe the design of the PM holders, which were validated by laser metrology. We analyse the effects of finite permeability, sensitivity to perturbations and magnetostatic forces. An exact analytic formula for the magnetic field from a finite-volume PM tower is presented to compute PM–PM forces and stress on the PM holder. Stellarator construction is complete and experiments are underway.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Fast correlation heating in moderately coupled electron–ion plasmas

If the electrons in a plasma are suddenly heated, the resulting change in Debye shielding causes the ion kinetic energy to quickly increase. For the first time, this correlation heating , which is much faster than collisional energy exchange, is rigorously derived for a moderately coupled, electron–ion plasma. The electron–ion mass ratio is taken to be the smallest parameter in the Bogoliubov–Born–Green–Kirkwood–Yvon hierarchy, smaller even than the reciprocal of the plasma parameter. This ordering differs from conventional kinetic theory by making the electron collision rates faster than the ion plasma frequency, which allows stronger coupling and makes the ion heating a function only of the total energy supplied to the electrons. The calculation uses known formulae for correlations in a two-temperature plasma, for which a new, elementary derivation is presented. Suprathermal ions may be created more rapidly by this mechanism than by ion–electron Coulomb collisions. This means that the use of a femtosecond laser pulse could potentially help to achieve ignition in certain fast ignition approaches to inertial confinement fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

An adjoint-based method for optimising MHD equilibria against the infinite- n , ideal ballooning mode

We demonstrate a fast adjoint-based method to optimise tokamak and stellarator equilibria against a pressure-driven instability known as the infinite-n ideal ballooning mode. We present three finite-β (the ratio of thermal to magnetic pressure) equilibria: one tokamak equilibrium and two stellarator equilibria that are unstable against the ballooning mode. Using the self-adjoint property of ideal magnetohydrodynamics, we construct a technique to rapidly calculate the change in the eigenvalue, a measure of ideal ballooning instability. Using the SIMSOPT optimisation framework, we then implement our fast adjoint gradient-based optimiser to minimise the eigenvalue and find stable equilibria for each of the three originally unstable equilibria.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Fast-ion transport in quasisymmetric equilibria in the presence of a resonant Alfvénic perturbation

Significant progress has been made in designing magnetic fields that provide excellent confinement of the guiding-centre trajectories of alpha particles using quasisymmetry (QS). Given the reduction in this transport channel, we assess the impact of resonant Alfvén eigenmodes (AEs) on the guiding-centre motion. The AE amplitudes are chosen to be consistent with experimental measurements and large-scale simulations. We evaluate the drift resonance condition, phase-space island width and island overlap criterion for quasisymmetric configurations. Kinetic Poincaré plots elucidate features of the transport, including stiff transport above a critical perturbation amplitude. Our analysis highlights key departures from the AE-driven transport in tokamaks, such as the avoidance of phase-space island overlap in quasihelical configurations and the enhanced transport due to wide phase-space islands in low magnetic shear configurations. In configurations that are closer to QS, with QS deviations $\delta B/B_0 \lesssim 10^{-3}$ , the transport is primarily driven by the AE, while configurations that are further from QS, $\delta B/B_0 \sim 10^{-2}$ , experience significant transport due to the QS-breaking fields in addition to the AE.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗