skip to main content


Title: Data-driven Radiative Magnetohydrodynamics Simulations with the MURaM Code
Abstract

We present a method of conducting data-driven simulations of solar active regions and flux emergence with the MURaM radiative magnetohydrodynamics (MHD) code. The horizontal electric field that is derived from the full velocity and magnetic vectors is implemented at the photospheric (bottom) boundary to drive the induction equation. The energy equation accounts for thermal conduction along magnetic fields, optically thin radiative loss, and heating of coronal plasma by viscous and resistive dissipation, which allows for a realistic presentation of the thermodynamic properties of coronal plasma that are key to predicting the observational features of solar active regions and eruptions. To validate this method, the photospheric data from a comprehensive radiative MHD simulation of solar eruption (the ground truth) are used to drive a series of numerical experiments. The data-driven simulation reproduces the accumulation of free magnetic energy over the course of flux emergence in the ground truth with an error of 3%. The onset time is approximately 8 minutes delayed compared to the ground truth. However, a precursor-like signature can be identified at the correct onset time. The data-driven simulation captures key eruption-related emission features and plasma dynamics of the ground truth flare over a wide temperature span, fromlog10T=4.5tolog10T>8. The evolution of the flare and coronal mass ejection as seen in synthetic extreme ultraviolet images is also reproduced with high fidelity. This method helps to understand the evolution of magnetic field in a more realistic coronal environment and to link the magnetic structures to observable diagnostics.

 
more » « less
NSF-PAR ID:
10419781
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. 118
Size(s):
["Article No. 118"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We present a comprehensive radiative magnetohydrodynamic simulation of the quiet Sun and large solar active regions. The 197 Mm wide simulation domain spans from 18(10) Mm beneath the photosphere to 113 Mm in the solar corona. Radiative transfer assuming local thermal equilibrium, optically thin radiative losses, and anisotropic conduction transport provide the necessary realism for synthesizing observables to compare with remote-sensing observations of the photosphere and corona. This model self-consistently reproduces observed features of the quiet Sun, emerging and developed active regions, and solar flares up to M class. Here, we report an overview of the first results. The surface magneto-convection yields an upward Poynting flux that is dissipated in the corona and heats the plasma to over 1 MK. The quiescent corona also presents ubiquitous propagating waves, jets, and bright points with sizes down to 2 Mm. Magnetic flux bundles emerge into the photosphere and give rise to strong and complex active regions with over 1023Mx magnetic flux. The coronal free magnetic energy, which is approximately 18% of the total magnetic energy, accumulates to approximately 1033erg. The coronal magnetic field is clearly non-force-free, as the Lorentz force needs to balance the pressure force and viscous stress as well as drive magnetic field evolution. The emission measure fromlog10T=4.5tolog10T>7provides a comprehensive view of the active region corona, such as coronal loops of various lengths and temperatures, mass circulation by evaporation and condensation, and eruptions from jets to large-scale mass ejections.

     
    more » « less
  2. Abstract

    We measure the thermal electron energization in 1D and 2D particle-in-cell simulations of quasi-perpendicular, low-beta (βp= 0.25) collisionless ion–electron shocks with mass ratiomi/me= 200, fast Mach numberMms=1–4, and upstream magnetic field angleθBn= 55°–85° from the shock normalnˆ. It is known that shock electron heating is described by an ambipolar,B-parallel electric potential jump, Δϕ, that scales roughly linearly with the electron temperature jump. Our simulations haveΔϕ/(0.5miush2)0.1–0.2 in units of the pre-shock ions’ bulk kinetic energy, in agreement with prior measurements and simulations. Different ways to measureϕ, including the use of de Hoffmann–Teller frame fields, agree to tens-of-percent accuracy. Neglecting off-diagonal electron pressure tensor terms can lead to a systematic underestimate ofϕin our low-βpshocks. We further focus on twoθBn= 65° shocks: aMs=4(MA=1.8) case with a long, 30diprecursor of whistler waves alongnˆ, and aMs=7(MA=3.2) case with a shorter, 5diprecursor of whistlers oblique to bothnˆandB;diis the ion skin depth. Within the precursors,ϕhas a secular rise toward the shock along multiple whistler wavelengths and also has localized spikes within magnetic troughs. In a 1D simulation of theMs=4,θBn= 65° case,ϕshows a weak dependence on the electron plasma-to-cyclotron frequency ratioωpece, andϕdecreases by a factor of 2 asmi/meis raised to the true proton–electron value of 1836.

     
    more » « less
  3. Abstract

    We present a theory based on the conventional two-term (i.e. Lorentzian) approximation to the exact solution of the Boltzmann equation in non-magnetized weakly ionized plasma to efficiently obtain the electron rate and transport coefficients in a magnetized plasma for an arbitrary magnitude and direction of applied electric fieldEand magnetic fieldB. The proposed transcendental method does not require the two-term solution of the Boltzmann equation in magnetized plasma, based on which the transport parameters vary as a function of the reduced electric fieldE/N, reduced electron cyclotron frequencyωce/N, and angleE,BbetweenEandBvectors, whereNis the density of neutrals. Comparisons between the coefficients derived from BOLSIG+’s solution (obtained via the two-term expansion whenB0) and coefficients of the presented method are illustrated for air, a mixture of molecular hydrogen (H2) and helium (He) representing the giant gas planets of the Solar System, and pure carbon dioxide (CO2). The new approach may be used in the modeling of magnetized plasma encountered in the context of transient luminous events, e.g. sprite streamers in the atmosphere of Earth and Jupiter, in modeling the propagation of lightning’s electromagnetic pulses in Earth’s ionosphere, and in various laboratory and industrial applications of nonthermal plasmas.

     
    more » « less
  4. Abstract

    Parker Solar Probe (PSP) observed sub-Alfvénic solar wind intervals during encounters 8–14, and low-frequency magnetohydrodynamic (MHD) turbulence in these regions may differ from that in super-Alfvénic wind. We apply a new mode decomposition analysis to the sub-Alfvénic flow observed by PSP on 2021 April 28, identifying and characterizing entropy, magnetic islands, forward and backward Alfvén waves, including weakly/nonpropagating Alfvén vortices, forward and backward fast and slow magnetosonic (MS) modes. Density fluctuations are primarily and almost equally entropy- and backward-propagating slow MS modes. The mode decomposition provides phase information (frequency and wavenumberk) for each mode. Entropy density fluctuations have a wavenumber anisotropy ofkk, whereas slow-mode density fluctuations havek>k. Magnetic field fluctuations are primarily magnetic island modes (δBi) with anO(1) smaller contribution from unidirectionally propagating Alfvén waves (δBA+) giving a variance anisotropy ofδBi2/δBA2=4.1. Incompressible magnetic fluctuations dominate compressible contributions from fast and slow MS modes. The magnetic island spectrum is Kolmogorov-likek1.6in perpendicular wavenumber, and the unidirectional Alfvén wave spectra arek1.6andk1.5. Fast MS modes propagate at essentially the Alfvén speed with anticorrelated transverse velocity and magnetic field fluctuations and are almost exclusively magnetic due toβp≪ 1. Transverse velocity fluctuations are the dominant velocity component in fast MS modes, and longitudinal fluctuations dominate in slow modes. Mode decomposition is an effective tool in identifying the basic building blocks of MHD turbulence and provides detailed phase information about each of the modes.

     
    more » « less
  5. Abstract

    The best upper limit for the electron electric dipole moment was recently set by the ACME collaboration. This experiment measures an electron spin-precession in a cold beam of ThO molecules in their metastableH(3Δ1)state. Improvement in the statistical and systematic uncertainties is possible with more efficient use of molecules from the source and better magnetometry in the experiment, respectively. Here, we report measurements of several relevant properties of the long-livedQ(3Δ2)state of ThO, and show that this state is a very useful resource for both these purposes. TheQstate lifetime is long enough that its decay during the time of flight in the ACME beam experiment is negligible. The large electric dipole moment measured for theQstate, giving rise to a large linear Stark shift, is ideal for an electrostatic lens that increases the fraction of molecules detected downstream. The measured magnetic moment of theQstate is also large enough to be used as a sensitive co-magnetometer in ACME. Finally, we show that theQstate has a large transition dipole moment to theC(1Π1)state, which allows for efficient population transfer between the ground stateX(1Σ+)and theQstate viaXCQStimulated Raman Adiabatic Passage (STIRAP). We demonstrate 90 % STIRAP transfer efficiency. In the course of these measurements, we also determine the magnetic moment ofCstate, theXCtransition dipole moment, and branching ratios of decays from theCstate.

     
    more » « less