Abstract:In this paper, we present Gauss's law-preserving spectral methods and their efficient solution algorithms for curl-curl source and eigenvalue problems in two and three dimensions arising from Maxwell's equations. Arbitrary order $H(curl)$-conforming spectral basis functions in two and three dimensions are firstly proposed using compact combination of Legendre polynomials. A mixed formulation involving a Lagrange multiplier is then adopted to preserve the Gauss's law in the weak sense. To overcome the bottleneck of computational efficiency caused by the saddle-point nature of the mixed scheme, we present highly efficient solution algorithms based on reordering and decoupling of the resultant linear algebraic system and numerical eigen-decomposition of one dimensional mass matrix. The proposed solution algorithms are direct methods requiring only several matrix-matrix or matrix-tensor products of $N$-by-$N$ matrices, where $N$ is the highest polynomial order in each direction. Compared with other direct methods, the computational complexities are reduced from $O(N^6)$ and $O(N^9)$ to $O(N^3)$ and $O(N^4)$ with small and constant pre-factors for 2D and 3D cases, respectively, and can further be accelerated to $O(N^{2.807})$ and $O(N^{3.807})$, when boosted with the Strassen's matrix multiplication algorithm. Moreover, these algorithms strictly obey the Helmholtz-Hodge decomposition, thus totally eliminate the spurious eigen-modes of non-physical zero eigenvalues. Extensions of the proposed methods and algorithms to problems in complex geometries with variable coefficients and inhomogeneous boundary conditions are discussed to deal with more general situations. Ample numerical examples for solving Maxwell's source and eigenvalue problems are presented to demonstrate the accuracy and efficiency of the proposed methods.

The bi-Lebedev scheme for the Maxwell eigenvalue problem with 3D bi-anisotropic complex media

Bifurcation Analysis of the Eigenstructure of the Discrete Single-curl Operator in Three-dimensional Maxwell's Equations with Pasteur Media

Eigenfunction expansion of dyadic Green's function in biisotropic medium and its applications

Green’s functions for transversely isotropic magnetoelectroelastic media

Fundamental Solutions for Transversely Isotropic Magneto-Electro-elastic Media and Boundary Integral Formulation

Green'S Functions for A Two-Phase Infinite Magneto-Electro-Elastic Plane

Bilateral Coplanar Waveguide and Periodic Microstrip Line in Bianisotropic Superstrate-Substrate Structures

A Novel Pseudospectral Time Domain Method for Complex EM Simulations

An Adaptive High-Order Transient Algorithm to Solve Large-Scale Anisotropic Maxwell’s Equations

Babich-Like Ansatz for Three-Dimensional Point-Source Maxwell's Equations in an Inhomogeneous Medium at High Frequencies.

Highly efficient Gauss's law-preserving spectral algorithms for Maxwell's double-curl source and eigenvalue problems based on eigen-decomposition

A mixed element scheme of Helmholtz transmission eigenvalue problem for anisotropic media

A new numerical method for scalar eigenvalue problems in heterogeneous, dispersive, sign-changing materials

A Novel Mixed Spectral Method and Error Estimates for Maxwell Transmission Eigenvalue Problems

A Novel Lanczos-Type Procedure for Computing Eigenelements of Maxwell and Helmholtz Problems

L-P Estimates For Maxwell'S Equations With Source Term

Three-Dimensional Finite-Element Analysis of Magnetotelluric Data Using Coulomb-Gauged Potentials in General Anisotropic Media

Energy norm error estimates and convergence analysis for a stabilized Maxwell's equations in conductive media

The General Solution of Three-Dimensional Problems in Magnetoelectroelastic Media

Eigenstates of the full Maxwell equations for a two-constituent composite medium and their application to a calculation of the local electric field of a time dependent point electric dipole in a flat-slabs microstructure