[论文解读] Adjoint methods for obstacle problems and weakly coupled systems of PDE

Filippo Cagnetti, Diogo A. Gomes+1 more|arXiv (Cornell University)|Mar 16, 2011

Numerical methods for differential equations参考文献 12被引用 5

一句话总结

本文将伴随方法应用于障碍问题及哈密顿-雅可比方程的弱耦合系统惩罚逼近的收敛速率分析。通过推导解关于正则化参数 ε 的导数的统一有界性，即使在哈密顿量不具凸性的情况下，也建立了 L∞ 范数下的 O(ε^{1/2}) 收敛速率，从而将伴随方法的应用范围扩展至非凸及退化 PDE。

ABSTRACT

The adjoint method, recently introduced by Evans, is used to study obstacle problems, weakly coupled systems, cell problems for weakly coupled systems of Hamilton-Jacobi equations, and weakly coupled systems of obstacle type. In particular, new results about the speed of convergence of common approximation procedures are derived.

研究动机与目标

• 分析障碍问题及哈密顿-雅可比方程弱耦合系统的惩罚逼近的收敛速率。
• 将伴随方法——此前仅用于凸哈密顿量与弱 KAM 理论——推广至非凸及退化 PDE。
• 推导解对正则化参数 ε 的敏感性的定量估计。
• 控制来自障碍约束、随机切换及弱耦合系统中最佳切换的奇异项。
• 建立一个研究具有切换动力学与非光滑约束系统的收敛性的通用框架。

提出的方法

• 将伴随方法应用于惩罚系统：在 U 内有 u^ε + H(x, Du^ε) + γ^ε(u^ε - ψ) = εΔu^ε，且在 ∂U 上有 u^ε = 0。
• 推导惩罚 PDE 的形式线性化，并构造相应的伴随方程，以估计对 ε 的敏感性。
• 通过分部积分及对伴随解 σ^j,ε 的估计，建立 ∂u^ε/∂ε 的统一有界性。
• 利用最大值原理控制涉及 (γ^ε)' 和 γ^ε_ε 的项，特别是靠近障碍界面处的项，推导关键估计。
• 通过分析具有切换哈密顿量和障碍型约束的耦合方程，将方法扩展至弱耦合系统。
• 证明策略依赖于将 ∂u^ε/∂ε 的 L∞ 范数有界于 C/ε^{1/2}，从而得出 O(ε^{1/2}) 的收敛速率。

实验结果

研究问题

• RQ1 当哈密顿量非凸时，惩罚解 u^ε 到障碍问题的粘性解 u 的收敛速率是什么？
• RQ2 如何将伴随方法调整以控制弱耦合系统中切换动力学引起的奇异项？
• RQ3 是否可以在不假设哈密顿量 H 在动量变量上凸性的情况下，建立 O(ε^{1/2}) 的收敛速率？
• RQ4 伴随解 σ^j,ε 在系统存在不连续性或非光滑切换时的行为如何？
• RQ5 惩罚参数 ε 在控制障碍型 PDE 解的正则性与敏感性方面起什么作用？

主要发现

• 本文在 L∞ 范数下建立了障碍问题惩罚逼近的 O(ε^{1/2}) 收敛速率，即使不假设哈密顿量的凸性。
• 推导出解对 ε 的导数的统一有界性：||∂u^ε/∂ε||_{L∞} ≤ C/ε^{1/2}，这直接意味着 O(ε^{1/2}) 的收敛速率。
• 伴随方法通过仔细分析伴随方程，成功控制了来自障碍约束的奇异项，如 (γ^ε)'(u^ε - ψ)。
• 该方法被扩展至具有障碍型约束的哈密顿-雅可比方程弱耦合系统，得到类似的收敛估计。
• 该分析提供了一个稳健的框架，用于研究具有非凸或非光滑哈密顿量的退化 PDE 的收敛性。
• 结果表明，伴随方法在处理非凸及切换动力学的 PDE 系统中非常有效，为弱 KAM 理论与最优控制提供了新工具。