INTRODUCTION TO PDEs In this course we are interested in the numerical approximation of PDEs using finite difference methods (FDM). We will use some simple prototype boundary value problems (BVP) and initial boundary value problems (IBVP) to study the methods. We begin by giving a brief overview of PDEs and introducing some terminology. We will concentrate mainly on second order PDEs. Note that a second order PDE can be classified as one of 3 types: elliptic, parabolic, or hyperbolic. 1 BVPs for elliptic equations Suppose we are given an open, bounded, connected domain in R 2 denoted by Ω with smooth boundary Ω. The prototype problem for second order elliptic problems is: 1 Seek a (real-valued) function u = u(x, y) defined on Ω = Ω Ω satisfying u = f(x, y) (1) where w = w xx + w yy with w xx denoting the second partial derivative of w with respect to x. (1) is called Poisson s equation, is called the Laplacian operator and (1) is called Laplace s equation in the special case f = 0. In general, we say the DE is homogeneous if the source term is zero. Note that in R 3 the Laplacian operator is just w = w xx + w yy + w zz For a specific f(x, y) we know that there are many such u s satisfying (1). To uniquely specify a u we must impose auxiliary conditions, in this case boundary values for u or its derivatives, or a combination of both. (We will see that we need to be careful here!) If we specify u(x, y) = g(x, y) (x, y) Ω (2) then we are imposing a Dirichlet boundary condition. If g = 0, we say it is a homogeneous Dirichlet boundary condition; otherwise we call the boundary condition inhomogeneous. Equations (1) and (2) form a Dirichlet BVP for u.
If we specify u = g(x, y) (x, y) Ω (3) where n is the unit outer normal and u = u n. Equations (1) and (3) form a Neumann BVP. Remark: constant. The Neumann BVP is not uniquely determined since if u satisfies Equations (1) and (3) then so does u + C where C is an arbitrary We can also specify a mixed or Robin boundary condition of the type αu(x, y) + β u (x, y) = g(x, y) (4) Here, if α = 0 then we have a Neumann boundary condition and if β = 0 we have a purely Dirichlet boundary condition. Note that we can write Ω = Ω 1 Ω 2 where Ω 1 Ω 2 = 0 and specify one type of boundary condition on Ω 1 and another on Ω 2. 2 Remark: For sufficiently smooth boundary (say piecewise C 1 ), f continuous in Ω and g continuous on Ω, the Dircihlet problem (1) -(2) has a unique solution. For cases where Ω is simple (such as a rectangle, circle, etc), we can use the technique of separation of variables to obtain an infinite series form of the exact solution. Example: Consider Poisson s equation with homogeneous Dirichlet boundary data posed on the unit square: u = f(x, y) = x 2 y in Ω = (0, 1) (0, 1) u = 0 on Ω whose exact solution is given by with γ n,m = u(x, y) = 4 n 2 (m 2 + n 2 ) n,m=1 1 1 0 γ m,n sin(nπx) sin(mπy) 0 f(x, y) sin(nπx) sin(mπy) dxdy When the domain is more general or we have variable coefficients for example, then we must turn to approximating the solution to the BVP. In fact, even though we have the exact solution to the above example, it is in terms of an infinite series and convergence is slow. In addition, we may have to use numerical integration to evaluate the coefficients γ n,m for certain choices of f(x, y).
Physical setting: Elliptic problems model states of equilibrium. For example, consider a membrane whose boundary is attached to a fixed frame. u is the displacement at (x, y, t). For elliptic equations we consider only the equilibrium shape of the membrane. We give a force f independent of t and displacement u at Ω. We seek u satisfying the given displacement on the boundary such that u = f. This equation often arises in electric potential theory, equilibrium flow, irrotational 2-D steady motion of an inviscid incompressible fluid. 2 IVP and IBVP for parabolic and hyperbolic equations Parabolic and hyperbolic equations model the time evolution of a physical phenomena from some know time (such as t = 0). The prototype equation for parabolic type second order equations is the heat (or diffusion) equation. We can consider either a purely initial value problem or a initial boundary value problem. In one space dimension we have the following. Seek u(x, t) satisfying the purely initial value problem 3 u t = νu xx < x <, t > 0 (5) u(x, 0) = u 0 (x) < x < (6) where ν > 0 is a given constant and u 0 is a given function of x. This is called the IVP (or Cauchy problem) for the heat equation in one dimension. We can also consider an IBVP for the heat equation when we have a bounded domain, say [a, b]. Here we choose to impose Dirichlet boundary data. Seek u(x, t) satisfying the IBVP u t = νu xx a < x < b, t > 0 (7) u(x, 0) = u 0 (x) a x b (8) u(a, t) = f(t) u(b, t) = g(t) t 0 (9) Of course, we must have that f(0) = u 0 (a) and g(0) = u 0 (b) for the problem to be well-defined. If the diffusion equation is defined on a domain in R 2 then we have the obvious extension. For example, for the IBVP defined on [a, b] [c, d] with Dirichlet boundary data we have:
Seek u(x, y, t) satisfying the IBVP u t = ν u = ν(u xx + u yy ) a < x < b, c < y < d t > 0 u(x, y, 0) = u 0 (x, y) a < x < b, c < y < d u(a, t) = f(t) u(b, t) = g(t) t > 0 Physical setting: temperature diffusion, diffusion of some concentration A second type of equation which models a time-dependent process is the hyperbolic type. The prototype equation for hyperbolics is the wave equation. An IBVP for the wave equation in one space dimension is given by the following. Seek u(x, t) satisfying the IBVP u tt = c 2 u xx a < x < b, t > 0 (10) u(x, 0) = u 0 (x), u t (x, 0) = v 0 (x) a x b (11) u(a, t) = f(t) u(b, t) = g(t) t 0 (12) 4 Note that in this case we need to impose two initial conditions since the time derivative appears as a second derivative. Of course we could also impose Neumann or mixed boundary conditions for these IBVPs. Physical setting: A physical example of the wave equation is a string where u is the displacement. One might ask how it is possible to determine if a particular second order linear PDE is elliptic, hyperbolic or parabolic. A general second order linear PDE in two variables (not necesarily spatial variables) can be written as au xx + bu xy + cu yy + du x + eu y + fu = g The discriminant, b 2 4ac, determines the classification of the equation as follows if b 2 4ac { = 0 the equation is parabolic < 0 the equation is elliptic > 0 the equation is hyperbolic For example, in Poisson s equation, a = c = 1, b = 0, (d = e = f = 0 also) so b 2 4ac = 4 < 0 so the equation is elliptic. For the heat equation we have a = 1, b = c = 0 so b 2 4ac = 0 and the equation is parabolic.
3 Other classifications All of the equations we have considered thus far are linear equations. A differential equation (DE) is linear if the unknown and its derivatives appearly linearly. Otherwise the equation is nonlinear.. A standard example (in one space dimension) of a second order nonlinear equation is Burger s equation u t + uu x = νu xx (13) Other examples include where g(u) is a nonlinear function of u such as u 3 or sin u. u + g(u) = f(x, y) We can also consider systems of differential equations. For example, equations modeling fluid flow are well known. The steady state Navier-Stokes equations are given by 5 ν u + u grad u + grad p = f div u = 0 plus appropriate boundary conditions. Of particular interest to us are first order (symmetric) hyperbolic systems. For example, if u = (u 1, u 2,... u n ) satisfies where A is an n n matrix. u t = A u x Remark: In the sequel we shall assume that the mathematical problem we want to approximate has a unique solution with sufficient continuity and differentiability properties (i.e, sufficient regularity ) and that the solution depends continuously on the data. Remark: We will see that to show a particular numerical method results in an approximation which converges to the exact solution at an optimal rate, we will have to require stronger continuity and differentiability conditions on the solution than was required for the proof of existence and uniqueness. This is not unique to finite difference methods (FDM).