Magnetic Drift Research Articles

Three-dimensional inhomogeneity of electron-temperature-gradient turbulence in the edge of tokamak plasmas

Nonlinear multiscale gyrokinetic simulations of a Joint European Torus edge pedestal are used to show that electron-temperature-gradient (ETG) turbulence has a rich three-dimensional structure, varying strongly according to the local magnetic-field configuration. In the plane normal to the magnetic field, the steep pedestal electron temperature gradient gives rise to anisotropic turbulence with a radial (normal) wavelength much shorter than in the binormal direction. In the parallel direction, the location and parallel extent of the turbulence are determined by the variation in the magnetic drifts and finite-Larmor-radius (FLR) effects. The magnetic drift and FLR topographies have a perpendicular-wavelength dependence, which permits turbulence intensity maxima near the flux-surface top and bottom at longer binormal scales, but constrains turbulence to the outboard midplane at shorter electron-gyroradius binormal scales. Our simulations show that long-wavelength ETG turbulence does not transport heat efficiently, and significantly decreases overall ETG transport—in our case by ∼40%—through multiscale interactions.

The 3D magnetic topology and plasma dynamics in open stochastic magnetic field lines

The thermal quench triggered by locked modes is known to be mainly due to open stochastic magnetic field lines connected to the wall boundary. It is essential to understand the 3D structure of open stochastic field lines since it determines the overall plasma dynamics in the system. In this study, we analyze the 3D magnetic topology for two key concepts, the connection length Lc and the effective magnetic mirror ratio Meff, and present a comprehensive picture of electron and ion dynamics related to the magnetic topology. The connection length determines the 3D structure of the ambipolar potential, and a sharp potential drop across distinct Lc regions induces the E × B transport and mixing across the field line. The confinement of electrons and ions along the field line is determined by the ambipolar potential and Meff configuration. Electron and ion temperatures in magnetic hills (Meff<1) are lower than in magnetic wells (Meff>1) because particles in magnetic hills are more likely to escape toward the wall boundary along the field line. The mixing between the magnetic wells and hills by E × B and magnetic drift motions results in collisionless detrapping of electrons and ions, which reduces their temperature efficiently. Numerical simulations of two different magnetic configurations demonstrate the importance of the collisionless detrapping mechanism, which could be the main cause of plasma temperature drop during the thermal quench.

Magnetic Field Drifts of Small HTS Dipole Magnet Under Repeated Excitation

Shielding-current-induced fields (SCIFs) in magnets wound with coated conductors deteriorate the field quality of the magnets. Particularly, in magnets that are excited repeatedly, and must generate time-varying magnetic fields, for example, in certain types of accelerator magnets, the influence of the temporal behavior of the SCIFs on magnetic field is quite complicated. We focused on the magnetic field drifts and conducted magnetic field measurements using a small dipole magnet wound with coated conductors cooled using a cryocooler. The magnet was operated under various patterns of excitation, and magnetic fields were measured using a rotating pickup coil, which enabled us to measure the dipole and sextupole components of the magnetic field. Numerical electromagnetic field analyses were conducted to examine the current distributions in conductors, which influenced the temporal behavior of the magnetic field.

Long Term Vortex Flow Evolution around a Magnetic Island in Tokamaks.

We present a gyrokinetic analysis of the vortex flow evolution in a magnetic island in collisionless tokamak plasmas. In a short term ω[over ¯]_{D}t<1, where ω[over ¯]_{D} is the secular magnetic drift of the orbit center, initial monopolar vortex flow approaches to its residual level determined by the neoclassical enhancement of polarization shielding after collisionless relaxation. The residual level depends on the location inside an island and is higher than the Rosenbluth-Hinton level [M.N. Rosenbluth and F.L. Hinton, Phys. Rev. Lett. 80, 724 (1998)PRLTAO0031-900710.1103/PhysRevLett.80.724] due to finite island width. In a long term ω[over ¯]_{D}t>1, the residual vortex flow evolves to a dipolar zonal-vortex flow mixture due to toroidicity-induced breaking of a helical symmetry. The mixture forms localized flow shear layers near the island separatrix away from X points. The deviation of the streamlines of the mixture flows from magnetic surfaces allows turbulence advection across the island. We expect a small island w≲qρ_{Ti}/s[over ^] provides a favorable condition for this mixture flow formation, while the monopolar vortex flow persists for a larger island.

Space Weather Effects from Observations by Moscow University Cubesat Constellation

Moscow State University is developing a project for a multi-satellite constellation intended for the monitoring of space radiation. A number of small satellites of CubeSat format were launched into selected orbits crossing the wide range of magnetic drift shells. The primary scope for the project is the operational monitoring of near-Earth’s radiation environment, i.e., fluxes of electrons and protons of Earth’s radiation belts and energetic particles of solar and galactic origin. To date, there are four CubeSat satellites operating in near-Earth orbits, which deliver scientific and telemetric data. Thus, for the first time, a unique multi-satellite constellation has been implemented, which makes it possible to simultaneously measure the particle and quantum fluxes at different areas in the near-Earth space using the same type of instruments. A special compact detector of gamma quanta and energetic charged particles (electrons and protons) DeCoR has been developed to carry out radiation monitoring by CubeSats. With their help, observations of various effects of space weather have been made. These effects include a variety of electron fluxes in the outer belt during geomagnetic activity in late November–early December 2021, filling of polar caps by solar energetic particles accelerated in flares occurring in late October–early November, and the existence of stable electron fluxes near the geomagnetic equator.

Analytical study on magnetic component of geodesic acoustic mode

To understand the role of geodesic acoustic modes (GAMs) in regulating the transport by interacting with background electromagnetic turbulence, a thorough theoretical understanding of the magnetic components of the mode is required. In this study, the magnetic components of GAM are analytically investigated under the gyrokinetic framework with both the and harmonics are considered, where is the poloidal mode number. With the quasi-neutrality condition and Ampère’s law, the amplitudes of various poloidal magnetic components are derived. It is shown that both and magnetic components exist and are dominated by the cosine and sine components, respectively. In addition, it is found that the amplitudes of all magnetic components increase with respect to the ratio of plasma pressure to magnetic pressure and safety factor Most importantly, the amplitude of magnetic component is significantly enhanced due to the coupling of magnetic drift frequency with the first and second harmonics of the distribution functions, thus it can be comparable to that of magnetic component under certain conditions. The various magnetic components derived in this study may provide some guidance for experimental measurements.

P312. Frequency Drift Correction in 31P Magnetic Resonance Spectroscopy

Phosphorus-31 bearing neurometabolites play an important role in sustaining normal brain activity. To measure in-vivo high energy phosphates, 31P MRS is widely applied. However, magnetic resonance imaging scanners are utilized by multiple research groups. Various clinical protocols are implemented sequentially. Some well-known pulse sequences, such as DTI and fMRI, induce magnetic field drift resulting from heated shimming elements. Frequency drift has been associated with a detrimental effect on 31P MRS spectral quality.

Tunable rheological behaviour of magnetized complex plasma

The shear viscosity (η) of a 3D liquid dusty plasma has been estimated as a function of magnetic field (B) and normalized ion flow velocity (M) from the simulation data using Green–Kubo formalism with the help of Langevin dynamics simulation. It has been shown that in the strongly correlated liquid state, complex plasma may exhibit sharp changes in viscosity with magnetic field and ion drift velocity. In presence of ion drift, an oscillatory and attractive wake potential develops among charged dust particles immersed in a plasma. The amplitude of this wake potential can be modulated by applying an external magnetic field. In this work, we explore how an external magnetic field influences the rheological properties of complex plasma via anisotropic wake potential. It is observed that the rheological property of such plasma depends on the dominant interaction operating among the particles and can be controlled by applying an external magnetic field. A novel regime of magnetic field is observed in which strongly correlated complex plasma liquid exhibits a sharp response to an external magnetic field. Due to this unique property, complex plasma may be used as a platform to study magneto-rheological characteristics of soft matter and there is a possibility of using dusty plasma as a magneto-rheological material in near future.

Prediction of divertor heat flux width for ITER pre-fusion power operation using BOUT++ transport code

Prediction of divertor heat flux width is performed for the first and the second pre-fusion power operation (PFPO) phases specified in the new ITER research plan using BOUT++ transport code (Li et al 2018 Comput. Phys. Commun. 228 69–82). The initial plasma profiles inside the separatrix are taken from CORSICA scenario studies. Transport coefficients in transport code are calculated by inverting the plasma profiles inside the separatrix and are assumed to be constants in the scrape-off-layer. An anomalous thermal diffusivity scan is performed with E × B and magnetic drifts. The results in two scenarios identifying two distinct regimes: a drift-dominant regime when diffusivity is smaller than the respective critical thermal diffusivity χ c and a turbulence-dominant regime when diffusivity is larger than it. The Goldston heuristic drift model and the ITPA multi-machine experimental scaling yield a lower limit of the width λ q . From transport simulations, we obtain the critical thermal diffusivity χ c = 0.5 m2 s−1 for the PFPO-1 scenario with toroidal magnetic field B = 1.77 T and plasma current I p = 5 MA, and χ c = 0.3 m2 s−1 for the PFPO-2 scenario with toroidal magnetic field B = 2.65 T and plasma current I p = 7.5 MA. Separatrix temperature and collisionality also have a significant impact on the heat flux width in the drift-dominant regime. The investigation clearly yields a scaling for critical thermal diffusivity using ITER scenarios with fixed safety factor q 95, major radius R, aspect ratio R/a, and the separatrix temperature T sep, and establishes the connection with CFETR and C-Mod discharges. This scaling implies that for a given tokamak device with q 95, R, R/a, and T sep fixed, a reduction of poloidal magnetic field by a factor of 3 leads to a 9 times higher critical value of thermal diffusivity χ c, possibly yielding a transition from turbulence- to drift-dominant regime.

Observational evidence of ring current in the magnetosphere of Mercury

The magnetic gradient and curvature drift of energetic ions can form a longitudinal electric current around a planet known as the ring current, that has been observed in the intrinsic magnetospheres of Earth, Jupiter, and Saturn. However, there is still a lack of observational evidence of ring current in Mercury’s magnetosphere, which has a significantly weaker dipole magnetic field. Under such conditions, charged particles are thought to be efficiently lost through magnetopause shadowing and/or directly impact the planetary surface. Here, we present the observational evidence of Mercury’s ring current by analysing particle measurements from MErcury Surface, Space Environment, GEochemistry, and Ranging (MESSENGER) spacecraft. The ring current is bifurcated because of the dayside off-equatorial magnetic minima. Test-particle simulation with Mercury’s dynamic magnetospheric magnetic field model (KT17 model) validates this morphology. The ring current energy exceeds 5times {10}^{10} J during active times, indicating that magnetic storms may also occur on Mercury.

Ion temperature anisotropy model with cross‐field drifts in the scrape‐off layer

AbstractAn ion temperature anisotropy model including cross‐field drifts has been implemented in the 2D transport code UEDGE. Two ion temperature equations are solved in the model to give separate ion temperatures perpendicular and parallel to the magnetic field, which results in an additional mirror force for ion parallel transport and also leads to an additional curvature drift in the perpendicular direction. In our previous work, an ion temperature anisotropy model was implemented and applied to the full 2D version of UEDGE including recycled neutrals. The simulation results showed that the mirror force due to anisotropy reduces ion parallel flow velocity, however, the simulations were presented without cross‐field electric and magnetic drifts. Nevertheless, in present medium‐size tokamaks or future compact fusion devices, cross‐fields are important. In this work, the anisotropy model is extended to be compatible with the inclusion of cross‐field drifts. Preliminary results on comparison between the model solutions with and without cross‐field drifts are presented.

Representation and modeling of charged particle distributions in tokamaks

Experimental diagnostics, analysis tools and simulations represent particle distributions in various forms and coordinates. Algorithms to manage these data are needed on platforms like the ITER Integrated Modeling & Analysis Suite (IMAS), performing tasks such as archiving, modeling, conversion and visualization. A method that accomplishes some of the required tasks for distributions of charged particles with arbitrarily large magnetic drifts in axisymmetric tokamak geometry is described here. Given a magnetic configuration, we first construct a database of guiding center orbits, which serves as a basis for representing particle distributions. The orbit database contains the geometric information needed to perform conversions between arbitrary coordinates, modeling tasks, and resonance analyses. Using that database, an imported or newly modeled distribution is mapped to an exact equilibrium, where the dimensionality is reduced to three constants of motion (CoM). The orbit weight is uniquely given when the input is a true distribution: one that measures the number of physical particles per unit of phase space volume. Less ideal inputs, such as distributions estimated without drifts, or models of particle sources, can also be processed. As an application example, we reconstruct the drift-induced features of a distribution of fusion-born alpha particles in a large tokamak, given only a birth profile, which is not a function of the alpha's CoM. Repeated back-and-forth transformations between CoM space and energy-pitch-cylinder coordinates are performed for verification and as a proof of principle for IMAS.

Normal‐ and Reversed‐Boomerang Stripes on Electron Pitch Angle Distributions: Solar Wind Dynamic Pressure Effect

AbstractIn the inner magnetosphere dominated by the dipole magnetic field, energetic electrons of 90° pitch angle drift faster than those of other pitch angles. Electron flux oscillations initiated by injections or drift resonances will be observed earlier at a 90° pitch angle, that is, 90°‐leading boomerang‐shaped stripes on pitch angle distributions(PADs) are expected. In this work, reversed‐boomerang stripes (90° observed last) are first reported from the observations of RBSP‐A at L‐shell ∼5.9 on 11 March 2016. The corresponding solar wind dynamic pressure is over 10 nPa, which suggests that strong compression on the magnetosphere changes the magnetic field and drift frequencies of energetic electrons. Test particles running in an image‐dipole magnetic field (off‐equatorial minima exists at dayside large L‐shell) reproduce the abnormal (90° slower) drift velocities at L‐shell &gt;5.8 but are normal at L‐shell &lt;4.5 dominated by a dipole field. Normal boomerang stripes on PADs are indeed observed by RBSP‐B at L‐shell ∼4.0, indicating the solar wind dynamic pressure effect.

Average Ring Current Response to Solar Wind Drivers: Statistical Analysis of 61 Days of ENA Images

AbstractThis paper presents the first determination of the average global and regional response of ring current (RC) flux intensity versus four key solar wind (SW) and interplanetary magnetic field (IMF) parameters: SW electric field (ESW), IMF polarity (BZ,IMF), SW pressure (PSW), and SW speed (VSW). We analyze energetic neutral atom (ENA) images from Two Wide‐angle Imaging Neutral‐atom Spectrometers (TWINS) to measure the RC response, and use 5‐min downstream‐propagated data from the Advanced Composition Explorer (ACE) and Wind to measure SW driving. Our imaging data comprises 61 events, including 1,838 global 2D equatorial maps of line‐of‐sight (LOS) averaged ion flux derived from background‐subtracted, 15 min TWINS ENA images spanning 1–30 keV. Average spatial and spectral distributions were organized into four driver strength groups for each SW parameter. ESW has the strongest correlation with RC ion flux (coefficients are [0.95, 0.94, 0.88] at [1, 16, 30] keV, respectively), confirming that magnetospheric convection is a primary driver. In case studies, the global RC responds rapidly (≤15 min) to changes in ESW ranging from strong to weak. PSW is correlated with RC intensity, though not as strongly as ESW, indicating the secondary importance of PSW‐driven compression in RC injections. Strong ESW increases lower‐energy ions farther eastward and higher‐energy ions closer to midnight, consistent with an energy‐dependent interplay between magnetic drift and convection. The strongest SW driving preferentially increases the least energetic ions, suggesting that quasi‐steady convection can load the nightside RC with cooler ions.

Suppression of Low-Frequency Magnetic Drift Based on Magnetic Field Sensitivity in K-Rb-21Ne Atomic Spin Comagnetometer

Suppression of Low-Frequency Magnetic Drift Based on Magnetic Field Sensitivity in K-Rb-²¹Ne Atomic Spin Comagnetometer

Formation of the radial electric field profile in the WEST tokamak

Sheared flows are known to reduce turbulent transport by decreasing the correlation length and/or intensity of turbulent structures. The transport barrier that takes place at the edge during improved regimes such as H mode, corresponds to the establishment of a large shear of the radial electric field. In this context, the radial shape of the radial electric field or more exactly of the perpendicular E × B velocity appears as a key element in accessing improved confinement regimes. In this paper, we present the radial profile of the perpendicular velocity measured using Doppler back-scattering system at the edge of the plasma, dominated by the E × B velocity, during the first campaigns of the WEST tokamak. It is found that the radial velocity profile is clearly more sheared in lower single null configuration (with the B × ∇B magnetic drift pointing toward the active X-point) than in upper single null configuration for ohmic and low current plasmas (B = 3.7 T and q 95 = 4.7), consistently with the expectation comparing respectively ‘favourable’ versus ‘unfavourable’ configuration. Interestingly, this tendency is sensitive to the plasma current and to the amount of additional heating power leading to plasma conditions in which the E × B velocity exhibits a deeper well in USN configuration. For example, while the velocity profile exhibits a clear and deep well just inside the separatrix concomitant with the formation of a density pedestal during L–H transitions observed in LSN configuration, deeper E r wells are observed in USN configuration during similar transitions with less pronounced density pedestal.

Fast simulations for large aspect ratio stellarators with the neoclassical code KNOSOS

In this work, a new version of KNOSOS is presented. KNOSOS is a low-collisionality radially-local, bounce-averaged neoclassical code that is extremely fast, and at the same time, includes physical effects often neglected by more standard codes: the component of the magnetic drift that is tangent to the flux-surface and the variation of the electrostatic potential on the flux-surface. An earlier version of the code could only describe configurations that were sufficiently optimized with respect to neoclassical transport. KNOSOS can now be applied to any large aspect ratio stellarator, and its performance is demonstrated by means of detailed simulations in the configuration space of Wendelstein 7-X.

Passive deconfinement of runaway electrons using an in-vessel helical coil

A helical coil designed to passively generate non-axisymmetric fields during a plasma disruption is shown (via electromagnetic analysis, linear MHD modeling, and relativistic drift orbit tracing) to be effective at deconfining runaway electrons (REs) on a time scale significantly faster than the plasma current quench. Magnetic equilibria from DIII-D RE-producing scenarios are used to calculate the toroidal electric field generated during the current quench phase of a disruption, which in turn drives current in the proposed n = 1 in-vessel helical coil, without the need for any external power supplies or disruption detection or prediction techniques. Simulations of the plasma evolution using the TokSys GS Evolve code predict the inductive coupling of coil currents up to 12% of the pre-disruption plasma current into the helical coil. The coil geometry is parametrically varied to maximize both the non-resonant and resonant components of the 3D magnetic perturbation, resulting in δB/B ≈ 10−2 and a vacuum island overlap width of up to 0.7ψ N . The REORBIT module of the MARS-F code is used to model the full non-axisymmetric magnetic field and trace RE drift orbits to determine the effect on RE deconfinement, with up to 70% of the RE orbits lost after 0.2 ms. A two-stage evolution of the RE orbit loss fraction is observed to be caused by resonant trapping between multiple magnetic island chains. Finally, electromagnetic and thermal stresses on the coil are calculated to be within operational limits for installation in DIII-D, and scale favorably to a reactor-size device. These findings motivate future experimental study of the helical coil concept in DIII-D or other tokamaks.

Development of a persistent-mode NMR magnet with superconducting joints between high-temperature superconductors

This paper describes the first persistent-mode medium magnetic field (400 MHz; 9.39 T) nuclear magnetic resonance (NMR) magnet which uses superconducting joints between high-temperature superconductors (HTSs). As the ultimate goal, we aim to develop a high-resolution 1.3 GHz (30.5 T) NMR magnet operated in the persistent-mode. The magnet requires superconducting joints between HTSs and those between an HTS and a low-temperature superconductor (LTS). Towards this goal, we have been developing persistent-mode HTS inner coils to be operated in a 400 MHz (9.39 T) NMR magnet and here we present the first prototype inner coil wound with a single piece (RE = rare earth)Ba2Cu3O7−x (REBCO) conductor. The coil and a REBCO persistent current switch are connected with intermediate grown superconducting joints with high critical currents in external magnetic fields. To evaluate the performance of the joints in an ultimately stable and homogeneous magnetic field, the coil is operated in the persistent-mode, generating 0.1 T, in a 9.3 T background magnetic field of a persistent-mode LTS outer coil. The magnetic field drift over two years of the 400 MHz LTS/REBCO NMR magnet is as small as ∼1 ppm, giving high-resolution NMR spectra. The magnetic field drift rate over the second year was 0.03 × 10−3 ppm h−1, which is more than three orders of magnitude smaller than that required for an NMR magnet, demonstrating that the superconducting joints function satisfactorily in a high-resolution NMR system. The corresponding joint resistance is inferred to be <10−14 Ω.

Feasibility study of fast ion loss diagnostics for CFQS by beam ion loss calculation on vacuum vessel

A feasibility study for installing fast-ion loss diagnostics (FILDs) is initiated to understand energetic ion loss process in a quasi-axisymmetric stellarator CFQS. The calculations of beam ion orbit, the guiding center orbit code including the slowing down by bulk plasma in the Boozer coordinates and the collisionless Lorentz orbit code in the Cartesian coordinates, are performed to estimate the toroidal/poloidal distribution of beam ion losses on a vacuum vessel. We observed that loss points on the vacuum vessel are mainly located at the upper side of torus according to gradient magnetic field drift direction. Based on the loss distribution on the vacuum vessel, a possible location of FILD for CFQS is investigated. Although the energy of lost beam ions is ∼10 keV, relatively high-energy (> 20 keV) components can be measured on the upper side of the torus in the low-density case. Furthermore, the radial or vertical dependence of beam ion loss flux at candidate movable scintillator-type FILD positions is obtained. A considerable number of particles reach the scintillator-type FILD position. The energy and pitch angle distribution of beam ions reaching the FILD position shows that scintillator-type FILDs can measure co-going transit beam ions with energy of ∼25 keV and pitch angle of ∼140 degrees and transition beam ions with energy of 10–15 keV and pitch angle of ∼80–100 degrees.