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

ひとこと要約

本稿では、p-双曲微分方程式および∞-双曲微分方程式の数値的近似のためのC⁰混合有限要素法を開発し、離散解がp-双曲微分方程式の弱解に収束することを証明するとともに、p→∞の極限を厳密に取り入れ、∞-双曲微分解を捉える。この手法により、∞-双曲微分関数が単純な境界条件のもとでも複雑な界面パターンを持つ区分的2次関数であることが明らかになった。

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.

研究の動機と目的

• 最小限の正則性仮定のもとで、4階の完全非線形PDEであるp-双曲微分方程式の安定的かつ収束性を示す数値スキームの開発。
• 混合有限要素法がp-双曲微分問題の弱解に収束することを確立し、追加の正則性のもとで最適収束率を達成すること。
• p-双曲微分方程式のp→∞極限が、発散型でない3階完全非線形PDEとしての∞-双曲微分方程式に正当化できることを示すこと。
• ∞-双曲微分関数の構造的性質、特に区分的正則性および界面の複雑さを数値的に調査すること。
• 古典的解が存在しないことを踏まえ、∞-双曲微分方程式に対する弱解概念（D解）を定義し、検証すること。

提案手法

• ヘッセ行列とラプラシアンの補助変数を導入することで、C¹有限要素ではなくC⁰有限要素を用いることが可能となる混合有限要素式を構築。
• p-双曲微分エネルギー汎関数を変分的定式化し、4階PDEを分解する混合形に書き直す。
• 先行研究にインspiredされたインフラ・サブ条件が安定性を保証し、より高い正則性のもとで最適収束率を達成する収束解析を可能にする。
• 離散系は、区分的多項式近似を用いたガレルキン法により解かれ、コンパクトネスおよびヤング測度技法を用いて収束性が証明される。
• p→∞極限は、L∞における変分法の理論を用いて分析され、p-双曲微分PDEの形式的極限として∞-双曲微分方程式が得られる。
• 境界データのmが異なる2次元領域における数値実験により、pが増加するにつれて、Δuが区分的定数関数に近づく様子が観察された。

実験結果

リサーチクエスチョン

• RQ1 C¹有限要素を必要としない安定的かつ収束性を示す有限要素法を、p-双曲微分方程式に対してどのように設計できるか？
• RQ2 数値解の収束挙動はどのようなものか？正則性仮定のもとで最適収束率を確立できるか？
• RQ3 p-双曲微分方程式のp→∞極限が、数値的および解析的に∞-双曲微分方程式として正当化できるか？
• RQ4 ∞-双曲微分関数はどのような構造的性質を示すか？特に、区分的正則性および界面の複雑さの面で。
• RQ5 古典的解が存在しない∞-双曲微分方程式に対して、弱解概念（D解）をどのように定義し、数値的に検証できるか？

主な発見

• 混合有限要素法は、最小限の正則性仮定のもとで、p-双曲微分方程式の弱解に収束し、コンパクトネスおよびヤング測度を用いた手法により収束性が証明された。
• 追加の正則性のもとで、収束率が最適となり、その値はpの値に依存する。pが大きいほど安定性および精度が向上する。
• pが増加するにつれて、p-双曲微分方程式の数値解は∞-双曲微分極限に近づき、Δuが区分的定数関数に近づく傾向を示し、解が右辺が区分的定数関数であるポアソン方程式を満たしていることが示唆された。
• 境界データのm=1,2,3における数値シミュレーションにより、∞-双曲微分関数が領域全体で区分的2次関数であり、極めて複雑で非自明な界面パターンを示すことが明らかになった。
• 本手法は、古典的解が存在しない∞-双曲微分方程式のD解概念を効果的に捉え、弱定式化の妥当性を検証した。
• 数値実験での観察された収束率は理論的予測を上回っており、準ノルム安定化による改善の可能性が示唆され、低階数問題における既知の挙動と整合的である。