The numerical solution of the one-dimensional nonlinear Klein-Gordon equation on an unbounded domain is studied in this paper. Split local absorbing boundary (SLAB) conditions are obtained by the operator splitting method, then the original problem is reduced to an initial boundary value problem on a bounded computational domain, which can be solved by the finite difference method. Several numerical examples are provided to show the advantages and effectiveness of the given method, and some interesting collision behaviors are also observed.

This is an extension of our earlier work \cite{DS} in which a high order stable method was constructed for solving hyperbolic conservation laws on arbitrarily distributed point clouds. %The initial condition is given on such a point cloud, %and the algorithm solves for point values of the solution at %later time also on this point cloud. An algorithm of building a suitable polygonal mesh based on the random points was given and the traditional discontinuous Galerkin (DG) method was adopted on the constructed polygonal mesh. Numerical results in \cite{DS} show that the current scheme will generate spurious numerical oscillations when dealing with solutions containing strong shocks. In this paper, we adapt a simple weighted essentially non-oscillatory (WENO) limiter, originally designed for DG schemes on two-dimensional unstructured triangular meshes \cite{ZZSQ}, to our high order method on polygonal meshes. The objective of this simple WENO limiter is to simultaneously maintain uniform high order accuracy of the original method in smooth regions and control spurious numerical oscillations near discontinuities. The WENO limiter we adopt is particularly simple to implement and will not harm the conservativeness and compactness of the original method. Moreover, we also extend the maximum-principle-satisfying limiter for the scalar case and the positivity-preserving limiter for the Euler system to our method. Numerical results for both scalar equations and Euler systems of compressible gas dynamics are provided to illustrate the good behavior of these limiters.

No abstract uploaded!

The main purpose of this paper is to give stability analysis and error estimates of the local discontinuous Galerkin (LDG) methods coupled with three specific implicit-explicit (IMEX) Runge-Kutta time discretization methods up to third order accuracy, for solving one-dimensional time-dependent linear fourth order partial differential equations. In the time discretization, all the lower order derivative terms are treated explicitly and the fourth order derivative term is treated implicitly. By the aid of energy analysis, we show that the IMEX-LDG schemes are unconditionally energy stable, in the sense that the time step $\dt$ is only required to be upper-bounded by a constant which is independent of the mesh size $h$. The optimal error estimate is also derived by the aid of the elliptic projection and the adjoint argument. Numerical experiments are given to verify that the corresponding IMEX-LDG schemes can achieve optimal error accuracy.

For a class of convection-diffusion equations with variable diffusivity, we construct third order accurate discontinuous Galerkin (DG) schemes on both one and two dimensional rectangular meshes. The DG method with an explicit time stepping can well be applied to nonlinear convection–diffusion equations. It is shown that under suitable time step restrictions, the scaling limiter proposed in Liu and Yu (2014) [23]when coupled with the present DG schemes preserves the solution bounds indicated by the initial data, i.e., the maximum principle, while maintaining uniform third order accuracy. These schemes can be extended to rectangular meshes in three dimension. The crucial for all model scenarios is that an effective test set can be identified to verify the desired bounds of numerical solutions. This is achieved mainly by taking advantage of the flexible form of the diffusive flux and the adaptable decomposition of weighted cell averages. Numerical results are presented to validate the numerical methods and demonstrate their effectiveness.