Reconstructing inclusions from Electrostatic Data
|
|
- Lily Bell
- 6 years ago
- Views:
Transcription
1 Reconstructing inclusions from Electrostatic Data Isaac Harris Texas A&M University, Department of Mathematics College Station, Texas Joint work with: W. Rundell Purdue University - November 2017 Reconstructing inclusions from Electrostatic Data 1 / 35
2 Outline Introduction Formulation of the Problem Solution of the Inverse Problem Solution to the Inverse Shape Problem Solution to the Inverse Impedance Problem Current/Future Research Reconstructing inclusions from Electrostatic Data 2 / 35
3 What is an Inverse Problem? Input = Physical Model = Output Input: Boundary Data Physical Model: Partial Differential Equation Output: Values of the Solution Forward Problem: Input + Physical Model = Output Inverse Problem: Output + Input = Physical Model Reconstructing inclusions from Electrostatic Data 3 / 35
4 An Inverse Scattering Problem The Physical Model: Acoustic or Electromagnetic wave scattering. Inverse Problem Reconstruct the support of the defective region in a known material from the measured scattered wave. Reconstructing inclusions from Electrostatic Data 4 / 35
5 Some Important Questions Assume your inverse problem is given by A(f ) = g Well-Posedness: Solvability, Injectivity and Stability Does A 1 exists as a continuous mapping? (Usually Not!!!) Regularization Techniques Approximate (point-wise) A 1 by a continuous mapping Inversion Algoritms Develop a stable and computable method to solve for f Reconstructing inclusions from Electrostatic Data 5 / 35
6 Formulation of the Problem I. Harris and W. Rundell A direct method for reconstructing inclusions and boundary conditions from electrostatic data Submitted (arxiv: ) Reconstructing inclusions from Electrostatic Data 6 / 35
7 Electrical Impedance Tomography (EIT) EIT is a non-destructive imaging method in which the physical parameters are recovered from surface electrode measurements. Figure: Picture taken from Reconstructing inclusions from Electrostatic Data 7 / 35
8 The Inverse Electrostatic Problem Inverse Problem Recover the inclusion Ω and the boundary condition B from the knowledge of the Dirichlet to Neumann operators Λf = ν u D and Λ 0 f = ν u 0 D. Reconstructing inclusions from Electrostatic Data 8 / 35
9 In this talk we focus on B(u 0 ) = ν u 0 + γ(x)u 0 Assumptions: 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 Ω The results in this talk still hold for the case when u is replaced by A(x) u + c(x)u with symmetric A(x) C 1 (D, C d d ), where the Re(A) is positive definite and Im(A) is non-positive. As well as c(x) L (D) such that Re(c) is non-negative. Reconstructing inclusions from Electrostatic Data 9 / 35
10 Some Previous Works Uniqueness for the inverse problem for 2 pairs of input/output. 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 10 / 35
11 Iterative Reconstruction Methods Iterative/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 on D. Λ 0 f j : ( Ω, γ) ν u (j) 0 for j = 1, 2 Then you would solve the system of non-linear (exponentially) ill-posed equation using a iterative scheme. Reconstructing inclusions from Electrostatic Data 11 / 35
12 Restrictions for Iterative Methods Good initial estimates are usually needed to insure convergence. This requires a priori 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 priori information? A lack of information cannot be remedied by any mathematical trickery -C. Lanczos Reconstructing inclusions from Electrostatic Data 12 / 35
13 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 13 / 35
14 Solution to the Inverse Shape Problem Reconstructing inclusions from Electrostatic Data 14 / 35
15 The Green s Function & Range Test 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 Th = ν w D where w = 0 in D \ Ω with w D = 0 and ν w Ω = h Theorem (I.Harris and W.Rundell) ν G(, z) D Range(T ) if and only if z Ω. Where the operator T : H 1/2 ( Ω) H 1/2 ( D). Main Idea: Using that G(x, z) as x z and the Range of T consists of smooth functions gives the result. Reconstructing inclusions from Electrostatic Data 15 / 35
16 Theorem (I.Harris and W.Rundell) The difference of the DtN operators (Λ Λ 0 ) : H 1/2 ( D) H 1/2 ( D) is compact, injective with a dense range (under some assumptions). Moreover, Range(Λ Λ 0 ) Range(T ). Techniques used in the proof 1. Compactness follows from elliptic regularity 2. Injectivity is a consequence on unique continuation 3. Green s 1st identity gives that (Λ Λ 0 ) is symmetric 4. Derive the factorization (Λ Λ 0 ) = TH Reconstructing inclusions from Electrostatic Data 16 / 35
17 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(T ) which is a contradiction provided that z / Ω. Reconstructing inclusions from Electrostatic Data 17 / 35
18 Determining Ω via Sampling Theorem (I.Harris and W.Rundell) 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 18 / 35
19 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 19 / 35
20 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 20 / 35
21 Solution to the Inverse Impedance Problem Reconstructing inclusions from Electrostatic Data 21 / 35
22 Theorem (I.Harris and W.Rundell) The DtN mapping Λ 0 : H 1/2 ( D) H 1/2 ( D) uniquely determines the impedance parameter γ(x) L + ( Ω). Techniques used in the proof 1. Assume γ 1 and γ 2 give that same DtN mapping 2. Unique continuation implies that (γ 1 γ 2 )u Ω 0 = 0 for all f H 1/2 ( D) 3. Then prove that the set { } u Ω 0 : for all f H 1/2 ( D) is a dense in L 2 ( Ω) Reconstructing inclusions from Electrostatic Data 22 / 35
23 Data Completion for Recovering γ(x) We wish to derive a non-iterative method for recovering the impedance parameter using Data Completion, ( f, Λ0 f ) D ( u f 0, ν u f 0) Ω. Using the Boundary Condition on Ω we have that (Direct Solution) γ(x n ) = νu f 0 (x n) u f 0 (x n) for n = 1,, N with x n = x(θ n ) Ω. Reconstructing inclusions from Electrostatic Data 23 / 35
24 Data Completion Algorithm I Let D = Γ m (measurements) and Ω = Γ i (impedance) 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 24 / 35
25 Data Completion Algorithm II We make the ansatz that u 0 (x) = (D m ϕ)(x) + ( D i ψ)(x) for x D \ Ω. Now assume that we have ( f, Λ 0 f ) 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 25 / 35
26 Data Completion Algorithm III Define the operator associated with the system of BIE [ I K A = mm K ] im K mi I + K ii Theorem (I.Harris and W.Rundell) The mapping A is a bounded linear operator from H 1/2 (Γ m ) H 1/2 (Γ i ) to itself and has a bounded inverse. Techniques used in the proof 1. Prove that A is Fredholm of index zero 2. Prove injectivity using the jump-relations Reconstructing inclusions from Electrostatic Data 26 / 35
27 Data Completion Algorithm IV Using the fact the we know the normal derivative on we have that [ I K Λ 0 f = [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. Notice: This reduces to solving a linear system of equation that is independent of the boundary condition being recovered. Reconstructing inclusions from Electrostatic Data 27 / 35
28 Reconstruction of γ ( x(θ) ) = 2 sin 4 (θ) Figure: Left: Reconstructed boundary of the inclusion. Right: Reconstructed impedance parameterγ ( x(θ) ) = 2 sin 4 (θ). Reconstructing inclusions from Electrostatic Data 28 / 35
29 Current/Future Research Reconstructing inclusions from Electrostatic Data 29 / 35
30 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. I. Harris Detecting inclusions with a generalized impedance condition from electrostatic data via sampling Submitted (arxiv: ) Reconstructing inclusions from Electrostatic Data 30 / 35
31 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(γ) γ min > 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 / 35
32 Reconstruction of an Inclusion with GIBC Let the coefficients be given by η ( x(θ) ) = 1 and γ ( x(θ) ) = sin 2 (θ). 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 / 35
33 Inverse Source Problem for the Diffusion ( t )u(x, t) = F (x) for (x, t) R d (0, ). with u(x, 0) = 0. We consider two types of sources F (x) = M c m δ(x z m ) and F (x) = m=1 M m=1 χ Dm with constants c m R where the points z m R d and D m R d are bounded open sets. Inverse Problem Recover the heat source F (x) such that Supp(F ) B R from the boundary measurements { u(x j, t) } j J where x j B R. Reconstructing inclusions from Electrostatic Data 33 / 35
34 Some Preliminary Results Three measurements { u(x j, t) } 3 uniquely determines the j=1 Point Sources F (x) = M c m δ(x z m ) for any M N m=1 Recovering the Center Points z m and D m is asymptotic to the point source problem for small regions Exact z m Recovered z m Exact D m Recovered D m (1, 1) (1.04, 0.98) (2, 2) (2.10, 2.10) Table: The reconstructed center points z m and area D m Reconstructing inclusions from Electrostatic Data 34 / 35
35 Figure: Questions? Reconstructing inclusions from Electrostatic Data 35 / 35
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationApproximation of inverse boundary value problems by phase-field methods
Approximation of inverse boundary value problems by phase-field methods Luca RONDI Università degli Studi di Trieste Dipartimento di Matematica e Informatica MSRI Berkeley, 3 December 2010 Luca RONDI (Università
More informationThe inverse conductivity problem with power densities in dimension n 2
The inverse conductivity problem with power densities in dimension n 2 François Monard Guillaume Bal Dept. of Applied Physics and Applied Mathematics, Columbia University. June 19th, 2012 UC Irvine Conference
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 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 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 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 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 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 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 informationA GLOBALLY CONVERGENT NUMERICAL METHOD FOR SOME COEFFICIENT INVERSE PROBLEMS WITH RESULTING SECOND ORDER ELLIPTIC EQUATIONS
A GLOBALLY CONVERGENT NUMERICAL METHOD FOR SOME COEFFICIENT INVERSE PROBLEMS WITH RESULTING SECOND ORDER ELLIPTIC EQUATIONS LARISA BEILINA AND MICHAEL V. KLIBANOV Abstract. A new globally convergent numerical
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 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 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 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 informationRecovering point sources in unknown environment with differential data
Recovering point sources in unknown environment with differential data Yimin Zhong University of Texas at Austin April 29, 2013 Introduction Introduction The field u satisfies u(x) + k 2 (1 + n(x))u(x)
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 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 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 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 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 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 informationA uniqueness result and image reconstruction of the orthotropic conductivity in magnetic resonance electrical impedance tomography
c de Gruyter 28 J. Inv. Ill-Posed Problems 16 28), 381 396 DOI 1.1515 / JIIP.28.21 A uniqueness result and image reconstruction of the orthotropic conductivity in magnetic resonance electrical impedance
More informationSolving the complete-electrode direct model of ERT using the boundary element method and the method of fundamental solutions
Solving the complete-electrode direct model of ERT using the boundary element method and the method of fundamental solutions T. E. Dyhoum 1,2, D. Lesnic 1, and R. G. Aykroyd 2 1 Department of Applied Mathematics,
More informationRegularizing the optimization-based solution of inverse coefficient problems with monotonicity constraints
Regularizing the optimization-based solution of inverse coefficient problems with monotonicity constraints Bastian von Harrach (joint work with Mach Nguyet Minh) http://numerical.solutions Institute of
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 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 informationThe X-ray transform for a non-abelian connection in two dimensions
The X-ray transform for a non-abelian connection in two dimensions David Finch Department of Mathematics Oregon State University Corvallis, OR, 97331, USA Gunther Uhlmann Department of Mathematics University
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 informationare harmonic functions so by superposition
J. Rauch Applied Complex Analysis The Dirichlet Problem Abstract. We solve, by simple formula, the Dirichlet Problem in a half space with step function boundary data. Uniqueness is proved by complex variable
More informationApplied Numerical Mathematics. High-order numerical schemes based on difference potentials for 2D elliptic problems with material interfaces
Applied Numerical Mathematics 111 (2017) 64 91 Contents lists available at ScienceDirect Applied Numerical Mathematics www.elsevier.com/locate/apnum High-order numerical schemes based on difference potentials
More informationApplications of the Maximum Principle
Jim Lambers MAT 606 Spring Semester 2015-16 Lecture 26 Notes These notes correspond to Sections 7.4-7.6 in the text. Applications of the Maximum Principle The maximum principle for Laplace s equation is
More informationInverse parameter identification problems
Inverse parameter identification problems Bastian von Harrach harrach@math.uni-stuttgart.de Chair of Optimization and Inverse Problems, University of Stuttgart, Germany ICP - Institute for Computational
More informationCalderón s inverse problem in 2D Electrical Impedance Tomography
Calderón s inverse problem in 2D Electrical Impedance Tomography Kari Astala (University of Helsinki) Joint work with: Matti Lassas, Lassi Päivärinta, Samuli Siltanen, Jennifer Mueller and Alan Perämäki
More informationFast evaluation of boundary integral operators arising from an eddy current problem. MPI MIS Leipzig
1 Fast evaluation of boundary integral operators arising from an eddy current problem Steffen Börm MPI MIS Leipzig Jörg Ostrowski ETH Zürich Induction heating 2 Induction heating: Time harmonic eddy current
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 informationM.Sc. in Meteorology. Numerical Weather Prediction
M.Sc. in Meteorology UCD Numerical Weather Prediction Prof Peter Lynch Meteorology & Climate Centre School of Mathematical Sciences University College Dublin Second Semester, 2005 2006. In this section
More informationApplications of the periodic unfolding method to multi-scale problems
Applications of the periodic unfolding method to multi-scale problems Doina Cioranescu Université Paris VI Santiago de Compostela December 14, 2009 Periodic unfolding method and multi-scale problems 1/56
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 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 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 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 informationRecent progress on the explicit inversion of geodesic X-ray transforms
Recent progress on the explicit inversion of geodesic X-ray transforms François Monard Department of Mathematics, University of Washington. Geometric Analysis and PDE seminar University of Cambridge, May
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 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 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 informationInverse problems for hyperbolic PDEs
Inverse problems for hyperbolic PDEs Lauri Oksanen University College London Example: inverse problem for the wave equation Let c be a smooth function on Ω R n and consider the wave equation t 2 u c 2
More informationPoincaré Duality Angles on Riemannian Manifolds with Boundary
Poincaré Duality Angles on Riemannian Manifolds with Boundary Clayton Shonkwiler Department of Mathematics University of Pennsylvania June 5, 2009 Realizing cohomology groups as spaces of differential
More informationGreen s Functions and Distributions
CHAPTER 9 Green s Functions and Distributions 9.1. Boundary Value Problems We would like to study, and solve if possible, boundary value problems such as the following: (1.1) u = f in U u = g on U, where
More informationIncreasing stability in an inverse problem for the acoustic equation
Increasing stability in an inverse problem for the acoustic equation Sei Nagayasu Gunther Uhlmann Jenn-Nan Wang Abstract In this work we study the inverse boundary value problem of determining the refractive
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 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 informationRecovery of anisotropic metrics from travel times
Purdue University The Lens Rigidity and the Boundary Rigidity Problems Let M be a bounded domain with boundary. Let g be a Riemannian metric on M. Define the scattering relation σ and the length (travel
More informationAn inverse elliptic problem of medical optics with experimental data
An inverse elliptic problem of medical optics with experimental data Jianzhong Su, Michael V. Klibanov, Yueming Liu, Zhijin Lin Natee Pantong and Hanli Liu Abstract A numerical method for an inverse problem
More informationNovel tomography techniques and parameter identification problems
Novel tomography techniques and parameter identification problems Bastian von Harrach harrach@ma.tum.de Department of Mathematics - M1, Technische Universität München Colloquium of the Institute of Biomathematics
More informationA Direct Sampling Method for Electrical Impedance Tomography
A Direct Sampling Method for Electrical Impedance Tomography Yat Tin Chow Kazufumi Ito Jun Zou Abstract This work investigates the electrical impedance tomography (EIT) problem in the case when only one
More informationFEM-FEM and FEM-BEM Coupling within the Dune Computational Software Environment
FEM-FEM and FEM-BEM Coupling within the Dune Computational Software Environment Alastair J. Radcliffe Andreas Dedner Timo Betcke Warwick University, Coventry University College of London (UCL) U.K. Radcliffe
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
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 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 informationResistor Networks and Optimal Grids for Electrical Impedance Tomography with Partial Boundary Measurements
Resistor Networks and Optimal Grids for Electrical Impedance Tomography with Partial Boundary Measurements Alexander Mamonov 1, Liliana Borcea 2, Vladimir Druskin 3, Fernando Guevara Vasquez 4 1 Institute
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 informationSOLVING ELLIPTIC PDES
university-logo SOLVING ELLIPTIC PDES School of Mathematics Semester 1 2008 OUTLINE 1 REVIEW 2 POISSON S EQUATION Equation and Boundary Conditions Solving the Model Problem 3 THE LINEAR ALGEBRA PROBLEM
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 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 informationUniversity of North Carolina at Charlotte Charlotte, USA Norwegian Institute of Science and Technology Trondheim, Norway
1 A globally convergent numerical method for some coefficient inverse problems Michael V. Klibanov, Larisa Beilina University of North Carolina at Charlotte Charlotte, USA Norwegian Institute of Science
More informationMath 46, Applied Math (Spring 2009): Final
Math 46, Applied Math (Spring 2009): Final 3 hours, 80 points total, 9 questions worth varying numbers of points 1. [8 points] Find an approximate solution to the following initial-value problem which
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 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 informationSolving Ill-Posed Cauchy Problems in Three Space Dimensions using Krylov Methods
Solving Ill-Posed Cauchy Problems in Three Space Dimensions using Krylov Methods Lars Eldén Department of Mathematics Linköping University, Sweden Joint work with Valeria Simoncini February 21 Lars Eldén
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 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 informationHomogenization of micro-resonances and localization of waves.
Homogenization of micro-resonances and localization of waves. Valery Smyshlyaev University College London, UK July 13, 2012 (joint work with Ilia Kamotski UCL, and Shane Cooper Bath/ Cardiff) Valery Smyshlyaev
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 informationOn some weighted fractional porous media equations
On some weighted fractional porous media equations Gabriele Grillo Politecnico di Milano September 16 th, 2015 Anacapri Joint works with M. Muratori and F. Punzo Gabriele Grillo Weighted Fractional PME
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 informationSimultaneous reconstruction of outer boundary shape and conductivity distribution in electrical impedance tomography
Simultaneous reconstruction of outer boundary shape and conductivity distribution in electrical impedance tomography Nuutti Hyvönen Aalto University nuutti.hyvonen@aalto.fi joint work with J. Dardé, A.
More informationElectric potentials with localized divergence properties
Electric potentials with localized divergence properties Bastian Gebauer bastian.gebauer@oeaw.ac.at Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy of Sciences,
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 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 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 informationDeforming conformal metrics with negative Bakry-Émery Ricci Tensor on manifolds with boundary
Deforming conformal metrics with negative Bakry-Émery Ricci Tensor on manifolds with boundary Weimin Sheng (Joint with Li-Xia Yuan) Zhejiang University IMS, NUS, 8-12 Dec 2014 1 / 50 Outline 1 Prescribing
More informationCauchy problem for Ultrasound Modulated EIT
Cauchy problem for Ultrasound Modulated EIT Guillaume Bal December 18, 2011 Abstract Ultrasound modulation of electrical or optical properties of materials offers the possibility to devise hybrid imaging
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 informationc 2006 Society for Industrial and Applied Mathematics
SIAM J. APPL. MATH. Vol. 66, No. 3, pp. 1002 1026 c 2006 Society for Industrial and Applied Mathematics ACOUSTIC SCATTERING BY MILDLY ROUGH UNBOUNDED SURFACES IN THREE DIMENSIONS SIMON N. CHANDLER-WILDE,
More informationWeek 2 Notes, Math 865, Tanveer
Week 2 Notes, Math 865, Tanveer 1. Incompressible constant density equations in different forms Recall we derived the Navier-Stokes equation for incompressible constant density, i.e. homogeneous flows:
More information