Abstract:We propose a boundary integral formulation for the dynamic problem of electromagnetic scattering and transmission by homogeneous dielectric obstacles. In the spirit of Costabel and Stephan, we use the transmission conditions to reduce the number of unknown densities and to formulate a system of coupled boundary integral equations describing the scattered and transmitted waves. The system is transformed into the Laplace domain where it is proven to be stable and uniquely solvable. The Laplace domain stability estimates are then used to establish the stability and unique solvability of the original time domain problem. Finally, we show how the bounds obtained in both Laplace and time domains can be used to derive error estimates for semi discrete Galerkin discretizations in space and for fully discrete numerical schemes that use Convolution Quadrature for time discretization and a conforming Galerkin method for discretization of the space variables.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the transient problem of scattering and transmission of electromagnetic waves on uniform dielectric obstacles in regions with Lipschitz boundaries. Specifically, the author proposes a framework of boundary integral equations (BIE) to describe the scattering and transmission phenomena of electromagnetic waves between different media. The main contributions of the paper are as follows:
1. **Proposal of boundary integral equations**: By using the transmission conditions to reduce the number of unknown densities, the author constructs a set of coupled boundary integral equations for describing the scattered and transmitted waves.
2. **Stability analysis in the Laplace domain**: The system is transformed into the Laplace domain, and its stability and unique solvability are proved. Using the stability estimates in the Laplace domain, the stability and unique solvability of the original time - domain problem are further established.
3. **Error estimates**: Through the bounds in the Laplace domain and the time domain, the spatial error estimates of the semi - discrete Galerkin discretization method are derived, as well as the time error estimates of the fully discrete numerical scheme. The latter uses the convolution quadrature method for time discretization and the Galerkin - conforming method for spatial variable discretization.
### Background of the paper
Consider a bounded open region \(\Omega^{-}\subset\mathbb{R}^3\) whose boundary \(\Gamma\) is Lipschitz continuous, and define the complement \(\Omega^{+}:=\mathbb{R}^3\setminus\Omega^{-}\). Assume that \(\Omega^{-}\) and \(\Omega^{+}\) are filled with uniform materials with linear and constant electromagnetic properties respectively, which are represented by their permittivity \(\epsilon^{+}\) and permeability \(\mu^{+}\) (in the outer region \(\Omega^{+}\)) and permittivity \(\epsilon^{-}\) and permeability \(\mu^{-}\) (in the inner region \(\Omega^{-}\)) respectively.
When an electromagnetic wave propagates from \(\Omega^{+}\) to \(\Omega^{-}\), due to the difference in electromagnetic properties between the two regions, part of the incident wave will be reflected back to \(\Omega^{+}\), and another part will pass through the interface \(\Gamma\) and be transmitted to \(\Omega^{-}\). The incident electromagnetic perturbations include the electric field \(E_{\text{inc}}\) and the magnetic field \(H_{\text{inc}}\), and the total field is the superposition of the incident field and the scattered fields \(E_{\text{scat}}\) and \(H_{\text{scat}}\) in \(\Omega^{+}\), or the superposition of the transmitted fields \(E\) and \(H\) in \(\Omega^{-}\).
### Mathematical model
The mathematical model proposed in the paper includes the following equations:
\[
\begin{aligned}
& \nabla\times\nabla\times E_{\text{scat}}+c_+^{-2}\frac{\partial^{2}E_{\text{scat}}}{\partial t^{2}} = 0\quad\text{in}\ \Omega^{+}\times(0, T),\\
& \nabla\times\nabla\times E + c_-^{-2}\frac{\partial^{2}E}{\partial t^{2}} = 0\quad\text{in}\ \Omega^{-}\times(0, T),\\
& \mathbf{n}\times(E - E_{\text{scat}})\times\mathbf{n}=\mathbf{n}\times E_{\text{inc}}\times\mathbf{n}\quad\text{in}\ \Gamma\times(0, T),\\