A high‐order implicit hybridizable discontinuous Galerkin method for the Benjamin‐Bona‐Mahony equation

Abstract

SummaryThis article introduces an efficient numerical scheme for the solution of the Benjamin‐Bona‐Mahony equation. The scheme is fully implicit, and combines high‐order time stepping and a hybridizable discontinuous Galerkin method in space. Implicit time stepping allows for larger time steps, while remaining robust for long‐time simulations. It is shown that the scheme is energy stable, and conserves mass. Computational experiments demonstrate that the scheme is high‐order accurate, robust, and possesses good conservation properties.