In mathematics, and in particular analysis, a **partial differential equation** (**PDE**) is an equation involving partial derivatives of an unknown function. The idea is to describe a function indirectly by a relation between itself and its partial derivatives, rather than writing down a function explicitly. A solution of the equation is any function satisfying this relation.

A PDE usually has several solutions; a problem often includes additional boundary conditions which constrain the solution set. Where ordinary differential equations have solutions that are families with each solution characterized by the values of some parameters, for a PDE it is more helpful to think that the parameters are function data (informally put, this means that the set of solutions is much larger). That is true fairly generally, unless the equations are heavily over-determined.

Partial differential equations are ubiquitous in science, as they describe phenomena such as fluid flow, gravitational fields, and electromagnetic fields. They are important in fields such as aircraft simulation, computer graphics, and weather prediction. The central equations of general relativity and quantum mechanics are also partial differential equations.

## Notation and examples

In PDEs, it is common to write the unknown function as *u* and its partial derivative with respect to the variable *x* as *u*_{x}, that is:

### Laplace's equation

A very important and basic PDE is Laplace's equation:-

*u*_{xx} + *u*_{yy} + *u*_{zz} = 0

for the unknown function *u*(*x*,*y*,*z*). Solutions to this equation, known as harmonic functions, serve as the potentials of vector fields in physics, such as the gravitational or electrostatic fields.

A generalization of Laplace's equation is Poisson's equation:-

*u*_{xx} + *u*_{yy} + *u*_{zz} = *f*

where *f*(*x*,*y*,*z*) is a given function. The solutions to this equation describe potentials of gravitational and electrostatic fields in the presence of masses or electrical charges, respectively.

### Wave equation

The wave equation is an equation for an unknown function *u*(*x*,*y*,*z*,*t*) (where we think of *t* as a time variable) which reads:-

*u*_{tt} = *c*^{2}(*u*_{xx} + *u*_{yy} + *u*_{zz})

Its solutions describe waves such as sound or light waves; *c* is a number which represents the speed of the wave. In lower dimensions, this equation describes the vibration of a string or drum. Solutions will typically be combinations of oscillating sine waves.

### Heat equation

The heat equation describes the temperature in a given region over time. It is:-

*u*_{t} = *k*(*u*_{xx} + *u*_{yy} + *u*_{zz})

Solutions will typically "even out" over time. The number *k* describes the thermal conductivity of the material.

### Euler-Tricomi equation

The Euler-Tricomi equation is used in the investigation of transonic flow. It is

*u*_{xx} = *x**u*_{yy}

### Advection equation

The advection equation describes the transport of a conserved scalar ψ in a velocity field . It is:

- ψ
_{t} + (*u*ψ)_{x} + (*v*ψ)_{y} + (*w*ψ)_{z} = 0.

If the velocity field is solenoidal (that is, ), then the equation may be simplified to

- ψ
_{t} + ψ.*u*_{x} + ψ.*u*_{y} + ψ.*w*_{z} = 0.

The one dimensional steady flow advection equation ψ_{t} + *u*.ψ_{x} = 0 (where *u* is constant) is commonly referred to as the pigpen problem . If *u* is not constant and equal to ψ the equation is referred to as Burgers' equation

### Other Examples

The Schrödinger equation is a PDE at the heart of non-relativistic quantum mechanics. In the WKB approximation it is the Hamilton-Jacobi equation.

Except for Burgers' equation, all the above equations are **linear** in the sense that they can be written in the form *Au* = *f* for a given linear operator *A* and a given function *f*. Other important non-linear equations include the Navier-Stokes equations describing the flow of fluids, and Einstein's field equations of general relativity.

## Methods to solve PDEs

Linear PDEs are generally solved, when possible, by decomposing the equation according to a set of basis functions, solving those individually and using superposition to find the solution corresponding to the boundary conditions. The method of separation of variables has many important particular applications.

There are no generally applicable methods to solve non-linear PDEs. Still, existence and uniqueness results (such as the Cauchy-Kovalevskaya theorem) are often possible, as are proofs of important qualitative and quantitative properties of solutions (getting these results is a major part of analysis).

Nevertheless, some techniques can be used for several types of equations. The h-principle is the most powerful method to solve underdetermined equations. The Riquier-Janet theory is an effective method for obtaining information about many analytic overdetermined systems.

The method of characteristics can be used in some very special cases to solve partial differential equations.

In some cases, a PDE can be solved via perturbation analysis in which the solution is considered to be a correction to an equation with a known solution. Alternatives are numerical analysis techniques from simple finite difference schemes to the more mature multigrid and finite element methods. Many interesting problems in science and engineering are solved in this way using computers, sometimes high performance supercomputers. However, most problems in science and engineering are tackled using scientific computing rather than numerical analysis, as usually it is not known whether the numerical methods used produce solutions close to the true ones.

## Classification

A single partial differential equation, or even a system of partial differential equations, can usually be classified as parabolic , hyperbolic or elliptic. Given a system of partial differential equations, there is an associated matrix. Properties of the eigensystem of the associated matrix determine the classification of the system.

### Examples

The matrix associated with the system

*u*_{t} + 2*u*_{x} = 0

*v*_{t} - *v*_{x} = 0

has coefficients,

The eigenvectors are (0,1) and (1,0) with eigenvalues 2 and -1. Thus, the system is hyperbolic.

## External links

Last updated: 10-21-2005 15:30:57