流動層シミュレーション
理論と物理
概要
先生、流動層シミュレーションって何をするものですか？
流動層は粒子充填層に下方からガスを吹き込み、粒子を浮遊・攪拌させる装置だ。石油精製のFCC（流動接触分解）、石炭ガス化、バイオマス燃焼、医薬品の造粒コーティングなど、化学工学の基幹技術だ。CFDで内部の粒子挙動やガス混合を予測する。
流動層のCFDにはどんな手法があるんですか？
大きく3つのアプローチがある。
| 手法 | 特徴 | 粒子数の目安 |
|---|---|---|
| Euler型粒体モデル（TFM） | 粒子を連続体として扱う | 無制限（粒子群） |
| DEM-CFD | 個々の粒子を追跡 | 〜$10^6$ 個 |
| CPFD法 | parcelで粒子群を代表 | $10^6$〜$10^{12}$相当 |
支配方程式
TFM（Two-Fluid Model）の方程式を教えてください。
気相と固相それぞれに連続の式と運動量方程式を解く。固相の応力にはKTGF（粒子運動論）を使う。
固相の圧力は Granular Temperature $\Theta_s$ から求まる。
流動化の鍵になるパラメータは何ですか？
最小流動化速度 $U_{mf}$ だ。これを超えるガス流速で粒子が浮遊し始める。Ergun式から推定できる。
流動化状態では圧力損失が層重量とバランスする。$\Delta p = (1-\varepsilon_{mf})(\rho_s - \rho_g) g L$ が流動化の判定基準だ。
Geldart分類
粒子の種類によって流動化の仕方が違うんですか？
Geldart（1973）の分類が基本だ。
| グループ | 粒径 | 流動化特性 | 例 |
|---|---|---|---|
| A | 20〜100 μm | 均一膨張後に気泡発生 | FCC触媒 |
| B | 100〜1000 μm | 直接気泡流動化 | 砂、ガラスビーズ |
| C | < 20 μm | 凝集性が強く流動化困難 | 小麦粉、タルク |
| D | > 1000 μm | 噴流（スパウト）形成 | 穀物、石炭塊 |
流動化の発見——FCCプロセスと流動層革命の前夜
流動層（Fluidized Bed）技術の産業応用が急拡大したのは1940年代で、Standard Oil（現ExxonMobil）が流動接触分解（FCC）プロセスを開発したことが発端です。砂粒が空気で浮揚して「液体のように」振る舞うこの現象は、20世紀前半の化学工学者には魔法に見えたと言われています。流動化の基礎理論であるErgun方程式（1952年）は今も流動層設計の根幹で、ε（ボイド率）とΔP（圧力損失）の関係を半経験的に結びつけます。CFDによる流動層シミュレーションは、このErgunモデルが粒子間力をどう表現するかに強く依存します。
各項の物理的意味
- 時間項 $\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の数値的なポイントを教えてください。
TFM（Euler型粒体モデル）での流動層シミュレーションには、いくつかの特有の課題がある。
メッシュとメソスケール構造
流動層では「クラスター」と呼ばれる粒子の密集構造が重要だ。クラスターのサイズは粒子径の10〜100倍程度で、これを解像するにはメッシュを十分に細かくする必要がある。
| メッシュ | 解像度 | 計算コスト | 精度 |
|---|---|---|---|
| Fine | $\Delta x \approx 5 d_p$ | 非常に高い | 高 |
| Standard | $\Delta x \approx 10$〜$20 d_p$ | 中程度 | 良好 |
| Coarse + Filter | $\Delta x > 50 d_p$ | 低い | フィルタモデル必要 |
粗いメッシュではクラスターが解像できないとどうなるんですか？
層膨張率を過大評価し、ガスバイパスを過小評価する。つまり実際よりも均一に流動化しているように見えてしまう。Filtered TFM（Igci et al., 2008; Ozel et al., 2013）で補正するか、十分に細かいメッシュを使う必要がある。
抗力モデルの選択
流動層で最も重要なクロージャモデルは気固間抗力だ。
| モデル | 特徴 | 推奨用途 |
|---|---|---|
| Gidaspow | Ergun + Wen-Yuの切替 | BFB（気泡流動層）標準 |
| Syamlal-O'Brien | 連続式、パラメータ調整可能 | 汎用 |
| EMMS | メソスケール構造を考慮 | CFB（循環流動層） |
| Koch-Hill | LBMデータベース | 高精度 |
EMMSモデルって何ですか？
Energy Minimization Multi-Scale（EMMS）モデルは、Li & Kwauk（中国科学院）が提案した手法で、クラスター構造によるガスバイパス効果を抗力に反映する。粗いメッシュでもある程度クラスターの影響を取り込めるため、工業スケールの循環流動層で広く使われている。
タイムステップと計算時間
流動層TFMは非定常計算が必須で、物理時間として数秒〜数十秒の計算が必要だ。
| パラメータ | 推奨値 | 備考 |
|---|---|---|
| $\Delta t$ | $10^{-4}$〜$10^{-3}$ s | Courant数 < 0.5 |
| 物理時間 | 5〜30 s | 統計的定常に達するまで |
| 平均化開始 | 2〜5 s後 | 初期過渡を除外 |
TFM vs DEM-CFD——流動層シミュレーションの二大潮流
流動層CFDには大別してTwo-Fluid Model（TFM/Euler-Euler）とDEM-CFD（Euler-Lagrange）の二つのアプローチがあります。TFMは粒子を連続体として扱うため100万以上の粒子系にスケールしやすい一方、個々の粒子接触は平均化されて失われます。DEM-CFDは個々の粒子を追跡するため物理的正確さでTFMを上回りますが、粒子数が10万を超えると計算コストが急増します。工業規模の流動層（直径3 m × 高さ10 m）のフルスケールCFDは2020年代でもTFMが現実的で、DEM-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次だが、流れの方向を正しく捕捉するため安定性が高い。
実践ガイド
実践ガイド
流動層シミュレーションの手順を教えてください。
Geldart B粒子のバブリング流動層（BFB）を例に説明しよう。
1. 形状: 円筒形リアクター（直径0.2 m、高さ1 m）、分散板
2. メッシュ: 六面体、セルサイズ 5 mm（粒径500 μmの10倍）
3. 初期充填: 下部0.3 mに$\alpha_s = 0.6$で固相を配置
4. ガス入口: 分散板から均一速度（$U/U_{mf} = 3$〜5）
5. 出口: 上部を圧力出口
6. KTGFパラメータ: Gidaspow粘性、Lun固相圧力、反発係数0.9
7. 非定常計算: $\Delta t = 5 \times 10^{-4}$ s で10秒計算
圧力損失の検証
まず何を検証すればいいですか？
流動層のCFDで最初に確認すべきは圧力損失だ。流動化状態では理論的に次が成り立つ。
CFDで計算した圧力損失がこの理論値と一致すれば、質量保存と力のバランスが正しいことの確認になる。
気泡径と気泡速度
気泡の挙動はどう評価するんですか？
Darton et al.（1977）の気泡径相関やDavidson & Harrison（1963）の気泡上昇速度理論と比較する。
CFDでの気泡は体積分率のiso-surface（$\alpha_g = 0.8$等）で抽出し、等価直径と上昇速度を計測する。2D計算と3D計算で気泡径が異なるので注意が必要だ。
よくある失敗
| 症状 | 原因 | 対策 |
|---|---|---|
| 流動化しない | ガス流速 < $U_{mf}$ | 速度を上げる |
| 過剰膨張 | 抗力モデルが不適切 | EMMSモデル or メッシュ細分化 |
| 固相が壁面に固着 | 壁面境界条件 | no-slip→free-slipに変更 |
| 非対称パターン | 2Dの限界 | 3D計算に切り替え |
石炭流動層燃焼炉——脱炭素時代の移行技術CFD
循環流動層（CFB）ボイラーは、石炭以外にもバイオマス・廃棄物固形燃料（SRF）を混焼できる柔軟性から、脱炭素移行期の重要技術として再評価されています。ライザー内の粒子循環速度・燃焼温度分布・SO2/NOx生成をCFDで予測することが新燃料の導入認証に直結します。三菱パワーや川崎重工がCFB炉のTFM-CFDを使った設計最適化事例を公開しており、ライザー断面形状の変更だけで燃焼効率を2%向上させた事例があります。CFDなしでは数百時間の実炉試験が必要な検証作業が、シミュレーションで1/10以下に短縮されています。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ CFDでも入口条件のミスは下流全体に波及します。
ソフトウェア比較
商用ツール比較
流動層シミュレーションに使えるツールを教えてください。
| ツール | 手法 | 抗力モデル | 反応 | 特徴 |
|---|---|---|---|---|
| Ansys Fluent | TFM (KTGF) | Gidaspow, Syamlal, EMMS(UDF) | Heterogeneous | 最も多機能 |
| STAR-CCM+ | TFM | 標準モデル | 対応 | ポリヘドラルメッシュ |
| Barracuda VR | CPFD | 独自モデル | 詳細反応 | 流動層専用、GPU対応 |
| OpenFOAM | TFM / DEM-CFD | Gidaspow, Koch-Hill | 基本対応 | 完全OSS |
| MFIX | TFM / DEM | 豊富 | 詳細反応 | DOE開発OSS、流動層特化 |
MFIXって何ですか？
米国エネルギー省NETLが開発したオープンソースの多相流コードで、流動層に特化している。TFMとDEMの両方に対応し、反応モデルも充実している。化学ルーピング燃焼やガス化の研究でよく使われるが、GUIの操作性は商用ツールに劣る。
用途別推奨
| 用途 | 推奨ツール | 理由 |
|---|---|---|
| FCC再生塔 | Barracuda VR | 大規模+反応、業界標準 |
| CFB（循環流動層） | Fluent + EMMS | EMMS抗力の有効性 |
| バイオマスガス化 | MFIX, Fluent | 反応モデルの充実 |
| 医薬品造粒 | Fluent + DEM | DEM-CFDで粒子レベルの解析 |
| 基礎研究 | OpenFOAM, MFIX | コード改変の自由度 |
Barracuda VRは流動層では定番なんですね。
石油メジャーやエンジニアリング会社で広く使われている。CPFD法はparcelで粒子群を代表するため、数十億粒子相当の工業スケールを現実的な計算時間で解ける。最近GPU対応も進み、さらに高速化されている。
MFiX vs Barracuda VR——流動層専用ツールの実力比較
流動層シミュレーションに特化したツールとして、DOEが開発したオープンソースのMFiXと商用のBarracuda VR（CPFD Software）が双璧をなします。MFiXはTFMとDEMの両方を実装し、学術研究の標準プラットフォームとして世界中の大学で使われています。Barracuda VRはMP-PIC（Multiphase Particle-in-Cell）法という独自手法で数千万粒子系をGPUで実用的な時間内に計算でき、石油・ガス産業の大型FCC装置設計で採用実績があります。商用ツールとオープンソースの差は「サポートの有無」と「大規模系の計算速度」に集約されます。
選定で最も重要な3つの問い
- 「何を解くか」：流動層シミュレーションに必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
先端技術と研究動向
流動層シミュレーションの最新研究にはどんなものがありますか？
いくつかの重要な方向性がある。
Filtered TFM
粗いメッシュでの精度劣化を補正するために、LESのサブグリッドモデルに相当するフィルタリング手法が研究されている。Igci & Sundaresan（2011）のフィルタリング抗力・応力モデルが代表的だ。
フィルタリングすると何が変わるんですか？
メソスケール構造（クラスター）によるガスバイパス効果を補正抗力として取り入れる。具体的には、細かいメッシュで計算したデータからフィルター幅依存のモデルを構築し、粗いメッシュに適用する。計算コストを1〜2桁削減できる可能性がある。
Coarse-Grained DEM
工業スケールの流動層をDEM-CFDで解くために、代表粒子（coarse particle）で実粒子群を表現する手法が進んでいる。Sakai（2016）のCG-DEM、Lu et al.のMP-PICなどがある。
化学ルーピング燃焼（CLC）のCFD
CLCってCO2回収に関係するんですか？
そうだ。金属酸化物粒子（キャリア）と燃料を流動層で反応させ、CO2を分離回収する次世代燃焼技術だ。2つの流動層（Air Reactor, Fuel Reactor）間で粒子を循環させるため、CFDでの粒子輸送と反応の連成が不可欠だ。EUのSUCCESSプロジェクト等で大規模CFDが行われている。
機械学習によるスケールアップ
流動層のスケールアップは伝統的に経験則に依存してきたが、CFDデータを教師データとして機械学習でスケーリング則を構築する研究が増えている。ラボスケールのCFD結果からパイロットスケールの挙動を予測するサロゲートモデルが期待されている。
Coffee Break よもやま話
化学ループ燃焼——酸素キャリア粒子のCFD最前線
化学ループ燃焼（CLC: Chemical Looping Combustion）は、金属酸化物粒子（Fe2O3、CuO等）を酸素キャリアとして燃料と空気を完全に分離することで、CO2を高濃度で回収できる次世代燃焼技術です。空気炉・燃料炉を循環する粒子の反応・酸化・還元をCFDで追跡するには、TFMに化学反応ソルバーと粒子物性の変化を組み込む必要があります。ダルムシュタット工科大学とチャルマース大学のグループが独立して開発したCLC-CFDモデルは、粒子転化率の実験値を10%以内で再現しており、CLCプロセスの商業化設計で参照されています。
工業スケールの流動層をDEM-CFDで解くために、代表粒子（coarse particle）で実粒子群を表現する手法が進んでいる。Sakai（2016）のCG-DEM、Lu et al.のMP-PICなどがある。
CLCってCO2回収に関係するんですか？
そうだ。金属酸化物粒子（キャリア）と燃料を流動層で反応させ、CO2を分離回収する次世代燃焼技術だ。2つの流動層（Air Reactor, Fuel Reactor）間で粒子を循環させるため、CFDでの粒子輸送と反応の連成が不可欠だ。EUのSUCCESSプロジェクト等で大規模CFDが行われている。
流動層のスケールアップは伝統的に経験則に依存してきたが、CFDデータを教師データとして機械学習でスケーリング則を構築する研究が増えている。ラボスケールのCFD結果からパイロットスケールの挙動を予測するサロゲートモデルが期待されている。
化学ループ燃焼——酸素キャリア粒子のCFD最前線
化学ループ燃焼（CLC: Chemical Looping Combustion）は、金属酸化物粒子（Fe2O3、CuO等）を酸素キャリアとして燃料と空気を完全に分離することで、CO2を高濃度で回収できる次世代燃焼技術です。空気炉・燃料炉を循環する粒子の反応・酸化・還元をCFDで追跡するには、TFMに化学反応ソルバーと粒子物性の変化を組み込む必要があります。ダルムシュタット工科大学とチャルマース大学のグループが独立して開発したCLC-CFDモデルは、粒子転化率の実験値を10%以内で再現しており、CLCプロセスの商業化設計で参照されています。
トラブルシューティング
トラブルシューティング
流動層シミュレーションでよくあるトラブルを教えてください。
順番に見ていこう。
1. 流動化せず固着する
症状: ガスを流しても粒子層が浮き上がらない。
対策:
- ガス流速が $U_{mf}$ を超えていることを確認（Ergun式で計算）
- 分散板の境界条件が velocity-inlet になっていることを確認
- 初期固相体積分率が高すぎないか確認（$\alpha_s = 0.55$〜$0.60$）
- 摩擦モデルの設定を確認（Frictional viscosity が過大だと動かない）
2. 層が過剰に膨張する
粒子が全部飛んでいってしまいます…
対策:
- メッシュが粗すぎないか確認（クラスター未解像→抗力過小→過膨張）
- 抗力モデルをGidaspowからEMMSに変更
- フリーボード領域を十分に確保
- 出口境界で固相体積分率のbackflow値を0に設定
3. 圧力損失が理論値と合わない
対策:
- 壁面境界条件を確認（固相のno-slip vs. free-slip）
- 粒子密度と粒子径が正しいか確認
- 重力の方向と大きさが正しいか確認
- 計算が定常に達してから測定しているか確認
4. 非対称な流動パターン
症状: 本来対称であるべき系で非対称なパターンが出現。
2D計算の場合は本質的な限界であり、3D計算に切り替える。3D計算でも発生する場合は、入口条件の微小な非対称性が原因の可能性があるので、入口プロファイルを確認する。
5. ツール固有の注意点
| ツール | 注意点 |
|---|---|
| Fluent | Granular Temperature方程式のPDE/Algebraic切替が結果に影響大 |
| STAR-CCM+ | 固相の壁面境界条件タイプの選択（partial-slip推奨） |
| Barracuda VR | Close-pack volume fractionとsmall number limitの設定 |
| MFIX | 旧バージョンと新バージョンで入力ファイル形式が異なるので注意 |
流動化しない——最小流動化速度の計算ミス対処法
流動層CFDで最初につまずくのが「最小流動化速度（Umf）の設定ミス」です。入口ガス速度がUmfを下回ると固定層として振る舞い、粒子が全く流動しないシミュレーションになります。Ergun方程式から理論的にUmfを計算できますが、粒子径分布・形状因子（sphericity）が実際と異なると計算値が2倍以上ずれることがあります。確認手順として「圧力損失 vs 速度」の掃引計算（少なくとも5点）でUmfをまず数値的に特定し、理論値との整合を取ることを強く推奨します。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——流動層シミュレーションの問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告