Development and analysis of Crank‐Nicolson scheme for metamaterial Maxwell's equations on nonuniform rectangular grids

In this paper, we develop a fully implicit scheme on staggered grids to solve the Maxwell's equations when Drude metamaterial is involved. Unconditional stability and optimal error estimate of the scheme are proved. Numerical results are provided to support the theoretical analysis, and used to demonstrate the applicability of the scheme to simulate the complicated backward wave propagation phenomenon occurring in metamaterials.