通过多尺度混合方法实现十亿网格油藏的高性能计算模拟

我们提出了一种基于非重叠区域分解技术的多尺度混合方法的三维并行实现，面向多核计算机，并通过数值实验评估了其计算性能。我们选取多尺度 Robin 耦合方法作为原型方法（许多其他方法可由此导出），并详细说明了其实现过程。基于 SPE10 基准，我们展示了在高度非均质各向异性岩层中，从数百万直至超过 20 亿计算单元的问题的数值结果。所提出的实现在局部问题和界面耦合系统中均依赖直接求解器。与基于代数多重网格预条件的最先进全局细网格求解器相比，我们在单相和两相流问题中发现了良好的弱可扩展性与强可扩展性。