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
36 3 Boundary Integral Equationswith the Calderon projection matrix 1C int := 2 I − K κV κD κ12 I + K∗ κFor the solution in an unbounded domain we have the formula.u = W κ γ 0,ext u − V κ γ 1,ext u in Ω ext .Applying the Dirichlet and Neumann trace operators, respectively, and using properties(3.12), (3.10) and the definitions of the boundary integral operators V κ and D κ we get 1γ 0,ext u =2 I + K κ γ 0,ext u − V κ γ 1,ext u on ∂Ω ext , (3.17)γ 1,ext u = −D κ γ 0,ext u + 12 I − K∗ κγ 1,ext u on ∂Ω ext . (3.18)Again, the boundary integral equations (3.17), (3.18) can be rewritten as γ 0,ext uγ 1,ext = C ext γ 0,ext uu γ 1,ext uwith the Calderon projection matrix 1C ext := 2 I + K κ−V κ1−D κ 2 I − K∗ κTheorem 3.16. The Calderon operators C int and C ext are projection operators, i.e.,.(C int ) 2 = C int (C ext ) 2 = C ext .Proof. The proof is analogous to the proof of Lemma 6.18 in [18].Corollary 3.17. For the boundary integral operators we have the following identities 1 1V κ D κ = κ2 I + K 2 I − K κ ,D κ V κ = 12 I + K∗ κ 12 I − K∗ κD κ K κ = K ∗ κD κ ,K κ V κ = V κ K ∗ κ.Proof. The relations follow directly from the comparison of the matrices (C int ) 2 and C int .,
37The following theorems are standard results of functional analysis and will be usedto prove solvability of boundary integral equations studied in the next sections. In thetheorems we assume that X denotes a Hilbert space.Theorem 3.18 (Lax-Milgram Lemma). Let A: X → X ∗ denote a linear, bounded andX-elliptic operator. Then for any f ∈ X ∗ there exists a unique element u ∈ X satisfyingMoreover, there holds the estimatewith the ellipticity constant c A 1 .Au = f.∥u∥ X ≤ 1c A ∥f∥ X ∗1Proof. For the proof of the Lax-Milgram Lemma see [18], proof following Theorem 3.4 orany standard functional analysis textbook.Theorem 3.19 (Fredholm Alternative). Let K : X → X denote a compact operator. Eitherthe homogeneous equation(I − K)u = 0has a nontrivial solution u ∈ X or the inhomogeneous problem(I − K)u = fhas a unique solution u ∈ X for all f ∈ X. In the latter case there exists a constant c ∈ R +such that∥u∥ X ≤ c∥f∥ X for all f ∈ X.Proof. For the proof of the Fredholm Alternative see, e.g., [7], Theorem 2.2.9.Theorem 3.20. Let A: X → X ∗ denote a linear, bounded and coercive operator and letA be injective, i.e., from Au = 0 it follows u = 0. Then for every f ∈ X ∗ there exists aunique solution to the equationAu = f. (3.19)Moreover, there exists a constant c ∈ R + such that∥u∥ X ≤ c∥f∥ X ∗ for all f ∈ X ∗ . (3.20)Proof. From coercivity of A we know that there exists a compact operator C such that thelinear operator D := A + C : X → X ∗ is X-elliptic. From the Lax-Milgram Lemma 3.18we get that there exists the inverse operator D −1 : X ∗ → X. Therefore, we obtain thatthe equation (3.19) is equivalent toBu = D −1 f
- 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 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: 35with the surface curl operatorcur
- 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
- Page 93 and 94: 81The approximate solution u h to (
- Page 95 and 96: 83with the unit direction vector d
- Page 97 and 98: 8540.540.430.430.220.320100.20.110
36 3 <strong>Boundary</strong> Integral <strong>Equation</strong>swith <strong>the</strong> Calderon projection matrix 1C int := 2 I − K κV κD κ12 I + K∗ κFor <strong>the</strong> solution in an unbounded domain we have <strong>the</strong> <strong>for</strong>mula.u = W κ γ 0,ext u − V κ γ 1,ext u in Ω ext .Applying <strong>the</strong> Dirichlet and Neumann trace operators, respectively, and using properties(3.12), (3.10) and <strong>the</strong> definitions of <strong>the</strong> boundary integral operators V κ and D κ we get 1γ 0,ext u =2 I + K κ γ 0,ext u − V κ γ 1,ext u on ∂Ω ext , (3.17)γ 1,ext u = −D κ γ 0,ext u + 12 I − K∗ κγ 1,ext u on ∂Ω ext . (3.18)Again, <strong>the</strong> boundary integral equations (3.17), (3.18) can be rewritten as γ 0,ext uγ 1,ext = C ext γ 0,ext uu γ 1,ext uwith <strong>the</strong> Calderon projection matrix 1C ext := 2 I + K κ−V κ1−D κ 2 I − K∗ κ<strong>The</strong>orem 3.16. <strong>The</strong> Calderon operators C int and C ext are projection operators, i.e.,.(C int ) 2 = C int (C ext ) 2 = C ext .Proof. <strong>The</strong> proof is analogous to <strong>the</strong> proof of Lemma 6.18 in [18].Corollary 3.17. For <strong>the</strong> boundary integral operators we have <strong>the</strong> following identities 1 1V κ D κ = κ2 I + K 2 I − K κ ,D κ V κ = 12 I + K∗ κ 12 I − K∗ κD κ K κ = K ∗ κD κ ,K κ V κ = V κ K ∗ κ.Proof. <strong>The</strong> relations follow directly from <strong>the</strong> comparison of <strong>the</strong> matrices (C int ) 2 and C int .,