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