An efficient Monte Carlo interior penalty discontinuous Galerkin method for the time-harmonic Maxwell's equations with random coefficients

Xiaobing Feng,Junshan Lin,Cody Lorton

DOI: https://doi.org/10.48550/arXiv.1806.05221

2018-06-14

Abstract:This paper develops an efficient Monte Carlo interior penalty discontinuous Galerkin method for electromagnetic wave propagation in random media. This method is based on a multi-modes expansion of the solution to the time-harmonic random Maxwell equations. It is shown that each mode function satisfies a Maxwell system with random sources defined recursively. An unconditionally stable IP-DG method is employed to discretize the nearly deterministic Maxwell system and the Monte Carlo method combined with an efficient acceleration strategy is proposed for computing the mode functions and the statistics of the electromagnetic wave. A complete error analysis is established for the proposed multi-modes Monte Carlo IP-DG method. It is proved that the proposed method converges with an optimal order for each of three levels of approximations. Numerical experiments are provided to validate the theoretical results and to gauge the performance of the proposed numerical method and approach.

Numerical Analysis

What problem does this paper attempt to address?

### Problems the paper attempts to solve

This paper aims to solve the problem of numerical simulation of electromagnetic wave propagation in random media. Specifically, the author has developed an efficient Monte Carlo Interior Penalty Discontinuous Galerkin (MCIP - DG) method for solving the time - harmonic Maxwell equations with random coefficients. The following are the main problems the paper attempts to solve:

1. **Propagation of electromagnetic waves in random media**:
- The propagation of electromagnetic waves in random media (such as atmospheric and biological media) has been a long - studied topic because of its wide applications in communication, remote sensing, detection, and imaging, etc.
- In these applications, describing the statistical characteristics of the electromagnetic wave field is a practical requirement.

2. **Limitations of existing methods**:
- Despite the rapid development of modern computing power, full - scale numerical modeling of three - dimensional random Maxwell equations remains a challenge. This is not only related to the large - scale and uncertainty of the problem, but also involves the modeling of multiple - scattering effects.
- Existing methods, such as the direct Monte Carlo sampling technique or the stochastic Galerkin method based on Karhunen - Loève expansion or Wiener chaos expansion, are still too computationally expensive in solving three - dimensional random vector Maxwell equations.

3. **Need for an efficient numerical method**:
- A method that can handle large - scale problems and effectively reduce computational complexity is required.
- The proposed MCIP - DG method decomposes the original problem into a series of almost deterministic Maxwell systems through multi - mode expansion, and combines the Monte Carlo method and acceleration strategies to efficiently calculate the mode functions and the statistical characteristics of electromagnetic waves.

### Overview of the solution

To address the above challenges, the paper proposes the following solutions:

- **Multi - mode expansion**: Use multi - mode expansion to represent the electric field \( E \) as a power series of a small perturbation parameter \( \varepsilon \). Each mode function satisfies an almost deterministic Maxwell system, and the source term is defined recursively.
- **Unconditionally stable IP - DG method**: Adopt an unconditionally stable interior penalty discontinuous Galerkin method to discretize the almost deterministic Maxwell system.
- **Monte Carlo method and acceleration strategy**: Combine the Monte Carlo method and LU decomposition acceleration strategy to quickly calculate each mode function and its statistical characteristics.
- **Error analysis**: Establish a complete error analysis and prove the optimal convergence order of the proposed method at three different approximation levels.

Through these methods, the paper provides an efficient and accurate numerical tool for the problem of electromagnetic wave propagation in random media.