Stable $C^1$-conforming finite element methods for a class of nonlinear fourth-order evolution equations

技術要約：一連の非線形4階発展方程式に対する安定な$C^1$適合有限要素法

問題定式化
本論文は、有界領域 $\Omega \subset \mathbb{R}^d$ ($d \in \{1, 2, 3\}$) において定義される、一般的なクラスの非線形4階発展偏微分方程式（PDE）の数値解析を取り扱う。支配方程式は以下のように与えられる：
$$\frac{\partial u}{\partial t} - \beta_1 \Delta u + \beta_2 \Delta^2 u = \beta_3(1 - |u|^2)u - \beta_4 u \times \Delta u + \beta_5 \Delta(|u|^2 u) + \beta_6 (j \cdot \nabla)u$$
ただし、ノイマン境界条件（$\frac{\partial u}{\partial n} = 0, \frac{\partial (\Delta u)}{\partial n} = 0$）および初期条件 $u(0) = u_0$ を満たすものとする。ここで、$u$ はスカラー値（$m=1$）またはベクトル値（$m=3$）であり得る。この方程式は、係数 $\beta_i$ および次元 $d$ の設定に応じて、以下のような様々な物理モデルを包含している：

• マイクロマグネティクス（微磁性学）： ランドー・リフシッツ・バリアクラー（LLB）方程式（スピントルクの有無を含む）、およびランドー・リフシッツ・ブロッホ方程式の正則化。
• パターン形成： スウィフト・ホーエンバーグ方程式。
• 相分離および拡散： ソース項を伴う対流型カーン・ヒリアン方程式、異常バイフラックス拡散モデル、および一般化された個体群分散モデル。

著者らは、特定の例（例：スカラーのスウィフト・ホーエンバーグや標準的なカーン・ヒリアン）に対しては数値的手法が存在するものの、すべての係数が潜在的に非ゼロである可能性のある完全なシステム (1.1) に対する統一的かつ一般的な数値スキームは、文献において欠如していると指摘している。

手法
著者らは、$C^1$適合要素を利用した完全離散有限要素法を提案し、分析している。この選択は、4階の空間微分（$\Delta^2 u$）によって必要とされる、有限要素空間が $H^2(\Omega)$ の部分空間であるという要件に基づいている。

• 空間離散化： 本手法は、次数 $r$ の区分的多項式を用いる（$d=1,2$ では $r=3$、$d=3$ では $C^1$ 連続性を維持するために $r=9$）。言及されている具体的な要素には、$C^1$-P3 ヘルミット要素（1D）、Hsieh–Clough–Tocher (HCT) マクロ要素（2D）、および Zenísek 要素（3D）がある。
• 時間離散化： 2つの線形化時間ステップスキームが提案されている：
  1. 半黙示型オイラー法（時間に関して1次）。
  2. 半黙示型 BDF2 法（時間に関して2次）。
     両方のスキームにおいて、非線形項は前の時刻ステップで評価されることで線形化され、その結果、各ステップで線形方程式の系が解かれる。これにより、ニュートン法による非線形系の求解に伴う計算コストを回避している。

主要な貢献および理論的結果
本論文は、提案されたスキームに関する厳密な誤差解析と安定性の証明を提供している：

1. 安定性：

   • これらのスキームは、次元 $d=1$ および $d=2$ において、$H^2$ ノルムで無条件に安定であることが証明されている。これは、メッシュサイズ $h$ に対する時間ステップサイズ $k$ の制限や、準一様（quasi-uniform）な三角形分割の仮定なしに安定性が保持されることを意味する。
   2. $d=3$ の場合、スキームは条件付き安定である。安定性は、初期データが十分に小さい場合、または時間ステップが初期エネルギーに対して特定の制限を満たす場合に保証される。著者らは、$C^1$ 適合要素（次数9の多項式を必要とする）は計算量的に過大なため、3Dにおいては非適合（non-conforming）手法の方が実用的であることを認めているが、理論的解析は完全性のために含まれている。

2. 誤差推定と収束：

   • 著者らは、$\beta = 0, 1, 2$ における $L^\infty(0, T; H^\beta)$ ノルムにおける最適な収束率を確立している。
   • 半離散スキームにおいて、誤差は $O(h^{r+1-\beta})$ で抑えられる。
   • 完全離散スキームについては：
     • 半黙示型オイラー法は $O(h^{r+1-\beta} + k)$ の誤差を達成する。
     • 半黙示型 BDF2 法は $O(h^{r+1-\beta} + k^2)$ の誤差を達成する。
   • この解析は、非線形問題に適応した特殊な楕円射影の開発、およびソボレフ埋め込みと逆不等式を用いた非線形項の慎重な推定に基づいている。

3. 数値実験：

   • 理論的な収束率は、2Dにおける HCT 要素を用いたオープンソースパッケージ FreeFEM を用いて検証された。
   • 実験は、ベクトル値のランドー・リフシッツ・バリアクラー方程式、およびソース項を伴うスカラー値のカーン・ヒリアン方程式に対して行われた。
   • 結果は、予測された収束次数を確認した。すなわち、3次多項式（$r=3$）の要素の使用と整合して、$L^2$ における4次、$H^1$ における3次、および $H^2$ における2次の空間離散化精度が得られた。

意義および主張
本論文は、一般的な非線形4階システム (1.1) の解を近似できる最初の数値スキームを提供していると主張している。強調されている主な利点は以下の通りである：

• 線形性： 非線形項を線形化することで、各時間ステップでの非線形系の求解に伴う複雑さや潜在的な不安定性を回避しており、微磁性シミュレーションのような用途において有利である。
• 安定性： $d=1,2$ における $H^2$ での無条件安定性は、厳格な時間ステップの制約なしに堅牢な性能を保証する。
• 精度： $C^1$ 適合要素の使用により、非適合手法と比較して、計算コストが増大するものの、解の勾配を捉えるための高次の精度が可能となる。

著者らは3Dのケースに関して控えめなトーンを維持しており、$C^1$ 適合要素の3Dにおける実用的な適用は、要求される多項式の次数が高いため、現在は計算量的に不可能であることを明示している。本研究は、物理学、材料科学、および生物学から生じる複雑な非線形4階発展方程式を解くための、基礎的な理論的枠組みとして機能する。