Small perturbations of an interface for elastostatic problems

We consider the Lamé system for an elastic medium consisting of an inclusion embedded in a homogeneous background medium. Based on the field expansion method and layer potential techniques, we rigorously derived the asymptotic expansion of the perturbed displacement field because of small perturbations in the interface of the inclusion. We extend these techniques to determine a relationship between traction‐displacement measurements and the shape of the object and derive an asymptotic expansion for the perturbation in the elastic moment tensors because of the presence of small changes in the interface of the inclusion. Copyright © 2016 John Wiley & Sons, Ltd.     

Electromagnetic scattering from perturbed surfaces

This paper studies the scattering of electromagnetic waves from a (local) perturbation of a fixed surface, the boundary of a given obstacle in ?³. The goal is to produce an algorithm for solving boundary value problems in the exterior of the perturbed domain solely based on the knowledge of the Green function for the original surface. This is done by solving a boundary integral equation which only involves the perturbed portion of the boundary. Copyright © 2006 John Wiley & Sons, Ltd.     

An asymptotic expansion for perturbations in the displacement field due to the presence of thin interfaces

We rigorously derive an asymptotic expansion for two-dimensional displacement field associated with thin elastic inclusion having no uniform thickness. Our approach is based on layer potential techniques through integral representation formulas of the fields. We extend these techniques to determine a relationship between traction–displacement measurements and the shape of the thin inclusion.     

Asymptotic expansions for currents caused by small interface changes of an electromagnetic inclusion

We consider solutions to the Helmholtz equation in two dimensions. The aim of this article is to advance the development of high-order asymptotic expansions for boundary perturbations of currents caused by small perturbations of the shape of an inhomogeneity with 𝒞²-boundary. The work represents a natural completion of Ammari et al. [H. Ammari, H. Kang, M. Lim, and H. Zribi, Conductivity interface problems. Part I: Small perturbations of an interface, Trans. Am. Math. Soc. 363 (2010), pp. 2901–2922], where the solution for the Helmholtz equation is represented by a system and the proof of our asymptotic expansion is radically different from Ammari et al. (2010). Our derivation is rigorous and is based on the field expansion method. Its proof relies on layer potential techniques. It plays a key role in developing effective algorithms to determine certain properties of the shape of an inhomogeneity based on boundary measurements.     

The steady two-dimensional flow over a rectangular obstacle lying on the bottom

We study a plane problem with mixed boundary conditions for a harmonic function in an unbounded Lipschitz domain contained in a strip. The problem is obtained by linearizing the hydrodynamic equations which describe the steady flow of a heavy ideal fluid over an obstacle lying on the flat bottom of a channel. In the case of obstacles of rectangular shape we prove unique solvability for all velocities of the (unperturbed) flow above a critical value depending on the obstacle depth. We also discuss regularity and asymptotic properties of the solutions.     

The singular sources method for an inverse problem with mixed boundary conditions

We use the singular sources method to detect the shape of the obstacle in a mixed boundary value problem. The basic idea of the method is based on the singular behavior of the scattered field of the incident point-sources on the boundary of the obstacle. Moreover we take advantage of the scattered field estimate by the backprojection operator. Also we give a uniqueness proof for the shape reconstruction.     

On the rate of convergence of solutions in domain with periodic multilevel oscillating boundary

In this paper we deal with the homogenization problem for the Poisson equation in a singularly perturbed domain with multilevel periodically oscillating boundary. This domain consists of the body, a large number of thin cylinders joining to the body through the thin transmission zone with rapidly oscillating boundary. Inhomogeneous Fourier boundary conditions with perturbed coefficients are set on the boundaries of the thin cylinders and on the boundary of the transmission zone. We prove the homogenization theorems and derive the estimates for the convergence of the solutions. Copyright © 2010 John Wiley & Sons, Ltd.     

Continuous Frechet differentiability with respect to a Lipschitz domain and a stability estimate for direct acoustic scattering problems

We consider direct acoustic scattering problems with eithera sound-soft or sound-hard obstacle, or lossy boundary conditions,and establish continuous Fréchet differentiability withrespect to the shape of the scatterer of the scattered fieldand its corresponding far-field pattern. Our proof is basedon the implicit function theorem, and assumes that the boundaryof the scatterer as well as the deformation are only Lipschitzcontinuous. From continuous Fréchet differentiability,we deduce a stability estimate governing the variation of thefar-field pattern with respect to the shape of the scatterer.We illustrate this estimate with numerical results obtainedfor a two-dimensional high-frequency acoustic scattering problem.     

Error bounds of singular boundary method for potential problems

The singular boundary method (SBM) is a recent strong‐form boundary collocation method free of integration, mesh, and fictitious boundary. Although an extensive study has been reported in the literature on improving its accuracy and stability as well as its applications to diverse problems, little, however, has been done to analyze its convergence mathematically. The main purpose of this paper is to derive the explicit error bounds of the SBM for potential problems as well as to explain the essential difference between the origin intensity factor (OIF) in the SBM and the singular integration in the boundary element method (BEM). In the process of derivation, we also illustrate the physical meaning of OIF and explain the reason why the OIF has the function to correct the discretization error on the boundary. Finally, several benchmark examples are given to verify the effectiveness of the conclusions obtained from this article, as well as to investigate the different convergence behaviors between the SBM and BEM. It can be found that the SBM has the explicit error bound and is mathematically a stable technique.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1987–2004, 2017     

Variational approach to contact problems in nonlinear elasticity

We use variational methods to study problems in nonlinear 3-dimensional elasticity where the deformation of the elastic body is restricted by a rigid obstacle. For an assigned variational problem we first verify the existence of constrained minimizers whereby we extend previous results. Then we rigorously derive the Euler-Lagrange equation as necessary condition for minimizers, which was possible before only under strong smoothness assumptions on the solution. The Lagrange multiplier corresponding to the obstacle constraint provides structural information about the nature of frictionless contact. In the case of contact with, e.g., a corner of the obstacle, we derive a qualitatively new contact condition taking into account the deformed shape of the elastic body. By our analysis it is shown here for the first time rigorously that energy minimizers really solve the mechanical contact problem. Received: 20 October 2000 / Accepted: 7 June 2001 / Published online: 5 September 2002     

Expansions of Dirichlet to Neumann operators under perturbations

We obtain an asymptotic expansion of the Dirichlet to Neumann operator (DNO) for the Dirichlet problem on perturbations of the unit disk. We write our result in terms of pseudodifferential operators which themselves have expansions in the perturbation parameter. For a given power of the perturbation parameter, m > 0, and a given order, n < 0, we give an algorithm which allows for the expansion of the symbol of the DNO up to mth power in the perturbation parameter, with error terms belonging to symbols of order n.     

The inverse scattering problem for partially penetrable obstacles

We consider the direct and inverse problems for the scattering of a partially penetrable obstacle. Here ‘partially penetrable obstacle’ means that the waves transmit into the obstacle just from partial boundary of the obstacle with the rest of the boundary touching a known perfect and thin scatterer. The solvability of the direct scattering problem is presented using the classical boundary integral equation method. An interesting interior transmission problem is investigated for the purpose of solving the inverse obstacle scattering problem. Then the linear sampling method is proposed to reconstruct the shape and location of the obstacle from near field measurements. We note that the inversion algorithm can be implemented by avoiding the use of background Green function as a test function due to a mixed reciprocal principle.     

On imaging obstacles inside inhomogeneous media

We show that an obstacle inside a known inhomogeneous medium can be determined from measurements of the scattering amplitude at one frequency, without a priori knowledge of the boundary condition. We also show that an obstacle inside a known inhomogeneous anisotropic conducting medium can be determined from electrostatic current and voltage measurements on the boundary of a domain containing the obstacle. Moreover, two obstacles with boundary measurements which are merely comparable as operators must be identical. The first part of the paper gives an extension of the factorization method which may be of independent interest and also yields a new reconstruction procedure.     

A technique to prove parameter-uniform convergence for a singularly perturbed convection–diffusion equation

A priori parameter explicit bounds on the solution of singularly perturbed elliptic problems of convection–diffusion type are established. Regular exponential boundary layers can appear in the solution. These bounds on the solutions and its derivatives are obtained using a suitable decomposition of the solution into regular and layer components. By introducing extensions of the coefficients to a larger domain, artificial compatibility conditions are not imposed in the derivation of these decompositions.     

Stabilization of the trial method for the Bernoulli problem in case of prescribed Dirichlet data

We apply the trial method for the solution of Bernoulli's free boundary problem when the Dirichlet boundary condition is imposed for the solution of the underlying Laplace equation, and the free boundary is updated according to the Neumann boundary condition. The Dirichlet boundary value problem for the Laplacian is solved by an exponentially convergent boundary element method. The update rule for the free boundary is derived from the linearization of the Neumann data around the actual free boundary. With the help of shape sensitivity analysis and Banach's fixed‐point theorem, we shed light on the convergence of the respective trial method. Especially, we derive a stabilized version of this trial method. Numerical examples validate the theoretical findings.Copyright © 2014 John Wiley & Sons, Ltd.     

Detecting inclusions with a generalized impedance condition from electrostatic data via sampling

In this paper, we derive a sampling method to solve the inverse shape problem of recovering an inclusion with a generalized impedance condition from electrostatic Cauchy data. The generalized impedance condition is a second order differential operator applied to the boundary of the inclusion. We assume that the Dirichlet‐to‐Neumann mapping is given from measuring the current on the outer boundary from an imposed voltage. A simple numerical example is given to show the effectiveness of the proposed inversion method for recovering the inclusion. We also consider the inverse impedance problem of determining the impedance parameters for a known material from the Dirichlet‐to‐Neumann mapping assuming the inclusion has been reconstructed where uniqueness for the reconstruction of the coefficients is proven.     

Shape calculus and finite element method in smooth domains

The use of finite elements in smooth domains leads naturally to polyhedral or piecewise polynomial approximations of the boundary. Hence the approximation error consists of two parts: the geometric part and the finite element part. We propose to exploit this decomposition in the error analysis by introducing an auxiliary problem defined in a polygonal domain approximating the original smooth domain. The finite element part of the error can be treated in the standard way. To estimate the geometric part of the error, we need quantitative estimates related to perturbation of the geometry. We derive such estimates using the techniques developed for shape sensitivity analysis.     

Decomposition in regular and singular parts (in domains with corners) and stability under perturbation of the geometry

We analyze the Dirichlet problem for the Laplacian in a polygonal domain where boundary and angles depend on a parameter. We use the boundary integral equation, localization and Mellin transformation techniques to show that the solution has a decomposition in regular and singular parts which blow up at certain exceptional angles. We derive a modified decomposition which depends continuously on the angle.     

Perturbed integral equation method on determining unknown geometry in fluid flow

We consider an inverse problem arising in fluid flow. An algorithm to find the shape of a body in uniform flow is proposed when the tangential velocity on its boundary is given a priori. The fluid flow is assumed to be inviscid, incompressible and irrotational.The essential idea to develop our algorithm is the boundary modification process toward the solution shape with the help of the perturbed integral equations. The perturbed integral equations are derived from the boundary perturbation. We also give examples exhibiting the reliability for our proposed algorithm.