A fourth-order compact difference scheme for Maxwell´s equations

A compact, central-difference approximation, in conjunction with the Yee grid, is used to compute the spatial derivatives in Maxwell´s equations. To advance the semi-discrete equations, the four-stage Runge-Kutta integrator is invoked. This combination of spatial and temporal differencing leads to a scheme that is fourth-order accurate, conditionally stable and highly efficient. Moreover, the use of compact differencing allows one to apply the compact operator in the vicinity of a perfect conductor-an attribute not found in other higher-order methods. Results are provided that quantify the spectral properties of the method. Simulations are conducted on problem spaces that span one and three dimensions and whose domains are of the open and closed type. Results from these simulations are compared with exact, closed-form solutions; the agreement between these results is consistent with numerical analysis.