ハーゲン・ポアズイユ流れ(円管内層流の速度分布)
十分に発達した円管内の層流では、速度分布が放物線になります。管内径と平均流速(または体積流量)から、最大流速 umax = 2U・各半径位置の流速・圧力損失・壁面せん断応力を求め、速度分布グラフを描画します。
入力
円管の内径 d(半径 R = d/2)。
平均流速 U と体積流量 Q は Q = U·A の関係。どちらで指定しても残りは自動換算します。
断面平均流速。最大流速はこの 2 倍(umax = 2U)になります。
体積流量。平均流速 U = Q/A に換算して計算します。
指定すると、その位置の流速 u(r) を求めます(r = 0 が中心、r = R が管壁)。「×R」を選ぶと r/R の比で入力できます。
入力すると壁面せん断応力 τw を、管長 L も入れると圧力損失 ΔP を計算します。水20°Cなら約 1.0 mPa·s。
粘度 μ と合わせて、ハーゲン・ポアズイユの式から圧力損失 ΔP を計算します。
粘度 μ と合わせてレイノルズ数 Re を計算し、層流かどうかを確認します。水なら約 1000。
計算結果
計算式
直径 d = 2R の円管内を流れが層流・定常で、十分に発達した状態(流速が軸方向に変化しない状態)になると、軸方向の流速 ux は中心からの距離 r だけで決まる放物線状の分布になります。これがハーゲン・ポアズイユ流れです。
速度分布
管壁(r = R)では粘着条件により ux(R) = 0、円管中央(r = 0)で最大流速 umax になります。断面平均流速 U とは次の関係があります。
体積流量と圧力勾配
体積流量 Q は軸方向の圧力勾配 ∂p/∂x に比例し、粘度 μ に反比例します(圧力勾配 ∂p/∂x の値は負)。
長さ L の管での圧力損失 ΔP、壁面せん断応力 τw は次のように書けます。
- $u_x(r)$:半径位置 r での軸方向流速[m/s]
- $u_{\max}$:最大流速(管中心 r = 0)[m/s]
- $U$:断面平均流速[m/s]
- $Q$:体積流量[m³/s]
- $d,\ R$:管内径・管半径(R = d/2)[m]
- $r$:管中心からの距離(0 ≤ r ≤ R)[m]
- $\mu$:流体の粘度[Pa·s]
- $\Delta P$:長さ L での圧力損失[Pa]
- $\tau_w$:壁面せん断応力[Pa]
適用範囲と注意事項
層流・定常で十分に発達した円管内の流れを前提にした計算です。
- 層流のみ:放物線状の速度分布は層流でのみ成り立ちます。円管の臨界レイノルズ数はおよそ Re = 2300 で、これを超えると乱流に遷移し、断面の速度分布は中央が平坦化します(密度 ρ を入力すると Re を確認できます)。
- 発達した流れ:「十分に発達した」流れを前提とします。入口付近の助走区間では速度分布が発達途中で、本ツールの分布とは異なります。
- 流体と管の前提:非圧縮性・ニュートン流体・一定断面の直管を前提としています。
- 最大流速と平均流速:最大流速は平均流速のちょうど 2 倍(umax = 2U)です。流速計測の代表点として、平均流速に等しくなるのは r/R ≒ 0.707 の位置です。
- 含まれない損失:圧力損失 $\Delta P$ は、直管部の粘性による摩擦分だけです。入口・出口、曲がり、弁、管径の変化による損失は含まれず、配管の高低差による圧力変化も含みません。
- 粘度の入力:粘度 $\mu$ は流れている温度での値を入力します。温度で粘度が大きく変わる流体で、管壁と中心の温度が違う場合は、式の前提(粘度が一定)から外れ、速度分布は放物線からずれます。
- 管が短い場合:助走区間の長さはレイノルズ数が大きいほど長くなります。管が短い場合は、全長が発達した流れにならず、$\Delta P$ の計算値と実際の値が一致しません。
参考資料
明解入門 流体力学(第2版)
Amazonで見る
化学工学計算支援ツール