Solar cycle variation of plasma mass density in the outer magnetosphere: Magnetoseismic analysis of toroidal standing Alfvén waves detected by Geotail
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.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Published observations on the toroidal oscillations of the earth are critically reviewed. A supplementary analysis of the record obtained by the Lamont strain seismometer is presented. Eleven toroidal modes are identified, and it is concluded that the periods are known to within 1 per cent. A perturbation scheme involving the ratio of the angular velocity of the earth to the resonant frequency is used in calculating the effects due to the rotation of the earth on the resonant frequency. The free oscillations are viewed as a superposition of traveling waves. In a nonrotating system two traveling waves combine to produce a stationary standing wave. In a rotating system, the rotation distinguishes between waves that travel in the direction of rotation and those that travel in the opposite direction. Rotation removes a degeneracy and results in a splitting of a spectral peak of order l into 2 times l plus 1 peaks. The fractional displacement in frequency for the lowest-order toroidal oscillations is 1/206 and of the same order as the Q of the peak, so that splitting will probably not be observed in the toroidal oscillations. Viewed locally, rotation causes a particle to precess about a direction parallel to the axis of rotation. This precession will cause a variation of amplitude with time if the motion is recorded by an instrument with an anisotropic response function. Care is therefore needed in studying the time decay of a given spectral peak. Rotation also couples the normal coordinates so that a motion that is initially purely horizontal will develop a vertical component. It is expected that vertical seismometers should record particle motion with the toroidal frequencies. The perturbations of the toroidal oscillations due to core-mantle interaction are treated in detail. An exact expression is obtained for the rate of energy dissipated by a finitely conducting plate oscillating across a magnetic field. The energy dissipated at the core-mantle boundary due to viscous and hydromagnetic coupling is shown to be insignificant as compared with the energy dissipated within the mantle. The toroidal magnetic field leaking into the lower mantle combines with the dipole field, resulting in a stress on the mantle, tending to stiffen the lower boundary. The stress is of sufficient magnitude to produce a displacement toward higher frequency in the lower-order toroidal oscillations. Observations on the (sub 0) T (sub 2) oscillations lead to an estimate of the toroidal magnetic field in the lower mantle. A calculation of elastic energy in the low-order oscillations suggests a value of 10 (sup 18) ergs per cycles per hour for the energy density at low frequencies in the Chilean earthquake. Each mode of oscillation has a characteristic radial distribution of elastic energy associated with it. This distribution determines which parts of the earth contribute most heavily in determining a particular resonant frequency. The distribution of energy for the lower 17 modes for a homogeneous and a Gutenberg model earth is calculated. The resonant frequencies for models of the earth based on the Gutenberg and Lehmann distribution of elastic properties are presented. It is shown that the Gutenberg model earth fits the observations more closely than the Lehmann model and that a slight alteration of the Gutenberg model gives a significantly better fit to the observations. The alteration involves a lower shear-wave velocity in the lower mantle while the Gutenberg velocity distribution is maintained in the upper mantle. Various studies of the earth's oscillations coupled with surface-wave investigations substantiate Gutenberg's hypothesis of a layer of low velocity in the upper mantle. The physical conditions required for the formation of a region of low velocity are examined in detail. The results confirm Birch's earlier statement that a temperature gradient in excess of 6 degrees to 7 degrees per kilometer is needed to produce a decrease in velocity. The low-velocity layer does not require that the temperature approach or exceed the melting temperature. If tile upper mantle is homogeneous, the region of lower velocity should commence at the base of the crust and extend to 150 kilometers under the oceans and about 100 kilometers under continental regions. The distribution of thermal conductivity and radioactivity consistent with the low-velocity layer is also considered.
Abstract We assess the toroidal magnetic field B t asymmetry in DIII-D due to a misalignment of the toroidal field coils with respect to the poloidal magnetic field coils and vacuum vessel. The peak-to-peak variation of the divertor strike point (SP) radius is measured to be 1 cm, with an n = 1 toroidal pattern. We use the centre of a narrow carbon deposition band on tungsten-coated divertor tiles just inside the outer strike point (OSP) as a proxy for the divertor SP location. The band occurred in a series of reverse B t discharges with the OSP positioned on the divertor inserts due to strong E × B drift transport of C from the inner to the outer SP through the private flux region. The variation in band radius (and hence the magnetic SP) is a (4.89 ± 0.31) mm shift toward (310 ± 4)° toroidal direction. These measurements agree well with previous measurements of the 3D magnetic field distribution (Luxon 2003 Nucl. Fusion 43 1813), simulations performed by the mafot field line integration code, and recent Langmuir probe measurements in the small-angle-slot (SAS) divertor (Watkins et al 2019 Nucl. Mater. Energy 18 46). Comparison of these measurements in the SAS divertor also indicates that there is the possibility of a tilt (in conjunction with the shift) of the B t coil field of (0.04 ± 0.07)° towards the toroidal angle of (215 ± 25)°. Previous measurements suggested a field misalignment of (4.6 ± 0.3) mm in the 270° toroidal direction, and a tilt of (0.06 ± 0.02)° toward the 114° toroidal direction, which is similar to the results reported here. These studies will be important for better understanding the radial variation of the toroidal strike line in DIII-D, for designing the new generation of SAS divertor, and for developing an understanding of the impact of error fields on tokamaks with tightly baffled slot divertors.
Abstract Ion cyclotron emission (ICE) is widely observed from spatially localised minority energetic ion populations in toroidal magnetically confined fusion (MCF) plasmas, both tokamaks and stellarators. Its spectral structure is typically regular with narrow suprathermal peaks, whose frequency separation matches a local energetic ion cyclotron frequency. Here we report the first computational (fully nonlinear self-consistent kinetic particle-in-cell code) and analytical (linear magnetoacoustic cyclotron instability (MCI)) studies of ICE observations from cylindrical plasmas contained in the Large Plasma Device (LAPD). Because LAPD is cylindrical, the plasma physics giving rise to the observed ICE spectrum necessarily excludes toroidal effects. Our approach, previously successful for toroidal plasma ICE, assumes slab geometry and hence is well adapted to LAPD. ICE from LAPD is strongly electrostatic, as distinct from electromagnetic, and is driven by 15 keV beam ions for which the ratio of perpendicular speed to the local Alfven speed, v ⊥ / v A , is 0.15, lower than in MCF plasmas from which beam-driven ICE has previously been observed. Our results are in good agreement with these observations. There is congruence between simulated ICE spectra, obtained in the saturated nonlinear regime of our computations, and observed ICE spectra; the underlying physics is essentially the same as in toroidal plasmas; and there is alignment with linear analytical theory where appropriate. The present work establishes a mapping from the cylindrical LAPD ICE observations to toroidal MCF ICE observations. The LAPD spectra are instances of beam-driven MCI-type ICE in its sub-Alfvenic, predominantly electrostatic manifestation, which has precedents in MCF stretching back to the 1990s. An interesting corollary is that, for many purposes, ICE in toroidal MCF plasmas ‘might as well’ be occurring in a cylinder.
A perturbation technique is used to compute the rotational corrections to the nonradial oscillation spectrum of a realistic neutron star model. The l = l(0) oscillations are coupled to l = l(0) + or - 1 oscillations by the Coriolis force. To calculate the rotational corrections to the quadrupole, l(0) = 2, oscillation modes, previous calculations of the nonradial oscillation spectrum of this neutron star model are extended to include l(0) = 1 and l(0) = 3 modes, as well as many higher overtones at each value of l(o). As an example of this technique, the rotational corrections to two quadrupole toroidal modes are computed. For these modes the coupling to the spheroidal class of oscillations produced by the Coriolis force introduces a nonzero radial component to the velocity field, as well as a nonvanishing Lagrangian density perturbation. This result is used to compute the neutrino damping rates for these toroidal modes of a rotating neutron star. The neutrino damping time can approach the gravitational radiation damping time in rotating neutron stars if the central temperature is high enough. The rotationally induced coupling of spheroidal oscillations to toroidal modes can also produce significant displacements at the stellar surface in some of the toroidal modes. The damping rates due to this process for the rotationally corrected toroidal modes are estimated. The coupling of surface motions to the toroidal modes also have interesting implications for channeling energy, e.g., that associated with a glitch in the crust, to the surface of the star. Perhaps this might produce observable effects in the pulsar emission process or a gamma-ray burst event.
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.
Optimizing the differential pair connection scheme (i.e. the set of pairs) of a toroidal array of magnetic sensors dedicated to measuring slowly rotating asymmetric fields can enhance the mode number detection capability and failure-resilience. In this work the condition number obtained from singular value decomposition (SVD) of the design matrix is used as a metric to evaluate the quality of a connection scheme. A large number of possible pair connections are usually available, so evaluating all of them may require extensive use of computational resources and can be very time consuming. Alternative methods to reduce the number of pairs evaluated without losing the capabilities of toroidal mode detection are presented in this paper. Here, three examples of the applications of such analysis for the 3D magnetic diagnostic system of DIII-D are also presented: the addition of two new toroidal arrays with n > 3 detection capabilities, the modification of an existing toroidal array in the low field side of the machine to accommodate the addition of a helicon antenna, and the design of changes in several toroidal arrays in the high field side to accommodate the addition of a lower hybrid current drive antenna on the center post.
Both near-field (<1 m away toroidally from the pellet) and far-field (>1 m away toroidally from the pellet) poloidal and toroidal impurity (carbon ion) flows are measured using visible imaging and fast bolometry during a polypropylene pellet-induced rapid plasma shutdown in the DIII-D tokamak. In the near field, the pellet appears to increase poloidal flow in the ion diamagnetic direction, possibly due to the strong radial temperature gradient caused by the pellet ablation. In the far field, the poloidal impurity flow typically appears slower and in the opposite direction. Toroidal impurity flow appears to be strongly influenced by the plasma initial toroidal rotation, especially in the far field. These results demonstrate that rapid shutdown impurity flows are not necessarily global in structure but can be quite different close to and far from the injected pellet.
Edge localized modes (ELMs) suppression by resonant magnetic perturbations (RMPs) was studied with the non-linear magneto-hydro-dynamic (MHD) code JOREK for the ITER H-mode scenarios at 15 MA, 12.5 MA, 10 MA/5.3 T. In this study, the main aim was to demonstrate that ELMs can be suppressed by RMPs while the divertor 3D footprints of heat and particle fluxes remain within divertor material limits. The unstable peeling–ballooning modes responsible for ELMs without RMPs were modelled first for each scenario using numerically accessible parameters for ITER. Then the stabilization of ELMs by RMPs was modelled with the same parameters. RMP spectra, optimized by the linear MHD MARS-F code, with main toroidal harmonics N = 2, N = 3, N = 4 have been used as boundary conditions of the computational domain of JOREK, including realistic RMP coils, main plasma, scrape off layer (SOL) divertor and realistic first wall. The model includes all relevant plasma flows: toroidal rotation, two fluid diamagnetic effects and neoclassical poloidal friction. With RMPs, the main toroidal harmonic and the non-linearly coupled harmonics remain dominant at the plasma edge, producing saturated modes and a continuous MHD turbulent transport thereby avoiding ELM crashes in all scenarios considered here. The threshold for ELM suppression was found at a maximum RMP coils current of 45 kAt–60 kAt compared to the coils maximum capability of 90 kAt. In the high beta poloidal steady-state 10 MA/5.3 T scenario, a rotating QH-mode without ELMs was observed even without RMPs. In this scenario with RMPs N = 3, N = 4 at 20 kAt maximum current in RMP coils, similar QH-mode behaviour was observed however with dominant edge harmonic corresponding to the main toroidal number of RMPs. The present MHD modelling was limited in time by few tens of ms after RMPs were switched on until the magnetic energy of the modes saturates. As a consequence the thermal energy was still evolving on this time scale, far from the ITER confinement time scale and hence only the form of 3D footprints on the divertor targets can be indicated within this set-up. Also note, that the divertor physics was missing in this model, so realistic values of fluxes are out of reach in this modelling. However the stationary 3D divertor and particle fluxes could be simply extrapolated from these results to the stationary situation considering that a large power fraction should be radiated in the core and SOL and only about 50 MW power is going to the divertor, which is an arbitrary, but reasonable number used here. The 3D footprints with RMPs show the characteristic splitting with the main RMP toroidal symmetry. The maximum radial extension of the footprints typically was ~20 cm in inner divertor and ~40 cm in outer divertor with stationary heat fluxes decreasing further out from the initial strike point from ~5 MW m ₋2 to ~1 MW m ₋2 assuming a total power in the divertor and walls is 50 MW. The heat fluxes remain within the divertor target and baffle areas, however with rather small margin in the outer divertor which could be an issue for the first wall especially in transient regimes when part of the plasma thermal energy is released due to switching on the RMP coils. This fact should be considered when RMPs are applied with a more favorable application before or soon after the L–H transition, although optimization is required to avoid increasing the L–H power threshold with RMPs.
We introduce a new approach to measure the magnetic pitch angle profile in tokamak plasmas with Doppler backscattering (DBS), a technique traditionally used for measuring flows and density fluctuations. The DBS signal is maximised when its probe beam's wavevector is perpendicular to the magnetic field at the cutoff location, independent of the density fluctuations [Hillesheim \emph{et al} 2022 \emph{Nucl. Fusion} \textbf{55} 073024]. Hence, if one could isolate this effect, DBS would then yield information about the magnetic pitch angle. By varying the toroidal launch angle, the DBS beam reaches cutoff with different angles with respect to the magnetic field, but with other properties remaining similar. Hence, the toroidal launch angle which gives maximum backscattered power is thus that which is matched to the pitch angle at the cutoff location, enabling inference of the magnetic pitch angle. We performed systematic scans of the DBS toroidal launch angle for repeated DIII-D tokamak discharges. Experimental DBS data from this scan were analysed and combined with Gaussian beam-tracing simulations using the Scotty code [Hall-Chen \emph{et al} 2022 \emph{Plasma Phys. Control. Fusion} \textbf{64} 095002]. The pitch-angle inferred from DBS is consistent with that from magnetics-only and motional-Stark-effect-constrained (MSE) equilibrium reconstruction in the edge. In the core, the pitch angles from DBS and magnetics-only reconstructions differ by one to two degrees, while simultaneous MSE measurements were not available. The uncertainty in these measurements was under a degree; we show that this uncertainty is primarily due to the error in toroidal steering, the number of toroidally separated measurements, and shot-to-shot repeatability. We find that the error of pitch-angle measurements can be reduced by optimising the poloidal launch angle and initial beam properties. Since DBS has high spatial and temporal resolutions, is non-perturbative, does not require neutral beams, and is likely robust to neutron damage of and debris on the first mirrors, using DBS to measure the pitch angle in future fusion energy systems is especially appealing.
ABSTRACT Strong magnetic fields play an important role in powering the emission of neutron stars. Nevertheless, a full understanding of the interior configuration of the field remains elusive. In this work, we present general relativistic magnetohydrodynamics (MHD) simulations of the magnetic field evolution in neutron stars lasting ${\sim } {880}\,$ms (∼6.5 Alfvén crossing periods) and up to resolutions of $0.1155\,$km using Athena++. We explore two different initial conditions, one with purely poloidal magnetic field and the other with a dominant toroidal component, and study the poloidal and toroidal field energies, the growth times of the various instability-driven oscillation modes, and turbulence. We find that the purely poloidal setup generates a toroidal field, which later decays exponentially reaching $1{{\ \rm per\ cent}}$ of the total magnetic energy, showing no evidence of reaching equilibrium. The initially stronger toroidal field setup, on the other hand, loses up to 20 per cent of toroidal energy and maintains this state till the end of our simulation. We also explore the hypothesis, drawn from previous MHD simulations, that turbulence plays an important role in the quasi-equilibrium state. An analysis of the spectra in our higher resolution setups reveals, however, that in most cases we are not observing turbulence at small scales, but rather a noisy velocity field inside the star. We also observe that the majority of the magnetic energy gets dissipated as heat increasing the internal energy of the star, while a small fraction gets radiated away as electromagnetic radiation.
This document serves as the closeout report of DOE Grant Award No. DE-FG02-99ER54563 with project period 1 April 2015 through 31 March 2020. The project comprised the international collaboration between MIT and EU scientists on the JET facility to improve our understanding of the physics of energetic particle-wave interactions by measuring the damping rates of stable Alfvén Eigenmodes (AEs) and unstable energetic particle driven modes. More specifically, this project involved the continued participation of MIT, the Culham Centre for Fusion Energy, the Swiss Plasma Center, and theorists from various European laboratories and from UC Irvine. Past contributions from the University of São Paulo are gratefully acknowledged. The Alfvén Eigenmode Active Diagnostic (AEAD) on JET discharges was successfully upgraded in this period and has extracted physically useful information, in primarily deuterium plasma discharges. It is expected that such experiments would be continued during the DT campaign in CY 2021 to assess the damping rates of similar modes in the presence of alpha particles, a product of fusing burning plasma. Only JET would carry out such experiments in the near term in the world. It is important to note that this collaboration was continued under DOE Grant Award No. DE-SC0014264 from 1 April 2020 through the present. In this grant period, the AEAD was upgraded with a set of individual amplifiers for each of a set of six antennas in two toroidally opposite locations. These new amplifiers allowed targeted selection of the antennas’ toroidal spectrum for AEs of interest with toroidal mode number | n | ≤ 20. The resonant detection and measurement of the damping rates of AEs was obtained from magnetic probes that could be compared with theory and simulations. This is a key topic of investigation for ITER and all other next-step fusion experiments, where such modes will interact with energetic particles produced by fusion reactions (alpha particles), Neutral Beam Injection (NBI) and Ion Cyclotron Resonance Heating (ICRH). The majority of the observations during this period were of Toroidal Alfvén Eigenmodes (TAEs) as these are most commonly observed on JET; a number of measurements were made during dedicated TAE experiments. However, a new set of lower frequency band filters were procured with the goal of studying Geodesic Acoustic Modes (GAMs), Beta (Acoustic) AEs (BAE/BAAEs), and Reverse Shear AEs (RSAEs), also predicted by theory. Commissioning and optimizations were successfully completed in this period, and the diagnostic has been in operation during the more recent JET campaigns. Initial results have been obtained in dedicated TAE experiments and successfully compared to drift-kinetic theory. A wide range of theoretical studies have been undertaken in support of the upcoming JET campaigns experiments. To supplement the ongoing use of ideal MHD codes, such as MISHKA during the studies of TAEs, gyrokinetic simulations of low frequency AEs were performed in collaboration with UC Irvine. The Gyrokinetic Toroidal Code (GTC) was used to determine the structure, frequencies and stability of AEs in JET plasmas. Thus, a solid scientific foundation has been laid for future DT campaigns in JET in the 2021 operation period when the damping rates of relevant modes in the presence of alpha particles could be assessed.
A high-efficiency, extreme ultraviolet (EUV) imaging spectrometer has been constructed and tested. The spectrometer employs a concave toroidal grating illuminated at normal incidence in a Rowland circle mounting and has only one reflecting surface. The toroidal grating has been fabricated by a new technique employing an elastically deformable submaster grating which is replicated in a spherical form and then mechanically distorted to produce the desired aspect ratio of the toroidal surface for stigmatic imaging over the selected wavelength range. The fixed toroidal grating used in the spectrometer is then replicated from this surface. Photographic tests and initial photoelectric tests with a two-dimensional, pulse-counting detector system have verified the image quality of the toroidal grating at wavelengths near 600 A. The basic designs of two instruments employing the spectrometer for astrophysical investigations in space are described, namely, a high-resolution EUV spectroheliometer for studies of the solar chromosphere, transition region, and corona; and an EUV spectroscopic telescope for studies of nonsolar objects.
A high-efficiency, extreme ultraviolet (EUV) imaging spectrometer is constructed and tested. The spectrometer employs a concave toroidal grating illuminated at normal incidence in a Rowland circle mounting and has only one reflecting surface. The toroidal grating has been fabricated by a new technique employing an elastically-deformable sub-master grating replicated in a spherical form and then mechanically distorted to produce the desired aspect ratio of the toroidal surface for stigmatic imaging over the selected wavelength range. The fixed toroidal grating used in the spectrometer is then replicated from this surface. Photographic tests and initial photoelectric tests with a two-dimensional, pulse-counting detector system verify the image quality of the toroidal grating at wavelengths near 600 A. The results of these tests and the basic designs of two instruments which could employ the imaging spectrometer for astrophysical investigations in space are described; i.e., a high-resolution EUV spectroheliometer for studies of the solar chromosphere, transition region, and corona; and an EUV spectroscopic telescope for studies of non-solar objects.
This paper presents results of extensive analysis of mode excitation observed during the operation of the Alfvén Eigenmode Active Diagnostic (AEAD) in the JET tokamak during the 2019-2020 deuterium campaign. Here, six of eight toroidally spaced antennas, each with independent power and phasing, were successful in actively exciting stable MHD modes in 479 plasmas. In total, magnetic resonances were detected with up to fourteen fast magnetic probes. In this work, we present the calculations of resonant frequencies $f_0$, damping rates γ < 0, and toroidal mode numbers $\textit{n}$, spanning the parameter ranges $f_0$ ≈ 30–250 kHz, –γ ≈ 0–13 kHz, and $\vert n \vert \leq {30}$. In general, good agreement is seen between the resonant and the calculated toroidal Alfvén Eigenmode frequencies, and between the toroidal mode numbers applied by the AEAD and estimated of the excited resonances. We note several trends in the database: the probability of resonance detection decreases with plasma current and external heating power; the normalized damping rate increases with edge safety factor but decreases with external heating. These results provide key information to prepare future experimental campaigns and to better understand the physics of excitation and damping of Alfvén Eigenmodes in the presence of alpha particles during the upcoming DT campaign, thereby extrapolating with confidence to future tokamaks.
Abstract The influence of negative triangularity (NT) of the plasma shape on the n = 1 (n is the toroidal mode number) tearing mode (TM) stability has been numerically investigated, with results compared to that of the positive triangularity (PT) counterpart. By matching the safety factor profile for a series of toroidal equilibria, several important plasma parameters, including the triangularity, the plasma equilibrium pressure, the plasma resistivity as well as the toroidal rotation, have been varied. The TM localized near the plasma edge is found to be more unstable in the NT plasmas as compared to the PT counterpart. The fundamental reason for this difference is the lack of favorable average curvature stabilization in NT configurations. Direct comparison of the Mercier index corroborates this conclusion. For the core-localized mode, where the difference in the local triangularity between NT and PT becomes small and the curvature stabilization is significantly reduced, larger Shafranov shift in the plasma core associated with the NT configuration results in more stable TM. The plasma toroidal flow generally stabilizes the TM in plasmas with both NTs and PTs. The flow stabilization is however weaker in the case of negative triangularity with finite plasma pressure.