“…b(u) is the non-negative diffusion mobility and Ψ(u) is the homogeneous free energy density. Based on the development of the model in Cahn-Hilliard equation, it is not surprising that many different numerical methods have been proposed for the Cahn-Hilliard equation, including the finite element methods [4][5][6][7][8][9][10][11][12][13], discontinuous Galerkin methods [14][15][16], local discontinous Galerkin methods [3,17,18], multi-grid method [19][20][21], finite difference methods [22][23][24] and so on. Second, the Camassa-Holm equations, introduced by Camassa and Holm in [25,26] as a model for the shallow water motion on the flat surface, are one of the basic models of peaked solitons.…”