Finite element analysis of concurrent radiation and conduction in participating media

A method is proposed to calculate temperature, conductive and radiative heat flux distributions in a participating medium. The method is based on the simultaneous solution of two non-linear and mutually conjugated equations describing distribution of both temperature and the so-called radiation function in the medium. In the case of isotropic scattering, the latter quantity, is proportional to the local energy density of radiation. The solution of the coupled non-linear equations is based on the finite element spatial discretization combined with the iterative technique.