一种先进的用于精细化扩散计算的均匀化方法.docVIP

外文文献复印件 一种先进的用于精细化扩散计算的均匀化方法 Yonatan Berman （Department of Physics, Nuclear Research Center-Negev, 84190 Beer-Sheva, Israel Department of Physics, Ben-Gurion University of the Negev, 84105 Beer-Sheva, Israel） 摘要：本文提出了一种用于精细化扩散计算的栅元均匀化方法。这种方法通过引入与空间变量弱相关的扩散系数来增加均匀化问题中的自由度。相对于传统的精细均匀化方法，如FVW或SPH，这种方法能使平均反应率、界面流和本征值得到更好的守恒。本文还提出了一种采用该方法的迭代算法：半均匀化迭代法（ISH），利用该方法对一组一维单群的基准题进行了计算并将计算结果与其他均匀化方法进行了比较。结果表明ISH在保持反应率和界面流的守恒上优于FVW和SPH法，在保持本征值的守恒上则优于FVW法。 关键词：精细均匀化；栅元计算 1 引言 在大多数反应堆堆芯计算中，栅元计算程序被用来得出均匀化截面和非连续因子（多数情况是对节块法而言）。栅元计算通常采用至少50个能群，对反应截面和其他均匀化参数并群处理得到一个少群问题。整个堆芯计算通常采用2-4个能群。而且，大多数栅元程序采用PN或SN近似求解中子输运方程，在堆芯计算中则通常采用扩散近似。这种栅元计算和堆芯计算之间的差异在预估中子通量密度（或功率）分布、反应率和有效增殖系数（或本征值）时可能会产生显著的误差。从这个角度来说，应该注意在标准的栅元计算程序中，均匀化截面并不能保证均匀化前后反应率的守恒（当没有采用非连续因子DFs时）。为了使上述差异引起的误差最小化，在堆芯计算中，均匀化过程应将上述差异考虑在内。 节块法在理论和数值计算中都被证实能够保证非均匀介质中不同参数的守恒，被广泛地用于反应堆分析和设计中。由于近几年计算机性能的发展，对详细的堆芯模型进行精确计算成为可能。三维细网的精细化计算便是其中的一种，该方法的均匀化只在单个栅元内进行，避免了对整个组件的空间均匀化。 最近关于精细均匀化技术的研究主要是在探讨一维平板基准题或二维方形栅元。这些技术适用于大多数现代压水堆，但并不适用于其它可能的几何形状。这些技术通常能很好地保证反应率的守恒，但并不能精确的保证界面流和增殖系数的守恒，特别是在具有很强吸收性的栅元中，比如控制棒栅元。通过改进这些技术使其计算更准确或减少网格数，从而能应用于三维堆芯的精细化计算中。 本文研究的目的是对标准均匀化方法提出一种改进，即半均匀化迭代法（ISH）。这种方法得到的均匀化反应截面和非均匀离散并与空间相关的扩散系数，在栅元均匀化问题中能维持几个重要物理量的守恒。 在第2部分写出了栅元均匀化公式并讨论了求解的困难，在第3部分提出了半均匀化迭代法（ISH），第4部分应用ISH法计算了一维基准题并与传统的体积-通量权重法（FVW）和SPH法进行了比较，最后在第5部分给出了结论。 2.均匀化理论 2.1.均匀化过程的数学公式 在均匀化过程中很难同时保持所有重要物理量在均匀化前后守恒，为了证明这一点，考虑中子输运方程的多群近似形式并假设我们能够得到该方程的解： （1） 式中：，，， ，，为反应堆本征值，为g群中子角通量，G是能群数，反应截面的符号采用常规形式。 在解出方程（1）后，在栅元均匀化过程中需保持守恒的物理量可以很容易计算得出。三个最重要的物理量分别是栅元平均反应率，平均界面流和本征值（或增殖系数）。与式（1）相似，可得到均匀化方程： （2） 式中：，，， ，，和是栅元均匀化后所得的与r无关的反应截面。通常，在整个区域中这些截面是非均匀的，即与r相关。 下面推导均匀化参数和非均匀参数之间的关系式，根据均匀化前后重要物理量的守恒： （3） （4） 式中，、分别是均匀化区域的表面积和体积。 2.2.均匀化中的数学悖论 由于输运方程中均匀化截面和均匀化通量间的强烈耦合，要计算出均匀化参数以使式（3）和式（4）成立并不容易。另一方面，均匀化过程中共需求解个均匀化参数（表示总能群数，表示反应类型数加上扩散系数），需保持个物理量守恒（个反应率和个均匀化区域各表面的界面流，其中表示在笛卡尔坐标系下问题的维数，），而均匀化问题的自由度不足。这一现实表明在栅元均匀化中要准确保持重要物理量的守恒是不可能的。等效均匀化理论和广义等效均匀化理论通过引入非连续因子绕开了这一数学上的难题。例如，SPH法中加入了个自由度（SPH因子），能适用于一维或多维且对称的问题，但并不适用于一般的二维或三维问题。要使所有重要物理量准确的守恒，就需要引入更多的自由度。组件结构的对称性可以有效的降低反应