連続: \(\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\)
流路形状を縮小ノズル・90°曲がり管・拡大管から選ぶと、連続の式・運動量方程式・ベルヌーイ式が同時に解かれ、流量・速度・圧力・反力がその場で更新されます。配管設計の検討に使える流体力学シミュレーター。
連続: \(\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\)
定常状態の連続の式(質量保存):
$$\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字管・バルブ部には運動量変化による力が加わります。この力を正確に計算してサポートブラケットを設計することで、配管の振動・破損を防ぎます。
ダム・スルースゲート:水門を通る流れの力や流量計算に応用されます。ゲート開度と水位から流量と水門にかかる力を計算し、安全設計に活用します。
レイノルズ輸送定理シミュレーターの物理モデルでは、縮小ノズル・90°曲がり管・拡大管の各要素に対し、固定されたコントロールボリュームを設定し、非定常項を無視した定常流れを仮定します。連続の式は質量保存則に基づき、入口と出口の断面積 \(A_1, A_2\) と平均流速 \(V_1, V_2\) の関係を \(A_1 V_1 = A_2 V_2\) で与えます。運動量方程式は、コントロールボリュームに作用する表面力(圧力と壁面反力)と体積力を考慮し、流出・流入する運動量流束の差と釣り合います。例えば曲がり管では、流れの方向変化により壁面に作用する反力 \(F\) が生じ、その成分は \(F_x = \rho Q (V_{2x} - V_{1x}) + (p_1 A_1 - p_2 A_2 \cos\theta)\) のように計算されます。ベルヌーイ方程式は粘性損失を無視した理想流れに適用し、ノズルや拡大管の圧力変化を \(p_1 + \frac{1}{2}\rho V_1^2 = p_2 + \frac{1}{2}\rho V_2^2\) で評価します。これらの式を連立することで、流量 \(Q\)、各断面の速度・圧力、および管壁が受ける反力をリアルタイムに出力します。
産業での実際の使用例
自動車業界では、エンジンの吸気マニホールドや排気系の触媒コンバーター前後における縮小・拡大管の解析に活用。例えばトヨタのハイブリッド車用エンジンでは、本シミュレーターで流量と圧力損失をリアルタイム計算し、過給機とのマッチングを最適化。また、化学プラントの配管設計(例:三菱ケミカルの反応器出口配管)では、90°曲がり管での流体反力を事前評価し、支持構造物の強度設計に反映している。
研究・教育での活用
大学の流体力学実験では、ベルヌーイの定理と運動量方程式の検証教材として使用。東京工業大学の学部講義では、縮小ノズル内の速度変化と圧力低下を学生が対話的に確認し、理論式とシミュレーション結果を比較する演習に採用。また、研究用途では、マイクロ流路内の非定常流れ解析の予備検討として、本ツールで流量と反力のオーダーを把握してから詳細CFDに進む手法が確立されている。
CAE解析との連携や実務での位置付け
本シミュレーターは、詳細な3D-CFD解析(例:ANSYS Fluent)の前段階として位置付けられる。設計初期段階でコントロールボリューム解析により基本諸元(流量・圧力・反力)を瞬時に算出し、設計変更の方向性を絞り込む。実務では、配管レイアウトの変更可否判断や、バルブ選定時の圧力損失見積もりに使用。CAE解析のインプット条件(境界条件や初期値)を本ツールで生成することで、解析時間を従来比で約30%短縮した事例もある。
「ベルヌーイの定理は流線に沿って成立する」という前提を忘れがちですが、曲がり管や急拡大・縮小部では流れが剥離・旋回し、エネルギー損失が生じるため、実際はベルヌーイ式をそのまま適用できません。特に90°曲がり管では遠心力による圧力分布の偏りが生じるため、入口と出口の圧力差だけから速度を推定すると誤差が大きくなります。また、「連続の式は非圧縮性流体なら常に断面積比で速度が決まる」と思い込みがちですが、実際には密度変化や温度変動がある場合や、キャビテーションが発生する条件下では体積流量が保存されず、計算が狂う点に注意が必要です。さらに、「運動量方程式で求めた反力は壁面にかかる力と常に等しい」と誤解されがちですが、実際にはコントロールボリューム内の流体の加速度項や、出入口での圧力・運動量流束のバランスを正しく考慮しないと、反力の向きや大きさを誤るため、特に曲がり管では壁面に作用する力の方向を直感的に逆に捉えないよう注意が必要です。
鋼製配管で入口径D1=50mm、出口径D2=30mm、入口流速V1=4m/s、入口圧力P1=150kPaの場合:連続の式により出口流速V2=(π×50²/4)/(π×30²/4)×4≈11.1m/sに加速。ベルヌーイ式P1+0.5ρV1²=P2+0.5ρV2²を適用すると、密度1000kg/m³の水流で圧力低下ΔP≈61kPa発生。出口圧力P2≈89kPaとなり、キャビテーション防止に注意が必要な領域です。