PHASE ERROR CONTROL FOR FD TD METHODS OF 2ND AND 4TH-ORDER ACCURACY

For FD-TD methods we determine the spatial resolution of the discretized domain in terms of the total computation time and the desired phase error. It is shown that the spatial step should vary as DELTAx is similar to g [e(phi)/t(c)]1/s in order to maintain a prescribed phase error level e(phi) throughout the computation time t(c), where s (=2 or 4) is the spatial order of accuracy of the scheme and g is a geometric factor. Significantly, we show that the rule of thumb of using 10-20 points per wavelength to determine the spatial cell size for the standard scheme is not optimal. Our results are verified by numerical simulations in two dimensions with the Yee scheme and a new fourth-order accurate FD-TD scheme.