GEO 325M/398M Numerical Modeling in the Geosciences

Jackson School of Geosciences - University of Texas at Austin

Project maintained by mhesse Hosted on GitHub Pages — Theme by mattgraham

Class rooom and time

Office hours

Additional course websites:

Matlab basics:

Here are some LiveScripts I prepared for the first class in 2018 that didn’t have a Matlab prerequisite. If you don’t have much Matlab experience, please look through them. Vectorized programming is a particularly important topic.

  1. demo_arrays.mlx [pdf]
  2. demo_functions.mlx [pdf]
  3. demo_control_flow.mlx [pdf]
  4. demo_matlab_functions.mlx [pdf]
  5. demo_plotting.mlx [pdf]
  6. demo_vectorized_programing.mlx [pdf]
  7. demo_odds_ends.mlx [pdf] (structures, logical indexing, anonymous functions)

Below are two files that I have sometimes used for the demos in class. If you put them into the folder with class files you should have no problem.

This years course project

In spring 2025 we will develop a model for the convection in high-pressure ice in icy ocean worlds that may controll the transfer of nutrients from the rocky interior into the ocean sandwiched between the layers of high and low pressure ice.

hi

For reference see [Vance at al. 2020]

Introduction

Lecture 1 (Jan 14): Course Project and Conservation Laws

Lecture 2 (Jan 16): Introduction to numerics

Lecture 3 (Jan 21): Conservative Finite Differences

1D Heat Flow

Lecture 4 (Jan 23): Discrete Operators

Lecture 5 (Jan 28): Boundary Conditions I: Dirichlet homogeneous

Lecture 6 (Jan 30): Boundary Conditions II: Dirichlet heterogenous

Lecture 7 (Feb 4): Boundary Conditons III: Neumann

Lecture 8 (Feb 6): Heterogenous coefficients and layered media

Lecture 9 (Feb 11): Non-linear problems and Netwon-Raphson method

Lecture 10 (Feb 13): Numerical Jscobian

Lecture 11 (Feb 18): Analytic Jacobian

Lecture 12 (Feb 20): Transient heat transport

1D ADVECTIVE-CONDUCTIVE HEAT TRANSFER

Lecture 13 (Feb 25): Advective heat transport

Lecture 14 (Feb 27): Time stepping ADE

DISCRETIZATION IN 2D

Lecture 15 (Mar 4): 2D Discrete operators - Part I

Lecture 16 (Mar 6): 2D Discrete operators - Part II

No class Mar 11!

Lecture 17 (Mar 13): 2D Advection matrix - Cooling of oceanic plates

STOKES FLOW

Lecture 18 (Mar 25): Derivation of Stokes equation

Lecture 19 (Mar 27): Staggered grid for Stress/Strain

no class Apr 1:

Lecture 20 (Apr 3): Discrete Stokes operators and BC’s

Lecture 21 (Apr 8): Streamfunction

Lecture 22 (Apr 10): Coupled Stokes flow and heat transport - Mid-ocean ridges

Lecture 23 (Apr 15): Rayleigh-Bernard convection

Lecture 24 (Apr 17): Implementing convection

Lecture 25 (Apr 22):

Lecture 26 (Apr 24): Class review & Discussion