pith. sign in

arxiv: 1906.08453 · v1 · pith:Z5L2YQ4Gnew · submitted 2019-06-20 · 🧮 math.NA · cs.NA

Wavelet-based Edge Multiscale Finite Element Method for Helmholtz problems in perforated domains

Shubin Fu , Guanglian Li , Richard Craster , Sebastien Guenneau This is my paper

Pith reviewed 2026-05-25 19:46 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords Helmholtz equationperforated domainsmultiscale finite element methodwaveletsconvergence analysisphotonic crystals
0
0 comments X

The pith

The wavelet-based edge multiscale finite element method achieves O(H) convergence for Helmholtz problems in perforated domains on regular coarse meshes.

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

The paper develops an algorithm for Helmholtz equations in domains containing holes that remains accurate for large wavenumbers without requiring meshes finer than the wavelength. It extends an existing wavelet-based edge multiscale finite element framework to these perforated settings and proves that the error stays proportional to the coarse mesh size H. The proof requires a resolution assumption to hold and the wavelet level parameter to be large enough. This matters because standard finite element approaches become prohibitively expensive when waves oscillate rapidly inside complex geometries such as those appearing in photonic crystals.

Core claim

For Helmholtz problems in perforated domains, the WEMsFEM algorithm achieves an O(H) error bound on a regular coarse mesh of size H when the resolution assumption holds and the level parameter is large enough. The method is tested numerically in two dimensions on problems motivated by photonic crystals.

What carries the argument

The Wavelet-based Edge Multiscale Finite Element Method (WEMsFEM) adapted to Helmholtz problems in perforated domains, which builds multiscale basis functions from wavelets along coarse-mesh edges to capture fine-scale perforation effects.

If this is right

  • Solutions can be computed accurately on meshes whose size H is independent of the wavenumber.
  • The approach applies to wave propagation through materials containing periodic perforations.
  • Two-dimensional numerical experiments confirm the predicted convergence rate on photonic-crystal-type geometries.
  • The scheme remains efficient when the wavenumber becomes large.

Where Pith is reading between the lines

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

  • Verifying or relaxing the resolution assumption would widen the range of perforated structures the method can treat.
  • Similar wavelet constructions on edges could be tested for other time-harmonic wave equations in domains with holes.
  • Adapting the edge-based basis functions to faces would be needed for a three-dimensional version.

Load-bearing premise

The resolution assumption holds and the original WEMsFEM framework extends directly to Helmholtz problems in perforated domains.

What would settle it

A numerical test on a perforated domain with a known exact solution in which the observed error fails to decrease proportionally to H once the level parameter is increased and the resolution assumption is satisfied.

Figures

Figures reproduced from arXiv: 1906.08453 by Guanglian Li, Richard Craster, Sebastien Guenneau, Shubin Fu.

Figure 1
Figure 1. Figure 1: Perforated domains. Two models of finite locally periodic photonic crystals, used later when they [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Illustration of a coarse neighborhood and coarse element with an overlapping constant [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Reference solution and multiscale solution for model 1 with centered source, [PITH_FULL_IMAGE:figures/full_fig_p015_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Reference solution and multiscale solution for model 1 with right source, [PITH_FULL_IMAGE:figures/full_fig_p015_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: shows the L 2 (Ω ) and H1 (Ω )-relative errors versus `, and both the L 2 (Ω ) and H1 (Ω )- relative errors decay rapidly as more wavelet basis functions are added. For example, for the case that the source lies at the center shown in [PITH_FULL_IMAGE:figures/full_fig_p016_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Reference solution and multiscale solution for model 2 with centered source, [PITH_FULL_IMAGE:figures/full_fig_p016_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Reference solution and multiscale solution for model 2 with right source, [PITH_FULL_IMAGE:figures/full_fig_p017_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Relative error against level, H = 1/10. 5.2 Performance of Algorithm 1: coarse-scale mesh size H Earlier we established theoretically that the error induced by Algorithm 1 can attain O(H) in Proposition 4.1, upon the condition on the coarse mesh size H and the level parameter `, cf. (4.11). We now test how the algorithm performs with respect to a different, finer, coarse-scale mesh, we take its size H := 1… view at source ↗
Figure 9
Figure 9. Figure 9: Relative error against level, H = 1/20. 0 1 2 3 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 level Relative error L 2 error H 1 error (a) Error for model 2, center source. 0 1 2 3 0 0.2 0.4 0.6 0.8 1 1.2 1.4 level Relative error L 2 error H 1 error (b) Error for model 2, right source [PITH_FULL_IMAGE:figures/full_fig_p018_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Relative error against level, H = 1/20. 18 [PITH_FULL_IMAGE:figures/full_fig_p018_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Perforated domains, models of locally non-periodic crystals [PITH_FULL_IMAGE:figures/full_fig_p019_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: The L 2 (Ω )-relative error is 3.70%; further decreasing the coarse mesh size H or increasing the level parameter ` improves the accuracy as shown in [PITH_FULL_IMAGE:figures/full_fig_p019_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Relative error against level, H = 1/20. Next, we study the performance of Algorithm 1 in the perforated domain of model 4. The perforations in model 4 cross the neighborhood boundary and we depict the resulting four coarse neighborhood ωi in [PITH_FULL_IMAGE:figures/full_fig_p020_13.png] view at source ↗
Figure 15
Figure 15. Figure 15: Reference solution and multiscale solution for model 4 with centered source, [PITH_FULL_IMAGE:figures/full_fig_p021_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: Error for model 4, H = 1/10, center source. 6 Conclusion We demonstrate that the Wavelet-based Edge Multiscale Finite Element Method (WEMsFEM) for Helmholtz problems in perforated domains, with possibly large wavenumbers, are an effective alternative to standard Finite Element methods with advantages for multiscale problems. Such problems have many applications in photonic and phononic crystals and having… view at source ↗
read the original abstract

We introduce a new efficient algorithm for Helmholtz problems in perforated domains with the design of the scheme allowing for possibly large wavenumbers. Our method is based upon the Wavelet-based Edge Multiscale Finite Element Method (WEMsFEM) as proposed recently in [14]. For a regular coarse mesh with mesh size H, we establish O(H) convergence of this algorithm under the resolution assumption, and with the level parameter being sufficiently large. The performance of the algorithm is demonstrated by extensive 2-dimensional numerical tests including those motivated by photonic crystals.

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 introduces a Wavelet-based Edge Multiscale Finite Element Method (WEMsFEM) for Helmholtz problems in perforated domains, extending the framework of reference [14]. For a regular coarse mesh of size H the central claim is an O(H) convergence rate under an unspecified resolution assumption together with a sufficiently large level parameter. The method is intended to accommodate possibly large wavenumbers, and its performance is illustrated by 2-D numerical experiments including photonic-crystal examples.

Significance. If the O(H) rate can be shown to hold with constants independent of the wavenumber k, the result would be a useful contribution to multiscale methods for high-frequency wave propagation in complex perforated geometries. Such problems arise in photonics and acoustics, where standard FEMs suffer from pollution and fine-mesh requirements; a coarse-mesh method with controlled error would therefore be of practical interest.

major comments (2)
  1. [Abstract] Abstract: the O(H) convergence statement is conditioned on an unspecified 'resolution assumption.' Because the central claim rests on this assumption, its precise formulation (including any implicit dependence on the wavenumber k) must be stated explicitly; without it the reader cannot verify whether the estimate controls the pollution terms that appear in standard Helmholtz a-priori bounds.
  2. [Analysis section (convergence theorem)] Analysis section (convergence theorem): the direct transfer of the WEMsFEM construction and error analysis from the coercive elliptic setting of [14] to the indefinite Helmholtz operator must be justified. Any hidden constants that grow with k would invalidate the advertised O(H) rate on a fixed coarse mesh; the proof must therefore exhibit the k-dependence (or k-independence) of all constants.
minor comments (1)
  1. [Abstract] The phrase 'level parameter being sufficiently large' is used without a quantitative bound; a concrete condition on this parameter in terms of H and k would improve readability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments on our manuscript. We address each major comment below and will revise the paper to improve clarity on the resolution assumption and the k-dependence in the analysis.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the O(H) convergence statement is conditioned on an unspecified 'resolution assumption.' Because the central claim rests on this assumption, its precise formulation (including any implicit dependence on the wavenumber k) must be stated explicitly; without it the reader cannot verify whether the estimate controls the pollution terms that appear in standard Helmholtz a-priori bounds.

    Authors: We agree that the resolution assumption should be stated explicitly in the abstract. The assumption (detailed in Assumption 3.1) requires that the coarse mesh size H is sufficiently small relative to the wavelength 1/k to ensure the error bound holds without pollution dominating. In the revised version we will add a concise formulation of this assumption to the abstract, including its dependence on k. revision: yes

  2. Referee: [Analysis section (convergence theorem)] Analysis section (convergence theorem): the direct transfer of the WEMsFEM construction and error analysis from the coercive elliptic setting of [14] to the indefinite Helmholtz operator must be justified. Any hidden constants that grow with k would invalidate the advertised O(H) rate on a fixed coarse mesh; the proof must therefore exhibit the k-dependence (or k-independence) of all constants.

    Authors: We acknowledge that the transfer requires explicit justification. The proof of Theorem 4.1 adapts the arguments from [14] using a Gårding-type inequality for the Helmholtz operator and invokes the resolution assumption to absorb the indefinite terms and control pollution. Under this assumption the constants in the O(H) estimate are independent of k. To make this transparent we will add a remark in the analysis section that explicitly tracks the k-dependence and confirms independence when the resolution assumption holds. revision: yes

Circularity Check

0 steps flagged

No circularity; convergence established by analysis under explicit assumptions

full rationale

The paper extends the WEMsFEM method from reference [14] and proves an O(H) convergence rate for the Helmholtz problem under a stated resolution assumption plus sufficiently large level parameter. The derivation chain consists of standard multiscale finite-element analysis applied to the perforated-domain Helmholtz operator; no equation reduces by construction to a fitted parameter, self-referential definition, or load-bearing self-citation whose validity is presupposed rather than re-proved. The result remains falsifiable by direct numerical counter-example on the coarse mesh and is therefore self-contained.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on an unspecified resolution assumption and on the prior WEMsFEM construction in reference [14]. No free parameters or invented entities are mentioned in the abstract.

axioms (1)
  • domain assumption Resolution assumption
    Invoked as a prerequisite for the O(H) convergence result.

pith-pipeline@v0.9.0 · 5620 in / 1089 out tokens · 23291 ms · 2026-05-25T19:46:21.650281+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

37 extracted references · 37 canonical work pages · 1 internal anchor

  1. [1]

    Berenger

    J.-P. Berenger. A perfectly matched layer for the absorption of electromagnetic waves. Journal of computational physics, 114(2):185–200, 1994

    work page 1994

  2. [2]

    Berlyand and H

    L. Berlyand and H. Owhadi. Flux norm approach to finite dimensional homogenization approximations with non-separated scales and high contrast. Arch. Ration. Mech. Anal., 198(2):677–721, 2010

    work page 2010

  3. [3]

    D. N. Chigrin, S. Enoch, C. M. S. Torres, and G. Tayeb. Self-guiding in two-dimensional photonic crystals. Optics Express, 11:1203–1211, 2003

    work page 2003

  4. [4]

    Collino and C

    F. Collino and C. Tsogka. Application of the perfectly matched absorbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media. Geophysics, 66(1):294–307, 2001

    work page 2001

  5. [5]

    www.comsol.com, 2012

    COMSOL. www.comsol.com, 2012

    work page 2012

  6. [6]

    R. V . Craster and S. Guenneau, editors. Acoustic Metamaterials. Springer-Verlag, 2012

    work page 2012

  7. [7]

    R. V . Craster, J. Kaplunov, and A. V . Pichugin. High frequency homogenization for periodic media. Proc R Soc Lond A, 466:2341–2362, 2010

    work page 2010

  8. [8]

    T. J. Cui, D. Smith, and R. Liu, editors. Metamaterials. Springer, Boston MA, 2010. 24

    work page 2010

  9. [9]

    Dubois, J

    M. Dubois, J. Perchoux, A. L. Vanel, C. Tronche, Y . Achaoui, G. Dupont, K. Bertling, A. D. Raki ´c, T. Antonakakis, S. Enoch, R. Abdeddaim, R. V . Craster, and S. Guenneau. Acoustic flat lensing using an indefinite medium. Phys. Rev. B, 99:100301(R), 2019

    work page 2019

  10. [10]

    W. E and B. Engquist. The heterogeneous multiscale methods. Commun. Math. Sci. , 1(1):87–132, 2003

    work page 2003

  11. [11]

    Efendiev, J

    Y . Efendiev, J. Galvis, and T. Hou. Generalized multiscale finite element methods. J. Comput. Phys., 251:116–135, 2013

    work page 2013

  12. [12]

    Efendiev, J

    Y . Efendiev, J. Galvis, and X.-H. Wu. Multiscale finite element methods for high-contrast problems using local spectral basis functions. J. Comput. Phys., 230(4):937–955, 2011

    work page 2011

  13. [13]

    Engquist and L

    B. Engquist and L. Ying. Sweeping preconditioner for the helmholtz equation: hierarchical matrix representation. Communications on pure and applied mathematics, 64(5):697–735, 2011

    work page 2011

  14. [14]

    S. Fu, E. Chung, and G. Li. Edge multiscale methods for elliptic problems with heterogeneous coeffi- cients. arXiv preprint arXiv:1810.10398, 2018

    work page arXiv 2018

  15. [15]

    Hiptmair, A

    R. Hiptmair, A. Moiola, and I. Perugia. A survey of Trefftz methods for the Helmholtz equation. In Building bridges: connections and challenges in modern approaches to numerical partial differential equations, volume 114 of Lect. Notes Comput. Sci. Eng., pages 237–278. Springer, [Cham], 2016

    work page 2016

  16. [16]

    Hou and X.-H

    T. Hou and X.-H. Wu. A multiscale finite element method for elliptic problems in composite materials and porous media. J. Comput. Phys., 134(1):169–189, 1997

    work page 1997

  17. [17]

    Hughes, G

    T. Hughes, G. Feij ´oo, L. Mazzei, and J.-B. Quincy. The variational multiscale method—a paradigm for computational mechanics. Comput. Methods Appl. Mech. Engrg., 166(1-2):3–24, 1998

    work page 1998

  18. [18]

    J. D. Joannopoulos, S. G. Johnson, R. D. Meade, and J. N. Winn. Photonic Crystals (2nd Edition) . Princeton University Press, Princeton, 2008

    work page 2008

  19. [19]

    S. G. Johnson and J. D. Joannopoulos. Block-iterative frequency-domain methods for Maxwell’s equations in a planewave basis. Optics Express, 8:173–190, 2001

    work page 2001

  20. [20]

    Lavrinenko, P

    A. Lavrinenko, P. I. Borel, L. H. Frandsen, M. Thorhauge, A. Harpøth, M. Kristensen, T. Niemi, and H. M. H. Chong. Comprehensive FDTD modelling of photonic crystal waveguide components.Optics Express, 12:234–248, 1 2004

    work page 2004

  21. [21]

    G. Li. On the convergence rates of GMsFEMs for heterogeneous elliptic problems without oversam- pling techniques. https://arxiv.org/abs/1802.08873, 2018

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  22. [22]

    G. Li, D. Peterseim, and M. Schedensack. Error analysis of a variational multiscale stabilization for convection-dominated diffusion equations in two dimensions. IMA J. Numer. Anal., 38(3):1229–1253, 2018

    work page 2018

  23. [23]

    L. Lu, J. D. Joannopoulos, and M. Soljacic. Topological photonics. Nature Photonics, 8:821–829, 2014

    work page 2014

  24. [24]

    www.lumerical.com, 2018

    Lumerical. www.lumerical.com, 2018. 25

    work page 2018

  25. [25]

    S. A. Maier, editor. World Scientific Handbook of Metamaterials and Plasmonics . World Scientific Series in Nanoscience and Nanotechnology, 2017

    work page 2017

  26. [26]

    Maling, D

    B. Maling, D. J. Colquitt, and R. V . Craster. Dynamic homogenisation of Maxwell’s equations with applications to photonic crystals and localised waveforms on gratings. Wave Motion, 69:35–49, 2017

    work page 2017

  27. [27]

    M ˚alqvist and D

    A. M ˚alqvist and D. Peterseim. Localization of elliptic multiscale problems. Math. Comp. , 83(290):2583–2603, 2014

    work page 2014

  28. [28]

    Melenk and I

    J. Melenk and I. Babu ˇska. The partition of unity finite element method: basic theory and applications. Comput. Methods Appl. Mech. Engrg., 139(1-4):289–314, 1996

    work page 1996

  29. [29]

    J. M. Melenk. On generalized finite element methods. PhD thesis, 1995

    work page 1995

  30. [30]

    J. M. Melenk, A. Parsania, and S. Sauter. General DG-methods for highly indefinite Helmholtz prob- lems. J. Sci. Comput., 57(3):536–581, 2013

    work page 2013

  31. [31]

    J. M. Melenk and S. Sauter. Wavenumber explicit convergence analysis for Galerkin discretizations of the Helmholtz equation. SIAM J. Numer. Anal., 49(3):1210–1243, 2011

    work page 2011

  32. [32]

    Peterseim

    D. Peterseim. Eliminating the pollution effect in Helmholtz problems by local subscale correction. Mathematics of Computation, 86(305):1005–1036, 2017

    work page 2017

  33. [33]

    Peterseim and B

    D. Peterseim and B. Verf ¨urth. Computational high frequency scattering from high contrast heteroge- neous media. arXiv preprint arXiv:1902.09935, 2019

    work page arXiv 1902

  34. [34]

    D. W. Prather, S. Shi, J. Murakowski, G. J. Schneider, A. Sharkawy, C. Chen, B. Miao, and R. Martin. Self-collimation in photonic crystal structures: a new paradigm for applications and device develop- ment. Journal of Physics D: Applied Physics, 40(9):2635, 2007

    work page 2007

  35. [35]

    S. A. Ramakrishna and T. Grzegorczyk. Physics and Applications of Negative Refractive Index Mate- rials. CRC Press and SPIE Press, 2008

    work page 2008

  36. [36]

    K. C. T. Strouboulis, I. Babu ˇska. The design and analysis of the generalized finite element method. Comput. Methods Appl. Mech. Engrg., 181:43–69, 2000

    work page 2000

  37. [37]

    Zolla, G

    F. Zolla, G. Renversez, A. Nicolet, B. Kuhlmey, S. Guenneau, and D. Felbacq.Foundations of photonic crystal fibres (2nd Edition). Imperial College Press, London, 2012. 26

    work page 2012