フィルタ流れ解析
理論と物理
概要
先生！ フィルタ内の流れをCFDで解析するって、どういう場面で使うんですか？
エアフィルタ、オイルフィルタ、排ガス触媒のDPF（Diesel Particulate Filter）、浄水フィルタなど、多孔質媒体を通過する流れの圧力損失予測や粒子捕集効率の評価に使うんだ。
支配方程式
フィルタって多孔質ですよね。Navier-Stokes方程式をそのまま使うんですか？
フィルタの細孔を直接解像する（Pore-Scale Simulation）アプローチと、体積平均化した多孔質メディアモデルの2つがある。実務ではほとんど後者だ。
多孔質メディアモデルでは、Darcy-Forchheimer式で流れ抵抗を表す。
第1項がDarcyの粘性抵抗、第2項がForchheimerの慣性抵抗ですね。
その通り。$\alpha$ [m²]は透過率（permeability）、$C_2$ [1/m]は慣性抵抗係数だ。低速流（Re < 1）ではDarcy項が支配的で、高速流では慣性項が重要になる。
充填層の場合はErgun式が使われる。
$\phi$ が空隙率で $d_p$ が充填粒子径ですね。Ergun式とDarcy-Forchheimerの関係はどうなりますか？
Ergun式から多孔質パラメータを逆算できる。
粒子捕集のモデリング
フィルタの粒子捕集効率はどうモデル化するんですか？
単繊維捕集理論（Single Fiber Theory）が基本だ。繊維1本に対する各種捕集メカニズムの効率を足し合わせる。
| 捕集メカニズム | 支配粒径 | 効率の式（概略） |
|---|---|---|
| さえぎり（Interception） | > 0.5 um | $E_R = \frac{R^2}{(1+R) Ku}$ |
| 慣性衝突（Impaction） | > 1 um | $E_I \propto Stk^2$ |
| ブラウン拡散（Diffusion） | < 0.3 um | $E_D \propto Pe^{-2/3}$ |
| 重力沈降 | > 5 um | $E_G = G/(Ku \cdot Re_f)$ |
| 静電気 | 全粒径 | 帯電量に依存 |
MPPS（Most Penetrating Particle Size）が0.1〜0.3 um付近にあるのは、拡散と慣性のギャップですよね。
その通り。HEPAフィルタの捕集効率曲線がV字型になるのはそのためだ。CFDではDPMで粒子軌道を追跡し、繊維層への到達を判定する方法が使われる。
実務上の注意点
フィルタCFDで特に気をつけるべきことは？
フィルタ捕集の「3つのメカニズム」——大きすぎても小さすぎても捕まる
フィルタ流れの理論で面白いのが、粒子捕集のメカニズムが粒径によって切り替わる点です。大きな粒子（数μm以上）は「慣性衝突」——流線をそれて繊維に激突する。中間サイズ（0.5〜2μm）は「さえぎり捕集」——流線に乗ったまま繊維に接触。超微細粒子（0.1μm以下）は「拡散捕集」——ブラウン運動でジグザグに動き偶然繊維に付着。つまり最も捕集しにくいのは「慣性も小さくブラウン運動も弱い」0.1〜0.3μmの過渡帯です。この3機構を組み合わせたSingle Fiber Theory（単一繊維理論）が、フィルタ設計の理論的基盤となっています。
各項の物理的意味
- 時間項 $\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の具体的な実装方法を教えてください。
多孔質メディアモデルの実装方法は2つある。Volumetric Porous Zone（体積全体に抵抗を分布）と、Porous Jump（薄い面に圧力ジャンプを設定）だ。
Porous Zone vs. Porous Jump
| 方式 | 適用 | メリット | デメリット |
|---|---|---|---|
| Porous Zone | 厚いフィルタ（充填層、触媒層） | フィルタ内部の流速分布が得られる | メッシュが必要 |
| Porous Jump | 薄いフィルタ（プリーツフィルタ） | メッシュ不要、計算が軽い | 内部流れの情報なし |
HEPAフィルタのような薄いフィルタはPorous Jumpでいいんですね。
そう。フィルタ厚さが数cm程度で、厚さ方向の流れの詳細が不要ならPorous Jumpが効率的だ。DPFのように厚さが100mm以上ある場合はPorous Zoneを使う。
Fluentでの多孔質設定手順
FluentでPorous Zoneを設定する具体的な手順を教えてください。
1. Cell Zoneでフィルタ領域を選択し、Porous Zoneを有効にする
2. Direction-1 Vectorをフィルタ法線方向に設定
3. Viscous Resistance（1/α [1/m²]）とInertial Resistance（C₂ [1/m]）を入力
4. Porosityを入力（デフォルト1.0は全空間が流路の意味で、物理的な空隙率に変更）
抵抗値の算出例を示す。HEPAフィルタで面風速0.5 m/sのとき圧損250 Pa、厚さ65 mmの場合。
低速域でDarcy項が支配的と仮定すると。
Viscous Resistanceが約4.3e8ですね。Inertial Resistanceは別の面風速でのデータから求めるんですか？
そう。2つ以上の面風速-圧損データがあれば、連立方程式でDarcy項とForchheimer項を分離できる。面風速0.3 m/sで圧損140 Pa、0.7 m/sで圧損380 Paのようなデータを使う。
DPMとの連成
フィルタでの粒子捕集をDPMでどうモデル化しますか？
Porous Zone内でDPM粒子の壁面条件をTrapに設定し、フィルタ繊維表面への衝突を捕集とみなす方法がある。ただしこれは簡略モデルで、より正確には捕集確率を粒径の関数としてUDFで実装する。
フィルタの総合捕集効率は次式で表される。
$\alpha_f$ がフィルタの充填率、$E_f$ が単繊維効率、$d_f$ が繊維径ですね。
そう。CFDでは各セルごとに局所の面風速から単繊維効率を計算し、DPM粒子の捕集確率として適用する。
メッシュの注意点
フィルタ領域のメッシュはどうしますか？
プリーツフィルタはひだの形状が複雑ですが、全部モデル化するんですか？
全プリーツを3Dモデル化するのは現実的でない場合が多い。1〜3プリーツの周期境界モデルで代表特性を求め、全体モデルでは等価な多孔質面として扱うのが実用的だ。
フィルタをCFDでモデル化する「多孔質媒体モデル」の中身
フィルタ流れ解析でよく使われる「多孔質媒体モデル（Porous Media Model）」は、フィルタ層内の詳細な繊維構造を無視して、体積全体に分布する圧力損失源として扱う手法です。Fluent・OpenFOAMともにDarcy則ベースの実装があり、圧力損失係数（粘性抵抗係数 $\alpha$ と慣性抵抗係数 $C_2$）を実測のΔP-流速曲線からフィッティングします。注意点は「フィルタの実際の厚み方向の不均一性」を無視していること——実使用で目詰まりが局所的に進行すると、モデルの精度が低下します。このため経時劣化を含む解析では粒子追跡（DPM）と連成させる必要があります。
風上差分（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: エアフィルタユニットの圧損最適化
HVACのエアハンドリングユニット内のフィルタですね。
フィルタバンクの前面にディフューザーがあり、フィルタ面への風速分布が不均一になるケースだ。CFDで評価して均一化を図る。
解析手順:
1. AHU内部の形状をCADでモデル化（ディフューザー、フィルタバンク、コイル）
2. フィルタをPorous Jumpで設定
3. 入口にファンBCまたはVelocity Inletを設定
4. フィルタ面上の風速分布を評価
5. フィルタ面の風速均一度指数を算出
均一度指数 $\gamma$ が1に近いほど均一なんですね。目標値はどのくらいですか？
$\gamma > 0.9$ が望ましい。0.8未満だとフィルタの局所寿命が短くなり、高風速部の捕集効率が低下する。
ケース2: 自動車DPFの圧損解析
ディーゼル車のDPFですね。
DPFはウォールフローモノリス構造で、隣接するチャネルの壁を通して排ガスを濾過する。CFDでは1チャネルの2D軸対称モデルで代表特性を求め、全体は1Dモデルと連成させる方法が一般的だ。
DPFの圧損は次式で表される。
| 圧損成分 | 典型値 [kPa] | 支配パラメータ |
|---|---|---|
| 入口チャネル摩擦 | 1〜3 | チャネル長、流速 |
| 壁面透過 | 3〜8 | 壁厚、透過率、煤堆積量 |
| 出口チャネル摩擦 | 0.5〜2 | チャネル長、流速 |
| 入口/出口損失 | 0.3〜1 | チャネル開口率 |
ケース3: 集塵機のバグフィルタ評価
工場の集塵機ですね。
バグフィルタ（布製フィルタ）の圧損はダスト層の厚さに大きく依存する。初期圧損とダスト堆積後の圧損を区別してモデル化する必要がある。
$W$ がダストの面積質量密度 [kg/m²] ですね。時間とともに増加する。
そう。逆洗（パルスジェット）後に $W$ がリセットされ、再び堆積が始まるサイクルをモデル化する場合はTransient解析+UDFが必要になる。
よくある失敗と対策
| 失敗パターン | 原因 | 対策 |
|---|---|---|
| 圧損が異常に大きい | Viscous Resistance値の単位ミス | 1/α [1/m²]であることを確認 |
| フィルタ面で流れが均一すぎる | 上流の偏流が多孔質抵抗で均される | 抵抗値が実際より高くないか確認 |
| DPM捕集率が0%または100% | 壁面BCの設定ミス | Trap/Escape/Reflectの使い分けを確認 |
| 計算が発散 | Porous Zoneの抵抗が急激すぎる | 初期から抵抗を入れずに流れ場を発達させてから投入 |
上流の偏流がフィルタの抵抗で均されてしまうのは、フィルタの抵抗が高すぎると起きやすいですよね。
その通り。フィルタの圧損が動圧に比べて十分大きい場合（$\Delta p_{filter} \gg \rho V^2/2$）、フィルタ面の流速分布は自然に均一化される。逆に言えば、低圧損のプレフィルタでは偏流が問題になりやすい。
プリーツフィルタの「プリーツ深さ」——多すぎると逆効果になる？
フィルタ設計の実践でよく誤解されるのが「プリーツ（ひだ折り）は深いほど有効面積が増えていい」という直感です。確かに浅いプリーツから深いプリーツにすると、有効ろ過面積が増えて初期圧損が下がります。ところがプリーツが深すぎると谷の部分に流れが届かなくなり、実効面積が理論値を大きく下回ります。CFDで可視化するとプリーツ谷底が「デッドゾーン」になっている様子が明確に分かり、最適なプリーツ深さ・間隔の比（アスペクト比）が存在することが実証されています。自動車エアフィルタの業界では、このCFD最適化で交換寿命を20〜30%延ばした事例もあります。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ CFDでも入口条件のミスは下流全体に波及します。
ソフトウェア比較
商用ツール比較
フィルタCFDに適したツールはどれですか？
多孔質メディアモデルの実装が充実しているツールが適している。
| ツール名 | 多孔質モデル | DPM/粒子追跡 | DPF専用機能 |
|---|---|---|---|
| Ansys Fluent | Porous Zone / Jump | DPM充実 | なし（UDFで対応） |
| STAR-CCM+ | Porous Region | Lagrangian | なし |
| Ansys CFX | Porous Domain | Particle Transport | なし |
| OpenFOAM | porousSimpleFoam | DPMFoam | コミュニティ実装あり |
| COMSOL | Darcy's Law / Brinkman | Particle Tracing | なし |
| AVL FIRE | Porous Media | Spray/DPM | DPFモジュールあり |
| Exothermia Suite | - | - | GPF/DPF専用 |
DPF専用のツールがあるんですね。
AVL FIREにはDPFモジュールが標準搭載されていて、煤の堆積・再生（酸化）をカップルで計算できる。Exothermia（旧Exothermia/axiom）はDPF/GPF専用の1D+2Dツールで、自動車業界で広く使われている。
用途別の推奨ツール
| 用途 | 推奨ツール | 理由 |
|---|---|---|
| HVAC/空調フィルタ | Fluent / STAR-CCM+ | Porous Jumpの簡便さ |
| 自動車DPF/GPF | AVL FIRE / Exothermia | 専用の煤堆積モデル |
| プロセスフィルタ（化学プラント） | Fluent / OpenFOAM | Ergun式の柔軟な設定 |
| 浄水・膜フィルタ | COMSOL | 膜透過モデル（Kedem-Katchalsky式） |
| 集塵機 | Fluent | DPMで粒径別捕集効率を評価 |
COMSOLが膜フィルタに強いのは、マルチフィジックスの得意分野ですね。
そう。COMSOLにはDarcy's Law、Brinkman Equations、Free and Porous Media Flowなど、多孔質流れのためのインタフェースが複数用意されていて、化学反応や物質輸送との連成も容易だ。
OpenFOAMでのフィルタ解析
OpenFOAMでフィルタ解析をする場合のソルバーは何を使いますか？
porousSimpleFoam（定常）またはsimpleFoam + fvOptionsで多孔質抵抗を追加する。設定はconstant/fvOptionsにDarcy-ForchheimerのパラメータをDictionaryで記述する。
設定例:
```
porosity
{
type DarcyForchheimer;
active true;
cellZone filter;
d (1e8 1e8 1e8); // 1/alpha [1/m^2]
f (100 100 100); // C2 [1/m]
}
```
d がViscous Resistance（1/alpha）、f がInertial Resistance（C2）ですね。方向ごとに異なる値を設定できると。
そう。異方性フィルタ（例えばプリーツフィルタの法線方向と平行方向で抵抗が異なる場合）にも対応できる。
フィルタ解析でOpenFOAMを使うときの「porousMediaモデル」の選び方
OpenFOAMでフィルタ流れを解析するとき、porous zoneの定義に使う抵抗係数をどう決めるかが初学者の最初の壁です。DarcyForchheimer則を使う場合、粘性抵抗係数 $d$（単位: m⁻²）と慣性抵抗係数 $f$（単位: m⁻¹）の2パラメータが必要で、これはフィルタ単体の実測Δp-u曲線から最小二乗フィッティングで求めます。商用ツール（ANSYS Fluent）は同様の設定をGUI上で行えるため習得が早い一方、OpenFOAMは物理モデルの透明性が高くカスタマイズ性に優れます。粒子追跡（DPM）との連成もOpenFOAMのほうが細かい制御が可能で、研究用途では強力な選択肢です。
選定で最も重要な3つの問い
- 「何を解くか」：フィルタ流れ解析に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
先端トピックと研究動向
フィルタCFDの最新研究はどんな方向に進んでいますか？
主に3つの方向がある。
1. Pore-Scaleシミュレーション
フィルタの細孔構造を直接解像する方法ですね。
X線CT（マイクロCT）でフィルタの3D微細構造を撮像し、それを直接CFDメッシュに変換してNavier-Stokesを解く。Lattice Boltzmann法（LBM）が良く使われる。
LBMはメッシュ生成が不要で、複雑な多孔質構造に適しているんですよね。
そう。CTのボクセルデータをそのまま計算格子として使えるのがLBMの大きな利点だ。Palabos（OSS）やPowerFLOW（商用）がLBMソルバーとして有名だ。
2. CFD-DEMによるダスト堆積の動的シミュレーション
CFD-DEM連成では、粒子の堆積過程を個々の粒子レベルで追跡する。粒子間の接触力、付着力（van der Waals力）を考慮して、ダストケーキの形成過程を時間発展的に計算できる。
ダストケーキの構造がフィルタ性能に影響するわけですね。
その通り。ダストケーキの空隙率は堆積条件によって変わり、圧損と捕集効率の両方に影響する。CFD-DEMではEDEM（商用）やLIGGGHTS（OSS）をFluentやOpenFOAMと連成させる。
3. フィルタ形状の最適化
フィルタの形状最適化ってどういうことですか？
プリーツフィルタのプリーツ深さ・間隔・枚数の最適化が典型例だ。プリーツが深いほど濾過面積が増えて面風速が下がるが、プリーツ間の通路が狭くなって入口損失が増える。このトレードオフをCFDパラメトリックスタディで最適化する。
最適プリーツ数の目安として、Chen & Pui (1995) の理論解がある。
$$ N_{opt} = \left(\frac{3 Q \mu}{2 W h_p (\Delta p_{media}/L)}\right)^{1/3} $$
$N_{opt}$ が最適プリーツ数で、$h_p$ がプリーツ深さ、$W$ がフィルタ幅ですね。
CFD-DEM連成では、粒子の堆積過程を個々の粒子レベルで追跡する。粒子間の接触力、付着力（van der Waals力）を考慮して、ダストケーキの形成過程を時間発展的に計算できる。
ダストケーキの構造がフィルタ性能に影響するわけですね。
その通り。ダストケーキの空隙率は堆積条件によって変わり、圧損と捕集効率の両方に影響する。CFD-DEMではEDEM（商用）やLIGGGHTS（OSS）をFluentやOpenFOAMと連成させる。
フィルタの形状最適化ってどういうことですか？
プリーツフィルタのプリーツ深さ・間隔・枚数の最適化が典型例だ。プリーツが深いほど濾過面積が増えて面風速が下がるが、プリーツ間の通路が狭くなって入口損失が増える。このトレードオフをCFDパラメトリックスタディで最適化する。
最適プリーツ数の目安として、Chen & Pui (1995) の理論解がある。
$N_{opt}$ が最適プリーツ数で、$h_p$ がプリーツ深さ、$W$ がフィルタ幅ですね。
トラブルシューティング
トラブルシューティング
フィルタCFDでよく遭遇する問題を教えてください。
パターン別に整理しよう。
1. Porous Zoneで速度が非物理的に大きい
症状: フィルタ領域内で速度が入口速度の何倍にもなる。
原因と対策:
- Porosityの設定を確認。デフォルトの1.0（=全空間が流路）のままだと、Superficial Velocityベースの計算では正しいが、Physical Velocityで表示すると空隙率で割った値になる
- FluentではSolver > GeneralのPorous Media Formulation（Superficial vs Physical Velocity）の設定を確認
Superficial Velocityは空のダクトを流れる場合の速度で、Physical Velocityは実際に細孔を通る速度ですよね。表示の問題だけですか？
その通り。Superficial Velocity formulation（デフォルト）ではDarcy-Forchheimerの式にSuperficial Velocityを使うから、抵抗パラメータもSuperficial基準で入力する必要がある。
2. 圧損が実測と合わない
チェックリスト:
| 確認項目 | よくある問題 |
|---|---|
| 抵抗パラメータの単位 | 1/α [1/m²] と α [m²] の取り違え |
| フィルタ厚さ | Porous Zoneの厚さとパラメータの整合性 |
| 面風速の範囲 | パラメータをフィッティングした速度範囲外で使用 |
| フィルタの目詰まり | 新品データで計算しているが実際は汚れている |
| 異方性 | 等方的に設定しているが実際は異方性がある |
3. Porous Zone周辺で残差が振動
対策:
- 多孔質領域と非多孔質領域の境界で急激な圧力変化が起きるため、メッシュ遷移を滑らかにする
- 初期条件で流れ場を予め発達させてから、段階的に多孔質抵抗を増やす
- Under-Relaxation Factorを下げる（Pressure: 0.2, Momentum: 0.5）
4. DPMの粒子がPorous Zoneを透過してしまう
本来捕集されるべき粒子がフィルタを素通りするケースですね。
原因: Porous ZoneはDPM粒子に対して壁として機能しない。粒子は多孔質領域の流れに沿って移動するだけだ。
対策:
- フィルタ下流面にDPM壁面条件（Trap）を設定し、UDFで粒径依存の捕集確率を実装
- あるいはフィルタ出口面にDPM Escape条件を設定し、出口で粒子数をカウントして透過率を算出
5. Porous Jump面の法線方向が逆
症状: フィルタを通過すると圧力が上昇する（通常は低下するはず）。
原因: Porous Jump面の法線ベクトルが流れと逆向き。FluentではFace Zone Directionを確認し、必要に応じてFlipする。
法線の向きは見落としやすいですね。圧損が負になったら真っ先にチェックすべきと。
その通り。Porous Jumpの設定後、必ずフィルタ前後の圧力差を確認して、物理的に正しい方向に圧損が発生しているかを検証すること。
「フィルタの端漏れ」——CFDで見えても実機で再現できない落とし穴
フィルタCFDのトラブルシューティングで悩ましいのが「シール（枠の密着）部分の端漏れ」です。CFDモデルでは理想的なシールを仮定するため端漏れは生じませんが、実機では枠の変形やガスケット不良で1〜2mm の隙間が生じ、フィルタ効率が設計値から大きく外れることがあります。大型産業フィルタの設置後に効率不良が判明するケースの相当数はこれが原因です。トラブル対応策としては、CFDに意図的なシール不良をモデル化し、どの程度の漏れで効率が何%低下するか「感度解析」を実施しておくと、現場での原因特定と改善策の議論が格段にスムーズになります。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——フィルタ流れ解析の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告