Testing the effect of resolution on gravitational fragmentation with Lagrangian hydrodynamic schemes

Abstract

ABSTRACT To study the resolution required for simulating gravitational fragmentation with newly developed Lagrangian hydrodynamic schemes, meshless finite-volume method (MFV) and meshless finite-mass method, we have performed a number of simulations of the Jeans test and compared the results with both the expected analytical solution and results from the more standard Lagrangian approach: smoothed particle hydrodynamics (SPH). We find that the different schemes converge to the analytical solution when the diameter of a fluid element is smaller than a quarter of the Jeans wavelength, λJ. Among the three schemes, SPH/MFV shows the fastest/slowest convergence to the analytical solution. Unlike the well-known behaviour of Eulerian schemes, none of the Lagrangian schemes investigated displays artificial fragmentation when the perturbation wavelength, λ, is shorter than λJ, even at low numerical resolution. For larger wavelengths (λ > λJ), the growth of the perturbation is delayed when it is not well resolved. Furthermore, with poor resolution, the fragmentation seen with the MFV scheme proceeds very differently compared to the converged solution. All these results suggest that, when unresolved, the ratio of the magnitude of hydrodynamic force to that of self-gravity at the sub-resolution scale is the largest/smallest in MFV/SPH, the reasons for which we have discussed in detail. These tests are repeated to investigate the effect of kernels of higher order than the fiducial cubic spline. Our results indicate that the standard deviation of the kernel is a more appropriate definition of the ‘size’ of a fluid element than its compact support radius.