Document Type




Publication Date




Source Publication

Journal of Quantitative Spectroscopy & Radiative Transfer

Source ISSN



In simulations of periodic or symmetric geometries, computational domains are reduced by imaginary boundaries that exploit the symmetry conditions. Two boundary conditions are proposed for Discrete Ordinate Methods to solve axisymmetric radiation problems. Firstly, a specularly reflective boundary condition similar to that is used in Photon Monte Carlo methods is developed for Discrete Ordinate Methods. Secondly, the rotational invariant formulation is revisited for axisymmetric wedge geometries. Correspondingly, a new rotationally invariant boundary condition specially designed for axisymmetric problems on wedge shape is proposed to enforce the rotational invariance properties possessed by the radiative transfer equation (RTE) but violated by three-dimensional conventional Discrete Ordinate Methods. Both boundary conditions have the advantage that the discretization and linear equation solution procedures of conventional three-dimensional DOM are not affected by changing to a reduced geometry. Consistency, accuracy and efficiency of the new boundary conditions are demonstrated by multiple numerical examples involving periodic symmetry and axisymmetry. A comparison between specularly reflective boundary conditions and the rotationally invariant formulation shows that the latter offers several advantages for wedge geometries. In other symmetry conditions, when the rotational invariant formulation is not applicable, specular reflective boundary conditions are still effective.


Accepted version. Journal of Quantitative Spectroscopy & Radiative Transfer, Vol. 182 (October 2016): 75-86. DOI. © 2016 Elsevier. Used with permission.

roy_9951acc.docx (571 kB)
ADA Accessible Version