撹拌槽CFD
撹拌槽CFDの理論基礎
概要
先生！ 撹拌槽のCFD解析ってどういう場面で使うんですか？
化学プラント、医薬品製造、食品加工、水処理などで使われる撹拌槽の流動パターン、混合時間、動力消費をCFDで予測する技術だ。インペラ（撹拌翼）の回転が作る複雑な3次元流れをNavier-Stokes方程式で解く。
支配方程式
撹拌槽の基本式を教えてください。
まず無次元数が重要だ。撹拌レイノルズ数とパワー数が基本になる。
$N$ が回転数 [rps]、$D$ がインペラ直径、$P$ が撹拌動力ですね。
そう。$Re_{imp} > 10^4$ で完全乱流、$Re_{imp} < 10$ で層流と判断する。遷移領域（10〜10,000）は最も解析が難しい。
混合時間 $\theta_m$ はトレーサー応答から定義する。
$T$ が槽径ですね。$\theta_m N$ が無次元混合時間で、完全乱流では定数（インペラ形状に依存）になると。
その通り。6枚平板タービン（Rushton Turbine）では $\theta_m N \approx 30$〜$50$ が典型値だ。
インペラの分類
| インペラ形状 | Np（乱流域） | 流れパターン | 用途 |
|---|---|---|---|
| Rushton Turbine (6枚平板) | 5.0〜5.5 | ラジアル | 気液混合、一般反応 |
| Pitched Blade Turbine (45°) | 1.2〜1.7 | アキシャル-ラジアル | 固液懸濁、混合 |
| Hydrofoil (A310, A320) | 0.3〜0.4 | アキシャル | 低剪断混合 |
| Anchor | 0.4〜0.8 (層流) | タンジェンシャル | 高粘度流体 |
| Helical Ribbon | 0.5〜1.0 (層流) | アキシャル+タンジェンシャル | 超高粘度 |
パワー数がインペラ形状で全然違いますね。Rushtonは5以上で、Hydrofoilは0.3程度ですか。
Rushtonは強い剪断場を作るから気液分散に適しているが、動力消費が大きい。Hydrofoilは軸流で効率的に液を循環させるが、気液分散能力は低い。用途に応じて使い分ける。
実務上の注意点
攪拌工学の父・Rushton——Rushtonタービンと無次元動力数の確立（1950年）
攪拌槽工学の礎を築いたのは米国のJ. H. Rushton（ラッシュトン）だ。1950年の論文「Power Characteristics of Mixing Impellers」で、インペラの無次元動力数Np = P/(ρN³D⁵)を定義し、Npが高Re域で一定値（完全乱流域でディスクタービン型は約5）に収束することを実験で証明した。この「ラッシュトンタービン」と動力数相関は以後70年にわたり攪拌設計のデファクトスタンダードとなった。現代のCFDでは彼の実験をベンチマークとして乱流モデルの検証が行われており、標準k-εモデルのNp予測誤差が10〜15%程度あることが多数の検証論文で確認されている。CFD精度の限界を知るために古典実験データの価値は今も変わらない。
撹拌槽CFDの数値計算手法
数値手法の詳細
撹拌槽でインペラが回転する流れをどうやって解くんですか？
インペラの回転をCFDでモデル化する方法は主に3つある。
回転モデルの選択
| 手法 | 概要 | 計算コスト | 精度 |
|---|---|---|---|
| MRF (Multiple Reference Frame) | 回転領域を定常的に扱う | 低（定常） | 中 |
| Sliding Mesh (SM) | 回転領域のメッシュを実際に回転 | 高（非定常） | 高 |
| Overset Mesh | 重合メッシュで回転 | 高（非定常） | 高 |
MRFとSliding Meshの使い分けはどうしますか？
MRFは定常解を得る手法で、時間平均的な流動パターンやパワー数の予測に使う。Sliding Meshは非定常解で、インペラとバッフルの干渉による周期的な力の変動（トルク変動）や、混合時間のトレーサー追跡に必要だ。
実務的にはまずMRFで大まかな流れ場を確認し、その後Sliding Meshで精密評価するのが効率的だ。
MRFの設定
FluentでMRFを設定する手順を教えてください。
1. メッシュでインペラ周辺に円筒形の回転ゾーンを作成
2. Cell Zone Conditions → 回転ゾーンにFrame Motion → Rotational Velocityを設定
3. 回転ゾーンの上下面はInterfaceで外部ゾーンと接続
4. 回転ゾーンにバッフルを含めないこと（バッフルは静止側）
回転ゾーンの寸法目安:
- 直径: インペラ直径の1.1〜1.3倍
- 高さ: インペラ高さの1.5〜2.0倍
- インペラとゾーン境界の距離: インペラ径の5〜15%
回転ゾーンの境界がインペラに近すぎるとどうなりますか？
インペラが生成する後流（ウェイク）が回転ゾーンの境界で不自然に切断され、パワー数やポンプ流量の予測精度が低下する。十分な余裕を持たせること。
メッシュ戦略
撹拌槽のメッシュで重要なポイント:
| 領域 | メッシュサイズ | 備考 |
|---|---|---|
| インペラ翼面 | D/100〜D/50 | 翼上面/下面で圧力差を解像 |
| インペラ翼端 | D/100 | 渦の発生点 |
| バッフル周辺 | T/100 | バッフル後方の渦 |
| 槽壁近傍 | T/50〜T/20 | 壁面境界層 |
| 液面付近 | 自由表面解析の場合細分化 | VOF使用時 |
全体のセル数の目安は？
標準的な単段インペラ+4枚バッフルの撹拌槽で100万〜500万セルが目安だ。Sliding Meshで長時間の混合シミュレーションを行う場合、計算時間はインペラ数十回転分（数百〜数千タイムステップ）が必要。
乱流モデル
完全乱流（$Re_{imp} > 10^4$）ではRealizable k-epsilon + Standard Wall Functionが撹拌槽の定番だ。Npの予測精度が高いことが多くの文献で検証されている。
ただしSST k-omegaの方がインペラ後流の渦構造をよく捉える場合があり、混合時間の予測にはSST k-omegaが良い結果を出すことがある。LESは研究用で、渦構造の詳細な可視化に使われる。
攪拌槽CFDのMRF法——インペラ回転の数値的扱いとその精度限界
攪拌槽CFDで最もよく使われる「MRF法（Multiple Reference Frame法）」は、インペラ周辺領域を回転座標系で解き、タンク本体は静止座標系で解くアプローチだ。定常計算が可能で計算が速い反面、インペラ-バッフル間の非定常干渉（Impeller-Baffle Interaction）を捉えられず、バッフル直後の局所流れ予測精度が下がる。より高精度な「スライディングメッシュ（SM）法」は回転領域と静止領域をリアルタイムで連結し非定常計算を行うため、MRFより精度は高いが計算コストが5〜10倍かかる。実務判断の目安は「バッフル周辺の詳細流れ・混合時間・ガス分散挙動 → SM法」「流量・圧力・全体的な流れパターン → MRF法」という使い分けだ。
撹拌槽CFDの実務適用
実践ガイド
撹拌槽CFDの実務的な評価手順を教えてください。
パワー数の検証
まずCFD結果の妥当性を確認するため、パワー数を実験相関式と比較する。
インペラのトルク $M$ から撹拌動力を計算:
Fluentではインペラ翼面でのMoment（モーメント）をReportから取得できる。
| インペラ | 文献のNp値 | CFD予測値（目安） | 許容誤差 |
|---|---|---|---|
| Rushton 6DT | 5.0〜5.5 | 4.5〜5.5 | ±10% |
| PBT 45° 4枚 | 1.2〜1.7 | 1.1〜1.8 | ±15% |
| A310 | 0.30〜0.35 | 0.28〜0.38 | ±15% |
±10〜15%の範囲に入っていれば妥当と考えていいですか？
MRFの定常計算では±15%程度の誤差は許容範囲だ。Sliding Meshで時間平均すると±10%以内に収まることが多い。
混合時間の評価
混合時間をCFDで求める方法を教えてください。
Sliding Meshの非定常計算でトレーサー（パッシブスカラー）を追跡する。
手順:
1. MRFで定常流れ場を収束させる
2. Sliding Meshに切り替え
3. 液面近くの特定セルにトレーサーを初期配置（質量分率 = 1.0）
4. Species Transportで拡散+対流によるトレーサーの混合を追跡
5. 槽内の複数モニタリング点でトレーサー濃度が最終値の±5%以内に到達する時間を $\theta_m$ とする
±5%基準で混合時間を定義するのが一般的ですよね。±2%や±1%もありますか？
ある。基準を厳しくすると混合時間は長くなる。±5%が最も一般的で、ISO 10932でも採用されている。
固液懸濁のモデリング
固体粒子の懸濁状態をCFDで評価できますか？
Eulerian-Eulerian多相流モデルで固体相と液相を両方連続体として解く。Zwietering (1958) の完全浮遊回転数 $N_{js}$ との比較で妥当性を検証する。
$S$ はインペラ形状定数、$X$ は固体濃度 [wt%] ですね。
よくある失敗と対策
| 失敗パターン | 原因 | 対策 |
|---|---|---|
| Npが文献値の半分以下 | 回転ゾーンの境界がインペラに近すぎる | ゾーン径をD×1.2以上に拡大 |
| 流れパターンが非対称 | メッシュの非対称性 | 周期モデル（90°セクター）で対称性を確保 |
| 混合時間が長すぎる | MRFの定常解で評価 | Sliding Meshの非定常解に切り替え |
| 液面が変形して発散 | VOFの時間ステップが大きすぎる | CFL < 0.5で時間ステップを制限 |
| トレーサーが拡散しない | Species Transportの拡散係数が小さすぎる | 乱流拡散が支配的なので乱流Sc数を確認 |
医薬品製造の攪拌スケールアップ——CFDが解いたバッチ再現性問題
医薬品製造における攪拌槽のスケールアップは、研究室スケール（1L）から製造スケール（10,000L）への拡大で混合時間・せん断力・泡立ちなどが大きく変化するため、再現性確保が難関だ。従来は「P/V（単位体積当たり動力）一定」のスケールアップ則を使っていたが、この則では槽スケールが大きくなるほど混合時間が長くなり、ゲル化や副反応が増える問題があった。CFD（インペラ回転をMRF法でモデル化）で槽内の混合時間分布を予測し、デッドゾーンを特定して槽形状（バッフル数・位置）を最適化することで、ある医薬品メーカーでは製造スケールでの収率を8%向上させた事例が製薬学会誌で報告されている。
撹拌槽CFDのソフトウェア比較
商用ツール比較
撹拌槽CFDに適したツールはどれですか？
| ツール | MRF | Sliding Mesh | 多相流 | 撹拌特化機能 |
|---|---|---|---|---|
| Ansys Fluent | あり | あり | Eulerian/VOF/DPM | MixSim（旧プラグイン） |
| STAR-CCM+ | あり | あり (Rigid Body Motion) | Eulerian/VOF | Macro Mixing Model |
| Ansys CFX | あり | あり | 多相Eulerian | ターボ機械ワークフロー |
| OpenFOAM | あり (MRF) | あり (AMI) | Eulerian/VOF | なし（手動設定） |
| COMSOL | あり | あり (Moving Mesh) | Phase Field/Level Set | Mixer Moduleあり |
| M-Star CFD | - | LBM (GPU) | VOF/DPM | 撹拌専用ツール |
M-Star CFDって何ですか？
M-Star CFDは撹拌槽に特化したGPUベースのLBM（Lattice Boltzmann Method）ソルバーだ。GPU上でLESを実行し、従来のRANSベースのCFDと比較して10〜100倍の高速化を実現している。撹拌槽のデモ解析が数分〜数十分で完了する。
LBMで撹拌槽を解くメリットは何ですか？
用途別推奨ツール
| 用途 | 推奨 | 理由 |
|---|---|---|
| パワー数・混合時間の予測 | Fluent / STAR-CCM+ | MRF + Sliding Meshの実績 |
| 気液混合（通気撹拌） | Fluent (Eulerian) | 気泡径分布モデルが充実 |
| 非ニュートン流体の混合 | Fluent / COMSOL | レオロジーモデルが豊富 |
| 迅速なLES評価 | M-Star CFD | GPU LESで高速 |
| スケールアップ検討 | Fluent + Workbench DOE | パラメトリック解析 |
| 研究・教育 | OpenFOAM | 無償、コミュニティ充実 |
OpenFOAMでの撹拌槽設定
OpenFOAMで撹拌槽を解く場合のソルバーは？
MRFの場合はsimpleFoam + MRF、Sliding Meshの場合はpimpleFoam + AMI（Arbitrary Mesh Interface）を使う。
MRF設定（constant/MRFProperties）:
- cellZone: rotatingZone
- origin: (0 0 0)
- axis: (0 0 1)（Z軸回転の場合）
- omega: 6.2832（1 rps = 2π rad/s）
omegaの単位がrad/sなのを忘れがちですね。
そう。rpmで指定する場合は変換が必要だ（300 rpm = 300/60 × 2π ≈ 31.4 rad/s）。
攪拌槽CFDのツール選択——ANSYS FluentのMRF実装とRUSHTONインペラの精度比較
攪拌槽CFDの商用ツールとしてはANSYS Fluent、STAR-CCM+、OpenFOAMが三大選択肢だ。FluentはMRF設定がGUI上で比較的容易でRushtonタービンの動力数Np予測精度の検証論文が豊富（Np誤差±10〜15%）。STAR-CCM+はポリヘドラルメッシュの自動生成でインペラ周辺の複雑形状メッシュ作成の手間が少なく、動的メッシュ（スライディング）機能も成熟している。OpenFOAMはMRFmixerFoamとrotatingWallVelocity境界条件の組み合わせが標準だが、チュートリアルが少なく初期設定に専門知識が必要。製薬・食品業界では規制当局（FDA）への検証データ提出が必要なため、商用ツールの実績・サポート体制が選定の重要因子となる。
撹拌槽CFDの先端研究
先端トピックと研究動向
撹拌槽CFDの最新トレンドを教えてください。
1. GPU LESの実用化
GPU（NVIDIA A100, H100）上でのLBM-LESが撹拌槽解析で実用化されつつある。M-Star CFDに加えて、Inducta/SimuliaPowerFLOWもGPU LBMを提供している。
計算時間の比較（標準的なRushton撹拌槽、100万セル、10回転分）:
| 手法 | ハードウェア | 計算時間 |
|---|---|---|
| RANS (MRF) | CPU 16コア | 30分 |
| RANS (Sliding Mesh) | CPU 16コア | 4〜8時間 |
| LES (FVM) | CPU 128コア | 1〜3日 |
| LES (LBM) | GPU 1枚 (A100) | 30分〜2時間 |
GPU LBMならLESがRANS Sliding Meshと同程度の時間で計算できるんですね。
2. Population Balance Model (PBM)
気液撹拌では気泡径分布が反応速度に影響する。PBMで気泡の合一と分裂を追跡し、局所的な気泡径分布を予測する。
$n(d,t)$ が気泡の数密度関数で、Breakup（分裂）とCoalescence（合一）のソース/シンク項があるんですね。
FluentにはPBM（Method of Moments、Discrete Method、QMOM）が実装されている。STAR-CCM+にもS-gamma modelがある。
3. 反応を伴う撹拌槽のCFD
CFDと化学反応動力学を連成させて、反応転化率の空間分布を予測する。特に急速反応（ミクロ混合が律速）の場合、混合のタイムスケールと反応のタイムスケールの比（Damkohler数）が重要になる。
Da >> 1 だと混合が律速で、Da << 1 だと反応が律速ですね。
そう。Da >> 1の場合、RANS平均濃度だけでは反応率を正確に予測できない。乱流-反応の相互作用をモデル化するために、PDF（Probability Density Function）法やEDC（Eddy Dissipation Concept）が使われる。
4. 機械学習によるスケールアップ予測
スケールアップは化学工学の永遠のテーマですからね。CFD+MLで解決の糸口が見えてきたと。
完全な予測はまだ難しいが、CFDで得られた局所的な情報（剪断速度分布、エネルギー散逸率分布）を使えば、従来の無次元相関式よりも精度の高いスケールアップが可能になりつつある。
攪拌槽のマイクロミキシング——分子スケール混合がCFDの次のフロンティア
攪拌槽のCFD解析はメカロスケール（インペラ周辺の流れ）の再現には成熟しているが、「マイクロミキシング（Micro-mixing）」と呼ばれる分子レベルの均一化はまだCFDの大きな課題だ。マイクロミキシング時間が反応時間より長い系（速い化学反応）では、CFDの濃度分布が平均値を与えるだけで局所的な反応率を過大評価する。IEM（Interaction by Exchange with the Mean）モデルやCFD-Population Balance Modelの組み合わせが研究フロンティアで、製薬業界の結晶化制御やナノ粒子合成プロセスでの実用化が進んでいる。マイクロミキシングの定量化には実験では蛍光消光法（Iodide-Iodate法）との比較が標準だ。
撹拌槽CFDのトラブル対応
トラブルシューティング
撹拌槽CFDでよくある問題を教えてください。
1. パワー数が文献値と合わない
チェックポイント:
- MRF回転ゾーンの境界位置（インペラに近すぎると過小評価）
- インペラ翼面のメッシュ解像度（翼面に少なくとも20セル）
- バッフルが回転ゾーンの外側にあるか確認（バッフルをMRFゾーンに含めるとNpが激減する）
- 壁面処理（Enhanced Wall Treatmentの場合y+ ≒ 1を確認）
- Moment参照点が回転軸と一致しているか
バッフルが回転ゾーンに入ってしまうのは初歩的ですが致命的なミスですね。
そう。バッフルが回転座標系で解かれると、相対速度がゼロに近くなり、バッフルの抵抗効果が消えてしまう。
2. MRFで非物理的なジェットが発生
症状: MRFゾーンの境界面で不自然な速度の不連続が見られる。
対策:
- MRFゾーンの境界面を、流れの大きな勾配がある場所から離す
- 境界面のメッシュサイズを内外で一致させる（大きな差があると補間誤差が大きくなる）
- Sliding Meshに切り替えて検証する
3. Sliding Meshの計算が発散
非定常計算が不安定になるケースですね。
対策:
- 時間ステップを小さくする（インペラが1ステップで回転する角度が1°〜5°以下）
- MRFの定常解を初期条件として使用（ゼロからスタートしない）
- AMI/Interface面のメッシュサイズを両側で一致させる
- PIMPLEアルゴリズム（OpenFOAM）またはPISO（Fluent）を使用
時間ステップの目安:
1 rps（60 rpm）の場合、$\Delta t < 2/(360 \times 1) \approx 0.0056$ s ですね。
4. 混合時間の予測が実験の2倍以上
チェックポイント:
- 乱流拡散が適切に計算されているか（乱流Schmidt数 $Sc_t$ のデフォルト値0.7は撹拌槽で妥当）
- トレーサーの初期配置が実験と同じ位置か
- 混合時間の判定基準が実験と同じか（±5% vs. ±2%）
- MRFの定常解で混合時間を評価していないか（Sliding Meshが必須）
5. 自由表面が過度に変形して発散
対策:
- VOFのCourant数を0.25以下に制限
- Geo-Reconstruct（Sharp Interface）の代わりにCICSAMを使う
- 液面変形が小さい場合はフラットフリーサーフェス（Symmetry BC）で近似
- 渦中心の空気巻き込み（Vortex Ingestion）が起きる条件では、VOF+LESが必要
自由表面を解かずにSymmetryで近似できるケースもあるんですね。
フルード数 $Fr = N^2 D / g$ が十分小さい（< 0.1程度）場合は液面変形が小さいのでSymmetry近似が妥当だ。高速撹拌でFr > 0.3になると渦による液面低下が無視できなくなる。
攪拌槽CFDの「渦の深さが合わない」——液面自由表面モデルの選択
攪拌槽CFDで「CFDでは渦が浅い／深すぎる」という実験との不一致は、自由表面の扱い方に原因があることが多い。多くの攪拌槽CFD解析では、液面を固定した「Flat Surface近似」（自由表面を剛壁として扱う）を用いるが、これでは液面変形（渦形成）を捉えられない。VOF（Volume of Fluid）法で液面を追跡すれば渦形状は再現できるが、計算コストが3〜5倍になる。実務的な判断基準は「Froude数（Fr = N²D/g）が0.1を超えたら液面変形を無視できない」——高速攪拌（Fr>0.1）の場合はVOF法が必須で、低速攪拌（Fr<0.1）ならFlat Surface近似で十分だ。設定前にFrを計算することで、過剰な計算コストを避けられる。
関連トピック
なった
詳しく
報告