轴对称介质内辐射源项反演算法

结合Abel变换和离散坐标法,提出了一种基于CCD相机采集的单幅辐射图像重建轴对称发射-吸收介质内辐射源项分布的反演算法。通过在求解辐射正问题得到的准确值的基础上,添加随机噪声模拟试验测量数据,分析了网格数目、辐射源项分布形式、吸收系数和测量误差对算法反演精度的影响。测试结果表明:该算法对测量误差不敏感,在有测量误差的情况下也能够准确的重建介质内的辐射源项分布。     

An inverse radiative transfer problem of simultaneously estimating profiles of temperature and radiative parameters from boundary intensity and temperature measurements

A parallel-plane space filled with absorbing, emitting, isotropically scattering, gray medium is studied in this paper. The boundary intensity and boundary temperature profiles are calculated for the inverse analysis. For the simultaneous estimation of temperature, absorption and scattering coefficient profiles in the medium, the sum of residuals of boundary intensity and temperature after being weighted by a balance factor is minimized through using a Newton-type iteration algorithm and the least-squares method. To avoid over-updating for the parameters, the relative updating magnitude during the iteration process is constrained not to be >0.5. It is shown that the boundary intensity measurement alone is not enough to estimate simultaneously the temperature (source) and the radiative properties (both absorption and scattering coefficients) when the measurement data contain sensitive random errors. The boundary temperature measurement can serve as a necessary supplementation to the boundary intensity to make this kind of inverse radiative transfer problem resolvable. It was shown that a compensation relationship between absorption and scattering coefficients makes it difficult to fix them accurately. Parabolic profiles for the three parameters are used to validate the estimation method. When the optical thickness approaches 4.0, the results for the radiative properties are not acceptable, although the result for temperature profile is reasonable. This means the method needs further improvements.     

A sensitivity function-based conjugate gradient method for optical tomography with the frequency-domain equation of radiative transfer

The Sensitivity Function-based Conjugate Gradient Method (SFCGM) is described. This method is used to solve the inverse problems of function estimation, such as the local maps of absorption and scattering coefficients, as applied to optical tomography for biomedical imaging. A highly scattering, absorbing, non-reflecting, non-emitting medium is considered here and simultaneous reconstructions of absorption and scattering coefficients inside the test medium are achieved with the proposed optimization technique, by using the exit intensity measured at boundary surfaces. The forward problem is solved with a discrete-ordinates finite-difference method on the framework of the frequency-domain full equation of radiative transfer. The modulation frequency is set to 600 MHz and the frequency data, obtained with the source modulation, is used as the input data. The inversion results demonstrate that the SFCGM can retrieve simultaneously the spatial distributions of optical properties inside the medium within a reasonable accuracy, by significantly reducing a cross-talk between inter-parameters. It is also observed that the closer-to-detector objects are better retrieved.     

The DRESOR method for radiative heat transfer in a one-dimensional medium with variable refractive index

This paper extends the DRESOR (Distribution of Ratios of Energy Scattered by the medium Or Reflected by the boundary surface) method to radiative transfer in a variable refractive index medium. In this method, the intensity is obtained from the source term along the curved integration paths determined only by the variable refractive index, and the DRESOR values are calculated by the Monte Carlo method in which the propagation of the energy bundles are affected by Snell's law. With given temperatures on the black boundaries of a one-dimensional medium, the temperature distribution inside the medium with a variable scattering property is calculated under the condition of radiative equilibrium. It is shown that the DRESOR method has a good accuracy in the cases studied. For an isotropic-scattering medium with the same optical thickness, the scattering albedo has no effect on the temperature distribution, which can be obtained from the general equations and can be seen as an extension of what exists for a constant refractive index; however, the different refractive index causes obvious changes in the temperatures inside the medium. The effect of anisotropic scattering on the temperature distribution cannot be ignored, although it is still weaker than the effect caused by variation in the refractive index.     

Temperature field of radiative equilibrium in a semitransparent slab with a linear refractive index and gray walls

The temperature field in a semitransparent slab of absorbing-emitting gray medium at radiative equilibrium is solved in this paper. The medium has a linear refractive index and the two boundaries are diffuse gray walls. A curved ray tracing technique is combined with a pseudo-source adding method to deduce the radiative intensities on the gray walls. And on the basis of the previous work done by Ben Abdallah and Le Dez, the discrete temperature field in the slab is deduced. The influences of refractive index distribution, boundary wall emissivities and optical thickness on the radiative equilibrium temperature field are examined. The results display the significant influences of the refractive index distribution and the boundary wall emissivities.     

Temperature field of radiative equilibrium in a two-dimensional graded index medium with gray boundaries

In this paper, a numerical method is presented for the study of the radiative transfer in a two-dimensional graded index semitransparent medium with diffuse gray boundaries. The numerical method is a combination of the linear refractive index bar model, the discrete curved ray-tracing technique and the pseudo source adding method (LRIB-CRTP). In the traditional ray-tracing technique, it is difficult to deal with the diffuse gray boundary while solving the radiative transfer. Using the pseudo source adding method, the diffuse gray boundary of the medium can be treated as a black boundary. We have also studied the radiative equilibrium temperature field of the medium and analyzed the influence of some parameters involved. The results show that the directional discrete number is important for the medium having a large absorption coefficient. The results also show that the refractive index distribution greatly influences the temperature field, whereas the linear absorption coefficient distribution has little influence on the temperature field.     

Efficient inverse radiation analysis in a cylindrical geometry using a combined method of hybrid genetic algorithm and finite-difference Newton method

An inverse radiation problem was considered to estimate boundary conditions such as temperature distribution and emissivity in axisymmetric absorbing, emitting, and scattering medium, given the measured incident radiative heat fluxes. The finite-volume method was employed to solve a direct radiative transfer equation for a two-dimensional axisymmetric geometry. Various parameter estimators, such as conjugate-gradient method, hybrid genetic algorithm, and finite-difference Newton method, were employed to solve the inverse problems, while discussing their performances in terms of estimation accuracy and computational efficiency. Based on this, we proposed, as a best inverse analysis tool, a new combined method that adopted the hybrid genetic algorithm as an initial value selector and used the finite-difference Newton method as a parameter estimator.     

吸收散射性三维矩形介质内辐射源项的反问题

提出了一种由边界出射辐射强度反演吸收散射性三维矩形介质内辐射源项分布的方法。该方法是在辐射传递方程离散坐标近似的基础上,用求目标函数极小值的共轭梯度法进行反演计算。通过对介质辐射特性、光学厚度等参数对反演精度影响的分析,结果表明,即使存在测量误差,本文所提出的方法可较精确地反演辐射源项。     

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.     

Radiative intensity solution and thermal emission analysis of a semitransparent medium layer with a sinusoidal refractive index

The radiative intensity in a sinusoidal refractive index semitransparent medium layer is solved by the curved ray-tracing method in combination with the pseudo-source adding method. One boundary of the medium layer is an opaque diffuse substrate wall. The other boundary is a semitransparent specular or diffuse surface, from which the medium thermal emission emerges. With considering a linear temperature distribution, the radiative intensity formulae are, respectively, deduced under the two boundary conditions. On the basis of the radiative intensity solutions, the directional and hemispherical emission of the medium layer with a specular surface as well as the hemispherical emission of that with a diffuse surface are calculated. The influences of the optical thickness, sinusoidal refractive index distribution and linear temperature distribution on the thermal emission are investigated. The results show that the effects of refractive index and temperature distribution are significant and are different under the two reflecting modes of the surface.     

An inverse analysis of a transient 2-D conduction-radiation problem using the lattice Boltzmann method and the finite volume method coupled with the genetic algorithm

This article deals with the simultaneous estimation of parameters in a 2-D transient conduction-radiation heat transfer problem. The homogeneous medium is assumed to be absorbing, emitting and scattering. The boundaries of the enclosure are diffuse gray. Three parameters, viz. the scattering albedo, the conduction-radiation parameter and the boundary emissivity, are simultaneously estimated by the inverse method involving the lattice Boltzmann method (LBM) and the finite volume method (FVM) in conjunction with the genetic algorithm (GA). In the direct method, the FVM is used for computing the radiative information while the LBM is used to solve the energy equation. The temperature field obtained in the direct method is used in the inverse method for simultaneous estimation of unknown parameters using the LBM-FVM and the GA. The LBM-FVM-GA combination has been found to accurately predict the unknown parameters.     

Two analytic methods applied to radiative transfer in the linear-anisotropic scattering medium with graded index and diffuse gray boundaries

The curved ray-tracing method is extended to radiative transfer in the graded index medium with diffuse gray boundary conditions instead of black boundary conditions and the pseudo-source adding method is extended to the case of the linear-anisotropic scattering medium with graded index from non-scattering medium. Furthermore, the equivalence of the two methods is verified by formulation derivation. As exact analytical solutions, both the methods have high accuracy and fast computational speed. The predicted temperature distributions and dimensionless radiative heat flux at radiative equilibrium are determined by the proposed methods, and the numerical results are compared with the data in references. The results show that the present methods have a good accuracy. Influences of various combinations of refractive index and boundary emissivities on the temperature distributions and dimensionless radiative heat flux are also investigated.     

A Polydisperse Sphere Model Describing the Propagation of Light in Biological Tissue

A polydisperse sphere model with the complex refractive index is employed to describe the propagation of light in biological tissue. The scattering coefficient, absorption coefficient and scattering phase function are calculated. At the same time, the inverse problem on retrieving the particles size distribution, imaginary part of the refractive index and number density of scatterers is investigated. The result shows that the retrieval scheme together with the Chahine algorithm is effective in dealing with such an inverse problem. It is also clarified that a group of parameters including the scattering coefficient, absorption coefficient and phase function are associated with another group including the refractive index, particle size distribution and number density of scatterers, which is a problem described in two different ways and the anisotropy factor is not an independent variable, but is determined by the phase function.     

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.     

Influence of the refractive index on the radiative source function of an isotropically scattering medium

The influence of refractive index on the radiative source function is presented for the case of emission from an isothermal, isotropically scattering medium. A closed-form, approximate solution is obtained and results are presented for finite and semi-infinite nonconservative cases. An increase in refractive index causes the source function to increase and that effect is more pronounced at higher scattering albedo and smaller optical depths.     

一维梯度折射率介质光学特性的反演

在单个方波脉冲入射情况下,利用共轭梯度法对一维梯度折射率介质的折射率、吸收系数以及散射系数进行了反演.正问题采用间断有限元法求解,反问题的解则在正问题的基础上通过共轭梯度法得到.研究结果表明,利用单个方波脉冲入射情况下的时域半球反射率以及时域半球透射率作为测试值能够有效地反演一维梯度折射率介质的光学特性及其分布情况.     

An inverse source problem in radiative transfer

The spherical-harmonics method is used to develop a solution to an inverse source problem in radiative transfer. It is assumed that, with the exception of the inhomogeneous source term, all aspects of the radiation-transport problem are known, and we seek to determine the inhomogeneous source term from specified angular distributions of radiation exiting the two surfaces of a homogeneous plane-parallel medium. Anisotropic scattering is included in the monochromatic radiative-transfer model and general reflecting boundary conditions are considered.     

Efficient inverse radiation analysis of temperature distribution in participating medium based on backward Monte Carlo method

An efficient numerical inverse radiation analysis based on the backward Monte Carlo (BMC) method is presented to determine the three-dimensional (3-D) temperature distribution in a large rectangular enclosure containing the participating medium, using radiative intensities in the visible range received by charge-coupled device (CCD) cameras. For large radiative sources and small detectors, when the radiation onto a small spot and onto a small direction cone is desired, the BMC method is more efficient than the forward Monte Carlo (FMC) method. Because the temperature reconstruction from the measured radiative intensities is an ill-posed inverse problem, least-square QR decomposition (LSQR) method is introduced to determine the meaningful temperature distribution. In order to gain insight into the effects on the accuracy of temperature distribution reconstruction, the detailed analyses are made using numerical simulations. The results show that the temperature distribution can be reconstructed accurately for the exact and noisy data.     

各向异性介质辐射特性参数联合反演

本文基于BP神经网络方法结合蒙特卡洛和BEER定律辐射传输模拟方法建立了联合反演各向异性散射介质的辐射特性参数模型。首先采用半球透射率结合半球反射率反演模型反演了各向同性介质的吸收系数和散射系数,在此基础上增加准直透射率,建立了联合反演各向异性介质的吸收系数、散射系数和散射不对称因子三参数联合反演模型。反演结果表明该模型能准确反演出介质辐射特性参数,具有实用意义。此外,为了检验测量误差对模型的反演准确性的影响,分别在不同程度测量误差情况下进行反演,结果显示测量误差对散射不对称因子反演值影响较大。     

Transmission of radiation through an aerosol medium

The properties of radiation through an aerosol medium have been achieved. This has been done by employing Mie scattering theory to calculate the radiation transfer scattering parameters in the form of extinction, absorption and scattering efficiencies. The equation of radiative transfer for the heat flux through a plane parallel atmosphere of aerosol has been solved. The aerosol size distributions are found in practical systems. Average efficiencies over size distribution for spherical particles of complex refractive index are calculated. Therefore, the radiative properties of stratospheric aerosols have been done. The obtained results found to be in a good agreement with the previous work.