H. Kuninaka and H. Hayakawa Phys. Rev. Lett. 93, 543 (4) M. Y. Louge & E. Adams Phys. Rev. E 65, 33 () Y.Tanaka et. al. Europhys. Lett. 63, 46 (3) K. Okumura et. al. Europhys. Lett. 6, 37 (3) R +X h R +X h (a) (b) J. W. Glasheen and T. A. McMahon, Nature, 38 (996) C. Clanet, F. Hersen and L. Bocquet. Nature, 47 (4)
Weight[Cwt.] distance [Yards] Number of ricochet 84 85 6 65 9 3 θ c = 8/ σ σ = c
G. Birkhoff et. al. (944) W. Johnson & S. R. Reid (975) I. M. Hutchings (976) 円柱衝突の理論 Birkhoffモデル 3 5 6 y O y dv m = dt w! pn ds x 4 E. G. Richardson, Proc. Phys. Soc. London, Sect. A 6 (948) surface π sin α pn = ρv 4 + π sin α (Rayleighの式) θc = 7.8/ σ 年月9日木曜日 Flow direction
v F = v gl v C. Clanet, F. Hersen and L. Bocquet. Nature, 47 (4)
6 6 5 5 入射角度 θ [deg.] Umin[m/sec] 石の水切りの実験 Skipping Stone 4 3 4 6 円盤の傾き角度 φ [deg.] v = 3.5[m/sec.] 4 3 Skipping stone 4 円盤の傾き角度 φ [deg.] Magic Angle φm 6 θ φ v C. Clanet, F. Hersen and L. Bocquet. Nature, 47 (4) 年月9日木曜日
SPH( Smoothed Particle Hydrodynamics) MPS( Moving Particle Semi-implicit ) CIP Cubic Interpolated Propagation Level-set
ρ (r) = mw (r r j ) j W(r j): Gaussian karnel. v j A (r) = j A (r) = j m A j W (r r j ) ρ j m A j W (r r j ) ρ j ρ j ρ (r j ) ( ) c ρ i ρ ρ i ρ ρ p i = = ρ i < ρ
SPH法 () Navier-Stokes方程式 u + (u )u = p + ν u + t ρ! " ζ + ν/3 u ρ! m!a"(r) = Aj W (r r j )!ρ" j j に基づく補間の手続き SPHの運動方程式! mb d!u"i = dt ρi ρj j " uij r ij pi + pj νξ rij # i Wij 斥力相互作用する粒子系の運動方程式 年月9日木曜日
Length of smoothing kernel vi Fluid particles fixed on solid wall (Wall particles)
r = v z g x Lz Lx
円柱と水面の衝突 ー衝突の様子と速度プロファイル 上図a-d; E. G. Richardson, Proc. Phys. Soc. London, Sect. A 6 (948)より引用 年月9日木曜日
斜め衝突の間に円柱が受ける力 粒子の初期配置の影響 粒子の配置を正方格子とランダム配置に選んだ二つの結果を比較 force v r.5 force v r Square Lattice Random.6 Square Lattice Random.5.4 fz.3 fz.5.. T=r/cで -.5 - 時間平均処理 fx -. fx -. -.3 -.4 3 4 t 5 6 7 8 3 4 5 t 6 7 8 r/v r/v 境界(底面)の影響 深さの異なる水槽 用いたシミュレーションの結果.4 force.3 v r fz Lz=7r Lz=4r fz fx. fz. fx -. fx -. -.3 3 4 5 t r/v 年月9日木曜日 6 7 8
θ c tank size lx=3r, lz=4r lx=4r, lz=8r θ c = 8/ σ Critical Angle [deg.] -.5. Specific gravity
v
5cm θ =.5, φ =, ω = 6[rounds/sec.]
5cm θ =.5, φ =, ω = [rounds/sec.]
v Minimum Velocity [m/sec.] 7 6 5 4 3 Experiment SPH Incident Angle [deg.] 6 5 4 3 Experiment SPH 3 4 5 6 Tilt Angle [deg.] 3 4 5 6
v 3 5 5 θ =8. φ = 6. 4 35 3 5 5 3 5 5 3 4
z φ = const. p ρ(v n). z' f R x = const. f = pnds = disk s face g C D ρ(v n) n/ds disk s face = C Dρ(v n) Sn d v Water x'
衝突のODEモデル 1 無次元化した運動方程式 x z = sin φ F = κs(z )z cos φ F! z (x,z) v F = gr f x R = const. CD Rρ κ= πdρ! g z' xz座標系 v d 円盤下側の角の位置 フロード数 #! "! " "!!!! z z π z + S(z! ) = arcsin + + sin φ sin φ sin φ 流体にひたっている面積 年月9日木曜日 x' Water
ODEモデルと SPHシミュレーションの比較 定数κの決定 κ=.94, (Fitting parameter) 円板の受ける力 円板の軌道.6.35 fz.3.5.4 fx.3.5. z/r....5 -. fx -.5 年月9日木曜日.5 -. SPH ODE -. -.5.5 time r/v.5 ODE SPH.5 fz 3 3.5 -.3 4 -.4 3 4 x/r 5 6 7 8
Minimum Velocity vmin [m/sec.] 7 6 5 4 3 Criterion A Experiment SPH Theory 3 4 5 6 Tilt Angle [deg.] Incident Angle [deg.] 6 5 4 3 Experiment SPH Theory 3 4 5 6 Tilt Angle [deg.]
衝突のODEモデル 解析解 反発 の定義 通常は 水面の高さを基準 にとる 位置条件 ここでは 重力方向の速度を基準にとる 速度条件 円盤の重力方向の速度が反転したら 反発が 起こったものと 見なす x (a) z (b) (c) (d) 年月9日木曜日 = sin φ F = κs(z )z cos φ F! 変曲点の存在条件をしらべる 球と水面の衝突の実験で得られた 衝突後の球の軌道 E. G. Richardson, Proc. Phys. Soc. London, Sect. A 6 (948)
v min = θ max = arccos gr cos(θ + φ) F { x sin φ + ( x sin φ + σd cos φ C D R sin φ } ) σd cos φ C D R sin φ φ Minimum Velocity vmin [m/sec.] 7 6 5 4 3 Experiment SPH Theory 3 4 5 6 Tilt Angle [deg.] Incident Angle [deg.] 6 5 4 3 Experiment SPH Theory 3 4 5 6 Tilt Angle [deg.]
35 3 σ. Angle [deg.] 5 5 5 Experiment SPH Theory 3 4 5 6 Incident Angle [deg.]
Incident speed [m/sec.] 3 5 5 5 3 4 n > 38 n > 3 n > n > n > 5 θ Angle of incidence [deg.] Angle of incidence [deg.] 4 3 n > 3 n > 38 v 5[m/sec.] n > n > n > 5 5 5 5 3 35 4 45 Tilt angle [deg.]
シミュレーション法の応用分野と 今後の課題 濡れ を伴う粉体系へのシミュレーション法の開発 地盤の液状化現象への応用 R+X 非常に大きな変形をともなう粘弾性体 のシミュレーション ゲルの衝突 水滴の衝突 R+X h h (a) (b) K. Okumura et. al. Europhys. Lett. 6, 37 (3) 高圧の物理学の実験 自由表面を持つ様々な現象への 数値的なアプローチ 円形跳水の問題 Clive Ellegaard et. al., Nonlinearity,, (999) Nature, 39, 3 (998) 年月9日木曜日