Abstract
Abstract We present a generalized algorithm based on a spherical harmonics expansion method for efficient computation of the three-dimensional gravitational potential on a multipatch grid in spherical geometry. Instead of solving for the gravitational potential by superposition of separate contributions from the mass density distribution on individual grid patches, our new algorithm directly computes the gravitational potential due to contributions from all grid patches in one computation step, thereby reducing the computational cost of the gravity solver. This is possible by considering a set of angular weights that are derived from rotations of spherical harmonics functions defined in a global coordinate system that is common for all grid patches. Additionally, our algorithm minimizes data communication between parallel computing tasks by eliminating its proportionality to the number of subdomains in the grid configuration, making it suitable for parallelized computation on a multipatch grid configuration with any number of subdomains. Test calculations of the gravitational potential of a triaxial ellipsoidal body with constant mass density on the Yin–Yang two-patch overset grid demonstrate that our method delivers the same level of accuracy as a previous method developed for the Yin–Yang grid while offering improved computation efficiency and parallel scaling behavior.
Sign-in/Register to access full text options 
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callDisclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
