Hyperbolic Equations

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


for x in Ω, where x represents a 2-D or 3-D point, with the initial conditions


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

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


with the initial conditions


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.

Was this topic helpful?