pith. sign in

arxiv: 2407.17130 · v2 · pith:UZDX3PRGnew · submitted 2024-07-24 · 🧮 math.NA · cs.NA

Multiscale modeling for a class of high-contrast heterogeneous sign-changing problems

Eric T. Chung , Patrick Ciarlet Jr. , Xingguang Jin , Changqing Ye This is my paper

Pith reviewed 2026-05-23 22:46 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords sign-changing problemsCEM-GMsFEMhigh-contrast coefficientsT-coercivityinf-sup stabilitya priori error estimatenegative-index metamaterialsmultiscale finite elements
0
0 comments X

The pith

A version of CEM-GMsFEM adapted for sign-changing coefficients yields inf-sup stable approximations and a priori error bounds for high-contrast heterogeneous elliptic problems.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a numerical method for linear second-order elliptic equations whose coefficient takes both positive and negative values in different subdomains. It modifies the constraint energy minimizing generalized multiscale finite element method by constructing auxiliary spaces suited to the sign-changing case. Numerical experiments illustrate that the resulting scheme remains effective for complex coefficient profiles and robust under large contrast ratios. Under several technical assumptions, T-coercivity theory supplies an inf-sup stability result together with an a priori error estimate. These properties matter because standard Galerkin methods lose stability when the coefficient changes sign, while high-contrast media require multiscale reduction to keep computations feasible.

Core claim

The proposed CEM-GMsFEM discretization, with auxiliary spaces tailored to accommodate sign-changing coefficients, is inf-sup stable and admits an a priori error estimate for high-contrast heterogeneous sign-changing problems; numerical tests confirm effectiveness on sophisticated profiles and robustness with respect to contrast ratios.

What carries the argument

The constraint energy minimizing generalized multiscale finite element method (CEM-GMsFEM) whose auxiliary-space construction is altered to respect sign-changing coefficients, thereby producing stable multiscale basis functions.

If this is right

  • The discrete problem remains uniquely solvable for any sufficiently fine multiscale mesh under the stated assumptions.
  • The approximation error is bounded by a constant times the sum of the fine-scale and coarse-scale contributions, with the constant independent of contrast under the robustness result.
  • The method applies directly to coefficient profiles arising from negative-index metamaterial inclusions embedded in a matrix or vice versa.
  • Standard multiscale basis construction without the sign-changing modification fails to guarantee stability, motivating the tailored auxiliary spaces.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same auxiliary-space modification could be tested on three-dimensional domains or on coefficients drawn from random fields to check whether the observed robustness persists beyond the reported examples.
  • If T-coercivity constants remain bounded uniformly in contrast, the approach might serve as a building block for time-harmonic Maxwell or acoustic problems that also feature sign changes.
  • An explicit dependence of the error constant on the number of auxiliary functions per coarse element would allow quantitative comparison with other multiscale techniques for indefinite problems.

Load-bearing premise

Several technical assumptions hold and the T-coercivity theory applies to the sign-changing setting.

What would settle it

Numerical computation of the discrete inf-sup constant that tends to zero as the contrast ratio increases, for a fixed mesh sequence satisfying the method's assumptions, would falsify the stability claim.

Figures

Figures reproduced from arXiv: 2407.17130 by Changqing Ye, Eric T. Chung, Patrick Ciarlet Jr., Xingguang Jin.

Figure 1
Figure 1. Figure 1: Illustration of the nested meshes Kh and KH. A fine element τ , two coarse elements Ki ′ and Ki ′′ , accompanied by their corresponding oversampling regions K2 i ′ and K2 i ′′ , are colored differently. In the original CEM-GMsFEM, a generalized eigenvalue problem, find λ ∈ R and v ∈ H1 (Ki) \ {0} s.t. Z Ki σ∇v · ∇w dx = λ Z Ki µmsh diam(Ki) −2σvw dx, (6) is solved on each coarse element Ki , where µmsh is … view at source ↗
Figure 2
Figure 2. Figure 2: (a) The coefficient profile and the marked coarse element. (b)–(d) The plot of the first/second/third eigenfunction corresponding to the marked coarse element. 4 Numerical experiments We conduct numerical experiments on a square domain Ω = (0, 1) × (0, 1). The fine mesh Kh is generated by dividing Ω into 400 × 400 squares. Consequently, the coefficient profile σ is represented by a 400 × 400 matrix/image, … view at source ↗
Figure 3
Figure 3. Figure 3: The subplots are marked as (x-y), where x can take a, b, or c, corresponding to the results for the first, second, or third eigenfunction, respectively. If y is 1, 2, or 3, the subplot displays the multiscale basis with m oversampling layers, m equal to y. Alternatively, if y is 4, the subplot shows the relative differences (y-axis) in the energy and L 2 norm of the multiscale bases between m = 8 and m = 1… view at source ↗
Figure 4
Figure 4. Figure 4: Numerical results for the flat interface model with [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Numerical results for the flat interface model with ( [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: (a) The smooth source term f that is constructed as the superposition of four 2D Gaussian functions centered at (0.25, 0.25), (0.75, 0.25), (0.25, 0.75), and (0.75, 0.75) with a variance of 0.01. (b)-(c) The reference solutions correspond to the 10 × 10 and 20 × 20 periodic configurations, respectively [PITH_FULL_IMAGE:figures/full_fig_p011_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Subplots (a) and (b) display visualizations for the coefficient σ used in the periodic cross-shaped inclusion model with 10 × 10 and 20 × 20 periodic configurations, respectively. Subplots (c) and (d) demon￾strate the reference solutions for 10 × 10 and 20 × 20 periodic configurations, respectively. The numerical results of the Q1 FEM and the proposed method with m ∈ {1, 2, 3, 4} are presented in [PITH_FU… view at source ↗
Figure 8
Figure 8. Figure 8: Numerical results for the periodic cross-shaped inclusion model with different periodic configura [PITH_FULL_IMAGE:figures/full_fig_p013_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: (a) The illustration of inclusions (dark regions) of the random inclusion model. (b) The plot of the reference solution by setting (σ + ∗ , σ− ∗ ) = (1, 10−3 ). (c) The plot of the reference solution by setting (σ + ∗ , σ− ∗ ) = (1, 103 ). increases, the growth rate is expected to slow down, leading to a deterioration of the approximation quality. Surprisingly, this deterioration commonly appears in severa… view at source ↗
Figure 10
Figure 10. Figure 10: Numerical results for the random inclusion model. Subplots [PITH_FULL_IMAGE:figures/full_fig_p015_10.png] view at source ↗
read the original 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.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The paper proposes a tailored CEM-GMsFEM for linear second-order sign-changing divergence-form PDEs with high-contrast heterogeneous coefficients, motivated by negative-index metamaterials. It asserts that, under several technical assumptions, T-coercivity yields inf-sup stability and an a priori error estimate; numerical experiments are said to demonstrate effectiveness on sophisticated profiles and robustness to contrast ratios.

Significance. If the technical assumptions hold and T-coercivity applies without additional restrictions on contrast or interface geometry, the result would supply the first multiscale method with provable stability for this class of sign-changing problems. The extension of the CEM-GMsFEM auxiliary-space construction is a natural technical step, and the reported numerical robustness would be a useful practical indicator.

major comments (2)
  1. [Abstract] Abstract (final paragraph): the inf-sup stability and a priori error estimate are asserted to follow from T-coercivity 'under several technical assumptions,' yet the manuscript provides no verification that these assumptions (on interface geometry, sign pattern, or contrast bounds) are satisfied by the high-contrast heterogeneous coefficients appearing in the numerical examples; this verification step is load-bearing for transferring the stability claim to the targeted regime.
  2. [Theory section] Theory development (T-coercivity application): the T-coercivity operator must map the trial space into the test space while controlling the contrast; without an explicit check or relaxation showing that the operator remains bounded independently of the contrast ratios used numerically, the stability result does not automatically extend to the high-contrast sign-changing setting the method is designed to address.
minor comments (1)
  1. [Numerical experiments] Numerical results are asserted to show robustness but lack tabulated error values, specific contrast ratios tested, or convergence tables that would allow quantitative assessment of the claimed effectiveness.

Simulated Author's Rebuttal

2 responses · 0 unresolved

Thank you for the careful review and constructive feedback. We address each major comment point by point below, indicating where revisions will be made to strengthen the connection between the theoretical assumptions and the numerical examples.

read point-by-point responses

  1. Referee: [Abstract] Abstract (final paragraph): the inf-sup stability and a priori error estimate are asserted to follow from T-coercivity 'under several technical assumptions,' yet the manuscript provides no verification that these assumptions (on interface geometry, sign pattern, or contrast bounds) are satisfied by the high-contrast heterogeneous coefficients appearing in the numerical examples; this verification step is load-bearing for transferring the stability claim to the targeted regime.

    Authors: We agree that an explicit link between the technical assumptions and the numerical test cases would improve clarity. The assumptions primarily concern the sign pattern (positive/negative subdomains with well-defined interfaces) and geometry, which are satisfied by construction in all reported examples. For contrast independence, the examples use ratios up to 10^6 that align with the regime where T-coercivity holds under the stated hypotheses. In the revision we will add a short paragraph after the abstract and a remark in Section 4 confirming that each numerical profile meets the interface and sign-pattern conditions, thereby justifying transfer of the stability result. revision: yes

  2. Referee: [Theory section] Theory development (T-coercivity application): the T-coercivity operator must map the trial space into the test space while controlling the contrast; without an explicit check or relaxation showing that the operator remains bounded independently of the contrast ratios used numerically, the stability result does not automatically extend to the high-contrast sign-changing setting the method is designed to address.

    Authors: The auxiliary-space construction in the CEM-GMsFEM is deliberately modified so that the T-coercivity operator T is chosen to satisfy ||T||_{V->V*} bounded by a constant independent of the contrast; this follows from the local corrector problems being posed with the sign-changing coefficient and the global constraint that enforces the inf-sup condition uniformly. The technical assumptions already encode the necessary control on the contrast through the T-coercivity constant. We will insert a clarifying sentence in the theory section (after the statement of Theorem 3.1) that explicitly notes the contrast-independent bound on T, together with a reference to the auxiliary-space design that achieves it. revision: partial

Circularity Check

0 steps flagged

No significant circularity; stability and error analysis invoke external T-coercivity under technical assumptions

full rationale

The paper adapts the existing CEM-GMsFEM framework to sign-changing coefficients and states that inf-sup stability plus a priori error estimates follow from T-coercivity theory once several technical assumptions hold. No equation or construction in the provided text reduces a claimed prediction or uniqueness result to a fitted parameter or self-citation by definition. The central theoretical claims are therefore not forced by the paper's own inputs; they rest on an external body of theory whose applicability is asserted rather than derived internally. This yields a low circularity score consistent with honest extension of prior non-circular work.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

Theoretical results rest on T-coercivity (standard in the field) plus unspecified technical assumptions required for the sign-changing case; no free parameters or new entities are mentioned.

axioms (2)
  • domain assumption T-coercivity theory applies to the sign-changing problem
    Invoked to establish inf-sup stability (abstract, final paragraph)
  • ad hoc to paper Several technical assumptions hold
    Required for both stability and a priori error estimate (abstract, final paragraph)

pith-pipeline@v0.9.0 · 5687 in / 1201 out tokens · 24618 ms · 2026-05-23T22:46:37.229177+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Multiscale Modeling for Time-harmonic Maxwell equations with impedance boundary conditions in highly heterogeneous media

    math.NA 2026-04 unverdicted novelty 7.0

    A multiscale method for Maxwell equations in high-contrast media achieves O(H) convergence independent of contrast by using an auxiliary spectral space that guarantees coercivity without divergence-free enforcement.

  2. Multiscale Methods for wave propagation in materials with sign-changing coefficients

    math.NA 2025-11 unverdicted novelty 5.0

    Adapts CEM-GMsFEM to sign-changing coefficients for EM waves, establishing inf-sup stability and a priori error estimates via T-coercivity and resolution conditions.

Reference graph

Works this paper leans on

52 extracted references · 52 canonical work pages · cited by 2 Pith papers

  1. [1]

    Abdulle, W

    A. Abdulle, W. E, B. Engquist, and E. Vanden-Eijnden , The heterogeneous multiscale method, Acta Numerica, 21 (2012), pp. 1–87, https://doi.org/10.1017/S0962492912000025

    work page doi:10.1017/s0962492912000025 2012

  2. [2]

    R. F. Almgern , An isotropic three dimensional structure with Poisson’s ratio = −1, Journal of Elasticity, 15 (1985), pp. 427–430, https://doi.org/10.1007/bf00042531

    work page doi:10.1007/bf00042531 1985

  3. [3]

    Altmann, P

    R. Altmann, P. Henning, and D. Peterseim , Numerical homogenization beyond scale separation, Acta Numerica, 30 (2021), pp. 1–86, https://doi.org/10.1017/S0962492921000015

    work page doi:10.1017/s0962492921000015 2021

  4. [4]

    Babuˇska and R

    I. Babuˇska and R. Lipton, Optimal local approximation spaces for generalized finite element methods with application to multiscale problems , Multiscale Modeling & Simulation. A SIAM Interdisciplinary Journal, 9 (2011), pp. 373–406, https://doi.org/10.1137/100791051

    work page doi:10.1137/100791051 2011

  5. [5]

    Babuˇska, R

    I. Babuˇska, R. Lipton, P. Sinz, and M. Stuebner , Multiscale-spectral GFEM and optimal over- sampling, Computer Methods in Applied Mechanics and Engineering, 364 (2020), pp. 112960, 28, https://doi.org/10.1016/j.cma.2020.112960

    work page doi:10.1016/j.cma.2020.112960 2020

  6. [6]

    Barr´e and P

    M. Barr´e and P. Ciarlet Jr. , The T-coercivity approach for mixed problems . In press in Comptes Rendus Mathematique, 2023, https://hal.science/hal-03820910

    work page 2023

  7. [7]

    Bensoussan, J.-L

    A. Bensoussan, J.-L. Lions, and G. Papanicolaou , Asymptotic analysis for periodic structures , AMS Chelsea Publishing, Providence, RI, 2011, https://doi.org/10.1090/chel/374. Corrected reprint of the 1978 original [MR0503330]

    work page doi:10.1090/chel/374 2011

  8. [8]

    Boffi, F

    D. Boffi, F. Brezzi, and M. Fortin , Mixed finite element methods and applications , Springer, Berlin, Germany and Heidelberg, Germany, 2013, https://doi.org/10.1007/978-3-642-36519-5 . 26

    work page doi:10.1007/978-3-642-36519-5 2013

  9. [9]

    Bonnet-BenDhia, C

    A.-S. Bonnet-BenDhia, C. Carvalho, and P. Ciarlet Jr. , Mesh requirements for the finite element approximation of problems with sign-changing coefficients, Numerische Mathematik, 138 (2018), pp. 801–838, https://doi.org/10.1007/s00211-017-0923-5

    work page doi:10.1007/s00211-017-0923-5 2018

  10. [10]

    Bonnet-BenDhia, L

    A.-S. Bonnet-BenDhia, L. Chesnel, and P. Ciarlet Jr. , T-coercivity for scalar interface prob- lems between dielectrics and metamaterials , ESAIM. Mathematical Modelling and Numerical Analysis, 46 (2012), pp. 1363–1387, https://doi.org/10.1051/m2an/2012006

    work page doi:10.1051/m2an/2012006 2012

  11. [11]

    Bonnet-BenDhia, L

    A.-S. Bonnet-BenDhia, L. Chesnel, and P. Ciarlet Jr. , T-coercivity for the Maxwell problem with sign-changing coefficients, Communications in Partial Differential Equations, 39 (2014), pp. 1007– 1031, https://doi.org/10.1080/03605302.2014.892128

    work page doi:10.1080/03605302.2014.892128 2014

  12. [12]

    Bonnet-BenDhia, L

    A.-S. Bonnet-BenDhia, L. Chesnel, and P. Ciarlet Jr. , Two-dimensional Maxwell’s equations with sign-changing coefficients, Applied Numerical Mathematics, 79 (2014), pp. 29–41, https://doi. org/10.1016/j.apnum.2013.04.006

    work page doi:10.1016/j.apnum.2013.04.006 2014

  13. [13]

    Bonnet-BenDhia, P

    A.-S. Bonnet-BenDhia, P. Ciarlet Jr., and C. M. Zw¨olf, Time harmonic wave diffraction prob- lems in materials with sign-shifting coefficients , Journal of Computational and Applied Mathematics, 234 (2010), pp. 1912–1919, https://doi.org/10.1016/j.cam.2009.08.041

    work page doi:10.1016/j.cam.2009.08.041 2010

  14. [14]

    Bunoiu, K

    R. Bunoiu, K. Karim, and C. Timofte , T-coercivity for the asymptotic analysis of scalar problems with sign-changing coefficients in thin periodic domains , Electronic Journal of Differential Equations, 2021 (2021), pp. 1–22, https://doi.org/10.58997/ejde.2021.59

    work page doi:10.58997/ejde.2021.59 2021

  15. [15]

    Bunoiu and K

    R. Bunoiu and K. Ramdani , Homogenization of materials with sign changing coefficients , Commu- nications in Mathematical Sciences, 14 (2016), pp. 1137–1154, https://doi.org/10.4310/cms.2016. v14.n4.a13

    work page doi:10.4310/cms.2016 2016

  16. [16]

    Bunoiu, K

    R. Bunoiu, K. Ramdani, and C. Timofte , T-coercivity for the homogenization of sign-changing coefficients scalar problems with extreme contrasts, Mathematical Reports, Romanian Academy of Sci- ences, 24(74) (2022), pp. 113–123

    work page 2022

  17. [17]

    Bunoiu, K

    R. Bunoiu, K. Ramdani, and C. Timofte , Homogenization of a transmission problem with sign- changing coefficients and interfacial flux jump , Communications in Mathematical Sciences, 21 (2023), pp. 2029–2049, https://doi.org/10.4310/cms.2023.v21.n7.a13

    work page doi:10.4310/cms.2023.v21.n7.a13 2023

  18. [18]

    Carvalho, L

    C. Carvalho, L. Chesnel, and P. Ciarlet Jr., Eigenvalue problems with sign-changing coefficients, Comptes Rendus Mathematique, 355 (2017), pp. 671–675, https://doi.org/10.1016/j.crma.2017. 05.002

    work page doi:10.1016/j.crma.2017 2017

  19. [19]

    Chaumont-Frelet and B

    T. Chaumont-Frelet and B. Verf¨urth, A generalized finite element method for problems with sign- changing coefficients, ESAIM. Mathematical Modelling and Numerical Analysis, 55 (2021), pp. 939–967, https://doi.org/10.1051/m2an/2021007

    work page doi:10.1051/m2an/2021007 2021

  20. [20]

    Chesnel and P

    L. Chesnel and P. Ciarlet Jr. , T-coercivity and continuous Galerkin methods: application to transmission problems with sign changing coefficients , Numerische Mathematik, 124 (2013), pp. 1–29, https://doi.org/10.1007/s00211-012-0510-8

    work page doi:10.1007/s00211-012-0510-8 2013

  21. [21]

    E. T. Chung and P. Ciarlet Jr. , A staggered discontinuous Galerkin method for wave propagation in media with dielectrics and meta-materials , Journal of Computational and Applied Mathematics, 239 (2013), pp. 189–207, https://doi.org/10.1016/j.cam.2012.09.033

    work page doi:10.1016/j.cam.2012.09.033 2013

  22. [22]

    E. T. Chung, Y. Efendiev, and W. T. Leung, Constraint energy minimizing generalized multiscale finite element method, Computer Methods in Applied Mechanics and Engineering, 339 (2018), pp. 298– 319, https://doi.org/10.1016/j.cma.2018.04.010. 27

    work page doi:10.1016/j.cma.2018.04.010 2018

  23. [23]

    E. T. Chung, Y. Efendiev, and G. Li , An adaptive GMsFEM for high-contrast flow problems , Journal of Computational Physics, 273 (2014), pp. 54–76, https://doi.org/10.1016/j.jcp.2014. 05.007

    work page doi:10.1016/j.jcp.2014 2014

  24. [24]

    Ciarlet Jr., T-coercivity: Application to the discretization of Helmholtz-like problems , Computers & Mathematics with Applications, 64 (2012), pp

    P. Ciarlet Jr., T-coercivity: Application to the discretization of Helmholtz-like problems , Computers & Mathematics with Applications, 64 (2012), pp. 22–34, https://doi.org/10.1016/j.camwa.2012. 02.034

    work page doi:10.1016/j.camwa.2012 2012

  25. [25]

    Ciarlet Jr

    P. Ciarlet Jr. and M. Vohral ´ık, Localization of global norms and robust a posteriori error con- trol for transmission problems with sign-changing coefficients , ESAIM. Mathematical Modelling and Numerical Analysis, 52 (2018), pp. 2037–2064, https://doi.org/10.1051/m2an/2018034

    work page doi:10.1051/m2an/2018034 2018

  26. [26]

    Cioranescu and P

    D. Cioranescu and P. Donato, An introduction to homogenization, vol. 17 of Oxford Lecture Series in Mathematics and its Applications, The Clarendon Press and Oxford University Press, New York, NY, 1999

    work page 1999

  27. [27]

    W. E and B. Engquist , The heterogeneous multiscale methods , Communications in Mathematical Sciences, 1 (2003), pp. 87–132

    work page 2003

  28. [28]

    Efendiev, J

    Y. Efendiev, J. Galvis, and T. Y. Hou, Generalized multiscale finite element methods (GMsFEM), Journal of Computational Physics, 251 (2013), pp. 116–135, https://doi.org/10.1016/j.jcp.2013. 04.045

    work page doi:10.1016/j.jcp.2013 2013

  29. [29]

    Efendiev and T

    Y. Efendiev and T. Y. Hou , Multiscale finite element methods: Theory and applications , vol. 4 of Surveys and Tutorials in the Applied Mathematical Sciences, Springer, New York, NY, 2009. Theory and applications

    work page 2009

  30. [30]

    Efendiev, T

    Y. Efendiev, T. Y. Hou, and X.-H. Wu, Convergence of a nonconforming multiscale finite element method, SIAM Journal on Numerical Analysis, 37 (2000), pp. 888–910, https://doi.org/10.1137/ s0036142997330329

    work page 2000

  31. [31]

    Gorb and Y

    Y. Gorb and Y. Kuznetsov , Asymptotic expansions for high-contrast scalar and vectorial PDEs , SIAM Journal on Applied Mathematics, 81 (2021), pp. 2246–2264, https://doi.org/10.1137/ 20m1357937

    work page 2021

  32. [32]

    T. Y. Hou and X.-H. Wu , A multiscale finite element method for elliptic problems in composite materials and porous media , Journal of Computational Physics, 134 (1997), pp. 169–189, https:// doi.org/10.1006/jcph.1997.5682

    work page doi:10.1006/jcph.1997.5682 1997

  33. [33]

    T. Y. Hou, X.-H. Wu, and Z. Cai , Convergence of a multiscale finite element method for elliptic problems with rapidly oscillating coefficients , Mathematics of Computation, 68 (1999), pp. 913–943, https://doi.org/10.1090/S0025-5718-99-01077-7

    work page doi:10.1090/s0025-5718-99-01077-7 1999

  34. [34]

    T. J. R. Hughes , Multiscale phenomena: Green’s functions, the Dirichlet-to-Neumann formulation, subgrid scale models, bubbles and the origins of stabilized methods , Computer Methods in Applied Mechanics and Engineering, 127 (1995), pp. 387–401, https://doi.org/10.1016/0045-7825(95) 00844-9

    work page doi:10.1016/0045-7825(95 1995

  35. [35]

    T. J. R. Hughes and G. Sangalli , Variational multiscale analysis: the fine-scale Green’s function, projection, optimization, localization, and stabilized methods , SIAM Journal on Numerical Analysis, 45 (2007), pp. 539–557, https://doi.org/10.1137/050645646

    work page doi:10.1137/050645646 2007

  36. [36]

    Ivrii, 100 years of Weyl’s law , Bulletin of Mathematical Sciences, 6 (2016), pp

    V. Ivrii, 100 years of Weyl’s law , Bulletin of Mathematical Sciences, 6 (2016), pp. 379–452, https: //doi.org/10.1007/s13373-016-0089-y . 28

    work page doi:10.1007/s13373-016-0089-y 2016

  37. [37]

    Kadic, T

    M. Kadic, T. B ¨uckmann, N. Stenger, M. Thiel, and M. Wegener , On the practicability of pentamode mechanical metamaterials , Applied Physics Letters, 100 (2012), https://doi.org/10. 1063/1.4709436

    work page 2012

  38. [38]

    R. S. Lakes, Negative-Poisson’s-ratio materials: Auxetic solids , Annual Review of Materials Research, 47 (2017), pp. 63–81, https://doi.org/10.1146/annurev-matsci-070616-124118

    work page doi:10.1146/annurev-matsci-070616-124118 2017

  39. [39]

    C. Ma, R. Scheichl, and T. Dodwell , Novel design and analysis of generalized finite element methods based on locally optimal spectral approximations , SIAM Journal on Numerical Analysis, 60 (2022), pp. 244–273, https://doi.org/10.1137/21m1406179

    work page doi:10.1137/21m1406179 2022

  40. [40]

    M˚alqvist and D

    A. M˚alqvist and D. Peterseim , Localization of elliptic multiscale problems , Mathematics of Com- putation, 83 (2014), pp. 2583–2603, https://doi.org/10.1090/S0025-5718-2014-02868-8

    work page doi:10.1090/s0025-5718-2014-02868-8 2014

  41. [41]

    M˚alqvist and D

    A. M˚alqvist and D. Peterseim , Numerical homogenization by localized orthogonal decomposition , vol. 5 of SIAM Spotlights, Society for Industrial & Applied Mathematics (SIAM), Philadelphia, PA, 2021

    work page 2021

  42. [42]

    G. W. Milton and A. V. Cherkaev, Which elasticity tensors are realizable?, Journal of Engineering Materials and Technology, 117 (1995), pp. 483–493, https://doi.org/10.1115/1.2804743

    work page doi:10.1115/1.2804743 1995

  43. [43]

    Ming and S

    P. Ming and S. Song, Error estimate of multiscale finite element method for periodic media revisited , Multiscale Modeling & Simulation. A SIAM Interdisciplinary Journal, 22 (2024), pp. 106–124, https: //doi.org/10.1137/22m1511060

    work page doi:10.1137/22m1511060 2024

  44. [44]

    Nicaise and J

    S. Nicaise and J. Venel , A posteriori error estimates for a finite element approximation of trans- mission problems with sign changing coefficients , Journal of Computational and Applied Mathematics, 235 (2011), pp. 4272–4282, https://doi.org/10.1016/j.cam.2011.03.028

    work page doi:10.1016/j.cam.2011.03.028 2011

  45. [45]

    L. A. Poveda, J. Galvis, and E. T. Chung, A second-order exponential integration constraint energy minimizing generalized multiscale method for parabolic problems, Journal of Computational Physics, 502 (2024), p. 112796, https://doi.org/10.1016/j.jcp.2024.112796

    work page doi:10.1016/j.jcp.2024.112796 2024

  46. [46]

    R. A. Shelby, D. R. Smith, and S. Schultz , Experimental verification of a negative index of refraction, Science, 292 (2001), pp. 77–79, https://doi.org/10.1126/science.1058847

    work page doi:10.1126/science.1058847 2001

  47. [47]

    Scale-Free Networks: Complex Webs in Nature and Technology

    R. Verf¨urth, A posteriori error estimation eechniques for finite element methods , Oxford University Press, Oxford, UK, Apr. 2013, https://doi.org/10.1093/acprof:oso/9780199679423.001.0001

    work page doi:10.1093/acprof:oso/9780199679423.001.0001 2013

  48. [48]

    V. G. Veselago, The electrodynamics of substances with simultaneously negative values of ϵ and µ, So- viet Physics Uspekhi, 10 (1968), pp. 509–514, https://doi.org/10.1070/pu1968v010n04abeh003699

    work page doi:10.1070/pu1968v010n04abeh003699 1968

  49. [49]

    Weyl, Ueber die asymptotische Verteilung der Eigenwerte , Nachrichten von der Gesellschaft der Wissenschaften zu G¨ ottingen, Mathematisch-Physikalische Klasse, 1911 (1911), pp

    H. Weyl, Ueber die asymptotische Verteilung der Eigenwerte , Nachrichten von der Gesellschaft der Wissenschaften zu G¨ ottingen, Mathematisch-Physikalische Klasse, 1911 (1911), pp. 110–117, http: //eudml.org/doc/58792

    work page 1911

  50. [50]

    Ye and E

    C. Ye and E. T. Chung , Constraint energy minimizing generalized multiscale finite element method for inhomogeneous boundary value problems with high contrast coefficients , Multiscale Modeling & Simulation. A SIAM Interdisciplinary Journal, 21 (2023), pp. 194–217, https://doi.org/10.1137/ 21m1459113

    work page 2023

  51. [51]

    C. Ye, H. Dong, and J. Cui , Convergence rate of multiscale finite element method for various boundary problems, Journal of Computational and Applied Mathematics, 374 (2020), p. 112754, https: //doi.org/10.1016/j.cam.2020.112754

    work page doi:10.1016/j.cam.2020.112754 2020

  52. [52]

    Zhao and E

    L. Zhao and E. T. Chung , An analysis of the NLMC upscaling method for high contrast problems , Journal of Computational and Applied Mathematics, 367 (2020), p. 112480, https://doi.org/10. 1016/j.cam.2019.112480. 29

    work page arXiv 2020