Search NASA⌕ Search

SEARCH · Search NASA

Results for “magnetohydrodynamics (MHD)”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

3D modeling of the Staged Z-pinch with the FLASH code

The Staged Z-pinch (SZP) fusion concept is a magneto-inertial compression scheme developed by Magneto-Inertial Fusion Technologies, Inc. (MIFTI), in which small amounts of fusion fuel are brought to fusion-relevant conditions by passing multi-million amperes strong current through a cylindrical shell of high atomic number material. One- and two-dimensional modeling performed by MIFTI with the MACH2 code, suggests that net fusion energy gain can be achieved when currents in the 10 million amperes range compress a 50%-50% mixture of deuterium and tritium gas. In this project we enhance and use the FLASH code to execute high-fidelity, simulations of various SZP configurations, a collaboration between the Flash Center for Computational Science at the University of Rochester and the Laboratory for Laser Energetics, and MIFTI. FLASH is a high-performance computing, multi-physics, radiation magnetohydrodynamic (MHD) code with extended physics capabilities, which is developed at the Flash Center. The goal of this project is to assess the shell/fuel stability of the pinch to two- and three-dimensional MHD instabilities and to utilize FLASH’s extended physics capabilities to understand how extended-MHD effects impact implosion dynamics and plasma conditions at stagnation. The project is a natural extension of an ongoing collaborative effort between MIFTI and the Flash Center through the U.S. DOE Advanced Research Projects Agency-Energy (ARPA-E) BETHE program and can provide MIFTI with simulation capabilities that are currently beyond their reach with MACH2. Ultimately, MIFTI wants FLASH to become one of their simulation workhorses for reliable, high-fidelity SZP platform design. During the project, the Flash Center team developed a suite of in one- and two-dimensional FLASH simulations to model different variants of the SZP platform, performed code-to-code comparisons with MACH2, and used the FLASH code to design and model SZP experiments like the Double Eagle experiments that MIFTI recently performed. Also, the Flash Center team developed the capabilities of FLASH to be able to do three-dimensional simulations of pulsed-power experiments with the code for the first time. Two publications from this effort are currently under review and two more are in preparation. MIFTI has also committed to the use of FLASH for future SZP simulation efforts.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Anisotropic Energy Transfer and Conversion in Magnetized Compressible Turbulence

We present a spatial filtering (or coarse-graining) analysis on 3D magnetized magnetohydrodynamic (MHD) turbulence simulations. The filtered compressible MHD formulae show transfer of kinetic and magnetic energies from large to small scales, as well as energy conversion between kinetic, magnetic, and thermal energies. The anisotropic filtering enables separate analyses of the energy flows perpendicular and parallel to the global mean magnetic field. Anisotropy in energy cascade is demonstrated by the larger perpendicular energy cascade rate and also the larger perpendicular wavenumbers associated with the peak energy transfer rate. We also find that the “inertial range” along the parallel (perpendicular) direction in the anisotropic energy cascade formulation is no longer strictly dissipation-free, because it includes the dissipation in the perpendicular (parallel) direction. A change in the driving force (kinetic only versus kinetic and magnetic) affects the energy conversion between kinetic and magnetic energies. While the compressibility of the driving force changes the partition of different channels of energy transfer and conversion, and also increases the total energy transfer rate, the global energy flow remains unaffected by compressibility qualitatively. Our analysis can be applied to multispacecraft observations of turbulence in the solar wind or a planetary magnetosphere.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The theory of kinetic effects on resistive wall mode stability in tokamaks

Tokamak fusion plasmas benefit from high pressures but are then susceptible to modes of instability. These magnetohydrodynamic (MHD) modes are macroscopic distortions of the plasma, but certain collective motions of individual particles can provide stabilizing effects opposing them. The presence of a resistive wall slows the mode growth, converting a kink to a resistive wall mode (RWM). A kinetic MHD model includes Maxwell's equations, ideal MHD constraints, and kinetic effects included through the pressure tensor, calculated with the perturbed drift-kinetic distribution function of the particles. The kinetic stabilizing effects on the RWM arise through resonances between the plasma rotation and particle drift motions: precession, bounce, and transit. A match between particle motions and the mode allows efficient transfer of energy that would otherwise drive the growth of the mode, thus damping the growth. The first approach to calculating RWM stability is to write a set of equations for the complex mode frequency in terms of known quantities and then to solve the system. The “energy principle” approach, which has the advantage of clarity in distinguishing the various stabilizing and destabilizing effects, is to change the force balance equation into an equation in terms of changes of kinetic and potential energies, and then to write a dispersion relation for the mode frequency in terms of those quantities. These methods have been used in various benchmarked codes to calculate kinetic effects on RWM stability. Importantly, the theory has illuminated the important roles of plasma rotation, energetic particles, and collisions in RWM stability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Grad–Shafranov equilibria via data-free physics informed neural networks

A large number of magnetohydrodynamic (MHD) equilibrium calculations are often required for uncertainty quantification, optimization, and real-time diagnostic information, making MHD equilibrium codes vital to the field of plasma physics. In this paper, we explore a method for solving the Grad–Shafranov equation by using physics-informed neural networks (PINNs). For PINNs, we optimize neural networks by directly minimizing the residual of the partial differential equation as a loss function. We show that PINNs can accurately and effectively solve the Grad–Shafranov equation with several different boundary conditions, making it more flexible than traditional solvers. This method is flexible as it does not require any mesh and basis choice, thereby streamlining the computational process. We also explore the parameter space by varying the size of the model, the learning rate, and boundary conditions to map various tradeoffs such as between reconstruction error and computational speed. Additionally, we introduce a parameterized PINN framework, expanding the input space to include variables such as pressure, aspect ratio, elongation, and triangularity in order to handle a broader range of plasma scenarios within a single network. Parameterized PINNs could be used in future work to solve inverse problems such as shape optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Three-dimensional modeling of a plasma in a strong azimuthal magnetic field

Three-dimensional magnetohydrodynamic (MHD) simulations are able to model the generation of disc-shaped plasma, driven by laser ablation from a current carrying rod in a pulsed-power machine producing azimuthal magnetic fields of 2–3 MG. The plasma at such extreme conditions is unique in that the parameter space for the plasma β and Hall parameter χ transition from below unity to greater than unity at different stages of the plasma generation. In simulations, the formation of the plasma disc in the azimuthal direction, is driven by heat flux from the laser spot and depends on the set of transport coefficients used in simulations. The most recent set of transport coefficients leads to the formation of plasma ejecta at the back end of the rod, that qualitatively matches experiments. Specifically, the cross-gradient Nernst effect, which twists magnetic field, is shown to have a large effect on the shape of the back end ejecta. In the direction along the axis of the rod there is propagation of perturbations from the disc as observed in experiments. In simulations the period of temperature perturbations is in good agreement with experimental results. Furthermore, an instability due to coupling of heat flux and the magnetic field advection provides a possible explanation for perturbation growth along the axis of the rod, and the instability growth rate is consistent with experimental results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A study of resistive peeling–ballooning modes across spherical tokamaks

We investigate how non-ideal-magnetohydrodynamics (MHD) effects, in particular plasma resistivity, impact the peeling–ballooning stability thresholds in spherical tokamaks. This analysis follows the discovery of resistive kink-peeling modes in ELMing National Spherical Torus Experiment (NSTX) discharges. In the present study we extend this modeling to ELMing pulses in the Mega Ampere Spherical Tokamak (MAST) and MAST—Upgrade (MAST-U), where we find a clear resistive scaling for peeling–ballooning modes. While in NSTX ideal-MHD predicts stability for ELMing discharges, in MAST-U we find that the plasma is slightly unstable to peeling–ballooning modes, but is fully stabilized once diamagnetic effects are considered in terms of a growth rate normalization. A resistive power law scaling is calculated for these modes on MAST-U, which lies in between that of tearing modes and resistive interchange modes. A comparison between M3D-C1 and NIMROD shows reasonable agreement for this scaling. Resistivity destabilizes the modes and the peeling–ballooning unstable domain is considerably expanded in both, MAST and MAST-U. In addition to the MAST/-U pulses we also analyze resistive PB stability in a NSTX-similarity discharge on DIII-D. While having a different aspect ratio from NSTX, this discharge uses NSTX-like shaping parameters, toroidal field and plasma current. By considering these discharges alongside NSTX cases, we identify conditions influencing the onset of resistive peeling–ballooning modes. Furthermore, our findings indicate that magnetic shear in the pedestal region is closely linked to the emergence of resistive edge modes.

MAST-U↗
DOE OSTI2025

Biermann-Battery-Driven Magnetized Collisionless Shock Precursors in Laser-Produced Plasmas

This Letter reports the first complete observation of magnetized collisionless shock precursors formed through the compression of Biermann-battery magnetic fields in laser produced plasmas. At OMEGA, lasers produce a supersonic CH plasma flow which is magnetized with Biermann-battery magnetic fields. The plasma flow collides with an unmagnetized hydrogen gas jet plasma to create a magnetized shock precursor. The situation where the flowing plasma carries the magnetic field is similar to the Venusian bow shock. Imaging 2⁢𝜔 Thomson scattering confirms that the interaction is collisionless and shows density and temperature jumps. Proton radiographs have regions of strong deflections and FLASH magnetohydrodynamic (MHD) simulations show the presence of Biermann fields in the Thomson scattering region. Electrons are accelerated to energies of up to 100 keV in a power-law spectrum. Here, OSIRIS particle-in-cell (PIC) simulations, initialized with measured parameters, show the formation of a magnetized shock precursor and corroborate the experimental observables.

High-energy-density plasmas↗
DOE OSTI2022

Turbulent Regimes in Collisions of 3D Alfvén-wave Packets

Using three-dimensional gyrofluid simulations, we revisit the problem of Alfvén-wave (AW) collisions as building blocks of the Alfvénic turbulent cascade and their interplay with magnetic reconnection at magnetohydrodynamic (MHD) scales. Depending on the large-scale value of the nonlinearity parameter χ 0 (the ratio between the AW linear propagation time and nonlinear turnover time), different regimes are observed. For strong nonlinearities (χ 0 ~ 1), turbulence is consistent with a dynamically aligned, critically balanced cascade—fluctuations exhibit a scale-dependent alignment $\sin {\theta }_{{k}_{\perp }}\propto {k}_{\perp }^{-1/4}$, resulting in a ${k}_{\perp }^{-3/2}$ spectrum and ${k}_{\parallel }\propto {k}_{\perp }^{1/2}$ spectral anisotropy. At weaker nonlinearities (small χ 0 ), a spectral break marking the transition between a large-scale weak regime and a small-scale ${k}_{\perp }^{-11/5}$ tearing-mediated range emerges, implying that dynamic alignment occurs also for weak nonlinearities. At χ 0 < 1 the alignment angle ${\theta }_{{k}_{\perp }}$ shows a stronger scale dependence than in the χ 0 ~ 1 regime, namely $\sin {\theta }_{{k}_{\perp }}\propto {k}_{\perp }^{-1/2}$ at χ 0 ~ 0.5, and $\sin {\theta }_{{k}_{\perp }}\propto {k}_{\perp }^{-1}$ at χ 0 ~ 0.1. Dynamic alignment in the weak regime also modifies the large-scale spectrum, scaling approximately as ${k}_{\perp }^{-3/2}$ for χ 0 ~ 0.5 and as ${k}_{\perp }^{-1}$ for χ 0 ~ 0.1. A phenomenological theory of dynamically aligned turbulence at weak nonlinearities that can explain these spectra and the transition to the tearing-mediated regime is provided; at small χ 0 , the strong scale dependence of the alignment angle combines with the increased lifetime of turbulent eddies to allow tearing to onset and mediate the cascade at scales that can be larger than those predicted for a critically balanced cascade by several orders of magnitude. Such a transition to tearing-mediated turbulence may even supplant the usual weak-to-strong transition.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Mallat Scattering Transformation based surrogate for Magnetohydrodynamics

Abstract A Machine and Deep Learning (MLDL) methodology is developed and applied to give a high fidelity, fast surrogate for 2D resistive MagnetoHydroDynamic (MHD) simulations of Magnetic Liner Inertial Fusion (MagLIF) implosions. The resistive MHD code is used to generate an ensemble of implosions with different liner aspect ratios, initial gas preheat temperatures (that is, different adiabats), and different liner perturbations. The liner density and magnetic field as functions of x , y , and z were generated. The Mallat Scattering Transformation (MST) is taken of the logarithm of both fields and a Principal Components Analysis (PCA) is done on the logarithm of the MST of both fields. The fields are projected onto the PCA vectors and a small number of these PCA vector components are kept. Singular Value Decompositions of the cross correlation of the input parameters to the output logarithm of the MST of the fields, and of the cross correlation of the SVD vector components to the PCA vector components are done. This allows the identification of the PCA vectors vis-a-vis the input parameters. Finally, a Multi Layer Perceptron (MLP) neural network with ReLU activation and a simple three layer encoder/decoder architecture is trained on this dataset to predict the PCA vector components of the fields as a function of time. Details of the implosion, stagnation, and the disassembly are well captured. Examination of the PCA vectors and a permutation importance analysis of the MLP show definitive evidence of an inverse turbulent cascade into a dipole emergent behavior. The orientation of the dipole is set by the initial liner perturbation. The analysis is repeated with a version of the MST which includes phase, called Wavelet Phase Harmonics (WPH). While WPH do not give the physical insight of the MST, they can and are inverted to give field configurations as a function of time, including field-to-field correlations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Free-Surface Liquid Lithium Flow Modeling and Stability Analysis for Fusion Applications

Liquid metal plasma facing components are considered an attractive design choice for fusion devices including pilot plants. Virtual prototyping of such devices includes modeling of free-surface flow of the electrically conductive liquid, which requires computational fluid dynamics (CFD) and magnetohydrodynamics (MHD) simulations. Numerical tools capable of simulating flows and heat transfer in the free-surface MHD flow were developed at PPPL based on the customized ANSYS CFX. Here, MHD is introduced using a magnetic vector potential approach. Free-surface flow capabilities are available in the code and were tested. Special stabilization procedures were derived and applied to improve convergence of the momentum equations with the source terms due to the Lorentz force and surface tension. Important characteristics of the fusion-relevant liquid metal flow is free surface smoothness and stability. Heat flux from the plasma impacts the liquid surface at a very acute angle, so any change of the free surface from axisymmetry can dramatically increase the local heat flux density and thus create excessive evaporation of liquid lithium into the plasma, which is detrimental to operations. Stability analysis of the liquid metal film flow was performed to determine applicable flow regimes. Thin film flow along horizontal wall is considered. Effects of gravity, magnetic field, and surface tension are included in the analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Simulation Models for Exploring Magnetic Reconnection

Simulations have played a critical role in the advancement of our knowledge of magnetic reconnection. However, due to the inherently multiscale nature of reconnection, it is impossible to simulate all physics at all scales. For this reason, a wide range of simulation methods have been crafted to study particular aspects and consequences of magnetic reconnection. This article reviews many of these methods, laying out critical assumptions, numerical techniques, and giving examples of scientific results. Plasma models described include magnetohydrodynamics (MHD), Hall MHD, Hybrid, kinetic particle-in-cell (PIC), kinetic Vlasov, Fluid models with embedded PIC, Fluid models with direct feedback from energetic populations, and the Rice Convection Model (RCM).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A Large Eddy Simulation Study of Flow Turbulence, Alumina Transport, and Bath Temperature Evolution in Conventional Aluminum-Smelting Cell Using OpenFOAM

In this study, a Large Eddy Simulation (LES) of the aluminum-smelting process is performed using OpenFOAM. To understand the coupled behavior of heat transfer, mass transfer, and flow of the smelting process, a multi-physics computational fluid dynamics (CFD) model based on the Eulerian–Eulerian multi-fluid approach is adopted. The model accounts for CO 2 bubble and magnetohydrodynamics (MHD)-driven flow, along with alumina dissolution, transport, and bath temperature evolution. The simulation predictions show small-scale turbulent vortical structures in the anode–cathode space caused by combined effect of MHD and CO 2 bubble-bath interactions and relatively large-scale asymmetric vortices in the inter-anode space caused by the CO 2 bubble-bath interactions. The vortex formation at the edges of the anodes evidently aids in transporting alumina from the central channel to the bottom of the anodes and prevents accumulation of gas bubbles in the periphery of the anode bottom. Symmetric bath cold spots are observed in the vicinity of the feeder. Cold spots are also observed in the anode–cathode distance space below the anode bottom due to the transport of undissolved solid to this region by the flow. The findings from the work are useful in developing and designing alumina-feeding strategy leading to reduced anode effects and smooth operation of the cell. The work also highlights the important flow structures in conventional aluminum-smelting cell.

36 MATERIALS SCIENCE↗
DOE OSTI2024

TokaMaker: An open-source time-dependent Grad-Shafranov tool for the design and modeling of axisymmetric fusion devices

Here, in this paper, we present a new static and time-dependent MagnetoHydroDynamic (MHD) equilibrium code, TokaMaker, for axisymmetric configurations of magnetized plasmas, based on the well-known Grad-Shafranov equation. This code utilizes finite element methods on an unstructured triangular grid to enable capturing accurate machine geometry and simple mesh generation from engineering-like descriptions of present and future devices. The new code is designed for ease of use without sacrificing capability and speed through a combination of Python, Fortran, and C/C++ components. A detailed description of the numerical methods of the code, including a novel formulation of the boundary conditions for free-boundary equilibria, and validation of the implementation of those methods using both analytic test cases and cross-code validation is shown. Results show expected convergence across tested polynomial degree for analytic and cross-code test cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

MFC 5.0: An exascale many-physics flow solver

Many problems of interest in engineering, medicine, and the fundamental sciences rely on high-fidelity flow simulation, making performant computational fluid dynamics solvers a mainstay of the open-source software community. Previous work MFC 3.0 was made a published, documented, and open-source solver via Bryngelson et al. Comp. Phys. Comm. (2021) with numerous physical features, numerical methods, and scalable infrastructure. MFC 5.0 is a significant update to MFC 3.0, featuring a broad set of well-established and novel physical models and numerical methods, as well as the introduction of GPU and APU (or superchip) acceleration. Here, we exhibit state-of-the-art performance and ideal scaling on the first two exascale supercomputers, OLCF Frontier and LLNL El Capitan. Combined with MFC’s single-accelerator performance, MFC achieves exascale computation in practice, and achieved the largest-to-date public CFD simulation at 200 trillion grid points as a 2025 ACM Gordon Bell Prize finalist. New physical features include the immersed boundary method, N-fluid phase change, Euler–Euler and Euler–Lagrange sub-grid bubble models, fluid-structure interaction, hypo- and hyper-elastic materials, chemically reacting flow, two-material surface tension, magnetohydrodynamics (MHD), and more. Numerical techniques now represent the current state-of-the-art, including general relaxation characteristic boundary conditions, WENO variants, Strang splitting for stiff sub-grid flow features, and low Mach number treatments. Weak scaling to tens of thousands of GPUs on OLCF Summit and Frontier and LLNL El Capitan achieves efficiencies within 5% of ideal to over 90% of their respective system sizes. Strong scaling results for a 16-times increase in device count show parallel efficiencies over 90% on OLCF Frontier. MFC’s software stack has undergone further improvements, including continuous integration, which ensures code resilience and correctness through over 300 regression tests; metaprogramming, which reduces code length while maintaining performance portability; and code generation for computing chemical reactions

Computational fluid dynamics↗
DOE OSTI2026

Fast solvers for tokamak fluid models with PETSc

Multigrid (MG) is widely recognized as a highly effective solver for the model problem, the Laplacian, but textbook MG fails on most problems of interest. MG methods have been applied to complex, real-world applications with careful consideration of the physical model and discretization. In this work we develop the first step in applying MG methods to science and engineering relevant magnetohydrodynamics (MHD) tokamak models in the M3D-C1 (https://m3dc1.pppl.gov) fusion energy science code. The semi-implicit time integrator in M3D-C1 is composed of many linear solves. The implicit advance of the momentum equation is the most challenging and is the focus of this work. The current production solver in M3D-C1 is a block Jacobi (BJ) preconditioner within a Krylov solver, where blocks group degrees of freedom on planes of constant toroidal coordinate. BJ convergence degrades as the number of planes increases due to the spectral properties of the matrix preconditioned with BJ. The partially magnetic field-aligned, regular toroidal grid structure in M3D-C1 is amenable to semi-coarsening geometric MG in the toroidal direction. This paper develops such a solver and demonstrates competitive performance on a runaway electron model of a SPARC (https://cfs.energy/technology/sparc) disruption, and superior robustness on a stellarator model on which the BJ solver fails to converge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Multi-device study of temporal characteristics of magnetohydrodynamic modes initiating disruptions

Disruptions in tokamaks are often preceded by magnetohydrodynamic (MHD) instabilities that can rotate or become locked to the wall. Measurements from magnetic diagnostics in the presence of MHD mode precursors to disruptions can yield potentially valuable input to the plasma control system, with a view to disruption avoidance, prediction and mitigation. This report presents an exploratory analysis of the growth of MHD modes and corresponding time scales on the basis of magnetic measurements in multiple tokamaks. To this end, a database was compiled using disruptive discharges from COMPASS, ASDEX Upgrade, DIII-D and JET, manually classified according to disruption root cause, and characterized by a great diversity of operational conditions and mode dynamics. The typical time during which a mode can be detected using saddle coils and the duration of the locked mode phase in the database both extend over several orders of magnitude, but generally the time scales increase with plasma size. Several additional factors are discussed that can influence these durations, including the disruption root cause. A scaling law for the locked phase duration was estimated, yielding predictions toward ITER of the order of hundreds of milliseconds or even seconds. In addition, a scaling law for the mode amplitude at the disruption onset, proposed earlier by de Vries et al. (2016), is applied to the database, and its predictive capabilities are assessed. Despite significant uncertainty on the predictions from both scaling laws, encouraging trends are observed of the fraction of disruptions that may be detected with sufficient warning time to allow mitigation or even avoidance, based solely on observations of MHD mode dynamics. When combined with similar analysis of measurements from diagnostics that are sensitive to other disruption precursors, our analysis methods and results may contribute to the reliability, robustness and generalization of disruption warning schemes for ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Nonlinear control of safety factor gradient in tokamaks using spatially variable electron cyclotron current drives

Active control of plasma properties may be necessary to achieve stable operation of next-generation tokamaks over large time scales. Such control algorithms can regulate the plasma properties to avoid the onset of magnetohydrodynamic (MHD) instabilities. For instance, the global and local properties of the safety factor profile are linked to the onset of neoclassical tearing modes (NTMs). This work proposes a model-based control approach for deterring/delaying NTMs through active modulation of one of the safety factor properties - the gradient of the safety factor profile at a particular rational safety factor surface. In particular, a novel control-oriented model for the local safety factor gradient is developed. The nonlinear control model is governed by a nonautonomous ordinary differential equation that accounts for a given rational safety factor surface’s spatial variation over time. Further, to improve the controllability of the spatially evolving parameter, the control model treats ECH&CD positions, along with noninductive powers, as controllable variables. A nonlinear control algorithm based on feedback linearization with optimization is synthesized to achieve the objective of regulating the safety factor gradient around a given target. The proposed algorithm allocates optimal ECH&CD positions, in addition to auxiliary powers, at each time instant as the rational safety factor surface drifts to locations with a low control authority. Stability guarantees of the proposed control law are also discussed in this work. The proposed algorithm is tested for a DIII-D tokamak scenario in nonlinear simulations carried out using the Control Oriented Transport SIMulator (COTSIM). Both fixed and moving ECH&CD cases are studied, and their outcomes are compared. Simulation results demonstrate that enthusiastic regulation of the safety factor gradient can be achieved during the ramp-up and flat-top phases of tokamak operation in both fixed and moving ECH&CD cases. However, real-time updates of ECH&CD positions can prevent the saturation of auxiliary powers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Nonlinear control of the minimum safety factor in tokamaks by optimal allocation of spatially moving electron cyclotron current drive

The minimum value of the safety factor profile is related to the magnetohydrodynamic (MHD) stability of the plasma confined in a tokamak. Therefore, active control of the minimum safety factor may mitigate MHD instabilities that can degrade or even terminate plasma confinement. Typically, in most tokamak scenarios, the minimum safety factor evolves spatially with time, i.e., the location at which the safety factor achieves the minimum value changes with time. In addition to the inherent nonlinearities in the minimum safety factor evolution, its spatial variation makes the control design challenging. In particular, complexity in control design may arise from the need for time-dependent nonlinear models that account for spatial variation of the minimum safety factor. Furthermore, the minimum safety factor may drift to locations where the actuator authority is low. The problem of minimum safety factor control with target location tracking and moving electron cyclotron current drive (ECCD) is addressed in this work. A nonlinear time-dependent model that incorporates the spatial variation of the minimum safety factor is presented. A nonlinear controller based on optimal feedback linearization is developed to track a target minimum safety factor. The proposed controller treats the ECCD position as a controllable variable. In other words, the controller prescribes the ECCD position (in addition to the non-inductive powers) in real time based on an optimal criterion that is defined a priori. This work also presents the steps necessary to integrate the minimum safety factor controller with a total energy controller to achieve multiple control objectives simultaneously. In conclusion, the proposed integrated control algorithm is tested using nonlinear simulations in the Control Oriented Transport SIMulator (COTSIM) for a DIII-D tokamak scenario.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗