20.1 Classification and characteristics

The general quasi-linear first order equation in two dimensions takes the form

\begin{displaymath}
\fbox{$\displaystyle
a u_x + b u_y = c
$} %
\end{displaymath} (20.1)

where $a$, $b$, and $c$ may depend on $x$, $y$, and $u$. All scalar first order equation are classified as hyperbolic equations.

The characteristics are defined by

\begin{displaymath}
\fbox{$\displaystyle
\frac{{\rm d}y}{{\rm d}x} = \frac{b}{a}
$} %
\end{displaymath} (20.2)

They will form a single family of lines in the $x,y$-plane. (In contrast, the characteristics of scalar second order hyperbolic differential equations form two intersecting families of lines.)

In general, the variation of a function $u$ of two variables along a line is given by the total differential of calculus,

\begin{displaymath}
\frac{{\rm d}u}{{\rm d}x} = u_x + u_y \frac{{\rm d}y}{{\rm d}x}
\end{displaymath}

but along the characteristic lines that becomes

\begin{displaymath}
\frac{{\rm d}u}{{\rm d}x} = u_x + u_y \frac{b}{a}
\end{displaymath}

and comparing with the partial differential equation, it is seen that the right hand side is $c/a$. So the variation of the solution $u$ along the characteristic lines is given by
\begin{displaymath}
\fbox{$\displaystyle
\frac{{\rm d}u}{{\rm d}x} = \frac{c}{a}
$} %
\end{displaymath} (20.3)

Note that frequently, you may have to solve the ordinary differential equations in a different form or order. For example, if $b/a$ depends on $u$, you will not be able to solve ${\rm d}{y}/{\rm d}{x}$ $\vphantom0\raisebox{1.5pt}{$=$}$ $b/a$ to find $y$ as a function of $x$ since $u$ in $b/a$ is still an unknown function of $x$. But maybe, say, $c/b$ is not a function of $x$, in which case you could solve ${\rm d}{u}/{\rm d}{y}$ $\vphantom0\raisebox{1.5pt}{$=$}$ $c/b$; then you could plug that solution for $u$ as a function of $y$ into ${\rm d}{x}/{\rm d}{y}$ $\vphantom0\raisebox{1.5pt}{$=$}$ $a/b$ to get an equation for $x$ that no longer involves $u$. The bottom line is that it is really best to write the characteristic equations as

\begin{displaymath}
\fbox{$\displaystyle
{\rm d}x : {\rm d}y : {\rm d}u = a : b : c
$} %
\end{displaymath} (20.4)

and pick from those proportionalities the ratio that is easiest to solve first.

In all the unsolved problems in the book, there is at least one solvable ratio. But if there is none, you may be forced to try to change variables, e.g. to polar, or eliminate one variable by, say, differentiating a ratio, hopefully producing a second order ordinary differential equation with one variable eliminated.


Example

Question: (5.30) Solve

\begin{displaymath}
y u_x + x u_y = cu
\end{displaymath}

Solution:

This example wants to solve the partial differential equation

\begin{displaymath}
y u_x + x u_y = cu
\end{displaymath}

For this equation, a ratio like ${\rm d}{u}/{\rm d}{x}$ $\vphantom0\raisebox{1.5pt}{$=$}$ $cu/y$ is not immediately solvable for $u$, since $y$ besides $u$ would be an unknown function of $x$. The only solvable ratio is in fact that between ${\rm d}{x}$ and ${\rm d}{y}$:

\begin{displaymath}
\frac{{\rm d}y}{{\rm d}x} = \frac{x}{y}
\qquad\Rightarro...
...y = x { \rm d}x
\qquad\Rightarrow\qquad
y^2 = x^2 + C_1
\end{displaymath}

where $C_1$ is the integration constant. These characteristic lines are hyperbola; they are sketched in figure 20.1.

Figure 20.1: Characteristics of the partial differential equation of problem 5.30.
\begin{figure}
\begin{center}
\leavevmode
{}
\setlength{\unitlength}{1p...
...(-1,135){\makebox(0,0)[r]{$y$}}
\end{picture}
\end{center}
\end{figure}

Now that $y$ is a known function of $x$, specifically $y$ $\vphantom0\raisebox{1.5pt}{$=$}$ $\sqrt{x^2+C_1}$ assuming it is positive, the ordinary differential equation for $u$ can be solved

\begin{displaymath}
\frac{{\rm d}u}{{\rm d}x} = \frac{cu}{y}
= \frac{cu}{\sqrt{x^2+C_1}}
\end{displaymath}

to give

\begin{displaymath}
\frac{{\rm d}u}{u} = c \frac{{\rm d}x}{\sqrt{x^2+C_1}}
\...
...t = c \ln\left(x + \sqrt{x^2+C_1}\right) + \ln\vert C_2\vert
\end{displaymath}

Taking exponentials and noting that the square root equals $y$, this simplifies to

\begin{displaymath}
u = C_2 (x+y)^c
\qquad\mbox{along a characteristic}\qquad
y^2 = x^2 + C_1
\end{displaymath}

For example, if it is given that $u$ $\vphantom0\raisebox{1.5pt}{$=$}$ 1 at the point $x$ $\vphantom0\raisebox{1.5pt}{$=$}$ $y$ $\vphantom0\raisebox{1.5pt}{$=$}$ 1 shown as a fat dot in figure 20.2, then it follows from the above general expressions that $C_1$ $\vphantom0\raisebox{1.5pt}{$=$}$ 0, so the characteristic line is the line $y$ $\vphantom0\raisebox{1.5pt}{$=$}$ $x$ shown in grey, and that $C_2$ $\vphantom0\raisebox{1.5pt}{$=$}$ $1/2^c$, so you would get $u(x,x)$ $\vphantom0\raisebox{1.5pt}{$=$}$ $x^c$ for $u$ on the grey line.

Figure 20.2: Given the value of $u$ at a single point on a characteristic line, $u$ can be found at every point on that line.
\begin{figure}
\begin{center}
\leavevmode
{}
\setlength{\unitlength}{1p...
...(-1,135){\makebox(0,0)[r]{$y$}}
\end{picture}
\end{center}
\end{figure}



Example

Question: (5.6) Solve the nasty example

\begin{displaymath}
x u_x + y u u_y = -xy.
\end{displaymath}

Solution:

In this case, none of the ratios

\begin{displaymath}
\frac{{\rm d}y}{{\rm d}x}=\frac{uy}{x}
\qquad
\frac{{\...
...rm d}x}=-y
\qquad
\frac{{\rm d}u}{{\rm d}y}=-\frac{x}{u}
\end{displaymath}

is a solvable ordinary differential equation; each involves three variables. You might try to take a derivative of an equation, like

\begin{displaymath}
\frac{{\rm d}^2 u}{{\rm d}x^2} =
- \frac{{\rm d}y}{{\rm ...
...}^2 u}{{\rm d}x^2} =
\frac{u}{x} \frac{{\rm d}u}{{\rm d}x}
\end{displaymath}

which is indeed an ordinary differential equation for $u(x)$ not involving the unknown $y$. But it is an awkward second order nonlinear equation.

The trick is to guess that the combination $xy$ can be found as a function of $u$:

\begin{displaymath}
y \frac{{\rm d}x}{{\rm d}u} + x \frac{{\rm d}y}{{\rm d}u} ...
...\Longrightarrow\qquad
\frac{{\rm d}xy}{{\rm d}u} = - 1 - u
\end{displaymath}

which produces

\begin{displaymath}
xy = - u - {\textstyle\frac12}u^2 + C_1.
\end{displaymath}

This can then be plugged into

\begin{displaymath}
\frac{{\rm d}x}{{\rm d}u} = - \frac{1}{y} = - \frac{x}{xy}
= \frac{x}{u + {\textstyle\frac12}u^2 - C_1}
\end{displaymath}

to get a separable equation giving $x$ as a function of $u$, with another integration constant $C_2$. However, that becomes a mess, involving either an arctan or logarithm, depending on the value of $C_1$.