[論文レビュー] Fixed-grid sharp-interface numerical solutions to the three-phase spherical Stefan problem

Yavkreet Swami, Jacob Barajas+1 more|arXiv (Cornell University)|Feb 9, 2026

Phase Change Materials Research被引用数 0

ひとこと要約

この論文は、有限サイズの粒子に対して3相球対流シュレーダー問題（ステファン問題）を固定格子のシャープインターフェース法へ拡張し、密度と運動エネルギーのジャンプを組み込み、二相結果と整合させ、ノイズ粒子スケールでは運動エネルギー効果が有意だがサイズが大きくなると低減することを示す。

ABSTRACT

Many metal manufacturing processes involve phase change phenomena, which include melting, boiling, and vaporization. These phenomena often occur concurrently. A prototypical 1D model for understanding the phase change phenomena is the Stefan problem. There is a large body of literature discussing the analytical solution to the two-phase Stefan problem that describes only the melting or boiling of phase change materials (PCMs) with one moving interface. Density-change effects that induce additional fluid flow during phase change are generally neglected in the literature to simplify the math of the Stefan problem. In our recent work [1], we provide analytical and numerical solutions to the three-phase Stefan problem with simultaneous occurrences of melting, solidification, boiling, and condensation in Cartesian coordinates. Our current work builds on our previous work to solve a more challenging problem: the three-phase Stefan problem in spherical coordinates for finite-sized particles. There are three moving interfaces in this system: the melt front, the boiling front, and the outer boundary which is in contact with the atmosphere. Although an analytical solution could not be found for this problem, we solved the governing equations using a fixed-grid sharp-interface method with second-order spatio-temporal accuracy. Using a small-time analytical solution, we predict a reasonably accurate estimate of temperature (in the three phases) and interface positions and velocities at the start of the simulation. Our numerical method is validated by reproducing the two-phase nanoparticle melting results of Font et al. [2]. Lastly, we solve the three-phase Stefan problems numerically to demonstrate the importance of kinetic energy terms during phase change of smaller (nano) particles. In contrast, these effects diminish for large particles (microns and larger).

研究の動機と目的

• 平面幾何学から有限サイズ粒子の球対称幾何へ、3相シュテファン問題のモデリングを拡張する。
• すべての界面において密度ジャンプと運動エネルギー跳躍条件を保持する。
• 球座標系での固定格子シャープインターフェースを用いた埋め込み境界法を開発・適用する。
• 界面の移動を正確に開始するため、初期時刻近傍の小時間解析解を用いて初期化する。
• 既知の二相結果と比較して手法を検証し、粒子サイズにわたる運動エネルギー項の役割を分析する。

提案手法

• 境界条件として密度・エネルギー・運動量を含む固体・液体・蒸気の支配方程式を、界面のジャンプ条件とディラックデルタ項とともに定式化する。
• 支配方程式を無次元化し、初期化のための小時間解析解を導出する。
• 球座標系で移動界面を埋め込み境界追跡する、二次精度の固定格子シャープインターフェース法を実装する。
• 低・高密度比極限下の小時間解析を用いて初期の界面位置・速度および温度場を得る。
• 有限サイズの球上で3相球対 Stefan問題を解き、界面の進化と密度・運動エネルギー項の影響を調べる。
• Font らの二相ナノ粒子溶融結果を再現することで数値的に検証する。

(a) Two-phase spherical Stefan problem

実験結果

リサーチクエスチョン

• RQ1 有限サイズ粒子における同時融解・沸騰境界がある3相球対 Stefan問題はどのように発展するか。
• RQ2 異なる粒子スケールで、密度ジャンプと運動エネルギー跳躍項が界面位置・速度・融解/沸騰時間に与える影響はどの程度か。
• RQ3 多相シュテファン問題に対して球座標系の固定格子シャープインターフェース法はどの程度正確か。
• RQ4 ノノ粒子尺度とミクロン粒子尺度で運動エネルギー項が融解時間にどの程度影響するか。

主な発見

• 運動エネルギー項は小さな（ナノ）粒子の融解時間に著しく影響する。
• より大きな粒子（ミクロンスケール以上）では運動エネルギー項の影響は低下する。
• 固定格子シャープインターフェース法は空間・時間で二次の精度を達成する。
• Font らによる二相ナノ粒子の融解結果を再現することで数値的検証に成功する。
• 三相シミュレーションは、同時融解・沸騰の過程で、相間の密度・運動エネルギー項の役割を定量化し、小さな粒子での重要性を強調する。

(b) Three-phase spherical Stefan problem