Stream function -- pressure virtual element methods for the Stokes--Darcy interface problem

技术摘要：用于 Stokes–Darcy 界面问题的流函数–压力虚单元方法

问题陈述
本文研究了耦合 Stokes–Darcy 界面问题的数值模拟，这是控制自由流区域（Stokes）与多孔介质（Darcy）之间相互作用的基础流体力学挑战。这种耦合通过法向速度连续性（质量守恒）、法向应力平衡以及 Beavers–Joseph–Saffman 切向滑移条件来实现。虽然存在各种有限元（FE）、间断 Galerkin（dG）和混合高阶（HHO）方法，且这些方法通常依赖于需要满足 inf–sup 稳定对的速度–压力表述，但本研究专注于在二维单连通区域中的一种特定重新表述。在此，通过引入标量流函数，在式中恒等满足不可压缩约束，从而从 Stokes 方程中消除了压力变量，并将系统简化为与二阶 Darcy 压力方程耦合的四阶双调和问题。

方法论
作者提出了一种针对该流函数–压力表述的新型虚单元方法（VEM）离散化方案。该方法利用了 VEM 处理通用多边形网格的能力，这对于处理复杂的界面几何形状和非匹配网格特别有利。

1. 表述： 连续问题被转化为原初–原初（primal-primal）形式。Stokes 区域使用流函数 $\chi \in H^2$，而 Darcy 区域使用压力势 $\phi \in H^1$。弱形式导致了一个扰动鞍点问题。
2. 离散化：
   • Stokes 区域： 对流函数采用 $C^1$ 一致最低阶 VEM 空间。这涉及二次多项式（$P_2$），并要求满足 $H^2$ 一致性，通过由节点值和节点梯度定义的增强虚拟空间来实现。
   • Darcy 区域： 使用基于线性多项式（$P_1$）的 $C^0$ 一致最低阶 VEM 空间进行压力计算。
   • 耦合： 界面条件在弱形式中自然实现。该方法无需重构网格即可处理界面处的非匹配网格。
3. 实现： 该方案在开源面向对象库 vem++ 中实现。利用了特定的 vem::vem2dDofFilter 类来管理不同子域之间的自由度（DoF），从而允许在单个网格上同时定义不同的 VEM 空间并基于单元标记进行过滤。
4. 稳定化： 离散双线性形式包含稳定项（$S_S^K$ 和 $S_D^K$），以确保在投影算子的核上满足强制性，遵循标准的 VEM 做法（例如对角线配方和 DOFI–DOFI）。

主要贡献
论文概述了五项主要贡献：

• 连续性分析： 开发了流函数–压力表述下 Stokes–Darcy 界面问题的连续性分析，这本质上避免了对 inf–sup 稳定性条件的需求。
• 离散方案： 设计了一种 $C^1$–$C^0$ 最低阶 VEM 离散化方案。作者强调了该方案的计算效率，指出其仅涉及 $3N_S + N_D + 4N_\Sigma$ 个自由度（其中 $N_S$、$N_D$ 和 $N_\Sigma$ 分别是 Stokes 区域内部顶点、Darcy 区域内部顶点以及界面上的顶点计数）。
• 误差分析： 推导了所提离散方案的先验误差估计，证明了在标准正则性假设下的最优收敛率。
• 软件实现： 将新功能集成到 vem++ 库中，通过网格过滤促进界面问题的简便实现。
• 数值验证： 通过数值实验验证了理论发现。

结果
数值实验验证了理论收敛率和该方法的鲁棒性：

• 制造解： 使用制造解进行的两个实验（一个具有零界面值，一个具有非零界面值），在各种网格类型（四边形、非凸、Voronci 以及扰动 Voronoi 网格）上展示了误差的渐近线性衰减（$O(h)$），这与最低阶方案的理论预测一致。
• 死端过滤： 四分之一环形死端过滤器的模拟成功恢复了以往研究的结果，证明了该方法处理弯曲界面和混合边界条件的能力。
• 生物人工器官： 在渗透系数变化的生物人工器官（封装式）中的血液流动模拟，展示了该方法模拟复杂生理流体以及渗透率对跨界面通量影响的能力。

意义与主张
论文将自身工作定位为填补了关于使用虚单元扩展流函数表述至耦合 Stokes–Darcy 问题的文献空白。作者声称，所提出的结构提供了一种“完全不受 inf–sup 限制的策略”，与传统的速度–压力方法相比，简化了离散空间的构建。利用 VEM 框架处理不规则界面几何形状和非匹配网格的能力被强调为一大关键优势。

作者对未来的挑战保持谦逊，明确指出将此流函数表述扩展到多连通区域（在此类区域中，流函数无法通过标准边界条件唯一确定）以及将该框架扩展到三维空间（需要矢量势和矢量值函数的高阶 $C^1$ 空间）仍是开放的挑战。此外，集成无稳定化技术以及推导鲁棒的后验误差估计器也被确定为未来的研究方向。