蒙特卡罗方法及其在中子输运问题中得应用技巧.docVIP

蒙特卡罗方法及其在中子输运问题中得应用技巧.doc

  1. 1、本文档共17页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  5. 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  6. 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  7. 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  8. 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
蒙特卡罗方法及其在中子输运问题中得应用 目录 蒙特卡罗方法及其在中子输运问题中得应用 1 1蒙特卡罗方法简介 3 1.1蒙特卡罗方法的基本原理 3 1.2 蒙特卡罗方法的误差 4 2 随机变量的抽样方法 4 2.1 直接抽样方法 4 2.1.1 离散型随机变量的抽样 4 2.1.2 连续型随机变量的抽样 4 2.2 挑选抽样法 5 2.3 复合抽样法 6 3 蒙特卡罗方法模拟中子输运过程 6 3.1 源抽样 6 3.2 输运距离的抽样 6 3.3 碰撞核素的抽样值 7 3.4 反应类型的抽样值 7 3.5 反应后中子状态的确定 7 3.5.1 弹性散射 7 3.5.2 非弹性散射 8 3.5.3 裂变反应 8 4 蒙特卡罗方法的减方差技巧 8 4.1 权 8 4.2 统计估计法 9 4.3 权窗 10 5 蒙特卡罗方法求解通量 10 5.1 通量的定义 10 5.2 点通量的计算 11 5.3 面通量的计算 11 5.3.1 统计估计法 11 5.3.2 加权法 12 5.4 体通量的计算 12 5.4.1 统计估计法 12 5.4.2 径迹长度法 13 5.4.3 碰撞密度法 13 5.4.4 几种体通量计算方法的比较 14 5.5 最终结果的统计 14 6 蒙特卡罗方法求解keff 15 6.1 有效增值因子keff的定义 15 6.2 蒙特卡罗方法求解keff 15 6.2.1 吸收估计法 15 6.2.2 碰撞估计法 15 6.2.3 径迹长度估计法 16 1蒙特卡罗方法简介 1.1蒙特卡罗方法的基本原理 蒙特卡罗方法(Mento Carlo Method)也叫统计模拟方法,是二十世纪四十年代由于计算机科学与技术发展和电子计算机的发明而提出来的一种基于概率论与数理统计的方法。蒙特卡罗方法广泛应用与金融工程、经济学、粒子输运模拟、热力学与统计物理学等领域。 为了说明蒙特卡罗方法的基本原理,先看两个例子。 例1 用蒲丰方法求解π 1977年,法国数学家蒲丰提出了一种计算圆周率的方法——随机投针法,及著名的蒲丰投针问题。这一方法如下: 取一张白纸,在上面画两条间距为2a的平行线; 取一根长度为2l(la)的针,随即的向平行线间投掷n次,观察针与平行线中得任一条相交的次数,记为m; 计算针与平行线相交的概率 可以证明这个概率的真实值为,也即当n很大时,。由此可以得到π的估计值 例2 求函数f(x)在[a,b]上的积分 f(x) [a,b]上的积分值I就是曲线y=f(x)、y轴、x=a以及x=b所围成的面积。使用蒙特卡罗方法求解该问题的方法就是:如右图所示,随机在y轴、y=M、x=a以及x=b围成的矩形均匀投掷n个点,如果落在y=f(x)下面的有m个点,则曲线y=f(x)、y轴、x=a以及x=b所围成的面积S 即积分值 由上面的两个例子可以看出,蒙特卡罗方法以一个“概率模型”为基础,将所求解的问题抽象为一个随机过程,使用已知分布抽样的方法求得试验结果的观察值,从而求得问题的近似解。 蒙特卡罗方法及其模拟粒子在物质中的输运过程 蒙特卡罗(Monte Carlo)方法又叫随机模拟法,是一种以概率论与数理统计为基础的方法。 例 用蒙特卡罗方法求解π 如下图所示,半径为1的圆外有一外接矩形,则圆与矩形的面积比为π/4。我们向矩形内均匀地投掷N个点,记录落在圆内部的点数M,则有 M/N ≈π/4 由此可以得到 π≈4M/N 由上例可以看出,蒙特卡罗方法即是根据所求的问题构造一个概率模型,使得所求问题的解等于该模型的某个变量,根据求解该变量而得到问题的解。 1.2 蒙特卡罗方法的误差 加入进行了N次模拟(对于例1,即投递了N次针,对于例2,即投掷了N个点),这N次模拟的标准差为σ,蒙特卡罗方法的误差为 2 随机变量的抽样方法 以下叙述中,f(x) 表示随机变量的概率密度函数,其分布函数为F(x),ξ为[0, 1)之间均匀分布的随机数,即ξ~U[0,1。 2.1.1 离散型随机变量的抽样 设离散型随机变量x的分布律如下: x x1 x2 x3 … xn pk p1 p2 p3 … pn 则x的抽样值为xi,其中i满足下述不等式 2.1.2 连续型随机变量的抽样 若F(x)的反函数F-1(x)存在,则随机变量x的抽样值为 例如,如果x~U[a,b],即x的分布函数为 F(x)的反函数,所以x的抽样值为 如果随机变量的概率分布函数的反函数不存在或很难求得,就应该采用其它的抽样方法。 2.2 挑选抽样法 若分布密

文档评论(0)

xuefei111 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档