Centre de documentation
scientifique

Numerical solution to the oblique derivative boundary value problem on non-uniform grids above the Earth topography / Matej Medl’a in Journal of geodesy, vol 92 n° 1 (January 2018)  

Public [article]  inJournal of geodesy > vol 92 n° 1 (January 2018) . - pp 1 - 19 Titre : Numerical solution to the oblique derivative boundary value problem on non-uniform grids above the Earth topography Type de document : Article/Communication Auteurs : Matej Medl’a, Auteur ; Karol Mikula, Auteur ; Robert Cunderlik, Auteur ; Marek Macák, Auteur Année de publication : 2018 Article en page(s) : pp 1 - 19 Note générale : Bibliographie Langues : Anglais (eng) Descripteur : [Vedettes matières IGN] Géodésie physique [Termes IGN] champ de pesanteur local [Termes IGN] discrétisation [Termes IGN] équation de Laplace [Termes IGN] méthode des éléments finis [Termes IGN] problème des valeurs limites Résumé : (Auteur) The paper presents a numerical solution of the oblique derivative boundary value problem on and above the Earth’s topography using the finite volume method (FVM). It introduces a novel method for constructing non-uniform hexahedron 3D grids above the Earth’s surface. It is based on an evolution of a surface, which approximates the Earth’s topography, by mean curvature. To obtain optimal shapes of non-uniform 3D grid, the proposed evolution is accompanied by a tangential redistribution of grid nodes. Afterwards, the Laplace equation is discretized using FVM developed for such a non-uniform grid. The oblique derivative boundary condition is treated as a stationary advection equation, and we derive a new upwind type discretization suitable for non-uniform 3D grids. The discretization of the Laplace equation together with the discretization of the oblique derivative boundary condition leads to a linear system of equations. The solution of this system gives the disturbing potential in the whole computational domain including the Earth’s surface. Numerical experiments aim to show properties and demonstrate efficiency of the developed FVM approach. The first experiments study an experimental order of convergence of the method. Then, a reconstruction of the harmonic function on the Earth’s topography, which is generated from the EGM2008 or EIGEN-6C4 global geopotential model, is presented. The obtained FVM solutions show that refining of the computational grid leads to more precise results. The last experiment deals with local gravity field modelling in Slovakia using terrestrial gravity data. The GNSS-levelling test shows accuracy of the obtained local quasigeoid model. Numéro de notice : A2018-011 Affiliation des auteurs : non IGN Thématique : POSITIONNEMENT Nature : Article nature-HAL : ArtAvecCL-RevueIntern DOI : 10.1007/s00190-017-1040-z Date de publication en ligne : 30/05/2017 En ligne : https://doi.org/10.1007/s00190-017-1040-z Format de la ressource électronique : URL Article Permalink : https://documentation.ensg.eu/index.php?lvl=notice_display&id=89054 [article]