A solution to the overdamping problem when simulating dust–gas mixtures with smoothed particle hydrodynamics

Daniel J. Price; Guillaume Laibe

doi:10.1093/MNRAS/STAA1366

A solution to the overdamping problem when simulating dust–gas mixtures with smoothed particle hydrodynamics

Daniel J. Price, Guillaume Laibe

School of Physics and Astronomy

Research output: Contribution to journal › Article › Research › peer-review

14 Citations (Scopus)

Abstract

We present a fix to the overdamping problem found by Laibe & Price when simulating strongly coupled dust–gas mixtures using two different sets of particles using smoothed particle hydrodynamics. Our solution is to compute the drag at the barycentre between gas and dust particle pairs when computing the drag force by reconstructing the velocity field, similar to the procedure in Godunov-type solvers. This fixes the overdamping problem at negligible computational cost, but with additional memory required to store velocity derivatives. We employ slope limiters to avoid spurious oscillations at shocks, finding the van Leer Monotonized Central limiter most effective.

Original language	English
Pages (from-to)	3929-3934
Number of pages	6
Journal	Monthly Notices of the Royal Astronomical Society
Volume	495
Issue number	4
DOIs	https://doi.org/10.1093/MNRAS/STAA1366
Publication status	Published - Jul 2020

Keywords

Dust, extinction
Hydrodynamics
Methods: numerical
Protoplanetary discs

Access to Document

10.1093/MNRAS/STAA1366

https://arxiv.org/abs/2005.06562Licence: Other

Cite this

@article{87f7332bdf7c4897b4dcee720519a382,

title = "A solution to the overdamping problem when simulating dust–gas mixtures with smoothed particle hydrodynamics",

abstract = "We present a fix to the overdamping problem found by Laibe & Price when simulating strongly coupled dust–gas mixtures using two different sets of particles using smoothed particle hydrodynamics. Our solution is to compute the drag at the barycentre between gas and dust particle pairs when computing the drag force by reconstructing the velocity field, similar to the procedure in Godunov-type solvers. This fixes the overdamping problem at negligible computational cost, but with additional memory required to store velocity derivatives. We employ slope limiters to avoid spurious oscillations at shocks, finding the van Leer Monotonized Central limiter most effective.",

keywords = "Dust, extinction, Hydrodynamics, Methods: numerical, Protoplanetary discs",

author = "Price, {Daniel J.} and Guillaume Laibe",

note = "Funding Information: We thank Pablo Loren-Aguilar, Matthew Bate, Christophe Pinte, Ugo Lebreuilly, Benoit Commerc¸on, Jim Stone, and James Wadsley for useful discussions, and the referee for helpful comments. DJP thanks Bernhard Mueller for useful lecture notes, his teaching load for inspiration, and is grateful for funding from the Australian Research Council via FT130100034 and DP180104235. We acknowledge computing time on Gadi via the Australian National Compute facility, and on Ozstar, funded by the Australian government and Swinburne University. GL acknowledges funding from PNP, PNPS, PCMI of CNRS/INSU, CEA and CNES, France, and via the IDEXLyon project (contract ANR-16-IDEX-0005) under University of Lyon. We acknowledge European Research Council (ERC) funding under H2020 grant 864965. We used SPLASH (Price 2007). Publisher Copyright: {\textcopyright} 2020 The Author(s)",

year = "2020",

month = jul,

doi = "10.1093/MNRAS/STAA1366",

language = "English",

volume = "495",

pages = "3929--3934",

journal = "Monthly Notices of the Royal Astronomical Society",

issn = "0035-8711",

publisher = "Oxford University Press",

number = "4",

}

TY - JOUR

T1 - A solution to the overdamping problem when simulating dust–gas mixtures with smoothed particle hydrodynamics

AU - Price, Daniel J.

AU - Laibe, Guillaume

N1 - Funding Information: We thank Pablo Loren-Aguilar, Matthew Bate, Christophe Pinte, Ugo Lebreuilly, Benoit Commerc¸on, Jim Stone, and James Wadsley for useful discussions, and the referee for helpful comments. DJP thanks Bernhard Mueller for useful lecture notes, his teaching load for inspiration, and is grateful for funding from the Australian Research Council via FT130100034 and DP180104235. We acknowledge computing time on Gadi via the Australian National Compute facility, and on Ozstar, funded by the Australian government and Swinburne University. GL acknowledges funding from PNP, PNPS, PCMI of CNRS/INSU, CEA and CNES, France, and via the IDEXLyon project (contract ANR-16-IDEX-0005) under University of Lyon. We acknowledge European Research Council (ERC) funding under H2020 grant 864965. We used SPLASH (Price 2007). Publisher Copyright: © 2020 The Author(s)

PY - 2020/7

Y1 - 2020/7

N2 - We present a fix to the overdamping problem found by Laibe & Price when simulating strongly coupled dust–gas mixtures using two different sets of particles using smoothed particle hydrodynamics. Our solution is to compute the drag at the barycentre between gas and dust particle pairs when computing the drag force by reconstructing the velocity field, similar to the procedure in Godunov-type solvers. This fixes the overdamping problem at negligible computational cost, but with additional memory required to store velocity derivatives. We employ slope limiters to avoid spurious oscillations at shocks, finding the van Leer Monotonized Central limiter most effective.

AB - We present a fix to the overdamping problem found by Laibe & Price when simulating strongly coupled dust–gas mixtures using two different sets of particles using smoothed particle hydrodynamics. Our solution is to compute the drag at the barycentre between gas and dust particle pairs when computing the drag force by reconstructing the velocity field, similar to the procedure in Godunov-type solvers. This fixes the overdamping problem at negligible computational cost, but with additional memory required to store velocity derivatives. We employ slope limiters to avoid spurious oscillations at shocks, finding the van Leer Monotonized Central limiter most effective.

KW - Dust, extinction

KW - Hydrodynamics

KW - Methods: numerical

KW - Protoplanetary discs

UR - http://www.scopus.com/inward/record.url?scp=85097041387&partnerID=8YFLogxK

U2 - 10.1093/MNRAS/STAA1366

DO - 10.1093/MNRAS/STAA1366

M3 - Article

AN - SCOPUS:85097041387

SN - 0035-8711

VL - 495

SP - 3929

EP - 3934

JO - Monthly Notices of the Royal Astronomical Society

JF - Monthly Notices of the Royal Astronomical Society

IS - 4

ER -

A solution to the overdamping problem when simulating dust–gas mixtures with smoothed particle hydrodynamics

Abstract

Keywords

Access to Document

Other files and links

DUSTBUSTERS: Dust and gas in planet forming discs

Gaps, rings and holes in protoplanetary discs - signs of newborn planets?

The birth of stars and planets

Cite this

A solution to the overdamping problem when simulating dust–gas mixtures with smoothed particle hydrodynamics

Abstract

Keywords

Access to Document

Other files and links

Projects

DUSTBUSTERS: Dust and gas in planet forming discs

Gaps, rings and holes in protoplanetary discs - signs of newborn planets?

The birth of stars and planets

Cite this