Application of Block Krylov Subspace Spectral Methods to Maxwell's Equations
- Stanford University, Department of Energy Resources Engineering, Stanford, CA 94305-2220 (United States)
Ever since its introduction by Kane Yee over forty years ago, the finite-difference time-domain (FDTD) method has been a widely-used technique for solving the time-dependent Maxwell's equations. This paper presents an alternative approach to these equations in the case of spatially-varying electric permittivity and/or magnetic permeability, based on Krylov subspace spectral (KSS) methods. These methods have previously been applied to the variable-coefficient heat equation and wave equation, and have demonstrated high-order accuracy, as well as stability characteristic of implicit time-stepping schemes, even though KSS methods are explicit. KSS methods for scalar equations compute each Fourier coefficient of the solution using techniques developed by Gene Golub and Gerard Meurant for approximating elements of functions of matrices by Gaussian quadrature in the spectral, rather than physical, domain. We show how they can be generalized to coupled systems of equations, such as Maxwell's equations, by choosing appropriate basis functions that, while induced by this coupling, still allow efficient and robust computation of the Fourier coefficients of each spatial component of the electric and magnetic fields. We also discuss the implementation of appropriate boundary conditions for simulation on infinite computational domains, and how discontinuous coefficients can be handled.
- OSTI ID:
- 21322933
- Journal Information:
- AIP Conference Proceedings, Vol. 1174, Issue 1; Conference: International multiconference of engineers and computer scientists 2009, Hong Kong (China), 28-30 Mar 2009; Other Information: DOI: 10.1063/1.3256248; (c) 2009 American Institute of Physics; Country of input: International Atomic Energy Agency (IAEA); ISSN 0094-243X
- Country of Publication:
- United States
- Language:
- English
Similar Records
Domain decomposed preconditioners with Krylov subspace methods as subdomain solvers
A CYLINDRICALLY SYMMETRIC UNIAXIAL PML MAXWELL SOLVER FOR TRANSIENT ATMOSPHERIC ELECTRICITY SIMULATIONS