Reverse Time Migration for Extended Obstacles: Acoustic Waves
|
|
- Sandra Copeland
- 6 years ago
- Views:
Transcription
1 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 Mathematics, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing 9, China Abstract. We consider the resolution of the single frequency reverse time migration (RTM) method for extended targets without the assumption of the validation of geometric optics approximation. The resolution analysis, which applies in both penetrable and non-penetrable obstacles with sound soft or impedance boundary condition on the boundary of the obstacle, implies that the imaginary part of the crosscorrelation imaging functional is always positive and thus may have better stability properties. Numerical experiments are included to illustrate the powerful imaging quality and to confirm our resolution results.
2 RTM for Inverse Scattering Problems. Introduction In this paper we propose and study an imaging algorithm to find the support of an unknown obstacle embedded in a known background medium from a knowledge of scattered acoustic waves measured on a given acquisition surface which is assumed to be far away from the obstacle. The algorithm does not require any a priori information of the physical properties of the obstacle such as penetrable or non-penetrable, and for non-penetrable obstacles, the type of boundary conditions on the boundary of the obstacle. Let the obstacle occupy a bounded Lipschitz domain R with ν the unit outer normal to its boundary. We assume the incident wave is emitted by a point source located at x s on a closed surface Γ s away from the obstacle. For penetrable obstacles, the measured wave u is the solution of the following acoustic scattering problem: u + k n(x)u = δ xs (x) in R, (.) u r r iku as r, r = x, (.) where k > is the wave number, n(x) L () is a positive scalar function supported in, δ xs is the irac source located at x s. The condition (.) is the well-known Sommerfeld radiation condition. For non-penetrable obstacles, the measured wave u is the solution of the following scattering problem: u + k u = δ xs (x) in R, (.3) u u = or + ikη(x)u = on, (.) u r r iku as r, r = x, (.5) where η(x) is a bounded function on. The irichlet boundary condition u = on corresponds to the sound soft obstacle. The second condition on in (.) is the impedance condition and it reduces to the sound hard obstacle when η(x) =. The existence and uniqueness of the problem (.)-(.) such that u s = u u i in Hloc (R ) and the problem (.3)-(.5) such that u s = u u i in Hloc (R \ ) is well-known [9, 7, ], where u i (x) = i H() (k x x s ) is the fundamental solution of the Helmholtz equation, H () (z) is the Hankle function of the first type and order zero. The direct methods for solving inverse scattering problems have drawn considerable interest in the literature in recent years. One example is the MUltiple SIgnal Classification (MUSIC) method which was first proposed in the signal processing in Schmidt []. It was used for imaging point scatterers under Born approximation for well-resolved targets by evaney [] and extended to locate small inclusions in Bruhl, Hanke, and Vogelious [5] and Ammari []. The key ingredient in the MUSIC algorithm is to construct a basis function for the signal space via singular value decomposition (SV) of the multi-static response matrix (MSR). For extended targets or cracks, the singular values of MSR matrix may decrease continuously without significant gap between signal space and null (noise) space.
3 RTM for Inverse Scattering Problems 3 The other class of direct methods for inverse scattering problems includes the linear sampling method of Colton and Kirsch [8], the factorization method of Kirsch in [], and the point source method of Potthast [8]. These methods are motivated by the uniqueness proof of inverse scattering problems and are based on the fact that the point source incident wave becomes singular as the source approaches to the boundary of scatterers (see Colton and Kress []). We refer to Potthast [9] and Kirsch and Grinberg [5] for more information on this class of methods. The reverse time migration (RTM) or the closely related prestack depth migration methods are nowadays widely used in exploration geophysics (see e.g., Berkhout [], Claerbout [7]). It is originated in the simple setting of exploding reflector model. For imaging the complex medium in practical applications, the analysis of the migration method is usually based on the high frequency assumption so that the geometric optics approximation can be used (see e.g., Bleistein, Cohen and Stockwell [3]). Our work is inspired by the recent study of RTM method in a noisy environment in Garnier [3] in which the Helmholtz-Kirchhoff identity is used to study the resolution of RTM for imaging small inclusions. Helmholtz-Kirchhoff identity plays an important role in the inverse source problem, see Bojarski []. The purpose of this paper is to provide a new mathematical understanding of the RTM method for extended obstacles without the assumption of geometric optics approximation. We study the resolution of the RTM method for both penetrable and non-penetrable obstacles. As the outcome of our resolution analysis we propose to use the imaginary part of the cross-correlation RTM functional. We show that this new imaging functional enjoys the nice feature that it is always positive and thus may have better stability properties. We will extend the results in this paper to study the RTM method for imaging extend targets using electromagnetic and elastic waves in forthcoming papers. The rest of this paper is organized as follows. In section we introduce the RTM algorithm. In section 3 we study the resolution of the imaging algorithm in section for both the penetrable and non-penetrable obstacles. In section we report extensive numerical experiments to show the competitive performance of the new RTM algorithm.. The reverse time migration method In this section we introduce the RTM method for inverse scattering problems. We assume that there are N s transducers uniformly distributed on Γ s = B s and N r transducers uniformly distributed on Γ r = B r. B s and B r are the circles of radius R s and R r, respectively. We denote by Ω the sampling domain in which the obstacle is sought. We assume the obstacle Ω and Ω is inside in B s, B r. Let G(x, y) be the fundamental solution of the Helmholtz equation G(x, y) + k G(x, y) = δ y (x) in R. Let u i (x, x s ) = G(x, x s ) be the incident wave and u s (x r, x s ) = u(x r, x s ) u i (x r, x s ) be the scattered field measured at x r, where u(x, x s ) is the solution of the problem either
4 RTM for Inverse Scattering Problems (.)-(.) or (.3)-(.5). The matrix (u s (x r, x s )) N r N s i,j= is called the multi-static response matrix (MSR) in the literature. Our RTM algorithm consists of two steps. The first step is the back-propagation in which we back-propagate the complex conjugated (time reversed) data u(x r, x s ) into the domain. The second step is the cross-correlation in which we compute the imaginary part of the cross-correlation of the incident field and the back-propagated field. Algorithm. (Reverse time migration algorithm) Given the data u s (x r, x s ) which is the measurement of the scattered field at x r when the source is emitted at x s, s =,..., N s and r =,..., N r. Back-propagation: For s =,..., N s, compute the solution v b of the following problem: v b (x, x s ) + k v b (x, x s ) = Γ N r r u s (x r, x s )δ xr (x) in R, (.) N r r= vb r r ikv b as r. (.7) Cross-correlation: For z Ω, compute { } I(z) = k Γ s N s Im u i (z, x s )v b (z, x s ). (.8) N s s= Taking the imaginary part of the correlation of the incidence field and the backpropagated field is motivated by the resolution analysis in the next section where we show that the imaginary part of the correlation functional is a positive function and thus is more stable than the real part of the correlation functional. By using the fundamental solution we can represent the solution v b of (.)-(.7) as which implies I(z) = k Im v b (z, x s ) = Γ r N r { Γ s Γ r N s N r N r r= N s N r s= r= G(z, x r )u s (x r, x s ), } G(z, x s )G(z, x r )u s (x r, x s ) z Ω. (.9) This formula is used in all our numerical experiments in section. Noticing that for z Ω which is a subdomain of Ω s, G(z, x s ) is a smooth function in x s Γ s. Similarly, G(z, x r ) is smooth in x r Γ r. We also know that since u s = u u i is the scattering solution of (.)-(.) or (.3)-(.), u s (x r, x s ) is also smooth in x r, x s. Therefore, the imaging functional I(z) in (.9) is a good trapezoid quadrature approximation of the following continuous functional: Î(z) = k Im G(z, x s )G(z, x r )u s (x r, x s )ds(x s )ds(x r ) z Ω. (.) Γ r Γ s This formula is the starting point of our resolution analysis in the next section.
5 RTM for Inverse Scattering Problems 5 3. The resolution analysis In this section we study the resolution of the imaging functional in (.). We start by recalling the the Helmholtz-Kirchhoff identity []. Lemma 3. Let be a bounded Lipschitz domain in R with ν being the unit outer normal to its boundary; we have ( ) G(ξ, y) G(x, ξ) G(x, ξ) G(ξ, y) ds(ξ) = i ImG(x, y) x, y. Proof. For the sake of completeness, we include a proof here. For any fixed x, since Im G(x, y) satisfies the Helmholtz equation, we obtain by the integral representation formula that for any y, Im G(x, ξ) G(ξ, y) Im G(x, y) = G(ξ, y) Im G(x, ξ) ds(ξ). By Im G(x, ξ) = (G(x, ξ) G(x, ξ)) we know the lemma follows if one can prove the i following identity: ( G(ξ, y) G(x, ξ) ) G(x, ξ) G(ξ, y) ds(ξ) = x, y. (3.) enote by B R a circle of radius R > that includes. Since x, y, G(x, ) and G(, y) satisfy Helmholtz equation in B R \. By integration by parts, we obtain G(ξ, y) G(x, ξ) G(x, ξ) G(ξ, y) ds(ξ) G(ξ, y) G(x, ξ) = G(x, ξ) G(ξ, y) ds(ξ) B R [ ] G(ξ, y) G(x, ξ) = ikg(ξ, y) G(x, ξ) ikg(x, ξ) G(ξ, y) ds(ξ). B R This shows (3.) by letting R since G(x, ξ) = O( ξ / ) and G(x,ξ) ikg(x, ξ) = O( ξ 3/ ) as ξ. This completes the proof. A direct consequence of the Helmholtz-Kirchhoff identity is the following lemma which plays a key role in our resolution analysis, see also [3]. Lemma 3. We have k G(x, x s )G(x s, z)ds(x s ) = Im G(x, z) + w s (x, z) x, z Ω, (3.) Γ s k G(x, x r )G(x r, z)ds(x r ) = Im G(x, z) + w r (x, z) x, z Ω, (3.3) Γ r where w s (x, z) + x w s (x, z) CRs, w r (x, z) + w r (x, z) CRr uniformly for any x, z Ω.
6 RTM for Inverse Scattering Problems Proof. The estimate for w s (x, z) in (3.) is a direct consequence of Lemma 3. and the following well-known asymptotic relations: G(x, x s ) = O(Rs / G(x, x s ) ), ikg(x, x s ) = O(Rs 3/ ) x Ω, x s Γ s. The estimate for x w s (x, z) follows again from Lemma 3. by using the following asymptotic relations: G(x, x s ) = O(Rs / G(x, xs ) ), ikg(x, x s ) = O(Rs 3/ ), x j x j for any x Ω, x s Γ s, j =,. Equation (3.3) can be proved similarly. This completes the proof. To proceed we recall the definition of the irichlet-to-neumann mapping Λ : H / ( ) H / ( ). For any g H / ( ), Λ (g) = w, where w Γ Hloc (R \ ) is the solution of the following scattering problem: w + k w = on R \, (3.) w w = g on, r r ikw as r. (3.5) The far field pattern w (ˆx) of the solution w to the scattering problem (3.)-(3.5) is defined by the asymptotic behavior { } w(x) = eik x w x / (ˆx) + O, x, (3.) x where ˆx = x/ x S := {x R : x = }. Lemma 3.3 Let g H / ( ) and w be the radiation solution satisfying (3.)-(3.5); then Im g, Λ (g) Γ = k w (ˆx) dˆx >, S where, Γ represents the duality pairing between H / ( ) and H / ( ). Proof. For the sake of completeness, we sketch a proof here. Let B R be a circle of radius R that includes. By integrating by parts one easily obtains g, Λ(g) Γ = w w ds = ( w k w )dx + w w B R \ Γ R r ds. Thus, by using the Sommerfeld radiation condition, we have Im g, Λ(g) Γ = lim Im R This completes the proof by using (3.). Γ R w w r ds = k lim R Γ R w ds. The following stability estimate of the forward scattering problem problem will be used in our resolution analysis. The proof can be found in Zhang [] by using the method of limiting absorption principle where a general transmission problem in two locally perturbed half-spaces is studied.
7 RTM for Inverse Scattering Problems 7 Lemma 3. Let n(x) L () be a positive scalar function supported in and the source f(x) L (R ) have compact support. Then the scattering problem U + k n(x)u = f(x) in R, U r r iku as r, r = x, admits a unique solution U H loc (R ). Moreover, there exists a constant C > such that U H () C f L (R ). We remark that the constant C in the lemma depends on the scatterer and the wave number k. The estimate of the explicit dependence of C on the wave number k when k is large draws considerable interests in the literature. We refer to the recent work of Chandler-Wilde and Monk [] and the references therein for results in that direction. In this paper we are interested in the case when the wavelength λ = π/k is of comparable size of the scatterer, that is, the case when k is not very large. Now we are in the position to show the resolution theorem for the RTM algorithm in this paper. We first consider the case of penetrable obstacles. Theorem 3. For any z Ω, let ψ(x, z) be the radiation solution of the Helmholtz equation with penetrable scatterer : ψ + k n(x)ψ = k (n(x) )Im G(x, z) in R. (3.7) Then if the measured field u s = u u i with u satisfying the problem (.)-(.), we have Î(z) = k ψ (ˆx, z) dˆx + wî(z) z Ω, S where wî L (Ω) C(Rs + Rr ). Proof. Since G(x, x s ) satisfies G(x, x s ) + k G(x, x s ) = δ xs (x) in R, we know that u s = u u i satisfies u s (x, x s ) + k u s (x, x s ) = k ( n(x))u(x, x s ), which implies u s satisfies the Lippmann-Schwinger equation u s (x, x s ) = k (n(ξ) )u(ξ, x s )G(x, ξ)dξ. Then, by Lemma 3., v b (z, x s ) = G(z, x r )u s (x r, x s )ds(x r ) Γ r = k (n(ξ) )u(ξ, x s ) G(x r, ξ)g(z, x r )dξds(x r ) Γ r [ ] = k (n(ξ) ) k (Im G(ξ, z) + w r(ξ, z)) u(ξ, x s )dξ, which yields after using (.) that Î(z) = k Im G(z, x s )v b (z, x s )ds(x s ) Γ s [ ] = k Im k (n(ξ) ) k (Im G(ξ, z) + w r(ξ, z)) v(ξ, z)dξ, (3.8)
8 RTM for Inverse Scattering Problems 8 where v(ξ, z) = k Γ s G(z, x s )u(ξ, x s )ds(x s ). Since by the Lippmann-Schwinger equation u(ξ, x s ) = G(ξ, x s ) + k (n(y) )u(y, x s )G(ξ, y)dy, we obtain by using Lemma 3. that v(ξ, z) = (Im G(ξ, z) + w s (ξ, z)) + k (n(y) )v(y, z)g(ξ, y)dy. Let w(ξ, z) = v(ξ, z) (Im G(ξ, z) + w s (ξ, z)); then w(ξ, z) = k (n(y) ) [w(y, z) + (Im G(y, z) + w s (y, z))] G(ξ, y)dy, and, consequently, w(ξ, z) = [ ] k (n(y) ) w(y, z) + (Im G(y, z) + w s (y, z)) G(ξ, y)dy. This implies that w(ξ, z) is the radiation solution of the following Helmholtz equation [ ] ξ w(ξ, z) + k w(ξ, z) = k (n(ξ) ) w(ξ, z) + (Im G(ξ, z) + w s (ξ, z)), which is equivalent to ξ w(ξ, z) + k n(ξ)w(ξ, z) = k (n(ξ) )(Im G(ξ, z) + w s (ξ, z)). Now by (3.7) we know that ζ(ξ, z) = w(ξ, z) ψ(ξ, z) satisfies ξ ζ(ξ, z) + k n(ξ)ζ(ξ, z) = k (n(ξ) )w s (ξ, z) in R, and the Sommerfeld radiation condition. By Lemma 3. and Lemma 3. we have ζ(, z) H () Ck (n( ) ) L () w s (, z) L () CR s, (3.9) uniformly for z Ω. This implies by using Lemma 3. again that v(ξ, z) = w(ξ, z) + (Im G(ξ, z) + w s (ξ, z)) = ψ(ξ, z) + ζ(ξ, z) + (Im G(ξ, z) + w s (ξ, z)), where ζ(, z) H () + w s (, z) L () CRs. Now by (3.8) we obtain Î(z) = Im k (n(ξ) )(Im G(ξ, z) + w r (ξ, z))v(ξ, z)dξ = Im k (n(ξ) )Im G(ξ, z)ψ(ξ, z)dξ + O(Rs + Rr ). Now by the equation satisfied by ψ(ξ, z) and integrating by parts we obtain [ Î(z) = Im ξ ψ(ξ, z) + k n(ξ)ψ(ξ, z) ] ψ(ξ, z)dξ + O(Rs + Rr ) ψ(ξ, z) = Im ψ(ξ, z)ds(ξ) + O(Rs + Rr ) = Im ψ(, z), Λ (ψ(, z) + O(R s This completes the proof by using Lemma R r ).
9 RTM for Inverse Scattering Problems 9 We remark that since ( + k n(x))im G(x, z) = k (n(x) )Im G(x, z), ψ(x, z) can be viewed as the scattering solution of the Helmholtz equation with the incident wave Im G(x, z). It is known that Im G(x, z) = Im i H() (k x z ) = J (k x z ) which peaks when x = z and decays as x z becomes large. Noticing that the source in (3.7) is supported in because of n(x) = outside. Thus the source in (3.7) becomes small when z moves away from outside the scatterer. On the other hand, the source in (3.7) will not be small if z is inside. Therefore we expect that the imaging functional will have a strong contrast at the boundary of the scatterer and decays outside the scatterer. This is indeed confirmed in our numerical experiments. Now we consider the resolution of the imaging functional in the case of nonpenetrable obstacles. We only prove the results for the case of impedance boundary condition. The case of irichlet boundary condition is similar and simpler and is left to the interested readers. We need the following result about the forward impedance scattering problem whose proof is similar to that for partially coated perfect conductor considered in []. It can also be proved by using the method of limiting principle, see e.g., Leis [, Remark.39]. Lemma 3.5 Let η(x) be bounded and f H / ( ). Then the scattering problem U + k U = in R \, U + ikη(x)u = g on, U r r iku as r, admits a unique solution U H loc (R ). Moreover, there exists a constant C such that U H / ( ) C g H / ( ). Theorem 3. For any z Ω, let ψ(x, z) be the radiation solution of the Helmholtz equation ψ(x, z) + k ψ(x, z) = in R \, (3.) [ ] ψ(x, z) + ikη(x)ψ(x, z) = + ikη(x) Im G(x, z) on. (3.) Then if the measured field u s = u u i with u satisfying the problem (.3)-(.5) with the impedance condition in (.), we have Î(z) = k ψ (ˆx, z) dˆx + k η(ξ) ψ(ξ, z) + Im G(ξ, z) ds(ξ) + wî(z) z Ω, S where wî L (Ω) C(R s + R r ). Proof. By the integral representation we know that ( u s (x r, x s ) = u s (ξ, x s ) G(x ) r, ξ) us (ξ, x s ) G(x r, ξ) ds(ξ). (ξ) (ξ)
10 RTM for Inverse Scattering Problems By using Lemma 3. we obtain that, for any z Ω, v b (z, x s ) = G(z, x r )u s (x r, x s )ds(x r ) Γ r = [ u s (ξ, x s ) k (Im G(ξ, z) + w r(ξ, z)) us (ξ, x s ) ] (Im G(ξ, z) + w r (ξ, z)) ds(ξ). By (.) we obtain then [ Î(z) = kim v s (ξ, z) (Im G(ξ, z) + w r(ξ, z)) v s(ξ, z) (Im G(ξ, z) + w r (ξ, z)) ] ds(ξ), (3.) where v s (ξ, z) = k Γ s G(z, x s )u s (ξ, x s )ds(x s ). By taking the complex conjugate, we have v s (ξ, z) = k G(z, x s )u s (ξ, x s )ds(x s ). Γ s Thus v s (ξ, z) is a weighted superposition of the scattered waves u s (ξ, x s ). Therefore, v s (ξ, z) is the radiation solution of the Helmholtz equation ξ v s (ξ, z) + k v s (ξ, z) = in R \ satisfying the impedance boundary condition (ξ) + ikη(ξ) v s (ξ, z) = k G(z, x s ) Γ s (ξ) + ikη(ξ) u s (ξ, x s )ds(x s ) = k G(z, x s ) Γ s (ξ) + ikη(ξ) G(ξ, x s )ds(x s ) = (ξ) + ikη(ξ) (Im G(ξ, z) + w s (ξ, z)) on, where we have used Lemma 3. in the last equality. This implies by using (3.)-(3.) that v s (ξ, z) = ψ(ξ, z) + ζ(ξ, z), where ζ satisfies the impedance scattering problem in Lemma 3.5 with g( ) = ( + ikη) w s (, z). By the Lemma 3. and Lemma 3.5, ζ satisfies ζ(, z) H / ( ) CRs uniformly for z Ω. Now from the boundary condition satisfied by ζ on we know also that ζ(,z) for z Ω. H / ( ) CR s Now substituting v s (ξ, z) = ψ(ξ, z) + ζ(ξ, z) into (3.) we obtain ( ) Im G(ξ, z) ψ(ξ, z) Î(z) = Im ψ(ξ, z) Im G(ξ, z) (ξ) (ξ) = Im By (3.) we have Im ( Im G(ξ, z) ψ(ξ, z) (ξ) ( Im G(ξ, z) ψ(ξ, z) (ξ) ψ(ξ, z) Im G(ξ, z) (ξ) ) ) ψ(ξ, z) Im G(ξ, z) ds(ξ) (ξ) ds(ξ) + O(R s ds(ξ) + O(R s uniformly + R r ) + R r ).
11 RTM for Inverse Scattering Problems [ Im G(ξ, z) = Im ψ(ξ, z) ikη(ξ)im G(ξ, z) (ξ) ψ(ξ, z) ] + ikη(ξ)ψ(ξ, z) Im G(ξ, z) + ikη(ξ)im G(ξ, z)ψ(ξ, z) ds(ξ) (ξ) ( ) [ ψ(ξ, z) = Im ψ(ξ, z) ikη(ξ)ψ(ξ, z) (ξ) Im G(ξ, z) ] + + ikη(ξ)im G(ξ, z) Im G(ξ, z) + ikη(ξ)im G(ξ, z)ψ(ξ, z) ds(ξ) (ξ) = Im ψ(, z), Λ (ψ(, z)) Γ + k η(ξ) ψ(ξ, z) + Im G(ξ, z) ds(ξ). This completes the proof by using Lemma 3.3. For the ease of later reference, we state the results for the sound soft non-penetrable obstacles in the following theorem. Theorem 3.3 For any z Ω, let ψ(x, z) be the radiation solution of the Helmholtz equation ψ(x, z) + k ψ(x, z) = ψ(x, z) = Im G(x, z) on. in R \, Then if the measured field u s = u u i with u satisfying the problem (.3)-(.5) with the irichlet condition in (.), we have Î(z) = k ψ (ˆx, z) dˆx + wî(z) z Ω, S where wî L (Ω) C(Rs + Rr ). We remark that for the non-penetrable obstacles, ψ(x, z) is again the scattering solution of the Helmholtz equation with the incident wave Im G(x, z). Similar to the remark after the proof of Theorem 3., we expect that the imaging functional will peak at the boundary of the scatterer and decay away from the scatterer.. Numerical experiments In this section, we show a variety of numerical experiments to illustrate the imaging power of the RTM algorithm proposed in this paper. To synthesize the scattering data we compute the solution u(x, x s ) of the scattering problem (.)-(.) or (.3)-(.5) by standard Nyström s methods [9]. The boundary integral equations on are solved on a uniform mesh over the boundary with ten points per probe wavelength. The boundaries of the obstacles used in our numerical experiments are parameterized as follows: Circle: x = ρ cos(θ), x = ρ sin(θ), θ (, π], Kite: x = cos(θ) +.5 cos(θ).5, x =.5 sin(θ), θ (, π], p-leaf: r(θ) = +. cos(pθ), θ (, π].
12 RTM for Inverse Scattering Problems Figure. Example : Non-penetrable obstacle with irichlet condition, probe wavelength λ =, and N s = N r =. The first row shows the surface plot and the contour plot of the imaging functional in Theorem 3.3 (ignoring wî). The second row shows the surface plot and the contour plot of the imaging functional in (.). Example. In this example we consider the imaging of a circle of radius ρ = centered at the origin. We compare the results of the imaging functional in (.9) with the corresponding theoretical results in Theorems 3. and 3.3. Let Γ r = Γ s be the circle of radius R =. Let Ω = ( 3, 3) ( 3, 3) be the search region and the imaging functional is computed at the nodal points of a uniform mesh. We test two probe wavelengths λ = and λ =.5, where λ = π/k. Figure and Figure show the comparison of the imaging functional for a nonpenetrable obstacle with irichlet condition on. Figure 3 and Figure show the results for a penetrable obstacle with n(x) =.5. We observe that the imaging functional (.9) agrees well with the theoretical results in Theorems 3. and 3.3. Figure 5 shows the comparison of the real part of the cross-correlation functional { } Ĩ(z) = k Γ s Γ r N s N r Re G(z, x s )G(z, x r )u N s N s (x r, x s ) z Ω. r s= r= and the imaginary part of the correlation functional in (.9). We observe that the imaging quality, i.e., the contrast near the boundary of the scatterer, of Ĩ(z) is worse than that of the imaginary part of the cross-correlation functional (.9). Figure depicts cross-sections of the imaging functional at x = for probe wavelengths λ =,,,.5. It shows clearly the resolution is improved with the increase of the wavenumber. Example. In this example we verify the diffraction limit in our resolution analysis
13 RTM for Inverse Scattering Problems 3 Figure. Example : Non-penetrable obstacle with irichlet condition, probe wavelength λ =.5, and N s = N r = 5. The first row shows the surface plot and the contour plot of the imaging functional in Theorem 3.3 (ignoring wî). The second row shows the surface plot and the contour plot of the imaging functional in (.). by considering two objects separated with a small fixed distance. The first model is a circle of radius ρ = 3 and a kite. The distance between two objects is about.5. We use the probe wavelength λ =,,.5 to image the objects. The search domain is Ω = (, ) (, ) with a sampling mesh. The results are shown in Figure 7 from which we observe that with the increase of the probe wave number, the gap between two objects is more and more visible. Figure 8 shows the imaging results when the number of emitters and receivers is reduced. The second model is a big circle of radius ρ = 5 and a small circle of radius ρ =.5 or ρ =.. The search domain is Ω = ( 7, 7) ( 7, 7) with a sampling 3 3 mesh. Figure 9 shows the imaging results. We observe that our imaging algorithm clearly locates the boundary of obstacles with different size as long as the high wave number content is available in our received data. Example 3. In this example we consider the stability of the imaging functional with respect to the additive Gaussian random noise. We introduce the additive Gaussian noise as follows: u noise = u s + ν noise, where u s is the synthesized data and ν noise is the Gaussian white noise with mean zero and standard deviation µ times the maximum of the data u s, i.e. ν noise N (, µ max u s ). For the fixed probe wavelength λ =, we choose one kite and one 8-leaf in our test.
14 RTM for Inverse Scattering Problems Figure 3. Example : Penetrable obstacles, probe wavelength λ =, and N s = N r =. The first row shows the surface plot and the contour plot of the imaging functional in Theorem 3.3 (ignoring wî). The second row shows the surface plot and the contour plot of the imaging functional in (.). The search domain is Ω = (, ) (, ) with a sampling mesh. Figure shows the imaging results with µ = %, %, %, % noise in the single frequency scattered data. The imaging quality can be improved by using multi-frequency data as illustrated in Figure in which we show the imaging results of summing the imaging functionals of probe wavelengths λ = /.8, /.9, /., /., /.. Example. In this example we consider the imaging of non-penetrable obstacles with different impedance conditions. We take the probe wavelength λ =.5. The search domain is Ω = (, ) (, ) with a sampling mesh. Figure shows the imaging results which indicate clearly the effectiveness of our imaging algorithm. Acknowledgments The work of J. Chen is supported in part by China NSF under the grant 5, 7, and that of Z. Chen is supported in part by National Basic Research Project under the grant CB397 and China NSF under the grant. The authors are also grateful to the referees for their constructive comments.
15 RTM for Inverse Scattering Problems 5 Figure. Example : Penetrable obstacles, probe wavelength λ =.5, and N s = N r = 5. The first row shows the surface plot and the contour plot of the imaging functional in Theorem 3.3 (ignoring wî). The second row shows the surface plot and the contour plot of the imaging functional in (.). References [] Ammari H and Kang H Reconstruction of Small Inhomogeneities from Boundary Measurements Lecture Notes in Mathematics vol 8 (Berlin: Springer Verlag) [] Berkhout A J 98 Seismic Migration: Imaging of Acoustic Energy by Wave Field Extrapolation (New York: Elsevier) [3] Bleistein N Cohen J and Stockwell J Mathematics of Multidimensional Seismic Imaging, Migration, and Inversion (New York: Springer) [] Bojarski N N 973 Inverse Scattering Naval Air Systems Command Report N9-73-C-3 (Washington.C.) [5] Bruhl M, Hanke M and Vogelius M 3 A direct impedance tomography algorithm for locating small inhomogeneities Numer. Math [] Cakoni F, Colton and Monk P The direct and inverse scattering problems for partially coated obstacles Inverse Problems [7] Claerbout J F 985 Imaging the Earth s Interior (Oxford: Blackwell Scientific Publication) [8] Colton and Kirsch A 99 A simple method for solving inverse scattering problems in the resonance region Inverse Problems [9] Colton and Kress R 998 Inverse Acoustic and Electromagnetic Scattering Problems (Heidelberg: Springer) [] Colton and Kress R Using fundamental solutions in inverse scattering Inverse Problems R9-R
16 RTM for Inverse Scattering Problems.5. Exact RTM Numerical RTM.9 Exact RTM Numerical RTM Exact RTM Numerical RTM.9.8 Exact RTM Numerical RTM Exact R TM Numerical R TM.. Exact RTM Numerical RTM Exact RTM Numerical RTM.. Exact RTM Numerical RTM Figure 5. Example : The comparison of the real (left) and imaginary (right) part of the cross-correlation functional. The first and second row is for the non-penetrable obstacle with irichlet condition and the third and fourth row is for the penetrable obstacle with n(x) =.5. The first and third row show the cross sections at x = of the real and imaginary part of the cross-correlation functional when λ = and N s = N r =. The second and fourth row show the cross sections at x = of the real and imaginary part of the cross-correlation functional when λ =.5 and N s = N r = 5.
17 RTM for Inverse Scattering Problems λ=.5 λ= λ= λ= Figure. Example : Image result comparisons for different probe wavelength: plus+line for λ = ; diamond+line for λ = ; star+line for λ = ; solid+line for λ = Figure 7. Example : The first picture is the exact obstacles. The other three pictures from left to right are imaging results using probe wavelengths λ =,,.5 and N s = N r = 8, respectively. [] Chandler-Wilde S N and Monk P 8 Wave-number explicit bounds in time-harmonic scattering SIAM J. Math. Anal [] evaney A J Super-resolution processing of multi-static data using time-reversal and MUSIC to appear in J. Acoust. Soc. Am. [3] Garnier J Sensor array imaging in a noisy environment National Institute for Mathematical Sciences NIMS Lecture Note Series TP3 (South Korea) [] Kirsch A 998 Characterization of the shape of a scattering obstacle using the spectral data of the far field operator Inverse Problems 89-5
18 RTM for Inverse Scattering Problems Figure 8. Example : The first and the second picture are imaging results using fixed probe wavelength λ = with reduced number of (N s, N r ) = (, ), (3, 8), respectively Figure 9. Example : The first and third picture show exact obstacles: one big circle with ρ = 5 and one small circle with radius ρ =.5 (first) or ρ =. (third). The second and the fourth picture are imaging results using λ =.5 and N s = N r = 38, Figure. Example 3: The imaging results using data added with additive Gaussian noise and µ = %, %, %, % from left to right, respectively. The probe wavelength λ = and N s = N r = 8. [5] Kirsch A and Grinberg N 8 The Factorization Method for Inverse Problems (Oxford: Oxford University Press) [] Leis R 98 Initial Boundary Value Problems in Mathematical Physics (Stuttgart: B.G. Teubner) [7] McLean W Strongly Elliptic Systems and Boundary Integral Equations (Cambridge: Cambridge University Press) [8] Potthast R 99 A fast new method to solve inverse scattering problems Inverse Problems 73-7 [9] Potthast R Point Sources and Multipoles in Inverse Scattering Theory (Boca Raton: Chapman & Hall/CRC)
19 RTM for Inverse Scattering Problems Figure. Example 3: The imaging results using multi-frequency data added with additive Gaussian noise and µ = %, %, %, % from left to right, respectively. The probe wavelengths λ = /.8, /.9, /., /., /. and N s = N r = Figure. Example : The first picture is the Neumann boundary condition η(x) =, the second picture is the impedance condition with η(x) =, and the third picture is the coated obstacle with the impedance η(x) = on the upper half part and η(x) = on the lower half part of the boundary of the obstacle. The probe wavelength λ =.5 and N s = N r = 8. [] Schmidt R 98 Multiple emitter location and signal parameter estimation IEEE Trans. Antennas. Propag [] Zhang B 99 On transmission problems for wave propagation in two locally perturbed half-spaces Math. Proc. Camb. Phil. Soc
An eigenvalue method using multiple frequency data for inverse scattering problems
An eigenvalue method using multiple frequency data for inverse scattering problems Jiguang Sun Abstract Dirichlet and transmission eigenvalues have important applications in qualitative methods in inverse
More 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 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 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 informationarxiv: v1 [math.na] 16 Jan 2015
PHASELESS IMAGING BY REVERSE TIME MIGRATION: ACOUSTIC WAVES ZHIMING CHEN AND GUANGHUI HUANG arxiv:5.957v [math.na] 6 Jan 5 Abstract. We propose a reliable direct imaging method based on the reverse time
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationA time domain probe method for threedimensional rough surface reconstructions
A time domain probe method for threedimensional rough surface reconstructions Article Accepted Version Burkard, C. and Potthast, R. (2009) A time domain probe method for three dimensional rough surface
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 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 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 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 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 informationThe Support of a Far Field
The Support of a Far Field John Sylvester University of Washington sylvest@math.washington.edu http://www.math.washington.edu/ sylvest/ Gratefully acknowledge research support from ONR and NSF Inverse
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 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 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 informationInverse scattering problem with underdetermined data.
Math. Methods in Natur. Phenom. (MMNP), 9, N5, (2014), 244-253. Inverse scattering problem with underdetermined data. A. G. Ramm Mathematics epartment, Kansas State University, Manhattan, KS 66506-2602,
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 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 informationHigh-Order Corrected Trapezoidal Quadrature Rules for Functions with a Logarithmic Singularity in 2-D
High-Order Corrected Trapezoidal Quadrature Rules for Functions with a Logarithmic Singularity in 2-D Juan C. Aguilar 1 Yu Chen 2 April 24, 2002 Abstract In this report we construct correction coefficients
More informationUniqueness in an inverse acoustic obstacle scattering problem for both sound-hard and sound-soft polyhedral scatterers
INSTITUTE OF PHYSICS PUBLISHING Inverse Problems 22 (2006) 515 524 INVERSE PROBLEMS doi:10.1088/0266-5611/22/2/008 Uniqueness in an inverse acoustic obstacle scattering problem for both sound-hard and
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 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 informationThe mathematics of scattering by unbounded, rough, inhomogeneous layers
The mathematics of scattering by unbounded, rough, inhomogeneous layers Simon N. Chandler-Wilde a Peter Monk b Martin Thomas a a Department of Mathematics, University of Reading, Whiteknights, PO Box 220
More informationELECTROMAGNETIC SCATTERING FROM PERTURBED SURFACES. Katharine Ott Advisor: Irina Mitrea Department of Mathematics University of Virginia
ELECTROMAGNETIC SCATTERING FROM PERTURBED SURFACES Katharine Ott Advisor: Irina Mitrea Department of Mathematics University of Virginia Abstract This paper is concerned with the study of scattering of
More informationVolume and surface integral equations for electromagnetic scattering by a dielectric body
Volume and surface integral equations for electromagnetic scattering by a dielectric body M. Costabel, E. Darrigrand, and E. H. Koné IRMAR, Université de Rennes 1,Campus de Beaulieu, 35042 Rennes, FRANCE
More 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 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 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 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 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 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 informationImage amplitudes in reverse time migration/inversion
Image amplitudes in reverse time migration/inversion Chris Stolk 1, Tim Op t Root 2, Maarten de Hoop 3 1 University of Amsterdam 2 University of Twente 3 Purdue University TRIP Seminar, October 6, 2010
More informationExtended Foldy Lax Approximation on Multiple Scattering
Mathematical Modelling and Analysis Publisher: Taylor&Francis and VGTU Volume 9 Number, February 04, 85 98 http://www.tandfonline.com/tmma http://dx.doi.org/0.3846/3969.04.893454 Print ISSN: 39-69 c Vilnius
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 informationThe Linear Sampling Method and the MUSIC Algorithm
CODEN:LUTEDX/(TEAT-7089)/1-6/(2000) The Linear Sampling Method and the MUSIC Algorithm Margaret Cheney Department of Electroscience Electromagnetic Theory Lund Institute of Technology Sweden Margaret Cheney
More informationThe Factorization Method for the Reconstruction of Inclusions
The Factorization Method for the Reconstruction of Inclusions Martin Hanke Institut für Mathematik Johannes Gutenberg-Universität Mainz hanke@math.uni-mainz.de January 2007 Overview Electrical Impedance
More informationThe Inside-Outside Duality for Scattering Problems by Inhomogeneous Media
The Inside-Outside uality for Scattering Problems by Inhomogeneous Media Andreas Kirsch epartment of Mathematics Karlsruhe Institute of Technology (KIT) 76131 Karlsruhe Germany and Armin Lechleiter Center
More informationEXISTENCE OF GUIDED MODES ON PERIODIC SLABS
SUBMITTED FOR: PROCEEDINGS OF THE FIFTH INTERNATIONAL CONFERENCE ON DYNAMICAL SYSTEMS AND DIFFERENTIAL EQUATIONS June 16 19, 2004, Pomona, CA, USA pp. 1 8 EXISTENCE OF GUIDED MODES ON PERIODIC SLABS Stephen
More informationarxiv: v2 [math.ap] 9 Sep 2014
Localized and complete resonance in plasmonic structures Hoai-Minh Nguyen and Loc Hoang Nguyen April 23, 2018 arxiv:1310.3633v2 [math.ap] 9 Sep 2014 Abstract This paper studies a possible connection between
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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationAnalysis of the anomalous localized resonance
Hyundae Lee(Inha University, Korea) Joint work with Habib Ammari, Giulio Ciraolo, Hyeonbae Kang, Graeme Milton. UCI, June 22, 2012 A conference on inverse problems in honor of Gunther Uhlmann Outline Introduction
More informationElectromagnetic wave scattering by many small bodies and creating materials with a desired refraction coefficient
168 Int. J. Computing Science and Mathematics, Vol. Electromagnetic wave scattering by many small bodies and creating materials with a desired refraction coefficient A.G. Ramm Mathematics epartment, Kansas
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 informationSeismic inverse scattering by reverse time migration
Seismic inverse scattering by reverse time migration Chris Stolk 1, Tim Op t Root 2, Maarten de Hoop 3 1 University of Amsterdam 2 University of Twente 3 Purdue University MSRI Inverse Problems Seminar,
More informationNumerical computation of the Green s function for two-dimensional finite-size photonic crystals of infinite length
Numerical computation of the Green s function for two-dimensional finite-size photonic crystals of infinite length F.Seydou 1,Omar M.Ramahi 2,Ramani Duraiswami 3 and T.Seppänen 1 1 University of Oulu,
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 informationUNCERTAINTY PRINCIPLES FOR INVERSE SOURCE PROBLEMS, FAR FIELD SPLITTING AND DATA COMPLETION
4 5 6 7 8 9 4 5 6 7 8 9 4 5 6 7 8 9 4 5 6 7 8 9 4 UNCERTAINTY PRINCIPLES FOR INVERSE SOURCE PROBLEMS, FAR FIELD SPLITTING AND DATA COMPLETION ROLAND GRIESMAIER AND JOHN SYLVESTER Abstract. Starting with
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 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 informationAnalysis of the Scattering by an Unbounded Rough Surface
Analysis of the Scattering by an Unbounded ough Surface Peijun Li and Jie Shen Abstract This paper is concerned with the mathematical analysis of the solution for the wave propagation from the scattering
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 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 informationOverview Avoiding Cycle Skipping: Model Extension MSWI: Space-time Extension Numerical Examples
Overview Avoiding Cycle Skipping: Model Extension MSWI: Space-time Extension Numerical Examples Guanghui Huang Education University of Chinese Academy of Sciences, Beijing, China Ph.D. in Computational
More informationProceedings of the 5th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15th July 2005
Proceedings of the 5th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15th July 2005 SOME INVERSE SCATTERING PROBLEMS FOR TWO-DIMENSIONAL SCHRÖDINGER
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationPropagation of Smallness and the Uniqueness of Solutions to Some Elliptic Equations in the Plane
Journal of Mathematical Analysis and Applications 267, 460 470 (2002) doi:10.1006/jmaa.2001.7769, available online at http://www.idealibrary.com on Propagation of Smallness and the Uniqueness of Solutions
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 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 informationImaging with Ambient Noise II
Imaging with Ambient Noise II George Papanicolaou Stanford University Michigan State University, Department of Mathematics Richard E. Phillips Lecture Series April 21, 2009 With J. Garnier, University
More informationIDENTIFICATION OF SMALL INHOMOGENEITIES: ASYMPTOTIC FACTORIZATION
IDENTIFICATION OF SMALL INHOMOGENEITIES: ASYMPTOTIC FACTORIZATION HABIB AMMARI, ROLAND GRIESMAIER, AND MARTIN HANKE Abstract. We consider the boundary value problem of calculating the electrostatic potential
More informationSeismic Waves and Earthquakes A Mathematical Overview
Seismic Waves and Earthquakes A Mathematical Overview The destruction caused by earthquakes is caused by the seismic waves propagating inside and in particular, near the surface of earth. We describe the
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 information