We present the first multidimensional, magnetohydrodynamic simulations that capture the initial formation and longterm 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 extend our earlier multidimensional, magnetohydrodynamic simulations of coronal rain occurring in magnetic arcades with higher resolution, grid-adaptive computations covering a much longer (> 6 hour) timespan. We quantify how in-situ forming blob-like condensations grow along and across field lines and show that rain showers can occur in limit cycles, here demonstrated for the first time in 2.5D setups. We discuss dynamical, multi-dimensional aspects of the rebound shocks generated by the siphon inflows and quantify the thermodynamics of a prominence-corona-transition-region like structure surrounding the blobs. We point out the correlation between condensation rates and the cross-sectional size of loop systems where catastrophic cooling takes place. We also study the variations of the typical number density, kinetic energy and temperature while blobs descend, impact and sink into the transition region. In addition, we explain the mechanisms leading to concurrent upflows while the blobs descend. As a result, there are plenty of shear flows generated with relative velocity difference around 80 km s −1 in our simulations. These shear flows are siphon flows set up by multiple blob dynamics and they in turn affect the deformation of the falling blobs. In particular, we show how shear flows can break apart blobs into smaller fragments, within minutes.

Aims. We intend to investigate the underlying physics for the coronal rain phenomenon in a representative bipolar magnetic field, including the formation and the dynamics of coronal rain blobs. Methods. With the MPI-AMRVAC code, we performed three dimensional radiative magnetohydrodynamic (MHD) simulation with strong heating localized on footpoints of magnetic loops after a relaxation to quiet solar atmosphere. Results. Progressive cooling and in-situ condensation starts at the loop top due to radiative thermal instability. The first large-scale condensation on the loop top suffers Rayleigh-Taylor instability and becomes fragmented into smaller blobs. The blobs fall vertically dragging magnetic loops until they reach low-β regions and start to fall along the loops from loop top to loop footpoints. A statistic study of the coronal rain blobs finds that small blobs with masses of less than 10 10 g dominate the population. When blobs fall to lower regions along the magnetic loops, they are stretched and develop a non-uniform velocity pattern with an anti-parallel shearing pattern seen to develop along the central axis of the blobs. Synthetic images of simulated coronal rain with Solar Dynamics Observatory Atmospheric Imaging Assembly well resemble real observations presenting dark falling clumps in hot channels and bright rain blobs in a cool channel. We also find density inhomogeneities during a coronal rain "shower", which reflects the observed multi-stranded nature of coronal rain.

We propose a model for the formation of loop-top hard X-ray (HXR) sources in solar flares through the Inverse Compton mechanism, scattering the surrounding soft X-ray (SXR) photons to higher energy HXR photons. We simulate the consequences of a flare-driven energy deposit in the upper chromosphere in the impulsive phase of single loop flares. The consequent chromosphere evaporation flows from both footpoints reach speeds up to hundreds of kilometers per second, and we demonstrate how this triggers Kelvin-Helmholtz instability (KHI) in the loop top, under asymmetric conditions. The KHI vortices further fragment the magnetic topology into multiple magnetic islands and current sheets, and the hot plasma within leads to a bright loop-top SXR source region. We argue that the magnetohydrodynamic turbulence that appears at the loop apex could be an efficient accelerator of non-thermal particles, which the island structures can trap at the loop-top. These accelerated non-thermal particles can upscatter the surrounding thermal SXR photons emitted by the extremely hot evaporated plasma to HXR photons.

Quasi-periodic propagating intensity disturbances have been observed in large coronal loops in EUV images over a decade, and are widely accepted to be slow magnetosonic waves. However, spectroscopic observations from Hinode/EIS revealed their association with persistent coronal upflows, making this interpretation debatable. We perform a 2.5D magnetohydrodynamic simulation to imitate the chromospheric evaporation and the following reflected patterns in a flare loop. Our model encompasses the corona, transition region, and chromosphere. We demonstrate that the quasi periodic propagating intensity variations captured by the synthesized Solar Dynamics Observatory/Atmospheric Imaging Assembly (AIA) 131, 94Å emission images match the previous observations well. With particle tracers in the simulation, we confirm that these quasi periodic propagating intensity variations consist of reflected slow mode waves and mass flows with an average speed of 310 km/s in an 80 Mm length loop with an average temperature of 9 MK. With the synthesized Doppler shift velocity and intensity maps of the Solar and Heliospheric Observatory/Solar Ultraviolet Measurement of Emitted Radiation (SUMER) Fe XIX line emission, we confirm that these reflected slow mode waves are propagating waves.