HOME

粘性抵抗と慣性抵抗

戻る

 地上で観測される色々な現象,例えば降雨やパラシュートによる安全な降下などを半定量的に理解するために、空気,即ち流体の抵抗を考える必要がある。そこで,ここでは空気抵抗の2つの原因を考察して、ニュートン力学が現実の世界を説明出来ている事を確認しよう。

粘性抵抗

 摩擦のない完全流体を除くと、流速の分布に伴う粘性抵抗が働き、粘性力は速さ v の勾配に比例する。
 図1のように、粘性を持つ川の水の流れを考えよう.経験的に、流体と固体(川岸)の接する面では流速が0となり、川の中央ほど速くなることが知られている.そのため、dy だけ離れた流線間の速さに dv = v(y + dy) – v(y) の差が生ずる。dy に対する流線間の速度差 dv の比(速度勾配)

に比例して水分子間の衝突による粘性(摩擦)が大きいと考えられる。川岸では常に速度が0なので、速度勾配も中央の速度 v に比例する。粘性抵抗 F と速度勾配の比例係数を粘性係数 η と定義すると、流線間の単位接触面積あたりの粘性力 τ は
(1)
と書き表せる.粘性係数の単位は、(1) 式から、


で与えられる.ここで、Pa = N/m2 は圧力の単位で、103 hPa = 105 Pa が約1気圧に相当する.
 一様な流速 v 中に置かれた半径 a の球の受ける粘性抵抗 F は、近似的に、球と流体が接する面積 S に粘性応力 τ を掛けて得られる:
=cv,
c=4πaη、正確な計算によるとaηv).
 (2)
ここで、球の表面に接する流体の流速 (v = 0) が、球からその半径 a 程度離れれば一様な流れの速度 v になると仮定し、球の表面の速度勾配を (v–0)/a で近似した。この場合、粘性抵抗は速さ v に比例し、20度C、1気圧の空気の粘性係数は約1.8x10-5 [Pa・s]で、水の粘性係数は、それより2桁近く大きく、0.80x10-3 [Pa・s]である。

慣性抵抗

 速度が遅い場合には、速度に比例する粘性抵抗が効果的であるが、ある程度速くなると(形状に依存)速度の2乗に比例する慣性抵抗が主になる.図2に示す速度 v で落下する半径 l のパラシュートを考えよう.このパラシュートが落下するにつれて、その下方の断面積 S (=πl2) の空気が排除される.この効果を、断面積 S の空気を下方に速度 v に加速すると理解出来る.即ち、単位時間の間に加速する空気の体積 ΔVtSv、その質量 ΔM は空気の密度をρとして ΔM/Δt= ρSv となる.結果として、空気は dt 時間当り ΔM = ρSvΔtv を掛けた運動量 dP/dt を得る。結局、運動方程式dP/dt= Fより、パラシュートの受ける慣性抵抗の大きさは F= ΔP/ΔtSv2 = Dv2 (D = ρS) になる。方向は「摩擦」なので、常に –v の方向。 この表式は近似式なので、正確にはこれに形状に依存した1程度の大きさの係数がかかる:D = kρS

粘性抵抗がある場合の運動

 速度に比例する粘性抵抗があるときには、運動方程式は速度に比例する減衰項を含み
  (3)
と書ける.t0 = m/c [Kg•m/N/s=s] は、慣性の大きさと摩擦による減速の比率で決まる時定数(速度が 1/e に減速されるに要する特性時間)である.慣性質量が大きければ摩擦に負けないため、減速までの時間が長くなる事に対応する。両式を時間で積分して、

が得られる.ここで、初期条件、t = 0 で v = v0 を代入すると積分定数が A = v0x, B = v0y + mg/c と決まり、

となる.これらの結果から、初速度は時定数 t0 で減衰し、落下速度は同じ時定数で終端速度に漸近する。

 時刻 t における質点の位置は、速度をもう一度時間で積分して、
 
  (5)
 
が得られる.初期条件、t = 0 x = y = 0 を代入し、 Cx1 = mv0x/c, Cy1 = mv0y/c + m2g/c2 から、
  (6)
となる.y の第1項は、初速度による到達距離、第2項は、摩擦のある場合の自由落下距離。

慣性抵抗がある場合の運動

 速度の2乗に比例する慣性抵抗がある場合の質点の運動方程式は、
  (7)
と書ける.(7) の両式は、vx, vy の両成分を含む連立微分方程式になる.ここでは、簡単のために (1) v = vx
(摩擦のある慣性運動、常に v > 0)或いは (2) v = vy(摩擦のある自由落下、常に v < 0)の場合の2通りを考えよう.

(8)
を時間で積分して、

(9a)
が得られる.(2) の場合は常に v < 0 で、v2mg/D(終端速度の2乗)が成立ち、

(9b)
となる.初期条件、t=0 で (1) の場合は、v=v0, (2) の場合は、v=0 を代入して得た C1=1/v0, A=1 を代入し、

(10)
を得る.また、(2) の場合の v の終端速度はで与えられ、(8) 式の(2) でdv/dt = 0と置いて得られる結果と等しいことが確認できる.この場合の特性時間は[秒]である.この時間の数倍でほぼ終端速度になると考えて良い.
 さて、到達距離 y は、(10) 式を更に時間で積分して、

(11)
が得られる.ここで、lnA、lnB は積分定数.ここで、初期条件、t=0 で (1) y=0, (2) y=h0 を代入して得た A=Dv0/m , B= exp(-Dh0/m) を代入し、

(12)
を得る.面白いことに、速度に比例する粘性抵抗の場合には、 x 成分は有限の距離までしか飛ばないが、速度の2乗に比例する慣性抵抗の (1) 場合には、速度が小さいところで限りなくゆっくり減衰するため極限値がないが、実際には粘性抵抗の寄与が効いてくる.また、(2) 場合には、(a) から分かるように、t→∞ の極限では y = h0vƒt + (vƒ2/g) ln 2 となるので、-vf t に比例して降下する.