arrow
Return

Gravity from surface triangulation: convergence acceleration with nested grids

delete2026-03-22
delete0
PRE
AI
DOI:10.1007/s10569-026-10286-2delete
deleteOriginal
deleteOriginal request for help
deleteShare
deleteSave
Abstract

Abstract

En 中文
The determination of the gravitational potential by the polyhedral method is revisited in the case where the surface of a body is composed of triangular facets. Based upon six test-shapes of astrophysical interest (sphere, spheroid, triaxial, lemon-shape, dumbell and torus) projected on nested grids, we verify that the convergence toward reference values is second-order in the step size of the grid, inside the body, at the surface and outside. We then show that the accuracy or computing time can be drastically enhanced by implementing the Repeated Richardson Extrapolation. This technique is especially efficient when the body's surface is smooth enough, and is therefore well adapted to the theory of figures (single and multi-layer fluids) and to dynamical studies (test-particle and mutual interactions), which require a large number of field evaluations. For real objects like asteroids that have very irregular terrains at small scales, the gain is modest. In that context, we estimate the discretization level beyond which the typical error in potential values due to altimetric uncertainties dominates over the contribution of sub-grid cavities and bumps. For bodies close to spherical, the criterion reads T greater than or similar to 64D3 lambda,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T \gtrsim \frac{64 D}{3 \lambda },$$\end{document} where D is the diameter of the body, lambda\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda $$\end{document} the typical shape error and T the number of triangular facets involved. The case of 433 Eros is considered as an example.
Keywords:
Analytical methods
Numerical methods
Gravitational potential
433 Eros

Journal

C
CELESTIAL MECHANICS & DYNAMICAL ASTRONOMY
IF:
1.4
Papers:
36
Citations:
0

Organization

C
centre national de la recherche scientifique (cnrs)
Scholars:
24.5W
Papers: 18.2W
Citations: 279