Effect of dust on Kelvin-Helmholtz instabilities

Dust is present in a large variety of astrophysical fluids, from tori around supermassive black holes to molecular clouds, protoplanetary discs, and cometary outflows. In many such fluids, shearing flows are present, leading to the formation of Kelvin-Helmholtz instabilities (KHI) and changing the properties and structures of the fluid through processes such as mixing and clumping of dust. We investigate how dust changes the growth rates of the KHI in 2D and 3D and how the it redistributes and clumps dust. We investigate if similarities can be found between the structures in 3D KHI and those seen in observations of molecular clouds. We do this by performing numerical hydrodynamical dust+gas simulations with in addition to gas a number of dust fluids. Each dust fluid represents a portion of the particle size-distribution. We study how dust-to-gas mass density ratios between 0.01 and 1 alter the growth rate in the linear phase of the KHI. We do this for a wide range of perturbation wavelengths, and compare these values to the analytical gas-only growth rates. As the formation of high-density dust structures is of interest in many astrophysical environments, we scale our simulations with physical quantities similar to values in molecular clouds. Large differences in dynamics are seen for different grain sizes. We demonstrate that high dust-to-gas ratios significantly reduce the growth rate of the KHI, especially for short wavelengths. We compare the dynamics in 2D and 3D simulations, where the latter demonstrates additional full 3D instabilities during the non-linear phase, leading to increased dust densities. We compare the structures formed by the KHI in 3D simulations with those in molecular clouds and see how the column density distribution of the simulation shares similarities with log-normal distributions with power-law tails sometimes seen in observations of molecular clouds.Comment: 14 pages, 20 figure

Radiative cooling in numerical astrophysics: the need for adaptive mesh refinement

Energy loss through optically thin radiative cooling plays an important part in the evolution of astrophysical gas dynamics and should therefore be considered a necessary element in any numerical simulation. Although the addition of this physical process to the equations of hydrodynamics is straightforward, it does create numerical challenges that have to be overcome in order to ensure the physical correctness of the simulation. First, the cooling has to be treated (semi-)implicitly, owing to the discrepancies between the cooling timescale and the typical timesteps of the simulation. Secondly, because of its dependence on a tabulated cooling curve, the introduction of radiative cooling creates the necessity for an interpolation scheme. In particular, we will argue that the addition of radiative cooling to a numerical simulation creates the need for extremely high resolution, which can only be fully met through the use of adaptive mesh refinement.Comment: 11 figures. Accepted for publication in Computers & Fluid

Solar flares and Kelvin-Helmholtz instabilities: A parameter survey

Hard X-ray (HXR) sources are frequently observed near the top of solar flare loops, and the emission is widely ascribed to bremsstrahlung. We here revisit an alternative scenario which stresses the importance of inverse Compton processes and the Kelvin- Helmholtz instability (KHI) proposed by Fang et al. (2016). This scenario adds a novel ingredient to the standard flare model, where evaporation flows from flare-impacted chromospheric foot-points interact with each other near the loop top and produce turbulence via KHI. The turbulence can act as a trapping region and as an efficient accelerator to provide energetic electrons, which scatter soft X-ray (SXR) photons to HXR photons via the inverse Compton mechanism. This paper focuses on the trigger of the KHI and the resulting turbulence in this new scenario. We perform a parameter survey to investigate the necessary ingredients to obtain KHI through interaction of chromospheric evaporation flows. When turbulence is produced in the loop apex, an index of -5/3 can be found in the spectra of velocity and magnetic field fluctuations. The KHI development and the generation of turbulence are controlled by the amount of energy deposited in the chromospheric foot-points and the time scale of its energy deposition, but typical values for M class flares show the KHI development routinely. Asymmetry of energy deposition determines the location where the turbulence is produced, and the synthesized SXR light curve shows a clear periodic signal related to the sloshing motion of the vortex pattern created by the KHI.Comment: 12 pages, 14 figure

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.Comment: 18 pages, 5 figure

Coronal rain in magnetic bipolar weak fields

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. 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. 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 beta 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.Comment: 8 figure