Pade: A code for protoplanetary disk turbulence based on Pade differencing

Abstract

The Pade code has been developed to treat hydrodynamic turbulence in protoplanetary disks. It solves the compressible equations of motion in cylindrical coordinates. Derivatives are computed using non-diffusive and conservative fourth-order Pade differencing, which has higher resolving power compared to both dissipative shock-capturing schemes used in most astrophysics codes, as well as non-diffusive central finite-difference schemes of the same order. The fourth-order Runge-Kutta method is used for time stepping. A previously reported error-corrected Fargo approach is used to reduce the time step constraint imposed by rapid Keplerian advection. Artificial bulk viscosity is used when shock-capturing is required. Tests for correctness and scaling with respect to the number of processors are presented. Finally, efforts to improve efficiency and accuracy are suggested.

References in corpus (5)

• PLUTO: a Numerical Code for Computational Astrophysics
• Athena: A New Code for Astrophysical MHD
• The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers
• A Vorticity-Preserving Hydrodynamical Scheme for Modeling Accretion Disk Flows
• Investigation of a Vorticity-preserving Scheme for the Euler Equations

Cited by in corpus (1)

• A high resolution simulation of protoplanetary disk turbulence driven by the vertical shear instability