蒙特卡罗算法
- 格式:pdf
- 大小:358.56 KB
- 文档页数:5
摘要:采用蒙特卡罗(Mc)计算方法,进行核反应粒子输运计算,是业内近十几年发展起来的一种有相对 优势的计算方法。本文针对并行计算过程中数据初始化、计算任务分配、计算结果归约、计算结果一致性、 程序功能一致性等关键问题进行了多种算法优劣的比较以及程序实现过程中的难点进行了分析,最终给出了 测试结果和研究结论。结果表明,Mc算法在粒子输运计算过程中有着相当好的计算效率。 关键词:蒙特卡罗算法;并行计算;计算效率 中图分类号:TL32 文献标识码:
第28卷2 O O 7 Nhomakorabea第4期 年8月
核动力工程
Nuclear Power Engineering
Vbl.28.No.4
Aug.2 O O 7
文章编号:0258-0926(2007)04—0020-05
蒙特卡罗算法并行计算研究
王家狮1,许川2,姚栋1
(1.中国核动力研究设计院核反应堆系统设计技术国家级重点实验室,成都。6l0041;2.中国核动力研究设计院,成都,6l0041)
A
1
引
言
配、计算归约、计算结果一致陛、程序功能一致
性进行了理论和实践的研究,给出了新的论述, 并实现了程序代码的编写,通过了计算测试和验
证。 2
蒙特卡罗(MC)算法是采用随机方法模拟物 理过程,应用数理统计获取计算结果的计算方法。 在核屏蔽设计计算和核反应堆临界安全分析采用 这种非确定论求解法,替代求解玻尔兹曼输运方 程,对于复杂几何结构中进行粒子输运计算,成 为一种主要趋势。这种求解方法与空间维数无关, 在求解动、静态与时间相关的二、三维复杂几何 结构粒子输运问题时,比其他确定论算法具有很 大的优势【lJ。 进行MC求解计算时,需要跟踪大量的粒子, 粒子数与计算精度成正比,这无疑需要花费相当 的求解时间。上世纪90年代以来,世界各国在高 性能并行计算机的研制、并行计算方法和软件研 究领域里,都取得了长足的发展。一般的科学研 究部门配备千亿次高性能计算机已成为现实【2】。 在这种背景下,采用并行计算方法成为解决MC 求解计算耗时问题的最佳途径。国外MC算法的 并行计算研究开展的比较早。美国oal(Ridge和
2.2需要解决的关键问题
Al锄os国家实验室都相继开发了基于MC算
法的粒子输运计算程序,成功地实现了并行计算。 国内的MC并行计算程序,已经有相当的可用性, 但在计算方法上依然存在着一些问题,可以改进 的余地很大。本文针对MC并行计算中的几个关 键技术问题:计算数据初始化、计算进程任务分
整个Mc算法并行计算研究工作需要解决下 列5个关键问题:①计算进程的数据初始化;② 计算进程任务的分配;③并行计算后数据的归约 处理;④并行计算和串行计算结果的一致性和并 行计算;⑤串行计算程序功能的一致性。
Los
MC算法并行计算研究
MC算法的并行计算研究,需要先确定算法
并行计算的可行性,再解决并行计算中的关键问 题,给出完整的解决方案,并完成程序代码的编 写和计算验证测试工作。
2.1
MC算法并行计算的可行性 Mc算法具备了良好的先天可并行化计算条
件,主要集中在以下4点:①计算数据相对独立, MC算法对每个粒子的跟踪计算,相互之问不需 要数据关联;②I/o通讯相对少,MC算法各进 程之间几乎没有I/o通讯;③循环粒度大,MC 算法的粒子计算一直处于循环状态中;④较少的 串行数据处理。MC算法在粒子跟踪计算过程外, 占用的时间不足5%。
收稿日期:200鲫6.08;修回日期:2006—08—22
万 方数据
万 方数据
万 方数据
万 方数据
万 方数据