删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

蒙特卡罗临界计算全局计数问题新策略研究

本站小编 Free考研考试/2021-12-29

摘要:基于一个裂变源分布对应香农熵序列的在线收敛性诊断方法, 为提高全局计数整体效率的均匀裂变点算法将在首次激活迭代步和首次判断收敛迭代步的最大值之后启动. 之后, 整体精度指标将每隔固定迭代步数计算一次, 一旦达到事先设定的精度标准, 整个临界迭代计算将提前终止. 这一判断过程将一直持续到事先设定的最大迭代步数为止. 通过一个反应堆基准问题C5G7模型的计算表明上述新策略在合适的参数下有助于提高蒙特卡罗临界计算全局计数问题的整体效率.
关键词: 临界问题/
全局计数/
蒙特卡罗方法/
香农熵

English Abstract


--> --> -->
计算核相关系统的临界性在核科学与工程领域是常见而重要的问题. 作为一种成熟的计算方法, 蒙特卡罗算法由于其几何、物理建模的精确性以及高度的可并行性等特点在临界计算中得到广泛的运用.
在最近的研究[1-4]中, 由于计算机硬件与模拟计算方法、软件的巨大进步, 使得过去由于能力不足而采取的近似手段有了可以进一步精确化处理的能力, 例如反应堆设备的pin-by-pin计算等. 这些趋势反过来对于计算方法和软件提出了更大的挑战. 对于临界问题的蒙特卡罗模拟而言, 由于核设备的精细化建模及多物理耦合计算需求, 在计算临界特性(即特征值)的同时进行大规模的全局计数成为一个难以高效解决的问题[5-7]. 困难一方面来自于计数规模的庞大; 另一方面来自于计算目标不能仅仅瞄准于最大程度地减少最小误差, 而是必须提高全局计数的整体效率, 这就使得单纯增加样本数是一种低效甚至无效的策略, 因为在功率较高、误差较小的区域更容易得到新增样本, 而功率较低、误差较大的区域反而得不到新增样本, 从而使得整体效率难以提高.
为提高临界计算全局计数问题的整体效率, 出现了一些算法上的研究. 这些研究以95/95准则[5] (即要求在一定时间内, 至少95%的栅元内计数的相对误差要以至少95%的置信度位于1%以下)为指导, 力图在总样本数一定的前提下获得更高的整体效率, 其中均匀裂变点(UFS)算法[5-7]是最早提出的一种. 该算法在假设裂变源本征分布已经达到的前提下, 通过统计已完成迭代步的裂变点数密度, 以此来偏倚当前迭代步的裂变次级平均粒子数, 然后通过纠偏的手段保持结果的无偏性, 最终的效果是在统计意义上使得计算误差较大的区域产生更多的低权重粒子而在计算误差较小的区域产生更少的高权重粒子, 从而使得整体效率得到较大的提高. 但是如果在裂变源分布还未收敛的情况下启动UFS算法, 由于统计得到的裂变点密度具有系统误差, 将会导致算法应用失败, 解决这一问题的方法一般是凭经验选择一个足够大的非激活迭代步数来保证裂变源分布的收敛. 同时, 由于只在完成所有的设定迭代步计算后才统计整体效率, 即使计算早已达到事先设定的整体效率标准, 也不得不进行多余的计算. 在后续的研究中还出现了效率更高的均匀计数密度算法[8,9].
本文提出了一种进一步提高临界计算全局计数问题整体效率的新策略. 这一策略在自主研发的定态蒙特卡罗粒子输运模拟软件JMCT上进行了验证. 本文第2节将介绍包含两个部分的新策略; 第3节给出了相应的数值结果; 第4节给出了结论.
这一新策略包含两个部分. 首先, 基于一个裂变源分布对应香农熵的实时收敛性诊断方法, 为提高全局计数整体效率而设计的UFS算法将在首次激活迭代步和首次判断收敛迭代步的最大值之后被激活, 这就为保证裂变源分布已经收敛提供了双重保障, 从而保证了UFS算法所使用的偏倚数据更加合理; 其次, 在UFS算法被启动后, 将定期监测全局计数的一个整体精度指标, 一旦这个指标小于事先约定的数值, 则认为全局计数已经很好地收敛, 整个计算将被终止, 而不必等到事先设定的最大迭代步数全部完成. 下面将详细介绍这两部分.
2
2.1.基于香农熵的裂变源分布实时收敛性诊断方法
-->香农熵最先引入蒙特卡罗临界计算是作为判断裂变源分布是否收敛的后验指标, 同时, 也存在另一些基于其他熵的实时收敛性诊断方法研究[10-12]. 基于香农熵的实时收敛性诊断方法在减少非定常输运问题的计算时间方面获得了成功的应用[13]. 将这一实时收敛性诊断方法应用在定态蒙特卡罗临界计算方面还属首次.
裂变源分布对应香农熵H的基本定义如下[10]:
$H = - \sum\limits_{i = 1}^N {{P_i}{{\log }_2}({P_i})} , $
式中, N为含裂变物质栅元的数目; ${P_i}$$\displaystyle\frac{{{M_i}}}{M}$, 其中M为所有裂变源中子的个数, ${M_i}$为裂变源中子其空间坐标位于第i个含裂变物质栅元中的总个数.
基于香农熵的实时收敛性判断准则依赖于多个香农熵值组合出的随机振子指标${K_n}$, 定义为[11]
${K_n}=\frac{{{H^n} - H_{\min }^{n, p}}}{{H_{\max }^{n, p} - H_{\min }^{n, p}}}.$
其中, ${H_n}$是当前第n个迭代步的香农熵值, $H_{\max }^{n, p}$$H_{\min }^{n, p}$分别是当前第n迭代步之前所有p个香农熵中的最大值和最小值. 在香农熵值序列已经收敛的前提下, ${K_n}$将在0.5附近做无规随机涨落. 实时判断香农熵是否收敛的准则在于判断当前第n迭代步及之前所有m个随机振子指标是否满足不等式[13]
$\left| {{K^{n - m + 1}} - 0.5} \right| < \varepsilon {\text{且}}\left| {\frac{1}{m}\sum\limits_{i = 0}^{m - 1} {{K^{n - i}}} - 0.5} \right| < \varepsilon, $
一旦满足, 则认为裂变源分布已经收敛, 可以在n和事先约定的非激活迭代步的最大值之后启动计数过程及UFS算法; $p,\; m,\; \varepsilon $按参考文献[11]的推荐取为20, 50, 0.1.
2
2.2.基于全局计数整体精度指标定期监测的提前结束计算方法
-->衡量全局计数整体效率的95/95准则的本质在于以尽可能少的计算时间获得尽可能高的整体精度. 整体精度指标其实有很多种, 都从一个方面反映了全局计数的整体收敛状况. 95/95准则使用的整体精度指标是Pre_95, 即所有计数的相对误差中至少95%的相对误差都不大于该值. 要在并行计算环境下定期监测这一整体精度指标, 就必须在激活UFS算法后每隔固定迭代步数计算一次所有计数的相对误差, 而通常来说, 并行蒙特卡罗程序需要进行规约之后才能得到计数的相对误差, 一般仅在所有迭代步完成之后进行规约并输出计数. 为解决这一问题, JMCT采用了对象序列化与在线反序列化技术(由于JMCT采用面向对象的编程方式实现, 程序运行时将产生多个对象实例, 其内存排布是无序的. 但重启动及续算功能需要将对象保存在文件中, 需要内存对象是有序且连续的. 在数据备份时对无序对象进行的操作就是序列化, 在恢复运行时对数据的操作就是反序列化), 以最小的I/O开销实现了并行环境的备份与恢复, 从而支持在任意迭代步中进行计数规约操作.
以C5G7模型(图1)作为验证上述策略的基准模型. 该模型是NEA发布的程序检验基准例题. 其中包含两种组件, 分别是铀氧化物UOX组件和铀钚混合氧化物MOX组件, 交叉2 × 2布置. 每个组件内有17 × 17个pin-cell, 其中又分为264个燃料pin、24个导向管pin和一个仪表管pin. UOX组件中, 包含一种富集度; MOX组件中包含三种富集度. 该模型大约包含2万个栅元, 其详细介绍可以参考文献[14].
图 1 C5G7模型图
Figure1. C5G7 model.

以全局体平均通量计数作为计算目标. 由于没有标准答案, 所以, 通过大规模1000迭代步的计算(每代5百万样本, 跳过前500代)获得基准解, 这一计算规模对于该问题是足够的. 图2所示裂变源分布对应香农熵的变化反映了这一点. 分四种情况进行了计算: 1)不启动策略和UFS算法(简称为no_stra_no_UFS); 2)不启动策略、启动UFS算法(简称为no_stra_with_UFS); 3)启动策略、不启动UFS算法(简称为with_stra_no_UFS); 4)启动策略和UFS算法(简称为with_stra_with_UFS). 最初的设置是1500迭代步, 跳过200代, 每代20万样本数. 图3显示了不启动策略时有无UFS算法裂变源分布所对应香农熵的变化, 可以看出, 200代的非激活迭代数目是不够的, 也就是UFS算法利用了不够精确的数据进行偏倚. 为考察全局计数的整体精度, 设计了整体精度指标E_global, 定义如下:
图 2 基准解裂变源分布对应香农熵的变化
Figure2. Shannon entropy of fission source distribution for benchmark result.

图 3 不启动策略时有无UFS算法的裂变源分布对应香农熵的变化
Figure3. Shannon entropy of fission source distribution for two cases (without strategy and with or without UFS algorithm).

$E{\rm{\_global }}=\sum\limits_{i = 1}^{{\rm{cell\_number}}} {{{\left(\frac{{{T^i} - T_{{\rm{benchmark}}}^i}}{{T_{{\rm{benchmark}}}^i}}\right)}^2}}, $
其中, $T_{{\rm{benchmark}}}^i$是第i个栅元内体平均通量的基准解, ${T^i}$是各种情况下第i个栅元内体平均通量的计算结果. 这个指标从一个方面度量了全局计数的整体精度高低, 越小表明整体精度越高.
当启动策略时, 原先约定的在第201步启动计数或UFS算法被推迟到第224步, 显示香农熵的实时诊断方法判断直到第224步裂变源分布才收敛得比较充分. 如果要求只要整体精度指标Pre_95不大于0.004就结束整个计算, 且每隔40步判断一次, 则不启动UFS算法时在第944步达到要求, 而启动UFS算法时在第824步就达到了要求. 从表1可以看出, 各种情况下计算得到的特征值keff基本没有什么变化, 说明各种方法得到的结果都是无偏的. 由于不启动策略时UFS算法启动得过早(在第201步启动), no_stra_with_UFS这种情况下的整体计算精度E_global是较差的, 即使其已经完成了全部1500个迭代步的计算. 对于with_stra_no_UFS情形, 虽然启动计数过程较晚, 裂变源分布收敛得已经比较充分, 但由于总共只有944个迭代步, 有效样本较no_stra_no_UFS情形少很多, 又没有UFS算法使全局计数的相对误差平均化, 所以整体精度指标E_global是最差的, 即使此时指标Pre_95已经小于0.004. 而对于with_stra_with_UFS情形, 虽然结束计算的时刻更早(在第824步), 从而有效样本数更少, 但由于启动UFS算法对全局计数误差的平均化效应, 其整体精度指标E_global是最好的.
方法keff计算时间/sE_globalPre_95
no_stra_no_UFS1.176608816.00.76890.00284
no_stra_with_UFS1.1765010189.21.25670.00259
with_stra_no_UFS (在第224启动计数, 在第944步结束)1.176575580.41.98000.00396
with_stra_with_UFS (在第224启动计数, 在第824步结束)1.176294831.10.56160.00396


表1结果比较
Table1.Comparison of results.

对于蒙特卡罗临界计算全局计数问题, 基于已有的UFS算法, 提出了一种进一步提高效率的新策略. 该策略通过基于裂变源分布对应香农熵的实时判断收敛准则和对全局计数整体精度指标的定期监控, 可进一步确保UFS算法获得质量更高的数据并且在保证精度的前提下减少冗余计算, 从而提高了整体效率. 对C5G7基准模型的计算表明本文提出的策略在参数合适时是高效的. 当然, 由于新策略中还有经验参数, 利用更智能化的方法消除这些参数是下一步工作的目标.
相关话题/计算 指标 统计 数据 程序

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 基于冗余图态的多人协作量子计算
    摘要:量子计算是一种基于量子力学基本原理设计的新型计算模型,在某些特定问题上表现出了远超经典计算机的处理能力.随着量子计算任务复杂度的提高,如何分配量子计算资源,实现多方协作的量子计算,将成为量子计算领域待解决的一个重要问题.本文在一次性量子计算的基础上,提出了基于冗余图态的多人协作量子计算方案.不 ...
    本站小编 Free考研考试 2021-12-29
  • 基于计算全息和<i>θ</i>调制的彩色图像加密方法
    摘要:提出了一种基于计算全息和θ调制的彩色图像光学加密新方法.该方法利用彩色三基色原理和计算全息编码技术,首先将彩色图像的红、绿、蓝三基色分量进行随机相位调制和菲涅耳衍射变换,然后经过θ调制后进行图像叠加并编码为计算全息图,即加密过程是将一幅彩色图像加密为一幅实值的二元计算全息图,得到单幅密文.解密 ...
    本站小编 Free考研考试 2021-12-29
  • <sup>7</sup>Li<sub>2</sub><sup>(0, ±1)</sup>分子体系基态振-转能级的全电子计算
    摘要:采用单参考与多参考耦合簇理论结合相关一致高斯基组计算研究了7Li2(0,±1)分子体系的电子基态的势能曲线,计算考虑了体系所有电子的关联效应与相对论效应,拟合得到了体系的光谱常数,并获得了电子基态的振动-转动能级信息.计算得到的中性与阳离子体系的光谱常数与实验值符合得很好;对于阴离子体系,平衡 ...
    本站小编 Free考研考试 2021-12-29
  • 三氨基三硝基苯基高聚物粘结炸药热力学性质的理论计算研究
    摘要:高聚物粘结炸药(PBX)的热力学性质是用于炸药结构响应、安全性评估、数值模拟分析等的重要参数.由于PBX结构的多尺度特性,完全采取实验方法精细表征这些参数存在巨大的挑战.本文运用第一性原理和分子动力学计算的方法,系统研究了三氨基三硝基苯(TATB)基高聚物粘结炸药的热力学参数和界面热传导性质. ...
    本站小编 Free考研考试 2021-12-29
  • 基于数据库进行乏燃料鉴别的多元统计分析研究
    摘要:近年来,随着国际核军控形势的变化,包含防扩散、防核恐及核安保的多边国际军控合作越来越受到重视.核取证技术作为防扩散、防核恐及核安保的一项核心技术,在对涉核非法活动的威慑、阻止以及响应方面具有重要作用,值得深入研究.目前针对核取证技术的研究较多,主要集中于材料的表征和数据的解读.其中解读作为核取 ...
    本站小编 Free考研考试 2021-12-29
  • 曲波域统计量自适应阈值探地雷达数据去噪技术
    摘要:非线性、非平稳探地雷达数据常掺杂各种复杂噪声源,其对精确提取弱反射波信号、识别绕射波双曲线同相轴特征具有严重影响,忽略噪声影响给探地雷达探测数据全波形偏移成像及后续解译造成较大误差.采用传统阈值函数的曲波变换去噪需要根据数据噪声水平人为确定合理阈值控制系数.对此,本文开展自适应阈值函数的曲波变 ...
    本站小编 Free考研考试 2021-12-29
  • 钫原子磁偶极超精细结构常数及其同位素的磁偶极矩的理论计算
    摘要:应用基于B样条基组的相对论耦合簇理论方法,计算了212Fr原子的nS(n=7—12),nP(n=7—12)和nD(n=6—11)态的磁偶极超精细结构常数.与精确实验值的比较说明这套理论方法能精确计算出磁偶极超精细结构常数,其中7P态的磁偶极超精细常数的理论值与实验值之间的差异小于1%.在忽略场 ...
    本站小编 Free考研考试 2021-12-29
  • 基于共心球透镜的多尺度广域高分辨率计算成像系统设计
    摘要:针对实时广域高分辨率成像需求,充分利用具有对称结构的多层共心球透镜视场大且各轴外视场成像效果一致性好的特点,设计基于共心球透镜的多尺度广域高分辨率计算成像系统.该系统基于计算成像原理,通过构建像差优化函数获得光学系统设计参数,结合球形分布的次级相机阵列进行全局性优化,提高系统性能的同时有效简化 ...
    本站小编 Free考研考试 2021-12-29
  • 惯性约束聚变装置中靶面光场特性的统计表征方法
    摘要:在激光驱动的惯性约束聚变装置中,常采用多种束匀滑手段对焦斑的时空特性进行调控.光传输链路中涉及的光学元件众多、传输变换复杂,往往导致光传输模型复杂,且在运用衍射光学方法分析焦斑形态和特征时面临大量的数据处理和计算,致使出现计算量大、计算效率低等问题,亟需寻求快速而简便的新方法来描述焦斑的统计特 ...
    本站小编 Free考研考试 2021-12-29
  • Ce-La-Th合金高压相变的第一性原理计算
    摘要:采用第一性原理计算对Ce0.8La0.1Th0.1在高压下fcc-bct的结构相变、弹性性质及热力学性质进行了研究讨论.通过对计算结果的分析,发现了合金在压力下的相变规律,压强升高到31.6GPa附近时fcc相开始向bct相转变,到34.9GPa时bct相趋于稳定.对弹性模量的计算结果从另一角 ...
    本站小编 Free考研考试 2021-12-29