scholarly journals Geodesic equations and their numerical solution in Cartesian coordinates on a triaxial ellipsoid

2019 ◽  
Vol 9 (1) ◽  
pp. 1-12 ◽  
Author(s):  
G. Panou ◽  
R. Korakitis

Abstract In this work, the geodesic equations and their numerical solution in Cartesian coordinates on an oblate spheroid, presented by Panou and Korakitis (2017), are generalized on a triaxial ellipsoid. A new exact analytical method and a new numerical method of converting Cartesian to ellipsoidal coordinates of a point on a triaxial ellipsoid are presented. An extensive test set for the coordinate conversion is used, in order to evaluate the performance of the two methods. The direct geodesic problem on a triaxial ellipsoid is described as an initial value problem and is solved numerically in Cartesian coordinates. The solution provides the Cartesian coordinates and the angle between the line of constant λ and the geodesic, at any point along the geodesic. Also, the Liouville constant is computed at any point along the geodesic, allowing to check the precision of the method. An extensive data set of geodesics is used, in order to demonstrate the validity of the numerical method for the geodesic problem. We conclude that a complete, stable and precise solution of the problem is accomplished.

2020 ◽  
Vol 14 (2) ◽  
pp. 205-213
Author(s):  
G. Panou ◽  
R. Korakitis

AbstractIn this work, the direct geodesic problem in Cartesian coordinates on a triaxial ellipsoid is solved by an approximate analytical method. The parametric coordinates are used and the parametric to Cartesian coordinates conversion and vice versa are presented. The geodesic equations on a triaxial ellipsoid in Cartesian coordinates are solved using a Taylor series expansion. The solution provides the Cartesian coordinates and the angle between the line of constant v and the geodesic at the end point. An extensive data set of geodesics, previously studied with a numerical method, is used in order to validate the presented analytical method in terms of stability, accuracy and execution time. We conclude that the presented method is suitable for a triaxial ellipsoid with small eccentricities and an accurate solution is obtained. At a similar accuracy level, this method is about thirty times faster than the corresponding numerical method. Finally, the presented method can also be applied in the degenerate case of an oblate spheroid, which is extensively used in geodesy.


2017 ◽  
Vol 7 (1) ◽  
Author(s):  
G. Panou ◽  
R. Korakitis

AbstractThe direct geodesic problem on an oblate spheroid is described as an initial value problem and is solved numerically using both geodetic and Cartesian coordinates. The geodesic equations are formulated by means of the theory of differential geometry. The initial value problem under consideration is reduced to a system of first-order ordinary differential equations, which is solved using a numerical method. The solution provides the coordinates and the azimuths at any point along the geodesic. The Clairaut constant is not used for the solution but it is computed, allowing to check the precision of the method. An extensive data set of geodesics is used, in order to evaluate the performance of the method in each coordinate system. The results for the direct geodesic problem are validated by comparison to Karney’s method. We conclude that a complete, stable, precise, accurate and fast solution of the problem in Cartesian coordinates is accomplished.


2021 ◽  
Vol 11 (1) ◽  
pp. 122-130
Author(s):  
C. Marx

Abstract The direct geodetic problem is regarded on the biaxial and triaxial ellipsoid. A known solution method suitable for low eccentricities, which uses differential equations in Cartesian coordinates and Taylor series expansions of these coordinates, is advanced in view of its practical application. According to previous works, this approach has the advantages that no singularities occur in the determination of the coordinates, its mathematical formulation is simple and it is not computationally intensive. The formulas of the solution method are simplified in the present contribution. A test of this method using an extensive test data set on a biaxial earth ellipsoid shows its accuracy and practicability for distances of any length. Based on the convergence behavior of the series of the test data set, a truncation criterion for the series expansions is compiled taking into account accuracy requirements of the coordinates. Furthermore, a procedure is shown which controls the truncation of the series expansions by accuracy requirements of the direction to be determined in the direct problem. The conducted tests demonstrate the correct functioning of the methods for the series truncation. However, the considered solution method turns out to be significantly slower than another current method for biaxial ellipsoids, which makes it more relevant for triaxial ellipsoids.


2021 ◽  
Vol 11 (1) ◽  
pp. 111-121
Author(s):  
G. Panou ◽  
R. Korakitis

Abstract In this work, two analytical and two numerical methods of converting Cartesian to ellipsoidal coordinates of a point in space are presented. After slightly modifying a well-known exact analytical method, a new exact analytical method is developed. Also, two well-known numerical methods, which were developed for points exactly on the surface of a triaxial ellipsoid, are generalized for points in space. The four methods are validated with numerical experiments using an extensive set of points for the case of the Earth. Then, a theoretical and a numerical comparative assessment of the four methods is made. Furthermore, the new exact analytical method is applied for an almost oblate spheroid and for the case of the Moon and the results are compared. We conclude that, the generalized Panou and Korakitis’ numerical method, starting with approximate values from the new exact analytical method, is the best choice in terms of accuracy of the resulting ellipsoidal coordinates.


2018 ◽  
Vol 2018 ◽  
pp. 1-7
Author(s):  
Yong-Hong Fan ◽  
Lin-Lin Wang

We propose a new algorithm for solving the terminal value problems on a q-difference equations. Through some transformations, the terminal value problems which contain the first- and second-order delta-derivatives have been changed into the corresponding initial value problems; then with the help of the methods developed by Liu and H. Jafari, the numerical solution has been obtained and the error estimate has also been considered for the terminal value problems. Some examples are given to illustrate the accuracy of the numerical methods we proposed. By comparing the exact solution with the numerical solution, we find that the convergence speed of this numerical method is very fast.


2012 ◽  
Vol 263-266 ◽  
pp. 1315-1318
Author(s):  
Kun Ming Yu ◽  
Ming Gong Lee

This paper is to discuss how Python can be used in designing a cluster parallel computation environment in numerical solution of some block predictor-corrector method for ordinary differential equations. In the parallel process, MPI-2(message passing interface) is used as a standard of MPICH2 to communicate between CPUs. The operation of data receiving and sending are operated and controlled by mpi4py which is based on Python. Implementation of a block predictor-corrector numerical method with one and two CPUs respectively is used to test the performance of some initial value problem. Minor speed up is obtained due to small size problems and few CPUs used in the scheme, though the establishment of this scheme by Python is valuable due to very few research has been carried in this kind of parallel structure under Python.


1990 ◽  
Vol 220 ◽  
pp. 339-354 ◽  
Author(s):  
S. D. Heister ◽  
J. M. Mcdonough ◽  
A. R. Karagozian ◽  
D. W. Jenkins

A numerical solution for the flow field associated with a compressible pair of counter-rotating vortices is developed. The compressible, two-dimensional potential equation is solved utilizing the numerical method of Osher et al. (1985) for flow regions in which a non-zero density exists. Close to the vortex centres, vacuum ‘cores’ develop owing to the existence of a maximum achievable flow speed in a compressible flow field. A special treatment is required to represent these vacuum cores. Typical streamline patterns and core boundaries are obtained for upstream Mach numbers as high as 0.3, and the formation of weak shocks, predicted by Moore & Pullin (1987), is observed.


1973 ◽  
Vol 187 (1) ◽  
pp. 71-78 ◽  
Author(s):  
B. R. Reason ◽  
D. Dyer

We present a numerical solution for the operating conditions of a hydrodynamic porous journal bearing. The numerical method allows for the possibility of variable porosity in the bearing matrix, but the solution has been achieved on the assumption of matrix homogeneity. The relation between the various bearing parameters have been shown for a variety of bearing geometries and permeabilities enabling the operating conditions for this type of bearing to be better appreciated. A comparison of the present solution with approximate solutions used by other authors has been made, which indicates the useful working range of the approximate solutions.


2013 ◽  
Vol 2013 ◽  
pp. 1-7 ◽  
Author(s):  
Berna Bülbül ◽  
Mehmet Sezer

A numerical method based on collocation points is developed to solve the nonlinear Klein-Gordon equations by using the Taylor matrix method. The method is applied to some test examples and the numerical results are compared with the exact solutions. The results reveal that the method is very effective, simple, and convenient. In addition, an error estimation of proposed method is presented.


Author(s):  
Qani Yalda

The main purpose of this paper is to obtain the real roots of an expression using the Numerical method, bisection method, Newton's method and secant method. Root analysis is calculated using specific, precise starting points and numerical methods and is represented by Maple. In this research, we used Maple software to analyze the roots of nonlinear equations by special methods, and by showing geometric diagrams, we examined the relevant examples. In this process, the Newton-Raphson method, the algorithm for root access, is fully illustrated by Maple. Also, the secant method and the bisection method were demonstrated by Maple by solving examples and drawing graphs related to each method.


Sign in / Sign up

Export Citation Format

Share Document