The Helmholtz Equation
|
|
- Martina Dean
- 5 years ago
- Views:
Transcription
1 The Helmholtz Equation Seminar BEM on Wave Scattering Rene Rühr ETH Zürich October 28, 2010
2 Outline Steklov-Poincare Operator Helmholtz Equation: From the Wave equation to Radiation condition Uniqueness of the Exterior Dirichlet Problem Boundary Integral Operators Representation formula Existence of the Exterior Dirichlet Problem Rene Rühr (ETH Zürich) The Helmholtz Equation October 28, / 34
3 Poincare - Steklov Operator Last time we saw the Calderon projection, enabling us to write the following system of BIE for the Cauchy data γ D u, γ N u ( ) ( γd u 1 = 2 I K ) ( ) 0 V 0 γd u 1 γ N u W 0 2 I + K 0 γ N u (here for for homogeneous case u = 0) where V 0 : H 1 2 (Γ) H 1 2 (Γ) K 0 : H 1 2 (Γ) H 1 2 (Γ) K 0 : H 1 2 (Γ) H 1 2 (Γ) W 0 : H 1 2 (Γ) H 1 2 (Γ) Rene Rühr (ETH Zürich) The Helmholtz Equation October 28, / 34
4 Steklov - Poincare Operator We proved (in 3d) V 0 : H 1 2 (Γ) H 1 2 (Γ) is elliptic invertible! Rewrite first equation: γ N u = V 1 0 ( 1 2 I + K 0)γ D u Define S 0 := V 1 0 ( 1 2 I + K 0) : H 1 2 (Γ) H 1 2 (Γ)...called Steklev-Poincare-Operator S 0 is a Dirichlet-to-Neumann map Rene Rühr (ETH Zürich) The Helmholtz Equation October 28, / 34
5 Steklov - Poincare Operator Second equation of the Calderon projection admits another representation: S 0 = W 0 + ( 1 2 I + K 0)V 1 0 ( 1 2 I + K 0) S 0 admits the same ellipticity estimates as W 0 : Sv, v Γ W 0 v, v Γ H 1 2 (Γ)-elliptic. Rene Rühr (ETH Zürich) The Helmholtz Equation October 28, / 34
6 Outline Steklov-Poincare Operator Helmholtz Equation: From the Wave equation to Radiation condition Uniqueness of the Exterior Dirichlet Problem Boundary Integral Operators Representation formula Existence of the Exterior Dirichlet Problem Rene Rühr (ETH Zürich) The Helmholtz Equation October 28, / 34
7 Wave equation reduces to Helmholtz Equation The wave equation is 2 t 2 Ψ c2 Ψ = 0 assume the solution to be time harmonic: Ψ(t, x) = e iwt u(x) then u will satisfy w 2 u c 2 u = 0 or u k 2 u = 0 with k = w c... called Helmholtz equation. Rene Rühr (ETH Zürich) The Helmholtz Equation October 28, / 34
8 Wave equation reduces to Helmholtz Equation Fundamental solution of u k 2 u = 0 in R 3 G k (x, y) = 1 e ik x y 4π x y x y also 1 e ik x y 4π x y is a fundamental solution - justify our choice later Just as for the Laplace, in n=2, things are a little bit different, Y 0 the second kind Bessel function of order zero, G k (x, y) = 1 2π Y 0(k x y ) x y whose singularity at x y behaves like the log s as s 0 we study the Exterior Dirichlet Problem (EDP) for a bounded domain Ω R 3 u k 2 u = 0 in Ω c and u = g in Ω = Γ Rene Rühr (ETH Zürich) The Helmholtz Equation October 28, / 34
9 Radiation condition Problem: in general won t be unique... for consider A eikr r Helmholtz equation and B e ikr r spherical waves will both satisfy the want our solution of the EDP also behaves like an outgoing wave! so let s look for a condition only u = A eikr r satisfies. r u iku = o( 1 r 2 ) as r Sommerfeld Radiation Condition! Rene Rühr (ETH Zürich) The Helmholtz Equation October 28, / 34
10 Radiation condition consider u k 2 u = f in R n. u = G k f is a solution.. that satifies, putting x = rω A(ω) eikr r + o(r 2 ) as r right choice of sign! Rene Rühr (ETH Zürich) The Helmholtz Equation October 28, / 34
11 Outline Steklov-Poincare Operator Helmholtz Equation: From the Wave equation to Radiation condition Uniqueness of the Exterior Dirichlet Problem Boundary Integral Operators Representation formula Existence of the Exterior Dirichlet Problem Rene Rühr (ETH Zürich) The Helmholtz Equation October 28, / 34
12 Uniqueness Uniqueness of the EDP! Theorem There is at most one radiating solution u H 1 loc (Ωc ) of where g H 1 2 (Γ). u k 2 u = 0 on Ω c γd c u = g on Γ, Rene Rühr (ETH Zürich) The Helmholtz Equation October 28, / 34
13 Uniqueness need the following Lemma established by Rellich: Lemma Let k > 0 and u a solution of u k 2 u = 0 on B ρ0 (0) c and suppose that then u = 0 on B ρ0 (0) c. lim u(x) 2 dσ = 0 ρ x =ρ Rene Rühr (ETH Zürich) The Helmholtz Equation October 28, / 34
14 Uniqueness Corollary and, with a little work, it follows that Let u H 1 loc (Ωc ) be a solution of the homogeneous exterior Helmholtz equation, u k 2 u = 0 on Ω c and satisfies in addition to the Sommerfeld radiation condition also ( ) I k (γ Nū)(γc c D u) dσ 0, then u = 0 on Ω c. Γ Rene Rühr (ETH Zürich) The Helmholtz Equation October 28, / 34
15 Uniqueness Proof: Let s bound our domain: Ω c ρ = Ω c B ρ Green s 1st Ω c ρ ū u k 2 ūu dx = ( k 2 )u, u L 2 (Ω c ρ) + γ N u, γ D u L 2 ( Ω c ρ ) first term vanishes by assumption. Multiply with k and take the imaginary part ( I(k) u 2 k 2 u 2 dx = I k ū ) ( B ρ ν u dσ I k ū ) Γ ν u dσ Ω c ρ Rene Rühr (ETH Zürich) The Helmholtz Equation October 28, / 34
16 Uniqueness last integral can be rewritten since u ρ iku 2 = u ρ I(k) u 2 + k 2 u 2 dx + 1 Ω c 2 ρ + u B ρ ρ iku I(k) > 0 Ω c ρ u 2 dx 0 I(k) = 0 Rellich applies 2 ( + k 2 u 2 + 2I k ū ) ν u B ρ 2 u 2 ρ + k 2 u 2 dσ = ( dσ I k ū ) ν u dσ. Γ Rene Rühr (ETH Zürich) The Helmholtz Equation October 28, / 34
17 Outline Steklov-Poincare Operator Helmholtz Equation: From the Wave equation to Radiation condition Uniqueness of the Exterior Dirichlet Problem Boundary Integral Operators Representation formula Existence of the Exterior Dirichlet Problem Rene Rühr (ETH Zürich) The Helmholtz Equation October 28, / 34
18 Integral operators and Boundary Integral Equations As for the Laplace, we define the boundary integral operators as (V k w)(x) = G k (x, y)w(y)ds y Γ (K k v)(x) = G k (x, y)v(y)ds y Γ n y (K k v)(x) = G k (x, y)v(y)ds y Γ n x (W k v)(x) = G k (x, y)v(y)ds y n x n y Γ Rene Rühr (ETH Zürich) The Helmholtz Equation October 28, / 34
19 Integral operators and Boundary Integral Equations with the same properties as for the Laplace Theorem For a bounded Lipschitz domain, the boundary integral operators are continous. Remark V k : H 1 2 (Γ) H 1 2 (Γ) K k : H 1 2 (Γ) H 1 2 (Γ) K k : H 1 2 (Γ) H 1 2 (Γ) W k : H 1 2 (Γ) H 1 2 (Γ) One has the equalities W k = K k Id and V k = K k Id. Rene Rühr (ETH Zürich) The Helmholtz Equation October 28, / 34
20 Integral operators and Boundary Integral Equations Lemma well... not exactly the same properties: V k no longer elliptic! but we still have V k : H 1 2 (Γ) H 1 2 (Γ) is coercive, i.e. there exists a compact operator C : H 1 2 (Γ) H 1 2 (Γ) such that the Gardings inequality V k w, w L 2 (Γ) + Cw, w L 2 (Γ) const w 2 H 1 2 (Γ) w H 1 2 (Γ) holds. Rene Rühr (ETH Zürich) The Helmholtz Equation October 28, / 34
21 Integral operators and Boundary Integral Equations Theorem the single and double layer potential give us radiating solution: If φ H 1 2 (Γ) then for u = Ψ SL φ satisfies the Sommerfeld radiation condition and ( k 2 )u = 0 in R n \ Γ Analogously for u = Ψ DL φ if φ H 1 2 (Γ). Rene Rühr (ETH Zürich) The Helmholtz Equation October 28, / 34
22 Integral operators and Boundary Integral Equations suppose u H 1 loc (Ωc ) with ( k 2 )u = 0 in Ω Ω c satisfies u = Ψ SL [γ N u] + Ψ DL [γ D u] in Ω Ω c Integral representation formula take traces to get the exterior Calderon projection ( ) ( γ c D u 1 γn c u = 2 I + K ) ( ) k V k γ c D u 1 W k 2 I K k γn c u different signs w.r.t to the interior Calderon Rene Rühr (ETH Zürich) The Helmholtz Equation October 28, / 34
23 Integral operators and Boundary Integral Equations boundary integral equations to the exterior Dirichlet problem: Direct Method - integral equality of the first kind: Let g H 1 2 (Γ). Find ψ H 1 2 (Γ) such that V k ψ, ϕ L 2 (Γ) = 1 2 g, ϕ L 2 (Γ) + K k g, ϕ L 2 (Γ) ϕ H 1 2 (Γ) Direct Method - integral equality of the second kind: Let g H 1 2 (Γ). Find ψ H 1 2 (Γ) such that 1 2 ψ, ϕ L 2 (Γ) + K kψ, ϕ L2 (Γ) = W k g, ϕ L2 (Γ) ϕ H 1 2 (Γ) Rene Rühr (ETH Zürich) The Helmholtz Equation October 28, / 34
24 Integral operators and Boundary Integral Equations Indirect Method - using the single layer potential: Let g H 1 2 (Γ). Find ψ H 1 2 (Γ) such that V k ψ, ϕ L2 (Γ) = g, ϕ L2 (Γ) ϕ H 1 2 (Γ) Indirect Method - using the double layer potential: Let g H 1 2 (Γ). Find ψ H 1 2 (Γ) such that 1 2 ψ, ϕ L 2 (Γ) + K k ψ, ϕ L2 (Γ) = g, ϕ L2 (Γ) ϕ H 1 2 (Γ) Rene Rühr (ETH Zürich) The Helmholtz Equation October 28, / 34
25 Outline Steklov-Poincare Operator Helmholtz Equation: From the Wave equation to Radiation condition Uniqueness of the Exterior Dirichlet Problem Boundary Integral Operators Representation formula Existence of the Exterior Dirichlet Problem Rene Rühr (ETH Zürich) The Helmholtz Equation October 28, / 34
26 Representation Formula for the exterior Calderon operator we assumed Greens representation formula.. does it hold? Theorem yes - thanks to the radiating property! Let g H 1 2 (Γ). And suppose u H 1 loc (Ωc ) is a radiating solution of u k 2 u = 0 on Ω c then u has the integral representation γd c u = g on Γ, u = Ψ DL g Ψ SL (γ c N u). Rene Rühr (ETH Zürich) The Helmholtz Equation October 28, / 34
27 Representation Formula Proof idea: Ω c ρ = Ω c B ρ Green u(x) = +Ψ SL (γn c u)(x) Ψ DL(γD c u)(x) G(x, y)γ N u(y) dσ + γ N G(x, y)γ D u(y) dσ B ρ B ρ Sommerfeld radiation condition two last terms vanish as ρ 0 Rene Rühr (ETH Zürich) The Helmholtz Equation October 28, / 34
28 Outline Steklov-Poincare Operator Helmholtz Equation: From the Wave equation to Radiation condition Uniqueness of the Exterior Dirichlet Problem Boundary Integral Operators Representation formula Existence of the Exterior Dirichlet Problem Rene Rühr (ETH Zürich) The Helmholtz Equation October 28, / 34
29 Existence Theorem we re ready to prove the existence!..by showing that we can solve the BIE. so we need Let g H 1 2 (Γ). And suppose u H 1 loc (Ωc ) is a radiating solution of u k 2 u = 0 on Ω c γ c Du = g on Γ, then φ = γ c N H 1 2 (Γ) is a solution of the boundary integral equation and u has the integral representation V k φ = ( Id /2 + K)g on Γ, (1) u = Ψ DL g Ψ SL φ. (2) Conversely, if φ H 1 2 (Γ) is a solution (1), then formula (2) defines a solution u H 1 loc (Ωc ) of the exterior Dirichlet problem. Rene Rühr (ETH Zürich) The Helmholtz Equation October 28, / 34
30 Existence Theorem the wave number k has influence on the BIE, namely V k is injective on H 1 2 (Γ) iff k 2 is not an eigenvalue of of the interior Dirichlet problem, i.e. u = k 2 u on Ω, γ D u = 0 u = 0 on Ω The kernel of V k is given by ker(v k ) = span [ γ N v : v = k 2 v on Ω and γ D v = 0 on Γ ] Rene Rühr (ETH Zürich) The Helmholtz Equation October 28, / 34
31 Existence Theorem to finalize, we also need the Fredholm Alternative Let A B(X, Y ) coercive. If Au = 0 only allows the trivial solution u = 0, then Au = f is uniquely solvable for all f Y. Else Au = f is solvable iff v, f = 0 for all v Y : A v = 0. Rene Rühr (ETH Zürich) The Helmholtz Equation October 28, / 34
32 Existence so at last Theorem Let Ω be a bounded Lipschitz domain with boundary Γ. Then for every g H 1 2 (Γ) the exterior Dirchlet problem u k 2 u = 0 on Ω c γ c D u = g on Γ has a unique solution u H 1 loc (Ωc ) that satisfies the Sommerfeld radiation condition. Rene Rühr (ETH Zürich) The Helmholtz Equation October 28, / 34
33 Existence need to solve V k φ = ( Id /2 + K k )g on Γ V k coercive, so apply Fredholm! k 2 not an eigenvalue Fredholm V k invertible Else: need to check solvability condition w, ( Id /2 + K k )g L 2 (Γ) = 0 for all w ker(v k ) What s ker(vk )? Can modifiy thm about ker(v k)! [ ] ker(v ) = span γ N v : v = k 2 v on Ω, γ D v = 0 on Γ Rene Rühr (ETH Zürich) The Helmholtz Equation October 28, / 34
34 Existence so take such a γ N v and apply Green s Identity γ N v, ( 1 2 I + K k)g L 2 (Γ) = γ N v, γ D (Ψ DL g) L 2 (Γ) = γ D v, γ N (Ψ DL g) L 2 (Γ) ( k 2 )v, Ψ DL g L 2 (Ω) + v, ( k 2 )Ψ DL g L 2 (Ω) = 0...since γ D v = 0 on Γ, ( k 2 )v = 0 and ( k 2 )Ψ DL g = 0 on Ω. Rene Rühr (ETH Zürich) The Helmholtz Equation October 28, / 34
Integral Representation Formula, Boundary Integral Operators and Calderón projection
Integral Representation Formula, Boundary Integral Operators and Calderón projection Seminar BEM on Wave Scattering Franziska Weber ETH Zürich October 22, 2010 Outline Integral Representation Formula Newton
More informationLECTURE 5 APPLICATIONS OF BDIE METHOD: ACOUSTIC SCATTERING BY INHOMOGENEOUS ANISOTROPIC OBSTACLES DAVID NATROSHVILI
LECTURE 5 APPLICATIONS OF BDIE METHOD: ACOUSTIC SCATTERING BY INHOMOGENEOUS ANISOTROPIC OBSTACLES DAVID NATROSHVILI Georgian Technical University Tbilisi, GEORGIA 0-0 1. Formulation of the corresponding
More informationA coupled BEM and FEM for the interior transmission problem
A coupled BEM and FEM for the interior transmission problem George C. Hsiao, Liwei Xu, Fengshan Liu, Jiguang Sun Abstract The interior transmission problem (ITP) is a boundary value problem arising in
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More information44 CHAPTER 3. CAVITY SCATTERING
44 CHAPTER 3. CAVITY SCATTERING For the TE polarization case, the incident wave has the electric field parallel to the x 3 -axis, which is also the infinite axis of the aperture. Since both the incident
More informationThe Factorization Method for Inverse Scattering Problems Part I
The Factorization Method for Inverse Scattering Problems Part I Andreas Kirsch Madrid 2011 Department of Mathematics KIT University of the State of Baden-Württemberg and National Large-scale Research Center
More informationCoercivity of high-frequency scattering problems
Coercivity of high-frequency scattering problems Valery Smyshlyaev Department of Mathematics, University College London Joint work with: Euan Spence (Bath), Ilia Kamotski (UCL); Comm Pure Appl Math 2015.
More informationFirst kind boundary integral formulation for the Hodge-Helmholtz equation
First kind boundary integral formulation for the Hodge-Helmholtz equation X.Claeys & R.Hiptmair LJLL UPMC / INRIA Alpines, SAM ETH Zürich Integral equation for low frequency Maxwell Ω = bounded Lipschitz,
More informationPreconditioned space-time boundary element methods for the heat equation
W I S S E N T E C H N I K L E I D E N S C H A F T Preconditioned space-time boundary element methods for the heat equation S. Dohr and O. Steinbach Institut für Numerische Mathematik Space-Time Methods
More informationThe double layer potential
The double layer potential In this project, our goal is to explain how the Dirichlet problem for a linear elliptic partial differential equation can be converted into an integral equation by representing
More informationTechnische Universität Graz
Technische Universität Graz Stability of the Laplace single layer boundary integral operator in Sobolev spaces O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2016/2 Technische
More informationSobolev spaces, Trace theorems and Green s functions.
Sobolev spaces, Trace theorems and Green s functions. Boundary Element Methods for Waves Scattering Numerical Analysis Seminar. Orane Jecker October 21, 2010 Plan Introduction 1 Useful definitions 2 Distributions
More informationIntroduction to the Boundary Element Method
Introduction to the Boundary Element Method Salim Meddahi University of Oviedo, Spain University of Trento, Trento April 27 - May 15, 2015 1 Syllabus The Laplace problem Potential theory: the classical
More informationTHE INTERIOR TRANSMISSION PROBLEM FOR REGIONS WITH CAVITIES
THE INTERIOR TRANSMISSION PROBLEM FOR REGIONS WITH CAVITIES FIORALBA CAKONI, AVI COLTON, AN HOUSSEM HAAR Abstract. We consider the interior transmission problem in the case when the inhomogeneous medium
More informationAnalysis in weighted spaces : preliminary version
Analysis in weighted spaces : preliminary version Frank Pacard To cite this version: Frank Pacard. Analysis in weighted spaces : preliminary version. 3rd cycle. Téhéran (Iran, 2006, pp.75.
More informationMATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Spring 2018 Professor: Jared Speck
MATH 8.52 COURSE NOTES - CLASS MEETING # 6 8.52 Introduction to PDEs, Spring 208 Professor: Jared Speck Class Meeting # 6: Laplace s and Poisson s Equations We will now study the Laplace and Poisson equations
More informationMATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Fall 2011 Professor: Jared Speck
MATH 8.52 COURSE NOTES - CLASS MEETING # 6 8.52 Introduction to PDEs, Fall 20 Professor: Jared Speck Class Meeting # 6: Laplace s and Poisson s Equations We will now study the Laplace and Poisson equations
More informationu xx + u yy = 0. (5.1)
Chapter 5 Laplace Equation The following equation is called Laplace equation in two independent variables x, y: The non-homogeneous problem u xx + u yy =. (5.1) u xx + u yy = F, (5.) where F is a function
More informationA class of domains with fractal boundaries: Functions spaces and numerical methods
A class of domains with fractal boundaries: Functions spaces and numerical methods Yves Achdou joint work with T. Deheuvels and N. Tchou Laboratoire J-L Lions, Université Paris Diderot École Centrale -
More informationBETI for acoustic and electromagnetic scattering
BETI for acoustic and electromagnetic scattering O. Steinbach, M. Windisch Institut für Numerische Mathematik Technische Universität Graz Oberwolfach 18. Februar 2010 FWF-Project: Data-sparse Boundary
More informationON THE EXISTENCE OF TRANSMISSION EIGENVALUES. Andreas Kirsch1
Manuscript submitted to AIMS Journals Volume 3, Number 2, May 29 Website: http://aimsciences.org pp. 1 XX ON THE EXISTENCE OF TRANSMISSION EIGENVALUES Andreas Kirsch1 University of Karlsruhe epartment
More informationTraces, extensions and co-normal derivatives for elliptic systems on Lipschitz domains
Traces, extensions and co-normal derivatives for elliptic systems on Lipschitz domains Sergey E. Mikhailov Brunel University West London, Department of Mathematics, Uxbridge, UB8 3PH, UK J. Math. Analysis
More informationAcoustic scattering : high frequency boundary element methods and unified transform methods
Acoustic scattering : high frequency boundary element methods and unified transform methods Book or Report Section Accepted Version Chandler Wilde, S.N. and Langdon, S. (2015) Acoustic scattering : high
More informationThe Interior Transmission Eigenvalue Problem for Maxwell s Equations
The Interior Transmission Eigenvalue Problem for Maxwell s Equations Andreas Kirsch MSRI 2010 epartment of Mathematics KIT University of the State of Baden-Württemberg and National Large-scale Research
More informationSPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS
SPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS TSOGTGEREL GANTUMUR Abstract. After establishing discrete spectra for a large class of elliptic operators, we present some fundamental spectral properties
More informationSecond Order Elliptic PDE
Second Order Elliptic PDE T. Muthukumar tmk@iitk.ac.in December 16, 2014 Contents 1 A Quick Introduction to PDE 1 2 Classification of Second Order PDE 3 3 Linear Second Order Elliptic Operators 4 4 Periodic
More informationInverse Obstacle Scattering
, Göttingen AIP 2011, Pre-Conference Workshop Texas A&M University, May 2011 Scattering theory Scattering theory is concerned with the effects that obstacles and inhomogenities have on the propagation
More informationTechnische Universität Graz
Technische Universität Graz Modified combined field integral equations for electromagnetic scattering O. Steinbach, M. Windisch Berichte aus dem Institut für Numerische Mathematik Bericht 2007/6 Technische
More informationChapter 1 Mathematical Foundations
Computational Electromagnetics; Chapter 1 1 Chapter 1 Mathematical Foundations 1.1 Maxwell s Equations Electromagnetic phenomena can be described by the electric field E, the electric induction D, the
More informationIntroduction to the boundary element method. A case study: the Helmholtz equation
Introduction to the boundary element method. A case study: the Helmholtz equation Francisco Javier Sayas January 2006 Contents 1 The Helmholtz equation 2 2 Scattering 4 3 Single layer acoustic potentials
More information[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,
269 C, Vese Practice problems [1] Write the differential equation u + u = f(x, y), (x, y) Ω u = 1 (x, y) Ω 1 n + u = x (x, y) Ω 2, Ω = {(x, y) x 2 + y 2 < 1}, Ω 1 = {(x, y) x 2 + y 2 = 1, x 0}, Ω 2 = {(x,
More informationDouble Layer Potentials on Polygons and Pseudodifferential Operators on Lie Groupoids
Double Layer Potentials on Polygons and Pseudodifferential Operators on Lie Groupoids joint work with Hengguang Li Yu Qiao School of Mathematics and Information Science Shaanxi Normal University Xi an,
More informationDarwin and higher order approximations to Maxwell s equations in R 3. Sebastian Bauer Universität Duisburg-Essen
Darwin and higher order approximations to Maxwell s equations in R 3 Sebastian Bauer Universität Duisburg-Essen in close collaboration with the Maxwell group around Dirk Pauly Universität Duisburg-Essen
More informationChapter 5 Fast Multipole Methods
Computational Electromagnetics; Chapter 1 1 Chapter 5 Fast Multipole Methods 5.1 Near-field and far-field expansions Like the panel clustering, the Fast Multipole Method (FMM) is a technique for the fast
More informationWhen all else fails, integrate by parts an overview of new and old variational formulations for linear elliptic PDEs
When all else fails, integrate by parts an overview of new and old variational formulations for linear elliptic PDEs E. A. Spence December 31, 2014 Abstract We give an overview of variational formulations
More informationThe Factorization Method for a Class of Inverse Elliptic Problems
1 The Factorization Method for a Class of Inverse Elliptic Problems Andreas Kirsch Mathematisches Institut II Universität Karlsruhe (TH), Germany email: kirsch@math.uni-karlsruhe.de Version of June 20,
More informationNavier-Stokes equations in thin domains with Navier friction boundary conditions
Navier-Stokes equations in thin domains with Navier friction boundary conditions Luan Thach Hoang Department of Mathematics and Statistics, Texas Tech University www.math.umn.edu/ lhoang/ luan.hoang@ttu.edu
More informationPoint estimates for Green s matrix to boundary value problems for second order elliptic systems in a polyhedral cone
Maz ya, V. G., Roßmann, J.: Estimates for Green s matrix 1 ZAMM Z. angew. Math. Mech. 00 2004 0, 1 30 Maz ya, V. G.; Roßmann, J. Point estimates for Green s matrix to boundary value problems for second
More informationChapter 7: Bounded Operators in Hilbert Spaces
Chapter 7: Bounded Operators in Hilbert Spaces I-Liang Chern Department of Applied Mathematics National Chiao Tung University and Department of Mathematics National Taiwan University Fall, 2013 1 / 84
More informationTHE STOKES SYSTEM R.E. SHOWALTER
THE STOKES SYSTEM R.E. SHOWALTER Contents 1. Stokes System 1 Stokes System 2 2. The Weak Solution of the Stokes System 3 3. The Strong Solution 4 4. The Normal Trace 6 5. The Mixed Problem 7 6. The Navier-Stokes
More informationON THE MATHEMATICAL BASIS OF THE LINEAR SAMPLING METHOD
Georgian Mathematical Journal Volume 10 (2003), Number 3, 411 425 ON THE MATHEMATICAL BASIS OF THE LINEAR SAMPLING METHOD FIORALBA CAKONI AND DAVID COLTON Dedicated to the memory of Professor Victor Kupradze
More informationPerron method for the Dirichlet problem.
Introduzione alle equazioni alle derivate parziali, Laurea Magistrale in Matematica Perron method for the Dirichlet problem. We approach the question of existence of solution u C (Ω) C(Ω) of the Dirichlet
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Sobolev Embedding Theorems Embedding Operators and the Sobolev Embedding Theorem
More informationFact Sheet Functional Analysis
Fact Sheet Functional Analysis Literature: Hackbusch, W.: Theorie und Numerik elliptischer Differentialgleichungen. Teubner, 986. Knabner, P., Angermann, L.: Numerik partieller Differentialgleichungen.
More informationSHARP BOUNDARY TRACE INEQUALITIES. 1. Introduction
SHARP BOUNDARY TRACE INEQUALITIES GILES AUCHMUTY Abstract. This paper describes sharp inequalities for the trace of Sobolev functions on the boundary of a bounded region R N. The inequalities bound (semi-)norms
More informationWhen is the error in the h BEM for solving the Helmholtz equation bounded independently of k?
BIT manuscript No. (will be inserted by the editor) When is the error in the h BEM for solving the Helmholtz equation bounded independently of k? I. G. Graham M. Löhndorf J. M. Melenk E. A. Spence Received:
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationTheory of PDE Homework 2
Theory of PDE Homework 2 Adrienne Sands April 18, 2017 In the following exercises we assume the coefficients of the various PDE are smooth and satisfy the uniform ellipticity condition. R n is always an
More informationVolume and surface integral equations for electromagnetic scattering by a dielectric body
Volume and surface integral equations for electromagnetic scattering by a dielectric body M. Costabel, E. Darrigrand, and E. H. Koné IRMAR, Université de Rennes 1,Campus de Beaulieu, 35042 Rennes, FRANCE
More informationSmall energy regularity for a fractional Ginzburg-Landau system
Small energy regularity for a fractional Ginzburg-Landau system Yannick Sire University Aix-Marseille Work in progress with Vincent Millot (Univ. Paris 7) The fractional Ginzburg-Landau system We are interest
More informationPartial regularity for fully nonlinear PDE
Partial regularity for fully nonlinear PDE Luis Silvestre University of Chicago Joint work with Scott Armstrong and Charles Smart Outline Introduction Intro Review of fully nonlinear elliptic PDE Our result
More informationA new method for the solution of scattering problems
A new method for the solution of scattering problems Thorsten Hohage, Frank Schmidt and Lin Zschiedrich Konrad-Zuse-Zentrum Berlin, hohage@zibde * after February 22: University of Göttingen Abstract We
More informationThe oblique derivative problem for general elliptic systems in Lipschitz domains
M. MITREA The oblique derivative problem for general elliptic systems in Lipschitz domains Let M be a smooth, oriented, connected, compact, boundaryless manifold of real dimension m, and let T M and T
More informationOblique derivative problems for elliptic and parabolic equations, Lecture II
of the for elliptic and parabolic equations, Lecture II Iowa State University July 22, 2011 of the 1 2 of the of the As a preliminary step in our further we now look at a special situation for elliptic.
More informationError analysis and fast solvers for high-frequency scattering problems
Error analysis and fast solvers for high-frequency scattering problems I.G. Graham (University of Bath) Woudschoten October 2014 High freq. problem for the Helmholtz equation Given an object Ω R d, with
More informationShort note on compact operators - Monday 24 th March, Sylvester Eriksson-Bique
Short note on compact operators - Monday 24 th March, 2014 Sylvester Eriksson-Bique 1 Introduction In this note I will give a short outline about the structure theory of compact operators. I restrict attention
More informationYour first day at work MATH 806 (Fall 2015)
Your first day at work MATH 806 (Fall 2015) 1. Let X be a set (with no particular algebraic structure). A function d : X X R is called a metric on X (and then X is called a metric space) when d satisfies
More informationON THE COUPLING OF BEM AND FEM FOR EXTERIOR PROBLEMS FOR THE HELMHOLTZ EQUATION
MATHEMATICS OF COMPUTATION Volume 68, Number 227, Pages 945 953 S 0025-5718(99)01064-9 Article electronically published on February 15, 1999 ON THE COUPLING OF BEM AND FEM FOR EXTERIOR PROBLEMS FOR THE
More informationInverse Scattering Theory
Chapter 1 Inverse Scattering Theory In this introductory chapter we provide an overview of the basic ideas of inverse scattering theory for the special case when the scattering object is an isotropic inhomogeneous
More informationBoundary problems for fractional Laplacians
Boundary problems for fractional Laplacians Gerd Grubb Copenhagen University Spectral Theory Workshop University of Kent April 14 17, 2014 Introduction The fractional Laplacian ( ) a, 0 < a < 1, has attracted
More informationAn inverse source problem in optical molecular imaging
An inverse source problem in optical molecular imaging Plamen Stefanov 1 Gunther Uhlmann 2 1 2 University of Washington Formulation Direct Problem Singular Operators Inverse Problem Proof Conclusion Figure:
More informationTechnische Universität Graz
Technische Universität Graz A note on the stable coupling of finite and boundary elements O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2009/4 Technische Universität Graz A
More informationTransmission Eigenvalues in Inverse Scattering Theory
Transmission Eigenvalues in Inverse Scattering Theory David Colton Department of Mathematical Sciences University of Delaware Newark, DE 19716, USA email: colton@math.udel.edu Research supported by a grant
More informationWeak Formulation of Elliptic BVP s
Weak Formulation of Elliptic BVP s There are a large number of problems of physical interest that can be formulated in the abstract setting in which the Lax-Milgram lemma is applied to an equation expressed
More informationGeometric bounds for Steklov eigenvalues
Geometric bounds for Steklov eigenvalues Luigi Provenzano École Polytechnique Fédérale de Lausanne, Switzerland Joint work with Joachim Stubbe June 20, 2017 luigi.provenzano@epfl.ch (EPFL) Steklov eigenvalues
More informationMixed exterior Laplace s problem
Mixed exterior Laplace s problem Chérif Amrouche, Florian Bonzom Laboratoire de mathématiques appliquées, CNRS UMR 5142, Université de Pau et des Pays de l Adour, IPRA, Avenue de l Université, 64000 Pau
More informationWave number explicit bounds in timeharmonic
Wave number explicit bounds in timeharmonic scattering Article Accepted Version Chandler Wilde, S.N. and Monk, P. 008) Wave numberexplicit bounds in time harmonic scattering. SIAM Journal on Mathematical
More informationINSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES. L q -solution of the Neumann, Robin and transmission problem for the scalar Oseen equation
INSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES L q -solution of the Neumann, Robin and transmission problem for the scalar Oseen equation Dagmar Medková Preprint No. 24-2017 PRAHA 2017 L q -SOLUTION
More informationA GALERKIN BOUNDARY ELEMENT METHOD FOR HIGH FREQUENCY SCATTERING BY CONVEX POLYGONS
A GALERKIN BOUNDARY ELEMENT METHOD FOR HIGH FREQUENCY SCATTERING BY CONVEX POLYGONS S N CHANDLER-WILDE AND S LANGDON Abstract. In this paper we consider the problem of time-harmonic acoustic scattering
More informationTransmission Eigenvalues in Inverse Scattering Theory
Transmission Eigenvalues in Inverse Scattering Theory Fioralba Cakoni Department of Mathematical Sciences University of Delaware Newark, DE 19716, USA email: cakoni@math.udel.edu Jointly with D. Colton,
More informationInverse wave scattering problems: fast algorithms, resonance and applications
Inverse wave scattering problems: fast algorithms, resonance and applications Wagner B. Muniz Department of Mathematics Federal University of Santa Catarina (UFSC) w.b.muniz@ufsc.br III Colóquio de Matemática
More informationSome Old and Some New Results in Inverse Obstacle Scattering
Some Old and Some New Results in Inverse Obstacle Scattering Rainer Kress Abstract We will survey on uniqueness, that is, identifiability and on reconstruction issues for inverse obstacle scattering for
More informationThe Gauss-Green Formula (And Elliptic Boundary Problems On Rough Domains) Joint Work with Steve Hofmann and Marius Mitrea
The Gauss-Green Formula And Elliptic Boundary Problems On Rough Domains) Joint Work with Steve Hofmann and Marius Mitrea Dirichlet Problem on open in a compact Riemannian manifold M), dimension n: ) Lu
More informationInverse Scattering Theory: Transmission Eigenvalues and Non-destructive Testing
Inverse Scattering Theory: Transmission Eigenvalues and Non-destructive Testing Isaac Harris Texas A & M University College Station, Texas 77843-3368 iharris@math.tamu.edu Joint work with: F. Cakoni, H.
More informationRegularity Theory a Fourth Order PDE with Delta Right Hand Side
Regularity Theory a Fourth Order PDE with Delta Right Hand Side Graham Hobbs Applied PDEs Seminar, 29th October 2013 Contents Problem and Weak Formulation Example - The Biharmonic Problem Regularity Theory
More informationClass Meeting # 9: Poisson s Formula, Harnack s Inequality, and Liouville s Theorem
MATH 8.52 COUSE NOTES - CLASS MEETING # 9 8.52 Introduction to PDEs, Fall 20 Professor: Jared Speck Class Meeting # 9: Poisson s Formula, Harnack s Inequality, and Liouville s Theorem. epresentation Formula
More informationMatrix construction: Singular integral contributions
Matrix construction: Singular integral contributions Seminar Boundary Element Methods for Wave Scattering Sophie Haug ETH Zurich November 2010 Outline 1 General concepts in singular integral computation
More informationarxiv: v3 [math.ap] 4 Jan 2017
Recovery of an embedded obstacle and its surrounding medium by formally-determined scattering data Hongyu Liu 1 and Xiaodong Liu arxiv:1610.05836v3 [math.ap] 4 Jan 017 1 Department of Mathematics, Hong
More informationFactorization method in inverse
Title: Name: Affil./Addr.: Factorization method in inverse scattering Armin Lechleiter University of Bremen Zentrum für Technomathematik Bibliothekstr. 1 28359 Bremen Germany Phone: +49 (421) 218-63891
More informationHeat Transfer in a Medium in Which Many Small Particles Are Embedded
Math. Model. Nat. Phenom. Vol. 8, No., 23, pp. 93 99 DOI:.5/mmnp/2384 Heat Transfer in a Medium in Which Many Small Particles Are Embedded A. G. Ramm Department of Mathematics Kansas State University,
More informationSTEKLOV REPRESENTATIONS OF GREEN S FUNCTIONS FOR LAPLACIAN BOUNDARY VALUE PROBLEMS. 1. Introduction
STEKLOV REPRESENTATIONS OF GREEN S FUNCTIONS FOR LAPLACIAN BOUNDARY VALUE PROBLEMS. GILES AUCHMUTY Abstract. This paper describes different representations for solution operators of Laplacian boundary
More informationi=1 α i. Given an m-times continuously
1 Fundamentals 1.1 Classification and characteristics Let Ω R d, d N, d 2, be an open set and α = (α 1,, α d ) T N d 0, N 0 := N {0}, a multiindex with α := d i=1 α i. Given an m-times continuously differentiable
More informationStability of an abstract wave equation with delay and a Kelvin Voigt damping
Stability of an abstract wave equation with delay and a Kelvin Voigt damping University of Monastir/UPSAY/LMV-UVSQ Joint work with Serge Nicaise and Cristina Pignotti Outline 1 Problem The idea Stability
More information6.3 Fundamental solutions in d
6.3. Fundamental solutions in d 5 6.3 Fundamental solutions in d Since Lapalce equation is invariant under translations, and rotations (see Exercise 6.4), we look for solutions to Laplace equation having
More informationAn introduction to the use of integral equation. for Harmonic Maxwell equation. J.C. NÉDÉLEC CMAP, École Polytechnique
An introduction to the use of integral equation for Harmonic Maxwell equation J.C. NÉDÉLEC CMAP, École Polytechnique 1. Introduction Sommaire 2. The plane waves solutions 3. Fundamental Solution and Radiation
More informationTrefftz-DG solution to the Helmholtz equation involving integral equations
Trefftz-DG solution to the Helmholtz equation involving integral equations H. Barucq, A. Bendali, M. Fares, V. Mattesi, S. Tordeux Magique 3D Inria Bordeaux Sud Ouest LMA UMR CNRS 5142 INSA of Toulouse
More informationClass Meeting # 9: Poisson s Formula, Harnack s Inequality, and Liouville s Theorem
MATH 85 COUSE NOTES - CLASS MEETING # 9 85 Introduction to PDEs, Fall 0 Professor: Jared Speck Class Meeting # 9: Poisson s Formula, Harnack s Inequality, and Liouville s Theorem epresentation Formula
More informationarxiv: v1 [math.ap] 21 Dec 2018
Uniqueness to Inverse Acoustic and Electromagnetic Scattering From Locally Perturbed Rough Surfaces Yu Zhao, Guanghui Hu, Baoqiang Yan arxiv:1812.09009v1 [math.ap] 21 Dec 2018 Abstract In this paper, we
More informationInverse scattering problem from an impedance obstacle
Inverse Inverse scattering problem from an impedance obstacle Department of Mathematics, NCKU 5 th Workshop on Boundary Element Methods, Integral Equations and Related Topics in Taiwan NSYSU, October 4,
More informationc Copyright by Qi Han August, 2012
c Copyright by Qi Han August, 2012 EXTERIOR REGULARIZED HARMONIC AND HARMONIC FUNCTIONS A Dissertation Presented to the Faculty of the Department of Mathematics University of Houston In Partial Fulfillment
More informationOn the Spectrum of Volume Integral Operators in Acoustic Scattering
11 On the Spectrum of Volume Integral Operators in Acoustic Scattering M. Costabel IRMAR, Université de Rennes 1, France; martin.costabel@univ-rennes1.fr 11.1 Volume Integral Equations in Acoustic Scattering
More informationLecture 8: Boundary Integral Equations
CBMS Conference on Fast Direct Solvers Dartmouth College June 23 June 27, 2014 Lecture 8: Boundary Integral Equations Gunnar Martinsson The University of Colorado at Boulder Research support by: Consider
More informationFrequency functions, monotonicity formulas, and the thin obstacle problem
Frequency functions, monotonicity formulas, and the thin obstacle problem IMA - University of Minnesota March 4, 2013 Thank you for the invitation! In this talk we will present an overview of the parabolic
More informationTHE DIRICHLET PROBLEM FOR THE LAPLACE OPERATOR
THE DIRICHLET PROBLEM FOR THE LAPLACE OPERATOR Stefano Meda Università di Milano-Bicocca c Stefano Meda 2013 ii A Francesco Contents I The Dirichlet problem via Perron s method 1 1 The classical Dirichlet
More informationNonlinear Analysis. Global solution curves for boundary value problems, with linear part at resonance
Nonlinear Analysis 71 (29) 2456 2467 Contents lists available at ScienceDirect Nonlinear Analysis journal homepage: www.elsevier.com/locate/na Global solution curves for boundary value problems, with linear
More informationTransmission eigenvalues with artificial background for explicit material index identification
Transmission eigenvalues with artificial background for explicit material index identification Lorenzo Audibert 1,, Lucas Chesnel, Houssem Haddar 1 Department STEP, EDF R&D, 6 quai Watier, 78401, Chatou
More informationFinite Elements. Colin Cotter. February 22, Colin Cotter FEM
Finite Elements February 22, 2019 In the previous sections, we introduced the concept of finite element spaces, which contain certain functions defined on a domain. Finite element spaces are examples of
More informationA modification of the factorization method for scatterers with different physical properties
A modification of the factorization method for scatterers with different physical properties Takashi FURUYA arxiv:1802.05404v2 [math.ap] 25 Oct 2018 Abstract We study an inverse acoustic scattering problem
More informationPartial Differential Equations, 2nd Edition, L.C.Evans The Calculus of Variations
Partial Differential Equations, 2nd Edition, L.C.Evans Chapter 8 The Calculus of Variations Yung-Hsiang Huang 2018.03.25 Notation: denotes a bounded smooth, open subset of R n. All given functions are
More informationA RELATIONSHIP BETWEEN THE DIRICHLET AND REGULARITY PROBLEMS FOR ELLIPTIC EQUATIONS. Zhongwei Shen
A RELATIONSHIP BETWEEN THE DIRICHLET AND REGULARITY PROBLEMS FOR ELLIPTIC EQUATIONS Zhongwei Shen Abstract. Let L = diva be a real, symmetric second order elliptic operator with bounded measurable coefficients.
More information