[論文レビュー] Benchmark Computation of Morphological Complexity in the Functionalized Cahn-Hilliard Gradient Flow

Andrew Christlieb, Keith Promislow+4 more|arXiv (Cornell University)|Jun 8, 2020

Solidification and crystal growth phenomena被引用数 4

ひとこと要約

本稿では、両親媒性系における形態的複雑性を研究するため、関数的化されたカーン・ヒリャード（FCH）勾配流れのベンチマーク計算を導入し、時間方向に2次精度の時間離散化とフーリエ擬スペクトル空間法を用いる。BDF2に基づくスケーム（IMEX、SAV、ETD、完全陰解法PSD）を比較した結果、SAVおよびIMEXが計算効率が最も高く、PSDは計算負荷が大きいものの最高の精度を達成した。

ABSTRACT

Reductions of the self-consistent mean field theory model of amphiphilic molecules in solvent can lead to a singular family of functionalized Cahn-Hilliard energies. We modify these energies, mollifying the singularities to stabilize the computation of the gradient flows and develop a series of benchmark problems that emulate the "morphological complexity" observed in experiments. These benchmarks investigate the delicate balance between the rate of absorption of amphiphilic material onto an interface and a least energy mechanism to disperse the arriving mass. The result is a trichotomy of responses in which two-dimensional interfaces either lengthen by a regularized motion against curvature, undergo pearling bifurcations, or split directly into networks of interfaces. We evaluate a number of schemes that use second order BDF2-type time stepping coupled with Fourier pseudo-spectral spatial discretization. The BDF2-type schemes are either based on a fully implicit time discretization with a PSD nonlinear solver, or upon IMEX, SAV, ETD approaches. All schemes use a fixed local truncation error target with adaptive time-stepping to achieve the error target. Each scheme requires proper "preconditioning" to achieve robust performance that can enhance efficiency by several orders of magnitude.

研究の動機と目的

• 相分離過程における両親媒性系で観察される形態的複雑性を再現する物理的根拠に基づいたベンチマーク問題の開発。
• 界面への両親媒性物質の供給速度とエネルギー最小化再分配機構との間の繊細なバランスの解明。
• 固定局所切捨て誤差を用いた適応的時間刻み法下での、剛性FCH勾配流れに対する複数の2次精度時間積分スケームの性能評価と比較。
• 目標とするグローバル誤差に到達するまでに必要な高速フーリエ変換（FFT）呼び出し回数を指標として、計算効率の公平な比較を確立。
• 極めて非線形な系ではエネルギー安定性のみでは誤った結果が隠蔽される可能性があるため、エネルギー減少と正確性の両立が不可欠であることを示す。

提案手法

• 勾配流れの数値計算を安定化させるために、関数的化されたカーン・ヒリャード（FCH）エネルギーの滑らかで特異性のないバージョンを定式化。
• 高精度およびスペクトル収束を得るために、2次後退差分（BDF2）時間離散化とフーリエ擬スペクトル空間離散化を組み合わせる。
• 4つの数値スケームを実装：(1) 事前条件付き最急降下法を用いた完全陰解法（PSD）、(2) 線形陰解法IMEX、(3) スカラー補助変数（SAV）、(4) 指数時間差分（ETD）を比較の目的で用いる。
• すべてのスケームに同一の固定局所切捨て誤差を目標とする適応的時間刻み法を適用し、一貫した精度を確保。
• 特に剛性非線形系においても、安定性と効率を向上させるために、すべてのスケームに事前条件付き処理を適用。
• SAVアプローチにより、不変エネルギー安定性が保証され、スケーム下でエネルギー減少が保証される。

実験結果

リサーチクエスチョン

• RQ1 適応的時間刻み法下でのFCH勾配流れの形態的複雑性を解明するにあたり、異なる2次精度時間積分スケームの性能はいかが。
• RQ2 剛性FCH系を解くにあたり、異なる数値スケームにおける精度と計算コスト（FFT呼び出し回数で測定）のトレードオフは。
• RQ3 エネルギー安定性のみが正確なシミュレーション結果を保証するのか、それとも複雑なフェーズフィールド系では正確性がエネルギー減少と併せて必要不可欠なのか。
• RQ4 両親媒性物質の供給速度は、曲率駆動運動、パール化、ネットワーク形成といった界面形態選択にどのように影響するか。
• RQ5 事前条件付き処理は、非線形および線形陰解法スケームのFCHシミュレーションにおける性能向上にどの程度寄与するか。

主な発見

• 完全陰解法BDF2-PSDスケームは、固定された局所切捨て誤差下で最小のグローバル離散化誤差を達成しており、テストされたスケームの中で最高の精度を示した。
• IMEXおよびSAVスケームは、同じグローバル誤差レベルに到達するまでにPSDよりもはるかに少ない高速フーリエ変換（FFT）呼び出し回数を要し、計算効率が最も高かった。
• SAVスケームの性能はIMEXに非常に近く、補助変数系の導入によりFFT呼び出し回数が1.4倍程度増加した。
• ETDスケームは比較の目的で含めたが、計算効率の点でIMEXおよびSAVに劣っており、性能が低いことが判明した。
• 事前条件付き処理は、すべてのスケームで安定した性能を実現するために不可欠であり、特に数個のオーダーの計算効率向上を実現した。
• ベンチマーク問題は、実験的に観察された形態的複雑性の三重性（曲率駆動運動、パール化分岐、曲線の分裂によるネットワーク形成）を成功裏に再現した。