辐射传输方程中相函数的研究.pdfVIP

辐射传输方程中相函数的研究 付娴娴杨春平吴健 (电子科技大学光电信息学院 成都610054) 摘要：精确、快速计算辐射传输方程一直是大气辐射学中的一个重要的课题。本文从辐射传输方程出发，讨论了辐射传输方 程中村j函数的计算方法，即H．G相函数展开、勒让德多项式展开和Dclta-M法展开，并指出了各种方法的优缺点。 关键词：辐射传输方程、相函数、Delta—M法 中文分类号：P422．3 文献标识符：A 1引言 由于大气辐射传输在气象学上的重要意义，也由于辐射传输问题本身的复杂性，国内外的许多学者一 直致力于研究辐射传输方程的数学求解，提出了许多算法，如球谐函数法【1】、间断坐标法121、高斯塞德迭代 法川等，通常都是对方程施以傅氏变换以得到一组与天顶余弦有关而与方位角无关的方程。在辐射传输计 算中相函数是一个重要的因子，如何处理高度不对称的相函数是一个很困难的任务。通常情况下我们利用 多项式展开的方式来求解，但是高度不对称的相函数很难用低阶多项式来近似。为了保证相函数的多项式 展开具有较高的精度，往往需要很大的展开项，甚至需要展开几百项，这在处理中所需要的计算机内存和 计算时间都要增加，而且为了避免数值计算的病态化，展开项数要求有一上刚制。常用的辐射传输精确计 算模型，比如DISORT，要求用户提供每一层介质的单次散射特性包括消光系数、单次散射反照率、不对称 因子和散射相函数。为了准确地在辐射传输模型中处理相函数，那么合适的展开技巧是很重要的。 2相函数的展开 忽略大气本身的热发射，在口I，见光、波段，均匀介质中辐射传输方程司表不为 ∥掣：，p，∥，≯) 一鲁r”删’鼬，P弘’，扎删)，㈨，≯’) 一掣P刊叩0，≯，‰丸) (1) 式中，1表示辐射强度，f为光学厚度，P为散射相函数，∥为天项角余弦，≯为方位角，∥o和九分别 为太阳天顶角余弦和方位角，090为单次散射反照率，z．F为大气上界太阳的辐照度。 2．1I卜G相函数 率相同。Henyey．Greenstein相函数，简称为H吒相函数，其表达式为 、 1——0,2 心’ ‰(cos0)2瓦g毒2Lgcos酽OU+‘一 J 其中口是散射角，2可以表示为 g=2zrP(eos和(cos功 f：lcos 渐减小。图l给出了g=0．70和0．85时的H-G相函数的变化情况 sOtilti rigAngle(degree) 图1 从图中我们可以看到H-G相函数单调地下降，前向散射也不是非常突出，而且不能描述后向散射的峰 明后向散射峰值。但是它定义的表达式中包含a、g．、92三个参数(这是它的缺陷)，增加了使用的复杂 性，于是Corette和Shanks提出了一个修正的H-G函数 嘣cos耻兰器蒜铬 ∽ 2．2用勒让德多项式(LegendrePolynomials)展开相函数 为了求维辐射传输方程，通常可以把相函数写成下面的表达式 N p(cose)=∑(2n+1)x。P．(cosO) (4) n=O 其中的系数X。是相函数关于勒让德多项式的矩，可以写成如下的表达式