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
30 3 Boundary Integral EquationsIntroducing the integral operatorsV κ : H −1/2 (∂Ω) → Hloc 1 (Ω), ( V κ s)(x) :=W κ : H 1/2 (∂Ω) → Hloc 1 (Ω), (W κt)(x) :=∂Ω∂Ωv κ (x, y)s(y) ds y , (3.3)∂v κ∂n y(x, y)t(y) ds y (3.4)with density functions s, t: ∂Ω → R allows us to rewrite the representation formulae (3.1)and (3.2) asu = V κ γ 1,int u − W κ γ 0,int u in Ω,u = − V κ γ 1,ext u + W κ γ 0,ext u in Ω ext .Remark 3.5. The operators V κ and W κ are usually called the single layer potential operatorand the double layer potential operator, respectively. However, there is a naming conflict.In the following sections we also use these terms for composite operators γ 0 Vκ and γ 0 W κ ,where γ 0 denotes the interior or exterior Dirichlet trace operator.Theorem 3.6. The operator V κ : H −1/2 (∂Ω) → Hloc 1 (Ω) is linear and continuous. Hence,for bounded domains there exists a constant c ∈ R + such that∥ V κ s∥ H 1 (Ω) ≤ c∥s∥ H −1/2 (∂Ω)for all s ∈ H −1/2 (∂Ω).Moreover, for all s ∈ H −1/2 (∂Ω) the function V κ s satisfies the Helmholtz equation in theweak sense (including the Sommerfeld radiation condition in the case of an unboundeddomain).Theorem 3.7. The operator W κ : H 1/2 (∂Ω) → Hloc 1 (Ω) is linear and continuous. Hence,for bounded domains there exists a constant c ∈ R + such that∥W κ t∥ H 1 (Ω) ≤ c∥t∥ H 1/2 (∂Ω)for all t ∈ H 1/2 (∂Ω).Moreover, for all t ∈ H 1/2 (∂Ω) the function W κ t satisfies the Helmholtz equation in theweak sense (including the Sommerfeld radiation condition in the case of an unboundeddomain).The preceding theorems allow us to seek the solution to the Helmholtz equation in theform u = V κ s or u = W κ t with unknown density functions s, t. This approach gives rise toindirect boundary element methods, which will be mentioned later.Properties of the above given potential operators, which will be discussed in the followingsections, can be found in [12] (see also [9] and [18]). These properties will be crucialfor the derivation of boundary integral equations (both direct and indirect), which will beused for the computation of the missing Cauchy data.
313.1 Single Layer Potential OperatorLet us first consider the operator V κ defined by (3.3). Recall, that for the Dirichlet traceoperator γ 0 we haveγ 0 : H 1 loc (Ω) → H1/2 (∂Ω).Combining this operator with V κ we obtain the single layer potential operatorV κ : H −1/2 (∂Ω) → H 1/2 (∂Ω), V κ := γ 0 Vκ .From linearity and continuity of γ 0 and V κ we get that the single layer potential operatoris linear and continuous, i.e., there exists a constant c ∈ R + such that∥V κ s∥ H 1/2 (∂Ω) ≤ c∥s∥ H −1/2 (∂Ω)for all s ∈ H −1/2 (∂Ω).Theorem 3.8. The single layer potential operator V κ : H −1/2 (∂Ω) → H 1/2 (∂Ω) is coercive.Proof. From Theorem 6.22 in [18] we have that he operator V 0 corresponding to the Laplaceequation, i.e., the Helmholtz equation with κ = 0, is H −1/2 (∂Ω)-elliptic. Moreover, theoperator C := V 0 − V κ : H −1/2 (∂Ω) → H 1/2 (∂Ω) is compact (see [18], Section 6.9). Thus,we have⟨(V κ + C)s, s⟩ = ⟨V 0 s, s⟩ ≥ c∥s∥ 2 H −1/2 (∂Ω)for all s ∈ H −1/2 (∂Ω),which completes the proof.Theorem 3.9. For s ∈ L ∞ (∂Ω) there holds the representation(V κ s)(x) = v κ (x, y)s(y) ds y for x ∈ ∂Ω.Proof. The proof is similar to the proof of Lemma 6.7 in [18].∂ΩMoreover, for the jump of the Dirichlet trace of the single layer potential V κ s on theboundary we have[γ 0 Vκ s] := γ 0,ext Vκ s − γ 0,int Vκ s = 0 for all s ∈ H −1/2 (∂Ω). (3.5)3.2 Adjoint Double Layer Potential OperatorIn Section 1.4 we introduced the Neumann trace operatorγ 1 : H 1 loc (Ω, ∆ + κ2 ) → H −1/2 (∂Ω).
- 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 38 and 39: 26 2 Helmholtz Equationand thus l
- Page 41: 293 Boundary Integral EquationsAt t
- 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
- Page 89 and 90: 77E N Err N Err N,p Err ϑ3962 1983
- Page 91 and 92: 792002000.81500.81500.60.41000.60.4
313.1 Single Layer Potential OperatorLet us first consider <strong>the</strong> operator V κ defined by (3.3). Recall, that <strong>for</strong> <strong>the</strong> Dirichlet traceoperator γ 0 we haveγ 0 : H 1 loc (Ω) → H1/2 (∂Ω).Combining this operator with V κ we obtain <strong>the</strong> single layer potential operatorV κ : H −1/2 (∂Ω) → H 1/2 (∂Ω), V κ := γ 0 Vκ .From linearity and continuity of γ 0 and V κ we get that <strong>the</strong> single layer potential operatoris linear and continuous, i.e., <strong>the</strong>re exists a constant c ∈ R + such that∥V κ s∥ H 1/2 (∂Ω) ≤ c∥s∥ H −1/2 (∂Ω)<strong>for</strong> all s ∈ H −1/2 (∂Ω).<strong>The</strong>orem 3.8. <strong>The</strong> single layer potential operator V κ : H −1/2 (∂Ω) → H 1/2 (∂Ω) is coercive.Proof. From <strong>The</strong>orem 6.22 in [18] we have that he operator V 0 corresponding to <strong>the</strong> Laplaceequation, i.e., <strong>the</strong> <strong>Helmholtz</strong> equation with κ = 0, is H −1/2 (∂Ω)-elliptic. Moreover, <strong>the</strong>operator C := V 0 − V κ : H −1/2 (∂Ω) → H 1/2 (∂Ω) is compact (see [18], Section 6.9). Thus,we have⟨(V κ + C)s, s⟩ = ⟨V 0 s, s⟩ ≥ c∥s∥ 2 H −1/2 (∂Ω)<strong>for</strong> all s ∈ H −1/2 (∂Ω),which completes <strong>the</strong> proof.<strong>The</strong>orem 3.9. For s ∈ L ∞ (∂Ω) <strong>the</strong>re holds <strong>the</strong> representation(V κ s)(x) = v κ (x, y)s(y) ds y <strong>for</strong> x ∈ ∂Ω.Proof. <strong>The</strong> proof is similar to <strong>the</strong> proof of Lemma 6.7 in [18].∂ΩMoreover, <strong>for</strong> <strong>the</strong> jump of <strong>the</strong> Dirichlet trace of <strong>the</strong> single layer potential V κ s on <strong>the</strong>boundary we have[γ 0 Vκ s] := γ 0,ext Vκ s − γ 0,int Vκ s = 0 <strong>for</strong> all s ∈ H −1/2 (∂Ω). (3.5)3.2 Adjoint Double Layer Potential OperatorIn Section 1.4 we introduced <strong>the</strong> Neumann trace operatorγ 1 : H 1 loc (Ω, ∆ + κ2 ) → H −1/2 (∂Ω).