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

蒙特卡罗方法是数值模拟中一个很常用的计算方法,应用范围很广泛,并常常作为数值仿真算法的基准.概率模型和伪随机数发生器是蒙特卡罗法中两个很重要的组成部分,它们决定了蒙特卡罗法的正确性和计算精度.本文提出了一种比较充分的检验方法-三维介质辐射对称性检验法,它可在不同条件下对各种伪随机数程序进行无限制的检验,在更精细的层次上区分伪随机数发生器的优劣,同时还可检验概率模型的正确性.     

微机上实用蒙特卡罗抽样库

描述一个实用蒙特卡罗抽样库(MCSL).它可以提供:伪随机数产生的优良算法和程序;重要分布的最佳随机抽样方法和程序;在粒子输运问题中位置、能量、方向的常用分布的抽样程序.还配有一个专门的抽样库检验系统.它在微机上运行,安装简单,使用方便,有较强的实用性和可移植性.     

蒙特卡罗与离散纵标耦合方法在压水堆堆腔漏束计算中的应用

为克服蒙特卡罗(MC)方法计算时间长和离散纵标(SN)方法复杂几何描述不精确的困难,采用SN-MC耦合计算流程,研究了基于蒙特卡罗方法和离散纵标方法的耦合方法。耦合方法的主要思想是根据离散纵标程序提供的中子角通量,利用接口程序计算出面源的分布概率,然后由修改后的源抽样子程序生成包括上下圆面源和圆柱面源的组合源,提供给蒙特卡罗程序进行计算分析。初步计算结果表明,该耦合方法是正确的,可用于压水堆堆腔漏束的计算分析。     

JMCT与子通道程序耦合方法研究及验证

蒙特卡罗与热工水力的耦合计算是目前反应堆数值模拟的重要研究方向,在蒙特卡罗方法连续能量点截面的基础上结合热工程序的温度反馈,反应堆中子计算的准确性得到大幅提高。为了提高计算精度,堆芯模型分辨率也需进一步提高,相比于组件均匀化模型,pin-by-pin的建模方式能够获得更好的结果。利用蒙特卡罗程序JMCT与子通道程序COBRA-EN实现了蒙特卡罗-热工的内耦合,内耦合方式通过内存进行数据传递,其计算效率及安全性均优于外耦合方法。随后利用NURISP项目迷你堆的pin-by-pin模型对耦合程序进行验证。计算结果与同类耦合程序相似,验证了程序的准确性。同时,对耦合过程的收敛性问题进行了初步分析。     

求解中子动力学方程的加权蒙特卡罗方法

为了实现基于蒙特卡罗方法的中子动力学计算,在传统的直接蒙特卡罗动力学方法的基础上,提出了一种加权蒙特卡罗动力学方法。该方法通过引入粒子权重的概念,隐式考虑中子俘获反应和裂变反应过程中中子数目的变化,避免了模拟粒子的数目随时间的变化,降低了统计偏差,消除了程序计算过程中粒子的存库操作,提高了计算精度。基于单能点堆模型,开发了中子动力学计算程序NECP-Dandi,进行了大量数值验证与分析,包括无缓发中子、单组缓发中子、六组缓发中子、正阶跃反应性引入、负阶跃反应性引入、正脉冲反应性、负脉冲反应性和正线性反应性引入等情况。数值结果表明,相比于直接蒙特卡罗动力学方法,加权蒙特卡罗动力学方法在计算结果的精度和计算效率上有较为明显的改进,程序结构更为简洁。     

热辐射输运问题的高效蒙特卡罗模拟方法

惯性约束聚变研究中,热辐射光子在介质中的输运以及热辐射光子与介质的相互作用是重要研究课题,蒙特卡罗方法是该类问题的重要研究手段之一.隐式蒙特卡罗方法虽然能正确地模拟热辐射在介质中的输运过程,但当模拟重介质(材料的吸收系数大)问题时,该方法花费的计算时间将变得很长,导致模拟效率很低.本文以离散扩散蒙特卡罗方法为基础,开发了"离散扩散蒙特卡罗方法辐射输运模拟程序",可以较好地解决重介质区的计算效率问题,但是离散扩散蒙卡罗方法在模拟轻介质区时精度不够高.辐射输运问题中通常既有轻介质也有重介质,为了能同时解决蒙特卡罗方法模拟的效率和精度问题,本文研究了离散扩散蒙特卡罗方法与隐式蒙特卡罗方法相结合的模拟方法,并提出了新的扩散区与输运区界面处理方法,研制了混合蒙特卡罗方法的辐射输运模拟程序.典型辐射输运问题模拟显示:在模拟重介质问题时,该程序能大幅缩短模拟时间,且能取得与隐式蒙特卡罗方法一致的结果;在模拟轻重介质均存在的问题时,与隐式蒙特卡罗方法相比,混合蒙特卡罗方法的模拟精度与其相当且计算效率同样能够得到显著提升.     

中国实验快堆燃料组件精细功率分布计算

在热工中,需要中子学计算给出燃料组件内各元件棒功率的相对分布。利用蒙特卡罗程序对中国实验快堆（CEFR）燃料组件内元件棒功率分布进行了理论计算分析,并保证计算结果相对统计误差小于0.8%。使用另一个基于六角形节块扩散理论的钠冷快堆中子学设计软件NAS计算得到的结果对蒙特卡罗程序计算结果进行了对比计算。结果表明,蒙特卡罗程序与NAS计算得到的元件棒相对功率分布结果的最大相对偏差小于3%。使用蒙特卡罗程序对CEFR燃料组件内精细功率分布的计算是可靠的,可用于设计计算当中。     

AMS02电磁量能器三维重建

研究了AMS02电磁量能器簇射三维重建, 发展出一个新的三维重建程序, 并用 蒙特卡罗模拟以及2002年束流测试数据对其进行了检验．该程序应用两次方向 重建、结合光纤衰减修正、横向以及纵向能量泄漏补偿, 完整地对事例进行了 重建, 得到了较为理想的能量分辨率和角度分辨率．     

多群输运-燃耗耦合蒙特卡罗计算

介绍了三维多群中子输运-燃耗耦合P3近似蒙特卡罗程序MCMG-BURN，该程序是在连续截面蒙特卡罗程序MCNP和反应堆栅元均匀化程序WIMS基础上发展而来的，使用多群截面为模拟临界实验堆和高通量工程实验堆，取得与MCNP和实验一致的结果，MCMG-BURN具有与MCNP相同的精度，但计算时间较MCNP要少得多。     

Fritiof程序的移植和开发

Fritiof程序是一个著名的粒子物理蒙特卡罗模拟程序,在国际上比较通用。本文简要介绍它的理论模型,功能及主要特点,描述了程序的移植、开发和应用。     

一个模拟光子照射人体的蒙特卡罗程序

本文分绍了一个用于辐射防护和医疗目的的蒙特卡罗程序。它能计算出受到光子照射时人体20种器官(组织)的吸收剂量和人体有效剂量当量。程序中使用了一种特殊的方法以加快计算结果的收敛。本文所介绍的方法也可用于其他复杂几何系统的剂量计算。     

粒子输运蒙特卡罗模拟现状概述

蒙特卡罗(MC)方法发展已有60多年历史,广泛应用于核科学及相关领域.MC方法超强的几何处理能力,使用精密的连续点截面参数,能够模拟各种复杂几何系统内的中子、光子、电子、α粒子、质子及其耦合输运问题.随着计算机的快速发展,通过大规模并行计算,MC方法及程序已成为模拟各种粒子输运问题的首选工具.     

聚变中子学的蒙特卡罗PVM并行模拟在ITER中子诊断设计中的应用

聚变能源很可能是人类文明得以维持发展的新型能源。未来的氘氚聚变堆的结构和工程设计很大程度上依赖于以聚变中子学为基础的计算。在过去的十余年中，很多的核数据库如FENDL和JENDL的检验工作围绕ITER设计而展开。聚变中子学计算包括中子和光子的输运计算。其计算目标是提供反应率和能谱等重要的信息。一维或二维的聚变中子学解析计算能提供一定精度的结果和高效率的优化设计，但对于一个三维的聚变托卡马克反应堆来说，只有蒙特卡罗方法能提供较精确的数值模拟结果。MCNP程序是由LANL实验室发展的用于中子和光子的蒙特卡罗计算的大型程序。PVM的并行计算环境能提高为MCNP程序的运行执行效率。     

纯中性事例事例顶点z坐标的最小二乘拟合

对于相互作用顶点z坐标未知的纯中性事例, 发展了一种运动学拟合的方法. 经蒙特卡罗样本检验, 拟合所用的程序是合理可靠的. 把该方法用于ψ′→γπ⁰→γγγ蒙特卡罗样本, 得到了更高的选择效率. 该方法可普遍应用于正负电子对撞实验中纯中性事例的分析.     

低能电子致DNA链断裂的纳剂量蒙特卡罗模拟

为了对DNA损伤进行纳剂量蒙特卡罗模拟,先根据D.Emfietzoglou给出的介电响应理论和Born修正模型计算了低能电子与液态水的非弹性反应截面,在Geant4-DNA的基础上建立了纳剂量蒙特卡罗程序的物理模块,同时参照国际纳剂量蒙特卡罗程序使用的模型参数建立了纳剂量蒙特卡罗程序的前化学模块、化学模块、DNA几何模块及DNA损伤模块。利用本程序,采用H. Nikjoo给出的对电子均匀照射DNA情形的模拟方法,还模拟计算了不同能量电子均匀照射DNA所导致的DNA链断裂产额,并与已公布的实验数据进行比对,二者符合良好。     

热辐射输运问题的隐式蒙特卡罗方法求解

热辐射与物质相互作用及辐射光子在物质中的传输是惯性约束聚变研究中的重要课题. 介绍了基于隐式蒙特卡罗方法的辐射输运方程,在该方程的积分-微分形式基础上,推导了利于蒙特卡罗方法模拟的等价的积分输运方程;基于积分方程设计数值模拟流程,编写三维蒙特卡罗数值模拟程序;针对热辐射输运典型问题及benchmark问题开展了数值实验,计算结果验证了方法的适应性及程序的正确性. 关键词： 热辐射 惯性约束聚变 输运方程 隐式蒙特卡罗     

MCNP-4C多粒子输运蒙特卡罗程序的MPI并行化

MCNP程序是由美国Los Alamos国家实验室研制的一个大型、多功能的粒子输运蒙特卡罗程序，可计算任意三维复杂几何系统内的中子、光子、电子或中子-光子-电子耦合输运问题，还可计算临界系统的多种本征值问题。MCNP程序的用户遍及全球，国内用户保守估计也在百家以上，过去主要应用在核科学领域，如今已推广到包括医学在内的许多领域。由于蒙特卡罗计算具有数据独立、循环粒度大、负载均衡的特性。因此，很适合作并行计算。虽然从MCNP-4A程序开始，MCNP程序具有了PVM并行计算功能，但并行系统的开发一直存在这样那样的问题，以致无法正常运行。由于MCNP程序有巨大的计算需求和计算量，我们每年使用MCNP程序完成的计算量就超过万CPU小时。许多问题采用串行计算，时间周期太长。因此，迫切需要研究缩短计算周期的并行程序。     

快堆中子学计算程序NECP-SARAX的改进与验证

NECP-SARAX是西安交通大学NECP团队开发的用于快中子反应堆的中子学程序系统。为准确处理快中子反应堆中中等质量核素散射共振以及非弹性散射导致的复杂的中子慢化效应,SARAX程序最初采用连续能量的蒙特卡罗方法计算中子能谱从而获得堆芯计算使用的有效多群截面。由于蒙特卡罗程序计算效率低,且在低能量段统计偏差较大,提出采用基于点截面的超细群方法计算中子慢化能谱,避免了蒙特卡罗方法产生参数时存在的缺陷。堆芯计算采用多群中子输运,通过优化简化几何建模,改进了程序的实用性。采用多种微扰方法计算堆芯各种反应性系数,提出了基于中子输运微扰理论的虚拟密度方法以计算堆内组件变形导致的反应性变化。在进行堆芯瞬态计算时,采用了点堆和改进准静态两种方法,可用于一般快堆和快谱ADS的典型事故分析。OECD发布的一系列快堆基准题测试表明,SARAX程序在快堆计算中具有良好的精度,达到了与国外著名快堆程序相当的水平。有效增殖因子与连续能量的蒙卡计算结果相比偏差在300 pcm以内。同时,由于引入了虚拟密度理论和三维时空动力学模型,程序功能更加完善,可以更好地满足快堆工程设计的需求。     

蒙特卡罗方法在解微分方程边值问题中的应用

介绍了蒙特卡罗方法的基本原理以及随机数的产生方法。基于蒙特卡罗方法的思想,结合有限差分方法,建立了求解微分方程边值问题的随机概率模型,并以第一类边界条件的拉普拉斯方程和一个给定初值及边界条件的非稳态热传导方程为数值算例,研究了蒙特卡罗方法在求解微分方程边值问题中的应用。结果表明:利用蒙特卡罗方法,不仅可以有效解决给定边界条件的微分方程,对于给定初值条件的微分方程,也可以从时域有限差分方程出发,采用蒙特卡罗方法进行求解。数值模拟和对误差的理论分析均表明,增加蒙特卡罗试验中的模拟粒子点数,可以提高计算结果的精度。     

激光在不同类型气溶胶中传输特性研究

激光在大气中的传输衰减特性是激光工程应用中需要考虑的一个重要问题.本文针对常用的1.06 μm和10.6μm激光,基于Mie散射理论计算了气溶胶粒子的单次散射参量;对于激光在气溶胶中多次散射传输衰减,建立了蒙特卡罗模拟计算模型,利用Matlab语言编制了相应的计算程序,计算分析了两种波长的激光分别在沙尘性、水溶性、海洋性和煤烟性四种不同类型气溶胶中透过率与传播距离、能见度的关系,并将蒙特卡罗方法和单次散射的计算结果进行了比较.结果表明,当能见度较低、气溶胶粒子反照率较高时,单次散射计算存在很大的误差,用蒙特卡罗方法更能揭示多重散射现象;煤烟性气溶胶对1.06 μm激光的传输衰减影响最大,沙尘性气溶胶对10.6 μm激光的传输衰减影响最大.