Subsections


21.2 Extension to finite regions

If $x$ is restricted by finite boundaries, the D'Alembert solution does not really apply. To use it anyway, we must somehow extend the problem to a doubly infinite $x$-range without boundaries. But our solution without boundaries should still satisfy the boundary conditions for the finite range. That is often possible by clever use of symmetry. An example can clarify that.


21.2.1 The physical problem

The problem is to find the pressure for sound wave propagation in a tube with one end closed and one end open:

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


21.2.2 The mathematical problem


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


21.2.3 Dealing with the boundary conditions

The D'Alembert solution applies to an infinite domain $-\infty$ $\raisebox{.3pt}{$<$}$ $x$ $\raisebox{.3pt}{$<$}$ $\infty$. So to use the D'Alembert solution, the given initial conditions, that are valid for 0 $\raisebox{.3pt}{$<$}$ $x$ $\raisebox{.3pt}{$<$}$ $\ell$ must be extended to all $x$. In other words, functions $f(x)$ and $g(x)$ must be converted into functions $\bar{f}(x)$ and $\bar{g}(x)$ that have values for all $x$. Of course, in the interval 0 $\raisebox{.3pt}{$<$}$ $x$ $\raisebox{.3pt}{$<$}$ $\ell$, they must stay the same as $f(x)$ and $g(x)$. Assume now for example that $f(x)$ looks as sketched below:

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

You might think that you could now simply take $\bar{f}(x)$ to be zero for all $x$ outside the range of the pipe. The corresponding D'Alembert solution will satisfy the wave equation everywhere, including inside the pipe 0 $\raisebox{.3pt}{$<$}$ $x$ $\raisebox{.3pt}{$<$}$ $\ell$. That is good, because the wave equation must indeed be satisfied. Unfortunately, the solution you get that way will not satisfy the boundary conditions at $x$ $\vphantom0\raisebox{1.5pt}{$=$}$ 0 and $x$ $\vphantom0\raisebox{1.5pt}{$=$}$ $\ell$. So it will still be wrong.

You must select the extension $\bar{f}(x)$ of $f(x)$ to all $x$ so that the correct boundary conditions become automatic.

The way to do it is as follows:

The process is shown for $f(x)$ below:

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


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


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


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


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

Create the extended function $\bar g(x)$ or $G$ the same way.

It is OK if you get kinks or discontinuities in your functions $\bar
f$ and $\bar g$ while creating (anti)symmetry. This happens when $f$ and/or $g$ does not satisfy the given boundary conditions. While then $u$ or $u_x$ may not have a unique value at the initial time, that problem will disappear when the time becomes greater than zero.


21.2.4 The final solution


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


\begin{displaymath}
u(x,t) = \frac{\bar f(x-at) + \bar f(x+at)}{2}
+ \frac{1}{2a} \int_{x-at}^{x+at} \bar g(\xi) { \rm d}\xi
\end{displaymath}

This is pretty easy to evaluate for simple functions $f$ and $g$. You will have fun doing it.


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

In the range 0 $\raisebox{-.3pt}{$\leqslant$}$ $x$ $\raisebox{-.3pt}{$\leqslant$}$ $\ell$, the found solution is exactly the same as for the finite pipe! The solution outside that range can simply be ignored.