不同燃耗计算模型对商用压水堆乏燃料组件核素成分的影响分析

燃耗计算精度对提高乏燃料贮存效率有着重要影响,在应用燃耗信用制时,燃耗计算得到的核素成分偏差决定了乏燃料贮存的临界安全裕量。不同燃耗计算模型所得到的核素成分偏差各不相同,为提高燃耗计算精度,提出了一种装载不同燃料富集度的多组件燃耗计算模型,并使用不同燃耗计算模型分别对TMI-1反应堆NJ07OG组件中的6个样本进行了计算、对比和分析。结果表明,相比其他模型,考虑不同燃料富集度的多组件模型得到的²³⁵U、²³⁸U和²³⁹Pu等核素平均相对偏差更接近于零且6个样本的相对偏差分布更为平均。     

组件计算程序PANDA研发及初步验证

PANDA是上海核工程研究设计院研发的压水堆组件计算程序。该程序采用基于特征线方法（MOC）的一步化计算流程,即在不引入能群压缩和栅元均匀化的情况下直接进行组件层面的两维非均匀输运计算。多群数据库采用基于ENDF/B-VI制作的70群结构中子数据库,基于ENDF/B-VII的新版数据库也正在开发中。共振自屏计算采用了空间相关丹可夫方法（SDDM）,既具备燃料芯块分区计算的能力,又保留了传统Stammler方法的计算效率。多群非均匀输运计算采用二维模块化MOC方法,并辅以双重粗网有限差分(CMFD)加速技术,具有良好的计算精度和效率。对传统线性子链解析（TTA）方法以及多种矩阵指数方法进行了研究,选取了适合PANDA程序燃耗链的燃耗方程求解技术。基于以上基本模型开发了PANDA程序,并从程序模块、总体集成和核设计程序系统确认等三个层面,初步验证了PANDA程序的计算性能,表明了PANDA程序的工程设计计算能力。     

基于线性核素链的燃耗算法与蒙特卡罗程序耦合计算

根据线性核素链原理和反应堆燃耗的特点,建立相应核数据库并采用回溯算法生成自适应核素链,完成核素的遍历和计算,形成多群点燃耗计算程序,能够独立进行燃耗计算并具有完整的输出结果形式。同时将所开发的燃耗计算程序与蒙特卡罗程序进行耦合,完成接口模块的设计,形成可用于研究堆和核电站的燃耗-输运耦合计算工具。将所开发的输运-耦合计算程序应用于中国实验快堆首炉堆芯燃耗的计算,将计算结果与现有设计数据进行比较。经过初步分析表明:新的燃耗计算程序能够精确计算锕系核素含量,包括易裂变核素和生成量很小的次锕系核素。对某些锕系核素如Pu241的计算结果还存在较大偏差,这需要对计算结果进一步分析,来确认偏差是来自计算过程还是相关的截面数据。整个燃耗-输运耦合计算系统对裂变产物的处理和反应性变化的计算也与现有的设计数据符合良好。     

先进中子学栅格程序KYLIN-2的开发与初步验证

针对先进核反应堆中结构复杂的燃料组件,中国核动力研究设计院开发了先进中子学栅格中子学程序KYLIN-2,该程序采用子群方法进行共振处理,采用特征线方法（MOC）进行复杂几何中子输运计算,并利用采用广义粗网格有限差分加速方法（GCMFD）来加速中子输运求解流程,采用基于改进预估校正临界-燃耗迭代方法（PPC）的切比雪夫方法求解复杂燃耗链,同时,为了方便用户使用,开发形成了支持复杂组件图形化建模工具和后处理显示工具。通过初步的数值检验,针对典型压水堆燃料组件,KYLIN-2具有较高的计算精度,满足工程使用需求。     

行波堆深燃耗的计算特点

行波堆属于新概念堆型,卸料燃耗深度可达400 GWd/tHM,是现有快堆的3~4倍、压水堆的6~8倍,较高的卸料燃耗深度对堆芯物理分析工具计算正确性提出挑战。基于此,以KYLIN-1程序为基础,从能谱、裂变产物核素重要性、燃耗计算误差累积等方面探究行波堆深燃耗计算特点。对典型行波堆六角形组件分析结果表明:低富集度铀组件寿期初、末能谱差别较大,采用单一权重谱制备的多群截面库用于其燃耗计算时,无限增殖系数偏差较大;为保证行波堆深燃耗计算的正确性,燃耗链应包含重要的70种裂变产物核素;行波堆深燃耗计算时,由于燃耗步增多累积的误差较小,无限增殖系数偏差每燃耗步约为0.001%。     

基于RMC的微观截面参数化北大核心CSCD

为了进行堆芯计算,需要通过组件计算提前构建少群截面参数库。传统确定论的组件截面参数化方法针对宏观截面进行截面参数化,但这种方式不仅需要考虑多种物理状态参数,而且需要考虑历史效应对截面的影响。提出了基于核素微观截面的蒙卡程序参数化方法,该方法可以消除燃耗历史的影响,且考虑的物理状态仅为燃耗深度以及材料温度。利用蒙卡程序产生组件截面参数库耦合堆芯程序进行堆芯计算,首先用蒙卡程序同时统计对应状态点下的核素密度以及核素少群微观截面,再利用核素微观截面进而获得宏观截面进行后续堆芯计算。为了验证方法正确性,构造了一个自定义的压水堆模型,计算结果与连续能量蒙卡计算结果符合良好。     

两种大规模燃耗链求解算法对比分析

为严格追踪裂变反应堆中核素成分随燃耗的变化，基于燃耗矩阵法求解燃耗方程，分别采用自主编写的Chebyshev有理近似方法(CRAM)程序和广泛应用的ORIGEN2程序进行大规模燃耗链的点燃耗计算，并对两种算法的相关参数进行对比分析。结果表明:在计算精度方面，CRAM与ORIGEN2程序获得的重要核素的核密度较为一致，个别核素相对误差较大；在计算效率方面，单步燃耗计算ORIGEN2略胜一筹，但CRAM耗时也非常短；在步长稳定性方面，CRAM具有显著优势，而ORIGEN2的统计结果受步长变化的影响较大。     

MCNP-FISPACT耦合燃耗计算程序开发与测试

反应堆堆芯的燃耗计算关系到堆芯的燃料管理,并直接影响堆芯的经济性评估,因此如何快速且准确地对堆芯进行燃耗计算一直是反应堆物理设计的研究重点之一。随着反应堆的发展,其几何结构和物理特性日渐复杂,现有的一维、二维耦合燃耗程序因其在几何处理上的限制,难以满足先进反应堆精细设计分析的要求。为对复杂反应堆堆芯的燃耗情况进行计算,结合粒子输运程序MCNP 处理复杂几何和燃耗程序FISPACT处理核素全面的特点,开发了接口程序耦合MCNP 和FISPACT来进行燃耗计算,并对耦合程序进行了计算验证。采用了IAEA 基准校核例题和CFETR中国聚变工程实验堆例题进行程序验证,经计算得出的有效倍增因子随燃耗的变化曲线和TBR等数据与标准例题的结果符合良好,其误差在可接受范围内。     

基于MCMG-Ⅱ和STEP1.0的输运-燃耗耦合系统架构设计

实现燃耗过程的精细化计算需要开发新型输运-燃耗耦合计算系统,系统中采用了三维多群中子输运蒙特卡罗程序MCMG-Ⅱ和基于回溯算法的多群点燃耗计算程序STEP1.0。燃耗与输运程序之间的耦合需要考虑能谱、总源强、反应率统计、初始核密度及其变化、裂变产物核素等信息。耦合系统采用直译式脚本语言Python来实现,可以充分利用其强大的文本处理功能和直译式的特点,以准确、便捷地将这些数据在程序之间进行自动转换,并根据总功率、辐照时间、步长选择等参数自动完成输运-燃耗耦合的分步计算和结果的图形化处理,从而实现整个系统的无缝连接。在耦合系统的软件架构设计中采用分层与封装策略,降低了开发难度,提高系统的可移植性和可扩展性。耦合系统能够更为精细地考虑几何形状、能谱的变化、辐照时间等因素,自动进行精细化模拟计算。     

MCADS程序的开发和ADS基准题计算

由于加速器驱动次临界堆存在外中子源,堆芯结构复杂,中子注量的各向异性严重,所以相关燃耗计算在次临界系统设计中起着重要作用。为实现次临界系统的燃耗计算,结合粒子输运程序MCNP处理复杂几何和燃耗程序LITAC处理核素全面的特点,开发了接口程序MCADS耦合MCNP和LITAC。然后选取IAEA-ADS基准题对耦合程序进行了验证计算。结果表明,燃耗、外源强度、空泡效应、初始功率分布等方面的计算结果和其他国家的计算结果相比有很好的一致性,证实了MCADS在次临界模式计算中的可靠性。     

基于VERA基准题的cosRMC程序验证

为了验证反应堆物理软件和方法的计算能力,美国CASL (Consortium for Advanced Simulation of LWRs) 项目提出了VERA (Virtual Environment for Reactor Application) 堆芯物理基准题。该基准题以Watts Bar初始堆芯为模型,涵盖从二维单栅元到三维全堆芯的燃耗及换料的十个基准问题。针对VERA基准题模型,利用COSINE软件包中的反应堆蒙特卡罗分析程序cosRMC进行临界计算,得到了有效增殖因子、组件功率分布、控制棒微积分价值和反应性系数等结果。通过与基准题中提供的KENO结果对比,两种蒙特卡罗程序的计算结果吻合良好。这表明cosRMC程序具有从组件到堆芯的计算能力,其临界计算精度基本与KENO程序相当。     

数学建模中的高维数据挖掘技术优化研究

高维数据挖掘由于特征空间占用开销较大,挖掘的复杂度较高,挖掘精度不高,为了提高对高维数据挖掘的准确性能,提出一种基于相空间重构和K-L变换特征压缩的高维数据挖掘数学建模方法。采用集成学习技术,对高维数据信息流进行相空间重构处理,考虑类间的数据不平衡性,求得高维数据的关联维特征参量,根据数据的链距离进行稀疏性融合,计算高维数据流模型的最大Lyapunove指数谱,根据谱分析方法实现数据聚类,对聚类后的数据采用K-L特征压缩方法进行降维处理,降低数据挖掘的内存及计算开销。仿真结果表明,采用该方法进行高维数据挖掘,数据挖掘的准确概率较高,占用内存消耗较少,计算开销较小。     

潜望式激光通信瞄准机构误差计算

本文提出了一种机构指向误差计算方法。首先,采用反射镜单位法向量和光束单位方向向量分别表示机构反射镜的姿态和激光束的传播方向,接着采用坐标变换和矩阵传递求出反射镜在总体坐标系中的单位法向量,通过向量对称计算出激光束经过两次镜面反射后的方向向量,最后将机构各参数的理论值和实际值代入上述算法,分别求出激光束经过两次镜面反射后的理论方向向量和实际方向向量,两向量夹角即为机构的指向误差。在得到指向误差和部件误差的关系式后,机构设计时即可通过机构指向误差反计算求出部件的误差值。结果表明:机构物理样机测试指向误差值为0.005 7°,将机构部件参数测试误差值代入误差正计算算式,得出机构的指向误差值为0.006 1°,以该指向误差通过反计算得出的参数平均误差值为0.002 25°,与测试结果相近,可以为机构的精度设计提供参考依据。     

网络计算机模型下海量大数据存储系统设计

对网络计算机模型下海量大数据进行安全稳定的存储,可以提高网络计算机的使用价值,增加其使用周期。但目前的海量大数据存储方法在存储过程中,无法做到对其进行灵活高效的存储,存在大数据存储分布密度较低,存储开销较大等问题。为此,以网络计算机模型体系结构为基础,提出了一种基于ARM的海量大数据存储系统设计方法。该设计方法先利用ARM芯片对网络计算机模型下海量大数据存储系统进行硬件构造,将网络海量大数据中的可利用与不可利用数据进行分类处理,采用VISA结构根据数据分类结果对大数据存储系统软件部分进行设计,依据大数据调度模型和存储相似度特征对大数据存储的时间,质量等方面进行计算,利用计算结果对大数据传输的阈值以及分布密度进行记录,最后根据循环分段的计算方式进行冗余大数据特性压缩,并对海量大数据的常规数据和冗余数据进行存储。实验仿真证明,所提方法提高了海量大数据存储的兼容性,增强了大数据存储的精确性和灵活性。     

基于随机抽样法的多群核数据不确定性影响分析

基于随机抽样方法,研究多群核数据不确定性对反应堆物理计算的影响。首先利用SCALE软件包中核数据协方差矩阵和自主开发的随机抽样模块SAMP,得到多群微观截面等核数据的抽样值,之后分别使用SCALE/TRITON和PARCS程序进行组件计算及堆芯稳态计算,最后通过统计分析得到组件和堆芯计算结果的不确定度。以Almaraz压水堆核电厂装载的燃料组件和首循环堆芯为对象,研究了不同燃耗下有效增殖因子、动力学参数、核素浓度和双群均匀化宏观截面等组件计算结果,以及堆芯功率分布等堆芯计算结果的不确定度。分析结果表明:组件计算结果不确定度多随燃耗变化,快群宏观截面不确定度总体高于热群;堆芯计算结果受核数据不确定性影响显著,其中稳态径向功率分布的最大不确定度为1.9%左右。     

近红外光谱用于植物样品中水溶性氯离子含量的测定

基于离散小波变换(DWT)和最小二乘支持向量回归(LSSVR)方法,建立了近红外光谱测定植物样品中水溶性氯离子的回归校正模型。以烟草样品中水溶性氯离子含量的测定为研究对象,首先采用DWT对近红外光谱进行数据压缩和背景扣除,再用LSSVR建立氯离子的校正模型。结果表明,与偏最小二乘回归(PLSR)和传统的LSSVR方法相比,作者所建模型的预测准确性具有一定优势。     

基于PCA-BP神经网络的光谱消光法颗粒粒径反演算法研究

光谱消光法广泛应用于颗粒粒径测量领域,在利用光谱消光法对颗粒粒径进行反演的过程中,由于颗粒的消光系数存在理论复杂、计算繁琐、收敛速度慢以及求解不稳定等问题,很大程度上影响了整个反演过程的快速性和准确性。且在众多波长的消光数据中,存在较多重复冗余的信息,也很大程度上增加了反演算法的时间。针对光谱消光法粒径反演算法计算繁琐、反演效率低的问题,提出了基于主成分分析（PCA）和BP神经网络的光谱消光颗粒粒径分析方法。基于Mie散射理论对不同粒径、不同波长下的光谱消光值进行了仿真计算,通过对光谱消光数据集的主成分分析及各个波长综合载荷系数的计算,实现了最优特征波长的选取,利用降维后的光谱消光数据训练了PCA-BP神经网络模型,并利用该网络模型计算了粒径颗粒分布。通过仿真计算,比较了PCA-BP神经网络模型与传统的BP神经网络模型的预测精度,并分析了波长数目对两种神经网络模型预测结果的影响。针对训练得到的PCA-BP神经网络模型开展光谱消光法粒径参数反演算法的验证实验,搭建了光谱消光法颗粒粒径参数测量实验系统,测量了粒径范围在0.5～9.7 μm内的6种不同粒径参数的聚苯乙烯标准颗粒。仿真和实验结果表明：基于主成分分析方法可确定各个波长向量之间的相关性,利用综合载荷系数选取最优特征波长对应的消光值对整体的光谱数据具有较好的代表性,可实现光谱数据的降维。相比传统的BP神经网络模型,基于PCA-BP神经网络模型的颗粒粒径分布的分析方法预测精度更高,对于较分散颗粒系的分布参数的预测有更加明显的优势。而且,被选取的波长数较少时,PCA-BP神经网络模型依然有较高的预测精度。利用训练好的PCA-BP神经网络模型对颗粒粒径参数进行实验验证,预测结果可瞬时输出,颗粒粒径分布误差在5%以内,验证了该算法的可行性。     

基于投影和重采样技术的目标旋转角度测量方法

精密电子组装中贴装元件的检测定位通常采用针对元件形状特征的分析方法。但是随着贴装元件形状的多样性,使得基于元件形状特征分析的检测定位变得困难。提出了一种新颖的基于投影的目标旋转角度测量方法。首先预先精确建立元件目标投影长度与旋转角度之间的关系,然后检测运动目标的投影长度,最后通过查找建立的投影长度与旋转角度之间的关系,即可确定角度旋转量。为了提高该方法的定位精度,采用图像重采样技术来提高角度检测精度。实验结果表明了该方法的有效性,旋转角度测量具有快速、高效和高精度的特点。