InthispaperwedevelopaconservativelocaldiscontinuousGalerkin(LDG) method for the Schro ̈ dinger-Korteweg-de Vries (Sch-KdV) system, which arises in var- ious physical contexts as a model for the interaction of long and short nonlinear waves. Conservative quantities in the discrete version of the number of plasmons, energy of the oscillations and the number of particles are proved for the LDG scheme of the Sch-KdV system. Semi-implicit time discretization is adopted to relax the time step constraint from the high order spatial derivatives. Numerical results for accuracy tests of stationary traveling soliton, and the collision of solitons are shown. Numerical ex- periments illustrate the accuracy and capability of the method.

In this paper, we develop, analyze and test the Fourier spectral methods for solving the Degasperis-Procesi(DP) equation which contains nonlinear high order derivatives, and possibly discontinuous or sharp transition solutions. The $L^2$ stability is obtained for general numerical solutions of the Fourier Galerkin method and Fourier collocation (pseudospectral) method.By applying the Gegenbauer reconstruction technique as a post-processing method to the Fourier spectral solution, we reduce the oscillations arising from the discontinuity successfully. The numerical simulation results for different types of solutions of the nonlinear Degasperis-Procesi equation are provided to illustrate the accuracy and capability of the methods.

In this paper, we develop an energy stable fully discrete discontinuous Galekin (DG) finite element method for the thin film epitaxy problem. Based on the method of lines, we construct and prove the energy stability of the spatial semi-discrete DG scheme firstly. To avoid the strict time step restriction of the explicit time integration method, the first order convex splitting method is used to get an unconditionally stable fully discrete DG method, which is a linearly implicit scheme for this nonlinear problem. The energy stability of the fully discrete convex splitting DG scheme is also proved. To improve the temporal accuracy, spectral deferred correction (SDC) method is adapted to achieve the high order accuracy in both time and space. Combining with the convex splitting method, the SDC method can be linearly solvable, high order accurate and stable in our numerical tests. These advantages ensure that the resulting fully discrete DG scheme is efficient to perform the long time simulation of the thin film epitaxy model. Numerical experiments of the accuracy and long time simulation show the capability and efficiency of the method.

In this paper, we develop the high order weighted essentially non-oscillatory (WENO) schemes for solving the Degasperis-Procesi (DP) equation, including finite volume (FV) and finite difference (FD) methods. The DP equation contains nonlinear high order derivatives, and possibly discontinuous or sharp transition solutions. The finite volume method is de- signed based on the total variation bounded property of the DP equation. And the finite difference method is constructed based on the L2 stability of the DP equation. Due to the adoption of the WENO reconstruction, both schemes are arbitrary high order accuracy and shock capturing. The numerical simulation results for different types of solutions of the nonlinear Degasperis-Procesi equation are provided to illustrate the accuracy and capability of the methods.

In this paper we consider a discontinuous Galerkin discretization of the ideal magnetohydrodynamics (MHD) equations on unstructured meshes, and the divergence free constraint ($\nabla \cdot \B = 0$) of its magnetic field $\B$. We first present two approaches for maintaining the divergence free constraint, namely the approach of a locally divergence free projection inspired by locally divergence free elements \cite{Li2005}, and another approach of the divergence cleaning technique given by Dedner et al. \cite{Dedner2002}. By combining these two approaches we obtain a efficient method at the almost same numerical cost. Finally, numerical experiments are performed to show the capacity and efficiency of the scheme.