We present a new approach to simulating mixtures of gas and dust in smoothed particle hydrodynamics (SPH). We show how the two-fluid equations can be rewritten to describe a single-fluid mixture moving with the barycentric velocity, with each particl
e carrying a dust fraction. We show how this formulation can be implemented in SPH while preserving the conservation properties (i.e. conservation of mass of each phase, momentum and energy). We also show that the method solves two key issues with the two fluid approach: it avoids over-damping of the mixture when the drag is strong and prevents a problem with dust particles becoming trapped below the resolution of the gas. We also show how the general one-fluid formulation can be simplified in the limit of strong drag (i.e. small grains) to the usual SPH equations plus a diffusion equation for the evolution of the dust fraction that can be evolved explicitly and does not require any implicit timestepping. We present tests of the simplified formulation showing that it is accurate in the small grain/strong drag limit. We discuss some of the issues we have had to solve while developing this method and finally present a preliminary application to dust settling in protoplanetary discs.
We describe a simple method for simulating the dynamics of small grains in a dusty gas, relevant to micron-sized grains in the interstellar medium and grains of centimetre size and smaller in protoplanetary discs. The method involves solving one extr
a diffusion equation for the dust fraction in addition to the usual equations of hydrodynamics. This diffusion approximation for dust is valid when the dust stopping time is smaller than the computational timestep. We present a numerical implementation using Smoothed Particle Hydrodynamics (SPH) that is conservative, accurate and fast. It does not require any implicit timestepping and can be straightforwardly ported into existing 3D codes.
We have performed Smoothed Particle Magnetohydrodynamics (SPMHD) simulations demonstrating the production of collimated jets during collapse of 1 solar mass molecular cloud cores to form the `first hydrostatic core in low mass star formation. Recentl
y a number of candidate first core objects have been observed, including L1448 IRS2E, L1451-mm and Per Bolo 58, although it is not yet clear that these are first hydrostatic cores. Recent observations of Per Bolo 58 in particular appear to show collimated, bipolar outflows which are inconsistent with previous theoretical expectations. We show that low mass first cores can indeed produce tightly collimated jets (opening angles <~ 10 degrees) with speeds of ~2-7 km/s, consistent with some of the observed candidates. We have also demonstrated, for the first time, that such phenomena can be successfully captured in SPMHD simulations.
Accounting for the Reynolds number is critical in numerical simulations of turbulence, particularly for subsonic flow. For Smoothed Particle Hydrodynamics (SPH) with constant artificial viscosity coefficient alpha, it is shown that the effective Reyn
olds number in the absence of explicit physical viscosity terms scales linearly with the Mach number - compared to mesh schemes, where the effective Reynolds number is largely independent of the flow velocity. As a result, SPH simulations with alpha=1 will have low Reynolds numbers in the subsonic regime compared to mesh codes, which may be insufficient to resolve turbulent flow. This explains the failure of Bauer and Springel (2011, arXiv:1109.4413v1) to find agreement between the moving-mesh code AREPO and the GADGET SPH code on simulations of driven, subsonic (v ~ 0.3 c_s) turbulence appropriate to the intergalactic/intracluster medium, where it was alleged that SPH is somehow fundamentally incapable of producing a Kolmogorov-like turbulent cascade. We show that turbulent flow with a Kolmogorov spectrum can be easily recovered by employing standard methods for reducing alpha away from shocks.
