基于任意三角形网格的中子扩散变分节块法

传统的基于矩形和六角形几何的堆芯计算程序已不适用于具有复杂几何的新型反应堆堆芯计算,本文开展了基于任意三角形网格的多群中子扩散变分节块方法研究。首先,采用ANSYS软件对计算区域进行三角形网格剖分,并利用坐标变换将任意三角形变换为正三角形;其次,采用Galerkin变分技术建立包含节块中子平衡方程的泛函,将三角形节块内变量利用正三角形内正交基函数进行展开;最后,利用变分原理,获得中子通量密度与节块边界上分中子流的响应关系,并基于传统的源迭代法对其进行求解。基于上述理论模型开发了程序TriVNM,并采用不同几何基准题进行了验证。结果表明,TriVNM计算的堆芯keff和归一化功率分布与参考解吻合较好,该计算方法适用于复杂几何堆芯扩散计算。