RAM: Rapid Advection Algorithm on Arbitrary Meshes

• Published in: The Astrophysical journal Vol. 952; no. 2; pp. 106 - 118
• Main Authors: Benítez-Llambay, Pablo, Krapp, Leonardo, Ramos, Ximena S., Kratter, Kaitlin M.
• Format: Journal Article
• Language: English
• Published: Philadelphia The American Astronomical Society 01.08.2023; IOP Publishing
• Subjects: Accretion; Advection; Algorithms; Astrophysics; Computational grids; Computing costs; Fluid mechanics; High Mach number; Hydrodynamics; Interpolation; Protoplanetary disks; Truncation errors; Velocity
• ISSN: 0004-637X; 1538-4357; 1538-4357
• DOI: 10.3847/1538-4357/acd698

The study of many astrophysical flows requires computational algorithms that can capture high Mach number flows, while resolving a large dynamic range in spatial and density scales. In this paper we present a novel method, RAM: Rapid Advection Algorithm on Arbitrary Meshes. RAM is a time-explicit method to solve the advection equation in problems with large bulk velocity on arbitrary computational grids. In comparison with standard upwind algorithms, RAM enables advection with larger time steps and lower truncation errors. Our method is based on the operator splitting technique and conservative interpolation. Depending on the bulk velocity and resolution, RAM can decrease the numerical cost of hydrodynamics by more than one order of magnitude. To quantify the truncation errors and speed-up with RAM, we perform one- and two-dimensional hydrodynamics tests. We find that the order of our method is given by the order of the conservative interpolation and that the effective speed-up is in agreement with the relative increment in time step. RAM will be especially useful for numerical studies of disk−satellite interaction, characterized by high bulk orbital velocities and nontrivial geometries. Our method dramatically lowers the computational cost of simulations that simultaneously resolve the global disk and potential well inside the Hill radius of the secondary companion.