This work presents an OpenFOAM implementation of the PN approximation for radiative heat transfer, including higher orders P3, P5, and P7. Also described is a procedure which enables the sequential numerical computations of the coupled partial differential equations (PDEs) by re-expressing the boundary conditions in matrix form so that individual boundary conditions can be associated with each PDE. The implementation of the software programs are verified with derived analytical solutions for 1-D slabs with constant and variable properties, and are also tested with various orientations in order to demonstrate the geometric invariance properties of the 3-dimensional PN formulation. A few examples taken from the literature are also considered in this work and could be taken as benchmark solutions for the PN approximations.

This content is only available via PDF.
You do not currently have access to this content.