本研究は乱流場を形成する粒子層に対し、平均化操作を施し巨視的方程式を導出することを目的としている。まずダクト内にアクリル製の円柱群を配置することにより充填層を単純化し、円柱群内での圧力損失、流速および乱流エネルギーの測定、さらにスモークワイヤー法による流れの可視化を行い、充填層内における乱流挙動を検討した。次に円柱群内の乱流場に対して、数値解析モデルの構築を行った。乱流モデルにはk-εモデルまたはk-lmモデルを用いて、円柱群を配置した場合および円柱群を空隙率によって表現した場合について、円柱群内における乱流挙動を検討した。以下に得られた成果をまとめる。 1.圧力損失はレイノルズ数の増加に伴い増大し、レイノルズ数が一定の場合は空隙率が小さいほど増大した。 2.乱流エネルギーは円柱群への流入と同時に急激に増大するが、下流にいくに伴い一定値に漸近する傾向を示した。またレイノルズ数の増加に伴い乱流エネルギーは増大した。 3.流れの可視化より、前列の円柱間からの噴流が円柱に衝突することにより乱れが生じ、下流に流れるに従い乱れが発達する様子が観察された。しかしながら乱流渦は円柱間の空隙により制限されるため、下流域においては乱れの発達は見られなかった。 4.円柱群を配置した数値解析では、円柱群に流入すると同時に乱流エネルギーが増大し、下流にいくに伴い一定値に漸近するという実験結果と同様の傾向を示した。 5.円柱群を空隙率によって表現した数値解析では、円柱群による乱れの生成を新たに見積もった。k-lmモデルでは実験結果と同様の傾向が得られたが、混合距離の評価によって大きく解析結果が変わってしまったため、混合距離の決定法が課題となる。また、高次モデルへの展開としてk-εモデルによる数値解析を試みたが、流入条件の設定などの課題があり、今後の研究で解明していく予定である。
|