Volume integral equations of the scattering problem of poroelasticity and their properties

Scattering of monochromatic waves on an isolated inhomogeneity (inclusion) in an infinite poroelastic medium is considered. Wave propagation in the medium and the inclusion are described by Biot's equations of poroelasticity. The problem is reduced to 3D‐integro‐differential equations for displacement and pressure fields in the region occupied by the inclusion. Properties of the integral operators in these equations are studied. Discontinuities of the fields on the inclusion boundary are indicated. The case of a thin inclusion with low permeability is considered. The corresponding scattering problem is reduced to a 2D integral equation on the middle surface of the inclusion. The unknown function in this equation is the pressure jump in the transverse direction to the inclusion middle surface. An inclusion with a thin layer of low permeability on its interface is considered. The appropriate boundary conditions on the inclusion interface are pointed out. Methods of numerical solution of the volume integral equations of the scattering problems of poroelasticity are discussed.     

Method of Integral Equations for the Three-Dimensional Problem of Wave Reflection from an Irregular Surface

The problem on the reflection of the field of a plane H-polarized three-dimensional electromagnetic wave from a perfectly conducting interface between media which contains a local perfectly conducting inhomogeneity is considered. To construct a numerical algorithm, the boundary value problem for the system of Maxwell equations in an infinite domain with irregular boundary is reduced to a system of singular integral equations, which is solved by the approximation–collocation method. The elements of the resulting complex matrix are calculated by a specially developed algorithm. The solution of the system of singular integral equations is used to obtain an integral representation for the reflected electromagnetic field and computational formulas for the directional diagram of the reflected electromagnetic field in the far region.     

The operator equations of Lippmann–Schwinger type for acoustic and electromagnetic scattering problems in L 2

This article is concerned with the scattering of acoustic and electromagnetic time harmonic plane waves by an inhomogeneous medium. These problems can be translated into volume integral equations of the second kind – the most prominent example is the Lippmann–Schwinger integral equation. In this work, we study a particular class of scattering problems where the integral operator in the corresponding operator equation of Lippmann–Schwinger type fails to be compact. Such integral equations typically arise if the modelling of the inhomogeneous medium necessitates space-dependent coefficients in the highest order terms of the underlying partial differential equation. The two examples treated here are acoustic scattering from a medium with a space-dependent material density and electromagnetic medium scattering where both the electric permittivity and the magnetic permeability vary. In these cases, Riesz theory is not applicable for the solution of the arising integral equations of Lippmann–Schwinger type. Therefore, we show that positivity assumptions on the relative material parameters allow to prove positivity of the arising volume potentials in tailor-made weighted spaces of square integrable functions. This result merely holds for imaginary wavenumber and we exploit a compactness argument to conclude that the arising integral equations are of Fredholm type, even if the integral operators themselves are not compact. Finally, we explain how the solution of the integral equations in L ² affects the notion of a solution of the scattering problem and illustrate why the order of convergence of a Galerkin scheme set up in L ² does not suffer from our L ² setting, compared to schemes in higher order Sobolev spaces.     

The reissner-sagoci problem for a multilayer base with a cylindrical cavity

We study the problem of the torsional oscillations of a plane disk-shaped die coupled with the upper boundary of a multilayer elastic base containing a vertical cylindrical cavity whose axis is perpendicular to the interface of the layers. The problem is stated as paired integral equations connected with the Weber integral transforms. To couple the solutions in the layers we use the method of initial parameters, which makes it possible to express the stress-strain state in any layer in terms of the solution of a Fredholm integral equation of second kind, to which the paired equations reduce. We exhibit an algorithm for numerical implementation of the problem. Translated fromDinamicheskie Sistemy, No. 13, 1994, pp. 55–61.     

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.     

Fast iteration method in the problem of waves interacting with a set of thin screens

A new iteration method is proposed for the wave equation describing the scattering of a harmonic wave from an arbitrary configuration in the form of an array of thin straight barriers. The problem is reduced to a system of boundary integral equations, which are discretized by applying the Belotserkovskii-Lifanov method. In discrete form, a finite number of systems with Toeplitz matrices (the number of systems is equal to the number of barriers) are solved at each iteration step by applying special fast methods. The algorithm is tested on several geometries, and its convergence in these cases is analyzed.     

Scattering of the electromagnetic field on a finite impedance section of an interface

The scattering of a three-dimensional electromagnetic field on a local impedance section of a wavy surface is considered. The boundary-value problem for the system of Maxwell equations is reduced to solving a system of hypersingular integral equations. A numerical algorithm is developed and a program is designed for computing the electrodynamic characteristics of the given scattering problem. The accuracy of the simulation results is investigated. __________ Translated from Prikladnaya Matematika i Informatika, No. 25, pp. 56–69, 2007.     

Compact gradient tracking in shape optimization

In the present paper we consider the minimization of gradient tracking functionals defined on a compact and fixed subdomain of the domain of interest. The underlying state is assumed to satisfy a Poisson equation with Dirichlet boundary conditions. We proof that, in contrast to the situation of gradient tracking on the whole domain, the shape Hessian is not strictly H ^1/2-coercive at the optimal domain which implies ill-posedness of the shape problem under consideration. Shape functional and gradient require only knowledge of the Cauchy data of the state and its adjoint on the boundaries of the domain and the subdomain. These data can be computed by means of boundary integral equations when reformulating the underlying differential equations as transmission problems. Thanks to fast boundary element techniques, we derive an efficient algorithm to solve the problem under consideration.     

An ultra-weak method for acoustic fluid–solid interaction

We introduce the ultra-weak variational formulation (UWVF) for fluid–solid vibration problems. In particular, we consider the scattering of time-harmonic acoustic pressure waves from solid, elastic objects. The problem is modeled using a coupled system of the Helmholtz and Navier equations. The transmission conditions on the fluid–solid interface are represented in an impedance-type form after which we can employ the well known ultra-weak formulations for the Helmholtz and Navier equations. The UWVF approximation for both equations is computed using a superposition of propagating plane waves. A condition number based criterion is used to define the plane wave basis dimension for each element. As a model problem we investigate the scattering of sound from an infinite elastic cylinder immersed in a fluid. A comparison of the UWVF approximation with the analytical solution shows that the method provides a means for solving wave problems on relatively coarse meshes. However, particular care is needed when the method is used for problems at frequencies near the resonance frequencies of the fluid–solid system.     

Interaction of a hollow cylinder of finite length and a plate with a cylindrical cavity with a rigid insert

Two problems of the interaction of a hollow circular cylinder with load-free ends and an unbounded plate with a cylindrical cavity and a symmetrically imbedded rigid insert are considered. Homogeneous solutions are found and the generalized orthogonality of these solutions is used when the modified boundary conditions are satisfied. As a result, we have a system of two integral equations in functions of the displacements of the outer and inner surfaces of the hollow cylinder. These functions are sought in the form of sums of a trigonometric series and a power function with a root singularity. The ill-posed infinite systems of linear algebraic equations obtained are regularized by the introduction of small positive parameters. Since the elements of the matrices of the systems as well as the contact stresses are defined by poorly converging numerical and functional series, an efficient method for calculating of the remainders of the above-mentioned series is developed. Formulae are found for the contact pressure distribution function and the integral characteristic. Examples of the calculation of the interaction of the cylinder and the plate with an insert are given.The method of solving contact problems described here has been used earlier1, 2 and the generalized orthogonality of the solutions found for bodies of finite dimensions, that is, for a rectangle and cylinders of finite length, is its basis. Problems for hollow cylinders with a band ² and an insert reduce to a system of two integral equations, and the problem for a rectangle¹ reduces to one integral equation. Solving these integral equations, ill-posed systems of linear algebraic equations are obtained which are subject to regularization³.     

Solution procedures for three-dimensional eddy current problems

The problem under consideration is that of the scattering of time periodic electromagnetic fields by metallic obstacles. A common approximation here is that in which the metal is assumed to have infinite conductivity. The resulting problem, called the perfect conductor problem, involves solving Maxwell's equations in the region exterior to the obstacle with the tangential component of the electric field zero on the obstacle surface. In the interface problem different sets of Maxwell equations must be solved in the obstacle and outside while the tangential components of both electric and magnetic fields are continuous across the obstacle surface. Solution procedures for this problem are given. There is an exact integral equation procedure for the interface problem and an asymptotic procedure for large conductivity. Both are based on a new integral equation procedure for the perfect conductor problem. The asymptotic procedure gives an approximate solution by solving a sequence of problems analogous to the one for perfect conductors.     

Acoustic scattering by a thin cylindrical screen with Dirichlet and impedance boundary conditions on opposite sides of the screen

A problem on scattering acoustic waves by a thin cylindrical screen is studied. In doing so, the Dirichlet condition is specified on one side of the screen, while the impedance boundary condition is specified on the other side of the screen. The solution of the problem is subject to the radiating condition at infinity and to the propagative Helmholtz equation. By using the potential theory the scattering problem is reduced to a system of singular integral equations with additional conditions. By regularization and subsequent transformations, this system is reduced to a vector Fredholm equation of the second kind and index zero. It is proved that the obtained vector Fredholm equation is uniquely solvable. Therefore the integral representation for a solution of the original scattering problem is obtained. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Integral equation methods in electromagnetic scattering from anisotropic media

We investigate scattering of time‐harmonic electromagnetic waves by an anisotropic inhomogeneous medium. The problem is equivalently transformed into a system of strongly singular integral equations. The uniqueness and existence of a solution is shown and we examine the regularity of the solution by means of integral equations. We also prove the analyticity of the scattered field with respect to the refractive matrix and give a characterization of the derivatives in terms of solutions to anisotropic scattering problems. Copyright © 2000 John Wiley & Sons, Ltd.     

Integral equations in the wave scattering problem at an irregular interface of domains

Modeling the scattering of electromagnetic waves at an interface of media with different characteristics, one encounters the conjugation problem. Using the method of boundary integral equations and the theory of generalized potentials, we prove the classical resolvability of this problem. The boundary is assumed to be irregular. This means that the plane is divided into two domains by a curve which coincides with a straight line, except for a finite part, producing the irregularity. We propose algorithms for the approximate solution of the conjugation problem based on the spline methods for the solution of integral equations. We theoretically substantiate the computational scheme, namely, we prove the convergence and estimate the convergence rate.     

Scattering of Plane Harmonic Waves in an Infinite Elastic Solid by an Arbitrary Configuration of Two-dimensional Obstacles

The problem of the multiple scattering of plane harmonic wavesin an infinite elastic solid by arbitrary configurations oftwo-dimensional obstacles is considered. Sets of integral equationsrelating the multiple scattering functions to the single scatteringfunctions are derived. The integral equations are replaced bya set of algebraic equations which are then solved by an iterativeprocedure. Some results are given for the scattering of a planeharmonic P or SV wave by an identical pair of parallel circularcylindrical cavities.     

Shape Derivatives of Boundary Integral Operators in Electromagnetic Scattering. Part II: Application to Scattering by a Homogeneous Dielectric Obstacle

We develop the shape derivative analysis of solutions to the problem of scattering of time-harmonic electromagnetic waves by a penetrable bounded obstacle. Since boundary integral equations are a classical tool to solve electromagnetic scattering problems, we study the shape differentiability properties of the standard electromagnetic boundary integral operators. The latter are typically bounded on the space of tangential vector fields of mixed regularity T H^-\frac12(div_G,G){\mathsf T \mathsf H^{-\frac{1}{2}}({\rm div}_{\Gamma},\Gamma)}. Using Helmholtz decomposition, we can base their analysis on the study of pseudo-differential integral operators in standard Sobolev spaces, but we then have to study the Gateaux differentiability of surface differential operators. We prove that the electromagnetic boundary integral operators are infinitely differentiable without loss of regularity. We also give a characterization of the first shape derivative of the solution of the dielectric scattering problem as a solution of a new electromagnetic scattering problem.     

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

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

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.     

A low velocity approximation for the relativistic Vlasov-Maxwell system

An approximation for the relativistic Vlasov-Maxwell (RVM) system of partial differential equations in the one-space, two-momenta case is proposed. The speed of light, c, appears as a parameter in this system. The approximation is obtained by modifying certain integral operators appearing in integral representations, due to Glassey and Strauss, of the electric and magnetic fields, and replaces the hyperbolic Maxwell system with one that is elliptic in nature (for each fixed t). Solutions of the modified problem are shown to converge in a pointwise sense to solutions of (RVM) at the asymptotic rate of 1/c² as c tends to infinity.     

Existence of solution in elastic wave scattering by unbounded rough surfaces

We consider the two‐dimensional problem of the scattering of a time‐harmonic wave, propagating in an homogeneous, isotropic elastic medium, by a rough surface on which the displacement is assumed to vanish. This surface is assumed to be given as the graph of a function ?∈C^1,1(?). Following up on earlier work establishing uniqueness of solution to this problem, existence of solution is studied via the boundary integral equation method. This requires a novel approach to the study of solvability of integral equations on the real line. The paper establishes the existence of a unique solution to the boundary integral equation formulation in the space of bounded and continuous functions as well as in all L^p spaces, p∈[1, ∞] and hence existence of solution to the elastic wave scattering problem. Copyright © 2002 John Wiley & Sons, Ltd.