Inverse scattering problem from an impedance obstacle
|
|
- Leslie Mason
- 5 years ago
- Views:
Transcription
1 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, 2014
2 Outline Inverse 1 Direct 2 Inverse 3 Ellipse Peanut Bean
3 Scattering problem Inverse Direct Object : time harmonic acoustic scattering Modelling : Obstacle Exterior boundary value problem for the Helmholtz equation B) A typical boundary value problem: Point Source Approximation for Obstacle Problems p.5/47
4 Scattering problem Inverse Direct Object : time harmonic acoustic scattering Modelling : Obstacle Exterior boundary value problem for the Helmholtz equation B) A typical boundary value problem: Point Source Approximation for Obstacle Problems p.5/47
5 Scattering problem Inverse Direct Object : time harmonic acoustic scattering Modelling : Obstacle Exterior boundary value problem for the Helmholtz equation B) A typical boundary value problem: Point Source Approximation for Obstacle Problems p.5/47
6 Direct problem Inverse Direct Definition 1 Find: u s C 2 (R 2 \ D) C(R 2 \ D) satisfies 1 the Helmholtz equation u s + k 2 u s = 0, in R 2 \ D 2 the impedance boundary condition u ν for the total field u := u i + u s + λu = 0 on D (1) 3 the Sommerfeld radiation condition(src) ( ) lim r r u s ν ikus = 0, r := x, ˆx := x x
7 Inverse Green s representation formula Direct u s (x) = D u s Φ(x, y) (y) ν(y) us (y) Φ(x, y)ds(y), ν(y) x IR2 \ D (2)
8 Solution ansatz Inverse Direct u s (x) = D Φ(x, y) u(y) ν(y) +λ(y)u(y)φ(x, y)ds(y), x IR2 \D (3)
9 Integral operators Inverse Direct Sϕ(x) := 2 D K ϕ(x) := 2 D Φ(x, y)ϕ(y)ds(y) (4) Φ(x, y) ϕ(y)ds(y) (5) ν(y)
10 Inverse Well-posedness of DP Direct Theorem 1 The direct problem has a unique solution given by u s (x) = D Φ(x, y) u(y) + λ(y)u(y)φ(x, y)ds(y), x IR 2 \ ν(y) D (6) where (the total field) u is the (unique) solution to the following boundary integral equation u Ku S(λu) = 2u i, on D (7)
11 Inverse Direct Far field pattern u The far field pattern or the scattering amplitude is given by ( )} u s (x) = {u eik x 1 (ˆx) + O x x x uniformly for all directions ˆx Ω := {x R 2 x = 1}. In our case u (ˆx) = (c 1 < ν(y), ˆx > +c 2 λ(y)) e ik<ˆx,y> u(y)ds(y) (8) where c 1 = 1 i 4 D k π, c 2 = 1+i 4 kπ
12 Inverse Direct Far field pattern u The far field pattern or the scattering amplitude is given by ( )} u s (x) = {u eik x 1 (ˆx) + O x x x uniformly for all directions ˆx Ω := {x R 2 x = 1}. In our case u (ˆx) = (c 1 < ν(y), ˆx > +c 2 λ(y)) e ik<ˆx,y> u(y)ds(y) (8) where c 1 = 1 i 4 D k π, c 2 = 1+i 4 kπ
13 Inverse Summary : Direct Problem Direct The direct problem can be understood as the process of calculating the far-field pattern from an impedance obstacle. Mathematically, it is equivalent to the solving of the system: { u Ku S(λu) = 2u i, on D u (ˆx) = F( D, λ, u), ˆx Ω (9)
14 Inverse Summary : Direct Problem Direct The direct problem can be understood as the process of calculating the far-field pattern from an impedance obstacle. Mathematically, it is equivalent to the solving of the system: { u Ku S(λu) = 2u i, on D u (ˆx) = F( D, λ, u), ˆx Ω (9)
15 Inverse Problem Inverse Definition 2 (IP) Determine both the scatterer D and the impedance λ if the far field pattern u (, d) is known for one incident direction d and one wave number k > 0.
16 Unique solvability Inverse? Uniqueness Not available Existence Not available
17 Unique solvability Inverse? Uniqueness Not available Existence Not available
18 Unique solvability Inverse? Uniqueness Not available Existence Not available
19 Unique solvability Inverse? Uniqueness Not available Existence Not available
20 Unique solvability Inverse? Uniqueness Not available Existence Not available
21 Inverse Comments on Existence Solving the inverse problem means to solve the far field equation F( D, λ, u) = u (10) However (10) is an equation of the first kind The operator F is compact F has no bounded inverse in general This means that equation (10) cannot be resonably solved!
22 Inverse Comments on Existence Solving the inverse problem means to solve the far field equation F( D, λ, u) = u (10) However (10) is an equation of the first kind The operator F is compact F has no bounded inverse in general This means that equation (10) cannot be resonably solved!
23 Inverse Comments on Existence Solving the inverse problem means to solve the far field equation F( D, λ, u) = u (10) However (10) is an equation of the first kind The operator F is compact F has no bounded inverse in general This means that equation (10) cannot be resonably solved!
24 Inverse Fredholm integral Equations 2. Kind ϕ(x) (x+1)e xy ϕ(y)dy = e x e (x+1), 0 x 1 Trapzoidal rule n x = 0 x = 0.5 x = Simpson s rule n x = 0 x = 0.5 x =
25 Inverse Fredholm integral Equations 1. Kind 1 0 (x + 1)e xy ϕ(y)dy = 1 e (x+1), 0 x 1 Trapzoidal rule n x = 0 x = 0.5 x = Simpson s rule n x = 0 x = 0.5 x =
26 Inverse Ill-Posed Problems : Definition 3 () Assume X, Y are normed spaces. Let the operator A : X Y be linear, bounded and injective. A family of bounded linear operators R α : Y X, α > 0 is called a regularization scheme for Aϕ = f, if it satisfies the following pointwise convergence lim R αaϕ = ϕ, for all ϕ X α 0 In this case, the parameter α is called the regularization parameter.
27 Inverse Ill-Posed Problems : Definition 3 () Assume X, Y are normed spaces. Let the operator A : X Y be linear, bounded and injective. A family of bounded linear operators R α : Y X, α > 0 is called a regularization scheme for Aϕ = f, if it satisfies the following pointwise convergence lim R αaϕ = ϕ, for all ϕ X α 0 In this case, the parameter α is called the regularization parameter.
28 Inverse : Error Find a stable approximation to the equation Aϕ = f The regularized approximation ϕ δ α := R α f δ The total approximation error ϕ δ α ϕ = R α f δ R α f + R α Aϕ ϕ We have ϕ δ α ϕ δ R α + R α Aϕ ϕ
29 Inverse : Error Find a stable approximation to the equation Aϕ = f The regularized approximation ϕ δ α := R α f δ The total approximation error ϕ δ α ϕ = R α f δ R α f + R α Aϕ ϕ We have ϕ δ α ϕ δ R α + R α Aϕ ϕ
30 Inverse : Parameter How to choose the regularization parameter α? 1 a priori choice based on some information of the solution. In general not available 2 a posteriori choice based on the data error level δ Discrepancy Principle of Morozov : AR α f δ f δ = γδ, γ 1
31 Inverse : Parameter How to choose the regularization parameter α? 1 a priori choice based on some information of the solution. In general not available 2 a posteriori choice based on the data error level δ Discrepancy Principle of Morozov : AR α f δ f δ = γδ, γ 1
32 Inverse : Parameter How to choose the regularization parameter α? 1 a priori choice based on some information of the solution. In general not available 2 a posteriori choice based on the data error level δ Discrepancy Principle of Morozov : AR α f δ f δ = γδ, γ 1
33 Inverse : Parameter How to choose the regularization parameter α? 1 a priori choice based on some information of the solution. In general not available 2 a posteriori choice based on the data error level δ Discrepancy Principle of Morozov : AR α f δ f δ = γδ, γ 1
34 Inverse : Example X, Y Hilbert spaces. Theorem 2 Assume A : X Y compact and linear. Then for every α > 0, the operator αi + A A : X X is bijective and has a bounded inverse. Furthermore, if the operator A is injective, then R α := (αi + A A) 1 A, α > 0 describes a regularization scheme with R α 1 2 α.
35 Inverse : Example X, Y Hilbert spaces. Theorem 2 Assume A : X Y compact and linear. Then for every α > 0, the operator αi + A A : X X is bijective and has a bounded inverse. Furthermore, if the operator A is injective, then R α := (αi + A A) 1 A, α > 0 describes a regularization scheme with R α 1 2 α.
36 Inverse Tikhonov Theorem 3 Let A : X Y be a linear and bounded operator. Assme α > 0. Then for each f Y there exists a unique ϕ α X such that { Aϕ α f + α ϕ α = inf ϕ X Aϕ f 2 + α ϕ 2} The minimizer ϕ α is given by the unique solution of the equation αϕ α + A Aϕ α = A f and depends continuously on f.
37 Inverse : Newton s method F( D, λ, u) = u 3 F Θ = u F (11) where αi βi γi + A A Θ = A (u F) (12) F D 0 0 A = F 0 λ 0 F 0 0 u
38 Inverse : Newton s method F( D, λ, u) = u 3 F Θ = u F (11) where αi βi γi + A A Θ = A (u F) (12) F D 0 0 A = F 0 λ 0 F 0 0 u
39 Inverse Modified Newton s Method Recall the solution of direct problem (9) { u Ku S(λu) = 2u i, on D u (ˆx) = F (γ, λ, u), ˆx Ω Split the inverse problem into two parts: 1 B(γ, λ)u = 2u i (13) 2 F (u)(γ, λ) = u (14) (13) is solved as a (well-posed) direct problem. (14) is solved as an ill-posed problem with two regularization parameters.
40 Inverse Modified Newton s Method Recall the solution of direct problem (9) { u Ku S(λu) = 2u i, on D u (ˆx) = F (γ, λ, u), ˆx Ω Split the inverse problem into two parts: 1 B(γ, λ)u = 2u i (13) 2 F (u)(γ, λ) = u (14) (13) is solved as a (well-posed) direct problem. (14) is solved as an ill-posed problem with two regularization parameters.
41 Inverse Modified Newton s Method Recall the solution of direct problem (9) { u Ku S(λu) = 2u i, on D u (ˆx) = F (γ, λ, u), ˆx Ω Split the inverse problem into two parts: 1 B(γ, λ)u = 2u i (13) 2 F (u)(γ, λ) = u (14) (13) is solved as a (well-posed) direct problem. (14) is solved as an ill-posed problem with two regularization parameters.
42 Inverse Modified Newton s Method Recall the solution of direct problem (9) { u Ku S(λu) = 2u i, on D u (ˆx) = F (γ, λ, u), ˆx Ω Split the inverse problem into two parts: 1 B(γ, λ)u = 2u i (13) 2 F (u)(γ, λ) = u (14) (13) is solved as a (well-posed) direct problem. (14) is solved as an ill-posed problem with two regularization parameters.
43 Inverse Modified Newton s Method Recall the solution of direct problem (9) { u Ku S(λu) = 2u i, on D u (ˆx) = F (γ, λ, u), ˆx Ω Split the inverse problem into two parts: 1 B(γ, λ)u = 2u i (13) 2 F (u)(γ, λ) = u (14) (13) is solved as a (well-posed) direct problem. (14) is solved as an ill-posed problem with two regularization parameters.
44 Inverse Modified Newton s Method Recall the solution of direct problem (9) { u Ku S(λu) = 2u i, on D u (ˆx) = F (γ, λ, u), ˆx Ω Split the inverse problem into two parts: 1 B(γ, λ)u = 2u i (13) 2 F (u)(γ, λ) = u (14) (13) is solved as a (well-posed) direct problem. (14) is solved as an ill-posed problem with two regularization parameters.
45 Inverse Modified Newton s Method : iterative scheme 1 Given a pair of initial guesses γ 0, λ 0. 2 Solve (13) for u. 3 Solve the regularized version of (14) for updates of γ, λ : ([ ] ) [ ] αi 0 + A χ A = A (u 0 βi q F) (15) where A = [ F γ 0 0 F λ ] 4 Set γ 0 = γ 0 + χ, λ 0 = λ 0 + q and repeat steps 2,3 until some criterion is fullfilled.
46 Inverse Modified Newton s Method : iterative scheme 1 Given a pair of initial guesses γ 0, λ 0. 2 Solve (13) for u. 3 Solve the regularized version of (14) for updates of γ, λ : ([ ] ) [ ] αi 0 + A χ A = A (u 0 βi q F) (15) where A = [ F γ 0 0 F λ ] 4 Set γ 0 = γ 0 + χ, λ 0 = λ 0 + q and repeat steps 2,3 until some criterion is fullfilled.
47 Inverse Modified Newton s Method : iterative scheme 1 Given a pair of initial guesses γ 0, λ 0. 2 Solve (13) for u. 3 Solve the regularized version of (14) for updates of γ, λ : ([ ] ) [ ] αi 0 + A χ A = A (u 0 βi q F) (15) where A = [ F γ 0 0 F λ ] 4 Set γ 0 = γ 0 + χ, λ 0 = λ 0 + q and repeat steps 2,3 until some criterion is fullfilled.
48 Inverse Modified Newton s Method : iterative scheme 1 Given a pair of initial guesses γ 0, λ 0. 2 Solve (13) for u. 3 Solve the regularized version of (14) for updates of γ, λ : ([ ] ) [ ] αi 0 + A χ A = A (u 0 βi q F) (15) where A = [ F γ 0 0 F λ ] 4 Set γ 0 = γ 0 + χ, λ 0 = λ 0 + q and repeat steps 2,3 until some criterion is fullfilled.
49 Inverse Modified Newton s Method : iterative scheme 1 Given a pair of initial guesses γ 0, λ 0. 2 Solve (13) for u. 3 Solve the regularized version of (14) for updates of γ, λ : ([ ] ) [ ] αi 0 + A χ A = A (u 0 βi q F) (15) where A = [ F γ 0 0 F λ ] 4 Set γ 0 = γ 0 + χ, λ 0 = λ 0 + q and repeat steps 2,3 until some criterion is fullfilled.
50 Major Advantages Inverse No extra equations Two smaller systems Fréchet derivatives are easily obtained The unique solvability of the numerical scheme followed straight forward
51 Inverse Numerical Settings: Solution spaces V m := span{1, cos t, cos 2t,..., cos mt; sin t,..., sin(m 1)t} Γ := (γ 1 (t), γ 2 (t)) V m V m λ V k Stopping criterion for the Newton s method: Discrepancy Principle: u,k u ɛ
52 Inverse Numerical Settings: Solution spaces V m := span{1, cos t, cos 2t,..., cos mt; sin t,..., sin(m 1)t} Γ := (γ 1 (t), γ 2 (t)) V m V m λ V k Stopping criterion for the Newton s method: Discrepancy Principle: u,k u ɛ
53 Inverse Numerical Settings: Solution spaces V m := span{1, cos t, cos 2t,..., cos mt; sin t,..., sin(m 1)t} Γ := (γ 1 (t), γ 2 (t)) V m V m λ V k Stopping criterion for the Newton s method: Discrepancy Principle: u,k u ɛ
54 Inverse Ellipse with exact data Ellipse Peanut Bean Γ = (0.4 cos t, 0.3 sin t)
55 Inverse Ellipse with exact data Ellipse Peanut Bean λ = cos t
56 Inverse Ellipse with 3% noise Ellipse Peanut Bean Γ = (0.4 cos t, 0.3 sin t)
57 Inverse Ellipse with 3% noise Ellipse Peanut Bean λ = cos t
58 Peanut A Inverse Ellipse Peanut Bean Γ = γ(t)(cos t, sin t), γ(t) = cos 2 t sin 2 t
59 Peanut A Inverse Ellipse Peanut Bean λ = sin t
60 Peanut B Inverse Ellipse Peanut Bean Γ = γ(t)(cos t, sin t), γ(t) = cos 2 t sin 2 t
61 Peanut B Inverse Ellipse Peanut Bean λ = 0.1 sin 2t cos t
62 Peanut C Inverse Ellipse Peanut Bean Γ = γ(t)(cos t, sin t), γ(t) = cos 2 t sin 2 t
63 Peanut C Inverse Ellipse Peanut Bean λ = 0.3e 0.25 cos3 t
64 Inverse Bean with exact data Ellipse Peanut Bean Γ = γ(t)(cos t, sin t), γ(t) = cos t sin(2t) cos t
65 Inverse Bean with exact data Ellipse Peanut Bean λ = sin t
66 Inverse Bean with 3% noise Ellipse Peanut Bean Γ = γ(t)(cos t, sin t), γ(t) = cos t sin(2t) cos t
67 Inverse Bean with 3% noise Ellipse Peanut Bean λ = sin t
Nonlinear 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 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 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 informationReconstructing 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 informationIterative regularization of nonlinear ill-posed problems in Banach space
Iterative regularization of nonlinear ill-posed problems in Banach space Barbara Kaltenbacher, University of Klagenfurt joint work with Bernd Hofmann, Technical University of Chemnitz, Frank Schöpfer and
More informationA 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 informationRegularization and Inverse Problems
Regularization and Inverse Problems Caroline Sieger Host Institution: Universität Bremen Home Institution: Clemson University August 5, 2009 Caroline Sieger (Bremen and Clemson) Regularization and Inverse
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 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 Chi-squared Distribution of the Regularized Least Squares Functional for Regularization Parameter Estimation
The Chi-squared Distribution of the Regularized Least Squares Functional for Regularization Parameter Estimation Rosemary Renaut DEPARTMENT OF MATHEMATICS AND STATISTICS Prague 2008 MATHEMATICS AND STATISTICS
More informationRECENT DEVELOPMENTS IN INVERSE ACOUSTIC SCATTERING THEORY
RECENT DEVELOPMENTS IN INVERSE ACOUSTIC SCATTERING THEORY DAVID COLTON, JOE COYLE, AND PETER MONK Abstract. We survey some of the highlights of inverse scattering theory as it has developed over the past
More informationDIRECT SAMPLING METHODS FOR INVERSE SCATTERING PROBLEMS
Michigan Technological University Digital Commons @ Michigan Tech Dissertations, Master's Theses and Master's Reports 2017 DIRECT SAMPLING METHODS FOR INVERSE SCATTERING PROBLEMS Ala Mahmood Nahar Al Zaalig
More informationRegularization in Banach Space
Regularization in Banach Space Barbara Kaltenbacher, Alpen-Adria-Universität Klagenfurt joint work with Uno Hämarik, University of Tartu Bernd Hofmann, Technical University of Chemnitz Urve Kangro, University
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 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 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 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 informationA Hybrid Method for Inverse Obstacle Scattering Problems
A Hybrid Method for Inverse Obstacle Scattering Problems Dissertation zur Erlangung des Doktorgrades der Mathematisch-Naturwissenschaftlichen Fakultäten der Georg-August-Universität zu Göttingen vorgelegt
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 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 informationRecent Developments in Inverse Acoustic Scattering Theory
SIAM REVIEW Vol. 42, No. 3, pp. 369 414 c 2000 Society for Industrial and Applied Mathematics Recent Developments in Inverse Acoustic Scattering Theory David Colton Joe Coyle Peter Monk Abstract. We survey
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 informationWhen is the error in the h BEM for solving the Helmholtz equation bounded independently of k?
BIT manuscript No. (will be inserted by the editor) When is the error in the h BEM for solving the Helmholtz equation bounded independently of k? I. G. Graham M. Löhndorf J. M. Melenk E. A. Spence Received:
More informationNonlinear error dynamics for cycled data assimilation methods
Nonlinear error dynamics for cycled data assimilation methods A J F Moodey 1, A S Lawless 1,2, P J van Leeuwen 2, R W E Potthast 1,3 1 Department of Mathematics and Statistics, University of Reading, UK.
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 informationInverse problems based on partial differential equations and integral equations
Inverse problems based on partial differential equations and integral equations Jijun Liu Department of Mathematics, Southeast University E-mail: jjliu@seu.edu.cn Korea, November 22, 214 Jijun Liu Inverse
More informationKirchhoff, 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 informationIterative Methods for Ill-Posed Problems
Iterative Methods for Ill-Posed Problems Based on joint work with: Serena Morigi Fiorella Sgallari Andriy Shyshkov Salt Lake City, May, 2007 Outline: Inverse and ill-posed problems Tikhonov regularization
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 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 informationConvergence rates of spectral methods for statistical inverse learning problems
Convergence rates of spectral methods for statistical inverse learning problems G. Blanchard Universtität Potsdam UCL/Gatsby unit, 04/11/2015 Joint work with N. Mücke (U. Potsdam); N. Krämer (U. München)
More informationA numerical algorithm for solving 3D inverse scattering problem with non-over-determined data
A numerical algorithm for solving 3D inverse scattering problem with non-over-determined data Alexander Ramm, Cong Van Department of Mathematics, Kansas State University, Manhattan, KS 66506, USA ramm@math.ksu.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 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 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 informationLinear Inverse Problems
Linear Inverse Problems Ajinkya Kadu Utrecht University, The Netherlands February 26, 2018 Outline Introduction Least-squares Reconstruction Methods Examples Summary Introduction 2 What are inverse problems?
More informationA Hybrid LSQR Regularization Parameter Estimation Algorithm for Large Scale Problems
A Hybrid LSQR Regularization Parameter Estimation Algorithm for Large Scale Problems Rosemary Renaut Joint work with Jodi Mead and Iveta Hnetynkova SIAM Annual Meeting July 10, 2009 National Science Foundation:
More informationHigh Frequency Scattering by Convex Polygons Stephen Langdon
Bath, October 28th 2005 1 High Frequency Scattering by Convex Polygons Stephen Langdon University of Reading, UK Joint work with: Simon Chandler-Wilde Steve Arden Funded by: Leverhulme Trust University
More informationAnalysis of the Hessian for inverse scattering problems: I. Inverse shape scattering of
Home Search Collections Journals About Contact us My IOPscience Analysis of the Hessian for inverse scattering problems: I. Inverse shape scattering of acoustic waves This article has been downloaded from
More informationREGULARIZATION PARAMETER SELECTION IN DISCRETE ILL POSED PROBLEMS THE USE OF THE U CURVE
Int. J. Appl. Math. Comput. Sci., 007, Vol. 17, No., 157 164 DOI: 10.478/v10006-007-0014-3 REGULARIZATION PARAMETER SELECTION IN DISCRETE ILL POSED PROBLEMS THE USE OF THE U CURVE DOROTA KRAWCZYK-STAŃDO,
More informationICES REPORT Analysis of the Hessian for Inverse Scattering Problems. Part I: Inverse Shape Scattering of Acoustic Waves
ICES REPORT 11-20 June 2011 Analysis of the Hessian for Inverse Scattering Problems. Part I: Inverse Shape Scattering of Acoustic Waves by Tan Bui-Thanh and Omar Ghattas The Institute for Computational
More informationA posteriori error estimates for the adaptivity technique for the Tikhonov functional and global convergence for a coefficient inverse problem
A posteriori error estimates for the adaptivity technique for the Tikhonov functional and global convergence for a coefficient inverse problem Larisa Beilina Michael V. Klibanov December 18, 29 Abstract
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 informationA far-field based T-matrix method for three dimensional acoustic scattering
ANZIAM J. 50 (CTAC2008) pp.c121 C136, 2008 C121 A far-field based T-matrix method for three dimensional acoustic scattering M. Ganesh 1 S. C. Hawkins 2 (Received 14 August 2008; revised 4 October 2008)
More informationTwo-parameter regularization method for determining the heat source
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 13, Number 8 (017), pp. 3937-3950 Research India Publications http://www.ripublication.com Two-parameter regularization method for
More informationNumerische Mathematik
Numer. Math. 1999 83: 139 159 Numerische Mathematik c Springer-Verlag 1999 On an a posteriori parameter choice strategy for Tikhonov regularization of nonlinear ill-posed problems Jin Qi-nian 1, Hou Zong-yi
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 information4.2 Green s representation theorem
4.2. REEN S REPRESENTATION THEOREM 57 i.e., the normal velocity on the boundary is proportional to the ecess pressure on the boundary. The coefficient χ is called the acoustic impedance of the obstacle
More informationThe Chi-squared Distribution of the Regularized Least Squares Functional for Regularization Parameter Estimation
The Chi-squared Distribution of the Regularized Least Squares Functional for Regularization Parameter Estimation Rosemary Renaut Collaborators: Jodi Mead and Iveta Hnetynkova DEPARTMENT OF MATHEMATICS
More informationA Galerkin boundary element method for high frequency scattering by convex polygons
A Galerkin boundary element method for high frequency scattering by convex polygons Article Published Version Chandler Wilde, S. N. and Langdon, S. (2007) A Galerkin boundary element method for high frequency
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 to the boundary element method. A case study: the Helmholtz equation
Introduction to the boundary element method. A case study: the Helmholtz equation Francisco Javier Sayas January 2006 Contents 1 The Helmholtz equation 2 2 Scattering 4 3 Single layer acoustic potentials
More 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 informationAcoustic scattering : high frequency boundary element methods and unified transform methods
Acoustic scattering : high frequency boundary element methods and unified transform methods Book or Report Section Accepted Version Chandler Wilde, S.N. and Langdon, S. (2015) Acoustic scattering : high
More informationA semismooth Newton method for L 1 data fitting with automatic choice of regularization parameters and noise calibration
A semismooth Newton method for L data fitting with automatic choice of regularization parameters and noise calibration Christian Clason Bangti Jin Karl Kunisch April 26, 200 This paper considers the numerical
More informationAdaptive discretization and first-order methods for nonsmooth inverse problems for PDEs
Adaptive discretization and first-order methods for nonsmooth inverse problems for PDEs Christian Clason Faculty of Mathematics, Universität Duisburg-Essen joint work with Barbara Kaltenbacher, Tuomo Valkonen,
More informationThe Calderon-Vaillancourt Theorem
The Calderon-Vaillancourt Theorem What follows is a completely self contained proof of the Calderon-Vaillancourt Theorem on the L 2 boundedness of pseudo-differential operators. 1 The result Definition
More informationConvergence Rates in Regularization for Nonlinear Ill-Posed Equations Involving m-accretive Mappings in Banach Spaces
Applied Mathematical Sciences, Vol. 6, 212, no. 63, 319-3117 Convergence Rates in Regularization for Nonlinear Ill-Posed Equations Involving m-accretive Mappings in Banach Spaces Nguyen Buong Vietnamese
More informationFredholm Theory. April 25, 2018
Fredholm Theory April 25, 208 Roughly speaking, Fredholm theory consists of the study of operators of the form I + A where A is compact. From this point on, we will also refer to I + A as Fredholm operators.
More informationConjugate Gradient algorithm. Storage: fixed, independent of number of steps.
Conjugate Gradient algorithm Need: A symmetric positive definite; Cost: 1 matrix-vector product per step; Storage: fixed, independent of number of steps. The CG method minimizes the A norm of the error,
More informationLecture 2: Tikhonov-Regularization
Lecture 2: Tikhonov-Regularization Bastian von Harrach harrach@math.uni-stuttgart.de Chair of Optimization and Inverse Problems, University of Stuttgart, Germany Advanced Instructional School on Theoretical
More informationPhD Course: Introduction to Inverse Problem. Salvatore Frandina Siena, August 19, 2012
PhD Course: to Inverse Problem salvatore.frandina@gmail.com theory Department of Information Engineering, Siena, Italy Siena, August 19, 2012 1 / 68 An overview of the - - - theory 2 / 68 Direct and Inverse
More informationMOSCO 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 informationInverse problems and medical imaging
Inverse problems and medical imaging Bastian von Harrach harrach@math.uni-frankfurt.de Institute of Mathematics, Goethe University Frankfurt, Germany Seminario di Calcolo delle Variazioni ed Equazioni
More informationA GALERKIN BOUNDARY ELEMENT METHOD FOR HIGH FREQUENCY SCATTERING BY CONVEX POLYGONS
A GALERKIN BOUNDARY ELEMENT METHOD FOR HIGH FREQUENCY SCATTERING BY CONVEX POLYGONS S N CHANDLER-WILDE AND S LANGDON Abstract. In this paper we consider the problem of time-harmonic acoustic scattering
More informationUnbiased Risk Estimation as Parameter Choice Rule for Filter-based Regularization Methods
Unbiased Risk Estimation as Parameter Choice Rule for Filter-based Regularization Methods Frank Werner 1 Statistical Inverse Problems in Biophysics Group Max Planck Institute for Biophysical Chemistry,
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 informationAn Iteratively Regularized Projection Method for Nonlinear Ill-posed Problems
Int. J. Contemp. Math. Sciences, Vol. 5, 2010, no. 52, 2547-2565 An Iteratively Regularized Projection Method for Nonlinear Ill-posed Problems Santhosh George Department of Mathematical and Computational
More informationRobust error estimates for regularization and discretization of bang-bang control problems
Robust error estimates for regularization and discretization of bang-bang control problems Daniel Wachsmuth September 2, 205 Abstract We investigate the simultaneous regularization and discretization of
More informationInitial Temperature Reconstruction for a Nonlinear Heat Equation: Application to Radiative Heat Transfer.
Proceedings of the 5th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15th July 2005 Initial Temperature Reconstruction for a Nonlinear Heat Equation:
More informationLevenberg-Marquardt method in Banach spaces with general convex regularization terms
Levenberg-Marquardt method in Banach spaces with general convex regularization terms Qinian Jin Hongqi Yang Abstract We propose a Levenberg-Marquardt method with general uniformly convex regularization
More informationA globally convergent numerical method and adaptivity for an inverse problem via Carleman estimates
A globally convergent numerical method and adaptivity for an inverse problem via Carleman estimates Larisa Beilina, Michael V. Klibanov Chalmers University of Technology and Gothenburg University, Gothenburg,
More informationFinite Element Analysis of Acoustic Scattering
Frank Ihlenburg Finite Element Analysis of Acoustic Scattering With 88 Illustrations Springer Contents Preface vii 1 The Governing Equations of Time-Harmonic Wave Propagation, 1 1.1 Acoustic Waves 1 1.1.1
More informationError analysis and fast solvers for high-frequency scattering problems
Error analysis and fast solvers for high-frequency scattering problems I.G. Graham (University of Bath) Woudschoten October 2014 High freq. problem for the Helmholtz equation Given an object Ω R d, with
More informationThe linear sampling method and energy conservation
The linear sampling method and energy conservation R Aramini, G Caviglia, A Massa and M Piana $ Dipartimento di Ingegneria e Scienza dell Informazione, Università di Trento, via Sommarive 4, 3823 Povo
More informationA model function method in total least squares
www.oeaw.ac.at A model function method in total least squares S. Lu, S. Pereverzyev, U. Tautenhahn RICAM-Report 2008-18 www.ricam.oeaw.ac.at A MODEL FUNCTION METHOD IN TOTAL LEAST SQUARES SHUAI LU, SERGEI
More informationStatistically-Based Regularization Parameter Estimation for Large Scale Problems
Statistically-Based Regularization Parameter Estimation for Large Scale Problems Rosemary Renaut Joint work with Jodi Mead and Iveta Hnetynkova March 1, 2010 National Science Foundation: Division of Computational
More informationON THE DYNAMICAL SYSTEMS METHOD FOR SOLVING NONLINEAR EQUATIONS WITH MONOTONE OPERATORS
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume, Number, Pages S -9939(XX- ON THE DYNAMICAL SYSTEMS METHOD FOR SOLVING NONLINEAR EQUATIONS WITH MONOTONE OPERATORS N. S. HOANG AND A. G. RAMM (Communicated
More informationIll-Posedness of Backward Heat Conduction Problem 1
Ill-Posedness of Backward Heat Conduction Problem 1 M.THAMBAN NAIR Department of Mathematics, IIT Madras Chennai-600 036, INDIA, E-Mail mtnair@iitm.ac.in 1. Ill-Posedness of Inverse Problems Problems that
More informationIntegral Representation Formula, Boundary Integral Operators and Calderón projection
Integral Representation Formula, Boundary Integral Operators and Calderón projection Seminar BEM on Wave Scattering Franziska Weber ETH Zürich October 22, 2010 Outline Integral Representation Formula Newton
More informationChapter 7: Bounded Operators in Hilbert Spaces
Chapter 7: Bounded Operators in Hilbert Spaces I-Liang Chern Department of Applied Mathematics National Chiao Tung University and Department of Mathematics National Taiwan University Fall, 2013 1 / 84
More informationSolutions: Problem Set 4 Math 201B, Winter 2007
Solutions: Problem Set 4 Math 2B, Winter 27 Problem. (a Define f : by { x /2 if < x
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 informationFunctionalanalytic tools and nonlinear equations
Functionalanalytic tools and nonlinear equations Johann Baumeister Goethe University, Frankfurt, Germany Rio de Janeiro / October 2017 Outline Fréchet differentiability of the (PtS) mapping. Nonlinear
More informationA DISCREPANCY PRINCIPLE FOR PARAMETER SELECTION IN LOCAL REGULARIZATION OF LINEAR VOLTERRA INVERSE PROBLEMS. Cara Dylyn Brooks A DISSERTATION
A DISCREPANCY PRINCIPLE FOR PARAMETER SELECTION IN LOCAL REGULARIZATION OF LINEAR VOLTERRA INVERSE PROBLEMS By Cara Dylyn Brooks A DISSERTATION Submitted to Michigan State University in partial fulfillment
More informationMeasure and Integration: Solutions of CW2
Measure and Integration: s of CW2 Fall 206 [G. Holzegel] December 9, 206 Problem of Sheet 5 a) Left (f n ) and (g n ) be sequences of integrable functions with f n (x) f (x) and g n (x) g (x) for almost
More informationarxiv: v1 [math.na] 11 Mar 2019 A high frequency boundary element method for scattering by a class of multiple obstacles
(2019) Page 1 of 39 arxiv:1903.04449v1 [math.na] 11 Mar 2019 A high frequency boundary element method for scattering by a class of multiple obstacles ANDREW GIBBS, DEPT. OF COMPUTER SCIENCE, K.U. LEUVEN,
More informationPerron method for the Dirichlet problem.
Introduzione alle equazioni alle derivate parziali, Laurea Magistrale in Matematica Perron method for the Dirichlet problem. We approach the question of existence of solution u C (Ω) C(Ω) of the Dirichlet
More informationGEOPHYSICAL INVERSE THEORY AND REGULARIZATION PROBLEMS
Methods in Geochemistry and Geophysics, 36 GEOPHYSICAL INVERSE THEORY AND REGULARIZATION PROBLEMS Michael S. ZHDANOV University of Utah Salt Lake City UTAH, U.S.A. 2OO2 ELSEVIER Amsterdam - Boston - London
More informationIterative Regularization Methods for Inverse Problems: Lecture 3
Iterative Regularization Methods for Inverse Problems: Lecture 3 Thorsten Hohage Institut für Numerische und Angewandte Mathematik Georg-August Universität Göttingen Madrid, April 11, 2011 Outline 1 regularization
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 informationGeneralized Local Regularization for Ill-Posed Problems
Generalized Local Regularization for Ill-Posed Problems Patricia K. Lamm Department of Mathematics Michigan State University AIP29 July 22, 29 Based on joint work with Cara Brooks, Zhewei Dai, and Xiaoyue
More informationInverse problems and medical imaging
Inverse problems and medical imaging Bastian von Harrach harrach@math.uni-frankfurt.de Institute of Mathematics, Goethe University Frankfurt, Germany Colloquium of the Department of Mathematics Saarland
More informationInverse Scattering Theory
Chapter 1 Inverse Scattering Theory In this introductory chapter we provide an overview of the basic ideas of inverse scattering theory for the special case when the scattering object is an isotropic inhomogeneous
More informationSparse Recovery in Inverse Problems
Radon Series Comp. Appl. Math XX, 1 63 de Gruyter 20YY Sparse Recovery in Inverse Problems Ronny Ramlau and Gerd Teschke Abstract. Within this chapter we present recent results on sparse recovery algorithms
More informationIdentification of Temperature Dependent Parameters in a Simplified Radiative Heat Transfer
Australian Journal of Basic and Applied Sciences, 5(): 7-4, 0 ISSN 99-878 Identification of Temperature Dependent Parameters in a Simplified Radiative Heat Transfer Oliver Tse, Renè Pinnau, Norbert Siedow
More informationLocating Multiple Multiscale Acoustic Scatterers
Locating Multiple Multiscale Acoustic Scatterers Jingzhi Li Hongyu Liu Jun Zou March, Abstract We develop three inverse scattering schemes for locating multiple multiscale acoustic scatterers in a very
More informationModel-aware Newton-type regularization in electrical impedance tomography
Model-aware Newton-type regularization in electrical impedance tomography Andreas Rieder Robert Winkler FAKULTÄT FÜR MATHEMATIK INSTITUT FÜR ANGEWANDTE UND NUMERISCHE MATHEMATIK SFB 1173 KIT University
More informationSolutions for Problem Set #5 due October 17, 2003 Dustin Cartwright and Dylan Thurston
Solutions for Problem Set #5 due October 17, 23 Dustin Cartwright and Dylan Thurston 1 (B&N 6.5) Suppose an analytic function f agrees with tan x, x 1. Show that f(z) = i has no solution. Could f be entire?
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 information