关于《快中子脉冲堆特性参数数值计算》的学术研究报告

本文旨在向同行研究者介绍邱东研究员于2000年9月发表在《原子能科学技术》（Atomic Energy Science and Technology）第34卷增刊上的一项原创性数值模拟研究。该研究聚焦于快中子脉冲反应堆（Fast Burst Reactor, FBR）的核心物理过程，通过创新的多物理场耦合方法，对脉冲爆发过程中的关键特性参数进行了高精度的理论计算，并成功指导了实验，显著提升了理论预测与实验数据的一致性。

一、 研究团队与发表信息 本项研究由来自中国工程物理研究院核物理与化学研究所（四川绵阳）的邱东副研究员独立完成。研究成果以题为《快中子脉冲堆特性参数数值计算》的论文形式，发表于2000年9月出版的《原子能科学技术》期刊。该期刊是中国核科学技术领域的重要学术刊物，此研究发表于其第34卷的增刊。

二、 学术背景与研究目的 本研究属于核反应堆物理与工程领域，具体涉及特种反应堆——快中子脉冲堆的瞬态动力学分析。快中子脉冲堆能够在极短时间内（微秒量级）释放出巨大的裂变能量脉冲，这种特性使其在科研和特定应用中具有独特价值。然而，脉冲爆发过程涉及中子学、热力学、结构力学等多个物理场的强烈耦合与瞬变，对其进行精确的理论预测极具挑战性。

在本文工作之前，国内外虽已开展相关理论研究，但计算结果与实验数据之间时常存在较大差异。一个被普遍忽视的关键因素是堆体外反照中子（reflected neutron） 的影响。传统理论通常认为缓发中子对快脉冲过程无关紧要，因而也未能细致考虑从堆体外部反射回来的中子（即反照中子）对脉冲特性（如脉冲宽度、峰值功率、总产额）的潜在影响。这些反照中子虽然返回堆芯的时间延迟很短，但其效应类似于具有特殊时间常数的“缓发”中子源，可能在快速瞬态过程中起到不可忽略的作用。

因此，本研究的主要目的包括：1) 建立一个更符合物理实际的理论模型，通过将反照中子作为“模拟缓发中子”纳入考虑，以改进对快脉冲堆瞬态行为的描述；2) 开发一套耦合中子输运、热弹性力学和反应堆动力学的数值计算程序，对脉冲爆发的主要特性参数进行精细化数值模拟；3) 将理论计算结果与XFBR脉冲堆的实验数据进行对比验证，并利用改进的理论指导完成更高产额的脉冲倍增实验。

三、 详细研究流程与方法 本研究本质上是理论建模与数值计算工作，其核心流程是构建物理模型、开发求解算法、进行数值模拟并与实验对比。研究并未涉及传统意义上的实验样本处理，但其“研究对象”可视为XFBR脉冲堆的物理系统及其在脉冲瞬态下的行为。整个工作流程包含以下几个关键步骤：

第一步骤：多物理场耦合理论模型的建立。 研究者针对快脉冲堆脉冲爆发这一特定物理过程，建立了三个相互耦合的核心方程： 1. 热弹性力学方程： 用于描述脉冲期间裂变能瞬间沉积导致堆芯材料受热、膨胀并产生应力的过程。由于脉冲持续时间极短，模型忽略了热传导效应。在球对称几何下，该方程描述了材料密度、径向应力与切向应力随时间和空间的变化关系，并设定了球心位移为零、外边界径向应力为零的边界条件。 2. 反应堆动力学方程： 采用空间平均近似，用于计算瞬态功率随时间的变化。本研究的创新点在于，在经典的点堆动力学方程中，除了常规的缓发中子先驱核组外，显式地加入了代表堆体外反照中子效应的附加项。这些项被处理为具有特殊时间常数的“模拟缓发中子”组（即公式中的 n 包括了真实缓发中子组和反照中子组）。模型中假设初始状态为缓发临界，并在t=0时刻引入一个阶跃反应性扰动，从而推导出功率满足的积分方程形式。 3. 中子输运方程： 采用多群、各向同性散射近似的球几何模型，用于计算任意时刻堆芯的中子通量分布和系统的有效中子倍增系数（keff）。在脉冲过程中，虽然材料密度和宏观截面因膨胀而变化，但研究者基于物理洞察，假设中子通量形状在短时间内变化不大，因此可以在每个时间步长内将时间视为参数，求解keff本征值问题来获得当前构型下的反应性。

这三个方程通过预加反应性这一实际操作参数作为共同输入和耦合纽带，形成了一个完整的闭合系统：动力学方程计算功率变化 -> 功率变化导致材料温升和膨胀（热弹性力学方程）-> 材料膨胀改变几何与密度，从而影响中子学性能（中子输运方程）-> 中子学性能变化（keff变化）反馈回动力学方程，影响下一步的功率计算。

第二步骤：数值求解算法的开发与实现。 针对上述耦合方程组，研究者开发了基于有限差分方法的数值求解方案，并编写了Fortran计算程序。具体方法如下： - 热弹性力学方程（式1）： 采用中心差分格式进行离散求解。为了处理XFBR脉冲堆多层球壳的实际结构，模型假设各层球壳之间均存在一个薄真空层，使得每层壳的内外表面在力学上都是自由的，简化了边界条件处理。空间网格步长的选择严格遵守数值稳定性条件（库朗条件）。 - 中子输运方程（式5）： 采用不连续坐标Sn差分格式进行求解。该格式具有绝对稳定的优点，空间步长划分主要依据精度需求。此方程求解可同时得到本征值keff、中子角通量/通量分布以及出射中子流等关键中子学参数。 - 反应堆动力学积分方程（式4）： 采用数值递推公式进行求解。首先根据缓发中子（含模拟反照中子）参数计算出一组常数 s_j 和 a_j，然后进行数值积分计算瞬态功率。

关键的耦合计算策略与性能优化： 考虑到三个方程计算代价差异巨大（中子输运计算最耗时），研究者设计了一套自适应多尺度时间推进算法来高效实现耦合。具体流程是：从一个初始状态开始，先求解动力学方程得到一个小时间步长后的功率；以此功率作为输入，求解热弹性力学方程得到材料的温升和位移；然后将新的材料构型输入中子输运方程，求解得到新的keff（即反应性）和通量分布。完成一个“大”时间步的耦合循环。 为了大幅缩短计算时间，算法根据快脉冲的物理机制进行了优化：每进行一次昂贵的中子输运计算后，可以进行多次力学计算；而每进行一次力学计算后，又可以多次推进动力学方程（功率计算）。这种基于物理过程时间尺度的解耦策略，在保证精度的前提下显著提升了计算效率。

第三步骤：计算执行与结果分析。 研究者将上述方法和程序应用于XFBR脉冲堆（文中提及为类球型结构）的模拟。计算中，预加反应性（δkp） 被设定为唯一的输入控制变量。程序系统地计算了不同预加反应性下，脉冲的峰值裂变率、脉冲半宽度（FWHM）、脉冲总产额、材料应力以及边界振动等关键特性参数随时间的变化关系。

四、 主要研究结果及其逻辑关联 研究获得了与实验高度吻合的系统性计算结果，揭示了多个重要物理现象：

1. 脉冲特性参数与预加反应性的关系： 图1展示了峰值裂变率、脉冲半宽度和脉冲产额随预加反应性变化的计算曲线，并与XFBR脉冲堆的实验数据进行了对比。结果显示，在理论模型中考虑反照中子效应后，计算值与实验值符合得非常好。这直接验证了将反照中子作为模拟缓发中子处理的合理性与必要性。理论计算成功预测了这些关键脉冲参数随反应性增加的变化趋势，为指导实验提供了定量依据。

2. 材料力学行为的转变阈值： 计算结果（图2及相关论述）揭示了一个重要的物理阈值：当预加反应性较小时（δkp < 6.2×10^-4），材料中的应力主要由热应力主导，且其大小随反应性变化较为平缓。然而，当预加反应性超过该临界值（δkp > 6.2×10^-4）时，动态效应开始显现并占据主导。此时，应力主要来源于动应力，它强烈依赖于温度升高的速率，而不再主要取决于堆芯的壳层结构。总应力（热应力与动应力之和）随预加反应性的增加而急剧上升。这一发现对评估脉冲堆的结构安全性和确定安全运行区间具有重要意义。

3. 预加反应性与实际反应性的“定向差”： 研究通过理论分析指出，在实际脉冲爆发过程中，测得的“实际反应性”（由中子波形前沿周期反推得到）与操作输入的“预加反应性”之间总是存在一个固定的差值，即“定向差”。本文指出，这一定向差正是由动力学效应和堆体周围常设反射体（导致反照中子存在） 等因素共同造成的不可避免的结果。本文的理论模型通过引入反照中子组，在物理上解释了这一差异的来源，从而在理论上弥合了“预加”与“实测”反应性之间的鸿沟。

结果之间的逻辑关系： 首先，改进的模型（引入反照中子）成功预测了宏观脉冲参数（图1），这证明了模型在整体动力学上的正确性。基于这一正确模型，进一步分析微观的力学行为，发现了应力机制转变的阈值（结果2），这为理解堆芯材料在强瞬态载荷下的响应提供了新见解。同时，模型对“定向差”的解释（结果3）从物理机理上统一了操作参数与物理参数，使理论框架更加自洽和完整。所有这些结果共同指向一个结论：反照中子效应是精确模拟快脉冲堆瞬态行为不可忽略的关键因素。

五、 研究结论与价值 本研究的主要结论是：通过将中子输运方程、热弹性力学方程和反应堆动力学方程进行耦合，并创新性地将堆体外反照中子效应作为额外的缓发中子项引入动力学方程，发展出了一套功能完整的快脉冲堆理论设计与计算程序。该程序的计算结果与XFBR脉冲堆的实验结果符合良好，显著改善了以往理论计算与实验数据差异较大的问题。

其科学价值在于：深化了对快中子脉冲堆瞬态复杂多物理场耦合过程的理解，特别是明确了反照中子在快速瞬态动力学中的作用，为快脉冲堆物理提供了更精确的理论模型。 其应用价值在于：所开发的数值计算程序具有较好的通用性和适应性，可用于快脉冲堆的设计优化、安全分析及运行预测。文中特别提到，在理论计算结果的指导下，成功完成了XFBR脉冲堆的脉冲产额倍增实验，这直接证明了该研究对实际工程实验具有重要的指导意义。

六、 研究亮点 1. 物理模型的创新： 首次在快脉冲堆瞬态理论模型中，系统性地考虑了堆体外反照中子效应，并将其处理为“模拟缓发中子”，这是提升计算精度的关键创新点。 2. 多物理场强耦合数值方法的实现： 成功实现了中子学-热工-力学三者在瞬态过程中的紧密耦合数值求解，并开发了高效的自适应时间推进策略，平衡了计算精度与效率。 3. 理论与实验的紧密结合： 研究不仅追求理论上的自洽与改进，更以实验数据为验证基准和最终目标，利用改进的理论成功指导了更高水平的实验（脉冲产额倍增），形成了“理论改进-实验验证-指导新实验”的完整闭环。 4. 揭示了重要的物理阈值： 发现了材料应力机制从热应力主导向动应力主导转变的预加反应性阈值，这对脉冲堆的安全运行范围界定具有实际参考价值。

七、 其他有价值的内容 本研究还展示了针对特定工程问题开展深入基础研究的典范。研究者对XFBR脉冲堆“类球型结构”的具体情况进行了针对性建模（如球对称假设、多层球壳间设真空层等），使得理论研究能够紧密贴合实际装置，增强了结果的可信度和实用性。此外，全文对耦合算法的描述细致，尤其是对不同方程采用不同差分格式以及为提升效率而设计的计算流程控制策略，对从事反应堆多物理场耦合计算的研究人员具有很好的方法论参考价值。