クリーンルーム気流解析
理論と物理
概要
先生！ クリーンルームの気流解析って、半導体工場とかで使うやつですよね？ どういう物理が関わってくるんですか？
クリーンルーム気流解析は、室内の清浄度を維持するために一方向流（ユニフロー）や乱流置換方式の気流パターンを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でパーティクル追跡を行う際は、この分布特性を理解した上で粒径範囲を設定することが重要です。
各項の物理的意味
- 時間項 $\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を実際に解くとき、有限体積法ですよね？ 具体的な離散化スキームはどう選ぶんですか？
クリーンルーム気流は低マッハ数の非圧縮性流れだから、圧力ベースソルバーを使う。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）を組み合わせれば、どの作業位置でコンタミが発生しやすいかも定量評価できます。
風上差分（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の実務的な解析フローを教えてください。
全体の流れは、(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の解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ 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」という使い分けが多く、どちらか一方で完結させようとすると非効率になりがちです。
選定で最も重要な3つの問い
- 「何を解くか」：クリーンルーム気流解析に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
先端トピックと研究動向
クリーンルーム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曲線）をモデルに組み込んだ連成解析を行うと、現場の不均一流れを再現できることが多いです。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——クリーンルーム気流解析の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告