A Direct Method for reconstructing inclusions from Electrostatic Data
|
|
- Pauline Elliott
- 5 years ago
- Views:
Transcription
1 A Direct Method for reconstructing inclusions from Electrostatic Data Isaac Harris Texas A&M University, Department of Mathematics College Station, Texas Joint work with: W. Rundell PDE Seminar, University of Houston - April 2017 Reconstructing inclusions from Electrostatic Data 1 / 38
2 Electrical Impedance Tomography Electrical impedance tomography (EIT) is a nondestructive type of imaging method in which the physical parameters of a material is recovered from surface electrode measurements. Figure: Picture taken from Reconstructing inclusions from Electrostatic Data 2 / 38
3 Formulation of the Problem Formulation of the Problem Reconstructing inclusions from Electrostatic Data 3 / 38
4 Formulation of the Problem The Inverse Electrostatic Problem Inverse Problem Reconstruct the support of the inclusion Ω and the boundary condition B from the knowledge of the voltage to current maps Λf = ν u D and Λ 0 f = ν u 0 D. Reconstructing inclusions from Electrostatic Data 4 / 38
5 Formulation of the Problem In this talk we focus on B(u 0 ) = ν u 0 + γ(x)u 0 Assumptions: We consider the Electrostatic Problem in R 2 and R 3 1. D R d be a bounded open region with D is class C 2 2. Ω D such that dist( D, Ω) > 0 with Ω is class C 2 3. The function γ(x) L + ( Ω) is the impedance L + ( Ω) := { γ L ( Ω) : γ(x) 0 } for all x Ω Remark: The results in section can be extended for the case when is replaced by A(x) with symmetric A(x) C 1 (D, C d d ), where the Re(A) is positive definite and Im(A) is non-positive. Reconstructing inclusions from Electrostatic Data 5 / 38
6 Formulation of the Problem Some Previous Works Uniqueness for the inverse problem for 2 Cauchy pairs. V. Bacchelli, Uniqueness for the determination of unknown boundary and impedance with the homogeneous Robin condition, Inverse Problems 25, (2009) An iterative method for one inclusion and γ(x) small. F. Ben Hassen, Y. Boukari and H. Haddar, Inverse impedance boundary problem via the conformal mapping method: the case of small impedances Revue ARIMA, 13 (2010), pp Newton s method for recovering the inclusion and γ(x). W. Rundell, Recovering an obstacle and its impedance from Cauchy data, Inverse Problems 24 (2008), Reconstructing inclusions from Electrostatic Data 6 / 38
7 Formulation of the Problem Iterative Reconstruction Methods Optimization Methods: Solve nonlinear model using iterative scheme to reconstruct the material parameters, which requires a priori information, can be computationally expensive. Now assume that we have f j = u (j) 0 and g j = ν u (j) 0 on D. Λ 0 f j : ( Ω, γ) g j for j = 1, 2 Then you would solve a non-linear and (exponentially) ill-posed posed problem using a iterative scheme. Reconstructing inclusions from Electrostatic Data 7 / 38
8 Formulation of the Problem Figure: Reconstruction of boundary and impedance via an iterative method based on non-linear BIE for the boundary and impedance. W. Rundell, Recovering an obstacle and its impedance from Cauchy data, Inverse Problems 24 (2008), Reconstructing inclusions from Electrostatic Data 8 / 38
9 Formulation of the Problem Restrictions for Iterative Methods Good initial guest is usually needed to insure Convergence This requires a prior information that may not be known: The number of inclusions. Knowledge of the type of Boundary Condition. For the impedance condition a initial estimate is needed. QUESTION: Can one develop an inversion algorithm that requires little to no a prior information? A lack of information cannot be remedied by any mathematical trickery -C. Lanczos Reconstructing inclusions from Electrostatic Data 9 / 38
10 Formulation of the Problem Alternative Reconstruction Method Direct (Qualitative) Methods: Using nonlinear model to reconstruct limited information with no a priori information, Usually by solving many linear integral equations. Our Reconstruction Method We split our inverse problem into 2 parts: 1. Reconstruct Ω via a Sampling Method 2. Once Ω is known reconstruct (non-negative) γ(x) Reconstructing inclusions from Electrostatic Data 10 / 38
11 Formulation of the Problem Useful Monographs for Elliptic IBVP: F. Cakoni and D. Colton, A Qualitative Approach to Inverse Scattering Theory. Applied Mathematical Sciences, Vol 188, Springer, Berlin A. Kirsch and N. Grinberg, The Factorization Method for Inverse Problems. Oxford University Press, Oxford Qualitative Methods in the Time Domain: H. Haddar, A. Lechleiter, S. Marmorat, An improved time domain linear sampling method for Robin and Neumann obstacles, App Analysis, 93 (2014) pp H. Heck, G. Nakamura and H. Wang, Linear sampling method for identifying cavities in a heat conductor, Inverse Problems 28 (2012) Reconstructing inclusions from Electrostatic Data 11 / 38
12 The Linear Sampling Method A Sampling Method to Recover the Boundary Reconstructing inclusions from Electrostatic Data 12 / 38
13 The Linear Sampling Method Let G(x, z) be the solution to the problem G(, z) = δ( z) in D and G(, z) = 0 on D. Now definite operator such that Gh = ν w D where w = 0 in D \ Ω with w D = 0 and ν w Ω = h Using that G(x, z) as x z and the Range of G consists of smooth functions we have that following result. Theorem (Range Test) ν G(, z) D Range(G) if and only if z Ω. Where the operator G : H 1/2 ( Ω) H 1/2 ( D). Reconstructing inclusions from Electrostatic Data 13 / 38
14 The Linear Sampling Method Theorem (Properties of the Data Operator) The difference of the voltage to current operators (Λ Λ 0 ) : H 1/2 ( D) H 1/2 ( D) is compact, injective with a dense range (under some assumptions). Moreover, satisfies that Range(Λ Λ 0 ) Range(G). Techniques used in the proof 1. Derive the factorization (Λ Λ 0 ) = G(L L 0 ) 2. Compactness follows from elliptic regularity 3. Injectivity is a consequence on unique continuation 4. Green s 1st identity gives that (Λ Λ 0 ) is symmetric Reconstructing inclusions from Electrostatic Data 14 / 38
15 The Linear Sampling Method Since (Λ Λ 0 ) has a dense range that there is a sequence such that (Λ Λ 0 )f z,ε ν G(, z) H 0 as ε 0. 1/2 ( D) Assume that f z,ε H 1/2 ( D) is bounded as ε 0. Since f z,ε f z,0 we have that (Λ Λ 0 )f z,ε (Λ Λ 0 )f z,0 as ε 0. Therefore, we have that ν G(, z) D Range(Λ Λ 0 ) Range(G) which is a contradiction provided that z / Ω. Reconstructing inclusions from Electrostatic Data 15 / 38
16 The Linear Sampling Method Determining Ω via Sampling Theorem (The Sampling Method) If the sequence { f z,ε }ε>0 H1/2 ( D) is such that (Λ Λ 0 )f z,ε ν G(, z) H 0 as ε 0 1/2 ( D) then f z,ε H 1/2 ( D) as ε 0 for all z / Ω. Reconstructing the boundary Ω via Sampling: 1. Choose a grid of points in D 2. For each grid point solve (Λ Λ 0 )f z,ε = ν G(, z) D 3. Compute the indicator W (z) = f z,ε 1 H 1/2 ( D) 4. Then Ω δ = { z D : W (z) = δ 1 } approximates Ω Reconstructing inclusions from Electrostatic Data 16 / 38
17 The Linear Sampling Method Numerical Experiments in the Unit Disk 1. Approximate the function u 0 via a System of BIE 2. Approximate the mappings Λ 0 e inθ and Λe inθ = ne inθ 3. Let f z n=0 n=0 f z n e inθ which gives the equation fn z (Λ Λ 0 )e inθ = 1 1 z 2 2π z z cos(θ θ z ) 4. Solve the above equation by a spectral cut-off 5. To visualize the inclusion we let W cutoff (z) = [ 31 n=0 f z n ] 1/2 2 and plot W (z) = W cutoff(z) W cutoff (z) Reconstructing inclusions from Electrostatic Data 17 / 38
18 The Linear Sampling Method Reconstruction of an Ellipse Ellipse: x(θ) = ( 0.5 cos(θ), 0.3 sin(θ) ) 1 Contour Plot Reconstruction 1 Level Curve Reconstruction W(z)=0.005 Reconstruction Inclusion Unit Circle Figure: Reconstruction via the Sampling Method with impedance parameter γ ( x(θ) ) = 2 sin 4 (θ) with cut-off parameter Reconstructing inclusions from Electrostatic Data 18 / 38
19 The Linear Sampling Method Reconstruction of an Cardioid Cardioid: x(θ) = cos(θ) sin(2θ) cos(θ) ( cos(θ), sin(θ) ) Contour Plot Reconstruction Level Curve Reconstruction W(z)=0.005 Reconstruction Inclusion Unit Circle Figure: Reconstruction via the Sampling Method with impedance parameter γ ( x(θ) ) = sin3 (θ) with cut-off parameter Reconstructing inclusions from Electrostatic Data 19 / 38
20 Recovering the BC via Integral Equations A Direct Method to Recover the Impedance Reconstructing inclusions from Electrostatic Data 20 / 38
21 Recovering the BC via Integral Equations Theorem (Uniqueness) The voltage to current mapping Λ 0 : H 1/2 ( D) H 1/2 ( D) uniquely determines the inclusion Ω (provided Ω = Lipshitz) and the impedance parameter γ(x) L + ( Ω). Techniques used in the proof 1. Using the fact that (Λ Λ 0 ) has a dense range gives uniqueness for the inclusion 2. Then prove that the set { } u Ω 0 : for all f H 1/2 ( D) is a dense in L 2 ( Ω) which gives uniqueness for the impedance Reconstructing inclusions from Electrostatic Data 21 / 38
22 Recovering the BC via Integral Equations Data Completion and Recovery of γ(x) Let D = Γ m (measurements) and Ω = Γ i (impedance) both be class C 2 boundaries and define the boundary integral operators (D m ϕ)(x) = 2 ϕ(y) ν(y) Φ(x, y) ds y for x R d \ Γ m Γ m and ( D i ψ)(x) = 2 [ ψ(y) ν(y) Φ(x, y) + x d 2] ds y for x R d \ Γ i Γ i where Φ(x, y) is given by Φ(x, y) = 1 2π ln x y in R2 and Φ(x, y) = 1 1 4π x y in R 3. Reconstructing inclusions from Electrostatic Data 22 / 38
23 Recovering the BC via Integral Equations Data Completion and Recovery of γ(x) We make the ansatz that u 0 (x) = (D m ϕ)(x) + ( D i ψ)(x) for x D \ Ω. Now assume that we have (f, g) = (u 0, ν u 0 ) on Γ m. Therefore (I K mm ) ϕ K im ψ = f on Γ m K mi ϕ + (I + K ii ) ψ = u Γi 0 on Γ i. Here the double layer operators are given by K pq ϕ = (D p ϕ)(x) for Kpq ψ = ( D p ψ)(x) for x Γ q with the index p,q =m,i. Reconstructing inclusions from Electrostatic Data 23 / 38
24 Recovering the BC via Integral Equations Data Completion and Recovery of γ(x) Define the operator [ I K A = mm K mi K im I + K ii ] which represents the integral operator associated with the system of BIE on the previous slide. Theorem (Inversion of A) The mapping A is a bounded linear operator from H 1/2 (Γ m ) H 1/2 (Γ i ) to itself and has a bounded inverse. R. Kress, Linear Integral Equations, Springer, New York, 3rd edition, Reconstructing inclusions from Electrostatic Data 24 / 38
25 Recovering the BC via Integral Equations Data Completion and Recovery of γ(x) Using the fact the we know the normal derivative on we have that [ I K g = [T mm Tim ] mm K mi K ] 1 [ ] im f I + K ii u Γi for x Γ m 0 where the hyper singular operators are given by T mm ϕ = ν(x) (D m ϕ)(x) and Tim ψ = ν(x) ( D i ψ)(x) for x Γ m. Solve for the impedance by using the boundary condition ν u 0 + γ(x)u 0 = 0 for x Γ i. Reconstructing inclusions from Electrostatic Data 25 / 38
26 Recovering the BC via Integral Equations Data Completion and Recovery of γ(x) One can consider two methods for solving for the impedance: (Naive) Direct Solution γ(x n ) = νu 0 (x n ) u 0 (x n ) for n = 1,, N with x n = x(θ n ) Ω and Least Squares Method min a,b N ν u 0 (x n ) + γ(x n )u 0 (x n ) n=1 2 with the impedance function given by γ ( x(θ) ) M a 0 + a m cos(mθ) + b m sin(mθ) m=1 Reconstructing inclusions from Electrostatic Data 26 / 38
27 Recovering the BC via Integral Equations Reconstruction of γ ( x(θ) ) = 5 4.8exp ( 10 sin 4 (θ 3.5) ) sin(θ) Perturbed Ellipse: x(θ) = ( cos(θ), sin(θ) ) (6, 7) 6 5 Exact Reconstructed Figure: Reconstruction for a perturbed ellipse from 10 Cauchy pairs via the Directly Solve Method. Reconstructing inclusions from Electrostatic Data 27 / 38
28 Recovering the BC via Integral Equations Reconstruction of γ ( x(θ) ) = 2 sin 4 (θ) Reconstructed Ellipse: x(θ) r(θ) ( cos(θ), sin(θ) ) r(θ) = Trig Polynomial determined by the level curve Exact Reconstructed Figure: Reconstruction for a reconstructed ellipse via the Least Squares Method with 16-Fourier Modes used. Reconstructing inclusions from Electrostatic Data 28 / 38
29 Current/Future Research An Extension to Other Inverse BVP Reconstructing inclusions from Electrostatic Data 29 / 38
30 Current/Future Research The IBVP for Generalized Impedance condition Reconstruct the support of the inclusion Ω and the boundary condition B from the knowledge of the voltage to current maps Λf = ν u D and Λ 0 f = ν u 0 D. Where the Generalized Impedance Boundary Condition (GIBC) B(u 0 ) = ν u 0 d ds η d ds u 0 + γu 0. F. Cakoni, Y. Hu and R. Kress Simultaneous reconstruction of shape and generalized impedance functions in electrostatic imaging Inverse Problems 30 (2014) Reconstructing inclusions from Electrostatic Data 30 / 38
31 Current/Future Research Some Current Work The direct problem is well-posed for η and γ L ( Ω) provided that the real-parts of the coefficients satisfy Re(η) η min > 0 and Re(γ) 0 where as the imaginary-parts satisfy Im(η) 0 and Im(γ) 0. Uniqueness for recovering Ω hold for a C 2 -boundary. Uniqueness holds for real-valued coefficients provided that η C( Ω) and γ L ( Ω). The Sampling Method is valid for recovering Ω. Reconstructing inclusions from Electrostatic Data 31 / 38
32 Current/Future Research Reconstruction of an Inclusion with GIBC Let the coefficients be given by η = 1/5 and γ = 2 + i where we use separation of variables to determine ( u 0 (r, θ) = a 0 + b 0 ln r + an r n + b n r n ) e inθ. n =1 1 Contour Plot Reconstruction 1 Level Curve W(z)=0.01 Reconstruction Reconstruction Inclusion Unit Circle Figure: Reconstruction for a inclusion with a GIBC via Sampling. Reconstructing inclusions from Electrostatic Data 32 / 38
33 Current/Future Research Reconstruction of (real-valued) Coefficients with GIBC Using Data Completion implies that the mapping ( f, Λ0 f ) ( Ω D u 0 (f ), ν u 0 (f )) is known. Multiplying the GIBC by u 0 and integrating given that Ω u 0 ν u 0 ds = Now assume that Ω η du 0 ds 2 + γ u 0 2 ds for all f H 1/2 ( D) η ( x(s) ) N n=1 η n Ψ (1) ( ) n x(s) and γ ( x(s) ) N n=1 γ n Ψ (2) ( ) n x(s) Notice that by taking f m for m = 1,, M the above equality gives a 2N M linear system of equations to recover η and γ. Reconstructing inclusions from Electrostatic Data 33 / 38
34 Current/Future Research The Inverse Scattering Problem The Physical Model: We consider the time harmonic Acoustic scattering in R 3 or Electromagnetic scattering in R 2 (scalar). Here we wish to use acoustic or electro-magnetic waves to recover physical parameters of a material. Reconstructing inclusions from Electrostatic Data 34 / 38
35 Current/Future Research The Direct Scattering Problem We consider the scattering by an inhomogeneous material with an interior impedance boundary. The total field u = u s + u i H 1 loc (R3 \ Ω) satisfies: u + k 2 (1 + q(x) ) u = 0 in R 3 \ Ω ν u + iγ(x)u = 0 on Ω r u s iku s = o ( r 1) as r. Let u i = e ikx ŷ then it is known that the scattering problem is well-posed for that radiating scattered fields u s (x, ŷ). Reconstructing inclusions from Electrostatic Data 35 / 38
36 Current/Future Research The Far Field Pattern The scattered field, has the following asymptotic expansion in R 3 (and similar in 2 dimensions) { ( )} u s (x, ŷ) = eik x 1 u (ˆx, ŷ) + O 4π x x as x. Using Green s Representation theorem one can show that the far field pattern is given by u (ˆx, ŷ) = u s (z, ŷ) ν e ik ˆx z ν u s (z, ŷ)e ik ˆx z ds(z) C where ˆx, ŷ S = {ˆp R 3 : ˆp = 1}. Reconstructing inclusions from Electrostatic Data 36 / 38
37 Current/Future Research Some Future Works Develop a Sampling Method to determine the inclusion. Uniqueness and Stability of determining Ω and γ(x). Derive Integral Equation to determine the Cauchy data on D. Derive a System of BIE to determine γ(x) on Ω. Some useful work for this problem: Y. Boukari and H. Haddar, A convergent data completion algorithm using surface integral equations, Inverse Problems 25, (2015) H. Qin and X. Liu, The linear sampling method for inhomogeneous media and buried objects from far field measurements, App. Num. Math, 105 (2016), pp Reconstructing inclusions from Electrostatic Data 37 / 38
38 Current/Future Research Figure: Questions? Reconstructing inclusions from Electrostatic Data 38 / 38
Reconstructing inclusions from Electrostatic Data
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: W. Rundell Purdue
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 informationAn 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 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 informationHomogenization 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 informationNonlinear Integral Equations for the Inverse Problem in Corrosion Detection from Partial Cauchy Data. Fioralba Cakoni
Nonlinear Integral Equations for the Inverse Problem in Corrosion Detection from Partial Cauchy Data Fioralba Cakoni Department of Mathematical Sciences, University of Delaware email: cakoni@math.udel.edu
More informationElectrostatic Imaging via Conformal Mapping. R. Kress. joint work with I. Akduman, Istanbul and H. Haddar, Paris
Electrostatic Imaging via Conformal Mapping R. Kress Göttingen joint work with I. Akduman, Istanbul and H. Haddar, Paris Or: A new solution method for inverse boundary value problems for the Laplace equation
More informationEstimation 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 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 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 informationThe 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 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 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 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 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 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 informationInverse Transport Problems and Applications. II. Optical Tomography and Clear Layers. Guillaume Bal
Inverse Transport Problems and Applications II. Optical Tomography and Clear Layers Guillaume Bal Department of Applied Physics & Applied Mathematics Columbia University http://www.columbia.edu/ gb23 gb23@columbia.edu
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 informationOrthogonality Sampling for Object Visualization
Orthogonality ampling for Object Visualization Roland Potthast October 31, 2007 Abstract The goal of this paper is to propose a new sampling algorithm denoted as orthogonality sampling for the detection
More informationNotes 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 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 informationReconstructing conductivities with boundary corrected D-bar method
Reconstructing conductivities with boundary corrected D-bar method Janne Tamminen June 24, 2011 Short introduction to EIT The Boundary correction procedure The D-bar method Simulation of measurement data,
More informationSampling methods for low-frequency electromagnetic imaging
Sampling methods for low-frequency electromagnetic imaging Bastian Gebauer gebauer@math.uni-mainz.de Institut für Mathematik, Joh. Gutenberg-Universität Mainz, Germany Joint work with Jin Keun Seo, Yonsei
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 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 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 informationThe 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 informationSurvey of Inverse Problems For Hyperbolic PDEs
Survey of Inverse Problems For Hyperbolic PDEs Rakesh Department of Mathematical Sciences University of Delaware Newark, DE 19716, USA Email: rakesh@math.udel.edu January 14, 2011 1 Problem Formulation
More informationThe 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 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 informationReverse Time Migration for Extended Obstacles: Acoustic Waves
Reverse Time Migration for Extended Obstacles: Acoustic Waves Junqing Chen, Zhiming Chen, Guanghui Huang epartment of Mathematical Sciences, Tsinghua University, Beijing 8, China LSEC, Institute of Computational
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 informationBEHAVIOR OF THE REGULARIZED SAMPLING INVERSE SCATTERING METHOD AT INTERNAL RESONANCE FREQUENCIES
Progress In Electromagnetics Research, PIER 38, 29 45, 2002 BEHAVIOR OF THE REGULARIZED SAMPLING INVERSE SCATTERING METHOD AT INTERNAL RESONANCE FREQUENCIES N. Shelton and K. F. Warnick Department of Electrical
More informationContributors and Resources
for Introduction to Numerical Methods for d-bar Problems Jennifer Mueller, presented by Andreas Stahel Colorado State University Bern University of Applied Sciences, Switzerland Lexington, Kentucky, May
More informationDetermination of thin elastic inclusions from boundary measurements.
Determination of thin elastic inclusions from boundary measurements. Elena Beretta in collaboration with E. Francini, S. Vessella, E. Kim and J. Lee September 7, 2010 E. Beretta (Università di Roma La
More informationImproved near-wall accuracy for solutions of the Helmholtz equation using the boundary element method
Center for Turbulence Research Annual Research Briefs 2006 313 Improved near-wall accuracy for solutions of the Helmholtz equation using the boundary element method By Y. Khalighi AND D. J. Bodony 1. Motivation
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 informationDiscreteness 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 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 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 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 informationUniqueness 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 informationNEW 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 informationWeierstrass Institute for Applied Analysis and Stochastics Direct and Inverse Elastic Scattering Problems for Diffraction Gratings
Weierstrass Institute for Applied Analysis and Stochastics Direct and Inverse Elastic Scattering Problems for Diffraction Gratings Johannes Elschner & Guanghui Hu Mohrenstrasse 39 10117 Berlin Germany
More informationSTEKLOFF 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 informationA probing method for the transmission eigenvalue problem
A probing method for the transmission eigenvalue problem Fang Zeng Jiguang Sun Liwe Xu Abstract In this paper, we consider an integral eigenvalue problem, which is a reformulation of the transmission eigenvalue
More informationStability and instability in inverse problems
Stability and instability in inverse problems Mikhail I. Isaev supervisor: Roman G. Novikov Centre de Mathématiques Appliquées, École Polytechnique November 27, 2013. Plan of the presentation The Gel fand
More informationTransmission 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 informationThe purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and engineering.
Lecture 16 Applications of Conformal Mapping MATH-GA 451.001 Complex Variables The purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and
More informationMonotonicity arguments in electrical impedance tomography
Monotonicity arguments in electrical impedance tomography Bastian Gebauer gebauer@math.uni-mainz.de Institut für Mathematik, Joh. Gutenberg-Universität Mainz, Germany NAM-Kolloquium, Georg-August-Universität
More informationMonotonicity-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 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 informationINTRODUCTION TO PDEs
INTRODUCTION TO PDEs In this course we are interested in the numerical approximation of PDEs using finite difference methods (FDM). We will use some simple prototype boundary value problems (BVP) and initial
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 informationLooking 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 informationElectrostatic 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 informationSome issues on Electrical Impedance Tomography with complex coefficient
Some issues on Electrical Impedance Tomography with complex coefficient Elisa Francini (Università di Firenze) in collaboration with Elena Beretta and Sergio Vessella E. Francini (Università di Firenze)
More informationEnhanced resolution in structured media
Enhanced resolution in structured media Eric Bonnetier w/ H. Ammari (Ecole Polytechnique), and Yves Capdeboscq (Oxford) Outline : 1. An experiment of super resolution 2. Small volume asymptotics 3. Periodicity
More informationA CALDERÓN PROBLEM WITH FREQUENCY-DIFFERENTIAL DATA IN DISPERSIVE MEDIA. has a unique solution. The Dirichlet-to-Neumann map
A CALDERÓN PROBLEM WITH FREQUENCY-DIFFERENTIAL DATA IN DISPERSIVE MEDIA SUNGWHAN KIM AND ALEXANDRU TAMASAN ABSTRACT. We consider the problem of identifying a complex valued coefficient γ(x, ω) in the conductivity
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 informationStability Estimates in the Inverse Transmission Scattering Problem
Stability Estimates in the Inverse Transmission Scattering Problem Michele Di Cristo Dipartimento di Matematica Politecnico di Milano michele.dicristo@polimi.it Abstract We consider the inverse transmission
More informationNEAR-FIELD IMAGING OF THE SURFACE DISPLACEMENT ON AN INFINITE GROUND PLANE. Gang Bao. Junshan Lin. (Communicated by Haomin Zhou)
Volume X, No X, X, X XX Web site: http://wwwaimsciencesorg NEAR-FIELD IMAGING OF THE SURFACE DISPLACEMENT ON AN INFINITE GROUND PLANE Gang Bao Department of Mathematics, Zhejiang University, Hangzhou,
More informationInverse 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 informationAsymptotic 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 informationThe 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 informationReconstruction Scheme for Active Thermography
Reconstruction Scheme for Active Thermography Gen Nakamura gnaka@math.sci.hokudai.ac.jp Department of Mathematics, Hokkaido University, Japan Newton Institute, Cambridge, Sept. 20, 2011 Contents.1.. Important
More informationComplex geometrical optics solutions for Lipschitz conductivities
Rev. Mat. Iberoamericana 19 (2003), 57 72 Complex geometrical optics solutions for Lipschitz conductivities Lassi Päivärinta, Alexander Panchenko and Gunther Uhlmann Abstract We prove the existence of
More informationA local estimate from Radon transform and stability of Inverse EIT with partial data
A local estimate from Radon transform and stability of Inverse EIT with partial data Alberto Ruiz Universidad Autónoma de Madrid U. California, Irvine.June 2012 w/ P. Caro (U. Helsinki) and D. Dos Santos
More informationCOMPLEX SPHERICAL WAVES AND INVERSE PROBLEMS IN UNBOUNDED DOMAINS
COMPLEX SPHERICAL WAVES AND INVERSE PROBLEMS IN UNBOUNDED DOMAINS MIKKO SALO AND JENN-NAN WANG Abstract. This work is motivated by the inverse conductivity problem of identifying an embedded object in
More informationFast shape-reconstruction in electrical impedance tomography
Fast shape-reconstruction in electrical impedance tomography Bastian von Harrach bastian.harrach@uni-wuerzburg.de (joint work with Marcel Ullrich) Institut für Mathematik - IX, Universität Würzburg The
More informationOn spherical-wave scattering by a spherical scatterer and related near-field inverse problems
IMA Journal of Applied Mathematics (2001) 66, 539 549 On spherical-wave scattering by a spherical scatterer and related near-field inverse problems C. ATHANASIADIS Department of Mathematics, University
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 informationRecent progress on the factorization method for electrical impedance tomography
Recent progress on the factorization method for electrical impedance tomography Bastian von Harrach harrach@math.uni-stuttgart.de Chair of Optimization and Inverse Problems, University of Stuttgart, Germany
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 informationModule 7: The Laplace Equation
Module 7: The Laplace Equation In this module, we shall study one of the most important partial differential equations in physics known as the Laplace equation 2 u = 0 in Ω R n, (1) where 2 u := n i=1
More informationThe Factorization method applied to cracks with impedance boundary conditions
The Factorization method applied to cracks with impedance boundary conditions Yosra Boukari, Houssem Haddar To cite this version: Yosra Boukari, Houssem Haddar. The Factorization method applied to cracks
More informationMath 126 Winter 2012: Rough lecture notes
Math 126 Winter 2012: Rough lecture notes Alex Barnett, Dartmouth College March 1, 2012 1 Introduction Numerical mathematics is at the intersection of analysis (devising and proving theorems), computation
More informationPhys.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 informationLecture Introduction
Lecture 1 1.1 Introduction The theory of Partial Differential Equations (PDEs) is central to mathematics, both pure and applied. The main difference between the theory of PDEs and the theory of Ordinary
More informationA note on the MUSIC algorithm for impedance tomography
A note on the MUSIC algorithm for impedance tomography Martin Hanke Institut für Mathematik, Johannes Gutenberg-Universität, 55099 Mainz, Germany E-mail: hanke@math.uni-mainz.de Abstract. We investigate
More informationSplitting methods with boundary corrections
Splitting methods with boundary corrections Alexander Ostermann University of Innsbruck, Austria Joint work with Lukas Einkemmer Verona, April/May 2017 Strang s paper, SIAM J. Numer. Anal., 1968 S (5)
More informationPut Paper Number Here
Proceedings of Proceedings of icipe : rd International Conference on Inverse Problems in Engineering June,, Port Ludlow, WA, USA Put Paper Number Here NEW RECONSTRUCTION ALGORITHMS IN ELECTRICAL IMPEDANCE
More informationA local estimate from Radon transform and stability of Inverse EIT with partial data
A local estimate from Radon transform and stability of Inverse EIT with partial data Alberto Ruiz Universidad Autónoma de Madrid ge Joint work with P. Caro (U. Helsinki) and D. Dos Santos Ferreira (Paris
More informationELLIPTIC PARTIAL DIFFERENTIAL EQUATIONS EXERCISES I (HARMONIC FUNCTIONS)
ELLIPTIC PARTIAL DIFFERENTIAL EQUATIONS EXERCISES I (HARMONIC FUNCTIONS) MATANIA BEN-ARTZI. BOOKS [CH] R. Courant and D. Hilbert, Methods of Mathematical Physics, Vol. Interscience Publ. 962. II, [E] L.
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 informationRecursive 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 informationDetecting Interfaces in a Parabolic-Elliptic Problem
Detecting Interfaces in a Parabolic-Elliptic Problem Bastian Gebauer bastian.gebauer@oeaw.ac.at Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy of Sciences, Linz,
More informationSensor and Simulation Notes. Note 384. November 1995
Sensor and Simulation Notes Note 384 November 1995 Optimization of the Feed Impedance of Impulse Radiating Antennas Part II: TEM Horns and Lens IRAs Everett G. Farr Farr Research 614 Paseo Del Mar NE Albuquerque,
More informationFast Multipole BEM for Structural Acoustics Simulation
Fast Boundary Element Methods in Industrial Applications Fast Multipole BEM for Structural Acoustics Simulation Matthias Fischer and Lothar Gaul Institut A für Mechanik, Universität Stuttgart, Germany
More informationScattering 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 informationOn 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 informationb n x n + b n 1 x n b 1 x + b 0
Math Partial Fractions Stewart 7.4 Integrating basic rational functions. For a function f(x), we have examined several algebraic methods for finding its indefinite integral (antiderivative) F (x) = f(x)
More informationA DIRECT IMAGING ALGORITHM FOR EXTENDED TARGETS
A DIRECT IMAGING ALGORITHM FOR EXTENDED TARGETS SONGMING HOU, KNUT SOLNA, AND HONGKAI ZHAO Abstract. We present a direct imaging algorithm for extended targets. The algorithm is based on a physical factorization
More informationDISTRIBUTION 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 informationTravelling bubbles in a moving boundary problem of Hele-Shaw type
Travelling bubbles in a moving boundary problem of Hele-Shaw type G. Prokert, TU Eindhoven Center for Analysis, Scientific Computing, and Applications (CASA) g.prokert@tue.nl Joint work with M. Günther,
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 informationThe linear sampling method for three-dimensional inverse scattering problems
ANZIAM J. 42 (E) ppc434 C46, 2 C434 The linear sampling method for three-dimensional inverse scattering problems David Colton Klaus Giebermann Peter Monk (Received 7 August 2) Abstract The inverse scattering
More informationIntroduction of Partial Differential Equations and Boundary Value Problems
Introduction of Partial Differential Equations and Boundary Value Problems 2009 Outline Definition Classification Where PDEs come from? Well-posed problem, solutions Initial Conditions and Boundary Conditions
More informationInverse obstacle scattering problems using multifrequency measurements
Inverse obstacle scattering problems using multifrequency measurements Nguyen Trung Thành Inverse Problems Group, RICAM Joint work with Mourad Sini *** Workshop 3 - RICAM special semester 2011 Nov 21-25
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 information