Search NASA⌕ Search

Engineering topics

Turnbull, Alan D.

Publications and source records attributed to Turnbull, Alan D..

Accuracy of kinetic equilibrium reconstruction of NSTX and NSTX-U plasmas and its impact on the transport and stability analysis

An accurate magnetohydrodynamic (MHD) equilibrium reconstruction is an essential starting point for stability and transport plasma analysis. Herein this work describes an approach for obtaining kinetic equilibrium reconstructions using the OMFIT framework, which has been applied for the first time to spherical tokamak data from NSTX and NSTX-U. The EFIT equilibrium solver is integrated with experimental data analysis procedures and subsequent TRANSP transport simulations to enhance the accuracy of the reconstruction, in particular, at the edge region, by adding constraints on the total pressure and current density profiles, based on the transport code solution. The accuracy of the equilibrium reconstruction depends on the uncertainty and number of constraints, as well as the choice of basis functions to represent the pressure and current density profiles. Improved fidelity of the equilibrium reconstruction is demonstrated by reducing the variability of the magnetic axis and boundary locations from several centimeters, for reconstructions based on magnetic and experimental pressure constraints, to only several millimeters, for kinetic reconstructions based on transport code constraints, when different representations of basis functions were tested. The variability of the safety factor on axis was reduced ten times in the same sensitivity study. The accuracy of the equilibrium reconstruction and subsequent mapping of the experimental kinetic profile data have a significant impact on the trapped gyro Landau fluid and linear CGYRO turbulence simulations, which predict different spectra of unstable modes and turbulent fluxes for cases with different numbers of constraints in the equilibrium reconstruction. Conversely, the stability analysis performed using the GATO code shows plasmas that are stable to n = 1 MHD modes in both equilibria using magnetic and experimental pressure constraints as well as the transport code constrained equilibrium. However, a scan of parameters away from these conditions shows considerable deviation in the threshold of unstable modes between these reconstructions. Therefore, for reliable plasma analysis and use in turbulence and stability calculations, a high-fidelity equilibrium reconstruction with accurate kinetic constraints based on transport code solutions is necessary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Resistive wall mode stability and resonant field amplification in MAST high beta plasma

The n=1 (n is the toroidal mode number) resistive wall mode (RWM) instability and the resonant field amplification (RFA) due to a stable RWM response are numerically investigated for a MAST high-pressure plasma scenario, utilizing the MARS-F (Liu Y.Q. et al 2000 Phys. Plasmas 7 3681) and MARS-K (Liu Y.Q. et al 2008 Phys. Plasmas 15 112503) codes, for the purpose of understanding the potential RWM behavior in MAST-U plasmas. Ideal kink stability analysis identifies a target plasma, with parameters similar to that of the reference equilibrium reconstructed from a MAST high-pressure discharge, that accesses the RWM regime. The unstable n=1 RWM for the target equilibrium is subject to strong damping by the plasma toroidal flow and/or the drift kinetic effects from thermal particles. As a result, the mode is found to be stable under the experimental flow conditions, or even without flow stabilization if drift kinetic stabilization is included. The stability prediction is robust against variation of the assumed resistive wall minor radius. Active magneto-hydrodynamic (MHD) spectroscopy modeling, using the magnetic coils designed for controlling the edge localized modes (ELMs) in MAST-U as the antenna, shows strong resonant field amplification due to a stable RWM response in the target plasma. Maximal amplification, of 6-7 times larger than the applied vacuum field within the plasma, is obtained assuming a coil phasing of 120-180 degrees between the upper and lower rows of the ELM control coils. The MHD-kinetic hybrid model predicts generally higher RFA than the fluid model, in particular near the low-field and high-field sides of the plasma boundary surface. Furthermore, the MARS-F/K modeling thus shows that the RWM in MAST-U is subject to strong flow and kinetic stabilization, but can nevertheless be detected by active MHD spectroscopy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Magnetic-Confinement Fusion—Plasma Theory: Tokamak Magnetohydrodynamic Equilibrium and Stability

Magnetohydrodynamics (MHD) provides a useful model to describe the crucial plasma macroscopic equilibrium and stability behaviors in toroidal tokamak devices by considering the plasma as a conducting fluid interacting with a surrounding confining electromagnetic field. MHD is the most basic plasma model, incorporating most large-scale phenomena, including plasma equilibrium and all major instabilities. MHD equations are obtained by taking moments of the Boltzmann equations for different plasma species. They provide a set of comprehensive physics constrains to compute and optimize the equilibrium plasma shape and pressure and current profiles that are critical to its stability and performance. In the ideal case, the equations have special properties that lead to efficient numerical calculation schemes, the most important of which is the ideal MHD energy principle for linear stability against small departures from equilibrium. In a tokamak plasma, equilibrium pressure is mostly destabilizing for MHD modes, whereas equilibrium current is also often a major driving force. Plasma resistivity creates new freedom for a MHD instability to grow, but there are also cases where the plasma resistivity plays a stabilizing role. Equilibrium toroidal flow and/or flow shear can affect MHD instabilities. Principal MHD instabilities include the internal kink mode, sawtooth, fishbone, external kink, resistive wall mode, resistive interchange, tearing and neoclassical tearing modes (NTMs), locked modes, toroidal Alfven eigenmodes (TAEs), and edge localized modes (ELMs). Fast-growing MHD instabilities can lead to an abrupt plasma disruption and termination that can potentially damage the device plasma facing components (PFCs) and in-vessel structures. Furthermore, an important MHD application is to develop robust techniques to mitigate and control MHD instabilities.

Dispersive pellet injection↗

Influence of up-down asymmetry in plasma shape on RMP response

Shaping effect on the plasma response to the externally applied resonant magnetic perturbation field is investigated for both DIII-D and MAST experiments, utilizing toroidal modeling. The plasma boundary shape is systematically varied ranging from single-null (SN) to double-null (DN) configurations, while other equilibrium quantities are kept largely unchanged. The relative amplitude of the computed plasma surface displacement, between the top/bottom of the torus and the outboard mid-plane, is identified as the most reliable indicator that distinguishes the plasma response between the SN and DN configurations. Here, the underlying physics is the weakening of the edge-peeling component in the plasma response with increasing up-down symmetry of the plasma boundary shape.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Enhanced radial energy transport induced by radially Curved Alfvén Eigenmode wavefronts

Radially curved Alfvén eigenmode wavefronts as observed on the Electron Cyclotron Emission Imaging diagnostic, indicate that the modes are associated with a net radial energy flow. This energy flow is due to a radial mismatch in the location of the drive and damping of the mode. Energy transport is created by an enhanced parallel fluctuating magnetic field that is generated by the radial wavefront curvature.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Progress in disruption prevention for ITER

Key plasma physics and real-time control elements needed for robustly stable operation of high fusion power discharges in ITER have been shown in recent research worldwide. Recent analysis has identified the current density profile as the main drive for disruptive instabilities in discharges simulating ITER's baseline scenario with high and low external torque. Ongoing development of model-based profile control and active control of magnetohydrodynamic instabilities is improving the stability of multiple scenarios. Significant advances have been made toward real-time physics-based prediction of instabilities, including path-oriented analysis, active sensing, and machine learning techniques for prediction that are beginning to go beyond simple disruption mitigation trigger applications. Active intervention adds to prevention of disruptions, including forced rotation of magnetic islands to prevent wall locking, and localized heating/current drive to shrink the islands. Stable discharge rampdowns have been achieved with the fastest ITER-like scaled current ramp rates, while maintaining an X-point configuration. Moreover, these elements are being integrated into stable operating scenarios and new event-handling systems for off-normal events in order to develop the physics basis and techniques for robust control in ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗