Simulation of electromagnetic wave propagation in metamaterials leads to more complicated time domain Maxwell's equations than the standard Maxwell's equations in free space. In this paper, we develop and analyze a non-dissipative discontinuous Galerkin (DG) method for solving the Maxwell's equations in Drude metamaterials. Previous discontinuous Galerkin methods in the literature for electromagnetic wave propagation in metamaterials were either non-dissipative but sub-optimal, or dissipative and optimal. Our method uses a different and simple choice of numerical fluxes, achieving provable non-dissipative stability and optimal error estimates simultaneously. We prove the stability and optimal error estimates for both semi- and fully discrete DG schemes, with the leap-frog time discretization for the fully discrete case. Numerical results are given to demonstrate that the DG method can solve metamaterial Maxwell's equations effectively.

In this paper, a new type of finite difference Hermite weighted essentially non-oscillatory (HWENO) schemes are constructed for solving Hamilton-Jacobi (HJ) equations. Point values of both the solution and its first derivatives are used in the HWENO reconstruction and evolved via time advancing. While the evolution of the solution is still through the classical numerical fluxes to ensure convergence to weak solutions, the evolution of the first derivatives of the solution is through a simple dimension-by-dimension non-conservative procedure to gain efficiency. The main advantages of this new scheme include its compactness in the spatial field and its simplicity in the reconstructions. Extensive numerical experiments in one and two dimensional cases are performed to verify the accuracy, high resolution and efficiency of this new scheme.

This paper is concerned with an interior penalty discontinuous Galerkin (IPDG) method based on a flexible type of non-polynomial local approximation space for the Helmholtz equation with varying wavenumber. The local approximation space consists of multiple polynomial-modulated phase functions which can be chosen according to the phase information of the solution. We obtain some {approximation} properties for this space and \textit{a prior} $L^2$ error estimates for the \textit{h}-convergence of the IPDG method using duality argument. We also provide ample numerical examples to show that, building phase information into the local spaces often gives more accurate results comparing to using the standard polynomial spaces.

In this paper, we analyze implicit-explicit (IMEX) Runge-Kutta (RK) time discretization methods for solving linear convection-diffusion equations. The diffusion operator is treated implicitly via the embedded discontinuous Galerkin (EDG) method and the convection operator explicitly via the upwinding discontinuous Galerkin method.

In this paper, we apply the local discontinuous Galerkin (LDG) method to 2D Keller--Segel (KS) chemotaxis model. We improve the results upon (Y. Epshteyn and A. Kurganov, SIAM Journal on Numerical Analysis, 47 (2008), 368-408) and give optimal rate of convergence under special finite element spaces before the blow-up occurs (the exact solutions are smooth). Moreover, to construct physically relevant numerical approximations, we consider $P^1$ LDG scheme and develop a positivity-preserving limiter to the scheme, extending the idea in (Y. Zhang, X. Zhang and C.-W. Shu, Journal of Computational Physics, 229 (2010), 8918-8934). With this limiter, we can prove the $L^1$-stability of the numerical scheme. Numerical experiments are performed to demonstrate the good performance of the positivity-preserving LDG scheme. Moreover, it is known that the chemotaxis model will yield blow-up solutions under certain initial conditions. We numerically demonstrate how to find the approximate blow-up time by using the $L^2$-norm of the $L^1$-stable numerical solution.