Primer on Partial Differential Equations
Table of Contents
A partial differential equation (PDE) is an equation that involves partial derivatives of a multivariable function. Typically, they describe how a quantity (like temperature, displacement, pressure) changes over space and time.
If ordinary differential equations (ODEs) describe evolution over one variable (often time), PDEs handle evolution over multiple variables.
The basic form of a PDE is:
F(x1,x2,…,xn,u,∂x1∂u,∂x2∂u,…,∂xi2∂2u,…)=0
where:
- u is the unknown function of the variables x1,x2,…,xn.
- F is a given function that relates the variables and their derivatives.
- The partial derivatives of u with respect to the variables are denoted as ∂xi∂u, and higher-order derivatives can be represented similarly.
Some examples of PDEs include:
The Heat Equation: Describes how heat diffuses through a medium.
∂t∂u=α∇2u
where α is the thermal diffusivity and ∇2 is the Laplacian operator.
The Wave Equation: Describes how waves propagate through a medium.
∂t2∂2u=c2∇2u
where c is the wave speed.
The Laplace Equation: A special case of the heat equation where there is no time dependence.
∇2u=0
This equation describes steady-state solutions.
The Poisson Equation: A generalization of the Laplace equation that includes a source term.
∇2u=f(x,y,z)
where f(x,y,z) is a known function.
The Navier-Stokes Equations: Describe the motion of fluid substances.
∂t∂u+(u⋅∇)u=−ρ1∇p+ν∇2u
where u is the velocity field, ρ is the fluid density, p is the pressure, and ν is the kinematic viscosity.
Methods of solving PDEs
There are several methods to solve PDEs, including:
- Separation of Variables: This method involves assuming that the solution can be written as a product of functions, each depending on a single variable. This reduces the PDE to a set of ODEs.
- Method of Characteristics: This method is used for first-order PDEs and involves transforming the PDE into a set of ODEs along certain curves called characteristics.
- Fourier Series and Transforms: These methods are used to solve linear PDEs by transforming them into algebraic equations in the frequency domain.
- Finite Difference and Finite Element Methods: These are numerical methods used to approximate solutions to PDEs by discretizing the domain and solving the resulting algebraic equations.
Example: Solving the Heat Equation with Separation of Variables
Consider the one-dimensional heat equation:
∂t∂u=α∂x2∂2u
with boundary conditions:
- u(0,t)=0
- u(L,t)=0
- u(x,0)=f(x)
where L is the length of the rod and f(x) is the initial temperature distribution.
To solve this using separation of variables, we assume a solution of the form:
u(x,t)=X(x)T(t)
Substituting this into the heat equation gives:
αT1dtdT=X1dx2d2X=−λ
where λ is a separation constant.
This leads to two ODEs:
- For T(t):
dtdT+αλT=0
- For X(x):
dx2d2X+λX=0
The boundary conditions lead to a Sturm-Liouville problem for X(x), which has solutions of the form:
Xn(x)=sin(Lnπx)
with eigenvalues λn=(Lnπ)2.
The solution for T(t) is:
Tn(t)=e−αλnt
The general solution is a sum of these modes:
u(x,t)=n=1∑∞Ane−αλntsin(Lnπx)
where An are coefficients determined by the initial condition u(x,0)=f(x) using Fourier series expansion.
The coefficients can be found using:
An=L2∫0Lf(x)sin(Lnπx)dx
This gives the complete solution to the heat equation with the specified boundary and initial conditions.