第三章：辐射传输方程方案.ppt

参考式： 对上式从0 到 τ0 积分： 即： 整理得I(0, Ω) 与 I(τ 0, Ω) 之间的关系： 对上式的解释： 位于τ=0处的辐射强度由两部分组成： τ= τ 0处的辐射强度穿过整层介质而经过衰减的值， 整层介质中的每个辐射源被衰减后到达τ=0处的辐射强度的总和。 I(0, Ω) 与 I(τ 0, Ω) 之间的关系也可以表述为： 请注意，此时μ0，若将其变为正数，则上式可变为： 对上式的解释： 位于τ= τ 0 处的辐射强度由两部分组成： τ= 0 处的辐射强度穿过整层介质而经过衰减的值， 整层介质中的每个辐射源被衰减后到达τ= τ 0处的辐射强度的总和。 源函数只考虑介质发射情况下的解 当源函数只考虑介质发射时，辐射传输方程相对考虑散射时要简单得多，因为它不需要考虑各方向散射辐射因素，而且J 与I 无关。此时的辐射传输方程可以写为： B(T)为普朗克函数，是物体亮温为T时发射的出射辐射亮度，它的强度与出射方向无关，即各向均一。 总结 辐射传输方程的求解是对 τ 的积分，而J 与I 是否有关决定了求解难易。 不考虑源函数的解为比尔-布格-朗伯定律，只考虑发射的解也相对简单。 注意辐射传输方程中单次散射项也与I 无关： 下一小节将重点解决该问题。 源函数J与待求强度I无关时的解 单次散射解与散射逐次计算法 二流 (two-stream) 近似 不考虑发射和多次散射，仅考虑源函数为单次散射情况时的传输方程为： 此时源函数与待求强度I 无关，可套用上一小节的解法，即上式可转换为： 其中 参照上一小节的解： 代入 即可求得仅考虑源函数为单次散射情况时的传输方程的解。 散射的逐次计算方法 散射的逐次计算方法是这样一种方法，我们单独对散射一次、二次、三次等的光子计算其强度，而总强度则为所有各次散射之和。 式中 n表示散射的次数。 注意到多次散射的源函数为： 由于二次散射是由一次散射引起的，因而从一次散射强度 I1(τ, Ω)即可求出二次散射源函数： 而二次散射强度是可以由其源函数计算出来的： 同样我们可以由二次散射强度推导出三次散射源函数，继而推出三次散射强度。 依此类推，我们可以得到任意次散射的强度，其递归关系式可以表示为： 注意是对τ 积分, 还是对Ω积分 利用递归关系式可以设计数值方法，逐级导出源函数和强度，进而根据： 求出包含多次散射的总强度。 总结 在辐射传输方程中，单次散射源函数J 与待求强度I 无关，可以求出解析解。 单次散射解中的第 1 项反映了比尔-布格-朗伯定律，有时也称为零次散射解，而将第 2项，即对源函数的积分结果称为单次散射解。 利用逐次计算方法可以依次得到各次散射的源函数和强度，进而求出考虑多次散射的方程解。 源函数J与待求强度I无关时的解 单次散射解与散射逐次计算法 二流 (two-stream) 近似 辐射传输与方位无关时的简化 观察我们已熟知的辐射传输方程(不考虑发射)： 当辐射传输与方位 φ无关，而仅与 μ有关时，注意到此时 ，则有： 注意 μ = cos θ 勒让德（Legendre）展开 散射相函数可以表征为散射角余弦的函数： 上式可以用有限N项的勒让德多项式进行展开： 其中 l 阶勒让德多项式： 前几阶的勒让德多项式为： 针对P(cosΘ)进行勒让德多项式展开的系数为： 前 2 阶的展开的系数为： 注意： P为散射相函数，Pl 为勒让德多项式的 l 阶展开， 二者符号差不多，不要搞混。 引入不对称因子： 对各向同性散射，g为零； 当相函数的衍射峰变得越来越尖锐时，g也随之增大； 若相函数峰值位于后向，g为负值； (1+g)/2 可以看作积分前向散射能量的百分比数； (1-g)/2 可以看作积分后向散射能量的百分比数。 不对称因子 g 观察与方位 φ无关时的辐射传输方程 对其进行勒让德多项式展开，有： -μ0中的负号用于表征方向 为了用解析方法求解上式，必须用有限个求和代替积分。业已发现，对于区间[-1，1]上的求积分，可用高斯公式展开，即对任何函数 f (μ)，有： 式中权重值 其中μj是偶阶勒让德多项式P2n (μ)的零点，并有： 考虑： 用高斯公式展开后得： 式中μi(-n,n)代表辐射流的方向。 取n=1，则得到两个辐射流，即j=-1和1，此时 N=1，并且相应的高斯点和权重值分别为 重排各项并以I↑表示I(τ, μ1)和I↓表示I(τ, -μ1)后，可导出两个联立方程，即 式中 上式即为二流近似的辐射传输方程，它可以得到解析解，这里不继续推导其求解过程，有兴趣者可以翻看相关参考资料 当取n=2时，即可得到四个辐射流，列出4个方程，称为四流近似。同样，我们可以采用六流近似、