Authors
Abstract
The PADÉ 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 nondiffusive and conservative fourth-order Padé differencing, which has higher resolving power compared to both dissipative shock-capturing schemes used in most astrophysics codes, as well as nondiffusive 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.Details
| Publication | The Astrophysical Journal Supplement Series, Volume 273, Issue 2, id.37, 16 pp. |
| Publication Date | August 2024 |
| DOI | |
| arXiv | arXiv:2406.15611 |
| Bibcode | 2024ApJS..273...37S |
| Collections | |
| Keywords | |
| UAT Keywords (generated)BETA | |
| E-Print Comment(s) | Astrophysical Journal Supplement Series, 273:37 (16pp), August 2024; doi:10.3847/1538-4365/ad5af3 |