利用多尺度预条件子加速多孔介质流动模拟的并行性能

《Mathematics and Computers in Simulation》：Speeding up parallel performance of porous media flow simulations using multiscale preconditioners

　　

本文采用多尺度Robin耦合方法（MRCM）开发了新的三维预条件算法；开展了不同预条件策略的对比研究；评估了大规模三维线性系统求解的配置；并在多进程并行计算环境中进行了实际应用。严格的测试程序验证了该方法的有效性与鲁棒性。

多孔介质中的流体流动模拟对于提高石油采收率等工程应用至关重要，此类模拟通常需要求解包含数十亿网格单元的大规模数值模型。尽管现代并行计算架构提供了算力支持，但由达西（Darcy）方程离散得到的大型稀疏线性系统仍是主要计算瓶颈。传统的预处理技术如代数多重网格（AMG）在对称正定问题中表现优异，但在渗透率对比度极高的复杂介质中，其收敛性和可扩展性可能退化。多尺度方法通过区域分解和低维界面问题可以在较低成本下获得近似解，然而精度随非均质性增强而下降。将多尺度方法作为预条件子，既能保留其可扩展性，又能保证求解精度。为此，本文首次系统研究了多尺度Robin耦合方法（MRCM）作为预处理共轭梯度法（PCG）预条件子的可行性，并针对三维高非均质多孔介质问题进行了大规模并行性能评估。该研究由Pericles R.P. Carvalho等研究人员完成，发表在《Mathematics and Computers in Simulation》。

在方法上，研究人员针对达西流方程建立椭圆边值问题，采用两通量近似（TPFA）有限体积法离散，并基于MRCM构造预条件子。MRCM通过在子域界面施加Robin型边界条件，结合多尺度基函数（MBF）和低维界面空间，形成代数系统。通过消去局部未知量，推导出关于压力变量的对称正定（SPD）预条件矩阵。由于原始MRCM解在子域界面上存在压力不连续，研究人员在预条件子前后引入平滑操作，包括Gauss-Seidel（GS）和逐次超松弛（SOR），并证明平滑后预条件子仍保持正定性。并行实现基于PETSc库，利用其PCSHELL接口将MRCM嵌入PCG迭代，并将MBF的计算限制在首次迭代中，后续迭代复用分解信息，界面系统由MUMPS直接求解器处理。数值实验使用了基于SPE10基准数据的各向同性（ISO）和各向异性（ANI）渗透率场，以及一个合成的十字形通道渗透率场，计算在巴西LNCC的Santos Dumont和CeMEAI的Euler超算上进行。

5.1 MRCM预条件子评估：通过一个简化的各向同性模型，改变区域分解、界面逼近阶次和Robin参数α，发现分段线性界面形状函数与α=10^-6的组合获得最少迭代次数，表明预条件子性能可通过参数调节，且对界面空间选择敏感。

5.2 各向同性模型（ISO）：在SPE10-derived各向同性场上，研究不同平滑策略（FB-GS、BF-GS、FB-SOR、BF-SOR）与AMG的性能差异。结果显示，SOR平滑优于GS，FB与BF差异很小；存在最优平滑迭代次数（约20~50次），过多会增加计算开销。弱扩展性分析中，MRCM的相对偏差σ/μ小于AMG，表明其扩展性更好。剖析各阶段耗时发现，MBF和界面求解占主导，而SOR能有效降低这两部分时间。在渗透率对比度变化实验中，当对比度降低时AMG耗时上升而MRCM耗时下降，表明MRCM更适应油藏典型对比度范围。

5.3 各向异性模型（ANI）：基于SPE10水平和垂直渗透率构建各向异性场。结果表明，在高对比度情形下AMG仍具优势，但MRCM在中等对比度下保持竞争性。可扩展性分析显示，MRCM的相对偏差随对比度增大而增大，但在适度对比度下仍表现良好。

5.4 三维十字形渗透率场的弱扩展性：该场由周期性十字形高渗通道组成，形成完全连通网络，对求解器极具挑战。弱扩展性测试表明，当核心分解使局部域尺寸与单元结构对齐（如15×15×15对应每个局部域8×8×8单元）时性能最佳。研究还考察了不同节点配置，发现增加节点数可降低内存竞争，提高性能。与AMG相比，MRCM的MBF时间扩展性极佳，总时间显著偏低。在更大规模（每核8×5³自由度）测试中，MRCM的弱扩展性与先前文献相当且相对偏差更小。

论文在讨论部分指出，MRCM预条件子结合优化的平滑策略能显著改善收敛行为和并行可扩展性。对于各向同性问题，MRCM在计算时间和扩展性上优于AMG，且SOR平滑和适当界面形状函数带来明显收益；在各向异性和高对比度场景中，AMG仍具竞争力，但MRCM在更适中的对比度下保持优势。十字形渗透率测试中，MRCM展现出优异的弱扩展性，甚至对超大规模问题也一致优于AMG。未来研究方向包括扩展到多相流、集成自适应和物理信息界面空间，以及进一步优化界面求解器以降低大规模核数下的成本。此外，采用近期提出的带过采样与平滑的MRCM（MRCM-OS）将大幅提高高度非均质渗透率场的计算效率。