On the use of a meshless method for solving radiative transfer with the discrete ordinates formulations

A meshless method is presented for solving the radiative transfer equation in the discrete ordinates approach. It is shown that the primitive variables formulation is unstable for low values of the absorption coefficient while the even parity formulation is always stable and accurate.     

Radiative transfer with internal reflection via the converged discrete ordinates method

In response to the challenge of establishing highly accurate solutions to the plane layer radiative transfer equation with the simplest of methods, the converged discrete ordinates method is presented. With this algorithm of only finite difference, quadrature and acceleration, we show how to obtain highly accurate intensities for radiative transfer in a finite layer with internal surface reflection. The method features angular smoothing and angular interpolation through “faux” quadrature. In addition, a manufactured solution demonstrates the high accuracy of the method for forward peaked scattering. We consider scattering in a heterogeneous medium as a final demonstration.     

The method of lines solution of discrete ordinates method for radiative heat transfer in cylindrical enclosures

A radiation code based on method of lines solution of discrete ordinates method for radiative heat transfer in axisymmetric cylindrical enclosures containing absorbing-emitting medium was developed and tested for predictive accuracy by applying it to (i) test problems with black and grey walls (ii) a gas turbine combustor simulator enclosing a non-homogeneous absorbing-emitting medium and benchmarking its steady-state predictions against exact solutions and measurements. Comparisons show that it provides accurate solutions for radiative heat fluxes and can be used with confidence in conjunction with CFD codes based on the same approach.     

Evaluation of quadrature schemes for the discrete ordinates method

After reviewing the presently available quadrature schemes for the discrete ordinates method, the accuracy of different schemes is analyzed and evaluated. It is shown from a comprehensive error analysis that the moment conditions have to satisfied not only for the principal coordinates directions, as it is mostly carried out, but for any arbitrary test direction. Among the schemes with approximately 50 discrete ordinates the DCT-020-2468 quadrature was found to give the best accuracy. The highest accuracy among all schemes is achieved by the LC-11 quadrature which requires 96 discrete ordinates. This scheme is rarely used up to date and deserves more attention for high accuracy predictions.     

Finite element formulation of the first- and second-order discrete ordinates equations for radiative heat transfer calculation in three-dimensional participating media

Two finite element methods (FEMs), FEDOM1 and FEDOM2 (standing for the first and the second finite element discrete ordinates methods, respectively), are formulated and numerically tested. The reference second-order discrete equation is modified in its scattering terms and is applied to the problems of absorbing/emitting and anisotropically scattering media by using the FEM. Numerical features of the developed FEMs are compared with one of the discrete ordinates interpolation method (DOIM), which uses a finite difference scheme. Prediction results of radiative heat transfer by these two FEMs are compared with reference solutions and verified in three-dimensional enclosures containing participating media. The results of FEDOM1 and FEDOM2 agree well with exact solutions for the problem of absorbing/emitting medium with various range of optical thickness. Generally, the two FEMs show more accurate results than DOIM. And FEDOM1 shows more accurate results than FEDOM2 in most of the test problems. Both of the developed FEMs show reasonable results compared with published Monte Carlo solutions for the tested absorbing/emitting and anisotropically scattering media. Although the FEDOM2 is not as accurate as the FEDOM1, it shows its own advantages that it reduces CPU time and memory space of dependent variable to half.     

The method of lines solution of discrete ordinates method for non-grey media

A radiation code based on method of lines (MOL) solution of discrete ordinates method (DOM) for radiative heat transfer in non-grey absorbing-emitting media was developed by incorporation of a gas spectral radiative property model, namely wide band correlated-k (WBCK) model, which is compatible with MOL solution of DOM. Predictive accuracy of the code was evaluated by applying it to 1-D parallel plate and 2-D axisymmetric cylindrical enclosure problems containing absorbing-emitting medium and benchmarking its predictions against line-by-line solutions available in the literature. Comparisons reveal that MOL solution of DOM with WBCK model produces accurate results for radiative heat fluxes and source terms and can be used with confidence in conjunction with computational fluid dynamics codes based on the same approach.     

Application of multidimensional scheme and the discrete ordinate method to radiative heat transfer in a two-dimensional enclosure with diffusely emitting and reflecting boundary walls

Radiative heat transfer is the dominant mode of heat transfer in many engineering problems, including combustion chambers, space, greenhouses, rocket plume sensing, among others. The aim of this study is to develop an efficient method capable of eliminating ray effects in complex 2D situations and to use the developed code for other problems including combined conduction and convection in connection with CFD codes. A complete genuinely multidimensional discretization in two-dimensional discrete ordinates method is formulated to solve radiative heat transfer in a rectangular enclosure composed of diffusely emitting and reflecting boundaries and containing homogeneous media that absorbs, emits and scatters radiation. A new genuinely multidimensional differencing scheme is used to solve the radiative transfer equation with S₄, S₆, S₈, T₆, T₇, T₈ and T₉ angular quadrature schemes. Different cases are analyzed and the results are compared when possible with those obtained by others researchers.     

Theoretical analysis of porous radiant burners under 2-D radiation field using discrete ordinates method

This paper describes a theoretical study to investigate the heat transfer characteristics of porous radiant burners (PRBs). In the present work, a 2-D rectangular model is used to solve the governing equations for porous medium and gas flow before the premixed flame to the exhaust gas. The gas and the solid phases are considered in non-local thermal equilibrium and combustion in the porous medium is modeled by considering a non-uniform heat generation zone. The homogeneous porous media, in addition to its convective heat exchange with the gas, may absorb, emit and scatter thermal radiation. The radiation effect in the gas flow is neglected but the conductive heat transfer is taken into account. In order to analyze the thermal characteristics of porous burners, the coupled energy equations for the gas and porous medium in steady condition are solved numerically and the discrete ordinates method (DOM) is used to obtain the distribution of radiative heat flux in the porous media. Finally, the effects of various parameters on the performance of porous radiant burners are examined. The present results are compared with some reported theoretical and experimental results by other investigators and good agreement is found.     

Discrete ordinates solution of coupled conductive radiative heat transfer in a two-layer slab with Fresnel interfaces subject to diffuse and obliquely collimated irradiation

The coupled conductive radiative heat transfer in a two-layer slab with Fresnel interfaces subject to diffuse and obliquely collimated irradiation is solved. The collimated and diffuse components problems are treated separately. The solution for diffuse radiation is obtained by using a composite discrete ordinates method and includes the development of adaptive directional quadratures to overcome the difficulties usually encountered at the interfaces. The complete radiation numerical model is validated against the predictions obtained by using the Monte Carlo method.     

The discrete ordinates method for anisotropic scattering in the neutron transport theory: The critical sphere problem

The spherical harmonics method for anisotropic scattering in the neutron transport theory related to the critical sphere problem was investigated by Yildiz [The spherical harmonics method for anisotropic scattering in neutron transport theory: the critical sphere problem. JQSRT 2001;71:25-37]. Some numerical results and figures that they provided are incorrect. The correct numerical results for the critical radius are obtained and tabulated for different scattering parameters by using the discrete ordinates method.     

Simultaneous derivation of intensities and weighting functions in a general pseudo-spherical discrete ordinate radiative transfer treatment

The retrieval of atmospheric constituents from measurements of backscattered light requires a radiative transfer forward model that can simulate both intensities and weighting functions (partial derivatives of intensity with respect to atmospheric parameters being retrieved). The radiative transfer equation is solved in a multi-layer multiply-scattering atmosphere using the discrete ordinate method. In an earlier paper dealing with the upwelling top-of-the-atmosphere radiation field, it was shown that a full internal perturbation analysis of the plane-parallel discrete ordinate solution leads in a natural way to the simultaneous generation of analytically-derived weighting functions with respect to a wide range of atmospheric variables. In the present paper, a more direct approach is used to evaluate explicitly all partial derivatives of the intensity field. A generalization of the post-processing function is developed for the derivation of weighting functions at arbitrary optical depth and stream angles for both upwelling and downwelling directions. Further, a complete treatment is given for the pseudo-spherical approximation of the direct beam attenuation; this is an important extension to the range of viewing geometries encountered in practical radiative transfer applications. The numerical model LIDORT developed for this work is able to generate intensities and weighting functions for a wide range of retrieval scenarios, in addition to the passive remote sensing application from space. We present a number of examples in an atmosphere with O₃ absorption in the UV, for satellite (upwelling radiation) and ground-based (downwelling radiation) applications. In particular, we examine the effect of various pseudo-spherical parameterizations on backscatter intensities and weighting functions with respect to O₃ volume mixing ratio. In addition, the use of layer-integrated multiple scatter output from the model is shown to be important for satellite instruments with wide-angle off-nadir viewing geometries.     

On radiative transfer in water spray curtains using the discrete ordinates method

Radiative transfer through water spray curtains has been presently addressed in conditions similar to devices used in fire protection systems. The radiation propagation from the heat source through the medium is simulated using a 2D Discrete Ordinates Method. The curtain is treated as an absorbing and anisotropically scattering medium, made of droplets injected in a mixing of air, water vapor and carbon dioxide. Such a participating medium requires a careful treatment of its spectral response in order to model the radiative transfer accurately. This particular problem is dealt with using a correlated-K method. Radiative properties for the droplets are calculated applying the Mie theory. Transmissivities under realistic conditions are then simulated after a validation thanks to comparisons with some experimental data available in the literature. Owing to promising results which are already observed in this case of uncoupled radiative problem, next step will be to combine the present study with a companion work dedicated to the careful treatment of the spray dynamics and of the induced heat transfer phenomena.     

Radiative heat transfer between two concentric spheres separated by a two-phase mixture of non-gray gas and particles using the modified discrete-ordinates method

The radiative heat transfer between two concentric spheres separated by a two-phase mixture of non-gray gas and a cloud of particles is investigated by using the combined finite-volume and discrete-ordinates method, named modified discrete-ordinates method (MDOM), which integrates the radiative transfer equation (RTE) over a control volume and a control angle simultaneously like in the finite-volume method (FVM) and treats the angular derivative terms due to spherical geometry as the conventional discrete-ordinates method (DOM). The radiative properties involving non-gray gas and particle behavior are modeled by using the extended weighted sum of gray gases model (WSGGM) with particles. Mathematical formulation and final discretization equations for the RTE are introduced by considering the behavior of a two-phase mixture of non-gray gas and particles in a spherically symmetric concentric enclosure. The present approach is validated by comparing with the results of previous works including gray and non-gray radiative heat transfer. Finally, a detailed investigation of the radiative heat transfer with non-gray gases and/or a two-phase mixture is conducted to examine the dependence of the radiative heat transfer upon temperature ratio between inner and outer spherical enclosure, particle concentration, and particle temperature.     

Assessment of the re-ordered wide band model for non-grey radiative transfer calculations in 3D enclosures

The objective of the present study is to evaluate variations of the re-ordered wide band model for non-grey radiative transfer calculations in 3D enclosures using the discrete ordinates method. First, the performance of various angular and spatial discretisation schemes of the discrete ordinates method is investigated. Then, several formulations, averaging procedures, and scaling methods of the re-ordered wide band model are tested, and the results are validated against those of a statistical narrow band model. The grey gases formulation using three optimised absorption coefficient is found to be the most efficient method.     

A theoretical study of radiation between two concentric spheres using a modified discrete ordinates method associated with Legendre transform

A modified discrete ordinates method (DOM) is used in spherical participating media. The radiative intensity is broken up into two components. One component is traced back to the enclosure's source. It is called direct intensity. The other component is rather traced back to the contribution of the medium itself. It is called diffuse intensity. Thus, the radiative transfer equation (RTE) is transformed into two simultaneous equations: a direct RTE and a diffuse RTE. The direct RTE is solved analytically. The diffuse RTE is solved numerically using the DOM. The streaming angular derivative term appearing in spherical geometry is modeled by making use of the Finite Legendre Transform. We study a pure radiation transfer problem between two concentric spheres. The medium is assumed to be gray and isotropically scattering. The limiting spheres are considered to be opaque, gray, diffusely emitting and diffusely reflecting with uniform emissivity over each surface. The obtained results are compared with available cases reported in the literature. In particular, relative importance of the direct radiation in optically thin media is studied.     

Discrete ordinates interpolation method for radiative heat transfer problems in three-dimensional enclosures filled with non-gray or scattering medium

The discrete ordinates interpolation method (DOIM) is applied to three groups of problems of radiative heat transfer in three-dimensional rectangular enclosures containing non-gray or scattering medium. The original DOIM is first extended to a gray gas model using a new geometric interpolation scheme. It is applied to participating media for different scattering phase functions and optical thicknesses. For the non-gray gas model, the DOIM coupled with the narrow band-based weighted-sum-of-gray-gases (WSGG) model is developed. A few test problems with real gases such as pure H₂O and a mixture of CO₂, H₂O and N₂ are taken. The wall heat flux is calculated and compared with the exact solutions or reference values. All results of test problems are found to be reliable in this study. The DOIM closely reproduces the Monte Carlo reference solutions for different scattering phase functions and optical thicknesses. The non-gray gas results are compared with reference calculations based on the statistical narrow band model and they also show good agreements. The DOIM shows a remarkable merit in the computation time and the grid compatibility, to prove its usefulness for engineering applications.     

VLIDORT: A linearized pseudo-spherical vector discrete ordinate radiative transfer code for forward model and retrieval studies in multilayer multiple scattering media

We describe a new vector discrete ordinate radiative transfer model with a full linearization facility. The VLIDORT model is designed to generate simultaneous output of Stokes vector light fields and their derivatives with respect to any atmospheric or surface property. We develop new implementations for the linearization of the vector radiative transfer solutions, and go on to show that the complete vector discrete ordinate solution is analytically differentiable for a stratified multilayer multiply scattering atmospheric medium. VLIDORT will generate all output at arbitrary viewing geometry and optical depth. The model has the ability to deal with attenuation of solar and line-of-sight paths in a curved atmosphere, and includes an exact treatment of the single scatter computation. VLIDORT also contains a linearized treatment for non-Lambertian surfaces. A number of performance enhancements have been implemented, including a facility for multiple solar zenith angle output. The model has been benchmarked against established results in the literature.     

Analysis of the radiative transfer equation with highly assymetric phase function

This paper considers a scalar radiative transfer problem with high scattering anisotropy. Two computational methods are presented based on decomposition of the diffuse light field into a regular and anisotropic part. The first algorithm (DOMAS) singles out the anisotropic radiance in the forward scattering peak using the Small-Angle Modification of RTE. The second algorithm (DOM2+) separates the single scattering radiance as an anisotropic part, which largely defines the fine detail of the total radiance in the backscattering directions. In both cases, the anisotropic part is represented analytically. With anisotropy subtraction, the regular part of the signal, which requires a numerical solution, is essentially smoothed as a function of angles. Further, the transport equation is obtained for the regular part that contains an additional source function from the anisotropic part of the signal. This equation is solved with the discrete ordinates method. A conducted numerical analysis of this work showed that algorithm DOMAS has a strong advantage as compared to the standard discrete ordinates method for simulation of the radiance transmission, and DOM2+ is the best of the three for the reflection computations. Both algorithms offer at least a factor of three acceleration of convergence of the azimuthal series for highly anisotropic phase functions.     

GPU accelerated simulations of 3D deterministic particle transport using discrete ordinates method

Graphics Processing Unit (GPU), originally developed for real-time, high-definition 3D graphics in computer games, now provides great faculty in solving scientific applications. The basis of particle transport simulation is the time-dependent, multi-group, inhomogeneous Boltzmann transport equation. The numerical solution to the Boltzmann equation involves the discrete ordinates (S_n) method and the procedure of source iteration. In this paper, we present a GPU accelerated simulation of one energy group time-independent deterministic discrete ordinates particle transport in 3D Cartesian geometry (Sweep3D). The performance of the GPU simulations are reported with the simulations of vacuum boundary condition. The discussion of the relative advantages and disadvantages of the GPU implementation, the simulation on multi GPUs, the programming effort and code portability are also reported. The results show that the overall performance speedup of one NVIDIA Tesla M2050 GPU ranges from 2.56 compared with one Intel Xeon X5670 chip to 8.14 compared with one Intel Core Q6600 chip for no flux fixup. The simulation with flux fixup on one M2050 is 1.23 times faster than on one X5670.