深穿透粒子输运问题的自适应抽样方法

屏蔽计算中的深穿透问题一直是蒙特卡罗计算的一个难题,研究了一种发射点作为驿站的随机游动机制,推导了相应的自适应抽样方法。其主要优势在于,在蒙特卡罗方法求解粒子输运的同时,利用已经获得的信息,自适应地控制各次抽样数,不断完善计算进程。通过对碰撞点引进重要性函数,实现发射点作为驿站的重要性抽样,并结合自适应控制达到最佳抽样状态。数值结果表明:基于发射点作为驿站的自适应抽样方法,在一定程度上克服了深穿透计算中估计值偏低现象。相应的重要函数抽样方法获得了满意的结果。     

蒙特卡罗粒子输运程序JMCT研制

介绍了具有自主知识产权的蒙特卡罗(MC)粒子输运程序JMCT的初步研制成果。JMCT基于三维组合几何支撑软件框架JCOGIN,采用模块化,分成多层管理结构,可处理多群碰撞和连续能量碰撞,可进行粒子并行和区域分解并行两种并行方法,并具有良好的可扩展性和高速通信技术,同时配有可视化建模前处理软件。介绍了JMCT采用的多群物质碰撞机制,展示了程序模拟计算测试模型的结果,与MCNP程序计算结果一致。JMCT串行计算速度相比MCNP提高了约3 倍;在20 480个处理器核上模拟2109样本,并行效率可达70%。     

蒙特卡罗粒子输运程序JMCT研制

介绍了具有自主知识产权的蒙特卡罗(MC)粒子输运程序JMCT的初步研制成果。JMCT基于三维组合几何支撑软件框架JCOGIN,采用模块化,分成多层管理结构,可处理多群碰撞和连续能量碰撞,可进行粒子并行和区域分解并行两种并行方法,并具有良好的可扩展性和高速通信技术,同时配有可视化建模前处理软件。介绍了JMCT采用的多群物质碰撞机制,展示了程序模拟计算测试模型的结果,与MCNP程序计算结果一致。JMCT串行计算速度相比MCNP提高了约3 倍;在20 480个处理器核上模拟2109样本,并行效率可达70%。     

多物理耦合计算中动态输运问题高效蒙特卡罗模拟方法

多物理耦合计算在众多领域都有重要应用.如果其包含粒子输运过程,用蒙特卡罗方法模拟粒子输运常占据大部分的计算时间,因此多物理耦合计算中动态输运问题的高效蒙特卡罗模拟方法意义重大,其不可避免地依赖于大规模并行.基于动态输运问题的特点,本文提出了两种新方法:一是针对输运燃耗耦合计算的新型计数规约算法;二是动态输运计算样本数自适应算法.两种算法都能在保持计算结果基本不变的前提下使计算时间大幅减少,从而提高了效率.     

蒙特卡罗粒子输运问题的全局降方差方法

精确获得中子通量的时间-空间-能量精细分布对于一类非定常粒子输运蒙特卡罗模拟至关重要.以零方差理论为基础,采取化整为零的策略,把蒙特卡罗方法与离散纵标方法相耦合,用近似重要函数指导蒙特卡罗模拟,给出一种实现蒙特卡罗全局降方差的计算方案.数值算例表明,该计算方案获得了全局降方差的效果.     

蒙特卡罗粒子输运软件JMCT抽样工具库设计

为通用型蒙特卡罗粒子输运程序JMCT设计了抽样工具库,通过两种技术途径提供各分布的抽样。一是针对各种常见分布提供特定抽样子程序;二是提供一个通用型的抽样子程序,可以实现任意离散分布和任意一维有限区间上连续分布的自动抽样。在设计任意一维有限区间上连续分布的自动抽样工具时利用了部分开源代码,利用其功能提供给用户最大的方便性。对抽样工具库的检验表明,其可以正确、方便地实现各种输运模拟中常见分布的抽样。     

非定常输运问题适应于消息传递并行编程环境的香农熵计算方法

在多计算步的非定常输运问题的蒙特卡罗模拟中,为自动调整每一步的样本数以获得较高的计算效率,可以有多种准则.一种可选的方法是在每一步每隔若干样本监测一次系统中未死亡粒子属性分布对应的香农熵的收敛情况以决定何时停止追加样本,此种方法需要在每一步频繁计算香农熵值.由于在MPI消息传递并行编程环境下香农熵的经典计算方法必须广播大量的数据,导致每一步的计算时间随香农熵计算频率的提高而快速增大,这显然是不能满足实际需求的.本文提出了一种适应于消息传递并行编程环境的香农熵计算新方法,该方法计算得到的香农熵值并不等价于经典方法,但二者之间的差别会随着样本数的增加而趋于零.新方法的最大优势是高频计算香农熵值的时间代价大为降低,为最终实现基于香农熵收敛判断的每步样本数的自动调整奠定了必要的基础.     

非定常输运模拟中基于粒子标识分类的源偏倚算法

在非定常输运问题的多步蒙特卡罗模拟中,根据粒子的不同属性进行标识分类可以得到非常细致的系统相关标识物理量.对于某些目标标识物理量,模拟的样本中造成有效贡献的样本相对较少,导致这些物理量模拟结果的涨落较大,而单靠增加总样本数不能高效地使有效样本达到一个合理的水平.本文提出一种基于标识分类的源偏倚算法,将增加的所有样本定向赋予目标类粒子,从而高效地降低目标标识物理量的统计误差且不影响非目标标识物理量的计算.以一维多层介质非定常输运模型验证上述结论.     

非均匀大气中γ射线输运的蒙特卡罗模拟方法

γ射线在高空大气中的输运问题常采用蒙特卡罗方法进行模拟,但由于大气是非均匀连续介质,通常需要根据高度上的密度变化将大气进行分层处理.高空大气密度变化非常剧烈,精细的分层会带来很大的计算量.针对非均匀连续介质中的粒子输运问题,利用质量距离抽样方法代替传统的步长抽样方法,提出一种可直接模拟粒子在非均匀连续介质中输运的蒙特卡罗方法,实现γ射线在非均匀大气中输运的蒙特卡罗模拟,并与MCNP程序的计算结果进行比较,结果符合较好.     

非定常粒子输运蒙特卡罗并行算法研究

蒙特卡罗方法应用于粒子输运求解已是众所周知的事，然而由于MC的误差与跟踪的样本数的平方根成反比，通常要达到误差要求，需要模拟大量粒子历史，从而花费大量的计算时间。特别对非定常粒子输运问题的MC模拟，需要模拟上万个时间步，每步需要跟踪上百万个粒子。由于粒子寿命远超过时间步步长，因此，上个时间步剩余粒子的广播和本时间步剩余粒子状态量的汇总，I/O需求量大，成为制约并行的瓶颈。     

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

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

解粒子输运问题的半随机模拟数值方法

分析了解粒子输运问题的随机模拟方法和确定性方法的优缺点,提出了零方差迭代格式和半随机模拟方法——把MC方法和确定性方法有机结合起来,达到扬长避短的目的。扼要地介绍了近年来在这方面的工作。     

粒子输运方程的确定论计算方法

针对实际应用中辐射和中子输运数值模拟,讨论球一维和柱二维几何粒子输运方程确定论计算方法的研究现状,包括离散纵标、球谐函数、迭代加速、并行计算等方法.重点讨论输运计算方法所取得的若干研究进展,包括离散纵标求积组、自适应时间离散格式、本征值迭代求解方法、简化球谐函数方法、修正的子网格隅角平衡方法、灰体综合加速方法、迭代初值选取方法、输运与扩散耦合方法、基于预估校正的并行格式等.简要介绍了相关输运计算程序的研制情况,并分析输运计算方法存在的难点,提出待开展研究的内容.     

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

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

非线性中子输运问题的蒙特卡罗模拟

针对考虑中子之间相互碰撞的非线性中子输运方程,提出一种线性化近似处理方法,导出相应的积分输运方程,利用蒙特卡罗方法求解此方程,数值实验表明算法的有效性,为研究超高能中子产生与输运问题提供了必要的模拟工具.     

SN共轭函数用于蒙特卡罗粒子输运自动减方差的研究

基于S_N输运计算平台ARES编制了三维共轭输运计算模块,根据一致性共轭驱动重要性抽样方法自动生成减方差参数,用于加速MCNP5计算.数值结果表明,自动生成的减方差参数可有效提高蒙特卡罗计算效率,并保证结果无偏.自动减方差技术利用S_N共轭函数可更经济准确的估计粒子重要性,避免手动估算减方差参数的复杂工作,对于复杂屏蔽问题的蒙特卡罗计算具有较好的应用前景.     

ECR微波等离子体离子输运的数值模拟

建立了ECR微波等离子体源离子输运的平板和圆柱模型,对离子历经的空间区域的输运过程进行了数值研究。采用Monte Carlo(M-C)方法模拟了存在外磁场情况下,离子离开开放电室后历经中性区、鞘层区、最后被加负偏压的工作表面吸收的全过程,考虑了离子与中性粒子的电荷交换碰撞玫弱性散性,统一处理了中性区和鞘层区电势的衔接,采用曲线拟合,电势自洽迭代方法把中性区和鞘层区衔接起来,得到了光滑自治的电势分布曲线和鞘层区不同位置处的速度分布、能量分布及角分布。