In this paper, a DG (Discontinuous Galerkin) method which has been widely employed in CFD (Computational Fluid Dynamics) is used to solve the two-dimensional time-domain Maxwell’s equations for complex geometries on unstructured mesh. The element interfaces on solid boundary are treated in both curved way and straight way. Numerical tests are performed for both benchmark problems and complex cases with varying orders on a series of grids, where the high-order convergence in accuracy can be observed. Both the curved and the straight solid boundary implementation can give accurate RCS (Radar Cross-Section) results with sufficiently small mesh size, but the curved solid boundary implementation can significantly improve the accuracy when using relatively large mesh size. More importantly, this CFD-based high-order DG method for the Maxwell’s equations is very suitable for complex geometries.

In this paper we present a fully discrete A-ø finite element method to solve Maxwell’s equations with a nonlinear degenerate boundary condition, which represents a generalization of the classical Silver-Müller condition for a non-perfect conductor. The relationship between the normal components of the electric field E and the magnetic field H obeys a power-law nonlinearity of the type H x n = n x (|E x n|α-1E x n) with α ∈ (0,1]. We prove the existence and uniqueness of the solutions of the proposed A-ø scheme and derive the error estimates. Finally, we present some numerical experiments to verify the theoretical result.

An implicit discontinuous Galerkin method is introduced to solve the time-domain Maxwell’s equations in metamaterials. The Maxwell’s equations in metamaterials are represented by integral-differential equations. Our scheme is based on discontinuous Galerkin method in spatial domain and Crank-Nicolson method in temporal domain. The fully discrete numerical scheme is proved to be unconditionally stable. When polynomial of degree at most p is used for spatial approximation, our scheme is verified to converge at a rate of O(τ2+hp+1/2). Numerical results in both 2D and 3D are provided to validate our theoretical prediction.

An algebraic multilevel method is proposed for efficiently simulating linear wave propagation using higher-order numerical schemes. This method is used in conjunction with the Finite Volume Time Domain (FVTD) technique for the numerical solution of the time-domain Maxwell’s equations in electromagnetic scattering problems. In the multilevel method the solution is cycled through spatial operators of varying orders of accuracy, while maintaining highest-order accuracy at coarser approximation levels through the use of the relative truncation error as a forcing function. Higher-order spatial accuracy can be enforced using the multilevel method at a fraction of the computational cost incurred in a conventional higher-order implementation. The multilevel method is targeted towards electromagnetic scattering problems at large electrical sizes which usually require long simulation times due to the use of very fine meshes dictated by point-per-wavelength requirements to accurately model wave propagation over long distances.

A fully discrete discontinuous Galerkin method is introduced for solving time-dependent Maxwell’s equations. Distinguished from the Runge-Kutta discontinuous Galerkin method (RKDG) and the finite element time domain method (FETD), in our scheme, discontinuous Galerkin methods are used to discretize not only the spatial domain but also the temporal domain. The proposed numerical scheme is proved to be unconditionally stable, and a convergent rate is established under the L2-norm when polynomials of degree atmost r and k are used for temporal and spatial approximation, respectively. Numerical results in both 2-D and 3-D are provided to validate the theoretical prediction. An ultra-convergence of order in time step is observed numerically for the numerical fluxes w.r.t. temporal variable at the grid points.

In this paper, the time-dependent Maxwell’s equations used to modeling wave propagation in dispersive lossy bi-isotropic media are investigated. Existence and uniqueness of the modeling equations are proved. Two fully discrete finite element schemes are proposed, and their practical implementation and stability are discussed.

In this paper an explicit finite-difference time-domain scheme for solving the Maxwell’s equations in non-staggered grids is presented. The proposed scheme for solving the Faraday’s and Ampere’s equations in a theoretical manner is aimed to preserve discrete zero-divergence for the electric and magnetic fields. The inherent local conservation laws in Maxwell’s equations are also preserved discretely all the time using the explicit second-order accurate symplectic partitioned Runge-Kutta scheme. The remaining spatial derivative terms in the semi-discretized Faraday’s and Ampere’s equations are then discretized to provide an accurate mathematical dispersion relation equation that governs the numerical angular frequency and the wavenumbers in two space dimensions. To achieve the goal of getting the best dispersive characteristics, we propose a fourth-order accurate space centered scheme which minimizes the difference between the exact and numerical dispersion relation equations. Through the computational exercises, the proposed dual-preserving solver is computationally demonstrated to be efficient for use to predict the long-term accurate Maxwell’s solutions.

In this paper, hp-adaptive finite element methods are studied for time-harmonic Maxwell’s equations. We propose the parallel hp-adaptive algorithms on conforming unstructured tetrahedral meshes based on residual-based a posteriori error estimates. Extensive numerical experiments are reported to investigate the efficiency of the hp-adaptive methods for point singularities, edge singularities, and an engineering benchmark problem of Maxwell’s equations. The hp-adaptive methods show much better performance than the h-adaptive method.

Inverse scattering problem is discussed for the Maxwell’s equations. A reduction of theMaxwell’s system to a new Fredholm second-kind integral equation with a scalarweakly singular kernel is given for electromagnetic (EM) wave scattering. Thisequation allows one to derive a formula for the scattering amplitude in which only ascalar function is present. If this function is small (an assumption that validates aBorn-type approximation), then formulas for the solution to the inverse problem areobtained from the scattering data: the complex permittivityϵ′(x) in a bounded regionD ⊂ R3 is found from the scattering amplitudeA(β,α,k) known for a fixed k = ω √ϵ0μ0 >0 and allβ,α ∈ S2, whereS2 is the unit sphere in R3,ϵ0 and μ0 are constantpermittivity and magnetic permeability in the exterior regionD′ = R3\D. The novel pointsin this paper include: i) A reduction of the inverse problem for vectorEM waves to a vector integral equation with scalar kernelwithout any symmetry assumptions on the scatterer, ii) A derivation of thescalar integral equation of the first kind for solving the inversescattering problem, and iii) Presenting formulas for solving this scalar integralequation. The problem of solving this integral equation is an ill-posed one. A method fora stable solution of this problem is given.

Since matrix compression has paved the way for discretizing the boundary integralequation formulations of electromagnetics scattering on very fine meshes, preconditionersfor the resulting linear systems have become key to efficient simulations. Operatorpreconditioning based on Calderón identities has proved to be a powerful device fordevising preconditioners. However, this is not possible for the usual first-kind boundaryformulations for electromagnetic scattering at general penetrable composite obstacles. Wepropose a new first-kind boundary integral equation formulation following the reasoningemployed in [X. Clayes and R. Hiptmair, Report 2011-45, SAM, ETH Zürich (2011)] foracoustic scattering. We call it multi-trace formulation, because itsunknowns are two pairs of traces on interfaces in the interior of the scatterer. We give acomprehensive analysis culminating in a proof of coercivity, and uniqueness and existenceof solution. We establish a Calderón identity for the multi-trace formulation, which formsthe foundation for operator preconditioning in the case of conforming Galerkin boundaryelement discretization.

The symmetric energy-conserved splitting FDTD scheme developed in is a very new and efficient scheme for computing the Maxwell’s equations. It is based on splitting the whole Maxwell’s equations and matching the x-direction and y-direction electric fields associated to the magnetic field symmetrically. In this paper, we make further study on the scheme for the 2D Maxwell’s equations with the PEC boundary condition. Two new energy-conserved identities of the symmetric EC-S-FDTD scheme in the discrete H1-norm are derived. It is then proved that the scheme is uncondi-tionally stable in the discrete H1-norm. By the new energy-conserved identities, the super-convergence of the symmetric EC-S-FDTD scheme is further proved that it is of second order convergence in both time and space steps in the discrete H1-norm. Numerical experiments are carried out and confirm our theoretical results.

In this paper, we study high order discretization methods for solving the Maxwell equations on hybrid triangle-quad meshes. We have developed high order finite edge element methods coupled with different high order time schemes and we compare results and efficiency for several schemes. We introduce in particular a class of simple high order low dissipation time schemes based on a modified Taylor expansion.

In this paper we present a rigorous derivation of the material parameters for both the cylinder and rectangle cloaking structures. Numerical results using these material parameters are presented to demonstrate the cloaking effect.

In this paper, we consider an interior penalty discontinuous Galerkin (DG) method for the time-dependent Maxwell’s equations in cold plasma. In Huang and Li (J. Sci. Comput., 42 (2009), 321-340), for both semi and fully discrete DG schemes, we proved error estimates which are optimal in the energy norm, but sub-optimal in the L2-norm. Here by filling this gap, we show that these schemes are optimally convergent in the L2-norm on quasi-uniform tetrahedral meshes if the solution is sufficiently smooth.

In this paper, we develop a multi-symplectic wavelet collocation method for three-dimensional (3-D) Maxwell’s equations. For the multi-symplectic formulation of the equations, wavelet collocation method based on autocorrelation functions is applied for spatial discretization and appropriate symplectic scheme is employed for time integration. Theoretical analysis shows that the proposed method is multi-symplectic, unconditionally stable and energy-preserving under periodic boundary conditions. The numerical dispersion relation is investigated. Combined with splitting scheme, an explicit splitting symplectic wavelet collocation method is also constructed. Numerical experiments illustrate that the proposed methods are efficient, have high spatial accuracy and can preserve energy conservation laws exactly.