Solving the fixed gravimetric boundary value problem by the finite element method using mapped infinite elements

The numerical approach for solving the fixed gravimetric boundary value problem (FGBVP) based on the finite element method (FEM) with mapped infinite elements is developed and implemented. In this approach, the 3D semi-infinite domain outside the Earth is bounded by the triangular discretization of...

Celý popis

Uložené v:
Podrobná bibliografia
Vydané v:Computational geosciences Ročník 27; číslo 4; s. 649 - 662
Hlavní autori: Macák, Marek, Minarechová, Zuzana, Tomek, Lukáš, Čunderlík, Róbert, Mikula, Karol
Médium: Journal Article
Jazyk:English
Vydavateľské údaje: Cham Springer International Publishing 01.08.2023
Springer Nature B.V
Predmet:
Boundary conditions
Boundary value problems
Computer applications
Dirichlet problem
Discretization
Earth
Earth and Environmental Science
Earth Sciences
Earth surface
Exact solutions
Finite element analysis
Finite element method
Geotechnical Engineering & Applied Earth Sciences
Harmonic functions
Hydrogeology
Infinite elements
Infinity
Laplace equation
Mathematical Modeling and Industrial Mathematics
Mathematical models
Nodes
Numerical methods
Original Paper
Prisms
Soil Science & Conservation
Finite element method
Oblique derivative boundary condition
3D semi-infinite domain
Fixed gravimetric boundary value problem
Reconstruction of the EGM2008
Mapped infinite element
Global gravity field modelling
ISSN:1420-0597, 1573-1499
On-line prístup:Získať plný text
Tagy: Pridať tag
Žiadne tagy, Buďte prvý, kto otaguje tento záznam!
Popis
Shrnutí:The numerical approach for solving the fixed gravimetric boundary value problem (FGBVP) based on the finite element method (FEM) with mapped infinite elements is developed and implemented. In this approach, the 3D semi-infinite domain outside the Earth is bounded by the triangular discretization of the whole Earth’s surface and extends to infinity. Then the FGBVP consists of the Laplace equation for unknown disturbing potential which holds in the domain, the oblique derivative boundary condition (BC) given directly at computational nodes on the Earth’s surface, and regularity of the disturbing potential at infinity. In this way, it differs from previous FEM approaches, since the numerical solution is not fixed by the Dirichlet BC on some part of the boundary of the computational domain. As a numerical method, the FEM with finite and mapped infinite triangular prisms has been derived and implemented. In experiments, at first, a convergence of the proposed numerical scheme to the exact solution is tested. Afterwards, a numerical study is focused on a reconstruction of the harmonic function (EGM2008) above the Earth’s topography. Here, a special discretization of the Earth’s surface which is able to fulfil the conditions that arise from correct geometrical properties of finite elements, and it is suitable for parallel computing is implemented. The obtained solutions at nodes on the Earth’s surface as well as nodes that lie approximately at the altitude of the GOCE satellite mission have been tested.
Bibliografia:ObjectType-Article-1
SourceType-Scholarly Journals-1
ObjectType-Feature-2
content type line 14
ISSN:1420-0597
1573-1499
DOI:10.1007/s10596-023-10224-3