Multiple scattering by particles embedded in an absorbing medium. 2. Radiative transfer equation

This paper continues a systematic theoretical analysis of electromagnetic scattering by a group of arbitrarily sized, shaped, and oriented particles embedded in an absorbing, homogeneous, isotropic, and unbounded medium. The previously developed microphysical approach is used to derive the generalized form of the radiative transfer equation (RTE) applicable to a large group of sparsely, randomly, and uniformly distributed particles. The derivation of the RTE directly from the macroscopic Maxwell equations yields unambiguous and definitive analytical expressions for the participating quantities and thereby fully resolves the lasting controversy caused by the conflicting outcomes of several phenomenological approaches.     

Radiative transfer in absorbing, emitting and linearly anisotropic-scattering inhomogeneous cylindrical medium

Radiative integral transfer equations for one-dimensional solid cylinder with absorbing, emitting and linearly anisotropic-scattering inhomogeneous medium were derived by Abulwafa et al. (JQSRT 62 (1999) 755). The anisotropic terms in the integral equations and their results for anisotropic benchmark problems (JQSRT 66 (2000) 487) are inaccurate. In this study, the integral equations for absorbing, emitting and linearly anisotropic-scattering medium are rederived, and the integral equations for one-dimensional solid cylindrical medium are solved. The results are compared with selected cases using the discrete ordinates S₁₆ and the exact solutions available in the literature.     

高斯光束中吸收双层球形微粒的横向光俘获

为了研究吸收双层球形微粒的横向光俘获,基于几何光学模型提出了双层带吸收球形微粒的光俘获模型,对TEM₀₀模式高斯光束照射下外层有光吸收的双层电介质球形微粒受到的横向光俘获力进行了数值模拟,取得了光俘获力特性的一系列结果.结果显示,双层球形微粒的外层吸收系数对包括稳态俘获位置,峰值强度,稳态俘获的刚度等光俘获特性有很大影响.此外,内外径的比率对吸收双层球形微粒的光俘获特性也有调制性的影响.在一定条件下,带吸收的双层球形微粒可以被俘获在光轴上,也可能被俘获在中心在光轴上的圆环上. 关键词： 光俘获 几何光学模型 高斯光束 吸收双层球     

Determination of optical properties of tissues

Recent results and used way of tissue optical properties measurements in vitro with laser goniophotometer are presented. The optical properties were obtained on the basis of an indirect technique using solution of the radiative transfer equation, which relate the optical parameters of tissue with the measured parameters of reflected and transmitted radiation. The angular patterns of reflected and transmitted radiation by optically thick saline-soaked samples of tumor tissue, namely, mammary fibroadenoma, were measured at the wavelength of 0.63 μm, and the measurement data were inverted to find the absorption and extinction coefficients and the parameters of the phase function of these tumor samples. Obtained values of the optical parameters have been found to be of the same order of magnitude as the earlier reported data for similar tumor samples.     

Geometrical optics approximation for light scattering by absorbing spherical particles

By means of geometrical optics, an approximation method is presented to compute the light scattering intensity of absorbing spherical particles illuminated by a plane wave. For absorbing particles, the effective refractive index and the effective refractive angle are related to the complex refractive index and incident angle. The formulas for calculation of the break of phases of reflection and refraction, which are different from the case of transparent particles, are exactly derived. Verification of the geometrical optics approximation (GOA) was performed by case studies and comparison of the present results with the Mie scattering. It is found that agreement between the GOA and the Mie theory is excellent in forward directions for weakly/moderately absorbing particles. Differently, for strongly absorbing particles, good agreement between the calculation methods is in the forward directions and large scattering angles. The agreement between the GOA and the Mie theory is better for larger particles.     

Noninvasive determination of tissue optical properties based on radiative transfer theory

We present a feasibility study of a new method for determining the tissue optical properties, including the absorption and scattering coefficients and the scattering asymmetry factor. A state-of-the-art radiative transfer model for the coupled air/tissue system, based on rigorous radiative transfer theory, is used in our forward modeling simulations. The concept of the effective photon penetration depth is introduced and used to help determine the depth below, which information about the tissue will not be available through noninvasive imaging of a biological tissue using reflected diffuse light. Simulation results show that for accurate determination of tissue optical properties, one can use radiative transfer theory in conjunction with measurements of reflected radiances as well as other existing techniques.     

Difference of optical properties between porous alumina and sapphire using two-substrate method at elevated temperature

We report on simultaneous measurements of the emittance, transmittance and reflectance of a porous alumina under same environment at 573 K, which has not been obtained. It is demonstrated that the correction of substrate radiation and multiple reflections is important for the accurate measurement of infrared (IR) optical properties of porous materials in a two-substrate method. We deduced the extinction, absorption and scattering coefficients of the alumina from the measured three quantities, based on the classical geometry model. The coefficients explain the difference of the IR optical properties between the alumina and the sapphire. The increasing reflectance of the alumina with thickness, unlike the sapphire, is attributed to the relatively high scattering coefficient at short wavelengths.     

Radiative transfer equation for graded index medium in cylindrical and spherical coordinate systems

In graded index medium, the ray goes along a curved path determined by Fermat principle, and the curved ray-tracing is very difficult and complex. To avoid the complicated and time-consuming computation of curved ray trajectory, the methods not based on ray-tracing technique need to be developed for the solution of radiative transfer in graded index medium. For this purpose, in this paper the streaming operator along a curved ray trajectory in original radiative transfer equation for graded index medium is transformed and expressed in spatial and angular ordinates and the radiative transfer equation for graded index medium in cylindrical and spherical coordinate systems are derived. The conservative and the non-conservative forms of radiative transfer equation for three-dimensional graded index medium are given, which can be used as base equations to develop the numerical simulation methods, such as finite volume method, discrete ordinates method, and finite element method, for radiative transfer in graded index medium in cylindrical and spherical coordinate systems.     

吸收媒质反射束Goos-Hnchen位移的直接计算公式

本文导出了从吸收媒质上反射的s和p两种偏振的Goos-Hnchen位移的一般表达式。若令消光系数等于零,它们就简化成McGuilk和Carniglia给出的仅用于非吸收媒质的公式。吸收媒质的位移,在任何不为零的入射角度时都有。在镓的表面上,位移有几个波长。     

杂质颗粒对受激布里渊散射介质光学击穿阈值影响的研究

建立了受激布里渊散射介质中热作用破坏的物理模型, 数值模拟了杂质颗粒的温度随其半径的变化曲线.结果显示,杂质颗粒存在一个最大热作用半径, 介质所含颗粒的尺寸在此半径附近时,介质最容易发生光学击穿现象, 其光学击穿阈值最低.在Continuum Nd: YAG种子注入式激光系统中,选取FC-3283, GF-180和HFE-7100介质, 通过不同孔径的过滤膜进行过滤,并研究了过滤前后的光学击穿阈值和能量反射率. 结果表明,随着过滤孔径的变小,介质光学击穿阈值逐渐提高, 且过滤之后介质的能量反射率有了明显的提高.介绍了一种利用He-Ne激光透射光光斑变化来判断是否发生光学击穿现象的方法,该方法具有方便、 准确的特点,可有效地减小由于肉眼观测引起的误差.     

Tau method approximation for radiative transfer problems in a slab medium

An approximate method for solving the radiative transfer equation in a slab medium with an isotropic scattering is proposed. The method is based upon constructing the double Legendre series to approximate the required solution using Legendre tau method. The differential and integral expressions which arise in the radiative transfer equation are converted into a system of linear algebraic equations which can be solved for the unknown coefficients. Numerical examples are included to demonstrate the validity and applicability of the method and a comparison is made with existing results.     

Curved ray tracing method for one-dimensional radiative transfer in the linear-anisotropic scattering medium with graded index

The curved ray tracing method (CRT) is extended to radiative transfer in the linear-anisotropic scattering medium with graded index from non-scattering medium. In this paper, the CRT is presented to solve one-dimensional radiative transfer in the linear-anisotropic scattering gray medium with a linear refractive index and two black boundaries. The predicted temperature distributions and radiative heat flux at radiative equilibrium are determined by the proposed method, and numerical results are compared with the data in references. The results show that the CRT has a good accuracy for radiative transfer in the linear-anisotropic scattering medium with graded index and the dimensionless emissive power and dimensionless radiative heat flux depend on the dimensionless refractive index gradient. It can also be seen that the dimensionless refractive index gradient has important effects on the temperature discontinuity at the boundaries.     

Stability in Debye series calculation for light scattering by absorbing particles and bubbles

The Debye-series decomposition is of importance for understanding of light scattering features and for the validity of the geometrical optics approximation to light scattering. The numerical stability and accuracy for calculating light scattering with Debye series is studied and an improved algorithm is proposed in this work. The ratios of the Riccati-Bessel functions and the logarithmic derivatives of the Riccati-Bessel functions are employed and calculated with proper recurrences. Exemplifying results are provided to show the improvement of the algorithm.     

Optical properties of single mixed-phase aerosol particles

We examine the light scattering from wet aerosol particles whose core morphology is a random Gaussian sphere coated by different amounts of water. We compare the results with those calculated from a concentric-sphere particle and from a homogeneous-sphere particle whose refractive index is determined from the Bruggeman mixing rule. In nearly all cases the differences between those of the Gaussian-core and the concentric-sphere particles are small. The most significant differences are seen in calculations of the asymmetry parameters and the intensity and polarization phase functions. The results of the homogeneous-sphere particles vary significantly from those of the Gaussian-core particles, typically an order of magnitude greater than for the concentric-sphere particles. It is not uncommon to see differences of 10% in the efficiencies of large, homogeneous-sphere particles; whereas, in the intensity and polarization phase functions, differences of several tens of percent are not uncommon.     

Three-dimensional radiative transfer in an anisotropically scattering, semi-infinite medium: generalized reflection function

Three-dimensional radiative transfer in an anisotropic scattering medium exposed to spatially varying, collimated radiation is studied. The generalized reflection function for a semi-infinite medium with a very general scattering phase function is the focus of this investigation. An integral transform is used to reduce the three-dimensional transport equation to a one-dimensional form, and a modified Ambarzumian's method is applied to formulate a nonlinear integral equation for the generalized reflection function. The integration is over both the polar and azimuthal angles; hence, the integral equation is said to be in the double-integral form. The double-integral, reflection function formulation can handle a variety of anisotropic phase functions and does not require an expansion of the phase function in a Legendre polynomial series. Complicated kernel transformations of previous single-integral studies are eliminated. Single and double scattering approximations are developed. Numerical results are presented for a Rayleigh phase function to illustrate the computational characteristics of the method and are compared to results obtained with the single-integral method. Agreement between the two approaches is excellent; however, as the transform variable increases beyond five the number of quadrature points required for the double-integral method to produce accurate solutions significantly increases. A new interpolation scheme produces accurate results when the transform variable is large.     

Three-dimensional radiative transfer in an isotropically scattering, semi-infinite medium: generalized reflection function

The focus of this study is the generalized reflection function of multidimensional radiative transfer. The physical situation considered is spatially varying, collimated radiation incident on the upper boundary of an isotropically scattering, semi-infinite 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 applied to formulate a nonlinear integral equation for the generalized reflection function. The resulting equation is said to be in double-integral form because the integration is over both angular variables. Computational issues associated with this generalized reflection function formulation are investigated. The source function and reflection function formulations are compared, and the relative merits of the two approaches are discussed.     

超短脉冲在光学介质中传输的B积分特性

从非线性薛定谔方程出发,利用分步傅里叶算法对其进行求解,数值模拟了在超短脉冲激光系统中几种光学介质中的B积分传输规律,并对其特性进行了简要分析。数值计算结果表明,初始脉冲的输入光强与光学介质的增益系数的增大会使B积分值增加。脉冲形状对B积分值有一定影响,与飞秒高斯脉冲比较,当输入为ps级的啁啾脉冲时,B积分较小。针对所选计算模型,适当的群速色散与高阶色散可以减小B积分值。这对超短脉冲放大系统的设计有参考意义。     

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.     

Non-Fourier heat conduction with radiation in an absorbing, emitting, and isotropically scattering medium

This article numerically analyses the combined conductive and radiative heat transfer in an absorbing, emitting, and isotropically scattering medium. The non-Fourier heat conduction equation, which includes the time lag between heat flux and the temperature gradient, is used to model the conductive heat transfer in the medium. It predicts that a temperature disturbance will propagate as a wave at finite speed. The radiative heat transfer is solved using the P₃ approximation method. In addition, the MacCormack's explicit predictor-corrector scheme is used to solve the non-Fourier problem. The effects of radiation including single scattering albedo, conduction-to-radiation parameter, and optical thickness of the medium on the transient and steady state temperature distributions are investigated in detail. Analysis results indicate that the internal radiation in the medium significantly influences the wave nature. The thermal wave nature in the combined non-Fourier heat conduction with radiation is more obvious for large values of conduction-to-radiation parameter, small values of optical thickness and higher scattering medium. The results from non-Fourier-effect equation are also compared to those obtained from the Fourier equation. Non-Fourier effect becomes insignificant as either time increases or the effect of radiation increases.