【機械設計マスターへの道】運動量の法則[流体力学の基礎知識⑤]
煩雑な数値計算によらずに比較的簡単に流れの現象を把握することができて、流体が管や板などの様々な形状の物体に及ぼす力を求めることができる「運動量の法則」とその適用例についてご紹介します。
1.ニュートンの法則
まず、力学の基礎となるニュートンの法則について再確認します。
第1法則:慣性の法則
外部から力が作用しない限り、動いている物体は動き続けようとし、止まっている物体は止まり続けようとする性質のこと。
第2法則:運動の法則
物体に力F[N]を加えて物体が運動するとき、
物体の質量をm[kg]、加速度をα[m/s2]とすると、
F = mα ・・・(1)
の関係になる。
第3法則:作用・反作用の法則
物体に力を加える(作用)と、逆方向に同じ大きさの力を受ける(反作用)こと。
2.運動量の法則
質量m[kg]の物体が速度v[m/s]で運動しているとき、運動量Mは、
M=mv [kgm/s] ・・・(2)
で定義されます。
速度の1回微分は加速度ですから、(1)(2)式より
F = m dv/dt =d/dt(mv) =dM/dt ・・・(3)
すなわち、運動量の時間的な変化割合は、その物体に及ぼす力に等しくなります。
この(3)式を運動量の法則といい、質点の集合体とみなすことができる連続体の流体要素にも適用することができます。
検査面と検査体積、運動量の式
簡単にするために、定常流の状態に限定して考えてみます。
図1で、管路を切断する断面①,②と管路内壁面で囲まれる仮想の境界面を「検査面」、検査面で囲まれる領域を「検査体積」、検査体積内の流体要素を「系」といいます。
検査面①を通って検査体積に入る運動量をMin、検査面②を通って検査体積から出る運動量をMoutとすると、系における運動量変化量⊿Mは、
⊿M=Mout-Min ・・・(4)
となります。
断面①,②の断面積、平均流速、密度をそれぞれ、A1,A2,v1,v2,ρ1,ρ2 とします。
断面①から入る流体は、微小時間⊿tの間に距離v1⊿t だけ移動し、流体質量m1=ρ1A1v1⊿tが平均流速v1で入りますから
Min=m1v1=(ρ1A1v1⊿t) v1 ・・・(5)
同様に、断面②から出る運動量は、
Mout= m2 v2=(ρ2A2v2⊿t) v2 ・・・(6)
となります。
ここで各断面を流れる流量をQ1,Q2とすれば、Q1= A1v1、Q2= A2v2であるので、式(3),(4),(5),(6)より、検査体積内の流体に作用する外力Fは、
F=dM/dt=(Mout-Min)/⊿t =ρ2Q2v2-ρ1Q1v1=(ρQv)out-(ρQv)in ・・・(7)
すなわち、定常流において、検査体積内の流体に及ぼす力は、検査面上の断面を通して単位時間に流出入する運動量の変化に等しくなります。
(7)式を運動量の式といい、粘性や圧縮性のある流体に対しても適用することができます。
非圧縮性流体であれば密度ρを一定として扱うことができます。
連続の式より非圧縮性流体の場合は Q= A1v1= A2v2=const. ですから、(7)式は、
F=ρQ(v2-v1) ・・・(8)
と表すことができます。
つまり、質量流量に速度差をかけた値の力が外力として検査体積内の流体に作用するということになります。
ニュートンの第3法則(作用反作用の法則)により、Fと値が同一で、方向が逆向きの力が物体に作用します。
図1で検査体積を囲む物体が管路であれば、管路に-Fの力が反作用としてかかることになります。
(8)式は、Q=A1V1=A2V2 の関係により、図のような A2<A1 の形状のノズルであれば V2>V1 になりますから、Fの方向は図の右向き(+方向)になります。
管路は、Fの反作用として‐Fの力を図の左向き(-方向)に受けます。
検査面の断面積や圧力が異なる場合
検査面の断面積あるいは、検査面に作用する圧力が異なる場合は、運動量による力の他に、圧力による力fpを考慮する必要があります。
図1の場合、fpは検査面上の断面①,②において圧力によって働く表面力で
fp= A1p1– A2p2 ・・・(9)
となります。
(9)式は、流体がもつ圧力により、検査面上の流入出面を通じて管路におよぼす力であり、入出面における圧力と断面積の積が圧力による力になります。
したがって、検査体積内の流体が管路に及ぼす力をfとすれば(8)、(9)式より
f=fp-F= A1p1– A2p2 + ρQ(v1-v2) ・・・(10)
となります。
管路が受ける力fは、圧力による力fpと、流体から受ける運動量による力の反作用‐Fの合計となりますので、(10)式の通り f=fp-F となります。
3.運動量の式の適用例
(1)消防ホースノズル
図2のような入口直径d1=40[mm], 出口直径d2=13[mm]の消防ノズル先端から、入口圧力p1=0.2[MPa]で水が大気中に放水される状態を考えます。
水がノズル内面に囲まれるような検査体積と検査面をとります。検査面の水の流入口を断面①、流出口(放出口=大気圧)を断面②とします。
流量をQ(m3/s)とすれば、「連続の式」(本連載コラム「連続の式とベルヌーイの定理」の回を参照)より
Q= A1v1= A2v2 したがって v1= (A2/ A1) v2 ・・・(11)
ノズル出口は大気圧ですので出口圧力p2=0となります。
ベルヌーイの式より、
v12/2+p1/ρ= v22/2 したがって p1=(ρ/2)( v22– v12) ・・・(12)
(11),(12)式よりv1を消去してv2について解けばv2=20.1[m/s]となります。
ただし、ρ=1000[kg/s](常温水)
A2=(π/4)(d2x10-3)2=1.33 x10-4[m2] A1=(π/4)(d1x10-3)2=1.26 x10-3[m2]
Q= A2v2=1.33 x10-4 x 20.1=2.67×10-3[m3/s](=160リッター毎分)
v1=Q/A1=2.67×10-3/((π/4) (d1x10-3)2 =2.12 m/s (d1=0.04[m])
(10)式より、ノズルが流出する水から受ける力fは、
f= A1p1 +ρQ(v1-v2)= 1.26 x10-3x0.2×106+1000×2.67×10-3x(2.12-20.1) =203[N]となります。
(2)噴流がタンクに及ぼす力とロケット推進の原理
図3のようなタンク側面の底部に小孔があって、水が噴出しているとします。
簡単のため水面高さhは一定に保持されていると仮定します。タンク内面を囲むように検査面を取ります。タンク断面積Aに比較して小孔の断面積aが十分小さい場合、噴流の速度vは、
v=√2gh ・・・(13)
で求めることができます。これを「トリチェリの定理」といいます。
噴流には縮流という現象が生じるため、噴流の断面積は小孔の断面積aよりも小さくなり、流量を求める際に流量係数をかける必要があります。
ここでは簡単のため、流量係数を1とすれば、
流量Q=av ・・・(14)
検査面にかかる圧力は各方向に対して一定ですから、運動量の式より、タンクが噴流から受ける力F(反作用力)は、
F=ρQv=ρav2=2gρah ・・・(15)
このタンクが噴流から受ける力Fを「推力」(スラスト)といいます。
ロケットの推進機構は、この原理に基づいています。
今回ご紹介した運動量の法則は、検査体積を適切にとることで、検査体積の境界面である検査面から流出入する流体に対して、各検査面の断面積、断面の圧力、密度、流量、流速を知ることができれば、(対象とする物体の内部状態がわからなくても)物体が流体から受ける力を求めることができる便利な手法です。
ぜひ設計等の実務でご活用ください。
次回の連載コラムでは、物体まわりの流れとカルマン渦について解説いたします。
(日本アイアール株式会社 特許調査部 S・Y)
<参考文献>
- 西海孝夫 著『図解 はじめて学ぶ 流体の力学』 日刊工業新聞社、2010
※本コラムで基礎を概説した流体力学についてさらに深く学びたい方に、おススメの書籍です。
- 社会人のための流体力学 [※社会人のための機械工学シリーズ 第3回](講師:技術士 近藤孝邦 氏)