On the integral solution of hyperbolic Kepler's equation

M. Calvo,A. Elipe,L. Rández

DOI: https://doi.org/10.1007/s10569-024-10184-5

2024-03-30

Celestial Mechanics and Dynamical Astronomy

Abstract:In a recent paper of Philcox, Goodman and Slepian, the solution of the elliptic Kepler's equation is given as a quotient of two contour integrals along a Jordan curve that contains in its interior the unique real solution of the elliptic Kepler's equation and does not include other complex zeroes. In this paper, we show that a similar explicit integral solution can be given for the hyperbolic Kepler's equation. With this purpose, we carry out a study of the complex zeros of the hyperbolic Kepler's equation in order to define suitable Jordan contours in the integrals. Even more, we show that appropriate elliptic Jordan contours can be defined for such integrals, which reduces the computing time. Moreover, using the ideas behind the fast Fourier transform (FFT) algorithm, these integrals can be approximated by the composite trapezoidal rule which gives an algorithm with spectral convergence as a function of the number of nodes. The results of some numerical experiments are presented to show that this implementation is a reliable and very accurate algorithm for solving the hyperbolic Kepler's equation.

astronomy & astrophysics,mathematics, interdisciplinary applications

What problem does this paper attempt to address?

This paper aims to solve the problem of integral solutions of the Hyperbolic Kepler's Equation (HKE). Specifically, the goal of the paper is to provide an explicit analytical solution through integral methods and study the effectiveness and accuracy of its numerical implementation. The following are the main contributions of the paper:

1. **Study of complex zeros**: The paper first conducts a comprehensive study of the complex zeros of HKE to define the appropriate contour (Jordan curve) in the integral. This step is crucial for ensuring the correctness and effectiveness of the integral solution.

2. **Definition of elliptical contour**: The paper proposes using an elliptical contour instead of a circular contour to define the path of the integral. Compared with the circular contour, the elliptical contour is more computationally efficient and can improve the accuracy of the solution.

3. **Approximation method of integral solution**: The paper shows how to use the Composite Trapezoidal Rule to approximate the integral solution. This method not only has spectral convergence but can also further simplify the calculation process by taking advantage of the symmetry of the function, thereby improving the computational efficiency.

4. **Numerical experiments**: The paper verifies the effectiveness and accuracy of the proposed method through a series of numerical experiments. The experimental results show that even in the case of approaching singular points, this method can achieve high accuracy, and as the number of nodes increases, the accuracy improves significantly.

5. **Conclusion**: The conclusion part of the paper summarizes the above research results, pointing out that this integral solution method is not only an elegant expression in theory but also can provide high - precision solutions in practical applications while maintaining a reasonable computational cost.

In summary, through theoretical analysis and numerical experiments, this paper successfully provides an efficient and accurate method for integral solutions of the hyperbolic Kepler's equation.