RECONSTRUCT LAMÉ PARAMETERS OF LINEAR ISOTROPIC ELASTICITY SYSTEM. Ru-Yu Lai

Transcription

1 RECONSTRUCT LAMÉ PARAMETERS OF LINEAR ISOTROPIC ELASTICITY SYSTEM Ru-Yu Lai arxiv: v3 [math.ap] 6 Aug 13 Abstract. Transient Elastography enables detection characterization of tissue abnormalities. In this paper we assume that the displacements are modeled by linear isotropic elasticity system the tissue displacement has been obtained by the first step in hybrid methods. Then we reconstruct the Lamé parameters of this system from knowledge of tissue displacement. We show that for a sufficiently large number of solutions of the elasticity system for an open set of the well-chosen boundary conditions, λ, can be uniquely stably reconstructed. The set of well-chosen boundary conditions is characterized in terms of appropriate complex geometrical optics solutions. 1. Introduction Medical imaging is the technique process used to create images of the human body for clinical purposes or medical science. Each available imagine method has its advantages disadvantages. Medical imaging modalities such as Computerized Tomography CT, Magnetic Resonance Imaging MRI Ultrasound Imaging UI are examples of modalities providing high resolution. In some situations, these modalities fail to exhibit a sufficient contrast between different types of tissues. For instance, in breast imaging ultrasound provides high resolution, while suffers from a low contrast. Other modalities, based on optical, elastic, or electrical properties of these tissues, display high contrast, such as Optical Tomography OT Electrical Impedance Tomography EIT. For example, some breast tumors on early stages might have no contrast with the healthy tissues with respect to ultrasound propagation, but a huge contrast in their optical electric properties. In order to obtain better image, the natural idea is to try to combine different imaging modalities. These are coupled-physics imaging methods, also called hybrid methods. It is to combine the high resolution modality with another high contrast modality. Examples of possible physical couplings include: Photo-Acoustic Tomography PAT Thermo- Acoustic Tomography TAT, Ultrasound modulated Optical Tomography UMOT elasticity with ultrasound in Transient Elastography TE. Reconstructions in hybrid inverse problems involve two steps. The first step is to solve the inverse problem concerning the high-resolution-low-contrast modality. For instance, in PAT TAT, this corresponds to the reconstructing the initial condition of a wave equation from boundary measures. In Transient Elastography, this is to solving an inverse scattering problem in a time-dependent wave equation. In this paper, we assume that this first step has been performed. We will focus on the second step. In the second step, we assume that the first step has been done we try to reconstruct the coefficients that display high contrast from the mappings obtained during the first step. In this paper, the modality we consider is Transient Elastography. TE is a non-invasive tool for measuring liver stiffness. The device creases high-resolution shear stiffness images of human tissue for diagnostic purposes. Shear stiffness is targeted because shear wave speed is larger in abnormal tissue than in normal tissue. In the experiment tissue initially 1991 Mathematics Subject Classification. Primary: 35R3 ; Secondary: 35Q74. Key words phrases. Inverse Problems, elasticity system, Lamé parameters. Department of Mathematics, University of Washington, Seattle, WA 98195, USA. rylai@uw.edu. 1

2 is excited with pulse at the boundary. This pulse creates the shear wave passing through the liver tissue. Then the tissue displacement is measured by using the ultrasound. The displacement is related to the tissue stiffness because the soft tissue has larger deformation than stiff tissue. When we have tissue displacement, we want to reconstruct shear modulus the first parameter λ. See [1] references there for more details. In TE, the high resolution modality is also ultrasound. The tissue displacement data can be obtained by the ultrasound in the first step. The second step is to recover the Lamé parameters from the knowledge of the tissue displacement. In the following paper, we will assume the first step has been performed. Let Ω R n, n =,3, be an open bounded domain with smooth boundary. Let u be the displacement satisfying the linear isotropic elasticity system { λ ui +S u+k u = in Ω, 1 u = g on Ω, where SA = A + A T / is the symmetric part of the matrix A. Here λ, are Lamé parameters k R is the frequency. Assume that k is not the eigenvalue of the elasticity system. The set of internal functions obtained by ultrasound in TE is given by : Hx = u j x 1 j J for some positive integer J. DenotethatP = {λ, C 7 Ω C 9 Ω; < m λ C7 Ω, C 9 Ω M λ, > }. Let H H be two sets of internal data of the elasticity system with parameters λ, λ,, respectively. Below is our main result: Theorem 1.1. Let Ω be an open bounded domain of R n with smooth boundary. Suppose that the Lamé parameters λ, λ, P Ω = Ω. Let u j ũ j be the solutions of the elasticity system with boundary data g j for parameters λ, λ,, respectively. Let H = u j 1 j J H = ũ j 1 j J be the corresponding internal data for λ, λ,, respectively for some integer J 3n+1. Then there is an open set of the boundary data g j 1 j J such that if H = H implies λ, = λ, in Ω. Moreover, we have the stability estimate in Ω CΩ + λ λ CΩ C H H C Ω. The remainder of this paper is organized as follows. In section, we introduce the CGO solutions of the elasticity system. Section 3 is devoted to constructing the Lamé parameters in two-dimensional case. The reconstruction of λ, in 3D is presented in section 4.. Complex geometric optics solutions of the elasticity system In this section, we will briefly introduce the complex geometric optics CGO solutions of the elasticity system. Based on [9], we can derive the following reduced system. Let w,f T satisfy w f f +V 1 x w w +V x f =. Here V contains the third derivative of 1/ 1 + 3/ k log V 1 x = Then the solution of the elasticity system 1 is u := 1/ w+ 1 f f 1. λ+ λ+ 1/.

3 3 Here g denotes the Hessian matrix g/ x i x j. Note that we will not need the explicit form of V in the construction of CGO solution. The construction of CGO solutions of with linear phase was first deduced by Nakamura Uhlmann in [13] [14], where they introduced the intertwining property in hling the first order term. Eskin Ralston [7] also gavea similar result in. Later Uhlmann Wang[17] used Carleman estimate to deduce the CGO solutions in two-dimensional case. In the following sections, we will use the CGO solutions constructed in Eskin Ralston s paper [7]. For the rest of this section, we will introduce the key lemma from [7]. To construct CGO solutions of, it is convenient to work on R n instead of Ω. Since Ω is bounded, we pick a ball B R for R >, such that Ω B R extend λ to R n by preserving its smoothness also suppλ, supp B R. Let α β be two orthogonal unit vectors in R n. Denote that = τα + iβ C n θ = α+iβ with τ >. Eskin Ralston [7] showed the following important result in three-dimensional case. For n =, it still holds. Lemma.1. Eskin-Ralston Consider the Schrödinger equation with external Yang-Mills potentials 3 Lu = u iax u+bxu =, x B R R n where Ax = A 1 x,...,a n x C n B R, n 6 with A j x Bx are n + 1 n+1 matrices. Then there are solutions of 3 of the form where C C n B R is the solution of u = e i x C x,θgθ x+oτ 1 iθ x C x,θ = θ AxC x,θ with detc is never zero, gz is an arbitrary polynomial in complex variables z, Oτ 1 is bounded by C1/τ in H l B R, l n. By the lemma above, the CGO solutions of can be written as follows w = e f i x r +Oτ s 1, where r, s T is C x,θgθ x. We can write w = e i x r + Oτ 1 f = e i x s + Oτ 1. Plugw,f T intou = 1/ w+ 1 f f 1. Thenwehavethe CGOsolutions of the elasticity system. Note that r,s T satisfies the equation 4 r θ = V s 1 x n n θ θ T r s Remark 1. Since C is invertible at every point in Ω, we can conclude that r, s T is not zero everywhere in Ω. This does not imply that both r s never vanish in Ω. However, for any point y Ω, there is a small neighborhood B y of y in Ω a CGO solution of w,f T such that both r s do not vanish in B y. 3. Reconstruction of Lamé parameters in two-dimensional case In the previous section, we already have the CGO solutions of the elasticity system. Now wewanttousethemtogiveareconstructionoffirst. Letu = u 1,u T bethedisplacement which satisfies the elasticity system 5 λ ui +S u+k u =. We will recover λ separately in the following two subsections..

4 Reconstruction of in D. From 5, we can deduce the following equation 6 Here we denote 7 u = 1 u u u u, F = u F +u G = k u. λ+ λ+ 1 λ+ λ+, u = a+b a b 1 a+b a b, G = with a = u u, b = 1 u 1 u, u = u 1 +u. The component u u are known since they are vectors which only depend on the internal data u. In order to recover, we want to eliminate the first term of 6 so that satisfies the transport equation. Obverse that the vector u has three different entries, we only need to construct three linearly independent vectors on some subdomain in Ω. With these three vectors, we can remove the first term of the left h side of the equation 6. More precisely, suppose that u j, for j =,1,, are three different solutions of 6 which satisfy 8 u j F +u j G = k u j. The notations Rf If are defined to be the real imaginary part of f, respectively. Now for j = 1,,3 =,,, we let u 1 = Rχxu, u = Iχxu u 1 = Rχ 3 xu 1, u = Iχ 3 xu 1, u 3 = R χ 1 xu 1 +χ xu, where χ j x is a nonzero function. Then we get 9 1 u l F +u l G = k u l, l = 1, u j F +u j G = k u j for j = 1,,3. Assume that { u 1,u,u 3 } are three linearly independent vectors in some subdomain of Ω, say Ω. Then there exist three functions Θ l 1,Θl, Θl 3 such that u l + 3 Θl j uj =. For l = 1,, multiplying 1 by Θ l j summing over j with equation 9, we have v l G = k u l + Θ l ju j, where v l = u l + 3 Θl j uj. Let β l = v l e 1 ẽ 1 +v l e ẽ γ l = e 3 +e 4 v l. Here e j R 4, ẽ j R with the j th entry is 1 others are zero. Then the above equation can be rewritten as β l +γ l = k u l + Θ l j uj Suppose that β 1 x β x are linearly independent for every x Ω. Then we can recover in Ω Ω for each frequency k independently. Lemma 3.1. Let u j for j be C solutions of the elasticity system with boundary conditions u j = g j on Ω. Let us define u = 1 u, u, u, u T assume that 1 { u 1,u,u 3 } are three linearly independent vectors in Ω, the neighborhood of x in Ω. {β 1 x,β x} are linearly independent in Ω.. 1

5 5 Then the reconstruction is stable in the sense that CΩ C x + x+ + H x H 11 x C Ω where x + Ω H x = u j j. Proof. Since {β 1 x,β x} are linearly independent in Ω, we can construct two vectorvalued functions Γx, Φx CΩ such that 1 +Γx = Φx in Ω. Since are solutions of 1 with coefficients Γ,Φ Γ, Φ, respectively, we have +Γx = Γx Γx + Φx Φx. Let x Ω, there exists a integral curve ψt with ψ = x + 13 ψt = ψe t Γψs ψ sds + t ψ1 = x. Thus Γψs Γψs +Φ Φψs ψ sds, t [,1]. Since Γ Φ only depend on Ω u j,j =,1,, u j are in the class of C Ω, the value of Γ Γ Φ Φ are bounded by the sum of α u j α ũ j for α. There is a constant C such that ψt C ψ +C H x H x C Ω, t [,1] for, P. Thus, for any x Ω, the value x is controlled by the internal data at the boundary point x Ω Global reconstruction of in D. Letλ, P,byLemma.1, itimpliesthatrs are in C 7 Oτ 1 H 5. Now we will show that how we can get three linear independent vectors of the form u. We plug the CGO solutions u = 1/ w + 1 f f 1 into u. Then we have the expression 14 u = e i x λ+ 1 r+o r+o i r i r by using the following equality which is part of the equation 4 15 s = λ+ λ+ 1/ r. +O1 Now we fix any point x Ω let = τ1,i = i C with τ >. Since, in Lemma.1, the matrix solution C x,θ is invertible, we can choose a constant vector g such that C x,θg = r,s T with the conditions s x = 1, s r in a neighborhood of x in B R, say U. Then we have the CGO solution of the elasticity system, that is, where u = 1/ w + 1 f f 1, w f = e i x r s +Oτ 1.

6 6 Let θ = /τ, θ = /τ. Let C 1 x,θ C x, θ satisfy the following two equations iθ x C 1x,θ = θ V 1 xc 1 x,θ, i θ x C x, θ = θ V 1 xc x, θ, respectively. Since = i, we can choose C x, θ = C 1 x,θ. Moreover, with suitable constant vector g, we can get that r = r 1 s = s 1 s l x =, r 1 x = 1, i = r x. By continuity of r l, we have r l in a neighborhood U of x. Then the CGO solutions are where w 1 f 1 = e i x r 1 s 1 u l = 1/ w l + 1 f l f l 1, +Oτ 1, w f = e i x r s So far we have three CGO solutions, that is, u, u 1, u. We consider e i x u 1 = 1 e i x 1 u 1 = 1 e i x 1 u = 1 Since r = r 1 s = s 1, we have e i x 1 u 1 We define +e i x 1 u = 1 λ+ r1 1 +ir 1 λ+ ir1 1 r 1 λ+ τr1 1 +ir 1 +O1 λ+ τir1 1 r 1 +O1 i λ+ r1 1 +ir 1 i λ+ r1 1 +ir 1 λ+ τ r 1 ir +O1 λ+ τ ir 1 +r +O1 i λ+ ir 1 r i λ+ ir 1 r +Oτ 1, +Oτ 1 +Oτ 1. +Oτ 1. O1 O1 i λ+ 1+ir ir 1 i λ+ 1+ir ir 1 u 1 = Re i x u ; u = Ie i x u ; =,, u 1 { Then u 1,u,u 3 is sufficiently large. u 3 = Re i x u 1 ; u = Ie i x u 1 ; = R e i x 1 u 1 +e i x 1 u. +Oτ 1. } are linearly independent in a small neighborhood U of x when τ

7 7 Therefore, for l = 1,, there exist functions Θ l j, j = 1,,3 such that 16 u l + Θ l j uj =, l = 1,. Since u, u 1, u are solutions of the equation 6, we have the following equations: u j u l F +u l G = k u l, l = 1,; F +u j G = k u j, j = 1,,3. Summing over the equations above using 16, we have two equations 17 β,l +γ,l = k u l + Θ l j uj, for l = 1,, where β,l = v l e 1 ẽ 1 +v l e ẽ γ,l = v l e 3 +v l e 4, Here we define v l = u l + 3 Θl j uj. Remark. By choosing suitable g, Θ l 3 x can be as small as we want. To show that, we choose a new constant vector ĝ, instead of the original g, such that ˆr,ŝ T = C x,θĝ θ x where ˆr x = r x /M, M >, ŝ x = s x ˆr,ŝ T satisfies the original assumption, that is, ŝ x = 1 ŝ ˆr in a neighborhood of x. Then Θ l 3x = P 3j 1 R ˆΘ l 3x = P 3j 1 R M λ+ 1,i r x λ+ i1,i r x λ+ 1,i r x λ+ i1,i r x +Oτ 1 p j +Oτ 1 p j, where P = p ij P 1 = P ij with p ij = p i p j. Note that P is a boundedly invertible symmetric matrix. Here p 1 = u 1, p = u, p 3 = u 3. From above, we obtain that ˆΘ l 3 is close to Θ l 3/M as τ is large. Therefore, the new ˆΘ l 3x is small when M τ is sufficiently large. Lemma 3.. Given any point x Ω, there exists a small neighborhood V of x such that β,j is not zero in V for j = 1,. Note that we denote by equalities up to terms which are asymptotically negligible as τ goes to infinity.

8 8 Proof. It is sufficiently to provethat β,1 doesnot vanish in someneighborhoodofx. Recall that 18 where β,1 = v 1 e 1 ẽ 1 +v 1 e ẽ, v 1 = u 1 + Θ 1 j uj. By 16, we have that Θ 1 j τ for j = 1,,3. Since s 1 x =, we get that 19 u 1 e 1 ẽ 1 + u 1 e ẽ x τ 1 u e 1 ẽ 1 + u e ẽ x τ 1. Since s k x =, r = r 1, s = s 1, 1, i r 1 x =, we obtain that u 3 e 1 ẽ 1 + u 3 e ẽ x 1 = 4 1 R 1 i 1 s 1 + s 1 1 s 1 +i s 1 +Oτ 1 x. Combining from 18 to 1, since Θ 1 3 x can be taken as small as we want See Remark in the construction of CGO solutions above, it follows that 4 β,1 x 1 +Θ 1 i 1 s 1 + s 1 3R 4 1 s 1 +i s 1 x. Similarly, β, x 1 4 +Θ i 1 s 1 + s 1 3R 4 1 s 1 +i s 1 x. By continuity of β,j, we complete the proof. Let g j = u j Ω be the given boundary data. By Lemma.1, since λ, P, we knew that u j H 4 B R. Let g j C 1,α Ω be the boundary data close to g j, that is, g j g j C 1,α Ω < ε for some ε >, then we can find solutions u j of the elasticity system with boundary data g j the existence of such solutions, see e.g. Ch.4 of [1]. By elliptic regularity theorem, we have u j u j C Ω < Cε for some constant C which is independent of λ,. Then we obtain that u j u j C Ω < Cε. Herethe notationu j is constructedin the samewayasu j with the CGOsolutionsu j replaced by the solutions u j. Therefore, { u 1,u,u 3 } are also linearly independent when ε is sufficiently small. We construct β j as in the equation 17 with u j replaced by u j,j =,1,. Therefore, by the definitions of β j β,j, it follows that β j β,j C 1 Ω is small when ε is small. Since β,j is not zero in V by Lemma 3., we can deduce that β j is also not zero in V if ε is small enough τ is sufficiently large. Moreover, with the suitable chosen CGO solutions u, {β,1,β, } are linearly independent in V as τ is sufficiently large

9 9 See the proof of Lemma 3.. When ε is sufficiently small, it implies that {β 1,β } are also linearly independent in V. Then we have the following equations: 3 β l +γ l = k u l + Θ l j uj, l = 1, with {β 1,β } a basis in R for every point x Ω. Here we denote Ω = U V Ω. Thus, there exists an invertible matrix A = a ij such that β l = a lk ẽ k with inverse of class C Ω. Thus, we have constructed two vector-valued functions Γx, Φx CΩ such that 3 can be rewritten as 4 Then we obtain the following result: +Γx = Φx in Ω. Theorem 3.3. Suppose that λ, λ, P. For any fixed x Ω, let u j be the corresponding CGO solutions for λ, u j constructed as above with internal data H x = u j j with ε sufficiently small. Let H x be constructed similarly with the parameters λ,. Assume that Ω = Ω. Then H x = H x implies that = in Ω, the neighborhood of x in Ω. Proof. Based on the construction above, the domain Ω can be taken as a small open ball with center x Ω B R. Since H x = H x, we have that solve the same equation +Γx = Φx in Ω where the functions Γ Φ depend on u j. Let x Ω denote ψt = 1 tx +tx, t [,1]. Restricted to this curve, we have { ψ t +Γψt ψ t = Φψt ψ t in Ω 5 x = x, The solution of 5 is given by ψt = ψe t Γψs ψ sds + t Φψs ψ sds, t [,1]. The solution x is given by the same formula since Ω = Ω so that = in Ω. We have constructed { u 1,u,u 3 } are linearly independent {β 1,β } forms a basis in R for every point x Ω when ε is sufficiently small τ is large. Applying Lemma 3.1, we have Theorem 3.4. Suppose that λ, λ, P. For any fixed x Ω, let u j be the corresponding CGO solutions for λ, u j constructed as above with internal data H x = u j j with ε sufficiently small. Let H x = ũ j j be constructed similarly for the parameters λ, with u j Ω = ũ j Ω. Then there exists an open neighborhood Ω of x in Ω such that CΩ C x x + H x H 6 x C Ω, x Ω. 7 CΩ C x + x+ + H x H x C Ω, x + Ω, x Ω. The global uniqueness stability result are stated as follows. Theorem 3.5. Global reconstruction of Let Ω be an open bounded domain of R with smooth boundary. Suppose that the Lamé parameters λ, λ, P Ω = Ω. Let u j ũ j be the solutions of the elasticity system with boundary data g j for parameters λ, λ,, respectively. Let H = u j 1 j J H = ũ j 1 j J be the corresponding internal data for λ, λ,, respectively for some integer J 3.

10 1 Then there is an open set of the boundary data g j 1 j J such that if H = H implies = in Ω. Moreover, we have the stability estimate CΩ C H H C Ω. Note that for the uniqueness of, we suppose that the two set of internal data are the same, that is, H = H. Since is uniquely reconstructed near a fixed point x Ω under the condition Ω = Ω, from the stability of in Ω, we can obtain that = in Ω. Proof. In section, we constructed CGO solutions in a ball B R which contains Ω. First, we consider any point x in Ω. Then we can find an open neighborhood B x B R of x. By Theorem 3.4, we have the estimate 8 CBx Ω C H x H x C Ω since Ω = Ω. Second, for any point y Ω, by Theorem 3.4, there exists an open neighborhood B y Ω of y with B y Ω = such that CBy C y + y + + H y H 9 y C Ω for some y + B y. Therefore, the compact set Ω is covered by x Ω B x. Then there exists finitely many B x, say, B x1,...,b xn, such that Ω N l=1 B x l. Nowforarbitrarypointz Ω, thereisb xj suchthatz B xj. SupposethatB xj Ω, this means that x j Ω. Then, by 8, we have 3 z z C H xj H xj C Ω. Otherwise, if B xj Ω is empty, then, by 9, we get that z z C x + j x+ j + H x j H 31 xj C Ω for x + j B x j. For the point x + j, since Ω is covered by finitely many subdomain B x l, after at most N 1 steps, we have 3 x + j x+ j C N Combining Then we have 33 z z C l j,l=1 H xl H xl C Ω. N H xl H xl C Ω. With 3 33, we have the global stability N 34 CΩ C H xl H xl C Ω. l=1 l=1 3.. Reconstruction of λ in D. The elasticity system can also be written in this form 35 where u = u u a+b a b, F = u F +u G = k u, 1 λ+ λ+ 1, u = 1 u u 1 a+b a b, G = λ+ λ+.

11 11 As in the reconstruction of, we will construct there linear independent vectors such that the first term of the equation 35 can be eliminated. Suppose that u j, for j =,1,,3, are three different solutions of 6 which satisfy 36 u j F +u j G = k u j. Now for j = 1,,3 =,,, we let u = Rχxu, u 1 = Rχ 3 xu 1, u = Iχ 3 xu 1, u 3 = R χ 1 xu +χ xu 3, where χ j x is a nonzero function. Then we get 37 u F +u G = k u 38 u j F +u j G = k u j for j = 1,,3. Assume that { u 1,u,u 3 } are three linearly independent vectors in some subdomain of Ω, say Ω. Then there exist three functions Θ 1,Θ, Θ 3 such that u + 3 Θ ju j =. Multiplying 38 by Θ j summing over j with equation 37, we have v G = k u + Θ j u j, where v = u + 3 Θ ju j. Let Then the above equation can be rewritten as κ = 1,1,, T v, σ = 1,1,1,1 T v. κλ = σ k u + Θ j u j. Suppose that κx does not vanish in Ω. Then we can deduce the following lemma. Lemma 3.6. Let u j for j 3 be C solutions of the elasticity system with boundary conditions u j = g j on Ω. Let u = u 1,u, a = u u b = 1 u 1 u. We define u = u, u, a+b, a b assume that 1 { u 1,u,u 3 } are three linearly independent vectors in Ω, the neighborhood of x in Ω. κx does not vanish in Ω. Then the reconstruction is stable in the sense that λ λ CΩ C CΩ + H x H 39 x C Ω where H x = u j j 3.

12 Global reconstruction of λ. Wewillfindthreelinearlyindependentvectors { u 1,u,u 3 } first. Then we can deduce stability of λ by using Lemma 3.6 Theorem 3.5. Plugging the CGO solution u = 1/ w + 1 f f 1 into u, we get λ+ r u = ie i x λ+ r 1 is ω 1 s+ 1/ ω 1 r +O1, 1 is 1 1 +ω s+ 1/ ω r where ω 1 = 1 +, 1 ω = 1, 1 +. For the same fixed point x Ω. Denote that = τe 1 +ie = i. We choose a constant vector g such that C x,θg = r,s with s, r, r x = 1 in a neighborhood of x, U. Then we get the CGO solution of the elasticity system with u = 1/ w + 1 f f 1 w f = e i x r s +Oτ 1. We choose another constant vector g 1 such that C 1 x,θg 1 = r 1,s 1 with s 1 in U r 1 x =. Then we get the CGO solution of the elasticity system with u 1 = 1/ w f 1 f 1 1 w 1 f 1 = e i x r 1 s 1 +Oτ 1. For l =,3, we choose a constant vector g l such that C l x,θg l = r l,s l with r l in U. Here we can choose r = r 3, s = s 3 by takingg = g 3 C x,θ = C 3 x, θ. Then we getthe CGOsolutionofthe elasticity system u l = 1/ w l + 1 f l f l 1 with w f = e i x r s +Oτ 1 w 3, f 3 = e i x r 3 s 3 +Oτ 1. We let e i x i 1 u 1 u,1 := e i x i 1 1 u = 1 = i 1 s 1+i i 1 +Oτ 1, λ+ 1,i r λ+ 1,i r 1 τ1+iis +ν 1 τi 1is +iν +Oτ 1

13 13 u, := e i x i 1 1 u 3 = 1 λ+ i, 1 r3 λ+ i, 1 r3 1 τ 1 iis 3 +ν 3 1 τ1 iis 3 +iν 3 +Oτ 1, where ν = 1 1+i,1 i s + 1/ 1+i,1 i r ν 3 = 1 1+i,i+ 1 s 3 + 1/ 1+i,i+1 r 3. To eliminate the higher order term of u,j, we consider the summation of two vectors : We define e i x i 1 1 u +e i x i 1 1 u 3 = 1 u 1 λ+ 1+i,i 1 r3 λ+ 1+i,i 1 r3 1 i, s + 1/ i, r 3 1,i s 3 + 1/,i r 3 u = Re i x i 1 u ; =,, +Oτ 1. = Re i x i 1 u 1 ; u = Ie i x i 1 u 1 ; u 3 = R e i x i 1 1 u +e i x i 1 1 u 3. { } Thus we have constructed three linear independent vectors u 1,u,u 3 in U as τ is sufficiently large. Therefore, there are three functions Θ j, j = 1,,3, such that 4 u + They also satisfy the following equations: u j F +u j Θ j u j =. G = k u j, j 3. Summing over j using 4, we get the following equation 41 v G = k u + Θ j u j, where 4 Then we obtain that where v = u + κλ = σ k u + Θ j u j. κ = 1,1,, T v, σ = 1,1,1,1 T v. Θ j u j Lemma 3.7. κ does not vanish in some neighborhood of x.,

14 14 Proof. Since Θ j τ r 1 x =, by using similar argument as in the proof of Lemma 3., we have that 43 Observe that 1,1,, T Θ 1 u 1 κ 4 := 1,1,, T u 3 1 e i x 1 u =R 1 i u The Oτ term of κ 4 is 1 i τ λ+ 1 + r τ 1, 1,1,, T Θ u τ 1. + e i x i 1 u 3 u 3. λ r 3 = since r = r 3. Thus the leading order term of κ 4 is O1, that is, where Since ϕ = ϕ +ϕ 3, ϕ l = 1 i 1/ i l τ 1r l 1 +il 1 r l +il 1 r l 1 +il 1r l +i r l + 1 i 1 s l + s l +i 1 + s l +i 1 1 l s l 11 1 i l 1 +l sl +i 1 l s l 1 i l 1 +l sl + 1 1/ i l rl 1 +il rl +il 1 rl 1 + 1/ i l 1 rl 1 +il 1 rl +il rl. κ 1 :=1,1,, T u x 1 ie i x 1 u = R 1 u 43, we get that κx =1,1,, T vx =κ 1 x +Θ 3 κ 4 x +Oτ 1 λ+ +Θ 3x ϕx. x = λ+ +Oτ 1 Since P, ϕx is some fixed number. Moreover, we can take s small such that Θ 3 x to be sufficiently small It can be done by following similar argument in Remark, then we can obtain that λ+ +Θ 3x ϕx. By continuity of κ, there exists a neighborhood V such that κ never vanishes in V. Let Ω = U V Ω. We have { u 1,u,u 3 } are linearly independent in Ω as τ is sufficiently large also κ does not vanish in Ω. Then it follows that 47 λ = σ κ k u + Θ j u j κ in Ω.

15 15 48 As in section 3.1, we also can find the solutions u j of the elasticity system such that u j u j C Ω < Cε, j =,1,,3. Now we let the internal data H x contains the three solutions we constructed in Theorem 3.3 the four solutions u j,j =,1,,3 in this section. Then we have the following result. Theorem 3.8. Suppose that λ, λ, P. For any fixed x Ω, let u j be the corresponding CGO solutions for λ, u j constructed in sections 3.1, 3. with internal data H x with ε sufficiently small. Let Hx be constructed similarly with the parameters λ,. Assume that Ω = Ω. Then H x = H x implies that λ = λ in Ω. Proof. Applying Theorem 3.3 equation 47, we have the uniqueness of λ near the point x. We deduce the following result by applying Lemma 3.6 Theorem 3.4. Theorem 3.9. Suppose that λ, λ, P. For any fixed x Ω, let u j be the corresponding CGO solutions for λ, u j constructed in sections 3.1, 3. with internal data H x with ε sufficiently small. Let H x be constructed similarly for the parameters λ, with u j Ω = ũ j Ω. Then we have the estimates λ λ CΩ C x x + H H 49 C Ω, x Ω. 5 λ λ CΩ C x + x+ + H H C Ω, x + Ω, x Ω. With 47 Theorem 3.9, we follow the same proof as in Theorem 3.5, then we can get the following global reconstruction of λ. Theorem 3.1. Global reconstruction of λ Let Ω be an open bounded domain of R with smooth boundary. Suppose that the Lamé parameters λ, λ, P Ω = Ω. Let u j ũ j be the solutions of the elasticity system with boundary data g j for parameters λ, λ,, respectively. Let H = u j 1 j J H = ũ j 1 j J be the corresponding internal data for λ, λ,, respectively for some integer J 7. Then there is an open set of the boundary data g j 1 j J such that if H = H implies λ = λ in Ω. Moreover, we have the stability estimate λ λ CΩ C H H C Ω. 4. Reconstruction of Lamé parameter in three-dimensional case The reconstruction of λ in R 3 basically follows the similar argument as in Section 3. In R 3, we need more CGO solutions to get linearly independent vectors locally Global reconstruction of in 3D. Let u = u 1,u,u 3 T be the displacement which satisfies the elasticity system 51 λ ui +S u+k u =.

16 16 Denote u = 1 u u 3 u u u u, F = λ+ λ+ λ+ 1 λ+ λ+ 3 λ+, u = b 3 b 13 b 1 1 b 3 b 13 3 b 1, G = where b ij = l u l i u i j u j + i u l + l u i + j u l + l u j with l,i,j {1,,3} are distinct numbers. From 51, we can deduce the following equation: u F +u G = k u. Here u = u 1 +u +u 3. In the following we will show that how we can get four linearly independent vectors of the form u on some subdomain of Ω. The key thing is to observe the behavior of u. We plug the CGO solutions u = 1/ w + 1 f f 1 into u. Then we get 5 u = e i x λ+ 1 r+o λ+ r+o λ+ 3 r+o i λ+ r i λ+ r i λ+ r +O1. Note that r = r 1,r,r 3 T = 1,, 3 T. Now we fix any point x Ω. Let = τ1,i, = τ1,,i with τ >. Let θ = /τ, θ = /τ. Since, in Lemma.1, the matrix solutions C x,θ C x, θ are invertible, we can choose two constant vectors g g such that C x,θg = r,s T C x, θ g = r, s T withs x = 1 = s x s, s r, r in a neighborhood of x, say U. Then we have the CGO solutions of the elasticity system, that is, u = 1/ w + 1 f f 1, ũ = 1/ w 1 + f f 1 with w f w f = e i x r s = e i x r s +Oτ 1, +Oτ 1. Let = τi, 1, θ = /τ. Let C 1 x,θ C x, θ satisfy that iθ x C 1x,θ = θ V 1 xc 1 x,θ, i θ x C x, θ = θ V 1 xc x, θ, respectively. Since = i, we can choose C x, θ = C 1 x,θ. Moreover, r = r 1 s = s 1. With suitable constant vector g, we can get that s l is zero at point x r 1 x = 1, i, = r x. By continuity of r l, we have r l in a neighborhood U 1 of x. Then the CGO solution is u l = 1/ w l + 1 f l f l 1 1 3,

17 17 with w 1 f 1 w f = e i x r 1 s 1 = e i x r s +Oτ 1, +Oτ 1. For, with a suitable constant vector g, we can get that s 3 is zero at point x r 3 x = 1,, i. By continuity of r 3, we have r 3 in a neighborhood U of x. Then the CGO solution is u 3 = 1/ w f 3 f 3 1 with w 3 f 3 = e i x r 3 s 3 +Oτ 1. Let U = l= U l. So far we have five CGO solutions, that is, u,ũ,u 1,u, u 3. Let r 1 = r 1 1,r1. We define,r1 3 u 1 := e i x u 1 = 1 u,1 := e i x 1 u 1 = 1 u, := e i x 1 u = 1 λ+ r1 1 +ir 1 λ+ ir1 1 r 1 +Oτ 1 λ+ τr1 1 +ir 1 +O1 λ+ τir1 1 r 1 +O1 i λ+ r1 1 +ir 1 +Oτ 1 i λ+ r1 1 +ir 1 i λ+ r1 1 +ir 1 λ+ τ r 1 ir +O1 λ+ τ ir 1 +r +O1 i λ+ ir 1 r +Oτ 1. i λ+ ir 1 r i λ+ ir 1 r Let u = u,1 +u,. Since r = r 1 s = s 1, we have u = 1 O1 O1 i λ+ 1+ir ir 1 i λ+ 1+ir ir 1 i λ+ 1+ir ir 1 +Oτ 1.

18 18 We also define We denote Then u 3 u 3 := e i x u 3 = 1 = R { u j : 1 j 4 λ+ r3 1 +ir 3 3 λ+ ir3 1 r 3 3 +Oτ 1 u 1 = Re i x u ; u = Ie i x u ; u 1 e i x 1 u 1 } u 3 = Re i x ũ ; =,, = Re i x u 1 ; u +e i x 1 u ; u 4 = Ie i x u 1 ; = Re i x u 3. are linearly independent in the neighborhood U of x as τ is sufficiently large. Therefore, for fixed l = 1,,3 there exist functions Θ l j, j = 1,,3,4, such that u l + 4 Θ l ju j =. As in Section 3.1, we summing over equations, then we have 4 β,l +γ,l = k u l + Θ l j uj for l = 1,,3, where β,l γ,l are functions which depend on,ω CGO solutions ũ,u j j =,...,3. Lemma 4.1. Given any point x Ω, there exists an open neighborhood V of x such that β,l is not zero in V, for l = 1,,3. Proof. Following the similar proof as in Lemma 3., we can prove this Lemma. Let Ω = U V Ω. Based on the lemma above, we may suppose that β,j are linearly independent in Ω as τ sufficiently large. Let g j = u j Ω g = ũ Ω be the given boundary data for j =,1,,3. Let g, g j C 1,α Ω be the boundary data close to g, g j, respectively, that is, g j g j C 1,α Ω < ε, g g C 1,α Ω < ε. Then we can find solutions ũ, u j of the elasticity system with boundary data g,g j, respectively. By regularity theorem, it follows that 53 u j u j C Ω < Cε, ũ ũ C Ω < Cε. From 53, we have u j u j CΩ Cε, 1 j 4. Therefore, { u j : 1 j 4 } are also linearly independent when ε is sufficiently small. We construct β l by replacing u j by u j. Then from 53, we can deduce that β l β,l C1 Ω for

19 19 is small when ε is sufficiently small. Then we have the following equation 4 54 β l +γ l = k u l + Θ l ju j, l = 1,,3, with {β l x} l=1,,3 a basis in R 3 for every point x Ω. There exists an invertible matrix A = a ij such that β l = a lk e k with inverse of class CΩ. Thus, we have constructed two vector-valued functions Γx, Φx CΩ such that 54 can be rewritten as 55 +Γx = Φx in Ω. Then we have the following uniqueness stability theorem. Theorem 4.. Global reconstruction of Let Ω be an open bounded domain of R 3 with smooth boundary. Suppose that the Lamé parameters λ, λ, P Ω = Ω. Let u j ũ j be the solutions of the elasticity system with boundary data g j for parameters λ, λ,, respectively. Let H = u j 1 j J H = ũ j 1 j J be the corresponding internal data for λ, λ,, respectively for some integer J 5. Then there is an open set of the boundary data g j 1 j J such that if H = H implies = in Ω. Moreover, we have the stability estimate CΩ C H H C Ω. Proof. The proof is similar to Theorem Global reconstruction of λ in 3D. The elasticity system can also be written in this form 56 where u = u u u b 3 b 13 b 1, F = u F +u G = k u, 1 λ+ λ+ 3 λ+ 1 3, u = 1 u u 3 u 1 b 3 b 13 3 b 1, G = λ+ λ+ λ+ As in the reconstruction of, we will construct CGO solutions such that the first term of 56 can be eliminated. Plug the CGO solution u = 1/ w + 1 f f 1 into u. Then we get u = e i x i λ+ r i λ+ r i λ+ r s+o s+o s+o +O1. Forthesamefixedpointx Ω. Wechooseaconstantvectorg suchthatc x,θg = r,s with s r in a neighborhood of x, say U, r x = 1. Then we get the CGO solution of the elasticity system u = 1/ w + 1 f f 1.

20 with w f = e i x r s +Oτ 1. We choose another constant vector g 1 such that C 1 x,θg 1 = r 1,s 1 with s 1 in a neighborhood of x, say U 1, r 1 x =. Then we get the CGO solution of the elasticity system u 1 = 1/ w f 1 f 1 1 with w 1 f 1 = e i x r 1 s 1 +Oτ 1. For l =,3, we choose a constant vector g l such that C l x,θg l = r l,s l with r l in a neighborhood of x, say U. Here we can choose r = r 3, s = s 3 by taking g = g 3 C x,θ = C 3 x, θ. Then we get the CGO solution of the elasticity system u l = 1/ w l + 1 f l f l 1 with w f w 3 f 3 = e i x r s = e i x r 3 s 3 +Oτ 1, +Oτ 1. For, we choose another constant vector g 4 such that C 4 x,θg 4 = r 4,s 4 with s 4 in a neighborhood of x, say U 3, r 4 x =. Then we get the CGO solution of the elasticity system u 4 = 1/ w f 4 f 4 1 with We define w 4 f 4 u 1 := e i x u 1 u,1 :=e i x 1 u = 1 = e i x r 4 s 4 = 1 s +Oτ 1. i λ+ r i λ+ r 1+i i 1 i λ+ r 1 τ1+is +O1 1 τi 1s +O1 O1 +Oτ 1 +Oτ 1

21 1 u, :=e i x 1 u 3 = 1 i λ+ r3 i λ+ r3 i λ+ r3 1 τ 1 is 3 +O1 1 τ1 is 3 +O1 O1 +Oτ 1. Let u = u,1 +u,, then the higher order is eliminated. Thus we have Also, we define that We denote u 3 u = 1 i λ+ + r3 + r3 i λ+ i λ+ u 3 := e i x u 4 u + r3 O1 O1 O1 = 1 s +Oτ 1. 1+i i 1 = Re i x u ; =,, u 1 = R e i x 1 u } = Re i x u 1 ; u +e i x 1 u 3 ; u 4 +Oτ 1. = Ie i x u 1 ; = Re i x u 4. { Then u j : 1 j 4 are linearly independent in the neighborhood U = 3 l= U l of x as τ is sufficiently large. Therefore, there exist functions Θ j, j = 1,,3,4, such that u + 4 Θ j u j =. Summing over, we get the following equation as in Section 3.: 4 57 v G = k u + Θ j u j, with v = u + 4 Θ ju j. We obtain that 58 κλ = σ k u + 4 Θ j u j,

22 where κ = 1,1,1,,, T v, σ = 1,1,1,1,1,1 T v. Lemma 4.3. κ does not vanish in some neighborhood of x. Proof. Similar argument as Lemma 3.7. Since 1 u λ+ u = e i x 1 r+o λ+ r+o 3 u λ+ r+o r 1 x = = r 4 x, we have 1,1,1,,, T Hence, we obtain that κx = 1,1,1,,, T Θ 1 u 1 +Θ u u λ+ +Θ 3u 3 x. +O1 +Θ 4 u 4 x τ 1. +Θ 3 u 3 x +Oτ 1 We can take s small enough such that Θ 3 x is small. Thus, κx. By continuity of κ, κ does not vanish in some neighborhood V of x. Let Ω = U V Ω. Since κ does not vanish in Ω, we have λ = σ 4 59 κ k u + Θ j u j in Ω. κ Applying the similar proof as in Theorem 3.5, we can deduce the following result. Theorem 4.4. Global reconstruction of λ Let Ω be an open bounded domain of R 3 with smooth boundary. Suppose that the Lamé parameters λ, λ, P Ω = Ω. Let u j ũ j be the solutions of the elasticity system with boundary data g j for parameters λ, λ,, respectively. Let H = u j 1 j J H = ũ j 1 j J be the corresponding internal data for λ, λ,, respectively for some integer J 1. Then there is an open set of the boundary data g j 1 j J such that if H = H implies λ = λ in Ω. Moreover, we have the stability estimate Acknowlegments λ λ CΩ C H H C Ω. The author are grateful to professor Gunther Uhlmann for his encouragements helpful discussions. The author also would like to thank Professor Jenn-Nan Wang for taking the time to discuss some properties of the elasticity system with her. The author is partially supported by NSF. References [1] G. Alessrini, A. Morassi E. Rosset, Detecting an inclusion in an elastic body by boundary measurements, SIAM J. Math. Anal., 336, [] G. Bal, Hybrid inverse problems internal functionals, Inside Out II, MSRI Publications, 6 1. [3] G. Bal G. Uhlmann, Reconstruction of coefficients in scalar second-order elliptic equations from knowledge of their solutions, To appear in Comm. on Pure Applied Math, 13. [4] G. Bal G. Uhlmann, Inverse diffusion theory for photoacoustics, Inverse Problems, 6 1.

23 [5] P. G. Ciarlet, Mathematical Elasticity, Studies in Math. its Appl., [6] G. Eskin, Global uniqueness in the inverse scatteing problem for the Schrödinger operator with external Yang-Mills potentials, Comm. Math. Phys., 1, [7] G. Eskin J. Ralston, On the inverse boundary value problem for linear isotropic elasticity, Inverse Problems, 18, [8] H. Heck, X. Li, J.-N. Wang, Identification of viscosity in an incompressible fluid, Indiana University Matheatics Journal, 56 6, [9] M. Ikehata, A remark on an inverse boundary value problem arising in elasticity, preprint. [1] W. McLean, Strongly Elliptic Systems Boundary Integral Equations, Cambridge University Press,. [11] J. McLaughlin D. Renzi, Shear wave speed recovery in transient elastography supersonic imaging using propagating fronts, Inverse Problems, 6, [1] J. R. McLaughlin, N. Zhang, A. Muca, Calculating tissue shear modulus pressure by D log-elastographic methods, Inverse Problems, 6 1, 857, 5. [13] G. Nakamura G. Uhlmann, Global uniequness for an inverse boundary problem arising in elasticity, Invent. Math., , [14] G. Nakamura G. Uhlmann, Erratum: Global uniqueness for an inverse boundary problem arising in elasticity, Invent. Math., 15 3, 5 7. [15] G. Nakamura G. Uhlmann, Complex geometrical optics solutions pseudoanalytic matrics, In: S.I. Kabanikhin, V.G. Romanov eds, Ill-Posed Inverse Problems,, [16] G. Nakamura, G. Uhlmann, J.-N. Wang, Oscillating-decaying solutions, Runge approximation property for the anisotropic elasticity system their applications to inverse problems, J. Math. Pures Appl., 84 5, [17] G. Uhlmann J.-N. Wang, Complex geometrical optics solutions reconstruction of discontinuities, SIAM J. Appl. Math., 68 8,