next up previous
Next: Additional Components in Standard Up: Numerical solutions of PDEs Previous: Crank-Nicolson method

Methods for Hyperbolic Equations

The FTCS-Method can also be applied on Hyperbolic equations of the form $u_{tt}=u_{xx}$. By using the approximation for the second derivatives (31) for both sides of the equations we obtain:
$\displaystyle \frac{u_l^{m-1}-2u_l^m+u_{l}^{m+1}}{k^2}$ $\textstyle =$ $\displaystyle \frac{u_{l-1}^m-2u_l^m+u_{l+1}^m}{h^2}$  
$\displaystyle u_l^{m+1}$ $\textstyle =$ $\displaystyle \frac{k^2}{h^2}(u_{l-1}^m-2u_l^m+u_{l+1}^m)
+2u_l^m-u_l^{m-1}$ (35)

To become a stable approximation, the constraint $\displaystyle{0<\frac{k^2}{h^2}\le1}$ has to be hold. By setting $\displaystyle{r=\frac{k^2}{h^2}=1}$, (35) becomes:
$\displaystyle u_l^{m+1}$ $\textstyle =$ $\displaystyle u_{l-1}^m+u_{l+1}^m-u_l^{m-1}.$ (36)


\begin{picture}
% latex2html id marker 680
(160,110)
%vertical lines
\color{cyan...
...put(110,60){\circle*{4}}
%unfilled circles
\put(85,85){\circle{6}}
\end{picture}

\begin{picture}
% latex2html id marker 704
(160,110)
%vertical lines
\color{cyan...
...put(110,60){\circle*{4}}
%unfilled circles
\put(85,85){\circle{6}}
\end{picture}



Heiko Enderling 2003-11-13