ベルヌーイの定理とベンチュリ効果

管を細くして、流速が上がり圧力が下がるのを観察しましょう。やがて水は、沸騰し始める限界点に達します。

インタラクティブシミュレーションを読み込んでいます...

レッスン

理論 — ベルヌーイの定理とベンチュリ効果

管に関する二つの基本原理。このページのすべては、そこからの帰結です。連続の式:入った分だけ出ていく。ゆえに A₁v₁ = A₂v₂ であり、管が細くなるほど水は速く流れます。ベルヌーイの定理:流線上では、単位体積あたりの圧力と運動エネルギーの和が一定に保たれます。つまり p + ½ρv² = constant となり、流速が上がれば圧力は下がります。何かが圧力を押し下げているわけではありません。エネルギーがすでに速度へと使われたからです。

各記号の意味

A
管の断面積
v
そこを流れる流体の速度
ρ
流体の密度
前提
単一の流線に沿う、非圧縮かつ非粘性の定常流。ここで最も重要になるのが粘性です。実際のベンチュリ管では下流で圧力が完全には回復せず、その差分こそが摩擦によって永久に失われた圧力となります。そしてそれこそが、このメーターを稼働させるための代償なのです。

流速は面積比の2乗で、圧力は流速の2乗で効いてくる

入口を 50 mm、スロート部を 25 mm に設定し、水が 2 m/s で流れているとします。スロート部の流速の行には 8.00 m/s、圧力変化の行には −29.95 kPa と表示されます。前者を決定しているのは連続の式です。四分の一の面積を同じ体積が通過するには四倍の速さで移動しなければならず、流速比の行は 50/25 の2乗である 4.00× を示しています。残りを引き受けるのがベルヌーイの定理です。その運動エネルギーの項には v² が含まれるため、四倍の速さは十六倍の動圧を意味します。

今度はスロート部の入力欄に 10 と入力してみてください。流速比は 25.00×、スロート部の流速は 50.00 m/s、圧力変化は −1245.75 kPa になります。入口はスロート部の二倍ではなく五倍になり、圧力低下は先ほど読み取った値の 41.6 倍に達します。ここからは計算を追うよりも、二つ目のパネルを見る方が視覚的に理解しやすいでしょう。そこでは入口と流体を固定し、スロート部だけを右に向かって細く変化させています。曲線は広いスロート部ではほぼ平坦ですが、あるところから急激に落ち込みます。

計算上は成り立っても、水はそこまでついていかない

その −1245.75 kPa は周囲圧の 101.325 kPa から引かれるため、スロート部の絶対圧力は −1144.43 kPa になります。いかなる流体であっても絶対圧力がマイナスになることはありませんし、水はゼロになるよりずっと手前で限界を迎えます。20 °C であれば 2339 Pa で沸騰するため、両方のパネルの下部にある網掛けの帯はまさにその境界線に引かれています。圧力の曲線がそこへ入り込み、赤く変わるのを確認してください。

物理的にあり得ないという印がついていても、行にはそのまま数値が表示されます。その数値こそが学びになるからです。また、このツールは逆側から境界を解き明かします。入口が 50 mm、水が 2 m/s、周囲圧が海面気圧のとき、このスロート径未満でキャビテーション発生 の行には 18.75 mm と表示されます。これ以上細く絞り込むと、液体は細くなる代わりにちぎれてしまいます。スロート部の入力欄に 19 と打ち込み、続いて 18 と入れてみてください。この二つの数字の間で、判定が切り替わることがわかります。

10 m の水柱は、ほぼ 1 気圧を使い果たす

高さを 10 m 上げるだけで大気圧のほぼすべてが失われます。管が細くなっていようがいまいが関係ありません。両方の直径を 50 mm に設定し、スロート部が入口からどれだけ上にあるかを決める高さのスライダーを +10.0 m まで動かしてみてください。圧力変化の行は −97.92 kPa になりますが、これはすべて ρgh の結果です。密度 998.2 kg/m³ の水、g = 9.81 m/s²、そして 10 m という高さです。大気がもたらす 101.325 kPa に対して、管は完全に真っ直ぐなままであるにもかかわらず、高さだけでその 96.6% を使い果たしてしまいました。スロート部の絶対圧力は 3.40 kPa を示し、これは蒸気圧のラインを 1 kPa ほど上回るだけです。

スライダーは 10 m で止まりますが、水も 10.11 m で止まります。残された圧力を ρgh が使い切ってしまうからです。吸引ポンプは自ら水を持ち上げているわけではありません。水面の上の圧力を下げることで、大気に水柱を押し上げさせているのです。ポンプが残せる最も低い圧力は、二つ前の段落で触れた 2339 Pa にすぎません。そこまでいくと水はそれ以上上がる代わりに沸騰してしまいます。つまり 10 m というのが与えられた予算のすべてであり、ポンプの性能が試されているわけではないのです。

全プロセスの詳細解説

  1. ベンチュリ管で流量を読み取る仕組みと、決して測れない流量 8 ステップ

    50 mm の水道管に、25 mm のスロート部を持つベンチュリ管が組み込まれています。標準的なベンチュリ管 を押して設定してください。二つの取り出し口をつなぐ差圧計は 29.95 kPa を示しています。現在の流量はいくらでしょうか。また、この計器が示せる最大の流量はどれくらいでしょうか。

    この問題をツールで開く

    1. 毎秒入口を通過する体積はすべてスロート部も通過します。そのため、面積が減少した分だけ流速が上がり、両者の積は一定に保たれます。面積は直径の2乗に比例するため、内径が半分になっても流速は二倍にはならず、四倍になります。流速比の行は 4.00×、スロート部の流速は 8.00 m/s を示しています。

      A1A2=(d1d2)2=(5025)2=4.00v2=v1(d1d2)2=2×4.00=8.00 m s−1\frac{A_1}{A_2}=\left(\frac{d_1}{d_2}\right)^{2}=\left(\frac{50}{25}\right)^{2}=4.00 \qquad v_2=v_1\left(\frac{d_1}{d_2}\right)^{2}=2\times4.00=8.00\ \mathrm{m\,s^{-1}}
    2. ベルヌーイの定理は、その流速の代償を圧力で支払います。管は水平なので二つの ρgh の項は打ち消し合い、運動エネルギーの項だけが残ります。運動エネルギーの項には v² が含まれるため、四倍の速さを生み出すには入口の動圧の十六倍が必要になります。そのうちの一倍分はすでに使われていたため、実際の圧力低下は ½ρv₁² の十五倍です。ここでの ½ρv₁² は 1996.4 Pa に当たります。圧力変化の行は −29.95 kPa を示しています。

      Δp=12ρ(v12−v22)=12(998.2)(22−82)=−29,946 Pa\Delta p=\tfrac{1}{2}\rho\left(v_1^{2}-v_2^{2}\right) =\tfrac{1}{2}(998.2)\left(2^{2}-8^{2}\right)=-29{,}946\ \mathrm{Pa}
    3. 入口側の取り出し口における大気圧からその値を引くと、スロート部は絶対圧力で 71.38 kPa になります。まだ余裕はあります。20 °C の水は 2339 Pa で沸騰するため、69 kPa ほどの余裕が残されているわけです。スロート部のマッハ数は 0.005 を示しており、ベルヌーイの定理に組み込まれたもう一つの仮定もまた、まったく無理がないことを報告しています。

      p2=pamb+Δp⇒101.325−29.946=71.379 kPaM=v2c=8.001481=0.005p_2=p_{\mathrm{amb}}+\Delta p\quad\Rightarrow\quad 101.325-29.946=71.379\ \mathrm{kPa} \qquad M=\frac{v_2}{c}=\frac{8.00}{1481}=0.005
    4. しかし、ベンチュリ管の働きはこうした計算とは異なります。計器には流速を決めるスライダーなどありません。あるのは二つの取り出し口と差圧計だけであり、流量こそが導き出すべき未知数なのです。そこで、今度は Δp が既知で v₁ を求めるものとして、同じ方程式を逆向きに解いてみましょう。差圧計の 29.95 kPa を代入し直すと、入口の流速は 2.000 m/s となります。これはまさに入口が設定されていた値です。

      Δp=−12ρv12[(d1d2)4−1]⇒v1=2∣Δp∣ρ[(d1/d2)4−1]=2(29,950)(998.2)(15)=4.0005=2.000 m s−1\Delta p=-\tfrac{1}{2}\rho v_1^{2}\Bigl[\Bigl(\frac{d_1}{d_2}\Bigr)^{4}-1\Bigr] \quad\Rightarrow\quad v_1=\sqrt{\frac{2\lvert\Delta p\rvert}{\rho\bigl[(d_1/d_2)^{4}-1\bigr]}} =\sqrt{\frac{2(29{,}950)}{(998.2)(15)}}=\sqrt{4.0005}=2.000\ \mathrm{m\,s^{-1}}
    5. これに入口の面積 π(0.050)²/4 = 1.9635×10⁻³ m² を掛ければ、流量の行にある 3.93 L/s は流速からではなく圧力の読み取り値から得られたことになります。計器のデータシートでは、これら二つのステップを直径比 d₂/d₁(ここでは 0.5)を表す β を用いて一つの方程式にまとめています。この装置が存在する理由こそがその方程式であり、このページ上でツールが決して計算しない唯一のものなのです。

      Q=v1A1=2.000×1.9635×10−3=3.927×10−3 m3 s−1Q=A22∣Δp∣ρ(1−β4)β=d2d1=0.5Q=v_1A_1=2.000\times1.9635\times10^{-3}=3.927\times10^{-3}\ \mathrm{m^{3}\,s^{-1}} \qquad Q=A_2\sqrt{\frac{2\lvert\Delta p\rvert}{\rho\left(1-\beta^{4}\right)}}\quad\beta=\frac{d_2}{d_1}=0.5
    6. どんな計器にも上限がありますが、この計器の上限は差圧計ではなく水そのものに依存しています。スロート部を絞っていくと、そこでの圧力は直径比の4乗で下落していき、やがて 2339 Pa に達します。その時点で液体は加速する代わりに沸騰を始めます。現在の流量でそれが起こるスロート部を逆算してみましょう。括弧内の値は 50.582 となり、その4乗根は 2.6669 です。このスロート径未満でキャビテーション発生 の行は 18.75 mm を示します。組み込まれた 25 mm のスロート部は、それよりも三分の一ほど広いことになります。

      d2∗=d1[1+2(pamb−pv)ρv12]−1/4=d1(50.582)−1/4502.6669=18.75 mmd_2^{\ast}=d_1\Bigl[1+\frac{2(p_{\mathrm{amb}}-p_v)}{\rho v_1^{2}}\Bigr]^{-1/4} =d_1(50.582)^{-1/4}\qquad\frac{50}{2.6669}=18.75\ \mathrm{mm}
    7. 今度はスロート部を 25 mm に固定して、同じ境界を逆側から見てみましょう。流体が消費できる圧力は周囲圧から蒸気圧を引いた 98,986 Pa であり、それを消費するのは ρ[(d₁/d₂)⁴ − 1] です。この流速以上でキャビテーション発生 の行は 3.64 m/s を示しています。

      v1max⁡=2(pamb−pv)ρ[(d1/d2)4−1]=2(98,986)(998.2)(15)=13.222=3.636 m s−1v_1^{\max}=\sqrt{\frac{2\left(p_{\mathrm{amb}}-p_v\right)}{\rho\bigl[(d_1/d_2)^{4}-1\bigr]}} =\sqrt{\frac{2(98{,}986)}{(998.2)(15)}}=\sqrt{13.222}=3.636\ \mathrm{m\,s^{-1}}
    8. その流速を流量に変換してくれる行はないため、自分で計算しましょう。50 mm の内径を 3.636 m/s で流れる場合、流量は 7.14×10⁻³ m³/s、すなわち 7.14 L/s となります。この計器のフルスケールは現在読み取っている流量の 1.82 倍であり、発生し得る最大の差圧は 98.99 kPa です。これは周囲圧から蒸気圧を引いたものであり、これを超えて利用できるものはありません。これ以上押し込むとスロート部は蒸気で満たされ、差圧は流量と連動しなくなり、差圧計は単に沸騰した空洞を測っているだけになります。

      Qmax⁡=v1max⁡A1=3.636×1.9635×10−3=7.14×10−3 m3 s−17.143.927=1.818Q_{\max}=v_1^{\max}A_1=3.636\times1.9635\times10^{-3}=7.14\times10^{-3}\ \mathrm{m^{3}\,s^{-1}} \qquad\frac{7.14}{3.927}=1.818

    解答

    現在は 3.93 L/s であり、最大で 7.14 L/s です。 ツールはそのどちらも表示しませんが、設計を左右するのは後者の値です。

    緩やかな絞り を押してスロート部を 40 mm に広げてみてください。この流速以上でキャビテーション発生 の値は 3.64 から 11.73 m/s(23.0 L/s に相当)へと上昇しますが、圧力変化は −29.95 kPa から −2.88 kPa へと激減します。どちらの数字も (d₁/d₂)⁴ − 1 という同じ括弧から導かれています。一方は分子に、もう一方は分母の平方根の中に位置しているため、測定範囲は手放した信号の逆数平方根としてしか広がりません。つまり、差圧を十分の一に減らすことで、3.2 倍の流量を測れるようになるわけです。β を選ぶということはまさにこのトレードオフであり、決してキャビテーションを起こさないほど広いスロート部というのは、差圧計にわずか数百パスカルの違いを読み取らせるようなものなのです。

    Qmax⁡=A12(pamb−pv)ρ[(d1/d2)4−1]Q_{\max}=A_1\sqrt{\frac{2\left(p_{\mathrm{amb}}-p_v\right)}{\rho\bigl[(d_1/d_2)^{4}-1\bigr]}}

学習の道すじ

流体の運動:粘性が奪い去るもの