Optimized Operator-Splitting Methods in  Numerical Integration of Maxwell's Equations

Optimized operator splitting methods for numerical integration of the time domain Maxwell's equations in Computational ElectroMagnetics (CEM) are proposed for the first time. The methods are based on splitting the time domain evolution operator of Maxwell's equations into suboperators, and corresponding time coefficients are obtained by reducing the norm of truncation terms to a minimum. The general high-order staggered finite difference is introduced for discretizing the three-dimensional curl operator in the spatial domain. The detail of the schemes and explicit iterated formulas are also included.