A new high order FDTD method based on wave equations

A (2M, 4) scheme of the high-order finite-difference time-domain (FDTD) method based on the second order wave equation is proposed. This scheme is of the fourth order accuracy in time domain using the symplectic integrator propagator (SIP), and is of the 2Mth order accuracy in space domain using the discrete singular convolution (DSC) method. This scheme is similar to the traditional (2M, 4) scheme where SIP and DSC method is used and the distinction between them are listed as follows. Firstly, this new scheme is based on the wave equation. Secondly, computational region is discretized by uniform mesh rather than Yee mesh. Thirdly, the new scheme can save more memory and computation time because only one field component is used for calculation. Numerical examples show that when compared wi