Light transfer problem in turbid media with surface reflectivity

Light transfer problems in turbid media with surface reflectivity satisfying Fresnel's law are formulated. The intensity of light is considered as a sum of collimated and diffuse radiances. The problem with a collimated source and boundary conditions with surface reflectivity is solved in terms of the corresponding source-free problem with simple boundary conditions. In order to solve the source-free problem, two-flux models in the differential and integral forms are obtained. For the differential form, weight functions are introduced in order to force the boundary conditions to be met, while in the integral form the boundary conditions are embedded. The integral form has the advantage over the differential form that it permits one to extend the validity of the model and to consider inhomogeneous media. Numerical results are given for albedos and partial fluxes depending upon the refractive index of the medium. The calculations are performed in both semi-infinite and finite media and are compared with the published calculations.     

Pomraning–Eddington approximation for radiative transfer in a homogeneous solid cylinder

Abstract

The radiative transfer in a solid cylinder containing a homogeneous turbid medium with anisotropic scattering is considered. The medium has a diffuse and specular reflecting boundary illuminated by an external incidence and contains an internal energy source. This general problem can be solved in terms of the solution of the corresponding source-free problem with a specular reflecting boundary and isotropic external incidence. The Pomraning–Eddington approximation is used to solve the source-free problem. Three different weight functions are used to verify the boundary condition to find the constants of the solution. The partial flux, the irradiance and the net flux at the boundary for the general problem are calculated.     

Stochastic radiative transfer in a finite plane-parallel medium with general boundary conditions

The time-independent radiative transfer problem in a scattering and absorbing planar random medium with general boundary conditions and internal energy source is considered. The medium is assumed to consist of two randomly mixed immiscible fluids, with the mixing statistics described as a two-state homogeneous Markov process. The problem is solved in terms of the solution of the corresponding free-source problem with simple boundary conditions which is solved using Pomraning-Eddington approximation in the deterministic case. A formalism, developed to treat radiative transfer in statistical mixtures, is used to obtain the ensemble-averaged solution. The average partial heat fluxes are calculated in terms of the albedoes of the source-free problem. Results are obtained for isotropic and anisotropic scattering for specular and diffused reflecting boundaries.     

基于P3近似的空间分辨漫反射研究

光源附近组织的空间分辨漫反射是近年来生物医学光子学领域的一个研究热点,其目的是发展一种能够测定活体生物组织光学参量的新技术。漫射近似理论研究光源附近组织的空间分辨漫反射具有很大局限性。P3近似理论考虑了相函数的三阶矩,能较准确地描述光源附近组织的光辐射分布。研究了基于P3近似的空间分辨漫反射,从输运理论的PN方程组出发,导出了P3近似方程组和P3近似的格林函数解;阐述了漫射近似与P1近似的关系,给出了外推边界条件下,准直光束近似后的P3近似漫反射率的完整表示,讨论了相函数二阶参量对P3近似漫反射的影响,并与漫射近似和蒙特卡罗模拟结果进行了比较,指出了P3近似的应用范围。     

Two-dimensional radiative transfer in a cylindrical geometry with anisotropic scattering

Exact integral equations are derived describing the source function and radiative flux in a two-dimensional, radially infinite cylindrical medium which scatters anisotropically. The problem is two-dimensional and cylindrical because of axisymmetric loading. Radially varying collimated radiation is incident normal to the upper surface while the lower boundary has no radiation incident upon it. The scattering phase function is represented by a spike in the forward direction plus a series of Legendre polynomials. The two-dimensional integral equations are reduced to a one-dimensional form by separating variables for the case when the radial variation of the incident radiation is a Bessel function. The one-dimensional form consists of a system of linear, singular Fredholm integral equations of second kind. Other more complex boundary conditions are shown to be solvable by a superposition of this basic Bessel function case. Diffusely incident radiation is also considered.     

Light scattering by multilayer ellipsoids in the Rayleigh approximation

A problem of light scattering by multilayer confocal ellipsoids is solved in the Rayleigh approximation. The electric field of a light wave is assumed constant and a set of Laplace equations with the corresponding boundary conditions is considered. The final expression for the polarizability of a particle is represented in the matrix form (2×2 matrices) in terms of parameters of a nucleus and subsequent layers. Numerical calculations of the scattering and absorption efficiencies of small multilayer spheres obtained using the exact (the generalization of the Mie theory) and approximate solutions well agree with each other.     

Relation between multidimensional radiative transfer in cylindrical and rectangular coordinates with anisotropic scattering

An exact integral equation is derived for the source function in a three-dimensional rectangular medium which scatters anisotropically. The upper boundary of the finite medium is exposed to collimated radiation, while the lower boundary has no radiation incident on it. The problem is multidimensional because the incident radiation varies spatially. The scattering phase function is represented by a series of Legendre polynomials. A double Fourier transform reduces the problem to a one-dimensional integral equation for the source function. The transformed equation is compared with the integral equation for a two-dimensional cylindrical medium which scatters anisotropically and is exposed to Bessel-varying collimated radiation. A simple relation is found between the two source functions which will greatly reduce the number of computations required for the three-dimensional case. The relation also illustrates the wide utility of the generalized one-dimensional source function. Simplification of the two-dimensional rectangular case to the generalized source function is also presented. The results are extended to problems with a strong anisotropic phase function which has a diffraction spike in the forward direction.     

Two-dimensional isotropic scattering in a semi-infinite medium

Exact results are presented for the source function, radiative flux, and intensity at the boundary of a two-dimensional, isotropically scattering, semi-infinite medium subjected to collimated or diffuse radiation. The spatial distributions of incident radiation considered are (1) cosine-varying, (2) semi-infinite step, (3) step at the origin and (4) finite strip. Two-dimensional effects are most pronounced at large albedos.     

Two-dimensional linearly anisotropic scattering in a semi-infinite cylindrical medium exposed to a laser beam

A modification of Ambarzumian's method is used to develop the integro-differential equations for the source function, flux, and intensity at the boundary of a two-dimensional, semi-infinite cylindrical medium which scatters linearly. The incident radiation is collimated, normal to the top surface of the medium, and is dependent only on the radial coordinate. The radial variation is assumed to be a Bessel function or a Gaussian distribution. The Gaussian boundary condition is used to simulate a laser beam. Numerical results are presented in graphical and tabular forms for both boundary conditions. Results for forward and backward scattering phase functions are compared with those for isotropic scattering. A method is presented for extending these results to the problem of a strongly anisotropic phase function which is made up of a spike in the forward direction superimposed on a linear phase function.     

Vavilov-Cherenkov radiation from a model cascade developing near the lunar surface

The model problem of Vavilov-Cherenkov radiation emergence from lunar regolith into vacuum under LORD experimental conditions is considered. The boundary problem on radiation emergence into vacuum is solved numerically in the given field approximation (Kirchhoff approximation) from the cascade near the boundary (near-field region of the radiation source). The results are of great importance to interpret future experimental data.     

Two-dimensional radiative transfer due to curved Dirac delta line sources

In this article, the two-dimensional radiative transport equation is considered for the curved Dirac delta line source. In order to account this source type, Green’s function of the radiative transport equation for the half-plane is derived in parts of the ballistic and diffuse contribution as well as under consideration of the Fresnel reflection at the boundary. The final results are verified with the Monte Carlo method for different collimated beams and several curved sources such as the elliptic and logarithmic spiral line source.     

Three-dimensional radiative transfer in an isotropically scattering, plane-parallel medium: generalized X- and Y-functions

The topic of this work is the generalized X- and Y-functions of multidimensional radiative transfer. The physical problem considered is spatially varying, collimated radiation incident on the upper boundary of an isotropically scattering, plane-parallel medium. An integral transform is used to reduce the three-dimensional transport equation to a one-dimensional form, and a modified Ambarzumian's method is used to derive coupled, integro-differential equations for the source functions at the boundaries of the medium. The resulting equations are said to be in double-integral form because the integration is over both angular variables. Numerical results are presented to illustrate the computational characteristics of the formulation.     

放射性同位素伽玛源准直照射辐射场模拟研究

针对各向同性伽玛源参考辐射场尺寸关键技术问题,GB/T 12162系列标准虽然进行了相关规定,但是该规定并未对准直照射状态下照射室尺寸提出具体要求。为减小用于辐射检测或监测类仪器仪表检定与量值校准时伽玛参考辐射场内散射影响,本文采用蒙特卡罗方法,研究了同位素放射源准直照射时,照射室尺寸变化对检验点处的剂量率值与能量分布的影响情况,获得了准直照射时伽玛辐射场照射室尺寸的边界条件,建立并完善了伽玛参考辐射场边界研究方法及相关标准细节,为准直照射状态下照射室尺寸设计提供了一种新方法或途径。     

Analysis of the method of generalized separation of variables in the problem of light scattering by small axisymmetric particles

In the problem of light scattering by small axisymmetric particles, we have constructed the Rayleigh approximation in which the polarizability of particles is determined by the generalized separation of variables method (GSVM). In this case, electric-field strengths are gradients of scalar potentials, which are represented as expansions in term of eigenfunctions of the Laplace operator in the spherical coordinate system. By virtue of the fact that the separation of variables in the boundary conditions is incomplete, the initial problem is reduced to infinite systems of linear algebraic equations (ISLAEs) with respect to unknown expansion coefficients. We have examined the asymptotic behavior of ISLAE elements at large values of indices. It has been shown that the necessary condition of the solvability of the ISLAE coincides with the condition of correct application of the extended boundary conditions method (ЕВСМ). We have performed numerical calculations for Chebyshev particles with one maximum (also known as Pascal’s snails or limaçons of Pascal). The obtained numerical results for the asymptotics of ISLAE elements and for the matrix support theoretical inferences. We have shown that the scattering and absorption cross sections of examined particles can be calculated in a wide range of variation of parameters with an error of about 1–2% using the spheroidal model. This model is also applicable in the case in which the solvability condition of the ISLAE for nonconvex particles is violated; in this case, the SVM should be considered as an approximate method, which frequently ensures obtaining results with an error less than 0.1–0.5%.     

Two-dimensional linearly anisotropic scattering in a finite-thick cylindrical medium exposed to a laser beam

Integro-differential equations are developed for the source function, flux and intensity at the boundaries of a two-dimensional finite-thick medium which scatters in a linear fashion. The incident radiation is collimated, normal to the upper surface of the medium and dependent only on the radial coordinate. Two radial distributions are investigated: (1) a Bessel function and (2) a Gaussian laser beam. The solution for the Gaussian beam is constructed from the Bessel solution. Numerical results are presented in graphical and tabular forms for both boundary conditions. Comparisons are made between forward and isotropic scattering and between the finite and semi-infinite cases.     

Effect of band or line shape on radiative transfer in a nongray two-dimensional medium

An exact formulation is presented for a nongray two-dimensional, finite, planar, absorbing-emitting medium in radiative equilibrium. The absorption coefficient consists of an array of equal intensity, nonoverlapping bands or lines. Rectangular, triangular, exponential, Doppler and Lorentz shapes are specifically considered. Exact expressions are obtained for a medium subjected to collimated and diffuse radiation. The integral equations are linearized by the narrow-band approximation. The solution for the cosine-varying, collimated, monochromatic radiation model is used to construct the solutions for other boundary conditions. The two-dimensional equations are reduced to one-dimensional equations by the method of separation of variables. Results for the diffuse case are presented for several spatial variations.     

Pure-triplet scattering for radiative transfer in binary discrete random finite media with reflective boundaries and internal source of radiation

The stochastic solution of the monoenergetic radiative transfer equation in a finite slab random medium with pure-triplet anisotropic scattering is considered. The random medium is assumed to consist of two randomly mixed immiscible fluids labelled by 1 and 2. The extinction function, the scattering kernel, and the internal source of radiation are treated as discrete random variables, which obey the same statistics. The theoretical model used here for stochastic media transport assumes Markovian processes and exponential chord length statistics. The boundaries of the medium under consideration are considered to have specular and diffuse reflectivities with an internal source of radiation inside the medium. The ensemble-average partial heat fluxes are obtained in terms of the average albedos of the corresponding source-free problem, whose solution is obtained by using the Pomraning-Eddington approximation. Numerical results are calculated for the average forward and backward partial heat fluxes for different values of the single scattering albedo with variation of the parameters that characterize the random medium. Compared to the results obtained by Adams et al. in the case of isotropic scattering based on the Monte Carlo technique, it can be demonstrated that we have good comparable data.     

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.     

On the uniqueness of the problem of acoustic diffraction by an infinite plate with local irregularities

The question concerning the uniqueness of the solution to the problem of the acoustic diffraction by an immersed and isolated thin infinite plate with a finite scatterer is studied. It is shown that, to provide the uniqueness of the solution, the conditions at the scatterer must lead to an energy inequality for a source-free field, which determines the absence of the energy-carrying field components at infinity. A formula that generalizes the Sommerfeld formula is obtained and is used to prove the uniqueness of the solution to the problem of diffraction by a plate immersed in an acoustic medium. For the problem of diffraction of a flexural wave by an irregularity of the plate, the uniqueness theorem is proved only for the case of a fixed or hinged edge. When boundary conditions of a general form are imposed on the scatterer in an isolated plate, the uniqueness of the solution is generally lost, which is also corroborated by an example.     

Vorticity conditioning in the computation of two-dimensional viscous flows

The problem of evaluating the boundary values of the vorticity in the calculation of two-dimensional viscous flows is considered. It is shown that the splitting of the fourth-order equation for the stream function into two second-order problems implies specific integral conditions which fix the abstract projection of the vorticity field with respect, to the linear manifold of the harmonic functions. These conditions are a direct consequence of the boundary conditions on the velocity, and ensure satisfaction of physically essential conservation laws for the vorticity. The discrete analogue of, the projection conditions produces as many algebraic equations as the number of boundary points and requires the solution of an equal number of Dirichlet problems. In the particular case of stationary linearized equations (Stokes equations) a direct, i.e., noniterative method of solution is obtained. Steady and unsteady computational schemes relying on the projection conditions on the vorticity are introduced and extensive numerical results of finite difference calculations of the driven-cavity model problem are reported and discussed.