Planet-disc interaction on a freely moving mesh

D. J. Muñoz, Kaitlin Kratter, V. Springel, L. Hernquist

Research output: Contribution to journalArticle

12 Citations (Scopus)

Abstract

General-purpose, moving-mesh schemes for hydrodynamics have opened the possibility of combining the accuracy of grid-based numerical methods with the flexibility and automatic resolution adaptivity of particle-based methods. Due to their supersonic nature, Keplerian accretion discs are in principle a very attractive system for applying such freely moving-mesh techniques. However, the high degree of symmetry of simple accretion disc models can be difficult to capture accurately by these methods, due to the generation of geometric grid noise and associated numerical diffusion, which is absent in polar grids. To explore these and other issues, in this work we study the idealized problem of two-dimensional planet-disc interaction with the moving-mesh code AREPO. We explore the hydrodynamic evolution of discs with planets through a series of numerical experiments that vary the planet mass, the disc viscosity and the mesh resolution, and compare the resulting surface density, vortensity field and tidal torque with results from the literature.We find that the performance of the moving-mesh code in this problem is in accordance with published results, showing good consistency with grid codes written in polar coordinates. We also conclude that grid noise and mesh distortions do not introduce excessive numerical diffusion. Finally, we show how the moving-mesh approach can help in resolving an outstanding challenge for polar-coordinate grid codes, namely the successful implementation of adaptive mesh refinement in regions of high density around planets and planetary wakes, while retaining the background flow at low resolution.

Original languageEnglish (US)
Pages (from-to)3475-3495
Number of pages21
JournalMonthly Notices of the Royal Astronomical Society
Volume445
Issue number4
DOIs
StatePublished - Sep 29 2014

Fingerprint

planets
mesh
planet
grids
hydrodynamics
accretion
interactions
polar coordinates
accretion disks
torque
numerical method
symmetry
viscosity
retaining
code
wakes
flexibility
experiment
method

Keywords

  • Hydrodynamics
  • Instabilities
  • Methods: numerical
  • Planet-disc interactions
  • Protoplanetary discs

ASJC Scopus subject areas

  • Space and Planetary Science
  • Astronomy and Astrophysics

Cite this

Planet-disc interaction on a freely moving mesh. / Muñoz, D. J.; Kratter, Kaitlin; Springel, V.; Hernquist, L.

In: Monthly Notices of the Royal Astronomical Society, Vol. 445, No. 4, 29.09.2014, p. 3475-3495.

Research output: Contribution to journalArticle

Muñoz, D. J. ; Kratter, Kaitlin ; Springel, V. ; Hernquist, L. / Planet-disc interaction on a freely moving mesh. In: Monthly Notices of the Royal Astronomical Society. 2014 ; Vol. 445, No. 4. pp. 3475-3495.
@article{e852938614a345fcb029b65f7d2d1a1c,
title = "Planet-disc interaction on a freely moving mesh",
abstract = "General-purpose, moving-mesh schemes for hydrodynamics have opened the possibility of combining the accuracy of grid-based numerical methods with the flexibility and automatic resolution adaptivity of particle-based methods. Due to their supersonic nature, Keplerian accretion discs are in principle a very attractive system for applying such freely moving-mesh techniques. However, the high degree of symmetry of simple accretion disc models can be difficult to capture accurately by these methods, due to the generation of geometric grid noise and associated numerical diffusion, which is absent in polar grids. To explore these and other issues, in this work we study the idealized problem of two-dimensional planet-disc interaction with the moving-mesh code AREPO. We explore the hydrodynamic evolution of discs with planets through a series of numerical experiments that vary the planet mass, the disc viscosity and the mesh resolution, and compare the resulting surface density, vortensity field and tidal torque with results from the literature.We find that the performance of the moving-mesh code in this problem is in accordance with published results, showing good consistency with grid codes written in polar coordinates. We also conclude that grid noise and mesh distortions do not introduce excessive numerical diffusion. Finally, we show how the moving-mesh approach can help in resolving an outstanding challenge for polar-coordinate grid codes, namely the successful implementation of adaptive mesh refinement in regions of high density around planets and planetary wakes, while retaining the background flow at low resolution.",
keywords = "Hydrodynamics, Instabilities, Methods: numerical, Planet-disc interactions, Protoplanetary discs",
author = "Mu{\~n}oz, {D. J.} and Kaitlin Kratter and V. Springel and L. Hernquist",
year = "2014",
month = "9",
day = "29",
doi = "10.1093/mnras/stu1918",
language = "English (US)",
volume = "445",
pages = "3475--3495",
journal = "Monthly Notices of the Royal Astronomical Society",
issn = "0035-8711",
publisher = "Oxford University Press",
number = "4",

}

TY - JOUR

T1 - Planet-disc interaction on a freely moving mesh

AU - Muñoz, D. J.

AU - Kratter, Kaitlin

AU - Springel, V.

AU - Hernquist, L.

PY - 2014/9/29

Y1 - 2014/9/29

N2 - General-purpose, moving-mesh schemes for hydrodynamics have opened the possibility of combining the accuracy of grid-based numerical methods with the flexibility and automatic resolution adaptivity of particle-based methods. Due to their supersonic nature, Keplerian accretion discs are in principle a very attractive system for applying such freely moving-mesh techniques. However, the high degree of symmetry of simple accretion disc models can be difficult to capture accurately by these methods, due to the generation of geometric grid noise and associated numerical diffusion, which is absent in polar grids. To explore these and other issues, in this work we study the idealized problem of two-dimensional planet-disc interaction with the moving-mesh code AREPO. We explore the hydrodynamic evolution of discs with planets through a series of numerical experiments that vary the planet mass, the disc viscosity and the mesh resolution, and compare the resulting surface density, vortensity field and tidal torque with results from the literature.We find that the performance of the moving-mesh code in this problem is in accordance with published results, showing good consistency with grid codes written in polar coordinates. We also conclude that grid noise and mesh distortions do not introduce excessive numerical diffusion. Finally, we show how the moving-mesh approach can help in resolving an outstanding challenge for polar-coordinate grid codes, namely the successful implementation of adaptive mesh refinement in regions of high density around planets and planetary wakes, while retaining the background flow at low resolution.

AB - General-purpose, moving-mesh schemes for hydrodynamics have opened the possibility of combining the accuracy of grid-based numerical methods with the flexibility and automatic resolution adaptivity of particle-based methods. Due to their supersonic nature, Keplerian accretion discs are in principle a very attractive system for applying such freely moving-mesh techniques. However, the high degree of symmetry of simple accretion disc models can be difficult to capture accurately by these methods, due to the generation of geometric grid noise and associated numerical diffusion, which is absent in polar grids. To explore these and other issues, in this work we study the idealized problem of two-dimensional planet-disc interaction with the moving-mesh code AREPO. We explore the hydrodynamic evolution of discs with planets through a series of numerical experiments that vary the planet mass, the disc viscosity and the mesh resolution, and compare the resulting surface density, vortensity field and tidal torque with results from the literature.We find that the performance of the moving-mesh code in this problem is in accordance with published results, showing good consistency with grid codes written in polar coordinates. We also conclude that grid noise and mesh distortions do not introduce excessive numerical diffusion. Finally, we show how the moving-mesh approach can help in resolving an outstanding challenge for polar-coordinate grid codes, namely the successful implementation of adaptive mesh refinement in regions of high density around planets and planetary wakes, while retaining the background flow at low resolution.

KW - Hydrodynamics

KW - Instabilities

KW - Methods: numerical

KW - Planet-disc interactions

KW - Protoplanetary discs

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

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

U2 - 10.1093/mnras/stu1918

DO - 10.1093/mnras/stu1918

M3 - Article

AN - SCOPUS:84923007644

VL - 445

SP - 3475

EP - 3495

JO - Monthly Notices of the Royal Astronomical Society

JF - Monthly Notices of the Royal Astronomical Society

SN - 0035-8711

IS - 4

ER -