李刚 作品数:94 被引量:186 H指数:9 供职机构: 北京应用物理与计算数学研究所 更多>> 发文基金: 国家自然科学基金 国家高技术研究发展计划 中国工程物理研究院基金 更多>> 相关领域: 理学 核科学技术 自动化与计算机技术 医药卫生 更多>>
蒙特卡罗粒子输运程序JMCT研制 被引量:46 2013年 介绍了具有自主知识产权的蒙特卡罗(MC)粒子输运程序JMCT的初步研制成果。JMCT基于三维组合几何支撑软件框架JCOGIN,采用模块化,分成多层管理结构,可处理多群碰撞和连续能量碰撞,可进行粒子并行和区域分解并行两种并行方法,并具有良好的可扩展性和高速通信技术,同时配有可视化建模前处理软件。介绍了JMCT采用的多群物质碰撞机制,展示了程序模拟计算测试模型的结果,与MCNP程序计算结果一致。JMCT串行计算速度相比MCNP提高了约3倍;在20 480个处理器核上模拟2×109样本,并行效率可达70%。 李刚 张宝印 邓力 胡泽华 马彦关键词:蒙特卡罗方法 粒子输运 三维中子-光子输运的蒙特卡罗程序MCMG 被引量:2 2013年 三维中子-光子输运蒙特卡罗程序MCMG发展了针对物质的碰撞机制,几何块、几何面动态可扩展,随机数周期进一步扩大到261。可进行多群-连续截面耦合计算,多群散射展开到P5,并考虑了中子上散射,程序配备了通用和专用多群截面库。MCMG模拟取得了与MCNP程序和实验一致的结果,串行计算速度较MCNP快2~4倍,可进行上万处理器核的并行计算。 邓力 胡泽华 李刚 李树 上官丹骅 姬志成BNCT放疗计划软件系统MDP的研制与初步验证 基于人体医学影像数据和通用蒙特卡罗程序MCNP,初步研制成功了可供BNCT研究使用的放疗计划系统的软件原型MDP.简要介绍该放疗计划系统的整体框架结构、主要功能以及所涉及到的主要技术,最后通过一个假设的治疗方案对该系统功... 陈朝斌 邓力 李刚 叶涛关键词:PHANTOM 文献传递 组合几何蒙特卡罗粒子输运支撑软件框架JCOGIN介绍 被引量:11 2013年 介绍了JCOGIN支撑软件框架层次式、模块化的体系结构及其核心数据结构。这些使得JCO-GIN软件框架能够支撑大规模复杂几何的蒙特卡罗(MC)粒子输运并行计算。同时介绍了在JCOGIN支撑软件框架下研发的JMCT程序,并进行了性能测试,在2万处理器核上模拟20亿粒子,并行效率达到70%。 张宝印 李刚 邓力关键词:并行计算 蒙特卡罗源偏倚抽样在辐射输运问题中的应用 在高温、高压、多介质、网格大变形的辐射输运计算中,粒予输运采用蒙特卡罗方法模拟,由于粒子分布不均匀,加之网格体积悬殊,导致符网格通量的统计误差涨落很大,误差积累可能致使计算结果失真。采用全局偏倚抽样,根据各网格的方差来调... 李刚 邓力 许海燕 李树关键词:蒙特卡罗方法 文献传递 基于能量密度分布的辐射源粒子空间抽样方法研究 被引量:2 2014年 利用隐式蒙特卡罗方法模拟热辐射光子在物质中的输运过程时,物质辐射源粒子是需要细致处理的物理量.传统的物质辐射源粒子抽样方法是体平均抽样方法,对于大多数问题,这样处理不会带来大的偏差.但是对于一些辐射吸收截面大、单一网格内温差显著的问题,体平均抽样方法的计算结果偏差较大.分析了产生偏差原因,提出一种基于辐射能量密度分布的辐射源粒子空间位置抽样方法,并推导了相应的抽样公式以解决此类问题.数值实验表明,新方法计算结果明显优于原方法且与解析结果基本一致. 李树 邓力 田东风 李刚关键词:输运 抽样方法 蒙特卡罗区域分解并行计算中确保串并行结果一致的伪随机数应用 被引量:5 2017年 物理建模的精细化和三维模拟给蒙特卡罗粒子输运计算的规模成千上万倍的增加,甚至超过单核内存的最大规模,仅仅依靠传统的粒子并行蒙特卡罗计算无法实现对模型模拟,区域分解并行是可能的解决方法之一.然而区域分解带来了粒子在各区域间进行迁移,导致现有的伪随机数应用方式无法确保串行计算和并行计算的结果一致.针对这种现象,本文提出赋予粒子随机数属性和动态派生次级粒子随机数的技巧来确保区域分解并行计算的串并行结果一致. 李刚 张宝印 邓力 上官丹骅 李瑞 马彦 付元光 胡小利关键词:伪随机数 蒙特卡罗方法 MCDB蒙特卡罗剂量计算系统及应用 被引量:1 2012年 硼中子俘获治疗(BNCT)蒙特卡罗剂量计算软件系统MCDB(Monte Carlo dosimetry code for brain)已经开发成功。它包括医学前处理、剂量计算和后处理。前处理把CT、MRI图像数据自动转化为剂量计算的输入文件,剂量计算基于蒙特卡罗(MC)方法,后处理是确定照射方向和照射时间。为了提高剂量计算的精度和缩短计算时间,MCDB发展了针对体素模型的快速粒子径迹算法,构造材料矩阵和计数矩阵,程序实现了MPI并行化。通过一个病例,MCDB完成了从CT、MRI提取数据、剂量计算和后处理的全过程。计算取得了与MCNP程序一致的结果,串行计算速度较MCNP提高3倍以上,并行效率可以达到90%,完全满足临床对计算精度和计算时间的要求。 邓力 李刚 陈朝斌 叶涛蒙特卡罗输运计算的组合几何区域分解算法 <正>针对核反应堆全堆芯pin-by-pin精细建模所带来的几何体数剧增以至于单核内存无法容纳的现象,在组合几何蒙特卡罗粒子输运支撑软件框架JCOGIN(J COmbinatorial Geometry MC trans... 李刚 张宝印 邓力 莫则尧文献传递 非齐次燃耗方程数值解法 被引量:1 2018年 非齐次燃耗方程常用于描述具有显著核素迁移效应的核能系统中核素含量随时间的变化规律.国内外许多燃耗计算程序无法求解方程非齐次项含时的情况.本文在方程非齐次项能够被有限阶关于时间的泰勒展开逼近这一前提下,研究了非齐次项含时情况下方程的两种解法.首先通过Laplace变换推导出了方程基于线性子链方法的解析解形式,然后使用Carathéodory-Fejér方法计算出了方程矩阵级数解的近最佳Chebyshev有理逼近式.将两种方法在燃耗计算程序JBURN中实现,并进行了数值计算,绝大部分计算结果符合很好,部分结果在较长有效数字内仍能保持一致,验证了方法的正确性和精度.同时为求解具有其他非齐次项形式的燃耗方程提供了一种思路. 付元光 邓力 邓力关键词:LAPLACE变换