A family of second order, linear, unconditionally stable methods for the Cahn-Hilliard-Navier-Stokes equations

技術要約：Cahn-Hilliard-Navier-Stokes方程式に対する、2次精度・線形・無条件安定な一族の手法

問題設定
複雑な多相流体力学、特にCahn-Hilliard-Navier-Stokes (CHNS) 方程式によってモデル化される現象の正確な数値シミュレーションは、重大な課題を伴う。これは、4階のCahn-Hilliard演算子の数学的な剛性と、非圧縮Navier-Stokes方程式との非線形な結合に起因する。標準的な陽解法（explicit time-stepping schemes）は、安定性を維持するために非常に厳しい時間ステップ制限（例：$\Delta t \sim O(\Delta x^4)$）に直面する。完全陰解法や凸分割法（convex-splitting methods）は無条件のエネルギー安定性を提供する一方で、各タイムステップで大規模かつ結合した非線形系を解く必要があり、計算コストが増大する場合が多い。したがって、物理的な保存特性を犠牲にすることなく、2次の時間精度、無条件安定性、および線形な可解性を組み合わせたアルゴリズムが求められている。

手法
著者らは、密度が一致した二相流のための、2次精度・線形・無条件安定なIMEX（implicit-explicit）法のファミリーを提案している。この半離散スキームは、主に以下の3つの要素で構成されている：

1. 非線形項の外挿： 非線形な対流項および結合項は、外挿された量（$H_{n+\theta}$ で表記）を用いて陽的に扱われる。これにより、非線形反復計算を回避している。
2. 補助変数定式化： 非線形な自由エネルギー項は、補助変数 $q$（二重井戸ポテンシャル $F(\phi)$ から導出）を介して処理され、問題の線形性を維持するように再定式化される。
3. 時間的曲率正則化： パラメータ $\epsilon \geq 0$ によって制御される安定化メカニメントが導入されている。これには、Navier-Stokesの正則化に関する先行研究から適応された概念である、解の離散的な時間的曲率を組み込んだ特定の補間演算子 $J^\epsilon_{n+\theta}$ が含まれる。

結果として得られるアルゴリズム（Algorithm 2.1）は、各タイムステップにおいて、定数係数を持つ線形系を解くことのみを要求する。このスキームは、パラメータ $\theta \in (1/2, 1]$ および $\epsilon \geq 0$ によってパラメータ化される。

主な貢献と理論的結果

• 無条件安定性： 本論文では、提案されたスキームが $\theta \in (1/2, 1]$ かつ $\epsilon \geq 0$ において無条件に長期安定であることを示す厳密な離散エネルギー評価を確立している。証明には、離散ノルムを定義するための対称正定値行列が用いられ、本スキームがタイムステップの制限を受けることなく、離散エネルギー散逸則を満たすことを示している。
• 線形可解性： 非線形ソルバーを必要とする可能性のある完全陰解法や凸分割法とは異なり、本手法はすべてのタイムステップにおいて線形系をもたらし、計算複雑性を大幅に軽減する。
• 2次精度： 理論的解析および数値実験により、本手法が2次の時間精度を保持することが確認されている。

数値結果
著者らは、一連のベンチマーク計算を通じて手法を検証している：

• 収束解析： 製造解法（Method of Manufactured Solutions, MMS）を用い、本スキームが異なる有限要素構成（$P2-P1-P2-P2$ および $P3-P2-P3-P3$）において、速度、圧力、および相場変数に対して近似的に2次の時間収束を示すことを実証している。
• ロバスト性と正則化： 二相流リッド駆動キャビティ問題において、$\theta$ と $\epsilon$ の相互作用を検討している。$\theta$ がCrank-Nicolsonの極限（例：$\theta=0.51$）に近く、$\epsilon=0$ の場合、弱く減衰する時間的振動によってスキームが不安定になる可能性があることが観察された。正の曲率正則化（$\epsilon > 0$）を導入することで、これらの振動を効果的に減衰させ、大きなタイムステップにおいても有界性とロバスト性を回復させることができる。
• 物理的ベンチマーク：
  • スピノーダル分解： 本手法は、長時間のシミュレーション（$T=100$）において、エネルギー散逸と質量保存を正確に捉え、効果的な相分離と粗大化を示している。
  • 液滴緩和： アルゴリズムは、表面張力駆動の界面運動を正常に解決し、質量を保持しながら正方形の液滴を円形の平衡状態へと遷移させている。
  • リッド駆動キャビティおよびレイリー・テイラー不安定性： 本スキームは、強いせん断、界面変形、および密度駆動の不安定性（低粘性領域におけるケルビン・ヘルムホルツのロールアップを含む）を、目立った偽の振動なしに処理できる。

意義
本論文は、提案された手法のファミリーが、CHNS系に対する実用的かつ効率的な代替案を提供すると主張している。外挿と補助変数によって非線形性をデカップリングし、曲率正則化によって時間離散化を安定化させることで、本手法は2次の精度、無条件の安定性、および線形の計算コストのバランスを実現している。著者らは、本スキームがエネルギー散逸や質量保存を含むモデルの主要な物理構造を保持しており、様々な領域における複雑な界面流のシミュレーションに適していることを強調している。