Enregistré dans:
| Auteurs principaux: | , |
|---|---|
| Format: | Preprint |
| Publié: |
2025
|
| Sujets: | |
| Accès en ligne: | https://arxiv.org/abs/2502.18794 |
| Tags: |
Ajouter un tag
Pas de tags, Soyez le premier à ajouter un tag!
|
| _version_ | 1866929732341202944 |
|---|---|
| author | Hanawa, Tomoyuki Mullen, Patrick D. |
| author_facet | Hanawa, Tomoyuki Mullen, Patrick D. |
| contents | High order algorithms have emerged in numerical astrophysics as a promising avenue to reduce truncation error (proportional to a power of the linear resolution $Δx$) with only a moderate increase to computational expense. Significant effort has been placed in the development of finite volume algorithms for (magneto)hydrodynamics, however, state-of-the-art astrophysical simulations tightly couple a plenitude of physics, additionally including gravity, photon transport, cosmic ray transport, chemistry, and/or diffusion, to name a few. Algorithms frequently operator split this additional physics (often a first order error in time) and/or adopt a model wherein their evaluation is limited to second order accuracy in space. In this work, we present a fourth order accurate finite volume scheme for self-gravitating hydrodynamics on a uniform Cartesian grid. The method supplies source terms for the gravitational acceleration ($ρ{\bf g}$) and gravitational energy release ($ρ{\bf v} \cdot {\bf g}$) associated with fourth-order accurate solutions to the Poisson equation. Our scheme (1) guarantees the conservation of total linear momentum, while (2) decreasing (in proportion to $Δx^4$) the effects of spurious heating and/or cooling associated with truncation error in the gravity. We demonstrate expected convergence rates for the algorithm by measuring errors in test problems evolving self-gravity modified linear waves and 3D polytropic equilibria. We test robustness of the algorithm by integrating an induced "inside-out" adiabatic collapse. We also discuss a method to smoothly downgrade the solution to second-order spatial accuracy to avoid spurious overshoots near steep density and/or pressure gradients. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2502_18794 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | Towards Higher Order Accuracy in Self-Gravitating Hydrodynamics Hanawa, Tomoyuki Mullen, Patrick D. Instrumentation and Methods for Astrophysics High order algorithms have emerged in numerical astrophysics as a promising avenue to reduce truncation error (proportional to a power of the linear resolution $Δx$) with only a moderate increase to computational expense. Significant effort has been placed in the development of finite volume algorithms for (magneto)hydrodynamics, however, state-of-the-art astrophysical simulations tightly couple a plenitude of physics, additionally including gravity, photon transport, cosmic ray transport, chemistry, and/or diffusion, to name a few. Algorithms frequently operator split this additional physics (often a first order error in time) and/or adopt a model wherein their evaluation is limited to second order accuracy in space. In this work, we present a fourth order accurate finite volume scheme for self-gravitating hydrodynamics on a uniform Cartesian grid. The method supplies source terms for the gravitational acceleration ($ρ{\bf g}$) and gravitational energy release ($ρ{\bf v} \cdot {\bf g}$) associated with fourth-order accurate solutions to the Poisson equation. Our scheme (1) guarantees the conservation of total linear momentum, while (2) decreasing (in proportion to $Δx^4$) the effects of spurious heating and/or cooling associated with truncation error in the gravity. We demonstrate expected convergence rates for the algorithm by measuring errors in test problems evolving self-gravity modified linear waves and 3D polytropic equilibria. We test robustness of the algorithm by integrating an induced "inside-out" adiabatic collapse. We also discuss a method to smoothly downgrade the solution to second-order spatial accuracy to avoid spurious overshoots near steep density and/or pressure gradients. |
| title | Towards Higher Order Accuracy in Self-Gravitating Hydrodynamics |
| topic | Instrumentation and Methods for Astrophysics |
| url | https://arxiv.org/abs/2502.18794 |