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

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

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

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

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

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

电磁线圈炮计算程序的收敛性与实验验证

根据电磁线圈炮中推力线圈与电枢之间的感应耦合方程和电枢动态响应方程,考虑发射过程中推力线圈和电枢的欧姆加热,开发了电磁线圈炮的计算程序,给出了程序的算法和流程。分析了程序的收敛性,认为最大时间步长和电枢最大网格均依赖于驱动脉冲波形的上升时间。通过模拟结果与单级线圈炮实验测量结果的比较,验证了程序的有效性。模拟的线圈电流波形与实验波形吻合较好,模拟的出口速度比测量速度大约6.5%。程序收敛性和有效性的验证表明该程序可用于电磁线圈炮系统的初步设计。     

快堆中子学计算程序NECP-SARAX的改进与验证

NECP-SARAX是西安交通大学NECP团队开发的用于快中子反应堆的中子学程序系统。为准确处理快中子反应堆中中等质量核素散射共振以及非弹性散射导致的复杂的中子慢化效应,SARAX程序最初采用连续能量的蒙特卡罗方法计算中子能谱从而获得堆芯计算使用的有效多群截面。由于蒙特卡罗程序计算效率低,且在低能量段统计偏差较大,提出采用基于点截面的超细群方法计算中子慢化能谱,避免了蒙特卡罗方法产生参数时存在的缺陷。堆芯计算采用多群中子输运,通过优化简化几何建模,改进了程序的实用性。采用多种微扰方法计算堆芯各种反应性系数,提出了基于中子输运微扰理论的虚拟密度方法以计算堆内组件变形导致的反应性变化。在进行堆芯瞬态计算时,采用了点堆和改进准静态两种方法,可用于一般快堆和快谱ADS的典型事故分析。OECD发布的一系列快堆基准题测试表明,SARAX程序在快堆计算中具有良好的精度,达到了与国外著名快堆程序相当的水平。有效增殖因子与连续能量的蒙卡计算结果相比偏差在300 pcm以内。同时,由于引入了虚拟密度理论和三维时空动力学模型,程序功能更加完善,可以更好地满足快堆工程设计的需求。     

快堆中子学计算程序NECP-SARAX的改进与验证

NECP-SARAX是西安交通大学NECP团队开发的用于快中子反应堆的中子学程序系统。为准确处理快中子反应堆中中等质量核素散射共振以及非弹性散射导致的复杂的中子慢化效应,SARAX程序最初采用连续能量的蒙特卡罗方法计算中子能谱从而获得堆芯计算使用的有效多群截面。由于蒙特卡罗程序计算效率低,且在低能量段统计偏差较大,提出采用基于点截面的超细群方法计算中子慢化能谱,避免了蒙特卡罗方法产生参数时存在的缺陷。堆芯计算采用多群中子输运,通过优化简化几何建模,改进了程序的实用性。采用多种微扰方法计算堆芯各种反应性系数,提出了基于中子输运微扰理论的虚拟密度方法以计算堆内组件变形导致的反应性变化。在进行堆芯瞬态计算时,采用了点堆和改进准静态两种方法,可用于一般快堆和快谱ADS的典型事故分析。OECD发布的一系列快堆基准题测试表明,SARAX程序在快堆计算中具有良好的精度,达到了与国外著名快堆程序相当的水平。有效增殖因子与连续能量的蒙卡计算结果相比偏差在300 pcm以内。同时,由于引入了虚拟密度理论和三维时空动力学模型,程序功能更加完善,可以更好地满足快堆工程设计的需求。     

MOC/SN耦合三维中子输运程序KYCORE开发与初步验证

KYCORE程序是中国核动力研究设计院开发的径向MOC（特征线方法）与轴向SN耦合三维中子输运程序。KYCORE将二维MOC与一维SN通过角通量实现高精度耦合,并通过粗网有限差分实现快速收敛,是目前可工程化应用于三维中子输运计算中精度最高的方法之一。介绍了2D/1D计算与加速理论,并通过与蒙特卡罗程序的计算对比,数值验证了KYCORE三维中子计算的准确性与高效性。     

MOC/SN耦合三维中子输运程序KYCORE开发与初步验证

KYCORE程序是中国核动力研究设计院开发的径向MOC（特征线方法）与轴向SN耦合三维中子输运程序。KYCORE将二维MOC与一维SN通过角通量实现高精度耦合,并通过粗网有限差分实现快速收敛,是目前可工程化应用于三维中子输运计算中精度最高的方法之一。介绍了2D/1D计算与加速理论,并通过与蒙特卡罗程序的计算对比,数值验证了KYCORE三维中子计算的准确性与高效性。     

钍铀循环关键核素的多群数据库加工与初步验证

钍基先进核能系统的发展对核数据提出了新的需求。钍铀循环关键核素的多群数据库是堆芯物理设计与分析计算的基础。因此,基于权威核数据处理程序NJOY对ENDF/B-VII.1评价核数据库进行处理得到WIMS格式多群截面数据,利用专用更新处理程序WILLIE得到WIMS格式数据库,并利用多群堆芯计算程序WIMSD-5B对数据库进行了一系列临界基准检验,结果表明:基于ENDF/B-VII.1加工的WIMS库与WIMS库更新计划（WLUP）的网站上发布的最新版本WIMS库的计算结果基本一致,并且,在16个钍铀循环基准题检验中,新加工的WIMS库计算平均误差要比WIMSD5B程序的自带WIMS库的计算结果小0.225 3%,精度更高,可靠性更好。     

钍铀循环关键核素的多群数据库加工与初步验证

钍基先进核能系统的发展对核数据提出了新的需求。钍铀循环关键核素的多群数据库是堆芯物理设计与分析计算的基础。因此,基于权威核数据处理程序NJOY对ENDF/B-VII.1评价核数据库进行处理得到WIMS格式多群截面数据,利用专用更新处理程序WILLIE得到WIMS格式数据库,并利用多群堆芯计算程序WIMSD-5B对数据库进行了一系列临界基准检验,结果表明：基于ENDF/B-VII.1加工的WIMS库与WIMS库更新计划（WLUP）的网站上发布的最新版本WIMS库的计算结果基本一致,并且,在16个钍铀循环基准题检验中,新加工的WIMS库计算平均误差要比WIMSD5B程序的自带WIMS库的计算结果小0.225 3%,精度更高,可靠性更好。     

A pellet tracking system for the PANDA experiment

Frozen microspheres of hydrogen (pellets) will be one of the target types for the future hadron physics experiment PANDA at FAIR (GSI, Darmstadt, Germany) [1]. Pellets with a diameter of 25– μm are generated about 3 meters above the interaction region, to which they travel with a velocity around 80 m/s inside a narrow pipe. The interaction region is defined by the overlap of the pellet stream and the accelerator beam and has a size of a few millimeters. One would like to know the interaction point more precisely, to have better possibilities to reconstruct particle tracks and events e.g. in charmonium decay studies. One would also like to suppress background events that do not originate in a pellet, but e.g. may occur in rest gas, that is present in the beam pipe. A solution is provided by the presented pellet tracking system together with a target operation mode that provides one and only one pellet in the interaction region most of the time. The goal is to track individual pellets in order to know their position with a resolution of a few tenths of a millimeter at the time of an interaction. The system must also be highly efficient and provide tracking information for essentially all pellets that pass the interaction region. Presented results from the design studies show that the goals can be fulfilled by this solution.     

一维磁化套筒惯性聚变模拟程序的设计与校验

在调研跟踪国外磁化套筒惯性聚变(MagLIF)研究领域最新进展的基础上,以一维磁流体力学方程为基础,结合氘氚燃料能量交换方程,建立一维磁化套筒惯性聚变物理模型,并使用Fortran语言编写完成一维数值模拟程序MagLIF-1D;通过与圣地亚实验室LASNEX,HYDRA等程序计算结果对比展示,完成程序校验工作,讨论后认为程序计算结果存在差异的主要原因可能来自不同程序对于材料状态方程库的选择;通过计算,MagLIF-1D程序可以直接获得内爆速度、燃料压力、燃料密度、聚变产额等关键物理量,这为后续更好地开展磁化套筒惯性聚变实验设计提供了有力工具。     

一维磁化套筒惯性聚变模拟程序的设计与校验

在调研跟踪国外磁化套筒惯性聚变(MagLIF)研究领域最新进展的基础上,以一维磁流体力学方程为基础,结合氘氚燃料能量交换方程,建立一维磁化套筒惯性聚变物理模型,并使用Fortran语言编写完成一维数值模拟程序MagLIF-1D;通过与圣地亚实验室LASNEX,HYDRA等程序计算结果对比展示,完成程序校验工作,讨论后认为程序计算结果存在差异的主要原因可能来自不同程序对于材料状态方程库的选择;通过计算,MagLIF-1D程序可以直接获得内爆速度、燃料压力、燃料密度、聚变产额等关键物理量,这为后续更好地开展磁化套筒惯性聚变实验设计提供了有力工具。     

熊猫型保偏光纤定轴仿真研究

应用光波导理论分析了保偏光纤的POL定轴方法的基本原理.根据熊猫型保偏光纤的折射率分布模型,采用光束传播方法对影响熊猫型保偏光纤POL曲线分布的各种因素进行了仿真研究,包括不同观测平面、不同折射率分布(匹配型或非匹配)、是否加匹配液等因素.通过仿真结果与实验观测结果以及保偏光纤与单模光纤仿真结果比较,表明了理论分析与仿真方法的可行性.     

BEAVRS基准模型的cosRMC精细建模与临界测试验证

BEAVRS基准题是来自麻省理工学院(MIT)的基于商用压水堆全堆芯精细建模的国际基准模型,具有详细的实测数据,是由多种富集度燃料及控制棒组件组成的大型压水堆分析验证基准题。使用我国自主研发的蒙卡程序cosRMC对BEAVRS基准题进行精细建模,主要计算了热零功率(HZP)状态下的临界本征值、全堆功率分布和控制棒价值,并与实测值以及国际知名蒙卡程序MCNP,OpenMC,MC21等结果进行对比。HZP状态下,cosRMC临界本征值结果与MCNP计算结果相差仅7.1 pcm,符合较好;不同控制棒组件插入情况下的临界本征值与理论值1.000的偏差小于0.74%,控制棒价值结果与实测值误差小于100 pcm,计算精度与同类软件相当;此外还对比分析了全堆功率分布与实测值结果的吻合程度及误差产生原因。初步验证了cosRMC程序对复杂堆芯精细建模计算的可行性和准确性,为程序以后的应用及完善打下基础。