Dusty gas with smoothed particle hydrodynamics - I. Algorithm and test suite

Guillaume Laibe, Daniel James Price

Research output: Contribution to journalArticleResearchpeer-review

84 Citations (Scopus)


We present a new algorithm for simulating two-fluid gas and dust mixtures in smoothed particle hydrodynamics (SPH), systematically addressing a number of key issues including the generalized SPH density estimate in multifluid systems, the consistent treatment of variable smoothing-length terms, finite particle size, timestep stability, thermal coupling terms and the choice of kernel and smoothing length used in the drag operator. We find that using double-hump-shaped kernels improves the accuracy of the drag interpolation by a factor of several hundred compared to the use of standard SPH bell-shaped kernels, at no additional computational expense. In order to benchmark our algorithm, we have developed a comprehensive suite of standardized, simple test problems for gas and dust mixtures: DUSTYBOX, DUSTYWAVE, DUSTYSHOCK, DUSTYSEDOV and DUSTYDISC, the first three of which have known analytic solutions.
Original languageEnglish
Pages (from-to)2345 - 2364
Number of pages20
JournalMonthly Notices of the Royal Astronomical Society
Issue number3
Publication statusPublished - 2012

Cite this