[论文解读] On the Numerical Approximation of $\infty$-Harmonic Mappings

Nikos Katzourakis, Tristan Pryer|arXiv (Cornell University)|Nov 4, 2015

Spectral Theory in Mathematical Physics参考文献 10被引用 5

一句话总结

本文提出了一套数值框架，用于近似求解向量型 ∞-Laplacian 系统，该系统是 L∞ 中极值变分法中出现的退化椭圆型偏微分方程。通过在单位正方形上使用有限元方法处理秩一边界数据，作者展示了 ∞-调和映射中的相分离与界面形成现象，揭示了当 p→∞ 时的结构转变，其结果与奇异解和非光滑梯度行为的理论预测一致。

ABSTRACT

Given a map $u : Ω\subseteq \mathbb{R}^n \longrightarrow \mathbb{R}^N$, the $\infty$-Laplacian is the system \[ \label{1} Δ_\infty u \, :=\, \Big( ext{D}u \otimes ext{D}u + | ext{D}u|^2 [ ext{D}u]^\bot \! \otimes I \Big) : ext{D}^2 u\, = \, 0. ag{1} \] \eqref{1} is the model system of vectorial Calculus of Variations in $L^\infty$ and arises as the "Euler-Lagrange" equation in relation to the supremal functional \[ \label{2} E_\infty(u,Ω)\, :=\, \| ext{D}u \|_{L^\infty(Ω)}. ag{2} \] The scalar case of \eqref{1} has been introduced by Aronsson in the 1960s and by now is relatively classical and well understood. The general system \eqref{1} has been discovered and studied by the first author in a series of recent papers. Supremal functionals are fundamental for applications because they provide more realistic models as opposed to conventional integral models. Herein we provide numerical approximations of solutions to the Dirichlet problem when $n=2$ and $N=2,3$ for certain carefully selected boundary data on the unit square. Our experiments demonstrate interesting and unexpected phenomena occurring in $L^\infty$ and provide insights on the structure of general solutions and the natural separation to phases they present.

研究动机与目标

• 开发并实现一种用于在二维和三维中求解向量型 ∞-Laplacian 系统的数值方法。
• 研究在非光滑、秩一边界数据下的 Dirichlet 问题中 ∞-调和映射解的结构。
• 通过数值实验探索 L∞-变分问题中相分离与界面的出现机制。
• 验证 p-Laplacian 解在 p→∞ 时收敛于 ∞-调和极限的性质。

提出的方法

• 在单位正方形区域 Ω⊂ℝ² 上，使用有限元法（FEM）求解向量型 p-Laplacian 系统的数值解。
• 对不断增加的 p 值（最高达 1024）求解 p-Laplacian 系统，以近似 ∞-调和映射（当 p→∞ 时）。
• 通过光滑过渡函数 K(t) 的积分变换构造边界数据，确保梯度具有秩一结构。
• 计算并可视化梯度矩阵的行列式 det(DU)，以检测秩降低与界面形成。
• 使用间隔为 0.05 的 det(DU) 等高线图，追踪不同 p 值下解结构的演化过程。
• 该方法利用粘性解理论，并处理了 ∞-Laplacian 算子中由正交投影 [Du]⊥ 引入的不连续系数。

实验结果

研究问题

• RQ1 在二维和三维中，向量型 p-Laplacian 的解如何随着 p→∞ 收敛于 ∞-调和映射？
• RQ2 在秩一边界条件下，∞-调和映射中会涌现出哪些结构性特征（如界面或相分离）？
• RQ3 在 ∞-调和解中，梯度矩阵的秩如何在空间上演化？这对其解的正则性有何含义？
• RQ4 数值实验能否揭示分析构造中预测的三重节点界面或矩形界面？
• RQ5 不连续的投影 [Du]⊥ 在 ∞-调和映射的数值逼近中起到何种作用？

主要发现

• 随着 p 增大，向量型 p-Laplacian 的解显现出梯度秩降低的区域，表明相分离结构正在形成。
• 当 p≥64 时，梯度矩阵的行列式 det(DU) 在域的大部分区域显著减小，表明秩降低并形成奇异相。
• 当 p=1024 时，数值解清晰显示出三重节点界面，与具有此类特征的显式 ∞-调和映射的分析示例完全一致。
• 随着 p 增大，矩形界面结构愈发显著，证实了 p-Laplacian 解向具有分段常数梯度结构的 ∞-调和解收敛。
• 数值实验确认，即使对于光滑映射，∞-调和映射的 PDE 中也可能出现不连续系数，这是由正交投影 [Du]⊥ 引起的。
• 结果表明，向量型 ∞-Laplacian 支持非光滑、奇异的解，具有明显的相分离结构，与粘性解理论及 L∞-变分问题的理论一致。