  This paper addresses the numerical solution of the Westervelt equation, which
arises as one of the model equations in nonlinear acoustics. The problem is
rewritten in a canonical form that allows the systematic discretization by
Galerkin approximation in space and time. Exact energy preserving methods of
formally arbitrary order are obtained and their efficient realization as well
as the relation to other frequently used methods is discussed.
