ボルダ・カルノー方程式

ボルダ・カルノー方程式（テンプレート:Lang-en-short）は、流体力学において流れの急拡大による流体のエネルギー損失を表す経験式である。損失によって全水頭がどの程度減少するかを記述する。これは、全水頭が流線に沿って一定である（不可逆損失のない）散逸のない流れに対するベルヌーイの式とは対照的である。この名称はジャン＝シャルル・ド・ボルダとラザール・カルノーにちなんで名付けられた。

この式は開水路の流れおよび配管内の流れの両方に適用できる。不可逆的なエネルギー損失が無視できる流れの部分ではベルヌーイの定理が代わりに使用される。

定式化

ボルダ・カルノー方程式は次式である：

$\mathrm{\Delta }E=\xi \frac{\rho }{2}(v_1-v_2{)}^2,$

ただし

• テンプレート:Math は流体のエネルギー損失
• テンプレート:Mvar は無次元の経験的な損失係数。テンプレート:Math の値をとる。
• テンプレート:Mvar は流体の密度
• テンプレート:Math は拡大前と拡大後の平均流速

大きく急拡大する場合、損失係数 テンプレート:Mvar は 1 になる。その他の場合、損失係数は他の手段、一般的には実験で得られたデータに基づく経験式で決定する必要がある。ボルダ・カルノー方程式は速度が減少する（テンプレート:Math）場合にのみ有効である。それ以外の場合、損失 テンプレート:Math は 0 である。つまり追加の外力による機械的仕事がなければ、流体のエネルギーの増加はあり得ない。

損失係数 テンプレート:Mvar は流線の影響を受けることがある。例えばディフューザーを使用して配管径を徐々に拡大するとエネルギー損失を減らすことができる。

全水頭とベルヌーイの定理との関係

ボルダ・カルノー方程式により、ベルヌーイの式の定数は減少する。非圧縮性流れの場合、結果は次のようになる：位置1の下流に位置2がある場合、流線に沿って

$p_1+\frac{1}{2}\rho v_1^2+\rho g{z}_1=p_2+\frac{1}{2}\rho v_2^2+\rho g{z}_2+\mathrm{\Delta }E,$

ただし

• テンプレート:Math は位置1, 2の圧力
• テンプレート:Math は（ある基準位置から測った）高さ
• テンプレート:Mvar は重力加速度

両辺にある最初の 3 つの項はそれぞれ圧力、流体の運動エネルギー密度、重力による位置エネルギー密度である。圧力は実質的にポテンシャルエネルギーの形態として作用することが分かる。

高圧配管の流れの場合、重力の影響を無視でき、テンプレート:Math は テンプレート:Math に等しくなる。すなわち

$\mathrm{\Delta }E=\mathrm{\Delta }(p+\frac{1}{2}\rho v^2).$

開水路流れの場合、テンプレート:Math は総損失水頭 テンプレート:Math と次のように関係する：

$\mathrm{\Delta }E=\rho g\mathrm{\Delta }H,$

ただし テンプレート:Mvar は全水頭

$H=h+\frac{v^2}{2g},$

テンプレート:Mvar は水頭（基準面から水面までの高さ）：$h=z+p/(\rho g)$。

例

急拡大する配管

急拡大する配管

ボルダ・カルノー方程式は、水平配管の急拡大部を通る流れに適用される。拡大部の上流にある断面1での平均流速を テンプレート:Math、圧力を テンプレート:Math、断面積を テンプレート:Math とする。断面2（拡大部およびテンプレート:仮リンクの領域よりかなり下流）での対応する量をそれぞれ テンプレート:Math および テンプレート:Math とする。拡大部では流れが剥離し、エネルギー損失を伴う乱流循環流領域が発生する。この急拡大による損失係数 テンプレート:Mvar はほぼ 1 に等しい：テンプレート:Math。流体の密度 テンプレート:Mvar を一定と仮定すると、質量保存則により断面1と断面2のそれぞれを通る体積流量は等しくなる：

$A_1{v}_1=A_2{v}_2$     すなわち     $v_2=\frac{A_1}{A_2}{v}_1.$

その結果、ボルダ・カルノー方程式により、この急拡大におけるエネルギー損失は次のようになる：

$\mathrm{\Delta }E=\frac{1}{2}\rho {(1-\frac{A_1}{A_2})}^2{v}_1^2.$

対応する全水頭損失 テンプレート:Math は次のようになる：

$\mathrm{\Delta }H=\frac{\mathrm{\Delta }E}{\rho g}=\frac{1}{2g}{(1-\frac{A_1}{A_2})}^2{v}_1^2.$

テンプレート:Math の場合、2か所の断面間の運動エネルギーの総変化が消散する。その結果、両方の断面間の圧力変化は次のようになる (重力の影響がない水平配管の場合)：

$\mathrm{\Delta }p=p_1-p_2=-\rho \frac{A_1}{A_2}(1-\frac{A_1}{A_2})v_1^2,$

そして水頭 テンプレート:Math の変化は次のようになる：

$\mathrm{\Delta }h=h_1-h_2=-\frac{1}{g}\frac{A_1}{A_2}(1-\frac{A_1}{A_2})v_1^2.$

ここで右辺の先頭にあるマイナス記号は圧力（および水頭）が配管の拡大後に大きくなることを意味する。配管拡大の直前と直後の圧力（および水頭）のこの変化がエネルギー損失に対応することは、ベルヌーイの定理の結果と比較すると明らかになる。散逸のないこの定理によれば、流速の低下はエネルギー損失を伴う今回のケースで見られるよりもはるかに大きな圧力の増加と関連している。

配管の急収縮

配管径の急収縮部を通る流れ。断面3付近に剥離泡が発生する。

配管径が急収縮する場合、流れは急に曲がる配管形状に沿いながら細い配管に流れ込むことができない。その結果、流れの剥離が生じ、細い配管の入口に再循環する剥離領域が形成される。主流は、剥離した流れ領域の間で収縮し、その後再び拡大して配管領域全体に広がる。

収縮前の断面1と主流が最も収縮する断面3（テンプレート:仮リンク部）との間では水頭損失はあまりない。しかし、断面3から断面2への流れの拡大では、かなりの損失がある。これらのヘッドロスは収縮係数 テンプレート:Mvar を使用してボルダ・カルノー方程式で表すことができる：

$\mu =\frac{A_3}{A_2},$

ただし テンプレート:Math は主流が最も収縮する位置3の断面積、テンプレート:Math は配管の狭い方の断面積である。テンプレート:Math なので、収縮係数は 1 未満である：テンプレート:Math。ここでも質量保存則が成り立つため、3か所の断面の体積流束は一定である（流体の密度 テンプレート:Mvar が一定の場合）：

$A_1{v}_1=A_2{v}_2=A_3{v}_3,$

ただし テンプレート:Math は対応する断面の平均流速である。ボルダ・カルノー方程式（損失係数 テンプレート:Math）により、流体の単位体積あたりのエネルギー損失 テンプレート:Math は配管の収縮によって次のように表される：

$\mathrm{\Delta }E=\frac{1}{2}\rho {(v_3-v_2)}^2=\frac{1}{2}\rho {(\frac{1}{\mu }-1)}^2{v}_2^2=\frac{1}{2}\rho {(\frac{1}{\mu }-1)}^2{\left(\frac{A_1}{A_2}\right)}^2{v}_1^2.$

対応する全水頭損失 テンプレート:Math は テンプレート:Math で計算できる。

ワイスバッハの測定によれば、鋭い角をもつ収縮部の収縮係数はおおよそ次の通りである：

$\mu =0.63+0.37{\left(\frac{A_2}{A_1}\right)}^3.$

急拡大に対する運動量バランスからの導出

配管の急拡大については、上図に示すように、質量と運動量の保存則からボルダ・カルノー方程式を導くことができる。断面積 テンプレート:Mvar を通過する（配管軸に平行な運動量成分に対する）運動量流束 テンプレート:Mvarは、オイラー方程式によれば次のようになる：

$S=A(\rho v^2+p).$

拡大部のすぐ上流の断面1と、流れが（拡大部で剥離した後）配管壁に再び付着する場所の下流の断面2、および配管壁によって囲まれたコントロールボリュームの質量と運動量の保存について考える。流入部によるコントロールボリュームの運動量増加 テンプレート:Math があり、流出部による運動量減少 テンプレート:Math がある。さらに拡大部の壁 (配管軸に垂直) によって流体に及ぼされる圧力による力

$F=(A_2-A_1)p_1$

の寄与もある。ここで圧力は近くの上流圧力 テンプレート:Math に等しいと仮定している。

これらの寄与を加えるとコントロールボリュームの運動量バランスは次のようになる：

$S_1-S_2+F=A_1(\rho v_1^2+p_1)-A_2(\rho v_2^2+p_2)+(A_2-A_1)p_1=0.$

その結果、質量保存則より テンプレート:Math となるため、

$\mathrm{\Delta }p=p_1-p_2=-\rho (\frac{A_1}{A_2}{v}_1^2-v_2^2)=-\rho \frac{A_1}{A_2}(1-\frac{A_1}{A_2})v_1^2,$

は上記の例で現れる圧力降下 テンプレート:Math と一致する。

エネルギー損失 テンプレート:Math は

$\begin{array}{cc}\mathrm{\Delta }E& =(p_1+\frac{1}{2}\rho v_1^2)-(p_2+\frac{1}{2}\rho v_2^2)\\ & =\mathrm{\Delta }p+\frac{1}{2}\rho (v_1^2-v_2^2)\\ & =-\rho v_2(v_1-v_2)+\frac{1}{2}\rho (v_1^2-v_2^2)\\ & =\frac{1}{2}\rho {(v_1-v_2)}^2\end{array}$

となる。これはボルダ・カルノー方程式（テンプレート:Math）である。

参考文献

• テンプレート:Citation, 634 pp.
• テンプレート:Citation, 634 pp.
• テンプレート:Citation, 706 pp.

関連項目

• ダルシー・ワイスバッハの式
• en:Prony equation