2-D isotropic finite difference time domain method

The numerical dispersion inherent in the conventional FDTD method causes the numerical phase speed to become a function of frequency and direction, which is the main source of error. Unlike the convention treatment to use 1-D finite difference scheme to approximate the spatial partial differential operator (PDO) in Maxwell´s equations, the use of the high-dimensional isotropic finite difference scheme in the FDTD method will greatly reduce the numerical anisotropy in the FDTD method, which demonstrates its superiority and applicability. In addition, the stability condition is strictly derived.