論文「任意多項式次数における離散最大値原理を持つ Cahn-Hilliard 方程式のための不連続ガラーキン法」の技術的概要

本論文は、Cahn-Hilliard (CH) 方程式、特に移動度（mobility）が退化する（degenerate）場合に対して、任意の多項式次数 $p$ において離散最大値原理（Discrete Maximum Principle, DMP）を満たす構造保存型の不連続ガラーキン（DG）法を提案するものです。著者らは、セル平均（cell averages）に基づいた移動度の評価と、Zhang-Shu スケーリングリミッターを組み合わせることで、厳密な退化性を保ちつつ、数値解が物理的な範囲 $[-1, 1]$ に収束することを証明しています。

以下に、問題設定、手法、主要な貢献、結果、および意義について詳細にまとめます。

1. 問題設定 (Problem)

Cahn-Hilliard 方程式は、二相流や相分離現象を記述する 4 階の非線形偏微分方程式です。本研究では、以下の系を扱います。

$$\begin{aligned} \partial_t \psi &= \text{Pe}^{-1} \nabla \cdot (M(\psi) \nabla \varphi), \\ \varphi &= W'(\psi) - \text{Cn}^2 \Delta \psi, \end{aligned}$$

ここで、$\psi$ は位相場（$[-1, 1]$ の範囲）、$\varphi$ は化学ポテンシャル、$M(\psi) = \max\{1-\psi^2, 0\}$ は退化移動度です。

• 課題: 4 階方程式は 2 階方程式（熱方程式など）とは異なり、強い最大値原理を持たないため、数値離散化において解が物理的に許容される範囲 $[-1, 1]$ を逸脱する（オーバーシュート/アンダーシュート）リスクがあります。
• 既存手法の限界: 従来の有限要素法（FEM）や対称内部ペナルティガラーキン（SIPG）法では、離散最大値原理が保証されず、数値実験でも範囲逸脱が観測されます。また、既存の制限付き（limited）手法（FEM-L, SIPG-L など）は、セル平均自体が最大値原理を満たすことを保証する理論的根拠が不足しており、特に移動度の評価方法（多項式評価）に正則化（regularization）を必要とするため、厳密な退化性を損なう可能性があります。

2. 手法 (Methodology)

著者らは、SWIPDP-L 法（Symmetric Weighted Interior Penalty Diffusive Projection with Scaling Limiter）を提案しました。

2.1. 空間離散化：SWIP 定式化

• セル平均に基づく移動度評価: 移動度 $M(\psi)$ を要素内の多項式として評価するのではなく、セル平均 $\Pi_0 \psi$ において評価します。
  • これにより、セル平均が $\pm 1$ に達した時点で移動度が厳密に 0 となり、そのセル境界でのフラックスが 0 になります。
  • この特性を利用することで、正則化（$\delta > 0$ を加えるなど）なしに、移動度の退化性を厳密に保持できます。
• 対称重み付き内部ペナルティ（SWIP）: 拡散項の離散化に SWIP 形式を採用し、ペナルティパラメータ $\eta_e$ を移動度に依存しない十分に大きな値に設定することで、強制性（coercivity）を証明しています。

2.2. 時間離散化と非線形処理

• Eyre スプリッティング: 非線形ポテンシャル $W'(\psi)$ を凸部分と凹部分に分割し、陰的・陽的に処理することで、エネルギー安定性を確保します。
• Leray-Schauder 不動点定理: 退化移動度を含む非線形方程式の解の存在を、正則化された問題から退化極限へ通過させることで証明しました。

2.3. 離散最大値原理の達成：Zhang-Shu スケーリングリミッター

• セル平均の有界性: 上記の SWIP 定式化により、セル平均 $\Pi_0 \psi$ が $[-1, 1]$ に収束することを証明（定理 4.5）。
• スケーリングリミッター: 要素内の解が $[-1, 1]$ を逸脱した場合、Zhang-Shu スケーリングリミッターを適用して点ごとの値を補正します。
  • セル平均が既に有界であることが保証されているため、このリミッターを適用することで、最終的な解 $\psi_h$ も $L^\infty$ 有界性を満たします。
  • このプロセスは、セル平均を調整する必要がある既存の 2 段階リミッター（例：[27]）に比べ、単一ステップで効率的です。

3. 主要な貢献 (Key Contributions)

1. 任意次数 $p$ における離散最大値原理の証明:
   • 従来の手法は $p=0$ に限定されるか、証明が不十分でしたが、本手法は任意の多項式次数 $p \ge 0$ でセル平均レベルでの有界性を証明し、リミッターを介して点ごとの有界性を保証します。
2. 正則化不要の退化移動度処理:
   • 移動度をセル平均で評価することで、正則化パラメータ $\delta$ を必要とせず、移動度の厳密な退化性（$M(\psi)=0$）を保持します。これにより、強制定数が移動度に依存せず、数値的安定性が向上します。
3. 解の存在証明:
   • 退化移動度を含む場合でも、Leray-Schauder 不動点定理を用いて離散解の存在を証明しました。
4. 保存則とエネルギー安定性:
   • 質量保存: 無制限版および制限版（リミッター適用後）の両方で質量保存が成立します（リミッターはセル平均を変化させないため）。
   • エネルギー減衰: $p=0$ の場合、単調なエネルギー減衰が理論的に証明されます。$p \ge 1$ の場合、リミッターの適用により理論的な単調性は保証されませんが、数値実験では常に減衰が観測されています。

4. 数値結果 (Numerical Results)

Dune-Fem ソフトウェアを用いた数値実験により、以下の点が確認されました。

• 収束性:
  • $p=0$ で $L^2$ ノルムにおいて 1 次収束、$p=1$ で $L^2$ ノルムにおいて 2 次収束、$H^1$ ノルムにおいて 1 次収束を示し、理論的な最適収束率と一致しました。
• 最大値原理の検証:
  • 数値計算における浮動小数点誤差や非線形ソルバの許容誤差（TOL）により、厳密な境界値 $1$をわずかに超えることがありますが、その誤差は機械精度（$\sim 10^{-14}$）のオーダーに収まり、理論的な構造を破るものではありません。
  • $p=0$ と $p=1$ の両方で、解が $[-1, 1]$ の範囲内に強く拘束され、振動しないことが確認されました。
• エネルギーと質量:
  • 時間発展に伴いエネルギーが単調に減少し、質量が保存されていることが確認されました。

5. 意義と展望 (Significance and Outlook)

• 理論的厳密性と実用性の両立: 構造保存（構造を壊さない）数値解法において、離散最大値原理を任意次数で保証する最初の DG 法のひとつであり、かつ実用的な効率性（単一リミッター適用）を持っています。
• 正則化の排除: 従来の手法が抱えていた「正則化パラメータの選択による誤差」や「強制定数のデータ依存性」といった問題を解決し、より物理的に忠実なシミュレーションを可能にします。
• 将来の展開:
  • $p \ge 1$ におけるエネルギー減衰の理論的証明（リミッターがエネルギーを減少させる条件の特定など）。
  • 対流項を含む Cahn-Hilliard 方程式への拡張（アップウィンド法との組み合わせ）。
  • Cahn-Hilliard-Navier-Stokes 連成系や、$hp$ 適応メッシュへの適用。

結論

本論文は、Cahn-Hilliard 方程式の数値解法において、**「任意次数での離散最大値原理」「質量保存」「エネルギー安定性」「正則化不要」**という、従来はトレードオフとされてきた複数の要件を同時に満たす革新的な DG 法を提案し、その数学的基礎と数値的有効性を確立した重要な研究です。