The Boundary Element Method for the Helmholtz Equation ... - FEI VÅ B
The Boundary Element Method for the Helmholtz Equation ... - FEI VÅ B The Boundary Element Method for the Helmholtz Equation ... - FEI VÅ B
26 2 Helmholtz Equationand thus lim∂u 2ε→∞ ∂B ε(0) ∂n (y)Since+ κ 2 |u(y)| 2 ds + 2κ Im u(y) ∂ū∂B ε(0) ∂n (y) ds = 0. (2.10)Now we apply the first Green’s identity (1.7) in Ω ε for functions u and ū to obtain∂ū∆ū(y)u(y) dy + ∇ū(y)∇u(y) dy = (y)u(y) ds. (2.11)Ω ε Ω ε∂n∆ū + κ 2 ū = 0in Ω εand ∂Ω ε = ∂Ω ∪ ∂B ε (0), we can rewrite the formula (2.11) as|∇u(y)| 2 dy − κ 2 |u(y)| 2 ∂ūdy =Ω ε Ω ε ∂B ε(0) ∂n∂Ω(y)u(y) ds − ∂ū(y)u(y) ds, (2.12)∂nwhere we had to take into account the opposite direction of n for y ∈ ∂Ω (see Figure 2.3).Because the left-hand side of (2.12) is real, we obtain∂ūIm∂n∂Ω(y)u(y) ds = Im ∂ū(y)u(y) ds < ∞. (2.13)∂n∂B ε(0)Inserting (2.13) into (2.10) we getand thus lim∂u 2ε→∞ ∂B ε(0) ∂n (y)limε→∞∂B ε(0) ∂u 2∂n (y)∂Ω ε+ κ 2 |u(y)| 2 ds + 2κ Im u(y) ∂ū∂Ω ∂n (y) ds = 0+ κ 2 |u(y)| 2 ds = −2κ Im∂Ωu(y) ∂ū (y) ds < ∞.∂nBecause both terms on the left-hand side are non-negative, they have to be bounded forε → ∞. Therefore, we have proved that|u(y)| 2 ds = O(1) for ε → ∞. (2.14)∂B ε(0)From Theorem 2.3 we know that the radiation condition and thus also (2.14) is validfor the fundamental solution v κ . Using the Hölder inequality we thus obtain ∂vκI 1 := u(y) (x, y) − iκv κ (x, y) ds y∂B ε(0) ∂n y 1/2 ≤ |u(y)| 2 ds y ∂v κ2 1/2 (x, y) − iκv κ (x, y)∂B ε(0)∂B ε(0) ∂n y ds y → 0 for ε → ∞(2.15)
27andI 2 :=∂B ε(0) ∂uv κ (x, y)∂n (y) − iκu(y) ds y 1/2 ≤ |v κ (x, y)| 2 ds y ∂u2 1/2∂B ε(0)∂B ε(0) ∂n (y) − iκu(y) ds y → 0 for ε → ∞.(2.16)Finally, for the bounded domain Ω ε we can apply the representation formula (2.7) toobtain∂uu(x) =∂Ω ε∂n (y)v κ(x, y) ds y − u(y) ∂v κ(x, y) ds y∂Ω ε∂n y∂u= −∂n (y)v κ(x, y) ds y + u(y) ∂v κ(x, y) ds y + I 2 − I 1∂n y∂Ωfor all x ∈ Ω ε . Letting ε → ∞, using (2.15) and (2.16), we eventually getu(x) = u(y) ∂v κ∂u(x, y) ds y −∂n y ∂n (y)v κ(x, y) ds y∂Ωfor all x ∈ Ω ext , which was to be proved.Theorems 2.6 and 2.7 provide representation of the solution to the Helmholtz equationby means of boundary integrals. The functions( V ∂v κκ s)(x) := v κ (x, y)s(y) ds y and (W κ t)(x) := (x, y)t(y) ds y∂n y∂Ωare called potentials with density functions s, t. Although we have only discussed the caseof smooth data, in the following section we will show that the potentials can also be definedfor more general density functions and domains. The properties of the integral operatorsV κ and W κ will play a crucial role in obtaining the missing Cauchy data.∂Ω∂Ω∂Ω
- Page 3: I hereby declare that this thesis i
- Page 7: AbstractIn this work we study the a
- Page 13 and 14: 1IntroductionThe boundary element m
- Page 15 and 16: 31 Function SpacesIn this section w
- Page 17 and 18: 5∂Ωy i 2ΩU + iU − iΓ iy i 1F
- Page 19 and 20: 7For a special choice of p = 2 we g
- Page 21 and 22: 9Note that the trace operator is no
- Page 23: 11Definition 1.8 (Weak Solution). C
- Page 27: 152 Helmholtz EquationLet us first
- Page 32 and 33: 20 2 Helmholtz Equation∂ΩΩ εε
- Page 34 and 35: 22 2 Helmholtz Equationthe normal v
- Page 36 and 37: 24 2 Helmholtz Equationand due to p
- Page 41 and 42: 293 Boundary Integral EquationsAt t
- Page 43 and 44: 313.1 Single Layer Potential Operat
- Page 45 and 46: 33For the jump of the Neumann trace
- Page 47 and 48: 35with the surface curl operatorcur
- Page 49 and 50: 37The following theorems are standa
- Page 51 and 52: 39u ∈ H 1 (Ω, ∆+κ 2 ) if and
- Page 53 and 54: 41Another approach to computing the
- Page 55 and 56: 433.5.3 Interior Mixed Boundary Val
- Page 57 and 58: 45respectively.The unique solvabili
- Page 59: 47with an unbounded domain Ω ext :
- Page 62 and 63: 50 4 Discretization and Numerical R
- Page 64 and 65: 52 4 Discretization and Numerical R
- Page 66 and 67: 54 4 Discretization and Numerical R
- Page 68 and 69: 56 4 Discretization and Numerical R
- Page 70 and 71: 58 4 Discretization and Numerical R
- Page 72 and 73: 60 4 Discretization and Numerical R
- Page 74 and 75: 62 4 Discretization and Numerical R
- Page 76 and 77: 64 4 Discretization and Numerical R
- Page 78 and 79: 66 4 Discretization and Numerical R
- Page 80 and 81: 68 4 Discretization and Numerical R
- Page 82 and 83: 70 4 Discretization and Numerical R
- Page 84 and 85: 72 4 Discretization and Numerical R
- Page 87 and 88: 755 Numerical ExperimentsIn this pa
27andI 2 :=∂B ε(0) ∂uv κ (x, y)∂n (y) − iκu(y) ds y 1/2 ≤ |v κ (x, y)| 2 ds y ∂u2 1/2∂B ε(0)∂B ε(0) ∂n (y) − iκu(y) ds y → 0 <strong>for</strong> ε → ∞.(2.16)Finally, <strong>for</strong> <strong>the</strong> bounded domain Ω ε we can apply <strong>the</strong> representation <strong>for</strong>mula (2.7) toobtain∂uu(x) =∂Ω ε∂n (y)v κ(x, y) ds y − u(y) ∂v κ(x, y) ds y∂Ω ε∂n y∂u= −∂n (y)v κ(x, y) ds y + u(y) ∂v κ(x, y) ds y + I 2 − I 1∂n y∂Ω<strong>for</strong> all x ∈ Ω ε . Letting ε → ∞, using (2.15) and (2.16), we eventually getu(x) = u(y) ∂v κ∂u(x, y) ds y −∂n y ∂n (y)v κ(x, y) ds y∂Ω<strong>for</strong> all x ∈ Ω ext , which was to be proved.<strong>The</strong>orems 2.6 and 2.7 provide representation of <strong>the</strong> solution to <strong>the</strong> <strong>Helmholtz</strong> equationby means of boundary integrals. <strong>The</strong> functions( V ∂v κκ s)(x) := v κ (x, y)s(y) ds y and (W κ t)(x) := (x, y)t(y) ds y∂n y∂Ωare called potentials with density functions s, t. Although we have only discussed <strong>the</strong> caseof smooth data, in <strong>the</strong> following section we will show that <strong>the</strong> potentials can also be defined<strong>for</strong> more general density functions and domains. <strong>The</strong> properties of <strong>the</strong> integral operatorsV κ and W κ will play a crucial role in obtaining <strong>the</strong> missing Cauchy data.∂Ω∂Ω∂Ω