ダクト内流れ
理論と物理
概要
先生！ ダクト内流れの解析って、空調配管やプラント配管で使うやつですよね？ 基礎から教えてください。
ダクト内流れのCFD解析は、配管やダクト系統の圧力損失予測、流量分配、偏流評価を目的とする。設計段階でDarcy-Weisbach式の手計算だけでは捉えきれない局所損失や二次流れをCFDで可視化するんだ。
支配方程式
圧力損失の基本式はDarcy-Weisbachですよね。
そう。直管部の摩擦損失はDarcy-Weisbach式で記述される。
ここで $f$ は管摩擦係数、$L$ は管長、$D_h$ は水力直径、$V$ は断面平均流速だ。層流の場合は $f = 64/Re$、乱流の場合はColebrookの式で求める。
Colebrookは陰的な式だから反復計算が必要ですね。実務ではSwamee-Jainの近似式を使うことも多いですか？
その通り。Swamee-Jainは陽的で実用上十分な精度がある。
局所損失（エルボ、分岐、拡大・縮小）は損失係数 $K$ で表す。
| 要素 | 損失係数 K（目安） |
|---|---|
| 90° エルボ（R/D=1.5） | 0.2〜0.3 |
| 90° マイター（ベーンなし） | 1.1〜1.3 |
| T字分岐（直進） | 0.3〜0.5 |
| T字分岐（分流） | 0.8〜1.3 |
| 急拡大 | $(1 - A_1/A_2)^2$ |
| 急縮小 | $0.5(1 - A_2/A_1)$ |
手計算での損失係数は文献値ですが、CFDだとジオメトリ固有の正確な値が出せるわけですね。
そう。特に角型ダクトのコーナーピースや複雑な分岐管は文献値がない場合が多いから、CFDで求める価値がある。
乱流モデルの選択
ダクト内流れに適した乱流モデルは何ですか？
| 乱流モデル | 推奨用途 | 備考 |
|---|---|---|
| Realizable k-epsilon | 直管・エルボ | 汎用、壁関数で高速計算 |
| SST k-omega | 剥離・急拡大 | 逆圧力勾配に強い |
| RSM (Reynolds Stress) | 旋回流・二次流れ | 精度高いが計算コスト大 |
角型ダクトでは二次流れ（コーナー渦）が発生しますが、k-epsilonで捉えられますか？
角型ダクトの二次流れはReynolds応力の非等方性に起因するから、厳密にはRSMが必要だ。ただし圧損予測が目的なら k-epsilon でも誤差は5%程度に収まる。
「助走区間」の理論——ダクト入口から何D進めば乱流は完成するか
ダクト内流れの理論で最初に習う重要概念が「水力学的助走区間」です。入口から流れが壁面の境界層に影響されながら、断面全体に渡る充分発達した乱流プロファイルになるまでの距離を指します。乱流の場合、おおよそ $x \approx 10 \sim 60 D$ （直径Dに対して）が必要とされています。実務のCFD解析でよくあるミスが「入口条件を一様流にしたまま解析領域をギリギリの長さにする」こと。助走区間を十分確保するか、実測の速度プロファイルを入口条件に設定しないと、ダウンストリームの圧力損失が過小評価されがちです。
各項の物理的意味
- 時間項 $\partial(\rho\phi)/\partial t$：蛇口をひねった瞬間を思い浮かべてください。最初は水がバタバタと不安定に出て、しばらくすると安定した流れになりますよね？ この「変化している最中」を記述するのが時間項です。心臓の拍動で血流が脈打つのも、エンジンのバルブが開閉するたびに流れが変動するのも、すべて非定常現象。では定常解析とは？ 「十分時間が経って流れが落ち着いた後」だけを見る——つまりこの項をゼロにする。計算コストが大幅に下がるため、まず定常で解いてみるのがCFDの基本戦略です。
- 対流項 $\nabla \cdot (\rho \mathbf{u} \phi)$：川に落ち葉を落としたらどうなりますか？ 流れに乗って下流に運ばれますよね。これが「対流」——流体の動きが物を運ぶ効果です。暖房の温風が部屋の端まで届くのも、空気という「運び屋」が熱を対流で輸送しているから。ここが面白いところ——この項は「速度×速度」を含むため非線形です。つまり、流れが速くなるとこの項が急激に強くなり、制御が難しくなる。これが乱流の根本原因です。よくある勘違い：「対流と伝導は同じようなもの」→ 全然違います！ 対流は流れが運ぶ、伝導は分子が伝える。桁違いの効率差があります。
- 拡散項 $\nabla \cdot (\Gamma \nabla \phi)$：コーヒーにミルクを入れて放置したことはありますか？ かき混ぜなくても、しばらく経つと自然に混ざりますよね。あれが分子拡散です。では次の質問——ハチミツとお水、どちらが流しやすいですか？ 当然お水ですよね。ハチミツは粘性（$\mu$）が高いから流れにくい。粘性が大きいと拡散項が強くなり、流体は「もったりした」動きになります。レイノルズ数が小さい流れ（ゆっくり、ドロドロ）では拡散が支配的。逆にRe数が大きい流れでは対流が圧倒し、拡散は脇役になります。
- 圧力項 $-\nabla p$：注射器のピストンを押すと、液体が針先から勢いよく出ますよね？ なぜでしょう？ ピストン側が高圧、針先が低圧——この圧力差が流体を押す力になるからです。ダムの放水も同じ原理。天気図で等圧線がギュッと密になっている場所では？ そう、強風が吹きます。「圧力差があるところに流れが生まれる」——これがナビエ-ストークス方程式の圧力項の物理的意味。ここでの勘違いポイント：CFDの「圧力」は絶対圧ではなくゲージ圧のことが多い。圧縮性解析に切り替えたとたんに結果がおかしくなる場合、絶対圧/ゲージ圧の混同が原因かもしれません。
- ソース項 $S_\phi$：暖められた空気が上に昇る——なぜでしょう？ 周囲より軽く（密度が低く）なったから、浮力で押し上げられるのです。この浮力はソース項として方程式に追加されます。他にも、ガスコンロの炎で化学反応熱が発生する、工場の電磁ポンプで金属溶湯にローレンツ力がかかる…これらはすべて「外部から流体にエネルギーや力を注入する」作用であり、ソース項で表現します。ソース項を忘れるとどうなるか？ 自然対流の解析で浮力を入れ忘れると、流体は一切動かない——冬の部屋で暖房をつけたのに暖かい空気が上に行かない、という物理的にありえない結果になります。
仮定条件と適用限界
- 連続体仮定：クヌッセン数 Kn < 0.01（分子平均自由行程 ≪ 代表長さ）で成立
- ニュートン流体仮定：せん断応力と歪み速度が線形関係（非ニュートン流体では粘度モデルが必要）
- 非圧縮性仮定（Ma < 0.3の場合）：密度を一定として扱う。マッハ数0.3以上では圧縮性効果を考慮
- ブシネスク近似（自然対流）：密度変化を浮力項のみで考慮し、他の項では一定密度を使用
- 適用外ケース：希薄気体（Kn > 0.1）、超音速・極超音速流れ（衝撃波捕捉が必要）、自由表面流れ（VOF/Level Set等が必要）
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 速度 $u$ | m/s | 入口条件で体積流量から換算する際、断面積の単位に注意 |
| 圧力 $p$ | Pa | ゲージ圧と絶対圧の区別。圧縮性解析では絶対圧を使用 |
| 密度 $\rho$ | kg/m³ | 空気: 約1.225 kg/m³@20°C、水: 約998 kg/m³@20°C |
| 粘性係数 $\mu$ | Pa·s | 動粘性係数 $\nu = \mu/\rho$ [m²/s] との混同に注意 |
| レイノルズ数 $Re$ | 無次元 | $Re = \rho u L / \mu$。層流/乱流遷移の判定指標 |
| CFL数 | 無次元 | $CFL = u \Delta t / \Delta x$。時間刻みの安定性に直結 |
数値解法と実装
数値手法の詳細
ダクト内流れをCFDで解くとき、メッシュや境界条件で気をつけることは何ですか？
まずメッシュから説明しよう。
メッシュ戦略
円管と角型ダクトでメッシュの作り方は変わりますか？
大きく変わる。円管にはO-gridトポロジー（蝶ネクタイ型）が推奨で、壁面に直交するプリズム層を確保しやすい。角型ダクトはスイープメッシュでプリズム層を入れる。
壁面第一層の高さは、使用する壁モデルに合わせる。
| 壁モデル | 必要な y+ | 第一層高さの目安（Re=10⁵, D=300mm） |
|---|---|---|
| Enhanced Wall Treatment | ≒ 1 | 約0.05 mm |
| Standard Wall Function | 30〜300 | 1〜10 mm |
| Scalable Wall Function | > 11.225 | > 0.4 mm |
y+ = 1にするとセル数がかなり増えますね。圧損精度の観点で壁関数でも十分ですか？
直管の摩擦損失だけなら壁関数で十分だ。ただし剥離を伴う急拡大やバルブ後方では、壁面分解（y+ ≒ 1）の方が精度が上がる。
境界条件の設定
入口・出口の境界条件はどう設定しますか？
典型的な設定パターンを示す。
| 境界 | 条件タイプ | 設定値 |
|---|---|---|
| ダクト入口 | Velocity Inlet | 設計風速 + 乱流強度5%、水力直径 |
| ダクト出口 | Pressure Outlet | ゲージ圧0 Pa |
| ファン位置 | Fan BC (Pressure Jump) | ファン特性曲線 |
| ダンパー | Porous Jump | 開度に応じた抵抗係数 |
| 壁面 | No-Slip Wall | 粗さ高さ（鋼管: 0.045 mm） |
壁面粗さをCFDに入れるんですね。材質ごとの粗さ高さはどこで調べますか？
ASHRAE Handbook FundamentalsやCrane TP-410に代表的な値が載っている。
| 材質 | 等価粗さ [mm] |
|---|---|
| 亜鉛鉄板ダクト | 0.09〜0.15 |
| 鋼管 | 0.045 |
| 塩ビ管 | 0.0015 |
| コンクリートダクト | 0.3〜3.0 |
| フレキシブルダクト | 1.0〜4.6 |
入口助走区間の処理
完全発達流を仮定する場合、助走区間はどう扱いますか？
乱流の助走区間は概ね $L_{entry} \approx 10 D_h$ だ。入口直後の圧損評価が目的でなければ、十分な助走区間を設けるか、Fully Developed Profileを入口条件に与える。Fluentでは入口にMapped条件（出口の速度プロファイルを入口にマッピングする周期条件）を使う方法もある。
ソルバー設定
具体的なソルバー設定の推奨値を教えてください。
| パラメータ | 推奨設定 |
|---|---|
| ソルバー | Pressure-Based, Steady |
| 圧力-速度連成 | SIMPLEC |
| 対流スキーム | Second Order Upwind |
| 圧力補間 | Second Order |
| 勾配 | Least Squares Cell-Based |
| 収束判定 | 残差 1e-4以下 + 出入口流量バランス < 0.1% |
出入口の流量バランスを収束判定に使うのは実用的ですね。残差だけだと見落とすことがありますか？
ある。残差が1e-4に下がっていても、出入口の質量流量差が1%以上あるケースがある。必ず物理量モニター（入口圧力、出口流量）の定常化も確認すること。
ダルシー・ワイスバッハ式——160年前の公式が今も現役な理由
ダクト内流れの圧力損失計算で今も使われるダルシー・ワイスバッハ式 $\Delta p = f \cdot (L/D) \cdot (\rho u^2/2)$ は、1850年代にヘンリー・ダルシーとジュリウス・ワイスバッハが独立に提唱したものです。160年以上経った現代でも、CFD解析の検証用ベンチマークや配管システムの概算設計で第一線で使われています。面白いのは「摩擦係数 f 」の部分で、層流ならば $f = 64/Re$ と解析的に決まりますが、乱流だとムーディー線図や各種実験式に頼ることになる。CFDはこの実験式が使えない複雑形状・多分岐系でこそ真価を発揮します。
風上差分（Upwind）
1次風上: 数値拡散が大きいが安定。2次風上: 精度向上するが振動のリスク。高レイノルズ数流れでは必須。
中心差分（Central Differencing）
2次精度だが、Pe数 > 2で数値振動が発生。低レイノルズ数の拡散支配流れに適する。
TVDスキーム（MUSCL、QUICK等）
リミッタ関数により数値振動を抑制しつつ高精度を維持。衝撃波や急勾配の捕捉に有効。
有限体積法 vs 有限要素法
FVM: 保存則を自然に満足。CFDの主流。FEM: 複雑形状・マルチフィジックスに有利。SPH等のメッシュフリー法も発展中。
CFL条件（クーラン数）
陽解法: CFL ≤ 1が安定条件。陰解法: CFL > 1でも安定だが、精度と反復回数に影響。LES: CFL ≈ 1を推奨。物理的意味: 1タイムステップで情報が1セル以上進まないこと。
残差モニタリング
連続の式・運動量・エネルギーの各残差が3〜4桁低下で収束と判断。質量保存の残差は特に重要。
緩和係数
圧力: 0.2〜0.3、速度: 0.5〜0.7が一般的な初期値。発散する場合は緩和係数を下げる。収束後は上げて加速。
非定常計算の内部反復
各タイムステップ内で定常解に収束するまで反復。内部反復数: 5〜20回が目安。残差がタイムステップ間で変動する場合は時間刻みを見直す。
SIMPLE法のたとえ
SIMPLE法は「交互に調整する」手法。まず速度を仮に求め（予測ステップ）、その速度で質量保存が満たされるよう圧力を補正し（補正ステップ）、補正された圧力で速度を修正する——このキャッチボールを繰り返して正解に近づく。2人で棚を水平にする作業に似ている：片方が高さを合わせ、もう片方がバランスを取り、これを交互に繰り返す。
風上差分のたとえ
風上差分は「川の流れに立って上流の情報を重視する」手法。川の中にいる人が下流を見ても水の出所は分からない——上流の情報が下流を決めるという物理を反映した離散化手法。精度は1次だが、流れの方向を正しく捕捉するため安定性が高い。
実践ガイド
実践ガイド
先生、実際のダクト系統CFDではどんなケースが多いですか？
典型的なケースを3つ紹介しよう。
ケース1: 分岐ダクトの流量分配
主ダクトから複数の分岐に均等に風量を分配したい場合ですね。
手計算ではT字分岐の損失係数を合算するが、連続する分岐では上流の偏流が下流に影響するため精度が落ちる。CFDなら全体を一体でモデル化して各分岐の流量を正確に予測できる。
設計検討の流れはこうだ。
1. 初期レイアウトでCFDを回して各分岐流量を確認
2. 流量が不均一な場合、ダンパー開度やガイドベーンの追加を検討
3. ダクトサイズの変更（断面縮小で動圧回復を利用）
4. 再計算で改善効果を検証
動圧回復ってStatic Regainのことですか？
そう。主ダクトの断面を下流に向かって段階的に縮小すると、流速増加で動圧が回復し、各分岐口での静圧が均一化される。Static Regain法は空調ダクト設計の定番手法だが、CFDで最適な縮小率を検証するのが効果的だ。
ケース2: エルボの圧損低減
エルボにガイドベーンを入れるとどのくらい圧損が減りますか？
90°マイターベンドの場合、ベーンなしでK≒1.2だが、シングルベーンでK≒0.5、ダブルベーンでK≒0.2まで低減できる。CFDではベーンの枚数・角度を最適化できる。
| ベーン構成 | 損失係数 K | 圧損低減率 |
|---|---|---|
| ベーンなし | 1.1〜1.3 | 基準 |
| シングルベーン | 0.4〜0.6 | 55% |
| ダブルベーン | 0.15〜0.25 | 82% |
| R/D=1.5 エルボ | 0.2〜0.3 | 78% |
ケース3: ファン出口の偏流評価
ファン直後のダクトで偏流が問題になるケースですね。
遠心ファンの出口は旋回成分を含む非一様な速度分布を持つ。これがダクトに接続されたとき、エルボまでの距離が短いと偏流が悪化してノイズや効率低下を引き起こす。ASHRAEではファン出口から2.5 Dh以上の直管区間を推奨している。
よくある失敗と対策
| 失敗パターン | 原因 | 対策 |
|---|---|---|
| 圧損が実測の2倍 | 壁面粗さの過大設定 | 材質に合った粗さ値を使用 |
| 分岐流量がアンバランス | 入口速度分布が一様すぎる | 上流のファンや曲がりも含めてモデル化 |
| エルボ後の剥離が消える | メッシュが粗い、y+が大きすぎる | エルボ外壁にプリズム層を追加 |
| ファンBCで逆流警告 | ファン前後の圧力差が特性曲線外 | ファンの運転点を確認 |
壁面粗さの影響って大きいんですね。設定を間違えるだけで圧損が倍になることもあると。
Re数が高い（10⁵以上）と粗さの影響が顕著になる。粗さ高さを0にした滑面と、コンクリートダクト（3mm）では摩擦係数が3倍以上違うことがある。
エルボー（曲がりダクト）の「ガイドベーン」は効果があるのか？
ダクトの90°エルボーで圧力損失を減らすために「ガイドベーン（整流翼）」を挿入するのはHVACの定番手法です。ガイドベーンがあると損失係数が0.9から0.1〜0.2程度まで激減することがあります。ところが実務でよくある落とし穴が「ガイドベーン間へのゴミ詰まり」。特に産業用ダクトでは粉塵や繊維が絡まりやすく、最初は低圧損でも時間経過で詰まり圧損が急増するケースがあります。CFDでは清浄状態しかモデル化できないため、メンテナンス性を含めた設計最適化は現場知識が不可欠です。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ CFDでも入口条件のミスは下流全体に波及します。
ソフトウェア比較
商用ツール比較
ダクト流れ解析に使えるツールを比較してもらえますか？
汎用CFDに加えて、ダクト系統専用の1D解析ツールも紹介しよう。
3D CFDツール
| ツール | 特徴 | ダクト解析での強み |
|---|---|---|
| Ansys Fluent | Fan BC、Porous Jumpが充実 | ファン特性曲線の直接入力 |
| STAR-CCM+ | パイプライン自動メッシュ | ポリヘドラルでエルボのメッシュが容易 |
| Ansys CFX | 結合型ソルバー | ターボ機械（ファン・ブロワ）との一体解析 |
| OpenFOAM | simpleFoam + pimpleFoam | 無償、スクリプト自動化 |
1Dネットワーク解析ツール
1D解析ツールって何ですか？ ダクト系統全体を1次元的に解くんですか？
そう。ダクト系統をノードとブランチのネットワークとして表現し、各要素の圧損係数を使って系統全体の流量配分を解く。3D CFDの前段階の概略設計や、系統全体のバランス調整に使われる。
| ツール | 開発元 | 特徴 |
|---|---|---|
| AFT Fathom | Applied Flow Technology | 液体配管ネットワーク解析 |
| AFT Arrow | Applied Flow Technology | 気体ダクト/配管ネットワーク解析 |
| Flownex | Flownex SE | 1Dシステムシミュレーション |
| Ductulator (ASHRAE) | ASHRAE | HVAC簡易ダクトサイジング |
1Dと3Dの使い分けはどう考えればいいですか？
Ansys FluentのFan BC設定
FluentのFan境界条件ってどうやって使うんですか？
Fan BCは面（Internal Face Zone）に設定し、ファンの圧力上昇-流量特性曲線を多項式または点テーブルで入力する。
圧力上昇は流速の関数として次のように定義する。
典型的な軸流ファン（風量3000 CMH、静圧200 Pa）の場合。
| 点 | 流量 [m³/h] | 静圧 [Pa] |
|---|---|---|
| 1 | 0 | 350 |
| 2 | 1000 | 320 |
| 3 | 2000 | 260 |
| 4 | 3000 | 200 |
| 5 | 4000 | 100 |
| 6 | 4500 | 0 |
ファンの運転点がCFDの結果として自動的に決まるわけですね。系統の圧損曲線との交点が運転点と。
その通り。系統抵抗をCFDで正確に計算するから、実際の運転点がファン特性曲線のどこにあるかが自動的に求まる。
ダクト系CFDでの「1Dネットワーク解析 vs フルCFD」の使い分け
ダクト内流れの商用ツール選択でよく議論になるのが「1Dネットワーク解析（例:AFT Fathom/Arrow）とフル3D CFDをどう使い分けるか」です。建物全体の数百本のダクト系を一括で最適化するなら1Dツールが圧倒的に速く、ASHRAE規格に準拠した設計書も自動生成できます。一方、特定のエルボー・分岐部など局所的な乱流・騒音問題はフルCFDでないと解析できません。実務の流れとしては「1Dで系全体のサイジング→問題箇所をフルCFDで深掘り」が効率的で、ANSYS FluentとSysweld、またはOpenFOAMを1D解析結果の境界条件として使う連成アプローチが増えています。
選定で最も重要な3つの問い
- 「何を解くか」：ダクト内流れに必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
先端トピックと研究動向
ダクト流れ解析の最新動向を教えてください。
最近のトレンドは3つある。
1. 空力騒音の予測（CAA連成）
ダクト系統の騒音をCFDで予測できるんですか？
ダンパーやベンド部で発生する空力騒音は、LESで渦構造を解像し、Ffowcs Williams-Hawkings（FW-H）の式で音源を評価する。
Ansys FluentではBroadband Noise Source Modelが利用でき、RANSベースの簡易騒音予測も可能だ。ただし精度はLES + FW-Hの方が遥かに高い。
RANSベースの騒音モデルはスクリーニング用で、最終評価はLESですね。
2. トポロジー最適化によるダクト形状設計
近年はトポロジー最適化をダクト形状設計に適用する研究が進んでいる。設計空間にDarcy抵抗を分布させ、流路として残す領域（抵抗ゼロ）と壁として機能する領域（抵抗無限大）を最適化する。
構造最適化と同じ考え方をCFDに適用するわけですね。
そうだ。目的関数を「圧力損失の最小化」として、制約条件に体積率を設定する。STAR-CCM+のAdjoint Solverや、COMSOL、TopOptのOpenFOAMプラグインで実装できる。
3. 1D-3Dハイブリッド解析
1Dと3Dを連成させる手法があるんですか？
FlownexやGT-SUITEなどの1Dシステムコードと3D CFDソルバーを連成させ、系統全体は1Dで解きつつ、問題の局所だけ3Dで詳細解析する手法だ。
大規模プラント配管（数百メートル）全体を3D CFDで解くのは非現実的だから、この方法は非常に実用的だ。Ansys TwinBuilderやSimcenter System Simulationと3D CFDの連成も商用で提供されている。
4. 機械学習による損失係数の予測
機械学習を使った研究もあるんですか？
非標準形状のダクト継手（3方向分岐、オフセット合流など）の損失係数を、CFDの大量計算結果からニューラルネットワークで学習させるアプローチが増えている。パラメータ（曲率比、断面比、分岐角度）を入力すると損失係数Kを瞬時に予測できる。
設計初期段階で大量のバリエーションを評価するのに便利ですね。
Surrogate ModelやPhysics-Informed Neural Network (PINN) と組み合わせると、CFDを一度も回さなくても概略設計ができるようになる。
ダクト形状の「角丸め」が圧力損失を劇的に変える
ダクト内流れの先端研究で面白いのが、矩形断面の角部曲率（コーナーR）が圧力損失に与える影響です。完全な直角コーナーがあると、そこで流れの剥離とコーナー渦が発生し、同じ断面積の円形ダクトと比べて摩擦係数が最大30%増加することがあります。LES（大渦シミュレーション）や DNS（直接数値シミュレーション）を用いた研究では、角部のコーナーRをダクト幅の5%確保するだけで剥離が大幅に抑制されることが示されています。空調・換気ダクトの設計では板金の曲げ加工コストとのバランスで「どこまでRを付けるか」が実務上の課題です。
トラブルシューティング
トラブルシューティング
ダクト内流れのCFDでよく遭遇する問題とその対策を教えてください。
パターン別に整理しよう。
1. 圧力損失が理論値と合わない
Darcy-Weisbachの計算値とCFD結果の差が大きいケースです。
チェックポイント:
- 壁面粗さの設定は正しいか（CFDのRoughness HeightとRoughness Constantが文献と整合しているか）
- Fluentの場合、Roughness Height $K_s$ はCFDの等価砂粗さで、Moody線図の $\varepsilon$ とは定義が異なる場合がある。Roughness Constant $C_s = 0.5$ のとき $K_s \approx 2\varepsilon$
- メッシュの壁面第一層がy+の要件を満たしているか
- 入口の乱流諸量（TI, $D_h$）が適切か
Roughness Heightの定義の違いは盲点ですね。
2. Fan BCでの収束不良
症状: Fan BCを設定すると残差が振動して収束しない。
対策:
- ファン特性曲線の点数を増やし、滑らかな曲線にする
- 初期条件として設計流量に近い速度場を与える
- Under-Relaxation Factorを下げる（Pressure: 0.2, Momentum: 0.5）
- Coupled Solverに切り替えるとロバストになる場合がある
3. 出口で逆流が発生
Pressure Outletで "Reversed flow on X faces" の警告が大量に出ます。
対策:
- 出口の断面を十分に下流に延長する（直管5Dh以上追加）
- Backflow条件に適切な乱流諸量を設定（デフォルトのままだと物理的でない逆流が入る）
- 出口をOutflow条件に変更する（ただし圧縮性流れでは使用不可）
4. T字分岐で流量比が実験と合わない
チェックポイント:
| 確認項目 | よくある問題 |
|---|---|
| 入口速度分布 | 一様プロファイルだと偏流を過小評価 |
| 乱流モデル | k-epsilonだと曲率効果が弱い |
| メッシュ | 分岐部の剥離領域が粗すぎる |
| 分岐後の管長 | 短すぎると出口BCの影響を受ける |
入口の速度分布を一様にしてしまうのは初心者がやりがちですね。
その通り。上流にエルボがある場合、そのエルボを含めてモデル化するか、実測プロファイルを入口条件に与えるべきだ。
5. 計算が非常に遅い
対策:
- 対称性がある場合は半分モデルを使用
- 長い直管部分は省略して、入口に完全発達プロファイルを与える
- 壁関数モデル（y+≒30）に切り替えてセル数を削減
- AMG（Algebraic Multigrid）の設定を確認（デフォルトで有効だが、Coupled Solverでは追加調整が有効）
直管部を省略するのは計算効率化の定番テクニックですね。助走区間さえ確保すれば結果に影響しないと。
そう。助走区間として10Dh程度の直管を入口に付ければ十分だ。
「流量が設計値の70%しか出ない」——ダクトトラブルあるある
ダクト系のトラブルで実務でよく遭遇するのが「シミュレーションでは計画流量が出るのに、実際に組み上げたら70〜80%しか出ない」パターンです。原因の多くはダクト接続部の「形状不連続」と「漏れ」。特にフランジ接合部のシール不良は、目に見えないわずかな隙間（0.5mm程度）でも系全体の圧損特性を大きく変えます。CFDトラブルシューティングでは、まず接続部の境界条件を見直し、次に分岐・合流部のローカルな圧損係数を実測値で補正するのが定石です。設計段階でCFDと配管設計ツールを連成させておくと、このギャップを事前に減らせます。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——ダクト内流れの問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告