クリーンルーム気流解析
クリーンルーム気流の理論基礎
概要
先生！ クリーンルームの気流解析って、半導体工場とかで使うやつですよね？ どういう物理が関わってくるんですか？
クリーンルーム気流解析は、室内の清浄度を維持するために一方向流（ユニフロー）や乱流置換方式の気流パターンをCFDで予測する技術だよ。ISO 14644-1で定義される清浄度クラス（Class 1〜Class 9）を達成するため、FFU（Fan Filter Unit）からの吹き出し気流がパーティクルをどう搬送・排出するかを解析するんだ。
なるほど、清浄度クラスを数値的に検証できるってことですね。
支配方程式
気流解析で使う方程式はNavier-Stokesですよね？ パーティクルの追跡はどうするんですか？
まず連続気相はRANSベースのNavier-Stokes方程式で解く。非圧縮性を仮定する場合が多い。
連続の式とNavier-Stokes方程式はこうなる。
FFUのフィルタ部分は多孔質メディアモデルで表現する。Darcy-Forchheimer抵抗が入る。
$\alpha$ が透過率で $C_2$ が慣性抵抗係数ですね。フィルタのカタログ値から逆算できますか？
HEPAフィルタの場合、面風速0.45 m/sで圧力損失が約250 Paというのが典型値だ。これとフィルタ厚さから $\alpha$ と $C_2$ を算出する。パーティクル追跡にはDPM（Discrete Phase Model）を使い、粒子の運動方程式を解く。
ブラウン力まで入れるんですか。サブミクロン粒子だとブラウン運動が効いてくるわけですね。
そう、0.1 um以下の粒子ではブラウン拡散が支配的になる。Cunningham補正係数 $C_c$ も必要だ。
乱流モデルの選択
クリーンルームの乱流モデルは何が適していますか？
一方向流クリーンルームでは層流に近い領域と乱流領域が共存するから、SST $k$-$\omega$ モデルが推奨される。低Re数領域の壁面処理が自然にできるからだ。
| 乱流モデル | 推奨度 | 特徴 |
|---|---|---|
| SST k-omega | 高 | 低Re壁面処理、剥離予測に強い |
| Realizable k-epsilon | 中 | 汎用的だが壁関数が必要 |
| RNG k-epsilon | 中 | 旋回流に若干良好 |
| LES (Smagorinsky) | 非常に高（計算コスト大） | 非定常渦構造を直接解像 |
半導体ファブの実案件ではSST k-omegaが多いんですね。LESは研究用途ですか？
その通り。ただし最近はクリーンルーム内の人体動作による乱れの非定常解析でLESが使われることも増えてきている。
実務上の注意点
現場で気をつけるべきポイントを教えてください。
人体の発塵モデルまでCFDに入れるんですね。クリーンルーム特有の話で勉強になります。
HEPAフィルタの「99.97%」という数字の正体
クリーンルーム気流の理論を学ぶとき、避けて通れないのがHEPAフィルタの捕集効率です。「0.3μmの粒子を99.97%以上捕集」というのは最悪値、つまりこのサイズが最も素通りしやすいという意味です。なぜ0.3μm？ それ以下の超微細粒子はブラウン運動が支配的で繊維に衝突しやすく、それ以上の粒子は慣性力で捕まりやすい。ちょうど0.3μm付近が「慣性も小さくブラウン運動も弱い」過渡領域となるため、最難関サイズになります。クリーンルームCFDでパーティクル追跡を行う際は、この分布特性を理解した上で粒径範囲を設定することが重要です。
クリーンルーム気流の数値計算手法
数値手法の詳細
クリーンルームCFDを実際に解くとき、有限体積法ですよね？ 具体的な離散化スキームはどう選ぶんですか？
クリーンルーム気流は低マッハ数の非圧縮性流れだから、圧力ベースソルバーを使う。SIMPLE系アルゴリズム（SIMPLE, SIMPLEC, PISO）で圧力-速度連成を解くんだ。
圧力-速度連成
SIMPLEとSIMPLECの使い分けはありますか？
定常解析ならSIMPLEC（圧力補正の緩和が不要で収束が速い）、非定常解析ならPISO（時間ステップ毎の反復が少ない）が推奨だ。Coupled Solverも選択肢だが、メモリ消費が大きい。
| アルゴリズム | 定常/非定常 | 特徴 |
|---|---|---|
| SIMPLE | 定常 | 基本手法、緩和係数の調整が必要 |
| SIMPLEC | 定常 | 収束が速い、クリーンルームに推奨 |
| PISO | 非定常 | 人体動作の非定常解析向き |
| Coupled | 両方 | ロバストだがメモリ2〜3倍 |
空間離散化
対流項のスキームはどれがいいですか？
クリーンルームは低速流れ（0.3〜0.5 m/s程度）なので数値拡散が問題になりやすい。Second Order Upwind以上を推奨する。
QUICKスキームは四面体メッシュでは使えないんですよね。
その通り。QUICKは構造格子か六面体メッシュ前提だ。ポリヘドラルメッシュの場合はSecond Order Upwindが無難だ。
DPMの実装詳細
パーティクル追跡の具体的な設定を教えてください。
DPMでは粒子軌道を時間積分で追跡する。クリーンルーム解析での典型設定はこうだ。
| パラメータ | 推奨値 | 備考 |
|---|---|---|
| 粒子径分布 | Rosin-Rammler (0.1〜10 um) | ISO 14644-1の対象粒径 |
| 粒子数 | 10,000以上/インジェクション面 | 統計的信頼性 |
| 積分手法 | Trapezoidal | 精度と速度のバランス |
| ブラウン力 | ON (dp < 1 um) | サブミクロン粒子必須 |
| Saffman揚力 | ON | 壁面近傍の挙動改善 |
| 壁面条件 | Trap/Reflect | 堆積 vs. 反発 |
粒子数が1万以上って結構多いですね。計算時間への影響は？
DPMは一方向連成（One-Way Coupling）なら気相の計算後にポスト的に追跡するだけだから、追加コストは全体の10〜20%程度だ。クリーンルームの粒子濃度は低いので一方向連成で十分。
メッシュ戦略
クリーンルームは大空間ですが、メッシュ数の目安はどのくらいですか？
典型的な半導体ファブ1ベイ（10m x 20m x 3m）で500万〜2000万セルが目安だ。FFU吹出面とウェハ周辺は局所細分化が必須で、最小セルサイズ5〜10 mm程度になる。
床下プレナムのメッシュも結構細かくする必要があるんですね。圧力損失に影響しますか？
床下プレナムは開口率が約25%のグレーチング床で大きな圧損が発生する。多孔質ジャンプ条件で簡略化する場合もあるが、局所的な偏流が問題になる場合はフルモデル化が必要だ。
クリーンルームCFDで鍵を握るのは「層流維持」の手法選び
クリーンルーム気流解析の数値手法でよく使われるのが、低レイノルズ数乱流モデル（Low-Re k-ε）と層流モデルの使い分けです。FFU（ファンフィルタユニット）からの吹き出し直下はRe数が低く層流に近い性質があり、ここを乱流モデルで解くと過剰な拡散が生じてパーティクル挙動が非現実的になります。実務では吹き出し口付近は層流扱いとし、作業エリアへの流れ込みを境に乱流モデルに切り替えるゾーン分割アプローチが有効です。ラグランジュ粒子追跡（DPM）を組み合わせれば、どの作業位置でコンタミが発生しやすいかも定量評価できます。
クリーンルーム気流の実務適用
実践ガイド
先生、クリーンルームCFDの実務的な解析フローを教えてください。
全体の流れは、(1)レイアウト情報の収集、(2)CADモデル構築、(3)メッシュ生成、(4)境界条件設定、(5)求解、(6)清浄度評価の順だ。
解析フロー
まずCADモデルはどこまで詳細にするんですか？
1. 形状簡略化のポイント
- FFUは多孔質面として吹出風速を与える（内部構造は不要）
- 製造装置はブロック形状で簡略化（排気口は正確に配置）
- 人体は簡易円柱モデル（直径0.4m、高さ1.7m）+ 発熱量75 W
- 配管・ケーブルラックは省略可能（流路閉塞が5%未満の場合）
2. 境界条件の設定
- FFU吹出面: Velocity Inlet（0.45 m/s、一様）+ 温度（22℃ typ.）
- リターン口/グレーチング床: Pressure Outlet（ゲージ圧0 Pa）
- 装置排気口: Velocity Inlet（負の値）またはOutflow
- 壁面: 断熱No-Slip（天井・壁）、温度指定（発熱装置表面）
FFUの吹出し温度って、空調設計値をそのまま使えばいいですか？
基本的にはそうだ。ただし装置発熱による局所的な温度上昇を検証するのがCFDの目的の一つだから、吹出温度は空調機出口温度（通常18〜20℃）を設定して、室内温度分布を計算結果として得る方がいい。
清浄度の評価手法
CFD結果からISO清浄度クラスをどうやって判定するんですか？
DPMの粒子追跡結果から、評価点における粒子濃度を算出する。具体的には次の手順だ。
1. 発塵源（人体、装置）から粒子をリリース
2. 評価点（ウェハ上方300mm等）を通過する粒子数をカウント
3. 通過粒子数と発塵率から個数濃度 [個/m³] を算出
4. ISO 14644-1の表と照合してクラス判定
| ISOクラス | 0.1 um [個/m³] | 0.5 um [個/m³] | 用途例 |
|---|---|---|---|
| Class 1 | 10 | - | 最先端半導体リソ |
| Class 3 | 1,000 | 35 | 半導体前工程 |
| Class 5 | 100,000 | 3,520 | 半導体後工程 |
| Class 7 | - | 352,000 | 医薬品製造 |
定量的な評価ができるんですね。ただ粒子数の統計誤差が心配です。
その通り。DPMの粒子追跡数を十分に確保することが重要だ。粒子数を2倍にして結果が変わらないことを確認する「粒子数独立性の検証」は必須だよ。
よくある失敗と対策
初心者がやりがちなミスを教えてください。
| 失敗パターン | 原因 | 対策 |
|---|---|---|
| 面風速が不均一すぎる | FFU直下のメッシュが粗い | FFU面に最低10x10セルを確保 |
| 温度が室内で均一になる | 数値拡散が大きい | Second Order以上を使用 |
| パーティクル評価が安定しない | DPM粒子数不足 | 最低10,000粒子/ソース |
| 床下リターンが偏る | グレーチング開口率の設定ミス | 多孔質ジャンプのK値を確認 |
| 非定常計算が発散 | 時間ステップが大きすぎる | CFL数 < 1を目安に設定 |
数値拡散でせっかくの温度分布がなまってしまうのは痛いですね。離散化スキームの選択が大事だと。
その通りだ。クリーンルームは低速流れだから、一次風上だと数値拡散で温度分布がほぼ均一になってしまう。必ずSecond Order Upwind以上を使うこと。
「人が歩くだけで0.3m/sの乱流」——クリーンルーム設計の最大の敵
クリーンルーム気流解析の実践で最も見落とされがちなのが「人体の発熱と動作が生み出す乱流」です。研究によれば、人が歩くだけで体の周囲に0.2〜0.5m/sの乱気流が発生し、FFUから吹き出す0.4〜0.5m/sの清浄ダウンフローを乱します。実務ベストプラクティスとして、作業員の人体形状（熱源33〜35W/人）をCFDモデルに組み込み、作業動線ごとに汚染拡散をシミュレーションする「動的クリーンルーム解析」が普及してきました。入退室の動線設計や防護服の義務化判断にも活用されます。
クリーンルーム気流のソフトウェア比較
商用ツール比較
クリーンルームCFDに適したツールってどれですか？ 専用ソフトもあると聞きましたが。
汎用CFDソルバーに加えて、クリーンルーム専用の簡易ツールも存在する。それぞれの特徴を整理しよう。
対応ツール一覧
| ツール名 | 開発元 | 特徴 | クリーンルーム適性 |
|---|---|---|---|
| Ansys Fluent | Ansys Inc. | DPMが充実、多孔質モデル完備 | 非常に高い |
| Simcenter STAR-CCM+ | Siemens | ポリヘドラルメッシュ、自動化 | 非常に高い |
| Ansys CFX | Ansys Inc. | 結合型ソルバー | 高い |
| OpenFOAM | OSS | icoFoam/simpleFoam + DPM | 高い（設定に知識必要） |
| COMSOL Multiphysics | COMSOL AB | Particle Tracingモジュール | 中（小規模向き） |
| FlowDesigner | アドバンスドナレッジ研究所 | 建築空調特化、日本製 | 高い（空調設計向き） |
| Stream | クレイドル（MSC） | 建築・空調CFD | 高い（国内実績豊富） |
FlowDesignerやStreamは日本のクリーンルーム設計で実績が多いんですか？
そうだ。FlowDesignerは空調設計者向けのUIが充実していて、FFUや空調機の設定が直感的にできる。Streamも国内の建設会社で広く使われている。ただし高度なDPM解析や多相流が必要な場合は、FluentやSTAR-CCM+の方が柔軟だ。
ツール選定の判断基準
どういう基準で選べばいいですか？
用途に応じた選定が重要だ。
| 用途 | 推奨ツール | 理由 |
|---|---|---|
| 半導体ファブ全体の気流検証 | Fluent / STAR-CCM+ | 大規模DPM、HPC対応 |
| 空調設備の設計検討 | FlowDesigner / Stream | 設備設計者向けUI |
| ミニエン（局所清浄空間）の詳細評価 | Fluent / CFX | 精密なDPM + Brownian |
| 研究・パラメトリックスタディ | OpenFOAM | ライセンス無料、自動化 |
| 簡易的なレイアウト検討 | COMSOL | マルチフィジックス連成 |
Ansys Fluentでの設定例
Fluentで設定する場合の具体的な手順を教えてください。
1. Viscous Model: SST k-omega、Low-Re Corrections ON
2. Energy Equation: ON（温度分布評価時）
3. Porous Media: FFU面にFace Porous Jumpを設定
4. DPM: Injection → Surface injection（発塵面から）
5. Species Transport: 必要に応じてトレーサーガスの拡散解析
FFUのPorous Jump設定値の目安はこうだ。
| パラメータ | HEPA (H13) | ULPA (U15) |
|---|---|---|
| Face Permeability [m²] | 5.0e-10 | 2.0e-10 |
| Pressure Jump Coefficient [1/m] | 200 | 500 |
| Medium Thickness [m] | 0.065 | 0.065 |
フィルタグレードによって透過率が全然違うんですね。カタログの圧損曲線からフィッティングするわけですか。
その通り。メーカーのカタログに記載されている「面風速 vs. 圧力損失」のデータから最小二乗法でDarcy項（粘性抵抗）とForchheimer項（慣性抵抗）を分離する。
クリーンルーム専用ツールか汎用CFDか——選定の決め手
クリーンルーム気流解析のツール選びでよく聞かれるのが「専用ソフトと汎用CFDのどっちがいい？」という問いです。Flomerics（現Mentor Simcenter）のFloVENTやIDEAS/AirPackはクリーンルーム・データセンター向けの専用設定が豊富で、非エキスパートでも短時間で解析できます。一方ANSYS FluentやOpenFOAMなどの汎用CFDは自由度が高い分、境界条件の作り込みに時間がかかります。実務では「新棟設計の初期検討→専用ツール」「問題調査・詳細最適化→汎用CFD」という使い分けが多く、どちらか一方で完結させようとすると非効率になりがちです。
クリーンルーム気流の先端研究
先端トピックと研究動向
クリーンルームCFDの最新トレンドってどんな感じですか？
半導体の微細化が進むにつれて要求される清浄度が上がり、従来の定常RANS解析では捉えきれない現象が問題になってきている。
LESによる非定常気流評価
LES（Large Eddy Simulation）をクリーンルームに適用する研究が増えていると聞きました。
その通り。作業者が動くときの発塵は本質的に非定常現象だから、RANSの時間平均では捉えきれない瞬間的な高濃度イベントが問題になる。LESではSmagorinskyモデルやWALE（Wall-Adapting Local Eddy-Viscosity）モデルを使って渦構造を直接解像する。
Smagorinsky定数 $C_s$ はいくつくらいが適切ですか？
室内気流では $C_s = 0.1$〜$0.12$ が一般的だ。Dynamic Smagorinskyモデルなら自動調整されるから、定数を気にしなくて済む。
AMC（Airborne Molecular Contamination）解析
パーティクルだけじゃなくて、分子レベルの汚染も評価するんですか？
最先端EUVリソグラフィでは、有機分子（アウトガス）がマスクやレンズに付着して歩留まりを低下させる。これをAMC解析と呼ぶ。Species Transport方程式で有機ガスの拡散・対流を解くんだ。
$Y_i$ は化学種の質量分率で、$D_i$ は拡散係数ですね。ソース項 $S_i$ はアウトガス発生量ですか？
そう。装置やシール材からのアウトガス発生率 [ng/(cm² hr)] をソース項として与える。活性炭フィルタの除去効率もモデル化する必要がある。
デジタルツインとリアルタイム監視
クリーンルームのデジタルツインって実用化されていますか？
一部の半導体ファブでは、パーティクルカウンターや温湿度センサーのリアルタイムデータとCFDモデルを連携させたデジタルツインが運用されている。
ROMってどうやって作るんですか？
POD（Proper Orthogonal Decomposition）で支配的なモードを抽出して、少数のモード係数で気流場を近似する方法が一般的だ。100ケース程度のCFD結果からROMを構築すれば、FFU風量や装置発熱を変えたときの気流をミリ秒単位で予測できる。
クリーンルーム気流のトラブル対応
トラブルシューティング
クリーンルームCFDで計算がうまくいかないとき、どう対処すればいいですか？
よくある問題をパターン別に整理しよう。
1. 収束しない・残差が振動する
特に定常計算で残差が下がらないケースが多いです。
考えられる原因と対策:
- Boussinesq浮力と圧力補間の不整合: 浮力項があるときはPRESTO!またはBody Force Weighted圧力補間を使う。Standard圧力補間だと残差が振動する
- 多孔質ジャンプの設定値が極端: FFUの抵抗値が大きすぎると圧力場が不安定になる。初期は小さい値から始めて段階的に増やす
- Under-Relaxation Factorの調整: 圧力を0.2、運動量を0.5に下げて試す
圧力補間がPRESTO!でないとダメなのはBoussinesq近似特有の問題なんですね。
2. DPMの粒子が壁に全部トラップされる
症状: 粒子がリリース直後にほぼ全て壁面に捕捉されて、評価点まで到達しない。
対策:
- 壁面のDPM境界条件がTrapになっていないか確認（ReflectまたはEscapeに変更）
- ブラウン力をONにしているか確認（微小粒子はブラウン運動で壁から離れる）
- 粒子の初速度が気流速度と一致しているか確認（0で放出すると重力落下する）
3. 温度分布が室内でほぼ均一
装置発熱を入れているのに温度差が出ないケースです。
対策:
- 離散化スキームを確認（First Order Upwindだと数値拡散で温度がなまる）
- メッシュが粗すぎないか確認（装置近傍で最低20mm以下のセル）
- 装置発熱量の単位を確認（W vs. W/m² の間違い）
4. 実測との乖離が大きい
実測の風速やパーティクルカウントとCFD結果が合わないときはどうすればいいですか？
確認ポイント:
| チェック項目 | よくある問題 | 対処 |
|---|---|---|
| FFU実風量 | カタログ値と実際の差 | 実測値でBC更新 |
| リーク | ダクト接続部からの漏れ | 漏れ風量を追加 |
| 装置排気量 | 実際の運転状態と異なる | 実測値確認 |
| 床下プレナム | モデル簡略化の影響 | 開口率・配管障害物の精度向上 |
| 人体モデル | 動的発塵を無視 | 非定常+DPMで再計算 |
FFUの実風量がカタログ値と違うことが多いんですか？
フィルタの目詰まり、ダクト圧損、FFUファンの経年劣化で実際の面風速は新品時から10〜20%低下していることがある。可能な限り実測値を使うべきだ。
Fluent特有のエラー
Fluent固有の注意点はありますか？
DPMのIncompleteは粒子が室内を延々と漂っている状態ですよね。最大ステップ数を増やすか、タイムアウトを設定するのがいいと。
そうだ。Physical Time Limitを設定すると、一定時間後に強制的に打ち切れるから、解析時間の予測がしやすくなる。クリーンルームの換気回数から滞留時間を見積もって、その3倍程度に設定するのが目安だ。
「FFUの圧力ばらつき」が引き起こす謎の不均一流れ
クリーンルームCFDのトラブルシューティングで頻繁に遭遇するのが、複数のFFU（ファンフィルタユニット）間での流量ばらつきです。各FFUは独立した電動ファンを持ちますが、隣接ユニット間の静圧バランスが崩れると一部のFFUが「逆流モード」に近い状態になることがあります。シミュレーション上では全FFUを均一流量条件で設定したのに実測値と合わない、という場合はこのばらつきを疑いましょう。各FFUの静圧—流量特性（P-Q曲線）をモデルに組み込んだ連成解析を行うと、現場の不均一流れを再現できることが多いです。
関連トピック
なった
詳しく
報告