Subsections


18.8 Two-Dimensional Coordinate Transforms

More powerful simplifications by changing coordinates are possible in 2D.

Assume that in terms of coordinates $x$ and $y$, we have a partial differential equation:

\begin{displaymath}
a u_{xx} + 2 b u_{xy} + c u_{yy} = d
\end{displaymath}

Then, if we transform to new coordinates, call them $\xi$ and $\eta$, we will get a new partial differential equation of the form: $\xi,\eta$

\begin{displaymath}
a' u_{\xi\xi} + 2 b' u_{\xi\eta} + c' u_{\eta\eta} = d'
\end{displaymath}

The idea is again to choose the new coordinates $\xi$ and $\eta$ so that the new partial differential equation is as simple as possible.

For example, for a hyperbolic equation, you may like coordinates $\xi$ and $\eta$ such that $a'$ and $c'$ are zero. To find out for what coordinates $\xi$ and $\eta$ that is the case, expressions for the new coefficients $a'$, $b'$, $c'$ and $d'$ in terms of the new coordinates are needed. These can be found by writing out the general transformation formulae from section 18.7.2 for the special case of two dimensions. You get, {D.5}:

\begin{displaymath}
\fbox{$\displaystyle
\begin{array}{l}
a' =
a \left(\...
... b \eta_{xy} + c \eta_{yy}\right) u_\eta
\end{array}
$} %
\end{displaymath} (18.22)


18.8.1 Characteristic Coordinates

Characteristic coordinates are coordinates so that the the $u_{\xi\xi}$ and $u_{\eta\eta}$ derivatives are eliminated. That leaves only the $u_{\xi\eta}$ derivative, greatly simplifying the partial differential equation. It reduces to the two-di­men­sion­al canonical form:

\begin{displaymath}
\fbox{$\displaystyle
2 b' u_{\xi\eta} = d'
$}
\end{displaymath} (18.23)

The first thing is to find out how this may be achieved. In terms of the coefficients of the transformed equation as discussed above, $a'$ and $c'$ must vanish. The condition $a'$ $\vphantom0\raisebox{1.5pt}{$=$}$ 0 requires, according to the given formulae:

\begin{displaymath}
a \left(\xi_x\right)^2 +
2b \left(\xi_x\right)\left(\xi_y\right) +
c \left(\xi_y\right)^2 = 0
\end{displaymath}

That can be considered to be a partial differential equation for $\xi$. A nonlinear first order equation, to be sure. Similarly for $c'$ to vanish,

\begin{displaymath}
a \left(\eta_x\right)^2 +
2b \left(\eta_x\right)\left(\eta_y\right) +
c \left(\eta_y\right)^2 = 0
\end{displaymath}

Note that $\xi$ and $\eta$ must satisfy the exact same equation, but they must be different solutions. Otherwise they are not valid independent coordinates.

To solve the equation for $\xi$ ($\eta$ goes the same way), divide by $\left(\xi_y\right)^2$:

\begin{displaymath}
a \left(-\frac{\xi_x}{\xi_y}\right)^2 -
2b \left(-\frac{\xi_x}{\xi_y}\right) +
c = 0
\end{displaymath}

and note that, from your calculus or thermo,

\begin{displaymath}
-\frac{\xi_x}{\xi_y} =
\left(\frac{{\rm d}y}{{\rm d}x}\right)_{\hbox{$\xi$ is constant}}
\end{displaymath}

So the lines of constant $\xi$ should satisfy the ordinary differential equation

\begin{displaymath}
\fbox{$\displaystyle
\frac{{\rm d}y}{{\rm d}x} = \frac{b \pm \sqrt{b^2-ac}}{a}
$} %
\end{displaymath}

We can achieve this by taking $\xi$ to be the integration constant in the solution of this ordinary differential equation! Integration constants are, like the word says, constant for solutions.

By taking the other sign for the square root, you can get a second independent coordinate $\eta$.

Bottom line, to get characteristic coordinates, solve the plus and minus sign ordinary differential equations above, and equate the integration constants to $\xi$ and $\eta$.

A couple of notes:

  1. Since integration constants are not unique, the characteristic coordinates are not. But the lines of constant $\xi$ and $\eta$ are unique, and are called characteristic lines or characteristics.
  2. Elliptic equations do not have characteristics, because the square root in the ordinary differential equation would be imaginary. The coordinates $\xi$ and $\eta$ must be real; you do not want to deal with partial differential equations in complex coordinates.
  3. Parabolic equations have only one family of characteristic lines. That is because the square root is zero, so taking the other root does not make a difference.


Example

Question: Use characteristic coordinates to reduce the wave equation in multi-di­men­sion­al canonical form

\begin{displaymath}
u_{tt} - a^2 u_{xx} = 0
\end{displaymath}

to its equivalent two-di­men­sion­al canonical form. Then solve it.

Solution:

First find the characteristics by solving the ordinary differential equation given above:

\begin{displaymath}
\frac{{\rm d}x}{{\rm d}t} = \frac{b \pm \sqrt{b^2-ac}}{a} = \pm a
\end{displaymath}

Note that the final $a$ is the wave propagation speed, not the coefficient $a$ in the generic second order equation.

The solution is simple:

\begin{displaymath}
x = at + \xi \qquad x = -at + \eta
\end{displaymath}

where $\xi$ and $\eta$ are the integration constants (as well as the characteristic coordinates). So the lines $x-at$ $\vphantom0\raisebox{1.5pt}{$=$}$ $\mbox{constant}$ are one set of characteristic lines, and the lines $x+at$ $\vphantom0\raisebox{1.5pt}{$=$}$ $\mbox{constant}$ are the other set.

Now find the coefficient $d'$. The coefficient $d$ was zero, and the second order derivatives of $\xi$ and $\eta$ in the formula for $d'$ are also zero, so $d'$ is zero too. So the wave equation in characteristic coordinates is

\begin{displaymath}
\fbox{$\displaystyle
u_{\xi\eta} = 0
$}
\end{displaymath} (18.24)

Note that $b'$ could be divided out, so there is no need to figure out what it is.

The wave equation can now easily be solved. Integration with respect to $\eta$ gives

\begin{displaymath}
u_\xi = f(\xi)
\end{displaymath}

where the integration constant $f$ can be any arbitrary function of $\xi$. Integrating with respect to $\xi$ gives the final solution:

\begin{displaymath}
u = f_1(\xi)+f_2(\eta)
\end{displaymath}

Here $f_1$ is an antiderivative of $f$, so it is arbitrary just like $f$. The additional integration constant $f_2$ is an arbitrary function of $\eta$.

However, you would surely want the solution in terms of the physical coordinates $x$ and $t$, rather than the mathematical characteristic coordinates. So substitute for $\xi$ and $\eta$ using the obtained equations for the characteristics. That gives the final solution:

\begin{displaymath}
\fbox{$\displaystyle
u = f_1(x-at)+f_2(x+at)
$}
\end{displaymath} (18.25)

That is the general solution of the wave equation. In order to solve a particular problem, you will still need to figure out what $f_1$ and $f_2$ are using whatever the initial and boundary conditions are. One special case, in which the $x$-range is doubly infinite, will be solved in detail later.



Example

Question: Find and sketch the characteristics of the equation

\begin{displaymath}
u_{xx} + y u_{yy} = 0
\end{displaymath}

Solution:

Figure out the coefficients in the the characteristic equation by looking at the partial differential equation:

\begin{displaymath}
\frac{{\rm d}y}{{\rm d}x} = \frac{b\pm\sqrt{b^2 - ac}}{a}
= \pm \sqrt{-y}
\end{displaymath}

Note that there are only characteristics for negative $y$. For positive $y$ the equation is elliptic. And for zero $y$ there will only be one direction for the characteristics, horizontal.

Use separation of variables to solve. In other words, take the $y$ factors to one side and the $x$-factors to the other side:

\begin{displaymath}
\frac{{\rm d}- y}{\sqrt{-y}} = \pm {\rm d}x \quad\quad\Rightarrow\quad\quad 2 \sqrt{-y} = \pm (x - C)
\end{displaymath}

Squaring both sides to get rid of the square root gives

\begin{displaymath}
y = -{\textstyle\frac14}(x - C)^2
\end{displaymath}

These are parabolae.

\begin{displaymath}
\hbox{\epsffile{chexamp1.eps}}
\end{displaymath}



Example

Question: Reduce the equation

\begin{displaymath}
e^y u_{xx} + 2 e^x u_{xy} - e^{2x-y} u_{yy} = 0
\end{displaymath}

to two-di­men­sion­al canonical form.

Solution:

Two-di­men­sion­al canonical form means characteristic form. Find the ordinary differential equation for the characteristics:

\begin{displaymath}
\frac{{\rm d}y}{{\rm d}x} = \frac{b\pm\sqrt{b^2 - ac}}{a}
= (1\pm\sqrt{2}) e^{x-y}
\end{displaymath}

Solve it using separation of variables:

\begin{displaymath}
e^y {\rm d}y = (1 \pm \sqrt2) e^x {\rm d}x \quad\quad\Rightarrow\quad\quad
e^y = (1 \pm \sqrt2) e^x + C
\end{displaymath}

The integration constants are the new coordinates:

\begin{displaymath}
\xi = (1 + \sqrt2) e^x - e^y \qquad
\eta = (1 - \sqrt2) e^x - e^y
\end{displaymath}

Work out the partial differential equation in these coordinates using the formulae given at the start of this section:

\begin{displaymath}
b ' =
a \left(\xi_x\right)\left(\eta_x\right) +
b \lef...
...) +
c \left(\xi_y\right)\left(\eta_y\right) = - 4 e^{2x+y}
\end{displaymath}


\begin{displaymath}
d' = d -
\left(a \xi_{xx} + 2 b \xi_{xy} + c \xi_{yy}\ri...
...left(a \eta_{xx} + 2 b \eta_{xy} + c \eta_{yy}\right) u_\eta
\end{displaymath}

so

\begin{displaymath}
d' = - \left[(1+\sqrt2)e^{x+y}+ e^{2x}\right] u_\xi
- \left[(1-\sqrt2)e^{x+y}+ e^{2x}\right] u_\eta
\end{displaymath}

The partial differential equation becomes

\begin{displaymath}
8 e^{x+y} u_{\xi\eta} =
\left[(1+\sqrt2)e^y+ e^x\right] u_\xi
\left[(1-\sqrt2)e^y+ e^x\right] u_\eta
\end{displaymath}

Get rid of $x$ and $y$ completely using the equations for the characteristics:

\begin{displaymath}
e^x = \frac{1}{2\sqrt2} (\xi-\eta) \quad
e^y = - \frac{1-\sqrt2}{2\sqrt2} \xi- \frac{1+\sqrt2}{2\sqrt2}\eta
\end{displaymath}

The resulting partial differential equation is

\begin{displaymath}
(\xi -\eta)[(1-\sqrt2)\xi + (1+\sqrt2)\eta]u_{\xi\eta}
= (1+\sqrt2)\eta u_\xi + (1-\sqrt2)\xi u_\eta
\end{displaymath}

It does not look easily solvable.



Example

Question: Find the characteristic coordinates of the equation

\begin{displaymath}
\sin^2 (x) u_{xx} + 2 \cos(x) u_{xy} - u_{yy} = 0
\end{displaymath}

Solution:

Find the ordinary differential equation for the characteristics:

\begin{displaymath}
\frac{{\rm d}y}{{\rm d}x} = \frac{b\pm\sqrt{b^2 - ac}}{a}
= \frac{\cos(x) \pm 1}{\sin^2(x)}
\end{displaymath}

Solve it:

\begin{displaymath}
y = - \frac1{\sin(x)} \pm \hbox{cotg}(x) + C
\end{displaymath}

The characteristic coordinates are the integration constants:

\begin{displaymath}
\xi = y + \frac1{\sin(x)} + \hbox{cotg}(x) \quad
\eta = y + \frac1{\sin(x)} - \hbox{cotg}(x)
\end{displaymath}

It does not look like the partial differential is going to be very simple.



18.8.2 Parabolic equations in two dimensions

In the parabolic case there is only one equation for the characteristics because the discriminant $b^2-ac$ is zero:

\begin{displaymath}
\frac{{\rm d}y}{{\rm d}x} = \frac{b}{a}
\end{displaymath}

So you can only find one characteristic coordinates, call it $\eta$.

You will need to take the other coordinate something else, say $\xi$ $\vphantom0\raisebox{1.5pt}{$=$}$ $x$. You want to take something simple, but it should be independent of the other coordinate.

The partial differential equations then simplifies to the two-di­men­sion­al canonical form

\begin{displaymath}
\fbox{$\displaystyle
a' u_{\xi\xi} = d'
$}
\end{displaymath} (18.26)

You may be surprised by that. In choosing $\eta$, all we did was make the coefficient $c'$ zero. We did not explicitly make $b'$ zero. But $b'$ is zero automatically. The reason is that the physical properties of partial differential equations do not change just because you use different coordinates. A parabolic equation should stay parabolic; there are fundamental differences between the physical behaviors of parabolic, elliptic, and hyperbolic equations. And the equation above would not be parabolic if $b'$ was nonzero.


Example

Question: Reduce the equation

\begin{displaymath}
x u_{xx} + 2 \sqrt{xy} u_{xy} + y u_{yy} - u_y = 0
\end{displaymath}

to two-di­men­sion­al canonical form.

Solution:

Write the equation for the characteristics

\begin{displaymath}
\frac{{\rm d}y}{{\rm d}x} = \frac{b\pm\sqrt{b^2 - ac}}{a}
= \sqrt{\frac yx}
\end{displaymath}

The square root is zero, so the equation is parabolic.

Solve the equation and call the integration constant $\eta$:

\begin{displaymath}
\frac{{\rm d}y}{\sqrt y} = \frac{{\rm d}x}{\sqrt x} \quad\quad\Rightarrow\quad\quad
\sqrt y = \sqrt x + \eta
\end{displaymath}

So take the new coordinates as

\begin{displaymath}
\xi = x \qquad \eta = \sqrt y - \sqrt x
\end{displaymath}

The final partial differential equation then becomes

\begin{displaymath}
- 4 \xi u_{\xi\xi} =
\left[\frac{1}{\sqrt{\xi}} + \frac{3}{\eta+\sqrt{\xi}}\right] u_{\eta}
\end{displaymath}



18.8.3 Elliptic equations in two dimensions

Characteristic lines are solutions to the ordinary differential equation

\begin{displaymath}
\frac{{\rm d}y}{{\rm d}x} = \frac{b\pm\sqrt{b^2 - ac}}{a}
\end{displaymath}

Elliptic equations have no real characteristics, because the square root is imaginary. However, elliptic equations can still be simplified, assuming that the above ordinary differential equation can be solved analytically.

Take either sign of the square root. Solve the equation and call the integration constant, say, $\tilde\xi$. Then write this integration constant in the form

\begin{displaymath}
\fbox{$\displaystyle
\tilde\xi = \xi + {\rm i}\eta
$}
\end{displaymath} (18.27)

where $\xi$ and $\eta$ are real and ${\rm i}$ $\vphantom0\raisebox{1.5pt}{$=$}$ $\sqrt{-1}$. In other words, take $\xi$ $\vphantom0\raisebox{1.5pt}{$=$}$ $\Re(\tilde\xi)$ and $\eta$ $\vphantom0\raisebox{1.5pt}{$=$}$ $\Im(\tilde\xi)$.

Using $\xi$ and $\eta$ as the new coordinates, it turns out that the partial differential equation takes the two-di­men­sion­al canonical form:

\begin{displaymath}
\fbox{$\displaystyle
a' u_{\xi\xi} + a' u_{\eta\eta} = d'
$}
\end{displaymath} (18.28)

You may note that this is quite similar to what you can get from rotating the coordinate system, as in the previous section. However, the above procedure works even if the coefficients $a$, $b$, and $c$ of the original partial differential equation are not constants.

There are significant limitations on this procedure, however, {D.6}


Example

Question: Reduce the equation

\begin{displaymath}
u_{xx} + (1+y)^2 u_{yy} = 0
\end{displaymath}

to two-di­men­sion­al canonical form.

Solution:

Write the equation for the characteristics

\begin{displaymath}
\frac{{\rm d}y}{{\rm d}x} = \frac{b\pm\sqrt{b^2 - ac}}{a}
= \pm {\rm i}(1+y)
\end{displaymath}

This is complex, so the equation is elliptic.

Solve using separation of variables

\begin{displaymath}
\frac{{\rm d}y}{1+y} = {\rm i}{ \rm d}x
\quad\quad\Rightarrow\quad\quad \ln\vert 1+y\vert - {\rm i}x = \tilde\xi
\end{displaymath}

The new coordinates can therefore be chosen as

\begin{displaymath}
\xi = \ln\vert 1+y\vert \qquad \eta = - x
\end{displaymath}

In terms of these coordinates, the equation becomes

\begin{displaymath}
a' u_{\xi\xi} + a' u_{\eta\eta} = d'
\end{displaymath}

The new partial differential equation becomes

\begin{displaymath}
u_{\xi\xi} + u_{\eta\eta} = - u_\xi
\end{displaymath}


18.8.3 Review Questions
1.

Convert the equation

\begin{displaymath}
10 u_{xx} + 6 u_{xy} + 2 u_{yy} = u_x + x + 1
\end{displaymath}

to two-di­men­sion­al canonical form.

Using rotation and stretching of the coordinates you would get

\begin{displaymath}
u_{\xi\xi} + u_{\eta\eta} = \frac{3}{\sqrt{110}} u_{\xi} - \...
...frac{3\sqrt{11}}{\sqrt{10}} \xi - \frac{1}{\sqrt{10}} \eta + 1
\end{displaymath}

Do you get the same equation? Should you? Comment.

Solution 2dcanel-a