Multifrequency compressional magnetic field oscillations and their relation to multiharmonic toroidal mode standing Alfvén waves
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
In this work we describe a simulation framework for testing the fast ion loss mechanism associated with the experimentally observed beam ion losses when an externally applied toroidally rotating perturbed magnetic field is used to control edge localized modes in the KSTAR tokamak. The simulations reproduce the key qualitative features of neutral beam injection (NBI) ion detection by a fast ion loss detector. The NBI ion losses in the simulation mainly occur for passing particles due to orbit stochastization, which is caused by orbit resonance with the 3D field perturbations. The relative toroidal angle of NBI ion deposition to the orbit island determines the radial path of the particles between confinement or loss. The fast ion loss quantity depends on the toroidal angle rotation of the 3D field with respect to the beam deposition position. The averaged transport of toroidal canonical angular momentum reveals that linear and nonlinear resonance of the NBI ions within the 3D field are the main factors determining fast ion transport and loss.
The Thomas Jefferson National Accelerator Facility (JLab) has developed a unique spectrometer system to study the weak interaction between electrons. The "Measurement of Lepton-Lepton Electroweak Reaction" (MOLLER) experiment, utilizing JLab's recent 12 GeV electron beam upgrade, is scheduled to operate for three years. Central to the MOLLER experiment are five water-cooled toroidal magnets, each with a unique geometry and seven-fold symmetry, designed to focus the particles. These magnets generate the magnetic field needed to separate incident beam electrons scattered from target electrons (Møller scattering) and protons (elastic e-p scattering) within a liquid hydrogen target. Here, this paper details the magnet field measuring technique developed to map all five MOLLER toroidal magnets at multiple locations inside and along the bore. It covers the design, mounting, and operation of the probe, along with the calibration procedure to determine the field and to prepare a field map for GEANT4 analysis. Additionally, the paper addresses the challenges of accurately measuring low magnetic fields.
In this paper results of a preliminary investigation into the existence of large scale magnetic fields in the sun are given. Using a kinematic model with prescribed internal rotation and a standard solar model, the poloidal and toroidal components of the magnetic field are calculated. The basic decay time is of the order of the age of the sun. In addition the fields are quite sensitive to slight variations in the internal rotation. The boundary condition at the solar surface also does not seem to influence the inner regions where large scale fields seem possible.
Abstract The safe operation of most tokamaks, especially the large ones, relies on the feedback control of vertical displacement events (VDEs). However, most of these feedback control systems are based on axisymmetric VDE models. In this study, we use NIMROD simulations to study the role of non-axisymmetric perturbations in free drift vertical displacement events on EAST. The high- n modes in the non-axisymmetric VDE grow first, which drives the formation of high- n magnetic island chains. Subsequently, the magnetic island chains grow and overlap with each other, leading to the destruction of the magnetic flux surface, which induces a minor disruption and accelerates the start of the following major disruption. The magnetic island and the stochastic magnetic field allow the toroidally non-axisymmetric poloidal plasma current to jet towards the hoop force direction, forming finger-like and filamentary structures. Such a plasma current non-axisymmetry strongly depends on the anisotropy in the thermal transport coefficients.
Differentially rotating stellar interiors nonaxisymmetric perturbation stability in toroidal and poloidal magnetic fields
Reasons for describing cosmic hydromagnetic phenomena by introducing electric currents are summarized demonstrating that change from the traditional magnetic field description to a current description provides new aspects of cosmic electrodynamics. It is shown that the kink-instability of electric currents is the basic physical phenomenon responsible for magnetic flux generation in cosmic physics. A model which starts with a poloidal magnetic field is discussed; differential motions change the field configuration so that kinetic energy is transferred into magnetic energy of induced toroidal fields; the currents associated with these become unstable when the toroidal energy exceeds the poloidal energy resulting in an amplification of the original poloidal field.
Large scale flows within the solar convection zone are the primary drivers of the Sun's magnetic activity cycle and play important roles in shaping the Sun's magnetic field. Differential rotation amplifies the magnetic field through its shearing action and converts poloidal field into toroidal field. Poleward meridional flow near the surface carries magnetic flux that reverses the magnetic poles at about the time of solar maximum. The deeper, equatorward meridional flow can carry magnetic flux back toward the lower latitudes where it erupts through the surface to form tilted active regions that convert toroidal fields into oppositely directed poloidal fields. These axisymmetric flows are themselves driven by large scale convective motions. The effects of the Sun's rotation on convection produce velocity correlations that can maintain both the differential rotation and the meridional circulation. These convective motions can also influence solar activity directly by shaping the magnetic field pattern. While considerable theoretical advances have been made toward understanding these large scale flows, outstanding problems in matching theory to observations still remain.
Following Greene et al. [Phys. Fluids 14, 671 (1971)] and Connor et al. [Phys. Plasmas 31, 577 (1988); Plasma Phys. Control. Fusion 34, 161 (1992); and Nucl. Fusion 33, 1533 (1993)], the Grad-Shafranov equation for an axisymmetric tokamak plasma equilibrium is solved via an expansion in the, supposedly small, inverse aspect-ratio of the plasma, ϵ. The displacements of equilibrium magnetic flux-surfaces due to plasma shaping are assumed to be $\mathcal{O}$(ϵ) smaller than the minor radii of the surfaces, but no other restriction is placed on the nature of the shaping. The solution of the Grad-Shafranov equation is matched to a vacuum solution that extends to infinity, and consists of an expansion in toroidal functions. The external poloidal magnetic field generated by a finite set of discrete external poloidal magnetic field-coils is calculated, and incorporated into the toroidal function expansion. In this manner, the shape of a large aspect-ratio tokamak plasma is directly related to the currents flowing in the external poloidal field-coils. Finally, a pedestal in the plasma pressure, and the associated spike in the bootstrap current, are incorporated into the model.
We solve the force-free equation J x B = 0 for fields which are toroidally symmetric. The technique utilizes an expansion about a cylindrical field and is therefore valid or tori with a large aspect ratio such as long, thin, coronal loops. The calculation is performed in spatial toroidal coordinates, rather than in the flux coordinates used by previous authors; this allows direct calculation of the loci of flux surfaces and of surfaces of constant magnetic pressure. Our solutions differ significantly from toroidal fields in laboratories, which are in general not force-free. They are characterized by field lines whose projections in the poloidal planes are circles with centers displaced by varying distances from the axis of the torus. In general, flux surfaces do not correspond to surfaces of constant magnetic pressure. We have examined solutions corresponding to simple analytic zero-order cylindrical fields. For moderate twists in the zero-order (cylindrical) field, the magnetic pressure is larger on the inner toroidal radius. However, this effect diminishes with twist angle and in fact, for extreme initial twists, the magnetic pressure can be larger on the outer radius. We compare our results with previous work utilizing flux coordinates.
In a recent publication (Toler et al., J. Plasma Phys., vol. 89, issue 2, 2023, p. 905890210), we demonstrated that for axisymmetric geometries, the Kapur–Rokhlin quadrature rule provided an efficient and high-order accurate method for computing the normal component, on the plasma surface, of the magnetic field due to the toroidal current flowing in the plasma, via the virtual-casing principle. The calculation was indirect, as it required the prior computation of the magnetic vector potential from the virtual-casing principle, followed by the computation of its tangential derivative by Fourier differentiation, to obtain the normal component of the magnetic field. Our approach did not provide the other components of the virtual-casing magnetic field. In this letter, we show that a more direct and more general approach is available for the computation of the virtual-casing magnetic field. The Kapur–Rokhlin quadrature rule accurately calculates the principal value integrals in the expression for all the components of the magnetic field on the plasma boundary, and the numerical error converges at a rate nearly as high as the indirect method we presented previously.
Field-aligned Birkeland currents and the angle of the magnetic line twist were calculated for an axially symmetric pole-on magnetosphere (assumed to be in MHD equilibrium). The angle of the field line twist was shown to have a strong radial dependence on the axisymmetric magnetotail as well as on the ionospheric conductivity and the amount of thermal plasma contained in closed magnetotail flux tubes. The field line twist results from the planetary rotation, which leads to the development of a toroidal magnetic B-sub-phi component and to differentially rotating magnetic field lines. It was shown that the time development of the toroidal magnetic B-sub-phi component and the rotation frequency are related through an induction equation.
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.
Abstract The scaling laws for the magnetic compression of a toroidally rotating field reversed configuration (FRC) have been investigated in this work. The magnetohydrodynamics (MHD) simulations of the magnetic compression on rotating FRCs employing the NIMROD code (Sovinec et al 2004 J. Comput. Phys. 195 355), are compared with the Spencer’s one-dimensional (1D) theory (Spencer et al 1983 Phys. Fluids 26 1564) for a wide range of initial flow speeds and profiles. The toroidal flow can influence the scalings directly through the alteration of the compressional work as also evidenced in the 1D adiabatic model, and indirectly by reshaping the initial equilibrium. However, in comparison to the static initial FRC equilibrium cases, the pressure and the radius scalings remain invariant for the magnetic compression ratio B w 2 / B w 1 up to 6 in presence of the initial equilibrium flow, suggesting a broader applicable regime of the Spencer scaling law for FRC magnetic compression. The invariant scaling has been proven a natural consequence of the conservation of angular momentum of both fluid and magnetic field during the dynamic compression process.
A simple model for the production and propagation of magnetized loops in the Galactic center is developed. It is found that a magnetized model can reproduce the filamentary, asymmetric structures on all scales from 100 pc to subparsec scales in the Galactic center. The radio and molecular gas morphology on all scales inside the Galactic center lobe may be explained in terms of the expansion of strongly magnetized loops and their interaction with the 2-5 pc torus and 30-50 pc portion of the center lobe. The magnetic loops have a toroidal field which explains the braided appearance of the radio continuum emission on the 2 pc scale. The barlike feature on the 2 pc radio continuum maps can be viewed as a piece of a magnetic loop with field strength about 0.01 G which is colliding and reconnecting with the ionized inner edge of the 2 pc molecular torus. The loop-torus interaction feeds energetic particles and shear Alfven waves into the torus, and deposits mass and energy.
Global gyrokinetic particle simulations of electrostatic ion temperature gradient (ITG) instability show that the most unstable eigenmode is localized to some magnetic fieldlines or discrete locations on the poloidal plane in the Wendelstein 7-X (W7-X) stellarator due to its mirror-like magnetic fields, which vary strongly in the toroidal direction and induce coupling of more toroidal harmonics (n) to form the linear eigenmode than in the Large Helical Device (LHD) stellarator. Nonlinear electrostatic simulation results show that self-generated zonal flows are the dominant saturation mechanism for the ITG instabilities in both the LHD and W7-X. Additionally, radial widths of the fluctuation intensity in both the LHD and W7-X are significantly broadened from the linear phase to the nonlinear phase due to turbulence spreading. Finally, nonlinear spectra in the W7-X are dominated by low-n harmonics, which can be generated both by nonlinear toroidal coupling of high-n harmonics and by linear toroidal coupling with large amplitude zonal flows due to the 3D equilibrium magnetic fields.
In this paper we model a geometrically thin accretion disk interacting with an externally imposed, uniform, vertical magnetic field. The accretion flow in the disk drags and distorts field lines, amplifying the magnetic field in the process. Inside the disk the radial component of the field is sheared into a toroidal component. The aim of this work is to establish the character of the resultant magnetic field and its dependence on the disk's parameters. We concentrate on alpha-disks driven by turbulent viscosity. Axisymmetric, two-dimensional solutions are obtained without taking into account the back-reaction of the magnetic field on the structure of the disk. The character of the magnetic field depends strongly on the magnitude of the magnetic Prandtl number, P . We present two illustrative examples of viscous disks: a so-called 'standard' steady state model of a disk around a compact star (e.g., cataclysmic variable), and a steady state model of a proto-planetary disk. In both cases, P = 1, P = 10(sup -1), and P = 10(sup -2) scenarios are calculated. Significant bending and magnification of the magnetic field is possible only for disks characterized by P of the order of 10(sup -2). In such a case, the field lines are bent sufficiently to allow the development of a centrifugally driven wind. Inside the disk the field is dominated by its toroidal component. We also investigate the dragging of the magnetic field by a nonviscous protoplanetary disk described by a phenomenological model. This scenario leads to large distortion and magnification of the magnetic field.
The magnetic helicity of fully nonlinear, multi-dimensional Alfven simple waves are investigated, by using relative helicity formulae and also by using an approach involving poloidal and toroidal decomposition of the magnetic field and magnetic vector potential. Different methods to calculate the magnetic vector potential are used, including the homotopy and Biot-Savart formulas. Two basic Alfven modes are identified: (a) the plane 1D Alfven simple wave given in standard texts, in which the Alfven wave propagates along the z-axis, with wave phase varphi=k_0(z-lambda t), where k_0 is the wave number and lambda is the group velocity of the wave, and (b)\ the generalized Barnes (1976) simple Alfven wave in which the wave normal {bf n} moves in a circle in the xy-plane perpendicular to the mean field, which is directed along the z-axis. The plane Alfven wave (a) is analogous to the slab Alfven mode and the generalized Barnes solution (b) is analogous to the 2D mode in Alfvenic, incompressible turbulence. The helicity characteristics of these two basic Alfven modes are distinct. The helicity characteristics of more general multi-dimensional simple Alfven waves are also investigated. Applications to nonlinear Aifvenic fluctuations and structures observed in the solar wind are discussed.