Deep-Learning-Enabled Assessment of Magnetic Confinement in Magnetized Liner Inertial Fusion.
Abstract not provided.
SEARCH · Search NASA
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.
Abstract not provided.
Abstract not provided.
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.
As magnetic confinement fusion energy gains traction internationally to enable abundant energy production, designing components for fusion systems is a pressing challenge. During the planned lifetime of a fusion device, components evolve in extreme environments and must withstand large, repeated thermal loads and bombardment by 14 MeV neutrons, plasma ions, and neutral particles (deuterium, tritium, and helium), corrosive conditions, etc. All these physical processes take place simultaneously, interact in intricate ways, and impose important constraints that can affect performance. Experimental data is rare and costly to obtain, making design particularly challenging. Predictive computational frameworks must be an integral part of an accelerated and cost-effective design process by modeling fusion system performance in simulated environments. To better understand component degradation and operational impacts on their performance, the Software for Advanced Large-scale Analysis of MAgnetic confinement for Numerical Design, Engineering & Research (SALAMANDER) is designed as an open-source, fully integrated, multiphysics, multiscale, NQA-1 compliant framework facilitating 3D, high-fidelity fusion system modeling. To that end, SALAMANDER is a MOOSE-based framework, and therefore leverages MOOSE upstream libraries such as PETSc and libMesh to deliver sophisticated finite element, finite volume, and nonlinear solver technology for fusion energy simulations. SALAMANDER couples MOOSE physics module capabilities—such as thermal hydraulics, heat conduction, Navier-Stokes, and thermomechanics—with tritium transport via TMAP8, neutronics via Cardinal, and nascent particle-in-cell capabilities. Direct simulation Monte Carlo methods will be used to address neutral transport near the walls. By coupling all these physics in an integrated application, SALAMANDER will enable high-fidelity modeling of irradiation levels and plasma exposure conditions of plasma facing components and their impact on heat and tritium distributions, as well as the resulting mechanical constraints experienced by the plasma facing components and performance of blanket systems. Furthermore, SALAMANDER will be particularly suited for engineering studies thanks to the stochastic tool module readily available in MOOSE, allowing for extended uncertainty quantification and risk analysis studies. It is also able to use computer-aided design (CAD) meshes to model complex geometries, which is indispensable for fusion systems. SALAMANDER therefore supports design, safety, engineering, and research projects for magnetic confinement fusion systems
We present a novel numerical method for solving the anisotropic diffusion equation in magnetic fields confined to a periodic box which is accurate and provably stable. We derive energy estimates of the solution of the continuous initial boundary value problem. A discrete formulation is presented using operator splitting in time with the summation by parts finite difference approximation of spatial derivatives for the perpendicular diffusion operator. Weak penalty procedures are derived for implementing both boundary conditions and parallel diffusion operator obtained by field line tracing. We prove that the fully-discrete approximation is unconditionally stable. Discrete energy estimates are shown to match the continuous energy estimate given the correct choice of penalty parameters. A nonlinear penalty parameter is shown to provide an effective method for tuning the parallel diffusion penalty and significantly minimises rounding errors. Several numerical experiments, using manufactured solutions, the “NIMROD benchmark” problem and a single island problem, are presented to verify numerical accuracy, convergence, and asymptotic preserving properties of the method. Finally, we present a magnetic field with chaotic regions and islands and show the contours of the anisotropic diffusion equation reproduce key features in the field.
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.
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.
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.
Plasma transport driven by turbulence ultimately determines the energy confinement performance of controlled fusion devices regardless of their confinement schemes and configurations. A large variety of plasma instabilities have been proposed for driving turbulence responsible for anomalous plasma transport beyond classical/neoclassical transport due to collisions. Although ion-scale turbulence usually dominates due to its large eddy size and saturation level, electron-scale turbulence has been recognized to be important in regions where ion-scale turbulence is suppressed (e.g., in internal transport barrier and in spherical tokamak H-mode plasmas) or is close to marginality. Electron-scale turbulence has been shown to nonlinearly interact with ion-scale turbulence, which modifies the dynamics of both and affects the resulting plasma transport, particularly when ion-scale instability is weakly driven. In this review paper, we focus on electron-scale turbulence that is believed to operate in magnetically confinement fusion devices and aim to provide a review of theoretical, numerical, and experimental developments in understanding electron-scale turbulence and its role in driving anomalous plasma turbulence. In particular, we focus on the electrostatic electron temperature gradient (ETG) mode which is the most widely recognized plasma instability underlying electron-scale turbulence observed in magnetically confined plasmas. Finally, we note that there are other less studied instabilities that might be responsible for observed electron-scale turbulence, most notably ubiquitous mode, and short-wavelength ion temperature gradient (SWITG) mode, which will be briefly touched on in this review.
Efforts are underway to magnetically confine electron–positron pair plasmas to study their unique behaviour, which is characterized by significant changes in plasma time and length scales, supported waves and unstable modes. However, use of conventional plasma diagnostics presents challenges with these low-density and annihilating matter–antimatter plasmas. To address this problem, we propose to develop techniques based on the distinct emission provided by annihilation. This emission exhibits two spatial correlations: the distance attenuation of isotropic sources and the back-to-back propagation of momentum-preserving 2 $\gamma$ annihilation. We present the results of our analysis of the $\gamma$ emission rate and the spatial profile of the annihilation in a magnetized pair plasma from direct pair collisions, from the formation and decay of positronium as well as from transport processes. In order to demonstrate the effectiveness of annihilation-based techniques, we tested them on annular $\gamma$ emission profiles produced by a $\beta ^+$ radioisotope on a rotating turntable. Direct and positronium-mediated annihilation result in overlapping volumetric $\gamma$ sources, and the 2 $\gamma$ emission from these volumetric sources can be tomographically reconstructed from coincident counts in multiple detectors. Transport processes result in localized annihilation where field lines intersect walls, limiters or internal magnets. These localized sources can be identified by the fractional $\gamma$ counts on spatially distributed detectors.
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.
Bolometers measure the plasma radiation and provide essential diagnostic information both for safe, sustainable operation of magnetic-confinement fusion (MCF) devices and for study of the relevant physics. The fusion community has been replying on resistive bolometers for plasma radiation measurement. However, their performance is greatly compromised by the strong electromagnetic interference present in MCF devices. Previously, we have demonstrated a fiber-optic bolometer that has the promise to overcome the challenges related to the strong EMI in the MCF environment. In this project, we study the performance of the fiber-optic bolometer in a simulated MCF environment and improved the sensitivity through a novel sensor design. Part of the work reported here was performed between 09/01-12/31/2017 when the PI was with University of Nebraska-Lincoln under the same project but with a different award number (DE-SC0018273).
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.
In this work, we present Aurora, an open-source package for particle transport, neutrals and radiation modeling in magnetic confinement fusion plasmas. Aurora's modern multi-language interface enables simulations of 1.5D impurity transport within high-performance computing frameworks, particularly for the inference of particle transport coefficients. A user-friendly Python library allows simple interaction with atomic rates from the Atomic Data and Atomic Structure database as well as other sources. This enables a range of radiation predictions, both for power balance and spectroscopic analysis. We discuss here the superstaging approximation for complex ions, as a way to group charge states and reduce computational cost, demonstrating its wide applicability within the Aurora forward model and beyond. Aurora also facilitates neutral particle analysis, both from experimental spectroscopic data and other simulation codes. Leveraging Aurora's capabilities to interface SOLPS-ITER results, we demonstrate that charge exchange is unlikely to affect the total radiated power from the ITER core during high performance operation. Finally, we describe the ImpRad module in the one modeling framework for integrated task framework, developed to enable experimental analysis and transport inferences on multiple devices using Aurora.
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
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.
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.
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.