In this paper we introduce a definition of the local conservation property for numerical methods solving time dependent conservation laws, which generalizes the classical local conservation definition. The motivation of our definition is the Lax-Wendroff theorem, and thus we prove it for locally conservative numerical schemes per our definition in one and two space dimensions. Several numerical methods, including continuous Galerkin methods and compact schemes, which do not fit the classical local conservation definition, are given as examples of locally conservative methods under our generalized definition.

Since 2000, the study of metamaterial has been a very hot topic due to its potential applications in many areas such as design of invisibility cloak and sub-wavelength imaging. Although several metamaterial models are often used by physicists and engineers, the study of their mathematical properties has lagged behind. In this paper, we initiate our investigation in the plasma-Lorentz model. More specifically, we first discuss the well-posedness of this model, then develop two fully-discrete finite element methods for solving it. Detailed stability and error analysis are carried out, and 3-D numerical results justifying our theoretical analysis are presented.

No abstract uploaded!

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.

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.