The inverse scattering problem by an elastic inclusion
|
|
- Bennett Hilary Lawrence
- 5 years ago
- Views:
Transcription
1 The inverse scattering problem by an elastic inclusion Roman Chapko, Drossos Gintides 2 and Leonidas Mindrinos 3 Faculty of Applied Mathematics and Informatics, Ivan Franko National University of Lviv, Ukraine. 2 Department of Mathematics, National Technical University of Athens, Greece. 3 Computational Science Center, University of Vienna, Austria. arxiv:6.7376v math.ap] 24 Oct 26 Abstract In this work we consider the inverse elastic scattering problem by an inclusion in two dimensions. The elastic inclusion is placed in an isotropic homogeneous elastic medium. The inverse problem, using the third Betti s formula (direct method), is equivalent to a system of four integral equations that are non linear with respect to the unknown boundary. Two equations are on the boundary and two on the unit circle where the far-field patterns of the scattered waves lie. We solve iteratively the system of integral equations by linearising only the far-field equations. Numerical results are presented that illustrate the feasibility of the proposed method. Keywords linear elasticity, inverse scattering problem, integral equation method Introduction The inverse scattering problem consists on finding the shape and the location of an obstacle by measuring the scattered wave, close or far from the scatterer. Depending on the kind of illumination (acoustic, electromagnetic or elastic) and the properties of the obstacle (soft, hard, penetrable or not) one faces different kind of problems regarding the unique solvability of the problem and the numerical scheme for approximating the solution. In this work we place the obstacle in a two-dimensional homogeneous and isotropic elastic medium and we assume that it is penetrable, a so-called inclusion, with different Lamé parameters from the exterior domain. We consider as incident wave an elastic longitudinal or transversal wave that after interacting with the boundary of the medium is split into an interior (transmitted) and a scattered wave, propagating in the inclusion and the exterior, respectively. The scattered wave is also decomposed into a longitudinal and a transversal wave with different wavenumbers that behave like spherical waves with different polarizations at infinity. Before considering the inverse problem, we should have a good knowledge of the direct problem, which is to find the scattered field and its far-field patterns from the knowledge of the obstacle and the incident wave. The direct problem is linear and well posed for smooth obstacles 28]. The inverse problem that we consider here can be seen as a continuation of ] where the inverse problem was examined for a rigid scatterer and a cavity. The problem of detecting an elastic inclusion has been also considered for given boundary measurements, 3], using the factorization chapko@lnu.edu.ua dgindi@math.ntua.gr leonidas.mindrinos@univie.ac.at
2 method 9], the linear sampling method 3, 32], a gradient descent method 25, 29] or the probing method 8]. Here, we solve this inverse problem by formulating an equivalent system of non-linear integral equations that has to be solved with a regularization iterative scheme due to its ill-posedness. To avoid an inverse crime we consider the direct method (Betti s formula) for the inverse problem and we keep the indirect approach as proposed in 28] for the direct problem. This method was introduced in 23] and then applied in many different problems, see for instance 6, 4, 5, 27, 3] for some recent applications. The system consists of four equations, two on the unknown boundary taking advantage of the boundary conditions and two on the unit circle assuming that we know the far-field pattern of the scattered fields for one or more incident waves. Even though the first two equations are well-posed because of the equivalence to the system of integral equations for the direct problem, the last two inherit the ill-posedness of the system due to the smoothness of the far-field operators. Following 2, 7] we apply a two-step method meaning, we first solve the well-posed subsystem to obtain the corresponding densities and then we solve the linearized (with respect to the boundary) ill-posed subsystem to update the initial approximation of the radial function. We consider Tikhonov regularization and the normal equations are solved by the conjugate gradient method. The paper is organized as follows: in Section 2 we formulate the problem in two dimensions and in Section 3 we present the direct scattering problem, the elastic potential and the equivalent system of integral equations. The inverse problem is stated in Section 4 where we construct an equivalent system of integral equation using the direct method. In Section 5 the two-step method for the parametrized form of the system and the necessary Fréchet derivatives of the operators are presented. In the last section, the numerical examples give satisfactory results and demonstrate the applicability of the proposed method. 2 Problem formulation We consider the scattering of time-harmonic elastic waves by an isotropic and homogeneous elastic inclusion D i R 2 with smooth boundary Γ described by the Lamé parameters λ i, µ i and the constant density ρ i. The exterior of D i described by D e = R 2 \ D i is filled with an isotropic and homogeneous elastic medium with Lamé constants λ e, µ e and density ρ e. Henceforth, j = i, e counts for the interior D i and the exterior domain D e, respectively. In addition, we assume that λ j + µ j >, µ j > and ρ j >. By ˆτ we define the unit tangent vector to Γ and by ˆn = Q ˆτ the unit normal vector directed on D e, where Q denotes the unitary matrix Q = ( ). The incident field is either a longitudinal plane wave or a transversal plane wave u inc p (x; ˆd) = ˆd e ikp,e ˆd x, u inc s (x; ˆd) = Q ˆd e iks,e ˆd x, where ˆd is the propagation vector and the wavenumbers are given by k 2 p,j := ρ jω 2 λ j + 2µ j, k 2 s,j := ρ jω 2 µ j, where ω > is the circular frequency. In the following, α = p, s counts for the longitudinal and the transversal waves, respectively. 2
3 The scattering of u inc by the inclusion generates the scattered field u e, x D e and the transmitted field u i, x D i. Both of them satisfy the Navier equation in their domains of definition j u j + ρ j ω 2 u j =, x D j, () with the Lamé operator defined by j := µ j j + (λ j + µ j ). If u j satisfies (), due to the Helmholtz decomposition, it can be written as a sum of a longitudinal and a transversal wave which are defined by u j = u j p + u j s, u j p := kp,j 2 u j, u j s := u j u j p. On the boundary we impose transmission conditions of the form where the boundary traction operator T j is given by u i = u e + u inc, on Γ, (2a) T i u i = T e (u e + u inc ), on Γ, (2b) T j u j := λ j ˆn( u j ) + 2µ j (ˆn ) u j + µ j (Q ˆn) (Q u j). The field u e is required to satisfy also the Kupradze radiation condition ( ) u e lim r α r r ik α,eu e α =, r = x, (3) uniformly in all directions. Then, the direct elastic scattering problem reads: Given D j (geometry and elastic parameters) and the incident field u inc α, solve the boundary value problem () - (3) to obtain u j. At this point we recall that any solution of () satisfying (3) has an asymptotic behaviour of the form ( )} u e α = eikα,er {u r α ( x) + O, r, r uniformly in all directions x = x/r S, where S denotes the unit circle. The pair (u p, u s ) is called the far-field patterns of the scattered field u e. 3 The direct elastic scattering problem To represent the solution of the direct and the inverse problem as a combination of an elastic singleand a double-layer potential we first introduce the fundamental solution of the Navier equation Φ j (x, y) = i H () 4µ (k s,j x y )I + i ] j 4ρ j ω 2 H () (k s,j x y ) H () (k p,j x y ) in terms of the the identity matrix I and the Hankel function H () of order zero and of the first kind. The Green s tensor can be transformed into Φ j (x, y) = Φ,j ( x y )I + Φ 2,j ( x y )J(x y), 3
4 where the functions Φ,j, Φ 2,j : R C are given by 2] Φ,j (t) = Φ 2,j (t) = i 4µ j H () i 4ρ j ω 2 (k s,jt) 2ks,j with H () = H () and t i ] 4ρ j ω 2 k s,j H () t (k s,jt) k p,j H () (k p,jt), (4a) ] H () (k s,jt) k 2 s,jh () (k s,jt) 2k p,j t J(x) = xx x 2, x H () (k p,jt) + k 2 p,jh () (k p,jt) in terms of a dyadic product of x with its transpose x. Then, for the vector density ϕ C,a (Γ)] 2, < a, we introduce the elastic single-layer potential (S j ϕ)(x) = Φ j (x, y) ϕ(y)ds(y), x D j \Γ, (5) Γ and the elastic double-layer potential (D j ϕ)(x) = T j y Φ j (x, y) ] ϕ(y)ds(y), x Dj \Γ. (6) Γ It is well known that S j and T j xd j are continuous in R 2 but both D j and T j xs j satisfy the following jump relations 24] ( D j ϕ = ± ) 2 I + K j ϕ, on Γ, (7a) ( T j xs j ϕ = ) 2 I + L j ϕ, on Γ, (7b) (4b) T j xd j ϕ = N j ϕ, on Γ, (7c) where the upper (lower) sign corresponds to the limit x Γ from D e (D i ), and the integral operators are defined by (K j ϕ)(x) = T j y Φ j (x, y) ] ϕ(y)ds(y), x Γ, (8a) Γ (L j ϕ)(x) = T j xφ j (x, y) ϕ(y)ds(y), x Γ, (8b) Γ (N j ϕ)(x) = T j x T j y Φ j (x, y) ] ϕ(y)ds(y), x Γ. (8c) Γ All the above integrals are well defined and in particular the operator S j for x Γ is weakly singular, the operators K j, L j are singular and N j admits a hypersingular kernel. From the asymptotic behaviour of the Hankel functions we can compute also the far-field patterns of the single- (5) and double-layer potential (6) 7, 2] (S α ϕ)( x) = β α J α ( x) ϕ(y) e ikα,e x y ds(y), x S, (9a) Γ (D α ϕ)( x) = γ α J α ( x) F( x, y) ϕ(y) e ikα,e x y ds(y), x S, (9b) Γ 4
5 with the coefficients β p = eiπ/4, λ e + 2µ e 8πkp,e β s = eiπ/4 µ e 8πks,e, γ p = e iπ/4 kp,e λ e + 2µ e 8π, and the matrices J p ( x) = J( x), J s ( x) = I J( x) and γ s = e iπ/4 µ e ks,e 8π, F( x, y) = λ e x ˆn(y) + µ eˆn(y) x + µ e (ˆn(y) x)i. Considering the indirect integral equation method, we search the solution of the direct scattering problem in the form u j (x) = (D j ϕ j )(x) + (S j ψ j )(x), x D j. () To simplify the above representation, we set and the formula () is reduced to ϕ j (x) = τ j ϕ(x), ψ j (x) = ψ(x), τ j = λ j + 2µ j µ j (λ j + µ j ) u j (x) = τ j (D j ϕ)(x) + (S j ψ)(x), x D j. () Using this representation, applying the boundary conditions (2) and the jump relations (7) we see that the densities ϕ, ψ satisfy the system of integral equations ( ( ) ( I + Li L e τ i N i τ e N e ψ T = S i S e τi+τe 2 I + τ i K i τ e K e) e ) u inc Γ ϕ u inc. (2) Γ The following result regarding uniqueness and existence was proved in 28]. Theorem 3. The system of integral equations (2) has precisely one solution (ϕ, ψ), with ϕ C,a (Γ)] 2 and ψ C,a (Γ)] 2. Moreover, the corresponding displacement fields () solve the direct scattering problem () - (3). Then, the solution of the direct problem () provides us with the far-field pattern (u p, u s ) given by u α ( x) = J α ( x) τ e γ α F( x, y) ϕ(y) + β α ψ(y)] e ikα,e x y ds(y), x S, (3) Γ where we have used the asymptotic forms (9) and that ϕ, ψ are the solutions of (2). Remark 3.2 The choice of τ j is not random since, as we are going to see later, the combination τ i N i τ e N e turns out to be a weakly singular operator, reducing the hypersingularity of N j. For the numerical implementation, we present also different representations of the solutions in order to distinguish from the formulas we will derive later for the solution of the inverse problem, even though, we are going to consider the direct method. We do not consider the solvability of the above systems. Let the solution of the direct problem be given by u j (x) = (S j ψ j )(x), x D j. (4) Then, the densities satisfy the system of equations ( ) ( ) ( ) Si S e ψi u 2 I + L i 2 I L = inc Γ e ψ e T e u inc Γ 5
6 and we obtain the far-field patterns u α ( x) = (S α ψ e )( x). If we consider the representations the densities satisfy ( 2 I + K i 2 I K ) ( ) e ψi = N i N e ψ e resulting to the far-field patterns u α ( x) = (D α ψ e )( x). u j (x) = (D j ψ j )(x), x D j, (5) ( ) u inc Γ T e u inc, Γ 4 The inverse elastic scattering problem Now we can state the inverse problem, which reads: Find the shape and the position of the inclusion D i (i.e. reconstruct the boundary) from the knowledge of the far-field patterns (u p ( x), u s ( x)) for all x S, for one incident plane wave u inc α with direction ˆd either longitudinal or transverse. In general, the unique solvability of the inverse problem for one or even for a finite number of incident waves is an open problem. Uniqueness for the transmission problem exists only for infinitely many incident waves 2]. There exist also results for a rigid scatterer, local uniqueness in R 2 ] and measuring only u s for a transversal incident plane wave and simple geometries in R 3 3]. 4. The integral equation method To solve numerically this problem, we consider the non-linear integral equation method, introduced in 23], but here we apply the direct method in contrast to the forward problem. We recall for v, w C 2 (D j )] 2 the third Betti s formula (v j w w j v)dx = (v T j w w T j v)ds(x). D j Γ Using the definitions (5) and (6), we consider the above formula once for the field u e and the tensor Φ e in D e, and then for u inc, Φ e in D e to obtain u e (x) = (D e u e )(x) (S e (T e u e ))(x), x D e, (6a) = (D e u inc )(x) (S e (T e u inc ))(x), x D e. (6b) We define u t := u e + u inc and by adding (6a) and (6b) we obtain u e (x) = (D e u t )(x) (S e (T e u t ))(x), x D e. (7) Similarly, for u i, Φ i in D i, the third Betti s formula results to u i (x) = (D i u i )(x) (S i (T i u i ))(x) = (D i u t )(x) (S i (T e u t ))(x), x D i, (8) where for the last equality we have used the transmission conditions (2). We set κ = u t Γ and µ = T e u t Γ and letting x Γ in the above representations, taking the traction and considering the jump relations (7), we get u j (x) = ( 2 I ± K j)κ(x) (S j µ)(x), x Γ, T j u j (x) = ±(N j κ)(x) + ( 2 I L j)µ(x), x Γ. 6
7 We consider (2) to obtain ( 2 I K e)κ + S e µ = u inc, on Γ, (2a) ( 2 I + K i)κ S i µ =, on Γ, (2b) N e κ + ( 2 I + L e)µ = T e u inc, on Γ, (2c) N i κ + ( 2 I L i)µ =, on Γ. (2d) In addition, given the far-field operators (9) and the representation (7) of the exterior field we observe that the unknown boundary Γ and the densities satisfy the (far-field) equation ( ) ( ) ( ) D p S p u p κ µ =, on S, or in compact form D s S s u s D κ S µ = U (2) where the right-hand side is the known far-field patterns from the direct problem. We observe that we have six equations (2) and (2) for the three unknowns Γ, κ and µ. In order to take advantage of the well-posedness of the direct problem, we consider the linear combinations (2a) + (2b) and τ e (2c) + τ i (2d) for the equations on the boundary and we keep the overdetermined far-field equation. Then, we can state the following theorem as a formal formulation of the inverse problem. Theorem 4. Given an incident field u inc α, α = p or s and the far-field patterns U, for all x S, if Γ and the vector densities κ, µ satisfy the system of integral equations then, Γ solves the inverse problem. (I + K i K e )κ + (S e S i )µ = u inc Γ, (22a) (τ i N i τ e N e )κ + ( τi+τe 2 I + τ e L e τ i L i )µ = τ e T e u inc Γ, (22b) D κ S µ = U, (22c) The integral operators involved in (22) are linear with respect to the densities but non-linear with respect to the boundary Γ. The subsystem (22a) - (22b) is equivalent to (2), thus well-posed as already proved 28]. The ill-posedness of the inverse problem is then due to the smooth kernels of the far-field operators in (22c). In general, there exist three different iterative methods to solve the system (22) by linearization: A. Given initial guesses for the boundary and the densities, we linearize all three equations in order to update all the unknowns. B. Given initial guess for the boundary, we solve the subsystem (22a) - (22b) to obtain the densities. Then, keeping the densities fixed we solve the linearized equation (22c) to obtain the update for the boundary. C. Given initial guesses for the densities, we solve the far-field equation (22c) to obtain Γ and then we solve the linearized form of (22a) - (22b) to obtain the densities. The linearization, using Fréchet derivatives of the operators, and the regularization of the ill-posed equations are needed in all methods. However, the iterative method A requires the calculation of the Fréchet derivatives of the operators with respect to all the unknowns and the selection of two regularization parameters at every step. Thus, we prefer to use one of the so-called two-step methods B or C. Between the two methods, it is obvious that the second method is preferable since we solve first a well-posed linear system and then we linearize only the far-field operators (operators with smooth and simple kernels). From now on, we focus on Method B, a method introduced in 7] and then applied in different problems, see for instance 2, 26] for some recent applications. 7
8 5 The two-step method To analyse further the Method B, we consider the following parametrization for the boundary Γ = {z(t) = r(t)(cos t, sin t) : t, 2π]}, where z : R R 2 is a C 2 -smooth, 2π-periodic parametrization. We assume in addition that z is injective in, 2π), that is z (t), for all t, 2π]. The non-negative function r represents the radial distance of Γ from the origin. Then, we define and the parametrized form of (22) is given by where with the parametrized operators (A (r; ξ))(t) = ξ(t) + (A 2 (r; ξ))(t) = 2π ξ(t) = κ(z(t)), ζ(t) = µ(z(t)), t, 2π] A A(r; ξ) + B(r; ζ) = C, (23) B A = A 2, B = B 2, C = C 2, A 3 B 3 C 3 2π ( τ i T i z(t) (A 3 (r; ξ))(t) = (D (r; ξ))(t), (B (r; ζ))(t) = 2π (B 2 (r; ζ))(t) = τ i + τ e ζ(t) + 2 (B 3 (r; ζ))(t) = (S (r; ζ))(t), where Φ j (t, τ) := Φ j (z(t), z(τ)), C T i z(τ)φ i (t, τ) T e z(τ)φ e (t, τ)] ξ(τ) z (τ) dτ, ] ) T i z(τ)φ i (t, τ)] τe T e z(t) T e z(τ)φ e (t, τ) ξ(τ) z (τ) dτ, Φ e (t, τ) Φ i (t, τ)] ζ(τ) z (τ) dτ, 2π ( ) τ e T e z(t)φ e (t, τ) τ i T i z(t)φ i (t, τ) ζ(τ) z (τ) dτ, 2π (D α (r; ξ))(t) = γ α J α ( x(t)) F( x(t), z(τ)) ξ(τ) e ikα,e x(t) z(τ) z (τ) dτ, 2π (S α (r; ζ))(t) = β α J α ( x(t)) ζ(τ) e ikα,e x(t) z(τ) z (τ) dτ and the right-hand side (C (r))(t) = u inc (z(t)), (C 2 (r))(t) = τ e (T e u inc )(z(t)), C 3 (t) = U ( x(t)). Remark 5. The operators A k, B k, k =, 2, 3 act on the densities and the first variable r shows the dependence on the unknown parametrization of the boundary. Only C 3 is independent of the radial function. The two-step method for the system (23) reads as follows: 8
9 Iterative Scheme 5.2 Initially, we give an approximation of the radial function r (). Then, in the kth iteration step: i. We assume that we know r (k ) and we solve the subsystem ( ) ( ) A (r (k ) B ; ξ) + (r (k ) ; ζ) = A 2 B 2 to obtain the densities ξ (k), ζ (k). ( C C 2 ) (r (k ) ), (24) ii. Then, keeping the densities fixed, we linearize the third equation of (23), namely A 3 (r (k ) ; ξ (k) ) + (A 3(r (k ) ; ξ (k) ))(q) + B 3 (r (k ) ; ζ (k) ) + (B 3(r (k ) ; ζ (k) )(q) = C 3. (25) We solve this equation for q and we update the radial function r (k) = r (k ) + q. The iteration stops when a suitable stopping criterion is satisfied. The function q stands for the radial function of the perturbed boundary Γ q = {q(t) = q(t)(cos t, sin t) : t, 2π]}, and the Fréchet derivatives of the operators are calculated by formally differentiating their kernels with respect to r 8] ( (A (D ) ( 3(r; ξ))(q) = p (r; ξ))(q), (B (S ) (r; ξ))(q) 3(r; ζ))(q) = p (r; ζ))(q) (r; ζ))(q) with where (D s 2π ((D α (r; ξ))(q))(t) = γ α J α ( x(t)) G α ( x(t), z(τ), q(τ)) ξ(τ) e ikα,e x(t) z(τ) dτ, 2π ((S α (r; ζ))(q))(t) = β α g α ( x(t), z(τ), q(τ))j α ( x(t)) ζ(τ) e ikα,e x(t) z(τ) dτ, for v(τ) := Q q (τ) and G α ( x(t), z(τ), q(τ)) = λ e x(t)v(τ) + µ e v(τ) x(t) + µ e (v(τ) x(t))i (S s ik α,e ( x(t) q(τ)) z (τ) F( x(t), z(τ)) g α ( x(t), z(τ), q(τ)) = ik α,e ( x(t) q(τ)) z (τ) + z (τ) q (τ) z. (τ) To show injectivity of the integral operators involved in (25), we consider a simplified linearization. Assuming that z is known, we linearize with respect to z only, viewing z as independent of z, resulting to 2π ((D α (r; ξ))(q))(t) = ik α,e γ α ( x(t) q(τ))j α ( x(t)) F( x(t), z(τ)) ξ(τ)e ikα,e x(t) z(τ) dτ, 2π ((S α (r; ζ))(q))(t) = ik α,e β α ( x(t) q(τ))j α ( x(t)) ζ(τ) e ikα,e x(t) z(τ) dτ, where now ξ(τ) := ξ(τ) z (τ) and ζ(τ) := ζ(τ) z (τ). In addition, we recall that for sufficiently small q, the perturbed boundary Γ q can be represented by q(t) = q(t)q z (t), t, 2π] 6]. Now we can state the following theorem considering the above formulas for the Fréchet derivatives and as unknown the function q. 9
10 Theorem 5.3 Let ξ, ζ solve (24) and let r be the radial function of the unperturbed boundary Γ. If q C 2, 2π] satisfies the homogeneous form of equation (25), meaning then q =. (A 3(r; ξ) + B 3(r; ζ))( q) =, (26) Proof: We follow the ideas presented in 6] for the Laplace operator. Equation (26) is equivalent to (D α (r; ξ) S α (r; ζ))( q) =, α = p, s. (27) We introduce the function V (x) = 2π 2π ˆn(z(τ)) T e z(τ)φ e (x, z(τ))] ξ(τ) q(τ) z (τ) dτ ˆn(z(τ)) Φ e(x, z(τ)) ζ(τ) q(τ) z (τ) dτ, x D e, that is a radiating solution of () in D e. The far-field patterns of V are given by V α ( x) = 2π 2π z(τ) Ψ α ( x, z(τ)) ξ(τ) q(τ) z (τ) dτ z(τ) Φ α ( x, z(τ)) ζ(τ) q(τ) z (τ) dτ, x S, where Ψ α ( x, z(τ)) = γ α J α ( x) F( x, z(τ)) e ikα,e x z(τ) and Φ α ( x, z(τ)) = β α J α ( x) e ikα,e x z(τ). We observe that Vα coincide with the left hand-side of (27) since also F is independent of z. Then, Vα, and by Rellich s Lemma we get that V (x) =, x D e. In this equation the first integral has a hypersingular kernel and the second one a kernel with lower singularity. Since the fundamental solution of the Navier equation has the same (logarithmic) singularity as the fundamental solution of the Laplace equation, we can show that ξ(t) q(t) = for almost every t, 2π] 6]. An application of unique continuation and Holmgren s theorem 9] results to q =, since ξ cannot be zero on Γ. 6 Numerical implementation In this section we firstly justify numerically the convergence of the proposed scheme using analytic solutions of the direct problem and then we investigate the applicability of the Iterative Scheme 5.2 for solving the inverse problem. We solve both integral equations systems using the Nyström method. To handle the singularities of the kernels we consider the usual quadrature rules based on trigonometric interpolation 2, 22]. For smooth kernels we use the trapezoidal rule. The exact forms of the parametrized kernels are presented in 5, 7]. Thus, here we only briefly present the form of the kernel of the operator τ i N i τ e N e appearing in (2) and in A 2. This combination of operators consists of two hypersingular terms but it turns out to be weakly singular, as discussed in Remark 3.2. We consider the following decomposition T j z(t) ] ( )] T j z(τ) Φ j(t, τ) = T j z(t) T j z(τ) Φ j (t, τ) Φ () j (t, τ) + T j z(t) T j z(τ) Φ() j (t, τ)], (28) where Φ () j denotes the fundamental solution of the static (ω = ) Navier equation. The first term is weakly singular and the second one preserves the hypersingularity. The advantage, of
11 z e.6.4 z e.5.2 z i.2 z i.5 z e z i Figure. The source points considered for the different boundary curves. this decomposition, is that the second term coming from the static case is easier to handle by a Maue-type expression 7], although it is not needed here. The integral operator with kernel the second term can be written as 5, Equation 2.6] (N () j (r; ξ))(t) = = 2π T j z(t) 2π c j 2π z (t) T j z(τ) Φ() j (t, τ)] ξ(τ) z (τ) dτ cot τ t ] 2 ξ (τ) + K(t, τ) ξ(τ) dτ, for a smooth kernel K independent of D j and c j = µ j (λ j + µ j )/(λ j + 2µ j ). Since c j = τ j, we see that τ i (N () i (r; ξ))(t) τ e (N () e (r; ξ))(t) =. Then, the combination τ i N i τ e N e presents only weakly singularity due to the first term in (28). The error and convergence analysis of the proposed numerical method can be carried out based on the theory of operator approximations and on estimates for trigonometrical interpolation in Sobolev spaces 22, Section 2.4]. This analysis shows that the applied method admits superalgebraic convergence and in the case of analytical data it convergences exponentially. In the following examples, we consider three different parametrizations of the boundary curves. A peanut-shaped boundary with radial function r(t) = (.5 cos 2 t +.5 sin 2 t) /2, t, 2π], an apple-shaped boundary with radial function r(t) = and a kite-shaped boundary with parametrization cos t. sin 2t, t, 2π], +.7 cos t z(t) = (cos t +.7 cos 2t,.2 sin t), t, 2π]. 6. Example with analytic solution We consider two arbitrary points z i D i and z e D e and we define the vector-valued boundary functions f = Φ i (x, z e )] Φ e (x, z i )], on Γ, g = T i xφ i (x, z e )] T e xφ e (x, z i )], on Γ,
12 n u p,( x()) u p,( x(π/4)) i i i i i i i i φ p,( x()) φ p,( x(π/4)) i i Table. The computed and the exact longitudinal far-field for the peanut-shaped boundary considering the representation (). where ] denotes the first column of the tensor. Then, the fields u i (x) = Φ i (x, z e )], x D i, u e (x) = Φ e (x, z i )], x D e, satisfy the Navier equations () and the transmission boundary conditions u i = u e + f, on Γ, T i u i = T e u e + g, on Γ. In addition, u e satisfies the Kupradze radiation condition (3). The exact values of the far-field patterns of u e considering the asymptotic behaviour of the Hankel function are given by φ α ( x, z i ) = β α e ikα,e x z i J α ( x)], α = p, s. To compute numerically the far-field patterns we consider the three different integral representations of the solution, meaning equations (2), (4) and (5) in order to show the efficiency of the numerical scheme. Then, the densities satisfy the corresponding systems of equations where we have to replace u inc Γ by f and T e u inc Γ by g. In all examples we choose the Lamé constants to be λ e =, µ e = and ρ e = in D e and λ i = 2, µ i = 2 and ρ i = in D i and ω = 8 circular frequency. We consider the source points z i = (,.2) and z e = (.4,.6) for the peanut-shaped and the apple-shaped boundary and the points z i = (.5,.5) and z e = (,.5) for the kite-shaped boundary, see Figure. The Tables, 2 and 3 show some numerical values of the components of the far-fields patterns at given directions. We consider different representations of the solution for the different boundary parametrizations to show that our approach is applicable in all cases. We see that the exponential convergence is clearly exhibited and we obtain the correct values related to the point source located in D i. 6.2 The inverse problem To avoid an inverse crime in the following examples, the simulated far-field data were obtained by solving numerically the direct problem, replacing () by (4) and considering double amount of collocation points. We approximate the radial function q by a trigonometric polynomial of the form m m q(t) a k cos kt + b k sin kt, k= k= 2 t, 2π],
13 n u s,2( x()) u s,2( x(π/4)) i i i i i i i i φ s,2( x()) φ s,2( x(π/4)). + i i Table 2. The computed and the exact transversal far-field for the apple-shaped boundary considering the representation (5). n u p,( x(π/4)) u s,( x(π/4)) i i i i i i i i φ p,( x(π/4)) φ s,( x(π/4)) i i Table 3. The computed and the exact far-fields for the kite-shaped boundary considering the representation (4). and we consider 2n equidistant points t j = jπ/n, j =,..., 2n. The subsystem (24) is well-posed and no special treatment is required. We solve the ill-posed linearized equation (25) by minimizing the Tikhonov functional of the corresponding discretized equation ATx b λ x p p, λ >. where x R (2m+) is the vector with the unknowns coefficients a,..., a m, b,..., b m of the radial function, and A C 8n 8n, b C 8n are given by A kj = M A 3 (t k, t j ) + M B 3 (t k, t j ), b k = C 3 (t k ) (M A3 ξ)(t k ) (M B3 ζ)(t k ), for k, j =,..., 2n, where M K denotes the matrix related to the discretized kernel of the operator K. The multiplication matrix T R (8n) (2m+) stands for the trigonometric functions of the approximated radial function. Here p defines the corresponding Sobolev norm. Since q has to be real valued we actually solve the following regularized equation (T ( R(A) R(A) + I(A) I(A) ) T + λ k I p ) x = T ( R(A) R(b) + I(A) I(b) ), (29) on the kth step, where the matrix I p R (2m+) (2m+) corresponds to the Sobolev H p penalty term. We solve (29) using the conjugate gradient method. We update the regularization parameter 3
14 in each iteration step k by λ k = λ ( 2 3) k, k =, 2,... for some given initial parameter λ >. To test the stability of the iterative method against noisy data, we add also noise to the far-field patterns with respect to the L 2 norm U δ = U + δ U 2 V 2 V, for a given noise level δ, where V = V + iv 2, for V, V 2 R 8n with components normally distributed random variables. Already in the acoustic regime 2], one incident wave does not provide satisfactory results, thus we have to generalize Iterative Scheme 5.2 also for multiple illuminations u inc l, l =,..., L. Iterative Scheme 6. (Multiple illuminations) Initially, we give an approximation of the radial function r (). Then, in the kth iteration step: i. We assume that we know r (k ) and we solve the L subsystems ( ) ( ) ( ) A (r (k ) B ; ξ l ) + (r (k ) C,l ; ζ l ) = (r (k ) ), l =,..., L (3) A 2 to obtain the densities ξ (k) l, ζ (k) l. B 2 ii. Then, keeping the densities fixed, we solve the overdetermined version of the linearized third equation of (23) C 2,l A 3(r (k ) ; ξ (k) ) + B 3(r (k ) ; ζ (k) ) A 3(r (k ) ; ξ (k) 2 ) + B 3(r (k ) ; ζ (k) 2 ). q = A 3(r (k ) ; ξ (k) L ) + B 3(r (k ) ; ζ (k) L ) for q and we update the radial function r (k) = r (k ) + q. The iteration stops when a suitable stopping criterion is satisfied. 6.3 Numerical results C 3, A 3 (r (k ) ; ξ (k) ) B 3(r (k ) ; ζ (k) ) 2 ) C 3,2 A 3 (r (k ) ; ξ (k) 2 ) B 3(r (k ) ; ζ (k). C 3,L A 3 (r (k ) ; ξ (k) L ) B 3(r (k ) ; ζ (k) L ) In the following examples we choose the incident field to be a longitudinal plane wave with different incident directions given by ˆd l = (cos 2πl L 2πl, sin L ), l =,..., L. We choose the Lamé constants to be λ e =, µ e = and ρ e = in D e and λ i = 2, µ i = 3 and ρ i = in D i and ω = 8 circular frequency. We set n = 64 collocation points for the direct problem and n = 32 for the inverse. The regularized equation (29) is solved for p =, meaning H penalty term and for initial regularization parameter λ =.8. We present reconstructions for different boundary curves, different number of incident directions and initial guesses for exact and perturbed far-field data. When, we refer to noisy data, we have considered δ = 5%. In all figures the initial guess is a circle with radius r, a green solid line, the exact curve is represented by a dashed red line and the reconstructed by a solid blue line. The arrows denote the directions of the incoming incident fields. 4
15 In the first example we consider the peanut-shaped boundary. The reconstructions for m = 3 coefficients, two incident fields and r =.5 initial radius are presented in Figure 2 after 4 iterations for the exact data and 25 iterations for the noisy. In Figure 3, we see that the reconstructions are not highly dependent on the initial guess. In the second example, the boundary to be reconstructed is the apple-shaped. Here, we set m = 4, and r =.5. The reconstructions for exact data and different number of incident fields are presented in Figure 4 for 8 iterations (one incident direction) and 4 iterations (three incident directions). Figure 5 shows the effect of the initial guess for noisy data and 4 iterations. In the last example, we choose the kite-shaped boundary. We consider m = 7 coefficients and r =.5. In Figure 6 we see the improvement with respect to the number of incident fields for exact data, iterations for three illuminations and 4 iterations for four illuminations. The dependence on the initial guess is shown in Figure 7, for r = we needed 4 iterations and 25 for r =.5, in doth cases we considered noisy data. All examples show the feasibility of the proposed method that is also reasonably stable against noise. The results are considerably improved if we consider more that one incident wave. One could also considered more sophisticated regularization techniques and methods to compute the regularization parameter that could improve the reconstructions but are out of the scope of this paper. References ] G. Alessandrini, A. Morassi, and E. Rosset. Detecting an inclusion in an elastic body by boundary measurements. SIAM Review, 46(3): , 24. 2] A. Altundag and R. Kress. On a two-dimensional inverse scattering problem for a dielectric. Appl. Analysis, 9(4):757 77, 22. 3] C. J. S. Alves and N. F..M Martins. The direct method of fundamental solutions and the inverse kirsch-kress method for the reconstruction of elastic inclusions or cavities. J. Integral Equations Appl., 2:53 78, 29. 4] C. Bellis, F. Cakoni, and B. B. Guzina. Nature of the transmission eigenvalue spectrum for elastic bodies. IMA Journal of Applied Mathematics, 78(5): , 23. 5] R. Chapko. On the numerical solution of a boundary value problem in the plane elasticity for a bouble-connected domain. Math. Comput. Simulat., 66: , 24. 6] R. Chapko, O. Ivanyshyn, and O. Protsyuk. On a nonlinear integral equation approach for the surface reconstruction in semi-infinite-layered domains. Inverse Problems in Science and Engineering, 2(3):547 56, 23. 7] R. Chapko, R. Kress, and L. Mönch. On the numerical solution of a hypersingular integral equation for elastic scattering from a planar crack. IMA J. Numer. Anal., 2:6 69, 2. 8] A. Charalambopoulos. On the fréchet differentiability of boundary integral operators in the inverse elastic scattering problem. Inv. Probl., :37 6, ] A. Charalambopoulos, A. Kirsch, K. Anagnostopoulos, D. Gintides, and K. Kiriaki. The factorization method in inverse elastic scattering from penetrable bodies. Inverse Problems, 23():27 5, 27. ] D. Colton and R. Kress. Inverse acoustic and electromagnetic scattering theory, volume 93 of Applied Mathematical Sciences. Springer-Verlag, Berlin, 2 edition, 998. ] D. Gintides and L. Midrinos. Inverse scattering problem for a rigid scatterer or a cavity in elastodynamics. ZAMM Z. Angew. Math. Mech., 9(4): , 2. 2] P. Hähner and G.C. Hsiao. Uniqueness theorems in inverse obstacle scattering of elastic waves. Inv. Probl., 9: , ] G. Hu, A. Kirsch, and M. Sini. Some inverse problems arising from elastic scattering by rigid obstacles. Inv. Probl., 29():59, 23. 5
16 4] O. Ivanyshyn and B. T. Johansson. Nonlinear integral equation methods for the reconstruction of an acoustically sound-soft obstacle. J. Integral Equations Appl., 9(3):289 38, 27. 5] O. Ivanyshyn and B. T. Johansson. Boundary integral equations for acoustical inverse sound-soft scattering. J. Inv. Ill-posed Problems, 6():65 78, 28. 6] O. Ivanyshyn and R. Kress. Nonlinear integral equations for solving inverse boundary value problems for inclusions and cracks. J. Integral Equations Appl., 8():3 38, 26. 7] B. T. Johansson and B.D. Sleeman. Reconstruction of an acoustically sound-soft obstacle from one incident field and the far-field pattern. IMA J. Appl. Math., 72:96 2, 27. 8] M. Kar and M. Sini. On the inverse elastic scattering by interfaces using one type of scattered waves. J. Elast., 8():5 38, 25. 9] R. J. Knops and L. E. Payne. Uniqueness Theorems in Linear Elasticity. Springer, Berlin, 97. 2] R. Kress. On the numerical solution of a hypersingular integral equation in scattering theory. J. Comput. Appl. Math., 6(3):345 36, ] R. Kress. Inverse elastic scattering from a crack. Inv. Probl., 2: , ] R. Kress. Linear Integral Equations. Springer, New York, 3 edition, ] R. Kress and W. Rundell. Nonlinear integral equations and the iterative solution for an inverse boundary value problem. Inv. Probl., 2:27 223, ] V.D. Kupradze. Three-Dimensional problems of the mathematical theory of elasticity and thermoelasticity. North-Holland Publishing Co., New York, ] F. Le Louër. A domain derivative-based method for solving elastodynamic inverse obstacle scattering problems. Inverse Problems, 3():56, ] K.M. Lee. Inverse scattering problem from an impedance crack via a composite method. Wave Motion, 56:43 5, ] J. Li and G. Sun. A nonlinear integral equation method for the inverse scattering problem by soundsoft rough surfaces. Inverse Problems in Science and Engineering, 23(4): , ] P.A. Martin. On the scattering of elastic waves by an elastic inclusion in two dimensions. Quart. J. Mech. and Appl. Math., 43(3):275 29, ] G. Pelekanos, R.E. Kleinman, and P.M. van den Berg. Inverse scattering in elasticity a modified gradient approach. Wave Motion, 32():57 65, 2. 3] G. Pelekanos and V. Sevroglou. Inverse scattering by penetrable objects in two-dimensional elastodynamics. Journal of Computational and Applied Mathematics, 5():29 4, 23. 3] H. H. Qin and F. Cakoni. Nonlinear integral equations for shape reconstruction in the inverse interior scattering problem. Inv. Probl., 27:355, 2. 32] V. Sevroglou. The far-field operator for penetrable and absorbing obstacles in 2d inverse elastic scattering. Inverse Problems, 2(2):77 738, ] V. Sevroglou and G. Pelekanos. An inversion algorithm in two-dimensional elasticity. Journal of Mathematical Analysis and Applications, 263(): , 2. 6
17 Figure 2. Reconstruction of a peanut-shaped boundary for exact (left) and noisy (right) data Figure 3. Reconstruction of a peanut-shaped boundary for initial guess r =.5 (left) and r = (right) and exact data. 7
18 Figure 4. Reconstruction of a apple-shaped boundary for one (left) and three (right) incident fields Figure 5. Reconstruction of a apple-shaped boundary for initial guess r =.5 (left) and r =.2 (right) and noisy data. 8
19 Figure 6. Reconstruction of a kite-shaped boundary for three (left) and four (right) incident fields Figure 7. Reconstruction of a kite-shaped boundary for initial guess r = (left) and r =.5 (right) and noisy data. 9
Inverse 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 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 informationNonlinear Integral Equations for the Inverse Problem in Corrosion Detection from Partial Cauchy Data. Fioralba Cakoni
Nonlinear Integral Equations for the Inverse Problem in Corrosion Detection from Partial Cauchy Data Fioralba Cakoni Department of Mathematical Sciences, University of Delaware email: cakoni@math.udel.edu
More 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 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 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 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 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 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 informationFast iterative BEM for high-frequency scattering
Fast iterative BEM for high-frequency scattering elastodynamic problems In 3D. Marion Darbas LAMFA UMR CNRS 7352, Amiens, France S. Chaillat (POEMS, CNRS-INRIA-ENSTA) and F. Le Louër (LMAC, UTC) Conference
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 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 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 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 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 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 informationSome Old and Some New Results in Inverse Obstacle Scattering
Some Old and Some New Results in Inverse Obstacle Scattering Rainer Kress Abstract We will survey on uniqueness, that is, identifiability and on reconstruction issues for inverse obstacle scattering for
More informationDept. of Mathematics, Dept. of Mathematics, NTUA, Tel.: , Fax :
Gintides Drossos Dept. of Mathematics, Dept. of Mathematics, NTUA, Tel.: +30-1-725 017, Fax : +30-1-772 1775 email : dgindi@math.ntua.gr EDUCATION: 1987 Diploma in Mathematics, University of Patras, Greece.
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 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 informationarxiv: v3 [math.ap] 4 Jan 2017
Recovery of an embedded obstacle and its surrounding medium by formally-determined scattering data Hongyu Liu 1 and Xiaodong Liu arxiv:1610.05836v3 [math.ap] 4 Jan 017 1 Department of Mathematics, Hong
More informationInverse 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 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 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 informationLooking Back on Inverse Scattering Theory
Looking Back on Inverse Scattering Theory David Colton and Rainer Kress History will be kind to me for I intend to write it Abstract Winston Churchill We present an essay on the mathematical development
More 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 informationOn the Numerical Solution of the Laplace Equation with Complete and Incomplete Cauchy Data Using Integral Equations
Copyright 4 Tech Science Press CMES, vol., no.5, pp.99-37, 4 On the Numerical Solution of the Laplace Equation with Complete and Incomplete Cauchy Data Using Integral Equations Christina Babenko, Roman
More informationInverse wave scattering problems: fast algorithms, resonance and applications
Inverse wave scattering problems: fast algorithms, resonance and applications Wagner B. Muniz Department of Mathematics Federal University of Santa Catarina (UFSC) w.b.muniz@ufsc.br III Colóquio de Matemática
More informationA 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 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 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 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 informationAsymptotic Behavior of Waves in a Nonuniform Medium
Available at http://pvamuedu/aam Appl Appl Math ISSN: 1932-9466 Vol 12, Issue 1 June 217, pp 217 229 Applications Applied Mathematics: An International Journal AAM Asymptotic Behavior of Waves in a Nonuniform
More 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 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 informationA conformal mapping method in inverse obstacle scattering
A conformal mapping method in inverse obstacle scattering Houssem Haddar, Rainer Kress To cite this version: Houssem Haddar, Rainer Kress. A conformal mapping method in inverse obstacle scattering. Complex
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 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 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 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 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 informationDetailed Program of the Workshop
Detailed Program of the Workshop Inverse Problems: Theoretical and Numerical Aspects Laboratoire de Mathématiques de Reims, December 17-19 2018 December 17 14h30 Opening of the Workhop 15h00-16h00 Mourad
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 informationTransmission Eigenvalues in Inverse Scattering Theory
Transmission Eigenvalues in Inverse Scattering Theory Fioralba Cakoni Department of Mathematical Sciences University of Delaware Newark, DE 19716, USA email: cakoni@math.udel.edu Jointly with D. Colton,
More 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 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 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 informationStability and instability in inverse problems
Stability and instability in inverse problems Mikhail I. Isaev supervisor: Roman G. Novikov Centre de Mathématiques Appliquées, École Polytechnique November 27, 2013. Plan of the presentation The Gel fand
More informationDomain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions
Domain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions Ernst P. Stephan 1, Matthias Maischak 2, and Thanh Tran 3 1 Institut für Angewandte Mathematik, Leibniz
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 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 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 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 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 informationResearch Statement. Tonatiuh Sánchez-Vizuet
{tonatiuh@cims.nyu.edu} http://cims.nyu.edu/~tonatiuh/ I am especially interested in devising, rigorously analyzing and implementing methods for the numerical simulation of transient wave phenomena. During
More informationNEAR FIELD REPRESENTATIONS OF THE ACOUSTIC GREEN S FUNCTION IN A SHALLOW OCEAN WITH FLUID-LIKE SEABED
Georgian Mathematical Journal Volume 4 2007, Number, 09 22 NEAR FIELD REPRESENTATIONS OF THE ACOUSTIC GREEN S FUNCTION IN A SHALLOW OCEAN WITH FLUID-LIKE SEAED ROERT GILERT AND MIAO-JUNG OU Abstract. In
More informationarxiv:submit/ [math.ap] 6 Sep 2017
INVERSE OBSTACLE SCATTERING FOR ELASTIC WAVES IN THREE DIMENSIONS PEIJUN LI AND XIAOKAI YUAN arxiv:submit/1998861 [math.ap] 6 Sep 2017 Abstract. Consider an exterior problem of the three-dimensional elastic
More informationDETERMINATION OF AN UNKNOWN SOURCE TERM IN A SPACE-TIME FRACTIONAL DIFFUSION EQUATION
Journal of Fractional Calculus and Applications, Vol. 6(1) Jan. 2015, pp. 83-90. ISSN: 2090-5858. http://fcag-egypt.com/journals/jfca/ DETERMINATION OF AN UNKNOWN SOURCE TERM IN A SPACE-TIME FRACTIONAL
More informationyou expect to encounter difficulties when trying to solve A x = b? 4. A composite quadrature rule has error associated with it in the following form
Qualifying exam for numerical analysis (Spring 2017) Show your work for full credit. If you are unable to solve some part, attempt the subsequent parts. 1. Consider the following finite difference: f (0)
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 informationThe Method of Fundamental Solutions applied to the numerical calculation of eigenfrequencies and eigenmodes for 3D simply connected domains
ECCOMAS Thematic Conference on Meshless Methods 2005 C34.1 The Method of Fundamental Solutions applied to the numerical calculation of eigenfrequencies and eigenmodes for 3D simply connected domains Carlos
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 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 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 informationGraduate School of Engineering, Kyoto University, Kyoto daigaku-katsura, Nishikyo-ku, Kyoto, Japan.
On relationship between contact surface rigidity and harmonic generation behavior in composite materials with mechanical nonlinearity at fiber-matrix interface (Singapore November 2017) N. Matsuda, K.
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 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 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 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 informationarxiv: v1 [math.ap] 11 Jan 2014
THE UNIFIED TRANSFORM FOR THE MODIFIED HELMHOLTZ EQUATION IN THE EXTERIOR OF A SQUARE A. S. FOKAS AND J. LENELLS arxiv:4.252v [math.ap] Jan 24 Abstract. The Unified Transform provides a novel method for
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationGEOPHYSICAL INVERSE THEORY AND REGULARIZATION PROBLEMS
Methods in Geochemistry and Geophysics, 36 GEOPHYSICAL INVERSE THEORY AND REGULARIZATION PROBLEMS Michael S. ZHDANOV University of Utah Salt Lake City UTAH, U.S.A. 2OO2 ELSEVIER Amsterdam - Boston - London
More 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 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 informationGreen s Dyadic for the Three-dimensional Linear Elasticity in Periodic Structures
BULLETIN OF THE GREEK MATHEMATICAL SOCIETY Volume 5, 006 35 47 Green s Dyadic for the Three-dimensional Linear Elasticity in Periodic Structures Constantinos Anestopoulos, Elias Argyropoulos, Drossos Gintides
More information4.2 Green s representation theorem
4.2. REEN S REPRESENTATION THEOREM 57 i.e., the normal velocity on the boundary is proportional to the ecess pressure on the boundary. The coefficient χ is called the acoustic impedance of the obstacle
More 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 information2 Basic Equations in Generalized Plane Strain
Boundary integral equations for plane orthotropic bodies and exterior regions G. Szeidl and J. Dudra University of Miskolc, Department of Mechanics 3515 Miskolc-Egyetemváros, Hungary Abstract Assuming
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 informationElectrostatic Backscattering by Insulating Obstacles
Electrostatic Backscattering by Insulating Obstacles M. Hanke a,, L. Warth a a Institut für Mathematik, Johannes Gutenberg-Universität Mainz, 5599 Mainz, Germany Abstract We introduce and analyze backscatter
More informationSINC PACK, and Separation of Variables
SINC PACK, and Separation of Variables Frank Stenger Abstract This talk consists of a proof of part of Stenger s SINC-PACK computer package (an approx. 400-page tutorial + about 250 Matlab programs) that
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 informationASYMPTOTIC THEORY FOR WEAKLY NON-LINEAR WAVE EQUATIONS IN SEMI-INFINITE DOMAINS
Electronic Journal of Differential Equations, Vol. 004(004), No. 07, pp. 8. ISSN: 07-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) ASYMPTOTIC
More informationA posteriori error estimates for the adaptivity technique for the Tikhonov functional and global convergence for a coefficient inverse problem
A posteriori error estimates for the adaptivity technique for the Tikhonov functional and global convergence for a coefficient inverse problem Larisa Beilina Michael V. Klibanov December 18, 29 Abstract
More 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 informationMETHODS OF THEORETICAL PHYSICS
METHODS OF THEORETICAL PHYSICS Philip M. Morse PROFESSOR OF PHYSICS MASSACHUSETTS INSTITUTE OF TECHNOLOGY Herman Feshbach PROFESSOR OF PHYSICS MASSACHUSETTS INSTITUTE OF TECHNOLOGY PART II: CHAPTERS 9
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 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 informationLinear Differential Equations. Problems
Chapter 1 Linear Differential Equations. Problems 1.1 Introduction 1.1.1 Show that the function ϕ : R R, given by the expression ϕ(t) = 2e 3t for all t R, is a solution of the Initial Value Problem x =
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 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 informationThe inverse transmission eigenvalue problem
The inverse transmission eigenvalue problem Drossos Gintides Department of Mathematics National Technical University of Athens Greece email: dgindi@math.ntua.gr. Transmission Eigenvalue Problem Find (w,
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 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 informationTrefftz type method for 2D problems of electromagnetic scattering from inhomogeneous bodies.
Trefftz type method for 2D problems of electromagnetic scattering from inhomogeneous bodies. S. Yu. Reutsiy Magnetohydrodynamic Laboratory, P. O. Box 136, Mosovsi av.,199, 61037, Kharov, Uraine. e mail:
More informationReverse Time Migration for Extended Obstacles: Acoustic Waves
Reverse Time Migration for Extended Obstacles: Acoustic Waves Junqing Chen, Zhiming Chen, Guanghui Huang epartment of Mathematical Sciences, Tsinghua University, Beijing 8, China LSEC, Institute of Computational
More informationWeak solutions to anti-plane boundary value problems in a linear theory of elasticity with microstructure
Arch. Mech., 59, 6, pp. 519 539, Warszawa 2007 Weak solutions to anti-plane boundary value problems in a linear theory of elasticity with microstructure E. SHMOYLOVA 1, S. POTAPENKO 2, A. DORFMANN 1 1
More informationFast and accurate methods for the discretization of singular integral operators given on surfaces
Fast and accurate methods for the discretization of singular integral operators given on surfaces James Bremer University of California, Davis March 15, 2018 This is joint work with Zydrunas Gimbutas (NIST
More informationDirect and inverse scattering problems for inhomogeneous impedance cylinders of arbitrary shape
RADIO SCIENCE, VOL. 38, NO. 3, 1055, doi:10.109/00rs00631, 003 Direct and inverse scattering problems for inhomogeneous impedance cylinders of arbitrary shape Ibrahim Akduman Istanbul Technical University,
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 information