[论文解读] On the numerical approximation of $p$-Biharmonic and $\infty$-Biharmonic functions

Nikos Katzourakis, Tristan Pryer|arXiv (Cornell University)|Jan 25, 2017

Advanced Mathematical Modeling in Engineering参考文献 3被引用 3

一句话总结

该论文提出了一种C⁰混合有限元方法，用于数值逼近p-双调和函数与∞-双调和函数，证明了离散解收敛于p-双调和方程的弱解，并严格地取p→∞的极限以捕捉∞-双调和解。该方法揭示，即使在简单的边界条件下，∞-双调和函数也是分片二次的，且具有复杂的界面模式。

ABSTRACT

In [KP16] (arXiv:1605.07880) the authors introduced a second-order variational problem in $L^{\infty}$. The associated equation, coined the $\infty$-Bilaplacian, is a \emph{third order} fully nonlinear PDE given by $Δ^2_\infty u\, := (Δu)^3 | D (Δu) |^2 = 0.$ In this work we build a numerical method aimed at quantifying the nature of solutions to this problem which we call $\infty$-Biharmonic functions. For fixed $p$ we design a mixed finite element scheme for the pre-limiting equation, the $p$-Bilaplacian $Δ^2_p u\, := Δ(| Δu |^{p-2} Δu) = 0.$ We prove convergence of the numerical solution to the weak solution of $Δ^2_p u = 0$ and show that we are able to pass to the limit $p o\infty$. We perform various tests aimed at understanding the nature of solutions of $Δ^2_\infty u$ and in 1-$d$ we prove convergence of our discretisation to an appropriate weak solution concept of this problem, that of $\mathcal D$-solutions.

研究动机与目标

• 在最小正则性假设下，为四阶完全非线性PDE——p-双调和方程，设计一种稳定且收敛的数值格式。
• 建立混合有限元方法对p-双调和问题弱解的收敛性，并在附加正则性条件下获得最优收敛速率。
• 证明p-双调和方程在p→∞时的极限可合理地解释为∞-双调和方程，后者为非散度形式的三阶完全非线性PDE。
• 通过数值方法研究∞-双调和函数的结构特征，特别是其分片二次性质与复杂的界面模式。
• 针对缺乏经典解的∞-双调和方程，定义并验证一种弱解概念——D-解。

提出的方法

• 通过引入Hessian与Laplacian的辅助变量，构建混合有限元格式，使得可使用C⁰有限元而非C¹有限元。
• 该方法基于最小化p-双调和能量泛函的变分形式，将其改写为混合形式以解耦四阶PDE。
• 受先前研究启发的inf-sup条件确保了稳定性，并在更高正则性下支持最优收敛速率的分析。
• 采用Galerkin方法求解离散系统，使用分片多项式逼近，通过紧致性与Young测度技术证明收敛性。
• 利用L∞中变分法理论分析p→∞的极限，得出∞-双调和方程为p-双调和PDE的正式极限。
• 在二维区域上进行数值实验，通过改变边界数据中的参数m，观察随着p增大，Δu趋近于分片常数函数。

实验结果

研究问题

• RQ1 如何在不依赖C¹有限元的条件下，为p-双调和方程设计一种稳定且收敛的有限元方法？
• RQ2 数值解对p-双调和方程的收敛行为如何？在正则性假设下能否建立最优收敛速率？
• RQ3 p-双调和方程在p→∞时的极限能否在数值与分析上均合理地解释为∞-双调和方程？
• RQ4 ∞-双调和函数表现出何种结构性质，特别是关于分片正则性与界面复杂性？
• RQ5 对于缺乏经典解的∞-双调和方程，如何定义并数值验证一种弱解概念（D-解）？

主要发现

• 在最小正则性假设下，混合有限元方法收敛于p-双调和方程的弱解，收敛性通过紧致性与Young测度技术得到证明。
• 在附加正则性条件下，该方法实现了依赖于p值的最优收敛速率，且随着p增大，稳定性与精度均得到提升。
• 随着p增大，p-双调和方程的数值解趋近于∞-双调和极限，其中Δu趋于分片常数函数，表明解满足具有分片常数右端项的Poisson方程。
• 在边界数据中m=1,2,3的数值模拟结果表明，∞-双调和函数在整个区域上为分片二次函数，且界面模式高度复杂且非平凡。
• 该方法成功捕捉了∞-双调和方程的D-解概念，验证了弱形式的有效性，尽管经典解不存在。
• 数值测试中观察到的收敛速率超过理论预测，提示可通过拟范数稳定化实现进一步改进，这与低阶问题中的已知行为一致。