時間変化する境界をもつ粘性流体の運動に対するSPHを用いた基礎研究

流体表面と物体の衝突過程に対する最近の研究として、L.Boquetらの石の水切りの実験[Nature, vol.427,29(2004)]に着目した。彼らはアルミニウムの円盤に回転を与えて水面に対して衝突させ、円盤の入射角度θと水面に対する傾き角度φを変えた時の反発の起こる条件を調べた。その結果、傾き角度φ=20°の時に反発が起こる為に必要な速度v_<min>が最小になる事を見いだした。これを再現する事、また傾き角度φについて最小値が生じるその物理的理由を明らかにする事を目的として、3次元の流体シミュレーションとモデル解析の双方から研究を行った。シミュレーションを行う場合の困難の一つは、大変形を起こす流体表面の境界条件の扱いである。ナビエ・ストークス方程式をラグランジュ描像で数値的に解くSPH(Smoothed Particle Hydrodynamics)法を用いた。SPH法は、流体を流体粒子の集まりとして記述する為、変形する表面を扱う際に格子生成などの煩雑な手続きをさける事が出来る点がかかる問題に対して有利な特徴である。表面張力を無視し、粘性の殆ど効かない領域で行ったSPHシミュレーションが、L.Boquetらによる実験を定性的に再現することを確認した。この結果から「石の水切り」のようなレイノルズ数が10^5のオーダーに達する衝突過程においては、物体が流体に運動量を与える事によって受ける反作用力が最も重要であると予想できる。これを考慮し流体表面の変形を無視した単純な微分方程式モデルを提案し、円盤の衝突に対して適用した。微分方程式を近似的に扱い、v_<min>とθ,φの関係について実験と定性的に一致する解析的な表式を得ることができた。

流体表面と物体の衝突過程に対する最近の研究として、L.Boquetらの石の水切りの実験[Nature, vol.427,29(2004)]に着目した。彼らはアルミニウムの円盤に回転を与えて水面に対して衝突させ、円盤の入射角度θと水面に対する傾き角度φを変えた時の反発の起こる条件を調べた。その結果、傾き角度φ=20°の時に反発が起こる為に必要な速度v_<min>が最小になる事を見いだした。これを再現する事、また傾き角度φについて最小値が生じるその物理的理由を明らかにする事を目的として、3次元の流体シミュレーションとモデル解析の双方から研究を行った。シミュレーションを行う場合の困難の一つは、大変形を起こす流体表面の境界条件の扱いである。ナビエ・ストークス方程式をラグランジュ描像で数値的に解くSPH(Smoothed Particle Hydrodynamics)法を用いた。SPH法は、流体を流体粒子の集まりとして記述する為、変形する表面を扱う際に格子生成などの煩雑な手続きをさける事が出来る点がかかる問題に対して有利な特徴である。表面張力を無視し、粘性の殆ど効かない領域で行ったSPHシミュレーションが、L.Boquetらによる実験を定性的に再現することを確認した。この結果から「石の水切り」のようなレイノルズ数が10^5のオーダーに達する衝突過程においては、物体が流体に運動量を与える事によって受ける反作用力が最も重要であると予想できる。これを考慮し流体表面の変形を無視した単純な微分方程式モデルを提案し、円盤の衝突に対して適用した。微分方程式を近似的に扱い、v_<min>とθ,φの関係について実験と定性的に一致する解析的な表式を得ることができた。