面向空间非均匀源与晶界分辨的裂变气体扩散广义谱算法及其在SCIANTIX中的实现与验证

《Nuclear Materials and Energy》：A generalized spectral algorithm for fission gas diffusion: Implementation and verification in SCIANTIX

　　

原文未提供独立摘要；以下依据论文开篇概述性亮点内容作专业翻译：在SCIANTIX（一种燃料性能程序）中引入裂变气体扩散（fission gas diffusion）的广义算法（generalized algorithm）。新算法放宽均匀源项（uniform source term）与均匀初始条件（uniform initial condition）假设。源项非均匀性（source heterogeneity）被证明会影响裂变气体滞留与释放（retention and release）。晶界（grain boundary, GB）分辨（resolution）通过局域源项（localized source term）表示。该框架可对局域亚晶尺度（sub-grain-scale）现象进行建模。

研究背景与问题：反应堆运行中铀裂变产生氙（xenon）和氪（krypton），累积裂变产额约为0.25–0.30 atoms per fission。由于其在氧化物燃料中溶解度极低，这些气体向晶界（grain boundary, GB）迁移，可在晶间气泡中积累或最终释放至燃料棒腔室；部分气体滞留于晶内，部分逃逸，从而影响肿胀和裂变气体释放（Fission Gas Release, FGR）。气体原子在晶内产生的位置显著影响其输运行为：近晶心产生的原子到达晶界所需时间长于近晶界产生者。因此，准确模拟裂变气体行为不仅需考虑扩散、俘获和分辨等机制，还需考虑裂变源的空间分布；这是预测性燃料性能评估和现代燃料性能程序（Fuel Performance Codes, FPCs）的核心功能。现有FPCs采用FORMAS、URGAS等算法；SCIANTIX采用谱扩散算法（spectral diffusion algorithm），其中Polypole-1和Polypole-2等依赖晶内均匀源假设，难以描述燃料样品气体注入、晶界分辨等局域现象。为此，研究人员在SCIANTIX中开发、实现并验证广义谱扩散求解器，使空间变化源项和非均匀初始条件可投影至谱基，同时保持谱方法计算效率；数值验证采用制造解方法（Method of Manufactured Solutions, MMS），重点展示源异质性和晶界分辨的物理影响。

关键技术方法（不超过250字）：研究人员将总气体浓度c_t(r,t)展开为空间本征函数ψ_n(r)与时间系数x_n(t)的线性组合；在球几何加权内积下对源项投影，采用分段线性源函数并对各空间区间贡献求和，同时以类似方式处理非均匀初始条件。数值验证采用MMS并与原始SCIANTIX求解器基准比较。案例采用Baker报道的分离效应实验和Zullo et al.模拟条件（1273 K、5500 h、恒定裂变率等），比较均匀源与时变随机非均匀源；晶界分辨将源分解为均匀裂变源S₁(r)与近晶界宽度λ的局域分辨源S₂(r)，比较λ=0、10、100 nm。

研究结果：
3.1 Case studies: Impact of source distribution：研究人员模拟1000多个受同一体积平均值约束的线性源函数，并选五个代表案例：均匀源S₀(r)为参考；S₁(r)、S₂(r)随半径线性减小，表示产气集中于晶心，S₁(r)梯度更陡；S₃(r)、S₄(r)随半径增大，表示产气集中于晶界。结果表明近晶心产气导致迁移延迟和更大晶内滞留；近晶界产气导致更快晶界积累和更早释放；晶界浓度曲线拐点对应FGR开始；相同总产气下源异质性本身即可显著改变行为，氪也呈类似趋势。进一步按Zullo et al.条件复现Baker实验：1273 K、5500 h、恒定裂变率、稳态温度与静水应力，采用Turnbull扩散率、异质形核、异质分辨、扩散俘获和Ainscough et al.晶粒长大模型；比较均匀源和时变随机非均匀源并保持体平均源相同。非均匀源下晶内氙浓度更高，表明晶内滞留增强；晶界库存早期增长较慢，随源分布演化近晶界产气比例增大而出现曲率变化；释放开始由均匀源约3300 h延迟至随机源约4700 h，最终释放库存更小。结论是源异质性即使体平均产率相同，也可显著影响晶内滞留、晶界库存和释放动力学。
3.2 Application: Grain boundary resolution：研究人员将总源分解为均匀裂变源S₁(r)和局限于晶界外侧宽度λ壳层的分辨源S₂(r)，并在基质输运方程中加入S₁(r)+S₂(r)，在晶界浓度方程中减去S₂(r)；晶界气体浓度c_GB的时间导数可为负。分辨源随晶界气体积累逐步发展，而非辐照开始即完全活跃；若一开始完全活跃，其影响会更大。模拟λ=0、10、100 nm，10 nm基线依据离子在物质中的阻止和射程（Stopping and Range of Ions in Matter, SRIM）数据和先前研究。随λ增加，氙从晶界转回晶内，晶内库存增加、晶界库存降低；释放开始约从无晶界分辨的3520 h变为λ=10 nm的3575 h和λ=100 nm的3740 h。Δc_t为正表示晶内滞留增强，Δc_GB和Δc_released为负表示晶界库存减少和释放延迟；偏差随λ单调增加，表明λ是控制分辨影响的关键参数。更长辐照和更高裂变率会增强该影响，趋势与L?s?nen一致。

讨论总结：广义谱求解器表明，空间源异质性即使体平均产率不变也能实质影响裂变气体输运，其影响取决于物理机制：Baker测试中时变非均匀源将释放开始由约3300 h延迟至4700 h；晶界分辨导致较温和但系统的延迟，释放开始由约3520 h变为3575 h和3740 h。晶界分辨源与裂变率成正比，提高辐照强度预期增强气体回注、增大晶内库存并更明显影响释放和肿胀；其重要性取决于辐照条件，在本文工况下主要影响释放时机，对总释放分数影响较温和。该方法价值不在于降低扩散系数、晶界行为、气泡演化等模型不确定性，而在于提供数值稳健框架以显式表示物理相关的空间源分布；低阶和低计算成本适合不确定性传播。当前假设包括晶粒内温度空间均匀但可随时间变化；燃耗和扩散系数由体平均裂变率评价并假设每时间步空间均匀；本征函数仅依赖几何和边界条件，与温度、燃耗和源强度无关且解析获得，无需本征函数库；源表示采用分段线性近似，高度局域剖面可能需更多区域和谱模；仍保留准稳态近似。方法学上可视为谱降阶模型，类似本征正交分解-Galerkin投影（Proper Orthogonal Decomposition with Galerkin Projection, POD-GP），但使用解析本征函数；当前限于SCIANTIX等效球近似，方法不限于球形，可扩展至柱状晶等。计算开销可忽略，适合OFFBEAT、TRANSURANUS、z3st等工程尺度FPCs；晶界分辨只是应用之一，还可用于非均匀气泡分布、晶界附近贫气泡区、注入剖面和异质微结构。

研究结论翻译：本研究开发并在SCIANTIX中实现了一种广义谱扩散算法，用于处理空间变化源项和非均匀初始条件，同时保持原谱公式的计算效率。所提框架能够表示在常规空间均匀产气假设下无法捕捉的局域亚晶尺度现象。代表性应用表明，仅源异质性即可显著改变裂变气体输运；特别是时变非均匀源显著改变了释放动力学，而随着特征分辨深度增加，晶界分辨逐步延迟了氙释放的开始。由于广义求解器仅引入可忽略的计算开销，其适合并入工程尺度燃料性能程序。除晶界分辨外，该方法为注入剖面、异质微结构、局域气泡群及其他空间依赖现象建模提供了灵活基础。因此，本工作扩展了谱扩散方法的适用性，同时保持常规燃料性能模拟所需效率。