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

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

基于JASMIN框架多物理耦合程序的性能优化及分析

基于并行应用支撑软件框架JASMIN的辐射流体与粒子输运耦合程序RHSn2D,采用最小邦元固定处理器数目的并行策略,计算实际模型的并行规模扩展至8192核,并行效率约为16%.集成程序时间分析,验证软件框架底层MPI并行环境聚合通信对于并行优化算法(尤其是辐射流体计算时间)的影响.     

用于全电磁粒子模拟的复杂建模及网格生成技术

 基于Open CASCADE开源的计算机辅助几何设计类库,给出了全电磁粒子模拟中复杂模型构建技术及网格生成技术。介绍了Open CASCADE软件的基本特点;给出了基于Open CASCADE进行自主研发的用于全电磁粒子模拟的复杂模型构建软件EasyEMModeling的设计思想、程序框架及已具备的功能;基于Open CASCADE中的射线与几何体的求交算法,给出了Yee网格中共形描述3维复杂模型的共形网格生成技术;最后,给出了验证实例,证实了共形网格生成方法的正确性及有效性。     

基于JASMIN并行框架的2.5维粒子模拟程序NEPTUNE2D的研制北大核心CSCD

介绍了2.5维自主研制的并行电磁粒子模拟程序NEPTUNE2D初步研发情况。该程序基于JASMIN并行自适应结构网格支撑框架研制,并行效能高,可扩展性强,且支持动态负载平衡;采用新型PIC算法替代传统算法,避免求解泊松方程修正电场,更适用于大规模并行计算;程序支持r-z坐标系下的器件仿真,可应用于高功率微波器件、电真空器件的快速模拟设计。该程序现已完成电磁场更新、粒子推进、电磁场注入/引出、粒子发射/吸收等基本物理功能模块的研制,并通过同轴线、圆波导、同轴二极管及无箔二极管算例模拟验证了模块的正确性。最后,应用NEPTUNE2D程序设计了一个高效同轴相对论返波管,给出了粒子模拟结果和并行性能测试结果。     

复杂产品虚拟试验支撑框架

VITA支撑框架作为虚拟试验验证的共性技术支撑平台,主要用于解决复杂航天产品研制试验过程中的异构性、实时性、互联互通等问题。本文提出了一套完整的VITA虚拟试验体系框架,讨论了它的研究背景及体系结构,研究了VITA运行中间件的技术架构,并进行了性能测试,介绍了基于中间件的扩展服务,介绍了VITA目前的应用情况,初步应用表明,VITA支撑框架为复杂产品研发过程提供了有力的虚拟试验技术支持,可在航空、航天等多个行业的复杂产品研发中进行应用。     

并行自适应结构网格应用支撑软件框架研制进展

JASMIN框架全称为“并行自适应结构网格应用支撑软件框架（Jparallel Adaptive Structured Mesh Applications Infrastructure）”。该框架面向科学与工程计算的（自适应）结构网格应用,支持应用领域专家在无需了解并行计算和共性数值算法细节的前提下,通过编写串行数值计算程序,     

基于钯/碳纳米粒子复合薄膜的双氧水电化学传感性能研究

本文采用团簇束流沉积方法制备了一种复合纳米粒子电化学催化剂,在碳纳米粒子支撑层上沉积钯纳米粒子薄膜,发现其在双氧水电化学传感中具有较高的灵敏度.碳纳米粒子的覆盖率对钯纳米粒子薄膜的双氧水电化学催化活性有明显的影响.当碳纳米粒子覆盖满一个单层的时候,钯/碳纳米粒子复合薄膜对双氧水的检测灵敏度达到了最高值,是没有碳纳米粒子支撑层时的两倍之多.     

车载光电稳定平台外框架优化设计

为使车载稳定平台在动态载荷状态下满足系统精度要求,对车载光电稳定平台外框架进行了优化设计。介绍了外框架材料的选择原则与外框架的结构形式;运用UG软件对外框架进行三维建模;采用理论分析与有限元计算相结合的手段,对外框架结构进行了分析与设计;并基于MSC软件进行了仿真校核。结果表明：本文设计的外框架能在满足刚度强度的前提下...     

不对称二聚物粒子在两嵌段共聚物形成的支撑膜上的自组装

结合自洽场理论和密度泛函理论,研究了不对称的二聚物粒子在AB两嵌段共聚物形成的支撑膜上的自组装行为. 不对称的二聚物粒子是由两个不同的球组成的两性分子. 其中一个球与A嵌段相亲,另外一个喜欢B嵌段. 不对称粒子能在支撑膜上形成一个双层结构. 由于衬底的存在,所形成的支撑膜上下两叶的对称性被破坏,导致了二聚物粒子在膜上形成的结构也变的不对称. 随着二聚物粒子浓度的增加,在膜上的二聚物粒子的结构将发生变化,从稀四方、六角、密四方再到圆柱结构. 在一个高浓度的密堆积下,二聚物粒子将形成弯曲的圆柱结构. 在支撑膜     

3维并行全电磁粒子模拟软件UNIPIC-3D

介绍了自行研制的3维并行全电磁粒子模拟软件UNIPIC-3D。在该软件中,电磁场量用二阶精度的时域有限差分方法迭代,粒子用相对论牛顿-洛仑兹力方程推进。该软件拥有复杂器件的几何建模和网格自动剖分的功能,具有模拟相对论返波管、虚阴极振荡器、磁绝缘线振荡器等高功率微波器件的能力。且该软件具有强大的后处理功能,可以显示电场、磁场、电流、电压、功率、频谱、粒子相空间等。在高性能并行计算机上对软件的并行效率进行了测试。通过与2.5维UNIPIC软件的结果进行比较,验证了UNIPIC-3D软件的正确性。     

BEAVRS基准模型热零功率状态的JMCT分析

使用JMCT (J Monte Carlo Transport Code)对来自MIT的全堆芯pin-by-pin精细建模的国际基准模型BEAVRS的热零功率(HZP)状态进行了模拟计算, 并与测试数据进行了对比和分析. 比较的物理量包括临界本征值、控制棒价值、反应性温度系数、轴向积分的全堆探测器测量值和不同位置四个组件轴向相对功率密度分布. HZP状态下不同控制棒位置插入和硼浓度的临界本征值计算, JMCT结果与理论值1.000 的误差小于0.2%, 控制棒价值计算结果与测量值符合. JMCT对轴向积分的探测器径向相对功率分布和四个组件的轴向归一化的探测器的计算结果与测量值进行了比较和分析, 计算结果与测量值一致, 同时清晰地展示了模型增加格架后, 轴向功率曲线在相应位置出现下凹的现象. 此外, JMCT给出了轴向积分的组件径向相对功率密度分布和轴向相对功率最大处(Z轴位置)的pin径向相对功率密度分布, 并与国际知名程序MC21结果进行了对比, 两个图像都符合得非常好. 随着计算机与并行计算的高速发展, 蒙特卡罗程序开始从传统的反应堆校验工具向反应堆设计工具转变.     

基于VENUS-III国际基准模型的JMCT程序验证

基于自主研制的三维中子-光子耦合输运蒙特卡罗通用程序JMCT（J Monte Carlo Transport Code）,采用连续点截面,对国际基准屏蔽VENUS-III模型开展精细建模和中子输运临界及屏蔽计算.临界计算得到系统k_eff、重要区域的通量及能谱.结果表明,JMCT和MCNP程序的重要区域体通量计数吻合较好,偏差均在1%以内.深穿透屏蔽计算采用外源模式,点探测器计数,JMCT计算值与实验测量值偏差在15%以内,满足屏蔽设计对误差的要求.初步验证了JMCT程序临界及屏蔽计算的可用性.     

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

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

蒙特卡罗输运模拟软件JMCT的深穿透屏蔽计算

核设施辐射屏蔽计算,由于其大规模计算及深穿透等特性,一直是蒙特卡罗方法工程应用的难点之一。采用我国自主研发的三维中子-光子蒙特卡罗粒子输运模拟软件JMCT,结合可视化建模工具JLAMT,对OECD国际基准例题Winfrith Iron/Water Benchmark Experiment(ASPIS)两例实验装置进行建模与计算分析, 并将计算结果与实验值及MCNP计算值进行对比。结果表明,JMCT计算值与MCNP计算值符合较好,其中Winfrith Iron Benchmark Experiment(ASPIS)最大偏差不超过7%,平均偏差1.3%;Winfrith Water Benchmark Experiment(ASPIS)最大偏差小于20%,平均偏差小于10%,证明了JMCT在屏蔽计算以及深穿透问题的可靠性与工程应用性。     

JMCT程序临界安全基准校验计算与分析

JMCT程序是北京计算数学与应用物理研究所自主开发的蒙特卡罗方法粒子输运计算程序。从国际临界安全实验数据库中选取棒栅类、U系、Pu系、233U系实验方案对JMCT程序进行验证,验证涵盖了常见裂变元素（如U,Pu等）从低富集度到高富集度的快、热能区,同时也包含了常见的中子毒物以及反射层材料。将JMCT程序的计算结果与基准实验值进行对比,并且与MCNP,MONK程序的计算结果进行了比较。结果表明,在检验的范围内,JMCT程序具有与国际通用蒙特卡罗方法粒子输运程序一样良好的计算精度。     

基于JLAMT可视化建模的JMCT模拟计算

JMCT是中国工程物理研究院高性能数值模拟软件中心粒子输运团队自主研发的三维蒙特卡罗模拟软件,JLAMT为其前处理可视化建模工具。使用JLAMT和JMCT程序对BW,KRITZ,BEAVRS等系列基准题进行了模拟,并对有效增殖系数及计数功能进行了对比分析。其中有效增殖系数计算结果与MCNP5的最大偏差为KRITZ2装载方案19基准题的89.1 pcm,除BEAVRS基准题外的计数结果与MCNP5的偏差基本小于2%,平均偏差在1%左右;BEAVRS基准题功率分布模拟结果与MCNP5及实测值最大偏差分别为7.06%,16.6%,控制棒价值计算与MCNP5及实测值均吻合较好。,     

JMCT优化查找策略的实现

分析连续能量点截面中子输运程序JMCT中区间查找的时间复杂度.针对主能量网格,裂变释放中子数和不可分辨共振网格进行区间映射加速,并以等差和等比序列网格为例,对比分析不同散列网格的质量.最后选择寿期初与寿期中燃料棒栅元的k_∞问题进行对比测试.结果显示,对于核素数目较多的输运问题,优化查找策略具有很好的加速效果.     

反应堆蒙特卡罗临界模拟中均匀裂变源算法的改进

在反应堆pin-by-pin精细建模及蒙特卡罗模拟计算研究中, 由于不同栅元的功率密度差异较大, 导致蒙特卡罗方法临界计算的样本在不同栅元之间的分配不均衡, 由此引起栅元内的各种计数的统计误差差异较大. 为使大部分栅元内计数的统计误差降至一个合理的水平, 单纯增加总样本已不是一个高效的解决方法. 通过在特定临界计算迭代算法的基础上改进并实现均匀裂变源算法的思想, 对大亚湾压水堆pin-by-pin模型取得了具有较高效率的数值结果. 本工作为具有自主知识产权的蒙特卡罗粒子输运模拟软件JMCT最终达到反应堆pin-by-pin模型(包括一系列国际基准模型)的模拟性能要求提供了一个有效的工具.     

蒙特卡罗粒子输运软件JMCT的网格计数功能设计与实现

基于蒙特卡罗输运的网格计数能够统计指定区域内密集的分布计数量.为实现JMCT统计局部计数量分布,设计并实现了JMCT网格计数功能,该功能支持三种正交几何（xyz直角几何、rθz柱几何、rθφ球几何）的均匀及不均匀剖分;介绍xyz直角几何网格计数的算法;基于大亚湾反应堆pin-by-pin模型、Venus模型以及ITER次临界包层模型等初步验证了网格计数功能;应用铀阵列模型,采用若干种网格剖分办法,在单机上对比测试了xyz直角几何下JMCT和MCNP5网格计数功能的串行性能,结果显示JMCT的计数时间较短,具有较高的效率.