RAM: Rapid Advection Algorithm on Arbitrary Meshes
arXiv:2305.05362 · doi:10.3847/1538-4357/acd698
Abstract
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 up-wind 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 non-trivial geometries. Our method dramatically lowers the computational cost of simulations that simultaneously resolve the global disk and well inside the Hill radius of the secondary companion.
arXiv comment15 pages, 7 figures. Submitted to ApJ. Comments are welcome
References in corpus (12)
- E pur si muove: Galiliean-invariant cosmological hydrodynamical simulations on a moving mesh
- FARGO: a fast eulerian transport algorithm for differentially rotating disks
- Zonal Flows and Long-Lived Axisymmetric Pressure Bumps in Magnetorotational Turbulence
- FARGO3D: A new GPU-oriented MHD code
- Implementation of the Shearing Box Approximation in Athena
- High-order conservative reconstruction schemes for finite volume methods in cylindrical and spherical coordinates
- Computational Eulerian Hydrodynamics and Galilean Invariance
- A conservative orbital advection scheme for simulations of magnetized shear flows with the PLUTO code
- A fourth-order accurate finite volume method for ideal MHD via upwind constrained transport
- DISCO: A 3D Moving-Mesh Magnetohydrodynamics Code Designed for the Study of Astrophysical Disks
- Orbital Advection by Interpolation: A Fast and Accurate Numerical Scheme for Super-Fast MHD Flows
- Orbital advection with magnetohydrodynamics and vector potential