Math 273A (Advanced Techniques in Computational Mathematics)

Course Topics: Mathematical and Computational Tools for Partial Differential Equations
Instructor: Prof. Michael Holst (5739 AP&M,; Regular Office Hours: Mon 11-11:50am)
Term: Fall 2018
Lecture: 10:00a-10:50a MWF, 2402 AP&M
TA: None
Discussion: None

Main Class webpage:

Textbook(s): I. Stakgold and M. Holst, Green's Functions and Boundary Value Problems.
                      John Wiley & Sons, New York, NY, Third Edition, 855 pages, 2011.
                      This book is free to UCSD students, and a PDF version of the book can be found [ here ].

Printable Syllabus: Can be found [ here ].

Models of physical systems, calculus of variations, principle of least action. Discretization techniques for variational problems, geometric integrators, advanced techniques in numerical discretization. Project-oriented; projects designed around problems of current interest in science, mathematics, and engineering. Prerequisites: consent of instructor.

GENERAL THEME OF THE COURSE: Math 273ABC was designed as part of an NSF CAREER Award to be a new type of project-oriented course on the mathematical foundations of computational mathematics and physics.

The main idea of the course is that each participant (including the instructor!) picks a topic of particular interest to his/her own research program, involving the analysis and/or numerical solution of an elliptic, parabolic, or hyperbolic partial differential equation. By attending the lectures and also pursuing their own projects, each participant will gain some mastery of the technical tools for analyzing their particular mathematical model using modern PDE methods (is their model well-posed?), for choosing a "good" discretization method (is the method guaranteed to give the right answer to some known precision?), and for designing algorithms for solving the resulting discrete problem on a workstation or possibly a large cluster (does the algorithm have reasonable complexity, and does it scale on a cluster?).

The instructor's particular project this quarter is to analyze and solve a coupled nonlinear elliptic system that arises in mathematical and numerical general relativity; the application is to gravitational wave simulation, and is relevant to NSF's LIGO project. (This problem is the main focus of the instructor's current NSF FRG Project.)

GRADES AND IMPORTANT DATES: There are no required homeworks or exams, although I will mention one or two suggested problems nearly every lecture as suggestions for self-study after the lecture. Your grade in the course is based (33% each) on (1) attending and participating in class, (2) on your work on your own project, and (3) your 20-minute project presentation during the time of the final. NOTE: Your project, and your short presentation at the end of the quarter, can be done individually or in small groups (2-4 people). By the fourth week of the quarter you need to tell me what you picked for your project; I will discuss the projects in class the first week of the quarter, and can answer any questions that you have.

Important dates:

First lecture: FRI 09/28
Last lecture: FRI 12/07
Thanksgiving (no lecture): FRI 11/23
Finals week: MON-FRI, 12/10-12/14
Final Exam (Project Presentations): TBA

SCHEDULE OF LECTURES: We will cover the three major topics listed below. The three major topics will each take about three weeks of the quarter. The individual sub-topics listed within each of the major topics will be covered in varying degree, depending on the available lecture time and the interests of the students. My aim in choosing the topics for the lectures was to provide a general set of mathematical and computational tools for partial differential equations, for mathematicians, scientists, and engineers who want to have a toolkit for doing state of the art computational science. (A secondary aim is to make all of us stronger overall applied mathematicians!) NOTE: My focus will be on tools (mathematical and computational) for elliptic and parabolic PDE; if you choose a project based around hyperbolic PDE, then that is fine, but you will have to do a little self-study as well.
  1. Linear Functional Analysis and Linear PDE:
    Dirichlet, Neumann, and Robin boundary-value problems for linear elliptic partial differential equations (PDE), with applications in mathematical physics. Classical, strong, weak solutions. Green's functions, Distributions, and Sobolev spaces. Vector spaces, norms, inner-products, Banach and Hilbert spaces. Linear operators, adjoints, four fundamental subspaces, operator norms, bounded, unbounded, closed, compact, projection operators, Fredholm alternative. Linear functionals, dual space, reflexive and separable spaces. Ellipticity, inf-sup conditions, energy estimates, Lax-Milgram Theorem, existence, uniqueness. Maximum principles, a priori L-infinity estimates, Harnack inequality. Petrov-Galerkin methods, discrete inf-sup conditions, Babuska-Brezzi-Ladyzenskaya Theorem, a priori error estimates, applications to linear elliptic equations. Finite volume, finite element, and related discretization methods.
  2. Nonlinear Analysis and Nonlinear PDE:
    Dirichlet, Neumann, and Robin boundary-value problems for nonlinear elliptic partial differential equations (PDE), with applications in mathematical physics. Nonlinear maps in Banach spaces, variations, Gateaux and Frechet derivatives, Taylor expansion, Newton's Method. Inverse and Implicit Function Theorems, Bifurcation. Weakly closed sets, weak-lower-semi-continuity, weak convergence. Stationarity of functionals, existence of minimizers. Best approximation in Banach spaces, a priori error estimates for nonlinear problems, applications to nonlinear elliptic equations. Finite volume, finite element, and related discretization methods.
  3. Iterative Methods for Linear and Nonlinear Problems:
    Stationary iterative methods for linear problems, convergence, complexity. Multigrid and multilevel iterative methods for linear problems, convergence, complexity. Domain decomposition methods for linear problems. convergence, complexity. Applications to discretized linear elliptic equations. Fixed-point iteration for nonlinear problems. Banach Fixed-Point Theorem, ordered Banach spaces, maximum principle, monotone increasing maps. Method of sub- and super-solutions, compatible barriers. Newton's method, convergence, inexactness, globalization. Multilevel and domain decomposition methods for nonlinear problems. Basic techniques of bifurcation theory, homotopy methods, and numerical continuation methods. Applications to continuous and discretized nonlinear elliptic equations.