Hyperbolic Equations

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

d2ut2(cu)+au=f,

for x in Ω, with the initial conditions

u(x,0)=u0(x)ut(x,0)=v0(x)

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 triangulation of Ω, the method of lines yields the second order ODE system

Md2Udt2+KU=F

with the initial conditions

Ui(0)=u0(xi)iddtUi(0)=v0(xi)i

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?