CSE-690

Numerical Methods for Atmospheric Modeling

Quick Info

Comprehensive air quality computer models are widely used investigation tools in environmental research; many physical and chemical processes are modeled and their integrated impacts on atmospheric pollutant concentrations studied. Air quality models are also important tools for regulatory and policy communities. This course will introduce the basic numerical and computational techniques used in the state-of-the-art air quality models. Programming assignments will underline the theory, such that at the end of the class each student will have developed his/her own mini-model. The following topics will be discussed: the continuity equations, treatment of chemical terms, solution of transport equations, aerosol dynamics, non-split methods, sensitivity analysis. We shall use the following texts:

Lecture notes after ``Numerical insights in atmospheric modeling'' by A. Sandu, G. Carmichael and F. Potra (main text);
``Fundamentals of Atmospheric Modeling'' by M.Z. Jacobson;
``Atmospheric Chemistry and Physics: from Air Pollution to Climate Change'' by J.H. Seinfeld and S. Pandis;

Course syllabus (pdf).

Class Notes

!---->

Week Who Everything PostScript

1 Amy Hamlin week1.tar.gz wk1.ps

2 Philip Wells week2.tar.gz wk2.ps
8 Jie Yang week8.tar.gz wk8.ps

9 Sourabh Shikhare wk9.tex wk9.ps

10 Zhang Zhang week10.tar.gz wk10.ps

Home Projects

Number Topic Due date Download

1 Basic Chemical Modeling. April 7 File

2 Advanced Time Integration. May 10 File

Grades (out of 10)

Student Id Project 1

7836 10.0

2253 10.0

0064

2855 10.0

5709

7417 10.0

Useful Links

KPP© - The Kinetic PreProcessor:
- Information: User's Manual,
- Code: kpp.tar.gz
- A chemical system has to be described in a specific language; KPP analyses the chemical mechanism and builds the derivative function and the Jacobian describing the chemical transformations. Kpp offers transparent support for treating sparsity, checks chemical equations for balance, and provides links to different numerical integrators; the output is C or Fortran code (ready to run).
- KPP is being succesfully used in the academic communities as well as in industry in many countries around the world.
Fast chemical solvers plus a comprehensive set of test problems from atmospheric chemistry.
- Information: README file
- Code: via ftp from Ode_benchmark.tar.gz.
LSODE - The Livermore solver for ordinary differential equations
- It solves stiff and nonstiff systems of the form dy/dt = f. In the stiff case, it treats the Jacobian matrix df/dy as either a full or a banded matrix, and as either user-supplied or internally approximated by difference quotients. It uses Adams methods (predictor-corrector) in the nonstiff case, and Backward Differentiation Formula (BDF) methods in the stiff case. The linear systems that arise are solved by direct methods (LU factor/solve). LSODE supersedes the older GEAR and GEARB packages, and reflects a complete redesign of the user interface and internal organization, with some algorithmic improvements.
- Code: Download Index.
RUNGE-KUTTA, ROSENBROCK and other time integration codes from E. Hairer, U. Geneva.
- Codes and Friends.
OMEGA, The Operational Multiscale Environment model with Grid Adaptivity.
- Home Page.
Miscellaneous
- Matlab files useful for project 1: example_f.m, example_jac.m, example_fjac.m, example_exact.m, example_beuler.m.
- Sample Fortran 90 implementation of backward Euler/project 1: project1.tar.gz
asandu@mtu.edu (Adrian Sandu)
http://www.cs.mtu.edu/~asandu/Courses/CSE690/CSE690.html

Week	Who	Everything	PostScript
1	Amy Hamlin	week1.tar.gz	wk1.ps
2	Philip Wells	week2.tar.gz	wk2.ps
8	Jie Yang	week8.tar.gz	wk8.ps
9	Sourabh Shikhare	wk9.tex	wk9.ps
10	Zhang Zhang	week10.tar.gz	wk10.ps

Number	Topic	Due date	Download
1	Basic Chemical Modeling.	April 7	File
2	Advanced Time Integration.	May 10	File

Student Id	Project 1
7836	10.0
2253	10.0
0064
2855	10.0
5709
7417	10.0