651-4241-00L  Numerical Modelling I and II: Theory and Applications

SemesterAutumn Semester 2016
LecturersT. Gerya
Periodicityyearly recurring course
Language of instructionEnglish



Courses

NumberTitleHoursLecturers
651-4009-00 GNumerical Modeling I: Theory24s hrs
Mon/108:15-12:00NO F 39 »
T. Gerya
651-4011-00 GNumerical Modeling II: Applications28s hrs
Mon/208:15-12:00NO F 39 »
T. Gerya

Catalogue data

AbstractIn this 13-week sequence, students learn how to write programs from scratch to solve partial differential equations that are useful for Earth science applications. Programming will be done in MATLAB and will use the finite-difference method and marker-in-cell technique. The course will emphasise a hands-on learning approach rather than extensive theory.
ObjectiveThe goal of this course is for students to learn how to program numerical applications from scratch. By the end of the course, students should be able to write state-of-the-art MATLAB codes that solve systems of partial-differential equations relevant to Earth and Planetary Science applications using finite-difference method and marker-in-cell technique. Applications include Poisson equation, buoyancy driven variable viscosity flow, heat diffusion and advection, and state-of-the-art thermomechanical code programming. The emphasis will be on commonality, i.e., using a similar approach to solve different applications, and modularity, i.e., re-use of code in different programs. The course will emphasise a hands-on learning approach rather than extensive theory, and will begin with an introduction to programming in MATLAB.
ContentA provisional week-by-week schedule (subject to change) is as follows:

Week 1: Introduction to the finite difference approximation to differential equations. Introduction to programming in Matlab. Solving of 1D Poisson equation.
Week 2: Direct and iterative methods for obtaining numerical solutions. Solving of 2D Poisson equation with direct method. Solving of 2D Poisson equation with Gauss-Seidel and Jacobi iterative methods.
Week 3: Solving momentum and continuity equations in case of constant viscosity with stream function/vorticity formulation.
Weeks 4: Staggered grid for formulating momentum and continuity equations. Indexing of unknowns. Solving momentum and continuity equations in case of constant viscosity using pressure-velocity formulation with staggered grid.
Weeks 5: Conservative finite differences for the momentum equation. "Free slip" and "no slip" boundary conditions. Solving momentum and continuity equations in case of variable viscosity using pressure-velocity formulation with staggered grid.
Week 6: Advection in 1-D. Eulerian methods. Marker-in-cell method. Comparison of different advection methods and their accuracy.
Week 7: Advection in 2-D with Marker-in-cell method. Combining flow calculation and advection for buoyancy driven flow.
Week 8: "Free surface" boundary condition and "sticky air" approach. Free surface stabilization. Runge-Kutta schemes.
Week 9: Solving 2D heat conservation equation in case of constant thermal conductivity with explicit and implicit approaches.
Week 10: Solving 2D heat conservation equation in case of variable thermal conductivity with implicit approach. Temperature advection with markers. Creating thermomechanical code by combining mechanical solution for 2D buoyancy driven flow with heat diffusion and advection based on marker-in-cell approach.
Week 11: Subgrid diffusion of temperature. Implementing subgrid diffusion to the thermomechanical code.
Week 12: Implementation of radioactive, adiabatic and shear heating to the thermomechanical code.
Week 13: Implementation of temperature-, pressure- and strain rate-dependent viscosity, temperature- and pressure-dependent density and temperature-dependent thermal conductivity to the thermomechanical code. Final project description.


GRADING will be based on weekly programming homeworks (50%) and a term project (50%) to develop an application of their choice to a more advanced level.
LiteratureTaras Gerya, Introduction to Numerical Geodynamic Modelling, Cambridge University Press 2010

Performance assessment

Performance assessment information (valid until the course unit is held again)
Performance assessment as a semester course
ECTS credits6 credits
ExaminersT. Gerya
Typegraded semester performance
Language of examinationEnglish
RepetitionRepetition possible without re-enrolling for the course unit.
Additional information on mode of examinationCourses Numerical modelling I and II are examined together at the end of the semester.

Learning materials

No public learning materials available.
Only public learning materials are listed.

Groups

No information on groups available.

Restrictions

There are no additional restrictions for the registration.

Offered in

ProgrammeSectionType
Earth Sciences MasterGeophysics: Methods IW+Information
Computational Science and Engineering BachelorGeophysics: Subject 2WInformation
Computational Science and Engineering MasterGeophysics: Subject 2WInformation