FFT based algorithm for efficient gravity field calculation: comparison with exact results for polyhedral shape models

Manuel Pérez Molina, Adriano Campo Bagatin, Nair Trógolo
2022 Zenodo  
In the frame of the H2020 NEO-MAPP project, we propose a Fast Fourier Transform (FFT) based algorithm for efficient computation of the gravity field created by a body with any arbitrary mass distribution defined by a mass density discretized with suitable resolution. Our algorithm starts by considering a primary three-dimensional cartesian grid that contains the considered mass distribution and has uniform point spacings \(\Delta x\), \(\Delta y\) and \(\Delta z\) in the directions parallel to
more » ... he \(x\), \(y\) and \(z\) axes. The density \(\rho(r)\) is then discretized at each primary grid cell, thus characterizing the mass distribution of the body, which creates the gravity field. Next, our algorithm considers a secondary three-dimensional cartesian grid that applies an arbitrary translation to the primary one and represents the space region where the gravity field will be computed. Once the primary and secondary grids have been defined, our algorithm computes efficiently the components of the gravity vector \(\vec{g}(r)\) created by the (discretized) body within the primary grid at all secondary grid points. From the computed gravity values at such points, a suitable interpolation allows extending the calculation to any point inside the region within the secondary grid. Our algorithm has been applied to bodies having regular and irregular shape as well as uniform and non-uniform density. Using 16 GB of RAM memory in different PCs, we set up primary and secondary grids of up to 241 points at each axis and computation times that do not exceed one minute for each secondary grid. The numerical results provided by our algorithm have been tested to be in very good agreement with exact analytical gravity fields created by homogeneous and non-homogenous mass distributions. In addition, we have applied our algorithm to compute the gravity field created by uniform polyhedral shape models like those of Didymos (65803) and compare the numerical results with the exact analytical ones for homogeneous polyhedrons. From s [...]
doi:10.5281/zenodo.7046033 fatcat:btrcy3bqq5abrmn6nwslz75ity