In this work, the model of rate-and-state friction, which can be viewed as central to the numerical simulation of earthquakes, is considered from a mathematical point of view.
First, a framework is presented through which a general class of such friction laws can be understood and analysed. A prototypical viscoelastic problem of earthquake rupture is then formulated, both in strong and in variational form.
Analysis of this problem is difficult, since the incorporation of rate-and-state friction leads to a coupling of variables. In a time-discrete setting, nonetheless, results on existence, uniqueness, and continuous parameter dependence of solutions can be obtained.
The principal idea is to reformulate the variable interdependence as a fixed point problem and to prove convergence for a corresponding iteration. With that in mind, next, a numerical algorithm is proposed that resolves the coupling through a fixed point iteration. Since it puts a state-of-the-art solver and adaptive time stepping to use, it is not only stable but also fast.
Its applicability to problems of interest is demonstrated in the penultimate chapter, which focuses on simulations of megathrust earthquakes that form at the base of a subduction zone.
The main assumptions made throughout this work are summarised and discussed in the last chapter.
1 Friction: An introduction
1.1 Rate-and-state friction
1.2 The laws of Dieterich and Ruina
1.3 A general framework
2 Continuum mechanics
2.1 Stress and strain
2.3 Conservation laws
2.4 The strong problem
2.5 The weak problem
2.6 Detour: Equivalence of both problems
3 Time discretisation and analysis
3.1 Function spaces
3.2 Solution operators
3.3 Existence and uniqueness of solutions
3.4 Existence in a more general setting
4 Spatial discretisation
5.1 Problem description
5.2 Mathematical modelling
5.4 Comparison of simulation and experiment
5.5 Stability of the algorithm
6 Discussion of assumptions
A Numerical details
A.1 The grid
A.3 Adaptive time stepping
B Miscellaneous theorems
C Gradient flows in one dimension
PDF-Datei von FUDISS_thesis_000000098145
Falls Ihr Browser eine Datei nicht öffnen kann, die Datei zuerst herunterladen und dann öffnen.