Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “toroidal magnetic fields”

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

Spatially dependent modeling and simulation of runaway electron mitigation in DIII-D

New simulations with the Kinetic Orbit Runaway electron (RE) Code (KORC) show RE deconfinement losses to the wall during plasma scrape off are the primary current dissipation mechanism in DIII-D experiments with high-Z impurity injection, and not collisional slowing down. The majority of simulations also exhibit an increase in the RE beam energy due to acceleration by the induced toroidal electric field, even while the RE beam current is decreasing. In this study, KORC integrates RE orbits using the relativistic guiding center equations of motion and incorporates time-sequenced, experimental reconstructions of the magnetic and electric fields and line integrated electron density to construct spatiotemporal models of electron and partially ionized impurity transport in the companion plasma. Comparisons of experimental current evolution and KORC results demonstrate the importance of including Coulomb collisions with partially ionized impurity physics, initial RE energy, pitch angle, and spatial distributions, and spatiotemporal electron and partially ionized impurity transport. The research presented here provides an initial quantification of the efficacy of RE mitigation via injected impurities and identification of the critical role played by loss of confinement due to plasma scrape off on the inner wall as compared to the relatively slow collisional damping.

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↗

Conceptual design of ELM control coils for the TCABR tokamak

An upgrade of the Tokamak à Chauffage Alfvén Brésilien (TCABR) is being designed to make it capable of creating a well controlled environment where the impact of resonant magnetic perturbation (RMP) fields on edge localised modes can be addressed over a wide range of (i) plasma shapes, (ii) divertor configurations, (iii) RMP coil geometries and (iv) perturbed magnetic field spectra. To this end, a unique set of in-vessel RMP coils is being designed and, in this work, their conceptual design is presented. This unique set of coils is composed of three toroidal arrays of coils on the low field side and three toroidal arrays of coils on the high field side. Each of these six toroidal arrays is composed of 18 coils thus allowing for the creation of RMP fields with toroidal mode numbers n ≤ 9 and with increased control of the poloidal mode number spectrum. To study dynamical effects of RMP fields of different toroidal mode numbers, all rotating simultaneously with different velocities, each of the 108 RMP coils will be powered independently by power supplies that can provide voltages of up to 4 kV and electric currents of up to 2 kA, with frequencies varying continuously from 0Hz up to 10kHz. A set of physical criteria were used to determine the optimal coil geometry and their respective number of turns to reduce the coil currents and voltages during operation with alternate current. Further, the conceptual design was carried out using both the vacuum approach (no plasma response) and the single-fluid response approach, which accounts for the response of a linear, single-fluid, visco-resistive plasma calculated using the M3D-C 1 code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effects of neoclassical tearing modes and toroidal field ripple on lost alpha power in the SPARC tokamak

Here, using the SPIRAL Monte Carlo, full particle-orbit simulation code, we investigate the effects of neoclassical tearing modes (NTMs) and toroidal field (TF) ripple on alpha power losses during steady-state operation of the SPARC primary reference discharge. Model perturbations for TF ripple and the $m/n = 2/1$ and $3/2$ NTMs with exaggerated widths selected based on an H-mode plasma approaching thermal quench are added to a simulated SPARC magnetic equilibrium through which marker particles are tracked. The 3/2 and 2/1 NTMs are located at $\rho_{\mathrm{pol}} \sim 0.76$ and $\rho_\mathrm{pol} \sim 0.86$ respectively, well positioned to increase alpha particle transport into and within an outer lossy region of the plasma beyond $\rho_\mathrm{pol} \sim 0.8$ where over 95% of lost alpha particles are born. Total alpha power losses are shown to increase modestly from 1.73% lost at a minimum to 2.34% lost at a maximum, and alpha particle surface power densities form localized hotspots on the first-wall near the low-field side midplane due to NTMs and TF ripple. We establish a conservative upper limit for first-wall alpha surface power densities on a toroidally symmetric wall for typical, flattop operation and motivate the consideration of NTMs in the design of three dimensional limiter surfaces for SPARC.

alpha power losses↗

Weakly quasisymmetric near-axis solutions to all orders

Here, we show that the equations satisfied by weakly quasisymmetric magnetic fields can be solved to arbitrarily high order in powers of the distance from the magnetic axis. This demonstration does not consider force balance. The existence of solutions requires an appropriate choice of parameters, most notably the toroidal current or rotational transform profiles. We do not prove that the expansion converges (it is likely divergent but asymptotic), and thus, the demonstration here should not be taken as definitive proof of the existence of global solutions. Instead, we provide a systematic construction of solutions to an arbitrarily high order.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗

Chasing the multi-modal plasma response in MAST-U

Abstract Achieving edge localized modes (ELMs) suppression in spherical tokamaks by applying resonant magnetic perturbations (RMPs) has proven challenging. The poloidal spectrum of the applied RMP is a key parameter that has an impact on the capability to mitigate and eventually suppress ELMs. In this work the resistive magnetohydrodynamic code MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) is used to evaluate the possibility of directly measuring the plasma response in MAST-U, and particularly its variation as function of the applied poloidal spectrum, in order to guide the experimental validation of the predicted best RMP configuration for ELM suppression. Toroidal mode number n = 2 RMP is considered to minimize the presence of sidebands, and to avoid the deleterious core coupling of n = 1. Singular Value Decomposition is used to highlight linearly independent structures in the simulated magnetic 3D fields and how those structures can be measured at the wall where the magnetic sensors are located. Alternative ways to measure the multimodal plasma response and how they can be used to infer the best RMP configuration to achieve ELM suppression are also presented, including the plasma displacement and the 3D footprints at the divertor plates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

First observation of RMP ELM mitigation on MAST Upgrade

Abstract The first experimental attempts at controlling edge localised modes (ELMs) via the application of resonant magnetic perturbations on the MAST Upgrade tokamak are reported. Using the linear MHD model MARS-F, the phase shift between the upper and lower coil rows ΔΦ was optimised for toroidal mode number n = 1 and n = 2 fields, to provide forward guidance to experiments. In low β N discharges, the application of n = 1 3D fields caused the ELM frequency f E L M to increase by over a factor 20 relative to the reference, and also induced a locked mode, which did not cause a plasma termination nor an H-L back transition. However when β N was raised, this induced locked mode caused plasma termination which precluded mitigation access. Initially, applying a numerically optimised n = 2 field had no effect. However applying a rigid toroidal shift to this field caused a locked mode disruption, demonstrating the presence of a substantial n = 2 error field. Coil current ramps were conducted with ΔΦ set at 6 different values, resulting in either locked mode disruptions or no effect, but mitigation with n = 2 fields was not established.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

3D radiated power analysis of JET SPI discharges using the Emis3D forward modeling tool

Abstract Precise values for radiated energy in tokamak disruption experiments are needed to validate disruption mitigation techniques for burning plasma tokamaks like ITER and SPARC. Control room analysis of radiated power ( P rad ) on JET assumes axisymmetry, since fitting 3D radiation structures with limited bolometry coverage is an under-determined problem. In mitigated disruptions, radiation is toroidally asymmetric and 3D, due to fast-growing 3D MHD modes and localized impurity sources. To address this problem, Emis3D adopts a physics motivated forward modeling (‘guess and check’) approach, comparing experimental bolometry data to synthetic data from user-defined radiation structures. Synthetic structures are observed with the Cherab modeling framework and a best fit chosen using a reduced χ 2 statistic. 2D tomographic inversion models are tested, as well as helical flux tubes and 3D MHD simulated structures from JOREK. Two nominally identical pure neon shattered pellet injection (SPI) mitigated discharges in JET are analyzed. 2D tomographic inversions with added toroidal freedom are the best fits in the thermal quench (TQ) and current quench (CQ). In the pre-TQ, 2D reconstructions are statistically the best fits, but are likely over-optimized and do not capture the 3D radiation structure seen in fast camera images. The next-best pre-TQ fits are helical structures that extend towards the high-field side, consistent with an impurity flow under the magnetic nozzle effect also observed in JOREK simulations. Whole-disruption radiated fractions of 0.98 + 0.03 / − 0.29 and 1.01 + 0.02 / − 0.17 are found, suggesting that the stored energy may have been fully mitigated by each SPI, although mitigation efficiencies well below ITER and SPARC requirements for high energy pulses are still within the large uncertainties. Emis3D is also used to validate JOREK SPI simulations, and confirms improvements in matching experiment from changes to impurity modeling. Time-dependent toroidal peaking factors are calculated and discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

MARS-F/K modeling of plasma response and fast ion losses due to RMP in KSTAR

Abstract The toroidal single-fluid magnetohydrodynamic (MHD) code MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) and the MHD-kinetic hybrid code MARS-K (Liu et al 2008 Phys. Plasmas 15 112503) are utilized to study the plasma response to the n = 1 ( n is the toroidal mode number) resonant magnetic perturbation (RMP), applied to suppress the type-I edge localized mode (ELM) in a KSTAR discharge. Both the resistive-rotating and ideal-static plasma models identify strong screening of the resonant radial field harmonics of the applied RMP due to the plasma response, and predict a strong edge-peeling response of the plasma which is consistent with the optimal ELM control coil current configuration adopted in experiment. The RMP-induced radial displacement of the plasma, computed by the resistive-rotating plasma model, agrees reasonably well with that reconstructed from the measured data in the plasma core. Taking into account the drift kinetic response of fast ions, MARS-K hybrid modeling also finds quantitative agreement of the plasma core fluid pressure perturbation with experiment. Based on the MARS-F computed plasma response, a guiding-center orbit-tracing simulation finds about 0.3% of fast ion losses due the n = 1 RMP in the KSTAR ELM control experiment considered. Most losses are associated with counter-current fast ions located near the plasma edge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

High Current Density Cables for Simpler HTS Magnets in Fusion Energy Systems

High current density cables are needed for the engineering design of potentially low cost, simpler geometry high temperature superconducting (HTS) magnets in the promising magnetic configurations as a fusion pilot plant (FPP) option. Significant technology maturation efforts are underway by privately funded startups with the goal to demonstrate mature HTS magnet technology. Test results, however, indicate critical engineering issues remain to be addressed to meet performance goals, and demonstrate HTS coil operation repeatability and reliability. To this end, exploring and enabling multiple viable conductor and cable options is vital. Partnering with a private fusion startup and manufacturers of superconducting strands and cables, Princeton Plasma Physics Laboratory (PPPL) is exploring and seeking to de-risk the aggressive high field approach presently targeted by others. Our main objective is to develop, test and calibrate novel high current density cables for a broad deployment of affordable and reliable coils using Bi-2212 conductors. If successful, such a project will provide technical feasibility for promising FPP configurations including spherical tokamaks (ST) and compact stellarators. Here, we aim at the state-of-the-art Bi-2212 cable technologies toward a current density of 100 A/mm 2 at 16 T and 4.2 K–10 K operation for low cost, simpler geometry toroidal field (TF) coils for compact stellarators developed by fusion startup companies on a timeline consistent with the FPP initiatives and beyond.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Global ion heating/transport during merging spherical tokamak formation

Here we report global ion heating/transport characteristics of magnetic reconnection during merging spherical tokamak formation experiment on TS-6 (TS-3U). Using the 96CH/320CH ultra high resolution ion Doppler tomography diagnostics, the full- 2 D imaging measurement clearly revealed that magnetic reconnection initially forms localized hot spots in the downstream region of outflow jet with inboard/outboard asymmetry (more deposition in the high field side) but the continuous accumulation of the heating coupled with transport process expands the high temperature region globally and forms characteristic poloidally ring-like structure aligned with field lines. The dynamic ion heating/transport process is also affected by the polarity of toroidal field and poloidally tilted/rotating global structure has experimentally been found both during and after merging. The characteristic poloidal asymmetry gets flipped when toroidal field direction is reversed and it was found that higher temperature appears in the positive potential side, which is opposite to the conventional understanding/prediction of guide field reconnection. Through the parallel acceleration process coupled with global heat transport, poloidally asymmetric non-classical feature has experimentally been found for the first time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A 3D helical filament surrogate model for 3D tokamak equilibria

A novel approach for efficient representation of three-dimensional (3D) tokamak equilibria is investigated, where a set of helical current filaments occupying the plasma region are employed to resolve deviations from the two-dimensional (2D) axi-symmetric state. A discrete set of 3D filaments, located at rational surfaces for a given toroidal mode number n and following the 2D equilibrium field lines (thus forming closed current loops), are found to provide a surrogate model of 3D equilibria with reasonable accuracy. Specifically, application of the filament model to 3D perturbed equilibria, due to the resonant magnetic perturbation (RMP) in DIII-D and MAST-U discharges, reveals that (1) a single helical filament per rational surface is sufficient; (2) 21 such helical filaments are capable of representing the n = 2 3D response field in MAST-U with less than 10% relative error as compared to that computed by a full magnetohydrodynamic code; (3) optimizing currents (both amplitude and phase) flowing in 3D filaments with fixed geometry, the highest accuracy fitting is found to depend on the characteristics of the 3D equilibria such as the coil current phasing of the RMP coils in our case studies. Here, whis filament approach is also applicable for generating surrogate models of other type of 3D tokamak equilibria, including those during the initial phase of the plasma disruption.

MARS-F↗

Magnetic flux emergence in binary neutron star remnants

Using high-resolution AthenaK simulations of a twisted toroidal flux tube, we study the flux emergence of magnetic structures in the shear layer of a hot massive neutron star typical of a binary neutron star remnant. High-resolution simulations demonstrate that magnetic buoyant instabilities allow for emergence only for extremely large magnetic fields significantly exceeding 10 17 G, and more typical fields around 10 16 G are instead dominated by hydrodynamic effects. Because merger remnants tend to be stable against hydrodynamic convection, our work places strong limitations on the mechanisms by which massive binary neutron star remnants can produce the magnetically-driven outflows needed to power jets.

79 ASTRONOMY AND ASTROPHYSICS↗

Initial electron cyclotron heating/current drive scoping study for the SMall Aspect Ratio Tokamak (SMART)

In this work an initial electron cyclotron heating/current drive (EC H&CD) scoping study is presented for the SMall Aspect Ratio Tokamak (SMART), constructed and operated by the Plasma Science and Fusion Technology (PSFT) Laboratory of the University of Seville. We consider two planned phases of operation with different magnitude of magnetic field (B 0.4, and 1 T). This EC H&CD scoping study is carried out by using the ray tracing code TRAVIS, which allows to quickly perform a series of scans on the main SMART plasma scenarios. A scan in the poloidal and toroidal launching angles is performed for the B 0.4 T plasma scenario assuming three launching locations (mid-plane, off-mid-plane, and top-launcher) and the extraordinary (X-) mode polarization. The absorption occurs at the second harmonic and, overall, the highest power absorption is obtained launching the EC beam from the top of the machine. A numerical analysis of the fundamental X-mode EC H&CD start-up regime for the B 1.0 T plasma scenario with two different wave frequencies (28 and 35 GHz) is presented showing a very high current drive efficiency.

Bertelli, Nicola [Princeton Plasma Physics Laborat↗

Design of radial interferometer–polarimeter for internal magnetic and density fluctuation measurements at multiple space–time scales in the National Spherical Torus Experiment-Upgrade (NSTX-U)

A Faraday-effect radial interferometer–polarimeter is designed for the National Spherical Torus Experiment-Upgrade (NSTX-U) to measure multiscale magnetic and density fluctuations critical to understanding fusion plasma confinement and stability, including those originating from magnetohydrodynamic instabilities, energetic particle-driven modes, and turbulence. The diagnostic will utilize the three-wave technique with 5 MHz bandwidth to simultaneously measure line-integrated magnetic and density fluctuations up to the ion-cyclotron frequency. Probe beams will be launched radially from the low-field side at the NSTX-U midplane, where the measured Faraday fluctuations mainly correspond to radial magnetic fluctuations that directly link to magnetic transport. A correlation technique will be employed to reduce the measurement noise to below 0.01° enabling detection of small amplitude fluctuations. Two toroidally displaced chords with 7° separation will be installed to measure toroidal mode numbers up to n = 25 for mode identification. Finally, solid-state microwave sources operating at 321 μm (935 GHz) will be used to minimize the impact of the Cotton–Mouton effect.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NSTXU Diagnostic Disruption Dynamic Loading Represented by Response Spectra

This article presents the results of transient dynamic simulations of loads due to disruption eddy currents on the NSTXU vacuum vessel. Dynamic loading at diagnostic mounting locations is expressed as response spectra derived from the time history results of the dynamic structural simulations of a variety of disruption scenarios. The disruption simulations draw on a history of the project assessments of worst case disruptions for specific components. Major efforts to assess disruption loading have included the vacuum vessel which is the major structural support for the machine, as well as the passive plates (PPs), high harmonic fast wave (HHFW) antenna, and centerstack casing. Each one of these efforts included transient electromagnetic simulations producing time-dependent eddy current Lorentz loads (and in some cases halo loads) which then were applied to time-dependent structural dynamic analyses intended to obtain the proper dynamic amplification factors. In some instances, the EM model and structural model were identical allowing direct transfer of EM forces to the structural model. In other cases, the EM and structural model were not identical and the vector potential (VP) transfer method was used. The results files from these analyses were available (or re-run) to post process in ANSYS Classic time history postprocessor. In conclusion, the ANSYS command is used to create response spectra from time history data at desired points on the vessel.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Coupled Experimental/Computational Investigation of the Dynamics of Interacting Magnetized Plasmas

The interaction, or interpenetration, of magnetized plasmas of different density and/or pressure occurs in a wide variety of natural and man-made systems. Such systems include extragalactic jets propagating into the intergalactic medium, solar coronal mass ejections into background solar wind, compact toroid (CT) fueling of magnetic fusion plasmas, and jets of capsule shell impurities into DT fusion fuel, which can lead to enhanced impurity mix in inertial fusion implosions. These plasmas may take the form of jets, with open, helical magnetic structure, or plasma “bubbles” with closed magnetic fields (B-fields), such as spheromaks or CT’s. Such structures, both open and closed B-field cases, can transport heat, particles and magnetic flux or magnetic helicity into background plasma regions. For example, the origin of extragalactic magnetic fields may be due, at least in part, to transport by astrophysical jets. The goal of this proposed work was to elucidate the detailed plasma and magnetic field dynamics of high-density plasma jets (open B-field) and bubbles (closed B-field) propagating into lower density background magnetized plasma through controlled laboratory experiments and closely coupled nonlinear MHD modeling. These experiments were conducted in the HelCat (Helicon-Cathode) linear plasma device at the University of New Mexico (UNM). Plasma jets and bubbles were launched via an existing compact coaxial plasma gun, mounted on the HelCat device. This gun produced plasmas tens of cm in scale and lasting tens of microseconds, thereby allowing detailed multipoint, space- and time-resolved measurements to be made routinely. The experiments were directly modeled using the extended magnetohydrodynamic (XMHD) PERSEUS code, developed at Cornell University [23,24]. Both experimental and numerical modeling work are ongoing. The main results to date are reported here. Additional supplemental funding for one year (8/1/2019 – 7/31/2020) supported numerical investigation of photoionization processes important in many low temperature plasmas, including the HelCat device. Initial results of this modeling work is also reported.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗