Xingguang Jin

Eric T. Chung, Patrick Ciarlet, Xingguang Jin et al.

From a mathematical perspective, the extraordinary properties of metamaterials are often reflected in the coefficients of the governing partial differential equations (PDEs). These coefficients may fall outside the assumptions of classical theory, particularly when the effective dielectric permittivity and/or magnetic permeability are negative. This situation can transform a coercive operator into a non-coercive one, potentially leading to ill-posedness. In this paper, we utilize the Constraint Energy Minimizing Generalized Multiscale Finite Element Method (CEM-GMsFEM), specifically designed for time-harmonic electromagnetic wave problems, where the construction of auxiliary spaces in the original CEM-GMsFEM is tailored to accommodate the sign-changing setting. Based on the framework of \texttt{T}-coercivity theory and resolution conditions, we establish the inf-sup stability and provide an a priori error estimate for the proposed method. The numerical results demonstrate the effectiveness and robustness of our approach in handling such sophisticated coefficient profiles.

Eric T. Chung, Patrick Ciarlet, Xingguang Jin et al.

The mathematical formulation of sign-changing problems involves a linear second-order partial differential equation in the divergence form, where the coefficient can assume positive and negative values in different subdomains. These problems find their physical background in negative-index metamaterials, either as inclusions embedded into common materials as the matrix or vice versa. In this paper, we propose a numerical method based on the constraint energy minimizing generalized multiscale finite element method (CEM-GMsFEM) specifically designed for sign-changing problems. The construction of auxiliary spaces in the original CEM-GMsFEM is tailored to accommodate the sign-changing setting. The numerical results demonstrate the effectiveness of the proposed method in handling sophisticated coefficient profiles and the robustness of coefficient contrast ratios. Under several technical assumptions and by applying the \texttt{T}-coercivity theory, we establish the inf-sup stability and provide an a priori error estimate for the proposed method.

Xiang Zhong, Eric T. Chung, Xingguang Jin

Modeling time-harmonic Maxwell problems in heterogeneous media presents significant mathematical and computational challenges. Due to the inherent non-elliptic structure and non-coercive nature of Maxwell equations, conventional methods face severe numerical instabilities, particularly in high-contrast media and at high wave numbers. These challenges often lead to ill-conditioned discrete systems and prohibitively high computational costs, limiting their practical applicability. To overcome these challenges, we introduce an efficient multiscale framework for time-harmonic Maxwell equations with impedance boundary conditions in high-contrast media. A major novelty of this study lies in circumventing the need for an explicit divergence-free constraint on multiscale basis functions. To achieve this, an auxiliary space is constructed via local spectral problems incorporating a mass term and a Silver-Müller-type boundary penalty. This novel design guarantees the coercivity of the corresponding bilinear form and automatically excludes the kernel of the curl operator from the leading eigenspaces. Building upon the auxiliary space, we then construct the multiscale space by using a distinct bilinear form. By exploiting a resolution condition and establishing key norm relationships, we rigorously prove the coercivity of this modified bilinear form a crucial property that underpins the whole analysis. Theoretical analysis shows that, with appropriate oversampling, the method achieves $O (H)$ convergence independent of the local contrast and the approximation error increases with the wave number $k$. Extensive numerical experiments are reported to validate the effectiveness of the proposed approach.