Chun Xia

MPI-AMRVAC FOR SOLAR AND ASTROPHYSICS

In this paper, we present an update to the open source MPI-AMRVAC simulation toolkit where we focus on solar and non-relativistic astrophysical magnetofluid dynamics. We highlight recent developments in terms of physics modules, such as hydrodynamics with dust coupling and the conservative implementation of Hall magnetohydrodynamics. A simple conservative high-order finite difference scheme that works in combination with all available physics modules is introduced and demonstrated with the example of monotonicity-preserving fifth-order reconstruction. Strong stability-preserving high-order Runge-Kutta time steppers are used to obtain stable evolutions in multi-dimensional applications, realizing up to fourth-order accuracy in space and time. With the new distinction between active and passive grid cells, MPI-AMRVAC is ideally suited to simulate evolutions where parts of the solution are controlled analytically or have a tendency to progress into or out of a stationary state. Typical test problems and representative applications are discussed with an outlook toward follow-up research. Finally, we discuss the parallel scaling of the code and demonstrate excellent weak scaling up to 30, 000 processors, allowing us to exploit modern peta-scale infrastructure.

THREE-DIMENSIONAL PROMINENCE-HOSTING MAGNETIC CONFIGURATIONS: CREATING A HELICAL MAGNETIC FLUX ROPE

The magnetic configuration hosting prominences and their surrounding coronal structure is a key research topic in solar physics. Recent theoretical and observational studies strongly suggest that a helical magnetic flux rope is an essential ingredient to fulfill most of the theoretical and observational requirements for hosting prominences. To understand flux rope formation details and obtain magnetic configurations suitable for future prominence formation studies, we here report on three-dimensional isothermal magnetohydrodynamic simulations including finite gas pressure and gravity. Starting from a magnetohydrostatic corona with a linear force-free bipolar magnetic field, we follow its evolution when introducing vortex flows around the main polarities and converging flows toward the polarity inversion line near the bottom of the corona. The converging flows bring the feet of different loops together at the polarity inversion line, where magnetic reconnection and flux cancellation happen. Inflow and outflow signatures of the magnetic reconnection process are identified, and thereby the newly formed helical loops wind around preexisting ones so that a complete flux rope grows and ascends. When a macroscopic flux rope is formed, we switch off the driving flows and find that the system relaxes to a stable state containing a helical magnetic flux rope embedded in an overlying arcade structure. A major part of the formed flux rope is threaded by dipped field lines that can stably support prominence matter, while the total mass of the flux rope is in the order of 4-5× 1014 g.

Simulating the in situ condensation process of solar prominences

Prominences in the solar corona are hundredfold cooler and denser than their surroundings, with a total mass of 1.e13 up to 1.e15 g. Here we report on the first comprehensive simulations of three-dimensional, thermally and gravitationally stratified magnetic flux ropes, where in situ condensation to a prominence happens due to radiative losses. After a gradual thermodynamic adjustment, we witness a phase where runaway cooling happens while counter-streaming shearing flows drain off mass along helical field lines. After this drainage, a prominence-like condensation resides in concave upward field regions, and this prominence retains its overall characteristics for more than two hours. While condensing, the prominence establishes a prominence-corona transition region, where magnetic field-aligned thermal conduction is operative during the runaway cooling. The prominence structure represents a force-balanced state in a helical flux rope. The simulated condensation demonstrates a right-bearing barb, as a remnant of the drainage. Synthetic images at extreme ultraviolet wavelengths follow the onset of the condensation, and confirm the appearance of horns and a three-part structure for the stable prominence state, as often seen in erupting prominences. This naturally explains recent Solar Dynamics Observatory views with the Atmospheric Imaging Assembly on prominences in coronal cavities demonstrating horns.

MULTIDIMENSIONAL MODELING OF CORONAL RAIN DYNAMICS

We present the first multidimensional, magnetohydrodynamic simulations that capture the initial formation and long-term sustainment of the enigmatic coronal rain phenomenon. We demonstrate how thermal instability can induce a spectacular display of in situ forming blob-like condensations which then start their intimate ballet on top of initially linear force-free arcades. Our magnetic arcades host a chromospheric, transition region, and coronal plasma. Following coronal rain dynamics for over 80 minutes of physical time, we collect enough statistics to quantify blob widths, lengths, velocity distributions, and other characteristics which directly match modern observational knowledge. Our virtual coronal rain displays the deformation of blobs into V-shaped features, interactions of blobs due to mostly pressure-mediated levitations, and gives the first views of blobs that evaporate in situ or are siphoned over the apex of the background arcade. Our simulations pave the way for systematic surveys of coronal rain showers in true multidimensional settings to connect parameterized heating prescriptions with rain statistics, ultimately allowing us to quantify the coronal heating input.We present the first multidimensional, magnetohydrodynamic simulations which capture the initial formation and the long-term sustainment of the enigmatic coronal rain phenomenon. We demonstrate how thermal instability can induce a spectacular display of in-situ forming blob-like condensations which then start their intimate ballet on top of initially linear force-free arcades. Our magnetic arcades host chromospheric, transition region, and coronal plasma. Following coronal rain dynamics for over 80 minutes physical time, we collect enough statistics to quantify blob widths, lengths, velocity distributions, and other characteristics which directly match with modern observational knowledge. Our virtual coronal rain displays the deformation of blobs into V -shaped like features, interactions of blobs due to mostly pressure-mediated levitations, and gives the first views on blobs which evaporate in situ, or get siphoned over the apex of the background arcade. Our simulations pave the way for systematic surveys of coronal rain showers in true multidimensional settings, to connect parametrized heating prescriptions with rain statistics, ultimately allowing to quantify the coronal heating input. Subject headings: magnetohydrodynamics(MHD) — Sun: corona — Sun: filaments, prominences

Parametric survey of longitudinal prominence oscillation simulations

Context. Longitudinal filament oscillations recently attracted increasing attention, while the restoring force and the damping mechanisms are still elusive.

FORMATION AND PLASMA CIRCULATION OF SOLAR PROMINENCES

Solar prominences are long-lived cool and dense plasma curtains in the hot and rarefied outer solar atmosphere or corona. The physical mechanism responsible for their formation and especially for their internal plasma circulation has been uncertain for decades. The observed ubiquitous down flows in quiescent prominences are difficult to interpret as plasma with high conductivity seems to move across horizontal magnetic field lines. Here we present three-dimensional numerical simulations of prominence formation and evolution in an elongated magnetic flux rope as a result of in-situ plasma condensations fueled by continuous plasma evaporation from the solar chromosphere. The prominence is born and maintained in a fragmented, highly dynamic state with continuous reappearance of multiple blobs and thread structures that move mainly downward dragging along mass-loaded field lines. The prominence plasma circulation is characterized by the dynamic balance between the drainage of prominence plasma back to the chromosphere and the formation of prominence plasma via continuous condensation. Plasma evaporates from the chromosphere, condenses into the prominence in the corona, and drains back to the chromosphere, establishing a stable chromosphere-corona plasma cycle. Synthetic images of the modeled prominence with the Solar Dynamics Observatory Atmospheric Imaging Assembly closely resemble actual observations, with many dynamical threads underlying an elliptical coronal cavity.

THE DYNAMICS OF FUNNEL PROMINENCES

We present numerical simulations in 2.5D settings where large-scale prominences form in situ out of coronal condensation in magnetic dips, in close agreement with early as well as recent reporting of funnel prominences. Our simulation uses full thermodynamic magnetohydrodynamics with anisotropic thermal conduction, optically thin radiative losses, and parameterized heating as main ingredients to establish a realistic arcade configuration from chromosphere to corona. The chromospheric evaporation, especially from transition region heights, ultimately causes thermal instability, and we witness the growth of a prominence suspended well above the transition region, continuously gaining mass and cross-sectional area. Several hours later, the condensation has grown into a structure connecting the prominence-corona transition region with the underlying transition region, and a continuous downward motion from the accumulated mass represents a drainage that matches observational findings. A more dynamic phase is found as well, with coronal rain, induced wave trains, and even a reconnection event when the core prominence plasma weighs down the field lines until a flux rope is formed. The upper part of the prominence is then trapped in a flux-rope structure, and we argue for its violent kink-unstable eruption as soon as the (ignored) length dimension would allow for ideal kink deformations.

SOLAR PROMINENCES: "DOUBLE, DOUBLE ... BOIL AND BUBBLE"

Observations revealed rich dynamics within prominences, the cool 10,000 K, macroscopic (sizes of order 100 Mm) clouds in the million degree solar corona. Even quiescent prominences are continuously perturbed by hot, rising bubbles. Since prominence matter is hundredfold denser than coronal plasma, this bubbling is related to Rayleigh-Taylor instabilities. Here we report on true macroscopic simulations well into this bubbling phase, adopting a magnetohydrodynamic description from chromospheric layers up to 30 Mm height. Our virtual prominences rapidly establish fully non-linear (magneto)convective motions where hot bubbles interplay with falling pillars, with dynamical details including upwelling pillars forming within bubbles. Our simulations show impacting Rayleigh-Taylor fingers reflecting on transition region plasma, ensuring that cool, dense chromospheric material gets mixed with prominence matter up to very large heights. This offers an explanation for the return mass cycle mystery for prominence material. Synthetic views at extreme ultraviolet wavelengths show remarkable agreement with observations, with clear indications of shear-flow induced fragmentations.

Simulating coronal condensation dynamics in 3D

We present numerical simulations in 3D settings where coronal rain phenomena take place in a magnetic configuration of a quadrupolar arcade system. Our simulation is a magnetohydrodynamic simulation including anisotropic thermal conduction, optically thin radiative losses, and parametrised heating as main thermodynamical features to construct a realistic arcade configuration from chromospheric to coronal heights. The plasma evaporation from chromospheric and transition region heights eventually causes localised runaway condensation events and we witness the formation of plasma blobs due to thermal instability, that evolve dynamically in the heated arcade part and move gradually downwards due to interchange type dynamics. Unlike earlier 2.5D simulations, in this case there is no large scale prominence formation observed, but a continuous coronal rain develops which shows clear indications of Rayleigh–Taylor or interchange instability, that causes the denser plasma located above the transition region to fall down, as the system moves towards a more stable state. Linear stability analysis is used in the non-linear regime for gaining insight and giving a prediction of the system’s evolution. After the plasma blobs descend through interchange, they follow the magnetic field topology more closely in the lower coronal regions, where they are guided by the magnetic dips.

MAGNETO-FRICTIONAL MODELING OF CORONAL NONLINEAR FORCE-FREE FIELDS. I. TESTING WITH ANALYTIC SOLUTIONS

We report our implementation of the magneto-frictional method in the Message Passing Interface Adaptive Mesh Refinement Versatile Advection Code (MPI-AMRVAC). The method aims at applications where local adaptive mesh refinement (AMR) is essential to make follow-up dynamical modeling affordable. We quantify its performance in both domain-decomposed uniform grids and block-adaptive AMR computations, using all frequently employed force-free, divergence-free, and other vector comparison metrics. As test cases, we revisit the semi-analytic solution of Low and Lou in both Cartesian and spherical geometries, along with the topologically challenging Titov-Demoulin model. We compare different combinations of spatial and temporal discretizations, and find that the fourth-order central difference with a local Lax-Friedrichs dissipation term in a single-step marching scheme is an optimal combination. The initial condition is provided by the potential field, which is the potential field source surface model in spherical geometry. Various boundary conditions are adopted, ranging from fully prescribed cases where all boundaries are assigned with the semi-analytic models, to solar-like cases where only the magnetic field at the bottom is known. Our results demonstrate that all the metrics compare favorably to previous works in both Cartesian and spherical coordinates. Cases with several AMR levels perform in accordance with their effective resolutions. The magneto-frictional method in MPI-AMRVAC allows us to model a region of interest with high spatial resolution and large field of view simultaneously, as required by observation-constrained extrapolations using vector data provided with modern instruments. The applications of the magneto-frictional method to observations are shown in an accompanying paper.