high frequency boundary element methods simon chandler-wilde
TRANSCRIPT
![Page 1: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/1.jpg)
High Frequency Boundary Element Methods
Simon Chandler-Wilde
University of Reading, UK
Joint work with: Steve Langdon
Roland Potthast, Eric Heinemeyer, Gottingen
Peter Monk, Delaware
PhD Students Chris Ross, Mizanur Rahman
Funded by: Leverhulme Trust, EPSRC
Dundee, June 2005
1
![Page 2: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/2.jpg)
Overview
• The scientific computing problem
• Two scattering problems and their integral equation formulation
• Review of high frequency boundary element methods and analysis
• Some of our own results on:
• Understanding solution behaviour at high frequency
• Approximating the solution efficiently
• Numerical computations
• Understanding dependence of conditioning on frequency
2
![Page 3: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/3.jpg)
The Scientific Computing Problem
Snap-shot of a component of a time harmonic electromagnetic field.
(Courtesy of Weng Cho Chew, Centre for Computational
Electromagnetics, Illinois, whose team develops fast multipole methods
to solve full systems with > 107 unknowns.)
3
![Page 4: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/4.jpg)
Acoustic Waves
Wave equation, wave speed c:
∆U =1c2∂2U
∂t2
Seek time harmonic solution U(x, t) = u(x)e−iωt. Then u satisfies the
Helmholtz equation
∆u+ k2u = 0
with k = ω/c = 2π × frequency/c > 0.
Simplest solution is the plane wave
u(x) = eikx1 ⇒ U(x, t) = ei(kx−ωt)
with wavelength λ = 2π/k.
4
![Page 5: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/5.jpg)
Focus of This Talk (and my Research)
For
∆u+ k2u = 0,
and numerical methods for its solution:
1. How does everything depend on k?
Everything = solution, error estimates, conditioning, complexity, ...
2. How can we remove or reduce this dependence?
5
![Page 6: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/6.jpg)
Scattering Problem 1. (2D/3D)
@@
@@@R ui, incident wave
∆u + k2u = 0
u = 0Γ
D
obstacle
-
6
x1
x2
6
![Page 7: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/7.jpg)
@@
@@@R ui, incident wave
∆u + k2u = 0
u = 0Γ
D
obstacle
-
6
x1
x2
Green’s representation theorem:
u(x) = ui(x)−∫
Γ
Φ(x, y)∂u
∂n(y)ds(y), x ∈ D,
where Φ(x, y) := i4H
(1)0 (k|x− y|) (2D), :=
14π
eik|x−y|
|x− y|(3D).
7
![Page 8: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/8.jpg)
@@
@@@R ui, incident wave
∆u + k2u = 0
u = 0Γ
D
obstacle
-
6
x1
x2
From Green’s representation theorem (Burton & Miller 1971):
12∂u
∂n(x) +
∫Γ
(∂Φ(x, y)∂n(x)
+ iηΦ(x, y))∂u
∂n(y)ds(y) = f(x), x ∈ Γ,
where
f(x) :=∂ui
∂n(x) + iηui(x), η > 0.
8
![Page 9: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/9.jpg)
@@
@@@R ui, incident wave
∆u + k2u = 0
u = 0Γ
D
2D polygon
-
6
x1
x2
From Green’s representation theorem:
12∂u
∂n(x) +
∫Γ
(∂Φ(x, y)∂n(x)
+ iηΦ(x, y))∂u
∂n(y)ds(y) = f(x), x ∈ Γ.
Theorem (follows from Burton & Miller 1971, Selepov 1969) If η ∈ R,
η 6= 0, then this integral equation is uniquely solvable in L2(Γ).
9
![Page 10: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/10.jpg)
@@
@@@R ui, incident wave
∆u + k2u = 0
u = 0Γ
D
obstacle
-
6
x1
x2
From Green’s representation theorem (Burton & Miller 1971):
12∂u
∂n(x) +
∫Γ
(∂Φ(x, y)∂n(x)
+ iηΦ(x, y))∂u
∂n(y)ds(y) = f(x), x ∈ Γ.
Open Problem 1: show well-posed in L2(Γ) for general Lipschitz Γ.
10
![Page 11: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/11.jpg)
12∂u
∂n(x) +
∫Γ
(∂Φ(x, y)∂n(x)
+ iηΦ(x, y))∂u
∂n(y)ds(y) = f(x), x ∈ Γ.
Conventional BEM: Apply a Galerkin method, approximating ∂u/∂n
by a piecewise polynomial of degree P , leading to a linear system to
solve with N degrees of freedom.
Open Problem 2: Prove stability for general Lipschitz scatterer. (More
complex formulation for which coercivity holds in Buffa & Hiptmair,
2003.)
11
![Page 12: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/12.jpg)
12∂u
∂n(x) +
∫Γ
(∂Φ(x, y)∂n(x)
+ iηΦ(x, y))∂u
∂n(y)ds(y) = f(x), x ∈ Γ.
Conventional BEM: Apply a Galerkin method, approximating ∂u/∂n
by a piecewise polynomial of degree P , leading to a linear system to
solve with N degrees of freedom.
Problem: N of order of (kL)d−1, where L is diameter, d = 2, 3 the
dimension, so cost is O(N2) to compute full matrix and apply iterative
solver ... or close to O(N) if a fast multipole method (e.g. Amini &
Profit 2003, Darve 2004) is used.
12
![Page 13: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/13.jpg)
12∂u
∂n(x) +
∫Γ
(∂Φ(x, y)∂n(x)
+ iηΦ(x, y))∂u
∂n(y)ds(y) = f(x), x ∈ Γ.
Conventional BEM: Apply a Galerkin method, approximating ∂u/∂n
by a piecewise polynomial of degree p, leading to a linear system to solve
with N degrees of freedom.
Problem: N of order of (kL)d−1, where L is diameter, d = 2, 3 the
dimension, so cost is O(N2) to compute full matrix and apply iterative
solver ... or close to O(N) if a fast multipole method (e.g. Amini &
Profit 2003, Darve 2004) is used.
This is fantastic but still infeasible as kL→∞.
13
![Page 14: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/14.jpg)
12∂u
∂n(x) +
∫Γ
(∂Φ(x, y)∂n(x)
+ iηΦ(x, y))∂u
∂n(y)ds(y) = f(x), x ∈ Γ.
Alternative: Reduce N by using new basis functions, e.g.
(i) approximate ∂u/∂n by taking a large number of plane waves and
multiplying these by conventional piecewise polynomial basis functions
(Perry-Debain et al. 2003, 2004). This is very successful (in 2D, 3D,
for acoustic/elastic waves and Neumann/impedance b.c.s),
reducing number of degrees of freedom per wavelength from e.g.
6-10 to close to 2. However N still increases proportional to (kL)d−1.
14
![Page 15: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/15.jpg)
12∂u
∂n(x) +
∫Γ
(∂Φ(x, y)∂n(x)
+ iηΦ(x, y))∂u
∂n(y)ds(y) = f(x), x ∈ Γ.
Alternative: Reduce N by using new basis functions, e.g.
(ii) for convex scatterers, remove some of the oscillation by factoring out
the oscillation of the incident wave, i.e. writing
∂u
∂n(y) =
∂ui
∂n(y)× F (y)
and approximating F by a conventional BEM. (E.g. RJ Uncles (Dundee
NA Report 1974), Abboud, Nedelec, Zhou 1994, Darrigrand 2002, Bruno
et al 2004).
15
![Page 16: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/16.jpg)
Alternative: Reduce N by using new basis functions, e.g.
(ii) for convex scatterers, remove some of the oscillation by factoring out
the oscillation of the incident wave, i.e. writing
∂u
∂n(y) =
∂ui
∂n(y)× F (y) (∗)
and approximating F by a conventional BEM.
For smooth obstacles this works well: equation (∗) holds with
F (y) ≈ 2 on the illuminated side (physical optics) and F (y) ≈ 0 in the
shadow zone.
An analysis will appear in Dominguez, Graham, & Smyshlyaev, in
preparation: for a circle/sphere guaranteed error bounds with cost
O(k1/9).
16
![Page 17: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/17.jpg)
(ii) for convex scatterers, remove some of the oscillation by factoring out
the oscillation of the incident wave, i.e. writing
∂u
∂n(y) =
∂ui
∂n(y)× F (y) (∗)
and approximating F by a conventional BEM. Not very effective for
non-smooth scatterers.
17
![Page 18: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/18.jpg)
Understanding solution behaviour in the case of a
2D convex polygon
18
![Page 19: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/19.jpg)
Ω Γ
Let
G(x, y) := Φ(x, y)− Φ(x, y′)
be the Dirichlet Green function for the left half-plane Ω. By Green’s
representation theorem,
u(x) = ui(x) + ur(x) +∫
∂Ω\Γ
∂G(x, y)∂n(y)
u(y)ds(y), x ∈ Ω,
19
![Page 20: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/20.jpg)
γ
In the left half-plane Ω,
u(x) = ui(x) + ur(x) +∫
∂Ω\Γ
∂G(x, y)∂n(y)
u(y)ds(y)
⇒ ∂u
∂n(x) = 2
∂ui
∂n(x)+2
∫∂Ω\Γ
∂2Φ(x, y)∂n(x)∂n(y)
u(y)ds(y), x ∈ γ = ∂Ω∩Γ.
20
![Page 21: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/21.jpg)
γ
Explicitly, where s is distance along γ, and
φ(s) and ψ(s) are k−1∂u/∂n and u, at distance s along γ,
φ(s) = P.O.+i2
[eiksv+(s) + e−iksv−(s)
]where
v+(s) := k
∫ 0
−∞F
(k(s− s0)
)e−iks0ψ(s0)ds0.
and F (z) := e−izH(1)1 (z)/z
21
![Page 22: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/22.jpg)
φ(s) = P.O.+i2
[eiksv+(s) + e−iksv−(s)
]where
v+(s) := k
∫ 0
−∞F
(k(s− s0)
)e−iks0ψ(s0)ds0.
Now F (z) := e−izH(1)1 (z)/z which is non-oscillatory, in that
F (n)(z) = O(z−3/2−n) as z →∞.
22
![Page 23: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/23.jpg)
φ(s) = P.O.+i2
[eiksv+(s) + e−iksv−(s)
]where
v+(s) := k
∫ 0
−∞F
(k(s− s0)
)e−iks0ψ(s0)ds0.
Now F (z) := e−izH(1)1 (z)/z which is non-oscillatory, in that
F (n)(z) = O(z−3/2−n) as z →∞.
⇒ v(n)+ (s) = O(kn(ks)−1/2−n) as ks→∞.
23
![Page 24: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/24.jpg)
φ(s) = P.O.+i2
[eiksv+(s) + e−iksv−(s)
]where
k−n|v(n)+ (s)| = O
((ks)−1/2−n
)as ks→∞
and (by separation of variables local to the corner),
k−n|v(n)+ (s)| = O
((ks)−α−n
)as ks→ 0,
where α < 1/2 depends on the corner angle.
24
![Page 25: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/25.jpg)
φ(s) = P.O.+i2
[eiksv+(s) + e−iksv−(s)
]where
k−n|v(n)+ (s)| =
O((ks)−1/2−n
)as ks→∞
O ((ks)−α−n) as ks→ 0,
where α < 1/2 depends on the corner angle.
Thus approximate
φ(s) ≈ P.O.+i2
[eiksV+(s) + e−iksV−(s)
],
where V+ and V− are piecewise polynomials on graded meshes.
25
![Page 26: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/26.jpg)
k−n|v(n)+ (s)| =
O((ks)−1/2−n
)as ks→∞
O ((ks)−α−n) as ks→ 0.
Thus approximate
φ(s) ≈ P.O.+i2
[eiksV+(s) + e−iksV−(s)
],
where V+ and V− are piecewise polynomials on graded meshes.
s = 0 tm =(
mN
)qλ
tm = crm
26
![Page 27: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/27.jpg)
Thus approximate
φ(s) ≈ P.O.+i2
[eiksV+(s) + e−iksV−(s)
],
where V+ and V− are piecewise polynomials on graded meshes.
Figure 1: Scattering by a square
27
![Page 28: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/28.jpg)
Thus approximate
φ(s) ≈ P.O.+i2
[eiksV+(s) + e−iksV−(s)
],
where V+ and V− are piecewise polynomials on graded meshes.
Figure 2: Scattering by a square
28
![Page 29: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/29.jpg)
Thus approximate
φ(s) ≈ P.O.+i2
[eiksV+(s) + e−iksV−(s)
],
where V+ and V− are piecewise polynomials on graded meshes.
Figure 3: Scattering by a square
29
![Page 30: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/30.jpg)
Thus approximate
φ(s) ≈ P.O.+i2
[eiksV+(s) + e−iksV−(s)
],
where V+ and V− are piecewise polynomials on graded meshes.
Theorem Where φN is the best approximation in L2(Γ) from the
approximation space, n is the number of sides, N the number of degrees
of freedom, and p the polynomial degree,
k1/2||φ− φN ||2 ≤ C supx∈D
|u(x)| [n(1 + log(kL/n))]p+3/2
Np+1
√logN.
30
![Page 31: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/31.jpg)
Numerical results
scattering by a square, k = 5
31
![Page 32: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/32.jpg)
Numerical results (scattering by a square)
Solution minus P.O. approximation;
32
![Page 33: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/33.jpg)
Numerical results (scattering by a square)
Correction to P.O. approximation;
33
![Page 34: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/34.jpg)
Numerical results (scattering by a square)
Correction to P.O. approximation;
34
![Page 35: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/35.jpg)
Numerical results (scattering by a square)
Correction to P.O. approximation;
35
![Page 36: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/36.jpg)
Numerical results (scattering by a square)
Correction to P.O. approximation;
36
![Page 37: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/37.jpg)
Numerical results (scattering by a square)
Correction to P.O. approximation, k = 5;
37
![Page 38: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/38.jpg)
Numerical results (scattering by a square)
Correction to P.O. approximation, k = 10;
38
![Page 39: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/39.jpg)
Numerical results (scattering by a square)
Correction to P.O. approximation, k = 20;
39
![Page 40: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/40.jpg)
Table 1: Relative errors, k = 10
k N dof ‖φ− φN‖2/‖φ‖2 EOC
10 2 24 1.1187×10+0 1.5
4 48 4.0499×10−1 0.7
8 88 2.5348×10−1 0.9
16 176 1.3979×10−1 1.3
32 360 5.5216×10−2 0.9
64 712 3.0358×10−2
40
![Page 41: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/41.jpg)
Table 2: Relative errors, k = 160
k N dof ‖φ− φN‖2/‖φ‖2 EOC
160 2 32 1.0350×10+0 1.3
4 56 4.2389×10−1 0.5
8 120 3.0406×10−1 0.6
16 240 2.0471×10−1 1.5
32 472 7.3763×10−2 1.0
64 944 3.6983×10−2
41
![Page 42: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/42.jpg)
What we are actually computing . . .
The difference between the exact solution and the leading order physical
optics/Kirchhoff approximation;
Figure 4: square, k = 5
42
![Page 43: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/43.jpg)
What we are actually computing . . .
The difference between the exact solution and the leading order physical
optics/Kirchhoff approximation;
Figure 5: square, k = 10
43
![Page 44: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/44.jpg)
What we are actually computing . . .
The difference between the exact solution and the leading order physical
optics/Kirchhoff approximation;
Figure 6: square, k = 20
44
![Page 45: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/45.jpg)
What we are actually computing . . .
The difference between the exact solution and the leading order physical
optics/Kirchhoff approximation;
Figure 7: square, k = 40
45
![Page 46: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/46.jpg)
The main remaining computational difficulty
Efficient evaluation of integrals of the form∫ ym+1
ym
∫ yj+1
yj
H(1)0 (k
√(s− a)2 + (t− b)2)eik(t+s) dtds,
∫ ym+1
ym
∫ yj+1
yj
(ct+ d)H(1)1 (k
√(s− a)2 + (t− b)2)√
(s− a)2 + (t− b)2)eik(t+s) dtds.
46
![Page 47: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/47.jpg)
The main remaining computational difficulty
... or more generally ...∫ ym+1
ym
∫ yj+1
yj
H(1)0 (k
√(s− a)2 + (t− b)2)pM (s)pN (t)eik(t+s) dtds,
Open Problem 3: Efficient evaluation of these oscillatory integrals.
(Darrigrand 2002, Ganesh, Langdon, Sloan, in preparation, and cf.
Iserles 2004, 2005.)
47
![Page 48: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/48.jpg)
The main theoretical difficulty: stability
48
![Page 49: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/49.jpg)
For 2D polygon no problem with ‘conventional’ stability
• Integral equation is well-posed (coercive), so ||(I −K)−1|| ≤ C
• Coercivity implies Galerkin method is stable as N →∞
49
![Page 50: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/50.jpg)
No problem with ‘conventional stability’... but
• Integral equation is well-posed (coercive), so ||(I −K)−1|| ≤ C(k)
• Coercivity implies Galerkin method is stable for N ≥ c(k)
50
![Page 51: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/51.jpg)
Work in progress: estimating C(k)
• Integral equation is well-posed (coercive), so ||(I −K)−1|| ≤ C(k)
For a star-like (in particular convex) domain we can use integration by
parts and trace theorems. Cf.
Melenk, PhD, 1995, Cummings & Feng to appear M3AS (inf-sup
constant O(k) for interior problems)
C-W & Monk to appear SIAM J Math Anal (inf-sup constant O(k3) for
variational formulation of exterior (rough surface) scattering problem)
Jerison & Kenig 1981, Vechota 1984
(Rellich-Payne-Weinberg-Necas-type identities)
51
![Page 52: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/52.jpg)
Work in progress: estimating C(k)
• Integral equation is well-posed (coercive), so ||(I −K)−1|| ≤ C(k)
For a star-like (in particular convex) domain we can use integration by
parts and trace theorems. Cf.
Melenk, PhD, 1995, Cummings and X. Feng to appear M3AS (inf-sup
constant O(k) for interior problems)
C-W & Monk to appear SIAM J Math Anal (inf-sup constant O(k3) for
variational formulation of exterior (rough surface) scattering problem)
Jerison & Kenig 1981, Vechota 1984
(Rellich-Payne-Weinberg-Necas-type identities)
For a circle/sphere, C(k) = 1 !! (Dominguez, Graham,
Smyshlyaev, in preparation)
Open Problem 4: Any bound on C(k) for general scatterer.
52
![Page 53: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/53.jpg)
Estimating c(k)??
• Coercivity implies Galerkin method is stable for N ≥ c(k)
||(I −KN )−1|| ≤ 2C(k), for N ≥ c(k).
Appears hopeless to estimate c(k), unless we apply least squares
method, i.e. apply Galerkin method to
(I −K∗)(I −K)φ = (I −K∗)ψ,
but poor conditioning (as k →∞), more difficult quadrature?
Open Problem 5: Estimate c(k) (skated over in Buffa & Sauter,
Preprint).
53
![Page 54: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/54.jpg)
Scattering Problem 2. (2D/3D)
∆u + k2u = g
u = 0
∂D
D ⊂ Rn
ZZ HHH
@@
@@@
PPP
~Support of g
-
6
6
x = (x1, ..., xn−1)
f(x)xn
54
![Page 55: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/55.jpg)
Scattering Problem 2. (2D/3D)
∆u + k2u = g
u = 0
∂D
D ⊂ Rn
ZZ HHH
@@
@@@
PPP
~Support of gxn = f+
xn = f−-
6
6
x = (x1, ..., xn−1)
f(x)xn
For variational formulation: (inf-sup constant)−1 ≤ 1 +√
2κ(κ+ 1)2
where κ = k(f+ − f−) (C-W & Monk, to appear)
55
![Page 56: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/56.jpg)
Scattering Problem 2. (2D/3D)
∆u + k2u = g
u = 0
∂D
D ⊂ Rn
ZZ H
HH@@
@@@PPP
~Support of g
-
6
6
x = (x1, ..., xn−1)
f(x)xn
For integral equation: ||(I −K)−1||L2(∂D) ≤ 1 + 10(1 + L2)3/4
where L is the Lipschitz constant of f . (C-W, Heinemeyer, Potthast,
preprints.)
56
![Page 57: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/57.jpg)
Summary and Conclusions
• Reviewed work on high frequency BEMs
• Described technique to understand behaviour of the field on the
boundary for scattering by a convex polygon (extends to convex
polyhedron in 3D)
• For a convex polygon, design of an optimal graded mesh for
piecewise polynomial approximation is then straightforward
• The number of degrees of freedom need only grow logarithmically
with the wavenumber to maintain a fixed accuracy
• General approximation idea (but not the simple analysis techniques)
extends to more general problems
• Mentioned results on estimating k dependence of norms of inverse
operators for variational/integral equation formulations for rough
surface scattering
57
![Page 58: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/58.jpg)
References
• C-W, Rahman, Ross Num Math 2002 (Conventional BEM for a specific
2D problem with dependence of everything on k specified.)
• C-W, Langdon, Ritter Phil Trans R Soc 2004, Langdon, C-W to appear
SIAM J Numer Anal
(High frequency BEM for same problem with DoF = O(1) as k →∞)
• C-W, Langdon, in preparation.
(Convex polygon analysis/numerics. DoF = O(log k))
• C-W, Monk, to appear SIAM J Math Anal
(Dependence of inf-sup constant on k for rough surface scattering.)
• C-W, Heinemeyer, Potthast, preprints. (Dependence of inverse
operator/condition number on k for Brakhage-Werner integral equation
for rough surface scattering.)
58
![Page 59: High Frequency Boundary Element Methods Simon Chandler-Wilde](https://reader033.vdocument.in/reader033/viewer/2022042914/626a4860ab04aa7aed18d4a8/html5/thumbnails/59.jpg)
Newton Institute Programme
Highly Oscillatory Problems
January-June 2007
Opportunity to address many of these open problems.
Workshop announcements in this area coming soon.
59