MOOSE框架中THM与子通道耦合技术实现解析

背景与意义

在核反应堆热工水力分析领域，系统级代码（如THM）与子通道级代码（如Pronghorn-SC）的耦合计算具有重要意义。MOOSE框架作为多物理场仿真平台，通过其多应用程序（multiAPP）功能实现了这一跨尺度耦合能力。本文以ORNL-19燃料组件为案例，深入剖析耦合实现的技术细节。

耦合架构设计

耦合系统采用双向数据交换模式：

1. 子通道→系统级
   子通道计算模块对燃料组件进行全三维求解，输出组件平均压降数据。该数据通过MOOSE的multiAPP接口传递至THM系统代码，THM据此动态调整核心区摩擦因子，确保系统级模型能准确反映子通道尺度的流动特性。

2. 系统级→子通道
   THM向子通道模块提供核心入口质量通量和出口压力边界条件，形成闭合的耦合循环。这种设计既保留了系统级计算的全局特性，又嵌入了子通道尺度的局部精细化模拟。

关键技术实现

1. 数据映射算法
   采用基于几何重叠域的变量传递方法，确保不同分辨率网格间的数据一致性。特别处理了轴向位移问题（案例中THM模型轴向偏移0.5m），这是可视化需求与计算精度平衡的典型实践。

2. 流体物性处理
   当前实现中采用临时方案处理钠冷却剂物性，后续需完善为正式物性库接口。这反映了多物理耦合中物性传递的通用挑战。

3. 流动方向假设
   基于快堆设计特点，采用单向流动假设简化耦合计算。该假设要求确保核心出口不出现反向流动，是工程简化与计算精度权衡的典型案例。

工程应用价值

该耦合方案已成功应用于：

• 系统级瞬态分析与子通道局部热工参数的协同预测
• 全堆芯尺度下的热点定位与DNBR分析
• 多尺度耦合中的不确定性量化研究

发展展望

未来改进方向包括：

1. 完善钠冷却剂物性模块的标准化接口
2. 开发反向流动处理算法
3. 优化耦合迭代收敛策略
4. 扩展应用于其他冷却剂类型的通用耦合框架

该技术的实现标志着MOOSE框架在核工程多尺度仿真能力上的重要突破，为先进反应堆设计提供了强有力的分析工具。