本研究では、衝撃波と呼ばれる圧縮性流体における不連続現象の数値解析において、異常解と呼ばれる非物理的解の発生について理解を深め、異常解の発生を克服する計算手法を構築し、工学の諸問題へと応用することを目的として研究を実施した。特に、これまでに知見が多くない空間を伝播する「移動衝撃波」を対象に解析、数値実験を実施し、得られた知見を基に数値解法を提案した。 まず、数値解析で生じる移動衝撃波の内部構造に注目し、衝撃波マッハ数や計算スキーム(空間再構築、セル境界流束の計算方法)を変化させ計算結果を整理した。弱い衝撃波では内部構造の拡散が顕著であることが明らかとなり、いくつかの手法では計算が不安定となった。こうした内部構造がどういった物理に支配されているのかを明らかにするため、衝撃波内部におけるエントロピー増加に着目し解析した。解析の結果弱い衝撃波ではエントロピー生成が過剰となっており、結果として散逸的な解となることが示唆された。この結果は、現状の計算手法は衝撃波の幅広い強度に対応できていないことを示しており、今後改善していくことが必要である。 数値解析手法については、セルの物理量を急峻な関数を使用して補完する手法を使用することによって微弱衝撃波をシャープに捉える手法を提案した。結果として、衝撃波のみではなくその他の流れに対しても解像度が向上した。これらについて学会発表(国内、国外)を行い、流体力学の国際雑誌に論文が掲載された。 実験による数値解析手法の検証や、実験研究や工学流れへの応用を見据えた研究も推進した。超音速風洞を整備し、実験結果と数値解析結果を光学可視化によって比較する環境を整えた。また粒子を含む流れに対しても数値解析手法の展開を実施した。
|