光子学报, 2011, 40 (7): 1117, 网络出版: 2011-08-10
辐射传输方程的三阶球谐展开(P3)近似的有限元法求解
Finite Element Method for ThreeOrder (P3) Approximation of Radiative Transfer Equation
辐射传输方程 球谐函数 扩散方程 有限元方法 蒙特卡洛模拟 Radiative transfer Spherical harmonics Diffusion equation Finite element method Monte Carlo simulation
摘要
从辐射传输方程出发,应用球谐函数方法对辐射率进行多项式展开,并利用球谐函数的正交性和递推性,在二维笛卡儿坐标下,导出了三阶展开的球谐函数微分方程组(P3近似),改进了以往成像文献中忽略各项异性因子的P3近似,并用有限元方法对二维圆域均匀和非均匀两种情况做了数值模拟.与漫射近似模型相比较, P3近似能更准确地描述光源附近及吸收较强情况下边界的光辐射分布情况.
Abstract
The flux takes the form of an expansion in spherical harmonics, and the twodimensional spherical harmonics equations to threeorder for anisotropic scattering is drived. The equations is solved using finite element method and compared the solutions with the firstorder diffusion equation and Monte Carlo simulation. The resolut shows that the developed threeorder model with high absorb coefficient is able to significantly improve the diffusion solution in circle geometry, and the radiance distribution close to light source is more accurate.
马文娟, 高峰, 朱苹苹, 易茜. 辐射传输方程的三阶球谐展开(P3)近似的有限元法求解[J]. 光子学报, 2011, 40(7): 1117. MA Wenjuan, GAO Feng, ZHU Pingping, YI Xi. Finite Element Method for ThreeOrder (P3) Approximation of Radiative Transfer Equation[J]. ACTA PHOTONICA SINICA, 2011, 40(7): 1117.