Skip to content

Commit

Permalink
Merge pull request #1 from numerical-mooc/master
Browse files Browse the repository at this point in the history
Update from original
  • Loading branch information
pcgit committed Dec 10, 2014
2 parents 200ad8c + a427490 commit c424f28
Show file tree
Hide file tree
Showing 90 changed files with 149,973 additions and 75 deletions.
2 changes: 1 addition & 1 deletion .gitignore
Original file line number Diff line number Diff line change
@@ -1 +1 @@
lessons/phugoid/.ipynb_checkpoints/*checkpoint.ipynb
.ipynb_checkpoints
36 changes: 35 additions & 1 deletion LICENSE
Original file line number Diff line number Diff line change
@@ -1,3 +1,37 @@
Instructional Material

All instructional material is made available under the Creative
Commons Attribution license. You are free:

* to Share---to copy, distribute and transmit the work
* to Remix---to adapt the work

Under the following conditions:

* Attribution---You must attribute the work using "Copyright (c)
Barbagroup" (but not in any way that suggests that we
endorse you or your use of the work). Where practical, you must
also include a hyperlink to https://github.com/numerical-mooc/numerical-mooc.

With the understanding that:

* Waiver---Any of the above conditions can be waived if you get
permission from the copyright holder.
* Other Rights---In no way are any of the following rights
affected by the license:
* Your fair dealing or fair use rights;
* The author's moral rights;
* Rights other persons may have either in the work itself or in
how the work is used, such as publicity or privacy rights. *
* Notice---For any reuse or distribution, you must make clear to
others the license terms of this work. The best way to do this is
with a link to http://creativecommons.org/licenses/by/3.0/.

For the full legal text of this license, please see:
http://creativecommons.org/licenses/by/3.0/legalcode

Software

The MIT License (MIT)

Copyright (c) 2014 Barba group
Expand All @@ -18,4 +52,4 @@ FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
SOFTWARE.
SOFTWARE.
28 changes: 17 additions & 11 deletions README.md
Original file line number Diff line number Diff line change
@@ -1,29 +1,35 @@
#Practical Numerical Methods in Python
#Practical Numerical Methods with Python

A multi-campus, connected course (plus MOOC) on numerical methods for differential equations in science and engineering. Collaboratively developed by:
- Lorena A. Barba, George Washington University, USA
- [Lorena A. Barba](http://lorenabarba.com), George Washington University, USA
- Ian Hawke, Southampton University, UK
- Carlos Jerez, Pontificia Universidad Catolica de Chile

**Note:** David Ketcheson, from King Abdullah University of Science and Technology (KAUST), Saudi Arabia was going to be our fourth partner, but unfortunately the local course at KAUST got cancelled due to low enrollment.

[**"Practical Numerical Methods with Python"**](http://openedx.seas.gwu.edu/courses/GW/MAE6286/2014_fall/about) is an open, online course hosted on an independent installation of the [Open edX](http://code.edx.org) software platform for MOOCs.
The MOOC (massive open online course) is ran by Prof. Barba at the George Washington University, while each instructor is running a local course, for credit at their institution. You can register for the MOOC at any time in the [GW Online Open edX](http://openedx.seas.gwu.edu/) platform to experience the complete course (including quizzes, examples and discussion board).

All content is open —really open, i.e., you can use, share, mod, remix— and most is available outside the course platform also (on GitHub and YouTube).

####Find the list of IPython Notebooks, with links to nbviewer, in the [Wiki](https://github.com/numerical-mooc/numerical-mooc/wiki).

##List of Modules
##List of Course Modules

1. **The phugoid model of glider flight.**
Described by a set of two nonlinear ordinary differential equations, the phugoid model motivates numerical time integration methods, and we build it up starting from one simple equation, so that the unit can include 3 or 4 lessons on initial value problems. This includes: a) Euler's method, 2nd-order RK, and leapfrog; b) consistency, convergence testing, local vs. global error; c) stability
1. [**The phugoid model of glider flight.**](https://github.com/numerical-mooc/numerical-mooc/tree/master/lessons/01_phugoid)
Described by a set of two nonlinear ordinary differential equations, the phugoid model motivates numerical time integration methods, and we build it up starting from one simple equation, so that the unit can include 3 or 4 lessons on initial value problems. This includes: a) Euler's method, 2nd-order RK, and leapfrog; b) consistency, convergence testing; c) stability
Computational techniques: array operations with NumPy; symbolic computing with SymPy; ODE integrators and libraries; writing and using functions.
2. **Space and Time—Introduction to finite-difference solutions of PDEs.**
2. [**Space and Time—Introduction to finite-difference solutions of PDEs.**](https://github.com/numerical-mooc/numerical-mooc/tree/master/lessons/02_spacetime)
Starting with the simplest model represented by a partial differential equation (PDE)—the linear convection equation in one dimension—, this module builds the foundation of using finite differencing in PDEs. (The module is based on the “CFD Python” collection, steps 1 through 4.) It also motivates CFL condition, numerical diffusion, accuracy of finite-difference approximations via Taylor series, consistency and stability, and the physical idea of conservation laws.
Computational techniques: more array operations with NumPy and symbolic computing with SymPy; getting high performance with Numba.
3. **Riding the wave: convection problems.**
Starting with the inviscid Burgers’ equation in conservation form and a 1D shock wave, cover a sampling of finite-difference convection schemes of various types: upwind, Lax-Friedrichs, Lax-Wendroff, MacCormack, then MUSCL (discussing limiters). Traffic-flow equation with MUSCL (from HyperPython). Reinforce concepts of numerical diffusion and stability, in the context of solutions with shocks. It will motivate spectral analysis of schemes, dispersion errors, Gibbs phenomenon, conservative schemes.
4. **Spreading out: Parabolic PDEs.**
Start with heat equation in 2D (first introduction of two-dimensional FD discretization). Introduce implicit methods: backward Euler, trapezoidal rule (Crank-Nicolson), backward-differentiation formula (BDF). Pattern formation models (reaction-diffusion). Theory content: A-stability (unconditional stability), L-stability (?). Fourier spectral methods and splitting.
Computational techniques: more array operations with NumPy and symbolic computing with SymPy; getting better performance with NumPy array operations.
3. [**Riding the wave: convection problems.**](https://github.com/numerical-mooc/numerical-mooc/tree/master/lessons/03_wave)
Starting with an overview of the concept of conservation laws, this module uses the traffic-flow model to study different solutions methods for problems with shocks: upwind, Lax-Friedrichs, Lax-Wendroff, MacCormack, then MUSCL (discussing limiters). Reinforces concepts of numerical diffusion and stability, in the context of solutions with shocks. It will motivate spectral analysis of schemes, dispersion errors, Gibbs phenomenon, conservative schemes.
4. [**Spreading out: diffusion problems.**](https://github.com/numerical-mooc/numerical-mooc/tree/master/lessons/04_spreadout)
This module deals with solutions to parabolic PDEs, exemplified by the diffusion (heat) equation. Starting with the 1D heat equation, we learn the details of implementing boundary conditions and are introduced to implicit schemes for the first time. Another first in this module is the solution of a two-dimensional problem. The 2D heat equation is solved with both explicit and implict schemes, each time taking special care with boundary conditions. The final lesson builds solutions with a Crank-Nicolson scheme.
5. **Relax and hold steady: elliptic problems.**
Laplace and Poisson equations (steps 9 and 10 of “CFD Python”), explained as systems relaxing under the influence of the boundary conditions and the Laplace operator; introducing the idea of pseudo-time and iterative methods. Linear solvers for PDEs : Jacobi’s method, slow convergence of low-frequency modes (matrix analysis of Jacobi), Jacobi as a smoother, Multigrid.
6. **Boundaries take over: the boundary element method (BEM).**
Weak and boundary integral formulation of elliptic partial differential equations; the free space Green's function. Boundary discretization: basis functions; collocation and Galerkin systems. The BEM stiffness matrix: dense versus sparse; matrix conditioning. Solving the BEM system: singular and near-singular integrals; Gauss quadrature integration.
7. **Tsunami: Shallow-water equation with finite volume method.**
1D first … 2D problem with HPC solution (Python parallel or CUDA Python) -- *optional*.

1 change: 1 addition & 0 deletions lessons/01_phugoid/.gitignore
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
*.pyc
4 changes: 2 additions & 2 deletions lessons/01_phugoid/01_01_Phugoid_Theory.ipynb
Original file line number Diff line number Diff line change
Expand Up @@ -451,7 +451,7 @@
"cell_type": "code",
"collapsed": false,
"input": [
"plot_flight_path(64,16,numpy.pi/2)"
"plot_flight_path(64,16,-numpy.pi/2)"
],
"language": "python",
"metadata": {},
Expand Down Expand Up @@ -711,4 +711,4 @@
"metadata": {}
}
]
}
}
60 changes: 29 additions & 31 deletions lessons/01_phugoid/01_02_Phugoid_Oscillation.ipynb

Large diffs are not rendered by default.

Loading

0 comments on commit c424f28

Please sign in to comment.