每层中的单次散射反照率和相函数为常数.ppt

DISORT模式使用 赵德龙 2007202140 主要内容： 1.Disort 算法简单介绍 2.Diosrt的功用 3.附录 在辐射传输计算中，散射大气中的辐射传输计算有其特殊的复杂性。由于散射源的存在，单色波的传输不再服从指数衰减律，而散射的非各向同性更增加了问题的复杂性。 辐射传输的离散坐标法的一般表述由Chandrasekhar(1950)和LioM(1980)给出，然而由于方法数值解的困难，难以广泛地应用于作为辐射传输的计算。 1988年stamnes等人(1988)解决了离散坐标法矩阵形式中的特征值和持征矢量求解问题以及积分常数的求解问题，同时公布了离散坐标法(Discrete Ordiate Methd)的辐射传输软件包——DISORT，这使得离散坐标法的广泛应用成为可能。 1.基本算法的介绍： 平面平行大气的单色辐射传输方程为： 源函数 其中，对满足局地热力学平衡的热源J0表示为： 直接太阳辐射的单次散射引起的增强： 这里的m0I0为入射流，得 整个方程表示为： 对相函数P(t，cosq)按2N阶Legendre多项式展开，辐射强度按Fourier阶的余弦项展开 辐射传输方程被分解为2N个独立的方程组，每个方程对应Fourier分量： 其中 这里 Pl(cosq)是Legendre多项式，Plm(m)是综合Legendre多项式。q为散射前后光线矢量方向的夹角即散射角。利用高斯求积公式可将对天顶角的积分化为求和，此时方程成为： ---1.1 wj mi 分别是求积系数和求积节点 由于单次散射反照率w和相函数P在垂直非均匀大气中是光学厚度t的函数，所以上式为2N对非常系数的微分方程组，不存在解析解。在求解，假设介质由L个相邻的均匀层组成，每层中的单次散射反照率和相函数为常数，其变化只存在于层与层之，热源项以t的多项式来近似，这样对某一均匀层，上式的矩阵形式为： 使用高斯求积公式，使得相函数满足： 这里采用双高斯法，即将高斯公式分别用于-1＜m＜0和0＜m＜1两个区域，求积节点和求积系数满足m-j＝-mj，和w-j＝-wj 矩阵方程是2N对常微分方程组，方程是线性的，可用现有的线性代数方法求解。齐次解形式为： 这是一个标准的2N×2N阶代数特征值问题，特征值：k，特征矢量：G。 由于上式的特殊构造，其待征值将成对出现(土k)所以代数特征值问题可按2的指数降阶。可得： (a-b)(G+-G-)=K(G++G-)(a+b)(G++G-) = K(G+-G-) 并有(a-b)(a+b)(G++G-)=k2 (G+-G-) 因此可得特征值k和特征矢量G 对入射光： 很容易证明1.1方程有特解。（省去角标m） 对各向同性的热辐射源，J0m =0，m0, J0≡ J00(t)=(1-w)B(t).每层的PLANK函数可近似为：B[T(t)]= ,可以找到特解： 系数Z0(ml)和Yl(ml)由线性代数方程确定。 方程1.1的解的一般形式： Cj是积分常数。j=0时： C0G0（ml）≡Z0(ml),k0≡1/m0 2. Disort功用 可求解垂直非均匀、各向异性并含热源的平面平行介质中的辐射传输问题，计算包括了热辐射、散射、吸收、下边界双向反射和发射等物理过程。也称为DISORT平面模式。 波谱可覆盖紫外（0.1mm）~微波段（5cm-10cm. 1991年DaNbackden和Stamnes(1991)介绍了DISORT软件包的伪球面和球面模式，这个基于扰动技术将球面问题简化成一系列平面问题的球面模式，对于研究曙暮光情况下的辐射传输问题非常有效。 1994年Kyllin8等人(1994)在此基础上完成了快速可靠的DISORT的二流模式，该模式的发展和持点是：1)增加了直射、散射辐射计算的伪球面模 2)增加了对随光学厚度而变化的内源的有效处理。与多流算法的对比结果表明：当太阳天顶角q＜60。地表反照率＜0.5时，地表紫外和可见光波段的辐射质量误差小于10％，而对高纬较大大阳天顶角和地表反照率的情况，误差将增至20％或更大。这个模式快速、精确、可靠，为从事辐射和粒子输运问题研究的科学家提供了一简化的软件工具。 1995年Kylling为DISORT软件包加入了大气模式参数，以及消光系数、相函数、光学厚度计算的多个子程序，使DISORT软件包真正成为用户使