Abstract:In this paper, we address the approximation of the coupling problem for the wave equation and Maxwell's equations of electromagnetism in the time domain in terms of electric field by means of a nodal linear finite element discretization in space, combined with a classical explicit finite difference scheme for time discretization. Our study applies to a particular case where the dielectric permittivity has a constant value outside a subdomain, whose closure does not intersect the boundary of the domain where the problem is defined. Inside this subdomain, Maxwell's equations hold. Outside this subdomain, the wave equation holds, which may correspond to Maxwell's equations with a constant permittivity under certain conditions. We consider as a model the case of first-order absorbing boundary conditions. First-order error estimates are proven in the sense of two norms involving first-order time and space derivatives under reasonable assumptions, among which lies a CFL condition for hyperbolic equations. The theoretical estimates are validated by numerical computations, which also show that the scheme is globally of the second order in the maximum norm in time and in the least-squares norm in space.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the numerical solution of the coupling problem between Maxwell's equations and the wave equation in the time domain using the explicit P1 finite element method. Specifically, the research focuses on the following points:
1. **Coupling problem**: In a specific geometric region, Maxwell's equations and the wave equation are respectively applicable. The paper studies how to establish an effective coupling relationship between these two regions and ensure the stability and accuracy of the overall problem.
2. **Absorbing boundary conditions**: In order to simulate wave propagation in an unbounded domain, the paper adopts the first - order absorbing boundary conditions (Absorbing Boundary Conditions, ABC). These conditions can reduce the influence of boundary reflection on the calculation results, thereby more realistically simulating actual physical phenomena.
3. **Explicit time - discretization scheme**: The paper proposes an explicit time - discretization scheme, which combines the classical explicit finite - difference method and the finite element method. The advantage of this scheme is high computational efficiency and is suitable for large - scale parallel computing.
4. **Error estimate**: The paper proves the first - order error estimate in time and space of this method under reasonable assumptions (such as the CFL condition), and verifies these theoretical results through numerical experiments.
### Specific problem description
The coupling problem considered in the paper can be stated as follows:
- **Computational domain**: Let \(\Omega\) be a bounded Lipschitz region with boundary \(\Gamma\). Divide \(\Omega\) into two sub - regions: the inner sub - region \(\Omega_{\text{in}}\) and the outer sub - region \(\Omega_{\text{out}}\), where \(\Omega=\bar{\Omega}_{\text{in}}\cup\Omega_{\text{out}}\).
- **Dielectric constant**: The dielectric constant \(\epsilon\) satisfies:
\[
\epsilon\equiv1\quad\text{in}\quad\Omega_{\text{out}},\quad\text{and}\quad\epsilon\geq1\quad\text{in}\quad\Omega_{\text{in}}
\]
- **Initial conditions**: Given the initial fields \(e_0\in\{H^1(\Omega)\}^N\) and \(e_1\in H(\text{div},\Omega)\), satisfying \(\nabla\cdot(\epsilon e_0)=\nabla\cdot(\epsilon e_1)\equiv0\), and the external force term \(f\in L^2([0, T]; H(\text{div},\Omega))\), satisfying \(\nabla\cdot f\equiv0\).
- **Boundary conditions**: On the boundary \(\Gamma\), the first - order absorbing boundary condition is adopted:
\[
\partial_n e+\partial_t e = 0\quad\text{in}\quad\Gamma\times(0, T)
\]
### Mathematical model
The coupling problem considered in the paper can be represented as:
\[
\begin{cases}
\epsilon\partial_{tt}e+\nabla\times\nabla\times e = f&\text{in}\quad\Omega_{\text{in}}\times(0, T)\\
\nabla\cdot(\epsilon e)=0&\text{in}\quad\Omega_{\text{in}}\times(0, T)\\
\partial_{tt}e-\Delta e = f&\text{in}\quad\Omega_