skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Stream-aligned Magnetohydrodynamics for Solar Wind Simulations
We present a reduced magnetohydrodynamic (MHD) mathematical model describing the dynamical behavior of highly conducting plasmas with frozen-in magnetic fields, constrained by the assumption that there exists a frame of reference, where the magnetic field vector, B, is aligned with the plasma velocity vector, u, at each point. We call this solution “stream-aligned MHD” (SA-MHD). Within the framework of this model, the electric field, E = − u × B ≡ 0, in the induction equation vanishes identically and so does the electromagnetic energy flux (Poynting flux), E × B ≡ 0, in the energy equation. At the same time, the force effect from the magnetic field on the plasma motion (the Ampère force) is fully taken into account in the momentum equation. Any steady-state solution of the proposed model is a legitimate solution of the full MHD system of equations. However, the converse statement is not true: in an arbitrary steady-state magnetic field, the electric field does not have to vanish identically (its curl has to, though). Specifically, realistic three-dimensional solutions for the steady-state (“ambient”) solar atmosphere in the form of so-called Parker spirals can be efficiently generated within the stream aligned MHD (SA-MHD) with no loss in generality.  more » « less
Award ID(s):
1836821
PAR ID:
10321308
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
The astrophysical journal
Volume:
926
ISSN:
2041-8213
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We used the stream-aligned magnetohydrodynamics (SA-MHD) model to simulate Carrington rotation 2210, which contains Parker Solar Probe’s (PSP) first perihelion at 36.5Ron 2018 November 6, to provide context to the in situ PSP observations by FIELDS and SWEAP. The SA-MHD model aligns the magnetic field with the velocity vector at each point, thereby allowing for clear connectivity between the spacecraft and the source regions on the Sun, without unphysical magnetic field structures. During this Carrington rotation, two stream interaction regions (SIRs) form, due to the deep solar minimum. We include the energy partitioning of the parallel and perpendicular ions and the isotropic electrons to investigate the temperature anisotropy through the compression regions to better understand the wave energy amplification and proton thermal energy partitioning in a global context. Overall, we found good agreement in all in situ plasma parameters between the SA-MHD results and the observations at PSP, STEREO-A, and Earth, including at PSP’s perihelion and through the compression region of the SIRs. In the typical solar wind, the parallel proton temperature is preferentially heated, except in the SIR, where there is an enhancement in the perpendicular proton temperature. This is further showcased in the ion cyclotron relaxation time, which shows a distinct decrease through the SIR compression regions. This work demonstrates the success of the Alfvén wave turbulence theory for predicting interplanetary magnetic turbulence levels, while self-consistently reproducing solar wind speeds, densities, and overall temperatures, including at small heliocentric distances and through SIR compression regions. 
    more » « less
  2. Abstract Substorm‐type evolution of the Earth's magnetosphere is investigated by mining more than two decades (1995–2017) of spaceborne magnetometer data from multiple missions including the first two years (2016‐2017) of the Magnetospheric MultiScale mission. This investigation reveals interesting features of plasma evolution distinct from ideal magnetohydrodynamics (MHD) behavior: X‐lines, thin current sheets, and regions with the tailward gradient of the equatorial magnetic fieldBz. X‐lines are found to form mainly beyond 20RE, but for strong driving, with the solar wind electric field exceeding ∼5mV/m, they may come closer. For substorms with weaker driving, X‐lines may be preceded by redistribution of the magnetic flux in the tailwardBzgradient regions, similar to the magnetic flux release instability discovered earlier in PIC and MHD simulations as a precursor mechanism of the reconnection onset. Current sheets in the growth phase may be as thin as 0.2RE, comparable to the thermal ions gyroradius, and at the same time, as long as 15RE. Such an aspect ratio is inconsistent with the isotropic force balance for observed magnetic field configurations. These findings can help resolve kinetic mechanisms of substorm dipolarizations and adjust kinetic generalizations of global MHD models of the magnetosphere. They can also guide and complement microscale analysis of nonideal effects. 
    more » « less
  3. Springer (Ed.)
    We analyze non-perturbatively the one-dimensional Schrödinger equation describing the emission of electrons from a model metal surface by a classical oscillating electric field. Placing the metal in the half-space x 0, the Schrödinger equation of the system is i∂t ψ = − 2 1 ∂x 2 ψ + (x)(U − E x cos ωt)ψ, t > 0, x ∈ R, where (x) is the Heaviside function and U > 0 is the effective confining potential (we choose units so that m = e = = 1). The amplitude E of the external electric field and the frequency ω are arbitrary. We prove existence and uniqueness of classical solutions of the Schrödinger equation for general initial conditions ψ(x, 0) = f (x), x ∈ R. When 2the initial condition is in L the evolution is unitary and the wave function goes to zero at any fixed x as t → ∞. To show this we prove a RAGE type theorem and show that the discrete spectrum of the quasienergy operator is empty. To obtain positive electron current we consider non-L 2 initial conditions containing an incoming beam from the left. The beam is partially reflected and partially transmitted for all t > 0. For these initial conditions we show that the solution approaches in the large t limit a periodic state that satisfies an infinite set of equations formally derived, under the assumption that the solution is periodic by Faisal et al. (Phys Rev A 72:023412, 2005). Due to a number of pathological features of the Hamiltonian (among which unboundedness in the physical as well as the spatial Fourier domain) the existing methods to prove such results do not apply, and we introduce new, more general ones. The actual solution exhibits a very complex behavior, as seen both analytically and numerically. It shows a steep increase in the current as the frequency passes a threshold value ω = ωc , with ωc depending on the strength of the electric field. For small E, ωc represents the threshold in the classical photoelectric effect, as described by Einstein’s theory. 
    more » « less
  4. null (Ed.)
    Variations of vertical atmospheric electric field E z have been attributed mainly to meteorological processes. On the other hand, the theory of electromagnetic waves in the atmosphere, between the bottom ionosphere and earth’s surface, predicts two modes, magnetic H (TE) and electric E (TH) modes, where the E-mode has a vertical electric field component, E z . Past attempts to find signatures of ULF (periods from fractions to tens of minutes) disturbances in E z gave contradictory results. Recently, study of ULF disturbances of atmospheric electric field became feasible thanks to project GLOCAEM, which united stations with 1 sec measurements of potential gradient. These data enable us to address the long-standing problem of the coupling between atmospheric electricity and space weather disturbances at ULF time scales. Also, we have reexamined results of earlier balloon-born electric field and ground magnetic field measurements in Antarctica. Transmission of storm sudden commencement (SSC) impulses to lower latitudes was often interpreted as excitation of the electric TH 0 mode, instantly propagating along the ionosphere–ground waveguide. According to this theoretical estimate, even a weak magnetic signature of the E-mode ∼1 nT must be accompanied by a burst of E z well exceeding the atmospheric potential gradient. We have examined simultaneous records of magnetometers and electric field-mills during >50 SSC events in 2007–2019 in search for signatures of E-mode. However, the observed E z disturbance never exceeded background fluctuations ∼10 V/m, much less than expected for the TH 0 mode. We constructed a model of the electromagnetic ULF response to an oscillating magnetospheric field-aligned current incident onto the realistic ionosphere and atmosphere. The model is based on numerical solution of the full-wave equations in the atmospheric-ionospheric collisional plasma, using parameters that were reconstructed using the IRI model. We have calculated the vertical and horizontal distributions of magnetic and electric fields of both H- and E-modes excited by magnetospheric field-aligned currents. The model predicts that the excitation rate of the E-mode by magnetospheric disturbances is low, so only a weak E z response with a magnitude of ∼several V/m will be produced by ∼100 nT geomagnetic disturbance. However, at balloon heights (∼30 km), electric field of the E-mode becomes dominating. Predicted amplitudes of horizontal electric field in the atmosphere induced by Pc5 pulsations and travelling convection vortices, about tens of mV/m, are in good agreement with balloon electric field and ground magnetometer observations. 
    more » « less
  5. The validity of conventional Ohm’s law is tested in the context of a rapidly evolving quark–gluon plasma produced in heavy-ion collisions. Here, we discuss the electromagnetic response using an analytical solution in kinetic theory. As conjectured previously, after switching on an electric field in a nonexpanding plasma, the time-dependent current is given by J(t)=(1−e−t/τ0)σ0E, where τ0 is the transport relaxation time and σ0 is the steady-state electrical conductivity. Such an incomplete electromagnetic response reduces the efficiency of the magnetic flux trapping in the quark–gluon plasma, and may prevent the observation of the chiral magnetic effect. Here, we extend the study to the case of a rapidly expanding plasma. We find that the decreasing temperature and the increasing transport relaxation time have opposite effects on the electromagnetic response. While the former suppresses the time-dependent conductivity, the latter enhances it. 
    more » « less