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

12.07.2015 Views

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 (∂Ω).

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 (∂Ω).

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!