**耦合Stokes–Darcy系统杂交可间断Galerkin离散的参数鲁棒预处理器：算子预条件框架下的构造、凝聚分析与数值验证**

《COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING》：Preconditioning of a hybridizable discontinuous Galerkin method for the coupled Stokes–Darcy system

　　

研究人员提出了一种预处理器，用于处理耦合Stokes–Darcy系统（coupled Stokes–Darcy system）杂交可间断Galerkin离散（hybridizable discontinuous Galerkin discretization，HDG离散）所产生的静态凝聚（static condensation）线性系统。该预处理器被设计为关于网格尺寸（mesh size）、运动黏度（viscosity）与渗透率（perme量）具有鲁棒性，其前提假设是ακ^?1/2（其中α为Beavers–Joseph–Saffman定律（Beavers–Joseph–Saffman law）中的常数，κ为渗透率常数）受某一常数界限制。该设计策略首先应用算子预条件框架（operator-preconditioning framework）为非凝聚离散构造预处理器，具体途径是证明该离散格式存在与网格尺寸、运动黏度及渗透率无关的有界性常数与稳定性常数。随后，研究人员利用所提出的杂交可预条件框架（hybridizable preconditioning framework），证明了所得凝聚预处理器应用于约化（reduced）线性系统的鲁棒性。数值算例表明，所提出预处理器具有鲁棒性。

**一、研究背景与问题提出**

耦合Stokes–Darcy系统用于刻画自由流区域与多孔介质区域之间的流体相互作用，在地下水文学与能源勘探等领域具有重要应用价值。该模型在界面Γ_I上需满足法向速度连续、法向应力连续以及Beavers–Joseph–Saffman界面切向滑移条件，其中滑移项由αμκ^?1/2刻画，α为实验确定的无量纲常数。为求解该模型，研究人员采用杂交可间断Galerkin（hybridizable discontinuous Galerkin，HDG）方法——一种允许静态凝聚（static condensation，即消去局部单元自由度、仅保留整体面未知量的技术）的间断Galerkin方法——进行离散。静态凝聚中最昂贵的步骤是求解面未知量对应的约化线性系统。

目前存在的问题是：尽管参数鲁棒预处理器已被成功设计用于HDG离散的Stokes问题与Darcy问题各自独立的情形，对于耦合Stokes–Darcy问题的HDG离散，如何构造同时关于网格尺寸h、运动黏度μ与渗透率κ鲁棒的预处理器，此前尚无解决方案。这正是本研究要开展的核心问题。

**二、研究内容、主要结论与意义**

在假设ακ^?1/2有界（即存在与μ、κ、α无关的常数c_α使得α≤c_ακ^1/2）的条件下，研究人员提出了针对HDG离散静态凝聚后约化线性系统的预处理器，并证明了其关于h、μ、κ的鲁棒性。研究分两步展开：首先应用算子预条件框架（operator preconditioning framework，即基于加权范数下适定性分析构造预处理器的框架）为非凝聚离散构造预处理器，途径是证明离散格式的有界性常数与稳定性常数与h、μ、κ无关；随后利用Henríquez等研究人员提出的杂交可预条件框架，证明了所得预处理器Schur补构成的凝聚预处理器对约化系统同样鲁棒。数值算例验证了理论结果。其重要意义在于：这是首个针对耦合Stokes–Darcy系统杂交离散的预处理器，且既不依赖分数算子，也无需在有限元多项空间中施加界面条件，从而简化了实现。

**三、关键技术与方法（约250字）**

研究人员主要采用以下关键方法：其一，算子预条件框架，即在依赖问题参数的加权范数下证明离散格式的一致适定性（uniform well-posedness，指离散问题适定性常数与网格及物理参数无关），由此诱导出块对角预处理器；其二，杂交可预条件框架，通过验证Schur补条件，将非凝聚预处理器的鲁棒性传递至静态凝聚后的约化系统；其三，加权内积（weighted inner products）构造，针对Stokes区速度、Darcy区速度与压力及界面迹量分别定义含μ、κ、α权重的离散内积；其四，局部求解器与静态凝聚技术，通过单元局部问题消去体积自由度并导出约化变分形式；其五，基于NGSolve的数值实现，以MINRES（最小残差法）与GMRES（广义最小残差法）迭代求解，不精确版本结合辅助空间预条件（auxiliary space preconditioner，ASP）与Hypre的BoomerAMG代数多重网格处理迹速度块与迹压力块。

**四、研究结果**

（一）HDG离散的适定性（Well-posedness）。研究人员构造了一族加权内积及相应加权范数，证明了HDG离散双线性形式a_h(·,·)的有界性（boundedness，Lemma 1）与inf-sup稳定性（inf-sup stability，即Ladyzhenskaya–Babu?ka–Brezzi型速度-压力耦合条件，见Lemma 2与Lemma 3），其常数与h、μ、κ无关。据此得到定理1：在条件α≤c_ακ^1/2下，HDG离散关于h、μ、κ一致适定（Brezzi条件满足）。

（二）预处理器构造（Preconditioners）。基于上述加权内积，研究人员定义了具有块结构的预处理器P（含Stokes区迹速度块、Stokes区迹压力块、Darcy区迹压力块等四个对角块），由适定性分析可知P是静态凝聚前HDG离散关于h、μ、κ的鲁棒预处理器；进而将其Schur补（Schur complement，即块矩阵消去自由度后得到的补系统算子）S_P作为约化系统的候选预处理器。为分析S_P，研究人员引入Stokes区速度辅助问题与Darcy区压力辅助问题的内部罚HDG（interior penalty HDG）离散，并证明了相应的椭圆 coercivity 性质与inf-sup条件。

（三）凝聚预处理器的鲁棒性（Robustness of the condensed preconditioner）。通过静态凝聚导出约化变分形式，研究人员在定理2中验证了关键条件：对任意迹未知量?_h∈X?_h，由局部问题恢复的体积自由度范数满足‖(?A₁₁^?1A₂₁^T?_h,?_h)‖_Xh≤c_l‖?_h‖_X?h，从而条件数cond(S_P^?1S_A)≤c_bc_l²/c_i与网格及物理参数无关，即S_P是关于h、μ、κ鲁棒的凝聚预处理器。

（四）替代预处理器（An alternative preconditioner）。研究人员进一步提出更接近HDG双线性形式的替代内积，定义替代预处理器P?及其Schur补S_P?。定理3证明：在α≤c_ακ^1/2假设下，S_P?同样关于h、μ、κ鲁棒，且数值实验表明其性能优于S_P（尤其在MINRES迭代中），故后续实验均采用S_P?^?1。

（五）制造解算例（Manufactured solutions）。在单位方形区域Ω=[0,1]²（Stokes区与Darcy区各占一半）上，构造制造解（manufactured solutions，即预设解反推源项与边界条件的验证方法）进行验证。h-鲁棒性实验显示：随网格从h=2^?4细化至2^?8，精确与不精确形式的预处理器迭代次数基本不增长，例如S_P?^?1 preconditioned MINRES迭代次数仅在113至116次之间。μ与κ鲁棒性实验显示：当log₁₀(μ)∈[?4,4]、log₁₀(κ)∈[0,?4]时，MINRES与GMRES迭代次数仅随参数呈现有限增长。参数扫描（α∈{0.1,1,10,100}，κ∈{1,10^?2,10^?4}，多项式次数k∈{2,3,4}）表明：MINRES在ακ^?1/2足够大时出现中等程度的性能退化，而GMRES在测试范围内相对不敏感。

（六）非均质渗透率算例（Heterogeneous permeability case）。针对空间变化渗透率κ(x₁,x₂)的非均质多孔介质场，在含32768个单纯形的固定网格上、k=2、μ=0.1、α=1的配置下，验证了S_P?^?1关于h（2^?4至2^?8）与μ（log₁₀(μ)从0降至?6）的鲁棒性，精确与不精确两种形式均表现出参数鲁棒性。

**五、讨论与结论**

讨论部分指出，与已有依赖界面Lagrange乘子或界面分数算子的预处理器相比，本文预处理器不要求将界面条件施加于有限元多项空间，实现更为简便；同时，理论分析未涉及逼近空间多项式次数k的鲁棒性。研究结论部分可概括翻译如下：研究人员介绍了针对带Beavers–Joseph–Saffman界面条件的耦合Stokes–Darcy模型HDG离散的预处理器，在α≤c_ακ^1/2条件下关于h、μ和κ具有鲁棒性。首先利用算子预条件框架为非凝聚线性系统给出了预处理器；随后利用杂交可预条件框架证明了：非凝聚问题预处理器的Schur补构成对静态凝聚HDG离散同样鲁棒的预处理器（条件同为α≤c_ακ^1/2）。所提出预处理器避免了分数算子的使用以及界面条件在有限元空间上的施加，从而简化了实现。数值实验结果支持了理论分析。