A boundary integral equation method for the transmission eigenvalue problem

We propose a new integral equation formulation to characterize and compute transmission eigenvalues for constant refractive index that play an important role in inverse scattering problems for penetrable media. As opposed to the recently developed approach by Cossonnière and Haddar [1,2] which relies on a two by two system of boundary integral equations our analysis is based on only one integral equation in terms of Dirichlet-to-Neumann or Robin-to-Dirichlet operators which results in a noticeable reduction of computational costs. We establish Fredholm properties of the integral operators and their analytic dependence on the wave number. Further we employ the numerical algorithm for analytic non-linear eigenvalue problems that was recently proposed by Beyn [3] for the numerical computation of transmission eigenvalues via this new integral equation.     

Interior transmission eigenvalue problem with refractive index having C2-transition to the background medium

The interior transmission eigenvalue problem for scalar acoustics is studied for a new class of refractive index. Existence of an infinite discrete set of transmission eigenvalues in the case that the acoustic properties of a domain D???? ⁿ are allowed to have a C ²-transition to the homogeneous background medium is established. It is shown that the transmission problem has a weak formulation on certain weighted Sobolev spaces for this class of refractive index. The weak formulation and the discreteness of the spectrum is justified by using the Hardy inequality to prove compact imbedding theorems. Existence of transmission eigenvalues is demonstrated by investigating a generalized eigenvalue problem associated with the weak formulation.     

Boundary integral equations for bending of thermoelastic plates with transmission boundary conditions

The solution of an initial‐boundary value problem for bending of a piecewise‐homogeneous thermoelastic plate with transverse shear deformation is represented as various combinations of single‐layer and double‐layer time‐dependent potentials. The unique solvability of the boundary integral equations generated by these representations is proved in spaces of distributions. Copyright © 2009 John Wiley & Sons, Ltd.     

A system of boundary integral equations for the transmission problem in acoustics

In this paper, we reduce the classical two-dimensional transmission problem in acoustic scattering to a system of coupled boundary integral equations (BIEs), and consider the weak formulation of the resulting equations. Uniqueness and existence results for the weak solution of corresponding variational equations are established. In contrast to the coupled system in Costabel and Stephan (1985) [4], we need to take into account exceptional frequencies to obtain the unique solvability. Boundary element methods (BEM) based on both the standard and a two-level fast multipole Galerkin schemes are employed to compute the solution of the variational equation. Numerical results are presented to verify the efficiency and accuracy of the numerical methods.     

On the existence of transmission eigenvalues in an inhomogeneous medium

We prove the existence of transmission eigenvalues corresponding to the inverse scattering problem for isotropic and anisotropic media for both the scalar problem and Maxwell's equations. Considering a generalized abstract eigenvalue problem, we are able to extend the ideas of Päivärinta and Sylvester [Transmission eigenvalues, SIAM J. Math. Anal. 40, (2008) pp. 783–753] to prove the existence of transmission eigenvalues for a larger class of interior transmission problems. Our analysis includes both the case of a medium with positive contrast and of a medium with negative contrast provided that the contrasts are large enough.     

Integral dispersion equation method to solve a nonlinear boundary eigenvalue problem

In this work a nonlinear eigenvalue problem for a nonlinear autonomous ordinary differential equation of the second order is considered. This problem describes the process of propagation of transverse-electric electromagnetic waves along a plane dielectric waveguide with nonlinear permittivity. We demonstrate, as far as we know, a new method that allows one to derive an equation w.r.t. spectral parameter (the dispersion equation) which contains all necessary information about the eigenvalues. The method is based on a simple idea that the distance between zeros of a periodic solution to the differential equation is the same for the adjacent zeros. This method has no connections with the perturbation theory or the notion of a bifurcation point. Theorem of equivalence between the eigenvalue problem and the dispersion equation is proved. Periodicity of the eigenfunctions is proved, a formula for the period is found, and zeros of the eigenfunctions are determined. The formula for the distance between adjacent zeros of any eigenfunction is given. Also theorems of existence and localization of the eigenvalues are proved.     

A conformal mapping algorithm for the Bernoulli free boundary value problem

We propose a new numerical method for the solution of the Bernoulli free boundary value problem for harmonic functions in a doubly connected domain D in where an unknown free boundary Γ₀ is determined by prescribed Cauchy data on Γ₀ in addition to a Dirichlet condition on the known boundary Γ₁. Our main idea is to involve the conformal mapping method as proposed and analyzed by Akduman, Haddar, and Kress for the solution of a related inverse boundary value problem. For this, we interpret the free boundary Γ₀ as the unknown boundary in the inverse problem to construct Γ₀ from the Dirichlet condition on Γ₀ and Cauchy data on the known boundary Γ₁. Our method for the Bernoulli problem iterates on the missing normal derivative on Γ₁ by alternating between the application of the conformal mapping method for the inverse problem and solving a mixed Dirichlet–Neumann boundary value problem in D. We present the mathematical foundations of our algorithm and prove a convergence result. Some numerical examples will serve as proof of concept of our approach. Copyright © 2015 John Wiley & Sons, Ltd.     

Symmetric boundary integral formulations for Helmholtz transmission problems

In this work we propose and analyze numerical methods for the approximation of the solution of Helmholtz transmission problems in two or three dimensions. This kind of problems arises in many applications related to scattering of acoustic, thermal and electromagnetic waves. Formulations based on boundary integral methods are powerful tools to deal with transmission problems in unbounded media. Different formulations using boundary integral equations can be found in the literature. We propose here new symmetric formulations based on a paper by Martin Costabel and Ernst P. Stephan (1985), that uses the Calderón projector for the interior and exterior problems to develop closed expressions for the interior and exterior Dirichlet-to-Neumann operators. These operators are then matched to obtain an integral system that is equivalent to the Helmholtz transmission problem and uses Cauchy data on the transmission boundary as unknowns. We show how to simplify the aspect and analysis of the method by employing an additional mortar unknown with respect to the ones used in the original paper, writing it in an appropriate way to devise Krylov type iterations based on the separate Dirichlet-to-Neumann operators.     

Some results on electromagnetic transmission eigenvalues

The electromagnetic interior transmission problem is a boundary value problem, which is neither elliptic nor self-adjoint. The associated transmission eigenvalue problem has important applications in the inverse electromagnetic scattering theory for inhomogeneous media. In this paper, we show that, in general, there do not exist purely imaginary electromagnetic transmission eigenvalues. For constant index of refraction, we prove that it is uniquely determined by the smallest (real) transmission eigenvalue. Finally, we show that complex transmission eigenvalues must lie in a certain region in the complex plane. The result is verified by examples. Copyright © 2014 John Wiley & Sons, Ltd.     

Combined boundary integral equations for acoustic scattering‐resonance problems

In this paper, boundary integral formulations for a time‐harmonic acoustic scattering‐resonance problem are analyzed. The eigenvalues of eigenvalue problems resulting from boundary integral formulations for scattering‐resonance problems split in general into two parts. One part consists of scattering‐resonances, and the other one corresponds to eigenvalues of some Laplacian eigenvalue problem for the interior of the scatterer. The proposed combined boundary integral formulations enable a better separation of the unwanted spectrum from the scattering‐resonances, which allows in practical computations a reliable and simple identification of the scattering‐resonances in particular for non‐convex domains. The convergence of conforming Galerkin boundary element approximations for the combined boundary integral formulations of the resonance problem is shown in canonical trace spaces. Numerical experiments confirm the theoretical results. Copyright © 2016 John Wiley & Sons, Ltd.     

The inverse eigenvalue problem for a weighted helmholtz equation

Given the m lowest eigenvalues, we seek to recover an approximation to the density function ρ in the weighted Helmholtz equation -Δ=λρu on a rectangle with Dirchlet boundary conditions. The density ρ is assumed to be symmetric with respect to the midlines of the rectangle. Projection of the boundary value problem and the unknown density function onto appropriate vector spaces leads to a matrix inverse problem. Solutions of the matrix inverse problem exist provided that the reciprocals of the prescribed eigenvalues are close to the reciprocals of the simple eigenvalues of the base problem with ρ = 1. The matrix inverse problem is solved by a fixed—point iterative method and a density function ρ^* is constructed which has the same m lowest eigenvalues as the unknown ρ. The algorithm can be modified when multiple base eigenvalues arise, although the success of the modification depends on the symmetry properties of the base eigenfunctions.     

The method of integral equation for scattering problem with a mixed crack

We consider the scattering of an electromagnetic time‐harmonic plane wave by an infinite cylinder having a mixed open crack (or arc) in R² as the cross section. The crack is made up of two parts, and one of the two parts is (possibly) coated by a material with surface impedance λ. We transform the scattering problem into a system of boundary integral equations by adopting a potential approach, and establish the existence and uniqueness of a weak solution to the system by the Fredholm theory. Copyright © 2011 John Wiley & Sons, Ltd.     

Inverse scattering for planar cracks via nonlinear integral equations

We present a Newton‐type method for reconstructing planar sound‐soft or perfectly conducting cracks from far‐field measurements for one time‐harmonic scattering with plane wave incidence. Our approach arises from a method suggested by Kress and Rundell (Inv. Probl. 2005; 21 (4):1207–1223) for an inverse boundary value problem for the Laplace equation. It was extended to inverse scattering problems for sound‐soft obstacles (Mathematical Methods in Scattering Theory and Biomedical Engineering. World Scientific: Singapore, 2006; 39–50) and for sound‐hard cracks (Inv. Probl. 2006; 22 (6)). In both cases it was shown that the method gives accurate reconstructions with reasonable stability against noisy data. The approach is based on a pair of nonlinear and ill‐posed integral equations for the unknown boundary. The integral equations are solved by linearization, i.e. by regularized Newton iterations. Numerical reconstructions illustrate the feasibility of the method. Copyright © 2007 John Wiley & Sons, Ltd.     

Modified boundary integral formulations for the Helmholtz equation

In this paper we describe some modified regularized boundary integral equations to solve the exterior boundary value problem for the Helmholtz equation with either Dirichlet or Neumann boundary conditions. We formulate combined boundary integral equations which are uniquely solvable for all wave numbers even for Lipschitz boundaries Γ=∂Ω. This approach extends and unifies existing regularized combined boundary integral formulations.     

A boundary integral formulation of the Stefan problem

This paper considers a boundary integral formulation of the Stefan problem for two spatial dimensions. This formulation has the advantage that its numerical implementation does not require the discretization of the Stefan condition. Furthermore, the formulation is capable of solving problems with complex boundaries. Several illustrative examples are given.     

An inverse coefficient problem for a quasilinear parabolic equation with periodic boundary and integral overdetermination condition

In this paper, the inverse problem of finding the time‐dependent coefficient of heat capacity together with the solution periodic boundary and integral overdetermination conditions is considered. Under some natural regularity and consistency conditions on the input data, the existence, uniqueness, and continuous dependence upon the data of the solution are shown. Some considerations on the numerical solution for this inverse problem are presented with an example. Copyright © 2014 John Wiley & Sons, Ltd.     

A convolution quadrature Galerkin boundary element method for the exterior Neumann problem of the wave equation

The numerical solution of the Neumann problem of the wave equation on unbounded three‐dimensional domains is calculated using the convolution quadrature method for the time discretization and a Galerkin boundary element method for the spatial discretization. The mathematical analysis that has been built up for the Dirichlet problem is extended and developed for the Neumann problem, which is important for many modelling applications. Numerical examples are then presented for one of these applications, modelling transient acoustic radiation. Copyright © 2009 John Wiley & Sons, Ltd.     

A fast numerical method for a natural boundary integral equation for the Helmholtz equation

A Neumann boundary value problem of the Helmholtz equation in the exterior circular domain is reduced into an equivalent natural boundary integral equation. Using our trigonometric wavelets and the Galerkin method, the obtained stiffness matrix is symmetrical and circulant, which lead us to a fast numerical method based on fast Fourier transform. Furthermore, we do not need to compute the entries of the stiffness matrix. Especially, our method is also efficient when the wave number k in the Helmholtz equation is very large.