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.     

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.     

Variational principles for scattering in the generator coordinate method

The generator coordinate method (GCM) wave function is used as a trial function in a Kohn type variational principle for scattering phase shifts. It is shown that a GCM trial function is a solution of the variational equations if the Hill-Wheeler integral equation is satisfied subject to an appropriate boundary condition. A new method for introducing the scattering boundary condition is presented. There is a uniqueness theorem for the phase shift.     

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.     

Remarks on boundary conditions for scalar scattering

The question is discussed whether potential scattering problems can be treated as boundary value problems associated with differential equations, as is sometimes suggested in the literature. We show that, except in some very special cases, this is not possible. The values of the wave function and its normal derivative on the boundary of a finite-range potential cannot be prescribed arbitrarily but are implicit in the integral equation of potential scattering. We derive two coupled singular integral equations for the boundary values for the case when the scattering potential is homogeneous.     

水下光场的迭代求解

水体对光线的散射是水下图像质量劣化的重要因素,为了定量分析在特定光源照射下水体散射的影响,建立了光线水下传输的散射模型,以此为基础推导出求解水下光场分布的Fredholm积分方程.在水中光线能量随距离的增大呈指数规律衰减,基于此,在水体体散射函数为常数的情况下,给出了有边界条件时该积分方程的数值迭代求解方法,得到高精度...     

Boundary integral equation method for analysis of light scattering by 2D nanoparticles

The scattering of a plane electromagnetic wave by 2D nanoparticles is simulated using the boundary integral equation method. The accuracy and smoothness of the solution are improved by applying Hermitean interpolants for discretization. The accuracy of the numerical approach used in this paper is comparable to the Mie classical solution for scattering by a circular cylinder. It is shown that laser radiation incident on silver nanoparticles generates plasma resonances, giving rise to peak values of the extinction cross section and considerably enhancing the local electromagnetic field on the surface.     

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.     

A fast directional algorithm for high-frequency electromagnetic scattering

This paper is concerned with the fast solution of high-frequency electromagnetic scattering problems using the boundary integral formulation. We extend the O(N log N) directional multilevel algorithm previously proposed for the acoustic scattering case to the vector electromagnetic case. We also detail how to incorporate the curl operator of the magnetic field integral equation into the algorithm. When combined with a standard iterative method, this results in an almost linear complexity solver for the combined field integral equations. In addition, the butterfly algorithm is utilized to compute the far field pattern and radar cross section with O(N log N) complexity.     

An iterative method to solve acoustic scattering problems using a boundary integral equation

In this work, a simple iterative method to solve the acoustic scattering/radiation problems using the boundary integral equation (BIE) formulation is presented. The operator equation obtained in the BIE formulation is converted into a matrix equation using the well-known method of moments solution procedure. The present method requires much fewer mathematical operations per iteration when compared to other available iterative methods. Further, the present iterative method can easily handle multiple incident fields, a highly desirable feature not available in any other iterative method, much the same way as direct solution techniques. Several numerical examples are presented to illustrate the efficiency and accuracy of the method.     

Efficient wave scattering analysis for damaged cylindrical waveguides

This paper addresses scattering effects at arbitrarily shaped defects in waveguides with finite cross-sections. The main subject is to predict scattered wavefields induced by an incident wavefield in order to localize and characterize damages. The numerical approach involves the solution of a boundary value problem in combination with a decomposition method of scattered wavefields. In order to investigate mode conversion phenomena at defects, a boundary element model of the damaged waveguide section is built up. Implementation of the Boundary Element Method based on the elastodynamic boundary integral equation and the Waveguide Finite Element Method allows for a numerically efficient calculation of scattering coefficients. For defect characterization and mode sensitivity analysis, various types of surface opening defects as well as shear fractures are considered. Numerical results are presented for cylindrical waveguides and are verified experimentally.     

A solution of radiative transfer in isotropic plane-parallel atmosphere by integrating Milne equation

In this paper we solve the inversion problem of the radiative transfer process in the isotropic plane-parallel atmosphere by iterative integrations of the Milne integral equation. As a result, we obtain the scattering function in the form of a cubic polynomial in optical thickness. The author has already solved the same problem by iterative integrations of Chandrasekhar's integral equation. In the Milne integral equation, both the cosines of the viewing angles and the optical thickness are integral variables, while in Chandrasekhar's integral equation the cosines of the viewing angles are variables but the optical thickness is not. We derive several series of exponential-like functions as intermediate derivations. Their convergences are evaluated by the author's previous work in the solution of Chandrasekhar's integral equation. The truncated scattering function up to the third order in optical thickness thus obtained is identical to that obtained from Chandrasekhar's integral equation, though their apparent forms are different. Chandrasekhar pointed out that the solution of Chandrasekhar's integral equation does not have a uniqueness of solution. The Milne equation, in contrast, has been proven to have a unique solution. We discuss the uniqueness of the solution by these two methods.     

Numerical analysis of the one-dimensional heat equation subject to a boundary integral specification

In this research a numerical technique is developed for the one-dimensional heat equation that combines classical and integral boundary conditions. New matrix formulation techniques with arbitrary polynomial bases are proposed for the numerical/analytical solution of this kind of partial differential equation. We give a simple and efficient algorithm based on an iterative process for numerical solution of the method.     

Demonstration of Inverse Scattering Transform in G-L-M Formalism for NLS Equation in Normal Dispersion with Non-vanishing Boundary

Using the integral representation of the Jost solution, we deduce some conditions as the kernel function N(x, y,t) if the Jost solution satisfies the two Lax equations. Then we verify the multi-soliton solution of NLS equation with non-vanishing boundary conditions if we prove that these conditions can be demonstrated by the GLM equation, which determines the kernel function N(x,y, t) in according to the inverse scattering method.     

A renormalized equation for the three-body system with short-range interactions

We study the three-body system with short-range interactions characterized by an unnaturally large two-body scattering length. We show that the off-shell scattering amplitude is cutoff independent up to power corrections. This allows us to derive an exact renormalization group equation for the three-body force. We also obtain a renormalized equation for the off-shell scattering amplitude. This equation is invariant under discrete scale transformations. The periodicity of the spectrum of bound states originally observed by Efimov is a consequence of this symmetry. The functional dependence of the three-body scattering length on the two-body scattering length can be obtained analytically using the asymptotic solution to the integral equation. An analogous formula for the three-body recombination coefficient is also obtained.     

Thermal analysis for transient radiative cooling of a conducting semitransparent layer of ceramic in high-temperature applications

A method is developed for obtaining transient temperature distribution in a cooling semitransparent layer of ceramic. The layer is emitting, absorbing, isotropically scattering and heat conducting with a refractive index ranging from 1 to 2. The solution involves solving simultaneously the energy equation and the integral equation for the radiative flux gradient. The energy equation is solved using an implicit finite volume scheme and the integral equation of radiative heat transfer is solved using the singularity technique and Gaussian integration. The effects of scattering are investigated. It is shown that scattering has a significant effect on the transient temperature distribution and the transient mean temperature of the layer.     

An application of the K-integrals for solving the radiative transfer equation with Fresnel boundary conditions

We introduce the reader to an approximate method of solving the transport equation which was developed in the context of neutron thermalisation by Kladnik and Kuscer in 1962 [Kladnik R, Kuscer I. Velocity dependent Milne's problem. Nucl Sci Eng 1962;13:149]. Essentially the method is based upon two special weighted integrals of the one-dimensional transport equation which are valid regardless of the boundary conditions, and any solution must satisfy these integral relationships which are called the K-integrals. To obtain an approximate solution to the transport equation we turn the argument around and insist that any approximate solution must also satisfy the K-integrals. These integrals are particularly useful when the problem under consideration cannot be solved easily by analytic methods. It also has the marked advantage of being applicable to problems where there is energy exchange in a collision and anisotropy of scattering. To establish the feasibility of the method we obtain a number of approximate solutions using the K-integral method for problems to which we have exact analytical solutions. This enables us to validate the method. It is then applied to a new problem that has not yet been solved; namely the calculation of the discontinuity in the scalar intensity at the boundary between two optically dissimilar materials.     

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.     

Non-LTE transfer with complete redistribution. Scaling laws for a slab

Scaling laws for conservative scattering in a finite slab are extracted from an asymptotic analysis of the integral equation for the source function. The solution is separated into an interior and a surface boundary layer part. The matching between the two parts provides a scaling law for the surface value of the source function. When expressed in terms of a mean number of scatterings, this scaling law is generalizable to non-conservative scattering and the empirical formula of Jones and Skumanich is recovered.