沸騰モデル — 数値解法と実装
数値解法の詳細
沸騰モデルをCFDに実装するとき、どういうフレームワークで解くんですか？
沸騰解析はEuler-Eulerの二流体モデルをベースに、RPIモデルを壁面境界条件として組み込む形が主流だ。液相と気相それぞれに連続の式と運動量方程式を解く。
気相の体積分率 $\alpha_v$ の輸送方程式には、蒸発・凝縮によるソース項 $\dot{m}$ が入る。
蒸発量はRPIモデルの $q_{evap}$ から決まるんですよね？
その通り。壁面での蒸発質量流束は次の通りだ。
バルク内での凝縮はRanz-Marshall相関から界面熱伝達係数を求め、気泡周囲のサブクール液との熱交換で計算する。
気泡力モデル
沸騰で生成した気泡はどうやって動くんですか？
気泡に作用する力のモデリングが重要だ。相間力として以下を考慮する。
| 力 | モデル | 役割 |
|---|---|---|
| 抗力 | Schiller-Naumann, Ishii-Zuber | 気泡の速度差を支配 |
| 揚力 | Tomiyama | 速度勾配による横方向力 |
| 壁面潤滑力 | Antal | 気泡を壁面から引き離す |
| 乱流分散力 | Lopez de Bertodano | 気泡の乱流拡散 |
| 仮想質量力 | Auton | 加速度効果 |
Tomiyamaの揚力って符号が変わることがあるんですか？
いい質問だ。気泡径がEötvös数 $Eo$ で臨界値を超えると、揚力の方向が反転する。小さな気泡は壁面に向かい、大きな気泡は管中心に移動する。これがボイド分布のwall-peakingとcore-peakingを決める重要な物理だ。
壁面関数の扱い
沸騰面で通常の壁関数は使えるんですか？
使えない。気泡の攪拌により壁面近傍の流れ構造が単相流と大きく異なるからだ。Fluent等ではboiling-specific wall functionが実装されていて、RPIモデルと整合した壁面温度計算を行う。
壁面メッシュは $y^+$ の制約よりも、気泡離脱径 $d_w$ に対して第一セル高さが適切であることが重要だ。一般に最初のセル高さが $d_w$ 以上であることが推奨される。
タイムステップと安定性
沸騰解析で発散しやすいのはなぜですか？
蒸発に伴う体積膨張が急激で、局所的に大きな体積ソース項が発生するからだ。対策としては以下が有効。
- 壁面過熱度を段階的に上げる（ランプアップ）
- 初期は単相定常解を求め、そこから沸騰を有効化
- タイムステップを $10^{-4}$ s以下に設定
- 体積分率方程式にunder-relaxation（0.3〜0.5）を適用
レイノルズの実験（1883年）——乱流発見の瞬間
オズボーン・レイノルズは、管内の水にインクを流す実験で「層流から乱流への遷移」を発見しました。流速を上げていくと、インクの線がある瞬間にグチャグチャに乱れる。この劇的な瞬間を、レイノルズは数学的に $Re = \rho uD/\mu$ という無次元数で表現した。100年以上経った今も、CFDエンジニアが最初に確認するのはこのレイノルズ数です。
離散化手法の詳細解説
空間離散化における手法選択が数値精度・安定性・計算コストに与える影響を詳述する。
風上差分（Upwind）
1次風上: 数値拡散が大きいが安定。2次風上: 精度向上するが振動のリスク。高レイノルズ数流れでは必須。
中心差分（Central Differencing）
2次精度だが、Pe数 > 2で数値振動が発生。低レイノルズ数の拡散支配流れに適する。
TVDスキーム（MUSCL、QUICK等）
リミッタ関数により数値振動を抑制しつつ高精度を維持。衝撃波や急勾配の捕捉に有効。
有限体積法 vs 有限要素法
FVM: 保存則を自然に満足。CFDの主流。FEM: 複雑形状・マルチフィジックスに有利。SPH等のメッシュフリー法も発展中。
マトリクスソルバーの選定指針
問題規模と特性に応じた最適なソルバー選択のガイドライン。
| ソルバー種別 | 詳細・推奨条件 |
|---|---|
| 圧力-速度連成（SIMPLE系） | SIMPLE: 標準的だが収束が遅い。SIMPLEC: 圧力補正の緩和が改善。PISO: 非定常問題に適する。 |
| 連立系ソルバー | AMG（代数的マルチグリッド）: 大規模問題の標準。ILU前処理: メモリ効率良好。ブロックGauss-Seidel: 連成系に有効。 |
| DOF別推奨 | 〜10⁵セル: SIMPLE+AMG、10⁵〜10⁷セル: SIMPLEC+AMG+並列、10⁷セル〜: 結合型ソルバー（Coupled Solver）を検討 |
時間積分法と収束判定
ソルバー内部の制御パラメータと収束判定基準について記述する。
CFL条件（クーラン数）
陽解法: CFL ≤ 1が安定条件。陰解法: CFL > 1でも安定だが、精度と反復回数に影響。LES: CFL ≈ 1を推奨。物理的意味: 1タイムステップで情報が1セル以上進まないこと。
残差モニタリング
連続の式・運動量・エネルギーの各残差が3〜4桁低下で収束と判断。質量保存の残差は特に重要。
緩和係数
圧力: 0.2〜0.3、速度: 0.5〜0.7が一般的な初期値。発散する場合は緩和係数を下げる。収束後は上げて加速。
非定常計算の内部反復
各タイムステップ内で定常解に収束するまで反復。内部反復数: 5〜20回が目安。残差がタイムステップ間で変動する場合は時間刻みを見直す。
数値解法の直感的理解
FVMのイメージ
有限体積法は「会計帳簿」に似ている。各セル（口座）について「入ってくる量」と「出ていく量」の収支を厳密に管理する。隣のセルに流れ出た量は、そのセルに流れ込む量と完全に一致する——これが「保存性」であり、流体解析で質量やエネルギーが勝手に増減しないことを保証する。
SIMPLE法のたとえ
SIMPLE法は「交互に調整する」手法。まず速度を仮に求め（予測ステップ）、その速度で質量保存が満たされるよう圧力を補正し（補正ステップ）、補正された圧力で速度を修正する——このキャッチボールを繰り返して正解に近づく。2人で棚を水平にする作業に似ている：片方が高さを合わせ、もう片方がバランスを取り、これを交互に繰り返す。
風上差分のたとえ
風上差分は「川の流れに立って上流の情報を重視する」手法。川の中にいる人が下流を見ても水の出所は分からない——上流の情報が下流を決めるという物理を反映した離散化手法。精度は1次だが、流れの方向を正しく捕捉するため安定性が高い。
CFDメッシュの品質管理や乱流モデルの選定に悩む時間を、もっと創造的な設計作業に使えたら。 — Project NovaSolverはそんな実務者の声から生まれました。
CAEの未来を、実務者と共に考える
Project NovaSolverは、沸騰モデルにおける実務課題の本質に向き合い、エンジニアリングの現場を支える道具づくりを目指す研究開発プロジェクトです。
プロジェクトの最新情報を見る →