Search NASA⌕ Search

SEARCH · Search NASA

Results for “global magnetospheric models”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

Saturn's Auroral Response to the Solar Wind: Centrifugal Instability Model

We describe a model initially presented by Sittler et al. [2006] which attempts to explain the global response of Saturn's magnetosphere and its corresponding auroral behavior to variations in the solar wind. The model was derived from published simultaneous Hubble Space Telescope (HST) auroral images and Cassini upstream measurements taken during the month of January 2004. These observations show a direct correlation between solar wind dynamic pressure and (1) auroral brightening toward dawn local time, (2) an increase of rotational movement of auroral features to as much as 75% of the corotation speed, (3) the movement of the auroral oval to higher latitudes and (4) an increase in the intensity of Saturn Kilometric Radiation (SKR). This model is an alternative to the reconnection model of Cowley et al. [2004a,b; 2005] which is more Earth-like while ours stresses rotation. If angular momentum is conserved in a global sense, then when compressed the magnetosphere will tend to spin up and when it expands will tend to spin down. With the plasma sheet outer boundary at L approximates 15 we argue this region to be the dominant source region for the precipitating particles. If radial transport is dominated by centrifugal driven flux tube interchange motions, then when the magnetosphere spins up, outward transport will increase, the precipitating particles will move radially outward and cause the auroral oval to move to higher latitudes as observed. The Kelvin-Helmholtz instability may contribute to the enhanced emission along the dawn meridian as observed by HST. We present this model in the context of presently published observations by Cassini.

Sittler, Edward C.↗
NASA NTRS1994

Solution of the Chapman-Ferraro problem with an arbitrary magnetopause

We present a global model of the magnetic field of the magnetosphere that includes the effects of the Chapman-Ferraro currents at the magnetopause. In contrast to ealier models, the magnetopause shape is arbitrary, thus allowing the use of more realistic geometries. The internal magnetospheric field model of Hilmer and Voigt (1993), is completely shielded within the magnetopause by solving the Laplace equation with Neumann boundary conditions using a finite difference method on a non-orthogonal, curvilinear grid. The resulting model magnetosphere is perfectly closed although the method can also be applied with more general boundary conditions, to generate a set of open models based on the approach of Toffoletto and Hill (1989, 1993). The purpose of this paper is to demonstrate the feasibility of a purely numerical approach to solving the Chapman-Ferraro problem with arbitrary magnetopause shape and boundary conditions.

Toffoletto, F. R.↗
NASA NTRS2016

Global Effects of Transmitted Shock Wave Propagation Through the Earth's Inner Magnetosphere: First Results from 3-D Hybrid Kinetic Modeling

We use a new hybrid kinetic model to simulate the response of ring current, outer radiation belt, and plasmaspheric particle populations to impulsive interplanetary shocks. Since particle distributions attending the interplanetary shock waves and in the ring current and radiation belts are non-Maxwellian, waveparticle interactions play a crucial role in energy transport within the inner magnetosphere. Finite gyroradius effects become important in mass loading the shock waves with the background plasma in the presence of higher energy ring current and radiation belt ions and electrons. Initial results show that shocks cause strong deformations in the global structure of the ring current, radiation belt, and plasmasphere. The ion velocity distribution functions at the shock front, in the ring current, and in the radiation belt help us determine energy transport through the Earth's inner magnetosphere.

Lipatov, A. S.↗
DOE OSTI2021

A boundary value “reservoir problem” and boundary conditions for multi-moment multifluid simulations of sheaths

Multifluid simulations of plasma sheaths are increasingly used to model a wide variety of problems in plasma physics ranging from global magnetospheric flows around celestial bodies to plasma–wall interactions in thrusters and fusion devices. For multifluid problems, accurate boundary conditions to model an absorbing wall that resolves a classical sheath remain an open research area. This work justifies the use of vacuum boundary conditions for absorbing walls to show comparable accuracy between a multifluid sheath and lower moments of a continuum-kinetic sheath.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2017

A Method to Predict Magnetopause Expansion in Radial IMF Events by MHD Simulations

This paper presents a method for taking into account changes of solar wind parameters in the foreshock using global MHD simulations. We simulate four events with very distant subsolar magnetopause crossings that occurred during quasi-radial interplanetary magnetic field (IMF) intervals lasting from one to several hours. Using previous statistical results, we suggest that the density and velocity in the foreshock cavity decrease to approx. 60% and approx. 94% of the ambient solar wind values when the IMF cone angle falls below 50 deg. This diminishes the solar wind dynamic pressure to 53% and causes a corresponding magnetospheric expansion. We change the upstream solar wind parameters in a global MHD model to take these foreshock effects into account. We demonstrate that the modified model predicts magnetopause distances during radial IMF intervals close to those observed by THEMIS. The strong total pressure decrease in the data seems to be a local, rather than a global, phenomenon. Although the simulations with decreased solar wind pressure generally reproduce the observed total pressure in the magnetosheath well, the total pressure in the magnetosphere often agrees better with results for nonmodified boundary conditions. The last result reveals a limitation of our method: we changed the boundary conditions along the whole inflow boundary, although a more correct approach would be to vary parameters only in the foreshock. A model with the suggested global modification of the boundary conditions better predicts the location of part of the magnetopause behind the foreshock but may fail in predicting the rest of the magnetopause.

Samsonov, A. A.↗
NASA NTRS2003

Dayside Proton Aurora: Comparisons between Global MHD Simulations and Image Observations

The IMAGE mission provides a unique opportunity to evaluate the accuracy of current global models of the solar wind interaction with the Earth's magnetosphere. In particular, images of proton auroras from the Far Ultraviolet Instrument (FUV) onboard the IMAGE spacecraft are well suited to support investigations of the response of the Earth's magnetosphere to interplanetary disturbances. Accordingly, we have modeled two events that occurred on June 8 and July 28, 2000, using plasma and magnetic field parameters measured upstream of the bow shock as input to three-dimensional magnetohydrodynamic (MHD) simulations. This paper begins with a discussion of images of proton auroras from the FUV SI-12 instrument in comparison with the simulation results. The comparison showed a very good agreement between intensifications in the auroral emissions measured by FUV SI-12 and the enhancement of plasma flows into the dayside ionosphere predicted by the global simulations. Subsequently, the IMAGE observations are analyzed in the context of the dayside magnetosphere's topological changes in magnetic field and plasma flows inferred from the simulation results. Finding include that the global dynamics of the auroral proton precipitation patterns observed by IMAGE are consistent with magnetic field reconnection occurring as a continuous process while the iMF changes in direction and the solar wind dynamic pressure varies. The global simulations also indicate that some of the transient patterns observed by IMAGE are consistent with sporadic reconnection processes. Global merging patterns found in the simulations agree with the antiparallel merging model. though locally component merging might broaden the merging region, especially in the region where shocked solar wind discontinuities first reach the magnetopause. Finally, the simulations predict the accretion of plasma near the bow shock in the regions threaded by newly open field lines on which plasma flows into the dayside ionosphere are enhanced. Overall the results of these initial comparisons between global MHD simulation results and IMAGE observations emphasize the interplay between reconnection and dynamic pressure processes at the dayside magnetopause. as well as the intricate connection between the bow shock and the auroral region.

Berchem, J.↗
NASA NTRS2007

Bifurcation and Hysteresis of the Magnetospheric Structure with a varying Southward IMF: Field Topology and Global Three-dimensional Full Particle Simulations

Using a three-dimensional full electromagnetic particle model (EMPM), we have performed global simulations of the interaction between the solar wind and the terrestrial magnetosphere, and have investigated its asymptotic stability. The distance between the dayside magnetopause subsolar point and the Earth center, R(sub mp) is measured, as the intensity of southward IMF |B(sub z)| is slowly varying. Based on the field topology theory, one analyzes the variation of R(sub mp) as a reference index of the dynamics of this interaction, when IMF |B(sub z)| successively increases and decreases to its original value. Two striking results are observed. First, as the IMF |B(sub z)| increases above a critical value, the variation of R(sub mp) suddenly changes (so called 'bifurcation' process in field topology). Above this critical value, the overall magnetic field topology changes drastically and is identified as being the signature of magnetic reconnection at the subsolar point on the magnetopause. Second, this subsolar point recovers its original location R(sub mp) by following different paths as the IMF |B(sub z)| value increases (from zero to a maximum fixed value) and decreases (from this maximum to zero) passing through some critical values. These different paths are the signature of 'hysteresis' effect, and are characteristic of the so-called 'subcritical-type' bifurcation. This hysteresis signature indicates that dissipation processes take place via an energy transfer from the solar wind to the magnetosphere by some irreversible way, which leads to a drastic change in the magnetospheric field topology. This hysteresis is interpreted herein as a consequence of the magnetic reconnection taking place at the dayside magnetopause. The field topology reveals to be a very powerful tool to analyze the signatures of three-dimensional magnetic reconnection without the obligation for determining the mechanisms responsible for, and the consequences of the reconnection on the overall magnetospheric dynamics.

Cai, DongSheng↗
DOE OSTI2023

Three-dimensional Dynamics of Strongly Twisted Magnetar Magnetospheres: Kinking Flux Tubes and Global Eruptions

Abstract The origins of the various outbursts of hard X-rays from magnetars (highly magnetized neutron stars) are still unknown. We identify instabilities in relativistic magnetospheres that can explain a range of X-ray flare luminosities. Crustal surface motions can twist the magnetar magnetosphere by shifting the frozen-in footpoints of magnetic field lines in current-carrying flux bundles. Axisymmetric (2D) magnetospheres exhibit strong eruptive dynamics, i.e., catastrophic lateral instabilities triggered by a critical footpoint displacement of ψ crit ≳ π . In contrast, our new three-dimensional (3D) twist models with finite surface extension capture important non-axisymmetric dynamics of twisted force-free flux bundles in dipolar magnetospheres. Besides the well-established global eruption resulting (as in 2D) from lateral instabilities, such 3D structures can develop helical, kink-like dynamics, and dissipate energy locally (confined eruptions). Up to 25% of the induced twist energy is dissipated and available to power X-ray flares in powerful global eruptions, with most of our models showing an energy release in the range of the most common X-ray outbursts, ≲10 43 erg. Such events occur when significant energy builds up while deeply buried in the dipole magnetosphere. Less energetic outbursts likely precede powerful flares, due to intermittent instabilities and confined eruptions of a continuously twisting flux tube. Upon reaching a critical state, global eruptions produce the necessary Poynting-flux-dominated outflows required by models prescribing the fast radio burst production in the magnetar wind—for example, via relativistic magnetic reconnection or shocks.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2007

Multi-Scale Modeling of Magnetospheric Reconnection

One of the major challenges in modeling the magnetospheric magnetic reconnection is to quantify the interaction between large-scale global magnetospheric dynamics and microphysical processes in diffusion regions near reconnection sites. There is still considerable debate as to what degree microphysical processes on kinetic scales affect the global evolution and how important it is to substitute numerical dissipation and/or ad hoc anomalous resistivity by a physically motivated model of dissipation. Comparative studies of magnetic reconnection in small scale geometries demonstrated that MHD simulations that included non-ideal processes in terms of a resistive term $\eta J$ did not produce the fast reconnection rates observed in kinetic simulations. For a broad range of physical parameters in collisionless magnetospheric plasma, the primary mechanism controlling the dissipation in the vicinity of the reconnection site is non-gyrotropic effects with spatial scales comparable with the particle Larmor radius. We utilize the global MHD code BATSRUS and incorporate nongyrotropic effects in diffusion regions in terms of corrections to the induction equation. We developed an algorithm to search for magnetotail reconnection sites, specifically where the magnetic field components perpendicular to the local current direction approaches zero and form an X-type configuration. Spatial scales of the diffusion region and magnitude of the reconnection electric field are calculated selfconsistently using MHD plasma and field parameters in the vicinity of the reconnection site. The location of the reconnection sites is updated during the simulations. To clarify the role of nongyrotropic effects in diffusion region on the global magnetospheric dynamic we perform simulations with steady southward IMF driving of the magnetosphere. Ideal MHD simulations with magnetic reconnection supported by numerical resistivity produce steady configuration with almost stationary near-earth neutral line (NENL). Simulations with non-gyrotropic corrections demonstrate dynamic quasi-periodic response to the steady driving condition. The loading/unloading cycle in non-gyrotropic MHD results has a non-stationary reconnection site in the magnetotail, with the retreating during the stretching phase and then a new NENL forming in the resulting thin plasma sheet. We expect that this model will lead to improved representations of space weather event in the magnetosphere.

Kuznetsova, M. M.↗
NASA NTRS2005

The Flow of Plasma in the Solar-Terrestrial Environment

The overall goal of our NASA theory research is to trace the flow of mass, momentum, and energy through the magnetosphere-ionosphere-atmosphere system taking into account the coupling, time delays, and feedback mechanisms that are characteristic of the system. Our approach is to model the magnetosphere-ionosphere-atmosphere (M-I-A) system in a self-consistent quantitative manner using unique global models that allow us to study the coupling between the different regions on a range of spatial and temporal scales. The uniqueness of our global models stems from their high spatial and temporal resolutions, the physical processes included, and the numerical techniques employed. Currently, we have time-dependent global models of the ionosphere, thermosphere, polar wind, plasmasphere, and electrodynamics. It is now becoming clear that a significant fraction of the flow of mass, momentum, and energy in the M-I-A system occurs on relatively small spatial scales. Therefore, an important aspect of our NASA Theory program concerns the effect that mesoscale (100-l000 km) density structures have on the macroscopic flows in the ionosphere, thermosphere, and polar wind. The structures can be created either by structured magnetospheric inputs (i.e., structured electric field, precipitation, or Birkeland current patterns) or by time variations of these inputs due to geomagnetic storms and substorms. Some of the mesoscale structures of interest include sun-aligned polar cap arcs, propagating plasma patches, traveling convection vortices, subauroral ion drift (SAID) channels, gravity waves, and the polar hole.

Schunk, Robert W.↗
NASA NTRS2020

A Double Disturbed Lunar Plasma Wake

Under nominal solar wind conditions, a tenuous wake forms downstream of the lunar nightside. However, the lunar plasma environment undergoes a transformation as the Moon passes through the Earth's magnetotail, with hot subsonic plasma causing the wake structure to disappear. We investigate the lunar wake response during a passing coronal mass ejection (CME) on March 8, 2012 while crossing the Earth's magnetotail using both a magnetohydrodynamic (MHD) model of the terrestrial magnetosphere and a three-dimensional hybrid plasma model of the lunar wake. The CME arrives at 1 AU around 10:30 UT and its impact is first detected inside the geomagnetic tail after 11:10 UT by the Acceleration, Reconnection, Turbulence and Electrodynamics of the Moon's Interaction with the Sun (THEMIS-ARTEMIS) satellites in lunar orbit. A global magnetospheric MHD simulation using Wind data for upstream conditions with the OpenGGCM model reveals the magnetosheath compression to the lunar position from 11:20–12:00 UT, accompanied by multiple flux rope or plasmoid-like features developing and propagating tailward. MHD results support plasma changes observed by the THEMIS-ARTEMIS satellites. Lunar-scale simulations using the Amitis hybrid code show a short and misaligned plasma wake during the Moon's brief entry into the magnetosheath at 11:20 UT, with plasma expansion into the void being aided by the higher plasma temperatures. Sharply accelerated flow speed and a compressed magnetic field lead to an enhanced electric field in the lunar wake capable of generating sudden changes to the nightside near-surface electric potential.

A. P. Rasca↗
NASA NTRS1987

MHD aspects of magnetotail dynamics

The equations of ideal MHD represent the conservation of mass, momentum, and energy; thus, MHD must, in some broad sense, provide a correct description of large-scale structure of the earth's magnetosphere. However, steady-state, ideal MHD (the simplest description) allows only a static interaction between the solar wind and the magnetosphere. Hence an MHD theory of magnetospheric dynamics must involve unsteady flow turbulence, such as might be generated by Kelvin-Helmholtz turbulence at the magnetopause, and dissipative transport, such as reconnection, which violates ideal MHD. This brief review summarizes recent theoretical efforts to explore and extend the early MHD magnetospheric models of Axford and Petschek. Analytic calculations show that steady convection in the magnetotail is nearly impossible. Global magnetospheric simulations exhibit substorm reconnection phenomena for southward IMF and predict a novel convection pattern for northward IMF. Reconnection simulations show plasmoid formation and field-aligned current flows in the tail. The Kelvin-Helmholtz instability can provide a strong anomalous viscosity which can drive the viscous convection system.

Coroniti, F. V.↗
NASA NTRS2017

The Substorm Cycle as Reproduced by Global MHD Models

Recently, Gordeev et al. (2015) suggested a method to test global MHD models against statistical empirical data. They showed that four community-available global MHD models supported by the Community Coordinated Modeling Center (CCMC) produce a reasonable agreement with reality for those key parameters (the magnetospheric size, magnetic field, and pressure) that are directly related to the large-scale equilibria in the outer magnetosphere. Based on the same set of simulation runs, here we investigate how the models reproduce the global loading-unloading cycle. We found that in terms of global magnetic flux transport, three examined CCMC models display systematically different response to idealized2 h north then 2 h south interplanetary magnetic field (IMF) Bz variation. The LFM model shows a depressed return convection and high loading rate during the growth phase as well as enhanced return convection and high unloading rate during the expansion phase, with the amount of loaded unloaded magnetotail flux and the growth phase duration being the closest to their observed empirical values during isolated substorms. Two other models exhibit drastically different behavior. In the BATS-R-US model the plasma sheet convection shows a smooth transition to the steady convection regime after the IMF southward turning. In the Open GGCM a weak plasma sheet convection has comparable intensities during both the growth phase and the following slow unloading phase. We also demonstrate potential technical problem in the publicly available simulations which is related to post processing interpolation and could affect the accuracy of magnetic field tracing and of other related procedures.

Gordeev, E.↗
NASA NTRS1981

A magnetospheric signature of some F layer positive storms

Calculations of electron density distributions in the global thermosphere-ionosphere system perturbed by high-latitude thermospheric heating are presented which indicate a link between the heating and magnetospheric plasma disturbances near the equator. The calculations were made using a self-consistent model of the global sunlit thermosphere-ionosphere system describing the evolution of equatorial plasma disturbances. The heat input is found to cause electron density enhancements that propagate along magnetic field lines from the F2 maximum over mid-latitudes to the equator in the magnetosphere and which correspond to the positive phase of an F layer storm. The positive phase is shown to be generated by the induction of equatorward winds that raise the mid-latitude F layer through momentum transfer from neutral atoms to ionospheric ions, which ions pull electrons with them. Model results are used to identify plasma signatures of equatorward winds and an intensified magnetospheric electric field in Explorer 45 and Arial 4 measurements taken during the positive phase of an F layer storm.

Miller, N. J.↗
NASA NTRS1987

Magnetic field and current structures in the magnetosphere of Uranus

Voyager 2 magnetic field and plasma data are compared with theoretical model calculations for the magnetosphere of Uranus to derive a global picture from the limited set of measurements. The results suggest that Voyager 2 entered the Uranian magnetosphere during a relaxation phase which followed a compression of the entire magnetosphere. The plasma beta values in the Uranian neutral sheet were found to be smallar by a factor of 3-8 compared to corresponding beta values in the average terrestrial neutral sheet. The excellent agreement found between observed and calculated magnetic tail lobe field strengths indicates that the Uranian magnetosphere reaches the state of quasi-static (i.e., slowly time-dependent MHD equilibrium).

Voigt, G.-H.↗
NASA NTRS2012

Solar Wind Charge Exchange During Geomagnetic Storms

On March 31st. 2001, a coronal mass ejection pushed the subsolar magnetopause to the vicinity of geosynchronous orbit at 6.6 RE. The NASA/GSFC Community Coordinated Modeling Center (CCMe) employed a global magnetohydrodynamic (MHD) model to simulate the solar wind-magnetosphere interaction during the peak of this geomagnetic storm. Robertson et aL then modeled the expected 50ft X-ray emission due to solar wind charge exchange with geocoronal neutrals in the dayside cusp and magnetosheath. The locations of the bow shock, magnetopause and cusps were clearly evident in their simulations. Another geomagnetic storm took place on July 14, 2000 (Bastille Day). We again modeled X-ray emission due to solar wind charge exchange, but this time as observed from a moving spacecraft. This paper discusses the impact of spacecraft location on observed X-ray emission and the degree to which the locations of the bow shock and magnetopause can be detected in images.

Robertson, Ina P.↗
NASA NTRS2008

Global Circulation and Impact of Plasmaspheric Plumes

We report results from the global circulation model of Lyon, Fedder, and Mobarry with an embedded model of the inner magnetosphere including the plasmasphere. The combination is used to initiate large numbers of representative protons on the geosynchronous orbit L shell, to assign particle weightings, to track their: subsequent trajectories in the 3D fields. This permits us to study the global circulation of plasmaspheric plumes and to compare these with Polar observations from the dayside magnetopause region . A range of events is studied from an isolated period of SBz in the solar wind,to a large storm sequence. We consider effects on circulating plasma reaching the dayside reconnection X-line, the population of the plasma sheet with ionospheric protons and the generation of ring current pressure from this source, compared with solar wind, polar wind, and auroral wind sources. We find that the transient plasmaspheric plume source is large in terms of total fluence, but of modest proportions in terms of contribution to the ring current. Implications of this and other results for improved space weather modeling and prediction will be discussed.

Moore, Thomas E.↗
NASA NTRS2010

Dynamics of Ring Current and Electric Fields in the Inner Magnetosphere During Disturbed Periods: CRCM-BATS-R-US Coupled Model

We present simulation results from a one-way coupled global MHD model (Block-Adaptive-Tree Solar-Wind Roe-Type Upwind Scheme, BATS-R-US) and kinetic ring current models (Comprehensive Ring Current Model, CRCM, and Fok Ring Current, FokRC). The BATS-R-US provides the CRCM/FokRC with magnetic field information and plasma density/temperature at the polar CRCM/FokRC boundary. The CRCM uses an electric potential from the BATS-R-US ionospheric solver at the polar CRCM boundary in order to calculate the electric field pattern consistent with the CRCM pressure distribution. The FokRC electric field potential is taken from BATS-R-US ionospheric solver everywhere in the modeled region, and the effect of Region II currents is neglected. We show that for an idealized case with southward-northward-southward Bz IMF turning, CRCM-BATS-R-US reproduces well known features of inner magnetosphere electrodynamics: strong/weak convection under the southward/northward Bz; electric field shielding/overshielding/penetration effects; an injection during the substorm development; Subauroral Ion Drift or Polarization Jet (SAID/PJ) signature in the dusk sector. Furthermore, we find for the idealized case that SAID/PJ forms during the substorm growth phase, and that substorm injection has its own structure of field-aligned currents which resembles a substorm current wedge. For an actual event (12 August 2000 storm), we calculate ENA emissions and compare with Imager for Magnetopause-to-Aurora Global Exploration/High Energy Neutral Atom data. The CRCM-BATS-R-US reproduces both the global morphology of ring current and the fine structure of ring current injection. The FokRC-BATS-R-US shows the effect of a realistic description of Region II currents in ring current-MHD coupled models.

Buzulukova, N.↗