Abstract
In this paper, we propose new finite volume element schemes to numerically solve the improved Boussinesq equation with Stokes damping. The new schemes can inherit characteristic properties of the conservation of mass and the decrease of total energy from the improved Boussinesq equation with Stokes damping. Numerical experiments illustrate that the proposed schemes are second-order accuracy in space and time.