《International Journal of Hydrogen Energy》：FESTIM v2.0: Upgraded framework for multi-species hydrogen transport and enhanced performance
编辑推荐：
为解决氢同位素在聚变能源材料中输运过程复杂、多物理场耦合模拟需求迫切的问题，研究人员开发了FESTIM v2.0开源有限元框架。该研究实现了多物种完全耦合输运、高级陷阱与反应网络建模、界面处理与外部工具（如OpenMC、OpenFOAM）集成，相比前代性能提升约15倍，为工程尺度模拟及更广泛的储氢技术研究提供了强大平台，对聚变燃料循环设计与安全评估具有重要意义。
开发聚变能源是实现清洁、高效能源的重要途径，而氢同位素（特别是氚T）的输运、滞留行为是聚变反应堆材料设计与安全运行的核心挑战。氚在聚变堆的第一壁、增殖区等部件材料中的扩散、捕获与释放过程，直接影响氚燃料循环效率、放射性废物总量和系统安全性。准确预测这些过程，对于确保满足法规要求、指导实验验证以及优化燃料循环部件设计都至关重要。然而，现实材料系统通常涉及多种氢同位素（如氢H、氘D、氚T）共存、复杂的缺陷捕获机制、不同材料间的界面输运以及多种物理效应（如热耦合）的相互作用，使得精确建模变得异常复杂。
为此，麻省理工学院等离子体科学与聚变中心的James Dark等研究人员，在《International Journal of Hydrogen Energy》上发表了关于开源有限元框架FESTIM的重大升级版本FESTIM v2.0的研究。这项研究的核心目标是扩展和提升该框架的物理建模能力与软件性能，使其能够应对上述复杂挑战。该工作是对原有FESTIM v1版本的重构与扩展。
研究人员开展了对FESTIM框架的全面重设计。他们基于新一代FEniCSx和DOLFINx软件栈进行了底层重构，这解决了v1版本依赖已停止维护的旧版软件所带来的可持续性、可扩展性和性能限制。为了实现更复杂的物理过程模拟，v2.0引入了全新的、模块化的多物种耦合输运模型。其最显著的进展在于，框架现在能够模拟多种氢同位素及其相互作用，包括多级占据陷阱、同位素交换、衰变和对流等复杂过程，这些在v1中由于架构限制无法实现。同时，针对多材料界面问题，除了保留v1的变量变换法以兼容历史研究外，还集成了间断伽辽金（Discontinuous Galerkin）和尼切（Nitsche）方法，这些成熟的数值方法极大地提升了多材料模拟的计算效率与精度。为了增强其工程适用性，FESTIM v2.0还专门设计了与外部模拟器（如蒙特卡罗中子输运代码OpenMC和计算流体力学软件OpenFOAM）的耦合接口，支持在更大规模的多物理场工作流中使用。在性能上，通过对代表性多材料扩散问题的基准测试，新版本实现了相对于FESTIM v1高达约15倍的加速，这使得大规模的工程级模拟成为可能。此外，研究团队为该版本建立了持续更新的验证与验证资源库，并通过活跃的开源社区模式（在GitHub上开发）确保了代码的长期可持续性和可复现性。
核心方法
作者为开展研究并实现FESTIM v2.0的功能，主要运用了以下关键技术方法：
- 1.
基于DOLFINx的模块化软件架构：完全重构了底层代码，从依赖旧版FEniCS迁移至现代化的DOLFINx，获得了更好的可扩展性、互操作性和长期维护能力。
- 2.
基于反应的灵活物理建模框架：引入了“反应”（Reaction）概念，将氢同位素的扩散、陷阱捕获与释放、衰变、同位素交换等物理过程统一描述为化学反应形式，从而通过组合反应构建复杂的反应网络。这为实现多物种耦合输运和高级陷阱模型奠定了基础。
- 3.
先进的数值界面处理方法：为了高效处理不同材料间溶解度差异导致的浓度不连续问题，FESTIM v2.0集成了数值分析方法，包括间断伽辽金/Nitsche方法以及罚方法，这些方法替代了原有的变量变换法，显著降低了多材料模拟的计算成本。
- 4.
外部求解器的松耦合集成：通过模块化设计，使得温度场、速度场等物理量可以从外部代码（如OpenMC用于中子辐照损伤，OpenFOAM用于流体对流）读入并耦合进FESTIM的输运模拟中，从而构建多物理场工作流。
研究结果
几何、网格与子域（Section 3: Geometry, meshes and subdomains）
该研究展示了FESTIM v2.0灵活的几何处理能力。它支持从简单1D线段到复杂3D结构（可通过GMSH、SALOME等软件生成）的多种网格格式，并通过子域（Subdomain）对象将物理属性（如材料参数）与几何区域解耦。这使得用户能够更自然地构建多材料问题。
控制方程与多物种框架（Section 4: Governing equations）
研究提出了一个通用的偏微分方程框架（Eq. 1），用于描述每个物种的浓度演化，包含扩散、源项、耦合反应和对流项。新框架的关键创新是支持多种“移动”和“固定”物种。例如，空陷阱位点可以被定义为“隐式”物种，其浓度可通过总位点浓度与占据位点浓度之差直接计算（Eq. 2）。这种设计使得模拟氢同位素与非氢物种（如空位、空位团、间隙原子）的输运成为可能。
反应网络与陷阱模型（Section 4.2: Reactions）
通过新引入的“反应”框架，FESTIM v2.0可以模拟复杂的相互作用。
- •
基本捕获（Trapping）：研究定义了经典的麦克纳布-福斯特（McNabb and Foster）模型，通过正向（捕获）和反向（脱陷）反应速率常数（通常遵循阿伦尼乌斯定律（Eq. 8, Eq. 9）来描述氢原子与陷阱位点的相互作用（Eq. 5）。
- •
多级捕获（Multi-level trapping）：该框架支持将多级（或称为多占据）捕获现象建模为一组串联反应（Eq. 10a-c）。例如，一个可容纳最多三个氢原子的缺陷，其填充分别由三个独立的、具有不同结合能的反应来描述。这种模型可以更真实地模拟辐照缺陷等复杂陷阱行为。
- •
同位素交换（Isotope swapping）：模拟了例如氚（T）被陷阱捕获后，可被移动的氢（H）原子替换下来的过程（Eq. 12），这对于预测混合同位素系统中的氚存量至关重要。
- •
衰变（Decay）：放射性衰变可被建模为一种无产物的单向反应（Eq. 13），以考虑氚的衰变效应。
对流效应（Section 4.3: Advection）
研究扩展了输运方程，加入了速度场驱动的对流项（Eq. 14），并展示了在扩散主导和对流主导两种不同佩克莱特数（Peclet number）情况下物种浓度分布的演变（）。这为模拟流体环境中氢的输运提供了可能。
界面条件与数值方法（Section 5: Interface conditions）
为解决多材料界面处氢浓度的不连续问题，FESTIM v2.0引入了新的数值处理方法。对于相同类型的材料（如均遵循西维特定律（Sievert's law）或亨利定律（Henry's law）），可以使用间断伽辽金/Nitsche方法（Eq. 17）。这种方法以弱形式在界面上强制执行浓度平衡条件，同时允许浓度场在跨材料界面时发生跳跃，是处理材料属性突变问题的稳健数值方案。
结论与讨论
FESTIM v2.0代表了氢同位素在材料中输运建模领域的一个重大进展。它从一个单一物种的简单扩散-捕获模型，演变为一个全面、灵活且高性能的多物理场、多物种模拟平台。该研究的核心结论在于，通过彻底重构软件架构，引入模块化的多物种与反应框架，并采用先进的数值界面处理技术，成功突破了前代软件在物理能力和计算效率上的局限。
新框架的意义深远：首先，它为研究聚变能源领域的关键问题——如不同氢同位素在等离子体面向材料中的竞争性捕获与交换行为、复杂缺陷网络的演化对氚存量的影响、以及跨材料界面的输运——提供了前所未有的建模能力。其次，高达15倍的性能提升使得对工程尺度（米级）部件（如增殖包层模块）进行高保真模拟成为现实，这对于设计和优化聚变燃料循环系统至关重要。再者，该框架并不局限于聚变领域，其模块化设计使其同样适用于核裂变系统的氢致失效分析、储氢材料优化等更广泛的氢技术研究。最后，通过与OpenMC、OpenFOAM等外部工具的集成，FESTIM v2.0可以被嵌入更大的系统级模拟工作流中，实现中子学-热工-氢输运等多物理场的强耦合分析。
这项研究不仅提供了一款强大的开源软件工具，更通过其清晰、可扩展的架构为未来模型的进一步扩展（如考虑应力辅助扩散、索雷效应等）铺平了道路。其活跃的社区开发模式、持续的验证与验证体系以及详尽的文档，确保了该平台能够作为未来氢同位素输运研究的长期可靠基础。
