- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
蒙卡特罗方课件5
第四章 蒙特卡罗方法解粒子输运问题 屏蔽问题模型 直接模拟方法 简单加权法 统计估计法 指数变换法 蒙特卡罗方法的效率 作 业 第四章 蒙特卡罗方法解辐射屏蔽问题 辐射(光子和中子)屏蔽问题是蒙特卡罗方法最早广泛应用的领域之一。本章主要从物理直观出发,说明蒙特卡罗方法解决这类粒子输运问题的基本方法和技巧。而这些方法和技巧对于诸如辐射传播、多次散射和通量计算等一般粒子输运问题都是适用的。 屏蔽问题模型 在反应堆工程和辐射的测量与应用中,常常要用一些吸收材料做成屏蔽物挡住光子或中子。我们所关心的是经过屏蔽后射线的强度及其能量分布,这就是屏蔽问题。 当屏蔽物的形状复杂,散射各向异性,材料介质不均匀 , 核反应截面与能量、位置有关时,难以用数值方法求解,用蒙特卡罗方法能够得到满意的结果。 粒子的输运问题带有明显的随机性质,粒子的输运过程是一个随机过程。粒子的运动规律是根据大量粒子的运动状况总结出来的,是一种统计规律。蒙特卡罗模拟,实际上就是模拟相当数量的粒子在介质中运动的状况,使粒子运动的统计规律得以重现。不过,这种模拟不是用实验方法,而是利用数值方法和技巧,即利用随机数来实现的。 为方便起见,选用平板屏蔽模型,在厚度为 a,长、宽无限的平板左侧放置一个强度已知,具有已知能量、方向分布的辐射源 S 。求粒子穿透屏蔽概率(穿透率)及其能量、方向分布。穿透率就是由源发出的平均一个粒子穿透屏蔽的数目。 同时,假定粒子在两次碰撞之间按直线运动 , 且粒子之间的相互作用可以忽略。 直接模拟方法 直接模拟方法就是直接从物理问题出发,模拟粒子的真实物理过程。 状态参数与状态序列 模拟运动过程 记录结果 状态参数与状态序列 粒子在介质中的运动的状态,可用一组参数来描述,称之为状态参数。它通常包括:粒子的空间位置 r, 能量 E 和运动方向Ω,以 S=( r , E ,Ω ) 表示。 有时还需要其他的参数,如粒子的 时间 t 和附带的权重W ,这时状态参数 为 S=( r , E ,Ω , t ,W ) 。 状态参数 通常要根据所求问题的类型和所用的方法来确定。 对于无限平板几何,取 S=( z , E , cosα) 其中 z 为粒子的位置坐标,α为粒子的运动方向与 Z 轴的夹角。 对于球对称几何 , 取 S=( r , E , cosθ) 其中 r 表示粒子所在位置到球心的距离,θ为粒子的运动方向与其所在位置的径向夹角。 粒子第 m 次碰撞后的状态参数为 或 它表示一个由源发出的粒子,在介质中经过 m 次碰撞后的状态,其中 rm :粒子在第 m 次碰撞点的位置 Em :粒子第 m 次碰撞后的能量 Ωm:粒子第 m 次碰撞后的运动方向 tm :粒子到第 m 次碰撞时所经历的时间 Wm :粒子第 m 次碰撞后的权重 有时,也可选为粒子进入第 m 次碰撞时的状态参数。 一个由源发出的粒子在介质中运动,经过若干次碰撞后,直到其运动历史结束(如逃出系统或被吸收等)。假定粒子在两次碰撞之间按直线运动,其运动方向与能量均不改变,则粒子在介质中的运动过程可用以下碰撞点的状态序列 描述: S0 ,S1 ,…,SM-1 ,SM 或者更详细些 , 用 来描述。这里 S0 为粒子由源出发的状态,称为初态,SM 为粒子的终止状态。M 称为粒子运动的链长。 这样的序列称为粒子随机运动的历史,模拟一个粒子的运动过程,就变成确定状态序列的问题。 模拟运动过程 为简单起见,这里以中子穿透均匀平板的模型来说明,这时状态参数 取 S=( z , E , cosα)。 模拟的步骤如下: (1) 确定初始状态 S0 : 确定粒子的初始状态,实际上就是要从中子源的空间位置、能量和方向分布中抽样。设源分布为 则分别从各自的分布中抽样确定初始状态。 对于平板情况, 抽样得到 z0=0。 (2) 确定下一个碰撞点 : 已知状态Sm-1,要确定状态Sm,首先要确定下一个碰撞点的位置 zm。在相邻两次碰撞之间,中子的输运长度 l 服从如下分布: 对于平板模型,l 服从分布: 其中,Σt 为介质的中子宏观总截面, 积分 称为粒子输运的自由程数, 系统的大小通常就是用系统的自由程数表示的。 显然,粒子输运的自由程数服从指数分布, 因此从 f ( l ) 中抽样确定 l,就是要从积分方程 中解出 l。 对于单一介质 则下一个碰撞点的位置 如果 zm≥a,则中子穿透屏
您可能关注的文档
- 英文文献翻译硼作显色剂在分光光度计法中的应用一项有关氨甲苯酸与硼酸为定作用的新颖的褪色分光光度计法.doc
- 英文温度湿度和剪口龄期对腐烂病和轮纹病菌自剪锯口侵染和发病的影响.docx
- 英文缩400个.doc
- 英文翻译一种新的计算方法用于实双工件在冲压排样设计中的最佳规划方法.doc
- 英文翻译中国五强溪船闸边坡计和施工的工程地质力学分析和监测控制.doc
- 英文翻译从难处理矿石中氰化浸提取黄金的新工艺.doc
- 英文翻译关于用PVT集热器的设计和实验研究.doc
- 英文翻译在板料成形过程实验的特征性采用维氏硬技术1.doc
- 英文翻译基于flah型单片机二进制模数转换器的温度计解码器.doc
- 英文翻译应力穿越节理时的透反射研究.docx
最近下载
- 三阶魔方F2L公式研究及记忆方法.pdf
- 装表接电工考试(技师)习题库(第1部分).pdf VIP
- 幼儿园-健康-危险的食品(预防食物中毒)-课件(互动版).pptx
- 人教版小学四年级音乐下册《第二单元 五十六朵花》大单元整体教学设计[2022课标].docx
- Q∕GDW 11316-2018 高压电缆线路试验规程.pdf
- 心房颤动患者左心耳封堵术后护理查房.pptx VIP
- 0723S03016-监察法学-理论课程教学大纲-申静、李芹.docx VIP
- 阜外心血管病医院手术分级目录2015版.docx VIP
- DB5328T 44-2025 珠芽黄魔芋栽培技术规程.docx VIP
- 2025年辽宁事业单位考试真题及答案.docx
文档评论(0)