Numerical methods for glacial isostatic adjustment models = Métodos numéricos para modelos de ajuste isostático glacial
Manifestación
- Autores
-
- Autor: Araújo Cabarcas, Juan Carlos
- Identificador
- 871599
- Fecha de publicación
- 2013
- Forma obra
- Tesis
- Lugar de producción
- 2013
- Idioma
- español
- Nota de edición
- Digitalización realizada por la Biblioteca Virtual del Banco de la República (Colombia)
- Materias
-
- Tecnología; Tecnología / Ingeniería y operaciones afines
- Notas
-
- Colfuturo
- © Derechos reservados del autor
- Nordic countries experience post-glacial rebound, a movement where geographical contours slowly change elevation with respect to the mean sea level. The glacial isostatic adjustment (GIA) model aims to explain the phenomena, which combined with seismic data allows geoscientists to reconstruct elastic coefficients and viscosities of the Earth's lithosphere and upper mantle. The use of standard commercial codes are not adequate for GIA simulations and result in significant errors in the displacement field.
This negative outcome suggests the development of GIA codes that include advection of pre-stress in the model. The problem set up consists on a solid 2D elastic layer under a flat Earth approximation, described by three different models suggested by current studies in geophysics.For space discretization the mixed finite element method (mFEM) is used and efficient preconditioners are built for the resulting algebraic system in saddle point form. A three level GMRES iterative solution strategy is proposed, based on Schur Complement preconditioners coupled with Multigrid techniques.
The implementation is presented as a ready-to-use toolbox that easily deals with problem parameters, geometries, compressible and fully incompressible materials and provides higher accuracy for the displacement field compared with the previously existent codes. It also can be easily extended to 3D geometries and allows the implementation of a viscoelastic mantle. The code is written in C++ using the deal.II library designed for FEM, permitting the use of readily-made software packages, such as Trilinos that are straightforwardly parallelizable. - Complemento de Schur; Finite element method; GIA; Método de elementos finitos; Precondicionadores; Preconditioners; AGI; Schur complement
- Enlace permanente
- https://www.cervantesvirtual.com/obra/numerical-methods-for-glacial-isostatic-adjustment-models-metodos-numericos-para-modelos-de-ajuste-isostatico-glacial-871599
- Enlaces