[論文レビュー] The pear-shaped fate of an ice melting front

James N. Hewett, Mathieu Sellier|arXiv (Cornell University)|May 7, 2017

Icing and De-icing Technologies参考文献 9被引用数 3

ひとこと要約

本研究では、水の密度逆転点（4.0 °C）付近の加熱された水平シリンダー周囲の氷の融解を、動的メッシュ法とノードシャッフル法を用いて移動する固体-液体界面を追跡することで数値的にシミュレートした。結果として、熱上昇流（T_w > 4.0 °C の場合に上向き、T_w < 4.0 °C の場合に下向き）に起因するピアス型の融解フロントが得られ、自然対流が融解速度および界面形状に顕著な影響を及ぼすことが確認された。形状および全体的な融解進行に関する実験データと良好な一致を示した。

ABSTRACT

A fluid-structure interaction problem with the melting of water around a heated horizontal circular cylinder is analysed with numerical simulations. Dynamic meshing was used for evolving the flow domain in time as the melting front extended radially outward from the cylinder; a node shuffle algorithm was used to retain mesh quality across the significant mesh deformation. We simulated one case above the density inversion point of water and one case below, yielding pear-shaped melting fronts due to thermal plumes either rising or falling from the cylinder, respectively. Results were compared with previous experimental studies and the melting front profiles matched reasonably well and melting rates were in agreement. We confirm that natural convection plays a significant role in the transport of energy as the melt zone increases, and needs to be considered for accurately modelling phase change under these conditions.

研究の動機と目的

• 加熱された水平シリンダー周囲の相変化における、水の密度逆転が融解フロントの形状および動態に与える影響を調査すること。
• 自然対流効果を明示的に追跡する界面追跡によるステファン問題の数値シミュレーションの妥当性を検証すること。
• 上逆転点および下逆転点温度領域における、シミュレートされた融解フロント形状および融解速度を実験データと比較すること。
• 融解領域が拡大するに従い、自然対流が熱伝達をどのように強化するかを評価すること。
• ノードシャッフルを用いた動的メッシュ法が、移動境界問題における大きな変形を扱う際の精度を評価すること。

提案手法

• バーナウリ近似を用いて浮力駆動対流を考慮した2次元非圧縮性層流モデルを採用した。
• 移動境界（融解フロント）は動的メッシュ法により明示的に追跡し、メッシュ品質を保つためにノードシャッフル法を用いた。
• 各時間ステップで、更新されたメッシュ上での流れ場を定常状態として解き、準定常条件を仮定した。
• 水の物性値は密度以外を一定と仮定したが、密度は温度に線形に依存させ、密度逆転効果を再現した。
• 界面の運動は局所的熱フラックスとステファン条件に基づき決定され、人工的な相変化源項は使用しなかった。
• 実験条件に一致させるために、2種類の壁温度（2.3 °C（逆転点以下）および14.1 °C（逆転点以上））でシミュレーションを実施した。

実験結果

リサーチクエスチョン

• RQ1 水の密度逆転点（約4.0 °C）が、加熱された水平シリンダー周囲の融解フロントの形状および進化にどのように影響を与えるか。
• RQ2 融解領域が拡大するに従い、自然対流が熱伝達をどれほど支配的にするか。また、逆転点の上下でその影響はどのように変化するか。
• RQ3 ノードシャッフルを用いた動的メッシュ法が、実験的に観察された複雑なピアス型界面形状を正確に捉えられるか。
• RQ4 なぜ、伝導支配から対流支配への遷移域（τ ≈ 0.3）において、シミュレーションと実験との間に融解速度のずれが生じるのか。
• RQ5 2次元定常状態の仮定が、3次元実験で観察された一時的不安定性および強化された熱伝達の予測にどのように影響を与えるか。

主な発見

• T_w = 14.1 °C の場合に上向き、T_w = 2.3 °C の場合に下向きの熱上昇流に起因するピアス型の融解フロントが両ケースで発生し、実験観察と一致した。
• レイノルズ数は逆転点の場合6700、上向きの場合7400であり、異なる流れ構造にもかかわらず、対流の強度はほぼ同等であった。
• 融解体積が増大するに従い、自然対流が熱伝達の支配的メカニズムに発展し、界面の進化に顕著な影響を与えた。
• シミュレーション全期間にわたり、実験的プロファイルと妥当な一致を示し、モデルの妥当性を確認した。
• 遷移域（τ ≈ 0.3）における融解速度のずれは、2次元定常モデルに存在しない非定常的3次元効果および一時的不安定性に起因するとされた。
• 上向きピアス型形状では、密度逆転点の包含に起因し、シリンダー上部に逆向きの渦が発生するなど、より複雑な流れ特徴を示した。