A CALDERÓN PROBLEM WITH FREQUENCY-DIFFERENTIAL DATA IN DISPERSIVE MEDIA. has a unique solution. The Dirichlet-to-Neumann map
|
|
- Dwain Summers
- 5 years ago
- Views:
Transcription
1 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 equation γ(, ω) u(, ω) = 0 from knowledge of the frequency differentials of the Dirichlet-to-Neumann map. For a frequency analytic γ(, ω) = k=0 (σ k + iϵ k )ω k, in three dimensions and higher, we show that d j dω j Λ γ(,ω) ω=0 for j = 0, 1,..., N recovers σ 0,...σ N and ϵ 1,,,, ϵ N. This problem arises in frequency differential electrical impedance tomography of dispersive media. 1. INTRODUCTION Let R n, n 2, be a bounded domain with Lipschitz boundary. For some ω 0 > 0 and c > 0, let γ C([0, ω 0 ]; L ()) be complex valued satisfying (1) Re(γ(x, ω)) c > 0, (x, ω) [0, ω 0 ]. For each fixed ω and f H 1/2 () the Dirichlet problem (2) γ(, ω) u(, ω) = 0 in, u(, ω) = f, has a unique solution. The Dirichlet-to-Neumann map is a bounded operator defined by Λ γ(,ω) : H 1/2 () H 1/2 () (3) Λ γ(,ω) (f) := n γ(, ω) u(, ω), where n is the outer normal to the boundary. Moreover the map ω Λ γ(,ω) C([0, ω 0 ]; L(H 1 2 (), H 1 2 ())), where L(H 1 2 (), H 1 2 ()) denotes the space of linear operators endowed with the strong topology. The Calderón problem refers to the determination of γ(, ω) from Λ γ(,ω), for each fixed frequency ω. At ω = 0 this is the original question in 2000 Mathematics Subject Classification. 35R30, 35J65, 65N21. Key words and phrases. Calderón problem, frequency differential electrical impedance tomography, complex geometrical optics. This author was supported in part by the NSF Grant DMS
2 2 SUNGWHAN KIM AND ALEXANDRU TAMASAN [6], which has been mostly settled in the affirmative: In dimensions n 3, of direct relevance here, we mention the breakthrough result in [22], where Λ γ(,0) is shown to uniquely determine γ(, 0), and the reconstruction method in [16]. Although not stated, these results extend to ω 0, case in which the coefficient γ(, ω) becomes complex valued. The analogous two dimensional problem (also at ω = 0) was cracked in [17], refined in [3] (with reconstruction in [14]) and completely settled in [1]. The breakthrough result in [5] settles in two dimensions the uniqueness part: Λ γ(,ω) uniquely determines γ(, ω) for each fixed ω (an earlier result [9] assumed γ = σ 0 (x) + iωϵ 0 (x) with sufficiently small ωϵ 0 ). For an understanding on the breadth of development of the Calderón problem we refer to the reviews [7], [4], [10], [2], and [23]. Recent biomedical advances [11, 12, 13, 18, 19] allow to measure discrepancies in boundary data at different frequencies, and produce images based on such differential data. Given that biological materials have frequency dependent electric-magnetic properties (they are dispersive) [8] and only one image is produced from multiple frequencies, the following questions arise: What can be determined from the frequency derivative d k dω Λ k γ(,ω) ω=ω0 for a fixed order k 0, or for several orders? Can one expect quantitative information from multi-frequency data? What do the images represent? We address these questions when the coefficient γ(, ω) depends analytically in ω in dimensions three or higher. In [21] the authors considered the case k = 1 and the ansatz γ = σ 0 (x) + iωϵ 1 (x) specific to media with frequency-independent electrical properties (i.e., non-dispersive). For some r, c > 0, let us define the class A r,c of functions γ(x, ω) such that (4) Re(γ(x, ω)) = σ n (x)ω n, Im(γ(x, ω)) = ϵ n (x)ω n, n=0 where σ 0 C 1,1 () is real valued and constant near the boundary with (5) 0 < c 1 σ 0 (x) c, and σ k, ϵ k C 1,1 0 () are real valued and compactly supported in with σ k r, ϵ k (6) r, k 1. σ 0 Note that ω γ(, ω) is analytic in ω < 1/r with values in L (), We prove the following: σ 0 n=1
3 A CALDERÓN PROBLEM FOR DISPERSIVE MEDIA 3 Theorem 1.1. Let n 3 and R n be open with C 1,1 boundary. Assume that γ A r,c, for some r, c > 0 and N 1 be fixed. Then dk ω=0 Λ dω k γ(,ω) for k = 0,..., N, uniquely determine σ 0, σ 1,..., σ N, and ϵ 1,..., ϵ N inside. Consequently, we also obtain: Corollary 1.1. Let n 3 and R n be open with C 1,1 boundary. Assume γ A r,c, for some r, c > 0. Then the map ω Λ γ(,ω) given in an arbitrarily small neighborhood of some fixed ω 0 [0, (1 + 2r) 1 ) recovers γ on [0, r 1 ). The proof of Theorem 1.1 is based on the complex geometrical optic solutions of Sylvester-Uhlmann in [22], whose existence is recalled below to establish notations. For k, η, l R n with k η = k l = k η = 0, and η 2 = k 2 + l 2, consider the vectors 4 ( ) ( ) k k (7) ξ 1 (η, k, l) := η i 2 + l, ξ 2 (η, k, l) := η i 2 l, ( which satisfy ξ 1 ξ 1 = ξ 2 ξ 2 = 0, ξ 1 2 = ξ 2 2 k = l ), 2 and 2 ξ 1 + ξ 2 = ik. The result is stated in the equivalent form: Existence of CGOs (Theorem 2.3, [22]): Let n 3, and σ 0 C 1,1 (R n ) be constant outside. For 1 < δ < 0 there are two constants R, C > 0 dependent only on δ, σ 0 / σ 0 L (), and such that the following holds: For ξ j C n as in (7) with l > R, there exist w j := w(, ξ j ) Hloc 1 (Rn ) solutions in R n of of the form σ 0 w j = 0 (8) w j (x) = e x ξ j σ 1/2 0 (1 + ψ j (x)), with (9) ψ j L 2 δ C, ξ j j = 1, 2, where f 2 := f(x) 2 (1 + x 2 ) δ dx. L 2 δ R n The assumptions on the coefficients may be relaxed as follows. The boundary values σ 0 and σ 0 can be recovered from Λ γ(,0) as shown in [16] for this regularity (and earlier in [15] and [22] for C -regular). Then σ 0 can be extended outside with preserved regularity while making it constant near the boundary of a neighborhood of the original domain. Similarly, the compact support assumption on σ n, ϵ n, for n = 1, 2,... can σ be replaced by the knowledge of the boundary values of σ n,ϵ n, n, and ϵ n. With the same proof as in [17], the Dirichlet-to-Neumann map
4 4 SUNGWHAN KIM AND ALEXANDRU TAMASAN Λ γ(,ω) transfers from to and we can consider the problem on the larger domain. 2. A FREQUENCY RANGE FOR ANALYTIC DEPENDENCE In this section we show analytic dependence of the map ω Λ γ(,ω) in the space of linear operators from H 1/2 () to H 1/2 () endowed with the strong topology, and introduce two operators needed in the proof of Theorem 1.1. Existence and uniqueness of solution for the conductivity problem (2) with a complex coefficient has been known for some time; e.g., in [20] the Lax-Milgram lemma is employed. Below we provide a solution which is frequency-analytic. Since the solution of (2) is unique, the result below also shows the analytic dependence of the solution of (2) for a frequencyanalytic coefficient. Theorem 2.1. Let R n, n 2, be a domain with Lipschitz boundary. Consider the Dirichlet problem (2) for f H 1/2 () real valued and (10) γ(x, ω) := σ 0 (x) + (σ k (x) + iϵ k (x))ω n, k=1 where σ 0 is real valued satisfying (5) for some c > 0, and σ k, ϵ k L (), k 1, are real valued satisfying (6) for some r > 0. Then, for (11) ω < r, the Dirichlet problem (2) has a unique solution u(, ω) H 1 () with the series representation (12) Re{u(x, ω)} = v k (x)ω k, Im{u(x, ω)} = h k (x)ω k k=0 absolutely convergent in H 1 (), where h 0 0, v 0 solves (13) σ 0 v 0 = 0 in, v 0 = f, and, for k = 1, 2,..., the following recurrence holds: (14) (15) σ 0 v k = σ 0 h k = k=0 (σ j v k j ϵ j h k j ) in, v k = 0, (σ j h k j + ϵ j v k j ) in, h k = 0.
5 A CALDERÓN PROBLEM FOR DISPERSIVE MEDIA 5 We stress here that the recurrence (14) and (15) uses the reality of the Dirichlet boundary data f. The proof of the theorem is a consequence of the estimate in Lemma 2.1 below. We start by a simple remark, which can be justified by induction. Remark 2.1. Let {a k } k=0 be a sequence of non-negatives and m > 0 such that a 1 a 0 and Then a k a 0 + m(a a k 1 ), k 2. a k (1 + m) k a 0, k 0. We denote by the L 2 ()-norm, and by 1 the H 1 ()-norm. Lemma 2.1. Let v 0 satisfy (13) and assume the sequences {v n } n=1, {h n } n=1 in H 1 () satisfy the corresponding problems (14) and (15). Then, for k = 1, 2,... we have (16) v k 1 + h k 1 c(1 + 2r) k v 0. Proof. Let k be a fixed positive index. Since v k = 0, by applying the divergence theorem to (14) we estimate σ 0 v k 2 = (σ j v k j ϵ j h k j ) v k dx r r σ 0 v k σ 0 v k j v k + σ 0 h k j v k dx ( σ 0 v k j + σ 0 h k j ), where the first inequality uses (6). Thus σ 0 v k r ( σ 0 v k j + σ 0 h k j ). Similarly, using h k = 0 and the divergence theorem in (15) we obtain σ 0 h k r ( σ 0 v k j + σ 0 h k j ). Now apply Remark 2.1 for a k = σ 0 v k + σ 0 h k and m = 2r to obtain
6 6 SUNGWHAN KIM AND ALEXANDRU TAMASAN (17) σ 0 v k + σ 0 h k (1 + 2r) k σ 0 v 0. Since v k, h k H 1 0() the estimate (16) follows from the lower and upper bounds on σ 0 in (5). Proof. [of Theorem 2.1] We seek the solution in the ansatz (18) u(x, ω) := v 0 (x) + (v n (x) + ih n (x))ω n n=1 and then the recurrence in (14) and (15) is obtained by identifying the same order coefficients in a formal series. It is in here that the reality of the Dirichlet boundary data f (hence of v 0 ) is used. We construct the sequences {v n }, {h n } inductively by solving the Dirichlet problems (14) and (15). Part of the inductive step, the right hand hand sides (in divergence form) of the equations in (14) and (15) belong to H 1 (). Thus the recurrence is well defined to construct v 0 H 1 () and the two sequences {v n } n=1 and {h n } n=1 in H 1 0(). The estimate (16) and ω < 1/(1 + 2r) shows (19) u 1 c v 0 1 n=0 (1 + 2r) k ω k <. Therefore u(, ω) defined in (18) is in H 1 (). A direct calculation involving (14) and (15) shows that it solves the Dirichlet problem (2) for γ(x, ω) := σ 0 (x) + (σ k (x) + iϵ k (x))ω n. k=1 As a consequence of the theorem above we can explicit the real and imaginary part of the operators dk dω k Λ γ(,ω) ω=0 (f) for a real valued f as follows. Corollary 2.1. Let γ(x, ω) be as in the Theorem 2.1 and f H 1/2 () be real valued. Then, with h 0 0, v 0 solution of (2), and v k and h k, k = 1, 2,... given by the recurrences (14) and (15), we obtain (20) (21) { } 1 d k k! Re dω Λ k γ(,ω) (f) = ω=0 { } 1 d k k! Im dω Λ k γ(,ω) (f) = ω=0 ( σ j v k j ( σ j h k j ϵ j + ϵ j h k j ), ) v k j.
7 A CALDERÓN PROBLEM FOR DISPERSIVE MEDIA 7 Proof. For ω < (1 + 2r) 1 < r 1 the series in (10) defining γ is absolutely convergent in L (), and the series (18) defining u is absolutely convergent in H 1 (). Therefore the two series below u(, ω) γ(, ω) = k=0 + i ω k k=0 ω k ( σ j v k j ϵ j ( σ j h k j ) h k j + ϵ j ) v k j are absolutely convergent in H 1/2 (). The formulas (20) and (21) now follow by differentiation in ω. The proof of Theorem 1.1 requires applications of dk dω k Λ γ(,ω) ω=0 to traces of the complex (valued) geometrical optics solutions. However, if the Dirichlet data in (13) is complex valued, the formulas (20) and (21) no longer hold for v k and h k obtained by the recurrence (14) and (15). In order to preserve this compatibility we extend the real and imaginary part of d k dω k Λ γ(,ω) ω=0 by complex linearity: For k 0 and real valued f, g H 1/2 () we define R (k), I (k) : H 1/2 () H 1/2 () by (22) R (k) (f + ig) := 1 { } d k k! Re Λ γ(,ω) dω k (f) ω=0 (23) I (k) (f + ig) := 1 { } d k k! Im Λ γ(,ω) dω k (f) ω=0 + i k! Re { d k Λ γ(,ω) dω k + i k! Im { d k Λ γ(,ω) dω k } (g), ω=0 } (g). ω=0 Proposition 2.1. Let R (k), I (k) : H 1/2 () H 1/2 () be as defined in (22) and (23). Then ( ) R (k) v k j (f + ig) = σ j ϵ h k j (24) j (25) I (k) (f + ig) = where h 0 0, v 0 solves σ 0 v 0 = 0, ( σ j h k j + ϵ j v 0 = f + ig ) v k j, and h k and v k are defined inductively by the recurrence (14) and (15) for k 1.
8 8 SUNGWHAN KIM AND ALEXANDRU TAMASAN Proof. The identities hold when g 0 by Corollary 2.1. The left hand sides of (24) and (25) are complex linear by definition. But the right hand side is also complex linear since the sequences{h k } and {v k } produced by the recurrence (14) and (15) are complex linear with the Dirichlet data. 3. PROOF OF THEOREM 1.1 AND ITS COROLLARY In this section we assume the hypotheses in Theorem 1.1. The proof is inductive and relies on the basic identities in the lemma below which relate the quantities (26) (27) Q r N := ( σ N σ N ln σ 0 ) σ 0, Q i N := ( ϵ N ϵ N ln σ 0 ) σ 0, for N 1, with the boundary operator R (N) in (22), respectively I (N) in (23), and some interior data previously determined in the inductive step. In the following lemma σ j, ϵ j for j = 1,..., N, need not be compactly supported in, and σ 0 need not be constant near the boundary, but rather its normal derivative be vanishing at the boundary. Lemma 3.1. Let w 1, w 2 H 1 () be σ 0 -harmonic. Then ( σ 1 σ 1 ln σ 0 )w 1 w 2 dx = 2 R (1) [w 1 ]w 2 ds (w 1 w 2 ) (σ 1 σ 1 (28) w 1w 2 )ds, and, for N 2, ( σ N σ N ln σ 0 )w 1 w 2 dx = 2 R (N) (w 1 w 2 ) [w 1 ]w 2 ds (σ N σ N w 1w 2 )ds N 1 (29) 2 (σ j v N j ϵ j h N j ) w 2 dx. Also (30) ( ϵ 1 ϵ 1 ln σ 0 )w 1 w 2 dx = 2 I (1) [w 1 ]w 2 ds (w 1 w 2 ) (ϵ 1 ϵ 1 w 1w 2 )ds,
9 A CALDERÓN PROBLEM FOR DISPERSIVE MEDIA 9 and, for N 2, ( ϵ N ϵ N ln σ 0 )w 1 w 2 dx = 2 I (N) [w 1 ]w 2 ds ( (w 1 w 2 ) ϵ N ϵ ) N w 1w 2 ds N 1 (31) 2 (σ j h N j + ϵ j v N j ) w 2 dx. Proof. From Proposition 2.1, (32) R (N) [w 1 ] = N ( σ j v N j ϵ j ) h N j where v 0 = w 1, h 0 = ϵ 0 = 0, and v j, h j, j = 1,..., N are defined recursively in (14) and (15). Upon multiplication of (32) by w 2 and integration by parts we obtain (33) N ( R (N) v N j [w 1 ]w 2 ds = σ j ϵ h N j j ( N ) = w 2 (σ j v N j ϵ j h N j ) dx + = = N + N 1 (σ j v N j ϵ j h N j ) w 2 dx σ 0 v N w 2 dx + σ N w 1 w 2 dx N 1 ) w 2 ds (σ j v N j ϵ j h N j ) w 2 dx (σ j v N j dx ϵ j h N j ) w 2 dx + σ N w 1 w 2 dx, where in the third equality we use the recurrence (14), and for the last equality we use the facts that w 2 is σ 0 -harmonic and v N H0(). 1 Note that for N = 1 the identity (33) reduces to (34) R (1) [w 1 ]w 2 ds = σ 1 w 1 w 2 dx.
10 10 SUNGWHAN KIM AND ALEXANDRU TAMASAN (35) Since σ 0 C 1,1 (), w j also solves w j = ln σ 0 w j, j = 1, 2. By using the general identity 2 w 1 w 2 = (w 1 w 2 ) w 1 w 2 w 2 w 1, and (35) we get 2 σ N w 1 w 2 dx = σ N ( (w 1 w 2 ) w 1 w 2 w 2 w 1 )dx = σ N ( (w 1 w 2 ) + ln σ 0 (w 1 w 2 ))dx ( (w 1 w 2 ) = σ N σ ) N w 1w 2 ds + ( σ N )w 1 w 2 dx + σ N ln σ 0 (w 1 w 2 )dx ( (w 1 w 2 ) = σ N σ ) N w ln σ 0 1w 2 ds + σ N w 1w 2 ds + ( σ N σ N ln σ 0 )w 1 w 2 dx ( (w 1 w 2 ) = σ N σ ) N w 1w 2 ds (36) + ( σ N σ N ln σ 0 )w 1 w 2 dx, where, for the last equality, the vanishing of the second boundary integral accounts for σ 0 being constant near the boundary. The use of (34), respectively (36), in (33) proves the identity (28), respectively (29). The identity (31) follows similarly: An integration by parts and use of (15) yield (37) I (N) [w 1 ]w 2 ds = N 1 (σ j h N j + ϵ j v N j ) w 2 dx + and the calculation in (36) with ϵ N replacing σ N shows ( (w 1 w 2 ) 2 ϵ N w 1 w 2 dx = ϵ N ϵ ) N w 1w 2 ds (38) + ( ϵ N ϵ N ln σ 0 )w 1 w 2 dx,. By using (38) in (37) we obtain (31). ϵ N w 1 w 2 dx,
11 A CALDERÓN PROBLEM FOR DISPERSIVE MEDIA 11 We prove Theorem 1.1 by induction on the order of derivative in frequency. We will employ the basic identities in the lemma above with w 1, w 2 being the CGO solutions in (8). Recall the dependence of w 1 and w 2 on the vectors η, k, l, which we do not made explicit to keep the formulas readable. Case N = 0: Assume that Λ γ(,0) is known, then we appeal to [16, Lemmae 2.7 and 2.12 b)] to conclude that Λ γ(,0) recovers the traces w j for j = 1, 2, and to [16, Theorem 5.1] to recover coefficient σ 0 in. Case N = 1: Assume that R (1) and I (1) are known. By using the CGO solutions w 1, w 2 in (30), and accounting for the compact support of σ 1, ϵ 1, we get (39) e ix k Q r 1(x)(1 + ψ 1 (x))(1 + ψ 2 (x))dx = R n R (1) [w 1 ]w 2 ds. Due to the decay estimate (9), the left hand side in (39) has a limit with l, thus so does the right hand side, and the Fourier transform F[Q r 1] is determined. The Fourier inversion yields [ Q r 1 = F 1 lim l ] R (1) [w 1 ]w 2 ds, Starting with (30) and proceeding similarly with the boundary data I (1), we find by Fourier inversion the L ()-function [ Q i 1 = F 1 lim l ] I (1) [w 1 ]w 2 ds. d In the case N = 1 we note that dω ω=0 Λ γ(,ω) uniquely determines the quantities Q r 1 and Q i 1, independent of knowledge of the DtN map Λ γ(,0). However, since σ 0 is known, we can also recover σ 1 and ϵ 1 as the unique solutions of the respective Dirichlet problems: (40) (41) σ 1 σ 1 ln σ 0 σ 1 ln σ 0 = σ 0 Q r 1 in, σ 1 = 0, ϵ 1 ϵ 1 ln σ 0 ϵ 1 ln σ 0 = σ 0 Q i 1 in, ϵ 1 = 0. The inductive step: Assume that σ 0, and σ j, ϵ j, j = 1,..., N 1 have been determined, and R (N), I (N) are known. Below we reconstruct σ N and ϵ N. By solving the recurrence (14) and (15) starting with v 0 = w 1, we can determine v j, h j H0() 1 for j = 1,..., N 1 (recall v 0 = w 1 and h 0 = 0).
12 12 SUNGWHAN KIM AND ALEXANDRU TAMASAN Since σ N is assumed compactly supported in, the identity (29) yields e ix k Q r N(x)(1 + ψ 1 (x))(1 + ψ 2 (x))dx R n N 1 (42) = 2 R (N) [w 1 ]w 2 ds 2 (σ j v N j ϵ j h N j ) w 2 dx. Note that the right hand side of (42) is known. Due to the decay estimate (9) the left hand side of (42) has a limit with l, thus so does the right hand side, and the Fourier transform F[Q r N ] is recovered from { F [Q r N] (k) = lim R (N) [w 1 ]w 2 ds l (43) 2 N 1 (σ j v N j ϵ j h N j ) w 2 dx By Fourier inversion in (42), we recovered Q r N L 0 (). By a similar reasoning starting from (31) with w 1, w 2 being the CGO solutions in (8), we determine Q i N L 0 (). Finally, the coefficients σ N and ϵ N, are recovered as the unique solution to the respective Dirichlet problems: (44) (45) σ N σ N ln σ 0 σ N ln σ 0 = σ 0 Q r N in, σ N = 0, ϵ N ϵ N ln σ 0 ϵ N ln σ 0 = σ 0 Q i N in, ϵ N = 0. This finishes the proof of Theorem 1.1. For the proof of the corollary, we note that the analyticity of the map ω Λ γ(,ω) allows for each derivative dk dω k Λ γ(,ω) at ω = 0 to be determined from the derivatives dj dω j Λ γ(,ω) ω0, j 0, taken at any other frequency ω 0 within the domain of analyticity [0, (1 + 2r) 1 ). 4. ACKNOWLEDGMENT The authors would like to thank the anonymous referees for the constructive criticism, which uncover a gap in a previous version. REFERENCES [1] K. Astala and L. Päivärinta, Calderón s inverse conductivity problem in the plane, Ann. of Math. 163(2006), [2] R.H. Bayford, Bioimpedance Tomography (Electrical impedance tomography), Annu. Rev. Biomed. Eng., 8(2006), [3] R. M. Brown and G. Uhlmann, Uniqueness in the inverse conductivity problem for nonsmooth conductivities in two dimensions. Comm. Partial Differential Equations 22 (1997), no. 5-6, }.
13 A CALDERÓN PROBLEM FOR DISPERSIVE MEDIA 13 [4] L. Borcea, EIT Electrical impedance tomography, Inverse Problems, 18(2002), R99 R136. [5] A.L. Bukhgeim, Recovering a potential from Cauchy data in the two-dimensional case, J. Inv. Ill-posed Problems 16(2008), [6] A.P. Calderón, On an inverse boundary value problem, in Seminar on Numerical Analysis and its Applications to Continuum Physics, Soc. Brasiliera de Matematica, Rio de Janeiro, 1980 [7] M. Cheney, D. Isaacson, and J. C. Newell, Electrical impedance tomography, SIAM Review, 41(1999), [8] K.R. Foster and H.P. Schwan, Dielectric properties of tissues and biological materials, a critical review, Crit. Rev. Biomed. Eng. 17 (1989), no. 1, [9] E. Francini, Recovering a complex coefficient in a planar domain from the Dirichletto-Neumann map, Inverse Problems , doi: / /16/1/309 [10] D. Holder, Electrical Impedance Tomography: Methods, History and Applications (Bristol, UK: IOP Publishing), [11] S. Kim, J. Lee, J. K. Seo, E. J. Woo and H. Zribi, Multifrequency trans-admittance scanner: mathematical framework and feasibility, SIAM J. Appl. Math., 69(2008), no. 1, [12] S. Kim, J.K. Seo, and T. Ha, A nondestructive evaluation method for concrete voids: frequency differential electrical impedance scanning,siam J. Appl. Math., 69(2009), no. 6, [13] S. Kim, E.J. Lee, E.J. Woo, and J.K. Seo, Asymptotic analysis of the membrane structure to sensitivity of frequency-difference electrical impedance tomography, Inverse problems, 28(2012), (17pp). [14] K. Knudsen and A. Tamasan, Reconstruction of less regular conductivities in the plane, Comm. Partial Differential Equations 29 (2004), no. 3-4, [15] R.V. Kohn and M. Vogelius, Determining conductivity by boundary measurements, Commun. Pure Appl. Math., 37(1984), [16] A. Nachman, Reconstructions from boundary measurements, Ann. Math., 128(1988), [17] A. Nachman, Global uniqueness for a two-dimensional inverse boundary value problem, Ann. Math. 143(1996), [18] J.K. Seo, J. Lee, S.W. Kim, H. Zribi, and E.J. Woo, Frequency-difference electrical impedance tomography (fdeit): algorithm development and feasibility study, Inverse Problems, 29(2008), [19] J.K. Seo, B. Harrach, and E.J. Woo, Recent progress on frequency difference electrical impedance tomography,esaim: PROCEEDINGS, 26(2009), [20] E. Somersalo, M. Cheney, and D. Isaacson, Existence and Uniqueness for an Electrode Model for ECCT, SIAM J. App. Math. 52(1992), [21] K. Sungwhan and A. Tamasan, On a Calderón problem in frequency differential electrical impedance tomography, SIAM J. Math. Anal. 45 (2013), no. 5, [22] J. Sylvester and G. Uhlmann, A global uniqueness theorem for an inverse boundary value problem, Ann. Math., 125(1987), [23] G. Uhlmann, Electrical impedance tomography and Calderón s problem, Inverse Problems, 25(2009), 1-39.
14 14 SUNGWHAN KIM AND ALEXANDRU TAMASAN DIVISION OF LIBERAL ARTS, HANBAT NATIONAL UNIVERSITY, KOREA address: DEPARTMENT OF MATHEMATICS, UNIVERSITY OF CENTRAL FLORIDA, ORLANDO, FLORIDA address:
is measured (n is the unit outward normal) at the boundary, to define the Dirichlet-
SIAM J MATH ANAL Vol 45, No 5, pp 700 709 c 013 Society for Industrial and Applied Mathematics ON A CALDERÓN PROBLEM IN FREQUENCY DIFFERENTIAL ELECTRICAL IMPEDANCE TOMOGRAPHY SUNGWHAN KIM AND ALEXANDRU
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 informationReconstructing conductivities with boundary corrected D-bar method
Reconstructing conductivities with boundary corrected D-bar method Janne Tamminen June 24, 2011 Short introduction to EIT The Boundary correction procedure The D-bar method Simulation of measurement data,
More informationOn uniqueness in the inverse conductivity problem with local data
On uniqueness in the inverse conductivity problem with local data Victor Isakov June 21, 2006 1 Introduction The inverse condictivity problem with many boundary measurements consists of recovery of conductivity
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 informationAALBORG UNIVERSITY. A new direct method for reconstructing isotropic conductivities in the plane. Department of Mathematical Sciences.
AALBORG UNIVERSITY A new direct method for reconstructing isotropic conductivities in the plane by Kim Knudsen January 23 R-23-1 Department of Mathematical Sciences Aalborg University Fredrik Bajers Vej
More informationAALBORG UNIVERSITY. The Calderón problem with partial data for less smooth conductivities. Kim Knudsen. Department of Mathematical Sciences
AALBORG UNIVERSITY The Calderón problem with partial data for less smooth conductivities by Kim Knudsen R-2005-07 February 2005 Department of Mathematical Sciences Aalborg University Fredrik Bajers Vej
More informationJUSTIFICATION OF POINT ELECTRODE MODELS IN ELECTRICAL IMPEDANCE TOMOGRAPHY
JUSTIFICATION OF POINT ELECTRODE MODELS IN ELECTRICAL IMPEDANCE TOMOGRAPHY MARTIN HANKE, BASTIAN HARRACH, AND NUUTTI HYVÖNEN Abstract. The most accurate model for real-life electrical impedance tomography
More informationInverse Electromagnetic Problems
Title: Name: Affil./Addr. 1: Inverse Electromagnetic Problems Gunther Uhlmann, Ting Zhou University of Washington and University of California Irvine /gunther@math.washington.edu Affil./Addr. 2: 340 Rowland
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 informationCalderón s inverse problem in 2D Electrical Impedance Tomography
Calderón s inverse problem in 2D Electrical Impedance Tomography Kari Astala (University of Helsinki) Joint work with: Matti Lassas, Lassi Päivärinta, Samuli Siltanen, Jennifer Mueller and Alan Perämäki
More informationStability and instability in inverse problems
Stability and instability in inverse problems Mikhail I. Isaev supervisor: Roman G. Novikov Centre de Mathématiques Appliquées, École Polytechnique November 27, 2013. Plan of the presentation The Gel fand
More informationCOMPLEX SPHERICAL WAVES AND INVERSE PROBLEMS IN UNBOUNDED DOMAINS
COMPLEX SPHERICAL WAVES AND INVERSE PROBLEMS IN UNBOUNDED DOMAINS MIKKO SALO AND JENN-NAN WANG Abstract. This work is motivated by the inverse conductivity problem of identifying an embedded object in
More informationPut Paper Number Here
Proceedings of Proceedings of icipe : rd International Conference on Inverse Problems in Engineering June,, Port Ludlow, WA, USA Put Paper Number Here NEW RECONSTRUCTION ALGORITHMS IN ELECTRICAL IMPEDANCE
More informationMapping Properties of the Nonlinear Fourier Transform in Dimension Two
Communications in Partial Differential Equations, 3: 59 60, 007 Copyright Taylor & Francis Group, LLC ISSN 0360-530 print/53-433 online DOI: 0.080/036053005005304 Mapping Properties of the Nonlinear Fourier
More informationMonotonicity arguments in electrical impedance tomography
Monotonicity arguments in electrical impedance tomography Bastian Gebauer gebauer@math.uni-mainz.de Institut für Mathematik, Joh. Gutenberg-Universität Mainz, Germany NAM-Kolloquium, Georg-August-Universität
More informationarxiv: v1 [math.ap] 21 Dec 2014
CONSTRUCTION OF INVISIBLE CONUCTIVITY PERTURBATIONS FOR THE POINT ELECTROE MOEL IN ELECTRICAL IMPEANCE TOMOGRAPHY LUCAS CHESNEL, NUUTTI HYVÖNEN, AN STRATOS STABOULIS Abstract. We explain how to build invisible
More informationNon-uniqueness result for a hybrid inverse problem
Non-uniqueness result for a hybrid inverse problem Guillaume Bal and Kui Ren Abstract. Hybrid inverse problems aim to combine two imaging modalities, one that displays large contrast and one that gives
More informationDETERMINING A FIRST ORDER PERTURBATION OF THE BIHARMONIC OPERATOR BY PARTIAL BOUNDARY MEASUREMENTS
DETERMINING A FIRST ORDER PERTURBATION OF THE BIHARMONIC OPERATOR BY PARTIAL BOUNDARY MEASUREMENTS KATSIARYNA KRUPCHYK, MATTI LASSAS, AND GUNTHER UHLMANN Abstract. We consider an operator 2 + A(x) D +
More informationComputational inversion based on complex geometrical optics solutions
Computational inversion based on complex geometrical optics solutions Samuli Siltanen Department of Mathematics and Statistics University of Helsinki, Finland samuli.siltanen@helsinki.fi http://www.siltanen-research.net
More informationLOCALIZED POTENTIALS IN ELECTRICAL IMPEDANCE TOMOGRAPHY. Bastian Gebauer. (Communicated by Otmar Scherzer)
Inverse Problems and Imaging Volume 2, No. 2, 2008, 251 269 Web site: http://www.aimsciences.org LOCALIZED POTENTIALS IN ELECTRICAL IMPEDANCE TOMOGRAPHY astian Gebauer Institut für Mathematik Johannes
More informationThe elliptic sinh-gordon equation in the half plane
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 8 25), 63 73 Research Article The elliptic sinh-gordon equation in the half plane Guenbo Hwang Department of Mathematics, Daegu University, Gyeongsan
More informationInverse parameter identification problems
Inverse parameter identification problems Bastian von Harrach harrach@math.uni-stuttgart.de Chair of Optimization and Inverse Problems, University of Stuttgart, Germany ICP - Institute for Computational
More informationNumerical Harmonic Analysis on the Hyperbolic Plane
Numerical Harmonic Analysis on the Hyperbolic Plane Buma Fridman, Peter Kuchment, Kirk Lancaster, Serguei Lissianoi, Mila Mogilevsky, Daowei Ma, Igor Ponomarev, and Vassilis Papanicolaou Mathematics and
More informationDetecting stochastic inclusions in electrical impedance tomography
Detecting stochastic inclusions in electrical impedance tomography Bastian von Harrach harrach@math.uni-frankfurt.de (joint work with A. Barth, N. Hyvönen and L. Mustonen) Institute of Mathematics, Goethe
More informationA SHARP STABILITY ESTIMATE IN TENSOR TOMOGRAPHY
A SHARP STABILITY ESTIMATE IN TENSOR TOMOGRAPHY PLAMEN STEFANOV 1. Introduction Let (M, g) be a compact Riemannian manifold with boundary. The geodesic ray transform I of symmetric 2-tensor fields f is
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 informationINVERSE BOUNDARY VALUE PROBLEMS FOR THE PERTURBED POLYHARMONIC OPERATOR
INVERSE BOUNDARY VALUE PROBLEMS FOR THE PERTURBED POLYHARMONIC OPERATOR KATSIARYNA KRUPCHYK, MATTI LASSAS, AND GUNTHER UHLMANN Abstract. We show that a first order perturbation A(x) D + q(x) of the polyharmonic
More informationNon-uniqueness result for a hybrid inverse problem
Non-uniqueness result for a hybrid inverse problem Guillaume Bal and Kui Ren Abstract. Hybrid inverse problems aim to combine two imaging modalities, one that displays large contrast and one that gives
More informationReconstructing conductivities with boundary corrected D-bar method
J. Inverse Ill-Posed Probl., Ahead of Print DOI 10.1515/jip-2013-0042 de Gruyter 2014 Reconstructing conductivities with boundary corrected D-bar method Samuli Siltanen and Janne P. Tamminen Abstract.
More informationDIRECT ELECTRICAL IMPEDANCE TOMOGRAPHY FOR NONSMOOTH CONDUCTIVITIES
DIRECT ELECTRICAL IMPEDANCE TOMOGRAPHY FOR NONSMOOTH CONDUCTIVITIES K ASTALA, J L MUELLER, A PERÄMÄKI, L PÄIVÄRINTA AND S SILTANEN Abstract. A new reconstruction algorithm is presented for eit in dimension
More informationarxiv: v1 [math.ap] 11 Jun 2007
Inverse Conductivity Problem for a Parabolic Equation using a Carleman Estimate with One Observation arxiv:0706.1422v1 [math.ap 11 Jun 2007 November 15, 2018 Patricia Gaitan Laboratoire d Analyse, Topologie,
More informationHyperbolic inverse problems and exact controllability
Hyperbolic inverse problems and exact controllability Lauri Oksanen University College London An inverse initial source problem Let M R n be a compact domain with smooth strictly convex boundary, and let
More informationSensitivity Analysis of 3D Magnetic Induction Tomography (MIT)
Sensitivity Analysis of 3D Magnetic Induction Tomography (MIT) W R B Lionheart 1, M Soleimani 1, A J Peyton 2 1 Department of Mathematics UMIST, Manchester, UK, Email: bill.lionheart@umist.ac.uk, 2 Department
More informationDetermination of thin elastic inclusions from boundary measurements.
Determination of thin elastic inclusions from boundary measurements. Elena Beretta in collaboration with E. Francini, S. Vessella, E. Kim and J. Lee September 7, 2010 E. Beretta (Università di Roma La
More 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 informationA RELATIONSHIP BETWEEN THE DIRICHLET AND REGULARITY PROBLEMS FOR ELLIPTIC EQUATIONS. Zhongwei Shen
A RELATIONSHIP BETWEEN THE DIRICHLET AND REGULARITY PROBLEMS FOR ELLIPTIC EQUATIONS Zhongwei Shen Abstract. Let L = diva be a real, symmetric second order elliptic operator with bounded measurable coefficients.
More informationNovel tomography techniques and parameter identification problems
Novel tomography techniques and parameter identification problems Bastian von Harrach harrach@ma.tum.de Department of Mathematics - M1, Technische Universität München Colloquium of the Institute of Biomathematics
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 informationLiouville-type theorem for the Lamé system with singular coefficients
Liouville-type theorem for the Lamé system with singular coefficients Blair Davey Ching-Lung Lin Jenn-Nan Wang Abstract In this paper, we study a Liouville-type theorem for the Lamé system with rough coefficients
More informationContributors and Resources
for Introduction to Numerical Methods for d-bar Problems Jennifer Mueller, presented by Andreas Stahel Colorado State University Bern University of Applied Sciences, Switzerland Lexington, Kentucky, May
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 informationFast shape-reconstruction in electrical impedance tomography
Fast shape-reconstruction in electrical impedance tomography Bastian von Harrach bastian.harrach@uni-wuerzburg.de (joint work with Marcel Ullrich) Institut für Mathematik - IX, Universität Würzburg The
More informationIntegral Representation Formula, Boundary Integral Operators and Calderón projection
Integral Representation Formula, Boundary Integral Operators and Calderón projection Seminar BEM on Wave Scattering Franziska Weber ETH Zürich October 22, 2010 Outline Integral Representation Formula Newton
More informationNONLINEAR FOURIER ANALYSIS FOR DISCONTINUOUS CONDUCTIVITIES: COMPUTATIONAL RESULTS
NONLINEAR FOURIER ANALYSIS FOR DISCONTINUOUS CONDUCTIVITIES: COMPUTATIONAL RESULTS K. ASTALA, L. PÄIVÄRINTA, J. M. REYES AND S. SILTANEN Contents. Introduction. The nonlinear Fourier transform 6 3. Inverse
More informationNumerical computation of complex geometrical optics solutions to the conductivity equation
Numerical computation of complex geometrical optics solutions to the conductivity equation Kari Astala a, Jennifer L Mueller b, Lassi Päivärinta a, Samuli Siltanen a a Department of Mathematics and Statistics,
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 informationA local estimate from Radon transform and stability of Inverse EIT with partial data
A local estimate from Radon transform and stability of Inverse EIT with partial data Alberto Ruiz Universidad Autónoma de Madrid U. California, Irvine.June 2012 w/ P. Caro (U. Helsinki) and D. Dos Santos
More informationAllan Greenleaf, Matti Lassas, and Gunther Uhlmann
Mathematical Research Letters 10, 685 693 (2003) ON NONUNIQUENESS FOR CALDERÓN S INVERSE PROBLEM Allan Greenleaf, Matti Lassas, and Gunther Uhlmann Abstract. We construct anisotropic conductivities with
More informationIncreasing stability in an inverse problem for the acoustic equation
Increasing stability in an inverse problem for the acoustic equation Sei Nagayasu Gunther Uhlmann Jenn-Nan Wang Abstract In this work we study the inverse boundary value problem of determining the refractive
More informationINVERSE BOUNDARY VALUE PROBLEMS FOR THE MAGNETIC SCHRÖDINGER EQUATION
INVERSE BOUNDARY VALUE PROBLEMS FOR THE MAGNETIC SCHRÖDINGER EQUATION MIKKO SALO Abstract. We survey recent results on inverse boundary value problems for the magnetic Schrödinger equation. 1. Introduction
More informationSTOKES PROBLEM WITH SEVERAL TYPES OF BOUNDARY CONDITIONS IN AN EXTERIOR DOMAIN
Electronic Journal of Differential Equations, Vol. 2013 2013, No. 196, pp. 1 28. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu STOKES PROBLEM
More informationA mathematical model and inversion procedure for Magneto-Acousto-Electric Tomography (MAET)
A mathematical model and inversion procedure for Magneto-Acousto-Electric Tomography (MAET) Leonid Kunyansky University of Arizona, Tucson, AZ Suppported in part by NSF grant DMS-0908243 NSF grant "we
More informationTraces, extensions and co-normal derivatives for elliptic systems on Lipschitz domains
Traces, extensions and co-normal derivatives for elliptic systems on Lipschitz domains Sergey E. Mikhailov Brunel University West London, Department of Mathematics, Uxbridge, UB8 3PH, UK J. Math. Analysis
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 informationOn stable inversion of the attenuated Radon transform with half data Jan Boman. We shall consider weighted Radon transforms of the form
On stable inversion of the attenuated Radon transform with half data Jan Boman We shall consider weighted Radon transforms of the form R ρ f(l) = f(x)ρ(x, L)ds, L where ρ is a given smooth, positive weight
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 informationA new inversion method for dissipating electromagnetic problem
A new inversion method for dissipating electromagnetic problem Elena Cherkaev Abstract The paper suggests a new technique for solution of the inverse problem for Maxwell s equations in a dissipating medium.
More informationFrom Electrical Impedance Tomography to Network Tomography
From Electrical Impedance Tomography to Network Tomography Carlos A. Berenstein Institute for Systems Research ISR University of Maryland College Park With the collaboration of Prof. John Baras and Franklin
More informationThe oblique derivative problem for general elliptic systems in Lipschitz domains
M. MITREA The oblique derivative problem for general elliptic systems in Lipschitz domains Let M be a smooth, oriented, connected, compact, boundaryless manifold of real dimension m, and let T M and T
More informationExercises - Chapter 1 - Chapter 2 (Correction)
Université de Nice Sophia-Antipolis Master MathMods - Finite Elements - 28/29 Exercises - Chapter 1 - Chapter 2 Correction) Exercise 1. a) Let I =], l[, l R. Show that Cl) >, u C Ī) Cl) u H 1 I), u DĪ).
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 informationSobolev regularity for the Monge-Ampère equation, with application to the semigeostrophic equations
Sobolev regularity for the Monge-Ampère equation, with application to the semigeostrophic equations Alessio Figalli Abstract In this note we review some recent results on the Sobolev regularity of solutions
More informationDownloaded 11/20/13 to Redistribution subject to SIAM license or copyright; see
SIAM J. SCI. COMPUT. Vol. 33, No. 1, pp. 328 341 c 2011 Society for Industrial and Applied Mathematics A SIMPLE NUMERICAL METHOD FOR COMPLEX GEOMETRICAL OPTICS SOLUTIONS TO THE CONDUCTIVITY EQUATION KUI
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 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 informationPotential Analysis meets Geometric Measure Theory
Potential Analysis meets Geometric Measure Theory T. Toro Abstract A central question in Potential Theory is the extend to which the geometry of a domain influences the boundary regularity of the solution
More informationReconstructing inclusions from Electrostatic Data
Reconstructing inclusions from Electrostatic Data Isaac Harris Texas A&M University, Department of Mathematics College Station, Texas 77843-3368 iharris@math.tamu.edu Joint work with: W. Rundell Purdue
More informationElectrical Impedance Tomography: 3D Reconstructions using Scattering Transforms
Downloaded from orbit.dtu.dk on: Nov 11, 2018 Electrical Impedance Tomography: 3D Reconstructions using Scattering Transforms Delbary, Fabrice; Hansen, Per Christian; Knudsen, Kim Published in: Applicable
More informationAN ELECTROMAGNETIC INVERSE PROBLEM IN CHIRAL MEDIA
AN ELECTROMAGNETIC INVERSE PROBLEM IN CHIRAL MEDIA STEPHEN R. MCDOWALL Abstract. We consider the inverse boundary value problem for Maxwell s equations that takes into account the chirality of a body in
More informationUniqueness in determining refractive indices by formally determined far-field data
Applicable Analysis, 2015 Vol. 94, No. 6, 1259 1269, http://dx.doi.org/10.1080/00036811.2014.924215 Uniqueness in determining refractive indices by formally determined far-field data Guanghui Hu a, Jingzhi
More informationRecent progress on the factorization method for electrical impedance tomography
Recent progress on the factorization method for electrical impedance tomography Bastian von Harrach harrach@math.uni-stuttgart.de Chair of Optimization and Inverse Problems, University of Stuttgart, Germany
More informationMath The Laplacian. 1 Green s Identities, Fundamental Solution
Math. 209 The Laplacian Green s Identities, Fundamental Solution Let be a bounded open set in R n, n 2, with smooth boundary. The fact that the boundary is smooth means that at each point x the external
More informationarxiv: v2 [math.ap] 17 Jun 2014
NONLINEAR FOURIER ANALYSIS FOR DISCONTINUOUS CONDUCTIVITIES: COMPUTATIONAL RESULTS K. ASTALA, L. PÄIVÄRINTA, J. M. REYES AND S. SILTANEN arxiv:30.705v [math.ap] 7 Jun 04 Contents. Introduction. The nonlinear
More informationu( x) = g( y) ds y ( 1 ) U solves u = 0 in U; u = 0 on U. ( 3)
M ath 5 2 7 Fall 2 0 0 9 L ecture 4 ( S ep. 6, 2 0 0 9 ) Properties and Estimates of Laplace s and Poisson s Equations In our last lecture we derived the formulas for the solutions of Poisson s equation
More informationANGKANA RÜLAND AND MIKKO SALO
QUANTITATIVE RUNGE APPROXIMATION AND INVERSE PROBLEMS ANGKANA RÜLAND AND MIKKO SALO Abstract. In this short note we provide a quantitative version of the classical Runge approximation property for second
More informationShort note on compact operators - Monday 24 th March, Sylvester Eriksson-Bique
Short note on compact operators - Monday 24 th March, 2014 Sylvester Eriksson-Bique 1 Introduction In this note I will give a short outline about the structure theory of compact operators. I restrict attention
More informationarxiv: v1 [math.ap] 6 Mar 2017
A Fidelity-embedded Regularization Method for Robust Electrical Impedance Tomography Kyounghun Lee, Eung Je Woo, and Jin Keun Seo arxiv:1703.01821v1 [math.ap] 6 Mar 2017 March 7, 2017 Abstract Electrical
More informationInverse problems and medical imaging
Inverse problems and medical imaging Bastian von Harrach harrach@math.uni-frankfurt.de Institute of Mathematics, Goethe University Frankfurt, Germany Seminario di Calcolo delle Variazioni ed Equazioni
More informationA class of domains with fractal boundaries: Functions spaces and numerical methods
A class of domains with fractal boundaries: Functions spaces and numerical methods Yves Achdou joint work with T. Deheuvels and N. Tchou Laboratoire J-L Lions, Université Paris Diderot École Centrale -
More informationModelling in photonic crystal structures
Modelling in photonic crystal structures Kersten Schmidt MATHEON Nachwuchsgruppe Multiscale Modelling and Scientific Computing with PDEs in collaboration with Dirk Klindworth (MATHEON, TU Berlin) Holger
More informationS chauder Theory. x 2. = log( x 1 + x 2 ) + 1 ( x 1 + x 2 ) 2. ( 5) x 1 + x 2 x 1 + x 2. 2 = 2 x 1. x 1 x 2. 1 x 1.
Sep. 1 9 Intuitively, the solution u to the Poisson equation S chauder Theory u = f 1 should have better regularity than the right hand side f. In particular one expects u to be twice more differentiable
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 informationEXISTENCE OF SOLUTIONS TO ASYMPTOTICALLY PERIODIC SCHRÖDINGER EQUATIONS
Electronic Journal of Differential Equations, Vol. 017 (017), No. 15, pp. 1 7. ISSN: 107-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE OF SOLUTIONS TO ASYMPTOTICALLY PERIODIC
More informationGAUSSIAN MEASURE OF SECTIONS OF DILATES AND TRANSLATIONS OF CONVEX BODIES. 2π) n
GAUSSIAN MEASURE OF SECTIONS OF DILATES AND TRANSLATIONS OF CONVEX BODIES. A. ZVAVITCH Abstract. In this paper we give a solution for the Gaussian version of the Busemann-Petty problem with additional
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 informationOn Friedrichs inequality, Helmholtz decomposition, vector potentials, and the div-curl lemma. Ben Schweizer 1
On Friedrichs inequality, Helmholtz decomposition, vector potentials, and the div-curl lemma Ben Schweizer 1 January 16, 2017 Abstract: We study connections between four different types of results that
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 informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationESTIMATES FOR ELLIPTIC HOMOGENIZATION PROBLEMS IN NONSMOOTH DOMAINS. Zhongwei Shen
W,p ESTIMATES FOR ELLIPTIC HOMOGENIZATION PROBLEMS IN NONSMOOTH DOMAINS Zhongwei Shen Abstract. Let L = div`a` x, > be a family of second order elliptic operators with real, symmetric coefficients on a
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 informationBOUNDARY VALUE PROBLEMS ON A HALF SIERPINSKI GASKET
BOUNDARY VALUE PROBLEMS ON A HALF SIERPINSKI GASKET WEILIN LI AND ROBERT S. STRICHARTZ Abstract. We study boundary value problems for the Laplacian on a domain Ω consisting of the left half of the Sierpinski
More informationFuzzy Statistical Limits
Fuzzy Statistical Limits Mark Burgin a and Oktay Duman b a Department of Mathematics, University of California, Los Angeles, California 90095-1555, USA b TOBB Economics and Technology University, Faculty
More informationThe mathematics of scattering by unbounded, rough, inhomogeneous layers
The mathematics of scattering by unbounded, rough, inhomogeneous layers Simon N. Chandler-Wilde a Peter Monk b Martin Thomas a a Department of Mathematics, University of Reading, Whiteknights, PO Box 220
More informationRegularity for Poisson Equation
Regularity for Poisson Equation OcMountain Daylight Time. 4, 20 Intuitively, the solution u to the Poisson equation u= f () should have better regularity than the right hand side f. In particular one expects
More informationWavelets and regularization of the Cauchy problem for the Laplace equation
J. Math. Anal. Appl. 338 008440 1447 www.elsevier.com/locate/jmaa Wavelets and regularization of the Cauchy problem for the Laplace equation Chun-Yu Qiu, Chu-Li Fu School of Mathematics and Statistics,
More informationMixed exterior Laplace s problem
Mixed exterior Laplace s problem Chérif Amrouche, Florian Bonzom Laboratoire de mathématiques appliquées, CNRS UMR 5142, Université de Pau et des Pays de l Adour, IPRA, Avenue de l Université, 64000 Pau
More informationA NUMERICAL ITERATIVE SCHEME FOR COMPUTING FINITE ORDER RANK-ONE CONVEX ENVELOPES
A NUMERICAL ITERATIVE SCHEME FOR COMPUTING FINITE ORDER RANK-ONE CONVEX ENVELOPES XIN WANG, ZHIPING LI LMAM & SCHOOL OF MATHEMATICAL SCIENCES, PEKING UNIVERSITY, BEIJING 87, P.R.CHINA Abstract. It is known
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationRegularizing the optimization-based solution of inverse coefficient problems with monotonicity constraints
Regularizing the optimization-based solution of inverse coefficient problems with monotonicity constraints Bastian von Harrach (joint work with Mach Nguyet Minh) http://numerical.solutions Institute of
More information