Difference between revisions of "NumMethodsPDEs"
| Line 10: | Line 10: | ||
:<math>u_x = \frac{\partial u}{\partial x}</math> | :<math>u_x = \frac{\partial u}{\partial x}</math> | ||
| − | =Looking at Nature= | + | =Looking at Nature (& a quick philosophical aside)= |
When we observe nature, certain patterns crop up again and again. For example, in the the field of electrostatics, if a function <math>f</math> describes a distribution of electric charge, then Poisson's equation: | When we observe nature, certain patterns crop up again and again. For example, in the the field of electrostatics, if a function <math>f</math> describes a distribution of electric charge, then Poisson's equation: | ||
| Line 18: | Line 18: | ||
gives the electric potential <math>\varphi</math>. | gives the electric potential <math>\varphi</math>. | ||
| + | However, Poisson's equation also describes the steady state temperature of a material when subjected to some heating. Poisson's equation also crops up when considering gravitational potentials. And there's more... | ||
| + | So what's going on here? Is Poisson's equation fundamental in some way? Well, yes I suppose it is. We can look at it another way, however, and say that Poisson's equation is the way we describe steady state phenomena which involve potentials using the language of maths. OK, this sounds a bit out there, but consider the following very general equation for a moment (it's a second-order linear equation in two variables): | ||
| + | <math>Au_{xx} + Bu_{xy} + Cu_{yy} + Du_x + Eu_y + Fu = G</math> | ||
| − | + | It turns out that we can categorise certain instances of this equation and then relate these categories to the kind of phenomena that they describe: | |
| − | |||
| − | |||
* '''Parabolic''': <math>B^2 - 4AC = 0</math>. This family of equations describe heat flow and diffusion processes. | * '''Parabolic''': <math>B^2 - 4AC = 0</math>. This family of equations describe heat flow and diffusion processes. | ||
| Line 29: | Line 30: | ||
* '''Elliptic''': <math>B^2 - 4AC < 0</math>. Steady-state phenomena. | * '''Elliptic''': <math>B^2 - 4AC < 0</math>. Steady-state phenomena. | ||
| + | That's pretty handy! | ||
| + | |||
| + | =Discritisation= | ||
Taylor Series: | Taylor Series: | ||
<math>f(x + h) = f(x) + \frac{f'(x)}{1!}h + \frac{f''(x)}{2!}h^2 + \cdots + \frac{f^{(n)}(x)}{n!}h^n + R_n(x)</math> | <math>f(x + h) = f(x) + \frac{f'(x)}{1!}h + \frac{f''(x)}{2!}h^2 + \cdots + \frac{f^{(n)}(x)}{n!}h^n + R_n(x)</math> | ||
Revision as of 10:39, 26 January 2011
Numerical Methods for PDEs: Solving PDEs on a computer
Introduction
In this tutorial, we'll take a look at how we might model aspects of the world around us on a computer. When we use mathematics to describe many of the phenomena that we see we end up using Partial Differential Equations (PDEs), and so we need ways to solve these numerically.
As a notation shorthand, we'll use [math]\displaystyle{ u_x }[/math] to represent the partial derivative of u with respect to x:
- [math]\displaystyle{ u_x = \frac{\partial u}{\partial x} }[/math]
Looking at Nature (& a quick philosophical aside)
When we observe nature, certain patterns crop up again and again. For example, in the the field of electrostatics, if a function [math]\displaystyle{ f }[/math] describes a distribution of electric charge, then Poisson's equation:
- [math]\displaystyle{ {\nabla}^2 \varphi = f. }[/math]
gives the electric potential [math]\displaystyle{ \varphi }[/math].
However, Poisson's equation also describes the steady state temperature of a material when subjected to some heating. Poisson's equation also crops up when considering gravitational potentials. And there's more...
So what's going on here? Is Poisson's equation fundamental in some way? Well, yes I suppose it is. We can look at it another way, however, and say that Poisson's equation is the way we describe steady state phenomena which involve potentials using the language of maths. OK, this sounds a bit out there, but consider the following very general equation for a moment (it's a second-order linear equation in two variables):
[math]\displaystyle{ Au_{xx} + Bu_{xy} + Cu_{yy} + Du_x + Eu_y + Fu = G }[/math]
It turns out that we can categorise certain instances of this equation and then relate these categories to the kind of phenomena that they describe:
- Parabolic: [math]\displaystyle{ B^2 - 4AC = 0 }[/math]. This family of equations describe heat flow and diffusion processes.
- Hyperbolic: [math]\displaystyle{ B^2 - 4AC \gt 0 }[/math]. Describe vibrating systems and wave motion.
- Elliptic: [math]\displaystyle{ B^2 - 4AC \lt 0 }[/math]. Steady-state phenomena.
That's pretty handy!
Discritisation
Taylor Series:
[math]\displaystyle{ f(x + h) = f(x) + \frac{f'(x)}{1!}h + \frac{f''(x)}{2!}h^2 + \cdots + \frac{f^{(n)}(x)}{n!}h^n + R_n(x) }[/math]