V. Mantle dynamics - A case study

Solid state convection in the rocky mantles is a key to understanding the thermochemical evolution and tectonics of terrestrial planets and moons. It is driven by internal heat and can be described by a system of coupled partial differential equations. There are no analytic solutions for realistic configurations and numerical models are an indispensable tool for researching mantle convection. After a brief general introduction, we introduce the basic equations that govern mantle convection and discuss some common approximations. The following case study is a contribution towards a self-consistent thermochemical evolution model of the Earth. A crude approximation for crustal differentiation is coupled to numerical models of global mantle convection, focussing on geometrical effects and the influence of rheology on stirring. We review Earth-specific geochemical and geophysical constraints, proposals for their reconciliation, and discuss the implications of our models for scenarios of the Earth's evolution. Specific aspects of this study include the use of passive Lagrangian tracers, highly variable viscosity in 3-d spherical geometry, phase boundaries in the mantle and a parameterised model of the core as boundary condition at the bottom of the mantle

Towards realistic simulations of lava dome growth using the level set method

The level set method has been implemented in a computational volcanology context. New techniques are presented to solve the advection equation and the reinitialisation equation. These techniques are based upon an algorithm developed in the finite difference context, but are modified to take advantage of the robustness of the finite element method. The resulting algorithm is tested on a well documented Rayleigh–Taylor instability benchmark [19], and on an axisymmetric problem where the analytical solution is known. Finally, the algorithm is applied to a basic study of lava dome growth

Micromechanical investigation of granular ratcheting using a discrete model of polygonal particles

We use a two-dimensional model of polygonal particles to investigate granular ratcheting. Ratcheting is a long-term response of granular materials under cyclic loading, where the same amount of permanent deformation is accumulated after each cycle. We report on ratcheting for low frequencies and extremely small loading amplitudes. The evolution of the sub-network of sliding contacts allows us to understand the micromechanics of ratcheting. We show that the contact network evolves almost periodically under cyclic loading as the sub-network of the sliding contacts reaches different stages of anisotropy in each cycle. Sliding contacts lead to a monotonic accumulation of permanent deformation per cycle in each particle. The distribution of these deformations appears to be correlated in form of vortices inside the granular assembly