We present the first $N$-body simulations that adapt the equations of smoothed particle hydrodynamics to capture the effect of dark matter self-interactions which are too frequent to be resolved explicitly. The relevant energy transfer equations are derived, the appropriate thermal conductivity is determined and the effects of different smoothing kernels are studied. We apply our framework to simulate the formation of isothermal cores in isolated dark matter haloes and determine the core growth rate as a function of the self-scattering cross section. Our approach may be combined with explicit simulations of rare scatterings in order to simulate accurately the effects of arbitrary dark matter self-interactions in future cosmological simulations.

Self-interacting dark matter may have striking astrophysical signatures, such as observable offsets between galaxies and dark matter in merging galaxy clusters. Numerical N-body simulations used to predict such observables typically treat the galaxies as collisionless test particles, a questionable assumption given that each galaxy is embedded in its own dark matter halo. To enable a more accurate treatment we develop an effective description of small dark matter haloes taking into account the two major effects due to dark matter self-scatterings: deceleration and evaporation. We point out that self-scatterings can have a sizeable impact on the trajectories of galaxies, diminishing the separation between galaxies and dark matter in merging clusters. This effect depends sensitively on the underlying particle physics, in particular the angular dependence of the self-scattering cross section, and cannot be predicted from the momentum transfer cross section alone.

Self-interactions of dark matter particles can potentially lead to an observable separation between the dark matter halo and the stars of a galaxy moving through a region of large dark matter density. Such a separation has recently been observed in a galaxy falling into the core of the galaxy cluster Abell 3827. We estimated the DM self-interaction cross section needed to reproduce the observed effects and find that the sensitivity of Abell 3827 has been significantly overestimated in a previous study. Our corrected estimate is $\tilde{\sigma}/m_\text{DM} \sim 3\:\text{cm}^2\:\text{g}^{-1}$ when self-interactions result in an effective drag force and $\sigma/m_\text{DM} \sim 1.5\:\text{cm}^2\:\text{g}^{-1}$ for the case of contact interactions, in some tension with previous upper bounds.