[論文レビュー] 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

ひとこと要約

本稿では、L∞における変分法の極値的計算から生じる非退化楕円型偏微分方程式であるベクトル値∞-ラプラシアン系の数値的近似フレームワークを提示する。単位正方形上での有限要素法を用い、ランク1の境界データを適用することで、∞-調和写像における相分離および界面形成を示し、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.

研究の動機と目的

• 2次元および3次元におけるベクトル値∞-ラプラシアン系を解くための数値的手法の開発および実装。
• 滑らかでないランク1の境界データを有する∞-調和写像のディリクレ問題の解の構造の調査。
• 数値実験を通じて、L∞-変分問題における相分離および界面の出現の解明。
• p→∞におけるp-ラプラシアン解の∞-調和極限への収束の検証。

提案手法

• 単位正方形領域 Ω⊂ℝ² におけるベクトル値p-ラプラシアンの有限要素法（FEM）を用いた数値的解法。
• pの値を増加させ（最大1024まで）、p→∞における∞-調和写像の近似を実行。
• 滑らかな遷移関数 K(t) の積分変換を用いて境界データを構築し、ランク1の勾配構造を保証。
• 勾配行列の行列式 det(DU) を計算および可視化し、ランクの低下および界面形成を検出。
• det(DU) の等高線図（0.05刻み）を用いて、p値の変化に伴う解の構造の進化を追跡。
• 粘性解理論を活用し、∞-ラプラシアン作用素に由来する不連続係数（直交射影 [Du]⊥）を処理。

実験結果

リサーチクエスチョン

• RQ1 2次元および3次元において、ベクトル値p-ラプラシアンの解はp→∞に伴いどのように∞-調和写像に収束するか？
• RQ2 ランク1の境界条件下で、∞-調和写像に現れる界面や相分離といった構造的特徴は何か？
• RQ3 ∞-調和解において、勾配行列のランクは空間的にどのように変化するか？これは解の正則性にどのような意味を持つのか？
• RQ4 数値実験により、解析的構成で予測された三重接点または長方形型界面が実際に観測されるか？
• RQ5 不連続な射影 [Du]⊥ が、∞-調和写像の数値的近似に果たす役割は何か？

主な発見

• pが増加するにつれ、ベクトル値p-ラプラシアンの解は勾配ランクが低い領域が拡大し、相分離構造の出現を示す。
• p≥64 において、領域の大部分で勾配行列の行列式 det(DU) が著しく低下し、ランクの低下および特異相の形成を示唆する。
• p=1024 の場合、数値解は明確に三重接点界面を示しており、このような特徴を持つ明示的∞-調和写像の解析的例と一致する。
• pの増加に伴い、長方形型界面構造が顕著に現れ、p-ラプラシアン解が勾配構造が区分的定数である∞-調和解に収束することを確認する。
• 数値実験により、滑らかな写像に対しても、∞-調和写像のPDEに不連続係数が現れることを確認した。
• 結果は、ベクトル値∞-ラプラシアンが滑らかでない特異解および明確な相を有することを示しており、粘性解理論およびL∞-変分問題の理論と整合的である。