多维梯度折射率介质内辐射传递的扩散近似

推导多维梯度折射率介质内稳态辐射传递的扩散近似方程.使用有限元法对扩散近似进行离散和求解,利用两个二维半透明介质的稳态辐射传递问题验证该扩散近似的精度及适用性.算例考虑介质为均匀折射率及梯度折射率两种情况.利用扩散近似分别求解辐射平衡时的边界热流、介质内温度场分布,并与辐射传递方程的求解结果进行对比分析.结果表明:介质折射率变化、散射特性、光学厚度及散射反照率均直接影响扩散近似的精度;在光学厚及强散射条件下,该扩散近似可以作为一种快速算法应用于梯度折射率介质稳态辐射传递的求解.     

Finite element approximation of the Fokker-Planck equation for diffuse optical tomography

In diffuse optical tomography, light transport theory is used to describe photon propagation inside turbid medium. A commonly used simplification for the radiative transport equation is the diffusion approximation due to computational feasibility. However, it is known that the diffusion approximation is not valid close to the sources and boundary and in low-scattering regions. Fokker-Planck equation describes light propagation when scattering is forward-peaked. In this article a numerical solution of the Fokker-Planck equation using finite element method is developed. Approach is validated against Monte Carlo simulation and compared with the diffusion approximation. The results show that the Fokker-Planck equation gives equal or better results than the diffusion approximation on the boundary of a homogeneous medium and in turbid medium containing a low-scattering region when scattering is forward-peaked.     

Image reconstruction in diffuse optical tomography using the coupled radiative transport-diffusion model

The coupled radiative transport-diffusion model can be used as light transport model in situations in which the diffusion equation is not a valid approximation everywhere in the domain. In the coupled model, light propagation is modelled with the radiative transport equation in sub-domains in which the approximations of the diffusion equation are not valid, such as within low-scattering regions, and the diffusion approximation is used elsewhere in the domain. In this paper, an image reconstruction method for diffuse optical tomography based on using the coupled radiative transport-diffusion model is developed. In the approach, absorption and scattering distributions are estimated by minimising a regularised least-squares error between the measured data and solution of the coupled model. The approach is tested with simulations. Reconstructions from different cases including domains with low-scattering regions are shown. The results show that the coupled radiative transport-diffusion model can be utilised in image reconstruction problem of diffuse optical tomography and that it produces as good quality reconstructions as the full radiative transport equation also in the presence of low-scattering regions.     

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.     

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

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

Generalized diffusion equation for multiple scattered elastic waves

We study multiple scattering of elastic waves using a generalized diffusion approximation of the radiative transfer equation applied to spherically symmetric scatterers. This generalized diffusion equation allows us to keep track of the two elastic wave types as well as the mode conversions. It describes the process towards equipartition and fills, as such, a gap between the radiative transfer equation and the conventional diffusion approximation. The effects of boundary conditions and dissipation on the energy partitioning are studied.     

Adjoint problem of radiative transfer for a pseudo-spherical atmosphere and general viewing geometries

We formulate the adjoint radiative transfer for a pseudo-spherical atmosphere and various retrieval scenarios. The single scattering radiance is computed in a spherical atmosphere by using the source integration technique, while for the multiple scattering radiance we formulate an one-dimensional adjoint radiative transfer equation in a plane-parallel atmosphere. The adjoint solution of the radiative transfer equation is obtained by employing the discrete ordinate method with matrix exponential. We provide an abbreviated derivation of our formalism as well as a discussion of the numerical implementation of the theory.     

Time-dependent equation for the intensity in the diffusion limit using a higher-order angular expansion

The first two terms in the spherical-harmonic expansion (the P(1) approximation) of the radiative transfer equation yield the diffusion equation. This approximation applies to multiple scattering and results in a solution for the energy density, the gradient of which is proportional to the light intensity. In this work a higher-order spherical-harmonic expansion of the radiative transfer equation is developed. This equation applies to the radiant intensity rather than the energy density. The equation can be decomposed into two terms: a propagator term obtained from the determinant of the coupled equations describing the individual components of the intensity, and a mixing matrix that describes the cross coupling between different orders of the expansion. Using the Fourier transform, an approximation based on expanding at small wave vectors k leads to an equation similar to the diffusion equation. The equation is expected to predict the intensity for multiple scattering at earlier times and shorter distances than the diffusion equation can. The notion of an equivalent wave field is introduced.     

Comparing Diffusion Approximation with Radiation Transfer Analysis for Light Transport in Tissues

The diffusion model that is an approximation of the equation of radiation transfer is typically used to describe photon migration in scattering-dominant media. In general biological tissue is highly scattering and very weakly absorbing against near-infrared light, yet it is heterogeneous and may contain relatively highly absorbing or low-scattering regions. Here applicability of the diffusion approximation over the radiative transfer theory for describing ultrafast laser transport in biological tissues is numerically studied and investigated over different kinds of tissue conditions and geometries. Tissues having tumors of different sizes, locations and nature as well as dual-tumor and low-scattering conditions are considered. Radiation transfer analysis is taken as a comparison objective and it is initially proved to be accurate in benchmark comparisons with Monte Carlo simulation. The results predict systematically about the compatible conditions where and when we can use the diffusion approximation and the conditions in which the diffusion approximation may provide misleading results.     

Decorrelation and phase-shift of coda waves induced by local changes: multiple scattering approach and numerical validation

We report on theoretical predictions of the decorrelation and phase-shift of coda waves induced by local changes in multiple scattering media. Using the multiple scattering formalism, we show that both expressions (decorrelation and phase-shift) involve a same sensitivity kernel based on the intensity transport in the medium. We compare the kernels based on the diffusion approximation with the ones based on the radiative transfer approximation, showing that the latter is more accurate at short times or for changes located close to the source or the receiver. We also perform a series of numerical simulations of wave propagation (finite differences) to validate our models in different configurations.     

Analysis of combined conduction and radiation heat transfer in presence of participating medium by the development of hybrid method

The current study addresses the mathematical modeling aspects of coupled conductive and radiative heat transfer in the presence of absorbing, emitting and isotropic scattering gray medium within two-dimensional square enclosure. A blended method where the concepts of modified differential approximation employed by combining discrete ordinate method and spherical harmonics method, has been developed for modeling the radiative transport equation. The gray participating medium is bounded by isothermal walls of two-dimensional enclosure which are considered to be opaque, diffuse and gray. The effect of various influencing parameters i.e., radiation-conduction parameter, surface emissivity, single scattering albedo and optical thickness has been illustrated. The adaptability of the present method has also been addressed.     

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.     

The surface cooling effect in radiative transfer

The halfspace radiative equilibrium problem with an incident beam of energy is analyzed by both the equation of transfer and the diffusion approximation. The spatial dependence of the source function is shown to exhibit (unexpected) surface cooling for all angles of incidence. For grazing incidence, the source function additionally exhibits a maximum within the halfspace. The diffusion results agree qualitatively with these transport findings for near normal incidence, but not for grazing incidence. The relationship of this problem to one in neutron transport theory is established, graphical results are presented, and the relevance of this analysis to understanding detailed radiation-hydrodynamic computer results for a certain class of problems is discussed.     

Free-Space Propagation of the Generalized Radiance Defined in Terms of the Coherent-Mode Representation of Cross-Spectral Density Function

The equation for free-space propagation of the generalized radiance defined formerly by the authors in terms of the coherent-mode representation of the cross-spectral density function is derived within the accuracy of the paraxial approximation. It is shown that, in the short-wavelength limit, this generalized radiance obeys in a good approximation the fundamental radiative transfer law of classical radiometry.     

Analysis of the radiative transfer equation with highly assymetric phase function

This paper considers a scalar radiative transfer problem with high scattering anisotropy. Two computational methods are presented based on decomposition of the diffuse light field into a regular and anisotropic part. The first algorithm (DOMAS) singles out the anisotropic radiance in the forward scattering peak using the Small-Angle Modification of RTE. The second algorithm (DOM2+) separates the single scattering radiance as an anisotropic part, which largely defines the fine detail of the total radiance in the backscattering directions. In both cases, the anisotropic part is represented analytically. With anisotropy subtraction, the regular part of the signal, which requires a numerical solution, is essentially smoothed as a function of angles. Further, the transport equation is obtained for the regular part that contains an additional source function from the anisotropic part of the signal. This equation is solved with the discrete ordinates method. A conducted numerical analysis of this work showed that algorithm DOMAS has a strong advantage as compared to the standard discrete ordinates method for simulation of the radiance transmission, and DOM2+ is the best of the three for the reflection computations. Both algorithms offer at least a factor of three acceleration of convergence of the azimuthal series for highly anisotropic phase functions.     

A single-scattering approximation for infrared radiative transfer in limb geometry in the Martian atmosphere

We present a single-scattering approximation for infrared radiative transfer in limb geometry in the Martian atmosphere. It is based on the assumption that the upwelling internal radiation field is dominated by a surface with a uniform brightness temperature. It allows the calculation of the scattering source function for individual aerosol types, mixtures of aerosol types, and mixtures of gas and aerosol. The approximation can be applied in a Curtis-Godson radiative transfer code and is used for operational retrievals from Mars Climate Sounder measurements. Radiance comparisons with a multiple scattering model show good agreement in the mid- and far-infrared although the approximate model tends to underestimate the radiances in realistic conditions of the Martian atmosphere. Relative radiance differences are found to be about 2% in the lowermost atmosphere, increasing to ∼10% in the middle atmosphere of Mars. The increasing differences with altitude are mostly due to the increasing contribution to limb radiance of scattering relative to emission at the colder, higher atmospheric levels. This effect becomes smaller toward longer wavelengths at typical Martian temperatures. The relative radiance differences are expected to produce systematic errors of similar magnitude in retrieved opacity profiles.     

短脉冲激光与介质相互作用的扩散模型

推导了扩散近似方程,通过半无限大均匀介质计算,用扩散理论分析解验证了数值方法的有效性.模拟了光在非均匀介质内的传输过程,给出了介质内光通量随时间变化的空间分布.结果表明,该基于扩散模型的数值方法能够模拟短脉冲光在强散射介质中的传播过程以及漫散射光的时间变化特性,并且借助于光通量空间分布能够准确模拟非均匀介质内内含物的位置.     

An Inverse Problem of Determining Coefficients in a One-Dimensional Radiative Transport Equation

We consider an inverse problem of determining unknown coefficients for a one-dimensional analogue of radiative transport equation. We show that some combination of the unknown coefficients can be uniquely determined by giving pulse-like inputs at the boundary and observing the corresponding outputs. Our result can be applied for determination of absorption and scattering properties of an optically turbid medium if the radiative transport equation is appropriate for describing the propagation of light in the medium.     

用双点源δ-P1近似光学模型反演生物组织的光学参量

根据空间分辨漫反射的双点源δ-P₁近似理论模型,采用非线性最小二乘法,从反射率的测量数据中反演得到了生物组织的吸收系数μ_a、有效散射系数μ' _s和二阶参量γ.研究表明,在光源与探测器之间距离大于一个输运平均自由程的情况下,双点源δ-P₁近似能较好地描述光源附近的光辐射分布,而且能够根据参量γ与μ' _s的关系得到组织的各向异性因子g.这些研究对于生物组织的光学性质测量以及漫反射光谱技术的应用具有重要意义. 关键词： 组织光学 P₁近似')" href="#">δ-P₁近似 光学参量 双点源     

空间分辨漫反射的一阶散射参量灵敏度

以输运理论的P3近似为基础,推导了空间分辨漫反射的一阶散射参量灵敏度的解析表示,并进行了数值分析,比较了散射参量对P3近似和漫射近似漫反射的影响.结果表明:在距光源两个输运平均自由程以内,该灵敏度与相应于漫射近似的灵敏度差别较为明显;距离光源约四个输运平均自由程附近,该灵敏度等于零,并且与光源之距与相应于漫射近似情况也...