Search NASA⌕ Search

SEARCH · Search NASA

Results for “plasma global 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 145 records · Page 8

On the penetration of the solar wind into the cometary ionosphere

Questions of interaction between the cometary ionosphere and the solar wind are considered, taking into account present steady-state models of the global flow pattern. The stability of the ionopause and the feasibility of its penetration by the solar wind magnetic field and plasma are found to be central factors in connection with a consideration of the flow patterns. An investigation is conducted regarding a magnetic field penetration of the ionopause. It is concluded that a 'magnetic barrier' is expected to form ahead of the cometary ionopause. The stability of the ionopause is studied in some detail, taking into consideration a four-component fluid (neutrals, ions, protons, and electrons). The present study leads to the conclusion that any deep penetration of the solar wind magnetic field into the cometary ionosphere should be due to the Kelvin-Helmholtz instability of the ionopause, rather than due to the continuous diffusion process recently proposed by Marochnik (1982).

Ershkovich, A. I.↗

Magnetofluid Simulations of the Global Solar Wind Including Pickup Ions and Turbulence Modeling

I will describe a three-dimensional magnetohydrodynamic model of the solar wind that takes into account turbulent heating of the wind by velocity and magnetic fluctuations as well as a variety of effects produced by interstellar pickup protons. In this report, the interstellar pickup protons are treated as one fluid and the protons and electrons are treated together as a second fluid. The model equations include a Reynolds decomposition of the plasma velocity and magnetic field into mean and fluctuating quantities, as well as energy transfer from interstellar pickup protons to solar wind protons that results in the deceleration of the solar wind. The model is used to simulate the global steady-state structure of the solar wind in the region from 0.3 to 100 AU. Where possible, the model is compared with Voyager data. Initial results from generalization to a three-fluid model is described elsewhere in this session.

Goldstein, Melvyn L.↗

Hot relativistic winds and the Crab Nebula

Efforts to formulate a self-consistent model of pulsar magnetospheres which links the particle source near the pulsar to the outflowing relativistic wind and couples the wind to the surrounding nebula are reviewed. The use of a relativistic MHD wind is recommended to account for global photon emission and the invisibility of the method of plasma transport. Consideration of a magnetic monopole relativistic wind due to an axially symmetric aligned rotator is combined with calculations of the initial velocity of the wind to show that the flow velocity in such a model will never exceed Mach 1. Extending the solution to the case of a hot relativistic wind at supersonic speeds is noted to yield results consistent with observations of the Crab Nebula

Fujimura, F. S.↗

Multifluid Simulations of the Global Solar Wind Including Pickup Ions and Turbulence Modeling

I will describe a three-dimensional magnetohydrodynamic model of the solar wind that takes into account turbulent heating of the wind by velocity and magnetic fluctuations as well as a variety of effects produced by interstellar pickup protons. The interstellar pickup protons are treated in the model as one fluid and the protons and electrons are treated together as a second fluid. The model equations include a Reynolds decomposition of the plasma velocity and magnetic field into mean and fluctuating quantities, as well as energy transfer from interstellar pickup protons to solar wind protons that results in the deceleration of the solar wind. The model is used to simulate the global steady-state structure of the solar wind in the region from 0.3 to 100 AU. The simulation assumes that the background magnetic field on the Sun is either a dipole (aligned or tilted with respect to the solar rotation axis) or one that is deduced from solar magnetograms.

Goldstein, Melvyn L.↗

Global electric field determination in the Earth's outer magnetosphere using charged particles

Although many properties of the Earth's magnetosphere have been measured and quantified in the past 30 years since it was discovered, one fundamental (for a zeroeth order magnetohydrodynamic (MHD) equilibrium) measurement was made infrequently and with poor spatial coverage: the global electric field. This oversight is in part due to the difficulty of measuring a plasma electric field, and in part due to the difficulty of measuring a plasma electric field, and in part due to the neglect of theorists. However, there is renewed interest in the convection electric field, since it has been realized that it is vital for understanding many aspects of the magnetosphere: the global MHD equilibrium, reconnection rates, Region 2 Birkeland currents, magnetosphere-ionosphere coupling, ring current and radiation belt transport, substorm injections, acceleration mechanisms, etc. Unfortunately the standard experimental methods have not been able to synthesize a global field (excepting the pioneering work of McIlwain's geostationary models), and we are left with an overly simplistic theoretical field, the Volland-Stern electric field mode. Again, single point measurements of the plasma pause were used to infer the appropriate amplitudes of the model, parameterized by Kp (Maynard & Chen, JGR 1975). Although this result was never intended to be the definitive electric field model, it has gone nearly unchanged for 15 years. However, the data sets being taken today require a great deal more accuracy than can be provided by the Volland-Stern model. Nor has the variability of the electric field shielding been properly addressed, although effects of penetrating magnetospheric electric fields has been seen in mid- and low-latitude ionospheric data sets. The growing interests in substorm dynamics also requires a much better assessment of the electric fields responsible for particle injections. Thus, we proposed and developed algorithms for extracting electric fields from particle data taken in the earth's magnetosphere. As a test of the effectiveness of these techniques, we analyzed data taken by the AMPTE/CCE spacecraft in equatorial orbit between 1984-1988. Some analytic tools had to be developed before construction of computer algorithms, and they are discussed.

Eastman, T.↗

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↗

Modelling of auroral electrodynamical processes: Magnetosphere to mesosphere

Research conducted on auroral electrodynamic coupling between the magnetosphere and ionosphere-atmosphere in support of the development of a global scale kinetic plasma theory is reviewed. Topics covered include electric potential structure in the evening sector; morning and dayside auroras; auroral plasma formation; electrodynamic coupling with the thermosphere; and auroral electron interaction with the atmosphere.

Chiu, Y. T.↗

Mshpy23: A User-Friendly, Parameterized Model of Magnetosheath Conditions

Lunar Environment heliospheric X-ray Imager (LEXI) and Solar wind – Magnetosphere - Ionosphere Link Explorer (SMILE) will observe magnetosheath and its boundary motion in soft X-rays for understanding magnetopause reconnection modes under varioussolar wind conditions after their respective launches in 2024 and 2025. Magnetosheath conditions, namely, plasma density, velocity, and temperature, are key parameters for predicting and analyzing soft X-ray images from the LEXI and SMILE missions. We developed a user-friendly model of magnetosheath that parameterizes number density, velocity, temperature, and magnetic field by utilizing the global Magnetohydrodynamics (MHD) model as well as the pre-existing gas-dynamic and analytic models. Using this parameterized magnetosheath model, scientists can easily reconstruct expected soft X-ray images and utilize them for analysis of observed images of LEXI and SMILE without simulating the complicated global magnetosphere models. First, we created an MHD based magnetosheath model by running a total of 14 OpenGGCM global MHD simulations under 7 solar wind densities (1, 5, 10, 15, 20, 25, and 30cm−3) and 2 interplanetary magnetic field BZ components (± 4nT), and then parameterizing the results in new magnetosheath conditions. We compared the magnetosheath model result with THEMIS statistical data and it showed good agreement with a weighted Pearson correlation coefficient greater than 0.77, especially for plasma density and plasma velocity. Second, we compiled a suite of magnetosheath models incorporating previous magnetosheath models (gas-dynamic, analytic), and did two case studies to test the performance. The MHD based model was comparable to or better than the previous models while providing self consistency among the magnetosheath parameters. Third, we constructed a tool to calculate a soft X-ray image from any given vantage point, which can support the planning and data analysis of the aforementioned LEXI and SMILE missions. A release of the code has been uploaded to a Github repository.

Magnetosheath↗

Modeling of the Convection and Interaction of Ring Current, Plasmaspheric and Plasma Sheet Plasmas in the Inner Magnetosphere

Distinctive sources of ions reside in the plasmasphere, plasmasheet, and ring current regions at discrete energies constitute the major plasma populations in the inner/middle magnetosphere. They contribute to the electrodynamics of the ionosphere-magnetosphere system as important carriers of the global current system, in triggering; geomagnetic storm and substorms, as well as critical components of plasma instabilities such as reconnection and Kelvin-Helmholtz instability at the magnetospheric boundaries. Our preliminary analysis of in-situ measurements shoves the complexity of the plasmas pitch angle distributions at particularly the cold and warm plasmas, vary dramatically at different local times and radial distances from the Earth in response to changes in solar wind condition and Dst index. Using an MHD-ring current coupled code, we model the convection and interaction of cold, warm and energetic ions of plasmaspheric, plasmasheet, and ring current origins in the inner magnetosphere. We compare our simulation results with in-situ and remotely sensed measurements from recent instrumentation on Geotail, Cluster, THEMIS, and TWINS spacecraft.

Fok, Mei-Ching↗

The Distant Magnetotail Under Long Duration, Very Northward IMF Conditions: October 22-24, 2003

A unique 32 hour interval of very northward Interplanetary Magnetic Field (IMF) on October 22-24, 2003 created a exceptionally thick cold dense magnetotail plasma sheet, a small polar cap and accompanying small tail lobe. These features were detected by the Cluster DMSP and FAST spacecraft and modeled by a global simulation as described in papers by Oieroset et al. (2005) and Li et al. (2005). During the same interval the Wind spacecraft was passing through the center of the magnetotail about 130 Re downstream of Earth. Wind results will be described that reveal a very unusual magnetotail characterized by (1) continual tailward flow of 200-400 km/s with densities in the range 0.2-3/cc, both of whch are clearly less than those expected in the magnetosheath, (2) a mostly northward Bz but with a predominant Bx field component with sign reversals indicating crossings between the two hemispheres of the tail, and (3) velocity waves superposed on the downstream flow with peak-to-peak amplitudes of 100 to 200 km/s, periods of 10 to 20 minutes and clockwise polarization. Low altitude DMSP and Fast measurements reveal an auroral oval with enhanced latitudinal thickness and a small polar cap filled with structured precipitzting electrons and few ions. A new global MHD simulation of the event exhibits a highly elliptical tail of diminished cross-section at 130 Re with major axis aligned with the northward IMF. The tail current sheet also tends to be aligned in a north-south direction with the two tail hemispheres to the east and west with their polarities depending on prior history of the IMF. The simulation appears to be consistent with many, but not all, of the observations. High latitude cusp reconnection and subsequent downtail flow of closed field lines may explain the tail structure, but the waves are more likely due to the Kelvin-Helmholtz instability often thought to occur during northward IMF conditions.

Fairfield, Donald H.↗

Inferring Lower Boundary Driving Conditions Using Vector Magnetic Field Observations

Low-beta coronal MHD simulations of realistic CME events require the detailed specification of the magnetic fields, velocities, densities, temperatures, etc., in the low corona. Presently, the most accurate estimates of solar vector magnetic fields are made in the high-beta photosphere. Several techniques have been developed that provide accurate estimates of the associated photospheric plasma velocities such as the Differential Affine Velocity Estimator for Vector Magnetograms and the Poloidal/Toroidal Decomposition. Nominally, these velocities are consistent with the evolution of the radial magnetic field. To evolve the tangential magnetic field radial gradients must be specified. In addition to estimating the photospheric vector magnetic and velocity fields, a further challenge involves incorporating these fields into an MHD simulation. The simulation boundary must be driven, consistent with the numerical boundary equations, with the goal of accurately reproducing the observed magnetic fields and estimated velocities at some height within the simulation. Even if this goal is achieved, many unanswered questions remain. How can the photospheric magnetic fields and velocities be propagated to the low corona through the transition region? At what cadence must we observe the photosphere to realistically simulate the corona? How do we model the magnetic fields and plasma velocities in the quiet Sun? How sensitive are the solutions to other unknowns that must be specified, such as the global solar magnetic field, and the photospheric temperature and density?

Schuck, Peter W.↗

Propagation of whistler waves driven by fine structured ion beams in the magnetotail

In a previous paper, which examined the propagation of low-frequency whistler waves generated by ion beams in the Earth's plasma sheet boundary layer (PSBL), it was found that whistler waves driven in the PSBL are focused toward the central plasma sheet due to the global magnetotail inhomogeneities; this finding may help explain the observations of magnetic noise bursts in the tail (Burinskaya et al., 1993). In this paper the same phenomenon is examined, but this time a much more realistic model is used for the ion beam in the PSBL. While the PSBL has been modeled as a solid, homogeneous ion beams with a width of one Earth radius, observations and theoretical considerations have shown that PSBL ion beams actually have a decreasing velocity profile toward the plasma sheet and that the density of the beams within the PSBL can vary locally. We consider again the propagation and generation of electromagnetic waves but in the presence of fine structured ion beams in the PSBL. Our results show that whistler waves, generated quasi-parallel to the background magnetic field, can be trapped locally within small spatial regions where the ion beam density is enhanced compared to the density of the adjacent PSBL region. Wave spectra and nonlinear saturation mechanisms are discussed.

Burinskaya, T.↗

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.↗

Multispacecraft Observations and Modeling of the 22/23 June 2015 Geomagnetic Storm

The magnetic storm of 22-23 June 2015 was one of the largest in the current solar cycle. We present in situ observations from the Magnetospheric Multiscale Mission (MMS) and the Van Allen Probes (VAP) in the magnetotail, field-aligned currents from AMPERE (Active Magnetosphere and Planetary Electrodynamics Response), and ionospheric flow data from Defense Meteorological Satellite Program (DMSP). Our real-time space weather alert system sent out a "red alert," correctly predicting Kp indices greater than 8. We show strong outflow of ionospheric oxygen, dipolarizations in the MMS magnetometer data, and dropouts in the particle fluxes seen by the MMS Fast Plasma Instrument suite. At ionospheric altitudes, the AMPERE data show highly variable currents exceeding 20 MA. We present numerical simulations with the Block Adaptive Tree-Solarwind - Roe - Upwind Scheme (BATS-R-US) global magnetohydrodynamic model linked with the Rice Convection Model. The model predicted the magnitude of the dipolarizations, and varying polar cap convection patterns, which were confirmed by DMSP measurements.

Reiff, P. H.↗

Global MHD Simulations of Space Plasma Environments: Heliosphere, Comets, Magnetospheres of Plants and Satellites

Magnetohydrodynamics (MHD) provides an approximate description of a great variety of processes in space physics. Accurate numerical solutions of the MHD equations are still a challenge, but in the past decade a number of robust methods have appeared. Once these techniques made the direct solution of MHD equations feasible, a number of global three-dimensional models were designed and applied to many space physics objects. The range of these objects is truly astonishing, including active galactic nuclei, the heliosphere, the solar corona, and the solar wind interaction with planets, satellites, and comets. Outside the realm of space physics, MHD theory has been applied to such diverse problems as laboratory plasmas and electromagnetic casting of liquid metals. In this paper we present a broad spectrum of models of different phenomena in space science developed in the recent years at the University of Michigan. Although the physical systems addressed by these models are different, they all use the MHD equations as a unifying basis.

Kabin, K.↗

Lunar electrical conductivity, permeability,and temperature from Apollo magnetometer experiments

Magnetometers were deployed at four Apollo sites on the moon to measure remanent and induced lunar magnetic fields. Measurements from this network of instruments were used to calculate the electrical conductivity, temperature, magnetic permeability, and iron abundance of the lunar interior. Global lunar fields due to eddy currents, induced in the lunar interior by magnetic transients, were analyzed to calculate and electrical conductivity profile for the moon, and those profiles were used to calculate the lunar temperature for an assumed lunar material of olivine. Simultaneous measurements by magnetometers on the lunar surface and in orbit around the moon were use to construct a whole-moon hysteresis curve, from which the global lunar magnetic permeability is determined. Total iron abundance (sum of iron in the ferromagnetic and paramagnetic states) was calculated for two assumed compositional models of the lunar interior. Other lunar models with an iron core and with a shallow iron-rich layer also discussed in light of the measured global lunar permeability. Simultaneous magnetic field and solar plasma pressure measurements show that the remanent fields at the Apollo 12 and 16 sites interact with, and are compressed by, the solar wind. Velocities and thicknesses of the earth's magnetopause and bow shock were also estimated from simultaneous magnetometer measurements.

Dyal, P.↗

Plasma IMS Composition Measurements for Europa and the Other Galilean Moons

NASA and ESA are planning the joint Europa Jupiter System Mission (EJSM) to the Jupiter system with specific emphasis to Europa and Ganymede, respectively. The Japanese Space Agency is also planning an orbiter mission to explore Jupiter's magnetosphere and the Galilean satellites. For NASA's Jupiter Europa Orbiter (JEO) we are developing the 3D Ion Mass Spectrometer (IMS) with two main goals which can also be applied to the other Galilean moons, 1) measure the plasma interaction between Europa and Jupiter's magnetosphere and 2) infer the 4 pi surface composition to trace elemental and significant isotopic levels. The first goal supports the magnetometer (MAG) measurements, primarily directed at detection of Europa's sub-surface ocean, while the second gives information about transfer of material between the Galilean moons, and between the moon surfaces and subsurface layers putatively including oceans. The measurement of the interactions for all the Galilean moons can be used to trace the in situ ion measurements of pickup ions back to either Europa's or Ganymede's surface from the respectively orbiting spacecraft. The IMS instrument, being developed under NASA's Astrobiology Instrument Development Program, would maximally achieve plasma measurement requirements for JEO and EJSM while moving forward our knowledge of Jupiter system composition and source processes to far higher levels than previously envisaged. The composition of the global surfaces of Europa and Ganymede can be inferred from the measurement of ejected neutrals and pick-up ions using at minimum an in situ payload including MAG and IMS also fully capable of meeting Level 1 mission requirements for ocean detection and survey. Elemental and isotopic analysis of potentially extruded oceanic materials at the moon surfaces would further support the ocean objectives. These measurements should be made from a polar orbiting spacecraft about Europa or Ganymede at height 100 km. The ejecta produced by sputtering of the surfaces of Europa and Ganymede has been shown to be representative of the surface composition. Level 2 science on surface geology and composition can then be further enhanced by addition of the following: 3D Ion Neutral Mass Spectrometer (INNS), 3D plasma electron spectrometer (ELS), and hot plasma energetic particle instrument. The measurement approach is to alternate between times measuring pickup ions and times measuring plasma and magnetic field parameters along the spacecraft trajectory. By measuring the pickup ion energy, arrival direction and mass-per-charge, the ion can be traced back along the ejection trajectory to the approximate area of origin if the 3-D electric field and magnetic field are known. In situ observations of plasma flow velocities and vector magnetic fields can be used to determine the local convective electric field (E = -VXB) along the spacecraft trajectory. By combining this information with models of the magnetospheric interaction with Europa, one can generate 3D maps of the electric and magnetic field and compute the trajectories of the pickup ions back to the surface or exospheric points of origin. In the case of Ganymede there is the additional complexity of its own internal dipole magnetic field, while Io's volcanic activity introduces the complexity of a highly structured denser atmosphere. Callisto with its less globally extended exosphere will have a simpler interaction than for Europa (i.e., more like our moon). We will discuss these differences in light of the above proposed technique. Finally, the INNS observations and neutral exosphere models are needed to estimate production rates of pickup ions. The hot plasma measurements are needed to correct for sputtering rates which can be time dependent and electron plasma observations for electron impact ionization rates. Instrument characteristics, field-of-view requirements, modes of operation and effects of radiation on instrument functionality will be discussed.

Sittler, Edward↗

A theoretical and empirical study of the response of the high latitude thermosphere to the sense of the 'Y' component of the interplanetary magnetic field

Patterns of magnetospheric energetic plasma precipitation as a function of the Y component of the Interplanetary Magnetic Field (IMF) are studied. The development of a three-dimensional, time-dependent global thermospheric model using a polar conversion electric field with a dependence on the Y component of the IMF to evaluate thermospheric wind circulation is examined. Thermospheric wind data from the ISEE-3 satellite, Dynamics Explorer-2 satellite, and a ground-based Fabry-Perot interferometer in Kiruna, Sweden, collected on December 1, 2, 6, 25, 1981 and February 12, 13, 1982 are described. The observed data and simulations of polar thermospheric winds are compared. In the Northern Hemisphere a strong antisunward ion flow on the dawn side of the geomagnetic polar cap is observed when the BY is positive, and the flow is detected on the dusk side when the BY is negative. It is concluded that the strength and direction of the IMF directly control the transfer of solar wind momentum and energy to the high latitude thermosphere.

Rees, D.↗