The Essentials
ID |
MATH761 |
Title |
Numerical ODE/PDE I |
Times |
TuTh, 12:30-1:45, PH301 |
Instructor |
Sorin Mitran |
Office Hours |
Mo 11-12, Tu 9-10, Th 3-4 |
Syllabus
- Basic equations of mathematical physics, other application domains
- Minimal residual formulation
- Numerical methods for ODE's
- Linear advection diffusion equations
- Parabolic equations
- Hyperbolic equations
- Equations of mixed type
- Coupled PDE-ODE systems
Motivation and Objectives
The physical world is described by a small set of laws. One of the most useful mathematical statements of these physical laws is in the form of partial differential equations (PDE's) describing local changes in physical parameters. Though it is usually straightforward to write down the particular PDE's that correspond to a given problem, finding a solution is much more difficult. For the vast majority of real-world problems approximate techniques must be used. One of the most productive approaches is to use a numerical approximation of the PDE's of interest. These techniques are used in studying chemical reactions, ecological models, astrophysical phenomena, aircraft design, financial models and in many other application domains.
Grading Policy
Grading shall be determined based on homework (HW, 48 points), supplementary reading (SR, 8 points), final project (FP, 16 points), midterm examination (ME, 14 points) and final examination (FE, 14 points). The number of points accumulated during classwork is mapped to a graduate grade as shown in the following table.
Clear Excellence |
H |
86-100 points |
Entirely Satisfactory |
P |
71-85 points |
Low Passing |
L |
56-70 points |
Failed |
F |
0-55 points |
Each interval is 15 points; H+/H/H- correspond to 5 point subintervals and so on.
Course Texts and Notes
One of the tenets of graduate study is the ability to critically examine multiple sources in order to verify theories and to enhance understanding. Lecture notes for the course shall be provided online in PDF format. These are meant to document the progress of material presented in class. It is expected and required that the student supplement the lecture notes with reading from other sources.
Lecture |
Topic |
Notes |
Lecture |
Topic |
Notes |
1 (8/19) |
Conservation Laws |
2 (8/21) |
Weighted Residual Method (WRM) |
||
3 (8/26) |
WRM examples - finite difference, finite volume |
4 (8/28) |
ODE systems |
||
5 (9/2) |
ODE methods |
6 (9/4) |
ODE method convergence |
||
7 (9/9) |
ODE LMM methods: stability, consistency |
8 (9/11) |
ODE LMM methods: boundary locus method, A-stable, L-stable |
||
9 (9/16) |
ODE LSM methods: stability regions |
10 (9/18) |
Fourier analysis of PDEs |
||
11 (9/23) |
1D heat equation |
12 (9/25) |
2D heat equation |
||
13 (9/30) |
Advection equation |
14 (10/2) |
Advection equation - numerical methods |
(see lect13.pdf) |
|
15 (10/7) |
Advection equation numerical schemes - Stability analysis |
16 (10/9) |
Modified equations - artificial diffusion, dispersion |
||
17 (10/14) |
Nonlinear hyperbolic equation |
18 (10/16) |
PDE systems |
||
19 (10/21) |
Finite volume methods |
20 (10/16) |
Godunov methods |
||
21 (10/28) |
Second-order finite volume methods |
22 (10/16) |
Limiters |
||
23 (11/4) |
Homework questions |
|
24 (11/6) |
Spectral methods |
|
25 (11/11) |
Spectral methods |
26 (11/6) |
Compact finite differences |
||
27 (11/18) |
Homework questions |
|
28 (11/20) |
|
|
General bibliography
The following texts serve as general background material. Students are encouraged to actively read from these texts material related to the current lecture. Homework and final examination questions from this material should be expected.
Finite Difference Schemes and Partial Differential Equations, John Strikwerda
Numerical Methods for Conservation Laws, Randall LeVeque (Cambridge)
Finite Difference Schemes for Computational Fluid Dynamics, Phillip Colella & Gerry Puckett
Topic bibliography
Texts providing more details on specific aspects we shall discuss shall be listed here as the course progresses.
ODE methods
Finite Difference Methods for Ordinary and Partial Differential Equations, Randall LeVeque (SIAM)
A First Course in the Numerical Analysis of Differential Equations, Arieh Iserles (Cambridge)
Homework and Projects
Posted date |
Due date |
Homework |
Solution |
9/2 |
9/16 |
||
9/23 |
10/6 |
||
10/7 |
10/21 |
Upload new attachment "hw3sol.pdf" Upload new attachment "hw3sol.tm" |
|
10/23 |
11/6 |
||
11/18 |
12/2 |
Useful C and Fortran coding techniques
Defining a function to be an argument to another procedure fEuler.c fEuler.f90
Computer codes
9/19 |
Finite difference operators |
Software
The course homework solutions use the following freely available software:
TeXmacs http://www.texmacs.org
Python with the SciPy packages http://www.python.org, http://www.scipy.org, especially the PyLab Matlab-like environment http://www.scipy.org/PyLab. An excellent Python distribution for the PC and Mac is offered by Enthought http://www.enthought.com/products/epd.php (free for academic use)
Gnuplot http://www.gnuplot.info/
Maxima (symbolic computation package) http://maxima.sourceforge.net/
Examinations
Mid-term examination solution: midterm.pdf
Final examination solution: final.pdf
Final examination is Friday, Dec. 12, noon, Ph301.
Course Bulletin Boards
Please choose the appropriate section for your question. Questions are answered daily, typically late at night.
