skip to main content

This content will become publicly available on January 1, 2023

Title: Numerical Modeling of Suprathermal Electron Transport in the Solar Wind: Effects of Whistler Turbulence with a Full Diffusion Tensor
Abstract The electron VDF in the solar wind consists of a Maxwellian core, a suprathermal halo, a field-aligned component strahl, and an energetic superhalo that deviates from the equilibrium. Whistler wave turbulence is thought to resonantly scatter the observed electron velocity distribution. Wave–particle interactions that contribute to Whistler wave turbulence are introduced into a Fokker–Planck kinetic transport equation that describes the interaction between the suprathermal electrons and the Whistler waves. A recent numerical approach for solving the Fokker–Planck kinetic transport equation has been extended to include a full diffusion tensor. Application of the extended numerical approach to the transport of solar wind suprathermal electrons influenced by Whistler wave turbulence is presented. Comparison and analysis of the numerical results with observations and diagonal-only model results are made. The off-diagonal terms in the diffusion tensor act to depress effects caused by the diagonal terms. The role of the diffusion coefficient on the electron heat flux is discussed.
; ;
Award ID(s):
Publication Date:
Journal Name:
The Astrophysical Journal
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract A self-consistent hybrid model of standing and moving striations was developed for low-current DC discharges in noble gases. We introduced the concept of surface diffusion in phase space (r,u) (where u denotes the electron kinetic energy) described by a tensor diffusion in the nonlocal Fokker-Planck kinetic equation for electrons in the collisional plasma. Electrons diffuse along surfaces of constant total energy ε=u-eφ(r) between energy jumps in inelastic collisions with atoms. Numerical solutions of the 1d1u kinetic equation for electrons were obtained by two methods and coupled to ion transport and Poisson solver. We studied the dynamics of striation formation in Townsend and glow discharges in Argon gas at low discharge currents using a two-level excitation-ionization model and a “full-chemistry” model, which includes stepwise and Penning ionization. Standing striations appeared in Townsend and glow discharges at low currents, and moving striations were obtained for the discharge currents exceeding a critical value. These waves originate at the anode and propagate towards the cathode. We have seen two types of moving striations with the 2-level and full-chemistry models, which resemble the s and p striations previously observed in the experiments. Simulations indicate that processes in the anode region could control moving striationsmore »in the positive column plasma. The developed model helps clarify the nature of standing and moving striations in DC discharges of noble gases at low discharge currents and low gas pressures.« less
  2. null (Ed.)
    ABSTRACT We present a kinetic stability analysis of the solar wind electron distribution function consisting of the Maxwellian core and the magnetic-field aligned strahl, a superthermal electron beam propagating away from the sun. We use an electron strahl distribution function obtained as a solution of a weakly collisional drift-kinetic equation, representative of a strahl affected by Coulomb collisions but unadulterated by possible broadening from turbulence. This distribution function is essentially non-Maxwellian and varies with the heliospheric distance. The stability analysis is performed with the Vlasov–Maxwell linear solver leopard. We find that depending on the heliospheric distance, the core-strahl electron distribution becomes unstable with respect to sunward-propagating kinetic-Alfvén, magnetosonic, and whistler modes, in a broad range of propagation angles. The wavenumbers of the unstable modes are close to the ion inertial scales, and the radial distances at which the instabilities first appear are on the order of 1 au. However, we have not detected any instabilities driven by resonant wave interactions with the superthermal strahl electrons. Instead, the observed instabilities are triggered by a relative drift between the electron and ion cores necessary to maintain zero electric current in the solar wind frame (ion frame). Contrary to strahl distributions modelled by shiftedmore »Maxwellians, the electron strahl obtained as a solution of the kinetic equation is stable. Our results are consistent with the previous studies based on a more restricted solution for the electron strahl.« less
  3. The relative velocities and positions of monodisperse high-inertia particle pairs in isotropic turbulence are studied using direct numerical simulations (DNS), as well as Langevin simulations (LS) based on a probability density function (PDF) kinetic model for pair relative motion. In a prior study (Rani et al. , J. Fluid Mech. , vol. 756, 2014, pp. 870–902), the authors developed a stochastic theory that involved deriving closures in the limit of high Stokes number for the diffusivity tensor in the PDF equation for monodisperse particle pairs. The diffusivity contained the time integral of the Eulerian two-time correlation of fluid relative velocities seen by pairs that are nearly stationary. The two-time correlation was analytically resolved through the approximation that the temporal change in the fluid relative velocities seen by a pair occurs principally due to the advection of smaller eddies past the pair by large-scale eddies. Accordingly, two diffusivity expressions were obtained based on whether the pair centre of mass remained fixed during flow time scales, or moved in response to integral-scale eddies. In the current study, a quantitative analysis of the (Rani et al. 2014) stochastic theory is performed through a comparison of the pair statistics obtained using LS with those from DNS.more »LS consist of evolving the Langevin equations for pair separation and relative velocity, which is statistically equivalent to solving the classical Fokker–Planck form of the pair PDF equation. Langevin simulations of particle-pair dispersion were performed using three closure forms of the diffusivity – i.e. the one containing the time integral of the Eulerian two-time correlation of the seen fluid relative velocities and the two analytical diffusivity expressions. In the first closure form, the two-time correlation was computed using DNS of forced isotropic turbulence laden with stationary particles. The two analytical closure forms have the advantage that they can be evaluated using a model for the turbulence energy spectrum that closely matched the DNS spectrum. The three diffusivities are analysed to quantify the effects of the approximations made in deriving them. Pair relative-motion statistics obtained from the three sets of Langevin simulations are compared with the results from the DNS of (moving) particle-laden forced isotropic turbulence for $St_{\unicode[STIX]{x1D702}}=10,20,40,80$ and $Re_{\unicode[STIX]{x1D706}}=76,131$ . Here, $St_{\unicode[STIX]{x1D702}}$ is the particle Stokes number based on the Kolmogorov time scale and $Re_{\unicode[STIX]{x1D706}}$  is the Taylor micro-scale Reynolds number. Statistics such as the radial distribution function (RDF), the variance and kurtosis of particle-pair relative velocities and the particle collision kernel were computed using both Langevin and DNS runs, and compared. The RDFs from the stochastic runs were in good agreement with those from the DNS. Also computed were the PDFs $\unicode[STIX]{x1D6FA}(U|r)$ and $\unicode[STIX]{x1D6FA}(U_{r}|r)$ of relative velocity $U$ and of the radial component of relative velocity $U_{r}$ respectively, both PDFs conditioned on separation $r$ . The first closure form, involving the Eulerian two-time correlation of fluid relative velocities, showed the best agreement with the DNS results for the PDFs.« less
  4. In this note we establish hypocoercivity and exponential relaxation to the Maxwellian for a class of kinetic Fokker-Planck-Alignment equations arising in the studies of collective behavior. Unlike previously known results in this direction that focus on convergence near Maxwellian, our result is global for hydrodynamically dense flocks, which has several consequences. In particular, if communication is long-range, the convergence is unconditional. If communication is local then all nearly aligned flocks quantified by smallness of the Fisher information relax to the Maxwellian. In the latter case the class of initial data is stable under the vanishing noise limit, i.e. it reduces to a non-trivial and natural class of traveling wave solutions to the noiseless Vlasov-Alignment equation.

    The main novelty in our approach is the adaptation of a mollified Favre filtration of the macroscopic momentum into the communication protocol. Such filtration has been used previously in large eddy simulations of compressible turbulence and its new variant appeared in the proof of the Onsager conjecture for inhomogeneous Navier-Stokes system. A rigorous treatment of well-posedness for smooth solutions is provided. Lastly, we prove that in the limit of strong noise and local alignment solutions to the Fokker-Planck-Alignment equation Maxwellialize to solutions of themore »macroscopic hydrodynamic system with the isothermal pressure.

    « less
  5. Over three decades of in-situ observations illustrate that the Kelvin–Helmholtz (KH) instability driven by the sheared flow between the magnetosheath and magnetospheric plasma often occurs on the magnetopause of Earth and other planets under various interplanetary magnetic field (IMF) conditions. It has been well demonstrated that the KH instability plays an important role for energy, momentum, and mass transport during the solar-wind-magnetosphere coupling process. Particularly, the KH instability is an important mechanism to trigger secondary small scale (i.e., often kinetic-scale) physical processes, such as magnetic reconnection, kinetic Alfvén waves, ion-acoustic waves, and turbulence, providing the bridge for the coupling of cross scale physical processes. From the simulation perspective, to fully investigate the role of the KH instability on the cross-scale process requires a numerical modeling that can describe the physical scales from a few Earth radii to a few ion (even electron) inertial lengths in three dimensions, which is often computationally expensive. Thus, different simulation methods are required to explore physical processes on different length scales, and cross validate the physical processes which occur on the overlapping length scales. Test particle simulation provides such a bridge to connect the MHD scale to the kinetic scale. This study applies different testmore »particle approaches and cross validates the different results against one another to investigate the behavior of different ion species (i.e., H+ and O+), which include particle distributions, mixing and heating. It shows that the ion transport rate is about 10 25  particles/s, and mixing diffusion coefficient is about 10 10  m 2  s −1 regardless of the ion species. Magnetic field lines change their topology via the magnetic reconnection process driven by the three-dimensional KH instability, connecting two flux tubes with different temperature, which eventually causes anisotropic temperature in the newly reconnected flux.« less