Elucidation of ice slurry flow with heat transfer by microscopic approach including solid-liquid phase change

固液相変化モデル（融解・凝固モデル）の構築とそれを伴う熱流動解析手法の開発を行った．本研究では，応力テンソルの不連続に基づく埋め込み境界－格子ボルツマン法に，温度の計算と融解・凝固による境界の移動計算を導入し，融解・凝固を伴う熱流動問題を解析するための手法の開発を行った．まず，実装した融解・凝固モデルの妥当性検証として，一次元半無限領域，二次元正方形領域，二次元円管領域，および三次元半無限領域における問題の計算を行った．その結果，固相が自由運動しない場合では既存の結果と良く一致することがわかった．一方，固相が自由運動する場合では，固相にはたらく流体力の差に起因する結果の相違が見られた．ただし，比較した参照データは，内部質量の影響を計算していないため，流体力を正しく評価できていないことが考えられることに注意が必要である．次に，構築した融解・凝固を伴う熱流動計算法を氷スラリー（微細な氷と液体の固液二相混合物）の熱流動問題に適用した．二次元正方領域内に氷粒子（一定の低温条件）を一個配置し，壁面となる境界には一定の高温条件，流れ方向の入口・出口間には圧力差を伴う周期境界条件を用いることで，発熱する二次元流路内を氷スラリーが流動する様子を再現し，計算を実行した．その結果，流路幅および空塔速度で無次元化した時刻が450の時，氷粒子の質量は初期に比べて約45%減少した．また，融解による固相の変形がない場合と比較して，氷粒子の重心は流路中心に向かうことがわかった．さらに，壁面上の空間平均ヌッセルト数は若干低下する結果となった．これは，氷粒子の表面積が減少したことや，熱源である壁面から遠ざかったことが原因と考えられる．

固液相変化モデル（融解・凝固モデル）の構築とそれを伴う熱流動解析手法の開発を行った．本研究では，応力テンソルの不連続に基づく埋め込み境界－格子ボルツマン法に，温度の計算と融解・凝固による境界の移動計算を導入し，融解・凝固を伴う熱流動問題を解析するための手法の開発を行った．まず，実装した融解・凝固モデルの妥当性検証として，一次元半無限領域，二次元正方形領域，二次元円管領域，および三次元半無限領域における問題の計算を行った．その結果，固相が自由運動しない場合では既存の結果と良く一致することがわかった．一方，固相が自由運動する場合では，固相にはたらく流体力の差に起因する結果の相違が見られた．ただし，比較した参照データは，内部質量の影響を計算していないため，流体力を正しく評価できていないことが考えられることに注意が必要である．次に，構築した融解・凝固を伴う熱流動計算法を氷スラリー（微細な氷と液体の固液二相混合物）の熱流動問題に適用した．二次元正方領域内に氷粒子（一定の低温条件）を一個配置し，壁面となる境界には一定の高温条件，流れ方向の入口・出口間には圧力差を伴う周期境界条件を用いることで，発熱する二次元流路内を氷スラリーが流動する様子を再現し，計算を実行した．その結果，流路幅および空塔速度で無次元化した時刻が450の時，氷粒子の質量は初期に比べて約45%減少した．また，融解による固相の変形がない場合と比較して，氷粒子の重心は流路中心に向かうことがわかった．さらに，壁面上の空間平均ヌッセルト数は若干低下する結果となった．これは，氷粒子の表面積が減少したことや，熱源である壁面から遠ざかったことが原因と考えられる．