Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “EARTH MAGNETOSPHERE”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Scaling the Ion Inertial Length and Its Implications for Modeling Reconnection in Global Simulations

Abstract We investigate the use of artificially increased ion and electron kinetic scales in global plasma simulations. We argue that as long as the global and ion inertial scales remain well separated, (1) the overall global solution is not strongly sensitive to the value of the ion inertial scale, while (2) the ion inertial scale dynamics will also be similar to the original system, but it occurs at a larger spatial scale, and (3) structures at intermediate scales, such as magnetic islands, grow in a self‐similar manner. To investigate the validity and limitations of our scaling hypotheses, we carry out many simulations of a two‐dimensional magnetosphere with the magnetohydrodynamics with embedded particle‐in‐cell (MHD‐EPIC) model. The PIC model covers the dayside reconnection site. The simulation results confirm that the hypotheses are true as long as the increased ion inertial length remains less than about 5% of the magnetopause standoff distance. Since the theoretical arguments are general, we expect these results to carry over to three dimensions. The computational cost is reduced by the third and fourth powers of the scaling factor in two‐ and three‐dimensional simulations, respectively, which can be many orders of magnitude. The present results suggest that global simulations that resolve kinetic scales for reconnection are feasible. This is a crucial step for applications to the magnetospheres of Earth, Saturn, and Jupiter and to the solar corona.

Tóth, Gábor↗
DOE OSTI2017

PIC simulations of wave-particle interactions with an initial electron velocity distribution from a kinetic ring current model

Whistler wave-particle interactions play an important role in the Earth inner magnetospheric dynamics and have been the subject of numerous investigations. By running a global kinetic ring current model (RAM-SCB) in a storm event occurred on Oct 23–24 2002, we obtain the ring current electron distribution at a selected location at MLT of 9 and L of 6 where the electron distribution is composed of a warm population in the form of a partial ring in the velocity space (with energy around 15 keV) in addition to a cool population with a Maxwellian-like distribution. The warm population is likely from the injected plasma sheet electrons during substorm injections that supply fresh source to the inner magnetosphere. These electron distributions are then used as input in an implicit particle-in-cell code (iPIC3D) to study whistler-wave generation and the subsequent wave-particle interactions. Here, we find that whistler waves are excited and propagate in the quasi-parallel direction along the background magnetic field. Several different wave modes are instantaneously generated with different growth rates and frequencies. The wave mode at the maximum growth rate has a frequency around 0.62ω ce , which corresponds to a parallel resonant energy of 2.5 keV. Linear theory analysis of wave growth is in excellent agreement with the simulation results. These waves grow initially due to the injected warm electrons and are later damped due to cyclotron absorption by electrons whose energy is close to the resonant energy and can effectively attenuate waves. The warm electron population overall experiences net energy loss and anisotropy drop while moving along the diffusion surfaces towards regions of lower phase space density, while the cool electron population undergoes heating when the waves grow, suggesting the cross-population interactions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Simulating the effects of warm O + ions on the growth of electromagnetic ion cyclotron (EMIC) waves

Electromagnetic ion cyclotron (EMIC) waves are believed to play a crucial role in the dynamics of the Earth’s magnetosphere. It has been widely accepted that plasma compositions can influence the growth rate of EMIC waves in the inner magnetosphere, but how warm O + ions change the wave growth rate is not well known. In this study, we investigate the impact of ring current O + concentration on the EMIC growth rate during a specific storm when the O + ion flux significantly increases. We calculate the growth rate of EMIC waves by using physics-based ion distributions output from the Ring current–Atmosphere interactions Model with Self-Consistent magnetic (B) and Electric (E) fields (RAM-SCBE) model. The percentage of warm O + is parametrically varied to examine how the maximum EMIC growth rate changes over time and how the global distribution of the maximum EMIC growth rate is affected. We found that the maximum growth rate of H-band appears in the dusk-to-midnight sector near the plasmapause, while O-band is excited at a slightly outer region. The maximum growth rate of He-band is closely related to the cold plasma density. With the increase of warm O + density, the maximum growth rate of H-band and He-band EMIC wave is reduced, while that of O-band EMIC wave is increased, and the region with this wave excitation is widened. Furthermore, such variation in the maximum EMIC growth rate implies a potential impact on the associated wave–particle interactions and change of the decay rate in the ring current.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Statistical description of coalescing magnetic islands via magnetic reconnection

The physical picture of interacting magnetic islands provides a useful paradigm for certain plasma dynamics in a variety of physical environments, such as the solar corona, the heliosheath and the Earth's magnetosphere. In this work, we derive an island kinetic equation to describe the evolution of the island distribution function (in area and in flux of islands) subject to a collisional integral designed to account for the role of magnetic reconnection during island mergers. This equation is used to study the inverse transfer of magnetic energy through the coalescence of magnetic islands in two dimensions. We solve our island kinetic equation numerically for three different types of initial distribution: Dirac delta, Gaussian and power-law distributions. The time evolution of several key quantities is found to agree well with our analytical predictions: magnetic energy decays as $\tilde {t}^{-1}$ , the number of islands decreases as $\tilde {t}^{-1}$ and the averaged area of islands grows as $\tilde {t}$ , where $\tilde {t}$ is the time normalised to the characteristic reconnection time scale of islands. General properties of the distribution function and the magnetic energy spectrum are also studied. Finally, we discuss the underlying connection of our island-merger models to the (self-similar) decay of magnetohydrodynamic turbulence.

Physics↗
DOE OSTI2018

Test of Ion Cyclotron Resonance Instability Using Proton Distributions Obtained From Van Allen Probe-A Observations

Anisotropic velocity distributions of protons have long been considered as free energy sources for exciting electromagnetic ion cyclotron (EMIC) waves in the Earth's magnetosphere. Here we rigorously calculated the proton anisotropy parameter using proton data obtained from Van Allen Probe–A observations. The calculations are performed for times during EMIC wave events (distinguishing the times immediately after and before EMIC wave onsets) and for times exhibiting no EMIC waves. We find that the anisotropy values are often larger immediately after EMIC wave onsets than the times just before EMIC wave onsets and the non–EMIC wave times. The increase in anisotropy immediately after the EMIC wave onsets is rather small but discernible, such that the average increase is by ~15% relative to the anisotropy values during the non–EMIC wave times and ~8% compared to those just before the EMIC wave onsets. Based on the calculated anisotropy values, we test the criterion for ion cyclotron instability suggested by Kennel and Petschek by applying it to the EMIC wave events. We find that despite the weak increase in anisotropy, the majority of the EMIC wave events satisfy the instability criterion. We suggest that the proton distributions often remain close to the marginal state to ion cyclotron instability, and consequently, the proton anisotropy values should often be observed near threshold values for ion cyclotron instability. In conclusion, we demonstrate the usefulness and limitation of the instability criteria expressed in the form of an inverse relation between the anisotropy and plasma beta.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Drift Instabilities in Thin Current Sheets Using a Two–Fluid Model With Pressure Tensor Effects

The integration of kinetic effects in fluid models is important for global simulations of the Earth's magnetosphere. We use a two–fluid 10–moment model, which includes the pressure tensor and has been used to study reconnection, to study the drift kink and lower hybrid drift instabilities. Using a nonlocal linear eigenmode analysis, we find that for the kink mode, the 10–moment model shows good agreement with kinetic calculations with the same closure model used in reconnection simulations, while the electromagnetic and electrostatic lower hybrid instabilities require modeling the effects of the ion resonance using a Landau fluid closure. As a result, comparisons with kinetic simulations and the implications of the results for global magnetospheric simulations are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
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 OSTI2014

High and low frequency instabilities driven by counter-streaming electron beams in space plasmas

A four-component plasma composed of a drifting (parallel to ambient magnetic field) population of warm electrons, drifting (anti-parallel to ambient magnetic field) cool electrons, stationary hot electrons, and thermal ions is studied in an attempt to further our understanding of the excitation mechanisms of broadband electrostatic noise (BEN) in the Earth's magnetospheric regions such as the magnetosheath, plasmasphere, and plasma sheet boundary layer (PSBL). Using kinetic theory, beam-driven electrostatic instabilities such as the ion-acoustic, electron-acoustic instabilities are found to be supported in our multi-component model. The dependence of the instability growth rates and real frequencies on various plasma parameters such as beam speed, number density, temperature, and temperature anisotropy of the counter-streaming (relative to ambient magnetic field) cool electron beam are investigated. It is found that the number density of the anti-field aligned cool electron beam and drift speed play a central role in determining which instability is excited. Using plasma parameters which are closely correlated with the measurements made by the Cluster satellites in the PSBL region, we find that the electron-acoustic and ion-acoustic instabilities could account for the generation of BEN in this region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Ion acoustic solitons and supersolitons in a magnetized plasma with nonthermal hot electrons and Boltzmann cool electrons

Arbitrary amplitude, ion acoustic solitons, and supersolitons are studied in a magnetized plasma with two distinct groups of electrons at different temperatures. The plasma consists of a cold ion fluid, cool Boltzmann electrons, and nonthermal energetic hot electrons. Using the Sagdeev pseudo-potential technique, the effect of nonthermal hot electrons on soliton structures with other plasma parameters is studied. Our numerical computation shows that negative potential ion-acoustic solitons and double layers can exist both in the subsonic and supersonic Mach number regimes, unlike the case of an unmagnetized plasma where they can only exist in the supersonic Mach number regime. For the first time, it is reported here that in addition to solitions and double layers, the ion-acoustic supersoliton solutions are also obtained for certain range of parameters in a magnetized three-component plasma model. The results show good agreement with Viking satellite observations of the solitary structures with density depletions in the auroral region of the Earth's magnetosphere.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Electron heating during magnetic reconnection: A simulation scaling study

Electron bulk heating during magnetic reconnection with symmetric inflow conditions is examined using kinetic particle-in-cell simulations. Inflowing plasma parameters are varied over a wide range of conditions, and the increase in electron temperature is measured in the exhaust well downstream of the x-line. The degree of electron heating is well correlated with the inflowing Alfvén speed c{sub Ar} based on the reconnecting magnetic field through the relation ΔT{sub e}=0.033 m{sub i} c{sub Ar}{sup 2}, where ΔT{sub e} is the increase in electron temperature. For the range of simulations performed, the heating shows almost no correlation with inflow total temperature T{sub tot}=T{sub i}+T{sub e} or plasma β. An out-of-plane (guide) magnetic field of similar magnitude to the reconnecting field does not affect the total heating, but it does quench perpendicular heating, with almost all heating being in the parallel direction. These results are qualitatively consistent with a recent statistical survey of electron heating in the dayside magnetopause (Phan et al., Geophys. Res. Lett. 40, 4475, 2013), which also found that ΔT{sub e} was proportional to the inflowing Alfvén speed. The net electron heating varies very little with distance downstream of the x-line. The simulations show at most a very weak dependence of electron heating on the ion to electron mass ratio. In the antiparallel reconnection case, the largely parallel heating is eventually isotropized downstream due a scattering mechanism, such as stochastic particle motion or instabilities. The simulation size is large enough to be directly relevant to reconnection in the Earth's magnetosphere, and the present findings may prove to be universal in nature with applications to the solar wind, the solar corona, and other astrophysical plasmas. The study highlights key properties that must be satisfied by an electron heating mechanism: (1) preferential heating in the parallel direction; (2) heating proportional to m{sub i} c{sub Ar}{sup 2}; (3) at most a weak dependence on electron mass; and (4) an exhaust electron temperature that varies little with distance from the x-line.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Demonstration of a viable quantitative theory for interplanetary type II radio bursts

Between 29 November and 1 December 2013 the two widely separated spacecraft STEREO A and B observed a long lasting, intermittent, type II radio burst for the extended frequency range ≈ 4 MHz to 30 kHz, including an intensification when the shock wave of the associated coronal mass ejection (CME) reached STEREO A. We demonstrate for the first time our ability to quantitatively and accurately simulate the fundamental (F) and harmonic (H) emission of type II bursts from the higher corona (near 11 solar radii) to 1 AU. Our modeling requires the combination of data-driven three-dimensional magnetohydrodynamic simulations for the CME and plasma background, carried out with the BATS-R-US code, with an analytic quantitative kinetic model for both F and H radio emission, including the electron reflection at the shock, growth of Langmuir waves and radio waves, and the radiations propagation to an arbitrary observer. The intensities and frequencies of the observed radio emissions vary hugely by factors ≈ 10{sup 6} and ≈ 10{sup 3}, respectively; the theoretical predictions are impressively accurate, being typically in error by less than a factor of 10 and 20 %, for both STEREO A and B. We also obtain accurate predictions for the timing and characteristics of the shock and local radio onsets at STEREO A, the lack of such onsets at STEREO B, and the z-component of the magnetic field at STEREO A ahead of the shock, and in the sheath. Very strong support is provided by these multiple agreements for the theory, the efficacy of the BATS-R-US code, and the vision of using type IIs and associated data-theory iterations to predict whether a CME will impact Earth’s magnetosphere and drive space weather events.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Processes setting the structure of the electron distribution function within the exhausts of anti-parallel reconnection

In situ spacecraft observations within the exhausts of magnetic reconnection document a large variation in the velocity space structure of the electron distribution function. Multiple mechanisms help govern the underlying electron dynamics, yielding a range of signatures for collisionless reconnection. These signatures include passing beams of electrons separated by well-defined boundaries from betatron heated/cooled trapped electrons. The present study emphasizes how localized regions of non-adiabatic electron dynamics can mix electrons across the trapped/passing boundaries and impact the form of the electron distributions in the full width of the exhaust. Here, while our study is based on 2D simulations, the described principles shaping the velocity space distributions also apply to 3D geometries making our findings relevant to spacecraft observation of reconnection in the Earth's magnetosphere.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Simulations of anti-parallel reconnection using a nonlocal heat flux closure

The integration of kinetic effects in fluid models is important for global simulations of the Earth's magnetosphere. In particular, it has been shown that ion kinetics play a crucial role in the dynamics of large reconnecting systems, and that higher-order fluid moment models can account for some of these effects. Here, we use a ten-moment model for electrons and ions, which includes the off diagonal elements of the pressure tensor that are important for magnetic reconnection. Kinetic effects are recovered by using a nonlocal heat flux closure, which approximates linear Landau damping in the fluid framework. Moreover, the closure is tested using the island coalescence problem, which is sensitive to ion dynamics. We also demonstrate that the nonlocal closure is able to self-consistently reproduce the structure of the ion diffusion region, pressure tensor, and ion velocity without the need for fine-tuning of relaxation coefficients present in earlier models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Generation of kinetic Alfvén waves in dayside magnetopause reconnection: A 3-D global-scale hybrid simulation

In this paper, we perform a three-dimensional (3-D) global-scale hybrid simulation to investigate the generation of kinetic Alfvén waves (KAWs) in reconnection at the dayside magnetopause. In the simulation, the magnetopause reconnection takes place due to a self-consistent interaction between the solar wind and Earth's magnetosphere, in which the solar wind is assumed to carry a southward interplanetary magnetic field. As multiple X-line reconnection occurs, shear Alfvénic perturbations are generated from reconnection, which propagate toward polar regions along the magnetic field lines. Our detailed analysis of the electromagnetic characteristics reveals the existence of kinetic Alfvén waves on the ion kinetic scales in these perturbations, in which the parallel electric field is generated. The 3-D structures and the Poynting fluxes in the shear Alfvén waves/KAWs in the global magnetopause are presented. It is shown that these KAWs carry electromagnetic energy into the cusp, which may lead to aurora brightening.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

An improved ten-moment closure for reconnection and instabilities

The integration of kinetic effects in fluid models is important for global simulations of Earth's magnetosphere. The use of the two-fluid ten-moment model, which includes the pressure tensor for both species, has had some success in simulating Ganymede and Mercury with a simple closure model. Here, we discuss a heat flux closure which accounts for some limitations of the earlier work while remaining computationally tractable. Comparisons with kinetic simulations for magnetic reconnection and lower-hybrid drift instabilities show good agreement with kinetic results and improvements on previous closure models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Nonlinear coupling of whistler waves to oblique electrostatic turbulence enabled by cold plasma

Kinetic simulations and theory demonstrate that whistler waves can excite oblique, short-wavelength fluctuations through secondary drift instabilities if a population of sufficiently cold plasma is present. The excited modes lead to heating of the cold populations and damping of the primary whistler waves. The instability threshold depends on the density and temperature of the cold population and can be relatively small if the temperature of the cold population is sufficiently low. Furthermore, this mechanism may thus play a significant role in controlling amplitude of whistlers in the regions of the Earth's magnetosphere where cold background plasma of sufficient density is present.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Scaling theory of three-dimensional magnetic reconnection spreading

We develop a first-principles scaling theory of the spreading of three-dimensional (3D) magnetic reconnection of finite extent in the out of plane direction. This theory addresses systems with or without an out of plane (guide) magnetic field, and with or without Hall physics. The theory reproduces known spreading speeds and directions with and without guide fields, unifying previous knowledge in a single theory. New results include: (1) Reconnection spreads in a particular direction if an x-line is induced at the interface between reconnecting and non-reconnecting regions, which is controlled by the out of plane gradient of the electric field in the outflow direction. (2) The spreading mechanism for anti-parallel collisionless reconnection is convection, as is known, but for guide field reconnection it is magnetic field bending. We confirm the theory using 3D two-fluid and resistive-magnetohydrodynamics simulations. (3) The theory explains why anti-parallel reconnection in resistive-magnetohydrodynamics does not spread. (4) The simulation domain aspect ratio, associated with the free magnetic energy, influences whether reconnection spreads or convects with a fixed x-line length. (5) We perform a simulation initiating anti-parallel collisionless reconnection with a pressure pulse instead of a magnetic perturbation, finding spreading is unchanged rather than spreading at the magnetosonic speed as previously suggested. The results provide a theoretical framework for understanding spreading beyond systems studied here, and are important for applications including two-ribbon solar flares and reconnection in Earth’s magnetosphere.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Implementation of a drive cylinder for low collisional experiments on magnetic reconnection

For most laboratory plasma experiments, Coulomb collisions between the particle species are sufficiently frequent that the particle distribution functions are relaxed to a near-Maxwellian form. This hampers the applicability of such experiments to phenomena observed in tenuous and near-collisionless space plasma. The Terrestrial Reconnection EXperiment (TREX) at the Wisconsin Plasma Physics Laboratory aims to study collisionless reconnection for parameters relevant to the Earth’s magnetosphere. To reduce the role of collisional effects, a reconnection Drive Cylinder has been developed, which increases both the effective system size of the TREX configuration and the rate at which reconnection can be driven. These two effects now permit TREX to reach a kinetic reconnection regime where collisional effects are minimized. The Drive Cylinder is comprised of 12 single loop drive-coils connected in parallel to a 10 kV capacitor bank. Insulated sheets of aluminum are applied to smooth the magnetic fields and enhance the drive efficiency. Following is a description of the technical details and performance of the Drive Cylinder.

Instruments & Instrumentation↗