After successful completion of the course, the student is expected to be able to:
Apply the equations of continuum mechanics to geodynamical problems.
Code one- and two-dimensional finite difference and spectral models of linear systems using MATLAB.
Describe the differences and limitations of different numerical techniques and choose the most appropriate method accordingly.
Explain the difference between Dirichlet and Neumann boundary conditions.
Use the commercial code Comsol Multiphysics to solve a variety of dynamical problems.
Produce visualisation of the output (graphs, contour plots, movies, etc.).
Finite difference and spectral methods. Introduction to Comsol Multiphysics (commercial finite element code). Comparison of analytical and numerical solutions: heat equation, flexure of thin and thick plates, fluid flow. Asthenospheric counter-flow model. Post-glacial rebound. Diapirism as Rayleigh-Taylor instability. Thermal convection. One-dimensional flow with constant and variable viscosity. Shear heating in Couette flow and thermal run-away instability. Faulting, friction, and simple earthquake models. Stress diffusion.
Lectures, homework, problem solving, computer exercises.