改进准静态格林函数节块法

提出了一种新的反应堆动力学方法－－改进准静态格林函数节块法，在处理反应堆三维瞬间问题时，这种方法在时间方面采用了改进准静态方法，通过对中子通量的因式分解，放大了三维空间计算的时间步长，减少了进行空间计算的次数。     

三维圆柱几何格林函数节块法

发展了中子扩散计算三维圆柱几何格林函数节块法。通过横向积分将中子扩散方程化为３个互相耦合的一维偏通量方程。对于径向偏通量方程，将径向扩散微分算符分解为平板几何的扩散微分算符和１个附加项之和，将附加项移到方程右端作为１个附加源项，这样，３个方程都化为平板几何一维方程的形式。再借助平板几何第二类边界条件格林函数，建立偏通量积分方程。方程推导中，对圆柱形曲面几何的线积分和横向积分均需对相应的广义线元作积分，对于修正源项，通过分部积分方法将偏通量导数项转化为对格林函数的求导，通过源迭代法求解方程。基准计算表明，该计算精度高、速度快、成为三维圆柱几何堆芯计算的有效方法。     

节块格林函数法的微扰计算

在反应堆物理设计和分析时,经常要进行微扰计算,以快速分析堆芯中子截面扰动下反应性的变化量。本文从微扰计算的普遍公式出发,给出了节块格林函数法（ＮＧＦＭ）下微扰计算的具体公式。通过对比验算,验证了ＮＧＦＭ下的微扰公式,并且证明微扰计算需要的是节块法的数学共轭解而不是物理共轭解。     

节块格林函数法的数学共轭方程

在反应堆物理计算中，扰动计算需要中子本征值方程的共轭方程。对于节块法、需要的是数学共轭方程、而不是物理共轭方程，但数学共轭方程求解比较困难。本文研究了节块格林函数法的数学共轭方程。并提出严格的求迭代公式，以和前向方程计算相当的计算量得到了其数学其轭方程的解。     

用于轻水堆扩散计算的格林函数节块展开法

提出一种用于轻水堆多维扩散计算的新方法─—格林函数节块展开法（ＧＮＥＭ），它是对节块展开法（ＮＥＭ）和节快格林函数方法（ＮＧＦＭ）的改进。参照ＮＥＭ把节块内偏中于通量用高阶多项式展开，利用格林函数求出节块表面偏中子通量与中子流的耦合关系，并用净中子流作为中间变量以简化计算。根据提出的理论编制了二维计算程序ＧＮＥＭ，对ＬＷＲ基准问题的数值计算表明，ＧＮＥＭ的计算精度与ＮＧＦＭ相当，而计算速度较ＮＧＦＭ和ＮＥＭ分别提高了１倍与１／３倍。应用ＧＮＥＭ取代ＰＳＵＩ－ＬＥＯＰＡＲＤ／ＡＤＭＡＲＣ中的扩散程序，对三里岛核电站（ＴＭＩ－１）的第１、第６循环进行了燃耗计算，其计算速度提高了３倍多，且具有更高的精度。     

两种圆柱几何格林函数节块法的比较

用格要函数节块法求解圆柱几何中子扩散方程时,可用平板几何格林函数代替复杂的圆柱几何格林函数。本文将基于平板格林函数的圆柱几何格林函数节块法（RNGFM）与基于贝塞尔函数的和何格林函数节块法（CNGFM）进行对比分析。分析结果表明：RNGFM的优点是公式简单,而CNGFM的迭代过程所需的计算时间则较少。     

格林函数节块法堆芯燃料管理程序NGCFM2D

本文简述了格林函数节块法堆芯燃料管理程序 NGCFM2D 的数学模型和程序结构,并分析和比较了用该程序计算的部分结果。     

基于保角变换的六角形中子扩散半解析节块法及数值结果

基于横向积分技术的半解析节块法,在求解方形中子扩散方程中取得了非常良好的效果。基于横向积分技术的六角形中子扩散问题会遇到奇异问题,而通过采用保角变换理论使该问题得到很好的处理,并能够转化为与方形问题很类似的形式。本文根据保角变换方法,将六角形中子扩散方程的求解整合到快速精确的半解析节块方法求解核心中;通过基准问题的数值检验,证明这种方式是相当有效和准确的。     

基于格林函数节块法的物理与热工-水力耦合方法研究

开发了三维物理与热工-水力耦合的PWR堆芯瞬态分析程序NGFMN-K/COBRA-Ⅳ/COBRA-Ⅳ(NCC)。少群时空中子动力学计算采用格林函数节块法程序NGFMN-K,隐式耦合子通道程序COBRA-Ⅳ实现瞬态计算。采用P10H8B功率重构方法给出热组件栅元功率分布,耦合另一个COBRA-Ⅳ程序模块,进行热组件子通道分析得到安全参数。对NEACRP-L-335 C1弹棒基准问题的计算表明,NCC程序的计算结果与参考结果符合很好,说明程序计算正确,可用于评估事故结果。     

格林函数法在快速计算热应力中的应用研究

通过对格林函数在热应力分析中的应用研究,得到一种结合有限元分析和数值计算进行热应力快速计算的格林函数方法。以核电厂主管道安全注入接管在升温和降温瞬态下的热应力分析为算例,验证该方法的合理性和准确性。同时对比分析格林函数法和有限元分析法的优缺点。     

基于格林函数的应力简化计算方法及试验验证研究

以格林函数为基础,开发了高温高压条件下管道应力影响函数,在实验室条件下选取管道壁面上的典型点,通过高温应变测量对该方法进行了验证。结果表明,基于格林函数的应力简化计算方法的计算结果与试验结果吻合较好,该方法可应用于核电厂关键设备和管道的强度及疲劳寿命的快速评价,亦可用于疲劳监测系统开发。     

Variant of Green's Function Nodal Method for Neutron Diffusion

A variant of the Green's function nodal method derived from the boundary integral form of the multigroup neutron diffusion equation in rectangular geometry is presented. As usual in the nodal methods, the multi-dimensional diffusion equation is integrated in the transverse direction. The resulting 1D diffusion equation is solved following the Boundary Element technique in one dimension. In this way a weighted residual method is obtained, with a Green's function for weighting, but with different boundary conditions than normally applied in the Green's function nodal methods. Mathematical formulation of the method is given and the iteration procedure is described. A computer program BINDIF has been designed, based on the new method. Its capabilities include the solution of the multigroup neutron diffusion equation of ID, 2D and 3D rectangular lattices. The BINDIF program has been checked against other methods used for global reactor calculations on benchmark problems, representative of realistic power reactor cores. The results indicate that the method is attractive to design highly efficient algorithms for a large mainframe, a personal computer or a parallel processor.     

应用不连续因子修正的六角形解析节块方法

在六角形解析节块方法中引入不连续因子,对该方法进行改进研究,并研制了分析程序HANDF和HANDF-C。应用该程序对UO2和MOX燃料组件组成的基准题进行了计算,结果表明,改进后的方法能够有效地提高反应堆堆芯(特别是非均匀性较强的堆芯)功率分布和有效增殖因子的计算精度。     

离散纵标六角形节块法及CMFD加速研究

离散纵标节块法是一种求解六角形中子输运方程的有效方法。本文基于六角形横向积分离散纵标方程,解析得到横向积分通量出射通量与入射通量的关系,并根据类似于扩散方法的六角形输运节块中子平衡方程形式,得到了一种离散纵标六角形节块法数值迭代策略。由于离散纵标法收敛速度较慢,本文根据粗网有限差分（CMFD）技术导出离散纵标六角形CMFD加速方法。数值计算结果表明,该CMFD加速技术能取得约16倍的加速效果。     

解析函数展开节块法的数值稳定性分析

解析函数展开节块法是求解六角形几何堆芯中子注量率的常用方法,其核心思想是利用严格满足中子扩散方程的解析基函数,直接把节块内的各群中子注量率近似展开。研究发现,此方法在求解一些特殊问题时,出现计算发散的情况。针对此类问题,提出了两种解决方法：截断近似方法和泰勒展开方法。通过对改造后的VVER440基准题验证表明：这两种方法在保证计算精度的基础上很好地解决了解析函数展开节块法的数值不稳定性问题。     

基于非线性迭代的节块方法

使用粗网有限差分方法 ,同时利用解析方法通过非线性迭代对差分近似进行修正 ,解两群扩散方程。数值计算结果证明 ,这种非线性方法具有较高的精度 ,并且较传统方法计算效率成倍提高     

格林函数时间卷积法在输运求解中的应用

利用瞬时源的格林函数卷积,可求解具有一定时间分布源粒子的输运问题。首先从输运方程的一般形式推导得到了格林函数时间卷积法的公式,并结合蒙特卡罗数值模拟算例,验证了格林函数时间卷积法在定态含时输运问题中的适用性,为相关问题的求解提供了新的思路。进一步数值分析表明,对于某些实际模型的蒙特卡罗求解,在同样计算样本条件下,格林函数卷积法还能实现较高的计算精度。     

快堆三维六角形节块法输运计算研究

对快堆堆芯的三维六角形节块输运理论计算方法进行了研究。在NAS程序基础上,使用节块输运理论,开发了三维六角形节块法输运程序NAS-T。表面通量采用双球谐DP1近似或DP3近似展开,在节块内采用一维离散纵坐标（S_N）方法求解。并根据节块法的特点,在迭代计算时采用了响应矩阵方法,节块内的通量分布用二阶多项式表示,简化了节块内部计算过程,大幅节省了计算时间。对中国实验快堆和BN-600的k_eff进行了计算,结果符合良好。