Distributed optimal control problems governed by poroelasticity equations

技术摘要：受孔隙弹性方程控制的分布式最优控制问题

问题陈述
本文研究了由 Biot 孔隙弹性模型控制的分布式最优控制问题的数值近似。该物理系统描述了孔隙流体与多孔介质变形之间的相互作用，这一模型在地球物理学、生物组织力学和碳封存等应用领域至关重要。特定的控制问题旨在最小化一个代价泛函，该泛函衡量固体位移 ($u$) 和流体压力 ($p$) 与目标值 ($u_C, p_C$) 的偏差，同时对控制变量的大小进行惩罚。控制变量被确定为分布在区域内的流体源/汇项 ($m$)。本研究侧重于忽略加速度效应的准静态固结过程，并专门针对“孔隙弹性锁定”（poroelastic locking）相关的数学与数值挑战，这种现象发生在特定储水系数 ($s_0$) 趋于零时，导致流体与弹性介质趋于不可压缩。

方法论
作者提出并分析了一种新型双场对称格式，其主要未知量为固体位移和流体压力。该方法论通过以下阶段进行：

1. 连续形式化与适定性：

   • 作者通过对动量平衡方程求时间导数来重新构建控制方程，将系统转化为一组微分方程组。
   • 他们定义了适当的函数空间（用于位移的 $V$ 和用于压力的 $Q$），并利用鞍点理论和 Ladyzhenskaya–Babuška–Brezzi (LBB) 条件建立了弱解的适定性。
   • 利用变分法中的直接法，利用代价函数的凸性和状态算子的线性，证明了最优控制的存在性与唯一性。
   • 通过伴随状态法推导出一阶必要且充分的最优性条件，从而得到关于控制变量的变分不等式。

2. 离散方案：

   • 空间离散化： 使用满足 inf-sup 条件（例如 MINI 元或 Taylor-Hood 元）的有限元 ($V_h \times Q_h$) 对定义域进行离散化。
   • 时间离散化： 采用后向欧拉格式（具体为时间上的 dG(0) 方法），从而得到全离散系统。
   • 控制离散化： 本文对控制变量采用了变分离散化方法（由 Hinze 提出）。这意味着控制空间并未被显式离散化；相反，离散控制是通过将离散伴随状态投影到容许控制集上获得的。这避免了因离散化控制空间而引入的近似误差。

3. 误差分析：

   • 作者为全离散方案推导了先验误差估计。
   • 分析将总误差分为离散误差（针对固定控制）以及由于连续最优控制与离散最优控制之间的差异导致的扰动误差。
   • 建立了针对前向状态系统和后向伴随系统的稳定性估计，以限制这些误差。

核心贡献
本文对文献做出了三个主要贡献：

• 新型对称格式： 引入了一种专门为 Biot 模型中的最优控制问题定制的双场对称格式（位移与压力），这与以往的四场或混合格式不同。
• 存在性与唯一性： 严格证明了新格式的适定性，并确立了最优控制的存在性与唯一性，同时推导了一阶最优性条件。
• 先验误差估计： 为全离散方案提供了完整的先验误差分析。值得注意的是，作者推导了位移和压力在 $L^\infty(0, T; H^1(\Omega))$ 和 $L^2(0, T; L^2(\Omega))$ 范数下的误差，以及控制变量在 $L^2(0, T; L^2(\Omega))$ 范数下的误差，并在标准正则性假设下实现了最优收敛阶。

结果

• 理论方面： 本文证明了数值方案的收敛阶为 $O(h^s + \Delta t)$，其中 $h$ 是空间网格尺寸，$\Delta t$ 是时间步长，$s$ 与解的正则性相关。具体而言，如果解具有足够的正则性（$s=1$），则在空间和时间上均实现线性收敛（$O(h + \Delta t)$）。
• 数值方面： 在单位正方形定义域上的数值实验验证了理论发现。作者使用开源库 FEniCS，通过前向-后向算法求解了最优性系统。结果表明，当使用最低阶 MINI 元和 Taylor-Hood 元时，所提方案在位移、压力和控制变量上均达到了预期的最优收敛阶。

意义
作者将这项工作定位为研究受孔隙弹性系统控制的最优控制问题的基础性步骤。他们指出，虽然针对流体在可变形多孔介质中流动的最优控制已有研究（例如 Bociu 等人在储水系数消失的假设下进行的研究），但针对一般 Biot 模型的此类问题的数值近似此前尚未被研究过。通过建立一种新对称格式的适定性并为全离散方案提供严谨的误差估计，本文填补了文献空白，为生物组织建模和地球物理模拟等应用中的孔隙弹性过程控制提供了稳健的理论框架。