Numerical simulation of an experimental analogue of a planetary magnetosphere
Not Available
SEARCH · Search NASA
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.
Not Available
Abstract not provided.
By adopting a 2D time-dependent wave code, we investigate how mode-converted waves at the Ion-Ion Hybrid (IIH) resonance and compressional waves propagate in 2D density structures with a wide range of field-aligned wavenumbers to background magnetic fields. The simulation results show that the mode-converted waves have continuous bands across the field line consistent with previous numerical studies. These waves also have harmonic structures in frequency domain and are localized in the field-aligned heavy ion density well. Lastly, our results thus emphasize the importance of a field-aligned heavy ion density structure for ultra-low frequency wave propagation, and suggest that IIH waves can be localized in different locations along the field line.
The Active Magnetosphere and Planetary Electrodynamics Response Experiment uses magnetic field data from the Iridium constellation to derive the global Birkeland current distribution every 10 min. We examine cases in which the interplanetary magnetic field (IMF) rotated from northward to southward resulting in onsets of the Birkeland currents. Dayside Region 1/2 currents, totaling ~25% of the final current, appear within 20 min of the IMF southward turning and remain steady. In the onset of nightside currents occurs 40 to 70 min after the dayside currents appear. Afterwards, the currents intensify at dawn, dusk, and on the dayside, yielding a fully formed Region 1/2 system ~30 min after the nightside onset. Our results imply that the dayside Birkeland currents are driven by magnetopause reconnection, and the remainder of the system forms as magnetospheric return flows start and progress sunward, ultimately closing the Dungey convection cycle.
Magnetic reconnection—the topological rearrangement of magnetic field—underlies many explosive phenomena across a wide range of natural and laboratory plasmas.3 It plays a pivotal role in electron and ion heating, particle acceleration to high energies, energy transport, and self-organization. Reconnection can have a complex relationship with turbulence at both large and small scales, leading to various effects that are only beginning to be understood. In heliophysics, magnetic reconnection plays a key role in solar flares, coronal mass ejections, coronal heating, solar wind dissipation, the interaction of interplanetary plasma with magnetospheres, dynamics of planetary magnetospheres such as magnetic substorms, and the heliospheric boundary with the interstellar medium. Here, the magnetic reconnection is integral to the solar and planetary dynamo processes.
Global interactions of energetic ions with magnetoplasmas and neutral gases lie at the core of many space and laboratory plasma phenomena ranging from solar wind entry into and transport within planetary magnetospheres and exospheres to fast-ion driven instabilities in fusion devices to astrophysics-in-lab experiments. The ability of computational models to properly account for physical effects that underlie such interactions, namely ion kinetic, ion cyclotron, Hall, collisional and ionization processes is important for the success and planning of experimental research in plasma physics. Understanding the physics of energetic ions, in particular their nonlinear resonance interactions with Alfvén waves, is central to improving the heating performance of magnetically confined plasmas for future energy generation. Fluid models are not adequate for high-beta plasmas as they cannot fully capture ion kinetic and cyclotron physics (e.g., ion behavior in the presence of magnetic nulls, shock structures, plasma interpenetration, etc.). Recent results from global reconnection simulations show that even in a MHD-like regime there may be significant differences between kinetic and MHD simulations. Therefore, kinetic modeling becomes essential for meeting modern day challenges in plasma physics. The hybrid approximation is an intermediate approximation between the fluid and fully kinetic approximations. It eliminates light waves, removes the electron inertial temporal and spatial scales from the problem and enables full-orbit ion kinetics. As a result, hybrid codes have become effective tools for exploring ion-scale driven phenomena associated with ion beams, shocks, reconnection and turbulence that control the large-scale behavior of laboratory and space magnetoplasmas. A number of numerical issues, however, make three-dimensional (3D) large-scale hybrid simulations of inhomogeneous magnetized plasmas prohibitively expensive or even impossible. To resolve these difficulties we have developed a novel Event-driven Multiscale Asynchronous Parallel Simulation (EMAPS) technology that replaces time stepping with self-adaptive update events. Local calculations are carried out only on an “as needed basis”. EMAPS (i) guarantees accurate and stable processing of physical variables in time accurate simulations, and (ii) eliminates unnecessary computation. Applying EMAPS to the hybrid model has resulted in the development of a unique parallel code, dimension-independent (compile-time-configurable) HYPERS (Hybrid Parallel Event-Resolved Simulator) that scales to hundreds of thousands of parallel processors. HYPERS advances electromagnetic fields and particles asynchronously on time scales determined by local physical laws and mesh properties. To achieve high computational accuracy in complex device geometries, HYPERS employs high-fidelity Cartesian grids with masked conductive cells. The HYPERS model includes multiple ion species, energy and momentum conserving ion-ion collisions, and provides a number of approximations for plasma resistivity and vacuum regions. Both local and periodic boundary conditions are allowed. The HYPERS solver preserves zero divergence of magnetic field. The project has demonstrated HYPERS capabilities on a number of applications of interest to fusion and astrophysical plasma physics applications listed below. 1. Theta-pinch formation of FRCs The formation, spontaneous spin-up, and stability of theta-pinch formed field-reversed configurations have been studied self-consistently in 3D. The end-to-end hybrid simulations reveal poloidal profiles of implosion-driven fast toroidal plasma rotation and demonstrate three discharge regimes as a function of experimental parameters: the decaying stable configuration, the tilt unstable configuration, and the nonlinear evolution of a fast growing tearing mode. 2. FRC collisions with magnetic mirrors Interactions of fast plasma streams and objects with magnetic obstacles (dipoles, mirrors, etc) lie at the core of many space and laboratory plasma phenomena ranging from magnetoshells and solar wind interactions with planetary magnetospheres to compact fusion plasmas. HYPERS simulations are compared with data from the MSX experiment (LANL) that focuses on the physics of magnetized collisionless shocks through the acceleration and subsequent stagnation of FRC plasmoids against a strong magnetic mirrors and flux-conserving boundaries. 3. Exploding magnetoplasmas Results from hybrid simulations of two experiments at the LAPD and Nevada Terawatt Facility are discussed where short-pulse lasers are used to ablate solid targets to produce plasmas that expand across external magnetic fields. The first simulation recreates flutelike density striations observed at the leading edge of a carbon plasma and predicts an early destruction of the magnetic cavity in agreement with experimental evidence. In the second simulation a polyethylene target is ablated into a mixture of protons and carbon ions. A mechanism is demonstrated that allows protons to penetrate the magnetic field in the form of a collimated flow. The results are compared to experimental data and single-fluid MHD simulations. The EMAPS framework has the potential for wide application in many other engineering and scientific fields, such as climate models, biological systems, electronic devices, seismic events, oil reservation simulators that all involve advancing solutions of partial differential equations in time where the rate of activity can be adapted widely over the spatial domain depending on locally space/time phenomena (“events”).
Bow shocks are ubiquitous astrophysical phenomena resulting from the supersonic passage of an object through a gas. Recently, pre-transit absorption in UV metal transitions of the hot Jupiter (HJ) exoplanets HD 189733b and WASP12-b have been interpreted as being caused by material compressed in a planetary bow shock. Here we present a robust detection of a time-resolved pre-transit, as well as in-transit absorption signature around the HJ exoplanet HD 189733b using high spectral resolution observations of several hydrogen Balmer lines. The line shape of the pre-transit feature and the shape of the timeseries absorption provide the strongest constraints on the morphology and physical characteristics of extended structures around an exoplanet. The in-transit measurements confirm the previous exospheric Hα detection, although the absorption depth measured here is ∼50% lower. The pre-transit absorption feature occurs 125 minutes before the predicted optical transit, a projected linear distance from the planet to the stellar disk of 7.2 R{sub p}. The absorption strength observed in the Balmer lines indicates an optically thick, but physically small, geometry. We model this signal as the early ingress of a planetary bow shock. If the bow shock is mediated by a planetary magnetosphere, the large standoff distance derived from the model suggests a large planetary magnetic field strength of B{sub eq} = 28 G. Better knowledge of exoplanet magnetic field strengths is crucial to understanding the role these fields play in planetary evolution and the potential development of life on planets in the habitable zone.
We investigate solitary magnetic structures in the foreshock of Mars and Earth. The structures exhibit pulse-like magnetic field and density enhancements along with plasma heating and local solar wind (SW) ion reflection. The structures at Mars resemble the foreshock structures developed from ultra-low-frequency electromagnetic waves gyro-resonant with SW ions at Earth, and propagate toward the magnetosphere. We perform fully kinetic simulations to reproduce the solitary structures and the nonlinear evolution of ion distribution functions, illustrating their resonance with SW ions. The structures present self-induced foreshock turbulence that can have space weather effects. Our results advance the fundamental understanding of how SW interacts with planetary magnetospheres, and have potential impact on the current picture of planet-origin ion escape.
Linear and nonlinear electrostatic ion acoustic waves in a weakly relativistic magnetorotating plasma in the presence of non-Maxwellian electrons and warm ions have been examined. The system under consideration has yielded two solutions, namely, the fast and slow acoustic modes which have been observed to depend on the streaming velocity, ion to electron temperature ratio, and the nonthermality parameter of the non-Maxwellian electrons. Using the multiple time scale analysis, we have derived the three dimensional nonlinear Zakharov–Kuznetsov equation and also presented its solution. Both compressive and rarefactive solitary structures have been found in consonance with the satellite observations. It has been observed that although the linear dispersion relation gives both fast and slow ion acoustic waves, the solitary structures form only for the fast acoustic mode. The dependence of the characteristics of the solitary structures on several plasma parameters has also been explored. The present investigation may be beneficial to understanding the rotating plasma environments such as those found in the planetary magnetospheres of Saturn and Jupiter.
We investigate magnetic reconnection in systems simultaneously containing asymmetric (anti-parallel) magnetic fields, asymmetric plasma densities and temperatures, and arbitrary in-plane bulk flow of plasma in the upstream regions. Such configurations are common in the high-latitudes of Earth's magnetopause and in tokamaks. We investigate the convection speed of the X-line, the scaling of the reconnection rate, and the condition for which the flow suppresses reconnection as a function of upstream flow speeds. We use two-dimensional particle-in-cell simulations to capture the mixing of plasma in the outflow regions better than is possible in fluid modeling. We perform simulations with asymmetric magnetic fields, simulations with asymmetric densities, and simulations with magnetopause-like parameters where both are asymmetric. For flow speeds below the predicted cutoff velocity, we find good scaling agreement with the theory presented in Doss et al. [J. Geophys. Res. 120, 7748 (2015)]. Applications to planetary magnetospheres, tokamaks, and the solar wind are discussed.
Investigation of the lunar crustal magnetic anomalies offers a comprehensive long-term data set of observations of small-scale magnetic fields and their interaction with the solar wind. In this paper a review of the observations of lunar mini-magnetospheres is compared quantifiably with theoretical kinetic-scale plasma physics and 3D particle-in-cell simulations. The aim of this paper is to provide a complete picture of all the aspects of the phenomena and to show how the observations from all the different and international missions interrelate. The analysis shows that the simulations are consistent with the formation of miniature (smaller than the ion Larmor orbit) collisionless shocks and miniature magnetospheric cavities, which has not been demonstrated previously. The simulations reproduce the finesse and form of the differential proton patterns that are believed to be responsible for the creation of both the “lunar swirls” and “dark lanes.” Using a mature plasma physics code like OSIRIS allows us, for the first time, to make a side-by-side comparison between model and space observations. This is shown for all of the key plasma parameters observed to date by spacecraft, including the spectral imaging data of the lunar swirls. The analysis of miniature magnetic structures offers insight into multi-scale mechanisms and kinetic-scale aspects of planetary magnetospheres.
The classic Bernstein waves may be intimately related to banded emissions detected in laboratory plasmas, terrestrial, and other planetary magnetospheres. However, the customary discussion of the Bernstein wave is based upon isotropic thermal velocity distribution function. In order to understand how such waves can be excited, one needs an emission mechanism, i.e., an instability. In non-relativistic collision-less plasmas, the only known Bernstein wave instability is that associated with a cold perpendicular velocity ring distribution function. However, cold ring distribution is highly idealized. The present Brief Communication generalizes the cold ring distribution model to include thermal spread, so that the Bernstein-ring instability is described by a more realistic electron distribution function, with which the stabilization by thermal spread associated with the ring distribution is demonstrated. The present findings imply that the excitation of Bernstein waves requires a sufficiently high perpendicular velocity gradient associated with the electron distribution function.
In this white paper, we emphasize the importance of in situ ion composition measurements provided by plasma mass spectrometers in planetary magnetospheres (Uranus and Neptune) aimed for future planetary missions. We present our current understanding on magnetospheres of Uranus and Neptune that are largely based on Plasma Science Investigation (PLS, faraday cups) and Low Energy Charged Particle (LECP) instrument on Voyager 2 and what measurements need to be made to improve it. We emphasize that the future missions equipped with plasma mass spectrometers must visit these magnetospheres again to advance our current knowledge. Finally, this white paper summarizes key science questions that needs to be addressed which all of them requires a direct ion composition measurements from a plasma spectrometer.
We study the magnetospheric structure and the ionospheric Joule Heating of planets orbiting M-dwarf stars in the habitable zone using a set of magnetohydrodynamic models. The stellar wind solution is used to drive a model for the planetary magnetosphere, which is coupled with a model for the planetary ionosphere. Our simulations reveal that the space environment around close-in habitable planets is extreme, and the stellar wind plasma conditions change from sub- to super-Alfvénic along the planetary orbit. As a result, the magnetospheric structure changes dramatically with a bow shock forming in the super-Alfvénic sectors, while no bow shock forms in the sub-Alfvénic sectors. The planets reside most of the time in the sub-Alfvénic sectors with poor atmospheric protection. A significant amount of Joule Heating is provided at the top of the atmosphere as a result of the intense stellar wind. For the steady-state solution, the heating is about 0.1%-3% of the total incoming stellar irradiation, and it is enhanced by 50% for the time-dependent case. The significant Joule Heating obtained here should be considered in models for the atmospheres of habitable planets in terms of the thickness of the atmosphere, the top-side temperature and density, the boundary conditions for the atmospheric pressure, and particle radiation and transport. Here we assume constant ionospheric Pedersen conductance similar to that of the Earth. The conductance could be greater due to the intense EUV radiation leading to smaller heating rates. We plan to quantify the ionospheric conductance in future study.
We report magnetospheres are a ubiquitous feature of magnetized bodies embedded in a plasma flow. While large planetary magnetospheres have been studied for decades by spacecraft, ion-scale “mini” magnetospheres can provide a unique environment to study kinetic-scale, collisionless plasma physics in the laboratory to help validate models of larger systems. In this work, we present preliminary experiments of ion-scale magnetospheres performed on a unique high-repetition-rate platform developed for the Large Plasma Device at the University of California, Los Angeles. The experiments utilize a high-repetition-rate laser to drive a fast plasma flow into a pulsed dipole magnetic field embedded in a uniform magnetized background plasma. 2D maps of the magnetic field with high spatial and temporal resolution are measured with magnetic flux probes to examine the evolution of magnetosphere and current density structures for a range of dipole and upstream parameters. The results are further compared to 2D particle-in-cell simulations to identify key observational signatures of the kinetic-scale structures and dynamics of the laser-driven plasma. We find that distinct 2D kinetic-scale magnetopause and diamagnetic current structures are formed at higher dipole moments, and their locations are consistent with predictions based on pressure balances and energy conservation.
Uranus possesses the most extreme magnetic and rotational geometry in the solar system, resulting in a uniquely dynamic and asymmetric interaction between its magnetosphere and the solar wind. Here we investigate the diurnal evolution of the Uranian bow shock (BS) at equinox using global multifluid magnetohydrodynamic simulations constrained by Voyager 2 observations. The BS structure and variability are quantified using the stand‐off distance, terminator distance, flaring parameter, and BS parameter, enabling a systematic assessment of its global asymmetry over one planetary rotation. Our results show that the BS exhibits pronounced diurnal asymmetry that is primarily controlled by planetary rotation. Even under steady upstream solar wind conditions, the BS undergoes periodic expansion and contraction, reflecting rotation‐driven reconfiguration of the planetary magnetospheric topology. This diurnal modulation represents a remarkably large asymmetry in planetary BS geometry within the solar system. These results establish planetary rotation as the intrinsic primary driver of Uranus' BS dynamics at equinox and provide essential constraints for the design and formulation of future space missions to the ice giants. These findings also offer a unique framework for understanding shock physics in the abundant population of ice‐giant exoplanets.
Abstract The rate of magnetic reconnection is of the utmost importance in a variety of processes because it controls, for example, the rate energy is released in solar flares, the speed of the Dungey convection cycle in Earth’s magnetosphere, and the energy release rate in harmful geomagnetic substorms. It is known from numerical simulations and satellite observations that the rate is approximately 0.1 in normalized units, but despite years of effort, a full theoretical prediction has not been obtained. Here, we present a first-principles theory for the reconnection rate in non-relativistic electron-ion collisionless plasmas, and show that the same prediction explains why Sweet-Parker reconnection is considerably slower. The key consideration of this analysis is the pressure at the reconnection site (i.e., the x-line). We show that the Hall electromagnetic fields in antiparallel reconnection cause an energy void, equivalently a pressure depletion, at the x-line, so the reconnection exhaust opens out, enabling the fast rate of 0.1. If the energy can reach the x-line to replenish the pressure, the exhaust does not open out. In addition to heliospheric applications, these results are expected to impact reconnection studies in planetary magnetospheres, magnetically confined fusion devices, and astrophysical plasmas.
In the Earth’s dayside reconnection boundary layer, whistler mode waves coincide with magnetic field openings and the formation of the resultant anisotropic electrons. Depending on the energy range of anisotropic electrons, whistlers can grow at frequencies in the upper and/or lower band. Observations show that whistler mode waves modulate Langmuir wave amplitude as they propagate toward the X line. Observations of whistler mode wave phase and Langmuir waves packets, as well as coincident electron measurements, reveal that whistler mode waves can accelerate electrons via Landau resonance at locations where E || is antiparallel to the wave propagation direction. The accelerated electrons produce localized beams, which subsequently drive the periodically modulated Langmuir waves. The close association of those two wave modes reveals the microscale electron dynamics in the exhaust region, and the proposed mechanism could potentially be applied to explain the modulation events observed in planetary magnetospheres and in the solar wind.