Method of characteristics

From Example Problems
(Redirected from Characteristic curve)
Jump to: navigation, search

In mathematics, the method of characteristics is a way of solving partial differential equations, and systems of PDEs.

For a first-order PDE, a characteristic is a line in phase space (comprising the independents, the dependent, and the partial derivatives of the dependent with respect to the independents) along which the PDE degenerates into an ordinary differential equation.

In two dimensions, even a non-linear first-order PDE can always be written in the form

F(x,y,u,u_{x},u_{y})=0,

where x and y are the independents, u(x,y) is the unknown solution, and the remaining arguments are the partial derivatives of the solution u. In the method of characteristics the contour map of F, comprising the level curves where F is constant, is used to recover the solution. This procedure is exact as long as the solution is smooth and differentiable.

While level curves do not intersect, characteristics sometimes do. When this happens the solution function is multivalued, the correct branch must be selected, and where they meet discontinuities arise in the form of shockwaves. Characteristics may also fail to cover part of the domain of the PDE: This is called a rarefaction, and a solution typically exists only in a weak, i.e., integral equation, sense.

The first order wave equation

c{\frac  {\partial u}{\partial x}}+{\frac  {\partial u}{\partial t}}=0\!

describes the movement of a wave in one direction with no change of shape. A solution is shown in Figure 1.3 below as a surface plot and a contour plot.

File:MOC.jpg

Example

Consider the one-dimensional scalar conservation equation

u_{t}+f_{x}(u)=0.

Here u and f are scalar, with u(x,t) a function of x and t.

By the chain rule this equation implies that

u_{t}+f_{u}u_{x}=0

Consider a curve in the x-t plane \eta (t) parameterized as a function of t so that \eta (0)=x_{0} for some x_{0} and \eta _{t}=f_{u}

{\frac  {d(u)}{dt}}\, =u_{t}+\eta _{t}u_{\eta }\,
=u_{t}+f_{u}u_{\eta }\,
=0\,
\int _{0}^{t}{\frac  {d(u)}{dt}}dt=\int _{0}^{t}0\,dt
\int _{0}^{t}\,du=0
u[\eta (t),t]-u[\eta (0),0]=0
u[\eta (t),t]=u(x_{0},0)

So anywhere along the curve \eta (t) (the characteristic curve) the function is the same as at the beginning.

A similar analysis shows that the curve is in fact a straight line for this case.

Bibliography

  • L.C. Evans, Partial Differential Equations, American Mathematical Society, Providence, 1998. ISBN 0-8218-0772-2
  • A. D. Polyanin, V. F. Zaitsev, and A. Moussiaux, Handbook of First Order Partial Differential Equations, Taylor & Francis, London, 2002. ISBN 0-415-27267-X
  • A. D. Polyanin, Handbook of Linear Partial Differential Equations for Engineers and Scientists, Chapman & Hall/CRC Press, Boca Raton, 2002. ISBN 1-58488-299-9