Integration of Chandrasekhar's integral equation

We solve Chandrasekhar's integration equation for radiative transfer in the plane-parallel atmosphere by iterative integration. The primary thrust in radiative transfer has been to solve the forward problem, i.e., to evaluate the radiance, given the optical thickness and the scattering phase function. In the area of satellite remote sensing, our problem is the inverse problem: to retrieve the surface reflectance and the optical thickness of the atmosphere from the radiance measured by satellites. In order to retrieve the optical thickness and the surface reflectance from the radiance at the top-of-the atmosphere (TOA), we should express the radiance at TOA “explicitly” in the optical thickness and the surface reflectance. Chandrasekhar formalized radiative transfer in the plane-parallel atmosphere in a simultaneous integral equation, and he obtained the second approximation. Since then no higher approximation has been reported. In this paper, we obtain the third approximation of the scattering function. We integrate functions derived from the second approximation in the integral interval from 1 to ∞ of the inverse of the cos of zenith angles. We can obtain the indefinite integral rather easily in the form of a series expansion. However, the integrals at the upper limit, ∞, are not yet known to us. We can assess the converged values of those series expansions at ∞ through calculus. For integration, we choose coupling pairs to avoid unnecessary terms in the outcome of integral and discover that the simultaneous integral equation can be deduced to the mere integral equation. Through algebraic calculation, we obtain the third approximation as a polynomial of the third degree in the atmospheric optical thickness.     

Knudsen Layer for Gas Mixtures

The Knudsen layer in rarefied gas dynamics is essentially described by a half-space boundary-value problem of the linearized Boltzmann equation, in which the incoming data are specified on the boundary and the solution is assumed to be bounded at infinity (Milne problem). This problem is considered for a binary mixture of hard-sphere gases, and the existence and uniqueness of the solution, as well as some asymptotic properties, are proved. The proof is an extension of that of the corresponding theorem for a single-component gas given by Bardos, Caflisch, and Nicolaenko [Comm. Pure Appl. Math. 39:323 (1986)]. Some estimates on the convergence of the solution in a finite slab to the solution of the Milne problem are also obtained.     

水下光场的迭代求解

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

Stochastic modeling of coherent phenomena in strongly inhomogeneous media

A procedure of numerical simulation for coherent phenomena in multiply scattering media is developed on the basis of the juxtaposition of a Monte Carlo stochastic method with an iterative approach to the solution of the Bethe-Salpeter equation. The time correlation function and the interference component of coherent backscattering are calculated for scalar and electromagnetic fields. The results of simulation are in good agreement with experimental results, as well as with theoretical results obtained by generalizing the Milne solution.     

Anisotropic surface relief diffraction gratings under arbitrary plane wave incidence

A modal method is used for the analysis under oblique incidence of a diffraction grating made of anisotropic material. The problem is studied viewing the structure as the cascade of junctions between periodic arrays of anisotropic slab waveguides with the same period and different heights. This diffraction problem is formulated in terms of an integral equation that enforces the continuity of the transverse magnetic field at the junction. The unknown is the transverse electric field at the junction. It is possible to use also another formulation, where the role of the two fields is exchanged. The kernels of these equations are the relevant Green's functions, which are expressed in terms of eigenfunction expansions. The determination of the modes of the various regions composed of arrays of anisotropic dielectric slabs has been carried out by the method of spectral elements, whereby the field components are represented in a polynomial basis and the original differential eigenvalue problem is converted into an algebraic one. The integral equation is solved numerically by the method of moments and each junction is characterized by its generalized scattering matrix (GSM). Finally, the diffraction efficiencies of the grating are obtained by combining the various GSM's.     

An approximate analytical solution of the radiative transfer equation in hydrometeors accounting for scattering effects

An iterative method for the radiative transfer equation solution is suggested for the scattering hydrometeors. The method allows to find the layer scattering indicatrix. The method is shown to be more general as compared with the four-flux theory or the perturbation method and differs strongly from these ones because even in the case of the first iteration it gives the results which are very close to ones obtained by numerical methods for all rain scattering angles.     

Mass operator for wave scattering from a slightly random surface

This paper deals with the mass operator representing multiple-scattering effects in the theory of wave scattering from a slightly random surface. By means of the stochastic-functional approach, a recurrence equation for the mass operator is obtained in the form of an iterative integral. However, its solution oscillates in a non-physical manner against the number of iterations. Next, the recurrence equation may be regarded as a nonlinear integral equation, when the number of iterations goes to infinity. An analytical solution of the nonlinear integral equation is presented for a special case in which the roughness spectrum is the Dirac delta function. Then, the nonlinear integral equation is solved numerically for the Gaussian roughness spectrum by iteration, starting from such an analytical solution. It is shown that only a few iterations are required to obtain the mass operator, even when the correlation distance is small. Effects of the mass operators on the coherent reflection coefficient and the incoherent scattering cross section are calculated and shown in figures.     

The Milne problem in a continuous stochastic planar medium with Gaussian statistics

The Milne problem of radiative transfer in a planar medium, with isotropic scattering is considered. The medium is assumed to be continuous stochastic medium, with fluctuations described as Gaussian field. Pomraning-Eddington method is used to obtain an explicit form for the radiation energy density in the deterministic case. It depends on two random variables, namely the optical space variable and the optical thickness of the medium. The Gaussian joint probability density function of these two random variables is defined and used to find the ensemble-averaged energy density and the linear extrapolation distance. It is shown that the statistical nature of the medium leads to two quite different solutions of the Milne problem. Numerical results are implemented for the sake of clarification.     

Random integral equation formulation of a generalized Langevin equation

In this paper, the generalized Langevin equation introduced by Kubo and Mori is formulated as a random integral equation. We consider (1) the existence and uniqueness of the solution, (2) moments of the solution process, (3) a comparison theorem for solution processes, and (4) the Cauchy polygonal approximation to the solution.     

Whittaker方程的场方法

用场方法来求解Whittaker方程.将一个场变量取作为其余场变量和时间的函数并对这个函数建立基本偏微方程.如能求得它的完全积分，那么Whittaker方程的解可由解代数方程来得到. 关键词： 场变量 基本偏微分方程 场方法 积分     

Dark solitons of the Biswas–Milovic equation by the first integral method

This paper studies the Biswas–Milovic equation that is a generalized version of the familiar nonlinear Schrodinger's equation describing the propagation of solitons through optical fibers for trans-continental and trans-oceanic distances with Kerr law nonlinearity by the aid of the first integral method. The dark 1-soliton solution is retrieved by the aid of this method and a couple of other singular periodic solutions are also obtained.     

Large energy behavior of the velocity distribution for the hard-sphere gas

The initial-value problem for the Boltzmann-Lorentz equation for hard spheres at zero temperature is shown to be ill defined, the general solution depending on an arbitrary function. The uniqueness of the solution can be obtained by imposing the conservation of the number of particles (Carleman's type of condition does not suffice). The linearized Boltzmann equation for hard spheres is then analyzed, as it occurs in Enskog's method for calculating transport coefficients. It is demonstrated that in the case of viscosity and diffusion it is necessary to add supplementary conditions to obtain the uniqueness of the solution. The nonuniform character of Enskog's expansion and violation of positivity in the large velocity region are exhibited.     

Simplification of an analytical solution of the fourth-moment equation

The equation for the fourth moment of a wave propagating in a multiply scattering random medium has been solved by various methods. When the analytical solutions are compared with numerical solutions of the equation it is found that the fundamental solution together with a first-order correction term agree very closely with the numerical results over a wide range of distances and scattering strengths. Unfortunately, the correction term involves multiple integrals and so is difficult to evaluate. This paper shows how some of these integrations can be carried out and the results combined in such a way that an analytical form similar to the fundamental solution is obtained involving only a single integral. This simplified combined solution also agrees very closely with the numerical results.     

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.     

各向异性介质三维电磁响应模拟的Ho-GEBA算法

本文基于积分方程法研究并建立了一种模拟横向同性介质中任意各向异性异常 体三维电磁响应的高阶广义扩展Born近似(Ho-GEBA)算法. 首先利用逐次迭代技术给出积分方程的广义级数展开解, 为保证其收敛性, 引入一种各向异性条件下满足压缩映射的迭代算子. 然后利用异常体区域分解技术, 并结合扩展Born近似原理, 得到各向异性介质三维电磁响应的Ho-GEBA解. 为提高效率, 计算过程中采用并矢Green函数的解析表达式. 最后通过数值计算实例对比验证了本文算法的有效性. 关键词： 高阶广义扩展Born近似 积分方程 电磁模拟 解析Green函数     

Coherent effects of multiple scattering for scalar and electromagnetic fields: Monte-Carlo simulation and Milne-like solutions

Based on an expansion of the Bethe-Salpeter equation in scattering orders a semi-analytic approach for simulation of coherent phenomena of multiple scattering in random media has been developed. We found that for scalar field the manifestation of these phenomena, observed as temporal field correlation function and coherent backscattering, are universal and well agreed with the results predicted by diffusion approximation. For the electromagnetic field the temporal correlation function and coherent backscattering are noticeably differ from those found for the scalar field, depending strongly on the scattering anisotropy. The obtained numerical results, for the first time to our knowledge, are compared directly with the known generalizations of the Milne solution.     

Coherent effects in multiple scattering of linearly polarized light

Comparing the stochastic Monte Carlo technique with the iteration procedure for solving the Bethe-Salpeter equation in the framework of numerical simulation, the time correlation function and the interference component of the coherent backscattering of a linearly polarized light wave in a multiply scattering medium are calculated. The results of the simulation agree well with theoretical results obtained by generalizing the Milne solution, as well as with experimental data.     

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.     

Inverse scattering transform for the time dependent Schrödinger equation with applications to the KPI equation

For the direct-inverse scattering transform of the time dependent Schrödinger equation, rigorous results are obtained based on an opertor-triangular-factorization approach. By viewing the equation as a first order operator equation, similar results as for the first ordern x n matrix system are obtained. The nonlocal Riemann-Hilbert problem for inverse scattering is shown to have solution.     

势问题的径向基函数——局部边界积分方程方法

将基于径向基函数构造的具有插值特性的近似函数和局部边界积分方程方法相结合，建立了求解势问题的径向基函数——局部边界积分方程方法，推导了相应离散方程.与其他边界积分方程的无网格方法相比，本文方法具有数值实现过程简单、计算量小、精度高的优点，并可直接施加边界条件.最后通过算例说明了该方法的有效性. 关键词： 径向基函数 无网格方法 局部边界积分方程 势问题