Engineering Papers⌕ Search

Engineering topics

Zank, G. P.

Publications and source records attributed to Zank, G. P..

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Flux Ropes, Turbulence, and Collisionless Perpendicular Shock Waves: High Plasma Beta Case

With the onset of solar maximum and the expected increased prevalence of interplanetary shock waves, Parker Solar Probe is likely to observe numerous shocks in the next few years. An outstanding question that has received surprisingly little attention has been how turbulence interacts with collisionless shock waves. Turbulence in the supersonic solar wind is described frequently as a superposition of a majority 2D and a minority slab component. Here we formulate a collisional perpendicular shock-turbulence transmission problem in a way that enables investigation of the interaction and transmission of quasi-perpendicular fluctuations such as magnetic flux ropes/islands and vortices as well as entropy and acoustic modes in the large plasma beta regime. We focus on the transmission of an upstream spectrum of these modes, finding that the downstream spectral amplitude is typically increased significantly (a factor of 10 or more), and that the upstream spectral index of the inertial range, and indeed the general spectral shape, is unchanged for the downstream magnetic variance, kinetic energy, and density variance. A comparison of the theoretically predicted downstream magnetic variance, kinetic energy, and density variance spectra with those observed at 1, 5, and 84 au by Wind, Ulysses, and Voyager 2 shows excellent agreement. The overall theoretically predicted characteristics of the transmission of turbulence across shocks observed in the solar wind appear to be largely consistent with recent observational studies by Pitňa et al. and Borovsky.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Kinetic entropy-based measures of distribution function non-Maxwellianity: theory and simulations

We investigate kinetic entropy-based measures of the non-Maxwellianity of distribution functions in plasmas, i.e. entropy-based measures of the departure of a local distribution function from an associated Maxwellian distribution function with the same density, bulk flow and temperature as the local distribution. First, we consider a form previously employed by Kaufmann & Paterson (J. Geophys. Res., vol. 114, 2009, A00D04), assessing its properties and deriving equivalent forms. To provide a quantitative understanding of it, we derive analytical expressions for three common non-Maxwellian plasma distribution functions. We show that there are undesirable features of this non-Maxwellianity measure including that it can diverge in various physical limits and elucidate the reason for the divergence. We then introduce a new kinetic entropy-based non-Maxwellianity measure based on the velocity-space kinetic entropy density, which has a meaningful physical interpretation and does not diverge. We use collisionless particle-in-cell simulations of two-dimensional anti-parallel magnetic reconnection to assess the kinetic entropy-based non-Maxwellianity measures. We show that regions of non-zero non-Maxwellianity are linked to kinetic processes occurring during magnetic reconnection. We also show the simulated non-Maxwellianity agrees reasonably well with predictions for distributions resembling those calculated analytically. These results can be important for applications, as non-Maxwellianity can be used to identify regions of kinetic-scale physics or increased dissipation in plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Interstellar Mapping and Acceleration Probe (IMAP): A New NASA Mission

The Interstellar Mapping and Acceleration Probe (IMAP) is a revolutionary mission that simultaneously investigates two of the most important overarching issues in Heliophysics today: the acceleration of energetic particles and interaction of the solar wind with the local interstellar medium. While seemingly disparate, these are intimately coupled because particles accelerated in the inner heliosphere play critical roles in the outer heliospheric interaction. Selected by NASA in 2018, IMAP is planned to launch in 2024. The IMAP spacecraft is a simple sun-pointed spinner in orbit about the Sun-Earth L1 point. IMAP’s ten instruments provide a complete and synergistic set of observations to simultaneously dissect the particle injection and acceleration processes at 1 AU while remotely probing the global heliospheric interaction and its response to particle populations generated by these processes. In situ at 1 AU, IMAP provides detailed observations of solar wind electrons and ions; suprathermal, pickup, and energetic ions; and the interplanetary magnetic field. For the outer heliosphere interaction, IMAP provides advanced global observations of the remote plasma and energetic ions over a broad energy range via energetic neutral atom imaging, and precise observations of interstellar neutral atoms penetrating the heliosphere. Complementary observations of interstellar dust and the ultraviolet glow of interstellar neutrals further deepen the physical understanding from IMAP. IMAP also continuously broadcasts vital real-time space weather observations. Finally, IMAP engages the broader Heliophysics community through a variety of innovative opportunities. This paper summarizes the IMAP mission at the start of Phase A development.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2018

Interstellar Mapping and Acceleration Probe (IMAP): A New NASA Mission

The Interstellar Mapping and Acceleration Probe (IMAP) is a revolutionary mission that simultaneously investigates two of the most important overarching issues in Heliophysics today: the acceleration of energetic particles and interaction of the solar wind with the local interstellar medium. While seemingly disparate, these are intimately coupled because particles accelerated in the inner heliosphere play critical roles in the outer heliospheric interaction. Selected by NASA in 2018, IMAP is planned to launch in 2024. The IMAP spacecraft is a simple sun-pointed spinner in orbit about the Sun-Earth L1 point. IMAP's ten instruments provide a complete and synergistic set of observations to simultaneously dissect the particle injection and acceleration processes at 1 AU while remotely probing the global heliospheric interaction and its response to particle populations generated by these processes. In situ at 1 AU, IMAP provides detailed observations of solar wind electrons and ions; suprathermal, pickup, and energetic ions; and the interplanetary magnetic field. For the outer heliosphere interaction, IMAP provides advanced global observations of the remote plasma and energetic ions over a broad energy range via energetic neutral atom imaging, and precise observations of interstellar neutral atoms penetrating the heliosphere. Complementary observations of interstellar dust and the ultraviolet glow of interstellar neutrals further deepen the physical understanding from IMAP. IMAP also continuously broadcasts vital real-time space weather observations. Finally, IMAP engages the broader Heliophysics community through a variety of innovative opportunities. This papersummarizes the IMAP mission at the start of Phase A development.

McComas, D. J.↗
DOE OSTI2017

The Origin of Compressible Magnetic Turbulence in the Very Local Interstellar Medium

Voyager 1 observed compressible magnetic turbulence in the very local interstellar medium (VLISM). We show that inner heliosheath (IHS) fast- and slow-mode waves incident on the heliopause (HP) generate VLISM fast-mode waves only that propagate into the VLISM. We suggest that this is the origin of compressible turbulence in the VLISM. We show that fast- and slow-mode waves transmitted across a tangential discontinuity such as the HP are strongly refracted on crossing the HP and subsequently propagate at highly oblique angles to the VLISM magnetic field. Thus, fast-mode waves in the VLISM contribute primarily to the compressible and not the transverse components of the VLISM fluctuating magnetic field variance 〈δ B-hat {sup 2}〉 since 〈δ B-hat {sub fz}{sup 2}〉≠0 and 〈δ B-hat {sub fx,fy}{sup 2}〉≃0. If the fast- and slow-mode waves in the IHS exhibit a Kolmogorov-like k{sup −5/3} power spectral density, as appears to be observed by Voyager 1, then the corresponding transmitted spectral density in the VLISM forms an amplified k{sup −5/3} power law. Consequently, the HP “radiates” fast-mode fluctuations into the VLISM, and the heliosphere therefore mediates the character of turbulence in the VLISM. In particular, we predict the form of the VLISM magnetic turbulence power spectral density to be a superposition of the background pristine interstellar turbulence spectrum and the fast-mode spectrum generated by the interaction of fast- and slow-mode IHS waves with the HP, i.e., a power law with an enhanced feature or “bump” corresponding to the contribution by fast-mode turbulence radiated by the HP. We note the applicability of these results to the generation of compressible turbulence in the local interstellar medium surrounding the asterospheres of stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

II. Transport of Nearly Incompressible Magnetohydrodynamic Turbulence from 1 to 75 au

The thermal plasma beta in the solar wind and the solar corona is of the order of β∼1 and β≪1. Zank et al. developed 2D and slab turbulence transport model equations of the order of β∼1 and β≪1 using nearly incompressible (NI) theory. We solve the Zank et al. NI MHD coupled turbulence transport equations for the inhomogeneous solar wind from 1 to 75 au, and compare the numerical solutions to Voyager 2 observations. We find that (1) the 2D turbulent energies are larger than the slab energies throughout the heliosphere; (2) the 2D turbulent energies decrease more slowly than the slab turbulent energies within ∼4 au, while the slab energies increase and the 2D energies flatten in the outer heliosphere; (3) the 2D normalized cross-helicity decreases faster than the slab normalized cross-helicity within ∼4 au; (4) the 2D normalized residual energy is more magnetically dominated than the slab; (5) the variance of density fluctuations decreases more rapidly than r{sup −4} within ∼10 au, and more slowly in the outer heliosphere; and (6) the observed variance in magnetic field fluctuations as a function of the thermal plasma beta is described by the two-component turbulence transport model. In summary, the NI MHD two-component Zank et al. turbulence transport model captures the behavior of the forward, backward, and total energies in the fluctuating Elsässer variables, the variance in the magnetic field, kinetic energy, and density fluctuations, the cross-helicities and residual energies, the thermal temperature and plasma beta, and the various correlation lengths.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Structure of Energetic Particle Mediated Shocks Revisited

The structure of collisionless shock waves is often modified by the presence of energetic particles that are not equilibrated with the thermal plasma (such as pickup ions [PUIs] and solar energetic particles [SEPs]). This is relevant to the inner and outer heliosphere and the Very Local Interstellar Medium (VLISM), where observations of shock waves (e.g., in the inner heliosphere) show that both the magnetic field and thermal gas pressure are less than the energetic particle component pressures. Voyager 2 observations revealed that the heliospheric termination shock (HTS) is very broad and mediated by energetic particles. PUIs and SEPs contribute both a collisionless heat flux and a higher-order viscosity. We show that the incorporation of both effects can completely determine the structure of collisionless shocks mediated by energetic ions. Since the reduced form of the PUI-mediated plasma model is structurally identical to the classical cosmic ray two-fluid model, we note that the presence of viscosity, at least formally, eliminates the need for a gas sub-shock in the classical two-fluid model, including in that regime where three are possible. By considering parameters upstream of the HTS, we show that the thermal gas remains relatively cold and the shock is mediated by PUIs. We determine the structure of the weak interstellar shock observed by Voyager 1 . We consider the inclusion of the thermal heat flux and viscosity to address the most general form of an energetic particle-thermal plasma two-fluid model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

On the Parallel and Oblique Firehose Instability in Fluid Models

A brief analysis of the proton parallel and oblique firehose instability is presented from a fluid perspective, and the results are compared with kinetic theory solutions obtained by the WHAMP code. It is shown that the classical CGL model very accurately describes the growth rate of these instabilities at sufficiently long spatial scales (small wavenumbers). The required stabilization of these instabilities at small spatial scales (high wavenumbers) naturally requires dispersive effects, and the stabilization is due to the Hall term and finite Larmor radius corrections to the pressure tensor. Even though the stabilization is not completely accurate, since at small spatial scales a relatively strong collisionless damping comes into effect, we find that the main concepts of the maximum growth rate and the stabilization of these instabilities is indeed present in the fluid description. However, there are differences that are quite pronounced when close to the firehose threshold, and that clarify the different profiles for marginally stable states with a prescribed maximum growth rate γ{sub max} in the simple fluid models considered here and the kinetic description.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Turbulent Transport in a Three-dimensional Solar Wind

Turbulence in the solar wind can play essential roles in the heating of coronal and solar wind plasma and the acceleration of the solar wind and energetic particles. Turbulence sources are not well understood and thought to be partly enhanced by interaction with the large-scale inhomogeneity of the solar wind and the interplanetary magnetic field and/or transported from the solar corona. To investigate the interaction with background inhomogeneity and the turbulence sources, we have developed a new 3D MHD model that includes the transport and dissipation of turbulence using the theoretical model of Zank et al. We solve for the temporal and spatial evolution of three moments or variables, the energy in the forward and backward fluctuating modes and the residual energy and their three corresponding correlation lengths. The transport model is coupled to our 3D model of the inhomogeneous solar wind. We present results of the coupled solar wind-turbulence model assuming a simple tilted dipole magnetic configuration that mimics solar minimum conditions, together with several comparative intermediate cases. By considering eight possible solar wind and turbulence source configurations, we show that the large-scale solar wind and IMF inhomogeneity and the strength of the turbulence sources significantly affect the distribution of turbulence in the heliosphere within 6 au. We compare the predicted turbulence distribution results from a complete solar minimum model with in situ measurements made by the Helios and Ulysses spacecraft, finding that the synthetic profiles of the turbulence intensities show reasonable agreement with observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Structure of the Heliotail from Interstellar Boundary Explorer Observations: Implications for the 11-year Solar Cycle and Pickup Ions in the Heliosheath

Interstellar Boundary Explorer ( IBEX ) measurements of energetic neutral atoms (ENAs) from the heliotail show a multi-lobe structure of ENA fluxes as a function of energy between ∼0.71 and 4.29 keV. Below ∼2 keV, there is a single structure of enhanced ENA fluxes centered near the downwind direction. Above ∼2 keV, this structure separates into two lobes, one north and one south of the solar equatorial plane. ENA flux from these two lobes can be interpreted as originating from the fast solar wind (SW) propagating through the inner heliosheath (IHS). Alternatively, a recently published model of the heliosphere suggests that the heliotail may split into a “croissant-like” shape, and that such a geometry could be responsible for the heliotail ENA feature. Here we present results from a time-dependent simulation of the heliosphere that produces a comet-like heliotail, and show that the 11-year solar cycle leads to the formation of ENA lobes with properties remarkably similar to those observed by IBEX . The ENA energy at which the north and south lobes appear suggests that the pickup ion (PUI) temperature in the slow SW of the IHS is ∼10{sup 7} K. Moreover, we demonstrate that the extinction of PUIs by charge-exchange is an essential process required to create the observed global ENA structure. While the shape and locations of the ENA lobes as a function of energy are well reproduced by PUIs that cross the termination shock, the results appear to be sensitive to the form of the distribution of PUIs injected in the IHS.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Theory and Transport of Nearly Incompressible Magnetohydrodynamic Turbulence

The theory of nearly incompressible magnetohydrodynamics (NI MHD) was developed largely in the early 1990s, together with an important extension to inhomogeneous flows in 2010. Much of the focus in the earlier work was to understand the apparent incompressibility of the solar wind and other plasma environments, and the relationship of density fluctuations to apparently incompressible manifestations of turbulence in the solar wind and interstellar medium. Further important predictions about the “dimensionality” of solar wind turbulence and its relationship to the plasma beta were made and subsequently confirmed observationally. However, despite the initial success of NI MHD in describing fluctuations in the solar wind, a detailed application to solar wind turbulence has not been undertaken. Here, we use the equations of NI MHD to describe solar wind turbulence, rewriting the NI MHD system in terms of Elsässer variables. Distinct descriptions of 2D and slab turbulence emerge naturally from the Elsässer formulation, as do the nonlinear couplings between 2D and slab components. For plasma beta order 1 or less regions, predictions for 2D and slab spectra result from the NI MHD description, and predictions for the spectral characteristics of density fluctuations can be made. We conclude by presenting a NI MHD formulation describing the transport of majority 2D and minority slab turbulence throughout the solar wind. A preliminary comparison of theory and observations is presented.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

THE INTERACTION OF TURBULENCE WITH PARALLEL AND PERPENDICULAR SHOCKS: THEORY AND OBSERVATIONS AT 1 au

Shocks are thought to be responsible for the amplification of turbulence as well as for generating turbulence throughout the heliosphere. We study the interaction of turbulence with parallel and perpendicular shock waves using the six-coupled-equation turbulence transport model of Zank et al. We model a 1D stationary shock wave using a hyperbolic tangent function and the Rankine–Hugoniot conditions for both a reduced model with four coupled equations and the full model. Eight quasi-parallel and five quasi-perpendicular events in the WIND spacecraft data sets are identified, and we compute the fluctuating magnetic and kinetic energy, the energy in forward and backward propagating modes, the total turbulent energy, the normalized residual energy, and the normalized cross helicity upstream and downstream of the observed shocks. We compare the observed fitted values upstream and downstream of the shock with numerical solutions to our model equations. The comparison shows that our theoretical results are in reasonable agreement with observations for both quasi-parallel and perpendicular shocks. We find that (1) the total turbulent energy, the energy in forward and backward propagating modes, and the normalized residual energy increase across the shock, (2) the normalized cross helicity increases or decreases across the shock, and (3) the correlation length increases upstream and downstream of the shock, and slightly flattens or decreases across the shock.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

LATITUDE, ENERGY, AND TIME VARIATIONS IN THE ENERGETIC NEUTRAL ATOM SPECTRAL INDICES MEASURED BY THE INTERSTELLAR BOUNDARY EXPLORER ( IBEX )

We investigate the latitude, energy, and temporal variations of the ∼0.5–6 keV energetic neutral atom (ENA) fluxes and spectral indices measured by the Interstellar Boundary Explorer ( IBEX ) outside of the Ribbon from 2009 to 2013. By combining all-sky maps for years 1–3 and years 4–5, we find that the >1.2 keV globally distributed (GD) ENA fluxes at all latitudes decrease by ∼25% from years 1–3 to years 4–5, but there is no change (<5%) in the corresponding spectral indices. We also show that the latitudinal profile and energy-dependence of the ≥1 keV ENA spectral indices outside the Ribbon exhibit no significant time differences between years 1–3 and 4–5, lending strong support for the notion that the highly organized, persistent energy dependence and latitudinal pattern of the ENA spectral indices during 2009–2014 are determined by that of the solar wind (SW) speed observed in the inner heliosphere during the deep solar minimum conditions of 2006–2010. In contrast, the ∼25% decrease in the >1.2 keV ENA fluxes at all latitudes from years 1–3 to years 4–5 occurs because the magnitude of the driving SW parameters (SW density or dynamic pressure) diminished significantly from 2006 to 2010. Based on the reconstructed latitudinal and temporal profiles of SW parameters from 2011 to 2014, i.e., during the rising phase of solar cycle 24, we suggest that the GD ENA fluxes in years 2014–2017, i.e., in maps 6–9, will either stabilize or increase and the latitudinal pattern and energy dependence of the corresponding spectral indices will be disrupted.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Pickup Ion Effect of the Solar Wind Interaction with the Local Interstellar Medium

Pickup ions are created when interstellar neutral atoms resonantly exchange charge with the solar wind (SW) ions, especially in the supersonic part of the wind, where they carry most of the plasma pressure. Here we present numerical simulation results of the 3D heliospheric interface treating pickup ions as a separate proton fluid. To satisfy the fundamental conservation laws, we solve the system of equations describing the flow of the mixture of electrons, thermal protons, and pickup ions. To find the density and pressure of pickup ions behind the termination shock, we employ simple boundary conditions that take into account the \emph{Voyager} observations that showed that the decrease in the kinetic energy of the mixture at the termination shock predominantly contributed to the increase in the pressure of pickup ions. We show that this model adequately describes the flow of the plasma mixture and results in a noticeable decrease in the heliosheath width.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

MODELING THE SOLAR WIND AT THE ULYSSES , VOYAGER , AND NEW HORIZONS SPACECRAFT

The outer heliosphere is a dynamic region shaped largely by the interaction between the solar wind and the interstellar medium. While interplanetary magnetic field and plasma observations by the Voyager spacecraft have significantly improved our understanding of this vast region, modeling the outer heliosphere still remains a challenge. We simulate the three-dimensional, time-dependent solar wind flow from 1 to 80 astronomical units (au), where the solar wind is assumed to be supersonic, using a two-fluid model in which protons and interstellar neutral hydrogen atoms are treated as separate fluids. We use 1 day averages of the solar wind parameters from the OMNI data set as inner boundary conditions to reproduce time-dependent effects in a simplified manner which involves interpolation in both space and time. Our model generally agrees with Ulysses data in the inner heliosphere and Voyager data in the outer heliosphere. Ultimately, we present the model solar wind parameters extracted along the trajectory of the New Horizons spacecraft. We compare our results with in situ plasma data taken between 11 and 33 au and at the closest approach to Pluto on 2015 July 14.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

COMBINING DIFFUSIVE SHOCK ACCELERATION WITH ACCELERATION BY CONTRACTING AND RECONNECTING SMALL-SCALE FLUX ROPES AT HELIOSPHERIC SHOCKS

Computational and observational evidence is accruing that heliospheric shocks, as emitters of vorticity, can produce downstream magnetic flux ropes and filaments. This led Zank et al. to investigate a new paradigm whereby energetic particle acceleration near shocks is a combination of diffusive shock acceleration (DSA) with downstream acceleration by many small-scale contracting and reconnecting (merging) flux ropes. Using a model where flux-rope acceleration involves a first-order Fermi mechanism due to the mean compression of numerous contracting flux ropes, Zank et al. provide theoretical support for observations that power-law spectra of energetic particles downstream of heliospheric shocks can be harder than predicted by DSA theory and that energetic particle intensities should peak behind shocks instead of at shocks as predicted by DSA theory. In this paper, a more extended formalism of kinetic transport theory developed by le Roux et al. is used to further explore this paradigm. We describe how second-order Fermi acceleration, related to the variance in the electromagnetic fields produced by downstream small-scale flux-rope dynamics, modifies the standard DSA model. The results show that (i) this approach can qualitatively reproduce observations of particle intensities peaking behind the shock, thus providing further support for the new paradigm, and (ii) stochastic acceleration by compressible flux ropes tends to be more efficient than incompressible flux ropes behind shocks in modifying the DSA spectrum of energetic particles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

CGL description revisited

Solar wind observational studies have emphasized that the solar wind plasma data is bounded by the mirror and firehose instabilities, and it is often believed that these instabilities are of a purely kinetic nature. The simplest fluid model that generalizes magnetohydrodynamics with anisotropic temperatures is the Chew-Goldberger-Low model (CGL). Here we briefly revisit the CGL description and discuss its (otherwise well-documented) linear firehose and mirror instability thresholds; namely that the firehose instability threshold is identical to the one found from linear kinetic theory and that the mirror threshold contains a factor of 6 error. We consider a simple higher-order fluid model with time dependent heat flux equations and show that the mirror instability threshold is correctly reproduced. We also present fully nonlinear three-dimensional simulations of freely decaying turbulence for the Hall-CGL model with isothermal electrons. The spatial resolution of these simulations is 512{sup 3} and the formation of a spectral break in magnetic and velocity field spectra around the proton inertial length is found.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗