摩擦係数f_Dってどうやって決めるんですか？

流れの状態で使う式が違う。層流（Re 4000）では管壁粗さe/Dの影響が出て、Moody線図（実験データ）から読み取るか、Colebrook-White方程式（陰関数）か、その近似のSwamee-Jain式（陽関数）を使う。実務ではMoodychartのExcelやPython実装（scipy）で即計算できる。Re = 10^5〜10^6、e/D = 0.0001の場合 f_D ≈ 0.018程度という感覚値は持っておくと設計概算に便利だよ。

曲管やバルブの圧力損失はどう計算するんですか？

局所損失として扱う。delta_p_local = K * (rho*U^2/2) という形で、損失係数Kを使う。エルボ（90°曲管）でK = 0.3〜0.9（r/D比で変わる）、グローブバルブで全開K = 7〜10、チェックバルブで2〜3——これらはHansen, Idelchick, Craneのハンドブックに載っている。配管系全体の設計では「直管損失（Darcy-Weisbach）＋局所損失（K値）」の和でポンプ揚程を決める。CFDはこのK値が不確かな形状——急拡大・合流部・複雑な弁形状——の精密評価に使う。

CFDでDarcy-Weisbach式との検証をする手順を教えてください。

標準的な検証フローは——①十分長い直管モデル（入口から完全発達するまでのL/D > 60）をメッシュ切り、②CFDで圧力勾配dP/dLをx方向に出力、③その値からf_D = 2*D*(dP/dL)/(rho*U^2)を計算、④Colebrook式やMoody線図の値と比較——だ。乱流ではy+ の設定が重要で、ウォール関数使用なら30

Darcy-Weisbach式 — CAE用語解説

Q: 先生、Darcy-Weisbach式って管路の圧力損失を計算する式ですよね。CFDでも使いますか？

CFDでは直接使うというより検証（V&V）に使うことが多い。Darcy-Weisbach式は delta_p = f_D * (L/D) * (rho*U^2/2) で、摩擦係数f_D（Darcy摩擦係数）、管長さL/径Dの比、動圧rho*U^2/2の積で圧力損失を求める。この解析解と直管CFDの結果を比較してソルバーが正しく摩擦損失を計算できているか確認する——これが配管系CFD解析の基本的な検証ステップだ。実用設計ではCFDより計算が速いから、まず相関式で概算してから形状が複雑な部分だけCFDを使う使い分けが多い。

カテゴリ: 用語集 | 2026-01-15

CAE visualization for darcy weisbach - technical simulation diagram

Darcy-Weisbach式

🧑‍🎓

先生、Darcy-Weisbach式って管路の圧力損失を計算する式ですよね。CFDでも使いますか？

理論と物理

Darcy-Weisbach式の基本概念

🧑‍🎓

Darcy-Weisbach式って、教科書では圧力損失を計算する式と書いてありますが、具体的に何の圧力損失なんですか？

🎓

管路やダクト内を流体が流れる時に生じる「摩擦による沿程損失」を計算するための式です。例えば、内径50mm、長さ100mの鋼管を水が流れる時、この式で管壁との摩擦で失われる圧力を求めます。式の基本形は

$$ h_f = f \frac{L}{D} \frac{V^2}{2g} $$

です。ここで

$$ h_f $$

は損失ヘッド[m]、

$$ f $$

は摩擦係数、

$$ L $$

は管長[m]、

$$ D $$

は管径[m]、

$$ V $$

は平均流速[m/s]、

$$ g $$

は重力加速度です。

🧑‍🎓

損失ヘッドを圧力損失[Pa]に直すにはどうすればいいですか？

🎓

流体の密度

$$ \rho $$

を使って

$$ \Delta P = \rho g h_f $$

で変換します。水（

$$ \rho \approx 1000 \, \text{kg/m}^3 $$

）の場合、損失ヘッド1mは約9.8kPaの圧力損失に相当します。実務ではこの圧力損失をポンプの選定や配管システムの圧力評価に使います。

🧑‍🎓

式の中の摩擦係数fはどうやって決めるんですか？教科書には「ムーディ線図」と書いてありますが。

🎓

その通り、ムーディ線図が基本です。これは相対粗度（管の内面の凹凸の高さ/管径）とレイノルズ数

$$ Re $$

の関数です。実用的には、Colebrook-Whiteの式

$$ \frac{1}{\sqrt{f}} = -2.0 \log_{10} \left( \frac{\epsilon/D}{3.7} + \frac{2.51}{Re\sqrt{f}} \right) $$

を暗黙的に解くか、近似式を使います。例えば、新しい塩ビ管の絶対粗さεは約0.0015mm、鋳鉄管だと0.26mm程度です。

🧑‍🎓

レイノルズ数が非常に小さい層流の時は、この式は成り立たないんですか？

🎓

成り立ちますが、摩擦係数fが単純な理論式で与えられます。層流（

$$ Re < 2300 $$

）では、

$$ f = \frac{64}{Re} $$

です。これをDarcy-Weisbach式に代入すると、ハーゲン・ポアズイユの式と一致し、圧力損失が流速に比例することがわかります。乱流では流速の2乗に比例するので、流れの状態で損失特性が大きく変わります。

数値解法と実装

CAEソフトウェアでの取り扱い

🧑‍🎓

CFDで配管の流れを全部メッシュ切って解くのと、Darcy-Weisbach式を使うのと、どう使い分けるべきですか？

🎓

システム全体の圧力損失やポンプ動力を見積もる「1次元ネットワーク解析」ではDarcy-Weisbach式が使われます。Ansys Fluentの「Porous Jump」条件や、専用のシステム解析ソフト（例えばFlowmaster）の根幹です。一方、弁の近傍の複雑な渦や局所的な剥離を詳細に知りたい時は3D CFDが必要です。計算コストは桁違いで、1Dネットワーク解析は数秒、3D CFDは数時間から数日かかります。

🧑‍🎓

CFDソフトで「等価粗さ」を設定する項目がありますが、あれはこの式のεと関係あるんですか？

🎓

直接的に関係します。Ansys FluentやSTAR-CCM+の壁面条件で設定する「Sand-grain Roughness」が、まさにColebrook-White式で使う絶対粗さεに相当します。ソルバーはこの値を用いて壁面近傍の乱流モデル（例えば、k-εモデルの壁面関数）を修正し、粗い壁面での摩擦抵抗増加を再現します。設定値は材質ごとのカタログ値、例えば機械加工面なら3.2μm（Ra値）などを参考にします。

🧑‍🎓

ネットワーク解析ソフトでは、式の中のfはどう計算しているんですか？毎回Colebrook式を解いているんですか？

🎓

実用的には、高速で精度の良い近似式が使われます。例えば、Swamee-Jainの式

$$ f = \frac{0.25}{ \left[ \log_{10} \left( \frac{\epsilon/D}{3.7} + \frac{5.74}{Re^{0.9}} \right) \right]^2 } $$

です。これはColebrook式の陽的近似で、レイノルズ数が

$$ 5 \times 10^3 $$

から

$$ 10^8 $$

、相対粗度が

$$ 10^{-6} $$

から

$$ 10^{-2} $$

の範囲で誤差2%以内と言われています。ソフト内部ではこのような式で瞬時にfを算出しています。

実践ガイド

設計・解析ワークフロー

🧑‍🎓

実際に配管システムの圧力損失を手計算で見積もる時、最初に何をすべきですか？

🎓

まずシステムの単線図（P&ID）を書き、全ての直管区間の長さLと内径D、流量Qをリストアップします。次に、使用する配管材質から「絶対粗さε」の値を決めます。JISやASMEの規格、メーカーカタログに参考値があります。例えば、JIS G 3452の炭素鋼管（SGP）の内面はε=0.046mm程度です。これで各区間の相対粗度ε/Dが決まります。

🧑‍🎓

流量Qから流速Vを求める時、内径は公称径でいいんですか？

🎓

ダメです。公称径（A呼称）は外径を基準にした呼び方です。流速計算には実際の内径を使わなければなりません。例えば、公称径100AのSGP管（スケジュール40）の実際の内径は約102.3mmです。カタログやJIS B 2201（配管用フランジの呼び圧力）などの規格表で確認する必要があります。ここを間違えると流速とレイノルズ数が大きくずれ、損失計算が台無しになります。

🧑‍🎓

エルボや弁などの部品による損失は、この式では計算できないですよね？

🎓

その通りで、それらは「形状抵抗」や「局部抵抗」と呼ばれ、別途「損失係数K」を用いた式

$$ h_f = K \frac{V^2}{2g} $$

で計算します。実務では、直管のDarcy損失と、全ての継手・弁のK値による損失を合算します。K値は「CraneのTechnical Paper No.410」や各メーカーのカタログに掲載されています。例えば、標準エルボ（90°）のK値は約0.3です。

🧑‍🎓

計算結果をどう検証すれば信頼性が高まりますか？

🎓

二つの方法があります。第一に、既存の類似設備の実測データ（ポンプ吐出し圧力、流量計指示値）と比較する。第二に、異なる情報源から粗さεや損失係数Kを選んで計算し、結果のバラツキ（感度分析）を見る。例えば、配管の粗さを「新品」と「経年使用」の両方の値で計算し、圧力損失がどの範囲に収まるか評価します。これが設計上の安全マージンになります。

ソフトウェア比較

各ツールでの実装と特徴

🧑‍🎓

Ansys Fluentで配管流れを解析する時、壁面関数とDarcy-Weisbach式はどう関係しているんですか？

🎓

Fluentが内部で行っているのは、NS方程式と乱流モデルを解くことで、Darcy-Weisbach式を「導出する」ようなものです。壁面関数（Standard Wall FunctionsやScalable Wall Functions）は、粗さの影響を考慮して壁面せん断応力を計算するモデルです。ユーザーが壁面に粗さを設定すると、計算結果として得られる壁面せん断応力から逆算される摩擦係数fは、実測やColebrook式に近い値になります。つまり、CFDはより物理的なアプローチで「f」を決定していると言えます。

🧑‍🎓

1Dシステム解析ソフト、例えばAmesimやFlowmasterでは、この式はどのように使われていますか？

🎓

これらのソフトの配管（パイプ）コンポーネントの根幹がDarcy-Weisbach式です。ユーザーは管長、内径、粗さ、流体特性を入力するだけで、コンポーネント内部で自動的に摩擦係数fを計算し、流量と圧力損失の関係をモデル化します。Flowmasterでは「Colebrook-White」「Hazen-Williams」など複数の損失式から選択可能です。油圧システムでは、作動油の粘度変化によるレイノルズ数変化も自動追従するので、温度変化を含む過渡解析が可能です。

🧑‍🎓

建築設備の分野でよく聞く「イギリス式」「アメリカ式」ってDarcy-Weisbach式と違うんですか？

🎓

建築の空調・給排水では、Hazen-Williamsの式がよく使われます。これは

$$ h_f = 10.67 \cdot C^{-1.852} \cdot D^{-4.871} \cdot L \cdot Q^{1.852} $$

という経験式で、摩擦係数fの代わりに「C値」（粗度係数）を使います。Darcy式が物理量に基づくのに対し、Hazen-Williamsは水に特化した簡便式です。C値は材質と状態で決まり、新しい鋼管でC=130、新しい鋳鉄管でC=130、古い管でC=80-100など。ソフトではどちらの式も選べるものが多いです。

🧑‍🎓

COMSOL Multiphysicsで配管の圧力損失を簡単にモデル化する方法はありますか？

🎓

「Pipe Flow Module」を使う方法と、「分散損失」を設定する方法があります。Pipe Flow Moduleは1Dネットワーク解析を可能にし、内部でDarcy-Weisbach式などを用います。もう一つの方法は、3D流路モデルに「集中損失」特徴を追加し、圧力ジャンプ条件として

$$ \Delta P = K \cdot \frac{1}{2} \rho V^2 $$

を設定するものです。このK値に、直管部分の損失相当分

$$ fL/D $$

を代入すれば、Darcy損失を簡易的に表現できます。詳細な3Dモデルを作りたくない長い直管部分の省略に有効です。

トラブルシューティング

よくあるエラーと対策

🧑‍🎓

手計算で求めた圧力損失と、CFDで計算した結果が大きく違いました。考えられる原因は？

🎓

まず疑うべきは「粗さの設定不一致」です。手計算で使ったεと、CFDの壁面条件で入力したSand-grain roughnessが同じ値ですか？次に「完全発達流れかどうか」です。Darcy式は入口から十分下流の発達流れを仮定していますが、CFDモデルが短すぎると入口影響が残り損失が過小評価されます。目安として、入口から

$$ 10D $$

以上離れた区間の損失を比較しましょう。また、CFDの乱流モデル（特に壁面関数）の選択も影響します。

🧑‍🎓

ネットワーク解析ソフトで、流量を増やすと圧力損失が急激に増加する非線形な結果が出ます。これは正しい動作ですか？

🎓

はい、正しい動作です。乱流域では摩擦係数fがほぼ一定とみなせるため、Darcy-Weisbach式は

$$ \Delta P \propto V^2 \propto Q^2 $$

となり、圧力損失は流量の2乗に比例します。流量が2倍になれば損失は4倍です。もし比例関係（線形）に見えるなら、層流計算になっているか、ソフトの設定で「線形損失」のような簡易モードが選択されている可能性があります。ポンプカーブとの連成解析では、この2乗特性がシステムの動作点を決定する重要な要素です。

🧑‍🎓

Colebrook式をExcelで解こうとすると、循環参照になってしまいます。どうすればいいですか？

🎓

その通り、fが両辺にある陰方程式なので直接代入では解けません。Excelなら「ゴールシーク」機能を使うか、前述のSwamee-Jainの陽的近似式を使うのが現実的です。VBAでニュートン・ラフソン法を実装する方法もあります。初期値として、乱流の粗い管ではf=0.02、滑らかな管ではf=0.012などを与えて反復計算させます。実務では、いちいち解かずに既存の近似式か、ソフトウェアの計算結果を信用するのが早いです。

🧑‍🎓

配管が水平じゃない場合、Darcy-Weisbach式はそのまま使えますか？重力の影響は？

🎓

式そのものは摩擦損失のみを表すので、傾きに関係なく使えます。ただし、システム全体の圧力バランスを考える時は、位置エネルギー（重力による圧力変化）を別途考慮する必要があります。ベルヌーイの式の形で考えると、二点間の全圧損失

$$ \Delta P_{total} $$

は、摩擦損失

$$ \Delta P_{friction} $$

と重力項

$$ \rho g \Delta z $$

の和（または差）になります。上り勾配では重力が抵抗として加わり、下り勾配では重力が推進力として働くので、計算では符号に注意が必要です。