Accelerating the pace of engineering and science

# Documentation

## Hyperbolic Equations

Using the same ideas as for the parabolic equation, hyperbolic implements the numerical solution of

$d\frac{{\partial }^{2}u}{\partial {t}^{2}}-\nabla \cdot \left(c\nabla u\right)+au=f,$

for x in Ω, with the initial conditions

$\begin{array}{c}u\left(x,0\right)={u}_{0}\left(x\right)\\ \frac{\partial u}{\partial t}\left(x,0\right)={v}_{0}\left(x\right)\end{array}$

for all x in Ω, and usual boundary conditions. In particular, solutions of the equation utt - cΔu = 0 are waves moving with speed $\sqrt{c}$.

Using a given triangulation of Ω, the method of lines yields the second order ODE system

$M\frac{{d}^{2}U}{d{t}^{2}}+KU=F$

with the initial conditions

$\begin{array}{c}{U}_{i}\left(0\right)={u}_{0}\left({x}_{i}\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\forall i\\ \frac{d}{dt}{U}_{i}\left(0\right)={v}_{0}\left({x}_{i}\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\forall i\end{array}$

after we eliminate the unknowns fixed by Dirichlet boundary conditions. As before, the stiffness matrix K and the mass matrix M are assembled with the aid of the function assempde from the problems

–∇ · (cu) + au = f and –∇ · (0∇u) + du = 0.