Parameter-robust well-posedness and discretisation for coupled Darcy--Forchheimer and advection-diffusion-reaction equations

技术摘要：耦合 Darcy–Forchheimer 与平流-扩散-反应方程的参数鲁棒性适定性与离散化

问题陈述
本研究探讨了由 Darcy–Forchheimer 方程与平流-扩散-反应（ADR）方程组成的完全耦合系统的数学分析与数值离散化。该系统模拟了多孔介质中的非线性流体流动与溶质运移的耦合过程，其中流速 $u$ 与压力 $p$ 满足包含 Forchheimer 阻力项的非线性动量平衡，而溶质浓度 $\phi$ 与总通量 $\zeta$ 则满足输运方程。本文解决的一个关键挑战是稳定性与收敛性对模型参数的依赖性，具体包括渗透率（导水系数）$\kappa$、Forchheimer 系数 $F$、扩散系数 $\alpha$ 以及反应系数 $\varrho_0$。作者旨在建立在参数无关（即鲁棒）意义下的适定性，并推导误差估计，这些估计在实际应用（如储层建模和地热能优化）中参数可能跨越多个数量级的情况下依然保持稳健。

方法论
本文采用泛函分析方法，在反射 Banach 空间框架内进行研究，超越了标准的 Hilbert 空间框架，以适应问题的非线性和特定的定标要求。

1. 抽象理论： 作者借鉴了近期关于 Banach 空间中扰动鞍点问题的理论。他们扩展了 Brezzi–Braess 理论，以建立参数无关的稳定性界限。这涉及定义适用于速度、压力、通量和浓度的加权 Banach 空间。这些空间的范数是通过由模型参数（例如 $\kappa^{-1/2}$、$F^{1/3}$、$\alpha^{-1/2}$）加权的各个空间之和与交集构建而成的。
2. 适定性分析： 利用不动点论证，将耦合系统解耦为 Darcy–Forchheimer 子问题和 ADR 子问题。作者证明了单个子问题是适定的，随后证明了在数据足够小的情况下，其解算子的复合映射在速度空间的一个闭球上构成收缩映射。这确立了具有参数无关稳定性常数的连续解的存在性与唯一性。
3. 有限元离散化： 采用混合有限元方法对问题进行离散化。作者利用 Raviart–Thomas 元处理速度和通量空间，并使用不连续 Galerkin（破碎 Sobolev）空间处理压力和浓度。至关重要的是，离散空间配备了与连续分析中相同的参数加权范数。
4. 误差分析： 在这些加权范数下推导先验误差估计。分析依赖于离散 inf-sup 条件和离散不动点论证，以证明收敛速率是优的且对 $\kappa$、$F$、$\alpha$ 和 $\varrho_0$ 是鲁棒的。
5. 预条件处理： 为了应对求解由此产生的非线性及非对称线性系统带来的计算挑战，作者设计了块对角算子预条件器。这些预条件器受 Schur 补结构和算子预条件框架的启发，并针对耦合系统的特定参数权重进行了定制。

核心贡献

• 向 Banach 空间的扩展： 本文将用于扰动鞍点问题的 Brezzi–Braas 理论从 Hilbert 空间扩展到了反射 Banach 空间，特别是在参数鲁棒性的背景下。
• 参数鲁棒稳定性： 作者证明了耦合的 Darcy–Forchheimer/ADR 系统是适定的，且其稳定性界限与渗透率、Forchheimer、扩散及反应系数无关。这是通过精心构建的加权 Banach 空间和范数实现的。
• 鲁棒离散分析： 本文提供了严谨的离散分析，表明混合有限元格式能够保持这些鲁棒稳定性特性。推导出的先验误差估计是准最优的，且与模型参数及网格尺寸 $h$ 无关。
• 算子预条件处理： 提出了可变块对角预条件器的设计。研究表明，这些预条件器对于模型参数（除平流占优机制下的扩散系数外）以及离散化参数都是鲁棒的。

数值结果
理论发现通过 2D 和 3D 的一系列数值实验得到了验证：

• 收敛性： 通过构造光滑解和异质渗透率进行的测试，证实了在标准范数和参数加权范数下均具有最优收敛阶（$O(h^{k+1})$）。结果表明，虽然标准范数在极端参数值下可能表现出敏感性，但加权范数恢复了最优收敛行为。
• 预条件器性能： 对线性化系统的数值测试表明，所提出的预条件器在 $\kappa$、$F$、$\alpha$ 和 $\varrho_0$ 剧烈变化时，其条件数和 MINRES 迭代次数保持稳定。作者指出，在扩散占优机制下观察到了鲁棒性；但在平流占优机制（高 Péclet 数）下，性能并非完全鲁棒，这与所选稳定化方法的局限性一致。
• 实际应用： 该方法被应用于顶盖驱动腔流问题和瞬态混相置换问题（Hele–Shaw 池中的粘性指进），展示了该方法处理复杂流体-输运耦合及非线性的能力。

意义与主张
本文声称其主要意义在于提供了一个数学严谨的框架，用于分析和模拟多孔介质中耦合的流体与输运过程，且无论物理参数的大小如何，该框架都能保持稳定与精确。通过转向 Banach 空间分析并利用参数加权范数，作者克服了标准 Hilbert 空间方法的局限性，后者在处理此类非线性耦合系统时往往无法提供统一的稳定性。这项工作为开发高效且鲁棒的数值求解器提供了理论基础，适用于参数（如渗透率和反应速率）变化显著的情景，这在地下工程和环境建模中是一个常见场景。作者也谦逊地承认，虽然不动点框架产生了鲁棒的理论界限，但耦合系统中观察到的数值鲁棒性表明，某些加权范数中分析性的鲁棒性损失可能是定标或不动点框架本身的产物，而非该方法根本性的局限。