Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userSelf-consistent 3D radiative MHD simulations of coronal rain formation and evolution
431
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Coronal rain consists of cool and dense plasma condensations formed in coronal loops as a result of thermal instability. Previous numerical simulations of thermal instability and coronal rain formation have relied on artificially adding a coronal heating term to the energy equation. To reproduce large-scale characteristics of the corona, using more realistic coronal heating prescription is necessary. We analyse coronal rain formation and evolution in a 3-dimensional radiative magnetohydrodynamic simulation spanning from convection zone to corona which is self-consistently heated by magnetic field braiding as a result of convective motions. We investigate the spatial and temporal evolution of energy dissipation along coronal loops which become thermally unstable. Ohmic dissipation in the model leads to the heating events capable of inducing sufficient chromospheric evaporation into the loop to trigger thermal instability and condensation formation. The cooling of the thermally unstable plasma occurs on timescales comparable to the duration of the individual impulsive heating events. The impulsive heating has sufficient duration to trigger thermal instability in the loop but does not last long enough to lead to coronal rain limit cycles. We show that condensations can either survive and fall into the chromosphere or be destroyed by strong bursts of Joule heating associated with a magnetic reconnection events. In addition, we find that condensations can also form along open magnetic field lines.
rate research
Read More
Determining the temperature distribution of coronal plasmas can provide stringent constraints on coronal heating. Current observations with the Extreme ultraviolet Imaging Spectrograph onboard Hinode and the Atmospheric Imaging Assembly onboard the Solar Dynamics Observatory provide diagnostics of the emission measure distribution (EMD) of the coronal plasma. Here we test the reliability of temperature diagnostics using 3D radiative MHD simulations. We produce synthetic observables from the models, and apply the Monte Carlo Markov chain EMD diagnostic. By comparing the derived EMDs with the true distributions from the model we assess the limitations of the diagnostics, as a function of the plasma parameters and of the signal-to-noise of the data. We find that EMDs derived from EIS synthetic data reproduce some general characteristics of the true distributions, but usually show differences from the true EMDs that are much larger than the estimated uncertainties suggest, especially when structures with significantly different density overlap along the line-of-sight. When using AIA synthetic data the derived EMDs reproduce the true EMDs much less accurately, especially for broad EMDs. The differences between the two instruments are due to the: (1) smaller number of constraints provided by AIA data, (2) broad temperature response function of the AIA channels which provide looser constraints to the temperature distribution. Our results suggest that EMDs derived from current observatories may often show significant discrepancies from the true EMDs, rendering their interpretation fraught with uncertainty. These inherent limitations to the method should be carefully considered when using these distributions to constrain coronal heating.
We study the magnetic field and current structure associated with a coronal loop. Through this we investigate to what extent the assumptions of a force-free magnetic field break down and where they might be justified. We analyse a 3D MHD model of the solar corona in an emerging active region with the focus on the structure of the forming coronal loops. The lower boundary of this simulation is taken from a model of an emerging active region. As a consequence of the emerging magnetic flux and the horizontal motions at the surface a coronal loop forms self-consistently. We investigate the current density along magnetic field inside (and outside) this loop and study the magnetic and plasma properties in and around it. We find that the total current along the loop changes its sign from being antiparallel to parallel to the magnetic field. This is caused by the inclination of the loop together with the footpoint motion. Around the loop the currents form a complex non-force-free helical structure. This is directly related to a bipolar current structure at the loop footpoints at the base of the corona and a local reduction of the background magnetic field (i.e. outside the loop) caused by the plasma flow into and along the loop. The locally reduced magnetic pressure in the loop allows the loop to sustain a higher density, which is crucial for the emission in extreme UV. The acting of the flow on the magnetic field hosting the loop turns out to be also responsible for the observed squashing of the loop. The complex magnetic field and current system surrounding it can be modeled only in 3D MHD models where the magnetic field has to balance the plasma pressure. A 1D coronal loop model or a force-free extrapolation can not capture the current system and the complex interaction of the plasma and the magnetic field in the coronal loop, despite the fact that the loop is under low-$beta$ conditions.
Acoustic waves excited in the photosphere and below might play an integral part in the heating of the solar chromosphere and corona. However, it is yet not fully clear how much of the initially acoustic wave flux reaches the corona and in what form. We investigate the wave propagation, damping, transmission, and conversion in the lower layers of the solar atmosphere using 3D numerical MHD simulations. A model of a gravitationally stratified expanding straight coronal loop, stretching from photosphere to photosphere, is perturbed at one footpoint by an acoustic driver with a period of 370 seconds. For this period acoustic cutoff regions are present below the transition region (TR). About 2% of the initial energy from the driver reach the corona. The shape of the cutoff regions and the height of the TR show a highly dynamic behavior. Taking only the driven waves into account, the waves have a propagating nature below and above the cutoff region, but are standing and evanescent within the cutoff region. Studying the driven waves together with the background motions in the model reveals standing waves between the cutoff region and the TR. These standing waves cause an oscillation of the TR height. In addition, fast or leaky sausage body-like waves might have been excited close to the base of the loop. These waves then possibly convert to fast or leaky sausage surface-like waves at the top of the main cutoff region, followed by a conversion to slow sausage body-like waves around the TR.
Flare-driven coronal rain can manifest from rapidly cooled plasma condensations near coronal loop-tops in thermally unstable post-flare arcades. We detect 5 phases that characterise the post-flare decay: heating, evaporation, conductive cooling dominance for ~120 s, radiative / enthalpy cooling dominance for ~4700 s and finally catastrophic cooling occurring within 35-124 s leading to rain strands with s periodicity of 55-70 s. We find an excellent agreement between the observations and model predictions of the dominant cooling timescales and the onset of catastrophic cooling. At the rain formation site we detect co-moving, multi-thermal rain clumps that undergo catastrophic cooling from ~1 MK to ~22000 K. During catastrophic cooling the plasma cools at a maximum rate of 22700 K s-1 in multiple loop-top sources. We calculated the density of the EUV plasma from the DEM of the multi-thermal source employing regularised inversion. Assuming a pressure balance, we estimate the density of the chromospheric component of rain to be 9.21x10^11 +-1.76x10^11 cm-3 which is comparable with quiescent coronal rain densities. With up to 8 parallel strands in the EUV loop cross section, we calculate the mass loss rate from the post-flare arcade to be as much as 1.98x10^12 +/-4.95x10^11 g s-1. Finally, we reveal a close proximity between the model predictions of 10^5.8 K and the observed properties between 10^5.9 K and 10^6.2 K, that defines the temperature onset of catastrophic cooling. The close correspondence between the observations and numerical models suggests that indeed acoustic waves (with a sound travel time of 68 s) could play an important role in redistributing energy and sustaining the enthalpy-based radiative cooling.
The reconstruction of the solar spectral irradiance (SSI) on various time scales is essential for the understanding of the Earths climate response to the SSI variability. The driver of the SSI variability is understood to be the intensity contrast of magnetic features present on the Sun with respect to the largely non-magnetic quiet Sun. However, different spectral synthesis codes lead to diverging projections of SSI variability. In this study we compare three different radiative transfer codes and carry out a detailed analysis of their performance. We perform the spectral synthesis at the continuum wavelength of 665 nm with the Code for Solar Irradiance (COSI), and the Rybicki-Hummer (RH), and Max Planck University of Chicago Radiative MHD (MURaM) codes for three 3D MHD simulations snapshots, a non-magnetic case, and MHD simulations with 100 G, and 200 G magnetic field strength. We determine the intensity distributions, the intensity differences and ratios for the spectral synthesis codes. We identify that the largest discrepancies originate in the intergranular lanes where the most field concentration occurs. Overall, the applied radiative transfer codes give consistent intensity distributions. Also, the intensity variation as a function of magnetic field strength for the particular 100 G and 200 G snapshots agree within the 2-3% range.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
