Artificial Intelligence 🤖
Differential equations
Introduction

Introduction

Classification

Order of PDE

The order is defined to be the order of the highest derivative in the equation. If the highest derivative is of order kk, then the equation is said to be of order kk.

A first-order partial differential equation involves two or more independent variables (say, xx and yy ) and a function u(x,y)u(x, y) and its first partial derivatives, uxu_{x} and uyu_{y}. Recall that the subscripts ' xx ' and ' yy ' denote partial differentiation with respect to xx and yy, respectively. For example, Eq. (2.2) is a first-order equation:

uxuy=0\frac{\partial u}{\partial x}-\frac{\partial u}{\partial y}=0

Our objective is to find the most general form of a differentiable function, u(x,y)u(x, y) which satisfies the PDE in some region of interest of the x,yx, y variables. This is the most general form of the solution which normally contains an arbitrary function. Compare this with the general solution of a first-order ordinary differential equation (ODE) which contains an arbitrary constant.

A second-order partial differential equation involves two or more independent variables (say, xx and yy ) and a function u(x,y)u(x, y), its first partial derivatives, uxu_{x} and uyu_{y} and second partial derivatives, uxx,uxyu_{x x}, u_{x y} and uyyu_{y y}. Equation (2.3) is a second-order equation:

2ux22uy2=u(x,y).\frac{\partial^{2} u}{\partial x^{2}}-\frac{\partial^{2} u}{\partial y^{2}}=u(x, y) .

The general solution to Eq. (2.3) contains two arbitrary functions.

Linear & nonlinear PDEs

An equation is called linear if in Eq. (2.1), FF is a linear function of the unknown function uu and its derivatives. For example, the equation for u(x,y)u(x, y),

x7ux+exyuy+sin(x2+y2)u=x3x^{7} u_{x}+e^{x} y u_{y}+\sin \left(x^{2}+y^{2}\right) u=x^{3}

is a linear equation, while the equation,

ux2+u2+y=1u_{x}^{2}+u^{2}+y=1

is a nonlinear equation. Nonlinear equations are often further classified according to the type of the nonlinearity. The following two PDEs are both nonlinear:

uxx+uyy=u3,uxx+uyy=u2u.\begin{array}{r} u_{x x}+u_{y y}=u^{3}, \\ u_{x x}+u_{y y}=|\nabla u|^{2} u . \end{array}

While Eq. (2.7) is nonlinear, it is said to be linear as a function of the highest-order derivative. Such a nonlinearity is called quasilinear. In Eq. (2.6), the nonlinearity is only in the unknown function, uu. Such equations are often called semilinear.

Systems of equations

A single PDE with just one unknown function is called a scalar equation. In contrast, a set of mm equations with ll unknown functions is called a system of mm equations (like, for example, Eqs. (2.16) defined below).

Examples in mathematical modelling

Before we start our discussion of solution techniques for PDEs, we give a short list of examples which occur in physical applications.

Equations representing transport phenomena can take the forms:

ut+ux=0,ut+xux=0.\begin{aligned} u_{t}+u_{x} & =0, \\ u_{t}+x u_{x} & =0 . \end{aligned}

Such equations can be used to model a substance (mass or energy) flowing in a region of space. Both equations in (2.8) are first-order and linear.

The nonlinear transport equation or shock-wave equation is:

ut+uux=0u_{t}+u u_{x}=0

Such equations play a crucial role in the modelling of gas dynamics, traffic flow, chromatography and chemical reactions (see also Subsec. 2.2.4). Equation (2.9) is also known as the Riemann equation or the (inviscid) Burgers' equation.

The one-dimensional wave equation is given by:

utt=c2uxxu_{t t}=c^{2} u_{x x}

where cc is a constant. Equation (2.10) may be used to describe sound waves in a pipe or vibrations of a string. The two-dimensional version is expressed as follows,

utt=c22u.u_{t t}=c^{2} \nabla^{2} u .

The diffusion equation

ut=kuxx,u_{t}=k u_{x x},

is a conservation equation used to model heat or mass transfer. Equation (2.12) is the 1D1 \mathrm{D} version of the diffusion equation and it is second-order and linear. The 2D2 \mathrm{D} version is given by:

ut=k2uu_{t}=k \nabla^{2} u \text {. }

In the absence of temporal effects, i.e when ut=0u_{t}=0, Eq. (2.13) reduces to Laplace's equation,

2u=0\nabla^{2} u=0

whose solutions are described by harmonic functions.

The Shrödinger equation is a linear PDE:

(22m2+V)Ψ=iΨt\left(-\frac{\hbar^{2}}{2 m} \nabla^{2}+V\right) \Psi=i \hbar \Psi_{t}

where Ψ\Psi is the complex wave function.

As a final example, consider the Navier-Stokes equations for a vector u\boldsymbol{u} where u=(u,v,w)\boldsymbol{u}=(u, v, w) are the x,yx-, y- and zz - components of the vector. You will be introduced to these PDEs in Fluid Mechanics II. The Navier-Stokes equations are Newton's second law of motion for a fluid, based on the conservation of momentum. The equations are given by:

ρ(ut+(u)u)=P+μ2u+f,\rho\left(\boldsymbol{u}_{t}+(\boldsymbol{u} \cdot \nabla) \boldsymbol{u}\right)=-\nabla P+\mu \nabla^{2} \boldsymbol{u}+\boldsymbol{f},

where f\boldsymbol{f} models external forces. Equation (2.16) is considered a system of PDEs since u\boldsymbol{u} is a vector and they are used to model the motion of a Newtonian fluid. Note that we have 4 unknown quantities in (2.16): these are the 3 scalar velocities (i.e. u,vu, v and ww ) and the pressure, PP. We therefore need an additional equation: this is given by the continuity equation. In the case of an incompressible fluid, i.e. the density, ρ\rho is constant, the continuity equation takes the form:

u=0\nabla \cdot \boldsymbol{u}=0

Prerequisites

To understand PDEs, we need to know all the basic facts about partial derivatives and multiple integrals. These include use of the chain rule, directional derivatives, change of variables in a double integral, Green's theorem and the divergence theorem for the integrals of derivatives. Further, knowledge of Fourier series expansions is needed for the solution of some linear PDES. It is useful to keep the following in mind when dealing with PDEs:

  • We assume that all derivatives exist and are continuous. Hence, mixed derivatives are equal, e.g. uxy=uyxu_{x y}=u_{y x}.
  • Derivatives are local. To calculate the derivative ux(x,y)u_{x}(x, y) at x=x0x=x_{0}, y=y0y=y_{0} we just need to know the values of u(x,y0)u\left(x, y_{0}\right) for xx near x0x_{0} since the derivative is given by the limit xx0x \rightarrow x_{0}. Finally, in solving certain classes of PDEs, we often reduce them to ODEs. A review of ODEs covered in Year I Mathematics is recommended prior to studying PDEs.