Partial Differential Equations
Introduction
While searching for a quantitative description of physical phenomena, the engineer or the physicist establishes generally a system of ordinary or partial differential equations valid in a certain region (or domain) and imposes on this system suitable boundary and initial conditions.
The general form of a second order PDE is :
a(.)Uxx+b(.)Uxy+c(.)Uyy+d(.)Ux+ e(.)Uy+f(.)U+g(.) = 0
where Ux means the partial derivative of U with respect to x, Uxx means the second derivative of U with respect to x2, etc.
a(.), b(.), c(.), ... are functions that can depend on x, y, TIME, etc.
Depending on the coefficients PDEs can be clasified in :
- Hyperbolic : when b2-ac > 0, a typical example is the wave equation.
- Parabolic : when b2-ac = 0, a typical example is the diffusion equation.
- Eliptic : when b2-ac < 0, a typical example is the Poisson equation.
The two firsts types (hyperbolic and parabolic) defines an initial value problem, the third defines a boundary problem.
Boundary conditions can be clasified in three kinds :
- Dirichlet : if we specify the function that we are searching in the boundary, that is : u = g(x,t).
- Neumann : if the derivative of the function is specified : un = g(x,t).
- Robin : if the condition is specified in the form: un+a(x,t)u = g(x,t).
The conditions are called homogeneous when g is zero.
Resolution methods
To solve numerically those equations it is necessary to recast the problem into a purely algebraic form, involving only the basic arithmetic operations. To achieve this, various forms of discretization of the continuum problem can be used. Some of the most popular methods are :
- The finite difference method : replaces each derivative by a discretization, usually truncated Taylor series. There are a lot of schemes, depending on the chosen discretization. The scheme can be explicit - if there's no need to solve a system of equations, just to walk the grid nodes - or implicit if we have to solve a system of equations for each row of the grid.
- The finite elements method : Consists in aproximating the function in small pieces of the domain called finite elements. Once the approximation for each element is determined, the solution for the whole domain is obtained assembling the solution for each element.
- Cellular automatons (CA): Cellular automatons are discrete systems of lattice sites having various initial values. These sites evolve in discrete time steps, each site assumes a new value based on the values of some local neighborhood of sites and a finite number of previous time steps.
There are a lot of methods to discretize the domain, we can clasify them into :
Examples
We have prepared several examples of the solution of PDEs :
Using MGEN to modify the PDE and the conditions
Last modified 2/1/2000 by Juan de Lara ( Juan.Lara@ii.uam.es, http://www.ii.uam.es/~jlara) need help for using this courses?.
The SODA code