DEM-CFD連成
理論と物理
概要
先生、DEM-CFD連成って何ですか？粒子と流体を一緒に計算するんですか？
その通り。DEM（Discrete Element Method：離散要素法）で個々の粒子の運動を追跡し、CFDで流体場を解く。両者を双方向にカップリングすることで、粒子-流体間の相互作用を再現する手法だ。
Lagrangian粒子追跡（DPM）とは何が違うんですか？
決定的な違いは粒子間接触の扱いだ。DPMでは粒子同士の衝突は簡易モデル（確率的な衝突）か無視だが、DEMでは粒子間の接触力を弾性ばね-ダッシュポットモデルで厳密に計算する。だから粉体・顆粒のように粒子間力が重要な系に向いている。
支配方程式
DEM側の方程式を教えてください。
各粒子 $i$ の並進運動と回転運動を追跡する。
$\mathbf{F}_{c,ij}$ は粒子 $j$ との接触力で、Hertz-Mindlinモデルが代表的だ。法線方向の接触力は次のように表される。
ここで $E^$ は等価ヤング率、$R^$ は等価半径、$\delta_n$ は重なり量、$\gamma_n$ は減衰係数だ。
CFD側はどうなりますか？
局所平均化したNavier-Stokes方程式を解く。粒子の存在によるボイド率 $\varepsilon_f$ を考慮する。
$\mathbf{S}_p$ は粒子から流体への反力（運動量ソース項）で、CFDセル内の全粒子に作用する流体力の和を体積で割ったものだ。
抗力モデル
粒子に作用する流体力はどうモデル化するんですか？
最も重要なのは抗力で、局所ボイド率に依存する抗力モデルを使う。
粒子に作用する流体力はどうモデル化するんですか？
最も重要なのは抗力で、局所ボイド率に依存する抗力モデルを使う。
| モデル | 適用範囲 | 特徴 |
|---|---|---|
| Ergun | $\varepsilon_f < 0.8$ | 充填層向け |
| Wen-Yu | $\varepsilon_f > 0.8$ | 希薄領域向け |
| Gidaspow | 全領域 | Ergun + Wen-Yuの切り替え |
| Di Felice | 全領域 | 連続的な遷移 |
| Koch-Hill | 全領域 | 格子ボルツマン法データベース |
抗力以外にも力はあるんですか？
圧力勾配力、仮想質量力、Saffman揚力、Magnus力なども考慮できるが、密度比 $\rho_p / \rho_f \gg 1$（粉体-空気系など）では抗力が支配的で、他の力は省略可能な場合が多い。
Cundallの革命——1979年、粒子は「ぶつかる」と定式化された
DEM（離散要素法）の創始者Peter Cundallは、岩盤の破壊力学を研究する中で「各粒子の接触力をバネ-ダッシュポット系でモデル化する」という発想を1979年に論文化しました。当初は岩盤ブロックの崩壊解析が目的でしたが、30年後にCFDと組み合わせたDEM-CFDが流動層・ミキサー・錠剤コーティング機の設計ツールとして製薬・化学業界に急速普及します。Cundall自身は「こんなに広く使われるとは思わなかった」と後年語ったとされ、1つのシンプルなモデルが業界を変えた好例です。
各項の物理的意味
- 時間項 $\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$。時間刻みの安定性に直結 |
数値解法と実装
数値解法の詳細
DEMとCFDの連成はどうやって実現するんですか？
基本的な連成スキームを説明しよう。
1. CFDタイムステップ $\Delta t_{CFD}$ で流体場を解く
2. 各粒子位置での流体速度・圧力を補間
3. 流体力（抗力等）を計算して各粒子に適用
4. DEMタイムステップ $\Delta t_{DEM}$ で粒子を更新（複数サブステップ）
5. 粒子位置からボイド率を再計算
6. 粒子→流体の反力をCFDソース項に反映
7. 次のCFDステップへ
DEMとCFDのタイムステップが違うんですか？
DEMのタイムステップは接触力計算のために非常に小さく、Rayleigh時間の20〜30%が目安だ。
典型的にはCFDの1ステップに対してDEMを100〜1000サブステップ実行する。これがDEM-CFDの計算コストの主因だ。
ボイド率の計算
ボイド率ってどうやって計算するんですか？
CFDセルに含まれる粒子の体積から計算する。粒子がセル境界をまたぐ場合の処理には以下のアプローチがある。
| 手法 | 概要 | 特徴 |
|---|---|---|
| Cell Centre | 粒子中心があるセルに全体積を帰属 | 簡単だが不連続 |
| Divided Volume | 粒子体積をセル間で分配 | より滑らか |
| Diffusion-based | カーネル関数で平滑化 | 最も滑らかだが計算コスト大 |
CFDセルサイズは粒子径の3〜5倍以上であることがunresolved DEM-CFDの前提条件だ。セルが粒子より小さいとボイド率の定義が破綻する。
resolved vs. unresolved DEM-CFD
resolvedとunresolvedの違いは何ですか？
unresolvedは粒子径がCFDメッシュより小さく、抗力モデルで相互作用を表現する。resolvedは粒子径がメッシュより大きく、粒子表面のflow fieldを直接解像する。Immersed Boundary法やOverset Meshで実現するが、粒子数は数百個程度に限られる。
主要ソフトウェア
DEM-CFD連成ができるツールにはどんなものがありますか？
resolvedとunresolvedの違いは何ですか？
unresolvedは粒子径がCFDメッシュより小さく、抗力モデルで相互作用を表現する。resolvedは粒子径がメッシュより大きく、粒子表面のflow fieldを直接解像する。Immersed Boundary法やOverset Meshで実現するが、粒子数は数百個程度に限られる。
DEM-CFD連成ができるツールにはどんなものがありますか？
| DEM側 | CFD側 | 連成方法 |
|---|---|---|
| EDEM (Altair) | Fluent, STAR-CCM+ | API連成 |
| LIGGGHTS (OSS) | OpenFOAM | CFDEMcoupling (OSS) |
| Rocky DEM (ESSS) | Fluent, CFX | 組み込み連成 |
| Fluent DEM | Fluent | ネイティブ実装 |
| STAR-CCM+ DEM | STAR-CCM+ | ネイティブ実装 |
ボイド率の計算——DEMとCFDの「翻訳問題」
DEM-CFD連成の最大の実装課題は、離散的な粒子位置情報から連続場であるCFDのボイド率分布への変換です。単純なセルカウント法はメッシュサイズが粒子径より小さいときに破綻します。より一般的なDiffusionスムージングやKernel Weighting法が実用化されていますが、カーネル半径の選択次第でドラッグ力計算が30%以上変化します。ESTOCHASTICやMFiX-DEMなど主要コードではこの「マッピング感度テスト」をユーザーが必ず行う必要があります。
風上差分（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次だが、流れの方向を正しく捕捉するため安定性が高い。
実践ガイド
実践ガイド
DEM-CFD連成解析の手順を教えてください。
流動層リアクターを例に説明しよう。
1. 粒子物性の決定: 粒径分布、密度、ヤング率、反発係数、摩擦係数
2. 形状作成: リアクター形状、ガス分散板、出口
3. CFDメッシュ: 粒子径の4〜5倍のセルサイズ
4. 粒子充填: DEM側で初期充填をシミュレーション（自然落下）
5. 流体供給開始: ガス速度を段階的にランプアップ
6. 定常状態の確認: 圧力損失や層高のモニタリング
粒子パラメータの設定
粒子の力学パラメータはどう決めるんですか？
Hertz-Mindlinモデルのパラメータ設定が重要だ。
| パラメータ | 典型値（ガラスビーズ） | 測定方法 |
|---|---|---|
| ヤング率 | $10^7$〜$10^8$ Pa（※低減値） | ナノインデンテーション |
| ポアソン比 | 0.25 | 文献値 |
| 反発係数 | 0.9 | 落下衝突実験 |
| 静摩擦係数 | 0.3〜0.5 | 滑り角実験 |
| 転がり摩擦係数 | 0.01〜0.05 | 安息角キャリブレーション |
ヤング率を低減するのはなぜですか？
実際の材料のヤング率（ガラス: 70 GPa）を使うとDEMタイムステップが極端に小さくなり、計算が非現実的な時間を要する。$10^7$〜$10^8$ Pa程度に下げても、粒子の流動挙動はほとんど変わらないことが多くの研究で示されている。ただし反発係数や接触時間が変わるので、キャリブレーションは必要だ。
安息角キャリブレーション
キャリブレーションって具体的にどうやるんですか？
最も一般的なのは安息角（angle of repose）のキャリブレーションだ。実験で粒子を筒から排出して形成される山の角度を測定し、DEMの摩擦係数を調整して合わせる。
安息角だけでは一意にパラメータが決まらないので、排出流量やドラム回転試験など複数の実験を組み合わせて検証するのが理想的だ。
計算コストの最適化
DEM-CFDは計算コストが高いと聞きましたが、どう対処すればいいですか？
| 手法 | 効果 | 注意点 |
|---|---|---|
| ヤング率の低減 | DEMステップを大きくできる | $10^6$ Pa以下は挙動が変化 |
| Coarse-graining | 代表粒子で複数粒子を表現 | スケーリング則の適切な選択 |
| GPU計算 | DEM計算を大幅高速化 | EDEM, Rocky DEMがGPU対応 |
| 対称性の利用 | 計算領域の縮小 | 周期境界条件の適用 |
錠剤コーティング——製薬GMP対応のDEM-CFD活用事例
医薬品製造において錠剤コーティングの均一性はGMP上の重要品質特性です。パンコーティング機内部ではDEM-CFDにより錠剤の運動軌跡・スプレー露出時間分布が計算でき、コーティング重量変動係数（CV）を実験なしに1%未満に抑えるパン形状最適化が可能になっています。AstraZenecaやNovartisなどの大手製薬メーカーがDEM-CFDをQbD（品質設計）ツールとして規制当局への申請資料に組み込む事例が2010年代後半から増えており、シミュレーション結果のデータインテグリティ管理が新たな課題となっています。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ CFDでも入口条件のミスは下流全体に波及します。
ソフトウェア比較
商用ツール比較
DEM-CFD連成に使えるツールを比較してください。
大きく分けて、ネイティブ実装と外部連成の2タイプがある。
| ツール構成 | DEMエンジン | CFDエンジン | 連成方式 | GPU対応 |
|---|---|---|---|---|
| EDEM + Fluent | EDEM (Altair) | Ansys Fluent | API連成 | DEM側のみ |
| EDEM + STAR-CCM+ | EDEM | STAR-CCM+ | API連成 | DEM側のみ |
| Rocky DEM + Fluent | Rocky (ESSS) | Ansys Fluent | API連成 | DEM側GPU |
| Fluent DEM | 内蔵DEM | Fluent | ネイティブ | 限定的 |
| STAR-CCM+ DEM | 内蔵DEM | STAR-CCM+ | ネイティブ | 限定的 |
| CFDEMcoupling | LIGGGHTS | OpenFOAM | オープンソース | 限定的 |
ネイティブ実装と外部連成はどちらがいいですか？
ネイティブ実装（Fluent DEM, STAR-CCM+ DEM）はセットアップが容易でデータ交換のオーバーヘッドが少ない。ただしDEM機能の成熟度はEDEMやRocky DEMには及ばない。
EDEMは粒子形状（ポリヘドラル、繊維状）の自由度が高く、Rocky DEMはGPUでの大規模計算に強い。数百万粒子の計算ならRocky DEMのGPU実装が圧倒的に速い。
オープンソースのCFDEMcouplingはどうですか？
LIGGGHTS + OpenFOAMの組み合わせで、学術分野で最も広く使われている。カスタマイズの自由度が高く、新しい抗力モデルや連成スキームの実装が容易だ。ただしGUIがなく、設定はすべてスクリプトベースになる。
コスト比較
| ツール | ライセンス | 概算年間コスト |
|---|---|---|
| EDEM | 商用 | CFDソルバーとは別途必要 |
| Rocky DEM | 商用 | CFDソルバーとは別途必要 |
| Fluent/STAR-CCM+ DEM | 基本パッケージに含む | 追加費用なし |
| CFDEMcoupling + LIGGGHTS | GPL/オープンソース | 無償 |
EDEM vs MFiX-DEM vs Rocky——DEM-CFD商用ツール三国志
DEM-CFD商用ツール市場は3強が分割しています。Altair EDEMはFluentやStar-CCM+との連成が成熟しており、鉱業・農業機械分野で長い実績があります。ANSYS MFiX-DEMはオープンソース由来の強みを持ち、学術研究と産業応用の橋渡しをしています。Rocky DEM（ESSS社、現ANSYS）は非球形粒子・粒子破砕のサポートが他を圧倒し、セメント・鉱石粉砕の分野で事実上の標準となっています。粒子数が100万を超える大規模系ではGPU対応の有無が決定的な差となり、Rocky DEMのGPUアクセラレーションは4倍以上の速度向上を実証しています。
選定で最も重要な3つの問い
- 「何を解くか」：DEM-CFD連成に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
先端技術と研究動向
DEM-CFDの最新研究にはどんなものがありますか？
いくつかの方向性を紹介しよう。
Coarse-graining法
粒子数を減らすため、1つの「代表粒子」で複数の実粒子を表現するcoarse-graining（粗視化）法が注目されている。代表粒子の径を実粒子の$k$倍にすると、粒子数は$k^3$分の1になる。
大きな粒子にしたら物理が変わりませんか？
適切なスケーリング則が必要だ。力のスケーリングとして、Sakai et al.（2009）のcoarse-grainingモデルでは抗力を$k^3$倍に補正する。Exact Scaling法、Statistical Scaling法など複数の流派がある。
非球形粒子
実際の粉体は球形ではない。非球形粒子の取り扱いが重要な研究テーマだ。
| 手法 | 概要 | 計算コスト |
|---|---|---|
| Multi-sphere | 複数の球でクラスターを構成 | 中 |
| Superquadrics | 楕円体・直方体の一般化形状 | 中〜高 |
| Polyhedra | 多面体で任意形状を表現 | 高 |
| Level Set DEM | 符号付き距離関数で形状を定義 | 高 |
EDEMの非球形粒子は何を使っているんですか？
EDEMはMulti-sphere法が標準で、複雑な形状をオーバーラップする球の集合体で近似する。Rocky DEMはPolyhedra法をサポートしており、より忠実な形状表現が可能だ。
機械学習による抗力モデル
DNSやLBM（格子ボルツマン法）で計算した詳細な流体力データを教師データとして、ニューラルネットワークで抗力モデルを構築する研究が進んでいる。Beetstra et al.やTang et al.のモデルが代表的だ。
従来の相関式より精度が上がるんですか？
特に高ボイド率の遷移領域や非球形粒子の抗力では、従来の相関式より大幅に精度が向上する。ただしDNSデータの生成自体が計算コストが高い課題がある。
機械学習×DEM——粒子接触モデルのデータ駆動型構築
従来のDEM接触モデル（Hertz-Mindlinモデル等）は球形粒子の単純接触を仮定しますが、実際の粒子は非球形・表面粗さ・湿気による付着力が複雑に絡み合います。2020年代に入り、分子動力学シミュレーションや実験データからニューラルネットワークで接触力を直接学習する「ML-DEM」研究が加速しています。ETH ZürichとDCSEの共同研究では、非球形砂粒子の摩擦係数を95%の精度で再現するモデルを開発し、GPU上での推論実装により100万粒子系でもリアルタイムに近い速度を達成しつつあります。
トラブルシューティング
トラブルシューティング
DEM-CFD連成でよくあるトラブルを教えてください。
順番に見ていこう。
1. 粒子がメッシュを突き抜ける
症状: 粒子が壁面やCFDドメインの外に出てしまう。
対策:
- DEMタイムステップが大きすぎる。Rayleigh時間の20%以下に設定
- 壁面のDEM境界がCFDメッシュと一致していることを確認
- 接触力のばね定数が適切か確認
2. 圧力損失が実験と合わない
流動層の圧力損失がずれる場合はどうすればいいですか？
対策:
- 抗力モデルの選択を見直す（Gidaspow vs. Koch-Hill）
- CFDセルサイズが粒子径の3〜5倍であることを確認
- ボイド率計算のスムージング方法を変更
- 粒子充填率が実験と一致しているか確認
3. 計算が異常に遅い
原因と対策:
- DEMタイムステップが小さすぎる: ヤング率を低減（$10^7$〜$10^8$ Pa）
- 粒子数が多すぎる: coarse-graining法を適用
- 負荷分散: DEMとCFDで異なる並列分割戦略を使い、各コアの負荷を均等化
4. 粒子が不自然に凝集する
症状: 粒子が固まったまま動かない。
対策:
- 摩擦係数が高すぎないか確認（静摩擦 < 0.5が一般的）
- 転がり摩擦が過大でないか確認
- 粘着力モデル（JKR等）が意図せず有効になっていないか確認
5. ツール固有の注意点
| ツール | 注意点 |
|---|---|
| EDEM + Fluent | 連成時間間隔がCFDタイムステップと一致していることを確認 |
| Rocky DEM | GPU計算時のメモリ制限に注意（粒子数×属性データ） |
| CFDEMcoupling | ボイド率平滑化の半径パラメータ（voidfractionModel）が解に影響 |
| Fluent DEM | ネイティブDEMの粒子形状は球形のみ（非球形はUDF必要） |
粒子が重なって発散する——DEM計算の典型的崩壊パターン
DEM-CFD計算が発散する最も多い原因は「粒子の過剰重なり（overlap）」です。初期配置で粒子同士が重なっていると、接触バネが無限大の斥力を生じて速度が爆発します。対策として粒子を低密度ガスとして拡散充填するrun-in計算が標準的なプロセスです。DEM安定条件はΔt < 0.1×√(m/k_n)（m:粒子質量、k_n:法線剛性）で与えられ、粒子径が1 mmから0.1 mmに変わるとΔtが1/10になって計算時間が10倍になります。小径粒子系では「coarse-graining（粗視化）」による仮想的な粒径拡大が実務上の回避策です。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——DEM-CFD連成の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告