Search NASA⌕ Search

SEARCH · Search NASA

Results for “Magnetic confinement fusion”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Use of intrinsic hysteresis for the active control of internal transport barriers in magnetically confined fusion plasmas

In magnetically confined fusion devices, control of internal transport barriers (ITBs) is important both to enhance and suppress the turbulent transport to improve confinement control. Barrier control should allow for the improvement of confinement to aid in achieving the needed fusion criteria while also permitting the degradation of confinement to control profiles and clean the device by moving out the impurities accumulated near the core. In this work, we present a novel control scenario that takes advantage of the hysteresis intrinsic to transport barriers to easily cycle through enhanced and degraded confinement regimes. The control scenario is illustrated using a five-field simplified transport model for an ITB using typical parameters of a neutral beam injection-heated DIII-D tokamak discharge. Pellets and ion cyclotron resonance frequency power are used as control knobs for this active control scenario. These knobs adequately modify at will the local gradients and, therefore, the growth rates and shearing rates, allowing for easy and efficient control of the barrier by taking advantage of the barrier hysteresis. The result is a control cycle that could be operated with a relatively small amount of power in high performance regimes which, nowadays, typically require large power to control. It may also have advantages to avoid, or at least ameliorate, the appearance of magnetohydrodynamic instabilities in the barrier region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Iso-nuclear tungsten dielectronic recombination rates for use in magnetically-confined fusion plasmas

Under the auspices of the IAEA Atomic and Molecular Data Center and the Korean Atomic Energy Research Institute, our assembled group of authors has reviewed the current state of dielectronic recombination (DR) rate coefficients for various ion stages of tungsten (W). Subsequent recommendations were based upon available experimental data, first-principle calculations carried out in support of this paper and from available recombination data within existing atomic databases. If a recommendation was possible, data were compiled, evaluated and fitted to a functional form with associated uncertainty information retained, where available. This paper also considers the variation of the W fractional abundance due to the underlying atomic data when employing different data sets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Doubling fusion power with volumetric optimization in magnetic confinement fusion devices

A technique, volumetric power optimization, is presented for enhancing the power output of magnetic confinement fusion devices. Applied to a tokamak, this approach involves shifting the burning plasma region to a larger plasma volume while introducing minimal perturbations to the plasma boundary shape. This edge perturbation—squareness—is analogous to pinching and stretching the edge boundary. Stability calculations confirm that this edge alteration is compatible with maintaining plasma stability. This optimization method for optimizing fusion power output could improve the performance of magnetic confinement fusion power plants.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Fast piezoelectric valve offering controlled gas injection in magnetically confined fusion plasmas for diagnostic and fuelling purposes

In magnetically confined fusion plasmas controlled gas injection is crucial for plasma fuelling as well as for various diagnostic applications such as active spectroscopy. In this paper, we present a new, versatile system for the injection of collimated thermal gas beams into a vacuum chamber. This system consists of a gas pressure chamber, sealed by a custom made piezo valve towards a small capillary for gas injection. The setup can directly be placed inside of the vacuum chamber of fusion devices as it is small and immune against high magnetic fields. This enables gas injection close to the plasma periphery with high duty cycles and fast switch on/off times ≲ 0.5 ms. In this work, we present the design details of this new injection system and a systematic characterization of the beam properties as well as the gas flowrates which can be accomplished. The thin and relatively short capillary yields a small divergence of the injected beam with a half opening angle of 20°. The gas box is designed for pre-fill pressures of 10 mbar up to 100 bars and makes a flowrate accessible from 10 18 part/s up to 10 23 part/s. It hence is a versatile system for both diagnostic as well as fuelling applications. Finally, the implementation of this system in ASDEX Upgrade will be described and its application for line ratio spectroscopy on helium will be demonstrated on a selected example.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Fiber-optic silicon Fabry–Perot interferometric bolometer with improved detection limit for magnetic confinement fusion

Fiber-optic bolometers (FOBs) intended for plasma radiation measurement in magnetically confined fusion have been previously developed using a silicon pillar that functions as both a Fabry–Perot interferometer (FPI) for temperature measurement and an absorber for the radiation. We report an FOB design that can significantly improve the detection sensitivity over earlier designs by engineering the absorber of the FOB. Our design uses the fact that, compared with a silicon pillar, a gold film with the same x-ray absorption thickness will show a much higher temperature rise from a given power density of the radiation. Therefore, the responsivity of an FOB can be improved by attaching a large gold disk to the silicon FPI as the absorber. We have developed a fabrication method for FOBs of such design and obtained an FOB with a 4-µm-thick, 0.6-mm-diameter gold disk attached to a 200-µm-diameter, 100-µm-thick silicon FPI. We have characterized the noise, responsivity, response time, and noise-equivalent power density (NEPD) and compared these with the earlier design where the absorber is mainly the silicon FPI itself. The experimental result suggests that the FOB with the gold disk achieves a responsivity of ∼2.8 mK/(W/m2) and a noise-equivalent-power-density of 0.13 W/m2, which are, respectively, more than nine times larger and six times smaller compared to the FOB using a previous design. Improved NEPD and good absorption over a broad frequency range will make the FOB more attractive for applications in magnetic-confinement fusion devices.

Uddin, Nezam↗
DOE OSTI2023

Principal Challenges in Toroidal Magnetic Confinement Fusion Systems

Research efforts toward a practical and economical magnetically confined fusion power source have made steady progress over the years. Today, the principal approaches to confining a hot plasma of Deuterium and Tritium ions long enough for a significant frac/on to undergo thermonuclear burn employ toroidal geometries. Such geometries are attractive because charged particles are tightly bound around the field lines, undergoing locally helical orbits as they travel along those lines. In the core regions of the plasma, the field lines have no connection to the material surfaces making up the interior of the vacuum vessel, so that plasma confinement times can be long enough for significant thermonuclear burnup via fusion reactions to take place. For good confinement, the geometry of the field lines must be such that nested field lines (which lie on nested flux surfaces) twist around the toroidal axis (a quasi-circular line that approximately lies in the center of the confined plasma) at rates that vary based on their distance from that axis. The rotational transform (or field line pitch), iota/2π , is defined as the number of poloidal transits per single toroidal transit of a field line on a toroidal flux surface.

42 ENGINEERING↗
DOE OSTI2023

Advancement of hybrid fluid-kinetic modeling for HEDP and ICF science

We report on the development progress of a hybrid fluid-kinetic code for simulating fluids and plasmas in a wide range of environments, such as laser–matter interactions, inertial confinement fusion, magnetic confinement fusion, and pulsed power. The suite of numerical tools under development utilizes heterogeneous computer architectures and leverages the benefits of particle–based simulation techniques. By working to combine the kinetic particle-in-cell (PIC) model with a particle-based fluid simulation technique, such as smoothed particle hydrodynamics, we are developing a flexible framework capable of accurately modeling complex flows within and between kinetic and fluid regimes. The TriForce code is under development as a C++ framework for parallel, 3D, particle-based, hybrid fluid-kinetic plasma simulations. The fluid half of TriForce will be based upon the meshless smoothed-particle-hydrodynamics (SPH) approach, well-suited for shear, mixing, and turbulence, whereas the kinetic half resembles a traditional particle-in-cell (PIC) code; other particle-based approaches to fluid modeling that do use a mesh are also possible to use and are under investigation. Maxwell’s electromagnetic field equations are solved either via explicit or implicit algorithms or approximated via resistive magnetohydrodynamics (MHD) using an Ohm’s law and resulting induction equation (extended MHD is under development). A primary goal of enabling direct comparisons, from the same code, between results from the variants of MHD and implicit electromagnetic solutions is to improve our fundamental understanding of systems with magnetic fields. The code is under development to recover results from both radiation-MHD and fully kinetic codes in those limits, and is continuing to be developed from other follow-on grants to operate in between where both descriptions may co-exist and interact. For certain applications, it is desired for a simulation to contain fluid ions and electrons as well as kinetic ions and electrons. Typically, it is too computationally intensive to model a full-scale ICF or HEDP experiment fully kinetically since many cycles are expended with very small time steps on modeling the fluid part of a material that is well treated by the fluid approximation. In this case, many traditional PIC particles can be replaced with a single fluid particle representing the thermal part of the distribution function, and there are fewer needed kinetic particles, which describe the non-thermal part and can be sub-cycled relative to the fluid particle advance. Furthermore, a pure fluid code may, depending on the problem, simply lack many physically important details that are beyond the scope of its reduced approximations and assumptions. In this report, we summarize the objectives achieved in the development of the collisional and kinetic half of the code, and the physics problems to which the code has been applied in the areas of advanced and innovative fusion concepts, pulsed power, and magneto-inertial fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

On numerical errors to the fields surrounding a relativistically moving particle in PIC codes

The particle-in-cell (PIC) method is widely used to model the self-consistent interaction between discrete particles and electromagnetic fields. It has been successfully applied to problems across plasma physics including plasma based acceleration, inertial confinement fusion, magnetically confined fusion, space physics, astrophysics, high energy density plasmas. In many cases the physics involves how relativistic particles (those with high relativistic γ factors) are generated and interact with plasmas. However, when relativistic particles stream across the grid, both in “vacuum” and in plasma, many numerical issues may arise which can lead to unphysical results. We present a detailed analysis of how discretized Maxwell solvers used in PIC codes can lead to numerical errors to the fields that surround particles that move at relativistic speeds across the grid. Expressions for the axial electric field as integrals in k space are presented that reveal two types of errors. The first arises from errors to the numerator of the integrand and leads to unphysical fields that are antisymmetric about the particle. Furthermore, the second arises from errors to the denominator of the integrand and lead to Cherenkov like radiation in “vacuum”. These fields are not anti-symmetric, extend behind the particle, and cause the particle to accelerate or decelerate depending on the solver and parameters. The unphysical fields are studied in detail for two representative solvers - the Yee solver and the FFT based solver. Although the Cherenkov fields are absent, the space charge fields are still present in the fundamental Brillouin zone for the FFT based solvers. In addition, the Cherenkov fields are present in higher order zones for the FFT based solvers. Comparison between the analytical solutions and PIC simulation results are presented. A solution for eliminating these unphysical fields by modifying the k operator in the axial direction is also presented. Using a customized finite difference solver, this solution was successfully implemented into OSIRIS [1]. Results from the customized solver are also presented. Additionally, this solution will be useful for a beam of particles that all move in one direction with a small angular divergence.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

On the synchrotron emission in kinetic simulations of runaway electrons in magnetic confinement fusion plasmas

Developing avoidance or mitigation strategies of runaway electrons (REs) in magnetic confinement fusion (MCF) plasmas is of crucial importance for the safe operation of ITER. In order to develop these strategies, an accurate diagnostic capability that allows good estimates of the RE distribution function in these plasmas is needed. Synchrotron radiation (SR) of RE in MCF, besides of being one of the main damping mechanisms for RE in the high energy relativistic regime, is routinely used in current MCF experiments to infer the parameters of RE energy and pitch angle distribution functions. In the present paper we address the long standing question about what are the relationships between different REs distribution functions and their corresponding synchrotron emission simultaneously including: full-orbit effects, information of the spectral and angular distribution of SR of each electron, and basic geometric optics of a camera. We study the spatial distribution of the SR on the poloidal plane, and the statistical properties of the expected value of the synchrotron spectra of REs. Here we observe a strong dependence of the synchrotron emission measured by the camera on the pitch angle distribution of runaways, namely we find that crescent shapes of the spatial distribution of the SR as measured by the camera relate to RE distributions with small pitch angles, while ellipse shapes relate to distributions of runaways with larger the pitch angles. A weak dependence of the synchrotron emission measured by the camera with the RE energy, value of the q-profile at the edge, and the chosen range of wavelengths is observed. Furthermore, we find that oversimplifying the angular dependence of the SR changes the shape of the synchrotron spectra, and overestimates its amplitude by approximately 20 times for avalanching runaways and by approximately 60 times for mono-energetic distributions of runaways.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

An Accurate SUPG-stabilized Continuous Galerkin Discretization for Anisotropic Heat Flux in Magnetic Confinement Fusion

We present a novel spatial discretization for the anisotropic heat conduction equation, aimed at improved accuracy at the high levels of anisotropy seen in a magnetized plasma, for example, for magnetic confinement fusion. The new discretization is based on a mixed formulation, introducing a form of the directional derivative along the magnetic field as an auxiliary variable and discretizing both the temperature and auxiliary fields in a continuous Galerkin (CG) space. Both the temperature and auxiliary variable equations are stabilized using the streamline upwind Petrov–Galerkin (SUPG) method, ensuring a better representation of the directional derivatives and therefore an overall more accurate solution. This approach can be seen as the CG-based version of our previous work (Wimmer, Southworth, Gregory, Tang, 2024), where we considered a mixed discontinuous Galerkin (DG) spatial discretization including DG-upwind stabilization. We prove consistency of the novel discretization, and demonstrate its improved accuracy over existing CG-based methods in test cases relevant to magnetic confinement fusion. This includes a long-run tokamak equilibrium sustainment scenario, demonstrating a 35% and 32% spurious heat loss for existing primal and mixed CG-based formulations versus 4% for our novel SUPG-stabilized discretization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Geodesic least squares regression for scaling studies in magnetic confinement fusion

In regression analyses for deriving scaling laws that occur in various scientific disciplines, usually standard regression methods have been applied, of which ordinary least squares (OLS) is the most popular. However, concerns have been raised with respect to several assumptions underlying OLS in its application to scaling laws. We here discuss a new regression method that is robust in the presence of significant uncertainty on both the data and the regression model. The method, which we call geodesic least squares regression (GLS), is based on minimization of the Rao geodesic distance on a probabilistic manifold. We demonstrate the superiority of the method using synthetic data and we present an application to the scaling law for the power threshold for the transition to the high confinement regime in magnetic confinement fusion devices.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

First measurement in a magnetic confinement fusion experiment of the H 3 + H 3 → He 5 + n intermediate two-body resonant reaction

We report on the first experimental measurements made at a magnetic confinement fusion device of the tritium(T)-tritium(T) reaction T + T → He 4 + 2 n indicating the presence of the intermediate two-body resonant reaction T + T → He 5 + n . During the second deuterium-tritium campaign (DTE2) at the Joint European Torus, measurements of fusion plasmas with high tritium concentrations, n T / ( n T + n D ) ≈ 0.99 , heated with tritium neutral beam injection, were performed using the neutron time-of-flight (TOF) spectrometer TOFOR. We detect a peak in the neutron emission TOF spectrum consistent with the two-body resonant reaction. The TT neutron emission energy spectrum is modeled using an R -matrix framework where the distributions of the most likely model parameters given our experimental TOF data are determined utilizing a Markov chain Monte Carlo approach. We compare our best estimate of the T + T neutron emission energy spectrum with results obtained at inertial confinement fusion experiments at the OMEGA facility and find a spectral shape that is consistent with the energy dependency in the neutron spectrum observed at OMEGA. Published by the American Physical Society 2024

Physics↗
DOE OSTI2016

Multi-energy SXR cameras for magnetically confined fusion plasmas (invited)

A compact multi-energy soft x-ray camera has been developed for time, energy and space-resolved measurements of the soft-x-ray emissivity in magnetically confined fusion plasmas. Multi-energy soft x-ray imaging provides a unique opportunity for measuring, simultaneously, a variety of important plasma properties (T e , n Z , ΔZ eff , and n e,fast ). The electron temperature can be obtained by modeling the slope of the continuum radiation from ratios of the available brightness and inverted radial emissivity profiles over multiple energy ranges. Impurity density measurements are also possible using the line-emission from medium- to high-Z impurities to separate the background as well as transient levels of metal contributions. As a result, this technique should be explored also as a burning plasma diagnostic in-view of its simplicity and robustness.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Scalable Multiphysics Block Preconditioning for Low Mach Number Compressible Resistive MHD with Application to Magnetic Confinement Fusion

This study investigates multiphysics block preconditioners that are critical in devising scalable Newton–Krylov iterative solvers for longer time-scale fully implicit fluid plasma models. The specific model of interest is the visco-resistive, low Mach number, compressible magnetohydrodynamics (MHD) model. This model describes the dynamics of conducting fluids in the presence of electromagnetic fields and can be used to study aspects of astrophysical phenomena, important science and technology applications, and basic plasma physics. The specific application of interest that motivates this study is the macroscopic simulation of longer time-scale stability and disruptions of magnetic confinement fusion devices, specifically the ITER Tokamak. The computational solution of the governing balance equations for mass, momentum, heat transfer, and magnetic induction for resistive MHD systems can be extremely challenging. These difficulties arise from both the strong nonlinear, nonsymmetric coupling of fluid and electromagnetic phenomena as well as the significant range of time and length scales that the interactions of these physical mechanisms produce. To handle the range of time and spatial scales of interest, a fully implicit unstructured variational multiscale finite element formulation is employed. For the scalable solution of the Newton linearized systems, fully coupled block preconditioners are designed to leverage algebraic multigrid subsolves. In conclusion, results are presented for the strong and weak scaling of the method as well as the robustness of these techniques for a large range of Lundquist numbers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Robust analysis of space-, time-, and energy-resolved soft x-ray measurements of magnetically confined fusion plasmas (invited)

A novel compact multi-energy soft x-ray (ME-SXR) diagnostic based on the PILATUS3 100K x-ray detector has been developed in collaboration between the Princeton Plasma Physics Laboratory and the University of Wisconsin-Madison and tested on the Madison Symmetric Torus (MST) reversed-field pinch. This solid-state photon-counting detector consists of a two-dimensional array of ~100 000 pixels for which the lower photon absorption cutoff energy can be independently set, allowing it to be configured for a unique combination of simultaneous spatial, spectral, and temporal resolution of ~1 cm, 100 eV, and 500 Hz, respectively. The diagnostic is highly versatile and can be readily adapted to diverse plasma operating conditions and scientific needs without any required downtime. New results from improved-confinement and quasi-single helicity plasmas in the MST demonstrate how the detector can be applied to study multiple aspects of the evolution of magnetically confined fusion-grade plasmas. These include observing the evolution of thermal emissivity, characterizing the energy of mid-Z excitation lines, extracting the T e profile, and observing the evolution of non-thermal populations. A technique for integrating the ME-SXR diagnostic into an integrated data analysis framework based on Bayesian inference is also presented. As a result, this allows ME-SXR measurements to be combined with data for complementary diagnostics in order to simultaneously infer Te and nZ from all available information.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

A tutorial on radio frequency sheath physics for magnetically confined fusion devices

Radio frequency (RF) sheaths occur under a wide variety of conditions when RF waves, material surfaces and plasma coexist. RF sheaths are of special importance in describing the interaction of ion cyclotron range of frequency (ICRF) waves with the boundary plasma in tokamaks, stellarators and other magnetic confinement devices. In this article the basic physics of RF sheaths is discussed in the context of magnetic fusion research. Techniques for modelling RF sheaths, their interaction with RF wave fields and the resulting consequences are highlighted. The article is intended as a guide for the early-career ICRF researcher, but it may equally well serve to provide an overview of basic RF sheath concepts and modelling directions for any interested fusion scientist.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

First Measurement of Drift-Alfvén Wave Polarization in Magnetically Confined Fusion Plasmas

Polarization of drift-Alfvén waves, defined as the ratio of electrostatic to electromagnetic fluctuations, has remained unmeasurable in fusion plasmas for decades, despite its pivotal role in understanding wave dynamics and their impact on plasmas. We report the first measurements of drift-Alfvén wave polarization in a hot, magnetically-confined plasma. Here, the breakthrough is enabled by a novel methodology developed from gyrokinetic theory, utilizing fluctuations of electron temperature and density. Analysis of data from the DIII-D tokamak reveals that the waves above the geodesic acoustic mode frequency exhibit dominant electromagnetic polarization, whereas lower frequency waves show a mix of electromagnetic and electrostatic polarization, indicating a strong coupling between shear Alfvén waves and drift-acoustic waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Multimachine analysis of intermittent fluctuations in the scrape-off layer of magnetically confined fusion plasmas

Far scrape-off layer fluctuations, also known as blobs/filaments, are compared for the first time across six tokamaks: Alcator C-Mod, DIII-D, TCV, KSTAR, MAST, and MAST-U, under the context of stochastic modeling in Ohmic and low confinement mode scenarios. In this model, which well describes these fluctuations, measurements from Langmuir probes are modeled as a superposition of uncorrelated, two-sided exponential temporal pulses. The key parameters of this model are the mean waiting times between consecutive pulses, the pulse duration times and the mean of the pulse amplitudes. The study examines the influence of line-averaged density, plasma current, magnetic field and aspect ratio. The mean amplitudes and pulse waiting times increase with the line-averaged density across devices, while the pulse duration times remain constant. Whereas the variation in waiting times is gradual across the density scan, the mean amplitudes increase abruptly as the density profile broadens and flattens for Alcator C-Mod. The blob waiting times decrease with increasing plasma current and the toroidal magnetic field, while they increase with aspect ratio for conventional tokamaks. The duration times follow the same trends as the waiting times in these parameters but vary more slowly, resulting in the trends in the waiting times dominating the intermittency level of the fluctuations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗