First-order partial differential equation

In mathematics, a first-order partial differential equation is a partial differential equation that involves only first derivatives of the unknown function of n variables. The equation takes the form

${\displaystyle F(x_{1},\ldots ,x_{n},u,u_{x_{1}},\ldots u_{x_{n}})=0.\,}$

Such equations arise in the construction of characteristic surfaces for hyperbolic partial differential equations, in the calculus of variations, in some geometrical problems, and in simple models for gas dynamics whose solution involves the method of characteristics. If a family of solutions of a single first-order partial differential equation can be found, then additional solutions may be obtained by forming envelopes of solutions in that family. In a related procedure, general solutions may be obtained by integrating families of ordinary differential equations.

General solution and complete integral

The general solution to the first order partial differential equation is a solution which contains an arbitrary function. But, the solution to the first order partial differential equations with as many arbitrary constants as the number of independent variables is called the complete integral. The following n-parameter family of solutions

${\displaystyle u=\phi (x_{1},x_{2},\dots ,x_{n},a_{1},a_{2},\dots ,a_{n})}$

is a complete integral if ${\displaystyle {\text{det}}|\phi _{x_{i}a_{j}}|\neq 0}$.[1]

Characteristic surfaces for the wave equation

Characteristic surfaces for the wave equation are level surfaces for solutions of the equation

${\displaystyle u_{t}^{2}=c^{2}\left(u_{x}^{2}+u_{y}^{2}+u_{z}^{2}\right).\,}$

There is little loss of generality if we set ${\displaystyle u_{t}=1}$: in that case u satisfies

${\displaystyle u_{x}^{2}+u_{y}^{2}+u_{z}^{2}={\frac {1}{c^{2}}}.\,}$

In vector notation, let

${\displaystyle {\vec {x}}=(x,y,z)\quad {\hbox{and}}\quad {\vec {p}}=(u_{x},u_{y},u_{z}).\,}$

A family of solutions with planes as level surfaces is given by

${\displaystyle u({\vec {x}})={\vec {p}}\cdot ({\vec {x}}-{\vec {x_{0}}}),\,}$

where

${\displaystyle |{\vec {p}}\,|={\frac {1}{c}},\quad {\text{and}}\quad {\vec {x_{0}}}\quad {\text{is arbitrary}}.\,}$

If x and x₀ are held fixed, the envelope of these solutions is obtained by finding a point on the sphere of radius 1/c where the value of u is stationary. This is true if ${\displaystyle {\vec {p}}}$ is parallel to ${\displaystyle {\vec {x}}-{\vec {x_{0}}}}$. Hence the envelope has equation

${\displaystyle u({\vec {x}})=\pm {\frac {1}{c}}|{\vec {x}}-{\vec {x_{0}}}\,|.}$

These solutions correspond to spheres whose radius grows or shrinks with velocity c. These are light cones in space-time.

The initial value problem for this equation consists in specifying a level surface S where u=0 for t=0. The solution is obtained by taking the envelope of all the spheres with centers on S, whose radii grow with velocity c. This envelope is obtained by requiring that

${\displaystyle {\frac {1}{c}}|{\vec {x}}-{\vec {x_{0}}}\,|\quad {\hbox{is stationary for}}\quad {\vec {x_{0}}}\in S.\,}$

This condition will be satisfied if ${\displaystyle |{\vec {x}}-{\vec {x_{0}}}\,|}$ is normal to S. Thus the envelope corresponds to motion with velocity c along each normal to S. This is the Huygens' construction of wave fronts: each point on S emits a spherical wave at time t=0, and the wave front at a later time t is the envelope of these spherical waves. The normals to S are the light rays.

Two-dimensional theory

The notation is relatively simple in two space dimensions, but the main ideas generalize to higher dimensions. A general first-order partial differential equation has the form

${\displaystyle F(x,y,u,p,q)=0,\,}$

where

${\displaystyle p=u_{x},\quad q=u_{y}.\,}$

A complete integral of this equation is a solution φ(x,y,u) that depends upon two parameters a and b. (There are n parameters required in the n-dimensional case.) An envelope of such solutions is obtained by choosing an arbitrary function w, setting b=w(a), and determining A(x,y,u) by requiring that the total derivative

${\displaystyle {\frac {d\varphi }{da}}=\varphi _{a}(x,y,u,A,w(A))+w'(A)\varphi _{b}(x,y,u,A,w(A))=0.\,}$

In that case, a solution ${\displaystyle u_{w}}$ is also given by

${\displaystyle u_{w}=\phi (x,y,u,A,w(A))\,}$

Each choice of the function w leads to a solution of the PDE. A similar process led to the construction of the light cone as a characteristic surface for the wave equation.

If a complete integral is not available, solutions may still be obtained by solving a system of ordinary equations. To obtain this system, first note that the PDE determines a cone (analogous to the light cone) at each point: if the PDE is linear in the derivatives of u (it is quasi-linear), then the cone degenerates into a line. In the general case, the pairs (p,q) that satisfy the equation determine a family of planes at a given point:

${\displaystyle u-u_{0}=p(x-x_{0})+q(y-y_{0}),\,}$

where

${\displaystyle F(x_{0},y_{0},u_{0},p,q)=0.\,}$

The envelope of these planes is a cone, or a line if the PDE is quasi-linear. The condition for an envelope is

${\displaystyle F_{p}\,dp+F_{q}\,dq=0,\,}$

where F is evaluated at ${\displaystyle (x_{0},y_{0},u_{0},p,q)}$, and dp and dq are increments of p and q that satisfy F=0. Hence the generator of the cone is a line with direction

${\displaystyle dx:dy:du=F_{p}:F_{q}:(pF_{p}+qF_{q}).\,}$

This direction corresponds to the light rays for the wave equation. To integrate differential equations along these directions, we require increments for p and q along the ray. This can be obtained by differentiating the PDE:

${\displaystyle F_{x}+F_{u}p+F_{p}p_{x}+F_{q}p_{y}=0,\,}$

${\displaystyle F_{y}+F_{u}q+F_{p}q_{x}+F_{q}q_{y}=0,\,}$

Therefore the ray direction in ${\displaystyle (x,y,u,p,q)}$ space is

${\displaystyle dx:dy:du:dp:dq=F_{p}:F_{q}:(pF_{p}+qF_{q}):(-F_{x}-F_{u}p):(-F_{y}-F_{u}q).\,}$

The integration of these equations leads to a ray conoid at each point ${\displaystyle (x_{0},y_{0},u_{0})}$. General solutions of the PDE can then be obtained from envelopes of such conoids.

Definitions of linear dependence for differential systems

This part can be referred to ${\displaystyle \S 1.2.3}$ of Courant's book.[2]

We assume that these ${\displaystyle h}$ equations are independent, i.e., that none of them can be deduced from the other by differentiation and elimination.

— Courant, R. & Hilbert, D. (1962), Methods of Mathematical Physics: Partial Differential Equations, II, p.15-18

An equivalent description is given. Two definitions of linear dependence are given for first-order linear partial differential equations.

${\displaystyle (*)\left\{{\begin{array}{*{20}{c}}\sum \limits _{ij}^{}{a_{ij}^{(1)}{\dfrac {\partial {y_{j}}}{\partial {x_{i}}}}}+{f_{1}}=0\\\vdots \\\sum \limits _{ij}^{}{a_{ij}^{(n)}{\dfrac {\partial {y_{j}}}{\partial {x_{i}}}}}+{f_{n}}=0\end{array}}\right.}$

Where ${\displaystyle x_{i}}$ are independent variables; ${\displaystyle y_{j}}$ are dependent unknowns; ${\displaystyle a_{ij}^{(k)}}$ are linear coefficients; and ${\displaystyle f_{k}}$ are non-homogeneous items. Let ${\displaystyle {Z_{k}}\equiv \sum _{ij}^{}{a_{ij}^{(k)}{\frac {\partial {y_{j}}}{\partial {x_{i}}}}}+{f_{k}}}$.

Definition I: Given a number field ${\displaystyle P}$, when there are coefficients (${\displaystyle c_{k}\in P}$), not all zero, such that ${\displaystyle \sum _{k}{{c_{k}}{Z_{k}}=0}}$; the Eqs.(*) are linear dependent.

Definition II (differential linear dependence): Given a number field ${\displaystyle P}$, when there are coefficients (${\displaystyle {c_{k}},d_{kl}\in P}$), not all zero, such that ${\displaystyle \sum _{k}{{c_{k}}{Z_{k}}}+\sum _{kl}{{d_{kl}}{\frac {\partial }{\partial {x_{l}}}}{Z_{k}}=0}}$, the Eqs.(*) are thought as differential linear dependent. If ${\displaystyle {d_{kl}}\equiv 0}$, this definition degenerates into the definition I.

The div-curl systems, Maxwell's equations, Einstein's equations (with four harmonic coordinates) and Yang-Mills equations (with gauge conditions) are well-determined in definition II, whereas are over-determined in definition I.

References

1. Garabedian, P. R. (1964). Partial Differential Equations. New York: Wiley. OCLC 527754.
2. Courant, R. & Hilbert, D. (1962). Methods of Mathematical Physics: Partial Differential Equations. II. New York: Wiley-Interscience.

Further reading

• Evans, L. C. (1998). Partial Differential Equations. Providence: American Mathematical Society. ISBN 0-8218-0772-2.
• Polyanin, A. D.; Zaitsev, V. F.; Moussiaux, A. (2002). Handbook of First Order Partial Differential Equations. London: Taylor & Francis. ISBN 0-415-27267-X.
• Polyanin, A. D. (2002). Handbook of Linear Partial Differential Equations for Engineers and Scientists. Boca Raton: Chapman & Hall/CRC Press. ISBN 1-58488-299-9.
• Sarra, Scott (2003). "The Method of Characteristics with applications to Conservation Laws". Journal of Online Mathematics and its Applications.