[論文レビュー] High efficiently numerical simulation of the TDGL equation with reticular free energy in hydrogel

Jun Han, Hui Zhang+1 more|arXiv (Cornell University)|May 29, 2017

Material Dynamics and Properties参考文献 13被引用数 4

ひとこと要約

本稿では、網状自由エネルギーを伴う時間発展ギンツブルグ＝ランダウ（TDGL）方程式を用いて、マクロマoleculeマイクロスフィア複合（MMC）ハイドロゲルの相分離をシミュレートするための、二つの$L^2$安定で2次精度の半implicit有限差分スキームを提案する。適応的時間刻み戦略が導入されており、特に自由エネルギーの低下が著しい初期段階では小さなステップを、変化が小さい後期段階では大きなステップを用いることで、計算時間（CPU時間）を最大80％まで削減しつつ、高い精度を維持する。

ABSTRACT

In this paper, we focus on the numerical simulation of phase separation about macromolecule microsphere composite (MMC) hydrogel. The model equation is based on Time-Dependent Ginzburg-Landau (TDGL) equation with reticular free energy. We have put forward two $L^2$ stable schemes to simulate simplified TDGL equation. In numerical experiments, we observe that simulating the whole process of phase separation requires a considerably long time. We also notice that the total free energy changes significantly in initial time and varies slightly in the following time. Based on these properties, we introduce an adaptive strategy based on one of stable scheme mentioned. It is found that the introduction of the time adaptivity cannot only resolve the dynamical changes of the solution accurately but also can significantly save CPU time for the long time simulation.

研究の動機と目的

• MMCハイドロゲルの相分離を支配するTDGL方程式と網状自由エネルギーを用いた、堅牢な数値スキームの開発。
• $L^2$ノルム安定スキームを用いて、数値的安定性とエネルギー減少性を保証する。
• 長時間シミュレーションに伴う高い計算コストに対処するため、適応的時間刻み戦略を導入する。
• 数値実験を通じて、適応スキームの効率性と精度を検証する。

提案手法

• MMCハイドロゲルの網状自由エネルギーに基づく変分原理を用いて、TDGL方程式を導出し簡略化する。
• 二つの2次精度の半implicit有限差分スキームを提案：一つは高次元線形方程式系の共役勾配法による解法を含み、もう一つは非線形方程式系に対するニュートン反復法を用いる。
• 自由エネルギーの変化率に基づいて時間刻みを調整する適応的時間刻み戦略を実装する—エネルギーの低下が著しいときには小さなステップ、変化が小さいときには大きなステップを用いる。
• スキームのエネルギー減少性を根拠として、後期シミュレーション段階で大きな時間ステップを安全に使用できることを裏付ける。
• 定常時間刻みと適応的時間刻みの比較を、CPU時間、$L^2$誤差、自由エネルギーの時間発展の観点から数値実験で行う。
• 細かい定常時間刻み（$\Delta t = 0.001$）との比較を通じて、適応戦略の妥当性を自由エネルギー曲線と$L^2$誤差の観点から検証する。

実験結果

リサーチクエスチョン

• RQ1 網状自由エネルギーを伴う簡略化されたTDGL方程式に対して、$L^2$安定で2次精度の半implicitスキームを構築可能か？
• RQ2 MMCハイドロゲルの相分離の長時間シミュレーションにおける計算コストを、精度を損なわず低減する方法は何か？
• RQ3 自由エネルギーのダイナミクスに基づく適応的時間刻み戦略は、「粗いから滑らかへ」変化を示す系のシミュレーション効率を向上させるか？
• RQ4 適応的時間刻みと定常時間刻みの比較において、CPU時間と精度の定量的向上はどの程度か？

主な発見

• $t = 500$における適応的時間刻み戦略は、$\Delta t = 0.001$と比較してCPU時間を約80％削減しており、2809.4秒（適応）vs. 14341.4秒（定常）である。
• $t = 200$における適応的戦略の$L^2$誤差は0.000325であり、定常ステップ（$\Delta t = 0.05$：0.0376、$\Delta t = 0.01$：0.0083）と比較して顕著に低い。
• 適応的戦略下での全自由エネルギーの時間発展は、$\Delta t = 0.001$の結果と非常に良く一致しており、後期段階で大きな時間ステップを用いても高い精度が保たれていることが確認された。
• スキームはエネルギー減少性と$L^2$安定性の両方を維持しており、後期シミュレーション段階で大きな時間ステップを安全に使用できる。
• 初期段階では相分離が急速に進行し、自由エネルギーの低下が著しいため、初期段階での小さな時間ステップが正当化される。
• 適応法は、後期段階で$\Delta t = 0.001$の100倍の時間ステップを用いることで、誤差の増加が最小限に抑えられ、高い効率性を達成している。