A remark on the single scattering preconditioner applied to boundary integral equations

This article deals with boundary integral equation preconditioning for the multiple scattering problem. The focus is put on the single scattering preconditioner, corresponding to the diagonal part of the integral operator, for which two results are proved. Indeed, after applying this geometric preconditioner, it appears that, firstly, every direct integral equations become identical to each other, and secondly, that the indirect integral equation of Brakhage–Werner becomes equal to the direct integral equations, up to a change of basis. These properties imply in particular that the convergence rate of a Krylov subspaces solver will be exactly the same for every preconditioned integral equations. To illustrate this, some numerical simulations are provided at the end of the paper.     

On radiation conditions for rough surface scattering problems

^** Email: arens{at}numathics.com^*** Email: hohage{at}math.uni-goettingen.de It is well known that Sommerfeld's radiation condition is nota valid characterization of outgoing waves for scattering problemsat rough surfaces. Instead, a radiation condition called upwardpropagating radiation condition (UPRC) is commonly used. Recently,a different radiation condition called the pole condition hasbeen investigated for scattering problems at bounded obstacles.In this paper we show the equivalence between the UPRC and thepole condition. In doing so, we give a rigorous interpretationof a formula called the angular spectrum representation forDirichlet data in the space of bounded continuous functions.     

Water wave scattering by finite arrays of circular structures

The scattering of small amplitude water waves by a finite arrayof locally axisymmetric structures is considered. Regions ofvarying quiescent depth are included and their axisymmetricnature, together with a mild-slope approximation, permits anadaptation of well-known interaction theory which ultimatelyreduces the problem to a simple numerical calculation. Numericalresults are given and effects due to regions of varying depthon wave loading and free-surface elevation are presented.     

求解二维多区域声波散射问题的边界元方法

对于多散射区域的声波散射问题的外Neumann边值问题,用单层位势来逼近每个散射域上的散射波,再利用位势理论的跳跃关系将问题转换为第二类边界积分方程组的求解问题,然后用Nystrom方法进行了求解.对多个随机散射区域的声波散射问题,数值例子体现了该求解方法的可行性和准确性.     

Boundary layer techniques for solving the Helmholtz equation in the presence of small inhomogeneities

We consider solutions to the Helmholtz equation in two and three dimensions. Based on layer potential techniques we provide for such solutions a rigorous systematic derivation of complete asymptotic expansions of perturbations resulting from the presence of diametrically small inhomogeneities with constitutive parameters different from those of the background medium. It is expected that our results will find important applications for developing effective algorithms for reconstructing small dielectric inhomogeneities from boundary measurements.     

Convergence in competition models with small diffusion coefficients

It is well known that for reaction-diffusion 2-species Lotka-Volterra competition models with spatially independent reaction terms, global stability of an equilibrium for the reaction system implies global stability for the reaction-diffusion system. This is not in general true for spatially inhomogeneous models. We show here that for an important range of such models, for small enough diffusion coefficients, global convergence to an equilibrium holds for the reaction-diffusion system, if for each point in space the reaction system has a globally attracting hyperbolic equilibrium. This work is planned as an initial step towards understanding the connection between the asymptotics of reaction-diffusion systems with small diffusion coefficients and that of the corresponding reaction systems.     

Decay and scattering of small solutions for Rosenau equations

We study the long-time behavior of small solutions of the Cauchy problem for a Rosenau equation. For a class of nonlinearity of the perturbation, the global small solution was obtained, and the decay and scattering for small amplitude solution are established.     

Analysis of a finite PML approximation for the three dimensional time-harmonic Maxwell and acoustic scattering problems

We consider the approximation of the frequency domain three-dimensional Maxwell scattering problem using a truncated domain perfectly matched layer (PML). We also treat the time-harmonic PML approximation to the acoustic scattering problem. Following work of Lassas and Somersalo in 1998, a transitional layer based on spherical geometry is defined, which results in a constant coefficient problem outside the transition. A truncated (computational) domain is then defined, which covers the transition region. The truncated domain need only have a minimally smooth outer boundary (e.g., Lipschitz continuous). We consider the truncated PML problem which results when a perfectly conducting boundary condition is imposed on the outer boundary of the truncated domain. The existence and uniqueness of solutions to the truncated PML problem will be shown provided that the truncated domain is sufficiently large, e.g., contains a sphere of radius . We also show exponential (in the parameter ) convergence of the truncated PML solution to the solution of the original scattering problem inside the transition layer.

Our results are important in that they are the first to show that the truncated PML problem can be posed on a domain with nonsmooth outer boundary. This allows the use of approximation based on polygonal meshes. In addition, even though the transition coefficients depend on spherical geometry, they can be made arbitrarily smooth and hence the resulting problems are amenable to numerical quadrature. Approximation schemes based on our analysis are the focus of future research.

     

The Kirsch-Kress method for inverse scattering by infinite locally rough interfaces

This paper is concerned with the inverse problem of reconstructing an infinite, locally rough interface from the scattered field measured on line segments above and below the interface in two dimensions. We extend the Kirsch-Kress method originally developed for inverse obstacle scattering problems to the above inverse transmission problem with unbounded interfaces. To this end, we reformulate our inverse problem as a nonlinear optimization problem with a Tikhonov regularization term. We prove the convergence of the optimization problem when the regularization parameter tends to zero. Finally, numerical experiments are carried out to show the validity of the inversion algorithm.     

High frequency scattering by convex curvilinear polygons

We consider the scattering of a time-harmonic acoustic incident plane wave by a sound soft convex curvilinear polygon with Lipschitz boundary. For standard boundary or finite element methods, with a piecewise polynomial approximation space, the number of degrees of freedom required to achieve a prescribed level of accuracy grows at least linearly with respect to the frequency of the incident wave. Here we propose a novel Galerkin boundary element method with a hybrid approximation space, consisting of the products of plane wave basis functions with piecewise polynomials supported on several overlapping meshes; a uniform mesh on illuminated sides, and graded meshes refined towards the corners of the polygon on illuminated and shadow sides. Numerical experiments suggest that the number of degrees of freedom required to achieve a prescribed level of accuracy need only grow logarithmically as the frequency of the incident wave increases.     

A stable well-conditioned integral equation for electromagnetism scattering

Finding a formulation for electromagnetic scattering of surfaces which is both well-posed and produces a well-conditioned linear system is still a challenging problem. We here propose one such formulation valid in the high-frequency regime. The mathematical analysis is provided and numerical results on rather complex geometries show the performance of the method.     

Scattering of Acoustoelectric Waves on a Cylindrical Inhomogeneity in a Transversely Isotropic Piezoelectric Medium

Scattering of acoustoelectric waves on a inhomogeneity is studied. The scatterer is a circular, continuous piezoelectric cylinder (fiber) embedded in a transversely isotropic piezoelectric medium. Expressions are found for the scattering amplitudes and total cross sections of three acoustoelectric waves propagating in the direction normal to the fiber axis. In the long-wave approximation, these expressions are obtained in an explicit form.     

Asymptotics of optimal control in the problem of harmonic wave scattering by an obstacle

Optimal impedance control for the Helmholtz equation in an unbounded domain is studied. Asymptotics of the optimal control with respect to a regularization parameter are constructed.     

Some norm estimates for the solutions of a nonlinear diffusion problem

Summary By means of Rellich's identity bounds for the spectrum of the nonlinear problem v+e ^v=0 are derived and certain norms for the solutions are estimated.

---

Zusammenfassung Mit Hilfe einer Rellichschen Identität werden Schranken für das Spektrum des nichtlinearen Problems v+e ^v=0 angegeben. Ferner werden gewisse Normen für die Lösungen abgeschätzt.

     

A Second-Order Method for the Electromagnetic Scattering from a Large Cavity

In this paper,we study the electromagnetic scattering from a two dimen- sional large rectangular open cavity embedded in an infinite ground plane,which is modelled by Helmholtz equations.By introducing nonlocal transparent boundary con- ditions,the problem in the open cavity is reduced to a bounded domain problem.A hypersingular integral operator and a weakly singular integral operator are involved in the TM and TE cases,respectively.A new second-order Toeplitz type approximation and a second-order finite difference scheme are proposed for approximating the hyper- singular integral operator on the aperture and the Helmholtz in the cavity,respectively. The existence and uniqueness of the numerical solution in the TE case are established for arbitrary wavenumbers.A fast algorithm for the second-order approximation is pro- posed for solving the cavity model with layered media.Numerical results show the second-order accuracy and efficiency of the fast algorithm.More important is that the algorithm is easy to implement as a preconditioner for cavity models with more general media.     

Energy norm a posteriori error estimates for mixed finite element methods

This paper deals with the a posteriori error analysis of mixed finite element methods for second order elliptic equations. It is shown that a reliable and efficient error estimator can be constructed using a postprocessed solution of the method. The analysis is performed in two different ways: under a saturation assumption and using a Helmholtz decomposition for vector fields.

     

The method of boundary integral equations for a mixed scattering problem

In this paper we consider the scattering of an electromagnetic time-harmonic plane wave by an infinite cylinder having an open arc and a bounded domain in R² as cross section. To this end, we solve a scattering problem for the Helmholtz equation in R² where the scattering object is a combination of a crack Γ and a bounded obstacle D, and we have Dirichlet-impedance type boundary condition on Γ and Dirichlet boundary condition on ∂D (∂D∈C²). Applying potential theory, the problem can be reformulated as a boundary integral system. We establish the existence and uniqueness of a solution to the system by using the Fredholm theory.     

Inverse problem of acoustic scattering in a space with a local inhomogeneity

An inclusion of finite size with a variable wave propagation velocity is contained in a homogeneous space. It is exposed to plane waves being propagated in all possible directions. The inverse problem is to restore the velocity by means of the scattering amplitude which is determined in terms of the scattered wave asymptotic for large x. A procedure for restoration is described and a uniqueness theorem is presented in the paper.Translated from Zapiski Nauchnykh Seminarov Leningradskogo Otdeleniya Matematicheskogo Instituta im. V. A. Steklova AN SSSR, Vol. 156, pp. 24–34, 1986.The authors are deeply grateful to the participants of the A. S. Blagoveschenskii seminar for valuable comments and discussion.     

Finite plane deformations of a three-phase circular inhomogeneity-matrix system

We consider finite plane deformations of a three-phase circular inhomogeneity-matrix system in which the inhomogeneity, the interphase layer and the matrix belong to the same class of compressible hyperelastic materials of harmonic-type but with each phase possessing its own distinct material properties. We obtain the complete solution when the system is subjected to general classes of remote (Piola) stress, specifically, remote stress distributions characterized by stress functions described by general polynomials of order n?1 in the corresponding complex variable z used to describe the matrix. As a particular case of the aforementioned analysis, we establish an Eshelby-type result namely that, for this class of harmonic materials, a three-phase circular inhomogeneity under uniform remote stress and eigenstrain, admits an internal uniform stress field when subjected to plane deformations.     

Closed-form solutions for a mode III radial matrix crack penetrating a circular inhomogeneity

In this research we address in detail a mode III radial matrix crack penetrating a circular inhomogeneity. One tip of the radial crack lies in the matrix, while the other tip of the radial crack lies in the circular inhomogeneity. In addition the two tips of the crack are mutually image points (or inverse points) with respect to the circular inhomogeneity-matrix interface. First we conformally map the crack onto a unit circle C_a in the new ζ-plane. Meanwhile the inhomogeneity-matrix interface is mapped onto C_b, a part of another circle in the ζ-plane. In addition C_a and C_b intersect at a vertex angle π/2. By using the method of image in the ζ-plane, closed-form solutions in terms of elementary functions are derived for three loading cases: (1) remote uniform antiplane shearing; (2) a screw dislocation located in the unbounded matrix; and (3) a radial Zener–Stroh crack.