| 研究課題/領域番号 |
21K03873
|
|---|
| 研究機関 | 信州大学 |
|---|
研究代表者 |
吉野 正人 信州大学, 学術研究院工学系, 教授 (00324228)
|
|---|
| 研究分担者 |
鈴木 康祐 信州大学, 学術研究院工学系, 准教授 (10735179)
浅岡 龍徳 信州大学, 学術研究院工学系, 准教授 (30508247)
|
|---|
| 研究期間 (年度) |
2021-04-01 – 2024-03-31
|
| キーワード | 格子ボルツマン法 / 氷スラリー / 固液二相流 / 融解・凝固 / 相変化 |
|---|
| 研究実績の概要 |
固液相変化モデル（融解・凝固モデル）の構築とそれを伴う熱流動解析手法の開発を行った．本研究では，応力テンソルの不連続に基づく埋め込み境界－格子ボルツマン法に，温度の計算と融解・凝固による境界の移動計算を導入し，融解・凝固を伴う熱流動問題を解析するための手法の開発を行った．
まず，実装した融解・凝固モデルの妥当性検証として，一次元半無限領域，二次元正方形領域，二次元円管領域，および三次元半無限領域における問題の計算を行った．その結果，固相が自由運動しない場合では既存の結果と良く一致することがわかった．一方，固相が自由運動する場合では，固相にはたらく流体力の差に起因する結果の相違が見られた．ただし，比較した参照データは，内部質量の影響を計算していないため，流体力を正しく評価できていないことが考えられることに注意が必要である．
次に，構築した融解・凝固を伴う熱流動計算法を氷スラリー（微細な氷と液体の固液二相混合物）の熱流動問題に適用した．二次元正方領域内に氷粒子（一定の低温条件）を一個配置し，壁面となる境界には一定の高温条件，流れ方向の入口・出口間には圧力差を伴う周期境界条件を用いることで，発熱する二次元流路内を氷スラリーが流動する様子を再現し，計算を実行した．その結果，流路幅および空塔速度で無次元化した時刻が450の時，氷粒子の質量は初期に比べて約45%減少した．また，融解による固相の変形がない場合と比較して，氷粒子の重心は流路中心に向かうことがわかった．さらに，壁面上の空間平均ヌッセルト数は若干低下する結果となった．これは，氷粒子の表面積が減少したことや，熱源である壁面から遠ざかったことが原因と考えられる．
|
| 現在までの達成度 (区分) |
現在までの達成度 (区分)
2: おおむね順調に進展している
理由
構築した融解・凝固モデルは，固相が自由運動する場合においてさらなる検討が必要であるが，概ね確立したと言える．
|
| 今後の研究の推進方策 |
構築した融解・凝固を伴う熱流動解析手法を三次元に拡張し，大規模並列計算用の固液二相系埋め込み境界－格子ボルツマン計算コードを作成する．また，氷スラリー流れの三次元熱流動シミュレーションを実施し，氷粒子が流れる空間的位置と熱伝達特性の関係を明らかにする．
|
| 次年度使用額が生じた理由 |
データ保存用のHDDの必要個数が減少したため，当初の金額と使用額に若干差が生じた．次年度は，研究成果報告のための学会参加費の一部として使用する予定である．
|
