The fourth order average vector field (AVF) method is applied to solve the “Good” Boussinesq equation. The semi-discrete system of the “good” Boussinesq equation obtained by the pseudo-spectral method in spatial variable, which is a classical finite dimensional Hamiltonian system, is discretizated by the fourth order average vector field method. Thus, a new high order energy conservation scheme of the “good” Boussinesq equation is obtained. Numerical experiments confirm that the new high order scheme can preserve the discrete energy of the “good” Boussinesq equation exactly and simulate evolution of different solitary waves well.
