高效偶阶离散纵标法研究及初步验证

采用TAKEDA基准题验证高效偶阶离散纵标法(HEPSN)的临界计算功能,数值结果表明：有效增殖因数和各区域平均中子通量密度值与基准值吻合良好,堆芯区域的平均中子通量密度误差在0.7%以内,其他区域较扩散理论误差显著减小。相较于传统离散纵标法,HEPSN计算效率更高,需要的存储更小;相较于扩散理论,HEPSN的计算精度更高。因而,高效偶阶离散纵标法在大规模输运计算中具有很好的应用前景。     

基于JSNT程序的压水堆屏蔽计算

反应堆屏蔽计算是评估核电站安全性能的基础,是指导电站设计、运行的重要手段之一。JSNT程序是中物院高性能数值模拟软件中心研制的大规模并行离散纵标输运程序,具有较高的计算精度和计算效率。利用JSNT程序对某压水堆进行了建模计算,给出了中子通量密度的分布结果,并与实验测量值进行了对比。对比结果显示：无论是采用S8计算还是S16计算,计算结果都能满足工程要求;相比S8而言,采用S16计算可以显著提高计算精度,能够将某些测点处的相对误差降至1%以内。     

基于JSNT程序的压水堆屏蔽计算

反应堆屏蔽计算是评估核电站安全性能的基础,是指导电站设计、运行的重要手段之一。JSNT程序是中物院高性能数值模拟软件中心研制的大规模并行离散纵标输运程序,具有较高的计算精度和计算效率。利用JSNT程序对某压水堆进行了建模计算,给出了中子通量密度的分布结果,并与实验测量值进行了对比。对比结果显示:无论是采用S8计算还是S16计算,计算结果都能满足工程要求;相比S8而言,采用S16计算可以显著提高计算精度,能够将某些测点处的相对误差降至1%以内。     

蒙特卡罗粒子输运问题的全局降方差方法

精确获得中子通量的时间-空间-能量精细分布对于一类非定常粒子输运蒙特卡罗模拟至关重要.以零方差理论为基础,采取化整为零的策略,把蒙特卡罗方法与离散纵标方法相耦合,用近似重要函数指导蒙特卡罗模拟,给出一种实现蒙特卡罗全局降方差的计算方案.数值算例表明,该计算方案获得了全局降方差的效果.     

粒子输运方程的确定论计算方法

针对实际应用中辐射和中子输运数值模拟,讨论球一维和柱二维几何粒子输运方程确定论计算方法的研究现状,包括离散纵标、球谐函数、迭代加速、并行计算等方法.重点讨论输运计算方法所取得的若干研究进展,包括离散纵标求积组、自适应时间离散格式、本征值迭代求解方法、简化球谐函数方法、修正的子网格隅角平衡方法、灰体综合加速方法、迭代初值选取方法、输运与扩散耦合方法、基于预估校正的并行格式等.简要介绍了相关输运计算程序的研制情况,并分析输运计算方法存在的难点,提出待开展研究的内容.     

全堆芯Pin by Pin中子扩散节块法CMFD加速

采用节块方法中有效的半解析方法求解中子扩散方程，在径向采用pin by pin的计算方式且中子通量基于两阶展开，在轴向采用基于4阶展开中子通量的节块法，并且进行全堆粗网有限差分（CMFD）加速，以达到计算精度与计算效率的平衡.通过IAEA 2D/3D基准问题和自定义3D pin by pin参考问题的数值计算，验证所提出的方法能够达到预期的精度，同时CMFD能够取得很好的加速效果.     

节块内嵌离散纵标(SN)方法的算法研究及程序开发

介绍节块内嵌离散纵标(S_N)方法求解三维堆芯中子输运/扩散方程的算法框架.在基于扩散理论的三维粗网节块展开方法(NEM)的算法体系中,用基于输运理论的径向二维细网节块离散纵标方法(NDOM)的内迭代过程,替代节块展开方法(NEM)内迭代的径向求解过程.该算法充分考虑了核电厂反应堆堆芯的三维结构特点,另一方面,也充分利用了已经成熟的三维粗网节块展开方法(NEM)和二维离散纵标方法(S_N)的研究成果,同时有效避免了利用离散纵标方法(S_N)求解三维中子输运方程所面临的计算内存和计算时间的瓶颈问题.编制开发二维多群节块离散纵标方法(NDOM)模块程序NS_NM和三维多群节块展开方法(NEM)模块程序MGNEM,并以此为基础编制开发节块内嵌S_N方法的模块程序HANWIND;其中,NS_NM为HANWIND求解两维问题的功能模块.针对OECD/NEA-2D C5G7MOX基准问题以及两环路核电厂三维堆芯的数值验算结果表明,节块内嵌S_N方法的算法开发及程序编制有效、切实可行.     

瞬发中子密度衰减法计算中子代时间

采用蒙特卡罗程序MCNP计算了西安脉冲堆中子代时间。使用MCNP程序模拟了反应堆瞬发中子通量密度衰减,基于忽略缓发中子项的点堆动力学方程计算出中子代时间。在微次临界下,研究了次临界度、源的分布、计数区域等对西安脉冲堆中子代时间计算结果的影响。计算分析表明：采用瞬发中子密度衰减法计算中子代时间时,微次临界度、源分布、计数区域等对计算结果影响都很小;误差产生的主要原因是忽略缓发中子项的点堆动力学方程并不能较好地反应瞬发中子通量密度的衰减规律。     

瞬发中子密度衰减法计算中子代时间

采用蒙特卡罗程序MCNP计算了西安脉冲堆中子代时间。使用MCNP程序模拟了反应堆瞬发中子通量密度衰减,基于忽略缓发中子项的点堆动力学方程计算出中子代时间。在微次临界下,研究了次临界度、源的分布、计数区域等对西安脉冲堆中子代时间计算结果的影响。计算分析表明:采用瞬发中子密度衰减法计算中子代时间时,微次临界度、源分布、计数区域等对计算结果影响都很小;误差产生的主要原因是忽略缓发中子项的点堆动力学方程并不能较好地反应瞬发中子通量密度的衰减规律。     

基于曲四面体单元剖分的CN E-H TDFEM方法

从采用Crank-Nicolson差分格式的基于麦克斯韦旋度方程的E-H时域有限元方法出发,将展开电场和磁场的叠层矢量基函数与曲四面体单元相结合.对金属球谐振腔及介质填充的圆柱形谐振腔的数值模拟表明：相较于规则四面体单元,在剖分单元数目相同的情况下,曲四面体单元离散表面弯曲结构可以获得更高的计算精度.同时,与0.5阶基函数结合曲四面体单元相比,1.0阶基函数与曲四面体单元结合可以用更少的单元数及未知量数目来获得更高的计算精度.     

三维中子输运分流再平衡加速及其并行算法

为加速基于离散纵标方法（SN）的中子输运计算,针对群内迭代进行研究并实现了分流再平衡（PCR）算法及基于区域分解的并行,并在加速方程求解时进一步耦合了算术多重网格算法。并行PCR算法已集成至JSNT软件,并针对临界和外源计算典型算例验证了算法的正确性和有效性。74亿自由度的国内某商业压水堆压力容器屏蔽计算中,JSNT迭代次数减少至原来的1/3以下,总体计算时间减少至原来的1/2以下,1024核并行效率达42.8%。     

三维中子输运分流再平衡加速及其并行算法

为加速基于离散纵标方法（SN）的中子输运计算,针对群内迭代进行研究并实现了分流再平衡（PCR）算法及基于区域分解的并行,并在加速方程求解时进一步耦合了算术多重网格算法。并行PCR算法已集成至JSNT软件,并针对临界和外源计算典型算例验证了算法的正确性和有效性。74亿自由度的国内某商业压水堆压力容器屏蔽计算中,JSNT迭代次数减少至原来的1/3以下,总体计算时间减少至原来的1/2以下,1024核并行效率达42.8%。     

Analysis of relaxation phenomena in heavy-ion collisions

The transfer of mass, as well as the dissipation of relative kinetic energy and relative angular momentum in deeply inelastic heavy-ion collisions are studied as functions of time. Based on the determination of a parametrized deflection function from the experimental angular distributions, a classical model for the calculation of the mean interaction time as a function of initial relative angular momentum is presented. The method allows to include also those processes which correspond to long interaction times. The model is applied to determine mass transport coefficients from experimental mass (or element) distributions. The resulting mass drift and diffusion coefficients are accurate within less than 30% and compare well with the systematics obtained from the microscopic transport theory. The energy loss as function of interaction time is consistent with the following picture: Fast dissipation of the radial part of the kinetic energy accompanied by the loss of angular momentum with a larger relaxation time.     

Electron transport and small angle collisions

Anisotropic Coulomb scattering of electrons necessitates high order Legendre cross-section expansions or restricted angular quadratures in numerical multigroup, discrete ordinates transport calculations. We investigate the monoenergetic (Spencer-Lewis) collision operator in a small scattering angle approximation and give comparative results for electrons. Equivalence between the small-angle collision operator and diffusion-in term of the Fokker-Planck expansion is also exhibited. Such procedure reduces the forward scattering singularity and the number of moments required in applications. An overview of multigroup discrete-ordinates methodology for electron scattering is briefly detailed. To span diffusive to transport phenomena, a diffusion synthetic technique is also employed. Agreement with exact analytic and Monte Carlo predictions is good.     

Evaluation of quadrature schemes for the discrete ordinates method

After reviewing the presently available quadrature schemes for the discrete ordinates method, the accuracy of different schemes is analyzed and evaluated. It is shown from a comprehensive error analysis that the moment conditions have to satisfied not only for the principal coordinates directions, as it is mostly carried out, but for any arbitrary test direction. Among the schemes with approximately 50 discrete ordinates the DCT-020-2468 quadrature was found to give the best accuracy. The highest accuracy among all schemes is achieved by the LC-11 quadrature which requires 96 discrete ordinates. This scheme is rarely used up to date and deserves more attention for high accuracy predictions.     

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

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

Assessment of the re-ordered wide band model for non-grey radiative transfer calculations in 3D enclosures

The objective of the present study is to evaluate variations of the re-ordered wide band model for non-grey radiative transfer calculations in 3D enclosures using the discrete ordinates method. First, the performance of various angular and spatial discretisation schemes of the discrete ordinates method is investigated. Then, several formulations, averaging procedures, and scaling methods of the re-ordered wide band model are tested, and the results are validated against those of a statistical narrow band model. The grey gases formulation using three optimised absorption coefficient is found to be the most efficient method.     

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

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