基于平面通量展开的三维Pin by Pin输运SP3计算

基于Pin by Pin输运SP₃计算是下一代反应堆物理计算方法中一个重要的研究方向。节块法一般采用高阶通量和中子源展开以适应较大的节块尺寸，因此对于Pin by Pin计算来说效率偏低。由于反应堆堆芯不均匀性更多发生在径向，本文提出一种径向基于二维平面通量展开结合轴向常规节块法的综合方法进行三维Pin by Pin输运SP₃计算，其三个方向的节块耦合迭代均采用基于中子净流耦合的格式。通过IAEA 2D/3D基准问题和典型输运基准问题计算，验证该计算方法对于堆芯扩散计算和输运计算的正确性和可行性。     

节块内嵌离散纵标(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方法的算法开发及程序编制有效、切实可行.     

基于非线性迭代的压水堆瞬态计算程序开发

采用目前工业上成熟的非线性迭代计算策略,基于两群粗网有限差分方法和多群UNM节块方法,开发了针对压水堆工况的三维瞬态扩散计算程序。UNM方法采用解析基函数作为基函数,通过方程变换解决了解析节块法在临界节块计算不稳定的问题,提高了计算精度。热工计算采用单通道模型和燃料棒一维导热模型,相比目前堆芯在用的经验关系式方法,该模型可以更加准确地计算燃料棒温度分布。采用基于横向积分方程的三节块方法,可以有效减轻控制棒尖齿效应对瞬态计算的影响。为测试程序性能,采用NEACRP等基准算例对程序进行了校验。数值结果表明,开发的程序计算结果正确,适用于压水堆堆芯瞬态过程的模拟。     

三维圆柱几何格林函数节块法中子扩散计算

发展了中子扩散计算三维圆柱几何格林函数节块法。首先通过横向积分将中子扩散方程化为三个互相耦合的一维偏通量方程。对于径向偏通量方程，将径向扩散微分算符分解为平板几何的扩散微分算符和一个修正项之和，将修正项移到方程右端作为修正源项。这样，三个方程都化为平板几何的一维方程莆式，再借助平板几何第二类边界条件格林函数，对主几相应体源作积分，建立偏通量积分方程，对于修正源项，通过分部积分方法将偏通量导数项转化     

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

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

改进的节块展开方法求解圆柱几何对流扩散方程

为高效求解球床高温气冷堆物理-热工耦合问题,发展改进节块展开法求解圆柱几何下的对流扩散方程.针对圆柱几何和对流扩散方程的特殊性,采用三阶多项式和指数函数作为r向横向积分方程的展开函数,在节块展开法的框架下高效求解对流扩散方程.数值验证表明,改进的节块展开方法具有固有的迎风特性,在使用粗网节块时依然能保持稳定性和较高的计算精度.     

基于不连续因子的沸水堆芯三维瞬态数值模拟

基于不连续因子校正的粗网格有限差分法是实现堆芯瞬态三维数值模拟的高效方法之一, 粗网节块的界面不连续因子与边界反照率的计算方法决定了实时数值模拟过程中的精度。在计算不连续因子的过程中, 省去了细网节块计算与粗网均匀化过程, 直接在粗网格划分情况下, 基于节块展开法和非线性迭代策略, 推导了粗网格界面不连续因子比率与边界反照率的计算公式, 并编制了相应的计算程序。沸水堆典型算例的三维瞬态模拟证实该方法可在空间域和时间域两方面, 使静态、瞬态精度均达到与先进节块法相等同的程度, 并且计算效率优于先进节块法, 为核电站全范围模拟机三维堆芯的实时仿真模型开发提供了一种切实可行的选择。     

中国实验快堆燃料组件精细功率分布计算

在热工中,需要中子学计算给出燃料组件内各元件棒功率的相对分布。利用蒙特卡罗程序对中国实验快堆（CEFR）燃料组件内元件棒功率分布进行了理论计算分析,并保证计算结果相对统计误差小于0.8%。使用另一个基于六角形节块扩散理论的钠冷快堆中子学设计软件NAS计算得到的结果对蒙特卡罗程序计算结果进行了对比计算。结果表明,蒙特卡罗程序与NAS计算得到的元件棒相对功率分布结果的最大相对偏差小于3%。使用蒙特卡罗程序对CEFR燃料组件内精细功率分布的计算是可靠的,可用于设计计算当中。     

高温气冷堆堆芯中子时空动力学模拟计算

研究开发三维圆柱几何堆芯多群中子时空动力学改进准静态方法模拟计算程序.对给定的模块式高温气冷堆模型进行模拟计算.初始状态下,计算结果与中子扩散程序CITATION吻合很好.动态情况下,模拟堆芯反应性、堆内各能群中子平均注量率和堆芯相对功率等物理量随时间变化,计算结果与理论分析一致.     

基于大规模并行计算的三维多群中子扩散方程有限差分方法

三维多群中子扩散方程的精确、高效求解是核动力堆芯设计及燃料管理的基础。应用有限差分方法求解该方程具有简便、精确、成熟的优点;然而,该方法的计算量和存储量均较大,极大地限制了它的计算规模和应用范围。本文基于大规模并行计算,研究三维多群中子扩散方程有限差分方法:采用中心有限差分格式离散中子扩散方程;基于MPI并行编程模型,采用空间区域分解的方式实现大规模并行计算;采用多群多区域耦合PGMRES算法进行并行加速。在集群服务器上开发了ParaFiDi程序,并采用IAEA3D,PHWR等多个基准题对该程序进行验证。数值结果表明,ParaFiDi程序具有较高的计算精度和计算效率。     

Fast 3D T(2)-weighted MRI with Hadamard encoding in the slice select direction

A fast method to obtain 3-dimensional (3D) magnetic resonance imaging with long repetition times is presented. It can be used to obtain fast 3D MRI with for example T(2) or diffusion weighted imaging. The method uses a 3D multiple thin slab sequence with radio frequency encoding, preferably Hadamard encoding, in the slice select direction. The point-spread function of the Hadamard-encoded slices is close to ideal even at low encoding numbers. This allows the acquisition of 3D data volumes with tolerable image quality up to four times faster than is possible using Fourier phase encoding. The scope of the method includes both longitudinal and transverse encoding. Longitudinal encoding provides a better point spread function than transverse encoding, at the expense of having to discard one slice per slab. The method is demonstrated experimentally for 4th order longitudinal Hadamard encoding to obtain 3D T(2)-weighted images.     

Simulation study on the boron-coated GEM neutron beam monitor

A Monte Carlo method of the 3D GEM simulation is introduced.The physical process of the neutron detection is described with the Geant4 code and the Garfield code.The results of the optimized electric-fields,the emitted ions spectrum,the electrons transverse diffusion and the signal width are presented.Moreover,the preliminary result with a CF-252 neutron source is reported.These are useful in designing detector structures and to provide an optimized option for the development of the boron-coated GEM neutron beam monitor.     

对称性及多群中子扩散方程数值解

在多群中子扩散方程解析解的基础上,利用方程及求解域的对称性建立了新的数值求解中子扩散方程的理论模型.该模型显著的优点是适用于各种对称区域(二维、三维区域)尤其是非正方形区域中子扩散方程的求解,它彻底避免了常规节块法应用于非正方形几何时所出现的奇异性问题,且所得的解在求解域内任意点上均满足扩散方程.以二、三维六角形几何为例建立了数学模型,并用基准问题校核了模型的正确性. 关键词： 中子扩散方程 对称群 数值解 解析