球形微粒子近场辐射换热研究

从偶极辐射理论出发,结合普朗克黑体辐射分布,推导出近场区热辐射分布。采用MC法模拟了两个纳米球形 村子的辐射换热,考察相天参数对近场辐射换热的影响。表明其辐射热流可比黑体辐射高8个数量级。其热流随间距的增 大急剧减小。小研究对于纳米粒子群的传热研究有积极意义。     

Monte Carlo ray-tracing simulation for radiative heat transfer in turbulent fluctuating media under the optically thin fluctuation approximation

The Monte Carlo ray-tracing method (MCRT) based on the concept of radiation distribution factor is extended to solve radiative heat transfer problem in turbulent fluctuating media under the optically thin fluctuation approximation. A one-dimensional non-scattering turbulent fluctuating media is considered, in which the mean temperature and absorption coefficient distribution are assumed and the shape of probability density function is given. The distribution of the time-averaged volume radiation heat source is solved by MCRT and direct integration method. It is shown that the results of MCRT based on the concept of radiation distribution factor agree with these of integration solution very well, but results of MCRT based on the concept of radiative transfer coefficient do not agree with these of integration solution. The solution of time-averaged radiative transfer equation by the concept of radiative transfer coefficient should be treated with caution.     

矩形介质内辐射换热的有限元法

利用有限单元法离散求解辐射传递方程和能量控制方程.分别计算了边界为黑体和灰体条件下矩形吸收、发射、各向同性散射介质内的平均入射强度和温度分布,并同蒙特卡罗法(M-C法)计算结果进行了比较.     

管内耦合换热的辐射热流特征及抽样模式比较

通过数值模拟,研究了圆管层流入口段辐射与对流祸合换热中的辐射热流分布特征,提出并比较了采用蒙特卡罗法计算辐射换热的两种抽样模式。结果表明,在圆管入口段的耦合换热中,同一截面的径向辐射热流密度与轴向辐射热流密度数量级相当,世分布形态差别很大。用蒙特卡罗方法计算辐射换热时,采取等权能束的抽样模式,可消除因抽样能束的能量不同而造成的计算误差、降低计算量、提高计算精度。     

Comparison of the standard weighted-sum-of-gray-gases with the absorption-line blackbody distribution function for the computation of radiative heat transfer in H2O/CO2 mixtures

This paper presents the computation of radiation heat transfer in a cylindrical enclosure in which the dimensions, the chemical species concentrations and the temperature fields make a realistic representation of an actual combustion chamber. Two gas models are applied and compared: the absorption-line blackbody distribution function (ALBDF), and the standard weighted-sum-of-gray-gases (WSGG) based on coefficients and correlations that are widely used in engineering. While the standard WSGG is restricted to the assumption of homogeneous gas mixture, the ALBDF can be applied to both homogeneous and non-homogeneous media. For the two gas models, the radiative exchanges are computed with the aid of the Monte Carlo method. The results show considerable discrepancies between the WSGG and the ALBDF models for the homogeneous medium. In addition, the importance of considering the non-homogeneity of the medium for an accurate computation of the radiative heat transfer is shown.     

Radiation symmetry test and uncertainty analysis of Monte Carlo method based on radiative exchange factor

Using Monte Carlo method, the paper investigates the radiative heat transfer in participating media. Based on the radiative exchange factor, an uncertainty analysis of Monte Carlo method is undertaken and the corresponding mathematical expressions are deduced to predict its accuracy. Furthermore, randomness properties of pseudorandom number generators are investigated, and a model to test radiation symmetry is adopted to validate the performance of some generators. The paper studies the effects of energy bundle numbers, discretization schemes, emission location, optical thicknesses, wall emissivity and CPU time on the numerical accuracy. In addition, the simulation results are proved to give a reference for using Monte Carlo method, which is applicable for calculation of the radiative exchange factor.     

太阳能吸热器玻璃窗辐射与导热计算

在太阳能塔式热发电空气吸热器工作过程中,吸热器玻璃窗受到聚焦太阳光、环境以及吸热面（体）的复合辐射。本文采用改进蒙特卡罗法、有限容积法,研究了太阳能空气吸热器石英玻璃窗半透明介质内部辐射与导热的耦合换热,研究了不同厚度半透明介质石英玻璃窗在工作面处于不同第三类边界条件下的稳态温度分布。研究表明石英玻璃的厚度和内侧空气温...     

Monte Carlo discrete curved ray-tracing method for radiative transfer in an absorbing-emitting semitransparent slab with variable spatial refractive index

A Monte Carlo discrete curved ray-tracing method is developed to analyze the radiative transfer in one-dimensional absorbing-emitting semitransparent slab with variable spatial refractive index, in which the Monte Carlo method is combined with the discrete curved ray-tracing method. A problem of radiative equilibrium with linear variable spatial refractive index is taken as an example to examine the accuracy of the proposed method. The temperature distributions and the dimensionless radiative heat flux are determined by the proposed method and compared with the data in references, which are obtained by other different methods. The results show that the Monte Carlo discrete curved ray-tracing method has a good accuracy in solving the radiative transfer in one-dimensional semitransparent slab with variable spatial refractive index.     

Discontinuous finite element method for radiative heat transfer in semitransparent graded index medium

To avoid the complicated and time-consuming computation of curved ray trajectories, a discontinuous finite element method based on discrete ordinate equation is extended to solve the radiative transfer problem in a multi-dimensional semitransparent graded index medium. Two cases of radiative heat transfer in two-dimensional rectangular gray semitransparent graded index medium enclosed by opaque boundary are examined to verify this discontinuous finite element method. Special layered and radial graded index distributions are considered. The predicted dimensionless net radiative heat fluxes and dimensionless temperature distributions are determined by the discontinuous finite element method and compared with the results obtained by the curved Monte Carlo method in references. The results show that the discontinuous finite element method has a good accuracy in solving the multi-dimensional radiative transfer problem in a semitransparent graded index medium.     

Infrared radiative transfer modelling in a 3D scattering cloudy atmosphere: Application to limb sounding measurements of cirrus

The Monte Carlo cloud scattering forward model (McClouds_FM) has been developed to simulate limb radiative transfer in the presence of cirrus clouds, for the purposes of simulating cloud contaminated measurements made by an infrared limb sounding instrument, e.g. the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS). A reverse method three-dimensional Monte Carlo transfer model is combined with a line-by-line model for radiative transfer through the non-cloudy atmosphere to explicitly account for the effects of multiple scattering by the clouds. The ice cloud microphysics are characterised by a size distribution of randomly oriented ice crystals, with the single scattering properties of the distribution determined by accurate calculations accounting for non-spherical habit.A comparison of McClouds_FM simulations and real MIPAS spectra of cirrus shows good agreement. Of particular interest are several noticeable spectral features (i.e. H₂O absorption lines) in the data that are replicated in the simulations: these can only be explained by upwelling tropospheric radiation scattered into the line-of-sight by the cloud ice particles.     

气动加热下高温陶瓷材料内的瞬态耦合换热

采用控制容积法结合蒙特卡罗法与谱带模型,数值模拟了气动加热下高温陶瓷内的辐射与导热瞬态耦合换热.分析了冷却边界条件、各向同性散射和表面特性对高温陶瓷内瞬态温度分布的影响.结果表明,高温下陶瓷内的辐射换热作用十分明显,陶瓷材料的表面特性对瞬态温度分布影响很大,各向同性散射也有明显影响.     

Discrete vs. continuum-scale simulation of radiative transfer in semitransparent two-phase media

The mathematical formulation of the continuum approach to radiative transfer modeling in two-phase semi-transparent media is numerically validated by comparing radiative fluxes computed by (i) direct, discrete-scale and (ii) continuum-scale approaches. The analysis is based on geometrical optics. The discrete-scale approach uses the Monte Carlo ray-tracing applied directly to real 3D geometry measured by computed tomography. The continuum-scale approach is based on a set of continuum-scale radiative transfer equations and associated radiative properties, and employs the Monte Carlo ray-tracing for computations of radiative fluxes and for computations of the radiative properties. The model two-phase media are reticulate porous ceramics and a particle packed bed, each composed of semitransparent solid and fluid phases. The results obtained by the two approaches are in good agreement within the limits of statistical uncertainty. The continuum-scale approach leads to a reduction in computational time by approximately one order of magnitude, and is therefore suited to treat radiative transfer problems in two-phase media in a wide range of engineering applications.     

Application of the spatial averaging theorem to radiative heat transfer in two-phase media

The spatial averaging theorem is applied to rigorously derive continuum-scale equations of radiative transfer in two-phase media consisting of arbitrary-type phases in the limit of geometrical optics. The derivations are based on the equations of radiative transfer and the corresponding boundary conditions applied at the discrete-scale to each phase, and on the discrete-scale radiative properties of each phase and the interface between the phases. The derivations confirm that radiative transfer in two-phase media consisting of arbitrary-type phases in the range of geometrical optics can be modeled by a set of two continuum-scale equations of radiative transfer describing the variation of the average intensities associated with each phase. Finally, a Monte Carlo based methodology for the determination of average radiative properties is discussed in the light of previous pertinent studies.     

Microwave radiative transfer intercomparison study for 3-D dichroic media

Three different numerical methods capable of solving the radiative transfer of microwave radiation within 3-D dichroic media are compared. A case study, represented by an intense rain shaft populated by perfectly oriented oblate raindrops, is analysed in detail, including a discussion of the behaviour of all four Stokes components.Results demonstrate an acceptable agreement between all Monte Carlo methods. The method based on a discrete ordinates scheme agrees only qualitatively with the Monte Carlo outputs. Because of its lower computational cost the backward Monte Carlo technique based on importance sampling represents the most efficient way to face passive microwave radiative transfer problems related to optically thick 3-D structured clouds including non-spherical preferentially oriented hydrometeors.     

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

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

Benchmark numerical solutions for radiative heat transfer in two-dimensional medium with graded index distribution

In graded index media, the ray goes along a curved path determined by Fermat principle. Generally, the curved ray trajectory in graded index media is a complex implicit function, and the curved ray tracing is very difficult and complex. Only for some special refractive index distributions, the curved ray trajectory can be expressed as a simple explicit function. Two important examples are the layered and the radial graded index distributions. In this paper, the radiative heat transfer problems in two-dimensional square semitransparent with layered and radial graded index distributions are analyzed. After deduction of the ray trajectory, the radiative heat transfer problems are solved by using the Monte Carlo curved ray-tracing method. Some numerical solutions of dimensionless net radiative heat flux and medium temperature are tabulated as the benchmark solutions for the future development of approximation techniques for multi-dimensional radiative heat transfer in graded index media.     

A 3-D radiation model for non-grey gases

A three-dimensional radiation code based on method of lines (MOL) solution of discrete ordinates method (DOM) coupled with spectral line-based weighted sum of grey gases (SLW) model for radiative heat transfer in non-grey absorbing-emitting media for use in conjunction with a computational fluid dynamics (CFD) code based on the same approach was developed. The code was applied to three test problems: two containing isothermal homogenous/non-homogenous water vapor and one non-isothermal water vapor/carbon dioxide mixture. Predictive accuracy of the code was evaluated by benchmarking its steady-state predictions against accurate results, calculated by ray tracing method with statistical narrow band model, available in the literature. Comparative testing with solutions of other methods is also provided. Comparisons reveal that MOL solution of DOM with SLW model provides accurate solutions for radiative heat fluxes and source terms and can be used with confidence in conjunction with CFD codes based on MOL.     

Discontinuous spectral element method for solving radiative heat transfer in multidimensional semitransparent media

A discontinuous spectral element method (DSEM) is presented to solve radiative heat transfer in multidimensional semitransparent media. This method is based on the general discontinuous Galerkin formulation. Chebyshev polynomial is used to build basis function on each element and both structured and unstructured elements are considered. The DSEM has properties such as hp-convergence, local conservation and its solutions are allowed to be discontinuous across interelement boundaries. The influences of different schemes for treatment of the interelement numerical flux on the performance of the DSEM are compared. The p-convergence characteristics of the DSEM are studied. Four various test problems are taken as examples to verify the performance of the DSEM, especially the performance to solve the problems with discontinuity in the angular distribution of radiative intensity. The predicted results by the DSEM agree well with the benchmark solutions. Numerical results show that the p-convergence rate of the DSEM follows exponential law, and the DSEM is stable, accurate and effective to solve multidimensional radiative transfer in semitransparent media.     

混合模拟方法分析三维封闭腔内辐射换热

１引言工程上对炉内三维辐射换热的计算已越来越重视,较为精确的计算方法是Ｈｏｔｔｅｌ的区域法［１,２］和ＭｏｎｔｅＣａｒｌｏ方法［３］,但两者都有一定的局限性而得不到广泛应用,原因如下：对区域法,由于计算区域间直接辐射交换面积时需求解多重积分,在计算时...     

Approximate analytical solution to normal emittance of semi-transparent layer of an absorbing, scattering, and refracting medium

A two-step approximate analytical solution for the normal emittance of a plane layer of an absorbing, scattering and refracting medium is derived analytically. The analysis is based on the transport approximation and the two-step solution method for radiative transfer. The high accuracy of the approximate solution, examined by comparing its results to those obtained independently by the discrete ordinates and Monte Carlo methods, makes it suitable for application in combined experimental-analytical studies to identify selected spectral radiative properties of dispersed media in the range of semi-transparency.