In this paper, we develop a nodal discontinuous Galerkin method for solving the time-dependent Maxwell’s equations when metamaterials are involved. Both semi- and fully-discrete schemes are constructed. Numerical stability and error estimate are proved for both schemes. Numerical results are presented to demonstrate that the method is not only efficient but also very effective in solving metamaterial Maxwell’s equatons.

The goal of this paper is to prove the well-posedness for the governing equations which are used for cylindrical cloaking simulation. A new time-domain finite element scheme is developed to solve the governing equations. Numerical results demonstrating the cloaking phenomenon with the cylindrical cloak are presented. We finally extend the analysis and simulation to an elliptical cloak model.

For decades, the widely used finite difference method on staggered grids, also known as the marker and cell (MAC) method, has been one of the simplest and most effective numerical schemes for solving the Stokes equations and Navier–Stokes equations. Its superconvergence on uniform meshes has been observed by Nicolaides (SIAM J Numer Anal 29(6):1579–1591, 1992), but the rigorous proof is never given. Its behavior on non-uniform grids is not well studied, since most publications only consider uniform grids. In this work, we develop the MAC scheme on non-uniform rectangular meshes, and for the first time we theoretically prove that the superconvergence phenomenon (i.e., second order convergence in the L2 norm for both velocity and pressure) holds true for the MAC method on non-uniform rectangular meshes. With a careful and accurate analysis of various sources of errors, we observe that even though the local truncation errors are only first order in terms of mesh size, the global errors after summation are second order due to the amazing cancellation of local errors. This observation leads to the elegant superconvergence analysis even with non-uniform meshes. Numerical results are given to verify our theoretical analysis.

In this paper we propose a finite element time-domain method for modeling the optical black holes (OBHs) coupled with the perfectly matched layer (PML) technique. Stability analysis is carried out for the proposed scheme. Simulations of cylindrical, elliptical and square black holes demonstrate that our method is quite effective in modeling OBHs in time domain. To our best knowledge, this is the first OBHs simulation realized by the finite element time-domain method.

Since the development of Yee scheme back in 1966, it has become one of the most popular simulation tools for modeling electromagnetic wave propagation in various situations.However, its rigorous error analysis on nonuniform rectangular type gridswas carried out until 1994 by Monk and Süli. They showed that theYee scheme is still second-order convergent on a nonuniform mesh even though the local truncation error is only of first order. In this paper, we extend their results to Maxwell’s equations in metamaterials by a simpler proof, and show the second-order superconvergence in space for the trueYee scheme instead of the only semi-discrete form discussed in Monk and Süli’s original work. Numerical results supporting our analysis are presented.