聚变液态金属增殖包层内多区域磁流体动力学流动中电势耦合的分离式和块耦合方法

《Fusion Engineering and Design》：Segregated and block-coupled approaches for electric potential coupling in multi-region magnetohydrodynamic flows within fusion liquid metal breeding blankets

　　

聚变反应堆中的液态金属增殖包层在强磁场下运行，这导致出现由电流驱动的复杂磁流体动力学（MHD）效应，该效应耦合了流体区域和固体区域。准确建模这种电磁耦合对于可靠预测聚变包层中的MHD流动至关重要。本研究比较了多区域MHD流动中电势耦合的两种数值策略：一种块耦合方法（block-coupled approach）和一种分离式方法（segregated approach）。块耦合方法同时求解所有区域中的控制方程，提供强界面耦合和快速收敛，但要求在流体-固体界面处具有严格的网格一致性，这在复杂几何形状中使用自动网格划分工具时难以实现。一种分离式多区域求解器依次求解控制方程，通过迭代界面耦合容忍非一致性网格。该分离式方法得到验证，并应用于ITER水冷锂铅测试包层模块（WCLL TBM）模型中的液态金属流动模拟。数值结果与实验观测吻合良好，并确认了增殖单元之间高度非均匀的液态金属流量分布。本研究强调了两种方法各自的优势和局限性，并为在ITER和DEMO包层系统的高保真模拟中选择合适的电势耦合策略提供了指导。

**论文解读文章**

**研究背景与问题**

在聚变反应堆中，液态金属增殖包层（如水冷锂铅包层，WCLL）是维持氚自持的关键部件。这些包层在强磁场（约数特斯拉）下运行，导致液态金属（如PbLi）流动时产生复杂的磁流体动力学（MHD）效应。具体而言，导电流体在磁场中运动会产生感应电流，这些电流在流体与固体结构（如加强板、隔板）之间耦合，形成洛伦兹力，显著改变流动特性，并导致巨大的MHD压降。准确建模这种多区域电磁耦合对于可靠预测包层内的流动分布、氚输运和热工水力性能至关重要。

然而，现有数值方法存在挑战。块耦合方法（block-coupled approach）虽能同时求解所有区域方程，提供强耦合和快速收敛，但要求流体-固体界面处的网格严格一致。在ITER和DEMO包层等复杂几何中，自动网格划分工具难以保证这种一致性，手动块结构化网格则大幅增加预处理工作量并限制几何灵活性。因此，需要一种更灵活、鲁棒的数值算法，以支持包层设计迭代中的高效模拟。本研究旨在比较两种数值策略，并重点开发一种分离式求解器，以克服块耦合方法的局限性。

**研究内容与意义**

研究人员开展了块耦合方法与分离式方法的对比研究，并在OpenFOAM框架下开发并验证了分离式多区域MHD求解器。该求解器通过迭代界面耦合容忍非一致网格，适用于复杂几何。研究将分离式方法应用于ITER水冷锂铅测试包层模块（WCLL TBM）模型的MHD流动模拟，揭示了增殖单元（BU）之间高度非均匀的流量分布，并与实验观测一致。该工作为选择合适的高保真模拟策略提供了指导，对未来的包层设计优化具有重要意义。论文发表在《Fusion Engineering and Design》。

**关键技术方法**

研究人员为开展研究用到以下主要关键技术方法：
1. **有限体积法（FVM）**：在OpenFOAM开源代码中实现，采用基于单元的有限体积法离散控制方程。
2. **块耦合方法**：基于OpenFOAM-extend的BlockLduMatrix类，构建全局块矩阵，同时求解流体和固体区域的电势泊松方程，实现强界面耦合。
3. **分离式方法**：基于OpenFOAM 7开发，通过迭代界面条件（如公式(9)-(10)）依次求解各区域，并利用任意网格界面（AMI）插值处理非一致网格。
4. **验证案例**：采用Sloan-Smith（SS66）解析解（厚壁导电Hartmann壁管道）和Hunt流（薄壁近似）进行验证，并与块耦合求解器结果对比；使用网格收敛指数（GCI）方法评估数值误差。
5. **实验数据来源**：模拟对象为KIT MEKKA实验室的WCLL TBM模型，该模型装有电势传感器和压力传感器，用于测量MHD流动下的压降和流速。

**研究结果**

**2. 数值方法（Numerical approaches for electric potential coupling in multi-region MHD flows）**
研究人员归纳了MHD流动的控制方程，包括流体与固体区域的动量守恒、质量守恒、欧姆定律及电荷守恒，并推导出电势泊松方程。通过界面处电势和法向电流连续性条件，定义了耦合权重因子ω（见公式(3)），该因子取决于流体与固体的电导率比和网格间距比。两种方法均基于此耦合公式：块耦合方法在矩阵层面耦合所有区域（公式(8)），而分离式方法通过迭代更新界面值（公式(9)-(10)）。

**2.1 块耦合方法（Block-coupled approach）**
块耦合方法将流体和固体区域的电势方程组成一个全局系统矩阵，同时求解，实现强隐式耦合，收敛速度快。但该方法要求界面网格严格一致，且内存需求高、模块化差，难以在复杂几何中应用。并行分解时，界面两侧单元必须位于同一处理器，限制了负载平衡。

**2.2 分离式方法（Segregated approach）**
分离式方法在每个时间步内独立求解各区域，通过迭代界面边界条件传递信息（显式处理）。它容忍非一致网格，可利用AMI插值，适用于自动网格生成。该方法内存需求低，易于实现，基于OpenFOAM Foundation分支，可访问自适应网格细化等先进功能。但需注意下松弛或PIMPLE算法以保证收敛。

**3. 验证与收敛性（Validation and convergence）**
研究人员通过三个案例验证分离式求解器：
- **案例I**：厚壁导电Hartmann壁管道中的MHD流动，与Sloan-Smith（SS66）解析解对比。网格收敛性研究表明，在Ha=500、tw/L=0.01、σ_S/σ=1时，600×600网格下压力梯度误差仅0.12%，最大速度误差0.04%。不同哈特曼数（Ha=500~10000）下误差均小于0.36%。
- **案例II和III**：分别为方形管道和圆形管道，与经严格验证的块耦合求解器结果对比，差异极小（案例II中压力梯度差异0.005%，案例III中0.16%），验证了分离式方法的准确性。
- 收敛行为比较：对于强耦合界面（ω=0.926），分离式方法所需执行时间比块耦合方法高约两个数量级；但对于弱耦合界面（ω=0.943），分离式方法的计算成本反而低约两个数量级。界面网格尺寸比δ_S/δ影响收敛速度，细化流体侧网格可改善收敛。

**4. 应用于TBM模型模拟（Application to simulation for a TBM mock-up）**
研究人员将分离式求解器应用于WCLL TBM模型（无冷却管）的MHD流动模拟（Ha=500，N=6.9×10⁴）。数值结果揭示了高度非均匀的流量分布：
- 流线分析显示，大部分流动（红色流线）直接进入BU1，流向第一壁后反转，进入大的排水歧管；另一部分（蓝色流线）沿给水歧管流至BU8，然后进入排水歧管。BU2~7中流量可忽略。
- 速度场表明，BU1和BU8中存在沿平行于磁场方向壁面的射流结构，给水歧管中的速度沿极向逐渐降低，排水歧管中速度增加，直至两者平衡。
- 定量流量分配（表6，图8）显示：约48%的流量通过BU1，10%通过BU2，BU3~6几乎无流量（各约0.3%~2%），BU7传递6.6%，BU8传递剩余的32%。该结果与实验测量（基于电势的间接流量估算）一致。

**总结讨论与结论**

**讨论部分总结**：研究指出，如此非均匀的流量分配不利于TBM中氚的有效提取，可能导致局部氚浓度过高或低效。建议通过修改给水和排水歧管的几何设计来优化，例如调整隔板沿极向的位置，以改善流量均匀性。

**研究结论翻译**：本研究的动机源于对高效、鲁棒数值算法的需求，以预测液态金属增殖包层中的多区域MHD流动。研究人员研究了两种电势耦合建模方法：块耦合方法和分离式方法。块耦合方法同时求解所有区域的控制方程，在简单几何中提供强隐式耦合和快速收敛，但要求界面网格严格一致，这限制了其在自动网格生成用于复杂包层几何时的适用性。为克服此局限，研究人员开发并验证了分离式求解器。通过独立求解各区域并采用迭代界面条件，分离式方法容忍流体-固体界面处的非一致网格，因此适用于复杂几何，可借助自动网格生成技术，大幅减少网格生成工作量。随后，将分离式求解器应用于WCLL TBM模型的MHD流动模拟，展示了其在复杂几何中的预测能力。结果揭示了增殖单元之间高度非均匀的液态金属流量分布，大部分流量通过BU1和BU8，且存在偏向BU1的显著不对称性。全面的数值结果及与实验数据的详细比较将在后续专门出版物中报道。两种方法的比较评估突出了它们在精度、鲁棒性和计算性能方面的各自优势与局限。基于这些发现，建议在简单几何的强耦合多区域MHD流动中采用块耦合方法，而在需要自动网格生成的复杂几何中，分离式方法更为合适。