Euler型粒体モデル
理論と物理
概要
先生、Euler型粒体モデルって何ですか？粒子を連続体として扱うんですか？
その通り。Euler型粒体モデル（Eulerian Granular Model）は、粉体や顆粒の集合体を「粒体相」という擬似連続体として扱い、気体相とともにEuler-Euler法の枠組みで解く手法だ。流動層、空気輸送、サイクロン、粉体混合などの固気二相流で広く使われている。
DEM-CFDとは何が違うんですか？
DEM-CFDは個々の粒子を離散的に追跡するが、粒子数が数百万〜数十億になる工業スケールでは計算コストが膨大になる。Euler型粒体モデルは粒子群を連続体として扱うため、大規模な系でも現実的な計算時間で解ける。ただし個々の粒子の接触力や形状効果は直接扱えない。
支配方程式
粒体相の方程式はどうなりますか？
粒体相の運動を記述するために、KTGF（Kinetic Theory of Granular Flow：粒体の運動論）を使う。固相の連続の式は通常のEuler-Eulerと同じだが、固相の応力テンソルにKTGFから導かれる構成則を用いるのが特徴だ。
KTGFって気体分子運動論の粒子版ということですか？
まさにその通り。粒子の速度揺動をGranular Temperature $\Theta_s$ で特徴づける。
Granular Temperatureの輸送方程式は次の通りだ。
右辺は順に、せん断による生成、拡散（$\kappa_s$は拡散係数）、非弾性衝突による散逸（$\gamma_s$）、気体-固体間のエネルギー交換（$\phi_{gs}$）だ。
固相の構成則
固相圧力とか粘性はどう決まるんですか？
固相圧力 $p_s$ はGranular Temperatureと体積分率から求まる。Lun et al.（1984）のモデルが代表的だ。
ここで $e_{ss}$ は粒子間の反発係数、$g_0$ は動径分布関数だ。$g_0$ は体積分率が最密充填率 $\alpha_{s,max}$（約0.63）に近づくと急増し、粒子が密集した状態での接触圧力を表現する。
| 構成則 | モデル例 | 物理量 |
|---|---|---|
| 固相圧力 | Lun, Syamlal-O'Brien | $p_s(\alpha_s, \Theta_s)$ |
| 固相粘性 | Gidaspow, Syamlal | $\mu_s(\alpha_s, \Theta_s)$ |
| 固相体積粘性 | Lun et al. | $\lambda_s$ |
| 摩擦応力 | Schaeffer, Johnson-Jackson | 密充填領域の応力 |
砂時計の物理——粒体は「流体」でも「固体」でもない
砂時計の砂が「流れる」ように見えるのに、砂山の斜面で静止するのはなぜか。この問いにEuler型粒体モデルが答える鍵があります。粒体は流体力学的な連続体方程式と固体力学的な弾性圧力の両方を必要とする「第三の状態」です。granular temperatureという概念は、分子の熱運動に倣い粒子速度の揺らぎを温度と見なす独創的なアイデアで、1980年代にJenkins & Richmanが気体分子運動論から導出しました。この理論なしにはCFBボイラーの粒子循環速度を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$。時間刻みの安定性に直結 |
数値解法と実装
数値解法の詳細
Euler型粒体モデルの数値的なポイントを教えてください。
最大の難しさは、固相体積分率が最密充填率 $\alpha_{s,max}$ に近づいたときの取り扱いだ。この領域では固相圧力が急激に増大し、数値的に不安定になりやすい。
摩擦応力（Frictional Stress）モデルが重要で、$\alpha_s > \alpha_{s,min}$（通常0.5）のとき、SchafferモデルやJohnson & Jacksonモデルで摩擦圧力と摩擦粘性を追加する。
抗力モデルの選択
固気二相流の抗力モデルはどう選べばいいですか？
Gidaspowモデルが最も一般的だ。これはErgun式（密充填領域）とWen-Yu式（希薄領域）を体積分率0.8で切り替える。
| 領域 | $\alpha_g$ | モデル | 式 | ||
|---|---|---|---|---|---|
| 密充填 | $< 0.8$ | Ergun | $\beta = 150 \frac{\alpha_s^2 \mu_g}{\alpha_g d_s^2} + 1.75 \frac{\alpha_s \rho_g \ | \mathbf{u}_g - \mathbf{u}_s\ | }{d_s}$ |
| 希薄 | $\geq 0.8$ | Wen-Yu | $\beta = \frac{3}{4} C_D \frac{\alpha_s \alpha_g \rho_g \ | \mathbf{u}_g - \mathbf{u}_s\ | }{d_s} \alpha_g^{-2.65}$ |
切り替えが不連続だと問題になりませんか？
その通りで、切り替え点での不連続がボイド率の振動を引き起こすことがある。Huilin-Gidaspowモデルはスムーズな遷移関数を導入して改善している。Syamlal-O'Brienモデルも全領域で連続な式を使うため安定性が高い。
OpenFOAMでの実装
OpenFOAMではどのソルバーを使いますか？
multiphaseEulerFoam がEulerian Granularに対応している。KTGFの各構成則は kineticTheoryModel クラスで選択できる。主な設定は constant/phaseProperties で行う。
Fluentでの設定
Ansys FluentではEulerian Multiphase Model内でGranular Phaseを有効にする。重要な設定項目は以下の通り。
| 設定 | 推奨 | 備考 |
|---|---|---|
| Granular Viscosity | Gidaspow | 標準的 |
| Granular Bulk Viscosity | Lun et al. | 体積粘性 |
| Frictional Viscosity | Schaeffer | 密充填領域 |
| Packing Limit | 0.63 | ランダム充填率 |
| Restitution Coefficient | 0.9 | ガラスビーズ典型値 |
KTGF収束の壁——「圧力が負になる」現象との闘い
Euler型粒体モデルの実装で最初につまずくのが、固相圧力が負値になって発散するという問題です。これはgranular temperatureが急激に0へ収束するときに起きる数値的不安定で、フロア値（最小値1e-10 m²/s²程度）の設定が事実上の必須対策となっています。ANSYSのドキュメントには「Partial Differential Equationとして解くかAlgebraic近似で解くか」の選択肢が示されていますが、高密度充填（α_s > 0.4）ではPDE法でないと精度が出ないことが多く、計算コストと精度のトレードオフを現場で判断する必要があります。
風上差分（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次だが、流れの方向を正しく捕捉するため安定性が高い。
実践ガイド
実践ガイド
Euler型粒体モデルの解析手順を教えてください。
流動層リアクターの解析を例に説明しよう。
1. 形状定義: リアクター形状、分散板、フリーボード領域を含める
2. メッシュ: 粒子径の10〜20倍のセルサイズ、六面体推奨
3. 初期条件: 固相を静止充填状態（$\alpha_s = 0.6$）で下部に配置
4. 気体入口: 分散板からの均一速度（最小流動化速度 $U_{mf}$ の2〜10倍）
5. 非定常計算: $\Delta t = 10^{-4}$〜$10^{-3}$ s で開始
6. モニタリング: 層高変動、圧力損失、ガスバイパス
最小流動化速度の確認
最小流動化速度ってどう求めるんですか？
Ergun式から推定できる。
Geldart分類でA粒子（$20 < d_p < 100$ μm）は均一膨張、B粒子（$100 < d_p < 1000$ μm）は気泡流動化を示す。CFDでもこの違いを再現できることが重要な検証ポイントだ。
2Dと3D
流動層は2Dで計算してもいいですか？
2D計算は定性的な挙動の確認やパラメータスタディには有用だが、定量的な予測には3Dが必要だ。特に気泡の合体頻度や飛び出し挙動は2Dと3Dで大きく異なる。ただし計算コストの制約から、まず2Dで方向性を確認し、最終評価を3Dで行うアプローチが現実的だ。
メッシュ感度
メッシュサイズの影響は大きいですか？
非常に大きい。Euler型粒体モデルはメッシュ依存性が強く、セルサイズが変わると気泡径や層膨張率が変化する。フィルタリングモデル（Filtered TFM）やcoarse-grained simulationの研究が進んでいるが、実務ではメッシュ感度解析が必須だ。
| メッシュ | 典型値 | 備考 |
|---|---|---|
| 微細 | $5 d_p$ | 研究用、計算コスト大 |
| 標準 | $10$〜$20 d_p$ | 実務的なバランス |
| 粗い | $> 30 d_p$ | メソスケール構造を解像できない |
FCC装置の粒子循環——石油精製を支える流動の匠
流動接触分解（FCC）装置は、ライザー管内で触媒粒子（平均径70 µm、密度1500 kg/m³）を毎秒10 m以上の速度で輸送しながら、原油を軽質油に変換します。世界の石油精製能力の約30%がFCC技術に依存しており、粒子の循環量制御はそのまま生産量に直結します。Euler型粒体モデルでFCCライザーをシミュレーションする場合、メッシュサイズを粒子径の10倍以下（≈1 mm）にするとコア-アニュラス流れ構造が再現できますが、実際の装置スケール（直径1 m × 高さ30 m）では計算量が爆発するため、クラスター分解効率（CRE）を組み込んだコースグレインモデルが産業界の標準になっています。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ CFDでも入口条件のミスは下流全体に波及します。
ソフトウェア比較
商用ツール比較
Euler型粒体モデルに対応しているツールを教えてください。
| ツール | KTGF実装 | 摩擦モデル | 多粒径対応 | 特徴 |
|---|---|---|---|---|
| Ansys Fluent | Lun, Syamlal, Gidaspow | Schaeffer, Johnson-Jackson | DQMOM, QMOM | 最も充実したGranularオプション |
| STAR-CCM+ | 標準KTGF | 摩擦粘性モデル | 複数固相 | ポリヘドラルメッシュ対応 |
| Ansys CFX | 基本KTGF | 限定的 | 複数固相 | 結合型ソルバーの安定性 |
| OpenFOAM | multiphaseEulerFoam | Schaeffer | 複数固相 | 完全カスタマイズ可能 |
| Barracuda Virtual Reactor | CPFD法 | 独自モデル | Parcelベース | 流動層専用、GPU対応 |
Barracuda Virtual Reactorって何ですか？
CPFD Software社（現Convergent Science傘下）が開発した流動層専用ツールだ。CPFD法（Computational Particle Fluid Dynamics）はEuler-Lagrangeの一種だが、粒子をparcelで代表して計算効率を上げている。流動層に特化しているため、FCC（流動接触分解）やCFB（循環流動層ボイラー）の設計に広く使われている。
用途別推奨
| 用途 | 推奨ツール | 理由 |
|---|---|---|
| FCC/CFBプロセス設計 | Barracuda, Fluent | 流動層特化 or 汎用性 |
| 粉体輸送 | Fluent, STAR-CCM+ | 摩耗・堆積モデル連携 |
| 学術・モデル開発 | OpenFOAM | ソースコードアクセス |
| 原子力（燃料デブリ） | CFX | 安全解析実績 |
コスト的にはどうですか？
Barracuda Virtual Reactorは流動層専用なので、流動層しかやらないなら効率的だ。汎用CFDも必要なら Fluent や STAR-CCM+ の Eulerian Granular で一本化するのがコスト合理的だ。OpenFOAMは無償だが、Granular モデルの検証事例は商用ツールに比べて少ない。
Fluent vs STAR-CCM+ vs OpenFOAM——粒体モデルの商用格差
Euler型粒体モデルの商用ソフト比較で最も議論になるのは「粒度分布（PSD）の扱い」です。ANSYS Fluentは複数固相を別々のEuler相として定義できますが、粒径ビン数が増えると計算コストが線形増加します。STAR-CCM+のPolydisperse Granular Flow（PGF）モデルは粒径分布を矩モデルで近似し、5ビン相当の計算を単一相で完結させる独自アーキテクチャを持ちます。一方、無料のOpenFOAMはKTGF実装が比較的シンプルで、カスタムモデルの追加に強みがありますが、ライセンスサポートなしでの産業応用には相応のCFDスキルが必要です。
選定で最も重要な3つの問い
- 「何を解くか」：Euler型粒体モデルに必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
先端技術と研究動向
Euler型粒体モデルの最新研究にはどんなものがありますか？
いくつかの重要な方向性がある。
Filtered TFM（フィルタリング二流体モデル）
粗いメッシュでもメソスケール構造（クラスターやストリーマー）の効果を取り込むために、フィルタリングの概念を導入する手法だ。LESにおけるサブグリッドモデルのTFM版と言える。
Igci et al.（2008）やMilioli et al.（2013）のフィルタリング抗力モデルが代表的で、粗いメッシュでも細かいメッシュと同等の巨視的挙動を予測できる。
Polydisperse（多粒径）モデル
実際の粉体は粒径がバラバラですよね？
その通り。多粒径を扱う方法は2つある。
| 手法 | 概要 | 計算コスト |
|---|---|---|
| Multi-solid model | 粒径ごとに固相を定義 | 相数×輸送方程式数 |
| DQMOM/QMOM | モーメント法で分布を追跡 | 効率的 |
粒径分布の偏析（サイズセグレゲーション）は流動層の性能に大きく影響する。大きな粒子は沈降、小さな粒子は飛散しやすい。これを正確に予測するには多粒径モデルが必須だ。
反応を伴う流動層
石炭ガス化、バイオマス燃焼、化学ルーピング（Chemical Looping Combustion）など、反応を伴う流動層のCFDが活発に研究されている。粒子の収縮・膨張（Shrinking Core Model）、揮発分放出、チャー燃焼などのサブモデルと連成する必要がある。
反応があると計算はさらに複雑になりそうですね。
確かに計算コストは増えるが、Fluent等ではHeterogeneous Reaction Modelが統合されていて、比較的簡単に設定できる。OpenFOAMの coalChemistryFoam も石炭燃焼用の枠組みを提供している。
機械学習による抗力補正
粒子形状の先端研究——球ではない粒体のCFD
実際の工業用粒子（砕砂、非球形触媒、生体細胞）のほとんどは球ではありません。非球形粒子の抗力は同体積球の1.2〜2.5倍に達し、球形近似では流動層高さを30%以上過小評価するケースが報告されています。2010年代から盛んになった「Euler-Euler with morphology tensor」アプローチは、粒子のアスペクト比と配向を場の量として輸送方程式に組み込み、圧延プロセスの繊維配向予測でLES実験値と5%以内の一致を達成しました。OpenFOAMのmultiphaseEulerFoamにはまだ標準実装されておらず、研究グループが独自パッチで対応しているのが現状です。
トラブルシューティング
トラブルシューティング
Euler型粒体モデルでよくあるトラブルを教えてください。
順番に見ていこう。
1. 固相体積分率が最密充填率を超える
症状: $\alpha_s > \alpha_{s,max}$ となり計算が発散。
対策:
- Packing Limitを正しく設定（球形粒子で0.63、実測値推奨）
- 摩擦圧力モデル（Frictional Pressure）が有効になっていることを確認
- タイムステップを小さくする（$10^{-4}$ s以下）
- 固相体積分率のunder-relaxationを下げる（0.2〜0.3）
2. 流動化しない / 圧力損失が合わない
ガスを流しても粒子が動かないんですが…
対策:
- ガス流速が $U_{mf}$ を超えていることを確認
- 抗力モデルの選択を見直す（Gidaspow, Syamlal-O'Brien を比較）
- 粒子径と密度が正しく設定されているか確認
- 初期の固相体積分率が高すぎないか確認（0.55〜0.60が適当）
- 分散板の境界条件が正しいか確認（均一速度入口）
3. 非物理的な粒子飛散
症状: 粒子がフリーボード上部に過剰に飛散。
対策:
- フリーボード領域を十分に長く取る
- 出口境界条件でbackflowの固相体積分率を0に設定
- メッシュ解像度がフリーボードで粗すぎないか確認
4. Granular Temperatureが異常値
対策:
- Granular Temperature方程式を代数近似（Algebraic）から偏微分方程式（PDE）に変更して安定性を確認
- 反発係数 $e_{ss}$ が0.5〜0.99の範囲にあることを確認
- 初期のGranular Temperatureを小さな正の値（$10^{-5}$）に設定
5. ツール固有の注意点
| ツール | 注意点 |
|---|---|
| Fluent | Granular Temperature方程式のAlgebraic近似は希薄流では不正確。Dense bedではPDE推奨 |
| STAR-CCM+ | 固相壁面境界条件でJohnson-Jacksonのスペキュラー係数設定に注意 |
| OpenFOAM | kineticTheoryModelの選択肢がバージョンで異なる。tutorialケースで動作確認推奨 |
| Barracuda | CPFD法独自のパラメータ（close-pack volume fraction等）が結果に敏感 |
「粒子が壁に張り付く」——壁境界条件の落とし穴
Euler型粒体モデルで流動層シミュレーションを行うと、固相体積分率が壁面付近で非現実的に高くなる（α_s→0.63以上）現象が頻発します。原因の大半は壁面でのgranular temperature境界条件の設定ミスで、Johnson-Jacksonモデルの鏡面反射係数（specularity coefficient）を0に設定すると壁面に粒子が堆積し続けます。実務では0.05〜0.25の範囲でキャリブレーションしますが、値が実験データなしに決められないことが設計段階での大きな課題です。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——Euler型粒体モデルの問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告