A boundary integral method to compute Green's functions for the radiative transport equation

The Green's function for the time-independent radiative transport equation in the whole space can be computed as an expansion in plane wave solutions. Plane wave solutions are a general class of solutions for the radiative transport equation. Because plane wave solutions are not known analytically in general, we calculate them numerically using the discrete ordinate method. We use the whole space Green's function to derive boundary integral equations. Through the solution of the boundary integral equations, we compute the Green's function for bounded domains. In particular we compute the Green's function for the half space, the slab, and the two-layered half space. The boundary conditions used here are in their most general form. Hence, this theory can be applied to boundaries with any kind of reflection and transmission law.     

有界空间中的非线性声学和收敛的积累解

在微扰近似下,拉格朗日体系下的一阶、二阶波动方程解是具势运动,应用拉格朗日变动参数法来寻求积累解。在一般情况下,二阶波的波动方程在半空间会出现各式各样的积累解,它们沿着3个坐标变量的方向都有积累,在理想介质中它们不满足辐射条件。本文的结果表明,在考虑到介质的非理想性之后,也只有沿着平面边界面法线方向有积累的积累解才满足辐射条件,因而是收敛的。      

半空间电大涂敷目标散射的高频分析方法

研究了半空间内(海面,地面)电大尺寸涂敷目标散射的高频求解方法.将半空间并矢格林函数引入物理光学方法中,对半空间复杂环境影响进行考虑,并利用阻抗边界条件考虑涂敷目标表面的复杂电磁散射,推导出半空间物理光学分析方法,同时结合图形电磁学,对半空间电大涂敷目标进行消隐判断,提取像素面元法矢量和深度缓存等有效信息,快速有效地计算了半空间电大涂敷目标的雷达散射截面.数值结果证明了方法的有效性和准确性. 关键词： 半空间物理光学方法 半空间并矢格林函数 图形电磁学 雷达散射截面     

The complex equivalent source method for sound propagation over an impedance plane

The sound field caused by a monopole source above an impedance plane can be calculated by using a superposition of equivalent point sources located along a line in the mirror space below the plane. Originally, such an approach for representing the half-space Green's function was described by Sommerfeld at the beginning of the last century, in order to treat half-space problems of heat conduction. However, the representation converges only for masslike impedances and cannot be used for the more important case of reflecting planes with springlike surface impedances. The singular part of the line integral can be transformed into a Hankel function, which shows that surface waves are contained in the whole solution. Unfortunately, this representation suffers from the lack of validity at certain receiver points and from restrictions on wave number and impedance range to ensure the necessary convergence. The main idea of the present method is to use also a superposition of equivalent point sources, but to allow that these sources can be located at complex source points. The corresponding form of the half-space Green's function is suitable for both masslike and springlike surface impedances, and can be used as a cornerstone for a boundary element method.     

Comparison between radiative transfer theory and the simplified spherical harmonics approximation for a semi-infinite geometry

In this study, the third-order simplified spherical harmonics equations (SP3), an approximation of the radiative transfer equation, are solved for a semi-infinite geometry considering the exact simplified spherical harmonics boundary conditions. The obtained Green's function is compared to radiative transfer calculations and the diffusion theory. In general, it is shown that the SP3 equations provide better results than the diffusion approximation in media with high absorption coefficient values but no improvement is found for small distances to the source.     

On the three-body Coulomb scattering problem

We describe the smoothness properties and the asymptotic form of the Green's function (in configuration space) for three charged particles. We also discuss the integral equations and the boundary value problems for the Coulomb wavefunctions and we show that they form a complete set. Finally, we study the singularities of the Coulomb scattering operator, and we investigate the connection between the Dollard wave operators and the Coulomb wavefunctions.     

三维润滑问题的数值分析

本文求解了轴承的三维润滑问题。使用格林函数方法,将基本方程化为求解区边界上的积分方程。采用空间曲边三角形单元离散边界,并用六结点等参数无数值求解积分方程,确定了速度场,压力场和边界上的应力。算例与解析解比较表明有较高的精度。其它一些典型问题的计算也给出满意的结果,表明本文方法对数值分析这类问题的有效性。     

TOMOGRAPHY FORMULA FOR BIOCHEMICAL IMAGING OF THIN TISSUE WITH DIFFUSE-PHOTON DENSITY WAVES

Using the transport theory to describe the near infrared light propagating in tissue with finite parallel-plane geometry, and taking the zero-boundary condition, we obtain the analytical expression of average photon density and Green's function incorporating the boundary effects in the homogeneous tissure. Making use of perturbation theory we also obtain the analytical expression of scattered wave induced by the heterogeneity, and present the 2-dimensional spatial transform of scattered wave with respect to transverse coordinate. If the information of heterogeneity on depth and thickness is available, diffraction tomography formula is presented to save the time of image reconstruction; if the information is unknown, we suggest to obtain the inhomogeneous function from the one-dimensional integral equation of 2-dimensional spatial transform of scattered wave applying the direct matrix method or iterative method for image reconstruction. This approach avoids directly solving three-dimensional integral equation of scattered wave. In our proposed approach the strong points of the direct matrix method, iterative method, and diffraction tomography are fully combined.     

雪层覆盖地面电磁散射的矩量法

为实现分层介质粗糙面电磁散射的矩量法研究,给出一种分层介质粗糙面电磁积分方程的区域分解方法.将格林定理应用于粗糙面所分的各子空间,结合波动方程和格林函数推导分层粗糙面的电磁积分方程,利用矩量法对其进行离散,数值计算得到雪层覆盖地面散射系数的角分布曲线,其中,粗糙表面由一维带限Weierstrass分形谱和Monte Carlo方法模拟.通过与时域有限差分法数值结果的比对,验证该方法的准确性,并分析散射系数随雪和地面参数、介质参数以及入射波参数的变化,获得了较完整的散射特征.     

Diffraction of sound by an elastic cylinder near the surface of an elastic halfspace

Diffraction of an acoustic wave by an elastic cylinder near the surface of an elastic halfspace is considered. The solution relies on a Helmholtz-type integral equation and uses the Green function of an elastic halfspace. The latter function is represented in the form of an integral over the Sommerfeld contour on the plane of a complex variable that has the meaning of the angle of the wave incidence on the halfspace boundary. An integral equation for the sound pressure distribution over the cylinder surface is derived. This equation is reduced to an infinite system of equations for the Fourier-series expansion coefficients of this distribution. The results obtained are valid for the diffraction of a cylindrical wave and a plane wave. They also describe the diffraction of a spherical wave when the transmitter and receiver are far from the cylinder and lie in one plane that is orthogonal to the cylinder axis.     

The Green's matrix and the boundary integral equations for analysis of time-harmonic dynamics of elastic helical springs

Helical springs serve as vibration isolators in virtually any suspension system. Various exact and approximate methods may be employed to determine the eigenfrequencies of vibrations of these structural elements and their dynamic transfer functions. The method of boundary integral equations is a meaningful alternative to obtain exact solutions of problems of the time-harmonic dynamics of elastic springs in the framework of Bernoulli-Euler beam theory. In this paper, the derivations of the Green's matrix, of the Somigliana's identities, and of the boundary integral equations are presented. The vibrational power transmission in an infinitely long spring is analyzed by means of the Green's matrix. The eigenfrequencies and the dynamic transfer functions are found by solving the boundary integral equations. In the course of analysis, the essential features and advantages of the method of boundary integral equations are highlighted. The reported analytical results may be used to study the time-harmonic motion in any wave guide governed by a system of linear differential equations in a single spatial coordinate along its axis.     

The use of the virtual source technique in computing scattering from periodic ocean surfaces

In this paper the virtual source technique is used to compute scattering of a plane wave from a periodic ocean surface. The virtual source technique is a method of imposing boundary conditions using virtual sources, with initially unknown complex amplitudes. These amplitudes are then determined by applying the boundary conditions. The fields due to these virtual sources are given by the environment Green's function. In principle, satisfying boundary conditions on an infinite surface requires an infinite number of sources. In this paper, the periodic nature of the surface is employed to populate a single period of the surface with virtual sources and m surface periods are added to obtain scattering from the entire surface. The use of an accelerated sum formula makes it possible to obtain a convergent sum with relatively small number of terms (～40). The accuracy of the technique is verified by comparing its results with those obtained using the integral equation technique.     

A BEM approach to validate a model for predicting sound propagation over non-flat terrain

A two-dimensional boundary element model for sound propagation in a homogeneous atmosphere above non-flat terrain has been constructed. An infinite impedance plane is taken into account in the Green's function in the underlying integral equation, so that only the non-flat parts of the terrain need to be discretised in the boundary element model. This Green's function is undefined for points below the impedance plane, and therefore valleys and hollows are taken into account by coupling the exterior domain above the ground with one or several interior domains below the ground, as suggested in a recent paper [J. Sound Vibrat. 223 (1999) 355]. The resulting BEM model, which can handle arbitrary combinations of barriers and hollows, has been used for validating a ray model for various difficult configurations, including combinations of valleys and barriers.     

Three-dimensional radiative transfer in an anisotropically scattering, plane-parallel medium: generalized reflection and transmission functions

The problem of spatially varying, collimated radiation incident on an anisotropically scattering, plane-parallel medium is considered. A very general phase function is allowed. An integral transform is used to reduce the three-dimensional radiative transport equation to a one-dimensional form, and a modified Ambarzumian's method is applied to derive nonlinear integral and integro-differential equations for the generalized reflection and transmission functions. The integration is over the polar and azimuthal angles—this formulation is referred to as the double-integral formulation. The integral equations are used to illustrate symmetry relationships and to obtain single- and double-scattering approximations. The generalized reflection and transmission functions are important in the construction of the solutions to many multidimensional problems. Coupled integral equations for the interior and emergent intensities are developed and, for the case of two identical homogeneous layers, used to formulate a doubling procedure. Results for an isotropic and Rayleigh scattering medium are presented to illustrate the computational characteristics of the formulation.     

Radiative Transport in a Periodic Structure

We derive radiative transport equations for solutions of a Schrödinger equation in a periodic structure with small random inhomogeneities. We use systematically the Wigner transform and the Bloch wave expansion. The streaming part of the radiative transport equations is determined entirely by the Bloch spectrum, and the scattering part by the random fluctuations.     

On the sound field of a shallow spherical shell in an infinite baffle

A method is presented for calculating the far field sound radiation from a shallow spherical shell in an acoustic medium. The shell has a concentrated ring mass boundary condition at its perimeter representing a loudspeaker voice coil and is excited by a concentrated ring force exerted by the end of the voice coil. A Green's function is developed for a shallow spherical shell, which is based upon Reissner's solution to the shell wave equation [Q. Appl. Math. 13, 279-290 (1955)]. The shell is then coupled to the surrounding acoustic medium using an eigenfunction expansion, with unknown coefficients, for its deflection. The resulting surface pressure distribution is solved using the King integral together with the free space Green's function in cylindrical coordinates. In order to eliminate the need for numerical integration, the radiation (coupling) integrals are solved analytically to yield fast converging expansions. Hence, a set of simultaneous equations is obtained which is solved for the coefficients of the eigenfunction expansion. These coefficients are finally used in formulas for the far field sound radiation.     

Boundary conditions in the radiative relaxation problem

The problem of the dissipation of temperature perturbations in a finite homogeneous atmosphere is solved for the situation in which the temperature at one boundary is maintained constant (that is, the temperature perturbation is zero for all times) while energy can be freely radiated to space through the other boundary. Exact solutions are shown for the exponential-sum fit to the kernel of the basic integral equation. These solutions constitute the set of radiative eigenfunctions. Also, approximate solutions in terms of the radiative eigenfunctions in the diffusion approximation (one exponential term in the expansion of the kernel) are obtained. These, in turn, are used in the solution of an initial value problem. The constant temperature boundary condition simulates the interface between two regions in one of which the relaxation processes are much more rapid than the purely radiative relaxation of the other.     

Focusing of electromagnetic waves into a uniaxial crystal

We derive integral representations suitable for studying the focusing of electromagnetic waves through a plane interface into a uniaxial crystal. To that end we start from existing exact solutions for the transmitted fields due to an arbitrary three-dimensional (3D) wave that is incident upon a plane interface separating two uniaxial crystals with arbitrary orientation of the optical axis in each medium. Then we specialize to the case in which the medium of the incident wave is isotropic and derive explicit expressions for the dyadic Green's functions associated with the transmitted fields as well as integral representations suitable for asymptotic analysis and efficient numerical evaluation. Relevant integral representations for focused 3D electromagnetic waves are also given. Next we consider the special case in which (i) the incident field is a two-dimensional (2D) TM wave and (ii) the optical axis in the crystal lies in the plane of incidence, implying that we have a 2D vectorial problem, and derive dyadic Green's functions, integral representations suitable for asymptotic and numerical treatment, and integral representations for focused TM fields. Numerical results for focused 2D TM fields based on these integral representations as well as corresponding experimental results will be presented in forthcoming papers.     

Traveling wave solutions for two nonlinear evolution equations with nonlinear terms of any order

In this paper, based on the known first integral method and the Riccati sub-ordinary differential equation (ODE) method, we try to seek the exact solutions of the general Gardner equation and the general Benjamin-Bona-Mahoney equation. As a result, some traveling wave solutions for the two nonlinear equations are established successfully. Also we make a comparison between the two methods. It turns out that the Riccati sub-ODE method is more effective than the first integral method in handling the proposed problems, and more general solutions are constructed by the Riccati sub-ODE method.     

Traveling wave solutions for two nonlinear evolution equations with nonlinear terms of any order

In this paper, based on the known first integral method and the Riccati sub-ordinary differential equation (ODE) method, we try to seek the exact solutions of the general Gardner equation and the general Benjamin-Bona-Mahoney equation. As a result, some traveling wave solutions for the two nonlinear equations are established successfully. Also we make a comparison between the two methods. It turns out that the Riccati sub-ODE method is more effective than the first integral method in handling the proposed problems, and more general solutions are constructed by the Riccati sub-ODE method.