核技术应用中的蒙卡计算问题研究.pptVIP

核技术应用研究中的蒙特卡罗计算问题 报告人：许淑艳 （中国原子能科学研究院） 目录 一．引言 二．蒙特卡罗方法简介 三．蒙特卡罗方法在探测器系列参数计算中的应用 四．蒙特卡罗方法在核辐射防护计算中的应用 五．蒙特卡罗方法在核辐射医学计算中的应用 六．蒙特卡罗方法在粒子辐射效应、抗辐射加固中的应用 七．蒙特卡罗模拟方法 八．常用蒙特卡罗模拟软件 九．应用实例 一．引言 在核技术应用领域的研究中，由于受实验条件的限制，有很多问题需要借助于理论计算来完成。对于几何结构复杂且粒子与核的反应机制复杂的问题，一般数值方法无法求解，而蒙特卡罗方法能够比较逼真地描述事物的特点及物理实验过程，解决数值方法难以解决的问题，因而该方法广泛地应用于核技术应用研究中。 本文重点介绍蒙特卡罗方法[1]在几个核技术领域中的应用，以及所使用的蒙特卡罗模拟方法；蒙特卡罗方法的软件；最后，给出三个应用实例。 二．蒙特卡罗方法简介 蒙特卡罗方法又称计算机模拟方法、随机抽样技巧或统计实验方法。半个多世纪以来，由于科学技术的发展和计算机的出现和发展，这种方法作为一种独立的方法被提出来，并首先在核武器的实验与研制中得到了应用。 蒙特卡罗方法是一种计算方法，但与一般数值计算方法有很大区别，它能够比较逼真地描述事物的特点及物理实验过程，解决数值方法难以解决的问题，因而该方法的应用领域是十分广泛的。 2.1 蒙特卡罗方法的计算原理 当所求问题的解为某个事件的概率，或者是某个随机变量的数学期望，或者是与概率、数学期望有关的量时，通过随机实验的方法，得到该事件发生的频率，或者该随机变量若干个具体观察值的算术平均值，以频率代替概率，或者算术平均值代替期望值得到问题的解。 蒙特卡罗方法能够比较逼真地描述具有随机性质的事物的特点及物理实验过程。它可以部分代替物理实验，甚至可以得到物理实验难以得到的结果。 用蒙特卡罗方法解决实际问题，可以直接从实际问题出发，而不必从物理方程或数学表达式出发。它有直观形象的特点。 蒙特卡罗方法受几何条件限制较小。它的收敛速度与问题的维数无关。 2.2 蒙特卡罗方法的特点 蒙特卡罗方法具有同时计算多个方案与多个未知量的能力。 蒙特卡罗方法的计算误差容易确定，它的程序结构简单、且易于实现。 同时，蒙特卡罗方法具有以下的缺点：收敛速度慢，误差具有概率性，在粒子输运问题中，计算结果与系统大小有关。 2.3 蒙特卡罗方法的主要应用范围 蒙特卡罗方法的主要应用范围包括：粒子输运问题，统计物理。在粒子输运问题中的应用，主要包括：实验核物理、反应堆物理以及高能物理等方面，以及真空技术，激光技术、医学、生物、探矿等方面。 三．蒙特卡罗方法在探测器系列参数计算中的应用 探测器系列参数的计算包括探测器对光子的探测效率、能量沉积谱、响应函数、全能蜂以及与它们有关的量。 三．蒙特卡罗方法在探测器系列参数计算中的应用 为了把实验测得的多能光子脉冲高度谱分解成单色光子脉冲高度谱，需要对探测器进行刻度。虽然可以通过实验的方法用光子源测得其能量响应，但是这样的放射性核素是有限的，特别是高能光子源更难以得到；对于低能光子源（小于2Mev）,实验测得的结果偏高。因此，要得到对任意能量的光子的响应函数，需要借助于理论计算来完成。 三．蒙特卡罗方法在探测器系列参数计算中的应用 光子源有点源、面源和体源。 由于光子与探测器中的物质发生反应后，产生次级光子和电子，而电子和正电子在输运过程中又可产生光子。因此，这一类计算问题需要考虑光子和电子的偶合输运过程。 四．蒙特卡罗方法在核辐射防护计算中的应用 一个核设施（例如：加速器，实验室等）的建立，辐射安全问题是必需考虑的问题。粒子（中子或光子）经过屏蔽层后，有多少泄露，屏蔽层的厚度为多少才能满足防护安全的要求，以及如何设计核设施中的结构，即所谓最佳设计方案问题，需要用蒙特卡罗方法模拟。这样可以节省大量的实验时间和经费。 五．蒙特卡罗方法在核辐射医学计算中的应用 核辐射医学在治疗疑难病症（如癌症）及重症病（如心脏病）有着非常重要的作用，在放射性治疗（外照射或内照射）过程中，病体所受剂量及其分布是治疗过程中必须考虑的问题，剂量小影响治疗效果，剂量大对病体造成伤害，因此，当射线进入人体后各部分受到的剂量必须要进行理论计算。 五．蒙特卡罗方法在核辐射医学计算中的应用 关于体内剂量场分布，主要有两个途径获得：一是基于各种测试数据和经验公式，使用近似插值方法；另一种方法是利用核物理知识和模型，进行理论计算。射线在人体内的反应机制是光子和电子偶合输运过程，非常复杂，一般数值方法难于求解。蒙特卡