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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Effect of resonant magnetic perturbations including toroidal sidebands on magnetic footprints and fast ion losses in HL-2M

Externally applied resonant magnetic perturbations (RMPs), generated by magnetic coils located outside the plasma (referred to as RMP coils), provide an effective way to control the edge localized mode (ELM) in tokamak devices. Due to the discrete nature of the toroidal distribution of these window-frame coils, toroidal sidebands always exist together with the fundamental harmonics designed for ELM control. In this work, the MARS-F code (Liu et al 2000 Phys. Plasmas7 3681) is applied to investigate the detailed features of the RMP spectra considering both the dominant harmonic (n = 2) and the associated sideband (n = 6), and the impact of the combined fields on magnetic footprints as well as on the fast ion losses for a reference double-null scenario in the HL-2M device. It is found that the sum of the n = 2 and n = 6 RMP fields splits the footprint and widens the footprint area, as compared to the single-n (n = 2) harmonic case. The resistive plasma response breaks the up–down symmetry of the footprint pattern on the outer divertor plates, which is otherwise symmetric assuming vacuum RMP fields. Considering fast ion losses, a threshold value exists for the initially launched radial position of test particles, as well as for the RMP coil current, before the loss occurs. When the threshold criterion is satisfied, the combined n = 2 and n = 6 RMP fields enhance the fast ion loss rate by , as compared to that of the n = 2 component alone. These results illustrate the important role of the sideband of RMP fields on the magnetic footprints and fast ion losses in tokamak plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

MARS-Q modeling of kink-peeling instabilities in DIII-D QH-mode plasma

In quiescent H-mode (QH-mode) regime, edge harmonic oscillations (EHOs) are believed to provide necessary radial transport to prevent occurrence of large edge localized modes. A systematic modeling study is performed here on the low- n EHOs in a DIII-D QH-mode plasma (Chen et al 2016 Nucl. Fusion 56 076011), by utilizing the MARS-Q code (Liu et al 2013 Phys. Plasmas 20 042503). Both the n = 1 and n = 2 instabilities are found to be strongly localized near the plasma edge, exhibiting the edge-peeling characteristics. The DIII-D resistive wall is found to have minor effects on these instabilities. The plasma resistivity is found to strongly modify the mode growth rate. Assuming the Spitzer model for the plasma resistivity, the computed mode growth rate scales as S sup>-1/3 with S being the Lundquist number. Toroidal flow of the plasma slightly stabilizes these edge localized kink-peeling modes. Drift kinetic effects all have a destabilization effect on these modes. Non-perturbative magneto-hydrodynamic-kinetic hybrid computations find that the drift kinetic effects associated with thermal particle species push the peak location of the eigenmode radially inward but still in the pedestal region. The modeled plasma temperature and density fluctuations in the plasma edge region, as well as the poloidal magnetic field perturbations along both the low and high field sides of the plasma surface, are in good agreement with experimental measurements. Finally, the quasi-linear initial value simulations find a strong non-linear interplay between the kink-peeling instability and the toroidal flow near the plasma edge. The combined effect of the damping of the flow amplitude and change of the edge flow shear is found to be the stabilizing factor for the kink-peeling mode, leading to the mode saturation and thus EHOs.

Physics↗
DOE OSTI2020

New heat flux model for non-axisymmetric divertor infrared structures

A convective heat flux model for perturbed plasmas, based on guiding center ion drift in vacuum elds (A. Wingen, et al., Phys. Plasmas 21 (2014) 012509), has been updated. The old model only considered ion heat flux, while here also electron heat flux is included. The updated model predicts divertor heat flux distributions in non-axisymmetric (3D) plasmas with applied Resonant Magnetic Perturbation (RMP) fields, and includes electric scalar potentials. It is found that a radial electric field in the near Scrape-o Layer (SOL) can considerably shift the footprints toroidally, leading to a smearing out eff ect of the incident heat flux, while a simple model for sheath potential has little impact on footprints. Various approaches to model electron heat flux are studied. A convective electron model, based on collisionless free streaming, is found to yield the best agreement with measurements, while a conductive model requires a at temperature gradient inside lobes to yield acceptable peak heat ux values. A heuristic heat flux layer approach, based on a fixed layer width also requires a limited heat flux inside the last closed flux surface (LCFS); by selecting various locations of the LCFS, the results of the conductive or convective model can be recovered respectively. The sum of ion and electron heat fluxes, both obtained by the convective model, is compared to experimental data for multiple time slices in DIII-D. Strike point splitting is observed with peak heat fluxes and layer widths that compare well to infrared camera (IR) measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Neoclassical toroidal viscosity torque prediction via deep learning

GPECnet is a densely connected neural network that has been trained on GPEC data, to predict the plasma stability, neoclassical toroidal viscosity (NTV) torque, and optimized 3D coil current distributions for desired NTV torque profiles. Using NTV torque, driven by non-axisymmetric field perturbations in a tokamak, can be vital in optimizing pedestal performance by controlling the rotation profile in both the core, to ensure tearing stability, and the edge, to avoid edge localized modes (ELMs). The generalized perturbed equilibrium code (GPEC) software package can be used to calculate the plasma stability to 3D perturbations and the NTV torque profile generated by applied 3D magnetic fields. These calculations, however, involve complex integrations over space and energy distributions, which takes time to compute. Initially, GPECnet has been trained solely on data representative of the quiescent H-mode (QH) scenario, in which neutral beams are often balanced and toroidal rotation is low across the plasma profile. Lastly, this work provides the foundation for active control of the rotation shear using a combination of beams and 3D fields for robust and high performance QH mode operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A magnetic diagnostic suite for the Pegasus-III experiment

Pegasus-III is an ultralow aspect ratio spherical tokamak providing a dedicated US experiment for comparative solenoid-free startup studies. A new magnetic diagnostic suite for equilibrium and low frequency (<200 kHz) magnetohydrodynamic mode analysis has been installed. These new diagnostics address the significant challenges of measuring magnetic field in a high noise environment with the majority constrained to fit in an 8 mm diagnostic gap on the high field side. Electrostatic switching noise generated by the 16 independent current feedback-controlled power supplies produces dV cm /dt ~ 1 kV/μs and volt level common mode noise on the magnetics. Immunity to this switching noise is accomplished through differential signal runs and signal processing, along with end-to-end electromagnetic interference shielding. The magnetic measurements are simultaneously digitized at 1 MHz and conditioned by precision 8 pole Butterworth filters with a corner frequency of 200 kHz to prevent aliasing down to the 16-bit level over the full passband. Ex-vessel calibrations of the B p coils were completed with a typical uncertainty of <0.5%. Stray toroidal field pickup from coil misalignment or positioning errors is corrected using a physics-based model. Comparisons of the corrected measurements to modeling agree to within 1.3% on average. This is within the 1.5% measurement uncertainty that a sensitivity analysis determined is needed for accurate fast boundary and equilibrium reconstruction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Effects of Injected Current Streams on MHD Equilibrium Reconstruction of Local Helicity Injection Plasmas in a Spherical Tokamak

Open field line currents are intrinsic to DC helicity injection plasma startup and pose a challenge for inferring the plasma equilibrium with standard reconstruction analysis. Local helicity injection (LHI) is a type of DC helicity injection which uses small, modular current sources to drive force-free current along helical field lines to produce tokamak plasmas. MHD modeling and magnetic measurements during LHI indicate the injected current streams remain coherent as helical structures on the outboard edge of a core toroidal plasma that is tokamak-like in a toroidally averaged sense. To extract core plasma equilibrium properties, external magnetic diagnostics corrected for contributions from the injected current streams are fitted by a standard Grad-Shafranov equilibrium code. An iterative approach for estimating and subtracting the stream contributions from the diagnostic signals is described and applied to a model equilibrium database to reduce systematic errors introduced by the streams. Convergence is usually attained with 2 to 4 iterations, with derived equilibrium parameters matching the prescribed axisymmetric core values to within estimated experimental uncertainties. Accurate recovery of core parameters occurs when the ratio of the net toroidal windup current from the streams to the core plasma current is less than 0.2, which is typically satisfied in most experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Stellarator equilibrium axis-expansion to all orders in distance from the axis for arbitrary plasma beta

A systematic theory of the asymptotic expansion of the magnetohydrostatics (MHS) equilibrium in the distance from the magnetic axis is developed to include arbitrary smooth currents near the magnetic axis. Compared with the vacuum and the force-free system, an additional magnetic differential equation must be solved to obtain the pressure-driven currents. It is shown that there exist variables in which the rest of the MHS system closely mimics the vacuum system. Thus, a unified treatment of MHS fields is possible. The mathematical structure of the near-axis expansions to arbitrary order is examined carefully to show that the double-periodicity of physical quantities in a toroidal domain can be satisfied order by order. The essential role played by the leading-order Birkhoff–Gustavson normal form in solving the magnetic differential equations is highlighted. Several explicit examples of vacuum, force-free and MHS equilibrium in different geometries are presented.

Physics↗
DOE OSTI2024

Two-stage crash process in resistive drift ballooning mode driven ELM crash

We report a two-stage crash process in edge localized mode (ELM) driven by resistive drift-ballooning modes (RDBMs) numerically simulated in a full annular torus domain with a scale-separated four-field reduced MHD (RMHD) model using the BOUT++ framework. In the early nonlinear phase, the small first crash is triggered by linearly unstable RDBMs, and m/n=2/1 magnetic islands are nonlinearly excited by nonlinear coupling of RDBMs as well as their higher harmonics. Here, m is the poloidal mode number, n is the toroidal mode number, the q = 2 rational surface exists near the pressure gradient peak, and q is the safety factor. Simultaneously, middle-n RDBM turbulence develops but is poloidally localized around X-points of the magnetic islands, leading to the small energy loss. The second large crash occurs in the late nonlinear phase. Higher harmonics of m/n=2/1 magnetic islands well develop around the q = 2 surface via nonlinear coupling and make the magnetic field stochastic by magnetic island overlapping. Turbulence heat transport develops at X-points of higher harmonics of m/n=2/1 magnetic islands, resulting in the turbulence spreading in the poloidal direction. The large second crash is triggered when the turbulence covers the whole poloidal region so that the magnetic island generation and magnetic field stochastization before the large crash can be interpreted as ELM precursors. It is concluded that the ELM trigger is attributed to the turbulent spreading in the poloidal direction in synchronization with the magnetic field stochastization and the crash is driven by E × B convection rather than the conventional Rechester–Rosenbluth anomalous electron heat transport.

Seto, H. (ORCID:0000000259278643)↗
DOE OSTI2024

Use of current-potential patches to obtain fundamental improvements to the coils of magnetic fusion devices

A central issue in the design of tokamaks or stellarators is the coils that produce the external magnetic fields. The freedom that remains unstudied in the design of coils is enormous. This freedom could be quickly studied computationally at low cost with high reliability. In particular, the space between toroidal field and modular coils that blocks access to the plasma chamber could be increased by a large factor. The concept of current-potential patches, which was developed in Todd Elder's thesis, provides a method for separating the study of the feasibility of coils with attractive physics properties from the engineering design of specific coils. In addition to enhanced accessibility, coils can be designed for increased plasma-coil separation, insensitivity to coil position errors, and plasma control.

Boozer, Allen H. [Columbia Univ., New York, NY (Un↗
DOE OSTI2023

Self-Consistent Simulation of the Excitation of Compressional Alfvén Eigenmodes and Runaway Electron Diffusion in Tokamak Disruptions

Alfvénic modes in the current quench (CQ) stage of the tokamak disruption have been observed in experiments. In DIII-D the excitation of these modes is associated with the presence of high-energy runaway electrons (REs), and a strong mode excitation is often associated with the failure of RE plateau formation. In this work we present results of self-consistent kinetic-MHD simulations of RE-driven compressional Alfvén eigenmodes (CAEs) in DIII-D disruption scenarios, providing an explanation of the CQ modes. Simulation results reveal that high energy trapped REs can have resonance with the Alfvén mode through their toroidal precession motion, and the resonance frequency is proportional to the energy of REs. The mode frequencies and their relationship with the RE energy are consistent with experimental observations. Further, the perturbed magnetic fields from the modes can lead to spatial diffusion of REs including the nonresonant passing ones, thus providing the theoretical basis for a potential approach for RE mitigation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Demonstration of reduced neoclassical energy transport in Wendelstein 7-X

Research on magnetic confinement of high-temperature plasmas has the ultimate goal of harnessing nuclear fusion for the production of electricity. Although the tokamak is the leading toroidal magnetic-confinement concept, it is not without shortcomings and the fusion community has therefore also pursued alternative concepts such as the stellarator. Unlike axisymmetric tokamaks, stellarators possess a three-dimensional (3D) magnetic field geometry. The availability of this additional dimension opens up an extensive configuration space for computational optimization of both the field geometry itself and the current-carrying coils that produce it. Such an optimization was undertaken in designing Wendelstein 7-X (W7-X), a large helical-axis advanced stellarator (HELIAS), which began operation in 2015 at Greifswald, Germany. A major drawback of 3D magnetic field geometry, however, is that it introduces a strong temperature dependence into the stellarator’s non-turbulent ‘neoclassical’ energy transport. Indeed, such energy losses will become prohibitive in high-temperature reactor plasmas unless a strong reduction of the geometrical factor associated with this transport can be achieved; such a reduction was therefore a principal goal of the design of W7-X. In spite of the modest heating power currently available, W7-X has already been able to achieve high-temperature plasma conditions during its 2017 and 2018 experimental campaigns, producing record values of the fusion triple product for such stellarator plasmas. The triple product of plasma density, ion temperature and energy confinement time is used in fusion research as a figure of merit, as it must attain a certain threshold value before net-energy-producing operation of a reactor becomes possible. Here we demonstrate that such record values provide evidence for reduced neoclassical energy transport in W7-X, as the plasma profiles that produced these results could not have been obtained in stellarators lacking a comparably high level of neoclassical optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Electric field and turbulence in global Braginskii simulations across the ASDEX Upgrade edge and scrape-off layer

Turbulence simulations in diverted geometry across the edge and scrape-off layer (SOL) of ASDEX Upgrade are performed with the GRILLIX code (Stegmeir et al 2019 Phys. Plasmas 26 052517). The underlying global (full-f) drift-reduced Braginskii model allows to concurrently study the self-consistent dynamics of the turbulence and the background as well as the evolution of toroidal and zonal flows. Different contributions to the radial electric field are identified. The dominant contribution on closed flux surfaces comes from the ion pressure gradient, due to the diamagnetic drift in the curved magnetic field. Large deviations can be induced, in particular, by the polarization particle flux, leading to zonal flows. The latter are driven by small-scale eddies, but do not exhibit much impact on the overall transport which is driven by ballooning modes at larger scales. Ion viscosity is found to be important in damping poloidal rotation through adjusting of the parallel velocity profile, but not via direct vorticity damping. The zonal flow drive peaks at the separatrix, where a strong shear layer forms due to the sheath-induced counter-propagating SOL flow, allowing for the formation of a transport barrier. The temperature profile across the separatrix is determined by the competition between cross-field transport and outflow in the SOL, the latter being largely controlled by the parallel heat conductivity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

A Low Loss, Fast Ramp HTS Solenoid Prototype for Compact Spherical Tokamaks

To explore AC loss characteristics of commercial Bi-2212 conductors made in the USA, a small ohmic heating (OH) prototype coil is being developed for compact tokamaks using high performance Bi-2212 wires newly manufactured by Bruker EST. Here, the goal is to build and test Bi-2212 prototype coils for fast ramp central solenoid (CS) operation of spherical tokamaks (ST). For a ST pilot plant, one of the key challenges is the tight space available for CS in a center column where in-board legs of the toroidal field (TF) and OH coils are assembled into a cryostat with a few millimeters clearance. A high current density Rutherford cable consisting of 17 legacy Bi-2212 wires, a contribution of the US Magnet Development Program (MDP), was fabricated at the Lawrence Berkeley National Laboratory (LBNL) for the development of fast ramp CS magnets in the ST Advanced Reactor (STAR) designed by PPPL. Characterization of the new and legacy HTS wires is presented, along with the coil design wound directly with new wire, and a Rutherford cable-wound solenoid using the legacy wire. A comparison of AC losses is made based on expected performance for a fast ramp ST operation. The HTS cable and cabled coil fabrication issues such as the heat treatment, coil mandrel design, assembly and 2212 quench characteristics are discussed, along with a test plan for CS design validation.

Zhai, Yuhu [Princeton Plasma Physics Laboratory (P↗
DOE OSTI2020

Error field correction strategies in preparation to MAST-U operation

In magnetic fusion devices, unwanted non-axisymmetric magnetic field perturbations, known as error fields (EF), can have detrimental effects on plasma stability and confinement. Such EFs may originate from several sources, such as axi-symmetric coil misalignments, coil feeds, 3D structures in the wall surrounding the plasma, presence of ferromagnetic materials near the plasma surface, blanket materials in future devices. To minimize their impact on plasma performance and on the available operational space, it is important to identify the EF sources and develop EF control strategies. MAST Upgrade (MAST-U) is a spherical tokamak which will operate in the near future after a series of enhancements from the previous MAST experiment. To deliver a machine with EF amplitude low enough to allow a high quality experimental programme, systematic analysis of the intrinsic EF sources has been carried out for poloidal eld (P) and divertor (D) coils. To deliver a machine with EF amplitude low enough to allow a high quality experimental programme, systematic analysis of the intrinsic EF sources has been carried out for poloidal field (P) and divertor (D) coils, whose magnetic eld measurements were available when writing this paper. Such measurements reveal that P and D coils are 3D deformed and thus are responsible for intrinsic EFs, with mainly n=1 and n=2 toroidal mode numbers. In preparation to MAST-U operation, both passive and active EF control strategies have been adopted for n=1 and n=2 EFs compensation. Passive EF control consisted of installing finely each D and P coil within MAST-U device so as to minimize the intrinsic n=1 EF amplitude. The optimal coil alignment has been determined based on magnetic eld measurements and the corresponding 3D electro-magnetic modelling, and envisaged coil shift and tilt of 3.2 mm and 0.7 mrad, respectively, in the case of P coil named P4. Conversely, active control will be use during MAST-U operation to reduce the n=2 EF which is associated mainly with P4 and P5, as well. Since these coils have been re-used from the MAST device, studies attempting n=2 EF control, based on MAST plasmas, have been modelled utilizing the MARS-F code to interpret experimental results in MAST and to give hints for future n=2 EF control studies in MAST-U. Here, a model-based control set for n=2 EF control has been identified which would allow for the minimization of rotation braking, of the resonant magnetic field at the q=2 and of the plasma displacement in MAST-U, simultaneously.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Coordinate parameterisation and spectral method optimisation for Beltrami field solver in stellarator geometry

The numerical solution of the stepped pressure equilibrium (Hudson et al 2012 Phys. Plasmas 19 112502) requires a fast and robust solver to obtain the Beltrami field in three-dimensional geometry such as stellarators. The spectral method implemented in the stepped pressure equilibrium code (SPEC) is efficient when the domain is a hollow torus, but ill-conditioning of the discretised linear equations occurs in the (solid) toroid due to the artificially singular coordinate parameterisation near the axis. Here, we propose an improved choice for the reference axis to prevent coordinates surfaces from overlapping. Then, we examine the parity and asymptotics of the magnetic vector potential near the axis and suggest the use of recombined and rescaled Zernike radial basis functions. The maximum relative error in the magnetic field of the Wendelstein 7-X geometry is shown to reach 10 –9 at high resolution in a series of convergence tests and benchmarks against the boundary integral equation solver for Taylor states. The new method is also reported to significantly improve the accuracy of multi-volume SPEC calculations. A comparison between free-boundary SPEC and the analytical Dommaschk potential is presented with higher-than-usual Fourier resolution. It is illustrated that we are able to resolve low amplitude current sheets when an interface is placed where there is no flux surface in the analytic solution. This was previously concealed because of insufficient numerical resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

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↗