Centre de documentation
scientifique

Fast calculation of gravitational effects using tesseroids with a polynomial density of arbitrary degree in depth / Fang Ouyang in Journal of geodesy, vol 96 n° 12 (December 2022)  

Public [article]  inJournal of geodesy > vol 96 n° 12 (December 2022) . - n° 97 Titre : Fast calculation of gravitational effects using tesseroids with a polynomial density of arbitrary degree in depth Type de document : Article/Communication Auteurs : Fang Ouyang, Auteur ; Long-wei Chen, Auteur ; Zhi-gang Shao, Auteur Année de publication : 2022 Article en page(s) : n° 97 Note générale : bibliographie Langues : Anglais (eng) Descripteur : [Vedettes matières IGN] Géodésie physique [Termes IGN] champ de gravitation [Termes IGN] coordonnées sphériques [Termes IGN] discrétisation [Termes IGN] intégrale de Newton [Termes IGN] inversion [Termes IGN] quadrature [Termes IGN] tesseroid [Termes IGN] transformation rapide de Fourier Résumé : (auteur) Fast and accurate calculation of gravitational effects on a regional or global scale with complex density environment is a critical issue in gravitational forward modelling. Most existing significant developments with tessroid-based modelling are limited to homogeneous density models or polynomial ones of a limited order. Moreover, the total gravitational effects of tesseroids are often calculated by pure summation in these methods, which makes the calculation extremely time-consuming. A new efficient and accurate method based on tesseroids with a polynomial density up to an arbitrary order in depth is developed for 3D large-scale gravitational forward modelling. The method divides the source region into a number of tesseroids, and the density in each tesseroid is assumed to be a polynomial function of arbitrary degree. To guarantee the computational accuracy and efficiency, two key points are involved: (1) the volume Newton’s integral is decomposed into a one-dimensional integral with a polynomial density in the radial direction, for which a simple analytical recursive formula is derived for efficient calculation, and a surface integral over the horizontal directions evaluated by the Gauss–Legendre quadrature (GLQ) combined with a 2D adaptive discretization strategy; (2) a fast and flexible discrete convolution algorithm based on 1D fast Fourier transform (FFT) and a general Toepritz form of weight coefficient matrices is adopted in the longitudinal dimension to speed up the computation of the cumulative contributions from all tesseroids. Numerical examples show that the gravitational fields predicted by the new method have a good agreement with the corresponding analytical solutions for spherical shell models with both polynomial and non-polynomial density variations in depth. Compared with the 3D GLQ methods, the new algorithm is computationally more accurate and efficient. The calculation time is significantly reduced by 3 orders of magnitude as compared with the traditional 3D GLQ methods. Application of the new algorithm in the global crustal CRUST1.0 model further verifies its reliability and practicability in real cases. The proposed method will provide a powerful numerical tool for large-scale gravity modelling and also an efficient forward engine for inversion and continuation problems. Numéro de notice : A2022-896 Affiliation des auteurs : non IGN Thématique : POSITIONNEMENT Nature : Article DOI : 10.1007/s00190-022-01688-9 Date de publication en ligne : 05/12/2022 En ligne : https://doi.org/10.1007/s00190-022-01688-9 Format de la ressource électronique : URL article Permalink : https://documentation.ensg.eu/index.php?lvl=notice_display&id=102248 [article]

Gravitational field modelling near irregularly shaped bodies using spherical harmonics: a case study for the asteroid (101955) Bennu / Blažej Bucha in Journal of geodesy, vol 95 n° 5 (May 2021)  

Public [article]  inJournal of geodesy > vol 95 n° 5 (May 2021) . - n° 56 Titre : Gravitational field modelling near irregularly shaped bodies using spherical harmonics: a case study for the asteroid (101955) Bennu Type de document : Article/Communication Auteurs : Blažej Bucha, Auteur ; Fernando Sanso, Auteur Année de publication : 2021 Article en page(s) : n° 56 Note générale : bibliographie Langues : Anglais (eng) Descripteur : [Vedettes matières IGN] Géodésie physique [Termes IGN] astéroïde [Termes IGN] champ de gravitation [Termes IGN] champ de pesanteur terrestre [Termes IGN] convergence [Termes IGN] harmonique sphérique [Termes IGN] intégrale de Newton Résumé : (auteur) We apply three spherical-harmonic-based techniques to deliver external gravitational field models of the asteroid (101955) Bennu within its circumscribing sphere. This region is known to be peculiar for external spherical harmonic expansions, because it may lead to a divergent series. The studied approaches are (i) spectral gravity forward modelling via external spherical harmonics, (ii) the least-squares estimation from surface gravitational data using external spherical harmonics and (iii) the combination of internal and external series expansions. While the first method diverges beyond any reasonable doubts, we show that the other two methods may ensure relative accuracy from ∼10−6 to 10−8 in the vicinity of Bennu. This is possible even with the second method, despite the fact that it relies on a single series of external spherical harmonics. Our main motivation was to study conceptual differences between spherical harmonic coefficients from satellite data (analogy to the first method) and from surface gravitational data (the second method). Such coefficients are available through the popular spherical-harmonic-based models of the Earth’s gravitational field and often are combined together. We show that the coefficients from terrestrial data may lead to a divergence effect of partial sums, though excellent accuracy can be achieved when such model is used in full. Under (presently) extreme but realistic conditions, the divergence effect of partial sums may affect many near-surface geoscientific applications, such as the geoid/quasigeoid computation or residual terrain modelling. Computer codes (Fortran, MATLAB) and data produced within the study are made freely available at http://edisk.cvt.stuba.sk/~xbuchab/. Numéro de notice : A2021-347 Affiliation des auteurs : non IGN Thématique : POSITIONNEMENT Nature : Article nature-HAL : ArtAvecCL-RevueIntern DOI : 10.1007/s00190-021-01493-w Date de publication en ligne : 22/04/2021 En ligne : https://doi.org/10.1007/s00190-021-01493-w Format de la ressource électronique : URL article Permalink : https://documentation.ensg.eu/index.php?lvl=notice_display&id=97591 [article]

Optimized formulas for the gravitational field of a tesseroid / Thomas Grombein in Journal of geodesy, vol 87 n° 7 (July 2013)  

Public [article]  inJournal of geodesy > vol 87 n° 7 (July 2013) . - pp 645 - 660 Titre : Optimized formulas for the gravitational field of a tesseroid Type de document : Article/Communication Auteurs : Thomas Grombein, Auteur ; Kurt Seitz, Auteur ; Bernard Heck, Auteur Année de publication : 2013 Article en page(s) : pp 645 - 660 Note générale : Bibliographie Langues : Anglais (eng) Descripteur : [Vedettes matières IGN] Géodésie physique [Termes IGN] champ de pesanteur terrestre [Termes IGN] coordonnées cartésiennes géocentriques [Termes IGN] coordonnées sphériques [Termes IGN] intégrale de Newton [Termes IGN] tesseroid Résumé : (Auteur) Various tasks in geodesy, geophysics, and related geosciences require precise information on the impact of mass distributions on gravity field-related quantities, such as the gravitational potential and its partial derivatives. Using forward modeling based on Newton’s integral, mass distributions are generally decomposed into regular elementary bodies. In classical approaches, prisms or point mass approximations are mostly utilized. Considering the effect of the sphericity of the Earth, alternative mass modeling methods based on tesseroid bodies (spherical prisms) should be taken into account, particularly in regional and global applications. Expressions for the gravitational field of a point mass are relatively simple when formulated in Cartesian coordinates. In the case of integrating over a tesseroid volume bounded by geocentric spherical coordinates, it will be shown that it is also beneficial to represent the integral kernel in terms of Cartesian coordinates. This considerably simplifies the determination of the tesseroid’s potential derivatives in comparison with previously published methodologies that make use of integral kernels expressed in spherical coordinates. Based on this idea, optimized formulas for the gravitational potential of a homogeneous tesseroid and its derivatives up to second-order are elaborated in this paper. These new formulas do not suffer from the polar singularity of the spherical coordinate system and can, therefore, be evaluated for any position on the globe. Since integrals over tesseroid volumes cannot be solved analytically, the numerical evaluation is achieved by means of expanding the integral kernel in a Taylor series with fourth-order error in the spatial coordinates of the integration point. As the structure of the Cartesian integral kernel is substantially simplified, Taylor coefficients can be represented in a compact and computationally attractive form. Thus, the use of the optimized tesseroid formulas particularly benefits from a significant decrease in computation time by about 45 % compared to previously used algorithms. In order to show the computational efficiency and to validate the mathematical derivations, the new tesseroid formulas are applied to two realistic numerical experiments and are compared to previously published tesseroid methods and the conventional prism approach. Numéro de notice : A2013-399 Affiliation des auteurs : non IGN Thématique : POSITIONNEMENT Nature : Article nature-HAL : ArtAvecCL-RevueIntern DOI : 10.1007/s00190-013-0636-1 Date de publication en ligne : 18/05/2013 En ligne : https://doi.org/10.1007/s00190-013-0636-1 Format de la ressource électronique : URL article Permalink : https://documentation.ensg.eu/index.php?lvl=notice_display&id=32537 [article]

Exemplaires(1)

Code-barres	Cote	Support	Localisation	Section	Disponibilité
266-2013071	SL	Revue	Centre de documentation	Revues en salle	Disponible