Revisiting the Numerical Evaluation and Visualization of the Gravity Fields of Asteroid 4769 Castalia Using Polyhedron and Harmonic Expansions Models

For the convenience of comparison with previous literature, the gravity calculations are revisited for the Asteroid 4769 Castalia, but with extensions on its surface and on intersecting planes and spheres around it, using the polyhedron and harmonic expansion methods with different order and degree for different cases, especially including the gravitational accelerations inside the asteroid, which did not appear at all before. In these evaluations, a few different facts of the these methods and results are revealed, such as the fact that gravity diverges when the position radius is less than the mean radius from harmonic-expansion method, and the maximum gravity is not at the deep valley and mountain top. For a surface that intersects the asteroid, the maximum gravity on it is at the intersection lines between the asteroid surface and the spheres or planes. This means that on the sphere and the plane, the gravities inside and outside the asteroid are smaller than the gravity on the intersection, i.e., on the surface. Some analyses of these conclusions are given with many examples with different radii of the sphere and with different order and degree harmonic expansion models for the above asteroid surface and surrounding spheres. It is interesting to note that very few researchers know that the polyhedral method can also be used to calculate the gravity inside an asteroid with just some modifications of the code. Some special gravity figures on surface and planes inside the asteroid Castalia are computed and made for the first time. The calculations also include tangential gravity, potential, and gravitational slope on surface. Specifically, we find that the overall mean gravitational slope could be one kind of indicator of the density of an asteroid. The minimum overall mean slope happens when the asteroid density is about 2.9 g/cm3, which is much larger than a usually assumed value between 1.7 and 2.5 for asteroid Castalia when its period is 4.07 h, since rotation period should be a more accurate parameter than its estimated density. These conclusions about this typical prolate-like asteroid could be a benchmark for analyzing other similar asteroids.