Changqing Ye,Xingguang Jin,Patrick Ciarlet Jr.,Eric T. Chung
Abstract: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.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how to effectively conduct numerical simulations in sign - changing problems with high contrast and heterogeneity. Specifically, these problems can be mathematically represented as a second - order linear partial differential equation (PDE), whose coefficients can take positive and negative values in different sub - domains and there are discontinuities on the interfaces between sub - domains. Such problems are commonly found in the applications of negative - refractive - index metamaterials in physical contexts, such as embedding negative - refractive - index materials in ordinary materials or vice versa.
### Main research objectives:
1. **Propose a new numerical method**: The paper proposes a new numerical method based on the Constrained Energy Minimization Generalized Multiscale Finite Element Method (CEM - GMsFEM), which is specifically used to deal with sign - changing problems.
2. **Improve the construction of the auxiliary space**: In the original CEM - GMsFEM method, the construction of the auxiliary space needs to be adjusted to adapt to the characteristics of sign - changing problems. The paper ensures the positivity of eigenvalues by modifying the coefficients in the generalized eigenvalue problem, thereby guaranteeing the approximation ability of the auxiliary space.
3. **Verify the effectiveness and robustness of the method**: Through numerical experiments, the effectiveness and robustness of the proposed method in dealing with complex coefficient distributions and high - contrast coefficient ratios are verified.
### Mathematical models and methods:
- **Model problem**: Consider a two - dimensional Lipschitz domain \(\Omega\), which can be divided into two non - overlapping sub - domains \(\Omega^{+}\) and \(\Omega^{-}\), with the interface \(\Gamma\). The coefficient \(\sigma\) is positive in \(\Omega^{+}\), negative in \(\Omega^{-}\), and discontinuous on \(\Gamma\).
\[
\sigma_{\text{max}}^{+}=\text{ess sup}_{x\in\Omega^{+}}\sigma(x),\quad\sigma_{\text{min}}^{+}=\text{ess inf}_{x\in\Omega^{+}}\sigma(x)
\]
\[
\sigma_{\text{max}}^{-}=\text{ess sup}_{x\in\Omega^{-}}|\sigma(x)|,\quad\sigma_{\text{min}}^{-}=\text{ess inf}_{x\in\Omega^{-}}|\sigma(x)|
\]
It is required that \(\sigma_{\text{max}}^{+}\geq\sigma_{\text{min}}^{+}>0\) and \(\sigma_{\text{max}}^{-}\geq\sigma_{\text{min}}^{-}>0\).
- **Variational form**: For \(f\in L^{2}(\Omega)\), the variational form of the model problem is:
\[
\text{Find }u\in V\text{ such that }\forall v\in V,\quad a(u, v)=\int_{\Omega}f v\,dx
\]
where
\[
a(v, w)=\int_{\Omega}\sigma\nabla v\cdot\nabla w\,dx
\]
- **T - coercivity theory**: To ensure the well - posedness of the model problem, the T - coercivity theory is introduced. If there exists a bijective mapping \(T: V\rightarrow V\) and a positive number \(\alpha\) such that for all \(v\in V\):
\[
|a(v, Tv)|\geq\alpha\|v\|_{\tilde{a}}^{2}
\]
then according to the Banach–Necas–Babuška theory, the model problem has a unique solution.
### Method steps:
1. **Construct local auxiliary spaces**: In each coarse - grid cell...