Krylov子间方法及先进重水堆在线通量拟构方法初步研究.pdfVIP

上海交通大学硕士论文 Krylov 子空间方法及先进重水堆在线通量拟构方法初步研究 摘 要 在反应堆物理领域，随着高阶谐波的广泛应用，对高阶谐波求解方 法的研究也越来越受到关注。近年来Krylov 子空间方法在矩阵特征值问 题和线性方程组的求解中得到了很大的研究和发展。在基于Krylov 子空 间投影的方法中，Rice 大学的D.C.Sorensen 等人提出的一种隐式重启动 的Arnoldi 方法（IRAM ）可以方便快速地同时求取矩阵的前几阶特征值 和特征向量，为此，本文部分的研究工作是对IRAM 方法在中子扩散方 程高阶谐波求解中的应用加以研究。研究结果表明：该方法在求解高阶 谐波时具有和常规的源迭代方法相同的计算精度；其求解速度与所选取 的子空间维数有关，选择适当的子空间维数可获得最快的计算速度；对 于Krylov 子空间方法，求解速度基本不受高阶谐波阶次的影响，而常规 源迭代法的求解速度随着高阶谐波阶次的升高而大大降低。因此对于需 要求解高阶次谐波的问题，Krylov 子空间方法的计算速度远比源迭代法 快；同时由于IRAM 计算过程中不断的正交化过程使得IRAM 方法能够 求解具有重根的本征值问题。本文还初步探讨了IRAM 方法在非线性本 征值问题求解中的应用。 本文第二部分工作是对先进重水堆（ACR ）的在线拟构方法进行研 究。对于重水堆来说，利用堆芯内探测器的测量信号来在线拟构全堆三 维中子通量密度空间分布在反应堆的运行管理中有着非常重要的作用。 由于本文研究的在线拟构方法，需在线地进行全堆三维扩散计算，因此 是否能够在线拟构取决于全堆芯扩散计算的时间。然而与压水堆相比， ACR 的几何尺寸要大得多，且必须进行全堆计算，因此通常的堆芯扩散 计算程序无法满足计算速度上的要求，因此开发一种快速堆芯扩散计算 方法是非常必要且迫切的。为此，本文探讨了多种方法来提高ACR 反应 堆全堆扩散计算的速度，包括非线性迭代法、全堆粗细网格交替扫描方 I 上海交通大学硕士论文 法以及粗细网格交替迭代过程中用反照率替代反射层加速的方法。通过 研究发现，反照率受堆芯状态的影响较小，因此本文最终采用反照率给 定条件下、全堆粗细网格交替迭代的方法，它能够满足ACR 反应堆在线 通量拟构对计算速度的要求。 最后，本文建立了一种将堆内探测器信号作为扩散计算时的内边界 条件、通过三维扩散计算在线拟构中子通量分布的方法，并通过数值验 证表明，该方法能够利用探测器的测量值，较有效地矫正中子通量分布 的计算值。 关键词： Krylov子空间，高阶谐波，非线性迭代，扩散计算加速，反照 率，通量拟构 II 上海交通大学硕士论文 Study on Krylov Subspace Method and Flux Mapping for Advanced CANDU Reactor ABSTRACT Calculation method of harmonics is being widely studied due to its importance in the reactor physics area. Recently, Krylov subspace method is rapidly developed for solving the large-scale linear equation and eigen-value equation. In this paper, Implicitly Restarted Arnoldi Method (IRAM), one of the Krylov subspace method, is applied to solving the higher order harmonics. The numerical results demonstrate that while its accuracy is comparable with the conventional