実在気体効果 — 実践ガイドとベストプラクティス
超臨界CO₂タービンの解析
超臨界CO₂って発電システムに使われるんですよね？ CFDではどう扱うんですか？
超臨界CO₂ブレイトンサイクルは次世代発電の有望技術だ。圧縮機入口が臨界点近傍（31°C, 7.4 MPa）で密度が液体並みに高いから、装置を小型化できる。CFDで扱う際の注意点を整理しよう。
| 項目 | 設定 | 備考 |
|---|---|---|
| EOS | Span-Wagner（CO₂専用高精度EOS） | PR-EOSでは精度不足の場合あり |
| 物性テーブル | REFPROP or CoolProp生成 | 温度200-800K, 圧力1-30 MPa |
| 乱流モデル | SST k-omega + 低Re補正 | 密度急変で乱流輸送が変化 |
| メッシュ | y⁺ < 1、擬臨界線周りに集中 | 壁面近傍の密度変化を解像 |
| ソルバー | Pressure-Based Coupled | 超臨界では密度ベースより安定 |
Span-Wagner EOSは何が違うんですか？
CO₂のためにフィッティングされたHelmholtz自由エネルギー型EOSで、臨界点近傍でもPR-EOSより格段に精度が高い。NISTのREFPROPに実装されている。FluentではNIST Real Gas Property（RGP）テーブル形式でインポートする。
LNGプロセスの実在気体効果
LNG（液化天然ガス）のバルブやパイプのCFDでも実在気体が重要ですか？
LNG条件（メタン主成分、-162°C、大気圧で液化）ではPR-EOS+混合則で多成分系を扱う。BOG（Boil Off Gas）の発生やフラッシュ蒸発の予測には気液平衡計算が必要で、通常のCFDの枠を超える。
FluentでLNGの気液二相流は計算できますか？
FluentのEvaporation-Condensationモデル（Lee model）とPR-EOSの組み合わせで概算できるが、多成分LNGの相平衡を正確に扱うにはCFXの方が強い。CFXはmulticomponent multiphaseで相平衡計算を内蔵しているからだ。
高温気体の検証事例
高温実在気体の検証にはどんな実験データを使いますか？
代表的なベンチマークを挙げよう。
- 正規衝撃波の温度: 衝撃波前後の温度・密度比を理想気体/平衡/非平衡モデルで比較（Andersonの教科書に詳細データあり）
- 衝撃波管実験: 既知のドライバー/テスト条件から衝撃波速度と後方状態を検証
- 円柱周りの弓状衝撃波: M=6-10でのスタンドオフ距離が気体モデルに強く依存
- NASA CEA: 化学平衡組成をCFDの結果と照合
理想気体と実在気体で結果がどれくらい違うか、定量的に見たいです。
例えばM=10の垂直衝撃波を考えよう。理想気体（$\gamma=1.4$）では $T_2/T_1 \approx 20.4$, $\rho_2/\rho_1 = 5.71$ だが、解離を含む実在気体では温度が下がり密度比は10以上になる。衝撃層が薄くなるのは密度が上がるからで、これは再突入体の衝撃波スタンドオフ距離に直結する。
実務チェックリスト
実在気体CFDを始めるときのチェックリストをまとめてもらえますか？
以下の項目を確認しよう。
- 作動流体の臨界点（$T_c, p_c$）と動作条件の距離は十分か → 臨界点近傍ならテーブル法推奨
- EOSの精度は検証済みか → 既知の状態点でREFPROPと比較
- 比熱の温度依存性は正しく入力されているか → NASAポリノミアル（7係数）の温度範囲を確認
- 輸送係数（粘性、熱伝導率）もRGPテーブルに含まれているか
- メッシュは物性急変領域で十分に細かいか
ライト兄弟は最初の「CFDエンジニア」だった？
ライト兄弟は1901年に自作の風洞で200以上の翼型を試験しました。当時のコンピュータは？ もちろん存在しません。彼らは手作業で揚力と抗力を測定し、最適な翼型を見つけ出した。現代のCFDエンジニアがFluent1発で計算する揚力係数を、ライト兄弟は何百回もの風洞実験で手に入れたのです。
実務者のための直感的理解
この解析分野のイメージ
CFDって、要は「デジタル風洞」です。自動車メーカーが巨大な風洞実験設備に何億円もかけるところを、PCの中で再現できる。でも1つ注意——風洞実験なら「風を当てれば結果が出る」けど、CFDでは「メッシュの品質」と「乱流モデルの選択」という見えない品質要因がある。ここを手抜きすると、きれいなコンター図が出ても中身はデタラメ…なんてことになりかねません。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ CFDでも入口条件のミスは下流全体に波及します。
CFDメッシュの品質管理や乱流モデルの選定に悩む時間を、もっと創造的な設計作業に使えたら。 — Project NovaSolverはそんな実務者の声から生まれました。
次世代CAEプロジェクト：開発者と実務者をつなぐ
Project NovaSolverは、実在気体効果を含む幅広い解析分野において、実務者の知見を最大限に活かせる環境の実現を探求しています。まだ道半ばですが、共に歩んでいただける方を募集しています。
開発パートナー登録 →