Integral Representation Formula, Boundary Integral Operators and Calderón projection

Similar documents
The Helmholtz Equation

Fact Sheet Functional Analysis

Technische Universität Graz

Sobolev spaces, Trace theorems and Green s functions.

Laplace s Equation. Chapter Mean Value Formulas

Traces, extensions and co-normal derivatives for elliptic systems on Lipschitz domains

Numerical Solutions to Partial Differential Equations

Partial Differential Equations

Preconditioned space-time boundary element methods for the heat equation

Theory of PDE Homework 2

Math The Laplacian. 1 Green s Identities, Fundamental Solution

Regularity Theory a Fourth Order PDE with Delta Right Hand Side

[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,

Introduction to the Boundary Element Method

CONVERGENCE THEORY. G. ALLAIRE CMAP, Ecole Polytechnique. 1. Maximum principle. 2. Oscillating test function. 3. Two-scale convergence

Scientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1

u xx + u yy = 0. (5.1)

When is the error in the h BEM for solving the Helmholtz equation bounded independently of k?

Your first day at work MATH 806 (Fall 2015)

Applied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.

Matrix construction: Singular integral contributions

Sobolev Spaces. Chapter 10

A PROPERTY OF SOBOLEV SPACES ON COMPLETE RIEMANNIAN MANIFOLDS

Stationary mean-field games Diogo A. Gomes

MATH 425, FINAL EXAM SOLUTIONS

LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,

The Factorization Method for Inverse Scattering Problems Part I

LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)

Friedrich symmetric systems

MATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Spring 2018 Professor: Jared Speck

MATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Fall 2011 Professor: Jared Speck

Kersten Schmidt 1 and Ralf Hiptmair 2

NEGATIVE NORM SOBOLEV SPACES AND APPLICATIONS

Green s Functions and Distributions

On some weighted fractional porous media equations

In this chapter we study elliptical PDEs. That is, PDEs of the form. 2 u = lots,

4 Sobolev spaces, trace theorem and normal derivative

Regularity of the obstacle problem for a fractional power of the laplace operator

Short note on compact operators - Monday 24 th March, Sylvester Eriksson-Bique

BOUNDARY VALUE PROBLEMS ON A HALF SIERPINSKI GASKET

Second Order Elliptic PDE

2 Two-Point Boundary Value Problems

Technische Universität Graz

Applications of the periodic unfolding method to multi-scale problems

Chapter 7: Bounded Operators in Hilbert Spaces

THE INVERSE FUNCTION THEOREM

LECTURE 5 APPLICATIONS OF BDIE METHOD: ACOUSTIC SCATTERING BY INHOMOGENEOUS ANISOTROPIC OBSTACLES DAVID NATROSHVILI

Finite Elements. Colin Cotter. February 22, Colin Cotter FEM

The Factorization Method for a Class of Inverse Elliptic Problems

Generalized Lax-Milgram theorem in Banach spaces and its application to the mathematical fluid mechanics.

SOLUTION OF THE DIRICHLET PROBLEM WITH L p BOUNDARY CONDITION. Dagmar Medková

Quantitative Homogenization of Elliptic Operators with Periodic Coefficients

Fundamental Solutions and Green s functions. Simulation Methods in Acoustics

On uniqueness in the inverse conductivity problem with local data

The Mortar Boundary Element Method

Acoustic scattering : high frequency boundary element methods and unified transform methods

Simple Examples on Rectangular Domains

Domain Perturbation for Linear and Semi-Linear Boundary Value Problems

Math 127C, Spring 2006 Final Exam Solutions. x 2 ), g(y 1, y 2 ) = ( y 1 y 2, y1 2 + y2) 2. (g f) (0) = g (f(0))f (0).

Applications of the Maximum Principle

SPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS

Adaptive Boundary Element Methods Part 2: ABEM. Dirk Praetorius

Double Layer Potentials on Polygons and Pseudodifferential Operators on Lie Groupoids

Partial Differential Equations and Sobolev Spaces MAT-INF4300 autumn Snorre H. Christiansen November 10, 2016

Fourier Transform & Sobolev Spaces

Math 5520 Homework 2 Solutions

Divergence Boundary Conditions for Vector Helmholtz Equations with Divergence Constraints

Volume and surface integral equations for electromagnetic scattering by a dielectric body

Méthodes d énergie pour le potentiel de double couche: Une inégalité de Poincaré oubliée.

Finite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product

Partial Differential Equations. William G. Faris

Exercises - Chapter 1 - Chapter 2 (Correction)

arxiv: v1 [math.ap] 25 Jul 2012

Elliptic Problems for Pseudo Differential Equations in a Polyhedral Cone

McGill University Department of Mathematics and Statistics. Ph.D. preliminary examination, PART A. PURE AND APPLIED MATHEMATICS Paper BETA

UNIVERSITY OF MANITOBA

Finite difference method for elliptic problems: I

THE STOKES SYSTEM R.E. SHOWALTER

4 Divergence theorem and its consequences

Chapter 2 Finite Element Spaces for Linear Saddle Point Problems

9. Boundary value problems in a constant-coefficient case

Finite Element Methods for Fourth Order Variational Inequalities

PIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED

Analysis in weighted spaces : preliminary version

The Dirichlet boundary problems for second order parabolic operators satisfying a Carleson condition

ANALYSIS OF A SCALAR PERIDYNAMIC MODEL WITH A SIGN CHANGING KERNEL. Tadele Mengesha. Qiang Du. (Communicated by the associate editor name)

Numerical Methods for Partial Differential Equations

Boundary integral equations on unbounded rough surfaces: Fredholmness and the finite section method

Boundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis

Algebras of singular integral operators with kernels controlled by multiple norms

Lecture No 1 Introduction to Diffusion equations The heat equat

c 2004 Society for Industrial and Applied Mathematics

THE GREEN FUNCTION. Contents

Regularity for Poisson Equation

A class of domains with fractal boundaries: Functions spaces and numerical methods

CONVERGENCE OF EXTERIOR SOLUTIONS TO RADIAL CAUCHY SOLUTIONS FOR 2 t U c 2 U = 0

i=1 α i. Given an m-times continuously

A RELATIONSHIP BETWEEN THE DIRICHLET AND REGULARITY PROBLEMS FOR ELLIPTIC EQUATIONS. Zhongwei Shen

2) Let X be a compact space. Prove that the space C(X) of continuous real-valued functions is a complete metric space.

Traces and Duality Lemma

Transcription:

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 Potential Single Layer Potential and Boundary Integral Operators associated to it Double Layer Potential and Boundary Integral Operators associated to it Calderón projector Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010projection 2 / 45

Newton Potential N 0 V 0 Casting Potentials Boundary Integral Operators Single Layer Potential Ψ 0 SL W 0 K 0 Double Layer Potential Ψ 0 DL K 0 Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010projection 3 / 45

Fundamental Solution of the Laplace Equation Consider the Poisson Equation for d = 2, 3 u(x) = f (x) for x Ω R d, (1) where we assume Ω to be a Lipschitz domain. The fundamental solution of the Laplace operator is given by G 0 (x, y) = { 1 1 4π 2π log x y for d = 2, 1 x y for d = 3. (2) How can this be derived? Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010projection 4 / 45

Fundamental Solution of the Laplace Equation Look for radial solution of homogeneous equation u(x) = 0 for x R d. (3) i.e. we assume the solution is of the form u(x) = v( x ) = v(r). Then for r > 0 which yields r = x i x i r and v(r) = v (r) x i x i r 2 ( 1 v(r) = v (r) x 2 i 2 x 2 i and r x 2 i r 3 r = 1 r x 2 i r 3 ) + v (r) x 2 i r 2. Therefore (3) is equivalent to (assuming u(x) = v( x ) = v(r)) v (d 1) (r) + v (r) = 0. (4) r Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010projection 5 / 45

Fundamental Solution of the Laplace Equation Multiply (4) by r d 1 we get 0 = r d 1 v (r) + r d 2 (d 1)v (r) = (r d 1 v (r)). (5) Integrate twice to obtain (for r > 0) v (r) = c r d 1 v(r) = { a log r + b for d = 2, c r d 2 for d 3. (6) Choosing the constants appropriately, we obtain (2). Notice that u(x) = v( x ) is smooth away from zero and we can write u(x) = δ 0 (x). Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010projection 6 / 45

Integral representation formula Recall Green s second identity (u v u v)dy = (γ D uγ N v γ N uγ D v)ds, (7) Ω where Γ = Ω and γ D, γ N denote the interior Dirichlet and the Neumann traces respectively. Γ Insert G 0 (x, y) for v(y), u(x) = f (x) (Poisson s equation) and use that G 0 (x, y) = δ 0 (x y) in the sense of distributions to obtain u(y)δ 0 (x y) dy u(y)g 0 (x, y) dy = Ω Ω (γ D u(y)γ N,y G 0 (x, y) γ N u(y)γ D G 0 (x, y)) ds(y), x Ω. (8) Γ Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010projection 7 / 45

Integral representation formula So, u(x) = (γ N u(y)γ D G 0 (x, y) γ D u(y)γ N,y G 0 (x, y)) ds(y) Γ + G 0 (x, y)f (y) dy, x Ω, (9) Ω Hence, if the integrals are all well defined and we know the Cauchy data (γ D u(x), γ N u(x)) for x Γ we have found a solution to Poisson s equation (1) in Ω. Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010projection 8 / 45

Integral representation formula Define for a compactly supported function f the volume or Newton potential by (N 0 f )(x) := G 0 (x, y)f (y) dy for x R d (10) R d and the Single Layer Potential for a function φ L 1 (Γ), (Ψ 0 SL φ)(x) := γ D G 0 (x, y)φ(y) ds(y), x R d \Γ. (11) Γ Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010projection 9 / 45

Integral representation formula For φ L 1 (Γ) the Double Layer Potential (Ψ 0 DL φ)(x) := γ N,y G 0 (x, y)φ(y) ds(y), x R d \Γ. (12) Γ For generalized functions we define them by Ψ 0 SL := N 0 γd Ψ0 DL := N 0 γn (13) Rewrite equation (9) as (assuming f is compactly supported in Ω) u(x) = (N 0 f )(x) Ψ 0 DL (γ Du)(x) + Ψ 0 SL (γ Nu)(x) for x R d \Γ. (14) Aim: Show well-definedness, boundedness and useful properties of the operators N 0, Ψ 0 SL and Ψ0 DL and hence derive boundary integral equations for u(x) in order to find the Cauchy data (γ D u, γ N u). Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 10 / 45

The Newton Potential G 0 (x, y) is not integrable in R d for d = 2, 3 and hence N 0 f might not be defined as the integral does not exist for arbitrary f (e.g. consider f 1 and d = 2). But for φ, ψ S(R d ) we have N 0 φ, ψ = ψ(x) G R d 0 (x, y)φ(y)dydx = φ, N 0 ψ (15) R d R d R d using Fubini and that G 0 (x, y) is symmetric, real valued. Hence, N 0 φ S(R d ) and we can define the Newton potential as an operator N 0 : S (R d ) S (R d ) by N 0 f, ψ R d := f, N 0 ψ R d ψ S(R d ) (16) Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 11 / 45

The Newton Potential Theorem We have The Newton Potential is continuous for all s R. N 0 : H s comp(r d ) H s+2 loc (Rd ) (17) Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 12 / 45

The Newton Potential Defining we have Theorem (N 0 f )(x) := Ω G 0 (x, y)f (y) dy for x R d. (18) For s [ 2, 0], the Newton Potential N 0 : H s (Ω) H s+2 (Ω) is a continuous map, i.e. N 0 f H s+2 (Ω) c f Hs (Ω) f H s (Ω) (19) Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 13 / 45

The Newton Potential Proof. 1 Show claim for φ C0 (Ω) and extend it to the completion with respect to the norm on H s (Ω) in the end. 2 Fouriertransform, cut-off functions, splitting of the integral and estimating terms separately, duality arguments Theorem N 0 f is a generalized solution of the partial differential equation { x (N 0 f f (x) for x Ω )(x) = f (x) = 0 for x R d \Ω. (20) Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 14 / 45

Proof. Let φ C0 (Rd ). Use integration by parts, Fubini s theorem and the symmetry of the fundamental solution to obtain ( x (N 0 f ) IBP )(x) φ(x)dx = (N 0 f ( )(x) x φ(x) ) dx R d R d Def. = G 0 (x, y) f (y)dy ( x φ(x) ) dx R d R d sym. = f (y) G 0 (y, x) ( x φ(x) ) dxdy R d R d IBP ( = f (y) x G 0 (y, x) ) φ(x)dxdy R d R d = f (y) δ 0 (x y)φ(x)dxdy R d R d = f (y)φ(y)dy R d Finally use density of C 0 (Rd ) restricted to Ω in H s (Ω). Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 15 / 45

The Newton Potential Applying the interior trace operator γ D (N 0 f )(x) = lim (N 0 f )( x) (21) Ω x x Γ defines a linear bounded operator i.e. it satisfies γ D N 0 : H 1 (Ω) H 1/2 (Γ) (22) γ D N 0 f H 1/2 (Γ) c f H 1 (Ω) f H 1 (Ω) (23) One can show Lemma For f L (Ω) it holds γ D (N 0 f )(x) = Ω G 0 (x, y)f (y)dy (24) Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 16 / 45

The Newton Potential Lemma For the application of the interior Neumann trace, we have The operator γ N N 0 : H 1 (Ω) H 1/2 (Γ) is bounded, i.e. γ N N 0 f H 1/2 (Γ) c f H 1 (Ω) f H 1 (Ω) (25) The proof uses the inverse trace theorem (i.e. that for a given w H 1/2 (Γ) there is a bounded extension Ew H 1 (Ω)), the application of Green s formula and the continuity of N 0 : H 1 (Ω) H 1 (Ω). Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 17 / 45

Outline Integral Representation Formula Newton Potential Single Layer Potential and Boundary Integral Operators associated to it Double Layer Potential and Boundary Integral Operators associated to it Calderón projector Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 18 / 45

The Single Layer Potential We have defined the Single Layer Potential (Ψ 0 SL φ)(x) := γ D G 0 (x, y)φ(y) ds(y), x R d \Γ (26) Theorem for functions φ L 1 (Γ). But the domain of definition can be extended: Γ The Single Layer Potential Ψ 0 SL : H 1/2+s (Γ) H 1+s (Ω) is bounded for s < 1/2, so, Ψ 0 SL φ H 1+s (Ω) c φ H 1/2+s (Γ) φ H 1/2+s (Γ) (27) Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 19 / 45

The Single Layer Potential Lemma Furthermore, The function u(x) := Ψ 0 SL φ(x), x Rd \Γ is a solution of the homogeneous partial differential equation Proof. u(x) = 0 x R d \Γ. (28) Observe that G 0 (x, y) is C for x R d \Γ and y Γ. Therefore we can exchange integration and differentiation to obtain x u(x) = x G 0 (x, y)φ(y)dy Γ = ( x G 0 (x, y))φ(y)dy = 0. Γ Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 20 / 45

The Weakly Singular Boundary Integral Operator V 0 Ψ 0 SL φ H1+s (Ω) for φ H 1/2+s (Γ), s < 1/2, so take the interior trace to obtain the weakly singular boundary integral operator (V 0 φ)(x) := γ D (Ψ 0 SLφ)(x) for x Γ. (29) V 0 : H 1/2+s (Γ) H 1/2+s (Γ) is continuous, i.e. it satisfies V 0 φ H 1/2+s (Γ) c φ H 1/2+s (Γ) φ H 1/2+s (Γ) (30) Lemma For φ L (Γ) we have the representation (V 0 φ)(x) = γ D (Ψ 0 SL φ)(x) = G 0 (x, y)φ(y)ds(y) x Γ (31) as a weakly singular surface integral. Γ Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 21 / 45

The Weakly Singular Boundary Integral Operator V 0 Similarly, we can take the exterior trace (V 0 φ)(x) = γd c (Ψ0 SLφ)(x) := lim Ω c x x Γ (Ψ0 SLφ)( x) for x Γ (32) We get the first jump relation of the single layer potential [γ D Ψ 0 SL φ] := γc D (Ψ0 SL φ)(x) γ D(Ψ 0 SLφ)(x) = 0 x Γ (33) for φ H 1/2+s. This follows from the properties of the trace operator and justifies that we denote the interior as well as the exterior trace of the single layer potential by the same symbol V 0. Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 22 / 45

Theorem Let φ H 1/2 (Γ) for d = 3 and φ H 1/2 (Γ) for d = 2. Then we have for a positive constant c. (H 1/2 (Γ) := { φ H 1/2 (Γ) : φ, 1 Γ = 0 }.) Proof. V 0 φ, φ Γ c φ 2 H 1/2 (Γ) (34) 1 application of interior and exterior Green s formulae a Ω (c)(u, u) = γ (c) N u, γ(c) D Γ u. 2 application of jump relations, triangle and Cauchy-Schwarz inequality 3 a property of the weak solution u H 1 (Ω) of the Dirichlet BVP ( u)(x) = 0, x Ω, γ D u(x) = g(x), x Γ, namely a(u, u) c γ N u 2 H 1/2 (Γ). Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 23 / 45

The Weakly Singular Boundary Integral Operator V 0 Remark: For diam(ω) < 1 one can show that V 0 is H 1/2 (Γ)-elliptic for d = 2 as well, i.e. V 0 φ, φ Γ c φ 2 H 1/2 (Γ) φ H 1/2 (Γ) (35) Notice that the H 1/2 (Γ)-ellipticity together with the H 1/2 (Γ)-continuity by the Lax-Milgram theorem imply that V 0 is invertible! Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 24 / 45

The Adjoint Double Layer Potential K 0 Apply the interior Neumann trace to the single layer potential of a function φ H 1/2+s (Γ) ( s < 1/2). Using (27) together with the properties of the Neumann trace operator γ N, we get a bounded operator γ N Ψ 0 SL : H 1/2+s (Γ) H 1/2+s (Γ) (36) Define the adjoint double layer potential K 0 by (K 0φ)(x) := lim γ N,x G 0 (x, y)φ(y)ds(y), x Γ (37) ɛ 0 and y Γ: y x ɛ 1 1 σ(x) := lim ɛ 0 2(d 1)π ɛ d 1 ds(y), x Γ (38) y Ω: y x =ɛ Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 25 / 45

The Adjoint Double Layer Potential K 0 Lemma Remark: If Γ is differentiable in a neighborhood of x Γ, we have σ(x) = 1 2. For φ H 1/2 (Γ), we have the representation γ N (Ψ 0 SL φ)(x) = σ(x)φ(x) + (K 0φ)(x), for x Γ (39) in the sense of H 1/2 (Γ), i.e. γn (Ψ 0 SL φ), ψ Γ = σφ + K 0φ, ψ Γ, ψ H1/2 (Γ) (40) Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 26 / 45

Proof: 1 Let u := Ψ 0 SL φ. Test with ψ C (Ω) and apply Green s first formula: γn (Ψ 0 SL φ), ψ Γ = x u(x) x ψ(x)dx Ω 2 Insert the definition of the single layer and write as a limit x u(x) x ψ(x)dx = x G 0 (x, y)φ(y)ds(y) x ψ(x)dx Ω Ω Γ ) = x (lim G 0 (x, y)φ(y)ds(y) x ψ(x)dx ɛ 0 Ω y Γ: x y ɛ (41) Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 27 / 45

Interchange the integral signs, integrate by parts and apply Green s formula again to the interior integral: x G 0 (x, y) x ψ(x)dx x Ω: x y ɛ = + x Γ: x y ɛ x Ω: x y =ɛ γ N,x G 0 (x, y)γ D ψ(x)ds(x) γ N,x G 0 (x, y)ψ(x)ds(x) (42) Estimate the obtained summands separately, by explicitely inserting the formula for γ N,x G 0 (x, y) in d = 2, 3. Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 28 / 45

The Adjoint Double Layer Potential K 0 If we go through the proof of the last lemma, we find for the exterior Neumann trace of the single layer potential the following representation: γ c N (Ψ0 SL φ)(x) = (σ(x) 1)φ(x) + (K 0φ)(x), for x Γ (43) in the sense of H 1/2 (Γ). We obtain the second jump relation of the single layer potential: [γ N Ψ 0 SL φ] := γc N (Ψ0 SL φ)(x) γ N(Ψ 0 SLφ)(x) = φ(x), for x Γ, (44) in the sense of H 1/2 (Γ). Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 29 / 45

Outline Integral Representation Formula Newton Potential Single Layer Potential and Boundary Integral Operators associated to it Double Layer Potential and Boundary Integral Operators associated to it Calderón projector Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 30 / 45

The Double Layer Potential The double layer potential (Ψ 0 DL φ)(x) := γ N,y G 0 (x, y)φ(y) ds(y), x R d \Γ (45) Γ defines a continuous operator Ψ 0 DL : H1/2+s (Γ) H 1+s (Ω) for s < 1/2, i.e. Ψ 0 DL φ H 1+s (Ω) c φ H 1/2+s (Γ) φ H 1/2+s (Γ). (46) As in the case of the single layer potential, u H 1 (Ω) defined by u(x) = (Ψ 0 DL φ)(x), x Rd \Γ for ψ H 1/2 (Γ) is a solution to the homogeneous partial differential equation x u(x) = 0, x R d \Γ. (47) Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 31 / 45

The Boundary Integral Operator K 0 Apply the interior trace operator to the double layer potential to obtain a bounded linear operator i.e. γ D Ψ 0 DL : H1/2+s (Γ) H 1/2+s (Γ), (48) γ D Ψ 0 DL φ H 1/2+s (Γ) c φ H 1/2+s (Γ) φ H 1/2+s (Γ) (49) Define for φ H 1/2+s (Γ) the boundary integral operator K 0 by (K 0 φ)(x) := lim (γ N,y G 0 (x, y))φ(y)ds(y). (50) ɛ 0 y Γ: y x ɛ Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 32 / 45

The Boundary Integral Operator K 0 Lemma For φ H 1/2 (Γ) we have the representation γ D (Ψ 0 DL φ)(x) = ( 1 + σ(x))φ(x) + (K 0φ)(x) for x Γ, (51) where σ is defined in (38). For the exterior trace we have the representation γ c D (Ψ0 DL φ)(x) = σ(x)φ(x) + (K 0φ)(x) for x Γ, (52) Hence, we obtain the first jump relation of the double layer potential: [γ D Ψ 0 DL φ] := γc D (Ψ0 DL φ)(x) γ D(Ψ 0 DLφ)(x) = φ(x) for x Γ. (53) Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 33 / 45

The Hypersingular Boundary Integral Operator W 0 For φ H 1/2+s (Γ), s < 1/2, we define for x Γ: (W 0 φ)(x) := γ N (Ψ 0 DLφ)(x) = lim n x x (Ψ 0 DLφ)( x) (54) Ω x x Γ W 0 is the hypersingular boundary integral operator. W 0 : H 1/2+s (Γ) H 1/2+s (Γ) is bounded for s < 1/2, so W 0 φ H 1/2+s (Γ) c φ H 1/2+s (Γ) φ H 1/2+s (Γ). (55) Associated to the hypersingular boundary integral operator is a bilinear form: W 0 v, w Γ = w(x)γ N,x γ N,y G 0 (x, y)v(y)ds(y)ds(x). (56) Γ Γ Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 34 / 45

The Hypersingular Boundary Integral Operator W 0 Insert u 0 (x) = 1 into the representation formula (9) to obtain 1 = γ N,y G 0 ( x, y)ds(y), x Ω. Γ Therefore x (Ψ 0 DL u 0)( x) = 0, x Ω. Taking the limit x x Γ (W 0 u 0 )(x) = 0, x Γ. This implies that we cannot get ellipticity of W 0 on H 1/2 (Γ). So we need to consider a subspace: Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 35 / 45

The Hypersingular Boundary Integral Operator W 0 Consider in d = 2, 3 the problem: Find (t, λ) H 1/2 (Γ) R such that V 0 t, τ Γ λ 1, τ Γ =0 τ H 1/2 (Γ) t, 1 Γ =1. (57) Using the ansatz t := t + 1/ Γ for arbitrary t H 1/2 (Γ), the second equation is already satisfied and the first one reads V0 t, τ Γ = 1 Γ V 01, τ Γ τ H 1/2 (Γ) (58) This variational problem is uniquely solvable which follows from the H 1/2 (Γ)-ellipticity of the operator V 0. The solution w eq := t + 1/ Γ is denoted as the natural density. We have λ = V 0 w eq, w eq Γ. Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 36 / 45

The Hypersingular Boundary Integral Operator W 0 Now we introduce the space H 1/2 (Γ) := { v H 1/2 (Γ) : v, w eq Γ = 0 }. (One can show then that V 0 : H 1/2 (Γ) H 1/2 (Γ) is an isomorphism.) Theorem The hypersingular operator is H 1/2 (Γ)-elliptic, i.e. W 0 φ, φ Γ c φ 2 H 1/2 (Γ) φ H 1/2 (Γ). (59) The proof uses the representation of the trace operators (51), (52), (47), Green s formulae, the jump relations, the norm equivalence theorem of Sobolev and the trace theorem. Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 37 / 45

The Hypersingular Boundary Integral Operator W 0 Finally, we have the second jump relation of the double layer potential: [γ N Ψ 0 DL φ] := γc N (Ψ0 DL φ)(x) γ N(Ψ 0 DLφ)(x) = 0, x Γ. (60) Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 38 / 45

Summary Summing up, we have Theorem Let Ω be a Lipschitz domain, Ω =: Γ its boundary. Then for all s [ 1 2, 1 2] the following boundary operators are bounded: V 0 : H 1/2+s (Γ) H 1/2+s (Γ), (61) K 0 : H 1/2+s (Γ) H 1/2+s (Γ), (62) K 0 : H 1/2+s (Γ) H 1/2+s (Γ), (63) W 0 : H 1/2+s (Γ) H 1/2+s (Γ), (64) and the jump relations for φ H 1/2 (Γ) and φ H 1/2 (Γ) [γ D Ψ 0 SL φ] = 0, [γ DΨ 0 DL ψ] = ψ, in H1/2 (Γ) (65) [γ N Ψ 0 SL φ] = φ, [γ NΨ 0 DL ψ] = 0, in H 1/2 (Γ) (66) Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 39 / 45

Outline Integral Representation Formula Newton Potential Single Layer Potential and Boundary Integral Operators associated to it Double Layer Potential and Boundary Integral Operators associated to it Calderón projector Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 40 / 45

Calderón projector Consider again representation formula (9). u(x) = (N 0 f )(x) Ψ 0 DL (γ Du)(x) + Ψ 0 SL (γ Nu)(x) for x R d \Γ. (67) Take Dirichlet and Neumann trace of the representation formula: γ D u(x) = (V 0 γ N u)(x) + (1 σ(x))γ D u(x) (K 0 γ D u)(x) + γ D N 0 f (x), (68) and γ N u(x) = σ(x)γ N u(x) + (K 0γ N u)(x) + (W 0 γ D u)(x) + γ N N 0 f (x), (69) for x Γ. Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 41 / 45

Calderón projector We have obtained a system of boundary integral equations where ( ) ( γd u 1 = 2 I K ) ( ) 0 V 0 γd u 1 γ N u W 0 2 I + K 0 + γ N u C := is the Calderón projector. ( 1 2 I K ) 0 V 0 1 W 0 2 I + K 0 ( γd N 0 ) f γ N N 0 f (70) (71) Lemma The operator C as defined in (71) is a projection, i.e. C = C 2. Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 42 / 45

Calderón projector Proof. 1 take (v, w) H 1/2 (Γ) H 1/2 (Γ). Then u( x) := (Ψ 0 SL v)( x) (Ψ0 DLw)( x) for x Ω is a solution of the homogeneous partial differential equation. 2 For the traces of u we find: γ D u(x) = (V 0 v)(x) + 1 2 w(x) (K 0w)(x), γ N u(x) = 1 2 v(x) + (K 0v)(x) + (W 0 w)(x). (72) 3 So we have found the Cauchy data (γ D u(x), γ N u(x)), which satify (68) and (69). This is equivalent to (70) with f 0. 4 Insert (72) for (γ D u(x), γ N u(x)) T to obtain the claim. Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 43 / 45

Half past six slide Corollary For the boundary integral operators, we have the relations ( )( ) 1 1 V 0 W 0 = 2 I + K 0 2 I K 0, (73) ( )( ) 1 1 W 0 V 0 = 2 I + K 0 2 I K 0, (74) V 0 K 0 =K 0 V 0, (75) K 0W 0 =W 0 K 0. (76) Since we have found that V 0 and W 0 are invertible operators (if we consider the correct subspaces of H 1/2 (Γ) and H 1/2 (Γ)) this allows us to find the complete Cauchy data if we have either γ D u(x) or γ N u(x) and hence find the solution to the Laplace equation. Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 44 / 45

Questions? Thanks you for listening to me. Aufwachen! Integral Representation Formula, Boundary Integral Operators October and Calderón 22, 2010 projection 45 / 45