ファノー流れ（摩擦付管路）シミュレーター

断熱・等断面のパイプの中を流れる圧縮性ガスが、壁面摩擦だけでマッハ 1 に近づき「摩擦チョーク」する現象を可視化するツールです。入口マッハ数・管径・管長・摩擦係数を変えて、出口マッハ数と圧力比、チョークまでの余裕をリアルタイムで確認できます。

パラメータ設定

入口マッハ数 M₁

亜音速流れのみ対象（0.05〜0.90）

配管直径 D

配管長さ L

ダルシー摩擦係数 f

ムーディ線図の値。商用鋼管で約0.015〜0.030

比熱比 γ

空気=1.40、天然ガス≈1.31、CO₂≈1.30

計算結果

—

入口 fL*/D

—

利用可能 fL/D

—

出口 fL*/D

—

出口マッハ数 M₂

—

圧力比 P₂/P₁

—

チョーク判定

—

摩擦付き管路の流れ — アニメーション

亜音速ガスが管に入り、壁面摩擦で加速しながらマッハ数が上昇します。チョーク時は管の途中に M=1 の壁が現れます。

Fanno関数 fL*/D vs マッハ数

出口マッハ数 vs 配管長さ

理論・主要公式

$$\frac{fL^{*}}{D}=\frac{1-M^{2}}{\gamma M^{2}}+\frac{\gamma+1}{2\gamma}\ln\!\frac{(\gamma+1)M^{2}}{2+(\gamma-1)M^{2}}$$

ファノー関数 fL*/D は、マッハ数 M の地点から摩擦だけで M=1（チョーク）に到達するまでに必要な無次元摩擦長さです。亜音速側では摩擦が M を 1 に向かって増加させ、超音速側では 1 に向かって減少させます。

$$\frac{fL_{1\to 2}}{D}=\left(\frac{fL^{*}}{D}\right)_{\!M_{1}}-\left(\frac{fL^{*}}{D}\right)_{\!M_{2}}$$

入口から出口までの実際の摩擦長さ fL/D は、入口と出口の Fanno 関数値の差に等しい。出口値が 0 以下になる管長では摩擦チョークが発生します。

$$\frac{P}{P^{*}}=\frac{1}{M}\sqrt{\frac{\gamma+1}{2+(\gamma-1)M^{2}}}$$

Fanno 線上の静圧比 P/P*。出口/入口の圧力比は P₂/P₁ = (P/P*)|_{M₂} / (P/P*)|_{M₁} で求められます。亜音速流れでは出口に向かって静圧が下がります。

ファノー流れ（摩擦付管路）とは

🙋

「ファノー流れ」って、長い配管を流れるガスの計算ですよね？普通のベルヌーイ式じゃダメなんですか？

🎓

いい質問だね。ベルヌーイ式は「非圧縮性・摩擦無視」が前提だから、短い水道管とかなら十分なんだ。でも、空気や天然ガスを数十メートル以上の細い管で流すと、密度が場所ごとに変わる「圧縮性」と、壁面摩擦による損失の両方が効いてくる。ファノー流れは「断熱（熱の出入りなし）＋等断面＋壁面摩擦あり」という3つの条件で、その2つを真面目に解いたモデルなんだよ。

🙋

面白いのは、左のマッハ数を 0.3 とか 0.5 にしてスライダーで管を長くしていくと、出口マッハ数がどんどん上がっていくところです。摩擦なのに「速くなる」って直感に反します。

🎓

そこがファノー流れの最大の見どころだね。非圧縮の感覚だと「摩擦＝抵抗＝遅くなる」だけど、ガスは違う。摩擦でエネルギーが熱に変わって静温度 T が下がる → 状態方程式から密度 ρ も下がる → 等断面では ρuA=一定だから速度 u が上がる、というドミノが起こるんだ。さらに M=u/√(γRT) は分子の u が上がって分母の √T が下がるから、マッハ数は二重に上昇する。これが亜音速流れで M が 1 に向かう理由だよ。

🙋

配管をさらに長くすると、画面で「チョーク」って赤くなりますよね。あれって何が起きているんですか？

🎓

「摩擦チョーク」と呼ばれる現象で、ファノー流れの一番有名な結果なんだ。入口の Fanno 関数値 fL*/D は「ここから M=1 に達するまでに必要な無次元摩擦長さ」を表していて、実際の管路の摩擦長さ fL/D がこの値に達した瞬間に管内のどこかで M=1 になる。それ以上長い管はそもそも物理的に成立しなくて、上流側の流量が自動的に絞られるんだ。だから長距離パイプラインを設計するときは「f·L/D が入口の fL*/D を超えていないか」が最重要チェックになる。

🙋

右の Fanno 関数のグラフ、M=1 付近で 0 に近づくのが見えます。これは「M が 1 に近いほどチョークが近い」ってことですか？

🎓

そう、まさにその通り。曲線は M=0 で発散して、M=1 でゼロになる単調減少関数なんだ。だから入口マッハ数が大きいほど fL*/D が小さい＝チョークまでの余裕が少ない。例えば M₁=0.3 なら fL*/D≈5.3 だけど、M₁=0.5 だと約 1.07、M₁=0.7 では約 0.21 まで急減する。だから「入口の流速を絞って流量を確保したい」ときは、配管径を太くするか、ブースターポンプで圧力を上げて入口マッハ数を下げる工夫をするんだよ。

🙋

実務では、出口の圧力比 P₂/P₁ も大事ですよね？このツールでは 0.84 とか出てます。

🎓

そう、配管設計の目的は「出口でいくらの圧力が残っているか」だから、P₂/P₁ は最終的な判断材料になる。ファノー流れでは亜音速の場合、出口で静圧 P は必ず下がる（圧力損失が出る）。今の条件だと約 16% の圧力降下だね。これに対し、レイリー流れ（加熱）では加熱量しだいで圧力が上がることも下がることもある。長距離天然ガスパイプラインは典型的なファノー問題で、コンプレッサーで定期的に再昇圧して P を回復させるんだ。

よくある質問

非圧縮性の感覚では摩擦は流れを遅くしますが、断熱・等断面の圧縮性管路では話が逆になります。摩擦でエネルギーが散逸すると静温度・静圧・密度が下がり、等断面では質量保存（ρuA=const）から速度 u が上がります。マッハ数 M は M=u/√(γRT) で速度上昇と温度低下の両方で増加するため、亜音速流れでは M が必ず M=1 に向かって増加します。逆に超音速流れでは摩擦が M を 1 に向かって減少させます。どちらも「終点は M=1」という点が摩擦チョークの本質です。

摩擦チョークとは、面積変化も熱の出入りもないのに、壁面摩擦だけで管路の途中で M=1 に達して詰まる現象です。判定指標は無次元摩擦長さ fL/D で、入口マッハ数 M1 における Fanno 関数 fL*/D（M1 から M=1 までに必要な摩擦長さ）と比較します。fL/D が fL*/D に達した時点で M=1 となり、そこから先の管路は物理的に存在できません。実務では fL/D ≤ 0.7·(fL*/D) を目安に余裕を持たせ、温度や流量の変動にも対応できる設計にします。

両者とも等断面の1次元圧縮性流れですが、駆動メカニズムが異なります。ファノー流れは「断熱＋摩擦」で、エントロピーが増えながら M=1 に近づきます。レイリー流れは「摩擦無視＋熱の出入り」で、亜音速では加熱で M↑、冷却で M↓となります。両方とも亜音速流れでは M=1 に向かって加速するように見えますが、ファノーは温度が下がりながら M↑、レイリー（加熱）は温度が上がりながら M↑という違いがあります。長距離パイプラインのように摩擦が支配的な場合はファノー、燃焼器のように熱付加が支配的な場合はレイリーで解析します。

f はムーディ線図またはコールブルック式から決めます。商用鋼管なら相対粗さ ε/D ≈ 5×10⁻⁵～5×10⁻⁴、レイノルズ数が 10⁵～10⁷ の完全乱流域で f ≈ 0.015～0.030 が目安です。長距離天然ガスパイプライン（D=500mm 程度）では f≈0.012、空調ダクト（D=300mm 程度）の鋼板で f≈0.020、ゴム/プラスチックホースは f≈0.030 を初期値に使うことが多いです。なお Fanning 摩擦係数 f_F とは f = 4·f_F の関係にあるので、文献の係数定義に注意してください。本ツールはダルシー摩擦係数（Darcy-Weisbach）を採用しています。

実世界での応用

長距離天然ガスパイプライン：シベリアや北米大陸を横断する数千km級のパイプラインでは、ファノー流れによる圧力損失と摩擦チョークの余裕設計が中核となります。数百kmごとにコンプレッサーステーションを設けて圧力を再回復させるのは、まさに「出口の Fanno 関数値が 0 に近づく前に上流に戻す」運用です。直径 600〜1200mm のパイプを摩擦係数 f≈0.012 で運用し、入口マッハ数を 0.05〜0.10 と低めに抑えることで、ステーション間の距離を稼いでいます。

建物の換気・空調ダクト設計：大規模ビルや工場の HVAC（空調換気）ダクトは、長さが 100m を超え、分岐や曲がりも多いため、摩擦による圧力損失が支配的になります。亜音速（M≈0.1〜0.2）での運用ですが、ダクト径と長さによっては摩擦チョークに至るケースもあり、ファノー解析でファンの必要静圧と最終吹出し風量を見積もります。大型データセンターの冷却ダクトでも、サーバーラックへの均等な風量分配にこの考え方が使われます。

ガスタービン翼の内部冷却通路：タービン入口温度 1500℃以上の現代の航空エンジンでは、翼の内部に細い冷却通路を貫通させて圧縮空気を流し、対流冷却します。直径 1〜3mm の細い通路を高速のガスが流れるため摩擦が大きく、ファノー流れに近い挙動を示します。設計者は通路出口でチョークしないように管径と長さの組み合わせを決め、必要な冷却空気量を確保しています。チョークすると翼の冷却が不足してクリープ破壊につながるため、極めて重要な解析です。

空気圧搬送（ニューマチック・コンベヤ）：セメント工場や食品工場で粉体を空気流で運ぶシステムは、長さ数百mの配管を高速ガスで流す典型例です。粉体添加で実効摩擦係数が増加するため、純空気よりさらに厳しいファノー条件下にあり、入口圧力を高めにとっておかないと粉体が途中で詰まる「ブロッキング」が発生します。設計時にファノー流れで「どの長さで M=1 に達するか」を逆算し、ブースターを配置する位置を決めます。

よくある誤解と注意点

まず最も多い誤解が、「摩擦は常に流れを遅くする」という非圧縮の直感をそのまま圧縮性に持ち込むこと。亜音速ファノー流れでは摩擦が速度を上げ、マッハ数を 1 に向かって増加させます。直感的には「ガスが温度低下でぎゅっと縮むから、同じ質量流量を保つために速くならざるを得ない」と理解すると間違えません。逆に超音速ファノー流れでは摩擦が M を 1 に向かって減少させますが、これも「M=1 が摩擦の終着点」という共通原理から導かれます。「亜音速も超音速も摩擦で M=1 に近づく」が要点で、「速くなる」「遅くなる」のどちらでも答えてしまうと半分しか合っていません。

次に、摩擦係数 f の定義の取り違えです。ファノー流れの文献には「Darcy 摩擦係数 f」と「Fanning 摩擦係数 f_F」の2種類が混在しており、両者は f = 4·f_F の関係にあります。式 fL*/D の f はダルシー（本ツールもダルシー基準）ですが、Fanning ベースで書かれた教科書では 4fL*/D と書かれることもあります。係数の取り違えは「実際の管路長より 4 倍短い場所でチョークすると誤判定する」という致命的なミスにつながります。ムーディ線図はダルシー、海外の化学工学ハンドブックは Fanning が多いので、出典を必ず確認してください。

最後に、「ファノー流れは厳密に断熱を仮定している」という前提を忘れること。実際の配管は外部と熱のやり取りがあり、特に長距離パイプラインでは地中温度に向かって冷却（または加熱）されます。本ツールの結果は「断熱・等断面・1次元・ガスは理想気体」という理想化下のものです。実務で 10% 以上の精度が必要なら、ファノー（摩擦）とレイリー（熱）の両方を組み合わせた一般化された圧縮性管路流れの数値解（区分積分や CFD）が必要です。とはいえ、ファノー解析は「摩擦支配かどうか」「チョークの余裕がどれくらいか」を一発で見るスクリーニングツールとして極めて有効で、本格的な CFD 解析の前段として今でも欠かせない手法です。