When and where:


Office hours: TBD 

Description: The modeling of many problems in the applied sciences and engineering is based on concepts of conservation of mass, momentum and energy, leading to systems of conservation laws. Prominent examples are the Maxwell equations of electromagnetics, the Euler and Navier-Stokes equations of fluid dynamics and equations of elasticity and the systems of magnetohydrodynamics of plasma physics.

In this course we shall develop, analyze and apply computational methods suitable for solving systems of conservation laws. We shall begin to discussing fundamental properties of conservation laws, including their ability to generate non-smooth solutions – shocks – from smooth initial conditions, leading to the introduction of weak solutions and entropy conditions.

After an initial discussion of finite difference methods for conservation laws, we introduce finite volume methods as the first major class of methods to study, including accuracy and stability of these methods. We discuss the importance of the numerical flux and approximate Riemann solvers and the extension of finite volume methods to general grids.

For larger and more complex problems, the ability to increase the order of the method is important and we discuss such extensions and the new challenges these introduce. This sets the stage for the development of nonlinear schemes, primarily of the type of essentially non-oscillatory (ENO) and weighted essentially non-oscillatory (WENO) methods. Higher order in time is achieved through the development of strongly stable Runge-Kutta methods (SSP-RK).

As an alternative but closely related techniques we use the last 1/3 of the class on the development and analysis of discontinuous Galerkin methods as a very general and robust high-order accurate extension of finite volume methods. We shall study these methods in some detail, including their mathematical properties and efficient implementation techniques. We emphasize a general approach that allows for the use of general unstructured grids to solve problems in complex geometries.

Time permitting, we extend the discussion to higher order problems such as advection-diffusion problems and, time permitting, to problems with global constraints such as the incompressible Navier-Stokes equations of fluid dynamics.

Throughout the course there will be an emphasis on mastering mathematical as well as computational aspects of the methods

Lecture notes: The class will be based on the text 

J.S. Hesthaven, 2017, Numerical Methods for Conservation Laws: From Analysis to Algorithms. 

Prerequisites: The class is taught as a Master level class and the expectations are that students are comfortable with analysis, differential equations and with some background in numerical methods for partial differential equations. Programming skills in C/Fortran/Matlab is assumed.

Exercises: There will be weekly exercises throughout the semester. The weekly exercises will mainly be theoretical but may also include computational problems.

Projects: There will be three projects during the class. These will included both theoretical and computational components and will require some programming. The projects have to be handed in and will be discussed as part of the oral examination.

Exams: The exam will be an oral examination with 30 min preparation and 30 min examination. Of this 20 min will focus on the prepared discussion of a given question and 10 min will be used to discuss the three project reports.

Grading: The grading of the oral examination will be based on a weighting of 50% of the three completed projects and 50% on the outcome of the discussion based on the prepared questions.