核技术应用中的蒙卡计算问题

核技术应用研究中的

蒙特卡罗计算(jìsuàn)问题

报告人：许淑艳（中国(zhōnɡɡuó)原子能科学研究院）精品资料目录(mùlù)一．引言二．蒙特卡罗方法简介三．蒙特卡罗方法在探测器系列参数计算中的应用四．蒙特卡罗方法在核辐射防护计算中的应用五．蒙特卡罗方法在核辐射医学计算中的应用六．蒙特卡罗方法在粒子辐射效应(xiàoyìng)、抗辐射加固中的应用七．蒙特卡罗模拟方法八．常用蒙特卡罗模拟软件九．应用实例精品资料一．引言(yǐnyán)在核技术应用领域的研究中，由于受实验条件的限制，有很多问题需要借助于理论计算来完成。对于几何结构复杂(fùzá)且粒子与核的反应机制复杂(fùzá)的问题，一般数值方法无法求解，而蒙特卡罗方法能够比较逼真地描述事物的特点及物理实验过程，解决数值方法难以解决的问题，因而该方法广泛地应用于核技术应用研究中。本文重点介绍蒙特卡罗方法[1]在几个核技术领域中的应用，以及所使用的蒙特卡罗模拟方法；蒙特卡罗方法的软件；最后，给出三个应用实例。

精品资料二．蒙特卡罗方法(fāngfǎ)简介蒙特卡罗方法又称计算机模拟方法、随机抽样技巧或统计实验方法。半个多世纪以来，由于科学技术的发展和计算机的出现和发展，这种方法作为一种独立的方法被提出来，并首先在核武器的实验与研制中得到了应用。蒙特卡罗方法是一种计算方法，但与一般数值计算方法有很大区别，它能够比较逼真地描述事物的特点及物理(wùlǐ)实验过程，解决数值方法难以解决的问题，因而该方法的应用领域是十分广泛的。精品资料2.1蒙特卡罗方法(fāngfǎ)的计算原理

当所求问题的解为某个事件的概率，或者是某个随机变量的数学期望，或者是与概率、数学期望有关的量时，通过随机实验(shíyàn)的方法，得到该事件发生的频率，或者该随机变量若干个具体观察值的算术平均值，以频率代替概率，或者算术平均值代替期望值得到问题的解。精品资料蒙特卡罗方法能够比较(bǐjiào)逼真地描述具有随机性质的事物的特点及物理实验过程。它可以部分代替物理实验，甚至可以得到物理实验难以得到的结果。用蒙特卡罗方法解决实际问题，可以直接从实际问题出发，而不必从物理方程或数学表达式出发。它有直观形象的特点。蒙特卡罗方法受几何条件限制较小。它的收敛速度与问题的维数无关。2.2蒙特卡罗方法(fāngfǎ)的特点精品资料2.2蒙特卡罗方法(fāngfǎ)的特点蒙特卡罗方法具有同时计算多个方案与多个未知量的能力。蒙特卡罗方法的计算误差容易确定，它的程序结构简单、且易于实现。同时，蒙特卡罗方法具有以下(yǐxià)的缺点：收敛速度慢，误差具有概率性，在粒子输运问题中，计算结果与系统大小有关。精品资料2.3蒙特卡罗方法的主要(zhǔyào)应用范围蒙特卡罗方法(fāngfǎ)的主要应用范围包括：粒子输运问题，统计物理。在粒子输运问题中的应用，主要包括：实验核物理、反应堆物理以及高能物理等方面，以及真空技术，激光技术、医学、生物、探矿等方面。精品资料探测器系列参数的计算包括(bāokuò)探测器对光子的探测效率、能量沉积谱、响应函数、全能蜂以及与它们有关的量。三．蒙特卡罗方法(fāngfǎ)在探测器系列参数计算中的应用精品资料三．蒙特卡罗方法在探测器系列参数(cānshù)计算中的应用为了把实验测得的多能光子脉冲高度谱分解成单色光子脉冲高度谱，需要对探测器进行刻度。虽然可以(kěyǐ)通过实验的方法用光子源测得其能量响应，但是这样的放射性核素是有限的，特别是高能光子源更难以得到；对于低能光子源（小于2Mev）,实验测得的结果偏高。因此，要得到对任意能量的光子的响应函数，需要借助于理论计算来完成。精品资料光子源有点源、面源和体源。由于光子与探测器中的物质发生反应后，产生次级光子和电子，而电子和正电子在输运(shūyùn)过程中又可产生光子。因此，这一类计算问题需要考虑光子和电子的偶合输运(shūyùn)过程。三．蒙特卡罗方法在探测器系列参数(cānshù)计算中的应用精品资料四．蒙特卡罗方法在核辐射防护(fánghù)计算中的应用一个核设施（例如：加速器，实验室等）的建立，辐射(fúshè)安全问题是必需考虑的问题。粒子（中子或光子）经过屏蔽层后，有多少泄露，屏蔽层的厚度为多少才能满足防护安全的要求，以及如何设计核设施中的结构，即所谓最佳设计方案问题，需要用蒙特卡罗方法模拟。这样可以节省大量的实验时间和经费。精品资料五．蒙特卡罗方法在核辐射医学计算(jìsuàn)中的应用核辐射医学在治疗疑难病症（如癌症）及重症病（如心脏病）有着非常重要的作用，在放射性治疗（外照射或内照射）过程中，病体所受剂量及其分布是治疗过程中必须考虑的问题，剂量小影响治疗效果，剂量大对病体造成伤害，因此，当射线(shèxiàn)进入人体后各部分受到的剂量必须要进行理论计算。精品资料关于体内剂量场分布，主要有两个途径获得：一是基于各种测试数据和经验公式，使用近似插值方法；另一种方法是利用核物理知识(zhīshi)和模型，进行理论计算。射线在人体内的反应机制是光子和电子偶合输运过程，非常复杂，一般数值方法难于求解。蒙特卡罗方法是解决该类问题十分有效的方法。五．蒙特卡罗方法在核辐射医学(yīxué)计算中的应用精品资料六．蒙特卡罗方法在粒子辐射(fúshè)效应、抗辐射(fúshè)加固中的应用在核聚变堆、空间卫星及飞行器中，其关键部分的元器件、晶体受到高能粒子（中子、光子、电子、质子等）辐射是造成各种元器件失灵、失控乃至事故的主要原因。因此，计算各种元器件材料的辐射损伤，对于材料的辐照筛选，设施的抗粒子辐射加固和安全运行(yùnxíng)有重要的意义。精品资料高能粒子辐射到晶体后，通过与核反应，产生多种粒子（如中子、光子、电子(diànzǐ)、质子、α粒子等）和反冲核，而其中的每一种粒子、反冲核通过各种反应道又可生成各种粒子与核。这种复杂的偶合输运问题不仅需要考虑粒子（不带电与带电）与核的反应，还需要考虑原子碰撞。六．蒙特卡罗方法在粒子(lìzǐ)辐射效应、抗辐射加固中的应用精品资料七．蒙特卡罗模拟(mónǐ)方法对于单个粒子（中子或光子），其模拟方法(fāngfǎ)通常有：直接模拟法，加权法，统计估计法、指数变换法等。精品资料2.两种粒子（比如(bǐrú)：光子与电子）偶合输运问题对于两种粒子（比如：光子与电子）偶合输运问题，使用“字典编辑式”的多分支方法。该方法是：在光子的输运过程(guòchéng)中，一旦有光子、电子产生，将光子存储，跟踪电子。模拟的电子历史结束后，依照后进先出的原则取出光子，进行模拟。直至所有的次级光子和电子的历史全部结束，一个由源出发的光子的历史结束。精品资料光子-电子的偶合输运多分支(fēnzhī)方法精品资料3.光子—电子—核—核偶合输运(shūyùn)的多分支方法该方法(fāngfǎ)是：在光子—电子偶合输运的多分支方法(fāngfǎ)的基础上，嵌入核—核偶合输运的多分支方法(fāngfǎ)。即：在光子—电子偶合输运的过程中，如有反冲核产生，即转向核—核偶合输运。当核—核偶合输运结束后，返回到光子—电子偶合输运的过程中去。精品资料光子-电子-核-核的多分支(fēnzhī)方法精品资料4.粒子—粒子—核—核偶合(ǒuhé)输运的多分支方法

该方法实现多种粒子的偶合输运(shūyùn)问题。粒子—粒子偶合输运(shūyùn)过程与光子—电子偶合输运(shūyùn)的多分支方法类似，由于粒子的种类不同，需要加入描述粒子类型的存储单元。这个方法是在粒子—粒子偶合输运(shūyùn)的基础上嵌入核—核偶合输运(shūyùn)过程。精品资料粒子(lìzǐ)-粒子(lìzǐ)-核-核偶合输运的分支方法精品资料八．几个常用(chánɡyònɡ)的蒙特卡罗模拟软件MCNP软件包EGS4软件包GEANT4程序包精品资料1.MCNP软件包

MCNP的全名是：MonteCarloNeutronandPhotoTransportCode.它是由美国LosAlamos实验室研制的一个大型的多功能的蒙特卡罗程序包。可用于计算中子、中子光子、光子电子及其组合的输运问题，以及临界(línjiè)（包括次临界(línjiè)和超临界(línjiè)）系统的本征值计算问题。是目前国内外普遍使用的程序包。精品资料MCNP—3B是1989年公开发表的版本，目前使用的版本大多为MCNP—4A，MCNP--4B，MCNP—4C版本。使用MCNP软件包需要(xūyào)填写输入卡片：INP。计算结果放在输出文件OUTP中。1.MCNP软件包

精品资料2．EGS4软件包

EGS（Electro一GammaShower）软件包是光子(guāngzǐ)一电子在任意几何中偶和输运的蒙特卡罗摸拟通用软件包。由美国斯坦福直线加速器中心编制，EGS4软件包是较新版本。EGS4软件包括四个部分精品资料第一部分(bùfen)：PEGS4一截面数据处理程序

它的功能是给出系统介质所需的截面及介质中各种元素的截面与各种反应的分支比。截面数据是使用分段线性拟合方法产生的。使用该程序时，首先需要按照规格要求填写输入卡片（输入文件(wénjiàn)），文件(wénjiàn)名为：USERINP。DAT，然后运行执行程序PEGS4。计算结果在输出文件(wénjiàn)FOR012.DAT中。运行过程在文件(wénjiàn)EGSOUT.DAT中可以看到。精品资料第二(dìèr)部分：EGS4软件包该软件包由两部分组成：SHOWER和HATCH。SHOWER是模拟光子一电子偶合输运过程的子程序，它在运行过程中又调用其他子程序和两个用户子程序HOWFAR（几何(jǐhé)输运程序）和AUSGAB（记录输出程序）；HATCH是输入介质数据子程序。精品资料第三部分：PRO—公共(gōnggòng)区文件

包括(bāokuò)八个公共区，可根据需要使用其中有关部分。精品资料第四部分(bùfen)：例题EGS4软件包给出了七个例题，供用户(yònghù)使用与程序检验。精品资料EGS4软件包的使用(shǐyòng)方法使用(shǐyòng)EGS4软件包，需要编写用户程序。其中包括主程序（及其所需要的若干子程序），和EGS4程序包所要求的两个子程序：HOWFAR子程序和AUSGAB子程序。精品资料3.GEANT4程序包Geant4是由欧洲核子中心主导开发的一套用于MonteCarlo模拟程序包.来自于美国、俄罗斯、日本、加拿大等国家10多个实验室的100多名科学家参与了Geant4程序的研制工作.它包括了实验装置的描述,粒子在材料和磁场(cíchǎng)中的运动,以及粒子与物质相互作用的物理过程模型等一整套工具包.精品资料3.GEANT4程序包它的一个(yīɡè)突出特点是它包含了非常丰富的物理模型,并且能够在非常大的能量范围内处理粒子与物质的相互作用.正是由于它的粒子种类多，物理模型全，能量范围大这些特点，使得他的应用领域越来越广泛，包括高能物理，核试验，加速器，医学，生物科学，辐射防护等多个领域。精品资料九．应用(yìngyòng)实例实例一、聚变堆第一壁材料SiC中电子辐射损伤实例二、高纯Ge探测器在不同几何上测量的“有效作用深度”的研究与计算实例三、I-125体内(tǐnèi)放射源剂量场分布计算精品资料实例(shílì)一、聚变堆第一壁材料SiC中电子辐射损失

计算了不同能量的入射电子在不同剂量下的位移(wèiyí)损伤强度，得到了入射能量为14Mev的电子，在不同的辐照温度，不同的损伤剂量及不同的高温弛豫时间，热扩散系数随测量温度的变化情况。

精品资料实例(shílì)一由计算结果可看出：当辐照温度低于700°C时，热扩散系数随损伤剂量的增加而减小，但对于(duìyú)损伤剂量为0.3dpa时，其热扩散系数比前者大，原因不仅是由于辐照温度高（770°C），而且由于高温弛豫时间长（13.3min）。当辐照温度800°C时，损伤剂量为0.1dpa时，其扩散系数却比第一、二种情况减少了许多，唯一的不同是高温弛豫时间短（6.3min），因此，可以得到一个重要结论：当高温弛豫时间大于6min时，能够减少对晶体的损伤。精品资料精品资料实例(shílì)一

1995年10月在莫斯科举行的“第七届国际聚变堆材料会议”上报告了该种实验设计、计算方法与计算结果，引起了较大的反响(fǎnxiǎng)，认为这一成果在国际上处于先进位置，实验、计算与计算方法首次使用。精品资料实例二、高纯Ge探测器在不同几何上测量的“有效作用深度”的研究