面向数据驱动计算力学中偏微分方程约束优化问题的变分多尺度方法

《COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING》：A variational multiscale approach to PDE-constrained optimization problems arising in data-driven computational mechanics

　　

研究针对数据驱动计算力学（DDCM, Data-Driven Computational Mechanics）中出现的偏微分方程（PDE, Partial Differential Equation）约束优化问题在反应-扩散情形下的原始形式与对偶形式最优性条件展开分析。首先，从连续框架出发，建立了这两类伴随表述的适定性（well-posedness）。随后，基于变分多尺度（VMS, Variational MultiScale）框架，针对其底层原始问题与对偶问题提出了稳定且相容的有限元逼近。对于拟一致（quasi-uniform）有限元剖分，研究分析了这些逼近的一般性质，并对两类典型子网格尺度（sub-grid scale）选择——即代数子网格尺度（ASGS, Algebraic Sub-Grid Scale）与正交子网格尺度（OSGS, Orthogonal Sub-Grid Scale）——建立了适定性。此外，在连续有限元函数情形下，仅通过改变稳定化参数的设计，便可在离散原始表述与离散对偶表述之间来回转换。最后，研究通过一系列复杂程度逐步提升的算例，对所提出的逼近格式进行了强度测试（stress-test），并给出了其数值性能的比较性与定性评估。

该文发表于《COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING》，研究主题聚焦于数据驱动计算力学（DDCM, Data-Driven Computational Mechanics）中反应-扩散问题所诱导的偏微分方程（PDE, Partial Differential Equation）约束优化子问题。文章的核心背景在于：DDCM作为连续介质力学中的新范式，不再预设显式本构关系，而是直接利用经验数据对场变量进行匹配。在反应-扩散问题中，目标是在满足守恒律、几何相容条件以及热力学第二定律的约束下，寻找与实验梯度-通量对最匹配的连续场。围绕这一框架，一个基础而关键的问题是：当离散数据场固定时，连续最优场是否存在、如何构造、以及采用何种数值格式才能稳定离散。既有工作虽然已在固体力学、结构动力学、电磁与接触等方向推动了DDCM的发展，但在扩散-反应情形下相关研究明显不足，且已有方法主要依赖交替方向法（ADM, Alternating Direction Method）反复求解优化子问题。本文正是针对这一连续优化子问题建立系统理论与稳定离散框架，从而弥补该方向的关键空白。

研究人员首先对原始的离散-连续非线性优化问题进行两层简化：其一，将热力学一致的辅助数据场 (e?, s?) 视为给定，从而把原问题化为带互补约束的数学规划；其二，在本文分析中暂不显式施加热力学第二定律约束，以便将问题线性化为鞍点问题，并集中考察其泛函分析结构。在此基础上，文章分别构造了连续层面的原始形式与对偶形式。原始形式以标量场 u、梯度场 e、通量场 s 及乘子场 λ, μ 为未知，建立混合鞍点变分问题；对偶形式则通过调整未知量所处函数空间，把标量变量置于 L₀²(Ω)，向量变量置于 H₀(div, Ω) 型空间中，形成另一组混合结构。文章证明：当无量纲物理参数 κζ²?_Ω⁴ 满足相应条件时，两种连续问题均满足混合鞍点问题的连续性、inf-sup稳定性与核上强制性，因此是适定的。原始形式的参数允许范围更宽，而对偶形式由于正则性转移，对参数区间要求更严格。

在离散层面，研究人员指出，若直接采用标准Galerkin有限元，则五场耦合系统需要满足复杂而空间相关的“小inf-sup条件”，不利于一般等阶连续有限元离散。为避免这种限制，文章引入变分多尺度（VMS）框架，将未知量分解为有限元可解析尺度与子网格尺度，并通过残差驱动方式对后者进行建模。围绕投影算子 P′ 的不同选取，研究建立了两类典型稳定化方法：若 P′ = I，得到ASGS格式；若 P′ = P^⊥，即投影到有限元空间的 L²(Ω) 正交补，则得到OSGS格式。文章分别为原始离散问题与对偶离散问题构造稳定双线性型，给出稳定化参数的尺度设计，并在拟一致网格与连续有限元空间下，严格证明了两类离散格式的适定性与相容性。尤其值得注意的是，由于连续有限元函数下原始与对偶双线性型本体一致，二者在离散层面的差别主要只体现在稳定化项及其参数设计上，因此可以仅通过改变稳定化参数结构，在原始与对偶离散表述之间切换，这为理解DDCM优化子问题的离散结构提供了统一视角。

本文采用的主要关键技术方法包括：第一，利用混合变分与鞍点理论，在连续泛函框架中分别分析原始形式与对偶形式的适定性；第二，基于VMS思想构造残差驱动的子网格尺度模型，系统推导ASGS与OSGS两类稳定化有限元格式；第三，在拟一致网格和连续等阶有限元空间上，通过逆不等式、Poincaré不等式与Young不等式建立离散inf-sup稳定性；第四，借助Firedrake有限元平台完成单体（monolithic）实现，并通过制造解算例与具有内部反应边界层的二维基准算例，对格式进行收敛性、稳定性与抗噪性能评估。数值样本均来自二维计算区域上的有限元离散解与其构造数据场，而非生物或临床样本队列。

以下结合论文主体结构，对研究结果进行分项解读。

2.1. Data-driven computational mechanics
本节首先系统回顾了DDCM在反应-扩散问题中的建模背景。研究人员将物理系统表示为浓度或温度 u、其梯度 e、通量 s 以及对偶乘子 λ, μ 的耦合体系，其中 λ 对应守恒律约束，μ 对应几何相容关系约束。通过引入经验数据集 D_m 中的梯度-通量对 (eⁿ, sⁿ)，构造辅助分片常值场 (e?, s?)，再以 L² 型距离度量连续场与数据场的偏差，建立原始的离散-连续非线性优化问题。本文进一步把问题简化为给定 (e?, s?) 与源项 q 时的PDE约束优化子问题，并指出：虽然文中为获得线性鞍点结构而暂不显式施加热力学第二定律约束，但该约束对于未来构造更稳健的优化算法具有重要意义。由此，本节得到全文研究对象，即一个线性化的PDE约束鞍点问题。

2.2. Primal form
本节推导并分析原始形式。研究人员选取 U=H₀¹(Ω)、E=[L²(Ω)]^d_Ω、S=[L²(Ω)]^d_Ω 等空间，将最优性条件写成混合系统，再进一步化简得到关于 u 与 λ 的最小不可约形式。研究表明，反应矩阵的特征值为 ±iκζ，说明在 ζ ≠ 0 时问题不能简单解耦，且具有四阶问题特征。随后文章定义原始问题的双线性型 a_p 与 b_p，建立对应的混合变分问题 P_p。通过调用既有混合问题理论并验证连续性和稳定性条件，研究证明：当 κζ²?_Ω⁴ ∈ (0, ∞) 时，原始连续问题适定，且解对数据 (e?, s?, q) 连续依赖。该结果说明原始形式在理论上具有较强鲁棒性。

2.3. Dual form
本节构造对偶形式，并将通量与乘子向量提升到 H₀(div, Ω) 空间。相对于原始形式，这一处理降低了标量场的正则性要求，但增强了向量场的散度控制能力。文章同样推导出其不可约强形式，指出该问题在本质上也表现为不可解耦的四阶型结构。随后，研究人员定义对偶双线性型 a_d 与 b_d，建立问题 P_d，并通过构造特定试探函数、利用散度算子的满射性、Young不等式与Poincaré不等式，证明对偶问题满足inf-sup条件与核上强制性。最终结论为：当 κζ²?_Ω⁴ ∈ (0, 2) 时，对偶连续问题适定。文章进一步指出，该更严格参数范围反映了原始-对偶转换中“标量正则性”与“向量正则性”之间的权衡。

3. Finite element approximation
本节进入有限元离散分析，是全文方法学核心。研究人员假定采用拟一致网格与连续有限元空间，并指出标准Galerkin格式会遭遇多重inf-sup兼容性困难，因此必须引入稳定化。VMS框架通过将未知量分解为有限元尺度与子网格尺度，并以方程残差显式建模后者，系统地产生稳定化项。文章给出了原始离散问题和对偶离散问题中五个子尺度分量 u′, e′, s′, λ′, μ′ 的统一表达，并据此构造稳定化双线性型和修正右端项。

3.1. Discretized primal problem
在原始离散问题中，研究人员为稳定化参数给出与网格尺度 h、物理参数 κ, ζ, ?_Ω 匹配的量纲一致设计，并分别对ASGS与OSGS情形证明适定性。对于ASGS，证明思路是选取特定测试函数组合，逐项利用Young不等式和Poincaré不等式控制各未知量及其导数；对于OSGS，则进一步借助投影分解与逆不等式，分别控制有限元空间内外的梯度与散度分量。由此，文章证明离散原始问题 P_p,h 适定且相容逼近连续问题。进一步地，研究还证明了若取特定参数选择，例如 k_u=k_λ=0 且 M_h=E_h，则可获得若干残量正交性，并推导出“最小稳定化”的原始格式；在这一情形下，ASGS与OSGS事实上等价。这一结果表明原始离散问题可在较简约的稳定结构下保持良好性质。

3.2. Discretized dual problem
对应地，文章为对偶离散问题构造了另一套稳定化参数设计，其中重点转向对 ?·s_h 与 ?·μ_h 的控制。与原始情形类似，研究分别证明ASGS与OSGS两种对偶离散格式均适定且相容。证明中除了控制散度项外，还需要通过额外的测试函数设计控制投影后的标量分量 P_Lh(u_h) 与 P_Uh(λ_h)。文章同样给出了特定参数下的残量正交性与最小稳定化结构，并说明在 κU_h=L_h 条件下，ASGS与OSGS也可视为等价。这部分结果从理论上确立了对偶离散方案的可行性。

3.3. Expected order of convergence
基于前述适定性与相容性，研究人员给出了误差行为的理论预期。对于原始格式，在定义的离散范数下，若所有未知量采用相同阶次 p 插值，则总体误差为 O(h^p)，其中 ?u 与 ?λ 具有 O(h^p) 精度，而 ?·s 与 ?·μ 为 O(h^p?1)。对偶格式则呈现相反特征：?u 与 ?λ 仅为 O(h^p?1)，但 ?·s 与 ?·μ 可达 O(h^p)。这清楚展示了从原始形式到对偶形式时，一个收敛阶从标量梯度转移到向量散度的现象，也是两类格式数值表现差异的理论根源。

4. Numerical results
本节通过数值实验系统检验理论。

4.1. Example 1: convergence to a manufactured solution
第一个算例采用制造解，在单位方形区域上构造光滑的标量场、梯度场、通量场及对偶场，并由此反推出数据场 e?, s? 与源项 q。研究分别在Dirichlet型与Neumann型边界条件下，对原始与对偶格式进行公平比较，并选取一组满足理论条件的稳定化参数。数值结果表明：两种格式在细网格上都能很好恢复全部场变量；但在粗网格和中等网格上，原始格式更擅长逼近标量场 u 与 λ，而对偶格式更擅长逼近向量乘子场 μ。在线性元情形，所有稳定格式都达到最优收敛率，原始格式在标量场误差上显著更优；在 H(div, Ω) 范数下，对偶格式对 μ 的收敛最为突出。二次元情形下，原始格式对标量场保持最优 H¹ 收敛，而对偶格式在 H(div, Ω) 范数下对向量场尤其是 μ 显示明显优势。与Natural格式、全稳定GLS格式和未稳定等阶Galerkin格式比较后，本文所提VMS稳定格式验证了其理论适定性与实现正确性。

4.2. Example 2: performance assessment in the presence of large gradients
第二个算例更具挑战性，考虑带圆形夹杂的二维反应-扩散问题，解具有明显内部反应边界层。研究先以标准Galerkin方法在同一网格上求得参考解，并将其按单元采样形成离散数据场，从而更接近真实DDCM应用场景。结果显示：在强梯度条件下，原始ASGS格式对标量势场 u 的逼近明显优于对偶格式，而对偶ASGS格式对通量场 s 的刻画更准确。沿对角线与不同半径圆周的剖面比较进一步确认了这一点，且误差分布图显示两类格式的优势区域与理论预期一致。与制造解算例相比，此时对偶格式在向量场上的优势不仅体现在 H(div, Ω) 范数中，而且在直接可视化结果上已十分明显。

Numerical results in the presence of noisy data
在同一反应边界层问题上，研究进一步向数据场 e? 与 s? 叠加均匀噪声，以模拟更实际的数据驱动环境。随着噪声水平 δ 增大，原始格式仍然在标量场恢复方面更稳健，对偶格式则持续在向量场恢复方面占优，即使在高噪声下这一趋势依然存在。与此同时，噪声破坏了 e 与 s 之间的热力学一致性，具体表现为部分区域出现 s·e > 0。文章强调，现有线性化框架本身并不阻止这类违反第二定律的现象，因此未来需要发展能够在噪声数据下兼顾误差过滤与热力学一致性的数值方法。研究还通过参数敏感性试验说明：在理论允许范围内，稳定化参数的适度变化对结果影响不大，表明所提格式对参数选择具有鲁棒性。

讨论部分可概括为：本文将DDCM中反应-扩散PDE约束优化子问题抽象为可分析的线性混合鞍点问题，并以VMS为基础构造了统一的原始/对偶稳定有限元框架。其理论价值在于同时建立了连续与离散层面的适定性，并阐明了原始与对偶表述之间的结构对应关系；其计算价值在于允许使用连续等阶有限元而不必满足复杂空间兼容条件。与此同时，文章也清楚指出研究的局限：为突出泛函分析主线，本文未显式纳入热力学第二定律与互补条件；此外，目前分析局限于连续有限元函数与拟一致网格。未来工作方向包括：将第二定律直接纳入优化格式，推广至非连续但相容的插值空间、非均匀网格、hp 收敛分析，以及更优的 L²(Ω) 误差估计。

研究结论部分可译为：本文提出了一种基于变分多尺度（VMS, Variational MultiScale）概念的有限元逼近方法，用于求解数据驱动计算力学（DDCM, Data-Driven Computational Mechanics）框架下反应-扩散问题中出现的偏微分方程（PDE, Partial Differential Equation）约束优化问题。研究证明了最优性条件的原始形式与对偶形式在连续层面均为适定问题，这是通过验证混合鞍点问题所需条件得以实现的。针对有限元逼近，为克服标准Galerkin方法对有限元空间复杂且问题依赖的inf-sup条件要求，研究采用VMS框架，通过系统建模子网格尺度，构造了稳定且相容的有限元格式，并进一步证明了代数子网格尺度（ASGS, Algebraic Sub-Grid Scale）与正交子网格尺度（OSGS, Orthogonal Sub-Grid Scale）两类方法的适定性。研究还表明，在VMS离散框架下，仅通过调整稳定化参数的设计，即可在问题的原始形式与对偶形式之间切换，从而为该类问题的离散结构提供统一认识。数值结果从定性与比较两个层面验证了所提稳定有限元格式的稳定性与相容性，显示出该方法在处理数据驱动模型固有挑战方面的有效性。尽管本文在应用于扩散-反应问题的DDCM方面取得了重要进展，研究仍承认若干简化假设：本文通过省略热力学第二定律的显式施加而对问题进行了线性化，但在计算后验上验证了解满足该规律；未来若能将其直接纳入优化格式，预计可获得更稳健的方法。在数值层面，当前工作聚焦于连续有限元函数与拟一致剖分，未来有必要进一步拓展到非连续空间、任意网格加密、hp 收敛以及改进的 L²(Ω) 收敛估计。