Radiative intensity with high directional and spatial resolutions can provide abundant useful information for combustion diagnosis systems based on radiative images. In this paper, an angular-spatial upwind element differential method (ASUEDM) is developed to discretize angular direction and spatial domain of radiative transfer equation (RTE) in a spherically participating medium. Because of the strong convection characteristic of radiative transfer equation, an upwind scheme is adopted to suppress the numerical oscillation. Meanwhile, Chebyshev-Gauss-Lobatto nodes are used to minimize the effect of the Runge phenomenon. Unlike conductive or convective boundary conditions, radiative boundary condition is unidirectional boundary condition, and a singularity node exists at the boundary. To deal with this singularity, we propose a discontinuous strategy. Three examples of radiative heat transfer in concentric spheres are chosen to test the capability of ASUEDM. Compared with benchmark solutions, ASUEDM can provide higher accuracy than discrete ordinates method or finite volume method. Besides, ASUEDM can flexibly provide hp convergence rate and achieve high-resolution characterization in angular direction and spatial domain.

具有高方向和空间分辨率的辐射强度可以为基于辐射图像的燃烧诊断系统提供丰富的有用信息。在本文中，开发了一种角空间逆风元微分法（ASUEDM）来离散球形参与介质中辐射传递方程（RTE）的角方向和空间域。由于辐射传递方程的强对流特性，采用迎风方案来抑制数值振荡。同时，Chebyshev-Gauss-Lobatto 节点用于最小化龙格现象的影响。与传导或对流边界条件不同，辐射边界条件是单向边界条件，在边界处存在奇异节点。为了应对这种奇点，我们提出了一个不连续的策略。选择了三个同心球中的辐射传热示例来测试 ASUEDM 的能力。与基准解决方案相比，ASUEDM 可以提供比离散坐标法或有限体积法更高的精度。此外，ASUEDM 可以灵活地提供hp收敛速度，实现角方向和空间域的高分辨率表征。