レイノルズ輸送定理シミュレーター 戻る
流体力学 / コントロールボリューム

レイノルズ輸送定理シミュレーター

縮小ノズル・90°曲がり管・拡大管のコントロールボリューム解析を実行。連続の式・運動量方程式・ベルヌーイ方程式から流量・速度・圧力・反力をリアルタイム計算します。

管路パラメータ

レイノルズ輸送定理(一般形)

$$\frac{dB_{sys}}{dt}= \frac{\partial}{\partial t}\!\int_{CV}\!\rho b\,dV + \!\int_{CS}\!\rho b(\vec{V}\!\cdot\!\hat{n})\,dA$$

連続: $\dot{m}=\rho A_1 V_1 = \rho A_2 V_2$
運動量: $\sum\vec{F}=\dot{m}(\vec{V}_2-\vec{V}_1)$
圧力(ベルヌーイ): $P_1+\tfrac{1}{2}\rho V_1^2 = P_2+\tfrac{1}{2}\rho V_2^2$

質量流量 ṁ
kg/s
出口流速 V₂
m/s
出口圧力 P₂
kPa
合計反力 |F|
N

レイノルズ輸送定理って何のために使うの?

🧑‍🎓
「レイノルズ輸送定理」って何ですか?流体力学の中でどこに使うのかイメージできなくて…
🎓
ざっくり言うと、「固定した箱(コントロールボリューム)」で流体の出入りを管理するための道具だよ。実際の流れを追うのではなく、箱の入口と出口だけを見ることで、力・エネルギー・質量をシンプルに計算できる。例えばホースの先を細くすると水が速く出るよね? その入口と出口だけを見て「どれだけ速くなるか」を計算するのが連続の式だ。シミュレーターで「縮小ノズル」を選んで、D₁を大きくすると出口速度V₂が上がるのが確認できるよ。
🧑‍🎓
なるほど!「90°曲がり管」にしたら、反力という値が出てきました。なぜ曲がっただけで力が生じるんですか?
🎓
流体は進む方向を変えるために「力が必要」だよね(ニュートンの第2法則)。x方向に流れていた水が90°曲がってy方向に出る場合、x方向の運動量が消えてy方向の運動量が生まれる。この「運動量の変化」を与えたのが管壁からの力で、その反作用が「管壁が受ける反力」なんだ。川が曲がる部分で土手が削れたり、配管の曲がり部がブラケットで固定されているのはこのためだよ。力成分解析タブで Fx と Fy の大きさを見てみて。
🧑‍🎓
「拡大管(ディフューザー)」にすると出口圧力P₂が入口P₁より高くなりましたね。速度が下がると圧力が上がるってことですか?
🎓
そう、それがベルヌーイの原理! $P + \frac{1}{2}\rho V^2 = const$ だから、速度Vが下がると圧力Pが上がる。拡大管(ディフューザー)はこの原理を使って速度エネルギーを圧力エネルギーに変換する装置だ。ターボ機械のディフューザー部分や、換気ダクトの出口部分でよく使われるよ。圧力・速度分布タブで入口・出口の値を見ると、この逆比例の関係が視覚的に確認できる。
🧑‍🎓
ベルヌーイ方程式とレイノルズ輸送定理は何が違うんですか?同じ流れを計算しているように見えるんですが…
🎓
ベルヌーイ方程式はレイノルズ輸送定理の「エネルギー方程式に非粘性・定常・非圧縮・単一流線」という特別な条件を仮定した場合の特殊解だ。このシミュレーターも実はベルヌーイ式でP₂を計算している。でも、レイノルズ輸送定理はもっと一般的で、粘性損失があるときは損失項を足せばいいし、複数の入出口がある場合も対応できる。複雑なポンプや圧縮機の設計では、この一般形が必要になるんだよ。

コントロールボリューム解析の数式

定常状態の連続の式(質量保存):

$$\dot{m} = \rho_1 A_1 V_1 = \rho_2 A_2 V_2$$

非圧縮流れ(ρ=const)では $A_1 V_1 = A_2 V_2$、出口速度 $V_2 = V_1 (A_1/A_2) = V_1 (D_1/D_2)^2$ となります。

定常状態の運動量方程式(x・y 各方向):

$$\sum F_x = \dot{m}(V_{2x} - V_{1x}) + P_2 A_2 n_{2x} + P_1 A_1 n_{1x}$$

90°曲がり管では $V_{1x}=V_1, V_{1y}=0, V_{2x}=0, V_{2y}=V_2$ として Fx・Fy を別々に求め、合力 $|F|=\sqrt{F_x^2+F_y^2}$ を計算します。

実際の工学への応用

ロケットエンジン設計:ノズルから高速で排出されるガスの運動量変化から推力を計算します。シミュレーターの縮小ノズルで「推力」を直感的に理解し、ロケット性能の基礎を学べます。

配管支持設計:工場の配管系統で、曲がり管・T字管・バルブ部には運動量変化による力が加わります。この力を正確に計算してサポートブラケットを設計することで、配管の振動・破損を防ぎます。

ダム・スルースゲート:水門を通る流れの力や流量計算に応用されます。ゲート開度と水位から流量と水門にかかる力を計算し、安全設計に活用します。

よくある質問

コントロールボリューム(CV)は解析対象の空間的な領域で、固定されていることも動くこともあります。コントロールサーフェス(CS)はそのCVの境界面です。レイノルズ輸送定理の体積積分はCV全体について、面積積分(フラックス項)はCS上で行います。このシミュレーターでは、ノズルや曲がり管の入口から出口まで囲んだ固定したCVを使っています。
ベルヌーイ方程式 $P + \frac{1}{2}\rho V^2 = const$ (全圧が一定)から説明できます。流量が保存されるため断面積が小さくなると流速が上がります。この速度の増加は運動エネルギーの増加を意味し、エネルギー保存の観点から圧力エネルギーが運動エネルギーに変換されます。このシミュレーターでD₂をD₁より小さくすると、V₂>V₁, P₂<P₁となることを確認できます。
このシミュレーターはベルヌーイ方程式(無損失)で計算しているため理論的な最大回復量を示しています。実際には、急拡大すると流体の剥離が起こり(境界層剥離)、乱流渦による損失が生じます。ディフューザーの半開き角度が通常7〜8°以下に制限されるのはこのためです。急角度では剥離により圧力回復が大幅に低下します。実際のCFD解析や設計では、この損失を考慮するために補正係数(圧力回復係数Cp)を使用します。
マッハ数(流速/音速)が0.3以下の場合、密度変化が約4%以下で非圧縮近似が有効です。水などの液体はほぼ常に非圧縮です。このシミュレーターは非圧縮流体を仮定しており、水(ρ≈1000 kg/m³)や低速の空気(ρ≈1.2 kg/m³)には正確です。高速ガス流(超音速ノズルなど)では圧縮性効果が重要になり、連続の式にも密度変化を含める必要があります。
一般に D₁=D₂(同径)の場合、x方向の力 Fx = ṁV₁ + P₁A₁ および y方向の力 Fy = ṁV₂ + P₂A₂ となります。入口と出口の径・速度・圧力が等しければ Fx=Fy で、合力方向は45°になります。径が変わる場合や一方の圧力が高い場合は非対称になります。このシミュレーターの力成分解析タブで D₁、D₂、P₁を変えたときの Fx、Fy の変化を確認できます。
レイノルズ輸送定理をエネルギー(B=E)に適用した定常エネルギー方程式は、入出口の全エネルギーの差がポンプ仕事・タービン仕事・損失ヘッドに等しいことを示します:$(P_1/\rho g + V_1^2/2g + z_1) = (P_2/\rho g + V_2^2/2g + z_2) + h_L - h_{pump}$。損失なし(hL=0)・ポンプなし(hpump=0)・水平管(z₁=z₂)の特殊ケースがベルヌーイ方程式です。このシミュレーターはベルヌーイ近似(損失なし)を使用しています。