求解辐射传递方程的离散坐标法

阐述了含吸收散射性介质三维空腔内辐射传递方程的离散坐标解法。讨论了入射散射项积分格式的选取,以及假散射和射线效应对解精度的影响。对三维矩形炉膛内辐射传递过程进行了数值模拟,并与区域法和离散传递法进行比较。比较结果表明离散坐标法具有较好的精度,是目前燃烧室内辐射传热过程数值模拟的一种较好的方法。     

Propagation of a light beam in an absorbing medium with large-scale inhomogeneities

A method of solving the radiative transfer equation is proposed; it enables one to take into account the influence of absorption on the angular and spatial distributions of radiation under conditions of sharply anisotropic multiple scattering. For phase functions that decrease with an increase in the scattering angle by the power law, the total flux attenuation and profiles of the angular and spatial distributions in a strongly absorbing medium are studied. The obtained analytical dependences exhibit a good agreement with results of numerical solution of the radiative transfer equation.     

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.     

Radiation transport and internal reflection in a sphere

We construct an integral equation for the flux intensity in a scattering and absorbing medium using the integro-differential form of the radiative transfer equation in a sphere. The sphere is uniformly irradiated by an external source of arbitrary angular distribution. The Fresnel boundary conditions, which incorporate reflection and refraction, are used. For the special cases of a non-scattering medium, and in the limit of an optically transparent medium, we obtain exact solutions for specular and diffuse refection. Some numerical examples are given which give qualitative agreement with some recent work of Tian and Chiu (JQSRT, 2005).     

Multidimensional radiative transfer in absorbing,emitting, and linearly anisotropic scattering cylindrical medium with space-dependent properties

The integral form of three-dimensional radiative transfer equation for an absorbing, emitting, and linear-anisotropic scattering medium with space-dependent properties is formulated. A product-integration method is subsequently applied to develop a numerical scheme for solving the corresponding integral transfer equations in a two-dimensional, axisymmetric and nonhomogeneous medium subjected to externally incident radiation or bounded by emitting and diffusely-reflecting walls. The numerical solutions for cases of constant, continuous, and stepwise variations of scattering albedo are presented to illustrate its accuracy and flexibility, and validated by comparing with results available in the literature.     

Radiative properties of densely packed spheres in semitransparent media: A new geometric optics approach

This contribution presents a new Ray-tracing method for calculating effective radiative properties of densely packed spheres in non-absorbing or semitransparent host medium. The method is restricted to the geometric optic objects and neglects the wave effects. The effective radiative properties such as the absorption and scattering coefficients, and phase function are retrieved from the calculation of mean-free paths of scattering and absorption, and the angular scattering probability of radiation propagating in the dispersed medium. The model accounts for the two geometric effects called here as non-point scattering and ray transportation effects. The successful comparison of the current model with data of radiative properties and transmittances of particle beds in a non-absorbing medium reported in the literature confirm its suitability. It is shown that: (i) for opaque or absorbing particles (not systematically opaque), the non-point scattering is the dominant geometric effects whereas both non-point scattering and ray transportation effects occur for weakly absorbing and transparent particles. In the later cases, these two geometric effects oppose and may cancel out. This may explain why the Independent scattering theory works well for packed of quasi-transparent particles; (ii) the non-point scattering and ray transportation effects can be captured through the scattering and absorption coefficients while using the classical form of phase function. This enables using the standard radiative transfer equation (RTE); (iii) the surrounding medium absorption can be accounted for without any homogenization rule. It contributes to increasing the effective absorption coefficient of the composite medium as expected but, at the same time, it reduces the particle extinction; and (iv) the current transfer calculation predicts remarkably the results of direct Monte Carlo (MC) simulation. This study tends therefore to confirm that the RTE can be applied to densely packed media by using effective radiative properties.     

漫反射各向异性散射介质内的温度热流场

本文采用射线踪迹结合节点分析法和谱带模型,研究了漫反射不透明边界下吸收、发射、各向异性散射介质内的热辐射传递过程。考虑介质辐射能的入射和散射方向,导出漫反射、不透明边界、各向异性散射介质的辐射传递系数。在辐射平衡的情况下,考察了表面发射率和散射反照率对介质内辐射热流和温度场的影响。研究表明,介质不透明边界处存在温度跃迁现象,而且,内界面发射率越大,相应界面温度跃迁越小。     

Non-gray radiative transfer in plane-parallel media with reflecting boundaries

A generalized equation of radiative transfer in the two-group picket-fence model is analyzed for a plane parallel, emitting, absorbing and isotropically scattering medium containing uniform heat sources and having boundary surfaces which are diffuse emitters and diffuse reflectors and are maintained at uniform but arbitrary temperatures. The solution of the general problem is expressed by the superposition of simpler problems which are solved by the application of the normal-mode-expansion technique. Highly accurate numerical results are presented for the temperature distribution and the radiative heat flux in the medium.     

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.     

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

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

DRESOR法对瞬态辐射传递问题的研究

用基于蒙特卡洛法(Monte Carlo Method,MCM)的DRESOR法(Distributions of Ratios of Energy Scattered by the medium Or Reflected by the boundary surface)求解入射辐射经过介质散射、壁面反射传递后辐射强度随时间变化的瞬态辐射传递方程(Transient RadiativeTransfer Equation,TRTE)问题。通过在系统内计算一单位瞬态入射辐射对介质的DRESOR数分布,就能计算任意时间内入射辐射在系统内时间响应特性,这样有效提高数值方法处理瞬态辐射问题的通用性。并且能够获得高方向分辨率的辐射强度随时间变化的结果,这是目前大多数数值处理方法比较难做到的,显示出了DRESOR法处理瞬态入射辐射问题的能力．     

Prediction of radiative heat transfer in 2D irregular geometries using the collocation spectral method based on body-fitted coordinates

Recently, an efficient numerical method, which is called the collocation spectral method (CSM), for radiative heat transfer problems, has been proposed by the present authors. In this numerical method there exists the exponential convergence rate, which can obtain a very high accuracy even using a small number of grids. In this article, the CSM based on body-fitted coordinates (BFC) is extended to simulate radiative heat transfer problems in participating medium confined in 2D complex geometries. This numerical method makes simultaneously the use of the merits of both the CSM and BFC. In this numerical approach, the radiative transfer equation (RTE) in orthogonal Cartesian coordinates should be transformed into the equation in body-fitted nonorthogonal curvilinear coordinates. In order to test the efficiency of the developed method, several 2D complex irregular enclosures with curved boundaries and containing an absorbing, emitting and scattering medium are examined. The results obtained by the CSM are assessed by comparing the predictions with those in references. These comparisons indicate that the CSM based on BFC can be recommended as a good option to solve radiative heat transfer problems in complex geometries.     

Discrete transfer method applied to radiative transfer in a variable refractive index semitransparent medium

Application of the discrete transfer method (DTM) has been extended to the analysis of radiative heat transfer in a variable refractive index participating medium. To validate the DTM formulation, radiative heat transfer in an absorbing, emitting and isotropically scattering planar medium was considered. The participating medium was assumed to be in radiative equilibrium. For both constant and variable refractive indices of the medium, the DTM results were compared with those available in the literature. The DTM was found to provide accurate results.     

Thermal analysis for transient radiative cooling of a conducting semitransparent layer of ceramic in high-temperature applications

A method is developed for obtaining transient temperature distribution in a cooling semitransparent layer of ceramic. The layer is emitting, absorbing, isotropically scattering and heat conducting with a refractive index ranging from 1 to 2. The solution involves solving simultaneously the energy equation and the integral equation for the radiative flux gradient. The energy equation is solved using an implicit finite volume scheme and the integral equation of radiative heat transfer is solved using the singularity technique and Gaussian integration. The effects of scattering are investigated. It is shown that scattering has a significant effect on the transient temperature distribution and the transient mean temperature of the layer.     

Radiative transfer theory with time delay for the effect of pulse entrapping in a resonant random medium: General transfer equation and point-like scatterer model

Abstract

A pulse propagation of a vector electromagnetic wave field in a discrete random medium under the condition of Mie resonant scattering is considered on the basis of the Bethe–Salpeter equation in the two-frequency domain in the form of an exact kinetic equation which takes into account the energy accumulation inside scatterers. The kinetic equation is simplified using the transverse field and far wave zone approximations which give a new general tensor radiative transfer equation with strong time delay by resonant scattering. This new general radiative transfer equation, being specified in terms of the low-density limit and the resonant point-like scatterer model, takes the form of a new tensor radiative transfer equation with three Lorentzian time-delay kernels by resonant scattering. In contrast to the known phenomenological scalar Sobolev equation with one Lorentzian time-delay kernel, the derived radiative transfer equation does take into account effects of (i) the radiation polarization, (ii) the energy accumulation inside scatterers, (iii) the time delay in three terms, namely in terms with the Rayleigh phase tensor, the extinction coefficient and a coefficient of the energy accumulation inside scatterers, respectively (i.e. not only in a term with the Rayleigh phase tensor). It is worth noting that the derived radiative transfer equation is coordinated with Poynting's theorem for non-stationary radiation, unlike the Sobolev equation. The derived radiative transfer equation is applied to study the Compton–Milne effect of a pulse entrapping by its diffuse reflection from the semi-infinite random medium when the pulse, while propagating in the medium, spends most of its time inside scatterers. This specific albedo problem for the derived radiative transfer equation is resolved in scalar approximation using a version of the time-dependent invariance principle. In fact, the scattering function of the diffusely reflected pulse is expressed in terms of a generalized time-dependent Chandrasekhar H-function which satisfies a governing nonlinear integral equation. Simple analytic asymptotics are obtained for the scattering function of the front and the back parts of the diffusely reflected Dirac delta function incident pulse, depending on time, the angle of reflection, the mean free time, the microscopic time delay and a parameter of the energy accumulation inside scatterers. These asymptotics show quantitatively how the rate of increase of the front part and the rate of decrease of the rear part of the diffusely reflected pulse become slower with transition from the regime of conventional radiative transfer to that of pulse entrapping in the resonant random medium.     

A parametric study of radiative transfer with anisotropic scattering in a one-dimensional system

Radiative heat transfer in an absorbing, emitting, anistropically-scattering, one-dimensional medium is analyzed. Unlike many of the existing works, the present analysis does not require a known temperature distribution within the medium. Assuming a model of linear anistropic scattering, the transfer equation and the energy equation are solved simultaneously by utilizing a recently developed successive approximation technique. Closed-form approximate solutions and accurate higher-order results are both presented. Calculations show that the relative importance of the anistropic scattering effect generally decreases with decreasing wall emissivity and decreasing optical thickness. For radiative equilibrium without internal heat generation, it is demonstrated that the anistropic-scattering heat-transfer results can be approximated quite adequately by the isotropic-scattering result with the introduction of the concept of an effective optical thickness. For media with internal heat generation, an interesting effect of the scattering albedo is observed. It is established that, in the limit of a large scattering albedo, the temperature of the medium approaches a constant value that is independent of anistropic-scattering effects and wall emissivity. The exact limiting expressions for the temperature and apparent emissivity of an isothermal slab are found.     

表面不透明三层漫反射散射性介质耦合换热

本文采用射线踪迹、节点分析法研究了三层吸收、各向同性散射性介质层内的一维辐射和导热瞬态耦合换热,复合层表面不透明漫反射,介质层交界面半透明漫反射,且半透明漫反射交界面的反射率采用Fresnel反射定律确定。采用一层和二层辐射能量传递模型跟踪辐射能量在三层介质内的传递,从而推导出辐射传递系数。运用辐射传递系数求解辐射源项,在辐射对流边界条件下、采用全隐格式求解瞬态能量方程,并从机理上研究了辐射和导热耦合换热过程。     

Apparent optical properties of spherical particles in absorbing medium

An apparent absorption efficiency for spherical particles in absorbing medium is introduced to take into account the non-exponential absorption of the near-field scattered radiation in the absorbing medium. The apparent extinction, which is the summation of the apparent scattering efficiency following previous studies and the apparent absorption efficiency, is the same as the actual extinction. These apparent optical properties are suited to radiative transfer equations.     

Exploiting the linearity of radiative transfer

A standard problem in radiative transfer is finding the external and internal radiative fields produced by uniform, parallel rays illuminating the top of a one-dimensional, scattering and absorbing medium of finite optical thickness. This problem has been solved in several ways with various physical restrictions. One approach is by finding the source function that represents the rate of production of scattered radiation per unit volume per unit solid angle at each point in the medium. The present paper develops and uses the idea that the standard source function is an influence function for a given medium. The linearity of radiative transfer is then used to find certain general source functions in terms of the standard one. The usefulness of the above concept is demonstrated by the following four problems: (1) derivation of Chandrasekhar's four principles of invariance from the radiative transfer equation, (2) derivation of the equations governing Chandrasekhar's X- and Y- functions without using the invariance principles or resolvent kernels, (3) finding the source function for a medium with a Lambert's-law bottom, and (4) finding the source function for a medium with a bottom that is a perfect specular reflector.     

Comparison of radiative transfer approximations for a highly forward scattering planar medium

We examine critically the accuracy of the two-flux, spherical harmonics and discrete ordinates methods for predicting radiative transfer in a planar, highly-forward scattering and absorbing medium. Numerical results for the radiative fluxes show that the two-flux and P₃-approximations yield accurate results compared to solutions based on the F_N-method. Indeed, these approximate methods are relatively simple and have potential for generalization to predict radiative transfer in multidimensional systems, as long as an appropriate simplification of the phase function is utilized.