The Interior Transmission Eigenvalue Problem for Maxwell s Equations

Size: px
Start display at page:

Download "The Interior Transmission Eigenvalue Problem for Maxwell s Equations"

Transcription

1 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 Center of the Helmholtz Association

2 The Interior Transmission Eigenvalue Problem for Maxwell s Equations 2/17 Outline of the talk The irect Scattering Problem - A Tutorial The Interior Transmission Eigenvalue Problem An Example iscreteness of the Spectrum Existence of Eigenvalues

3 The Interior Transmission Eigenvalue Problem for Maxwell s Equations 3/17 The irect Scattering Problem - A Tutorial Formulation of the Scattering Problem: Given: ε = ε 0 ε r (x), µ = µ 0 and σ = 0 (i.e. non-magnetic and non-conducting media) such that ε r 1 has compact support. Given: Incident waves E inc, H inc which satisfy Maxwell system in frequency domain: etermine: curl E inc iωµ 0 H inc = 0 in R 3, curl H inc + iωε 0 E inc = 0 in R 3. Total fields E and H such that curl E iωµ 0 H = 0 in R 3, curl H + iωε 0 ε r E = 0 in R 3. E s = E E inc, H s = H H inc satisfy Silver-Müller radiation condition H s (x) x x E s (x) = O(1/ x ) as x.

4 The Interior Transmission Eigenvalue Problem for Maxwell s Equations 4/17 Formulation of the Scattering Problem Assumption: ε r L (R 3 ) with ε r = 1 outside some bounded domain and ε r c 0 > 0 on. Elimination of E: curl [ ε 1 r curl H ] k 2 H = 0 in R 3 with wave number k = ω ε 0 µ 0. The solution has to be understood in the weak sense, i.e. H H loc (curl, R 3 ) and [ (ε 1 r curl H) curl ψ k 2 H ψ ] dx = 0 R 3 for all ψ H(curl, R 3 ) with compact support. Remember: ψ H(curl, Ω) ψ, curl ψ L 2 (Ω) 3, ψ H loc (curl, R 3 ) ψ B H(curl, B) for all balls B

5 Volume Potential Fundamental solution of scalar Helmholtz equation u + k 2 u = 0: Φ k (x, y) = exp(ik x y ) 4π x y, x y. Theorem: efine volume potential with density a L 2 () 3 by V (x) = curl Φ k (x, y) a(y) dy = x Φ k (x, y) a(y) dy, x R 3. Then V H loc (curl, R 3 ) is the only radiating solution of in the weak sense, i.e. curl curl V k 2 V = curl a R 3 [ curl V curl ψ k 2 V ψ ] dx = for all ψ H(curl, R 3 ) with compact support. a curl ψ dx The Interior Transmission Eigenvalue Problem for Maxwell s Equations 5/17

6 The Interior Transmission Eigenvalue Problem for Maxwell s Equations 6/17 Integro-ifferential Equation Recall: Rewrite as: curl [ ε 1 r curl H ] k 2 H = 0 curl curl H k 2 H = curl [ (1 ε 1 r ) curl H ] }{{} = q Subtraction of curl curl H inc k 2 H inc = 0 yields curl curl H s k 2 H s = curl [ q curl H ], thus H(x) H inc (x) = H s (x) = curl Φ k (x, y) q(y) curl H(y) dy This is an integro-differential equation for H in H(curl, ). efine T k : H(curl, ) H(curl, ) by T k H = H curl Φ k (, y) q(y) curl H(y) dy

7 Fredholm Property Recall: T k H = H curl Φ k (, y) q(y) curl H(y) dy Then T i : H(curl, ) H(curl, ) is isomorphism onto. Indeed, T i H = V is equivalent to H = V + H and H satisfies H = curl Φ i (, y) q(y) [ curl H(y) + curl V (y) ] dy, thus curl curl H + H = curl [ q(curl H + curl V ) ], and in variational form R 3 [ (ε 1 r curl H) curl ψ + H ψ ] dx = for all ψ H(curl, R 3 ). (q curl V ) curl ψ dx Also T k T i is compact in H(curl, ). The Interior Transmission Eigenvalue Problem for Maxwell s Equations 7/17

8 The Interior Transmission Eigenvalue Problem for Maxwell s Equations 8/17 The irect Scattering Problem Thus: T k H = H inc T i H + (T k T i )H = H inc Theorem: (a) Uniqueness existence, continuous dependence (b) Uniqueness, existence, continuous dependence if ε r C 1,α () Asymptotic behavior from H s = curl Φ k(, y) q(y) curl H(y) dy: H s (x) = exp(ik x ) 4π x uniformly with respect to ˆx = x/ x S 2. H (ˆx) + O(1/ x 2 ), x, Far field pattern H satisfies H (ˆx) ˆx = 0 for all ˆx S 2.

9 The Interior Transmission Eigenvalue Problem for Maxwell s Equations 9/17 The Far Field Operator Incident field: H inc (x) = S 2 g(ˆθ) e ik ˆθ x ds(ˆθ), x R 3, with g L 2 t (S2 ) = {g L 2 (S 2 ) 3 : g(ˆθ) ˆθ = 0 on S 2 }. Corresponding far field pattern H defines far field operator F : L 2 t (S2 ) L 2 t (S2 ), g H. F is compact and normal, scattering matrix S = I + ik F is unitary. 8π 2 Theorem: If F is one-to-one and q q 0 > 0 on then the Factorization Method is applicable, i.e. for any fixed p R 3 with p 0 and any z R 3 define φ z L 2 t (S2 ) by φ z (ˆx) = (ˆx p) e ik ˆx z, ˆx S 2. Then: z φ z R ( (F F) 1/4)

10 The Interior Transmission Eigenvalue Problem for Maxwell s Equations 10/17 Interior Transmission Eigenvalue Problem Question: When is F one-to-one? Fg = 0 = H = 0 = H s = 0 outside of = H = H inc outside of = Cauchy data of H and H inc coincide on. Thus H, H inc satisfy: curl [ ε 1 r curl H ] k 2 H = 0 in, curl curl H inc k 2 H inc = 0 in, ν H = ν H inc, ν [ ε 1 r curl H ] = ν curl H inc on. efinition: λ > 0 is interior transmission eigenvalue if there exists non-trivial solution (u, w) H(curl, ) H(curl, ) with curl[p curl u] λ u = 0 in, curl curl w λ w = 0 in, ν u = ν w, ν [p curl u] = ν curl w on.

11 An Example curl [ p curl u ] λ u = 0, curl curl w λ w = 0 in, ν u = ν w, ν [p curl u] = ν curl w on. Example: unit ball p > 0 constant, λ = k 2. Ansatz (ρ = 1/ p): w(x) = α curl [ j n (kr)y n (ˆx) x ], u(x) = β curl [ j n (ρkr)y n (ˆx) x ] Boundary values: ν w = α j n (k) Grad Y n, ν curl w = α [ k j n(k) + j n (k) ] ˆx Grad Y n. Boundary conditions imply: ( j n (k) j n (ρk) k j n(k) + j n (k) p [ (ρk) j n(ρk) + j n (ρk) ] ) ( α β ) = ( 0 0 Asymptotic form of determinant yields existence of inifinitly many real eigenvalues λ j = k 2 j, and they converge to infinity. The Interior Transmission Eigenvalue Problem for Maxwell s Equations 11/17 ).

12 Elimination of w curl [ p curl u ] λ u = 0, curl curl w λ w = 0 in, ν u = ν w, ν [p curl u] = ν curl w on. Set v = w u and q = p 1. Assumption: q q 0 > 0 on curl [ p curl v ] λ v = curl [ q curl w ] in, ν v = 0, ν [p curl v] = ν [q curl w] on. Multiplication with ϕ for any ϕ H 1 () and integration yields: curl [ p curl v ] ϕ dx λ v ϕ dx = curl [ q curl w ] ϕ dx. Thus: div v = 0 in, ν v = 0 on, i.e. v X 0 where X 0 = { v H 0 (curl, ) : v ϕ dx = 0 ϕ H 1 () }. Also: X = { w H(curl, ) : The Interior Transmission Eigenvalue Problem for Maxwell s Equations 12/17 w ϕ dx = 0 ϕ H 1 () }.

13 Elimination of w, cont. Given v X 0, find w X with curl [ q curl w ] = curl [ p curl v ] λ v in, ν [q curl w] = ν [p curl v] on. This defines A λ : X 0 X0 by [ ] A λ v, ψ = curl w curl ψ λ w ψ dx, v, ψ X0. Theorem: (a) Let (u, w) solve interior transmission eigenvalue problem for some λ > 0. Then A λ v = 0 for v = w u. (b) Let v X 0 satisfy A λ v = 0. efine w X according to above. Then there exists φ H 1 () such that (u, w + φ) solves interior transmission eigenvalue problem where u = w + φ v. The Interior Transmission Eigenvalue Problem for Maxwell s Equations 13/17

14 iscreteness of Spectrum [ ] A λ v, ψ = curl w curl ψ λ w ψ dx, v, ψ X0. [ ] (q curl w) curl ψ dx = (p curl v) curl ψ λv ψ dx for all ψ X. Theorem: A λ : X 0 X 0 has the form A λ = A 0 + λc 1 + λ 2 C 2 where A 0 is self-adjoint and coercive and C 1, C 2 are self-adjoint and compact. Furthermore, C 2 is non-negative. Idea of proof: v j w j, j = 1, 2 (recall: p = q + 1) [ A λ v 1, ψ = (p curl w1 ) curl ψ λ w ] 1 ψ dx (q curl w 1 ) curl ψ dx Set u = w v : A 0 v, v = = curl v 2 dx + curl v 2 dx + curl u curl v dx (q curl u) curl u dx curl v 2 L 2 () The Interior Transmission Eigenvalue Problem for Maxwell s Equations 14/17

15 The Interior Transmission Eigenvalue Problem for Maxwell s Equations 15/17 iscreteness of Spectrum, cont. Quadratic eigenvalue problem: Find λ > 0 and v 0 with A 0 v + λc 1 v + λ 2 C 2 v = 0! Rewrite as (set v 1 = A 1/2 0 v): v 1 + λ A 1/2 0 C 1 A 1/2 0 v 1 + λ }{{} = K 1 2 A 1/2 0 } C 2 A 1/2 0 v 1 {{} = 0. = K2 2 With v 2 = λ K 2 v 1 this becomes ( ) ( v1 K1 K = λ 2 K 2 0 v 2 ) ( v1 v 2 ) This is non-self adjoint eigenvalue problem for compact operator. Therefore, spectrum is discrete, eigenvalues (if they exist) tend to infinity.

16 The Interior Transmission Eigenvalue Problem for Maxwell s Equations 16/17 Existence of Eigenvalues Goal: Find λ > 0 such that A λ fails to be one-to-one, i.e. find λ > 0 such that µ = 0 is eigenvalue of A λ. Let µ min (λ) be the smallest eigenvalue of A λ. Then µ min (0) > 0 since A 0 is coercive. Idea (Päivärinta, Sylvester 2008): Find ˆλ > 0 with µ min (ˆλ) 0. Then there exists λ (0, ˆλ) with µ min (λ) = 0. Idea (Cakoni, Gintides 2010): Choose ball B and eigenvalue ˆλ > 0 with eigenfunctions û, ŵ H 1 (B) corresponding to the ball and constant p 0 = 1 + q 0. Set ˆv = ŵ û. Then ˆv X 0 (B) and ˆλˆv = 0. Extend ˆv by zero into, then ˆv X 0 () and Aˆλˆv, ˆv... ˆλˆv, ˆv = 0. Therefore, µ min (ˆλ) 0. Refined arguments: There exist infinitely many real eigenvalues λ j.

17 Remarks Some References: Colton, Kirsch, Päivärinta: Far field patterns for acoustic waves in an inhomogeneous medium, 1989 u + k 2 pu = 0, w + k 2 w = 0 in, u = w, u ν = w ν Päivärinta, Sylvester: Transmission eigenvalues, 2008 on. Cakoni, Gintides, Haddar: The existence of an infinite discrete set of transmission eigenvalues, 2010 Cakoni, Colton, Haddar: The interior transmission problem for regions with cavity, 2010 Cakoni, Kirsch: On the interior transmission eigenvalue problem, 2010 Open problems: u + k 2 pu = 0 in \ B, u = 0 on B Existence of complex eigenvalues? rop assumption p > 1? The Interior Transmission Eigenvalue Problem for Maxwell s Equations 17/17

ON THE EXISTENCE OF TRANSMISSION EIGENVALUES. Andreas Kirsch1

ON 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 information

The Factorization Method for Inverse Scattering Problems Part I

The 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 information

The Factorization Method for Maxwell s Equations

The Factorization Method for Maxwell s Equations The Factorization Method for Maxwell s Equations Andreas Kirsch University of Karlsruhe Department of Mathematics 76128 Karlsruhe Germany December 15, 2004 Abstract: The factorization method can be applied

More information

Transmission Eigenvalues in Inverse Scattering Theory

Transmission 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 information

The Inside-Outside Duality for Scattering Problems by Inhomogeneous Media

The Inside-Outside Duality for Scattering Problems by Inhomogeneous Media The Inside-Outside uality for Scattering Problems by Inhomogeneous Media Andreas Kirsch epartment of Mathematics Karlsruhe Institute of Technology (KIT) 76131 Karlsruhe Germany and Armin Lechleiter Center

More information

Transmission Eigenvalues in Inverse Scattering Theory

Transmission 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 information

NEW RESULTS ON TRANSMISSION EIGENVALUES. Fioralba Cakoni. Drossos Gintides

NEW RESULTS ON TRANSMISSION EIGENVALUES. Fioralba Cakoni. Drossos Gintides Inverse Problems and Imaging Volume 0, No. 0, 0, 0 Web site: http://www.aimsciences.org NEW RESULTS ON TRANSMISSION EIGENVALUES Fioralba Cakoni epartment of Mathematical Sciences University of elaware

More information

A coupled BEM and FEM for the interior transmission problem

A 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 information

Inverse Scattering Theory: Transmission Eigenvalues and Non-destructive Testing

Inverse 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 information

The Imaging of Anisotropic Media in Inverse Electromagnetic Scattering

The Imaging of Anisotropic Media in Inverse Electromagnetic Scattering The Imaging of Anisotropic Media in Inverse Electromagnetic Scattering Fioralba Cakoni Department of Mathematical Sciences University of Delaware Newark, DE 19716, USA email: cakoni@math.udel.edu Research

More information

Discreteness of Transmission Eigenvalues via Upper Triangular Compact Operators

Discreteness of Transmission Eigenvalues via Upper Triangular Compact Operators Discreteness of Transmission Eigenvalues via Upper Triangular Compact Operators John Sylvester Department of Mathematics University of Washington Seattle, Washington 98195 U.S.A. June 3, 2011 This research

More information

Completeness of the generalized transmission eigenstates

Completeness of the generalized transmission eigenstates Completeness of the generalized transmission eigenstates Eemeli Blåsten Department of Mathematics and Statistics, University of Helsinki Research carried out at the Mittag-Leffler Institute Joint work

More information

Homogenization of the Transmission Eigenvalue Problem for a Periodic Media

Homogenization of the Transmission Eigenvalue Problem for a Periodic Media Homogenization of the Transmission Eigenvalue Problem for a Periodic Media Isaac Harris Texas A&M University, Department of Mathematics College Station, Texas 77843-3368 iharris@math.tamu.edu Joint work

More information

Factorization method in inverse

Factorization 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 information

Chapter 1 Mathematical Foundations

Chapter 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 information

THE INTERIOR TRANSMISSION PROBLEM FOR REGIONS WITH CAVITIES

THE 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 information

The Factorization Method for a Class of Inverse Elliptic Problems

The 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 information

Some Old and Some New Results in Inverse Obstacle Scattering

Some 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 information

Transmission eigenvalues with artificial background for explicit material index identification

Transmission 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 information

Inverse Scattering Theory and Transmission Eigenvalues

Inverse Scattering Theory and Transmission Eigenvalues Inverse Scattering Theory and Transmission Eigenvalues David Colton Department of Mathematical Sciences University of Delaware Newark, DE 19716, USA email: colton@udel.edu Research supported a grant from

More information

Transmission eigenvalues and non-destructive testing of anisotropic magnetic materials with voids

Transmission eigenvalues and non-destructive testing of anisotropic magnetic materials with voids Transmission eigenvalues and non-destructive testing of anisotropic magnetic materials with voids Isaac Harris 1, Fioralba Cakoni 1 and Jiguang Sun 2 1 epartment of Mathematical Sciences, University of

More information

STEKLOFF EIGENVALUES AND INVERSE SCATTERING THEORY

STEKLOFF EIGENVALUES AND INVERSE SCATTERING THEORY STEKLOFF EIGENVALUES AND INVERSE SCATTERING THEORY David Colton, Shixu Meng, Peter Monk University of Delaware Fioralba Cakoni Rutgers University Research supported by AFOSR Grant FA 9550-13-1-0199 Scattering

More information

Transmission Eigenvalues in Inverse Scattering Theory

Transmission Eigenvalues in Inverse Scattering Theory Transmission Eigenvalues in Inverse Scattering Theory Fioralba Cakoni, Houssem Haddar To cite this version: Fioralba Cakoni, Houssem Haddar. Transmission Eigenvalues in Inverse Scattering Theory. Gunther

More information

Computation of Maxwell s transmission eigenvalues and its applications in inverse medium problems

Computation of Maxwell s transmission eigenvalues and its applications in inverse medium problems Computation of Maxwell s transmission eigenvalues and its applications in inverse medium problems Jiguang Sun 1 and Liwei Xu 2,3 1 Department of Mathematical Sciences, Michigan Technological University,

More information

Uniqueness in determining refractive indices by formally determined far-field data

Uniqueness in determining refractive indices by formally determined far-field data Applicable Analysis, 2015 Vol. 94, No. 6, 1259 1269, http://dx.doi.org/10.1080/00036811.2014.924215 Uniqueness in determining refractive indices by formally determined far-field data Guanghui Hu a, Jingzhi

More information

Homogenization of the transmission eigenvalue problem for periodic media and application to the inverse problem

Homogenization of the transmission eigenvalue problem for periodic media and application to the inverse problem Homogenization of the transmission eigenvalue problem for periodic media and application to the inverse problem Fioralba Cakoni, Houssem Haddar, Isaac Harris To cite this version: Fioralba Cakoni, Houssem

More information

Estimation of transmission eigenvalues and the index of refraction from Cauchy data

Estimation of transmission eigenvalues and the index of refraction from Cauchy data Estimation of transmission eigenvalues and the index of refraction from Cauchy data Jiguang Sun Abstract Recently the transmission eigenvalue problem has come to play an important role and received a lot

More information

Fioralba Cakoni 1, Houssem Haddar 2, Thi-Phong Nguyen 1 1 Department of Mathematics, Rutgers University, 110 Frelinghuysen Road, 1.

Fioralba Cakoni 1, Houssem Haddar 2, Thi-Phong Nguyen 1 1 Department of Mathematics, Rutgers University, 110 Frelinghuysen Road, 1. New Interior Transmission Problem Applied to a Single Floquet-Bloch Mode Imaging of Local Perturbations in Periodic Media arxiv:1807.11828v1 [math-ph] 31 Jul 2018 Fioralba Cakoni 1, Houssem Haddar 2, Thi-Phong

More information

An eigenvalue method using multiple frequency data for inverse scattering problems

An eigenvalue method using multiple frequency data for inverse scattering problems An eigenvalue method using multiple frequency data for inverse scattering problems Jiguang Sun Abstract Dirichlet and transmission eigenvalues have important applications in qualitative methods in inverse

More information

Notes on Transmission Eigenvalues

Notes on Transmission Eigenvalues Notes on Transmission Eigenvalues Cédric Bellis December 28, 2011 Contents 1 Scattering by inhomogeneous medium 1 2 Inverse scattering via the linear sampling method 2 2.1 Relationship with the solution

More information

The Asymptotic of Transmission Eigenvalues for a Domain with a Thin Coating

The Asymptotic of Transmission Eigenvalues for a Domain with a Thin Coating The Asymptotic of Transmission Eigenvalues for a Domain with a Thin Coating Hanen Boujlida, Houssem Haddar, Moez Khenissi To cite this version: Hanen Boujlida, Houssem Haddar, Moez Khenissi. The Asymptotic

More information

On Approximate Cloaking by Nonsingular Transformation Media

On Approximate Cloaking by Nonsingular Transformation Media On Approximate Cloaking by Nonsingular Transformation Media TING ZHOU MIT Geometric Analysis on Euclidean and Homogeneous Spaces January 9, 2012 Invisibility and Cloaking From Pendry et al s paper J. B.

More information

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

Volume 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 information

The Helmholtz Equation

The Helmholtz Equation The Helmholtz Equation Seminar BEM on Wave Scattering Rene Rühr ETH Zürich October 28, 2010 Outline Steklov-Poincare Operator Helmholtz Equation: From the Wave equation to Radiation condition Uniqueness

More information

arxiv: v1 [math.ap] 21 Dec 2018

arxiv: 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 information

A Direct Method for reconstructing inclusions from Electrostatic Data

A Direct Method for reconstructing inclusions from Electrostatic Data A Direct Method for reconstructing inclusions from Electrostatic Data Isaac Harris Texas A&M University, Department of Mathematics College Station, Texas 77843-3368 iharris@math.tamu.edu Joint work with:

More information

Darwin 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 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 information

LECTURE 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 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 information

The Mathematical Theory of Maxwell s Equations

The Mathematical Theory of Maxwell s Equations The Mathematical Theory of Maxwell s Equations Andreas Kirsch and Frank Hettlich epartment of Mathematics Karlsruhe Institute of Technology (KIT) Karlsruhe, Germany c January 4, 203 2 Contents Introduction

More information

A computer-assisted Band-Gap Proof for 3D Photonic Crystals

A computer-assisted Band-Gap Proof for 3D Photonic Crystals A computer-assisted Band-Gap Proof for 3D Photonic Crystals Henning Behnke Vu Hoang 2 Michael Plum 2 ) TU Clausthal, Institute for Mathematics, 38678 Clausthal, Germany 2) Univ. of Karlsruhe (TH), Faculty

More information

SPECTRAL APPROXIMATION TO A TRANSMISSION EIGENVALUE PROBLEM AND ITS APPLICATIONS TO AN INVERSE PROBLEM

SPECTRAL APPROXIMATION TO A TRANSMISSION EIGENVALUE PROBLEM AND ITS APPLICATIONS TO AN INVERSE PROBLEM SPECTRAL APPROXIMATIO TO A TRASMISSIO EIGEVALUE PROBLEM AD ITS APPLICATIOS TO A IVERSE PROBLEM JIG A JIE SHE 2,3 Abstract. We first develop an efficient spectral-galerkin method and an rigorous error analysis

More information

The KPP minimal speed within large drift in two dimensions

The KPP minimal speed within large drift in two dimensions The KPP minimal speed within large drift in two dimensions Mohammad El Smaily Joint work with Stéphane Kirsch University of British Columbia & Pacific Institute for the Mathematical Sciences Banff, March-2010

More information

ITERATIVE METHODS FOR TRANSMISSION EIGENVALUES

ITERATIVE METHODS FOR TRANSMISSION EIGENVALUES ITERATIVE METHODS FOR TRANSMISSION EIGENVALUES JIGUANG SUN Abstract. Transmission eigenvalues have important applications in inverse scattering theory. They can be used to obtain useful information of

More information

Transmission eigenvalues for the electromagnetic scattering problem in pseudo-chiral media and a practical reconstruction method

Transmission eigenvalues for the electromagnetic scattering problem in pseudo-chiral media and a practical reconstruction method Transmission eigenvalues for the electromagnetic scattering problem in pseudo-chiral media and a practical reconstruction method Tiexiang Li School of Mathematics, Southeast University, Nanjing 289, People

More information

ELECTROMAGNETIC SCATTERING FROM PERTURBED SURFACES. Katharine Ott Advisor: Irina Mitrea Department of Mathematics University of Virginia

ELECTROMAGNETIC SCATTERING FROM PERTURBED SURFACES. Katharine Ott Advisor: Irina Mitrea Department of Mathematics University of Virginia ELECTROMAGNETIC SCATTERING FROM PERTURBED SURFACES Katharine Ott Advisor: Irina Mitrea Department of Mathematics University of Virginia Abstract This paper is concerned with the study of scattering of

More information

Transmission eigenvalues in inverse scattering theory

Transmission eigenvalues in inverse scattering theory Inside Out II MSRI Publications Volume 6, 22 Transmission eigenvalues in inverse scattering theory FIORALBA CAKONI AN HOUSSEM HAAR In the past few years transmission eigenvalues have become an important

More information

Inverse Obstacle Scattering

Inverse 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 information

Kirchhoff, Fresnel, Fraunhofer, Born approximation and more

Kirchhoff, Fresnel, Fraunhofer, Born approximation and more Kirchhoff, Fresnel, Fraunhofer, Born approximation and more Oberseminar, May 2008 Maxwell equations Or: X-ray wave fields X-rays are electromagnetic waves with wave length from 10 nm to 1 pm, i.e., 10

More information

Inverse wave scattering problems: fast algorithms, resonance and applications

Inverse 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 information

Chapter 5 Fast Multipole Methods

Chapter 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 information

ON THE MATHEMATICAL BASIS OF THE LINEAR SAMPLING METHOD

ON 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 information

Factorization Method for Electromagnetic Inverse Scattering from Biperiodic Structures

Factorization Method for Electromagnetic Inverse Scattering from Biperiodic Structures Factorization Method for Electromagnetic Inverse Scattering from Biperiodic Structures Armin Lecheiter, inh Liem Nguyen To cite this version: Armin Lecheiter, inh Liem Nguyen. Factorization Method for

More information

Spectral Properties of Elliptic Operators In previous work we have replaced the strong version of an elliptic boundary value problem

Spectral Properties of Elliptic Operators In previous work we have replaced the strong version of an elliptic boundary value problem Spectral Properties of Elliptic Operators In previous work we have replaced the strong version of an elliptic boundary value problem L u x f x BC u x g x with the weak problem find u V such that B u,v

More information

Coercivity of high-frequency scattering problems

Coercivity 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 information

The inverse transmission eigenvalue problem

The inverse transmission eigenvalue problem The inverse transmission eigenvalue problem Drossos Gintides Department of Mathematics National Technical University of Athens Greece email: dgindi@math.ntua.gr. Transmission Eigenvalue Problem Find (w,

More information

Monotonicity-based inverse scattering

Monotonicity-based inverse scattering Monotonicity-based inverse scattering Bastian von Harrach http://numerical.solutions Institute of Mathematics, Goethe University Frankfurt, Germany (joint work with M. Salo and V. Pohjola, University of

More information

Interpolation in h-version finite element spaces

Interpolation in h-version finite element spaces Interpolation in h-version finite element spaces Thomas Apel Institut für Mathematik und Bauinformatik Fakultät für Bauingenieur- und Vermessungswesen Universität der Bundeswehr München Chemnitzer Seminar

More information

The inverse transmission eigenvalue problem

The inverse transmission eigenvalue problem The inverse transmission eigenvalue problem Drossos Gintides Department of Mathematics National Technical University of Athens Greece email: dgindi@math.ntua.gr. Transmission Eigenvalue Problem Find (w,

More information

Asymptotic Behavior of Waves in a Nonuniform Medium

Asymptotic Behavior of Waves in a Nonuniform Medium Available at http://pvamuedu/aam Appl Appl Math ISSN: 1932-9466 Vol 12, Issue 1 June 217, pp 217 229 Applications Applied Mathematics: An International Journal AAM Asymptotic Behavior of Waves in a Nonuniform

More information

REAL ANALYSIS II HOMEWORK 3. Conway, Page 49

REAL ANALYSIS II HOMEWORK 3. Conway, Page 49 REAL ANALYSIS II HOMEWORK 3 CİHAN BAHRAN Conway, Page 49 3. Let K and k be as in Proposition 4.7 and suppose that k(x, y) k(y, x). Show that K is self-adjoint and if {µ n } are the eigenvalues of K, each

More information

Hilbert Spaces. Hilbert space is a vector space with some extra structure. We start with formal (axiomatic) definition of a vector space.

Hilbert Spaces. Hilbert space is a vector space with some extra structure. We start with formal (axiomatic) definition of a vector space. Hilbert Spaces Hilbert space is a vector space with some extra structure. We start with formal (axiomatic) definition of a vector space. Vector Space. Vector space, ν, over the field of complex numbers,

More information

KPP Pulsating Traveling Fronts within Large Drift

KPP Pulsating Traveling Fronts within Large Drift KPP Pulsating Traveling Fronts within Large Drift Mohammad El Smaily Joint work with Stéphane Kirsch University of British olumbia & Pacific Institute for the Mathematical Sciences September 17, 2009 PIMS

More information

Laplace s Equation. Chapter Mean Value Formulas

Laplace 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 information

arxiv:submit/ [math.ap] 6 Sep 2017

arxiv:submit/ [math.ap] 6 Sep 2017 INVERSE OBSTACLE SCATTERING FOR ELASTIC WAVES IN THREE DIMENSIONS PEIJUN LI AND XIAOKAI YUAN arxiv:submit/1998861 [math.ap] 6 Sep 2017 Abstract. Consider an exterior problem of the three-dimensional elastic

More information

44 CHAPTER 3. CAVITY SCATTERING

44 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 information

On Transformation Optics based Cloaking

On Transformation Optics based Cloaking On Transformation Optics based Cloaking TING ZHOU Northeastern University CIRM Outline Inverse Problem and Invisibility 1 Inverse Problem and Invisibility 2 Singular cloaking medium 3 Regularization and

More information

A new method for the solution of scattering problems

A 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 information

MOSCO CONVERGENCE FOR H(curl) SPACES, HIGHER INTEGRABILITY FOR MAXWELL S EQUATIONS, AND STABILITY IN DIRECT AND INVERSE EM SCATTERING PROBLEMS

MOSCO CONVERGENCE FOR H(curl) SPACES, HIGHER INTEGRABILITY FOR MAXWELL S EQUATIONS, AND STABILITY IN DIRECT AND INVERSE EM SCATTERING PROBLEMS MOSCO CONVERGENCE FOR H(curl) SPACES, HIGHER INTEGRABILITY FOR MAXWELL S EQUATIONS, AND STABILITY IN DIRECT AND INVERSE EM SCATTERING PROBLEMS HONGYU LIU, LUCA RONDI, AND JINGNI XIAO Abstract. This paper

More information

Finite Element Methods for Maxwell Equations

Finite Element Methods for Maxwell Equations CHAPTER 8 Finite Element Methods for Maxwell Equations The Maxwell equations comprise four first-order partial differential equations linking the fundamental electromagnetic quantities, the electric field

More information

Preparatory Material for the European Intensive Program in Bydgoszcz 2011 Analytical and computer assisted methods in mathematical models

Preparatory Material for the European Intensive Program in Bydgoszcz 2011 Analytical and computer assisted methods in mathematical models Preparatory Material for the European Intensive Program in Bydgoszcz 2011 Analytical and computer assisted methods in mathematical models September 4{18 Basics on the Lebesgue integral and the divergence

More information

Factorization Method for Electromagnetic Inverse Scattering from Biperiodic Structures

Factorization Method for Electromagnetic Inverse Scattering from Biperiodic Structures Factorization Method for Electromagnetic Inverse Scattering from Biperiodic Structures Armin Lechleiter inh-liem Nguyen ecember 23, 2012 Abstract We investigate the Factorization method as an analytical

More information

Vector Spaces. Vector space, ν, over the field of complex numbers, C, is a set of elements a, b,..., satisfying the following axioms.

Vector Spaces. Vector space, ν, over the field of complex numbers, C, is a set of elements a, b,..., satisfying the following axioms. Vector Spaces Vector space, ν, over the field of complex numbers, C, is a set of elements a, b,..., satisfying the following axioms. For each two vectors a, b ν there exists a summation procedure: a +

More information

Trefftz-Discontinuous Galerkin Methods for Maxwell s Equations

Trefftz-Discontinuous Galerkin Methods for Maxwell s Equations Trefftz-Discontinuous Galerkin Methods for Maxwell s Equations Ilaria Perugia Faculty of Mathematics, University of Vienna Vienna P D E Ralf Hiptmair (ETH Zürich), Andrea Moiola (University of Reading)

More information

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

MATH 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 information

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

MATH 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 information

The factorization method for a cavity in an inhomogeneous medium

The factorization method for a cavity in an inhomogeneous medium Home Search Collections Journals About Contact us My IOPscience The factorization method for a cavity in an inhomogeneous medium This content has been downloaded from IOPscience. Please scroll down to

More information

An 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 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 information

TRANSPORT IN POROUS MEDIA

TRANSPORT IN POROUS MEDIA 1 TRANSPORT IN POROUS MEDIA G. ALLAIRE CMAP, Ecole Polytechnique 1. Introduction 2. Main result in an unbounded domain 3. Asymptotic expansions with drift 4. Two-scale convergence with drift 5. The case

More information

A modification of the factorization method for scatterers with different physical properties

A 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 information

Looking Back on Inverse Scattering Theory

Looking Back on Inverse Scattering Theory Looking Back on Inverse Scattering Theory David Colton and Rainer Kress History will be kind to me for I intend to write it Abstract Winston Churchill We present an essay on the mathematical development

More information

On the Spectrum of Volume Integral Operators in Acoustic Scattering

On 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 information

Electrostatic Backscattering by Insulating Obstacles

Electrostatic Backscattering by Insulating Obstacles Electrostatic Backscattering by Insulating Obstacles M. Hanke a,, L. Warth a a Institut für Mathematik, Johannes Gutenberg-Universität Mainz, 5599 Mainz, Germany Abstract We introduce and analyze backscatter

More information

Wave equation on manifolds and finite speed of propagation

Wave equation on manifolds and finite speed of propagation Wave equation on manifolds and finite speed of propagation Ethan Y. Jaffe Let M be a Riemannian manifold (without boundary), and let be the (negative of) the Laplace-Beltrami operator. In this note, we

More information

Composite Membranes. 8th August 2006

Composite Membranes. 8th August 2006 Composite Membranes Aldo Lopes UFMG, aldoelopes@hotmail.com Russell E. Howes, Brigham Young University rhowes@byu.edu Cynthia Shepherd, California State University - Northridge Cynthia.Shepherd.87@csun.edu

More information

Finite Element Methods for Optical Device Design

Finite Element Methods for Optical Device Design DFG Research Center Matheon mathematics for key technologies and Zuse Institute Berlin Finite Element Methods for Optical Device Design Frank Schmidt Sven Burger, Roland Klose, Achim Schädle, Lin Zschiedrich

More information

Inverse Scattering Theory

Inverse 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 information

Phys.Let A. 360, N1, (2006),

Phys.Let A. 360, N1, (2006), Phys.Let A. 360, N1, (006), -5. 1 Completeness of the set of scattering amplitudes A.G. Ramm Mathematics epartment, Kansas State University, Manhattan, KS 66506-60, USA ramm@math.ksu.edu Abstract Let f

More information

The eigenvalue problem for a bianisotropic cavity

The eigenvalue problem for a bianisotropic cavity The eigenvalue problem for a bianisotropic cavity Eftychia Argyropoulou, Andreas D.Ioannidis National and Kapodistrian University of Athens, Greece Linnaeus University, Sweden PIERS 2013, Stockholm, Sweden

More information

Time-Varying Systems; Maxwell s Equations

Time-Varying Systems; Maxwell s Equations Time-Varying Systems; Maxwell s Equations 1. Faraday s law in differential form 2. Scalar and vector potentials; the Lorenz condition 3. Ampere s law with displacement current 4. Maxwell s equations 5.

More information

Recursive integral method for transmission eigenvalues

Recursive integral method for transmission eigenvalues Recursive integral method for transmission eigenvalues Ruihao Huang Allan A. Struthers Jiguang Sun Ruming Zhang Abstract Recently, a new eigenvalue problem, called the transmission eigenvalue problem,

More information

Comparison of a One and Two Parameter Family of Transmission Conditions for Maxwell s Equations with Damping

Comparison of a One and Two Parameter Family of Transmission Conditions for Maxwell s Equations with Damping Comparison of a One and Two Parameter Family of Transmission Conditions for Maxwell s Equations with Damping M. El Bouajaji 1, V. Dolean 2, M. J. Gander 3 and S. Lanteri 1 1 Introduction Transmission conditions

More information

Corner Singularities

Corner Singularities Corner Singularities Martin Costabel IRMAR, Université de Rennes 1 Analysis and Numerics of Acoustic and Electromagnetic Problems Linz, 10 15 October 2016 Martin Costabel (Rennes) Corner Singularities

More information

Electromagnetic wave propagation. ELEC 041-Modeling and design of electromagnetic systems

Electromagnetic wave propagation. ELEC 041-Modeling and design of electromagnetic systems Electromagnetic wave propagation ELEC 041-Modeling and design of electromagnetic systems EM wave propagation In general, open problems with a computation domain extending (in theory) to infinity not bounded

More information

Travelling fronts for the thermodiffusive system with arbitrary Lewis numbers

Travelling fronts for the thermodiffusive system with arbitrary Lewis numbers Travelling fronts for the thermodiffusive system with arbitrary Lewis numbers François Hamel Lenya Ryzhik Abstract We consider KPP-type systems in a cylinder with an arbitrary Lewis number (the ratio of

More information

Outgoing wave conditions in photonic crystals and transmission properties at interfaces

Outgoing wave conditions in photonic crystals and transmission properties at interfaces Outgoing wave conditions in photonic crystals and transmission properties at interfaces A. Lamacz, B. Schweizer December 19, 216 Abstract We analyze the propagation of waves in unbounded photonic crystals.

More information

Integral Equations in Electromagnetics

Integral Equations in Electromagnetics Integral Equations in Electromagnetics Massachusetts Institute of Technology 6.635 lecture notes Most integral equations do not have a closed form solution. However, they can often be discretized and solved

More information

New Helmholtz-Weyl decomposition in L r and its applications to the mathematical fluid mechanics

New Helmholtz-Weyl decomposition in L r and its applications to the mathematical fluid mechanics New Helmholtz-Weyl decomposition in L r and its applications to the mathematical fluid mechanics Hideo Kozono Mathematical Institute Tohoku University Sendai 980-8578 Japan Taku Yanagisawa Department of

More information

DISTRIBUTION A: Distribution approved for public release.

DISTRIBUTION A: Distribution approved for public release. AFRL-AFOSR-VA-TR-2016-0371 Nondestructive Testing and Target Identification David Colton UNIVERSITY OF DELAWARE 12/20/2016 Final Report DISTRIBUTION A: Distribution approved for public release. Air Force

More information

Lecture 8: Boundary Integral Equations

Lecture 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 information

Scattering of electromagnetic waves by thin high contrast dielectrics II: asymptotics of the electric field and a method for inversion.

Scattering of electromagnetic waves by thin high contrast dielectrics II: asymptotics of the electric field and a method for inversion. Scattering of electromagnetic waves by thin high contrast dielectrics II: asymptotics of the electric field and a method for inversion. David M. Ambrose Jay Gopalakrishnan Shari Moskow Scott Rome June

More information