Local solutions to high-frequency 2D scattering problems

We consider the solution of high-frequency scattering problems in two dimensions, modeled by an integral equation on the boundary of a smooth scattering object. We devise a numerical method to obtain solutions on only parts of the boundary with little computational effort. The method incorporates asymptotic properties of the solution and can therefore attain particularly good results for high frequencies. We show that the integral equation in this approach reduces to an ordinary differential equation.     

Transient scattering by a circular cylinder

The scattering of transient plane waves by a circular cylinder is studied by using the Kirchhoff time-retarded potential boundary integral equation method. Two distinct problems are solved: (i) surface velocity potentials (or pressures) are found for rigid cylinders scattering ramp, ramp-step and Gaussian incident potential (or pressure) waves and (ii) surface velocities are found for free boundary (pressure release) cylinders scattering ramp, ramp-step and Gaussian incident velocity waves. The numerical schemes for both boundary conditions are very similar; since the same influence coefficients are used they differ only by a sign in the final formulation. Numerical results are readily obtained for the first few transit times. This approach is complementary to the usual modal approach in that it is best suited to early time values where the modal solutions converge most slowly.     

开口/封闭薄壳体声辐射和散射的统一边界积分方程解法

推导了在二维和三维空间下开口和封闭薄壳体在任意阻抗边界条件下声辐射和散射的统一边界积分方程.相对于以前的求解方法,该方程求解声辐射和散射问题具有相同的影响矩阵,能够同时求解薄壳体气动和振动噪声的辐射和散射现象,以及分析壳体声阻抗对声波传播的影响.推导的方程可以应用于叶轮机械、管道等噪声和消声器消声性能的预测等方面.在此方程基础上,可以进一步考虑运动边界和运动介质对声辐射和散射的影响. 关键词： 薄壳体 声阻抗 积分方程 边界元方法     

Unified boundary integral equation for the scattering of elastic and acoustic waves: solution by the method of moments

A unified boundary integral equation (BIE) is developed for the scattering of elastic and acoustic waves. Traditionally, the elastic and acoustic wave problems are solved separately with different BIEs. The elastic wave case is represented in a vector BIE with the traction and displacement vectors as unknowns whereas the acoustic wave case is governed by a scalar BIE with velocity potential or pressure as unknowns. Although these two waves can be unified in the form of a partial differential equation, the unified form in its BIE counterpart has not been reported. In this work, we derive the unified BIE for these two waves and then show that the acoustic wave case can be derived from this BIE by introducing a shielding loss for small shear modulus approximation; hence only one code needs to be maintained for both elastic and acoustic wave scattering. We also derive the asymptotic Green's tensor for zero shear modulus and solve the corresponding vector equation. We employ the method of moments, which has been widely used in electromagnetics, as a numerical tool to solve the BIEs involved. Our numerical experiments show that it can also be used robustly in elastodynamics and acoustics.     

Analysis of mutiple-shepers radiation and scattering problems by using a null-field integral equation approach

A systematic approach using the null-field integral equation in conjunction with the degenerate kernel is employed to solve the multiple radiation and scattering problems. Our approach can avoid calculating the principal values of singular and hypersingular integrals. Although we use the idea of null-field integral equation, we can locate the point on the real boundary thanks to the degenerate kernel. The proposed approach is seen as one kind of semi-analytical methods, since the error is attributed from the truncation of spherical harmonics. Finally, the numerical examples including one and two spheres are given to verify the validity of proposed approach.     

Unified boundary integral equation for the scattering of elastic and acoustic waves: solution by the method of moments

A unified boundary integral equation (BIE) is developed for the scattering of elastic and acoustic waves. Traditionally, the elastic and acoustic wave problems are solved separately with different BIEs. The elastic wave case is represented in a vector BIE with the traction and displacement vectors as unknowns whereas the acoustic wave case is governed by a scalar BIE with velocity potential or pressure as unknowns. Although these two waves can be unified in the form of a partial differential equation, the unified form in its BIE counterpart has not been reported. In this work, we derive the unified BIE for these two waves and then show that the acoustic wave case can be derived from this BIE by introducing a shielding loss for small shear modulus approximation; hence only one code needs to be maintained for both elastic and acoustic wave scattering. We also derive the asymptotic Green's tensor for zero shear modulus and solve the corresponding vector equation. We employ the method of moments, which has been widely used in electromagnetics, as a numerical tool to solve the BIEs involved. Our numerical experiments show that it can also be used robustly in elastodynamics and acoustics.     

Two-dimensional isotropic scattering in a semi-infinite cylindrical medium

Beginning with the integral equation for the source function, the solutions for the source function, flux and intensity at the boundary of a two-dimensional, isotropically scattering cylindrical medium are found. The incident radiation is collimated and normal to the surface of the medium and depends only on the radial coordinate. For a Bessel function boundary condition, separation of variables is used to reduce the source function integral equation to a one-dimensional equation. The resulting integral equation is shown to be the same as that for the two-dimensional planar case. Solutions for other boundary conditions are then shown to be superpositions of the Bessel function solution. Numerical results are presented for a Gaussian distribution of incident radiation which closely models a laser beam. These multiple scattering results are compared to the single scattering approximation. Also, the solution for a strongly anisotropic phase function which is made up of a spike in the forward direction superimposed on an otherwise isotropic phase function is expressed in terms of the isotropic results.     

Riemann?CHilbert Approach to the Elastodynamic Equation: Part I

We develop the Riemann?CHilbert (RH) approach to scattering problems in elastic media. The approach is based on the RH method introduced in the 1990s by Fokas (A unified approach to boundary value problems, CBMS-SIAM, 2008) for studying boundary problems for linear and integrable nonlinear PDEs. A suitable Lax pair formulation of the elastodynamic equation is obtained. The integral representations derived from this Lax pair are applied to Rayleigh wave propagation in an elastic half space and quarter space. The latter problem is reduced to the analysis of a certain underdetermined RH problem. We show that the problem can be re-formulated as a well-determined vector Riemann?CHilbert problem with a shift posed on a torus.     

Path of Momentum Integral in the Skorniakov-Ter-Martirosian Equation

The Skorniakov-Ter-Martirosian (STM) integral equation is widely used for the quantum three-body problems of low-energy particles (e.g., ultracold atom gases). With this equation these three-body problems can be efficiently solved in the momentum space. In this approach the boundary condition for the case that all the three particles are gathered together is described by the upper limit of the momentum integral, i.e., the momentum cutoff. On the other hand, in realistic systems, the three-body recombination (TBR) process can occur when all these three particles are close to each other. In this process two particles form a deep dimer and the other particle can gain high kinetic energy and then escape from the low-energy system. In the presence of the TBR process, the momentum-cutoff in the STM equation would include a non-zero imaginary part. As a result, the momentum integral in the STM equation should be done in the complex-momentum plane. In this case the result of the integral depends on the choice of the integral path. Obviously, only one integral path can lead to the correct result. In this paper we consider how to correctly choose the integral path for the STM equation. We take the atom-dimer scattering problem in a specific ultracold atom gas as an example, and show the results given by different integral paths. Based on the result for this case we explore the reasonable integral paths for general case.     

A new integral equation method for direct electromagnetic scattering in homogeneous media and its numerical confirmation

In this paper, we derive a new integral equation method for direct electromagnetic scattering in homogeneous media and present a numerical confirmation of the new method via a computer simulation. The new integral equation method is based on a paper written by DeSanto [1], originally for scattering from an infinite rough surface separating homogeneous dielectric half-spaces. Here, it is applied to a bounded scatterer, which can be an ohmic conductor or a dielectric, with some simplification of the continuity conditions for the fields. The new integral equation method is developed by choosing the electric field and its normal derivative as boundary unknowns, which are not the usual boundary unknowns. The new integral equation method may provide significant computational advantages over the standard Stratton-Chu method [2] because it leads to a 50% sparse, rather than 100% dense, impedance (collocation) matrix. Our theoretical development of the new integral equation method is exact.     

A fast direct solver for the integral equations of scattering theory on planar curves with corners

We describe an approach to the numerical solution of the integral equations of scattering theory on planar curves with corners. It is rather comprehensive in that it applies to a wide variety of boundary value problems; here, we treat the Neumann and Dirichlet problems as well as the boundary value problem arising from acoustic scattering at the interface of two fluids. It achieves high accuracy, is applicable to large-scale problems and, perhaps most importantly, does not require asymptotic estimates for solutions. Instead, the singularities of solutions are resolved numerically. The approach is efficient, however, only in the low- and mid-frequency regimes. Once the scatterer becomes more than several hundred wavelengths in size, the performance of the algorithm of this paper deteriorates significantly. We illustrate our method with several numerical experiments, including the solution of a Neumann problem for the Helmholtz equation given on a domain with nearly 10000 corner points.     

On some boundary value problems on a strip in the complex plane

The Sommerfeld integral inversion method for the Helmholtz equation in an angular region with different boundary values leads to boundary value problems in an infinite strip of the complex plane. We investigate a generic system for such boundary value problems and give the existence and uniqueness results with optimal growth estimates on the solution. We also give the solutions of the Dirichlet problem in a strip when the boundary functions grow exponentially.     

Simulation of electromagnetic scattering by random discrete particles using a hybrid FE-BI-CBFM technique

A hybrid finite element–boundary integral–characteristic basis function method (FE-BI-CBFM) is proposed for an efficient simulation of electromagnetic scattering by random discrete particles. Specifically, the finite element method (FEM) is used to obtain the solution of the vector wave equation inside each particle and the boundary integral equation (BIE) using Green's functions is applied on the surfaces of all the particles as a global boundary condition. The coupling system of equations is solved by employing the characteristic basis function method (CBFM) based on the use of macro-basis functions constructed according to the Foldy–Lax multiple scattering equations. Due to the flexibility of FEM, the proposed hybrid technique can easily deal with the problems of multiple scattering by randomly distributed inhomogeneous particles that are often beyond the scope of traditional numerical methods. Some numerical examples are presented to demonstrate the validity and capability of the proposed method.     

Light scattering by needle-type and disk-type particles

A powerful tool to analyze light scattering by 3D arbitrary-shaped homogeneous or inhomogeneous obstacles located in free space is based on volume integral equation. In this paper we apply a weak form of volume integral equation to simulate light scattering by needle- and disk-type particles such as straight and curved cylinders, cylindrical plate and hexagonal prism with high aspect ratio and low and high values of refractive indexes. For problems where discrete sources method could be applied, we calculated differential scattering cross-section using both methods and got excellent agreement in results.     

Numerical solution for curved crack problem in elastic half-plane using hypersingular integral equation

A hypersingular integral equation for the curved crack problems of an elastic half-plane is introduced. Formulation of the equation is based on the usage of a modified complex potential. The potential is generally expressed in the form of a Cauchy-type integral. The modified complex potential is composed of the principal part and the complementary part. The principal part of the complex potential is actually equivalent to the original complex potential for the curved crack in an infinite plate. The role of the complementary part is to eliminate the boundary traction along the boundary of the half-plane caused by the principal part. From the assumed boundary traction condition, a hypersingular integral equation is obtained for the curved crack problems of an elastic half-plane. The curve length coordinate method is used to obtain a final solution. Several numerical examples are presented that prove the efficiency of the suggested method.     

Radiation transport and internal reflection in a two region, turbid sphere

We construct an integral equation for the flux intensity in a scattering and absorbing two-region turbid spherical medium using the integro-differential form of the radiative transfer equation. The sphere is uniformly irradiated by an external source of arbitrary angular distribution and contains a distributed volume source. Anisotropic scattering is accounted for by the transport approximation. The Fresnel boundary conditions, which incorporate reflection and refraction, are used at the outer surface and at the interface between the two regions. In this respect, some new interfacial boundary conditions are introduced. For the special case of a non-scattering medium, we obtain exact solutions for specular reflection. Some numerical examples are given which show qualitative agreement with some recent work of other authors. Of particular interest are the emergent angular distribution and the albedo of the surface as a function of the refractive index and the radii of the two regions. We also draw attention to the fact that the boundary conditions at the interface differ according to the relative values of the refractive indices in the two regions. The interfacial boundary conditions for use in diffusion theory are derived and compared with those of Aronson [Boundary conditions for diffusion of light. J opt Soc Am 1995;12:2532]. In appendix B, we show how diffusion theory may be used to include scattering into the problem in a simple way.     

NUMERICAL SIMULATION OF OPEN WAVEGUIDE CONVERTERS USING FDTD METHOD

We study 3-dimensional asymmetric diffraction problems for waveguide-based electro-dynamic systems, radiating to infinite free space. For calculations we utilize the FDTD (Finite Difference Time Domain) numerical simulation method with the UPML (Unsplit Perfectly Matched Layer) absorbing boundary conditions. This paper states that the FDTD method, in spite of its relatively low calculation speed, has an approved ability of solving certain problems that cannot be solved by the other traditional numerical simulation methods (the method of integral equation, the method of scattering matrix).     

The Elliptic Sinh-Gordon Equation in the Quarter Plane

We study the elliptic sinh-Gordon equation formulated in the quarter plane by using the so-called Fokas method, which is a signi?cant extension of the inverse scattering transform for the boundary value problems. The method is based on the simultaneous spectral analysis for both parts of the Lax pair and the global algebraic relation that involves all boundary values. In this paper, we address the existence theorem for the elliptic sinh-Gordon equation posed in the quarter plane under the assumption that the boundary values satisfy the global relation. We also present the formal representation of the solution in terms of the unique solution of the matrix Riemann- Hilbert problem de?ned by the spectral functions.     

Generator coordinate amplitude for scattering

In the generator coordinate method for scattering the proper boundary condition is accomplished by requiring the GC amplitude to satisfy an integral equation of the first kind. Attempts to solve this problem are first reviewed and then an improved approximation is proposed which is applicable to a wider class of scattering problems in addition to the Coulomb scattering.A better approximation is obtained in the asymptotic region, where the generator coordinate, i.e., the distance between two shell-model wells of the fragments, is larger than the touching distance of the colliding nuclei, by deriving partial differential equations of first order for the terms of an asymptotic series in $\text{1}\text{E}$, where E is the scattering energy.Extracting the information on the GC amplitude for small values of the generator parameter from the integral equation of the first kind is an ill-posed problem. It is shown that the method of statistical regularization offers a powerful and controllable procedure to uncover the GC amplitude. The unknown GC amplitude is treated as a random function with an a priori distribution of probability which is based on the assumption that the amplitude is bounded and that the errors in the input are random with zero expectation value. A useful procedure is found for fixing parameters of the a priori distribution. The solution for small values of the GC parameter is expressed in the form of a Dini series.The method is applied to the calculation of the GC amplitude for scattering of two α-particles at 15 MeV c.m. energy. The measure of the accuracy is the difference between the input wave function of relative motion and the result of folding of the GC amplitude with the kernel of the integral equation. The prescribed accuracy is reached with this method on a much larger interval than with any previously proposed method.     

轴对称回转体瞬态温度场边界元分析

本文从三维瞬态势问题的边界积分方程及基本解出发,推导出轴对称势问题的边界积分方程及基本解,然后离散形成边界元方程,在此基础上对若干瞬态温度场数值例进行了分析计算,结果表明推导出的方法是可行的,并具有较高的精度和稳定性,能应用于工程中复杂的轴对称回转体的瞬态温度场分析。