monte-carlo方法在大气辐射传输中的应用肖辉2008202452.ppt

Monte-Carlo 方法  在大气辐射传输中的应用 肖辉 2008202452 1、蒙特卡罗方法介绍 2、蒙特卡罗方法在积分计算上的应用 3、蒙特卡罗方法在辐射传输中的应用 4、参考文献和参考网址 1、蒙特卡罗方法介绍 蒙特卡罗方法（Monte-Carlo Method）也称统计模拟方法，是二十世纪四十年代中期由于科学技术的发展和电子计算机的发明，而被提出的一种以概率统计理论为指导的一类非常重要的数值计算方法。是指使用随机数来解决很多计算问题的方法。 积分计算上应用 3、蒙特卡罗方法在辐射传输中的应用 Monte Carlo方法用于计算辐射传输时，即通过随机变量的统计试验求其近似解。 主要有两个特点： （1） Monte Carlo方法不受平面分层大气的限制，它可以解水平非均匀的问题，也可以解球面大气中的问题。 （2）Monte Carlo方法是模拟辐射在大气中传输的过程。它考察一个个光子在大气中传输时发生的散射、吸收，从统计角度来推算最终应有的辐射强度分布。 蒙特卡罗方法模拟的基本过程 光子以一个特殊的方向进入介质,确定发生碰撞时光子运动的距离,接着就能确定碰撞时是否发生了散射和吸收。若发生散射,则由适当的散射相位函数就能选取散射后新的运动方向。这个过程重复进行,直到光子被吸收或者从边界逸出为止。 (1)光子的产生 光子的位置用一个三位坐标(X,Y,Z)表示,光子的起始坐标为(0,0,0),模拟中假定为点光源。 (2)光子运动方向的确定 光子自由路径长度： 式中 上均匀分布的随机数； 消光系数： 散射方向： 其中g是非对称因子,表征散射函数的不对称性，反映了散射的各向异性程度。当g越接近-1,后向散射性越强烈；g越接近于0,散射的介质越接近于各向同性；g越接近于1,介质的前向散射越强烈。 (3)光子的新坐标 1、考虑一个电源向四周发射100000个光子 2、计算半径20cm内光子的传输，将20cm分成1000档，每档长为0.02cm 3、过程中，考虑了吸收、散射和折射 不同吸收系数 不同散射系数 不同的不对称因子 某个光子的轨迹图 某个光子的轨迹图 参考文献和参考网址 书目： 1、计算机随机模拟 裴鹿成等著 湖南科学技术出版社 1989 2、FORTRAN常用算法程序集 徐士良编 清华大学出版社 1992 文献： 1、用蒙特卡罗方法研究生物组织中的光分布 2、蒙特卡罗方法模拟光在云雾中的传输 3、蒙特卡罗方法及应用 4、水果组织中光子传输的蒙特卡罗模拟 网址： 1、/software/mc/ 2、/forum/index.php 3、.au/amiller/random.html * * 设随机点总数为N，这些点随机地落入D中任何一处，于是，落入曲边梯形内点的数目m与N之比反映了曲边梯形面积与正方形D的面积之比。 由此计算曲边梯形面积近似值 N分别取一千、两万和二十万做了10次试验 rm,rm+1 在xyz系中的方向余弦值为 ， *