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: Multifluid Simulation of Solar Chromospheric Turbulence and Heating Due to Thermal Farley–Buneman Instability
Abstract Models fail to reproduce observations of the coldest parts of the Sun’s atmosphere, where interactions between multiple ionized and neutral species prevent an accurate MHD representation. This paper argues that a meter-scale electrostatic plasma instability develops in these regions and causes heating. We refer to this instability as the Thermal Farley–Buneman Instability (TFBI). Using parameters from a 2.5D radiative MHD Bifrost simulation, we show that the TFBI develops in many of the colder regions in the chromosphere. This paper also presents the first multifluid simulation of the TFBI and validates this new result by demonstrating close agreement with theory during the linear regime. The simulation eventually develops turbulence, and we characterize the resulting wave-driven heating, plasma transport, and turbulent motions. These results all contend that the effects of the TFBI contribute to the discrepancies between solar observations and radiative MHD models.  more » « less
Award ID(s):
1903416 2149781
PAR ID:
10486194
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
949
Issue:
2
ISSN:
0004-637X
Format(s):
Medium: X Size: Article No. 59
Size(s):
Article No. 59
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The transport of waves and turbulence beyond the photosphere is central to the coronal heating problem. Turbulence in the quiet solar corona has been modeled on the basis of the nearly incompressible magnetohydrodynamic (NI MHD) theory to describe the transport of low-frequency turbulence in open magnetic field regions. It describes the evolution of the coupled majority quasi-2D and minority slab component, driven by the magnetic carpet and advected by a subsonic, sub-Alfvénic flow from the lower corona. In this paper, we couple the NI MHD turbulence transport model with an MHD model of the solar corona to study the heating problem in a coronal loop. In a realistic benchmark coronal loop problem, we find that a loop can be heated to ∼1.5 million K by transport and dissipation of MHD turbulence described by the NI MHD model. We also find that the majority 2D component is as important as the minority slab component in the heating of the coronal loop. We compare our coupled MHD/NI MHD model results with a reduced MHD (RMHD) model. An important distinction between these models is that RMHD solves for small-scale velocity and magnetic field fluctuations and obtains the actual viscous/resistive dissipation associated with their evolution whereas NI MHD evolves scalar moments of the fluctuating velocity and magnetic fields and approximates dissipation using an MHD turbulence phenomenology. Despite the basic differences between the models, their simulation results match remarkably well, yielding almost identical heating rates inside the corona. 
    more » « less
  2. Abstract The formation, development, and impact of slow shocks in the upstream regions of reconnecting current layers are explored. Slow shocks have been documented in the upstream regions of magnetohydrodynamic (MHD) simulations of magnetic reconnection as well as in similar simulations with thekglobalkinetic macroscale simulation model. They are therefore a candidate mechanism for preheating the plasma that is injected into the current layers that facilitate magnetic energy release in solar flares. Of particular interest is their potential role in producing the hot thermal component of electrons in flares. During multi-island reconnection, the formation and merging of flux ropes in the reconnecting current layer drives plasma flows and pressure disturbances in the upstream region. These pressure disturbances steepen into slow shocks that propagate along the reconnecting component of the magnetic field and satisfy the expected Rankine–Hugoniot jump conditions. Plasma heating arises from both compression across the shock and the parallel electric field that develops to maintain charge neutrality in a kinetic system. Shocks are weaker at lower plasmaβ, where shock steepening is slow. While these upstream slow shocks are intrinsic to the dynamics of multi-island reconnection, their contribution to electron heating remains relatively minor compared with that from Fermi reflection and the parallel electric fields that bound the reconnection outflow. 
    more » « less
  3. This paper develops a unified linear theory of cross field plasma instabilities, including the Farley–Buneman, electron thermal, and ion thermal instabilities, in spatially uniform collisional plasmas with partially unmagnetized multi-species ions. Collisional plasma instabilities in weakly ionized, highly dissipative, weakly magnetized plasmas play an important role in the lower Earth's ionosphere and may be of importance in other planetary ionospheres, stellar atmospheres, cometary tails, molecular clouds, accretion disks, etc. In the Earth's ionosphere, these collisional plasma instabilities cause intense electron heating. In the solar chromosphere, they can do the same—an effect originally suggested from spectroscopic observations and modeling. Based on a simplified 5-moment multi-fluid model, the theoretical analysis presented in this paper produces the linear dispersion relation for the combined Thermal Farley–Buneman Instability with an important long-wavelength limit analyzed in detail. This limit provides an easy interpretation of different instability drivers and wave dissipation. This analysis of instability, combined with simulations, will enable us to better understand plasma waves and turbulence in these commonly occurring collisional space plasmas. 
    more » « less
  4. Models currently fail to reproduce observations of the coldest regions in the Sun’s atmosphere, though recent work suggests the thermal Farley–Buneman instability (TFBI) may play a critical role. This meter-scale, electrostatic, multifluid plasma instability causes turbulence and heating in the coldest regions of the solar chromosphere. This paper describes how TFBI turbulence and heating varies across multifluid 2D, kinetic 2D, and kinetic 3D simulations. It also presents the first 3D simulations of the TFBI. We find that multifluid and kinetic 2D simulations produce similar results overall, despite using vastly different approaches. Additionally, our kinetic 3D simulations produce a similar or somewhat larger amount of heating compared to 2D, as contributions from the parallel electric field account for only (13 ± 2.5)% of the total turbulent heating in 3D. 
    more » « less
  5. Abstract The Magnetospheric Multiscale (MMS) mission has presented a new opportunity to study the fine scale structures and phenomena of the Earth’s magnetosphere, including cross scale processes associated with the Kelvin–Helmholtz Instability (KHI), but such studies of the KHI and its secondary processes will require a database of MMS encounters with Kelvin–Helmholtz (KH) waves. Here, we present an overview of 45 MMS observations of the KHI from September 2015 to March 2020. Growth rates and unstable solid angles for each of the 45 events were calculated using a new technique to automatically detect plasma regions on either side of the magnetopause boundary. There was no apparent correlation between solar wind conditions during the KHI and its growth rate and unstable solid angle, which is not surprising as KH waves were observed downstream of their source region. We note all KHI were observed for solar wind flow speeds between 295 and 610 km/s, possibly due to a filtering effect of the instability onset criteria and plasma compressibility. Two‐dimensional Magnetohydrodynamic (2D MHD) simulations were compared with two of the observed MMS events. Comparison of the observations with the 2D MHD simulations indicates that the new region sorting method is reliable and robust. The ability to automatically detect separate plasma regions on either side of a moving boundary and determine the KHI growth rate may prove useful for future work identifying and studying secondary processes associated with the KHI. 
    more » « less