Development of high-order PN models for radiative heat transfer in special geometries and boundary conditions

• Published in: Journal of quantitative spectroscopy & radiative transfer Vol. 172; pp. 98 - 109
• Main Authors: Ge, Wenjun, Modest, Michael F., Roy, Somesh P.
• Format: Journal Article
• Language: English
• Published: Elsevier Ltd 01-03-2016
• Subjects: High-order spherical harmonics; Partially diffuse and partially specular; Radiative transfer; RTE solver; Specified heat flux at the wall; Specular reflection; Partially diffuse and partially specular; Radiative transfer; Specified heat flux at the wall; Specular reflection; High-order spherical harmonics; RTE solver
• ISSN: 0022-4073; 1879-1352
• DOI: 10.1016/j.jqsrt.2015.09.001

The high-order spherical harmonics (PN) method for 2-D Cartesian domains is extracted from the 3-D formulation. The number of equations and intensity coefficients reduces to (N+1)2/4 in the 2-D Cartesian formulation compared with N(N+1)/2 for the general 3-D PN formulation. The Marshak boundary conditions are extended to solve problems with nonblack and mixed diffuse-specular surfaces. Additional boundary conditions for specified radiative wall flux, for symmetry/specular reflection boundaries have also been developed. The mathematical details of the formulations and their implementation in the OpenFOAM finite volume based CFD software platform are presented. The accuracy and computational cost of the 2-D Cartesian PN are compared with that of the 3-D PN solver and a Photon Monte Carlo solver for a square enclosure, as well as a 45° wedge geometry with variable radiative properties. The new boundary conditions have been applied for both test cases, and the boundary condition for mixed diffuse-specular surfaces is further illustrated by numerical examples of a rectangular geometry enclosed by walls with different surface characteristics.
•Derivation of the 2-D Cartesian formulation for the high order spherical harmonics (PN) methods.•Derivation of the boundary condition of PN for the mixed diffuse-specular surfaces.•Derivation of the boundary condition of PN for specified radiative heat flux at the wall.•Derivation of the boundary condition of PN for symmetry/specular boundaries.