基于球谐离散坐标法的大气辐射传输.pdfVIP

中国工程热物理学会 传热传质 学术会议论文 编号：093504 基于球谐离散坐标法的大气辐射传输 王希影，阮立明，齐宏，史萌 （哈尔滨工业大学 能源科学与工程学院，黑龙江省 哈尔滨市 150001 ） (Tel: 0451 Email: wxy102435@163.com) 摘要：本文首先应用球谐离散坐标法计算一维辐射传输模型，并将计算结果与 BMC 法进行比较，然 后选取层云为例，根据其散射特性应用球谐离散坐标法（SHDOM）计算二维算例的大气辐射传输。计 算结果包括计算域顶面出射辐射强度、上行热流密度、下行热流密度、净热流密度。结果表明：SHDOM 方法与 BMC 方法符合情况很好，二维大气辐射传输计算得到的顶面出射辐射强度的分布存在多个峰 值，热流密度的分布显示出云对地面的冷却作用和加热作用的效果。 关键词：球谐函数-离散坐标法；大气辐射；散射性质 0 前言 辐射传输过程是大气中最重要的物理过程之一，不论在民用领域还是在军用领域， 大气辐射过程的研究都具有重要的意义。在民用领域，大气辐射过程的研究可以帮助我 们预报天气，预防即将发生的气候灾难；在军用领域，大气辐射过程的研究可以帮助我 们识别目标，预警导弹。 辐射传输方程为微分—积分方程，对该方程的数值求解，包含空间坐标和空间立体 角两个离散过程。求解辐射传输模型的难点是求解内散射部分，内散射部分是由在视线 探测方向上的粒子对辐射向各个方向上的散射所形成的。不同的处理空间离散和立体角 离散的方法会产生不同的求解方法，例如离散坐标法、有限元法、热流法、离散传递法、 球谐函数法、蒙特卡洛法及有限体积法等。自 1998 年 Evens 发表球谐离散坐标法 （SHDOM ）的理论，近年来陆续有人将其应用在大气辐射传输等多个领域[1-5]。 本文首先应用球谐离散坐标法计算一维辐射传输模型，并将计算结果与 BMC 法进 行比较，然后选取层云为例，根据其散射特性应用球谐离散坐标法（SHDOM ）计算二 维算例的大气辐射传输。 1 球谐离散坐标法 辐射传输方程： dI s k ( ) s ′ ′ ′ k I s k I I s d ( ) ( , ) ( , ) ds =? e + a b + 4π ∫ ? Φ ? ? ? (1) 4π k 其中， s ′ ′ ′ ?k I (s) 表示衰减的能量，k I 表示本身辐射的能量， I (s, ? )Φ(?, ? )d ? 表 e a b 4π ∫ 4π k 示空间各个方向引起的散射源项，k I s I s ′ ′ d ′ a b + 4π ∫ ( , ? )Φ(?, ? )