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.     

Continuum radiative heat transfer modeling in media consisting of optically distinct components in the limit of geometrical optics

Continuum-scale equations of radiative transfer and corresponding boundary conditions are derived for a general case of a multi-component medium consisting of arbitrary-type, non-isothermal and non-uniform components in the limit of geometrical optics. The link between the discrete and continuum scales is established by volume averaging of the discrete-scale equations of radiative transfer by applying the spatial averaging theorem. Precise definitions of the continuum-scale radiative properties are formulated while accounting for the radiative interactions between the components at their interfaces. Possible applications and simplifications of the presented general equations are discussed.     

A Monte Carlo method for 3D thermal infrared radiative transfer

A 3D Monte Carlo model for specific application to the broadband thermal radiative transfer has been developed in which the emissivities for gases and cloud particles are parameterized by using a single cubic element as the building block in 3D space. For spectral integration in the thermal infrared, the correlated k-distribution method has been used for the sorting of gaseous absorption lines in multiple-scattering atmospheres involving 3D clouds. To check the Monte-Carlo simulation, we compare a variety of 1D broadband atmospheric fluxes and heating rates to those computed from the conventional plane-parallel (PP) model and demonstrate excellent agreement between the two. Comparisons of the Monte Carlo results for broadband thermal cooling rates in 3D clouds to those computed from the delta-diffusion approximation for 3D radiative transfer and the independent pixel-by-pixel approximation are subsequently carried out to understand the relative merits of these approaches.     

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.     

An efficient diffusion approximation for 3D radiative transfer parameterization: application to cloudy atmospheres

The three-dimensional (3D) diffusion radiative transfer equation, which utilizes a four-term spherical harmonics expansion for the scattering phase function and intensity, has been efficiently solved by using the full multigrid numerical method. This approach can simulate the transfer of solar and thermal infrared radiation in inhomogeneous cloudy conditions with different boundary conditions and sharp boundary discontinuity. The correlated k-distribution method is used in this model for incorporation of the gaseous absorption in multiple-scattering atmospheres for the calculation of broadband fluxes and heating rates in the solar and infrared spectra. Comparison of the results computed from this approach with those computed from plane-parallel and 3D Monte Carlo models shows excellent agreement. This 3D radiative transfer approach is well suited for radiation parameterization involving 3D and inhomogeneous clouds in climate models.     

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.     

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 net-exchange Monte Carlo approach to radiation in optically thick systems

A Monte Carlo approach to radiative transfer in participating media is described and tested. It solves to a large extent the well known problem of Monte Carlo simulation of optically thick absorption configurations. The approach which is based on a net-exchange formulation and on adapted optical path sampling procedures is carefully designed to insure satisfactory convergence for all types of optical thicknesses. The need for such adapted algorithms is mainly related to the problem of gaseous line spectra representation in which extremely large ranges of optical thicknesses may be simultaneously encountered. The algorithm is tested against various band average computations for simple geometries using the Malkmus statistical narrow band model.     

The role of cloud-scale resolution on radiative properties of oceanic cumulus clouds

Both individual and combined effects of the horizontal and vertical variability of cumulus clouds on solar radiative transfer are investigated using a two-dimensional (x- and z-directions) cloud radar dataset. This high-resolution dataset of typical fair-weather marine cumulus is derived from ground-based cloud radar observations. The domain-averaged (along x-direction) radiative properties are computed by a Monte Carlo method. It is shown that (i) different cloud-scale resolutions can be used for accurate calculations of the mean absorption, upward and downward fluxes; (ii) the resolution effects can depend strongly on the solar zenith angle; and (iii) a few cloud statistics can be successfully applied for calculating the averaged radiative properties.     

Application of the modified differential approximation for radiative transfer to arbitrary geometry

The first-order spherical harmonics method (or P₁ approximation) has found prolific usage for approximate solution of the radiative transfer equation (RTE) in participating media. However, the accuracy of the P₁ approximation deteriorates as the optical thickness of the medium is decreased. The modified differential approximation (MDA) was originally proposed to remove the shortcomings of the P₁ approximation in optically thin situations. This article presents algorithms to apply the MDA to arbitrary geometry—in particular, geometry with obstructions, and inhomogeneous media. The wall-emitted component of the intensity was computed using a combined view-factor and ray-tracing approach. The Helmholtz equation, arising out of the medium-emitted component, was solved using an unstructured finite-volume procedure. The general procedure was validated for both two-dimensional (2D) and three-dimensional (3D) geometries against benchmark Monte Carlo results. The accuracy of MDA was found to be superior to the P₁ approximation for all optical thicknesses. Its accuracy, when compared with the discrete ordinates method (both S₆ and S₈), was found to be clearly superior in optically thin situations, but problem dependent in optically intermediate and thick situations. For 3D geometries, calculation and storage of the view-factor matrix was found to be a major shortcoming of the MDA. In addition, for inhomogeneous media, calculation of optical distances requires a ray-tracing procedure, which was found to be a bottleneck from a computational efficiency standpoint. Several strategies to reduce both memory and computational time are discussed and demonstrated.     

The Monte Carlo atmospheric radiative transfer model McArtim: Introduction and validation of Jacobians and 3D features

A new Monte Carlo atmospheric radiative transfer model is presented which is designed to support the interpretation of UV/vis/near-IR spectroscopic measurements of scattered Sun light in the atmosphere. The integro differential equation describing the underlying transport process and its formal solution are discussed. A stochastic approach to solve the differential equation, the Monte Carlo method, is deduced and its application to the formal solution is demonstrated. It is shown how model photon trajectories of the resulting ray tracing algorithm are used to estimate functionals of the radiation field such as radiances, actinic fluxes and light path integrals. In addition, Jacobians of the former quantities with respect to optical parameters of the atmosphere are analyzed. Model output quantities are validated against measurements, by self-consistency tests and through inter comparisons with other radiative transfer models.     

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.     

高倍聚光太阳模拟器的设计

研究了一种高倍聚光、高均匀性太阳模拟器,它主要由短弧氙灯、椭球聚光镜、高速快门和积分器构成,辐射功率超过5 kW,辐射峰值超过1 800个太阳常数（标准AM1.5∶1太阳常数=1 kW/m2）。介绍了该模拟器的光学设计,主要工程特点和仿真结果。应用蒙特卡洛光线追迹技术优化几何配置,以使光源到目标面的辐射能量转换效率最大化。仿真结果显示：在直径20 mm的圆形目标区域其平均辐射超过1 000 kW/m2,相应的辐照不均匀度小于6%。利用该装置模拟了高倍聚光太阳系统的辐射特性,为研究高效率太阳能电池的热化学过程和测试先进的高温材料提供了实验平台。     

An alternative Monte Carlo approach to the thermal radiative transfer problem

The usual Monte Carlo approach to the thermal radiative transfer problem is to view Monte Carlo as a solution technique for the nonlinear thermal radiative transfer equations. The equations contain time derivatives which are approximated by introducing small time steps. An alternative approach avoids time steps by using Monte Carlo to directly sample the time at which the next event occurs. That is, the time is advanced on a natural event-by-event basis rather than by introducing an artificial time step.     

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.     

基于蒙特卡罗方法的卫星光学波段闪电辐射观测模拟

我国下一代静止气象卫星风云四号将首次搭载光学闪电成像仪上天,对闪电和与之相联系的强对流天气进行实时、连续观测。对未来风云四号卫星闪电观测资料的特性进行预先研究,一个重要的途径就是建立光学波段的闪电辐射传输模拟算法,对闪电辐射的物理过程及卫星观测到的闪电辐射光学特征进行模拟研究。结合风云四号卫星闪电成像仪观测的几何和波长参数,在构建简化的闪电光源模型和雷暴云模型的基础上,利用蒙特卡罗随机模拟和光子追踪方法,模拟得到了风云四号卫星闪电成像仪观测到的云顶闪电信号,并对其辐射特征进行了初步研究。结果表明,基于蒙特卡罗方法的闪电辐射传输模拟,可以从理论上揭示卫星光学波段闪电辐射观测与最重要影响参数之间的联系,这将为未来星载闪电资料的应用提供非常有价值的信息。     

A diffusion-based approximate model for radiation heat transfer in a solar thermochemical reactor

An approximate numerical method for fast calculations of the radiation heat transfer in a solar thermochemical reactor cavity is formulated based on the separate treatment of the solar and thermal radiative exchange by the diffusion approach. The usual P₁ approximation is generalized by applying an equivalent radiation diffusion coefficient for the optically thin central part of the cavity. The resulting boundary-value problems are solved using the finite element algorithm. The accuracy of the model is assessed by comparing the results to those obtained by a pathlength-based Monte Carlo simulation. The applicability of the proposed model is demonstrated by performing calculations for an example problem, which incorporates a range of parameters typical for a solar chemical reactor and the spectral radiative properties of polydisperse zinc oxide particles.     

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.     

Comparison of ray-tracing methods for semitransparent slab with variable spatial refractive index

The ray-tracing technique has the main difficulty in solving radiative transfer in the medium with variable spatial refractive index. Recently, three methods have been developed for the application of the ray-tracing technique in those medium. To compare and discuss the numerical characteristics of those methods, a semitransparent slab with variable spatial refractive index is taken as an example, and the reflectivity and the transmissivity of the slab are computed by the curved ray-tracing method, the multi-layer approach, and the discrete curved ray-tracing method, respectively. As the result, it is shown that, the discrete curved ray-tracing method gives the result with good accuracy and convergence characteristics than the multi-layer approach. Due to accounting physically inexistent reflection on the interface between sublayers, the multi-layer approach converges slowly.     

An adaptive emission model for Monte Carlo simulations in highly inhomogeneous media represented by stochastic particle fields

Traditional Monte Carlo ray-tracing (MCRT) methods for continuous participating media are not applicable in media represented by point masses (or stochastic particles) frequently encountered in combustion modeling. In the authors’ previous work several ray models and particle models have been proposed for radiation simulations in such media. In the present paper an efficient emission scheme is developed for MCRT in highly inhomogeneous media represented by particle fields. Ray energies are limited to a narrow range to reduce statistical error, by having particles emit numbers of photons proportional to their emissive power (including combination of weak particles). A method to evaluate the radiative heat source, required by the overall energy equation, is also developed. A particle field representing the highly inhomogeneous medium in a turbulent jet flame is employed to test the proposed methods.