Based on the Padé approximation and multistep method, we propose an implicit high-order unconditionally stable complex envelope algorithm to solve the time-dependent Maxwell’s equations. Unconditional numerical stability can be achieved simultaneously with a high-order accuracy in time. As we adopt the complex envelope Maxwell’s equations, numerical dispersion and dissipation are very small even at comparatively large time steps. To verify the capability of our algorithm, we compare the results of the proposed method with the exact solutions.
© 2008 Optical Society of America
Equations on this page are rendered with MathJax. Learn more.