Numerical solution methods for glacial rebound models
نویسندگان
چکیده
We consider the finite element discretization of the system of partial differential equations describing the stress field and the displacements in a (visco)elastic inhomogeneous layered media in response to a surface load. The underlying physical phenomenon, which is modelled, is glacial advance and recession, and the resulting crustal stress state. We analyse the elastic case in more detail and present discretization error estimates. The so-obtained linear system of equations is solved by an iterative solution method with suitable preconditioning and numerical experiments are presented.
منابع مشابه
The Influence of 5000 year-old and Younger Glacial Mass Variability on Present-day Crustal Rebound in the Antarctic Peninsula
Assessment of Antarctic rebound is complicated by two issues: (1) The total ice volume at Last Glacial Maximum is contentious, with estimates ranging from just a few meters to several tens of meters of equivalent eustatic sea level rise. (2) The late Holocene mass budget is also uncertain. Space-based geodesy may provide important data in the coming years for estimating the recent ice mass bala...
متن کاملApplication of Soft Computing Methods for the Estimation of Roadheader Performance from Schmidt Hammer Rebound Values
Estimation of roadheader performance is one of the main topics in determining the economics of underground excavation projects. The poor performance estimation of roadheader scan leads to costly contractual claims. In this paper, the application of soft computing methods for data analysis called adaptive neuro-fuzzy inference system- subtractive clustering method (ANFIS-SCM) and artificial neu...
متن کاملSymplectic and symmetric methods for the numerical solution of some mathematical models of celestial objects
In the last years, the theory of numerical methods for system of non-stiff and stiff ordinary differential equations has reached a certain maturity. So, there are many excellent codes which are based on Runge–Kutta methods, linear multistep methods, Obreshkov methods, hybrid methods or general linear methods. Although these methods have good accuracy and desirable stability properties such as A...
متن کاملAre post-glacial rebound model predictions consistent with the global space-geodetic secular velocity field?
On global scales the present-day secular velocity field of the Earth’s surface is dominated by two processes, namely plate tectonics and post-glacial isostatic rebound. Over the last 15 years, various empirical plate motion models have been derived and constrained by modern space-geodetic observations. Most of these models describe the global horizontal motions for a set of rigid plates distrib...
متن کاملRobust Preconditioned Iterative So- lution Methods for Large-scale Non- symmetric Problems
We study robust, preconditioned, iterative solution methods for largescale linear systems of equations, arising from different applications in geophysics and geotechnics. The first type of linear systems studied here, which are dense, arise from a boundary element type of discretization of crack propagation in brittle material. Numerical experiment show that simple algebraic preconditioning str...
متن کامل