arxiv: v3 [math.ap] 28 Mar 2018

Transcription

1 Electric field concentration in the presence of an inclusion with eccentric core-shell geometry Junbeom Kim Mikyoung Lim arxiv: v3 [math.ap] 28 Mar 218 May 3, 218 Abstract In this paper we analyze the gradient blow-up of the solution to the conductivity problem in two dimensions in the presence of an inclusion with eccentric core-shell geometry. Assuming that the core and shell have circular boundaries that are nearly touching, we derive an asymptotic formula for the solution in terms of the single and double layer potentials with image line charges. We also deduce an integral formula with image line charges for the problem relating to two nearly touching separated conductors. Mathematics subject classification MSC2): 35J25, 73C4 Key words. geometry, Image charge Contents Gradient blow-up, Conductivity equation, Bipolar coordinates, Anti-plane elasticity, Core-shell 1 Introduction 2 2 The bipolar coordinate system Geometry of the core-shell structure The bipolar coordinate system associated with the eccentric core-shell structure. 5 3 Decomposition of the solution in singular and regular terms Layer potential formulation Series expansions for the density functions in bipolar coordinates Blow-up contribution term in H and decomposition of the solution Asymptotics of the singular term in bipolar coordinates Series solution by the separation of variables Integral expression of the singular term in bipolar coordinates Optimal bounds This work is supported by the Korean Ministry of Science, ICT and Future Planning through NRF grant No. 216R1A2B41453 to M.L. and J.K.). Department of Mathematical Sciences, Korea Advanced Institute of Science and Technology, Daejeon 34141, Korea kjb2474@kaist.ac.kr, mklim@kaist.ac.kr). 1

2 5 Image line charge formula Line charge formula of the Lerch transcendent function Asymptotic formula for the core-shell geometry Asymptotic formula for the case of two separated disks Conclusion 27 1 Introduction We consider the conductivity problem in the presence of an inclusion with eccentric core-shell geometry which consists of a core and a surrounding shell, where the core and shell have circular boundaries that are nearly touching. Specifically, we let B i and B e be two disks such that B e contains B i, as drawn in Figure 2.1, where the core B i and shell B e \ B i have the constant conductivities k i and k e satisfying < k i, k e 1 <. The aim of the paper is to derive an asymptotic formula for the solution to the conductivity problem { σ u = in R 2, ux) Hx) = O x 1 1.1) ) as x, as the distance between B i and B e degenerates to zero. Here, σ denotes the conductivity profile σ = k i χb i ) + k e χb e \ B i ) + χr 2 \ B e ), H is an entire harmonic function, and χd) represents the characteristic function of a set D. The electrostatic interaction problem has been studied for a long time, including [37]. For two separated spheres, the series expression by separation of variables in bispherical coordinates and the method of image charges were developed [19, 2, 45]. The electric field may blow up as the distance tends toward zero [26], and it causes difficulty in the numerical computation. There have been many attempts to overcome this difficulty and to compute the asymptotic behavior of the electric field [11, 26, 38, 39, 4, 41]. For two dimensions, this problem can be also considered in relation to the stress concentration in composite materials. In composite structures, which consist of a matrix and inclusions, it is common for inclusions to be closely located. Having extreme material parameters is a necessary condition for inclusions to induce the electric field or the stress tensor blow-up. Li and Vogelius showed that the electric field is bounded if the inclusions conductivities are finite and strictly positive [31], and Li and Nirenberg extended this result to elliptic systems [32]. In recent years, the gradient blow-up feature of the electric field in the small area between the inclusions boundaries has attracted much attention due to its application to various imaging modalities, for example [28, 42]. In particular, the electric field concentration in the presence of two nearly touching separated inclusions has been intensively studied. We denote by ɛ the distance between the inclusions. The generic blow-up rate for conductors is ɛ 1/2 in two dimensions [9, 1, 12, 18, 26, 46, 47], and ɛ ln ɛ 1 in three dimensions [13, 14, 36]. In two dimensions, Ammari et al. deduced an optimal bound for two disks with arbitrary constant conductivities [7, 9], and Yun obtained the blow-up rate for perfect conductors of general shapes [46, 47]. The asymptotic behavior of the gradient blow-up solution was studied for perfect conductors of circular shapes by Kang et al. [22], and of arbitrary shapes by Ammari et al. [2] and Kang et al. [21]. For disks with arbitrary constant conductivities, Lim and Yu obtained an asymptotic formula for the electric field [34]. The blow-up feature for perfect conductors of a bow-tie 2

3 structure was investigated by Kang and Yun in [25]. In three dimensions, the blow-up rate for perfect conductors of general shapes was obtained by Bao et al. [13, 14]. For perfect conductors of spherical shapes, Ammari et al. showed that the potential difference in the conductors is necessary for the gradient blow-up [3], an optimal bound for the electric field concentration was obtained by Lim and Yun [36], and the asymptotic feature of the blow-up solution was investigated by Kang et al. [23] and by Lim and Yu [33]. Lekner calculated the field enhancement rate by using the bispherical coordinates [29]. Recently, Ammari and Yu developed a numerical method to compute the blow-up solution [44]. The blow-up rate for insulating inclusions is the same as perfectly conducting ones in two dimensions, but is different in three dimensions. Yun showed that the rate in three dimensions is 1/ɛ [48]. For the elasticity problem, Bao et al. obtained an upper bound in two dimensions [16], and this result was generalized to higher dimensions [17]. Recently, the gradient blow-up term of the solution in two dimensions was verified for convex hard inclusions by Kang and Yu [24] and for circular holes by Lim and Yu [35]. An inclusion with core-shell geometry can also induce the gradient blow-up for the solution to the conductivity problem as the distance between the boundaries of the core and shell degenerates to zero. In two dimensions, Ammari et al. obtained the blow-up rate ɛ 1/2 for core-shell geometry with circular boundaries [7]. Recently, Li et al. obtained an optimal bound for inclusions of general shapes in two dimensions and higher [3], and Bao et al. extended the result to the linear elasticity problem [15]. It is worth mentioning that a plasmonic structure of core-shell geometry induces the shielding effect on a region away from the structure [43]. Our goal in this paper is to characterize the electric field concentration for core-shell geometry with circular boundaries. We first show that only the linear term in Hx) contributes to the gradient blow-up. We then show that the gradient blow-up happens only in the shell area B e \B i. We finally express u by the single and double layer potentials with image line charges as follows: for x B e \ B i, ux) Hx) [ = C xr 2 ln x s ϕ 1 s) dσs) 2,),p 1 ] [ + C yr 2 x s, e 2 2 x s 2 ψ 1s) dσs),),p 1 ] [p 2,c i ] [p 2,c i ] ln x s ϕ 2 s) dσs) ] x s, e 2 x s 2 ψ 2s) dσs) + rx), where C x, C y, r 2 are some constants depending on the background potential function H and the geometry of the inclusion, and r is uniformly bounded independently of ɛ, x. The points p 1, p 2 and c i are defined according to B i and B e ; see Theorem 5.2 for further details. In the derivation, we use the bipolar coordinate system due to its convenience to solve the conductivity problem in the presence of two circular interfaces. We also deduce an integral expression with image line charges for the problem relating to two separated disks. The paper is organized as follows. In section 2, we explain the core-shell geometry and the associated bipolar coordinate system. We decompose the solution u into a singular term and a regular term in section 3 and derive an asymptotic formula for the solution in terms of bipolar coordinates in section 4. Section 5 is about the image line charge formula, and we conclude with some discussion. ] 3

4 2 The bipolar coordinate system 2.1 Geometry of the core-shell structure We let B i and B e be the disks centered at c i and c e with the radii r i and r e, respectively. As stated before, we assume that B e contains B i and that the distance between the boundary circles, say ɛ = dist B i, B e ), is small. After appropriate shifting and rotation, the disk centers are c i = c i, ) and c e = c e, ) with c i = r2 e ri 2 r e r i ɛ) 2, c e = c i + r e r i ɛ. 2.1) 2r e r i ɛ) Note that r e > r i > and c e > c i > as shown in Figure 2.1. The closest points on B i and B e between B i and B e are, respectively, x i = x i, ) and x e = x e, ), with Note that x i, x e = Oɛ). x i = c i r i = 2ɛr e ɛ 2 2r e r i ɛ), x e = c e r e = 2ɛr i + ɛ 2 2r e r i ɛ). Figure 2.1: An inclusion with eccentric core-shell geometry left) and the coordinate level curves of associated bipolar coordinates right). For a disk D centered at c with radius r, the reflection of x = x, y) R 2 across D is R D x) = r2 x c) x c 2 + c. We also define R D z) C for z = x + iy such that R D x) = R{R D z)}, I{R D z)}). Here the symbols R{w} and I{w} indicate the real and imaginary parts of a complex number w, respectively. We then define the reflection of a function f across D as and the combined reflections as R D [f]x) = f R D x)) for x R 2, R D1 R D2 R Dn [f]x) = f R Dn R Dn 1 R D1 )x) ). 4

5 Let R i and R e denote the reflections across the circles B i and B e, respectively. Then, one can easily show that the combined reflections R i R e and R e R i have the same fixed points p 1 = α, ) and p 2 = α, ) with α = ɛ2ri + ɛ)2r e ɛ)2r e 2r i ɛ). 2.2) 2r e r i ɛ) One can also easily show that α = r ɛ + Oɛ ɛ), where r = 2ri r e r e r i. 2.3) The boundaries of B i and B e are then the coordinate level curves {ξ = ξ i } and {ξ = ξ e }, respectively, in the bipolar coordinate system ξ, θ) that will be defined below, with ξ i = ln c i + α lnr i ), ξ e = ln c e + α lnr e ). 2.4) For later use we set { ξ i,k = 2kξ i ξ e ) + ξ i, ξ e,k = 2kξ i ξ e ) + ξ e 2.5) for each k N. One can show so that ξ i,k ξ i and ξ e,k+1 ξ e. ξ i > ξ e >, 2.2 The bipolar coordinate system associated with the eccentric core-shell structure We set the bipolar coordinate system ξ, θ) R π, π] with respect to the two poles located at p 1 = α, ) and p 2 = α, ): for any point x = x, y) in the Cartesian coordinate system, the coordinates ξ, θ) are defined such that e ξ+iθ = α + z α z with z = x + iy. 2.6) In other words z = α eξ+iθ 1 e ξ+iθ + 1 = α sinh ξ cosh ξ + cos θ + iα sin θ cosh ξ + cos θ. 2.7) We write either x = xξ, θ) or z = zξ, θ) to specify bipolar coordinates if necessary. For a function u = ux), we may write uξ, θ) to indicate uxξ, θ)) for notational simplicity. Figure 2.1 illustrates the coordinate level curves of the bipolar coordinate system. Note that the two level curves {θ = π 2 } and {θ = π 2 } form one circle that is centered at the origin with radius α. Therefore, we have that zξ, θ) α if θ π 2, zξ, θ) > α if θ > π ) 5

6 One can easily see that the two scale factors coincide, i.e., x ξ = x θ = 1 hξ, θ), where In particular, the components of x/ ξ x ξ satisfy that x, y ξ ξ hξ, θ) = = α 1 + cos θ cosh ξ cosh ξ + cos θ) 2, 1 hξ,θ) u = hξ, θ) u ξ êξ + hξ, θ) u θ êθ, cosh ξ + cos θ. α y ξ = α sin θ sinh ξ cosh ξ + cos θ) 2 2.9). The gradient for a function ux, y) can be written as with ê ξ = x/ ξ x/ ξ, ê θ = x/ θ x/ θ. 2.1) For any ξ, the level curve {ξ = ξ } is the circle centered at c = c, ) with radius r such that c = α e2ξ + 1 e 2ξ, r = α. 2.11) 1 sinh ξ From 2.7), the center c has the bipolar coordinates 2ξ, π). On this circle, the normal and tangential derivatives of a function u satisfy that u = sgnξ )hξ, θ) u, 2.12) ν ξ=ξ ξ ξ=ξ u = sgnξ )hξ, θ) u, T ξ=ξ θ ξ=ξ where ν is the outward unit normal vector and T is the rotation of ν by π 2 -radian angle. We may simply write u ν ξ, θ) for u ν ξ=ξ ξ, θ). The reflection across the circle {ξ = ξ } is actually a reflection with respect to ξ-variable, as follows. Lemma 2.1. For a disk D whose boundary is the coordinate level curve {ξ = ξ } with some ξ, we have R D xξ, θ)) = x2ξ ξ, θ), for all ξ, θ) 2ξ, π). Proof. Assume that D has the center c = c, ) and the radius r. From 2.11), we deduce that c = x2ξ, π) and c + α)c α) = α2 sinh 2 = r 2, ξ For z = x + iy c + i, we have α + R D z) α R D z) = α + r 2 z c) z c 2 α r2 z c) z c 2 = c + α c α + c c = c + α c α e 2ξ = c + α c α. z c)c α) z c z c)c+α) z c c α) + z c) c + α) + z c) = c + α c α α z α + z = e2ξ ξ+iθ. We prove the lemma by use of 2.6). 6

7 3 Decomposition of the solution in singular and regular terms The gradient blow-up of the solution to 1.1) may or may not blow up depending on the applied field H. In this section, we show that only the linear term of H contributes to the blow-up feature. This property naturally leads to decomposition of the solution into a singular term that blows up as ɛ tends to zero) and a regular term that remains bounded independently of ɛ); see Theorem 3.7 for further details. 3.1 Layer potential formulation Let D be a bounded domain in R 2 with a Lipschitz boundary D. We define the single layer potential S D and the Neumann-Poincaré operator K D as S D [ϕ]x) = 1 ln x y ϕy) dσy), x R 2, 2π D K D [ϕ]x) = 1 2π p.v. x y, ν x x y 2 ϕy) dσy), x D, D where ϕ L 2 D). Here p.v. means the Cauchy principal value. They satisfy the relation ν S D[ϕ] ± = ± 1 ) 2 + K D ϕ on D, where ν is the outward normal vector on D. For a C 1,α -domain D, the operator K D is compact on L 2 D) and the eigenvalues of K D on L2 D) lie in 1/2, 1/2]; see [27]. The operator λi K D is invertible on L2 Ω) for λ 1/2 [5]. For more properties of the Neumann- Poincaré, see [4, 6]. The solution u to the conductivity equation 1.1) satisfies the transmission conditions and it can be expressed as u + u = k e on B e, ν ν u k e + u = k i on B i, ν ν where ϕ i, ϕ e ) are mean-zero functions satisfying with ux) = Hx) + S Bi [ϕ i ]x) + S Be [ϕ e ]x), 3.1) λ i I K B i S B i ν e λ i = S B e ν i λ e I K B e ϕi ϕ e ) H = ν H ν Bi Be 3.2) k i + k e 2k i k e ), λ e = k e + 1 2k e 1). 3.3) Here, ν i and ν e denote the outward unit normal vectors on B i and B e, respectively. Note that λ i, λ e > 1 2 due to the positive assumption on the conductivities k i, k e. The problem 3.1) is 7

8 solvable for λ i, λ e > 1 2 as shown in [1]. Since B i and B e are disks, K B i and K B e are zero operators on L 2 B i) and L 2 B e), respectively. So, the system of boundary integral equations 3.2) becomes λ i I S B i ν e S B e ν i λ e I ϕi ϕ e ) = H ν H ν Bi Be. 3.4) We can relate the tangential component and the normal component of the electric field on B i and B e. The following lemma will be used later. One can prove the lemma by slightly modifying the proof of [9, Theorem 2]. Lemma 3.1. Let u be the solution to 1.1) with an entire harmonic function H. Let H be the harmonic conjugate of H and σ = 1 k i χb i ) + 1 k e χb e \ B i ) + χr 2 \ B e ). Then, the solution v to the problem { σ v = in R 2, satisfies that vx) Hx) = O x 1 ) as x u T = v + u on B e, ν T = 1 v + on B i. k e ν 3.2 Series expansions for the density functions in bipolar coordinates The single layer potential associated with a disk can be written in terms of the reflection with respect to the disk s boundary. Lemma 3.2 [9]). Let D be a disk centered at c and ν be the outward unit normal vector on D. If v is harmonic in D and continuous on D, then we have [ v ] S D 1 vc) vx) + for x D, 2 2 x) = ν D 1 2 R D[v]x) + vc) for x R 2 \ D. 2 If v is harmonic in R 2 \ D, continuous on R 2 \ D and vx) as x, then we have [ v S D ν + D ] x) = 1 2 R D[v]x) for x D, 1 2 vx) for x R2 \ D. We now express the density functions ϕ i and ϕ e in 3.1) in series form by use of Lemma 2.1. Lemma 3.3. The density functions in 3.1) can be expressed as λ i ϕ i θ) = 1 1 ) 1 ) m hξ i, θ) H 2λ e 4λ m= i λ e hξ i,m, θ) ν ξ i,m, θ), λ e ϕ e θ) = H ν ξ e, θ) ) 1 ) m hξ e, θ) H 2λ i 2λ e 4λ i λ e hξ e,m+1, θ) ν ξ e,m+1, θ). m= 3.5) 8

9 Proof. First, we formally expand the solution ϕ i, ϕ e ) to 3.4) as follows uniform and absolute convergence will be proved later): λ i ϕ i = 1 1 ) 1 ) m Re R i ) m H ) Bi, 2λ e 4λ m= i λ e ν i λ e ϕ e = H Be ) 1 ) m Ri R e ) m R i H ) 3.6) Be. ν e 2λ i 2λ e 4λ i λ e ν e Using Lemma 3.2, we compute S Bi [ϕ i ] ν e and S Be [ϕ e ] = 1 S Be ν i λ e ν i = ) λ i 2λ e m= = λ i 2λ e = λ e ϕ e H ν e Be, 1 2λ i λ e [ H ] Be ν e 1 1 2λ e ) ) m= = 1 2λ e H ν i Bi 1 4λ i λ e = λ i ϕ i H ν i Bi. m= 1 m= ) m S Bi 4λ i λ e ν e 1 4λ i λ e λ e m= [ Re R i ) m H ) ] Bi ν i ) m ν e Ri R e R i ) m H ) Bi ) m [ S Be Ri R e ) m R i H ) ] Be 4λ i λ e ν i ν e ) 1 ) m Ri R e ) m+1 H ) Bi 4λ i λ e ν i Therefore the series solutions in 3.6) satisfy 3.4). For the entire harmonic function H, we have from Lemma 2.1 that R e R i Hξ, θ) = HR i R e ξ, θ)) = H2ξ i ξ e ) + ξ, θ), R i R e R i Hξ, θ) = HR i R e R i ξ, θ)) = H2ξ i ξ e ) + 2ξ i ξ, θ). Repeating the reflections, it follows that R e R i ) m Hξ, θ) = H2mξ i ξ e ) + ξ, θ), R i R e ) m R i Hξ, θ) = H2mξ i ξ e ) + 2ξ i ξ, θ). Recall that ξ i,m and ξ e,m+1 are given as 2.5). By use of 2.12) we deduce that Re R i ) m H ) Bi = hξ i, θ) ν i ξ H 2mξ i ξ e ) + ξ, θ)) ξ=ξi = hξ i, θ) H hξ i,m, θ) ν ξ i,m, θ), Ri R e ) m R i H ) Be = hξ e, θ) ν e ξ H 2mξ i ξ e ) + 2ξ i ξ, θ)) = hξ e, θ) H hξ e,m+1, θ) ν ξ e,m+1, θ). ξ=ξe 9

10 Therefore, the equation 3.6) becomes 3.5). Now, we show uniform and absolute convergence for the series expansions. ξ i,m ξ i and ξ e,m+1 ξ e. Since xξ i,m, θ), xξ e,m+1, θ) B e for all m and Recall that hξ i, θ) hξ i,m, θ) 1, hξ e, θ) hξ e,m+1, θ) 1, we have hξ i, θ) H hξ i,m, θ) ν ξ i,m, θ) Hx) L B e), hξ e, θ) H hξ e,m+1, θ) ν ξ e,m+1, θ) Hx) L B e). Note that 1 4λ i λ e < 1 due to λ i, λ e > 1 2. Therefore, each term in 3.6) in other words, each term in 3.5)) is bounded 1 m by Hx) L B e) 4λ i λ e. This implies convergence and so the proof is completed. Corollary 3.4. If k i and k e are bounded but not small, then u H) is uniformly bounded in R 2 independently of ɛ. Proof. From Lemma 3.3 and the analysis in its proof, we have ϕ i L B i ), ϕ e L B e) 6 H L B e) λ i λ e The jump formula for the single layer potential and 3.4) implies u ± L = ν λ i ± 1 B i ) 2 ϕ i L B i ). 3.7) Therefore we have uniform boundedness for normal derivatives of u on B i and B e. We can also prove uniform boundedness for tangential derivatives by using a harmonic conjugate of u, in view of Lemma 3.1. Since u H)x) satisfies the decay condition as x, u H) is uniformly bounded in R 2 independently of ɛ by the maximum principle. 3.3 Blow-up contribution term in H and decomposition of the solution We prove the following upper bounds of two series sums related with the scale factor function h, which is essential for showing that only the linear term of H contributes to the gradient blow-up. 1

11 Lemma 3.5. Let N ɛ = r ɛ, then we have N ɛ m= N ɛ m= hξ i, θ) hξ i,m, θ) hξ e, θ) hξ e,m+1, θ) where C is a constant independent of ɛ, θ. C zξ i, θ) C zξ e, θ) for all θ > π 2, for all θ > π 2, Proof. We prove the first inequality only. The second one can be shown in the same way. Note that hξ i, θ) hξ i,m, θ) = cosh ξ i + cos θ cosh ξ i,m + cos θ = cosh ξ i + cos θ) cosh ξ i,m cos θ) cosh ξ i cos θ) cosh ξ i cos θ) cosh ξ i,m + cos θ) cosh ξ i,m cos θ) cosh ξ i + cos θ) α 2 cosh ξ i,m cos θ) α 2 = zξ i,m, θ) 2 cosh ξ i cos θ) cosh ξ i,m + cos θ) zξ i, θ). 3.8) We use 2.6), to estimate the last term in the above equation, and obtain with γ = e 2ξ i ξ e). This implies α + zξ i,m, θ) α zξ i,m, θ) = eξ i,m+iθ = e 2mξ i ξ e) e ξi+iθ = γ m α + zξ i, θ) α zξ i, θ) zξ i,m, θ) = α zξ i, θ)γ m + 1) + αγ m 1) zξ i, θ)γ m 1) + αγ m + 1) = 1 + α γ m 1 zξ i, θ) γ m γ m. 1 α γ m zξ i, θ) From 2.3) and 2.11), we deduce ξ i = r r i ɛ + Oɛ ɛ) and ξe = r r e ɛ + Oɛ ɛ) so that ξ i ξ e = 2 r ɛ + Oɛ ɛ) 3.9) and r ɛ < γ = e 2ξ i ξ e) < r ɛ. Hence, by using 2.8) as well, we have for 1 m < r ɛ and θ > π 2 that and 1 + 2m ) ɛ < γ m < 1 + 5r r ɛ ɛ < e 5 r zξ i,m, θ) 1 + α e 5 1 zξ i,θ) 2 2m ɛ/r αe 5 +1) + 1 zξ i,θ) C 1 C 2 m + 1 zξ i,θ) with some positive constants C 1 and C 2 independent of ɛ, θ and m. Hence, we have zξ i,m, θ) zξ i, θ) C for 1 m < r, θ > π m zξ i, θ) + 1 ɛ ) 11

12 Note that for any t >, the following holds: m= 1 mt + 1) xt + 1) 2 dx = t. This completes the proof because of 3.8), 3.1) and the boundedness of zξ i, θ). Now, we show that u H) is bounded regardless of ɛ if Hx) = O x ) near the origin. Lemma 3.6. If Hx) = O x ) in a bounded region containing B e, then we have u H) L R 2 ) C for some constant C independent of k i, k e and ɛ. Proof. Let H denote the gradient with respect to the Cartesian coordinates, and ξ i,m, ξ e,m+1 are defined as 2.5). We first show that hξ i, θ) Hξ i,m, θ) C. 3.11) hξ i,m, θ) m= Here and after, C denotes a constant independent of ɛ, θ and m. For notational simplicity we decompose the series in 3.11) as hξ i, θ) Hξ i,m, θ) = hξ i,m, θ) m= + m r ɛ m> r ɛ := I + II and separately prove the uniform boundedness. assumption Hx) = O x ) in B e that For θ π 2, we have from 2.8) and the I C r ɛ zξ i,m, θ) C r ɛ α = O1). For θ > π 2, we have from Lemma 3.5 and 3.1) { C I zξi,m, zξ i, θ) max θ) } r : m = O1). ɛ We now estimate II. For m > r ɛ, we have ξ i,m > 1 due to 3.9) so that and hξ i, θ) hξ i,m, θ) = Therefore, it follows that cosh ξ i + cos θ cosh ξ i,m + cos θ m> r zξ i,m, θ) < Cα cosh ξ i 1 2 cosh ξ i,m II 4e 2mξi ξe) Czξ i,m, θ) ɛ 4e ξ i ξ i,m = 4e 2mξ i ξ e). Cα ξ i ξ e = O1). 12

13 In the same way, one can easily show that m= hξ i, θ) hξ e,m+1, θ) Hξ e,m+1, θ) C. 3.12) From 3.11), 3.12) and Lemma 3.3 we obtain ϕ i < C 1 1 ), ϕ e < C λ i 2λ e λ e and, thus, u ν = λ i ± 1 ± 2 ϕ i L B i ) C. L B i ) The same result holds on B e. We can also prove uniform boundedness for tangential derivatives by using a harmonic conjugate for the harmonic function; see Lemma 3.1. Following the same argument in the proof of cor 3.4, we have uniform boundedness for u H) in R 2 independently of ɛ. As a direct consequence of Lemma 3.6, we can decompose the solution into a singular term and a regular term as follows. Theorem 3.7. Let H be an arbitrary entire harmonic function. Set H 1 x) = H) x and H 2 x) = Hx) H 1 x). Then, the solution u to the transmission problem 1.1) has the decomposition ux) = u 1 x) + u 2 x), x R 2, 3.13) where u 1 and u 2 are the solutions to 1.1) with H 1 and H 2 in the place of H, respectively, and u 2 H 2 )x) L R 2 ) < C for some constant C independent of ɛ. Corollary 3.8. If λ i λ e > 1 4 equivalently, either k i > k e > 1 or k i < k e < 1), then for any entire function H we have u H) L R 2 ) C with some positive constant C independent of ɛ. Proof. Thanks to Theorem 3.7, we may consider only the linear function for H, say Hx, y) = ax + by. Then, from 2.12) we have 1 ) m hξ i, θ) H 4λ m= i λ e hξ i,m, θ) ν ξ i,m, θ) = 1 ) m hξ i, θ) a x ) 4λ i λ e ξ + b y ξ m= ξ=ξ i,m. x Note that from 2.9), ξ and their derivatives in ξ have at most a finite number of changes in sign when θ is fixed. Moreover, hξ i, θ) x ξ ξ i,m, θ), hξ i, θ) y ξ ξ i,m, θ) hξ i, θ) hξ i,m, θ) 1. ξ, y 13

14 Therefore, from the convergence property of alternating series we derive that 1 ) m hξ i, θ) H 4λ i λ e hξ i,m, θ) ν ξ i,m, θ) C H) and, in the same way, m= m= 1 ) m hξ e, θ) H 4λ i λ e hξ e,m+1, θ) ν ξ e,m+1, θ) C H). Hence, ϕ i, ϕ e are bounded independent of ɛ. Following the same argument in cor 3.4 we have uniform boundedness for u H). For the case of separated disks, a boundedness result comparable to Lemma 3.6 was derived in [8], based on the series expression in the Cartesian coordinates that is similar to 3.6). We emphasize that in this paper, we provide a much simpler proof thanks to Lemma Asymptotics of the singular term in bipolar coordinates From the analysis in the previous section, the gradient blow-up does not occur either when Hx) = O x ) or when λ i λ e > 1 4 equivalently, either k i > k e > 1 or k i < k e < 1). In this section, we first derive an asymptotic expansion for the solution to 1.1) assuming that H is a linear function and λ i λ e < 1 4 equivalently, either k i > k e, 1 > k e or k i < k e, 1 < k e ). We then characterize the gradient blow-up. We follow the derivation in [34], where the asymptotic behavior for the electric field was derived in the presence of two separated circular conductors. The gradient blow-up feature is essentially related to the following quantities: and τ i = k i k e k i + k e, τ e = 1 k e, τ = τ i τ e, β = r ln τ) 1 + k e 4, 4.1) ɛ C x = H) e 1 ), C y = H) e 2 ), 4.2) where e 1 = 1, ) and e 2 =, 1). Note that τ i = 1 2λ i, τ e = 1 2λ e and β >. 4.1 Series solution by the separation of variables In bipolar coordinates, harmonic functions admit the general solution expansion fξ, θ) = a + b ξ + c θ + n=1 [ ] a n e nξ + b n e nξ ) cos nθ + c n e nξ + d n e nξ ) sin nθ, where a n, b n, c n and d n are some constants. The linear functions x and y can be expanded as [ ] x = sgnξ)α ) n e n ξ cos nθ, y = 2α n=1 1) n e n ξ sin nθ. n=1 14

15 Lemma 4.1. We define A n e nξ iθ 2ξi) + B n e nξ iθ 2ξe)) for ξ < ξ e, with n=1 Ux) := C + A n e nξ iθ 2ξi) + B n e n ξ iθ)) for ξ e < ξ < ξ i, n=1 A n e n ξ iθ) + B n e n ξ iθ)) n=1 2α 1) n A n = τ 1 e 2nξ i ξ e) 1 2α 1) n B n = τ 1 e 2nξ i ξ e) 1 ) 1 τ 1 e e 2nξ i ξ e), 1 + τ 1 i e 2nξ i ξ e) ), for ξ > ξ i and C = n=1 A ne n 2ξ i iπ) + B n e inπ. Then, the functions x + R{Ux)} and y + I{Ux)} are the solutions to 1.1) for Hx) = x and Hx) = y, respectively. Proof. One can easily show that x+r{ux)} and y+i{ux)} satisfy the transmission condition in the equation 1.1) and demonstrate uniform convergence. We can prove the decay property in the same way as in Lemma 3.1 of [34]. Definition 1. The Lerch transcendent function is defined as We also define Lz; β) = ze β+1)t dt for z C, z < ) 1 + ze t P z; β) = z) z Lz; β) = ze β+1)t 1 + ze t ) 2 dt. Let us derive some properties on P for later use. By the integration by parts it follows that for all s >, P e s iθ ; β) = One can easily show that P e s+iθ ) = e βt e s t iθ 1 + e s t iθ ) 2 dt βt 1 + coshs + t) cos θ i sinhs + t) sin θ e 2coshs + t) + cos θ) 2 dt. 4.4) e βt 1 2coshs + t) + cos θ) dt 1 2βcosh s + cos θ). 4.5) For any s 2 > s 1 >, β > and θ π, the following holds: P e s2+iθ ; β) P e s1+iθ s 2 s 2 ; β) 2cosh s 1 + cos θ). 4.6) The following integral approximation for an infinite series sum is crucial to obtaining the gradient blow-up term of u. 15

16 Lemma 4.2 [34]). Let < a < a and < τ < 1. For arbitrary θ π, π], we have a τ m 1 e ma ) +a+iθ 1 + e ma+a+iθ ) 2 P e a a)+iθ ; ln τ ) a 4a cosha a) + cos θ. m=1 4.2 Integral expression of the singular term in bipolar coordinates We verify the asymptotic behavior of u in bipolar coordinates as follows. We will later express the gradient blow-up term of u in terms of the Cartesian coordinates in section 5. Theorem 4.3. Let u be the solution to 1.1). Then we have the following: a) The solution u admits that ux) = Hx) + C x r 2 R { q x; β, τ i, τ e )} + Cy r 2 I { q x; β, τ i, τ e )} + rx), x R 2, 4.7) where the singular function q is qx; β, τ i, τ e ) := 1 ) τ + τ i Le ξ iθ 2ξ i ; β) + ) τ + τ e Le ξ iθ 2ξ e ; β) in R 2 \ B e, ) τ + τ 2 i Le ξ iθ 2ξ i ; β) + ) τ + τ e Le ξ iθ ; β) in B e \ B i, ) τ + τ i Le ξ iθ ; β) + ) τ + τ e Le ξ iθ ; β) in B i, and rx) L R 2 ) < C for some constant C independent of ɛ. b) For any k i, k e, we have that u H) L R 2 \B e\b i C for some constant C independent of ɛ. The gradient blow-up occurs only when τ i, τ e are similar to 1 equivalently, )) < k e 1, k i ). For such a case, we have that in B e \ B i, u H)x) = C xr τ i + τ e + 2τ) 2 ɛ + C yr τ i + τ e + 2τ) 2 ɛ ) { } cosh ξ + cos θ R P e ξ iθ ; β) ê ξ ) { } cosh ξ + cos θ I P e ξ iθ ; β) ê ξ + O1), where O1) is uniformly bounded independently of ɛ and x. Proof. In view of Theorem 3.7, we may consider only the linear function for H. In other words, Hx) = H) 1, ))x + H), 1))y. Moreover, we can assume τ = τ i τ e > since the blow-up happens only in such a case. One can derive a series expansion for Ux) by differentiating the formula in Lemma

17 From the definition of P, we have q ξ = 1 τ + τ i )P e ξ iθ 2ξ i; β) τ + τ e )P e ξ iθ 2ξe ; β) in R 2 \ B e, τ + τ 2 i )P e ξ iθ 2ξ i; β) + τ + τ e )P e ξ iθ ; β) in B e \ B i, τ + τ i )P e ξ iθ ; β) + τ + τ e )P e ξ iθ ; β) in B i, q θ = i τ + τ i )P e ξ iθ 2ξ i; β) + τ + τ e )P e ξ iθ 2ξe ; β) in R 2 \ B e, τ + τ 2 i )P e ξ iθ 2ξ i; β) + τ + τ e )P e ξ iθ ; β) in B e \ B i, τ + τ i )P e ξ iθ ; β) + τ + τ e )P e ξ iθ ; β) in B i. Thanks to Lemma 4.2, one can show that Ux) = r 2 qx; β, τ i, τ e ) + O1), x R 2, 4.8) in a similar way as in the proof of [34, Lemma 5.1]. Since the functions x + R{Ux)} and y + I{Ux)} are, respectively, the solutions to 1.1) with Hx) = x and Hx) = y, we complete the proof of a). Now we prove b). From 4.8) it follows that u H)x) =C x r 2 R{ qx; β, τ i, τ e )} ê ξ )ê ξ + C x r 2 R{ qx; β, τ i, τ e )} ê θ )ê θ +C y r 2 I{ qx; β, τ i, τ e )} ê ξ )ê ξ + C y r 2 I{ qx; β, τ i, τ e )} ê θ )ê θ + O1). Using 2.1) and 4.6), we can show in a similar way as in [34, section 4.1] that q ê ξ = O1) + 1 τ i τ e )P e ξe iθ ; β) in R 2 \ B e, hξ, θ) τ 2 i + τ e + 2τ)P e ξ iθ ; β) in B e \ B i, τ e τ i )P e ξ iθ ; β) in B i, q ê θ = O1) + i 2 hξ, θ)τ e τ i )P e ξ iθ ; β) in R 2. By use of 4.5), we obtain that hξ, θ)p e ξ iθ ; β) 1 αβ Note that we have C ln τ. τ i τ e τ i τ e ln τ i ln τ e ln τ owing to τ > and τ i, τ e < 1. For τ i, τ e <, it follows that q C in R 2, so that the blow-up does not happen. For τ i, τ e >, we have that q ê ξ C in R 2 \ B e ) B i, q ê θ C in R 2. From 4.9) and the definition of h, α, we can derive the formula in b) and complete the proof. Figure 4.1 shows the graphs of the normal and tangential components of u on B e. As stated in Theorem 4.3b), only the normal component blows up as ɛ tends to zero. 4.9) 17

18 a) ɛ =.1, Hx) = x b) u H) ê ξ, ɛ =.1 c) u H) ê θ, ɛ =.1 d) ɛ =.1, Hx) = y e) u H) ê ξ, ɛ =.1 f) u H) ê θ, ɛ =.1 Figure 4.1: Electric field for the core-shell structure with r i = 1, r e = 2, k i = 1, and k e =.1, where ɛ =.1 or.1. The figures are contour plots of u H) the first column) and the graphs of u H) ê ξ and u H) ê θ on B e the second and third column, solid curve). Dashed curves are computed by using the singular term in Theorem 4.3a). The background potential is Hx) = x the first row) or Hx) = y the second row). 4.3 Optimal bounds Assume that k i = 1 and k e 1. Then, by using 4.4), one can easily deduce the upper and lower bounds for cosh ξ i + cos θ)p e ξi iθ ; β). Specifically, we have that cosh ξi + cos θ ) { } R P e ξi iθ ; β) C 1 for θ =, 4.1) cosh ξi + cos θ ) I { P e ξ i iθ ; β)} β + 1 C 1 β + 1)β + 3) for θ = π 2 with some positive constant C 1 independent of ɛ, β. Indeed, for θ = π 2 we have {P I e ξ i π 2 i ; β)} = e βt sinhξ i + t) 2 cosh 2 ξ i + t) dt = eβξ i e βt e t e t ξ i e t + e t ) 2 dt e βξ i e βt et e t ξ i 2e t ) 2 dt 1 e ξ i ) 4 β + 1 e 3ξi β + 3 C 1 β + 1)β + 3). 4.11) 18

19 Also, we can show that cosh ξ i + cos θ)p e ξi iθ ; β) C 2 β + 1 for some positive constant C 2 independent of ɛ, β, θ. Note that we have owing to β = 1 ɛβ + 1) 1 r k e + ɛ r 4 1 τ + O1 τ) 2 ) ) r k e. ɛ ɛ 4.12) In total, assuming k i = 1 and k e 1, the following holds with some positive constant C 1, C 2 independent of k e, k i, ɛ: If C x and C y =, then we have C 1 C x ɛ k e + If C x =, C y and r ke ɛ then we have C 1 C x u H) L B e\b i ) C 2 r k e + ɛ. r < C for some constant C independent of ɛ i.e., β is bounded), C y ɛ k e + 5 Image line charge formula C y u H) L B e\b i ) C 2 r k e + We obtained an asymptotic formula for the electric field in bipolar coordinates in the previous section, where the singular term is a linear combination of the Lerch transcendent functions. In this section, we rewrite the formula in terms of the Cartesian coordinates. It turns out that the singular term is actually the single and double layer potentials with line charges on the x-axis. Also, we derive an asymptotic formula for the case of two separated disks. As well as the core-shell geometry, the singular term for two separated disks can be expressed as the single and double layer potentials with line charges on the x-axis. 5.1 Line charge formula of the Lerch transcendent function We first express each component of the singular function q, that is, a linear combination of the Lerch transcendent function, in the Cartesian coordinates. Before stating the formulas, let us define some notations. Here and after, we denote by [x 1, x 2 ] the line segment connecting the two points x 1, x 2. We also write fs, ) = fs) for a function f defined on R. Lemma 5.1. Assume ξ. We denote by c = c, ) and r the center and radius of the level curve {ξ = ξ } for ξ. We denote by c = +, ) and r = for ξ =. We set the charge distribution functions as ϕ + s) = 2αβe 2βξ s α)β 1 s + α) β+1, ψ +s) = e 2βξ ϕ s) = 2αβe 2βξ s + α)β 1 s α) β+1, ψ s) = e 2βξ ɛ. r ) s α β for α < s < c, s + α ) s + α β for c < s < α. s α 19

20 Then, the following holds: a) For ξ < 2ξ, we have Le ξ+iθ 2ξ ; β) = [p 2, c ] x s, e 2 ln x s ϕ + s) dσs) + r + x) i [p 2, c ] x s 2 ψ + s) dσs), where the remainder term satisfies that r + x) 1 r for all x if ξ. b) For ξ > 2ξ, we have Le ξ iθ 2ξ ; β) = [ c, p 1 ] x s, e 2 ln x s ϕ s) dσs)+r x) i [ c, p 1 ] x s 2 ψ s) dσs), where the remainder term satisfies that r x) 1 r for all x if ξ. Proof. By integration by parts we have Lz; β) = ze t e βt 1 + ze dt t = ln1 + z) + β ln1 + ze t )e βt dt ) 1 + ze t = β ln e βt dt. 1 + z Here, we used the fact that β > to evaluate the integrand at t =. Owing to Lz; β) = Lz; β), it follows that R {Lz; β)} = 1 2 Lz; β) + Lz; β)) = β ln 1 + ze t 1 + z e βt dt. The imaginary part of Lz; β) becomes I { Lz; β) } = 1 Lz; β) Lz; β)) 2i = 1 ze β+1)t 2i 1 + ze 1 + ze = 1 2i ze β+1)t dt t t z z)e β+1)t 1 + ze t 2 dt = I{z}e β+1)t 1 + ze t 2 dt. 5.1) To prove a), we now assume ξ < 2ξ and set z = e ξ+iθ 2ξ. We first consider R {Lz; β)}. Note that z 1. We compute 1 + ze t = 1 + e ξ+iθ e 2ξ +t) ) = 1 e 2ξ+t) )1 + e ξ+iθ ) 1 e 2ξ+t) eξ+iθ 1 e 2ξ +t) = 1 e 2ξ+t) ) ) α zξ, θ) = 1 e 2ξ +t) α zξ, θ) 2α α zξ, θ) 2α zξ, θ) + α 1 + e 2ξ +t) 1 e 2ξ +t) ) + e 2ξ +t) 1 e 2ξ +t) ). 5.2) 2

21 For ease of estimation, we set In other words, From 5.2) it follows that st) = α 1 + e 2ξ +t) 1 e 2ξ +t). e t = e 2ξ st) α st) + α. 1 + ze t = 1 e 2ξ +t) zξ, θ) + st)). α zξ, θ) It is immediately clear that st) is the center of circle {ξ = ξ + t 2 } owing to 2.11). In particular, s) = c. Hence we have { } R Le ξ+iθ 2ξ ; β) = β ln 1 + ze t 1 + ze e βt dt 1 e 2ξ+t) ) zξ, θ) st)) = β ln 1 e 2ξ ) zξ, θ) s)) e βt dt = β Due to the fact that ln zξ, θ) st) e βt dt β ds = we have { } R Le ξ+iθ 2ξ ; β) = β with One can easily see that ln zξ, θ) s) e βt dt + L e 2ξ ; β). 2αe 2ξ +t) ) 1 e 2ξ +t) 2 dt = 1 s + α)s α)dt, 2α = 2αβe 2βξ ln zξ, θ) st) e βt dt ln zξ, θ) c + L e 2ξ ; β) c α ln zξ, θ) s s α)β 1 s + α) β+1 ds + r +zξ, θ)) r + zξ, θ)) = ln zξ, θ) c + L e 2ξ ; β). rzξ, θ)) 1 r for all ξ < 2ξ. We now consider I {Lz; β)}. Note that I { e ξ+iθ} { } α + zξ, θ) = I α zξ, θ) { } α + zξ, θ))α zξ, θ)) = I α zξ, θ) 2 2α = α zξ, θ) 2 I{ zξ, θ) }. 21

22 Using 5.1) and 5.2), we have { } I Le ξ+iθ 2ξ ; β) = = = e 2βξ I{e ξ+iθ 2ξ }e β+1)t 1 + e ξ+iθ 2ξ t 2 dt I { zξ, θ) } e βt zξ, θ) st) 2αe 2ξ+t) 2 ) 1 e 2ξ +t) 2 dt c α I { zξ, θ) } zξ, θ) s 2 ) s α β ds. s + α Note that I { zξ, θ) } = x s, e 2. Hence we complete the proof of a). We can prove b) in a similar way, where the remainder term is r x) = ln zξ, θ) + c + L e 2ξ ; β). 5.2 Asymptotic formula for the core-shell geometry We now express the singular part of u in terms of the single and double layer potential with image line charges on the x-axis. Theorem 5.2. Let u be the solution to 1.1) corresponding to an entire harmonic function H. For any < k i, k e 1 <, we have that u H) L R 2 \B e\b i C for some constant C )) independent of ɛ. The gradient blow-up for u may occur only when τ i, τ e are similar to 1 that is, equivalently, < k e 1, k i ). For k i, k e such that < k e < 1 and k e < k i, we have that in B e \ B i, [ ux) = Hx) + C xr 2 ] ln x s ϕ 1 s) dσs) ln x s ϕ 2 s) dσs) 2,), p 1 ] [p 2, c i ] [ + C yr 2 ] x s, e 2 2 x s 2 ψ x s, e 2 1s) dσs) x s 2 ψ 2s) dσs) + rx),,), p 1 ] where the charge distributions are [p 2, c i ] s + α)β 1 ϕ 1 s) = τ + τ e ) 2αβ s α) β+1 χ, α]s), ) s + α β ψ 1 s) = τ + τ e ) χ s α, α]s), ϕ 2 s) = τ + τ i ) 2αβe 2βξ s α)β 1 i s + α) β+1 χ [α, c i ]s), ) ψ 2 s) = τ + τ i ) e 2βξ i s α β χ s + α [α, ci]s), and rx) is uniformly bounded independently of ɛ and x. Proof. Applying Lemma 5.1 to Theorem 4.3 a) we prove the theorem. 22

23 Figure 5.1: Bold intervals indicate the support of line charges for the core-shell structure. a) k i = 1, k e =.3, ɛ =.1 β = , τ =.8869) b) k i = 1, k e =.158, ɛ =.1 β = 1, τ =.9387) c) k i = 1, k e =.1, ɛ =.1 β =.6325, τ =.968) Figure 5.2: Image line charges for the core-shell structure with r i = 1, r e = 2, k i = 1, ɛ =.1 and various k e -values. The left column shows 1 2 ϕ 1 ϕ 2 ), and the right column shows 1 2 ψ 1 ψ 2 ) in Theorem 5.2. As k e decreases, β approaches to. The density 1 2 ϕ 1 ϕ 2 ) is bounded if β 1, but diverges to near p 1 and p 2 if β < 1. 23

24 a) k i = 1, k e =.3, ɛ =.1 β =.62, τ =.8869) b) k i = 1, k e =.3, ɛ =.36 β = 1, τ =.8869) c) k i = 1, k e =.3, ɛ =.1 β = , τ =.8869) Figure 5.3: Image line charges for the core-shell structure with r i = 1, r e = 2, k i = 1, k e =.3 and various ɛ-values. The left column shows 1 2 ϕ 1 ϕ 2 ), and the right column shows 1 2 ψ 1 ψ 2 ) in Theorem 5.2. The support of line charges for the core-shell structure are drawn in Figure 5.1. As ɛ decreases, β increases so that 1 2 ϕ 1 ϕ 2 ) becomes bounded. 24

25 Figure 5.1 shows the support of the density functions in Theorem 5.2. One can easily see that ψ js) s = ϕ j s) for j = 1, 2 and from 2.11), ϕ 1 s) ds = τ + τ e, ϕ 2 s) ds = τ + τ i. The density functions show a different feature depending on β. The functions ϕ 1 and ϕ 2 are bounded if β 1, but they diverge to ± near p 1 or p 2 if β < 1. See Figure 5.2 and Figure 5.3. Remark 1. When β = equivalently τ i = τ e = 1), which is the limiting case of ln τ ɛ, we have that Le ξ+iθ 2ξ x s, e 2 ; ) = ln x p 2 + r + x) i [p 2, c ] x s 2 ds, ξ < 2ξ, Le ξ iθ 2ξ x s, e 2 ; ) = ln x p 1 + r x) + i x s 2 ds, ξ > 2ξ, [ c, p 1 ] where r + x) and r x) are bounded independently of ɛ, x. Hence, the singular term of u in 4.7) satisfies C x r R 2 { q x;, 1, 1 )} + C y r I 2 { q x;, 1, 1 )} = C x r 2 ln x p 1 x p 2 + C yr 2 [p 2, c ] [ c, p 1 ] where rx) are bounded independently of ɛ, x. x s, e 2 x s 2 ds + rx), x B e \ B i, 5.3) Recall that we denote by x i = x i, ) and x e = x e, ) the closest points on B i and B e between B i and B e, respectively. The gradient of the function in 5.3) blows up at x e as ɛ tends zero, since the distance between p 1 and p 2 is O ɛ). However, the solution to the following problem does not blow up: u = in R 2 \ B e, u ν = on B e, u x) Hx) = O x 1 ) as x. These facts imply that ux e ) for the solution u to 1.1) does not converge to u x e ) when k e and ɛ, satisfying that k i is fixed and ln τ ɛ. See [34, section 5.3] for a non-convergence result for u in the case of separated disks. 5.3 Asymptotic formula for the case of two separated disks We now consider the gradient blow-up feature for the solution u to 1.1) in the presence of two separated disks, say B 1 and B 2. In other words, the conductivity profile is given by 5.4) σ = k 1 χb 1 ) + k 2 χb 2 ) + χ R 2 \ B 1 B 2 ) ). 5.5) As before we assume that < k i, k e 1 <. For two separated disks with arbitrary conductivities, the asymptotic formula for the potential function in bipolar coordinates was derived in 25

26 [34]. In this section, we express the singular term by the single and double layer potential with line charges. We set B j to be the disk centered at c j, ) with radius r j, j = 1, 2. We denote by ɛ the distance between B 1 and B 2 ; see Figure 5.4 for the geometry of the two separated disks. We set R j to be the reflection with respect to B j, j = 1, 2. We can assume that, by appropriate shifting and rotation, the fixed points of the combined reflection R 1 R 2 are p 1 = α, ), p 2 = α, ) with ɛ2r1 + ɛ)2r 2 + ɛ)2r 1 + 2r 2 + ɛ) α =. 2r 1 + r 2 + ɛ) As for the core-shell geometry, we denote by ξ, θ) the bipolar coordinates defined with respect to the poles p 1 and p 2. We now set 2r1 r 2 r = r 1 + r 2 and τ 1 = k 1 1 k 1 + 1, τ 2 = k 2 1 k 2 + 1, τ = τ 1τ 2, β = r ln τ) 4. ɛ Lemma 5.3 [34]). Let u be the solution to 1.1) with σ given by 5.5). Then u satisfies that where ux) = Hx) + C x r 2 R { q x; β, τ i, τ e )} + Cy r 2 I { q x; β, τ i, τ e )} + rx), x R 2, 5.6) qx; β, τ 1, τ 2 ) := 1 τ + τ 1 )Le ξ+iθ) 2ξ 1 ; β) τ + τ 2 )Le ξ+iθ) 2ξ 2 ; β) in R 2 \ B 1 B 2 ), τ + τ 2 1 )Le ξ iθ ; β) τ + τ 2 )Le ξ+iθ) 2ξ 2 ; β) in B 1, τ + τ 1 )Le ξ+iθ) 2ξ 1 ; β) τ + τ 2 )Le ξ+iθ ; β) in B 2, and rx) L R 2 ) < C for some constant C independent of ɛ. Using the line charge formula of the Lerch transcendent function in section 5.1, one can easily show the following. Theorem 5.4. Let u be the solution to 1.1) with σ as given by 5.5). The gradient blow-up for u may occur only when τ is similar to 1 equivalently, k 1, k 2 1 or < k 1, k 1 1). Moreover, u satisfies the following integral formulas in R 2 \ B 1 B 2 ): a) For 1 < k 1, k 2 <, we have [ ux) = Hx) + C xr 2 ln x s ϕ 1 s) dσs) 2 [c 1,p 1 ] where the charge distributions are ϕ 1 s) = τ + τ 1 )2αβe 2βξ 1 ϕ 2 s) = τ + τ 2 )2αβe 2βξ 2 [p 2,c 2 ] s + α)β 1 s α) β+1 χ [c 1, α]s), s α)β 1 s + α) β+1 χ [α, c 2 ]s), and rx) L R 2 ) < C for some constant C independent of ɛ. 26 ln x s ϕ 2 s) dσs) ] + rx),

27 b) For < k 1, k 2 < 1, we have [ ux) = Hx) + C yr 2 2 [c 1,p 1 ] ] x s, ê y x s 2 ψ x s, ê y 1s) dσs) + [p 2,c 2 ] x s 2 ψ 2s) dσs) + rx), where the charge distributions are ) ψ 1 s) = τ + τ 1 )e 2βξ 1 s + α β χ s α [c1, α]s), ) ψ 2 s) = τ + τ 2 )e 2βξ 2 s α β χ s + α [α, c2]s), and rx) L R 2 ) < C for some constant C independent of ɛ. One can derive a similar line charge formula for u in B 1 B 2 by use of the line charge formula of the Lerch transcendent function in section Conclusion We have analyzed the gradient blow-up of the solution to the conductivity problem in the presence of an inclusion with core-shell geometry. We showed that the gradient blow-up happens only when the conductivities of the core and shell satisfy the conditions k e < 1 and k e < k i. We then derived an asymptotic formula for the potential function in terms of the bipolar coordinates. We finally expressed the asymptotic formula in terms of the single and double layer potentials with line charges. It will be of interest to generalize the result to core-shell geometry with arbitrary shapes, or to the elasticity problem. References [1] Habib Ammari, Giulio Ciraolo, Hyeonbae Kang, Hyundae Lee, and Graeme W Milton. Spectral theory of a Neumann Poincaré-type operator and analysis of cloaking due to anomalous localized resonance. Archive for Rational Mechanics and Analysis, 282): , 213. [2] Habib Ammari, Giulio Ciraolo, Hyeonbae Kang, Hyundae Lee, and Kihyun Yun. Spectral analysis of the Neumann Poincaré operator and characterization of the stress concentration in anti-plane elasticity. Archive for Rational Mechanics and Analysis, 281):275 34, 213. [3] Habib Ammari, George Dassios, Hyeonbae Kang, and Mikyoung Lim. Estimates for the electric field in the presence of adjacent perfectly conducting spheres. Quarterly of applied mathematics, 652): , 27. [4] Habib Ammari, Josselin Garnier, Wenjia Jing, Hyeonbae Kang, Mikyoung Lim, Knut Sølna, and Han Wang. Mathematical and statistical methods for multistatic imaging Vol. 298) Springer, 213. [5] Habib Ammari and Hyeonbae Kang. Reconstruction of small inhomogeneities from boundary measurements. Springer,

28 Figure 5.4: Geometry of two separated disks. Bold intervals indicate the support of line charges for the core-shell structure. a) k 1 = 1, k 2 = 15 left) and k 1 = 1/1, k 2 = 1/15 right). ɛ =.5 β = , τ =.7159) b) k 1 = , k 2 = left) and k 1 = 1/ , k 2 = 1/ right). ɛ =.5 β = 1, τ =.857) c) k 1 = 2, k 2 = 3 left) and k 1 = 1/2, k 2 = 1/3 right). ɛ =.5 β =.772, τ =.8464) Figure 5.5: Image line charges for the two separated disks with r 1 = 1.5, r 2 = 2, ɛ =.5 and various k 1, k 2 -values. The left column shows 1 2 ϕ 1 ϕ 2 ), and the right column shows 1 2 ψ 1 ψ 2 ) in Theorem

29 a) ɛ =.1 β =.5459, τ =.8464) b) ɛ =.3 β = 1, τ =.8464) c) ɛ =.1 β = , τ =.8464) Figure 5.6: Image line charges for the two separated disks with r 1 = 1.5, r 2 = 2 and various ɛ-values, where k 1 = 2, k 2 = 3 left column) and k 1 = 1/2, k 2 = 1/3 right column). The left column shows 1 2 ϕ 1 ϕ 2 ), and the right column shows 1 2 ψ 1 ψ 2 ) in Theorem 5.4. The support of line charges for the core-shell structure are drawn in Figure 5.4. As ɛ decreases, β increases so that 1 2 ϕ 1 ϕ 2 ) becomes bounded. 29

30 [6] Habib Ammari and Hyeonbae Kang. Polarization and Moment Tensors: With Applications to Inverse Problems and Effective Medium Theory, volume 162. Springer Science & Business Media, 27. [7] Habib Ammari, Hyeonbae Kang, Hyundae Lee, Jungwook Lee, and Mikyoung Lim. Optimal estimates for the electric field in two dimensions. Journal de mathématiques pures et appliquées, 884):37 324, 27. [8] Habib Ammari, Hyeonbae Kang, Hyundae Lee, Mikyoung Lim, and Habib Zribi. Decomposition theorems and fine estimates for electrical fields in the presence of closely located circular inclusions. Journal of Differential Equations, 24711): , 29. [9] Habib Ammari, Hyeonbae Kang, and Mikyoung Lim. Gradient estimates for solutions to the conductivity problem. Mathematische Annalen, 3322): , 25. [1] I. Babu ska, B. Andersson, P. Smith and K. Levin, Damage analysis of fiber composites. I. Statistical analysis on fiber scale, Comput. Methods Appl. Mech. Engrg ), [11] G. K. Batchelor and R. W. O Brien, Thermal or Electrical Conduction Through a Granular Material, Proc. Roy. Soc. A ), [12] B. Budiansky and G. F. Carrier, High shear stresses in stiff fiber composites, Jour. Appl. Mech ), [13] Ellen Shiting Bao, Yan Yan Li, and Biao Yin. Gradient estimates for the perfect conductivity problem. Archive for rational mechanics and analysis, 1931): , 29. [14] Ellen Shiting Bao, Yan Yan Li, and Biao Yin. Gradient estimates for the perfect and insulated conductivity problems with multiple inclusions. Communications in Partial Differential Equations, 3511): , 21. [15] JiGuang Bao, Hongjie Ju, and Haigang Li. Optimal boundary gradient estimates for Lamé systems with partially infinite coefficients. Advances in Mathematics, 314: , 217. [16] JiGuang Bao, HaiGang Li, and YanYan Li. Gradient estimates for solutions of the Lamé system with partially infinite coefficients. Archive for Rational Mechanics and Analysis, 2151):37 351, 215. [17] JiGuang Bao, HaiGang Li, and YanYan Li. Gradient estimates for solutions of the Lamé system with partially infinite coefficients in dimensions greater than two. Advances in Mathematics, 35: , 217. [18] Eric Bonnetier and Faouzi Triki. Pointwise bounds on the gradient and the spectrum of the neumann poincaré operator: the case of 2 discs. Contemp. Math, 577:81 92, 212. [19] Milford Harold Davis. Two charged spherical conductors in a uniform electric field: Forces and field strength. The Quarterly Journal of Mechanics and Applied Mathematics, 174): , [2] GB Jeffery. On a form of the solution of Laplace s equation suitable for problems relating to two spheres. Proceedings of the Royal Society of London. Series A, Containing Papers of a Mathematical and Physical Character, 87593):19 12,

31 [21] Hyeonbae Kang, Hyundae Lee, and KiHyun Yun. Optimal estimates and asymptotics for the stress concentration between closely located stiff inclusions. Mathematische Annalen, ): , 215. [22] Hyeonbae Kang, Mikyoung Lim, and KiHyun Yun. Asymptotics and computation of the solution to the conductivity equation in the presence of adjacent inclusions with extreme conductivities. Journal de Mathématiques Pures et Appliquées, 992): , 213. [23] Hyeonbae Kang, Mikyoung Lim, and KiHyun Yun. Characterization of the electric field concentration between two adjacent spherical perfect conductors. SIAM Journal on Applied Mathematics, 741): , 214. [24] Hyeonbae Kang and Sanghyeon Yu. Quantitative characterization of stress concentration in the presence of closely spaced hard inclusions in two-dimensional linear elasticity. arxiv preprint arxiv: , 217. [25] Hyeonbae Kang and KiHyun Yun. Optimal estimates of the field enhancement in presence of a bow-tie structure of perfectly conducting inclusions in two dimensions. arxiv preprint arxiv:177.98, 217. [26] Joseph B Keller. Conductivity of a medium containing a dense array of perfectly conducting spheres or cylinders or nonconducting cylinders. Journal of Applied Physics, 344): , [27] Oliver Dimon Kellogg. Foundations of potential theory, volume 31. Springer Science & Business Media, 212. [28] J. B. Lassiter, J. Aizpurua, L. I. Hernandez, D. W. Brandl, I. Romero, S. Lal, J. H. Hafner, P. Nordlander and N. J. Halas. Close encounters between two nanoshells. Nano letters, 84): , 28. [29] John Lekner. Near approach of two conducting spheres: enhancement of external electric field. Journal of Electrostatics, 696): , 211. [3] Haigang Li and Longjuan Xu. Optimal estimates for the perfect conductivity problem with inclusions close to the boundary. arxiv preprint arxiv: , 217. [31] Yan Yan Li and Michael Vogelius. Gradient estimates for solutions to divergence form elliptic equations with discontinuous coefficients. Archive for Rational Mechanics and Analysis, 1532):91 151, 2. [32] Yanyan Li and Louis Nirenberg. Estimates for elliptic systems from composite material. Communications on pure and applied mathematics, 567): , 23. [33] Mikyoung Lim and Sanghyeon Yu. Asymptotic analysis for superfocusing of the electric field in between two nearly touching metallic spheres. arxiv preprint arxiv: , 214. [34] Mikyoung Lim and Sanghyeon Yu. Asymptotics of the solution to the conductivity equation in the presence of adjacent circular inclusions with finite conductivities. Journal of Mathematical Analysis and Applications, 4211): ,