Numerical homogenization for indefinite time-harmonic Maxwell equations

Abstract: We propose a novel numerical homogenization method based on the edge multiscale approach for solving indefinite time-harmonic Maxwell equations in heterogeneous media with large wavenumber. Numerical methods for these equations in homogeneous media with high wavenumber are particularly challenging due to the so-called pollution effect: the mesh size must be significantly smaller than the reciprocal of the wavenumber to achieve a desired accuracy. This challenge is amplified in heterogeneous media, which frequently occur in practical applications such as metamaterial simulations, since resolving the heterogeneity is necessary for obtaining reliable solutions. Our approach overcomes this difficulty by avoiding explicit resolution of the heterogeneity, while employing a mesh size that depends almost linearly on the reciprocal of the wavenumber. The approximation properties and stability of the method rely critically on the development and rigorous analysis of a novel, nonstandard variational formulation, which constitutes the main innovation of this work. Extensive numerical experiments are provided to validate our theoretical findings.

• The paper introduces an edge-oriented multiscale homogenization method that mitigates the pollution effect in high-frequency Maxwell simulations in heterogeneous media.
• The paper employs a local-global decomposition with Haar wavelet-based trace approximations to build a coarse-scale system with near-linear complexity.
• The paper provides rigorous a priori error analysis with explicit wavenumber-dependent estimates, supported by extensive 3D numerical experiments.

Numerical Homogenization for Indefinite Time-Harmonic Maxwell Equations: Theory, Algorithms, and Error Analysis

Introduction and Problem Statement

This work presents a new numerical homogenization methodology targeting time-harmonic Maxwell equations in highly heterogeneous media at elevated wavenumbers. The primary challenge addressed is the so-called “pollution effect” in high-frequency regimes: as the wavenumber $k$ increases, conventional edge-based finite element methods require the mesh size $h$ to be significantly smaller than $k^{-1}$ to maintain accuracy, resulting in prohibitive computational cost. The complication is further aggravated in the presence of strong heterogeneities, such as those encountered in metamaterial and composite applications, which impose additional finescale mesh constraints.

The authors propose a multiscale framework based on an edge-oriented strategy that enables the use of coarse meshes which do not resolve the heterogeneity explicitly, while retaining accuracy even for large $k$. The approach promises complexity control that is almost linear in $k^{-1}$ and independent of the heterogeneity scale.

Methodological Innovations

Variational Setting and Regularity

The authors start by rigorously formulating the impedance Maxwell problem in a three-dimensional Lipschitz domain $D\subset\mathbb{R}^3$ with rapidly varying, uniformly positive permittivity $\epsilon(x)$. They establish regularity assumptions on the solution $u\in \mathbf{H}_{imp}(\text{curl}; D)$, invoking duality-based regularity results suitable for low-regularity right-hand sides and relating higher-order Sobolev spaces to boundary trace regularity. They further introduce a nonstandard variational formulation—central to their analysis—by leveraging a transposition method inspired by Lions and Magenes to facilitate approximation error estimates for solutions with minimal regularity.

Local-Global Edge Multiscale Decomposition

The core methodology utilizes a local-global decomposition. The domain is covered by overlapping coarse neighborhoods $\{\omega_i\}$ (associated with a coarse mesh), and the global solution is written as a sum of local components supported on these neighborhoods. The key technical step is to represent each local field via its tangential traces on the boundary $\partial\omega_i$, and then to approximate these traces through multilevel Haar wavelet expansions. This constructs a multiscale ansatz space with hierarchical enrichment levels, where only a limited number of wavelet levels are typically required for accuracy.

A partition of unity is employed to guarantee global conformity and to combine local solutions into a globally admissible vector field in $h$0. The construction is compatible with edge degrees of freedom, ensuring the method's suitability for indefinite Maxwell problems.

Numerical Galerkin and Homogenization

Within the multiscale ansatz space, a global Galerkin formulation is advanced. The ansatz functions are built from local solutions to auxiliary Maxwell problems corresponding to Haar wavelet-based traces on $h$1. The method therefore bypasses direct resolution of heterogeneity within the bulk, encoding essential information at the coarse-grid skeleton via localized boundary problems. The resulting coarse-scale system is assembled and solved for the global field.

Error Analysis and Wavenumber Explicit Estimates

A potent contribution of the paper is a mathematically sharp, wavenumber-explicit a priori error analysis. The following points are of particular technical interest:

• Nonstandard Variational Technique: The error analysis is established using an auxiliary variational formulation that allows control of local approximation errors in very weak norms, accommodating both high-contrast coefficients and minimal regularity solutions.
• Wavelet Approximation Rates: The local Haar wavelet projections provide explicit convergence rates in the $h$2 and energy norms, exploiting the low-regularity trace behavior at subdomain interfaces.
• Quasi-optimal Global Error Bounds: The Galerkin orthogonality, together with regularity results for the dual Maxwell problem and the established local error bounds, culminates in a global error bound. The error in the energy norm decays at a rate governed by the wavelet enrichment level $h$3, the regularity exponent $h$4, and the coarse mesh size $h$5, with

$h$6

under reasonable assumptions on $h$7 and $h$8 relative to $h$9.

This result is significant: the required mesh resolution is nearly optimal (almost linear in $k^{-1}$0) and the method's accuracy is proved to be robust to the coefficient heterogeneity scale.

Numerical Verification

The authors conduct extensive 3D numerical experiments to support their theoretical claims:

• Homogeneous and Heterogeneous Cases: Simulations include homogeneous media as well as inclusions with moderate to high contrast and media with highly oscillatory coefficients.
• Robustness to Large Wavenumbers: The method remains stable and accurate up to $k^{-1}$1 on coarse meshes that would be otherwise insufficient for conventional discretizations.
• Efficiency of Wavelet Enrichment: A single Haar wavelet enrichment drastically reduces the error, even in highly heterogeneous configurations. For example, relative $k^{-1}$2 errors decrease from over $k^{-1}$3 to under $k^{-1}$4 in some strong contrast cases, with corresponding reductions in $k^{-1}$5 error.
• Complexity Reduction: The size of the global linear system is reduced by more than an order of magnitude compared to full fine-mesh discretization.

The data clearly demonstrate that this multiscale edge-based method equipped with wavelet enrichment levels achieves accuracy and robustness unattainable by traditional coarse-scale approaches.

Implications, Connections, and Future Directions

This work addresses long-standing bottlenecks for high-frequency electromagnetic simulations in complex media. It provides a rigorous and implementable recipe for constructing multiscale spaces in which the pollution effect is systematically controlled, overcoming the typical scaling issues in classical methods. The approach is also generic, suggesting applicability to other indefinite and multiscale problems, including those with random or parametric coefficients, and aligns with ongoing trends in domain decomposition and model reduction.

The theoretical framework opens several avenues for further exploration:

• Extension to Larger $k^{-1}$6: The stability and approximation theory established here is a promising foundation for addressing even higher-frequency regimes, although system size becomes a limiting factor unless complemented by iterative solvers tailored to the multiscale ansatz space.
• Divergence-Free Multiscale Bases: The partition of unity construction currently loses the exact divergence-free property. Constructing globally divergence-conforming multiscale spaces without sacrificing locality or computational tractability remains an open theoretical and practical challenge.
• Integration with Fast Solvers: Coupling this numerical homogenization with domain decomposition or multilevel preconditioners—using the multiscale ansatz as the coarse space—could yield solvers with both fast convergence and nearly optimal complexity.

Conclusion

This paper introduces a mathematically rigorous and practically validated edge-multiscale numerical homogenization method for indefinite time-harmonic Maxwell equations in highly heterogeneous media with large wavenumber (2604.22502). Its use of local wavelet-based trace approximations and nonstandard variational theory yields robust, quasi-optimal global error estimates, significantly reducing the computational burden associated with high-frequency wave propagation in complex materials. The work represents a substantial advance in multiscale electromagnetics, providing tools and analysis that are expected to inform both future theoretical developments and large-scale simulation practice.