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.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Global MHD Simulations of the Time-dependent Corona
Abstract We describe, test, and apply a technique to incorporate full-Sun, surface flux evolution into an MHD model of the global solar corona. Requiring only maps of the evolving surface flux, our method is similar to that of Lionello et al., but we introduce two ways to correct the electric field at the lower boundary to mitigate spurious currents. We verify the accuracy of our procedures by comparing to a reference simulation, driven with known flows and electric fields. We then present a thermodynamic MHD calculation lasting one solar rotation driven by maps from the magnetic flux evolution model of Schrijver & DeRosa. The dynamic, time-dependent nature of the model corona is illustrated by examining the evolution of the open flux boundaries and forward-modeled EUV emission, which evolve in response to surface flows and the emergence and cancellation flux. Although our main goal is to present the method, we briefly investigate the relevance of this evolution to properties of the slow solar wind, examining the mapping of dipped field lines to the topological signatures of the “S-Web” and comparing charge state ratios computed in the time-dependently driven run to a steady-state equivalent. Interestingly, we find that driving on its own does not significantly improve the charge state ratios, at least in this modest resolution run that injects minimal helicity. Still, many aspects of the time-dependently driven model cannot be captured with traditional steady-state methods, and such a technique may be particularly relevant for the next generation of solar wind and coronal mass ejection models.  more » « less
Award ID(s):
1854790
PAR ID:
10523486
Author(s) / Creator(s):
; ; ; ; ; ; ;
Publisher / Repository:
IOP
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
959
Issue:
2
ISSN:
0004-637X
Page Range / eLocation ID:
77
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; ; (, The Astrophysical Journal Letters)
    Abstract We present in this Letter the first global comparison between traditional line-tied steady-state magnetohydrodynamic models and a new, fully time-dependent thermodynamic magnetohydrodynamic simulation of the global corona. To approximate surface magnetic field distributions and magnitudes around solar minimum, we use the Lockheed Evolving Surface-Flux Assimilation Model to obtain input maps that incorporate flux emergence and surface flows over a full solar rotation, including differential rotation and meridional flows. Each time step evolves the previous state of the plasma with a new magnetic field input boundary condition, mimicking photospheric driving on the Sun. We find that this method produces a qualitatively different corona compared to steady-state models. The magnetic energy levels are higher in the time-dependent model, and coronal holes evolve more along the following edge than they do in steady-state models. Coronal changes, as illustrated with forward-modeled emission maps, evolve on longer timescales with time-dependent driving. We discuss implications for active and quiet Sun scenarios, solar wind formation, and widely used steady-state assumptions like potential field source surface calculations. 
    more » « less
  2. ; ; ; ; ; ; ; ; ; ; et al (, The Astrophysical Journal)
    Abstract This paper addresses the first direct investigation of the energy budget in the solar corona. Exploiting joint observations of the same coronal plasma by Parker Solar Probe and the Metis coronagraph aboard Solar Orbiter and the conserved equations for mass, magnetic flux, and wave action, we estimate the values of all terms comprising the total energy flux of the proton component of the slow solar wind from 6.3 to 13.3 R ⊙ . For distances from the Sun to less than 7 R ⊙ , we find that the primary source of solar wind energy is magnetic fluctuations including Alfvén waves. As the plasma flows away from the low corona, magnetic energy is gradually converted into kinetic energy, which dominates the total energy flux at heights above 7 R ⊙ . It is found too that the electric potential energy flux plays an important role in accelerating the solar wind only at altitudes below 6 R ⊙ , while enthalpy and heat fluxes only become important at even lower heights. The results finally show that energy equipartition does not exist in the solar corona. 
    more » « less
  3. ; ; (, The astrophysical journal)
    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
  4. ; ; ; ; (, Journal of Geophysical Research: Space Physics)
    Abstract During magnetospheric substorms, plasma from magnetic reconnection in the magnetotail is thought to reach the inner magnetosphere and form a partial ring current. We simulate this process using a fully kinetic 3D particle‐in‐cell (PIC) numerical code along with a global magnetohydrodynamics (MHD) model. The PIC simulation extends from the solar wind outside the bow shock to beyond the reconnection region in the tail, while the MHD code extends much further and is run for nominal solar wind parameters and a southward interplanetary magnetic field. By the end of the PIC calculation, ions and electrons from the tail reconnection reach the inner magnetosphere and form a partial ring current and diamagnetic current. The primary source of particles to the inner magnetosphere is bursty bulk flows (BBFs) that originate from a complex pattern of reconnection in the near‐Earth magnetotail at to . Most ion acceleration occurs in this region, gaining from 10 to 50 keV as they traverse the sites of active reconnection. Electrons jet away from the reconnection region much faster than the ions, setting up an ambipolar electric field allowing the ions to catch up after approximately 10 ion inertial lengths. The initial energy flux in the BBFs is mainly kinetic energy flux from the ions, but as they move earthward, the energy flux changes to enthalpy flux at the ring current. The power delivered from the tail reconnection in the simulation to the inner magnetosphere is  W, which is consistent with observations. 
    more » « less
  5. ; ; (, The Astrophysical Journal)
    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
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email