| 研究実績の概要 |
多相連成災害の素過程を解明する計算力学手法の構築するため、本研究では、水流と土粒子を想定した複数の固体粒子の連成現象を数値解析により予測する手法を構築し、基本的な実験結果を利用して解法の検証を行った。今年度および研究期間全体の具体的な検討内容を以下に示す。
底面から流入する上昇水流による飽和粒子層の内部流動化と破壊に至る過程を詳細に把握するため、粒子層厚を120 mm と増加させた実験を行い、粒子の挙動を高速ビデオカメラで撮影するとともに、間隙水圧計により粒子間の圧力の計測を行った。実験では、平均粒径約 4 mm の礫粒子と、粒径 4 mm のガラスビーズを利用して、粒子形状の違いによる現象の相違の確認も行った。一方、これらの実験を対象とする粒子スケールの数値解析を行い、実験結果と比較した。数値計算では、粒子周辺の流れも解析する分解能の流体計算を行った。四面体要素から構成される礫粒子モデルの総数は14,341で、流体計算セル数は約5千3百万とし、京都大学のスーパーコンピュータのコア数1,088を利用する並列計算により粒子スケールの数値計算を行った。粒子間の接触は個別要素法により評価している。数値計算の結果、水流が底面から流入すると、時間の経過とともにキャビティとよばれる空洞が徐々に大きくなり、それが粒子層上面に到達して水流が貫通し、粒子層の破壊に至る過程が観察されたが、これと同様の現象が再現されることが確認できた。粒子形状の相違による粒子運動のパターンについても同様の傾向が数値計算により得られることが示された。また、間隙水圧の時間的・空間的な変化についても、数値計算により概ね再現できることが確認された。
以上のように、固体粒子に作用する流体力を、粒子周辺の流体計算結果から評価するという、素過程を再現した数値解析手法を本研究により構築し、その有効性を確認することができた。
|
