c 2005 Society for Industrial and Applied Mathematics
|
|
- Primrose May
- 5 years ago
- Views:
Transcription
1 MULTISCALE MODEL. SIMUL. Vol. 4, No. 4, pp c 5 Society for Industrial and Applied Mathematics DIRECT ALGORITHMS FOR THERMAL IMAGING OF SMALL INCLUSIONS HABIB AMMARI, EKATERINA IAKOVLEVA, HYEONBAE KANG, AND KYOUNGSUN KIM Abstract. The goal of this paper is to reconstruct a collection of small inclusions inside a homogeneous object by applying a heat flux and measuring the induced temperature on its boundary. Taking advantage of the smallness of the inclusions, we design efficient noniterative algorithms for locating the inclusions from boundary measurements of the temperature. We illustrate the feasibility and the viability of our algorithms by numerical examples. Key words. thermal imaging, small inclusions, direct inversion algorithms AMS subject classifications. 35R, 35B3 DOI..37/4666. Introduction. Thermal imaging is a technique of wide utility in nondestructive testing and evaluation. To apply this procedure, one uses a heat source to apply a heat flux to the boundary of an object and then observes the resulting temperature response on the object s surface. From this information one attempts to determine internal thermal conductivity of the object. Thermal imaging has been significantly investigated as a method for damage and corrosion detection [, 8, 7]. For mathematical considerations of thermal imaging we refer to [7, 4, 5, 6,, ] and the references therein. In its most general form thermal imaging is severely ill-posed and nonlinear [7]. If, however, in advance we have additional structural information about the object, then we may be able to determine specific features about the internal thermal conductivity of the object with a good resolution. One such type of knowledge could be that the object consists of a smooth background containing a number of unknown small inclusions with significantly different thermal conductivities. In this case thermal imaging seeks to recover the unknown inclusions. Due to the smallness of the inclusions the associated temperatures measured on the surface of the object are very close to the temperatures corresponding to the object without inclusions. So unless one knows exactly what patterns to look for, noise will largely dominate the information contained in the measured temperature. Furthermore, in applications it is often not necessary to reconstruct precisely the thermal properties or the geometry of the inclusions. The information of real interest is their positions. Our goal in this paper is to develop, in a mathematically rigorous way, a numerical method to detect the location of the inclusions from boundary measurements of the temperature. We find an asymptotic expansion of solutions to the transmission problem for the heat equation in terms of the geometry of the inclusion. This formula Received by the editors December 6, 4; accepted for publication (in revised form June 6, 5; published electronically October 7, 5. This research was partly supported by ACI NIM 7, CNRS-KOSEF grant 4889, Korea Science and Engineering Foundation grant R-3-- -, and BK at Seoul National University. Centre de Mathématiques Appliquées, CNRS UMR 764 and Ecole Polytechnique, 98 Palaiseau Cedex, France (ammari@cmapx.polytechnique.fr, iakov@cmapx.polytechnique.fr. School of Mathematical Sciences, Seoul National University, Seoul 5-747, Republic of Korea (hkang@math.snu.ac.kr, kgsun@math.snu.ac.kr. 6
2 DIRECT ALGORITHMS FOR THERMAL IMAGING 7 describes the perturbation of the solution caused by the presence of an anomaly of small diameter. Based on this asymptotic expansion we derive formulae to find the location of the inclusions with good accuracy. The formulae are explicit and can be easily implemented numerically. The general approach we will take is in the spirit of that developed for reconstructing small conductivity inclusions from boundary measurements; see [6, 9, 4,,, 8, ] and the references therein. The paper is organized as follows. In the next section we formulate the problem and state the asymptotic expansion of the boundary perturbations that are due to a collection of conductivity inclusions. The rigorous derivation of this formula is explained in sections 3 and 4. In section 5 we use our asymptotic expansion to design direct reconstruction algorithms for locating the inclusions. Numerical examples given in this section illustrate the feasibility and the viability of our direct procedures.. Problem formulation. Let Ω be a bounded Lipschitz domain in R d, d =, 3. Assume that Ω contains a finite number m of small inclusions {D s } m, each of the form D s = εb s + z s, where B s, s =,...,m, is a bounded Lipschitz domain in R d containing the origin. We denote the collection of inclusions by D = m D s. We assume that the inclusions D s, s =,...,m, are well separated from each other and from Ω. More precisely, we assume that there exists a constant d > such that (. inf dist(x, Ω >d and dist(d i,d j >d, i j. x D s We also assume that the background is homogeneous with thermal conductivity γ and that the inclusion D s has thermal conductivity γ s, where γ s and γ are positive constants, γ s γ, s =,...,m. The thermal conductivity profile γ ε ofωisthe piecewise constant function { γ for x Ω \ D, γ ε = γ s for x D s,,...,m. In this paper we consider the following transmission problem for the heat equation: t u (γ ε u = in Ω (,T, (. u(x, = θ(x for x Ω, u γ = g on Ω (,T, where the Neumann boundary data g L (,T; L ( Ω and the initial data θ(x C (Ω are subject to the compatibility condition (.3 γ θ = g(, on Ω. Here and throughout this paper, the vector ν denotes the unit outward normal to Ω. We will also use the following notation: Ω T =Ω (,T, D st = D s (,T, Ω T = Ω (,T, D st = D s (,T. Let U be the background solution defined as the solution of t U γ ΔU = in Ω T, U(x, = θ(x for x Ω, U γ = g on Ω T.
3 8 H. AMMARI, E. IAKOVLEVA, H. KANG, AND K. KIM Introduce ψ (s j for j =,...,d and s =,...,m to be the unique solution of Δψ (s j (x = in R d \ B s, ψ (s j + (x ψ (s j (x =x j on B s, (.4 ψ (s j ψ (s γ j (x γ s x j + (x =γ on B s, ψ (s j (x =O( x d as x. Here ν denotes the outward unit normal to B s ; superscripts + and indicate the limiting values as we approach B s from outside B s and from inside B s, respectively. Let (.5 Ψ (s j (x := { ψ (s j (x, x R d \ B s, ψ (s j (x+x j, x B s. Then it follows from (.4 and (.5 that Ψ (s j + Ψ (s Ψ (s j Ψ (s j =, γ j γ s + =(γ γ s x j (.6 Define the d d matrix M (s =(m (s ij for s =,...,m by m (s ij = (s j ψ y i B s (y dσ y, on B s. where y i is the ith component of y. The matrix M (s is the Pólya Szegö polarization tensor associated with the domain B s and the conductivity γ s [3]. Its properties were extensively studied in []. Recall that if B s is a disk, then the polarization tensor M (s has the following explicit form [9]: M (s dγ B s (.7 = (γ s +(d γ I d, where I d is the d d identity matrix. The following asymptotic expansion holds. Theorem.. Let Φ C (Ω T satisfy ( t + γ ΔΦ(x, t = in Ω T with Φ(x, T =for x Ω. Define the weighted boundary measurements I Φ (T := γ (u U(x, t Φ Ω (x, t dσ x dt. Then the following asymptotic expansion for I Φ (T holds as ε : (.8 I Φ (T = ε d m (γ γ s { O(ε 4 U(z s,t M (s if d =3, Φ(z s,t dt + O(ε 3 log ε if d =. By making appropriate choices of test functions Φ and background solutions U we will develop from the asymptotic formula (.8 efficient location search algorithms for detecting the inclusions D s, s =,...,m, from measurements of u U on Ω T.
4 DIRECT ALGORITHMS FOR THERMAL IMAGING 9 3. Preliminary results. For the transmission problem for the heat equation there are convenient functional spaces, which we recall below. The first is formed of all functions v(x, t onω T with finite norm ( / v L (,T ;H(Ω = v(,t H(Ω dt, where H(Ω = L (Ω or H (Ω, and the second is formed of all functions v(x, t onω T with finite norm where v ΩT = sup v(,t,ω + v,ωt, t T ( / ( v(,t,ω := v(x, t dx, v =,ΩT Ω Ω v(x, t dx dt /. We also recall that the Lipschitz character of Ω is the L -norm of the gradient of the function parameterizing the boundary Ω. See, for example, []. It is clear that the Lipschitz character of the domain tω, t>, is independent of t. The following regularity result on the transmission problem for the heat equation will be of use to us. Lemma 3.. Let v be the solution of the following transmission problem for the heat equation: t v(x, t γ Δv(x, t =F (x, t in Ω \ D (,T, t v(x, t γ s Δv(x, t =F (x, t in D st, s =,...,m, v + (x, t v (x, t = on D st, v γ v (x, t γ s + (x, t =f(x, t on D st, v(x, = v (x for x Ω, v γ (x, t =h(x, t on Ω T. Then the following regularity estimate holds: (3. v ΩT [ C v,ω + T /( F L (,T ;L (Ω + h L (,T ;L ( Ω + f L (,T ;L ( D ], where the constant C depends only on the Lipschitz characters of Ω and B s, s =,...,m, and the constants γ s, s =,...,m. Proof. We first recall the trace estimate u L ( Ω C u H (Ω, where the constant C depends only on the Lipschitz character of Ω. See [5] and []. Through integrations by parts over Ω (,t for t T and straightforward
5 H. AMMARI, E. IAKOVLEVA, H. KANG, AND K. KIM calculations, we have t γ ε v = Ω = = t Ω\D t t Ω t Ω γ v v + t m v γ v + ( t v F v Ω t hv fv D t D s t γ s v v D s ( v γ s v γ v + Ω ( t v F v. It then follows from the above trace estimate that t t t t Ω t(v +γ ε v = Fv+ hv fv Ω Ω D t [ F (,τ L (Ω v(,τ L (Ω + h(,τ L ( Ω v(,τ L ( Ω + f(,τ L ( D v(,τ L ( D] dτ F L (,T ;L (Ω v L (,T ;L (Ω + h L (,T ;L ( Ω v L (,T ;L ( Ω + f L (,T ;L ( D v L (,T ;L ( D F L (,T ;L (Ω v L (,T ;L (Ω + C h L (,T ;L ( Ω v L (,T ;H (Ω + C f L (,T ;L ( D v L (,T ;H (D CT /( F L (,T ;L (Ω + h L (,T ;L ( Ω + f L (,T ;L ( D v ΩT for some constant C depending only on the Lipschitz characters of Ω and B s, s =,...,m. Let { } (3. μ = min, min γ s. s m Then, by using Young s inequality, we obtain that μ v Ω T v,ω + CT ( F δ L (,T ;L (Ω + h L (,T ;L ( Ω + f L (,T ;L ( D + 3δ v Ω T, which gives the desired estimate by taking δ small enough. This completes the proof. We start the derivation of the asymptotic formula (.8 for u with the following estimate of u U. Lemma 3.. There exists a constant C such that u U ΩT Cε d, where C = μ T / sup U(x, t γ γ s B s (x,t Ω T and μ is given by (3..
6 DIRECT ALGORITHMS FOR THERMAL IMAGING Proof. Let W = u U. Then W is the solution of [ ] U t W (γ ε W = (γ s γ χ Ds ΔU δ Ds (3.3 W (x, = for x Ω, W γ = on Ω T. in Ω T, Here χ Ds is the characteristic function of the domain D s, and δ Ds is the Dirac-delta function on the boundary of D s. The function W also satisfies W + W = on D st, W γ W γ s + = (γ s γ U on D st. Let f = [ (γ s γ χ Ds ΔU δ Ds U Multiplying the first equation in (3.3 by W, integrating by parts over Ω (,t, and using the divergence theorem and (3.3, we obtain t t ( fw dx dτ = t W (γ ε W Wdxdτ Ω Ω t t = Ω t(w W dx dτ γ ε Ω Wdσ x dτ t + γ ε W Wdxdτ Ω = t W (,t,ω + γ ε W dx dτ. On the other hand, making use of the divergence theorem and Hölder s inequality, we readily get t fw dx dτ = m [ t t (γ s γ ΔUW + (γ γ s U ] Ω D s D s W t = (γ γ s U Wdxdτ D s (γ γ s U(,t,Ds W (,t,ds dτ (γ γ s U,DsT W,DsT Ω ]. γ γ s sup (x,t Ω T U(x, t D s / T / W ΩT = cε d W ΩT,
7 H. AMMARI, E. IAKOVLEVA, H. KANG, AND K. KIM where c = T / sup (x,t ΩT U(x, t m γ γ s B s. Combining the above estimates, we obtain that and hence sup W (,t,ω + t T Ω μ W Ω T cε d W ΩT, γ ε W dx dτ cε d W ΩT, which completes the proof. Lemma 3.3. Let u, U, and Φ be the functions defined in Theorem.. Then the following identity holds: (3.4 I Φ (T = (γ γ s u(x, t Φ(x, t dx dt. D s Proof. Set W = u U. Then integrating by parts and using the conditions Φ t=t = and W t= =, we obtain that ( [ ] t W Φ+γ W Φ = W (x, T Φ(x, T W (x, Φ(x, dx Ω Ω W t Φ+γ W Φ Ω Ω γ W ΔΦ = γ Ω Ω W Φ. On the other hand, we have ( u m t W Φ+γ W Φ = γ Ω (Ω\D Φ+ u γ s D s Φ U + (γ γ s u Φ γ D s Ω Φ = (γ γ s u Φ. D s The proof of identity (3.4 is now complete. 4. Proof of the asymptotic expansion formula. In this section we prove Theorem.. Recall that D = m D s = m (εb s + z s. Let (4. V (x, t =u(x, t U(x, t+ε j= j U(z s,tψ (s j ( x zs ε, where Ψ (s j is defined by (.5.
8 DIRECT ALGORITHMS FOR THERMAL IMAGING 3 Lemma 4.. The following estimate holds: (4. V,ΩT C { ε 5/ if d =3, ε log ε if d = for some constant C independent of ε. Proof. We suppose that there is only a single inclusion D = εb (z = with conductivity γ. The general case of multiple inclusion follows by iteration of the arguments that we will present for the case m =. In other words, we may develop asymptotic formulae involving the difference between the fields u and U with s inclusions and those with s inclusions, s = m,...,, and then at the end essentially form the sum of those m formulae. Observe that V satisfies t V γ ΔV = ε ( x t j U(,tΨ j := F (x, t in (Ω\ D (,T, ε j= t V γδv =(γ ε γ χ D ΔU + ε We also have γ V ( x t j U(,tΨ j := F (x, t in D (,T. ε j= V + V = on D (,T, γ V [ U + =(γ γ j= j U(,t x ] j := G(x, t on D (,T and V (x, = ε V γ = γ ( x j θ(ψ j ε j= j= j U(,t Ψ j in Ω, ( x := H(x, t on Ω (,T. ε Then, Let Ω =ε Ω and let v(x, t :=V (εx, ε t, x Ω, t T/ε. V,Ω T /ε = ε +d v(,t,ω dt, and hence it suffices to show that (4.3 /ε v(,t,ω dt C for some constant C independent of ε.
9 4 H. AMMARI, E. IAKOVLEVA, H. KANG, AND K. KIM Observe that v satisfies t v γ Δv = ε F (εx, ε t:=f(x, t in ( Ω \ B (,T/ε, t v γδv = ε F (εx, ε t:=f(x, t in B (,T/ε, v + v = on B (,T/ε, v γ γ v + = εg(εx, ε t:=g(x, t on B (,T/ε, v(x, = ε j θ(ψ j (x in Ω, j= v γ = εh(εx, ε t:=h(x, t on Ω (,T/ε. Since Ψ j (x =O( x d as x and diam( Ω = O(ε, we have (4.4 Ψ j (x dx = Ω { O(log ε if d =, O( if d =3, and therefore { ε v(x, if d =3, dx C Ω ε log ε if d =. Moreover, we claim that (4.5 /ε f(,t, Ω dt Cε, /ε /ε h(,t, Ω dt Cε, g(,t, B dt Cε, and hence (4.3 follows from Lemma 3.. We now prove (4.5. By definition of f and F,weget /ε /ε f(,t, Ω dt Cε4 ΔU(εx, ε t dx dt + Cε 6 j= B /ε t j U(,ε t dt Ψ j (x dx. Ω Since U is a smooth function, the first inequality in (4.5 follows from (4.4. To prove the second inequality in (4.5, we note that U (x, t j U(,t x j = O( x j= for all (x, t D (,T. Therefore, g(x, t =εg(εx, ε t=o(ε, and hence the second inequality in (4.5 holds.
10 DIRECT ALGORITHMS FOR THERMAL IMAGING 5 For the third inequality, we have /ε h(,t dt ε, Ω j= /ε j U(,t dt Ω Ψ j (x dx. Since Ψ j (x = O( x d as x, we have Ψ j (x = O(ε d for all x Ω. Thus we get /ε h(,t, Ω dt Cεd+. This completes the proof of (4.5 and Lemma 4.. We are now ready to prove Theorem.. Since Ψ (s j (x =O( x d as x and the inclusions are well separated, it follows from (3.4 and (4. that [ ( ] I Φ(T (γ γ s U(x, t j U(z s,t Ψ (s x zs j D s ε j= Φ(x, t dx dt = (γ γ s V (x, t Φ(x, t dx dt D s { O(ε 4 Cε d/ ifd =3, V,ΩT = O(ε 3 log ε ifd =. Since U(x, t = d j= ju(z s,t x j + O(ε and Φ(x, t = d j=i iφ(z s,t x i + O(ε for all x D s, we get the following modulo O(ε 4 ifd = 3 and O(ε 3 log ε if d =: I Φ (T =ε d as desired. m (γ γ s m = ε d (γ γ s m = ε d (γ γ s j= i= j= i= j U(z s,t i Φ(z s,t j U(z s,t i Φ(z s,t U(z s,t M (s Φ(z s,t dt, (x j Ψ (s j B s (s j ψ x i B s (x x i dx dt dσ dt 5. Location search algorithms. In this section we apply formula (.8 in Theorem. (with an appropriate choice of test functions Φ and background solutions U for the purpose of identifying the location of the inclusions D s, s =,...,m. We first design a real-time location search algorithm with good resolution and accuracy for detecting a single inclusion (m =. This algorithm is related to the constant current projection algorithm developed in []. See also [3]. It makes use of constant heat flux and, not surprisingly, it is limited in its ability to effectively locate multiple small inclusions.
11 6 H. AMMARI, E. IAKOVLEVA, H. KANG, AND K. KIM Using many heat sources, we then describe an efficient method to locate multiple inclusions (m > and illustrate its feasibility. Our approach is related to the linear sampling method [9] and the MUSIC algorithm [3, ]. For the sake of simplicity we consider only the two-dimensional case. 5.. Detection of a single inclusion. For y R \ Ω, let (5. Φ(x, t =Φ y (x, t := 4πγ (T t e x y 4γ (T t. The function Φ satisfies ( t + γ ΔΦ = in Ω T and the final condition Φ t=t = in Ω. Suppose that there is only one inclusion D = z + εb with conductivity γ. For simplicity assume that B is a disk. Then it follows from (.7 that the polarization tensor associated with B and γ is given by M = (γ B /(γ + γ I, where I is the identity matrix. Choose the background solution U(x, t (= θ(x to be a harmonic (time-independent function in Ω T. We compute and But and therefore y z Φ y (z,t = 8πγ e (T t z y 4γ (T t, M Φ y (z,t = B z y γ + γ 4πγ (T t e M Φ y (z,t dt = d z y dt e 4γ (T t = B z y T γ + γ 4πγ z y 4γ z y 4γ (T t, e z y 4γ (T t e z y 4γ (T t (T t M Φ y (z,t dt = B z y γ + γ π z y e Then the asymptotic expansion (.8 in Theorem. yields (T t dt. z y 4γ (T t. (5. I Φ (T (y = ε (γ γ θ(z M Φ y (z,t dt + O(ε 3 log ε = ε γ γ γ + γ B θ(z (y z π y z e y z 4γ T + O(ε 3 log ε. Now we are in a position to present our location search algorithms for detecting a single inclusion. First algorithm (with one measurement. We prescribe the initial condition θ(x =a x for some fixed unit constant vector a and choose g = γ a ν as an applied time-independent heat flux on Ω T. Throughout this section, a is
12 DIRECT ALGORITHMS FOR THERMAL IMAGING 7 taken to be a coordinate unit vector. Take two observation lines Σ and Σ contained in R \ Ω such that Σ := a line parallel to a, Σ := a line normal to a. Next we find two points P i Σ i (i =, so that I Φ (T (P =and min I Φ (T (x if γ γ<, x Σ I Φ (T (P = max I Φ (T (x if γ γ>. x Σ Finally, we draw the corresponding lines Π (P and Π (P given by Π (P :={x a (x P =}, Π (P :={x (x P is parallel to a}. Then the intersecting point P of Π (P Π (P is close to the inclusion D. The proof of the following lemma follows from []. We sketch it for the reader s convenience. Lemma 5.. There exists a positive constant C independent of ε and z such that P z Cε log ε for ε small enough. Proof. Fix, for example, γ >γ. From (5. it follows that there exists a constant C, independent of y =(y,, z =(z,z, and ε, such that ( I Φ (T (y ε γ γ B (y z y z γ + γ π y z e 4γ T + Cε log ε and ( I Φ (T (y ε γ γ B (y z y z γ + γ π y z e 4γ T Cε log ε for all y Σ. Therefore, for y =(y, Σ satisfying y z y z y z e 4γ T C(γ + γ ε log ε (γ γ B we have I Φ (T (y. On the other hand, for y Σ satisfying y z y z y z e 4γ T C(γ + γ ε log ε (γ γ B we similarly have I Φ (T (y. This implies that the quantity I Φ (T (y changes sign on Σ for y satisfying C(γ + γ ε log ε (γ γ B y z y z y z e 4γ T C(γ + γ ε log ε (γ γ B for some positive constant C independent of ε. Therefore, the zero point P satisfies (P z a C ε log ε for some constant C depending only on Ω, B, T, γ, and γ. Let y =(,y Σ. The point y =(,z Σ is the unique point where ε γ γ (y z B θ(z γ + γ π y z e y z 4γ T
13 8 H. AMMARI, E. IAKOVLEVA, H. KANG, AND K. KIM attains its maximum. Following [], we can easily see that the point P Σ where I Φ (T (y attains its maximum belongs to an ε-neighborhood of y and we conclude that P =Π (P Π (P satisfies P z Cε log ε for some constant C depending only on Ω, B, T, γ, and γ. Second algorithm (with two measurements. Consider two measurements corresponding to the initial conditions θ (x =a x and θ (x =a x for some fixed constant vector a, where a is a unit vector orthogonal to a, and choose g = γ a ν and g = γ a ν as applied time-independent heat fluxes on Ω T. Here Σ and Σ are as in the first algorithm. Denote by IΦ i (T the function I Φ (T corresponding to θ = θ i, i =,. Then we find two points P i Σ i (i =, so that IΦ i (T (P i =. We draw the corresponding lines Π (P and Π (P given by Π (P :={x a (x P =}, Π (P :={x (x P is parallel to a}. From the proof of Lemma 5. it immediately follows that the intersecting point P of Π (P Π (P is close to the inclusion D: There exists a positive constant C independent of ε and z such that P z Cε log ε for ε small enough. 5.. Detection of multiple inclusions. Consider m inclusions D s = εb s + z s, s =,...,m. Choose U(x, t =U y (x, t := x y 4πγ t e 4γ t for y R \ Ω or, equivalently, g to be the heat flux corresponding to a heat source placed at the point source y and the initial condition θ(x = in Ω, to obtain that I Φ (T = ε m (γ γ s 64π γ 4 (y z s M (s (y z s ( t (T t exp y z s + O(ε 3 log ε, 4γ (T t y z s 4γ t dt where Φ is given by (5.. Suppose for the sake of simplicity that all the domains B s are disks. Then it follows from (.7 that M (s = m (s I, where m (s = γ B s /(γ s + γ. Let the source points y l R \ Ω for l N. We assume that the countable set {y l } l N has the property that any analytic function which vanishes in {y l } l N vanishes identically. Our location search algorithm for detecting multiple inclusions is as follows. For n N sufficiently large, define the matrix A =[A ll ] n l,l = by m A ll := ε (γ γ s 64π γ 4 m (s (y l z s (y l z s ( t (T t exp y l z s 4γ (T t y l z s 4γ t dt.
14 DIRECT ALGORITHMS FOR THERMAL IMAGING 9 For z Ω, we decompose the symmetric real matrix C defined by [ C := ( t (T t exp y l z 4γ (T t y z ] l dt 4γ t as follows: C = p v l (zv l (z l= l,l =,...,n for some p n, where v l R n and vl denotes the transpose of v l. Define the vector g z (l R n for z Ωby ( g z (l = (y zv l (z,...,(y n zv ln (z, l =,...,p. Here v l,...,v ln are the components of the vector v l, l =,...,p. Let y l =(y lx,y ly for l =,...,n, z =(z x,z y, and z s =(z sx,z sy. We also introduce ( g zx (l = (y x z x v l (z,...,(y nx z x v ln (z and ( g zy (l = (y y z y v l (z,...,(y ny z y v ln (z. Lemma 5.. The following characterization of the location of the inclusions in terms of the range of the matrix A holds: (5.3 g (l zx and g (l zy Range(A l {,...,p} iff z {z,...,z m }. Proof. For z Ω suppose that g zx (l and g zy (l Thus, g (l z = for some constants β s (l, which implies that y l z ( t (T t exp y l z = α s Range(A for all l {,...,p}. ( β s (l (y z s v l (z s,...,(y n z s v ln (z s 4γ (T t y l z 4γ t ( y l z s y l z s t (T t exp dt 4γ (T t y l z s 4γ t dt for all l =,...,n, for some constants {α s },...,m. Therefore, provided that the countable set of sources {y l } l N has the property that any analytic function which vanishes in {y l } l N vanishes identically, it follows that y z ( t (T t exp y z y z dt 4γ (T t 4γ t = α s y z s t (T t exp ( y z s 4γ (T t y z s 4γ t dt
15 3 H. AMMARI, E. IAKOVLEVA, H. KANG, AND K. KIM for all y z, y/ {z,...,z m }, which proves that U y (z,t Φ y (z,t dt = α s U y (z s,t Φ y (z s,t dt y z, y / {z,...,z m }. Next, taking the Laplacian in y yields Δ y U y (z,t Φ y (z,t dt = Δ y U y (z s,t Φ y (z s,t dt y z, y / {z,...,z m }. The singularity of Δ y U y(z,t Φ y (z,t dt at the point z, Δ y U y (z,t Φ y (z,t dt = δ y=z + smoother function, gives that z {z,...,z m }, as desired. Let the singular value decomposition (SVD of the matrix A be defined by A = UΣV. Suppose that the vectors {g zx,g (l zy} (l l=,...,p are linearly independent. The rank of A is then equal to p; the first p columns of U, {u,...,u p }, provide a basis for the column space of A, which is denoted by U S ; and the rest of the matrix U, {u p+,u p+,...,u n }, provides a basis for the left null space of A, denoted by U N. From Lemma 5. we have that a test point z coincides with one of the positions z s if and only if P (g z (l c =, for any vector c R \{}, where P = I (U S US is the orthogonal projection onto the null space of A. Thus we can form an image of the locations {z s } m by plotting, at each point z Ω, the quantity W c (l (z := g z (l c / P (g z (l c for l =,...,p. The resulting plot will have large peaks at the locations of z s, s =,...,m. Note that the smallest number n which is sufficient to efficiently recover the inclusions depends on the (unknown number m. This is the main reason to take n sufficiently large Numerical experiments Detection of a single inclusion. In this subsection we present numerical experiments for locating a single inclusion D. We consider Ω as a unit disk centered at the origin in R. The applied time-independent heat flux is given by g = e ν (and g = e ν, g = e ν for the second algorithm on Ω, and the inclusion D is taken as a homogeneous disk of diameter ε =., centered at z =(.43,.7 with conductivity γ = 3. The thermal conductivity of the background γ =. The observation lines are given by Σ = {(x,, x R} and Σ = {(,x, x R}. Figures 5. and 5. show the results of our algorithms. Observing the function I Φ (T (x on the lines Σ and Σ, we obtain the zero point P =(.43, of I Φ (T on Σ and the minimum point (or the zero point for the second algorithm P = (,.7 of I Φ (T onσ, noting that γ γ is negative. Then the intersecting point P of Π (P Π (P in the location search algorithm is given by P =(.43,.7.
16 DIRECT ALGORITHMS FOR THERMAL IMAGING 3 5 x 4.8 x I Φ (x I Φ (x.6.8 P. P x Σ (a (b x Σ.5 x 3.5 I Φ (x.5 P (c Fig. 5.. Patterns of I Φ (T (x calculated on the observation lines Σ = {(x,, x R} and Σ = {(,x, x R}. (a and (b correspond to the first algorithm; (a and (c correspond to the second algorithm. x Σ.8.6 Π.4 Π. Y axis X axis Fig. 5.. Location detection intersecting the lines Π (P and Π (P. The intersecting point is P =(P,P.
17 3 H. AMMARI, E. IAKOVLEVA, H. KANG, AND K. KIM Detection of multiple inclusions. In this subsection we present a numerical experiment for locating multiple inclusions. Let again Ω be a unit disk in R centered at the origin. We assume that we know the values A ll of the pattern I Φ (T for some positive T and for a small finite number of equidistantly distributed source points ( y l = ρ cos πl πl, sin n n, l =,...,n. Here ρ denotes the (minimal distance between the source points y l, l =,...,n, and the origin. We have two cases of interest: ρ>t and ρ T. In the first numerical experiment we take two homogeneous circular disks of diameter ε =. denoted as D and D and respectively centered at z =(.,.3 and z =(.43,.35, to be retrieved using n = source points. Corresponding conductivities γ s, s =,, are equated to 3. The thermal conductivity of the background γ =. Within the above setting, the retrieval of the inclusions involves the calculation of the SVD A = UΣV of the matrix A =[A ll ] R n n and the decomposition of the matrix C. It was observed numerically that the rank of C is equal to one (p =. Therefore, we expect to see four nonzero singular values of the matrix A. Denote by {e i } i= the orthonormal basis in R. Then, for each discrete location z Ω (the sampling step henceforth is h =.3, the identifier of interest is W c (z (:= W c (, c = {e,e,e + e }, calculated within Ω. Plots of z W c (z illustrate the result achieved, sharp peaks being expected to occur at the locations of the inclusions, z s, s =,. Other accompanying results displayed consist of the singular values of A, using a standard log scale. First we take ρ =.5 and T = (first case of interest. The singular values are displayed in Figure 5.3; the identifier W c (z, for c = e, e, and e + e, is displayed in Figure 5.4. Singular Values of A ( 3 4 log (σ j Singular Value Number, σ j Fig ρ>t: distribution of the singular values of A for n =sources. The results obtained are easy to interpret. Four singular values emerge from noise, and we note that two singular values are associated with one specific inclusion. The inclusions, being that they are observed via different W c (z, are clearly discriminated from the background, the visual aspect depending upon the choice of c. In the second numerical example, we take ρ =.5 and T = 3 (the second case of interest and next we take ρ = 4 and T = 3 (the first case of interest.
18 DIRECT ALGORITHMS FOR THERMAL IMAGING Y axis.. Y axis X axis X axis Y axis X axis Fig ρ>t: gray-level (or color map of W c(z, c = e, e, and e + e, ordered from top left to bottom, for all points z in Ω. Figure 5.5 shows the distribution of the singular values of A for two cases. Figure 5.6 shows the identifier W c (z, c = e, for two cases, where the images of W c (z are obtained by using the 5 and 4 largest singular vectors (associated with the 5 and 4 largest singular values of A, respectively in the second case of interest and the 4 largest singular vectors (associated with the 4 largest singular values of A in the first case of interest. The results obtained are evidently less easy to interpret than before. Five singular values emerge from noise in the second case, and four (as in the numerical example above emerge in the first case of interest. The inclusions, being that they are observed via W c (z, c = e, are clearly discriminated from the background (Figure 5.6(b and (c. Notice that in the second case the positions of the inclusions are moving away from the actual positions. The second numerical example involves one more inclusion than the previous one. We add one more homogeneous circular disk of the same diameter ε =. denoted as D 3 and centered at z 3 =(.4,.5, to be retrieved using n = 4 source points, with conductivity γ 3 =3. The result for ρ =4,T = is shown in Figure 5.7 (the distribution of the singular values of A for 4 source points and the identifier W c (z for c = e.
19 34 H. AMMARI, E. IAKOVLEVA, H. KANG, AND K. KIM Singular Values of A ( 5 Singular Values of A ( log (σ j log (σ j Singular Value Number, σ j Singular Value Number, σ j ρ =.5, T =3 ρ =4,T =3 Fig Distribution of the singular values of A for n =sources Y axis.. Y axis X axis (a X axis (b Y axis X axis (c Fig Gray-level (or color map of W c(z, c = e, for all points z in Ω. (a ρ =.5, T =3: the map of W c(z is obtained by using the first 4 singular vectors of A; (b ρ =.5, T =3: the map of W c(z is obtained by using the first 5 singular vectors of A; (c ρ =5, T =3: the map of W c(z is obtained by using the first 4 singular vectors of A.
20 DIRECT ALGORITHMS FOR THERMAL IMAGING 35 8 Singular Values of A ( log (σ j 3 4 Y axis Singular Value Number, σ j (a X axis (b Fig ρ>t: (a distribution of the singular values of A for n =4source points; (b graylevel (or color map of W c(z, c = e, for all points z in Ω. As in the first numerical example, the results obtained are easy to interpret. Six singular values emerge from noise, and we note that two singular values are associated with one specific inclusion. The inclusions, being that they are observed via W c (z, c = e, are clearly discriminated from the background. REFERENCES [] H. Ammari and H. Kang, High-order terms in the asymptotic expansions of the steady-state voltage potentials in the presence of conductivity inhomogeneities of small diameter, SIAM J. Math. Anal., 34 (3, pp [] H. Ammari and H. Kang, Reconstruction of Small Inhomogeneities from Boundary Measurements, Lecture Notes in Math. 846, Springer-Verlag, Berlin, 4. [3] H. Ammari and J. K. Seo, An accurate formula for the reconstruction of conductivity inhomogeneities, Adv. in Appl. Math., 3 (3, pp [4] M. Brühl, M. Hanke, and M. S. Vogelius, A direct impedance tomography algorithm for locating small inhomogeneities, Numer. Math., 93 (3, pp [5] K. Bryan and L. F. Caudill, Jr., An inverse problem in thermal imaging, SIAM J. Appl. Math., 56 (996, pp [6] K. Bryan and L. F. Caudill, Stability and reconstruction for an inverse problem for the heat equation, Inverse Problems, 4 (998, pp [7] W. J. Cantwell and J. Morton, The significance of damage and defects and their detection in composite materials: A review, J. Strain Anal., 7 (99, pp [8] Y. Capdeboscq and M. Vogelius, A general representation formula for boundary voltage perturbations caused by internal conductivity inhomogeneities of low volume fraction, MAN Math. Model. Numer. Anal., 37 (3, pp [9] D. J. Cedio-Fengya, S. Moskow, and M. S. Vogelius, Identification of conductivity imperfections of small diameter by boundary measurements: Continuous dependence and computational reconstruction, Inverse Problems, 4 (998, pp [] R. Chapko, R. Kress, and J. R. Yoon, On the numerical solution of an inverse boundary value problem for the heat equation, Inverse Problems, 4 (998, pp [] R. Chapko, R. Kress, and J. R. Yoon, An inverse boundary value problem for the heat equation: The Neumann condition, Inverse Problems, 5 (999, pp [] M. Cheney, The linear sampling method and the MUSIC algorithm, Inverse Problems, 7 (, pp [3] A. J. Devaney, Super-resolution processing of multi-static data using time reversal and MUSIC, J. Acoust. Soc. Amer., to appear. [4] A. Elayyan and V. Isakov, On uniqueness of recovery of the discontinuous conductivity coefficient of a parabolic equation, SIAM J. Math. Anal., 8 (997, pp
21 36 H. AMMARI, E. IAKOVLEVA, H. KANG, AND K. KIM [5] L. C. Evans and R. E. Gariepy, Measure Theory and Fine Properties of Functions, Stud. Adv. Math., CRC Press, Boca Raton, FL, 99. [6] A. Friedman and M. S. Vogelius, Identification of small inhomogeneities of extreme conductivity by boundary measurements: A theorem on continuous dependence, Arch. Rational Mech. Anal., 5 (989, pp [7] V. Isakov, Inverse Problems for Partial Differential Equations, Springer-Verlag, New York, 998. [8] A. Killey and J. P. Sargent, Analysis of thermal non-destructive testing, J. Phys. D: Appl. Phys., (989, pp [9] A. Kirsch, The MUSIC algorithm and the factorization method in inverse scattering theory for inhomogeneous media, Inverse Problems, 8 (, pp [] O. Kwon, J. K. Seo, and J. R. Yoon, A real-time algorithm for the location search of discontinuous conductivities with one measurement, Comm. Pure Appl. Math., 55 (, pp. 9. [] O. A. Ladyzhenshya, The Boundary Value Problems of Mathematical Physics, Appl. Math. Sci., Springer-Verlag, New York, 985. [] P. M. Patel, S. K. Lau, and D. P. Almond, A review of image analysis techniques applied in transient thermographic nondestructive testing, Nondestructive Testing Eval., 6 (99, pp [3] G. Pólya and G. Szegö, Isoperimetric Inequalities in Mathematical Physics, Ann. of Math. Stud. 7, Princeton University Press, Princeton, NJ, 95.
Determination 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 informationEnhanced resolution in structured media
Enhanced resolution in structured media Eric Bonnetier w/ H. Ammari (Ecole Polytechnique), and Yves Capdeboscq (Oxford) Outline : 1. An experiment of super resolution 2. Small volume asymptotics 3. Periodicity
More informationA note on the MUSIC algorithm for impedance tomography
A note on the MUSIC algorithm for impedance tomography Martin Hanke Institut für Mathematik, Johannes Gutenberg-Universität, 55099 Mainz, Germany E-mail: hanke@math.uni-mainz.de Abstract. We investigate
More informationIDENTIFICATION OF SMALL INHOMOGENEITIES: ASYMPTOTIC FACTORIZATION
IDENTIFICATION OF SMALL INHOMOGENEITIES: ASYMPTOTIC FACTORIZATION HABIB AMMARI, ROLAND GRIESMAIER, AND MARTIN HANKE Abstract. We consider the boundary value problem of calculating the electrostatic potential
More informationSome issues on Electrical Impedance Tomography with complex coefficient
Some issues on Electrical Impedance Tomography with complex coefficient Elisa Francini (Università di Firenze) in collaboration with Elena Beretta and Sergio Vessella E. Francini (Università di Firenze)
More informationEstimation of transmission eigenvalues and the index of refraction from Cauchy data
Estimation of transmission eigenvalues and the index of refraction from Cauchy data Jiguang Sun Abstract Recently the transmission eigenvalue problem has come to play an important role and received a lot
More informationComplex geometrical optics solutions for Lipschitz conductivities
Rev. Mat. Iberoamericana 19 (2003), 57 72 Complex geometrical optics solutions for Lipschitz conductivities Lassi Päivärinta, Alexander Panchenko and Gunther Uhlmann Abstract We prove the existence of
More informationInverse Transport Problems and Applications. II. Optical Tomography and Clear Layers. Guillaume Bal
Inverse Transport Problems and Applications II. Optical Tomography and Clear Layers Guillaume Bal Department of Applied Physics & Applied Mathematics Columbia University http://www.columbia.edu/ gb23 gb23@columbia.edu
More informationand finally, any second order divergence form elliptic operator
Supporting Information: Mathematical proofs Preliminaries Let be an arbitrary bounded open set in R n and let L be any elliptic differential operator associated to a symmetric positive bilinear form B
More 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 informationPiecewise Smooth Solutions to the Burgers-Hilbert Equation
Piecewise Smooth Solutions to the Burgers-Hilbert Equation Alberto Bressan and Tianyou Zhang Department of Mathematics, Penn State University, University Park, Pa 68, USA e-mails: bressan@mathpsuedu, zhang
More informationMultivariable Calculus
2 Multivariable Calculus 2.1 Limits and Continuity Problem 2.1.1 (Fa94) Let the function f : R n R n satisfy the following two conditions: (i) f (K ) is compact whenever K is a compact subset of R n. (ii)
More informationNonlinear stabilization via a linear observability
via a linear observability Kaïs Ammari Department of Mathematics University of Monastir Joint work with Fathia Alabau-Boussouira Collocated feedback stabilization Outline 1 Introduction and main result
More informationA global solution curve for a class of free boundary value problems arising in plasma physics
A global solution curve for a class of free boundary value problems arising in plasma physics Philip Korman epartment of Mathematical Sciences University of Cincinnati Cincinnati Ohio 4522-0025 Abstract
More informationLocal null controllability of the N-dimensional Navier-Stokes system with N-1 scalar controls in an arbitrary control domain
Local null controllability of the N-dimensional Navier-Stokes system with N-1 scalar controls in an arbitrary control domain Nicolás Carreño Université Pierre et Marie Curie-Paris 6 UMR 7598 Laboratoire
More informationA MUSIC-type Algorithm for Detecting Internal Corrosion from Electrostatic Boundary Measurements
A MUSIC-type Algorithm for Detecting Internal Corrosion from Electrostatic Boundary Measurements Habib Ammari Hyeonbae Kang Eunjoo Kim Kaouthar Louati Michael S. Vogelius October 24, 27 Abstract We establish
More informationSome recent results on controllability of coupled parabolic systems: Towards a Kalman condition
Some recent results on controllability of coupled parabolic systems: Towards a Kalman condition F. Ammar Khodja Clermont-Ferrand, June 2011 GOAL: 1 Show the important differences between scalar and non
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More information1 Math 241A-B Homework Problem List for F2015 and W2016
1 Math 241A-B Homework Problem List for F2015 W2016 1.1 Homework 1. Due Wednesday, October 7, 2015 Notation 1.1 Let U be any set, g be a positive function on U, Y be a normed space. For any f : U Y let
More informationThe first order quasi-linear PDEs
Chapter 2 The first order quasi-linear PDEs The first order quasi-linear PDEs have the following general form: F (x, u, Du) = 0, (2.1) where x = (x 1, x 2,, x 3 ) R n, u = u(x), Du is the gradient of u.
More informationSHARP BOUNDARY TRACE INEQUALITIES. 1. Introduction
SHARP BOUNDARY TRACE INEQUALITIES GILES AUCHMUTY Abstract. This paper describes sharp inequalities for the trace of Sobolev functions on the boundary of a bounded region R N. The inequalities bound (semi-)norms
More informationRelaxation methods and finite element schemes for the equations of visco-elastodynamics. Chiara Simeoni
Relaxation methods and finite element schemes for the equations of visco-elastodynamics Chiara Simeoni Department of Information Engineering, Computer Science and Mathematics University of L Aquila (Italy)
More informationA CALDERÓN PROBLEM WITH FREQUENCY-DIFFERENTIAL DATA IN DISPERSIVE MEDIA. has a unique solution. The Dirichlet-to-Neumann map
A CALDERÓN PROBLEM WITH FREQUENCY-DIFFERENTIAL DATA IN DISPERSIVE MEDIA SUNGWHAN KIM AND ALEXANDRU TAMASAN ABSTRACT. We consider the problem of identifying a complex valued coefficient γ(x, ω) in the conductivity
More informationDifferential equations, comprehensive exam topics and sample questions
Differential equations, comprehensive exam topics and sample questions Topics covered ODE s: Chapters -5, 7, from Elementary Differential Equations by Edwards and Penney, 6th edition.. Exact solutions
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 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 informationThe Method of Small Volume Expansions for Emerging Medical Imaging
Vienna p. 1/5 The Method of Small Volume Expansions for Emerging Medical Imaging Habib Ammari CNRS & Ecole Polytechnique Vienna p. 2/5 Motivation and Principles of the MSVE Inverse problems in medical
More informationREMARKS ON THE VANISHING OBSTACLE LIMIT FOR A 3D VISCOUS INCOMPRESSIBLE FLUID
REMARKS ON THE VANISHING OBSTACLE LIMIT FOR A 3D VISCOUS INCOMPRESSIBLE FLUID DRAGOŞ IFTIMIE AND JAMES P. KELLIHER Abstract. In [Math. Ann. 336 (2006), 449-489] the authors consider the two dimensional
More informationin Bounded Domains Ariane Trescases CMLA, ENS Cachan
CMLA, ENS Cachan Joint work with Yan GUO, Chanwoo KIM and Daniela TONON International Conference on Nonlinear Analysis: Boundary Phenomena for Evolutionnary PDE Academia Sinica December 21, 214 Outline
More informationA LOWER BOUND ON BLOWUP RATES FOR THE 3D INCOMPRESSIBLE EULER EQUATION AND A SINGLE EXPONENTIAL BEALE-KATO-MAJDA ESTIMATE. 1.
A LOWER BOUND ON BLOWUP RATES FOR THE 3D INCOMPRESSIBLE EULER EQUATION AND A SINGLE EXPONENTIAL BEALE-KATO-MAJDA ESTIMATE THOMAS CHEN AND NATAŠA PAVLOVIĆ Abstract. We prove a Beale-Kato-Majda criterion
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 informationAnalysis of the anomalous localized resonance
Hyundae Lee(Inha University, Korea) Joint work with Habib Ammari, Giulio Ciraolo, Hyeonbae Kang, Graeme Milton. UCI, June 22, 2012 A conference on inverse problems in honor of Gunther Uhlmann Outline Introduction
More informationLocal semiconvexity of Kantorovich potentials on non-compact manifolds
Local semiconvexity of Kantorovich potentials on non-compact manifolds Alessio Figalli, Nicola Gigli Abstract We prove that any Kantorovich potential for the cost function c = d / on a Riemannian manifold
More informationSome asymptotic properties of solutions for Burgers equation in L p (R)
ARMA manuscript No. (will be inserted by the editor) Some asymptotic properties of solutions for Burgers equation in L p (R) PAULO R. ZINGANO Abstract We discuss time asymptotic properties of solutions
More informationClass Meeting # 1: Introduction to PDEs
MATH 18.152 COURSE NOTES - CLASS MEETING # 1 18.152 Introduction to PDEs, Spring 2017 Professor: Jared Speck Class Meeting # 1: Introduction to PDEs 1. What is a PDE? We will be studying functions u =
More informationHOMEOMORPHISMS OF BOUNDED VARIATION
HOMEOMORPHISMS OF BOUNDED VARIATION STANISLAV HENCL, PEKKA KOSKELA AND JANI ONNINEN Abstract. We show that the inverse of a planar homeomorphism of bounded variation is also of bounded variation. In higher
More information1. Introduction. Consider the following parameterized optimization problem:
SIAM J. OPTIM. c 1998 Society for Industrial and Applied Mathematics Vol. 8, No. 4, pp. 940 946, November 1998 004 NONDEGENERACY AND QUANTITATIVE STABILITY OF PARAMETERIZED OPTIMIZATION PROBLEMS WITH MULTIPLE
More informationu xx + u yy = 0. (5.1)
Chapter 5 Laplace Equation The following equation is called Laplace equation in two independent variables x, y: The non-homogeneous problem u xx + u yy =. (5.1) u xx + u yy = F, (5.) where F is a function
More informationBIHARMONIC WAVE MAPS INTO SPHERES
BIHARMONIC WAVE MAPS INTO SPHERES SEBASTIAN HERR, TOBIAS LAMM, AND ROLAND SCHNAUBELT Abstract. A global weak solution of the biharmonic wave map equation in the energy space for spherical targets is constructed.
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 information{ } is an asymptotic sequence.
AMS B Perturbation Methods Lecture 3 Copyright by Hongyun Wang, UCSC Recap Iterative method for finding asymptotic series requirement on the iteration formula to make it work Singular perturbation use
More informationThe Factorization Method for the Reconstruction of Inclusions
The Factorization Method for the Reconstruction of Inclusions Martin Hanke Institut für Mathematik Johannes Gutenberg-Universität Mainz hanke@math.uni-mainz.de January 2007 Overview Electrical Impedance
More informationThe incompressible Navier-Stokes equations in vacuum
The incompressible, Université Paris-Est Créteil Piotr Bogus law Mucha, Warsaw University Journées Jeunes EDPistes 218, Institut Elie Cartan, Université de Lorraine March 23th, 218 Incompressible Navier-Stokes
More informationON BOUNDEDNESS OF MAXIMAL FUNCTIONS IN SOBOLEV SPACES
Annales Academiæ Scientiarum Fennicæ Mathematica Volumen 29, 2004, 167 176 ON BOUNDEDNESS OF MAXIMAL FUNCTIONS IN SOBOLEV SPACES Piotr Haj lasz and Jani Onninen Warsaw University, Institute of Mathematics
More informationBoundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis
Boundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis T. Tran Q. T. Le Gia I. H. Sloan E. P. Stephan Abstract Radial basis functions are used to define approximate solutions
More informationfor all subintervals I J. If the same is true for the dyadic subintervals I D J only, we will write ϕ BMO d (J). In fact, the following is true
3 ohn Nirenberg inequality, Part I A function ϕ L () belongs to the space BMO() if sup ϕ(s) ϕ I I I < for all subintervals I If the same is true for the dyadic subintervals I D only, we will write ϕ BMO
More informationASYMPTOTIC STRUCTURE FOR SOLUTIONS OF THE NAVIER STOKES EQUATIONS. Tian Ma. Shouhong Wang
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS Volume 11, Number 1, July 004 pp. 189 04 ASYMPTOTIC STRUCTURE FOR SOLUTIONS OF THE NAVIER STOKES EQUATIONS Tian Ma Department of
More informationReconstruction Scheme for Active Thermography
Reconstruction Scheme for Active Thermography Gen Nakamura gnaka@math.sci.hokudai.ac.jp Department of Mathematics, Hokkaido University, Japan Newton Institute, Cambridge, Sept. 20, 2011 Contents.1.. Important
More informationGlobal well-posedness of the primitive equations of oceanic and atmospheric dynamics
Global well-posedness of the primitive equations of oceanic and atmospheric dynamics Jinkai Li Department of Mathematics The Chinese University of Hong Kong Dynamics of Small Scales in Fluids ICERM, Feb
More informationMATH 425, FINAL EXAM SOLUTIONS
MATH 425, FINAL EXAM SOLUTIONS Each exercise is worth 50 points. Exercise. a The operator L is defined on smooth functions of (x, y by: Is the operator L linear? Prove your answer. L (u := arctan(xy u
More informationEXISTENCE AND UNIQUENESS FOR A THREE DIMENSIONAL MODEL OF FERROMAGNETISM
1 EXISTENCE AND UNIQUENESS FOR A THREE DIMENSIONAL MODEL OF FERROMAGNETISM V. BERTI and M. FABRIZIO Dipartimento di Matematica, Università degli Studi di Bologna, P.zza di Porta S. Donato 5, I-4126, Bologna,
More informationTHE HOT SPOTS CONJECTURE FOR NEARLY CIRCULAR PLANAR CONVEX DOMAINS
THE HOT SPOTS CONJECTURE FOR NEARLY CIRCULAR PLANAR CONVEX DOMAINS YASUHITO MIYAMOTO Abstract. We prove the hot spots conjecture of J. Rauch in the case that the domain Ω is a planar convex domain satisfying
More informationHOW TO APPROXIMATE THE HEAT EQUATION WITH NEUMANN BOUNDARY CONDITIONS BY NONLOCAL DIFFUSION PROBLEMS. 1. Introduction
HOW TO APPROXIMATE THE HEAT EQUATION WITH NEUMANN BOUNDARY CONDITIONS BY NONLOCAL DIFFUSION PROBLEMS CARMEN CORTAZAR, MANUEL ELGUETA, ULIO D. ROSSI, AND NOEMI WOLANSKI Abstract. We present a model for
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 informationBLOWUP THEORY FOR THE CRITICAL NONLINEAR SCHRÖDINGER EQUATIONS REVISITED
BLOWUP THEORY FOR THE CRITICAL NONLINEAR SCHRÖDINGER EQUATIONS REVISITED TAOUFIK HMIDI AND SAHBI KERAANI Abstract. In this note we prove a refined version of compactness lemma adapted to the blowup analysis
More informationMAT 257, Handout 13: December 5-7, 2011.
MAT 257, Handout 13: December 5-7, 2011. The Change of Variables Theorem. In these notes, I try to make more explicit some parts of Spivak s proof of the Change of Variable Theorem, and to supply most
More informationWELL-POSEDNESS OF THE TWO-DIMENSIONAL GENERALIZED BENJAMIN-BONA-MAHONY EQUATION ON THE UPPER HALF PLANE
WELL-POSEDNESS OF THE TWO-DIMENSIONAL GENERALIZED BENJAMIN-BONA-MAHONY EQUATION ON THE UPPER HALF PLANE YING-CHIEH LIN, C. H. ARTHUR CHENG, JOHN M. HONG, JIAHONG WU, AND JUAN-MING YUAN Abstract. This paper
More informationThe heat equation in time dependent domains with Neumann boundary conditions
The heat equation in time dependent domains with Neumann boundary conditions Chris Burdzy Zhen-Qing Chen John Sylvester Abstract We study the heat equation in domains in R n with insulated fast moving
More informationHerz (cf. [H], and also [BS]) proved that the reverse inequality is also true, that is,
REARRANGEMENT OF HARDY-LITTLEWOOD MAXIMAL FUNCTIONS IN LORENTZ SPACES. Jesús Bastero*, Mario Milman and Francisco J. Ruiz** Abstract. For the classical Hardy-Littlewood maximal function M f, a well known
More informationA Singular Integral Transform for the Gibbs-Wilbraham Effect in Inverse Fourier Transforms
Universal Journal of Integral Equations 4 (2016), 54-62 www.papersciences.com A Singular Integral Transform for the Gibbs-Wilbraham Effect in Inverse Fourier Transforms N. H. S. Haidar CRAMS: Center for
More informationDS-GA 1002 Lecture notes 10 November 23, Linear models
DS-GA 2 Lecture notes November 23, 2 Linear functions Linear models A linear model encodes the assumption that two quantities are linearly related. Mathematically, this is characterized using linear functions.
More informationTHE L 2 -HODGE THEORY AND REPRESENTATION ON R n
THE L 2 -HODGE THEORY AND REPRESENTATION ON R n BAISHENG YAN Abstract. We present an elementary L 2 -Hodge theory on whole R n based on the minimization principle of the calculus of variations and some
More informationInverse Brascamp-Lieb inequalities along the Heat equation
Inverse Brascamp-Lieb inequalities along the Heat equation Franck Barthe and Dario Cordero-Erausquin October 8, 003 Abstract Adapting Borell s proof of Ehrhard s inequality for general sets, we provide
More informationare harmonic functions so by superposition
J. Rauch Applied Complex Analysis The Dirichlet Problem Abstract. We solve, by simple formula, the Dirichlet Problem in a half space with step function boundary data. Uniqueness is proved by complex variable
More informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More informationOBSERVABILITY INEQUALITY AND DECAY RATE FOR WAVE EQUATIONS WITH NONLINEAR BOUNDARY CONDITIONS
Electronic Journal of Differential Equations, Vol. 27 (27, No. 6, pp. 2. ISSN: 72-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu OBSERVABILITY INEQUALITY AND DECAY RATE FOR WAVE EQUATIONS
More informationFinite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product
Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )
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 informationLecture notes: Introduction to Partial Differential Equations
Lecture notes: Introduction to Partial Differential Equations Sergei V. Shabanov Department of Mathematics, University of Florida, Gainesville, FL 32611 USA CHAPTER 1 Classification of Partial Differential
More informationClassical solutions for the quasi-stationary Stefan problem with surface tension
Classical solutions for the quasi-stationary Stefan problem with surface tension Joachim Escher, Gieri Simonett We show that the quasi-stationary two-phase Stefan problem with surface tension has a unique
More informationEXISTENCE OF SOLUTIONS TO THE CAHN-HILLIARD/ALLEN-CAHN EQUATION WITH DEGENERATE MOBILITY
Electronic Journal of Differential Equations, Vol. 216 216), No. 329, pp. 1 22. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE OF SOLUTIONS TO THE CAHN-HILLIARD/ALLEN-CAHN
More informationAsymptotic Spectral Imaging
Rennes p. /4 Asymptotic Spectral Imaging Habib Ammari CNRS & Ecole Polytechnique Rennes p. 2/4 Motivation Image internal defects from eigenvalue and eigenvector boundary measurements (or spectral measurements).
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 informationStability of an abstract wave equation with delay and a Kelvin Voigt damping
Stability of an abstract wave equation with delay and a Kelvin Voigt damping University of Monastir/UPSAY/LMV-UVSQ Joint work with Serge Nicaise and Cristina Pignotti Outline 1 Problem The idea Stability
More informationHomogenization and error estimates of free boundary velocities in periodic media
Homogenization and error estimates of free boundary velocities in periodic media Inwon C. Kim October 7, 2011 Abstract In this note I describe a recent result ([14]-[15]) on homogenization and error estimates
More informationarxiv:math/ v1 [math.ap] 28 Oct 2005
arxiv:math/050643v [math.ap] 28 Oct 2005 A remark on asymptotic completeness for the critical nonlinear Klein-Gordon equation Hans Lindblad and Avy Soffer University of California at San Diego and Rutgers
More information13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs)
13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs) A prototypical problem we will discuss in detail is the 1D diffusion equation u t = Du xx < x < l, t > finite-length rod u(x,
More informationA TWO PARAMETERS AMBROSETTI PRODI PROBLEM*
PORTUGALIAE MATHEMATICA Vol. 53 Fasc. 3 1996 A TWO PARAMETERS AMBROSETTI PRODI PROBLEM* C. De Coster** and P. Habets 1 Introduction The study of the Ambrosetti Prodi problem has started with the paper
More informationOn Smoothness of Suitable Weak Solutions to the Navier-Stokes Equations
On Smoothness of Suitable Weak Solutions to the Navier-Stokes Equations G. Seregin, V. Šverák Dedicated to Vsevolod Alexeevich Solonnikov Abstract We prove two sufficient conditions for local regularity
More informationREMARKS ON THE MEMBRANE AND BUCKLING EIGENVALUES FOR PLANAR DOMAINS. Leonid Friedlander
REMARKS ON THE MEMBRANE AND BUCKLING EIGENVALUES FOR PLANAR DOMAINS Leonid Friedlander Abstract. I present a counter-example to the conjecture that the first eigenvalue of the clamped buckling problem
More informationA uniqueness result and image reconstruction of the orthotropic conductivity in magnetic resonance electrical impedance tomography
c de Gruyter 28 J. Inv. Ill-Posed Problems 16 28), 381 396 DOI 1.1515 / JIIP.28.21 A uniqueness result and image reconstruction of the orthotropic conductivity in magnetic resonance electrical impedance
More informationSébastien Chaumont a a Institut Élie Cartan, Université Henri Poincaré Nancy I, B. P. 239, Vandoeuvre-lès-Nancy Cedex, France. 1.
A strong comparison result for viscosity solutions to Hamilton-Jacobi-Bellman equations with Dirichlet condition on a non-smooth boundary and application to parabolic problems Sébastien Chaumont a a Institut
More informationOn the infinity Laplace operator
On the infinity Laplace operator Petri Juutinen Köln, July 2008 The infinity Laplace equation Gunnar Aronsson (1960 s): variational problems of the form S(u, Ω) = ess sup H (x, u(x), Du(x)). (1) x Ω The
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 information1. INTRODUCTION 2. PROBLEM FORMULATION ROMAI J., 6, 2(2010), 1 13
Contents 1 A product formula approach to an inverse problem governed by nonlinear phase-field transition system. Case 1D Tommaso Benincasa, Costică Moroşanu 1 v ROMAI J., 6, 2(21), 1 13 A PRODUCT FORMULA
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 informationAuthor(s) Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 41(1)
Title On the stability of contact Navier-Stokes equations with discont free b Authors Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 4 Issue 4-3 Date Text Version publisher URL
More informationOn the p-laplacian and p-fluids
LMU Munich, Germany Lars Diening On the p-laplacian and p-fluids Lars Diening On the p-laplacian and p-fluids 1/50 p-laplacian Part I p-laplace and basic properties Lars Diening On the p-laplacian and
More informationTopological Derivative
Outline Theory & Outline About the lecturer 2009 Ph.D. in Mathematics at the Faculty of Mathematics, University of Lodz, Poland 2010 Ph.D. in Applied Mathematics at the Elie Cartan Institute, Université
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 15 Heat with a source So far we considered homogeneous wave and heat equations and the associated initial value problems on the whole line, as
More informationarxiv: v2 [math.ap] 9 Sep 2014
Localized and complete resonance in plasmonic structures Hoai-Minh Nguyen and Loc Hoang Nguyen April 23, 2018 arxiv:1310.3633v2 [math.ap] 9 Sep 2014 Abstract This paper studies a possible connection between
More informationarxiv:math/ v1 [math.dg] 7 Jun 2004
arxiv:math/46v [math.dg] 7 Jun 4 The First Dirichlet Eigenvalue and Li s Conjecture Jun LING Abstract We give a new estimate on the lower bound for the first Dirichlet eigenvalue for the compact manifolds
More informationAnomalous transport of particles in Plasma physics
Anomalous transport of particles in Plasma physics L. Cesbron a, A. Mellet b,1, K. Trivisa b, a École Normale Supérieure de Cachan Campus de Ker Lann 35170 Bruz rance. b Department of Mathematics, University
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 informationFactorization method in inverse
Title: Name: Affil./Addr.: Factorization method in inverse scattering Armin Lechleiter University of Bremen Zentrum für Technomathematik Bibliothekstr. 1 28359 Bremen Germany Phone: +49 (421) 218-63891
More informationA Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions
A Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions Zhiqiang Cai Seokchan Kim Sangdong Kim Sooryun Kong Abstract In [7], we proposed a new finite element method
More informationExistence of minimizers for the pure displacement problem in nonlinear elasticity
Existence of minimizers for the pure displacement problem in nonlinear elasticity Cristinel Mardare Université Pierre et Marie Curie - Paris 6, Laboratoire Jacques-Louis Lions, Paris, F-75005 France Abstract
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 informationIn this chapter we study elliptical PDEs. That is, PDEs of the form. 2 u = lots,
Chapter 8 Elliptic PDEs In this chapter we study elliptical PDEs. That is, PDEs of the form 2 u = lots, where lots means lower-order terms (u x, u y,..., u, f). Here are some ways to think about the physical
More informationOn an uniqueness theorem for characteristic functions
ISSN 392-53 Nonlinear Analysis: Modelling and Control, 207, Vol. 22, No. 3, 42 420 https://doi.org/0.5388/na.207.3.9 On an uniqueness theorem for characteristic functions Saulius Norvidas Institute of
More information