![Sample our Engineering & Technology journals, sign in here to start your access, latest two full volumes FREE to you for 14 days]({"type":"image" , "src" : "/sda/5933/TOC-Subject-Sample-2021-Engineering-Tech.jpg"})
Abstract
The original radiative transfer equation is a first-order integrodifferential equation, which can be taken as a convection-dominated equation. The presence of the convection term may cause nonphysical oscillation of solutions. This type of instability can occur in many numerical methods, including the finite-difference method and the finite-element method, if no special stability treatment is used. To overcome this problem, a second-order radiative transfer equation is derived, which is a diffusion-type equation similar to the heat conduction equation for an anisotropic medium. The consistency of the second-order radiative transfer equation with the original radiative transfer equation is demonstrated. The perturbation characteristics of error are analyzed and compared for both the first- and second-order equations. Good numerical properties are found for the second-order radiative transfer equation. To show the properties of the numerical solution, the standard Galerkin finite-element method is employed to solve the second-order radiative transfer equation. Four test problems are taken as examples to check the numerical properties of the second-order radiative transfer equation. The results show that the standard Galerkin finite-element solution of the second-order radiative transfer equation is numerically stable, efficient, and accurate.
Acknowledgments
The support of this work by the National Nature Science Foundation of China (50425619,50336010) is gratefully acknowledged.