• The Neumann-Poincar\'e operator on the sphere has $\frac{1}{2(2k+1)}$, $k=0,1,2,\ldots$, as its eigenvalues and the corresponding multiplicity is $2k+1$. We consider the bifurcation of eigenvalues under deformation of domains, and show that Frech\'et derivative of the sum of the bifurcations is zero. We then discuss the connection of this result with some conjectures regarding the Neumann-Poincar\'e operator.
  • Upon insertion of an inclusion into a medium with the uniform field, if the field is not perturbed at all outside the inclusion, then it is called a neutral inclusion. It is called a weakly neutral inclusion if the field is perturbed weakly. The inclusions neutral to multiple uniform fields are of circular shape if the medium is isotropic, and any other shape cannot be neutral. We consider in this paper the problem of constructing inclusions of general shape which are weakly neutral to multiple fields in two dimensions. We show that a simply connected domain satisfying a certain geometric condition can be realized as a weakly neutral inclusion to multiple fields by introducing an imperfect interface parameter on the boundary. The geometric condition on the domain and the imperfect interface parameter are determined by the first coefficient of the conformal mapping from the exterior of the unit disk onto the exterior of the domain. We provide some numerical examples to compare field perturbations by weakly neutral inclusions and perfectly bonding interfaces. They clearly show that the field perturbation by weakly neutral inclusions is much weaker.
  • In the region between close-to-touching hard inclusions, the stress may be arbitrarily large as the inclusions get closer. The stress is represented by the gradient of a solution to the Lam\'e system of linear elasticity. We consider the problem of characterizing the gradient blow-up of the solution in the narrow region between two inclusions and estimating its magnitude. We introduce singular functions which are constructed in terms of nuclei of strain and hence are solutions of the Lam\'{e} system, and then show that the singular behavior of the gradient in the narrow region can be precisely captured by singular functions. As a consequence of the characterization, we are able to regain the existing upper bound on the blow-up rate of the gradient, namely, $\epsilon^{-1/2}$ where $\epsilon$ is the distance between two inclusions. We then show that it is in fact an optimal bound by showing that there are cases where $\epsilon^{-1/2}$ is also a lower bound. This work is the first to completely reveal the singular nature of the gradient blow-up in the context of the Lam\'{e} system with hard inclusions. The singular functions introduced in this paper play essential roles to overcome the difficulties in applying the methods of previous works. Main tools of this paper are the layer potential techniques and the variational principle. The variational principle can be applied because the singular functions of this paper are solutions of the Lam\'{e} system.
  • This paper deals with the field enhancement, that is, the gradient blow-up, due to presence of a bow-tie structure of perfectly conducting inclusions in two dimensions. The bow-tie structure consists of two disjoint bounded domains which have corners with possibly different aperture angles. The domains are parts of cones near the vertices, and they are nearly touching to each other. We characterize the field enhancement using explicit functions and, as consequences, derive optimal estimates of the gradient in terms of the distance between two inclusions and aperture angles of the corners. The estimates show that the field is enhanced beyond the corner singularities due to the interaction between two inclusions.
  • We prove in a mathematically rigorous way the asymptotic formula of Flaherty and Keller on the effective property of densely packed periodic elastic composites with hard inclusions. The proof is based on the primal-dual variational principle, where the upper bound is derived by using the Keller-type test functions and the lower bound by singular functions made of nuclei of strain. Singular functions are solutions of the Lam\'e system and capture precisely singular behavior of the stress in the narrow region between two adjacent hard inclusions.
  • We prove that the elastic Neumann--Poincar\'e operator defined on the smooth boundary of a bounded domain in three dimensions, which is known to be non-compact, is in fact polynomially compact. As a consequence, we prove that the spectrum of the elastic Neumann-Poincar\'e operator consists of three non-empty sequences of eigenvalues accumulating to certain numbers determined by Lam\'e parameters. These results are proved using the surface Riesz transform, calculus of pseudo-differential operators and the spectral mapping theorem.
  • We investigate anomalous localized resonance on the circular coated structure and cloaking related to it in the context of elasto-static systems. The structure consists of the circular core with constant Lam\'e parameters and the circular shell of negative Lam\'e parameters proportional to those of the core. We show that the eigenvalues of the Neumann-Poincar\'e operator corresponding to the structure converges to certain non-zero numbers determined by Lam\'e parameters and derive precise asymptotics of the convergence. We then show with estimates that cloaking by anomalous localized resonance takes place if and only if the dipole type source lies inside critical radii determined by the radii of the core and the shell.
  • We show that the eigenvalues of the Neumann-Poincar\'e operator on analytic boundaries of simply connected bounded planar domains tend to zero exponentially fast, and the exponential convergence rate is determined by the maximal Grauert radius of the boundary. We present a few examples of boundaries to show that the estimate is optimal.
  • We study spectral properties of the Neumann-Poincar\'e operator on planar domains with corners with particular emphasis on existence of continuous spectrum and pure point spectrum. We show that the rate of resonance at continuous spectrum is different from that at eigenvalues, and then derive a method to distinguish continuous spectrum from eigenvalues. We perform computational experiments using the method to see whether continuous spectrum and pure point spectrum appear on domains with corners. For the computations we use a modification of the Nystr\"om method which makes it possible to construct high-order convergent discretizations of the Neumann-Poincar\'e operator on domains with corners. The results of experiments show that all three possible spectra, absolutely continuous spectrum, singularly continuous spectrum, and pure point spectrum, may appear depending on domains. We also prove rigorously two properties of spectrum which are suggested by numerical experiments: symmetry of spectrum (including continuous spectrum), and existence of eigenvalues on rectangles of high aspect ratio.
  • The neutral inclusion problem in two dimensional isotropic elasticity is considered. The neutral inclusion, when inserted in a matrix having a uniform applied field, does not disturb the field outside the inclusion. The inclusion consists of the core and shell of arbitrary shapes, and their elasticity tensors are isotropic. We show that if the coated inclusion is neutral to a uniform bulk field, then the core and shell must be concentric disks, provided that the shear and bulk moduli satisfy certain conditions.
  • We consider well-posedness of the boundary value problem in presence of an inclusion with complex conductivity $k$. We first consider the transmission problem in $\mathbb{R}^d$ and characterize solvability of the problem in terms of the spectrum of the Neumann-Poincar\'e operator. We then deal with the boundary value problem and show that the solution is bounded in its $H^1$-norm uniformly in $k$ as long as $k$ is at some distance from a closed interval in the negative real axis. We then show with an estimate that the solution depends on $k$ in its $H^1$-norm Lipschitz continuously. We finally show that the boundary perturbation formula in presence of a diametrically small inclusion is valid uniformly in $k$ away from the closed interval mentioned before. The results for the single inclusion case are extended to the case when there are multiple inclusions with different complex conductivities: We first obtain a complete characterization of solvability when inclusions consist of two disjoint disks and then prove solvability and uniform estimates when imaginary parts of conductivities have the same signs. The results are obtained using the spectral property of the associated Neumann-Poincar\'e operator and the spectral resolution.
  • We first investigate spectral properties of the Neumann-Poincar\'e (NP) operator for the Lam\'e system of elasto-statics. We show that the elasto-static NP operator can be symmetrized in the same way as that for Laplace operator. We then show that even if elasto-static NP operator is not compact even on smooth domains, its spectrum consists of eigenvalues which accumulates to two numbers determined by Lam\'e constants. We then derive explicitly eigenvalues and eigenfunctions on disks and ellipses. We then investigate resonance occurring at eigenvalues and anomalous localized resonance at accumulation points of eigenvalues. We show on ellipses that cloaking by anomalous localized resonance takes place at accumulation points of eigenvalues.
  • We study resonance for the Helmholz equation with a finite frequency in a plasmonic material of negative dielectric constant in two and three dimensions. We show that the quasi-static approximation is valid for diametrically small inclusions. In fact, we quantitatively prove that if the diameter of a inclusion is small compared to the loss parameter, then resonance occurs exactly at eigenvalues of the Neumann-Poincar\'e operator associated with the inclusion.
  • The purpose of this paper is to propose a non-iterative method for the inverse conductivity problem of recovering multiple small anomalies from the boundary measurements. When small anomalies are buried in a conducting object, the electric potential values inside the object can be expressed by integrals of densities with a common sparse support on the location of anomalies. Based on this integral expression, we formulate the reconstruction problem of small anomalies as a joint sparse recovery and present an efficient non-iterative recovery algorithm of small anomalies. Furthermore, we also provide a slightly modified algorithm to reconstruct an extended anomaly. We validate the effectiveness of the proposed algorithm over the linearized method and the MUSIC algorithm by numerical simulations.
  • We study an elliptic equation with measurable coefficients arising from photo-acoustic imaging in inhomogeneous media. We establish Holder continuity of weak solutions and obtain pointwise bounds for Green's functions subject to Dirichlet or Neumann condition.
  • We consider the neutral inclusion problem in three dimensions which is to prove if a coated structure consisting of a core and a shell is neutral to all uniform fields, then the core and the shell must be concentric balls if the matrix is isotropic and confocal ellipsoids if the matrix is anisotropic. We first derive an over-determined boundary value problem in the shell of the neutral inclusion, and then prove in the isotropic case that if the over- determined problem admits a solution, then the core and the shell must be concentric balls. As a consequence it is proved that the structure is neutral to all uniform fields if and only if it consists of concentric balls provided that the coefficient of the core is larger than that of the shell.
  • The purpose of this paper is to investigate the spectral nature of the Neumann-Poincar\'e operator on the intersecting disks, which is a domain with the Lipschitz boundary. The complete spectral resolution of the operator is derived, which shows in particular that it admits only the absolutely continuous spectrum, no singularly continuous spectrum and no pure point spectrum. We then quantitatively analyze using the spectral resolution the plasmon resonance at the absolutely continuous spectrum.
  • We investigate in a quantitative way the plasmon resonance at eigenvalues and the essential spectrum (the accumulation point of eigenvalues) of the Neumann-Poincar\'e operator on smooth domains. We first extend the symmetrization principle so that the single layer potential becomes a unitary operator from $H^{-1/2}$ onto $H^{1/2}$. We then show that the resonance at the essential spectrum is weaker than that at eigenvalues. It is shown that anomalous localized resonance occurs at the essential spectrum on ellipses, but cloaking does not occur on ellipses unlike the core-shell structure considered in [20]. It is shown that resonance does not occur at the essential spectrum on three dimensional balls.
  • If stiff inclusions are closely located, then the stress, which is the gradient of the solution, may become arbitrarily large as the distance between two inclusions tends to zero. In this paper we investigate the asymptotic behavior of the stress concentration factor, which is the normalized magnitude of the stress concentration, as the distance between two inclusions tends to zero. For that purpose we show that the gradient of the solution to the case when two inclusions are touching decays exponentially fast near the touching point. We also prove a similar result when two inclusions are closely located and there is no potential difference on boundaries of two inclusions. We then use these facts to show that the stress concentration factor converges to a certain integral of the solution to the touching case as the distance between two inclusions tends to zero. We then present an efficient way to compute this integral.
  • We present a new systematic method to construct the conformal mapping from outside the unit disc to outside of a simply connected domain using the generalized polarization tensors. We also present some numerical results to validate effectiveness of the method.
  • The size estimation problem in electrical impedance tomography is considered when the conductivity is a complex number and the body is two-dimensional. Upper and lower bounds on the volume fraction of the unknown inclusion embedded in the body are derived in terms of two pairs of voltage and current data measured on the boundary of the body. These bounds are derived using the translation method. We also provide numerical examples to show that these bounds are quite tight and stable under measurement noise.
  • This paper extends the concept of generalized polarization tensors (GPTs), which was previously defined for inclusions with homogeneous conductivities, to inhomogeneous conductivity inclusions. We begin by giving two slightly different but equivalent definitions of the GPTs for inhomogeneous inclusions. We then show that, as in the homogeneous case, the GPTs are the basic building blocks for the far-field expansion of the voltage in the presence of the conductivity inclusion. Relating the GPTs to the Neumann-to-Dirichlet (NtD) map, it follows that the full knowledge of the GPTs allows unique determination of the conductivity distribution. Furthermore, we show important properties of the the GPTs, such as symmetry and positivity, and derive bounds satisfied by their harmonic sums. We also compute the sensitivity of the GPTs with respect to changes in the conductivity distribution and propose an algorithm for reconstructing conductivity distributions from their GPTs. This provides a new strategy for solving the highly nonlinear and ill-posed inverse conductivity problem. We demonstrate the viability of the proposed algorithm by preforming a sensitivity analysis and giving some numerical examples.
  • If a core of dielectric material is coated by a plasmonic structure of negative dielectric material with non-zero loss parameter, then anomalous localized resonance may occur as the loss parameter tends to zero and the source outside the structure can be cloaked. It has been proved that the cloaking due to anomalous localized resonance (CALR) takes place for structures of concentric disks and the critical radius inside which the sources are cloaked has been computed. In this paper, it is proved that CALR takes place for structures of confocal ellipses and the critical elliptic radii are computed. The method of this paper uses the spectral analysis of the Neumann-Poincar\'e type operator associated with two interfaces (the boundaries of the core and the shell).
  • When two perfectly conducting inclusions are located closely to each other, the electric field concentrates in a narrow region in between two inclusions, and becomes arbitrarily large as the distance between two inclusions tends to zero. The purpose of this paper is to derive an asymptotic formula of the concentration which completely characterizes the singular behavior of the electric field, when inclusions are balls of the same radii in three dimensions.
  • We consider the problem of neutral inclusions for two-dimensional conductivity and anti-plane elasticity. The neutral inclusion, when inserted in a matrix having a uniform field, does not disturb the field outside the inclusion. The inclusion consists of a core and a shell. We show that if the inclusion is neutral to two linearly independent fields, then the core and the shell are confocal ellipses.