Geodynamic modeling with staggered finite differences and marker in cell: theory, teaching and examples
Taras Gerya, ETH Zurich
Numerical modeling of geodynamic processes is an essential approach in both science and industry with ever- growing demand and high efficiency/cost ratio. Current trend in geodynamic modeling is to develop universal approaches with potentially unlimited number of applications.
One simple and flexible method is based on staggered finite differences and marker in cell techniques (SFD-MIC), which demonstrated superior performance in several branches of modern quantitative Earth sciences. It is suitable for modeling various long-term and short-term thermomechanical processes involving large 3D deformation of rheologically complex materials. Recently, potential applicability of this method to technological processes (material science) and natural processes of industrial significance (geo-hydro-mechanics, waste deposits) has also been demonstrated.
This webinar gives a short theory of the SFD-MIC method, discuses Matlab-based teaching approach and presents modeling examples of natural and technological significance. [slides]
Gabriele Morra, University of Louisiana at Lafayette
Students and young researchers who want to learn to use computational tools for geodynamic modeling have the option to choose among a wide range of numerical tools. I will show how Python and its libraries represent an easy-to-use platform for self-learning, with performance close to compiled codes. I will present (1) how to visualize and run vectorial calculations, (2) examples from classical Mechanics like particles trajectories in 2D-3D, (3) a detailed description of how to write Lagrangian, Eulerian and Particles in Cell codes for solving linear and non-linear continuum mechanics problems and (4) advanced techniques like tree-codes, Boundary Elements, Lattice Boltzmann Method, as well as use Jupyter Notebooks for creating and distributing content. The goal is to encourage professional and students to learn by experimenting and experiencing, like children who learn by playing. [slides][Jupyter Notebooks]
Where have all the dimensions gone? Hands on methods for introducing students to non-dimensional numbers in laboratory and numerical modeling
Eric Mittelstaedt, University of Idaho
Experienced modelers are familiar with how non-dimensionalizing mathematical systems can help improve numerical stability, reduce the number of free variables needed to explain a physical system, capture the essential driving processes of a problem of interest, and scale laboratory experiments to the Earth. However, when first introduced to non-dimensional numbers, students often have difficulty understanding how the mantle can have a depth of 1, or how numbers such as the Nusselt number or the Rayleigh number are derived. In this webinar, I will discuss a hands-on, in-class experiment involving a simple oscillator (mass on a spring) that I have used to introduce students to non-dimensionalizing equations, deriving non-dimensional numbers, and scaling experimental results. The mass-spring system is familiar to many students from their introductory physics classes and the mathematical system is simpler than many problems of interest in geodynamics. The combination of familiarity, basic mathematics, and a hands-on experiment facilitate student comprehension and future application of non-dimensional numbers. These methods are aimed at an introductory graduate course or senior level undergraduate course on modeling.
New developments in AxiSEM/Instaseis for seismic wave propagation on local scales
Lion Krischer, Simon Staehler, and Martin van Driel, ETH Zurich; Tarje Nissen-Meyer, Oxford University
Instaseis (http://instaseis.net) is a Python tool to quickly extract high-frequency seismograms for any source-receiver geometry from databases generated with AxiSEM (http://axisem.info/). So far it has been limited to global datasets, which in turn limited the maximum frequencies due to exploding storage requirements; a number of global databases are hosted by the IRIS DMC (http://ds.iris.edu/ds/products/syngine/). We recently expanded the AxiSEM/Instaseis combination to additionally handle local scale simulations as well as databases, which enables much higher frequencies.
In this webinar, we present a short overview of the theory behind AxiSEM/Instaseis, highlight the new extensions for local scale simulations and databases, and show a short practical tutorial on how to use it.
ASPECT 2.0: Improved architecture, new features
Rene Gassmöller, Juliane Dannberg and John Naliboff, UC Davis