NEAR-FIELD IMAGING OF THE SURFACE DISPLACEMENT ON AN INFINITE GROUND PLANE. Gang Bao. Junshan Lin. (Communicated by Haomin Zhou)
|
|
- Bethanie Gregory
- 5 years ago
- Views:
Transcription
1 Volume X, No X, X, X XX Web site: NEAR-FIELD IMAGING OF THE SURFACE DISPLACEMENT ON AN INFINITE GROUND PLANE Gang Bao Department of Mathematics, Zhejiang University, Hangzhou, China; and Department of Mathematics, Michigan State University, East Lansing, MI 4884, USA Junshan Lin Institute for Mathematics and its Applications, University of Minnesota, Minneapolis, MN 55455, USA (Communicated by Haomin Zhou) Abstract This paper is concerned with the inverse diffraction problem for an unbounded obstacle which is a ground plane with some local disturbance The data is collected in the near-field regime with a distance above the surface displacement that is smaller than the wavelength In this regime, the evanescent modes carried by the scattered wave are significant, which makes it different from the far-field measurement We formulate explicitly the connection between the evanescent wave modes and the high frequency components of the surface displacement, and present a new numerical scheme to reconstruct the surface displacement from the boundary measurements By extracting the information carried by the evanescent modes effectively, it is shown that the resolution of the reconstructed image is significantly improved in the near field Numerical examples show that images with a resolution of λ/ are obtained Introduction There has been much interest in the near-field optics in the past two decades, in particular motivated by applications in the near-field microscopy The study of the near-field optics has its origin in an effort to break the diffraction limit imposed by the far-field imaging, where only the propagating wave components with spatial frequency below the wavenumber are available, and the resolution of the image is approximately λ/ [4, 5, 6] In the near field, however, the bandwidth of the spatial frequency may be expanded by taking account of the evanescent (exponentially decayed) waves We refer the reader to [3, 5, 8, 8] and the references therein for detailed discussions when single scattering (or Born approximation) is assumed In this paper, we consider the full scattering of the time harmonic electromagnetic wave that impinges on an infinite ground (the x x 3 ) plane with some local disturbance We focus the study on the TM polarization or E-parallel case, and assume that the surface displacement is invariant along the x 3 direction Consequently, the Maxwell equations are reduced to the two dimensional Helmholtz equation Mathematics Subject Classification: 35J5, 35R3, 65N Key words and phrases: inverse scattering,near-field imaging, Helmholtz equation G Bao s research was supported in part by the NSF grants DMS-9835, CCF-836, EAR-7457, DMS-96836, DMS-9, the ONR grant N4---39, a Key Project of the Major Research Plan of NSFC (No 934), and a special research grant from Zhejiang University c X AIMSciences
2 Gang Bao and Junshan Lin Figure Setup of the problem Before introducing the forward scattering model, we describe the geometry of the obstacle shown in Figure Let γ R be bounded, open, and γ be the boundary of γ Denote the closure of γ as γ The ground plane Γ := R\γ, and the local surface displacement is represented by Γ := {x = (x, x ) x γ, x = f(x )}, where the function f is defined on γ: f(x ) > for x γ, f(x ) = for x γ By requiring f > on γ, the surface displacement is directed upward Clearly, D := Γ Γ is the boundary of the whole unbounded obstacle on which the electromagnetic wave impinges The domain above D is denoted as D ( R ) The incident wave field u i = e ikq x is a plane wave that propagates along the direction q = (sin θ, cos θ) T, where θ is the incident angle and k = ω c is the wavenumber Here ω is the angular frequency, and c is the speed of the wave propagating in the vacuum We also denote the wavelength by λ If the obstacle is a flat perfect conductor, then the reflected field u r = e ikq x produced by the flat surface is a plane wave propagating along the direction q = (sin θ, cos θ) T In general, the total field u t from the scattering by Γ Γ consists of three parts: the incident wave u i, the reflected wave u r, and the scattered field The scattered field u satisfies the Helmholtz equation () u + k u = in D Assuming that the obstacle is a perfect conductor, the total field vanishes on the boundary Hence () u = (u i + u r ) on D It is easily seen that u = on Γ Moreover, the scattered field satisfies the Sommerfield radiation condition: ( ) u (3) lim r r r iku =, r = x The forward scattering problem ()-(3) admits a unique solution u C (D) C( D) if D is C and boundary data u D is continuous [] Clearly here the plane wave u i, u r C (R ) It is worth mentioning that there are many results on related scattering problems in the literature The well-posedness of the scattering problem for an obstacle with locally downward surface displacement (f(x ) < for x γ) was studied in [, ] There are also general studies on the scattering by a non-local perturbed half plane See for example [7, ] and the references therein Our goal of this paper is to study the inverse problem, more specifically near-field imaging of the local surface displacement In our framework, data is collected on
3 Near-field imaging of the surface displacement on an infinite ground plane 3 the line x = d above the surface displacement with a distance that is smaller than the wavelength λ (near-field regime) To be more precise, it is required that < d max x γ f(x ) < λ The inverse problem is to reconstruct f from the scattered field u(, d) collected on the line x = d Our work is originally motivated by the recent paper [9], in which a linearized model has been introduced for the nonlinear inverse scattering problem by the single scattering assumption The authors also proposed a broadband imaging strategy for denoising and improving the resolution of the image However, this linearized model is valid only if f λ and the modulus of its derivative f simultaneously Here we investigate the more general case by considering the full scattering model, for which the linearized model in [9] is no longer valid This imaging problem shares many of the well-known difficulties with other inverse boundary value problems, particularly nonlinearity and ill-posedness However, by collecting data in the near-field regime, the evanescent wave modes which are not accessible in the far-field regime (d max f λ) become significant This crucial fact may be confirmed by the analysis of the scattered field in Section Evanescent wave modes make it possible to break the diffraction limit It is shown that such exponentially decayed modes of the scattered wave contain exactly the high spatial frequency information (fine features) of the profile f Our study in this paper is to analyze the scattered field carefully, and design a numerical method that makes use of the evanescent modes effectively, thus to improve the resolution of the image Numerical examples confirm that a resolution of λ/ is obtained in the near field It should be pointed out that our work is different from the reconstruction of the star-like local disturbance by the far-field pattern in [3] The main focus in the near-field imaging is to break the diffraction limit In this paper the evanescent wave modes in the near-field regime are studied carefully, and a clear characterization of the connection between the high frequency components of the surface displacement and the evanescent wave modes is given explicitly ( (3) in Section ) The analysis leads naturally to a reconstruction method which extracts the information carried by the evanescent modes effectively and is very different from the far-field reconstruction in [3] Moreover, we do not impose any restriction on the parameterization of the local surface displacement The rest of the paper is organized as follows Section begins with a brief discussion on the layer potential and boundary integral equations The scattered field is studied in the near-field regime In particular, it is shown that the evanescent wave modes are significant in this regime and contain the high spatial frequency information of the profile f Based on the analysis in Section, a reconstruction method is proposed in Section 3 Several numerical examples are presented in Section 4 to demonstrate that the resolution via near-field imaging may be significantly improved compared to the far-field case The paper is concluded with some general remarks in Section 5 Analysis of the scattered wave Layer potential and boundary integral equations The boundary integral equation methods have become very popular in recent decades for the simulation of the scattering from an obstacle, and the theory is now very well developed [6, 7] We begins with a brief discussion in this section before examining the scattered wave in the near-field regime
4 4 Gang Bao and Junshan Lin Introduce Green s function G(x, y) := Φ(x, y) Φ(x r, y), where Φ(x, y) = i 4 H() (k x y ) is the fundamental solution for the Helmholtz equation in R and x r is the reflection of x by the x axis, ie, x r = (x, x ) Denote Γ = Γ γ, Γ r := { (x, x ) x Γ }, Γ Γ r = Γ Γ r γ, and D r := { (x, x ) x D } For a function ψ C(Γ), we define the single layer potential: u(x) = G(x, y)ψ(y)ds y, x R \Γ Γ r Γ The following lemma is concerned with the limit of the normal derivative of the single layer potential, when it is extended from above and below the boundary Γ The limit for the case when Γ is the smooth boundary of a bounded obstacle is well known [6, 7] Here Γ is not a closed curve For completeness, the proof of this lemma is provided in the Appendix A Lemma Assume that Γ is C For the single layer potential with continuous density ψ, the following holds: ( ) u G(x, y) (x) = ψ(y)ds y ψ(x), x Γ, ν ± Γ ( ) u where ν is the unit normal directed into D, (x) := lim h + ν(x) u(x ± ν ± hν(x)) The following representation result may serve as a starting point for the analysis of the scattered field Lemma If D is C 3, then there exists ψ C(Γ) such that the solution to ()-(3) can be expressed as the single layer potential () u(x) = G(x, y)ψ(y)ds y x D Γ Proof Let u be the solution of ()-(3) Note u = on Γ, by Green s theorem and the radiation condition, the scattered field takes the following form G(x, y) () u(x) = u(y) G(x, y) u(y) ds y, x D ν y ν y Γ Denote the domain bounded by Γ Γ r as D Let u := u i +u r Then u C (R ) satisfies the Helmholtz equation in D By Green s theorem, we have (3) ( + Γ Γ r ) ( Φ(x, y) ν y u (y) Φ(x, y) u ) (y) ds y =, x D, ν y where ν y is the unit normal directed to D for y Γ, ν y is the unit normal directed to D r for y Γ r
5 Near-field imaging of the surface displacement on an infinite ground plane 5 By noting that u (y) = u (y r ), and u (y) = u (y r ) for y Γ, we get ν y ν yr Φ(x, y) Φ(x r, y) u (y)ds y = ( u (y))ds y, Γ r ν y Γ ν y Φ(x, y) u (y) ds y = Φ(x r, y)( u (y) (4) )ds y Γ r ν y Γ ν y (5) Substituting (4) into (3) yields Γ G(x, y) u (y) G(x, y) u (y) ds y =, x D ν y ν y A combination of () and (5) leads to G(x, y) u(x) = (u + u )(y) G(x, y) Γ ν y ( u = G(x, y) + u Γ ν y ν y ( u Let ψ = ( u + u ) (y)ds y, ν y ν y ) (y)ds y x D ) ψ C(Γ) follows by the standard regularization theory + u ν y ν y for second-order elliptic equations and the Sobolev imbedding theorems [] The proof is now complete Define the integral operator K : C(Γ) C(Γ) by letting (6) (Kψ)(x) := G(x, y)ψ(y)ds y, x Γ Γ Denote g = (u i +u r ) Γ Then the existence of the solution to the integral equation Kψ = g follows from Lemma Note that since the single layer potential can be continuously extended to Γ, the density function ψ defined in Lemma is a solution of the integral equation Kψ = g Regarding the uniqueness, the following result holds Proposition If D is C 3, then Kψ = g is uniquely solvable if k is not the eigenvalue of in D for the Dirichlet problem Here D is the domain bounded by Γ Γ r For completeness, we sketch the proof here If Kψ = for some ψ C(Γ), then the single layer potential u(x) = G(x, y)ψ(y)ds y solves the exterior problem u + k u = in D, u (7) = ( ) on D, u lim r r r iku = Γ
6 6 Gang Bao and Junshan Lin and the interior problem (8) { u + k u = in D, u = on Γ Γ r If k is not the eigenvalue of in D for the Dirichlet problem, then (7) and (8) attains a unique solution respectively Hence u v u = By the limit of the v single layer potential in Lemma we have ψ = u v u = on Γ Note that v ψ C(Γ), thus ψ on Γ Scattered wave in the near-field regime Now we are ready to examine the scattered field in the near-field regime It is shown that the exponentially decayed (evanescent) wave modes which are localized to the surface the local disturbance are significant in the near field Furthermore, the connection between the evanescent wave modes and the high frequency components for the profile of the local disturbance is formulated explicitly For convenience, the Green s function G(x, y) and the fundamental solution Φ(x, y) are written with an explicit dependence on the first variable: G(x, y) = G(x y ; x, y ) := Φ(x y, x y ) Φ(x y, x + y ), where Φ(x, x ) = i 4 H() (k x + x ) For the free space fundamental solution Φ(x y, x y ), we have the following plane wave (Weyl) decomposition (Appendix B): (9) Φ(x y, x y ) = i ei(x y) κ e ik(κ) x y dκ, 4π k (κ) where k κ κ < k (propagating modes), () k (κ) = i κ k κ > k (evanescent modes) R The decomposition may be viewed as the sum of the plane waves that consist of propagating and evanescent modes It is clear that all the wave modes propagate along the x direction Along the x direction, when the magnitude of the spatial frequency κ is below k, the wave mode also propagates; otherwise, it decays exponentially along the x direction and is denoted as the evanescent mode A similar plane wave decomposition holds for G(x y ; x, y ) evaluated at x = d: () G(x y ; d, y ) = i 4π R k (κ) ei(x y) κ ( e ik(κ) d y e ik(κ) d+y ) dκ Thus for the scattered field at x = (x, d), by noting the single layer potential (), some simple calculations yield u(x, d) = i ( ei(x y) κ e ik(κ) d f(y) e ik(κ) d+f(y) ) ψ(y, f(y )) J dy dκ, 4π R γ k (κ) where J = + f This implies that the measured scattered field u on the line x = d can also be viewed as the superposition of the propagating and evanescent wave modes It is important to note that the evanescent modes with spatial
7 Near-field imaging of the surface displacement on an infinite ground plane 7 frequency beyond the wavenumber k decay exponentially along the x direction, and are localized to the surface of the obstacle within one wavelength Therefore, in the far-field regime, such evanescent modes carried by the scattered field are lost However, in the near-field regime, the evanescent modes are significant, and the measured scattered field carries more information for the profile of the local disturbance to be reconstructed On the other hand, following from the Taylor expansion of G(x y ; d, y ) at y =, the scattered field can also be expanded as u(x, d) = [ G(x y ; d, ) + G(x y ; d, ) f(y ) + G(x y ; d, ) (f(y )) y where γ +O(f 3 )] ψ(y, f(y )) J dy (a) G(x y ; d, ) = G(x y ; d, ) = by the symmetry property of G; y (b) for the high order term, by a direct calculation, asymptotically ( O(f 3 ) = O ( f ) G(x y ; d, ) λ ) f y ( ) f We assume that, and denote ϕ = ψ + f λ Then the scattered field is simplified as u(x, d) γ G(x y ; d, ) y f(y )ϕ(y )dy = G(x y ; d, ) () f(y )ϕ(y )dy R y The equality holds since the surface displacement is supported on γ Here we implicitly extend the definition of f and ϕ to R by setting f as for x R \ γ On the other hand, from () a simple calculation yields G(x y ; d, ) y = π R e i(x y) κ e ik(κ)d dκ Therefore, by taking the Fourier transform of (), we arrive at (3) û(κ, d) e ik(κ)d (fϕ)(κ), κ R, where denotes the Fourier transform Remark The expression (3) formulates explicitly the connection between the evanescent wave modes and the high frequency components of the profile f It indicates that a high spatial frequency mode of the scattered field û(κ, d) carries high spatial frequency information (fine features) of f to be reconstructed Remark If f λ and the modulus of its derivative f, then (3) is reduced automatically to the linear model discussed in [9] Now we distinguish the far-field and the near-field cases based on (3) When the spatial frequency κ > k, e ik(κ)d decreases exponentially with respect to the distance d and its value vanishes when d exceeds one wavelength (see Figure ) Thus in the far-field regime (d λ), û(κ, d) for κ > k, ie, the high spatial y
8 8 Gang Bao and Junshan Lin 8 κ= k κ=5 k d/λ Figure e ik(κ)d for d (, 5λ) when the spatial frequency κ = k and 5k respectively frequency information of f is lost in the far-field measurement In the context of imaging, this implies that it is impossible to recover f with very high resolution when any noise is present However, in the near-field regime with d < λ, e ik(κ)d is not close to and the exponentially decayed modes are still significant in the scattered field Therefore, the higher spatial frequency components of f can still be retrieved by inverting the evanescent modes of the scattered field 3 Near-field imaging 3 Inversion method Assume that the measurement u(, d) is polluted with some additive noise n(x ), which takes the following form n(x ) = σ rand(x ) u(x, d) σ is the noise level, and rand(x ) is a uniformly distributed random variable in [, ] for each x R Moreover, rand(x ) is mutually independent for different values of x Based on the previous analysis, we present a reconstruction method From (3), we introduce the pseudo-inverse operator I d as follows: { e ik (κ)d κ k () I d (κ) = c, κ > k c, where I d (κ) is a cut-off regularized operator, and k c is a regularization parameter For other regularization techniques such as the classical Tikhonov regularization and the Landweber iteration, we refer the reader to [7,,, 4, 9] for detailed discussions In the far-field case, as we discuss in the previous section, only the propagating modes can be used for imaging if noise is present, thus the cutoff frequency k c = k In the near-field regime, the bandwidth of the spatial frequency is expanded beyond the wavenumber k by taking account of the evanescent waves Note that e ik(κ)d is an exponentially increasing function with respect to κ when κ > k Hence, the noise may be exponentially amplified for large κ For fixed distance d, the cutoff frequency k c depends on the noise level (or signal-to-noise ratio) Here, following [9], we choose k c in such a way that () e ik(kc)d = e k c k d = σ
9 Near-field imaging of the surface displacement on an infinite ground plane k c /k d/λ Figure 3 k c /k versus the distance varying from λ/ to λ That is, the spatial frequency with the transfer function e ik(κ)d below the noise level σ is cut off More explicitly, ( ( log ) ) / (3) k c = k σ + d In view of () or (3), the pseudo-inverse () offers a regularization strategy for the inverse problem We plot the function k c /k for various distance d in Figure 3 at 5% noise level It is easily seen that at the fixed noise level, the cutoff frequency k c k when d < λ, ie, the bandwidth of the spatial frequency in the near field is much larger than in the far field This guarantees better resolution for the final reconstruction in the near-field regime, since the higher spatial frequency components of f are recovered Denote (4) ĥ(κ) = I d (κ)û(κ, d) To compute h, the FFT may be applied to compute the inverse Fourier transform, where h is an approximation of fϕ To reconstruct f from h, we need to take into account of the boundary data on Γ, which turns out to be a (well-posed) nonlinear problem We next introduce some notations for representing the surface displacement f Denote by C, (γ) the set of Lipschitz continuous functions on γ Introduce the Banach space C, (γ) := {f f C, (γ), f = on γ} with the usual norm f(x ) f(y ) f, = f + sup x,y γ,x y x y For a fixed small number δ, we define γ δ := {x γ dist(x, γ) < δ} For x γ δ, let x b γ such that x x b = dist(x, γ) Assume that { x, x, x x N } is a set of grid points defined on γ We represent f by a piecewise linear function, where f is linear on [x j, x j ] for j =,,, N, and is continuous on γ globally Moreover, f is strictly greater than for the interior grid points x, x x N, and is on the boundary x, x N We denote the set of all such functions by P It is clear that P is a subset of (γ), which is defined as follows: C,
10 Gang Bao and Junshan Lin C, (γ) := {f C, (γ) f(x) > for x γ; ɛ, δ >, st f(x) ɛ x x b for x γ δ } For fixed grid points { x, x, x x N }, δ = min{ x x, x N x N }, ɛ = min{ f(x ) x x, N f(x ) } On the other hand, C, x N xn (γ) is an open subset of C, (γ) We rewrite the integral operator (6) on γ by introducing the operator K f defined as (5) ( K f ϕ)(x ) := G(x y ; f(x ), f(y ))ϕ(y )dy, x γ, γ where ϕ = ψ + f The kernel of the integral operator K f has explicit dependence on f Here we adopt the subscript for K f to emphasize its dependence on f Similarly, define g f (x ) := g(x, f(x )), where g = (u i + u r ) Γ is the boundary data A natural way of separating f from h computed in (4) is to minimize the functional ( ) (6) min Kf ϕ g f + fϕ h f P L L (γ) (γ) However, this minimization problem is difficult to solve in practice, since f and ϕ are both unknowns One alternative is to solve min fϕ h L (γ) subject to K f ϕ = g f f P It should be pointed out that in general ϕ is not Fréchet differentiable with respect to f since the operator K f is compact Therefore, existing numerical methods such as Newton s method can not be applied directly We introduce a new function h such that h is Lipschitz continuous on γ, moreover { h(x ) x (7) h (x ) = γ\γ δ, for sufficiently small δ; x γ By choosing a small δ (usually the length of two neighboring grid points), h h L (γ) is small In practice, h may be chosen as the following piecewise linear function: { h (x j h(x j ) = ) j =,,, N ; j =, N Now we match the data on the boundary by solving the minimization problem with a constraint: (8) min f P Kf ϕ g f L (γ) where ϕ C(γ) satisfies fϕ = h Remark 3 The minimization problem (8) is a special case of (6) fϕ = h, (6) becomes min K L f ϕ g f + h h L (γ) (γ), and h h L (γ) is small by the definition of h By letting
11 Near-field imaging of the surface displacement on an infinite ground plane Remark 4 For f P well defined Moreover C, (γ), the function ϕ C(γ) that satisfies fϕ = h is ( min x γ\γ δ f(x ) + /ɛ) h, The problem (8) can be solved by Newton s method Since ( f λ ), the iteration is expected to converge fast to the real solution, which is confirmed by our numerical examples To linearize the problem, we require the mapping F (f) := K f ϕ g f be Fréchet differentiable with respect to f P 3 Differentiability of the mapping F (f) := K f ϕ g f Here and thereafter, M and M stand for some generic positive constants, whose values may vary from step to step but should be clear from the contexts Let L(C(γ), C(γ)) be the set of all bounded linear operators that map the functional space C(γ) to itself C, Lemma 3 If f (γ) C, (γ), then the mapping f K f is Fréchet differentiable from C, (γ) to L(C(γ), C(γ)) Moreover, the Fréchet derivative is the linear mapping δ f ( K f ) (δ f ) for δ f C, (γ), where ( K f ) (δ f ) L(C(γ), C(γ)) is defined as [( K f ) [ (δ f )]ϕ(x ) = Φ(x y, f(x ) f(y )) (δ f (x ) δ f (y )) x for ϕ C(γ) γ Φ(x y, f(x ) + f(y )) x (δ f (x ) + δ f (y )) ] ϕ(y )dy Proof It is clear that the mapping δ f ( K f ) (δ f ) is linear Step: the mapping δ f ( K f ) (δ f ) is bounded Denote r = (x y ) + (f(x ) f(y )), r = (x y ) + (f(x ) + f(y )) We first estimate the first part: Φ(x y, f(x ) f(y )) x = ik 4 ( H () (kr ) ) f(x ) f(y ) k ( 4 r H () (kr )) For a small fixed constant τ, if y x τ (away from the singularity), then ( H () )) (kr M It follows that γ\{ y x <τ } Φ(x y, f(x ) f(y )) (δ f (x ) δ f (y ))ϕ(y )dy x M δ f,
12 Gang Bao and Junshan Lin On the other hand, if y x < τ, small We have { y x <τ } { y x <τ } ( H () (kr )) O( r ) for τ sufficiently Φ(x y, f(x ) f(y )) (δ f (x ) δ f (y ))ϕ(y )dy x M r x y dy δ f, M δ f, For the second part, Φ(x y, f(x ) + f(y )) (δ f (x ) + δ f (y ))ϕ(y )dy γ x = Φ(x y, f(x ) + f(y )) [δ f (y ) δ f (x )]ϕ(y )dy γ x + Φ(x y, f(x ) + f(y )) ϕ(y )dy δ f (x ) x γ =: A + A From the estimate of the first term, the following inequality also holds: For x γ\γ δ, A M δ f, A M δ δ f,, where M δ = k For x γ δ, ( A k M ( M ( + γ\[x δ/,x +δ/] H () (kδ ) H () (kδ ) min x γ\γ δ ) ( [x δ/,x +δ/] x+δ/ ( H () (kf(x ))) H () (kr )) dy δ f (x ) (y x ) + f(x ) dy δ f (x ) δ f, + M x δ/ ) δ + ln ln( ) x x b δ f, δ + (ɛ/) x x b The last inequality follows from the fact that there exists an ɛ such that f(x ) ɛ x x b for x γ δ, and δ f (x ) δ f, x x b Therefore A M(ɛ, δ) δ f,, where M(ɛ, δ) = M ( ) (H () (kδ ) + ln inf x γ δ/ is a positive constant Therefore [( Kf ) (δ f )]ϕ M(ɛ, δ, f) δ f, δ ln( ) x x b δ + (ɛ/) x x b
13 Near-field imaging of the surface displacement on an infinite ground plane 3 for any ϕ C(γ), ie, ( K f ) (δ f ) M(ɛ, δ, f) δ f, L(C(γ),C(γ)) Step: An estimate of the remainder term K f+δf K f ( K f ) (δ f ) For any ϕ C(γ), by Taylor s formula, [ K f+δf K f ( K f ) (δ f )]ϕ(x ) Φ(x y, f(x ) f(y ) + t(δ f (x ) δ f (y )) ) = γ x dt[δ f (x ) δ f (y )] ϕ(y )dy Φ(x y, f(x ) + f(y ) + t(δ f (x ) + δ f (y )) ) γ x dt[δ f (x ) + δ f (y )] ϕ(y )dy = B + B By a similar argument as in Step, we can estimate B Assume that δ f, is sufficiently small For a fixed small constant τ, if y x τ (away from singularity), then Φ(x y, f(x ) f(y ) + t(δ f (x ) δ f (y )) ) γ\{ y x <τ } x dt[δ f (x ) δ f (y )] ϕ(y )dy M δ f, Φ If y x < τ, x O( r ) for τ sufficiently small Thus Φ(x y, f(x ) f(y ) + t(δ f (x ) δ f (y )) ) { y x <τ } x dt[δ f (x ) δ f (y )] ϕ(y )dy M { y x <τ } r x y dy δ f, M δ f, Next, the term B can also be split into two parts B and B: B Φ(x y, f(x ) + f(y ) + t(δ f (x ) + δ f (y )) ) = γ x dt[δ f (x ) δ f (y )] ϕ(y )dy and B = 4 γ Φ(x y, f(x ) + f(y ) + t(δ f (x ) + δ f (y )) ) x dt[δ f (x )δ f (y )]ϕ(y )dy It suffices to estimate B For x γ\γ δ, B Mδ δ f,, where M δ = min if δ f, is sufficiently small x γ\γ δ k ( ) H () ( k f(x ) ) + k ( ) H () f(x ) ( k f(x ) )
14 4 Gang Bao and Junshan Lin B For x γ δ, ( + γ\[x δ/,x +δ/] M δ δ f, + M [x δ/,x +δ/] x+δ/ x δ/ ) ( ) Φ dt dy δ f (x ) δ f, x M δ δ f, + M ɛ/ x x b δ f (x ) δ f, M(ɛ, δ) δ f, Therefore Kf+δf K f ( K f ) (δ f ) L(C(γ),C(γ)) = O( δ f, ), for sufficiently small δ f,, which completes the proof (y x ) + (ɛ/) x x b dy δ f (x ) δ f, Let h C, ( γ) be defined as in (7) For f P (γ), let ϕ C(γ) satisfy fϕ = h The following lemma concerns the Fréchet derivative of the mapping f ϕ ( C, (γ) C(γ) ) C, Lemma 3 If f (γ) C, (γ), then the mapping f ϕ defined as the above is Fréchet differentiable from C, (γ) to C(γ) Moreover, its Fréchet derivative is the linear mapping δ f ϕ for δ f C, (γ), where ϕ C(γ) and Proof It is easy to show that where M δ = min x γ\γ δ f(x ) C, ϕ (x ) = ϕ(x ) f(x ) δf(x ) for x γ ϕ (x ) (M δ + ɛ ) δ f, for x γ, is a constant Therefore ϕ (M δ + ɛ ) δ f, and the mapping δ f ϕ is bounded from C, (γ) C(γ) For a perturbation of f with δ f C, (γ), a perturbation of ϕ satisfies (ϕ + δ ϕ )(f + δ f ) = h If δ f, is sufficiently small, for x γ, the following estimate for the high order term holds: (ϕ + δ ϕ )(x ) ϕ(x ) ϕ (x ) = (f(x ) + δ f (x ))f(x ) ϕ(x ) δ f (x ) (M δ + ɛ ) δ f, Thus for sufficiently small δ f, (ϕ + δ ϕ ) ϕ ϕ = O( δ f, ),
15 Near-field imaging of the surface displacement on an infinite ground plane 5 Figure 4 Real part (left) and the imaginary part (right) of the scattered field By Taylor s expansion, it is easily seen that the mapping f g f (:= g(x, f(x )) ) is also Fréchet differentiable from C, (γ) to C(γ) We denote its Fréchet derivative as the mapping δ f g f Combining Lemma 3 and Lemma 3 and using the product rule, we have the following theorem:, Theorem 33 If f C (γ) C, (γ), F (f) := K f ϕ g f is Fréchet differentiable from C, (γ) to C(γ) Moreover, the Fréchet derivative maps δ f to DF (δf) = [( K) f (δ f )]ϕ + K f (ϕ ) g f 4 Numerical examples First, let us consider the solution of the forward scattering problem By Proposition, if k is not an eigenvalue of in D, to solve the forward scattering problem ()-(3) efficiently in our numerical simulation, we can firstly solve the integral equation Kψ = g and substitute ψ into () If k is an eigenvalue, then the forward scattering problem can be solved by introducing the artificial boundary (eg half circle ) and solving the problem in a bounded domain Since our focus is on inverse problem, without loss of generality, we assume that k is not an eigenvalue in our numerical examples In the following examples, an incident wave u i = e ikq x / with normal incident direction impinges on the obstacle The wavenumber k =, λ 68 cm, q = (, ) T In all the figures, the plots are rescaled with respect to the wavelength λ Example 4 The real surface displacement is represented by two bumps, each one with the size of order λ The two bumps are close to each other, and separated with distance λ/ The scattered field in the region [ 3λ, 3λ] [, 3λ] is plotted in Figure 4 Data are collected above the obstacle with distance d = λ/5 (near-field) We also assume that 5% noise is added to the simulation data It follows from (3) that k c 6k The near-field image f n and the real image f are plotted in Figure 5 (left) Though two bumps are close to each other (λ/), they are clear distinguishable Therefore, super-resolution is achieved with near-field measurements To confirm the convergence of Newton s method for solving the minimization problem (8), the relative error is shown with respect to the iteration number in Figure 6 Here the
16 6 Gang Bao and Junshan Lin y/λ 5 y/λ x/λ 6 6 x/λ Figure 5 Near-field (left) and far-field (right) images The solid line represents the real image, and the dotted line is the reconstruction Relative error Figure 6 3 Iterations Relative error with respect to the Newton iterations y/λ 5 y/λ x/λ x/λ Figure 7 Near-field (left) and far-field (right) images The solid line represents the real image, and the dotted line is the reconstruction relative error is defined as ( N j= ) f(x j ) f n(x j ) / ( N j= f(xj ) ) / The reconstruction converges fast, which leads to the real surface displacement after the first iterations
17 Near-field imaging of the surface displacement on an infinite ground plane 7 Y/λ X/λ Figure 8 Comparison of the reconstruction (dotted line) with the real profile (solid line) for the near-field case To compare with the far-field image, we collect the data again with d = 5λ and 5% noise It is obvious that k c = k in the far-field case The image f n is shown in Figure 5 (right) It is clear that the two bumps can not be distinguished, which is due to the fact that the high spatial frequency information of the surface displacement is lost Example 4 We consider a non-smooth profile in this example Two bumps are also separated with distance λ/ The measurement is polluted with 5% noise Figure 7 is the reconstructed near-field and far-field images when d = λ/5 and 5λ respectively Super-resolution is still achieved via near-field imaging, and the accuracy of the reconstruction is deteriorated when data is collected at far field Example 43 The real surface displacement is a long periodic structure (see Figure 8, solid line) Each period is a bump with size of order λ, and two neighboring bumps are separated with distance λ/ The measurement distance d again is λ/5, where 5% noise is added to the simulation data We compare the near-field image and the real image in Figure 8 The periodic structure is also accurately reconstructed with super-resolution 5 Summary and discussion We study the near-field imaging of the local surface displacement on an infinite ground plane An analysis of the scattered field is carried out, which indicates that in the far-field regime, the high spatial frequency information of the surface displacement is lost due to inaccessibility of the evanescent modes However, in the near-field, the evanescent modes become significant, for which high spatial frequency modes of the surface displacement may be recovered A numerical method is proposed in the near-field regime that yields super-resolution for the reconstructed image Numerically, the Newton iteration for the minimization problem (8) converges fast to the real solution However, no rigorous theoretical convergence analysis is presently available Another issue is denoising in the near-field regime Noise will be exponentially amplified in near-field regime The cutoff wavenumber k c, which is the critical parameter for the resolution of the reconstructed image, strongly depends on the noise level A denosing technique based on the broadband signal has recently been proposed in [9] for the linearized model For the nonlinear imaging, the problem becomes much more challenging and is completely open
18 8 Gang Bao and Junshan Lin Acknowledgments We thank Professor George C Papanicolaou for directing our attention to their paper [9] which motivates our work We also thank the referees for useful comments and suggestions Appendix A Proof of Lemma For x Γ, ɛ >, denote Γ x,ɛ := { y Γ, y x < ɛ } We first show that lim lim ds y = ɛ h + Γ x,ɛ Let Ω (the boundary of some bounded connected domain Ω) be a C closed curve such that Γ Ω Moreover, for x Γ, the unit outward normal of x on Ω coincides with ν x Let ψ on Ω, from the classical results in [6, 7], ( ) u Φ(x, y) (9) (x) = ds y ν + Ω, x Γ On the other hand, for any small fixed number ɛ >, lim h + Ω ds y = lim h + Now letting ɛ, we get () lim h + Ω From (9), (), we have () lim = ds y = Ω lim ɛ h + [ Ω\Γ x,ɛ + Φ(x, y) Ω\Γ x,ɛ Φ(x, y) Γ x,ɛ ds y + lim ] ds y ds y + lim h + lim ɛ h + Γ x,ɛ ds y = Γ x,ɛ Φ(x, y) ds y Γ x,ɛ Φ(x, y) ds y For the integral on Γ, for any fixed small number ɛ >, lim ψ(y)ds y h + Γ = lim ψ(y)ds y + lim (ψ(y) ψ(x))ds y h + Γ\Γ x,ɛ h + Γ x,ɛ + lim ds y ψ(x) h + Γ x,ɛ Φ(x, y) = ψ(y)ds y + lim (ψ(y) ψ(x))ds y Γ\Γ x,ɛ h + Γ x,ɛ + lim ds y ψ(x) h + Γ x,ɛ On the other hand, there exists some constant M that ds y Γ x,ɛ M
19 Near-field imaging of the surface displacement on an infinite ground plane 9 and Γ ψ(y) ψ(x), as ɛ for y Γ x,ɛ Therefore, by letting ɛ and noting (), we obtain Φ(x, y) lim ψ(y)ds y = ψ(y) h + ψ(x) For x Γ, ν ( ) u The proof for ν Φ(x r + hν x, y) lim ψ(y)ds y = h + Γ Γ By combining the above, we have ( ) u (x) = + Γ Γ Φ(x r, y) ψ(y) G(x, y) ψ(y)ds y ψ(x), x Γ can be carried out in the same fashion Appendix B Plane wave decomposition The fundamental solution Φ(x, x ) = i 4 H() (k x + x ) solves the equation Φ + k Φ = δ(x) in the sense of distribution, where δ is the Dirac mass at the origin Taking the Fourier transform with respect to x yields The solution ˆ Φ(κ, x ) x + (k κ )ˆ Φ(κ, x ) = δ(x ) i ˆ Φ(κ, x ) = e k i κ k κ x by noting the radiation condition The plane wave decomposition is obtained by applying the inverse Fourier transform References [] H Ammari, G Bao, and AW Wood, An integral equation method for the electromagnetic scattering from cavities, Math Meth Appl Sci 3 (), 57-7 [] H Ammari, G Bao, and AW Wood, Analysis of the electromagnetic scattering from a cavity, Japan J Indust Appl Math 9 (), 3-3 [3] H Ammari, J Garnier, and K Sølna, Resolution and stability analysis in full-aperture, linearized conductivity and wave imaging, Proc Amer Math Soc, to appear [4] M Born and E Wolf, Principles of Optics (6th ed), Cambridge University Press, 98 [5] P Carney and J Schotland, Inverse scattering for near-field microscopy, Appl Phys Lett 77 (), [6] D Colton and R Kress, Integral Equation Methods in Scattering Theory, Pure and Applied Mathematics, Wiley, New York (983) [7] D Colton and R Kress, Inverse Acoustic and Electromagnetic Scattering Theory, Applied Mathematical Sciences, Vol 93, Springer-Verlag, Berlin (998) [8] D Courjon and C Bainier, Near field microscopy and near field optics, Rep Prog Phys 57 (994), [9] G Derveaux, G Papanicolaou and C Tsogka, Resolution and denoising in near-field imaging, Inverse Problems (6),
20 Gang Bao and Junshan Lin [] H W Engl, M Hanke, and A Neubauer, Regularization of Inverse Problems, Mathematics and Its Application, Kluwer Academic Pubishers, New York (996) [] L C Evans, Partial Differential Equations, Graduate Studies in Mathematics, Vol 9, American Mathematical Society (997) [] A Kirsch, An Introduction to the Mathematical Theory of Inverse Problems, Applied Mathematical Sciences, Vol, Springer-Verlag, New York, 996 [3] R Kress and T Tran, Inverse scattering for a locally perturbed half-plane, Inverse Problems 6 (), [4] L Landweber, An iteration formula for Fredholm integral equations of the first kind, Am J Math 73 (95), [5] L Novotny and B Hecht, Principles of Nano-Optics, Cambridge University Press (6) [6] L Rayleigh, On the theory of optical images with special reference to the optical microscope, Phil Mag 5 (896), [7] F Reitich and C Turc, High-order solutions of three-dimensional roughsurface scattering problems at high-frequencies I: the scalar case, Waves Random and Complex Media 5 (5), -6 [8] J Sun, P Carney, and J Schotland, Near-field scanning optical tomography: a nondestructive method for three-dimensional nanoscale imaging, IEEE J Sel Top Quant (6), 7-8 [9] A V Tikhonov, On the solution of incorrectly formulated problems and the regularization method, Soviet Math Doklady, 4 (963), [] A Willers, The Helmholtz equation in disturbed half-spaces, Math Meth Appl Sci 9 (987), 3-33 [] B Zhang and S N Chandler-Wilde, Integral equation methods for scattering by infinite rough surfaces, Math Meth Appl Sci 6 (3), Received November ; revised June address: bao@mathmsuedu address: linxa@imaumnedu
arxiv: v1 [math.ap] 21 Dec 2018
Uniqueness to Inverse Acoustic and Electromagnetic Scattering From Locally Perturbed Rough Surfaces Yu Zhao, Guanghui Hu, Baoqiang Yan arxiv:1812.09009v1 [math.ap] 21 Dec 2018 Abstract In this paper, we
More 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationThe oblique derivative problem for general elliptic systems in Lipschitz domains
M. MITREA The oblique derivative problem for general elliptic systems in Lipschitz domains Let M be a smooth, oriented, connected, compact, boundaryless manifold of real dimension m, and let T M and T
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 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 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 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 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 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 information1 Lyapunov theory of stability
M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability
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 informationINVERSE MEDIUM SCATTERING PROBLEMS IN NEAR-FIELD OPTICS *1) 1. Introduction
Journal of Computational Mathematics, Vol.25, No.3, 2007, 252 265. INVERSE MEDIUM SCATTERING PROBLEMS IN NEAR-FIELD OPTICS *1) Gang Bao (Department of Mathematics, Michigan State University, East Lansing,
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 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 informationMath The Laplacian. 1 Green s Identities, Fundamental Solution
Math. 209 The Laplacian Green s Identities, Fundamental Solution Let be a bounded open set in R n, n 2, with smooth boundary. The fact that the boundary is smooth means that at each point x the external
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 informationPiecewise Smooth Solutions to the Burgers-Hilbert Equation
Piecewise Smooth Solutions to the Burgers-Hilbert Equation Alberto Bressan and Tianyou Zhang Department of Mathematics, Penn State University, University Park, Pa 68, USA e-mails: bressan@mathpsuedu, zhang
More informationViscosity Iterative Approximating the Common Fixed Points of Non-expansive Semigroups in Banach Spaces
Viscosity Iterative Approximating the Common Fixed Points of Non-expansive Semigroups in Banach Spaces YUAN-HENG WANG Zhejiang Normal University Department of Mathematics Yingbing Road 688, 321004 Jinhua
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 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 informationA method for creating materials with a desired refraction coefficient
This is the author s final, peer-reviewed manuscript as accepted for publication. The publisher-formatted version may be available through the publisher s web site or your institution s library. A method
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 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 informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
More informationGeneralized optical theorem for reflection, transmission, and extinction of power for electromagnetic fields
PHYSICAL REVIEW E 71, 5661 5 Generalized optical theorem for reflection, transmission, and extinction of power for electromagnetic fields D. R. Lytle II Department of Electrical and Computer Engineering,
More informationSHARP BOUNDARY TRACE INEQUALITIES. 1. Introduction
SHARP BOUNDARY TRACE INEQUALITIES GILES AUCHMUTY Abstract. This paper describes sharp inequalities for the trace of Sobolev functions on the boundary of a bounded region R N. The inequalities bound (semi-)norms
More information444/,/,/,A.G.Ramm, On a new notion of regularizer, J.Phys A, 36, (2003),
444/,/,/,A.G.Ramm, On a new notion of regularizer, J.Phys A, 36, (2003), 2191-2195 1 On a new notion of regularizer A.G. Ramm LMA/CNRS, 31 Chemin Joseph Aiguier, Marseille 13402, France and Mathematics
More informationIn this chapter we study elliptical PDEs. That is, PDEs of the form. 2 u = lots,
Chapter 8 Elliptic PDEs In this chapter we study elliptical PDEs. That is, PDEs of the form 2 u = lots, where lots means lower-order terms (u x, u y,..., u, f). Here are some ways to think about the physical
More informationAbstract. 1. Introduction
Journal of Computational Mathematics Vol.28, No.2, 2010, 273 288. http://www.global-sci.org/jcm doi:10.4208/jcm.2009.10-m2870 UNIFORM SUPERCONVERGENCE OF GALERKIN METHODS FOR SINGULARLY PERTURBED PROBLEMS
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 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 informationA new inversion method for dissipating electromagnetic problem
A new inversion method for dissipating electromagnetic problem Elena Cherkaev Abstract The paper suggests a new technique for solution of the inverse problem for Maxwell s equations in a dissipating medium.
More informationNon-uniqueness result for a hybrid inverse problem
Non-uniqueness result for a hybrid inverse problem Guillaume Bal and Kui Ren Abstract. Hybrid inverse problems aim to combine two imaging modalities, one that displays large contrast and one that gives
More information44 CHAPTER 3. CAVITY SCATTERING
44 CHAPTER 3. CAVITY SCATTERING For the TE polarization case, the incident wave has the electric field parallel to the x 3 -axis, which is also the infinite axis of the aperture. Since both the incident
More informationOn nonexpansive and accretive operators in Banach spaces
Available online at www.isr-publications.com/jnsa J. Nonlinear Sci. Appl., 10 (2017), 3437 3446 Research Article Journal Homepage: www.tjnsa.com - www.isr-publications.com/jnsa On nonexpansive and accretive
More informationTHE L 2 -HODGE THEORY AND REPRESENTATION ON R n
THE L 2 -HODGE THEORY AND REPRESENTATION ON R n BAISHENG YAN Abstract. We present an elementary L 2 -Hodge theory on whole R n based on the minimization principle of the calculus of variations and some
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 informationTwo-level multiplicative domain decomposition algorithm for recovering the Lamé coefficient in biological tissues
Two-level multiplicative domain decomposition algorithm for recovering the Lamé coefficient in biological tissues Si Liu and Xiao-Chuan Cai Department of Applied Mathematics, University of Colorado at
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 informationDouble Layer Potentials on Polygons and Pseudodifferential Operators on Lie Groupoids
Double Layer Potentials on Polygons and Pseudodifferential Operators on Lie Groupoids joint work with Hengguang Li Yu Qiao School of Mathematics and Information Science Shaanxi Normal University Xi an,
More informationON ILL-POSEDNESS OF NONPARAMETRIC INSTRUMENTAL VARIABLE REGRESSION WITH CONVEXITY CONSTRAINTS
ON ILL-POSEDNESS OF NONPARAMETRIC INSTRUMENTAL VARIABLE REGRESSION WITH CONVEXITY CONSTRAINTS Olivier Scaillet a * This draft: July 2016. Abstract This note shows that adding monotonicity or convexity
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 informationSobolev Spaces. Chapter 10
Chapter 1 Sobolev Spaces We now define spaces H 1,p (R n ), known as Sobolev spaces. For u to belong to H 1,p (R n ), we require that u L p (R n ) and that u have weak derivatives of first order in L p
More informationAPPLICATION OF THE MAGNETIC FIELD INTEGRAL EQUATION TO DIFFRACTION AND REFLECTION BY A CONDUCTING SHEET
In: International Journal of Theoretical Physics, Group Theory... ISSN: 1525-4674 Volume 14, Issue 3 pp. 1 12 2011 Nova Science Publishers, Inc. APPLICATION OF THE MAGNETIC FIELD INTEGRAL EQUATION TO DIFFRACTION
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 informationSome negative results on the use of Helmholtz integral equations for rough-surface scattering
In: Mathematical Methods in Scattering Theory and Biomedical Technology (ed. G. Dassios, D. I. Fotiadis, K. Kiriaki and C. V. Massalas), Pitman Research Notes in Mathematics 390, Addison Wesley Longman,
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 informationEuler Equations: local existence
Euler Equations: local existence Mat 529, Lesson 2. 1 Active scalars formulation We start with a lemma. Lemma 1. Assume that w is a magnetization variable, i.e. t w + u w + ( u) w = 0. If u = Pw then u
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 informationNonlinear error dynamics for cycled data assimilation methods
Nonlinear error dynamics for cycled data assimilation methods A J F Moodey 1, A S Lawless 1,2, P J van Leeuwen 2, R W E Potthast 1,3 1 Department of Mathematics and Statistics, University of Reading, UK.
More 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 informationCreating materials with a desired refraction coefficient: numerical experiments
Creating materials with a desired refraction coefficient: numerical experiments Sapto W. Indratno and Alexander G. Ramm Department of Mathematics Kansas State University, Manhattan, KS 66506-2602, USA
More informationNon-uniqueness result for a hybrid inverse problem
Non-uniqueness result for a hybrid inverse problem Guillaume Bal and Kui Ren Abstract. Hybrid inverse problems aim to combine two imaging modalities, one that displays large contrast and one that gives
More informationElectromagnetic fields and waves
Electromagnetic fields and waves Maxwell s rainbow Outline Maxwell s equations Plane waves Pulses and group velocity Polarization of light Transmission and reflection at an interface Macroscopic Maxwell
More informationOn the simplest expression of the perturbed Moore Penrose metric generalized inverse
Annals of the University of Bucharest (mathematical series) 4 (LXII) (2013), 433 446 On the simplest expression of the perturbed Moore Penrose metric generalized inverse Jianbing Cao and Yifeng Xue Communicated
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 informationSOLUTION OF THE DIRICHLET PROBLEM WITH L p BOUNDARY CONDITION. Dagmar Medková
29 Kragujevac J. Math. 31 (2008) 29 42. SOLUTION OF THE DIRICHLET PROBLEM WITH L p BOUNDARY CONDITION Dagmar Medková Czech Technical University, Faculty of Mechanical Engineering, Department of Technical
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 informationA higher-order Boussinesq equation in locally non-linear theory of one-dimensional non-local elasticity
IMA Journal of Applied Mathematics 29 74, 97 16 doi:1.193/imamat/hxn2 Advance Access publication on July 27, 28 A higher-order Boussinesq equation in locally non-linear theory of one-dimensional non-local
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 informationAn inverse problem for eddy current equations
An inverse problem for eddy current equations Alberto Valli Department of Mathematics, University of Trento, Italy In collaboration with: Ana Alonso Rodríguez, University of Trento, Italy Jessika Camaño,
More informationTakens embedding theorem for infinite-dimensional dynamical systems
Takens embedding theorem for infinite-dimensional dynamical systems James C. Robinson Mathematics Institute, University of Warwick, Coventry, CV4 7AL, U.K. E-mail: jcr@maths.warwick.ac.uk Abstract. Takens
More informationNonlinear stabilization via a linear observability
via a linear observability Kaïs Ammari Department of Mathematics University of Monastir Joint work with Fathia Alabau-Boussouira Collocated feedback stabilization Outline 1 Introduction and main result
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 informationFrequency functions, monotonicity formulas, and the thin obstacle problem
Frequency functions, monotonicity formulas, and the thin obstacle problem IMA - University of Minnesota March 4, 2013 Thank you for the invitation! In this talk we will present an overview of the parabolic
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 informationLARGE DEVIATIONS OF TYPICAL LINEAR FUNCTIONALS ON A CONVEX BODY WITH UNCONDITIONAL BASIS. S. G. Bobkov and F. L. Nazarov. September 25, 2011
LARGE DEVIATIONS OF TYPICAL LINEAR FUNCTIONALS ON A CONVEX BODY WITH UNCONDITIONAL BASIS S. G. Bobkov and F. L. Nazarov September 25, 20 Abstract We study large deviations of linear functionals on an isotropic
More informationand finally, any second order divergence form elliptic operator
Supporting Information: Mathematical proofs Preliminaries Let be an arbitrary bounded open set in R n and let L be any elliptic differential operator associated to a symmetric positive bilinear form B
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 15 Heat with a source So far we considered homogeneous wave and heat equations and the associated initial value problems on the whole line, as
More informationINSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES. Note on the fast decay property of steady Navier-Stokes flows in the whole space
INSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES Note on the fast decay property of stea Navier-Stokes flows in the whole space Tomoyuki Nakatsuka Preprint No. 15-017 PRAHA 017 Note on the fast
More informationITERATIVE APPROXIMATION OF SOLUTIONS OF GENERALIZED EQUATIONS OF HAMMERSTEIN TYPE
Fixed Point Theory, 15(014), No., 47-440 http://www.math.ubbcluj.ro/ nodeacj/sfptcj.html ITERATIVE APPROXIMATION OF SOLUTIONS OF GENERALIZED EQUATIONS OF HAMMERSTEIN TYPE C.E. CHIDUME AND Y. SHEHU Mathematics
More informationThe Dirichlet s P rinciple. In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation:
Oct. 1 The Dirichlet s P rinciple In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation: 1. Dirichlet s Principle. u = in, u = g on. ( 1 ) If we multiply
More informationXiyou Cheng Zhitao Zhang. 1. Introduction
Topological Methods in Nonlinear Analysis Journal of the Juliusz Schauder Center Volume 34, 2009, 267 277 EXISTENCE OF POSITIVE SOLUTIONS TO SYSTEMS OF NONLINEAR INTEGRAL OR DIFFERENTIAL EQUATIONS Xiyou
More informationGlobal well-posedness for the critical 2D dissipative quasi-geostrophic equation
Invent. math. 167, 445 453 (7) DOI: 1.17/s-6--3 Global well-posedness for the critical D dissipative quasi-geostrophic equation A. Kiselev 1, F. Nazarov, A. Volberg 1 Department of Mathematics, University
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 informationUncertainty principles for far field patterns and applications to inverse source problems
for far field patterns and applications to inverse source problems roland.griesmaier@uni-wuerzburg.de (joint work with J. Sylvester) Paris, September 17 Outline Source problems and far field patterns AregularizedPicardcriterion
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 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 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 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 information