On the application of the coupled finite-infinite element method to geodetic boundary-value problem

Our aim is to introduce the Coupled Finite-Infinite Element Method (CFIEM) as a new alternative approach to the Earth's gravity field modelling. We show that if the computational domain is large enough in radial direction, one can obtain the qualitatively and quantitatively comparable solution to the solution by the Finite Element Method (FEM). We study the influence of the size of the computational domain on the final CFIEM solution as well as the successive refinement of the discretization and its convergence to the exact solution. As an input data we use the synthetic boundary conditions computed from a Synthetic Earth Gravity Model (SEGM) and we test the CFIEM solution by the data generated directly from SEGM and the solution by the FEM.