安 萍，姚 栋

(中国核动力研究设计院 核反应堆系统设计技术重点实验室，四川 成都 610041)

近二十多年来，求解多维中子扩散方程的粗网格计算方法得到广泛发展，其中矩形先进节块方法的发展和应用尤为成功。Lawrence等[1-2]成功发展了矩形节块格林函数方法(NGFM)，并通过大量例题验证了理论及程序的正确性和高效性。随着采用六角形几何的快堆的发展和俄罗斯VVER型压水堆的应用，六角形节块法的研究也日益得到重视，将先进矩形节块法推广到六角形节块有非常重要的意义。

Chao等[3]研制的六角形节块扩散程序ANC-H是目前公认的国际上精度最高的六角形节块程序。程序采用保角变换思想将六角形几何变换为矩形几何。本文采用该保角变换将六角形节块变换为矩形节块，对变换后的矩形节块扩散方程进行横向积分。并应用第二类边界条件的格林函数法，建立六角形节块各表面净中子流为未知量的空间耦合方程，采用剩余权重源迭代方法进行求解。研制三维多群六角形格林函数堆芯程序NACK，采用大量基准题对其进行验证，并与国际上相应程序的计算结果进行比较。

1 保角变换

保角变换是一种特殊映射，它保持任一点邻域内变换前后的伸缩性和旋转角不变，并保持角的定向。存在六角形到矩形的保角变换因子f[3]。图1示出变换前后对应点的关系，其中a和b为矩形长和宽。拉普拉斯算子在保角变换下形式不变，即：

(1)

图1 六角形节块到矩形节块的保角变换