Search NASA⌕ Search

SEARCH · Search NASA

Results for “MAGNETIC EQUATOR”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Curvature of the chiral phase transition line from the magnetic equation of state of ( 2 + 1 )-flavor QCD

We analyze the dependence of the chiral phase transition temperature on baryon number and strangeness chemical potentials by calculating the leading order curvature coefficients in the light and strange quark flavor basis as well as in the conserved charge ( B , S ) basis. Making use of scaling properties of the magnetic equation of state (MEoS) and including diagonal as well as off-diagonal contributions in the expansion of the energylike scaling variable that enters the parametrization of the MEoS, allows to explore the variation of T c ( μ B , μ S ) = T c ( 1 − ( κ 2 B μ ^ B 2 + κ 2 S μ ^ S 2 + 2 κ 11 B S μ ^ B μ ^ S ) ) along different lines in the ( μ B , μ S ) plane. On lattices with fixed cutoff in units of temperature, a T = 1 / 8 , we find κ 2 B = 0.015 ( 1 ) , κ 2 S = 0.0124 ( 5 ) and κ 11 B S = − 0.0050 ( 7 ) . We show that the chemical potential dependence along the line of vanishing strangeness chemical potential is about 10% larger than along the strangeness neutral line. The latter differs only by about 3% from the curvature on a line of vanishing strange quark chemical potential, μ s = 0 . We also show that close to the chiral limit the strange quark mass contributes like an energylike variable in scaling relations for pseudocritical temperatures. The chiral phase transition temperature decreases with decreasing strange quark mass, T c ( m s ) = T c ( m s phy ) ( 1 − 0.097 ( 2 ) ( m s − m s phys ) / m s phy + O ( ( Δ m s ) 2 ) . Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

A robust fourth-order finite-difference discretization for the strongly anisotropic transport equation in magnetized plasmas

We propose a second-order temporally implicit, fourth-order-accurate spatial discretization scheme for the strongly anisotropic heat transport equation characteristic of hot, fusion-grade plasmas. Following Du Toit et al. (2018), the scheme transforms mixed-derivative diffusion fluxes (which are responsible for the lack of a discrete maximum principle) into nonlinear advective fluxes, amenable to nonlinear-solver-friendly monotonicity-preserving limiters. The scheme enables accurate multi-dimensional heat transport simulations with up to seven orders of magnitude of heat-transport-coefficient anisotropies with low cross-field numerical error pollution and excellent algorithmic performance, with the number of linear iterations scaling very weakly with grid resolution and grid anisotropy, and scaling with the square-root of the implicit timestep. We propose a multigrid preconditioning strategy based on a lower-order approximation that renders the scheme efficient and scalable under grid refinement. Several numerical tests are presented that display the expected spatial convergence rates and strong algorithmic performance, including fully nonlinear magnetohydrodynamics simulations of kink instabilities in a Bennett pinch in 2D helical geometry and of ITER in 3D toroidal geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A provably stable numerical method for the anisotropic diffusion equation in confined magnetic fields

We present a novel numerical method for solving the anisotropic diffusion equation in magnetic fields confined to a periodic box which is accurate and provably stable. We derive energy estimates of the solution of the continuous initial boundary value problem. A discrete formulation is presented using operator splitting in time with the summation by parts finite difference approximation of spatial derivatives for the perpendicular diffusion operator. Weak penalty procedures are derived for implementing both boundary conditions and parallel diffusion operator obtained by field line tracing. We prove that the fully-discrete approximation is unconditionally stable. Discrete energy estimates are shown to match the continuous energy estimate given the correct choice of penalty parameters. A nonlinear penalty parameter is shown to provide an effective method for tuning the parallel diffusion penalty and significantly minimises rounding errors. Several numerical experiments, using manufactured solutions, the “NIMROD benchmark” problem and a single island problem, are presented to verify numerical accuracy, convergence, and asymptotic preserving properties of the method. Finally, we present a magnetic field with chaotic regions and islands and show the contours of the anisotropic diffusion equation reproduce key features in the field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Prediction of the Wave Normal Angle of Proton‐Band EMIC Waves Near Geosynchronous Orbit

We investigate how the wave normal angle (WNA) and polarization of proton-band electromagnetic ion cyclotron (EMIC) waves change as they travel from their source to Earth. This paper marks a significant milestone as the first full-wave simulation of proton-band EMIC waves reflecting from the ionosphere. Our findings show that the WNA can change rapidly during propagation, primarily due to plasma inhomogeneities, such as variations in the Alfvén speed. The wave polarization is strongly related to the WNA, consistent with theory. Newly generated EMIC waves near the equator propagate with a WNA of 0° , then the WNA gradually shifts to 90° as they move toward Earth. In contrast, reflecting waves having 90° of WNA at Earth maintain a relatively larger WNA even near the magnetic equator. As a result, only the newly generated waves close to the source, where the magnetic latitude is less than approximately 20° , show left-handed polarization, while linear polarization remains dominant throughout the rest of the propagation.

EMIC wave↗
DOE OSTI2023

Code-verification techniques for the method-of-moments implementation of the magnetic-field integral equation

For computational physics simulations, code verification plays a major role in establishing the credibility of the results by assessing the correctness of the implementation of the underlying numerical methods. In computational electromagnetics, surface integral equations, such as the method-of-moments implementation of the magnetic-field integral equation, are frequently used to solve Maxwell's equations on the surfaces of electromagnetic scatterers. These electromagnetic surface integral equations yield many code-verification challenges due to the various sources of numerical error and their possible interactions. In this paper, we provide approaches to separately measure the numerical errors arising from these different error sources. We demonstrate the effectiveness of these approaches for cases with and without coding errors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Unraveling the Formation Region and Frequency of Chorus Spectral Gaps

Abstract The present study addresses two basic questions related to banded chorus waves in the Earth’s magnetosphere: 1) are chorus spectral gaps formed near the equatorial source region or during propagation away from the equator? and 2) why are chorus spectral gaps usually located below 0.5 f ce ( f ce : electron gyro‐frequency)? By analyzing Van Allen Probes data, we demonstrate that chorus spectral gaps are observed in the source region where chorus waves propagate both in the parallel and anti‐parallel directions to the magnetic field. Chorus spectral gaps below 0.5 f ce are associated with electron parallel acceleration at energies above the equatorial Landau resonant energies. We explain that initially generated chorus waves quickly isotropize the electron distribution through Landau resonant acceleration, and the isotropization occurs for higher energies at higher latitudes. The isotropized population, after returning to the magnetic equator, leads to a chorus gap typically below 0.5 f ce by suppressing wave excitation.

Li, Jinxing↗
DOE OSTI2022

An Eastward Current Encircling Mercury

Mercury has a terrestrial-like magnetosphere which is usually taken as a scaled-down-version of Earth's magnetosphere with a similar current system. We examine Mercury's magnetospheric current system based on a survey of Mercury's magnetic field measured by the Mercury Surface, Space Environment, Geochemistry, and Ranging spacecraft as well as computer simulations. We show that there is no significant Earth-like ring current flowing westward around Mercury, instead, we find, for the first time, an eastward current (EC) encircling the planet near the night-side magnetic equator with an altitude of ~500–1,000 km. The EC is closed with the dayside magnetopause current and could be driven by the gradient of plasma pressure as a diamagnetic current. Thus, Mercury's magnetosphere is not a scaled-down Earth magnetosphere, but a unique natural space plasma laboratory. Our findings offer fresh insights to analyze data from the BepiColombo mission, which is expected to orbit Mercury in 2025.

58 GEOSCIENCES↗
DOE OSTI2023

Statistical Study of the Nonthermal Continuum Radiation Beaming Angle Measured by the High Frequency Receiver on Van Allen Probes-A

The nonthermal continuum (NTC) radiation beaming angle is computed over the entire Van Allen Probes-A mission when the spacecraft was in the dawn sector. The conditions in the dawn sector are favorable for the wave vector to lie near/in the spacecraft's spin plane allowing a favorable estimate of the beaming angle, and the dawn sector is also advantageous in that previous studies show NTC occurrence to peak in this sector. We found that scatter plots, over the entire mission, of beaming angle versus magnetic latitude form a distinct inverted-V pattern, with the apex at/near the magnetic equator. This pattern was sharpest for frequencies ( f ) ≲ 100 kHz. Using the NTC beaming formula from the linear mode conversion theory (LMCT), we show that such an inverted-V pattern is expected due to the large variation in the plasmapause location over the entire mission. The theoretical derived pattern qualitatively reproduces the observed pattern but not quantitatively. The lack of quantitative agreement is discussed and is attributed to several factors, one factor is off-centered emissions from the radio window. The qualitative agreement strongly supports LMCT as being the dominant mechanism generating NTC for f ≲ 100 kHz. For f ≳ 100 kHz, the inverted-V pattern becomes less distinct, and strong near-equatorial beaming is observed. Finally, after considering contamination of our selections by left-handed polarized AKR, our study suggests that besides LMCT another unidentified NTC generation mechanism becomes important for f ≳ 100 kHz.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Observation and Numerical Simulation of Cold Ions Energized by EMIC Waves

This is the first report of significant energization (up to 7,000 eV) of low-energy He + ions, which occurred simultaneously with H-band electromagnetic ion cyclotron (EMIC) wave activity, in a direction mostly perpendicular to the ambient magnetic field. The event was detected by the Arase satellite in the dayside plasmatrough region off the magnetic equator on 15 May 2019. The peak energy of the He + flux enhancements is mostly above 1,000 eV. At some interval, the He + ions are energized up to ~7,000 eV. The H-band waves are excited in a frequency band between the local crossover and helium gyrofrequencies and are close to a linear polarization state with weakly left-handed or right-handed polarization. The normal angle of the waves exhibits significant variation between 0° and 80°, indicating a non-parallel propagation. Here, we run a hybrid code with parameters estimated from the Arase observations to examine the He + energization. The simulations show that cold He + ions are energized up to more than 1,000 eV, similar to the spacecraft observations. From the analysis of the simulated wave fields and cold plasma motions, we found that the ratio of the wave frequency to He + gyrofrequency is a primary factor for transverse energization of cold He + ions. As a consequence of the numerical analysis, we suggest that the significant transverse energization of He + ions observed by Arase is attributed to H-band EMIC waves excited near the local helium gyrofrequency.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

O-Mode Instability for Non-Thermal Continuum Radiation Generated Near the Plasmapause

Analysis of the ordinary mode (O-mode) instability is performed to comprehend the nonthermal continuum (NTC) radiation near the plasmapause, taking into account the relativistic wave-electron resonance effect. The energy source is the anisotropy in the velocity of the minority suprathermal electron population. Numerical solutions demonstrate that the O-mode can be unstable with multiple narrow frequency bands located close to harmonics of the electron cyclotron frequency above the local electron plasma frequency. These waves have narrow beaming angle bands of nearly 90° relative to the ambient magnetic field. Our findings indicate that NTC radiation generated by this wave-electron resonance instability near the plasmapause can propagate nearer to the magnetic equator with multiple harmonics, which is in agreement with a recent statistical study using Van Allen Probes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

An asymptotic-preserving semi-Lagrangian algorithm for the anisotropic heat transport equation with arbitrary magnetic fields

Here, we extend the recently proposed semi-Lagrangian algorithm for the extremely anisotropic heat transport equation [Chacón et al., J. Comput. Phys ., 272 (2014)] to deal with arbitrary magnetic field topologies. The original scheme (which showed remarkable numerical properties) was valid for the so-called tokamak-ordering regime, in which the magnetic field magnitude was not allowed to vary much along field lines. The proposed extension maintains the attractive features of the original scheme (including the analytical Green's function, which is critical for tractability) with minor modifications, while allowing for completely general magnetic fields. The accuracy and generality of the approach are demonstrated by numerical experiment with an analytical manufactured solution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Guiding center equations for the magnetic dipole

Since the discovery of Van Allen radiation belts in the 1960s, observations of energetic ions trapped in the Earth's dipole magnetic field have illustrated the remarkable confinement properties of this configuration. As such, it has been used for confining a hot plasma for nuclear fusion studies, starting from the pioneering work of Bo Lehnert and Akira Hasegawa, in the Levitated Dipole Experiment (LDX) at MIT until 2011 and in the RT-1 experiment at the University of Tokyo. More recently, the dipole has been subject to a renewed interest for fusion studies by a couple of startups and for smaller applications as a cold plasma source. While the equilibrium and magneto-hydrodynamic stability of the dipole have been investigated quite in detail, neoclassical properties of the dipole are comparatively much less known: the dipole is more known in geophysics than in fusion science. For this reason, in this paper, we propose a set of Hamiltonian, guiding-center equations to describe the motion of electrons and ions in a magnetic dipole configuration. We also developed a code, and we show the main features of particle motion, benchmarking our results with the analytical solutions for the bounce and precession motion, which are well documented in the literature. We also draw some general conclusions for the neoclassical transport in usual toroidal confinement schemes, such as the tokamak and the stellarator, pointing out the unique advantages of the dipole in confining energetic particles.

Hamiltonian mechanics↗
DOE OSTI2024

Structure preserving transport stabilized compatible finite element methods for magnetohydrodynamics

Here, we present compatible finite element space discretizations for the ideal compressible magnetohydrodynamic equations. The magnetic field is considered both in div- and curl-conforming spaces, leading to a strongly or weakly preserved zero-divergence condition, respectively. The equations are discretized in space such that transfers between the kinetic, internal, and magnetic energies are consistent, leading to a preserved total energy. We also discuss further adjustments to the discretization required to additionally achieve magnetic helicity preservation. Finally, we describe new transport stabilization methods for the magnetic field equation which maintain the zero-divergence and energy conservation properties, including one method which also preserves magnetic helicity. The methods' preservation and improved stability properties are confirmed numerically using a steady state and a magnetic dynamo test case.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

A magnetic analog of pressure–strain interaction

We study the evolution equation for magnetic energy density for a non-relativistic magnetized plasma in the (Lagrangian) reference frame comoving with the electron bulk velocity. Analyzing the terms that arise due to the ideal electric field, namely, perpendicular electron compression and magnetic field line bending, we recast them to reveal a quantity with a functional form analogous to the often-studied pressure–strain interaction term that describes one piece of internal energy density evolution of the species in a plasma, except with the species pressure tensor replaced by the magnetic stress tensor. We dub it the “magnetic stress–strain interaction.” We discuss decompositions of the magnetic stress–strain interaction analogous to those used for pressure–strain interaction. These analogies facilitate the interpretation of the evolution of the various forms of energy in magnetized plasmas and should be useful for a wide array of applications, including magnetic reconnection, turbulence, collisionless shocks, and wave–particle interactions. We display and analyze all the terms that can change magnetic energy density in the Lagrangian reference frame of the electrons using a particle-in-cell simulation of magnetic reconnection.

Barbhuiya, M. Hasan (ORCID:0000000163301650)↗
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 OSTI2025

Impact of Magnetic-field-driven Anisotropies on the Equation of State Probed in Neutron Star Mergers

Binary neutron star mergers can produce extreme magnetic fields, some of which can lead to strong magnetar-like remnants. While strong magnetic fields have been shown to affect the dynamics of outflows and angular momentum transport in the remnant, they can also crucially alter the properties of nuclear matter probed in the merger. In this work, we provide a first assessment of the latter, determining the strength of the pressure anisotropy caused by Landau-level quantization and the anomalous magnetic moment. To this end, we perform the first numerical relativity simulation with a magnetic polarization tensor and a magnetic-field-dependent equation of state using a new algorithm we present here, which also incorporates a mean-field dynamo model to control the magnetic field strength present in the merger remnant. Our results show that—in the most optimistic case—corrections to the anisotropy can be in excess of 10% and are potentially largest in the outer layers of the remnant. This work paves the way for a systematic investigation of these effects.

General relativity↗