Toward an efficient second-order method for computing the surface gravitational potential on spherical-polar meshes
Fuente:
arXiv
Guardado en:
| Autores principales: | , |
|---|---|
| Formato: | Preprint |
| Publicado: |
2024
|
| Materias: | |
| Acceso en línea: | |
| Etiquetas: |
Agregar Etiqueta
Sin Etiquetas, Sea el primero en etiquetar este registro!
|
| _version_ | 1866914872701222912 |
|---|---|
| author | Gressel, Oliver Ziegler, Udo |
| author_facet | Gressel, Oliver Ziegler, Udo |
| contents | Astrophysical accretion discs that carry a significant mass compared with their central object are subject to the effect of self-gravity. In the context of circumstellar discs, this can, for instance, cause fragmentation of the disc gas, and -- under suitable conditions -- lead to the direct formation of gas-giant planets. If one wants to study these phenomena, the disc's gravitational potential needs to be obtained by solving the Poisson equation. This requires to specify suitable boundary conditions. In the case of a spherical-polar computational mesh, a standard multipole expansion for obtaining boundary values is not practicable. We hence compare two alternative methods for overcoming this limitation. The first method is based on a known Green's function expansion (termed "CCGF") of the potential, while the second (termed "James' method") uses a surface screening mass approach with a suitable discrete Green's function. We demonstrate second-order convergence for both methods and test the weak scaling behaviour when using thousands of computational cores. Overall, James' method is found superior owing to its favourable algorithmic complexity of $\sim \mathcal{O}(n^3)$ compared with the $\sim\mathcal{O}(n^4)$ scaling of the CCGF method. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2407_11648 |
| institution | arXiv |
| publishDate | 2024 |
| record_format | arxiv |
| spellingShingle | Toward an efficient second-order method for computing the surface gravitational potential on spherical-polar meshes Gressel, Oliver Ziegler, Udo Earth and Planetary Astrophysics Instrumentation and Methods for Astrophysics Solar and Stellar Astrophysics Astrophysical accretion discs that carry a significant mass compared with their central object are subject to the effect of self-gravity. In the context of circumstellar discs, this can, for instance, cause fragmentation of the disc gas, and -- under suitable conditions -- lead to the direct formation of gas-giant planets. If one wants to study these phenomena, the disc's gravitational potential needs to be obtained by solving the Poisson equation. This requires to specify suitable boundary conditions. In the case of a spherical-polar computational mesh, a standard multipole expansion for obtaining boundary values is not practicable. We hence compare two alternative methods for overcoming this limitation. The first method is based on a known Green's function expansion (termed "CCGF") of the potential, while the second (termed "James' method") uses a surface screening mass approach with a suitable discrete Green's function. We demonstrate second-order convergence for both methods and test the weak scaling behaviour when using thousands of computational cores. Overall, James' method is found superior owing to its favourable algorithmic complexity of $\sim \mathcal{O}(n^3)$ compared with the $\sim\mathcal{O}(n^4)$ scaling of the CCGF method. |
| title | Toward an efficient second-order method for computing the surface gravitational potential on spherical-polar meshes |
| topic | Earth and Planetary Astrophysics Instrumentation and Methods for Astrophysics Solar and Stellar Astrophysics |
| url | https://arxiv.org/abs/2407.11648 |