[論文レビュー] Three-dimensional phase-field simulations of water freezing and thawing at pore-scale

Pavel Strachota|arXiv (Cornell University)|Jan 3, 2024

Climate change and permafrostEarth and Planetary Sciences被引用数 3

ひとこと要約

本研究では、ガラスビーズで構成されるポーラス媒体内の氷の凍結・融解を3次元スケールでシミュレートするためのフェーズフィールドシミュレーションフレームワークを提示する。新しい粒子動力学モデルを用いてポーラス構造を生成し、表面張力あり・なしの両状態での相転移をシミュレートするハイブリッド並列ソルバを採用した。主な貢献は、表面張力に起因する事前融解効果が、定性的な比較と計算効率の分析を通じて、ポーラススケールで数値的に成功裏に示されたことである。

ABSTRACT

This work deals with numerical simulation of water freezing and thawing in a complex three-dimensional geometry of a porous medium. The porous structure is represented by a virtual container filled with glass beads. Phase transition modeling is approached at both macro-scale and micro-scale, combining heat transfer in a heterogeneous medium and a phase-field approximation of the Gibbs-Thomson relation by means of the Allen-Cahn equation. The formulation of the model contains novel components tailored for the given purpose. At the macro-scale, surface tension effects are negligible and phase transition focusing based on temperature can replace the Allen-Cahn equation. In contrast to that, simulations of equilibrium states at the micro-scale allow to eliminate the heat equation by assuming constant supercooling. For numerical solution, an efficient hybrid parallel algorithm based on the finite volume method and the Runge-Kutta-Merson solver with adaptive time stepping are employed. The results of different model variants at different scales are discussed. In a parametric study, the full phase-field model is demonstrated to deliver consistent results across a wide range of surface tension values, exhibiting curvature-induced premelting if surface tension is artificially exaggerated. As surface tension tends to the realistic values, the results of the phase-field approach those of the simplifed temperature-driven phase transition model. In addition, micro-scale simulations of water freezing at different supercooling values aim to predict the unfrozen water content and compare the results with data from literature. Numerical stability, accuracy, and computational costs are also discussed.

研究の動機と目的

• ポーラス媒体内の水の凍結・融解を3次元スケールで計算的に効率的に行う数値フレームワークの開発。
• 跳ね返り係数を調整可能なシンプルで正確な粒子動力学アルゴリズムを用いて、球状のガラスビーズによる現実的なポーラス床の形成をモデル化。
• 表面張力を含む・含まないフェーズフィールドモデルを比較し、表面張力が事前融解行動に与える影響を調査。
• 有限体積法と適応時間刻みルンゲ・クッタ法を組み合わせたハイブリッド並列ソルバを実装し、効率的なシミュレーションを実現。
• 今後の実験データとの検証の基盤を提供するため、シミュレーション設定を実験的配置と整合させる。

提案手法

• 球体間の指数関数的減衰力に基づく粒子動力学モデルを開発し、重力と相互衝突を伴う条件下で200個のガラスビーズが容器内に沈降する過程をシミュレート。
• このモデルにより、0（完全に非弾性）から1（完全に弾性）までの跳ね返り係数を有する正確な衝突シミュレーションが可能となり、不自然な振動やジャイタリーな状態を伴わない安定した静止状態が得られる。
• 熱伝導と相転移をモデル化するためにフェーズフィールド法を適用し、表面張力を含む（事前融解を捉える）と含まない（無視）の2種類のバージョンを用意。
• 不均質なポーラス幾何形状に適した新しい要素を組み込んだフェーズフィールド定式化を採用し、温度依存の相転移進行を組み込んだ。
• 偏微分方程式系の効率的かつ安定した解法を実現するため、有限体積法と適応時間刻みルンゲ・クッタ法（Runge-Kutta-Merson）を組み合わせたハイブリッド並列アルゴリズムを実装。
• 数値ソルバはMATLABおよびC++で実装され、MITライセンスのもとで公開されたGitHubリポジトリを通じて完全に再現可能である。

Figure 1: Bounce test with the coefficient of restititution $e=1$ and $e=\sqrt{2}/2$ . The ode15s solver was used with options given by odeset(’RelTol’,1e-10,’AbsTol’,1e-8) . Increasing the accuracy further did not result in any visible difference.

実験結果

リサーチクエスチョン

• RQ1 3次元ポーラス媒体内におけるポーラススケールでの凍結・融解過程において、表面張力は液体・固体界面の形成と進化にどのように影響を与えるか？
• RQ2 指数関数的減衰力に基づくシンプルな粒子動力学モデルは、現実的なパッケージングを持つポーラス床への球状ビーズの自己組織化を正確に再現できるか？
• RQ3 フェーズフィールドモデルに表面張力を組み込むことで、高曲率のポーラス領域における未凍結水含量（事前融解）の予測にどのような影響が生じるか？
• RQ4 3次元ポーラススケールでの凍結・融解シミュレーションにおいて、グリッド解像度とモデルの複雑さが計算コストにどのように影響するか？
• RQ5 提示された数値フレームワークは、実験データとの今後の検証を支えるのに十分に安定的で正確かつ効率的と言えるか？

主な発見

• 粒子動力学モデルにより、重力による沈降と衝突処理を経て、不自然な振動やジャイタリーな状態を伴わず、安定的かつ現実的なガラスビーズのポーラス床が生成された。
• 表面張力を含むフェーズフィールドモデルは、物理的期待と整合するように、高曲率の氷表面に未凍結水の薄い層を正しく再現した。
• 表面張力を無視したモデルでは、より鋭く急激な相転移が得られ、界面の物理的挙動を的確に捉えるために表面張力の重要性が浮き彫りになった。
• 3次元ポーラススケールシミュレーションにおいて計算コストが実用的であることが判明し、384コアで100³グリッドでは約1時間、400³グリッドでは約90時間の合計シミュレーション時間がかかった。
• ハイブリッド並列ソルバは良好なスケーラビリティと数値的安定性を示し、適応時間刻みにより急激な相転移時でも精度が保たれた。
• 結果は既知の物理的現象と定性的に整合しており、今後の実験的未凍結水含量および界面形状の検証へのモデルの可能性を裏付けた。

(a) Initial configuration