关于U-Mo整体燃料肿胀行为预测模型研究的学术报告
本文旨在向中文科研界介绍并解读一篇发表于《Journal of Nuclear Materials》第544卷（2021年）的重要研究论文，题为“Swelling of U-Mo monolithic fuel: Developing a predictive swelling correlation under research reactor conditions”。该研究由爱达荷国家实验室（Idaho National Laboratory）的A.B. Robinson、W.J. Williams、W.A. Hanson、B.H. Rabin、N.J. Lybeck和M.K. Meyer团队完成，并于2020年12月5日在线发表。本研究属于类型a，即对一项原创性研究的详细报告。
一、 研究团队与发表信息 本研究的主要作者为A.B. Robinson及其在爱达荷国家实验室的同事。研究成果发表于材料科学与核工程领域的权威期刊《Journal of Nuclear Materials》2021年第544卷。该研究直接服务于美国高性能研究堆（High Performance Research Reactors, HPRRs）从高浓铀（Highly Enriched Uranium, HEU）向低浓铀（Low-Enriched Uranium, LEU）转换的国家级重大工程需求。
二、 学术背景与研究目标 本研究的主要科学领域是核燃料材料学，具体聚焦于研究堆用铀-钼（U-Mo）合金整体型板状燃料的辐照肿胀行为。理解并准确预测燃料在辐照下的尺寸稳定性，是确保反应堆安全运行和完成燃料资格认证（Qualification）的核心前提。
研究背景：为降低核扩散风险，全球致力于将研究堆燃料从高浓铀转换为低浓铀。U-10Mo（铀-10wt.%钼）整体型燃料因其高铀密度和良好的辐照稳定性，被选为美国HPRRs转换的燃料类型。在这种板状燃料设计中，U-10Mo燃料芯体表面覆盖一层锆（Zr）扩散屏障，再与6061铝合金（AA6061）包壳通过热等静压工艺结合。在辐照过程中，燃料的尺寸变化（肿胀）受多种机制共同影响，包括固体和气体裂变产物的积累导致的肿胀、辐照促进蠕变（Irradiation-Assisted Creep）、燃料相变以及燃料与包壳间相互作用层（Interaction-Layer）的形成。其中，裂变气体（如Xe和Kr）形成的气泡形态演变是导致高燃耗下肿胀行为复杂化的关键。在较低裂变密度（< ~2.5×10²¹ fissions/cm³）下，肿胀主要由固体裂变产物的体积膨胀主导，行为线性且稳定。当裂变密度超过约3.0×10²¹ fissions/cm³时，燃料发生辐照诱导的晶粒细化（或称再结晶），为裂变气体气泡在晶界上聚集和长大提供了场所，导致气体肿胀贡献加剧，肿胀速率可能发生变化。
研究目标：尽管已有多个肿胀模型被提出（如Kim-Hofman、Leenaers、Rest、Perez-Robinson和Wachs等人的模型），但这些模型在较高裂变密度区域（> ~5×10²¹ fissions/cm³）的预测与实验数据存在差异，且未能充分涵盖广泛的辐照条件。为了降低机理不确定性并为燃料性能分析提供更可靠的工程模型，本研究旨在利用大量实验数据，建立一个基于统计学的、能够预测U-Mo整体燃料在宽范围裂变密度下肿胀行为的经验关联式（Predictive Swelling Correlation），并提供相应的置信区间和预测区间。
三、 详细研究流程与方法 本研究的工作流程系统而严谨，主要包括数据收集、数据处理与“分箱”、模型拟合与比较、以及模型验证四个主要步骤。
1. 数据收集与实验对象： 研究团队收集了来自7个不同辐照实验的74个U-10Mo整体燃料试验板的辐照后数据。这些实验在先进试验堆（Advanced Test Reactor, ATR）中进行，覆盖了不同的燃料尺寸、功率密度和裂变密度范围。具体实验包括：AFIP-3（1个全尺寸板）、AFIP-4（6个子尺寸板）、RERTR-9（2个微型板）、RERTR-10（4个微型板）、RERTR-12（52个微型板，剔除了4个发生“鼓包”的板）、AFIP-6 MkII（高功率大型板）和AFIP-7（弯曲板元件组装）。总计分析了超过18,000个局部厚度测量数据点，构成了一个极其庞大的数据集。
2. 局部燃料肿胀计算流程： * 测量与修正：对每个燃料板进行辐照后轮廓测定（Profilometry），获得局部板厚度（t_p）。对于氧化层较厚的实验（AFIP-6MkII和AFIP-7），使用涡流技术测量氧化层厚度（t_ox），并从总厚度中扣除因氧化导致的体积增加部分（乘以系数1/1.975）。从制造记录中获取辐照前局部板厚度（t_p0）和燃料箔厚度（t_f），并已知锆层厚度（t_Zr）。 * 肿胀计算：采用公式 燃料肿胀 (%) = [ (t_p - t_ox/1.975) - t_p0 ] / (t_f - 2*t_Zr) * 100% 计算局部燃料肿胀。该公式假设所有厚度增加均源自燃料体积膨胀，但实际上也包含了因辐照蠕变导致的燃料重定位（Fuel Relocation）贡献，这是分析中的一个挑战。 * 裂变密度关联：使用蒙特卡洛中子粒子（MCNP）代码计算每个板的局部与平均比值（Local-to-Average Ratio, L2AR）裂变率梯度，再结合平均裂变密度，估算出每个厚度测量点对应的局部裂变密度（f_d）。将肿胀数据与裂变密度数据在统一的网格上进行对齐（通过插值或降采样）。
3. 数据处理与“分箱”方法： 由于原始数据点（局部肿胀 vs. 局部裂变密度）存在显著散射（主要源于燃料重定位、边缘约束效应以及测量和计算的不确定性），直接拟合难以揭示内在趋势。为此，研究采用了统计“分箱”（Binning）技术来减少散射。 * 方法比较：团队比较了四种分箱方法：按固定裂变密度步长（0.1, 0.25, 0.5 ×10²¹ fissions/cm³）分箱，以及按固定数据点数量（每100个点）分箱。 * 方法选择：通过线性回归分析发现，不同分箱方法得到的趋势模型其置信区间重叠，表明方法选择不影响统计趋势。最终选择“每100个数据点一箱”的方法，其优点是每个箱具有相同的统计权重，无需根据箱内数据点数量进行加权，且在高裂变密度区域（数据点较少）能提供更平衡的代表性。
4. 模型拟合与比较： 对分箱后的数据（包含所有实验，含AFIP-6MkII）尝试了四种函数进行拟合：线性、带可变拐点的分段线性、二次和三次函数。 * 分段线性拟合：虽能捕捉低裂变密度（固体肿胀主导）向高裂变密度（气体肿胀主导）的转变，但其拐点位置高度依赖于AFIP-6MkII实验在再结晶区域（~3-4.5×10²¹ fissions/cm³）的数据。若剔除AFIP-6MkII数据，拐点会从~3.3移至~4.5 ×10²¹ fissions/cm³，显示出不稳定性。且物理上再结晶是一个渐变过程，非尖锐转折，因此该模型被放弃。 * 多项式拟合：二次和三次函数的拟合优度（R²值）均超过96%，且统计上无显著区别。为简化模型，选择使用二次函数。值得注意的是，包含或排除AFIP-6MkII数据对二次拟合的稳定性影响较小。
5. 模型验证： 为了验证从局部厚度数据推导出的肿胀关联式的有效性，研究团队将其与RERTR-12实验微型板的板平均燃料肿胀数据进行了对比。板平均肿胀是通过辐照前后浸没法（Archimedes method）测量板体积变化直接计算得到的，该方法避免了局部厚度测量中燃料重定位等因素的干扰，能更纯粹地反映燃料的总体积膨胀。
四、 主要研究结果 1. 数据处理结果：原始局部肿胀数据与裂变密度的关系图显示了大量的数据散射，但在统计分箱后，数据清晰地呈现出一个随裂变密度增加而单调上升的趋势。分箱有效地平滑了由非体积膨胀机制（如蠕变重定位）引入的噪声，揭示了燃料肿胀的内在物理趋势。
2. 模型拟合结果： * 推荐肿胀关联式：对分箱数据（每100点一箱）进行强制过原点的二次拟合，得到最终的预测模型： 肿胀 (%) = 6.13×10⁻⁴³ * f_d² + 4.00×10⁻²¹ * f_d （方程17） 其中f_d为裂变密度（fissions/cm³）。该模型表明，在研究的裂变密度范围内，U-Mo整体燃料的肿胀行为可以用一个简单的二次多项式很好地描述。 * 置信区间与预测区间：研究提供了基于分箱数据的95%预测区间，适用于反应堆安全分析： * 下限：5.87×10⁻⁴³ * f_d² + 4.14×10⁻²¹ * f_d - 5.28 （方程18） * 上限：6.39×10⁻⁴³ * f_d² + 3.86×10⁻²¹ * f_d + 5.28 （方程19） 此外，也提供了基于原始数据的更宽的95%预测区间（方程20, 21），适用于实验人员比较单个测量值。
3. 模型验证结果：如图11所示，由局部数据推导出的推荐肿胀模型及其预测区间，与RERTR-12实验通过浸没法获得的板平均燃料肿胀数据高度吻合。这一关键验证表明，尽管局部厚度数据包含复杂因素，但通过大规模的统计分析和适当的数据处理方法，所建立的模型能够准确预测燃料的整体平均行为。模型预测在低浓铀燃料可达的最高裂变密度（约7.8×10²¹ fissions/cm³，对应100% LEU燃耗）以下，燃料肿胀是稳定且可预测的，未出现“突破性肿胀”（Breakaway Swelling）的迹象。
4. 与历史模型的比较：如图2所示，本研究推荐的二次模型与Kim-Hofman、Leenaers (U-7Mo弥散燃料)、Rest等人的模型在中等裂变密度下较为接近，但不同于Perez-Robinson和Wachs模型中预测的更高肿胀。本模型基于更广泛、更系统的数据集，并提供了统计意义上的不确定性量化，因此具有更高的可靠性和工程应用价值。
五、 研究结论与价值 本研究成功建立了一个基于大量实验数据的、用于预测U-Mo整体燃料在研究堆条件下肿胀行为的经验关联式。该模型采用二次多项式形式，简洁实用，并附有统计导出的预测区间，为燃料性能建模和反应堆安全分析提供了关键输入。
科学价值：研究证实了在宽裂变密度范围（直至近8×10²¹ fissions/cm³）内，U-Mo整体燃料的肿胀行为是连续、稳定且可预测的，未观察到突变的不稳定区域。这加深了对U-Mo燃料在辐照下（特别是经历再结晶和气体肿胀加剧阶段）宏观尺寸变化规律的理解。
工程与应用价值：该肿胀关联式及其预测区间可直接用于U-Mo整体燃料元件的设计、安全评估和运行限值制定，是美国HPRR向LEU转换燃料资格认证过程中的重要工具。它为反应堆设计师和监管机构提供了量化燃料行为及其不确定性的依据，有助于确保转换后研究堆的安全、可靠运行。
六、 研究亮点 1. 数据规模空前：研究分析了来自74个辐照燃料板、超过18,000个数据点，是迄今为止关于U-Mo整体燃料肿胀行为最全面、最系统的数据集分析。 2. 创新的数据处理方法：针对局部厚度数据中强烈的散射问题，创造性并系统地应用了统计“分箱”技术，有效剥离了燃料重定位等非肿胀机制的影响，成功提取出表征燃料本征体积膨胀的宏观趋势。 3. 稳健的模型构建与验证：通过比较多种函数形式，选择了物理意义合理且统计稳健的二次模型。尤为重要的是，利用独立的板平均肿胀数据（浸没法）对模型进行了验证，增强了结论的可信度。 4. 提供量化不确定性：不仅给出了最佳拟合曲线，还计算了95%置信区间和预测区间，为工程应用提供了至关重要的不确定性信息，体现了严谨的科学研究态度。 5. 明确的工程导向：研究始终围绕解决HPRR燃料工程认证中的实际难题——建立可靠的肿胀预测工具，具有明确的应用目标和价值。
七、 其他有价值内容 研究还指出，AFIP-6MkII实验的数据在再结晶区域的肿胀值似乎略高于其他实验，作者提出这可能与该实验较高的功率密度/温度、钼含量变化或杂质有关，并指出这将是未来研究的课题。同时，论文也提到，目前所有U-Mo燃料肿胀模型均未考虑功率或温度依赖性，如果未来实验发现这种依赖性，则需要开发适用于不同条件的额外肿胀关联式。这为后续研究指明了方向。