# 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.

## 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 x0 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.