蒙特卡罗模拟中三维介质辐射对称性检验研究.pdfVIP

中闰r秤热物耻学会 传热传质学 学术会议论文 88 编0：063I 蒙特卡罗模拟中三维介质辐射对称性检验 帅永，谈和平 (蚧尔滨T．ql,人学舰窄航灭热物理∥f究所，黑龙江哈尔滨150001) Tel：045l一8641 2308·Email：shuaiyong@hit．edu．cn 摘要：蒙特卡罗厅法足数值模拟中一个很常用的计算方法，心用北…很广泛．并常常作为数值仿真鳕： 法的堆准。概牢模型和伪随秽I数发生器足蔹特卡罗法中两个很重要的组成部分．它们决定r蒙特{}罗 注的jF确性羁l汁算精度。本文挺…J，一种比较充分的榆验方法——三维介质辐射对称性检验法，它口J． iJilq．．1条件下对再种伪随机数程序j2[行无限制的榆验，4：更精细的层次，J：区分伪随机数发生器的优劣， Jrd时还n¨惫验概半模型的正确性。 关键词：蒙特卡罗法：介质辎射：概牢模型；伪随机数 0引言 蒙特g-罗(MonteCarlo)法也叫做随机仿真方法或随机抽样方法，是一种基于“随机 数”的计算方法l|l。蒙特譬罗方法是i：卜。股数值计算方法有本质区别，属丁．试验数学的 ～个分支，起源丁．甲．期的心儿率近似概率的数学思想，它利刚随机数学进行统计试验， 以求得的统计特征值作为待解问题的数值解。蒙特墨罗方法已经在原子弹．1：科的科学研 究中发挥了极其重要的作J{j，并止住日黼J．泛地麻_Hj丁．物理、l：群、经济、金融的各个 方面。它的基本思想是：为了求解问题，首先建立一个概率模型或随机过祥，使它的参 数等丁．问题的解：然后通过对模型或过程的观察或抽样试验来计算所求随机参数的统计 特征，最后给山所求解的近似值。 对J：辐射传输数值仿真，将热辐射的传输过程分解为发射、透射、反射、吸收和散 射等1一系列独立的子过稗，并建立每个一r过祥的概率模刑。令每个单元(面元和体元) 发射一定域的光束(能术、光子)，跟踪、统计每束光束的门宿(被介质和界面吸收，或 从系统中透射出或逸出)，从而得剑该单元辐射能鬣分配的统计结果。蒙特R罗方法以概 率统计理沦为其主要理论基础，以随机抽样(随机变鼙的抽样)为其主要手段。概率模 喇的止确性和伪随机数发生器的随机均匀性决定了蒙特p罗法模拟计笄的止确性和运算 精度。本文±占合介质辐射传输特性，任¨二维介质辐射埘称性枪验法P—l基础上提i13了…种 ·叮以对概；簪模7啦和伪随机数发乍器进{丁检验的方法——二i维介质辐舅t埘称性枪验法，它 ¨，亿／1：nd条仆卜．时符种伪随机数程序进行无限制的检验，侄．更精细flqE-；k—t-Ⅸ分伪随机 数发生器的优劣，同时还T叮检验蒙特g-罗概率模型的止确性。 基金项目：Iq家[J然科学旗金tti项tt(No资ⅡJJ 526 l三维介质辐射对称。眭检验模型 二三维介质辐射对称性检验法的 基本思路赴取一个立方体腔体结构， X 网格数为～×NN，如图1所示．腔 体，}，『11】允满j，参j0性介质，从中心f眭 置点(图中阴影标识部分。取其坐标 为fa以“】)通过伪随机数程序向空间 随机发射一定数目(Ⅳ㈨I)的辐射能 求，而后统计经过一次及多次传递后 箨朋格吸收的能求射线数目(N¨，t)。 符介质的衰减、(吸收、散射)系数均 图1+二维介质辐射对称性检验法模刑 匀，n四刷壤面的发射率相等，则理 论上儿何位置对称网格吸收的能求数应相等。为简便起见，本文中不考虑散射，则衰减 系数等『．吸收系数。 (1) √(c，一‘．，，★)2+(口一咒，，，^)2+(口一刁，，。★)2=DoD5石-a，1．i,j，kN 式(1)农示。