No Arabic abstract
Transport coefficients in highly ionised plasmas like the intra-cluster medium (ICM) are still ill-constrained. They influence various processes, among them the mixing at shear flow interfaces due to the Kelvin-Helmholtz instability (KHI). The observed structure of potential mixing layers can be used to infer the transport coefficients, but the data interpretation requires a detailed knowledge of the long-term evolution of the KHI under different conditions. Here we present the first systematic numerical study of the effect of constant and temperature-dependent isotropic viscosity over the full range of possible values. We show that moderate viscosities slow down the growth of the KHI and reduce the height of the KHI rolls and their rolling-up. Viscosities above a critical value suppress the KHI. The effect can be quantified in terms of the Reynolds number Re = U{lambda}/{ u}, where U is the shear velocity, {lambda} the perturbation length, and { u} the kinematic viscosity. We derive the critical Re for constant and temperature dependent, Spitzer-like viscosities, an empirical relation for the viscous KHI growth time as a function of Re and density contrast, and describe special behaviours for Spitzer-like viscosities and high density contrasts. Finally, we briefly discuss several astrophysical situations where the viscous KHI could play a role, i.e., sloshing cold fronts, gas stripping from galaxies, buoyant cavities, ICM turbulence, and high velocity clouds.
In this paper we investigate whether Smoothed Particle Hydrodynamics (SPH), equipped with artificial conductivity, is able to capture the physics of density/energy discontinuities in the case of the so-called shearing layers test, a test for examining Kelvin-Helmholtz (KH) instabilities. We can trace back each failure of SPH to show KH rolls to two causes: i) shock waves travelling in the simulation box and ii) particle clumping, or more generally, particle noise. The probable cause of shock waves is the Local Mixing Instability (LMI), previously identified in the literature. Particle noise on the other hand is a problem because it introduces a large error in the SPH momentum equation. We also investigate the role of artificial conductivity (AC). Including AC is necessary for the long-term behavior of the simulation (e.g. to get $lambda=1/2, 1$ KH rolls). In sensitive hydrodynamical simulations great care is however needed in selecting the AC signal velocity, with the default formulation leading to too much energy diffusion. We present new signal velocities that lead to less diffusion. The effects of the shock waves and of particle disorder become less important as the time-scale of the physical problem (for the shearing layers problem: lower density contrast and higher Mach numbers) decreases. At the resolution of current galaxy formation simulations mixing is probably not important. However, mixing could become crucial for next-generation simulations.
Sloshing cold fronts (CFs) arise from minor merger triggered gas sloshing. Their detailed structure depends on the properties of the intra-cluster medium (ICM): hydrodynamical simulations predict the CFs to be distorted by Kelvin-Helmholtz instabilities (KHIs), but aligned magnetic fields, viscosity, or thermal conduction can suppress the KHIs. Thus, observing the detailed structure of sloshing CFs can be used to constrain these ICM properties. Both smooth and distorted sloshing CFs have been observed, indicating that the KHI is suppressed in some clusters, but not in all. Consequently, we need to address at least some sloshing clusters individually before drawing general conclusions about the ICM properties. We present the first detailed attempt to constrain the ICM properties in a specific cluster from the structure of its sloshing CF. Proximity and brightness make the Virgo cluster an ideal target. We combine observations and Virgo-specific hydrodynamical sloshing simulations. Here we focus on a Spitzer-like temperature dependent viscosity as a mechanism to suppress the KHI, but discuss the alternative mechanisms in detail. We identify the CF at 90 kpc north and north-east of the Virgo center as the best location in the cluster to observe a possible KHI suppression. For viscosities $gtrsim$ 10% of the Spitzer value KHIs at this CF are suppressed. We describe in detail the observable signatures at low and high viscosities, i.e. in the presence or absence of KHIs. We find indications for a low ICM viscosity in archival XMM-Newton data and demonstrate the detectability of the predicted features in deep Chandra observations.
A 2D linear theory of the instability of Sweet-Parker (SP) current sheets is developed in the framework of Reduced MHD. A local analysis is performed taking into account the dependence of a generic equilibrium profile on the outflow coordinate. The plasmoid instability [Loureiro et al, Phys. Plasmas {bf 14}, 100703 (2007)] is recovered, i.e., current sheets are unstable to the formation of a large-wave-number chain of plasmoids ($k_{rm max}Lsheet sim S^{3/8}$, where $k_{rm max}$ is the wave-number of fastest growing mode, $S=Lsheet V_A/eta$ is the Lundquist number, $Lsheet$ is the length of the sheet, $V_A$ is the Alfven speed and $eta$ is the plasma resistivity), which grows super-Alfvenically fast ($gmaxtau_Asim S^{1/4}$, where $gmax$ is the maximum growth rate, and $tau_A=Lsheet/V_A$). For typical background profiles, the growth rate and the wave-number are found to {it increase} in the outflow direction. This is due to the presence of another mode, the Kelvin-Helmholtz (KH) instability, which is triggered at the periphery of the layer, where the outflow velocity exceeds the Alfven speed associated with the upstream magnetic field. The KH instability grows even faster than the plasmoid instability, $gmax tau_A sim k_{rm max} Lsheetsim S^{1/2}$. The effect of viscosity ($ u$) on the plasmoid instability is also addressed. In the limit of large magnetic Prandtl numbers, $Pm= u/eta$, it is found that $gmaxsim S^{1/4}Pm^{-5/8}$ and $k_{rm max} Lsheetsim S^{3/8}Pm^{-3/16}$, leading to the prediction that the critical Lundquist number for plasmoid instability in the $Pmgg1$ regime is $Scritsim 10^4Pm^{1/2}$. These results are verified via direct numerical simulation of the linearized equations, using a new, analytical 2D SP equilibrium solution.
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.
We present results from two sim30 ks Chandra observations of the hot atmospheres of the merging galaxy groups centered around NGC 7618 and UGC 12491. Our images show the presence of arc-like sloshing cold fronts wrapped around each group center and sim100 kpc long spiral tails in both groups. Most interestingly, the cold fronts are highly distorted in both groups, exhibiting wings along the fronts. These features resemble the structures predicted from non-viscous hydrodynamic simulations of gas sloshing, where Kelvin-Helmholtz instabilities (KHIs) distort the cold fronts. This is in contrast to the structure seen in many other sloshing and merger cold fronts, which are smooth and featureless at the current observational resolution. Both magnetic fields and viscosity have been invoked to explain the absence of KHIs in these smooth cold fronts, but the NGC 7618/UGC 12491 pair are two in a growing number of both sloshing and merger cold fronts that appear distorted. Magnetic fields and/or viscosity may be able to suppress the growth of KHIs at the cold fronts in some clusters and groups, but clearly not in all. We propose that the presence or absence of KHI-distortions in cold fronts can be used as a measure of the effective viscosity and/or magnetic field strengths in the ICM.