個数密度関数法（PBM）
理論と物理
概要
先生、個数密度関数法（PBM）って何ですか？
PBM（Population Balance Model）は、気泡・液滴・粒子のサイズ分布の時間空間変化を追跡するモデルだ。合体（coalescence）、分裂（breakup）、核生成（nucleation）、成長（growth）によってサイズ分布が変化する過程を記述する。
どんな場面で使うんですか？
気泡塔での気泡径分布、乳化プロセスでの液滴径分布、晶析での結晶サイズ分布、エアロゾルの粒径変化など、分散相のサイズ分布が重要な問題全般で使う。Euler-Euler法と組み合わせて使うのが一般的だ。
支配方程式
PBMの方程式を教えてください。
体積 $V$ を持つ粒子（気泡・液滴）の個数密度 $n(V, \mathbf{x}, t)$ の輸送方程式は次の通りだ。
右辺の各項は以下を表す。
- $B_{coal}$: 合体による生成（2つの小さな粒子が合体してVのサイズになる）
- $D_{coal}$: 合体による消滅（サイズVの粒子が他と合体してより大きくなる）
- $B_{break}$: 分裂による生成（大きな粒子が分裂してVのサイズが生まれる）
- $D_{break}$: 分裂による消滅（サイズVの粒子が分裂して小さくなる）
合体と分裂の速度はどうモデル化するんですか？
代表的なクロージャモデルを示そう。
| プロセス | モデル例 | 駆動メカニズム |
|---|---|---|
| 合体頻度 | Prince & Blanch (1990) | 乱流衝突 + 膜排出 |
| 合体頻度 | Luo (1993) | 乱流エネルギー |
| 分裂頻度 | Luo & Svendsen (1996) | 乱流渦による分裂 |
| 分裂頻度 | Martínez-Bazán (2010) | 慣性-表面張力バランス |
Sauter平均径
サイズ分布から代表径を求めるにはどうするんですか？
Sauter平均径 $d_{32}$ が最もよく使われる。
$m_k = \int_0^\infty V^{k/3} n(V) dV$ は分布の $k$ 次モーメントだ。$d_{32}$ は体積-表面積比に対応し、物質移動や反応速度の評価に適している。
個体群バランス——工学の「進化論」としての多相流理論
Population Balance Equation（PBE）は、サイズ・年齢・組成などの「内部座標」を持つ個体群の時空間変化を記述する汎用的な枠組みです。生物の個体群動態（Lotka-Volterra方程式）と数学的に同型であり、気泡分布・結晶粒径分布・細胞濃度分布まで統一的に扱えます。化学工学でのPBE応用はRandolph & Larson（1971年）の結晶化論文が草分けで、50年後の今はCFDと連成したC-PBE（Coupled PBE）として医薬品製造・気泡塔・液液抽出の設計ツールに進化しています。
各項の物理的意味
- 時間項 $\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$。時間刻みの安定性に直結 |
数値解法と実装
数値解法の詳細
PBMの数値的な解き方を教えてください。
PBEは体積（サイズ）を追加の独立変数とする積分微分方程式なので、直接解くのは困難だ。以下の離散化手法が使われる。
| 手法 | 概要 | 計算コスト | 精度 |
|---|---|---|---|
| MUSIG (Multi-Size Group) | サイズ空間をビンに分割 | 高（ビン数に比例） | ビン数依存 |
| QMOM (Quadrature MOM) | モーメント法 + 求積公式 | 低（6〜8モーメント） | 良好 |
| DQMOM | 直接求積モーメント法 | 低 | 良好 |
| S-Gamma | 2パラメータ分布仮定 | 最低 | 分布形状に制約 |
MUSIGとQMOMの違いを教えてください。
MUSIGはサイズ空間を10〜30のビン（size group）に分割し、各ビンの個数密度を輸送方程式として解く。精度が高いが、ビン数だけ追加の輸送方程式を解くため計算コストが大きい。
QMOM（McGraw, 1997）はサイズ分布を直接解かず、最初の数個のモーメント（$m_0, m_1, ..., m_5$等）の輸送方程式を解く。Product-Difference Algorithm（PDA）やWheeler法で求積公式の節点と重みを決定し、合体・分裂のソース項を計算する。
Fluentでの実装
Ansys FluentではEulerian Multiphase + Population Balance Modelで利用可能だ。
| 手法 | Fluent名称 | ビン/モーメント数 |
|---|---|---|
| Discrete Method | Size Group | 10〜30ビン |
| QMOM | Quadrature MOM | 6モーメント |
| DQMOM | Direct QMOM | 2〜4環境 |
| SMM | Standard MOM | 6モーメント |
STAR-CCM+での実装
STAR-CCM+では S-Gamma モデル（Lo, 2000）が標準で、2パラメータ（平均径と分散）で分布を追跡する。計算コストが最も低いが、分布形状がガンマ分布に制約される。MUSIGも利用可能だ。
OpenFOAMでの実装
OpenFOAMでは populationBalanceModel クラスがv2006以降で利用可能だ。MUSIG（iMUSIG含む）とQMOMが実装されている。設定は constant/phaseProperties 内の populationBalanceCoeffs で行う。
QMOM vs DQMOM——モーメント法のトレードオフ
CFD-PBEの計算コストを現実的に抑えるためにモーメント法が使われます。QMOM（Quadrature Method of Moments）は分布関数を求積点で近似し、モーメント方程式のみを解くことで分布全体の時間発展を追跡します。しかし分裂・合体が複雑な系では求積点が交差する「moment inversion問題」が発生します。DQMOM（Direct QMOM）は各求積点の位置と重みを直接輸送方程式で解くことでこの問題を回避しますが、ソース項の非線形性が収束を困難にします。気泡塔CFDのベンチマークでは、QMOMとDQMOMの気泡径分布予測が10〜20%の差を持つケースが多く、モデル選択が結果に無視できない影響を与えます。
風上差分（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次だが、流れの方向を正しく捕捉するため安定性が高い。
実践ガイド
実践ガイド
PBM解析の手順を教えてください。
気泡塔の気泡径分布解析を例に説明しよう。
1. Euler-Euler設定: 液相（連続相）+ 気相（分散相）
2. PBM有効化: 気相にPopulation Balanceを設定
3. サイズ範囲: 最小1 mm〜最大20 mm（10〜20ビン）
4. 合体モデル: Prince-Blanch or Luo
5. 分裂モデル: Luo-Svendsen
6. 抗力: サイズ依存（各ビンの径でIshii-Zuberを適用）
7. 初期分布: 均一径 or 実験分布で初期化
8. 後処理: $d_{32}$、局所気泡径分布の時間空間変化
サイズ範囲とビン数の設定
ビン数はどのくらい必要ですか？
| 手法 | 推奨ビン/モーメント数 | 計算コスト増加 |
|---|---|---|
| MUSIG | 15〜25ビン | 輸送方程式×15〜25追加 |
| iMUSIG | 3〜5速度グループ×5〜10サイズ | 中程度 |
| QMOM | 6モーメント | 輸送方程式×6追加 |
| S-Gamma | 2パラメータ | 最小限 |
iMUSIGって何ですか？
inhomogeneous MUSIG の略で、CFX が最初に実装した手法だ。サイズグループをさらに速度グループに分け、大きな気泡と小さな気泡に異なる速度場を持たせる。Tomiyama揚力の符号反転（小気泡は壁面へ、大気泡は中心へ）を再現するのに不可欠だ。
検証のポイント
PBMの結果をどう検証しますか？
以下の実験量との比較が標準だ。
| 計測量 | 計測手法 | 備考 |
|---|---|---|
| 局所気泡径分布 | 光ファイバープローブ | 確率密度関数 |
| Sauter平均径 $d_{32}$ | 位相ドップラー法 | 径と速度同時計測 |
| ガスホールドアップ | 差圧計、ワイヤーメッシュセンサー | 断面平均・局所 |
| 気泡上昇速度 | 高速カメラ + 画像処理 | サイズ-速度相関 |
医薬品結晶化——粒径分布を設計するCFD-PBE
医薬品原薬の製造プロセスで最重要な単位操作の一つが「結晶化（Crystallization）」です。結晶の粒径分布（PSD）は溶解度・ろ過性・生体内溶解速度を決定し、最終製品の効き目に直結します。バッチ晶析槽のCFD-PBEにより、攪拌速度・温度プロファイル・種結晶量を最適化してCVが5%以内のPSDを実現した事例が報告されています。AstraZenecaは2015年以降、主要化合物の晶析条件をCFD-PBEなしに決定しない方針を社内で決定したとされており、規制申請書にもシミュレーション結果が添付されています。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ CFDでも入口条件のミスは下流全体に波及します。
ソフトウェア比較
商用ツール比較
PBMに対応しているツールを比較してください。
| ツール | PBM手法 | 合体モデル | 分裂モデル | 特徴 |
|---|---|---|---|---|
| Ansys Fluent | MUSIG, QMOM, DQMOM, SMM | Luo, Prince-Blanch | Luo-Svendsen, Laakkonen | 最多の手法選択肢 |
| STAR-CCM+ | S-Gamma, MUSIG | 標準モデル | 標準モデル | S-Gammaが軽量 |
| Ansys CFX | iMUSIG, MUSIG | Prince-Blanch | Luo-Svendsen | iMUSIGの元祖 |
| OpenFOAM | MUSIG, QMOM | 基本モデル | 基本モデル | カスタマイズ自由 |
用途別推奨
| 用途 | 推奨ツール | 理由 |
|---|---|---|
| 気泡塔リアクター | Fluent, CFX | PBM + Euler-Eulerの成熟度 |
| 乳化・分散 | Fluent | QMOM + 液液分散 |
| 晶析プロセス | Fluent (DQMOM) | 核生成 + 成長 |
| 二相流配管（ボイド分布） | CFX (iMUSIG) | 速度グループ分割 |
| 学術研究 | OpenFOAM | 新しいクロージャモデル実装 |
晶析でもPBMを使うんですね。
晶析では核生成、結晶成長、凝集、破砕がサイズ分布を決める。PBMはまさにこれらの過程を記述するための枠組みだ。Fluentの DQMOM は結晶の成長速度がサイズ依存する場合にも対応できる。
化学工学のプロセスシミュレーションでは、PBMの結果をAspen Plus等のプロセスシミュレータにフィードバックしてプラント全体の最適化を行うワークフローも確立されつつある。
gPROMS vs Fluent PBM——PBE連成シミュレーションの産業標準
Population Balance Modelingの商用エコシステムは、プロセスシミュレータ側とCFD側の二つの軸で発展しています。Process Systems EnterpriseのgPROMSはバッチ晶析・乳化プロセスのPBEを高度に実装しており、製薬業界で広く使われています。CFD側ではANSYS Fluent のPBMモジュールがQMOM/DQMOMを標準実装し、気泡塔・液液分散系のCFD-PBEが産業案件として受注されています。OpenFOAMのpopulationBalanceFoamはオープンソースとして急速に機能拡張されており、学術から産業への橋渡し役を担っています。
選定で最も重要な3つの問い
- 「何を解くか」：個数密度関数法（PBM）に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
先端技術と研究動向
PBMの最新研究にはどんなものがありますか？
いくつかの方向性を見ていこう。
DNSからのクロージャモデル構築
気泡の合体・分裂過程をVOF法のDNSで直接計算し、合体効率や分裂頻度のクロージャ相関を改良する研究が進んでいる。Tryggvason et al.（Notre Dame）やLu & Tryggvason（2013）のFront-Tracking DNSが先駆的だ。
多変数PBM
サイズだけでなく他の変数も追跡できますか？
2Dや多次元のPBMとして、サイズと組成、サイズと温度、サイズと形状など複数の内部変数を同時に追跡する研究がある。晶析での結晶形状（多形）制御や、噴霧での液滴温度-サイズ結合などに応用される。ただし次元の呪いで計算コストが急増するため、QMOMやCQMOMが効率的な解法として研究されている。
機械学習による合体・分裂カーネル
DNSデータを教師データとして、合体効率や分裂頻度をニューラルネットワークで予測するモデルが提案されている。従来の経験相関では表現できない複雑なパラメータ依存性を学習できる利点がある。
LES + PBM
RANSベースのPBMではなく、LESと組み合わせて乱流の大規模構造がサイズ分布に与える影響を直接捕捉する研究が進んでいる。気泡塔の過渡的な気泡径変化や混合性能の予測精度が向上する。
LESとPBMの組み合わせは計算コストが高そうですね。
その通り。MUSIG+LESは非常に重いので、S-GammaやQMOMのような低コスト手法との組み合わせが現実的だ。
Coffee Break よもやま話
Moments Closure——PBEのクロージャー問題との格闘
PBEを解析的に扱う上での根本的困難は「クロージャー問題」です。気泡分裂頻度・合体速度をサイズの関数として定式化する際、二体以上の相関を表す高次モーメントが低次モーメント方程式に現れ、閉じた方程式系が得られません。このクロージャーには物理的仮定（Homogeneous Turbulence近似等）が必要で、その妥当性が予測精度の上限を決めます。2020年代にはデータ駆動型クロージャー（機械学習で実験データからクロージャーモデルを直接学習）の研究が急増しており、従来の半経験式に比べて予測精度が30%向上したという報告もあります。
RANSベースのPBMではなく、LESと組み合わせて乱流の大規模構造がサイズ分布に与える影響を直接捕捉する研究が進んでいる。気泡塔の過渡的な気泡径変化や混合性能の予測精度が向上する。
LESとPBMの組み合わせは計算コストが高そうですね。
その通り。MUSIG+LESは非常に重いので、S-GammaやQMOMのような低コスト手法との組み合わせが現実的だ。
Moments Closure——PBEのクロージャー問題との格闘
PBEを解析的に扱う上での根本的困難は「クロージャー問題」です。気泡分裂頻度・合体速度をサイズの関数として定式化する際、二体以上の相関を表す高次モーメントが低次モーメント方程式に現れ、閉じた方程式系が得られません。このクロージャーには物理的仮定（Homogeneous Turbulence近似等）が必要で、その妥当性が予測精度の上限を決めます。2020年代にはデータ駆動型クロージャー（機械学習で実験データからクロージャーモデルを直接学習）の研究が急増しており、従来の半経験式に比べて予測精度が30%向上したという報告もあります。
トラブルシューティング
トラブルシューティング
PBMでよくあるトラブルを教えてください。
順番に見ていこう。
1. 気泡径が実験と大きくずれる
対策:
- 合体・分裂モデルの係数を確認（デフォルト値が適切でない場合がある）
- 表面張力の値が正しいか確認（分裂頻度は表面張力に強く依存）
- 乱流散逸率 $\varepsilon$ の計算精度を確認（分裂頻度は$\varepsilon$の関数）
- メッシュ依存性を確認（粗いメッシュでは$\varepsilon$が不正確になる）
2. サイズ分布が両端に偏る
最小ビンと最大ビンに集中してしまいます…
対策:
- サイズ範囲を広げる（最小・最大が実際の分布範囲をカバーしているか確認）
- ビン数を増やしてサイズ解像度を上げる
- 合体/分裂のバランスを確認（一方が極端に強いと偏る）
3. PBM方程式の残差が収束しない
対策:
- PBMの under-relaxation を下げる（0.3〜0.5）
- まずPBMなしでEuler-Euler部分を収束させてからPBMを有効化
- 初期の気泡径分布を物理的に妥当な値で設定
- タイムステップを小さくする
4. QMOMでモーメントが非物理的になる
症状: 負のモーメントや求積公式の破綻。
対策:
- モーメントのリアライザビリティ条件を確認
- QMOMからDQMOMに切り替えて安定性を改善
- モーメントの上限・下限クリッピングを有効化
5. ツール固有の注意点
| ツール | 注意点 |
|---|---|
| Fluent | MUSIGのビン数が多いとメモリ消費が急増。QMOMの方がスケーラブル |
| CFX | iMUSIG の速度グループ数は3〜5が推奨。多すぎると不安定 |
| STAR-CCM+ | S-Gamma はガンマ分布仮定なので、二峰性分布は表現できない |
| OpenFOAM | PBMの実装はバージョンで変更が多いので、最新ドキュメントを確認 |
粒径が発散する——PBEの数値安定性問題の診断
CFD-PBEで最も厄介なのが「計算途中で気泡径や結晶サイズが非現実的な値に発散する」現象です。多くの場合、分裂・合体のカーネル関数がサイズに対して高次の依存性を持ち、大径粒子が出現すると合体速度が爆発的に増加するポジティブフィードバックが起きます。対策として粒径の上限制限子（Dmax）を設定することが応急処置ですが、この値の物理的根拠が必要です。より根本的にはモーメント法の離散化スキームを上流差分から高次精度スキームに変え、発散の「タネ」となる数値振動を抑えることが効果的です。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——個数密度関数法（PBM）の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告