平面团簇稳定结构的蒙特卡罗树搜索

以平面团簇为例提出了一种结合结构识别和蒙特卡罗树技术搜索稳定结构的新方法.体系原子之间的相互作用由两类模型势能函数来描述:Lennard-Jones二体势函数与基于Lennard-Jones势的三体势函数.考虑可能的三角晶格碎片作为候选结构,引入编号策略对结构进行快速识别,并运用蒙特卡罗树搜索研究稳定结构随着原子数增大的演化过程;对于能量较低的候选结构,进一步采取局域优化来获得对应体系的稳定结构.计算表明,Lennard-Jones二体势函数对应的三角晶格团簇更稳定;在特定的参数下,三体势函数对应的六角晶格团簇更稳定.结合结构识别和蒙特卡罗树搜索可以对候选结构空间进行高效扫描,在较短时间内更容易搜索到稳定的团簇结构,并可以与第一原理计算结合实现材料的结构预测.     

细水雾遮蔽红外辐射的蒙特卡洛模拟及透射场分析

基于蒙特卡洛法,对水雾遮蔽热辐射衰减建立了光子状态序列的追踪模型.利用该模型对水雾粒子的多重散射进行模拟,分析了影响水雾消光性能的消光系数、水雾浓度、不对称因子等参数,发现不对称因子也是影响水雾隐身性能的重要参数.研究表明:对于消光系数相同的介质,不对称因子越小越有利于衰减目标热辐射强度,而采用朗伯比尔定律则不能获得此结论;在雾粒子半径远大于红外波长的情况下,提高水雾的浓度并减小粒子半径可有效提高水雾的衰减作用.     

微波器件微放电阈值计算的蒙特卡罗方法研究

为了精准快速地计算微波器件中微放电效应的阈值,在传统蒙特卡罗方法的基础上,提出了三种不同的蒙特卡罗方法,分别对二次电子的初始能量、出射角度和初始相位等参数进行随机,结合四阶龙格-库塔法和Furman模型计算电子的运动轨迹和单次碰撞产生的二次电子发射系数,然后应用不同的方法计算有效二次电子发射系数作为微放电效应的判据.以平板传输线TEM模式为研究对象,采用四种不同的蒙特卡罗方法计算微放电阈值,并与统计模型结果进行对比.结果表明单电子-多碰撞蒙特卡罗方法误差最小,而且稳定性最好.     

圆窗片表面次级电子倍增效应的数值模拟

基于Monte Carlo模拟算法, 建立了粒子输运模型, 通过对盒形窗内圆窗片表面次级电子倍增现象进行数值仿真, 获得了TE₁₁模非均匀分布电场作用下次级电子倍增的规律. 结果表明: 在微波输入端, 指向窗片表面的磁场力起到了维持次级电子与窗片碰撞的作用, 在电场强度较高的区域倍增剧烈, 有质动力对倍增无贡献; 在微波输出端, 受背离窗片表面磁场力的影响, 在表面静电场较弱的情况下, 次级电子倍增不能发生; 当表面静电场足以维持单面倍增的发生, 随着传输功率的增大, 电子渡越时间增长, 有质动力使得倍增强烈的区域由强电场区逐渐转移到弱电场区域. 对利用外静电场抑制微波输入端次级电子倍增效应的方法进行了数值模拟验证. 关键词： 圆窗片 11模')" href="#">TE₁₁模 次级电子倍增 MonteCarlo模拟     

遮光罩内壁蜂窝结构单元的双向反射特性

内壁蜂窝结构是决定遮光罩杂散光抑制性能的主要因素,掌握其各向异性辐射传输特性对遮光罩的杂散光分析和优化设计具有重要意义.针对蜂窝结构单元,建立辐射传输模型,导出了表面双向反射分布函数的离散表达式;通过蒙特卡罗法模拟蜂窝结构单元的辐射传输过程,分析了几何参量和涂层反射率对等效面反射特性的影响;根据模拟所得双向反射分布函数数据库,建立了用于遮光罩杂散光分析的蜂窝结构等效面的反射特性概率模型.分析结果表明,蜂窝结构对杂散光呈现很强的后向散射特征,等效面的反射特性概率模型与直接模拟结果符合很好,可用于遮光罩杂散光分析和设计.     

中子诱发232Th裂变初始碎片质量及动能分布Monte-Carlo研究

随着第四代反应堆以及先进核能利用系统的发展,对中子核数据提出了高精度、多核素、宽能区的新要求.目前,中国核数据评价库(CENDL库)中相关核裂变的数据较缺失,不足以满足当前核能发展的需求.因此,建立面向中子核数据需求的可靠计算方法和工具变得极为重要.本文基于Monte-Carlo方法建立了裂变碎片质量动能计算模型,研究了中低能中子诱发²³²Th(n,f)反应发射中子前裂变碎片的分布特性.对于裂变碎片质量分布,本模型计算结果与实验值最大偏差约1%,与GEF, TALYS程序计算结果 (与实验值最大偏差约2%)相比具有一定优势.对于发射中子前裂变碎片动能分布,本模型计算结果与实验数据一致.结果表明,所发展的计算模型能够较好地预测²³²Th(n,f)反应发射中子前裂变碎片数据,为中子诱发锕系核裂变反应计算提供一种新思路.     

压水堆核电站燃料组件弯曲对堆芯中子学的影响分析计算

在压水堆核电站中,由于燃料组件装配的压紧力、冷却剂流动、辐射蠕变、燃耗等因素会导致燃料组件的弯曲,燃料组件的弯曲对组件间的水隙分布产生影响,从而影响中子的慢化行为及堆芯的传热性能,进而对反应堆堆芯的运行参数造成影响。本文分析了组件弯曲的成因及机理、影响及后果（包括对堆芯功率分布、径向功率倾斜、焓升因子、热点因子等参数的影响）,并使用蒙特卡罗软件JMCT,对组件弯曲的确定论计算程序的正确性进行了验证。最后通过确定论的计算程序模块,对CPR1000核电站的组件弯曲情况进行了模拟分析,计算结果表明:在某一燃耗下,随着水隙增加或减小,燃料组件功率会随之增加或减小,使堆芯的功率分布发生倾斜,影响核电站的安全运行。     

基于VHTRC的棱柱式高温气冷堆核设计程序验证北大核心CSCD

高温气冷堆是国际公认的固有安全性高的反应堆堆型。针对高温气冷堆包覆颗粒燃料引入的燃料组件的双重非均匀性以及棱柱式堆芯布置的非均匀性和强空间耦合效应,提出基于蒙特卡罗均匀化-确定论输运方法的RMC-SaraGR程序系统作为棱柱式高温气冷堆的核设计程序。基于日本棱柱式高温气冷堆临界实验装置VHTRC基准题,针对此套核设计程序系统开展了均匀化模型研究和初步验证。研究结果表明,基于蒙特卡罗均匀化方法,采用全堆模型、合适的能群结构和分区方式产生组件群常数,并经过超级等效均匀化方法进行等效均匀化修正,可以保证堆芯多群均匀计算具有较高的计算精度。     

磁绝缘线振荡器内气体碰撞电离数值模拟

针对气体碰撞电离过程,介绍了蒙特卡罗碰撞（MMC）的处理方法,利用MMC方法编写了气体碰撞电离模块,将其移植到3维全电磁粒子模拟程序NEPTUNE之中,模拟了充有He气的磁绝缘线振荡器（MILO）。模拟结果表明:当He气密度较低时,电离的正离子由于较重无法自由移动,形成了正离子通道,可以有效中和电子束空间电荷场,有利于电子束传输和群聚,提高了束波互作用效率,微波输出功率得到了明显提高,起振时间也有所缩短;当进一步增加He气密度时,电离碰撞增强,电子和离子数目会雪崩式增长,电子束由于碰撞增强而导致能散度增大,其负效应已经远大于中和空间电荷场的正效应,反而不利于电子束的群聚和共振,从而导致输出微波功率降低乃至截断,起振时间缩短是由于其在非雪崩阶段的正效应积累所致,但是随着负效应的增强起振功率不能得以维持,二极管最终将闭合。另外,还模拟了MILO填充空气、水蒸气及二氧化碳等多原子、多组分气体的碰撞电离物理过程。模拟结果显示,同压强情况下,填充空气、水蒸气及二氧化碳的脉冲缩短现象要比填充He气等较低原子序数气体的情况严重得多。     

1+1维Wolf-Villain模型奇异标度行为的数值模拟

采用Kinetic Monte Carlo(KMC)方法对描述分子束外延生长(MBE)的1+1维Wolf-Villain模型进行大尺寸和长生长时间的数值模拟研究,以消除渡越行为的影响.计算得到整体和局域标度指数.结果显示,在所模拟的空间和时间尺度范围内,1+1维Wolf-Villain模型仍呈现出固有奇异标度行为.这一结论与López等人最近的理论分析结果不一致.