任意入射辐射条件下介质特性对瞬态辐射传递的影响

用基于Monte Carlo法的DRESOR法在平行平板系统内具有吸收、无发射介质中研究不同波形入射、壁面反射、介质散射率、光学厚度、各向异性散射等条件对瞬态辐射传递的影响.任意连续波形入射辐射是目前大多数数值方法很难处理的瞬态辐射问题,而DRESOR法通过在系统内计算一单位入射辐射能对介质的DRESOR数分布,就能计算任意连续波形入射辐射条件下高方向分辨率的瞬态辐射强度结果.DRESOR法和Monte Carlo法计算的结果进行了比较验证,两者吻合较好,证明了DRESOR法处理瞬态入射辐射问题的正确性和有效性.     

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法处理瞬态入射辐射问题的能力．     

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.     

The Iterative-DRESOR method to solve radiative transfer in a plane-parallel, anisotropic scattering medium with specular-diffuse boundaries

Using the intensity with high directional resolution obtained by the Basic-DRESOR method as an initial guess, which is substituted into the integrated radiative transfer equation (IRTE), an iterative algorithm is proposed, called the Iterative-DRESOR method. This method can reduce the error levels of the intensity from several percent using the Basic-DRESOR method to a level of less than 1.0×10⁻⁶ with acceptable computation costs. The method is also validated against the exact heat flux in literature in some cases. It further clarifies some uncertain results for the reflectance in a pure, linearly anisotropic scattering medium with specular-diffuse boundaries. The directional distributions of intensity are obviously influenced by the reflecting modes of the boundary, especially in the zone near the boundary. The reflecting mode of an emitting boundary has little effect on the transmittance or reflectance. The reflecting mode of a non-emitting boundary also has little effect on the transmittance, but it obviously influences the reflectance. The difference between the reflectance for specular and diffuse boundaries increases at first, and then decreases, as the optical thickness of the medium increases. The difference will decrease as the scattering albedo of the medium increases, and it is negligible when the medium is pure scattering. The effect of the scattering phase function of the medium on the difference can also not be ignored. The Iterative-DRESOR method is expected to strengthen the capability of the Monte Carlo method to produce accurate results and to validate the results of other methods to solve RTE.     

DRESOR法对平行入射辐射问题的研究

本采用一种基于蒙特卡洛法(Monte Carlo Method,MCM)求解辐射传递方程(Radiative Transfer Equation, RTE)的快捷、有效的方法-DRESOR法(Distributions of Ratios of Energy Scattered Or Reflected)在一维充满吸收、各向同性散射介质平行平板中,外部有平行入射条件下,求解计算空间点的辐射强度沿空间方向角的分布,而不需要辐射平衡和在空间位置坐标和方向角度坐标上同时离散辐射传递方程进行迭代求解。     

Radiative heat transfer in a planar medium with anisotropic scattering and directional boundaries

Radiation heat transfer in an absorbing, emitting and scattering medium has been the subject of many previous investigations. Most solutions are numerically complex and the existing analytical solutions are restricted in application by the simplifying assumptions involved. A plane-parallel medium is considered which scatters anisotropically. The boundaries are considered to be specular reflectors, as predicted by Fresnel's relations, while the diffusely incident radiation is refracted according to Snell's law. The emission is restricted to a medium with a uniform temperature distribution. Approximate closed-form solutions for the radiative heat flux and incident intensity are presented for dielectric layers and linear anisotropic scattering. Numerical results are also presented and show that the effects of directional boundaries, anisotropic scattering, scattering albedo and optical depth are accurately predicted by the approximate solution.     

Analysis of collimated radiation in participating media using the discrete transfer method

A general formulation of the discrete transfer method is provided to analyze radiative heat transfer problems in a participating medium subjected to collimated radiation. The formulation is validated by considering 1-D planar absorbing, emitting and anisotropically scattering gray medium in radiative equilibrium. Anisotropy of the medium is approximated by linear anisotropic phase function. For the purpose of comparison, the problem is also solved analytically. Results are obtained for different angles of incidence of the collimated radiation. At a given angle of incidence, results are obtained for forward, isotropic and backward scattering situations. Heat flux results are compared over a wide range of values of the extinction coefficient. Emissive power distributions in the medium are also obtained for some cases. The discrete transfer method results are found to compare very well with the analytic results.     

半透明平板边界入射辐射热流密度的反问题

对一维半透明平板内辐射、导热及边界对流耦合换热过程进行了研究。提出了一种由一侧边界出射辐射强度反演另一侧边界入射辐射热流密度的方法。通过对各向异性散射、吸收系数、散射系数、边界外侧来流温度、对流换热系数、半透明平板的导热系数和平板厚度等参数对反演精度影响的分析表明,方法是可行的。     

Time-dependent radiative transfer using Chapman-Enskog maximum entropy method

The time-dependent problems of radiative transfer involve a coupling between radiation and material energy fields and are nonlinear because of proposed temperature dependence of the medium characteristics in semi-infinite medium with Rayleigh anisotropic scattering. By means of the limited flux, Chapman-Enskog and maximum entropy technique the time-dependent radiative transfer equation has been solved explicitly. The maximum entropy method is used to solve the resulting differential equation for radiative energy density. The calculations are carried out for temperature (normalized dimensionless) Θ(x,τ), radiative energy density and net flux with Rayleigh and anisotropic scattering for different space at different times.     

Direct numerical simulation of turbulence/radiation interaction in premixed combustion systems

An important fundamental issue in chemically reacting turbulent flows is turbulence/radiation interaction (TRI); TRI arises from highly nonlinear coupling between temperature and composition fluctuations. Here, a photon Monte Carlo method for the solution of the radiative transfer equation has been integrated into a turbulent combustion direct numerical simulation (DNS) code. DNS has been used to investigate TRI in a canonical configuration with systematic variations in optical thickness. The formulation allows for nongray gas properties, scattering, and general boundary treatments, although in this study, attention has been limited to gray radiation properties, no scattering, and black boundaries. Individual contributions to emission and absorption TRI have been isolated and quantified. Of particular interest are intermediate values of optical thickness where, for example, the smallest hydrodynamic and chemical scales are optically thin while the largest turbulence scales approach an optically thick behavior. In the configuration investigated, the temperature self-correlation contribution (emission) is primarily a function of the ratio of burned-gas temperature to unburned-gas temperature, and is the dominant contribution to TRI only in the optically thin limit. Even in the most optically thin case considered, the absorption coefficient–Planck function correlation and absorption coefficient–intensity correlation are not negligible. At intermediate values of optical thickness, contributions from all three correlations are significant.     

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

Graphical and tabular results are presented for the back-scattered intensity from a finite two-dimensional cylindrical medium exposed to a Gaussian beam of radiation. Also, results for the source function and flux at the boundaries are presented. The influence of optical thickness and albedo are most pronounced at large optical radii. The semi-infinite results can be used to approximate the finite case for small optical radii. Ranges for single, double, and multiple scattering are discussed. For locations far from the incident beam, the results can be expressed in terms of universal functions independent of beam size. A method is presented for extending the isotropic results to the anisotropic case where the phase function is made up of a spike superimposed on an otherwise isotropic phase function.     

Radiative transfer in a two-dimensional rectangular medium exposed to diffuse radiation

The integral equation for radiative transfer in a two-dimensional rectangular scattering medium exposed to diffuse radiation is solved numerically by removing the singularity. This method yielded accurate results except at very large optical thicknesses. Graphical and tabular results for the source function, flux, and intensity are presented. The source function is also calculated using the first term of a Taylor series expansion. The Taylor series is fairly accurate for small optical thicknesses and columnar geometries. A method is presented for extending these results to the problem of a strongly anisotropic scattering phase function which is made up of a spike in the forward direction superimposed on an isotropic phase function.     

吸收-透射界面下正弦折射率介质内辐射换热

对具有吸收-透射性边界面的梯度折射率半透明介质层,建立了介质内热辐射传递与边界面辐射换热的数理模型,并采用数值弯曲光线跟踪法求解介质内的热辐射传递。通过数值模拟,分析了正弦折射率下,边界面的反射特性、吸收率以及介质层光学厚度对介质内热辐射平衡温度场及热流分布的影响。结果表明,边界面的反射特性与吸收率对介质内辐射换热均有重要影响,吸收率的影响与边界面反射特性、介质层光学厚度及环境条件相关,呈现特征不同的作用。     

球谐函数法求解辐射传输方程的假散射和射线效应

假散射和射线效应是辐射传输方程近似解法中出现的特有误差.从辐射传输方程的近似求解过程出发,在定性分析的基础上构造物理模型,通过数值模拟研究球谐函数法(P₁和P₃近似法)的假散射和射线效应.构造激光平行入射和倾斜入射二维半透明介质的物理模型,通过内部温度场的分布特征研究假散射.构造顶部侧面保持高温而其余侧面保持低温的二维半透明介质方案,通过对比底面边界净热流密度分析射线效应.计算结果表明球谐函数法中同时存在假散射和射线效应,P₃近似比P₁近似减小了射线效应.同时,球谐函数法的射线效应随光学厚度的增加而减小.     

Transient radiative transfer in participating media with pulse-laser irradiation—an approximate Galerkin solution

An assessment is made of the Galerkin technique as an effective method of solution for transient radiative transfer problems in participating media. A one-dimensional absorbing and isotropically scattering plane-parallel gray medium irradiated with a short-pulse laser on one of its boundaries is considered for the application of the method. The medium is non-emitting and the boundaries are non-reflecting and non-refracting. In the integral formulation of the problem for the source function, the time-wise variation of the radiation intensity at any point and in any direction in the medium is assumed to be the same as the time-wise variation of the average intensity at the same point as an approximation for the application of the method. The transient transmittance and reflectance of the medium are evaluated for various values of the optical thickness, scattering albedo and pulse duration. The results are in agreement with those available in the literature. It is demonstrated that the method is relatively simple to implement and yields accurate results.     

散射相函数对一维介质内辐射传递的影响规律

采用有限体积法研究了一维线性各向异性散射介质内散射相函数对辐射传递的影响规律.经与理论解、辐射元法、蒙特卡洛法计算结果比较表明,有限体积法的计算结果更可靠,且不同散射相函数的辐射换热系统中,其无因次热流之比与光学厚度之间存在某种单调变化的函数关系,利用该函数关系可以检验模型的准确度.     

Transient,combined conduction and radiation in an absorbing,emitting, and isotropically scattering solid sphere

Transient, combined conduction and radiation is solved in an absorbing, emitting, and isotropically scattering solid sphere with a black boundary initially at a uniform temperature and for times t > 0 subjected to a constant temperature at the spherical surface. The collocation method is used to solve the radiation part of this problem and the implicit finite difference scheme is used to solve the conduction part. The effects of the conduction-to-radiation parameter, the single scattering albedo, the optical thickness of the medium on the temperature distribution and the heat flux in the medium are examined.     

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.     

Temperature response in participating media with anisotropic scattering caused by pulsed lasers

It is not by isotropic scattering but by anisotropic scattering that radiant energy is redistributed in some materials containing real particles, fibers, or impurities. In some instances, great difference can be caused in transient thermal behavior between isotropic scattering and anisotropic scattering media. Ray tracing method combined with Hottel's zonal method is introduced to deduce thermal radiative source term for various optical boundary conditions induced by collimated incidence passing through translucent boundary. Temperature response caused by laser pulse at non-incident side of participating and anisotropic scattering media is examined. We investigate effects of scattering albedo, scattering phase function, initial temperature of media and thickness of media on temperature response. Results obtained for anisotropic scattering media are compared with those for isotropic scattering one and show that anisotropic scattering must be considered in the simulating measurement of thermophysical properties by the laser flash method for some materials with big scattering albedo which behave anisotropically, or big error will be introduced; forward scattering can increase excess temperature and backward scattering can decrease it at non-incident side of the considered sample irradiated by laser pulse.     

Radiation transfer in an isotropically scattering homogeneous solid sphere

Radiative heat transfer in an absorbing, emitting, isotropically scattering sphere with a diffusively reflecting and emitting boundary is solved by the Galerkin method. By using the expressions given in this work, the radiation heat flux, the angular distribution of radiation intensity and the divergence of the radiation heat flux anywhere in the sphere can be determined with a high degree of accuracy for all values of the single scattering albedo, from small to moderately large optical radii. Results are presented for representative cases to illustrate the accuracy and the convergence of the method.