Abstract
This paper is concerned with the design of two different classes of Galerkin boundary element methods for the solution of high-frequency sound-hard scattering problems in the exterior of two-dimensional smooth convex scatterers. We prove in this paper that both methods require a small increase (in the order of \(k^\epsilon \) for any \(\epsilon > 0\)) in the number of degrees of freedom to guarantee frequency independent precisions with increasing wavenumber k. In addition, the accuracy of the numerical solutions are independent of frequency provided sufficiently many terms in the asymptotic expansion are incorporated into the integral equation formulation. Numerical results validating \(\mathcal {O}(k^\epsilon )\) algorithms are presented.









Similar content being viewed by others
References
Abramowitz, M., Stegun, I.A.: Handbook of mathematical functions with formulas, graphs, and mathematical tables, National Bureau of Standards Applied Mathematics Series, vol. 55. U.S. Government Printing Office, Washington, D.C, For sale by the Superintendent of Documents (1964)
Amini, S., Profit, A.: Multi-level fast multipole solution of the scattering problem. Eng. Anal. Bound. Elem. 27(5), 547–564 (2003)
Anand, A., Boubendir, Y., Ecevit, F., Reitich, F.: Analysis of multiple scattering iterations for high-frequency scattering problems. II. The three-dimensional scalar case. Numer. Math. 114(3), 373–427 (2010)
Banjai, L., Hackbusch, W.: Hierarchical matrix techniques for low- and high-frequency Helmholtz problems. IMA J. Numer. Anal. 28(1), 46–79 (2008)
Boffi, D.: Finite element approximation of eigenvalue problems. Acta Numer. 19, 1–120 (2010)
Boubendir, Y., Turc, C.: Wave-number estimates for regularized combined field boundary integral operators in acoustic scattering problems with Neumann boundary conditions. IMA J. Numer. Anal. 33(4), 1176–1225 (2013)
Bruno, O., Geuzaine, C., Reitich, F.: On the \(O(1)\) solution of multiple-scattering problems. IEEE Trans. Magn. 41(5), 1488–1491 (2005)
Bruno, O.P., Domínguez, V., Sayas, F.J.: Convergence analysis of a high-order Nyström integral-equation method for surface scattering problems. Numer. Math. 124(4), 603–645 (2013)
Bruno, O.P., Geuzaine, C.A.: An \(O(1)\) integration scheme for three-dimensional surface scattering problems. J. Comput. Appl. Math. 204(2), 463–476 (2007)
Bruno, O..P., Geuzaine, C..A., Monro, J..A., Jr., Reitich, F.: Prescribed error tolerances within fixed computational times for scattering problems of arbitrarily high frequency: the convex case. Philos. Trans. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci. 362(1816), 629–645 (2004)
Bruno, O.P., Kunyansky, L.A.: A fast, high-order algorithm for the solution of surface scattering problems: basic implementation, tests, and applications. J. Comput. Phys. 169(1), 80–110 (2001)
Céa, J.: Approximation variationnelle des problèmes aux limites. Ann. Inst. Fourier (Grenoble) 14(fasc. 2), 345–444 (1964)
Chandler-Wilde, S.N., Graham, I.G.: Boundary integral methods in high frequency scattering. In: Highly oscillatory problems, London Math. Soc. Lecture Note Ser., vol. 366, pp. 154–193. Cambridge University Press, Cambridge (2009)
Chandler-Wilde, S.N., Graham, I.G., Langdon, S., Spence, E.A.: Numerical-asymptotic boundary integral methods in high-frequency acoustic scattering. Acta Numer. 21, 89–305 (2012)
Chandler-Wilde, S.N., Hewett, D.P.: Wavenumber-explicit continuity and coercivity estimates in acoustic scattering by planar screens. Integr. Equa. Oper. Theory 82(3), 423–449 (2015)
Chandler-Wilde, S.N., Hewett, D.P., Langdon, S., Twigger, A.: A high frequency boundary element method for scattering by a class of nonconvex obstacles. Numer. Math. 129(4), 647–689 (2015)
Chandler-Wilde, S.N., Langdon, S.: A Galerkin boundary element method for high frequency scattering by convex polygons. SIAM J. Numer. Anal. 45(2), 610–640 (2007)
Chandler-Wilde, S.N., Langdon, S., Mokgolele, M.: A high frequency boundary element method for scattering by convex polygons with impedance boundary conditions. Commun. Comput. Phys. 11(2), 573–593 (2012)
Chandler-Wilde, S.N., Spence, E.A., Gibbs, A., Smyshlyaev, V.P.: High-frequency bounds for the Helmholtz equation under parabolic trapping and applications in numerical analysis. SIAM J. Math. Anal. 52(1), 845–893 (2020)
Colton, D., Kress, R.: Inverse Acoustic and Electromagnetic Scattering Theory, Applied Mathematical Sciences, vol. 93. Springer, Berlin (1992)
Colton, D.L., Kress, R.: Integral Equation Methods in Scattering Theory. Pure and Applied Mathematics (New York). Wiley, New York (1983)
Davies, R.W., Morgan, K., Hassan, O.: A high order hybrid finite element method applied to the solution of electromagnetic wave scattering problems in the time domain. Comput. Mech. 44(3), 321–331 (2009)
Domínguez, V., Graham, I.G., Smyshlyaev, V.P.: A hybrid numerical-asymptotic boundary integral method for high-frequency acoustic scattering. Numer. Math. 106(3), 471–510 (2007)
Ecevit, F.: Frequency independent solvability of surface scattering problems. Turkish J. Math. 42(2), 407–422 (2018)
Ecevit, F., Eruslu, H.H.: A Galerkin BEM for high-frequency scattering problems based on frequency-dependent changes of variables. IMA J. Numer. Anal. 39(2), 893–923 (2019)
Ecevit, F., Özen, H.Ç.: Frequency-adapted Galerkin boundary element methods for convex scattering problems. Numer. Math. 135(1), 27–71 (2017)
Ecevit, F., Reitich, F.: Analysis of multiple scattering iterations for high-frequency scattering problems. I. The two-dimensional case. Numer. Math. 114(2), 271–354 (2009)
Fedoryuk, M.V.: The stationary phase method and pseudodifferential operators. Russ. Math. Surv. 26(1), 65 (1971)
Galkowski, J.: Distribution of resonances in scattering by thin barriers. Mem. Am. Math. Soc. 259(1248), ix+152 (2019)
Galkowski, J., Müller, E.H., Spence, E.A.: Wavenumber-explicit analysis for the Helmholtz \(h\)-BEM: error estimates and iteration counts for the Dirichlet problem. Numer. Math. 142(2), 329–357 (2019)
Giladi, E., Keller, J.B.: An asymptotically derived boundary element method for the helmholtz equations. Appl. Comput. Electromagn. Soc. (2004)
Groth, S., Hewett, D., Langdon, S.: A hybrid numerical-asymptotic boundary element method for high frequency scattering by penetrable convex polygons. Wave Motion 78, 32–53 (2018)
Groth, S.P., Hewett, D.P., Langdon, S.: Hybrid numerical-asymptotic approximation for high-frequency scattering by penetrable convex polygons. IMA J. Appl. Math. 80(2), 324–353 (2013)
Hesthaven, J., Warburton, T.: High-order accurate methods for time-domain electromagnetics. CMES Comput. Model. Eng. Sci. 5(5), 395–407 (2004)
Hewett, D.P.: Shadow boundary effects in hybrid numerical-asymptotic methods for high-frequency scattering. Eur. J. Appl. Math. 26(5), 773–793 (2015)
Hewett, D.P., Langdon, S., Chandler-Wilde, S.N.: A frequency-independent boundary element method for scattering by two-dimensional screens and apertures. IMA J. Numer. Anal. 35(4), 1698–1728 (2015)
Hewett, D.P., Langdon, S., Melenk, J.M.: A high frequency \$hp\$ boundary element method for scattering by convex polygons. SIAM J. Numer. Anal. 51(1), 629–653 (2013)
Huybrechs, D., Vandewalle, S.: A sparse discretization for integral equation formulations of high frequency scattering problems. SIAM J. Sci. Comput. 29(6), 2305–2328 (2007)
Johnson, W.P.: The curious history of Faà di Bruno’s formula. Am. Math. Monthly 109(3), 217–234 (2002)
Langdon, S., Mokgolele, M., Chandler-Wilde, S.: High frequency scattering by convex curvilinear polygons. J. Comput. Appl. Math. 234(6), 2020–2026 (2010)
Lazergui, S., Boubendir, Y.: Asymptotic expansions of the Helmholtz equation solutions using approximations of the Dirichlet to Neumann operator. J. Math. Anal. Appl. 456(2), 767–786 (2017)
Martinez, A.: An Introduction to Semiclassical and Microlocal Analysis. Springer, New York (2002)
Melrose, R.B.: Local Fourier-Airy integral operators. Duke Math. J. 42(4), 583–604 (1975)
Melrose, R.B.: Airy operators. Commun. Partial. Differ. Equ. 3(1), 1–76 (1978)
Melrose, R.B., Taylor, M.E.: Near peak scattering and the corrected Kirchhoff approximation for a convex obstacle. Adv. Math. 55(3), 242–315 (1985)
Namburu, R., Mark, E., Clarke, J.: Scalable electromagnetic simulation environment. CMES Comput. Model. Eng. Sci. 5(5), 443–453 (2004)
Schwab, C.: \(p\)- and \(hp\)-Finite Element Methods. Numerical Mathematics and Scientific Computation. The Clarendon Press, Oxford University Press, New York (1998)
Sherer, S.E., Visbal, M.R.: Time-domain scattering simulations using a high-order overset-grid approach. In: Workshop on Computational Electromagnetics in Time-Domain, 2005. CEM-TD 2005., pp. 44–47 (2005)
Taflove, A., Hagness, S.C.: Computational electrodynamics: The finite-difference time-domain method, 2nd edn. Artech House Inc, Boston (2000)
Tong, M.S., Chew, W.C.: Multilevel fast multipole acceleration in the Nyström discretization of surface electromagnetic integral equations for composite objects. IEEE Trans. Antennas Propagat. 58(10), 3411–3416 (2010)
Acknowledgements
We thank the anonymous referees for their constructive comments as this helped in improving the quality of the paper.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The authors declare that they have no conflict of interest.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
F. Ecevit is supported by the Scientific and Technological Research Council of Turkey through grant TÜBİTAK-1001-117F056. Y. Boubendir’s work was supported by the NSF through Grants DMS-1720014 and DMS-2011843. A. Anand gratefully acknowledges the support by Science & Engineering Research Board through File No MTR/2017/000643.
Appendices
Appendix
A Asymptotic expansion of the total field
The ansatz representing the asymptotic expansion of the total field \(\eta \) for the Neumann boundary value problem used in the present paper is given by Melrose and Taylor in [45]. However, the authors did not present all the mathematical steps needed in its derivation. In this section, we provide the missing details of this analysis.
Let \(K \subset {\mathbb {R}}^{n+1}\) be a compact strictly convex obstacle such that \(B= \partial K \subset {\mathbb {R}}^{n+1}\) is a smooth hyper-surface, and consider the Neumann-to-Dirichlet operator
where \({\mathcal {D}}'(B \times {\mathbb {R}})\) and \({\mathcal {E}}'({{\mathbb {R}}\times B})\) are the spaces of distributions and compactly supported distributions respectively, and v is the solution to the wave problem
wherein \(\varOmega = {\mathbb {R}}^{n+1} \setminus K\) is the exterior domain, and \(\nu \) is the outward unit normal.
In what follows, for an incident field \(v^{i}(t,x) = \delta (t - \alpha \cdot x)\) (\(\delta \) is the Dirac function) with direction \(\alpha \in {\mathcal {S}}^n\), we denote the solution of the wave problem (58) associated with \(f(t,x) = (\alpha \, \cdot \, \nu (x)) \partial _t v^{i}(t,x)\) by v. In this case, the total field \(v^t := v + v^{i}\) can be expressed on the boundary as [45, p.296]
where I is the identity operator. Using the same notation and procedure in [45], we define the Kirchhoff operator (see [45, Equation 8.26])
where F is the Fourier integral operator [41, Equation 9]
with kernel \(\kappa _{F}(t,\alpha ,x):=\delta (t-\alpha \cdot x)\).
As shown in [45], the asymptotic behavior of the total field \(\eta \) is determined by the kernel \(\kappa _{Q_N}\) of the Kirchhoff operator \(Q_N\).
Lemma 6
[45, Lemma 9.1] The asymptotic behavior as \(k \rightarrow \pm \infty \) of the total field \(\eta (\alpha ,k,x)\) obtained by inverse Fourier transformation of the kernel \(\kappa _{Q_N}(\alpha ,t,x)\)
is determined by the singularities of the kernel \(\kappa _{Q_N}\) modulo rapidly decreasing terms.
The main goal is therefore the study of the kernel \(\kappa _{Q_N}\). To this end, one first utilizes the general theory of Fourier integral operators with folding canonical relations [45, §5 and §6] to decompose the operator \(Q_N\). In what follows, we use the notation in [45] except that \(OPS^\theta _{\xi ,\zeta }\) denotes the collection of pseudodifferential operators with symbols in the Hörmander class \(S^\theta _{\xi ,\zeta }\), and we write \(OPS^\theta \) for \(OPS^\theta _{1,0}\).
Theorem 6
[45, Theorem 8.30] The Kirchhoff operator \(Q_N\) can be expressed as
where \(J_1\) and \(J_2\) are elliptic Fourier integral operators of order zero, \(D\in OPS^{-\frac{n}{2}-\frac{1}{6}}_{\frac{1}{3},0}\) has an asymptotic expansion
with \(A_{r,\ell } \in S^{-\frac{n}{2}-\frac{\ell }{3}-\frac{r}{3}+(\ell +1)_-}_{cl}\) and \(\varPhi ^{r,\ell }(k^{-\frac{1}{3}}\xi _1)\in S^{(\frac{\ell }{3}-\frac{2r}{3})_{+}}_{\frac{1}{3},0} ({\mathbb {R}})\) (see [44, p.11-12]) so that
where \((k,\xi )\) are variables dual to (t, x), and \(a_{r,\ell } \in S^{-\frac{n}{2}-\frac{r}{3}-\frac{\ell }{3}+(\ell +1)_-}_{1,0}\) admits an asymptotic expansion
wherein \(a_{q,r,\ell }\) are \(C^{\infty }\) functions uniformly bounded together with all their derivatives (cf. [42, Definition 2.5.6]), and \({\mathcal {A}}^{-1}\) is the convolution operator defined by Fourier transformation as [45, Equation 1.36]
where \(A_+(z) := Ai(e^{\frac{2\pi i}{3}}z)\) and Ai is the Airy function [43].
The Fourier integral operator \(J_2\) is given by
where \(a_{J_2} \in S^0_{1,0}\) (\(a_{J_2}\) does not depend on \(\xi \) and k because it is a symbol of order 0). Applying the Dirac function at the base point \((0,{\bar{\alpha }})\) (see [41, 45]) yields
where the operator \(P\in OPS^{-\frac{n}{2}+\frac{1}{6}}\) is specified by
Accordingly, use of (66) in (62) implies
Finally, since \(P\delta _0 = P^{\#}\delta _0 \mod OPS^{-\infty }\) and \(P^{\#}\) commutes with \({\mathcal {A}}^{-1}\) [45, p.295], (67) can be rewritten as
modulo rapidly decreasing terms.
In what follows, we briefly explain how representation (68) can be used to express the amplitude associated with the kernel \(\kappa _N\) as an asymptotic series of oscillatory integrals each of which is amenable to an application of the stationary phase method [28].
To this end, we first use (63) to deduce for the composition \(DP^{\#}\) of the pseudo-differential operators D and \(P^{\#}\)
with
where \(p^{\#} \in S^0_{1,0}\), and \((\tau ,\eta )\) is the dual variable to \((t',y)\). Using (70) in (69), we therefore get
with
As for the composition of the operator \(DP^{\#}\) with the Fourier integral operator \(J_1\) appearing in (62), let us first note that
where \((\tau ,\eta )\) are variables dual to (t, x), \(a_{J_1} \in S^0_{1,0}\), and the phase function \(\psi _1\) is defined in a neighborhood of the base point as [45, Equations 7.11 and 7.13]
In (73), we have used the notation \(x'=(x_2,\ldots ,x_n)\) for \(x = (x_1,\ldots ,x_n) \in B\) and similarly for \(\eta \in {\mathbb {R}}^n\). Combining (71) with (72), we obtain
and we rewrite this as
where
In light of (68), substituting [45, p.295]
for v in (74), we get
where \(q_{N_{r,\ell }}({\bar{\alpha }},x,k,\xi ):= a_J\# b_{r,\ell }({\bar{\alpha }},x,k,\xi )\).
For the kernel \(\kappa _{Q_N}\), we accordingly have
and, in virtue of Lemma 6, we deduce (see [45, p.298])
Making the change of variable \(\xi =k \zeta \) in (77), we therefore find [45, p.298]
where \(\psi _2(\alpha , x,\zeta ) = \psi _1(x,1,\zeta )-{\bar{\alpha }} \cdot \zeta \), and with a slight abuse of notation we have written \(q_{N_{r,\ell }}({\bar{\alpha }},x,k,\zeta )\) for \(k^n q_{N_{r,\ell }}({\bar{\alpha }},x,k,\zeta )\). In this case, \(q_{N_{r,\ell }}\) is a symbol of order \(\frac{n}{2}+\frac{1}{6}-\frac{r}{3}-\frac{\ell }{3}+(\ell +1)_{-}\) so that
Using (79) in (78), we arrive at
In order to further simplify (80), we introduce the function
The symbolic behavior of \(\varPsi ^{r,\ell }\) is as follows.
Lemma 7
[45, Lemma 9.34] The function \(\varPsi ^{r,\ell }\) defined in (81) belongs to \(S^{1+\ell -2r}({\mathbb {R}})\), admits an asymptotic expansion
as \(\tau \rightarrow +\infty \), and is rapidly decreasing in the sense of Schwarz as \(\tau \rightarrow -\infty \).
Rewriting (81) as
and using \({\mathcal {F}}^{-1}\) to denote the inverse Fourier transform, we obtain
Using (83) in (80), we finally conclude
where
with the phase function given by
As shown in [45], the integrals \(I_{q,r,\ell }(\alpha ,k,x)\) can be treated using the stationary phase method which results in
and this leads into the following for the envelope \(\eta ^{\mathrm{slow}}(\alpha ,k,x) = e^{-ik \alpha \cdot x} \eta (\alpha ,k,x)\).
Theorem 7
[45, Theorem 9.36] In a vicinity of the shadow boundary \(\{ x \in B : \alpha \cdot \nu (x) = 0 \}\), \(\eta ^{\mathrm{slow}}(\alpha ,k,x)\) belongs to the Hörmander class \(S^{0}_{\frac{2}{3},\frac{1}{3}}\) and admits an asymptotic expansion
with
where \(b_{p,q,r,\ell }\) are complex-valued \(C^{\infty }\) functions, Z is a real-valued \(C^{\infty }\) function that is positive on the illuminated region \(\{ x \in B : \alpha \cdot \nu (x) < 0 \}\), negative on the shadow region \(\{ x \in B : \alpha \cdot \nu (x) >0 \}\), and that vanishes precisely to first order at the shadow boundary.
Under certain assumptions, Theorem 7 is in fact valid over the entire boundary B. This is given in the next theorem where we use the notation \(B^{\epsilon }_{\gtrless } = \{ x \in B: \alpha \cdot \nu (x) \gtrless \epsilon \}\).
Theorem 8
Assume there exists \(\epsilon \in (0,1)\) such that on \(B^{\epsilon }_{<}\) the envelope \(\eta ^{\mathrm{slow}}\) belongs to \(S^0_{1,0}(B^{\epsilon }_{<} \times (0,\infty ))\) and admits an asymptotic expansion
and it is rapidly decreasing in the sense of Schwarz on \(B^{\epsilon }_{>}\) as \(k \rightarrow -\infty \). Then \(\eta ^{\mathrm{slow}} \in S^{0}_{\frac{2}{3},\frac{1}{3}}(B \times (0,\infty ))\) and the asymptotic expansion (88) is valid over the entire boundary B.
The proof of Theorem 8 follows the same lines as in the proof of [23, Corollary 5.3] (see also [27, Theorem 3.1 and Corollary 2.1]) and is based on the standard matching of asymptotic expansions technique (see e.g. [23] and the references therein). The expansion (89) related to Neumann problem is similar to the one given in [45, Equation 1.15] for the Dirichlet case. Furthermore, using the references provided in [23] (see the proof of Corollary 5.3) we can deduce that, for the two-dimensional Neumann boundary value problem, \(\eta ^{\mathrm{slow}}\) decays exponentially in \(B^{\epsilon }_{>}\) as \(k \rightarrow -\infty \) which implies the assumption of its rapid decay in the sense of Schwarz in Theorem 8.
B Auxiliary results
Here we provide auxiliary results used in the proofs.
Lemma 8
[24, Lemma 14] Let \(a(s,k) = k^{\theta } \, b(s) \, \varphi (k^{\omega }\varUpsilon (s))\) where \(b,\varphi \) and \(\varUpsilon \) are smooth functions, b and \(\varUpsilon \) are periodic, and \(\theta \in {\mathbb {R}} \backslash {\mathbb {N}}\) and \(\omega \in {\mathbb {R}} \backslash {\mathbb {Z}}_{+}\). Then
for all \(n,m \in {\mathbb {Z}}_{+}\) and all \(k >0\).
Theorem 9
[47, Corollary 3.12] Given a function \(f \in C^{\infty }([a,b])\) and \(n \in {\mathbb {Z}}_+\), there exists a constant \(C_n > 0\) such that
for all \(d \in {\mathbb {N}}\) with \(d+1 \ge n\).
Lemma 9
[24, Lemma 14] Suppose that either \([\alpha ,\beta ] \subseteq [t_{1},t_{2}] \subseteq (c,d)\) or \([\alpha ,\beta ] \cap (t_{1},t_{2}) = \varnothing \) and \([c,d] \subseteq (t_1,t_2)\). Then, for any \(a,b \in {\mathbb {R}}\), \(n \in {\mathbb {N}} \cup \{ 0 \}\), \(m \in {\mathbb {N}}\), there holds
Here we have
when \(2n-(p+q+j) = -1\), and
when \(2n-(p+q+j) \ne -1\).
Rights and permissions
About this article
Cite this article
Ecevit, F., Boubendir, Y., Anand, A. et al. Spectral Galerkin boundary element methods for high-frequency sound-hard scattering problems. Numer. Math. 150, 803–847 (2022). https://doi.org/10.1007/s00211-022-01269-0
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00211-022-01269-0