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.

In this paper we develop an adaptive edge finite element method based on a reliable and efficient recovery type a posteriori error estimator for time-harmonic Maxwell equations. The asymptotically exact a posteriori error estimator is based on the superconvergence result proved for the lowest-order edge element on triangular grids, where most pairs of triangles sharing a common edge form approximate parallelograms. The efficiency and robustness of the proposed method is demonstrated by extensive numerical experiments for electromagnetic cloaking problems with highly anisotropic permittivity and permeability.

In this work, we investigate wave transmission property through a zero index metamaterial (ZIM) waveguide embedded with rectangular dielectric defects. We show that total reflection and total transmission (cloaking) can be achieved by adjusting the geometric sizes and/or permittivities of the defects. Our work provides another possibility of manipulating wave propagation through ZIM in addition to the widely studied dielectric defects with cylindrical geometries.

In this paper, we study a time-domain spherical cloaking model recently introduced by Zhao and Hao (2009). This model is quite complicated and is composed of four coupled differential equations. Here we first prove the existence and uniqueness of a solution for this model. Then we obtain a stability result, which analysis is quite involved due to the coupling between the four variables. To our best knowledge, this is the first well-posedness study carried out for the time-domain cloaking model with metamaterials.

In this paper, a mathematical analysis of the popular carpet cloak model is carried out. The well-posedness of the model is first proved, and then a finite element time-domain method is developed for solving this model. Three carpet cloak simulations are demonstrated to show the effectiveness of our model and the developed algorithm. To the best of our knowledge, this is the first mathematical analysis carried out for the carpet cloak model. The numerical simulation using edge elements for the carpet cloak model is also original.