沸騰モデル — 実践ガイドとベストプラクティス
実践ガイド
沸騰解析を実務でやるときの手順を教えてください。
典型的なサブクール核沸騰(管内上向き流)の解析フローを示そう。
1. 形状作成: 管径、加熱長、入口助走区間を含めたモデル
2. メッシュ: 壁面近傍プリズム層(第一セル高さ ≈ 気泡離脱径)、管断面方向に20セル以上
3. 物性値設定: 飽和温度、$h_{fg}$、飽和圧力での液相・気相物性
4. 単相助走計算: まずヒーターなしで完全発達流を得る
5. 熱流束印加: 壁面均一熱流束を段階的に増加
6. 沸騰モデル有効化: RPIモデルの各パラメータを設定
7. モニタリング: 壁面温度、ボイド率分布、出口気泡径
メッシュ設計
沸騰解析のメッシュで特に気をつけるべきことは?
一般的なCFDメッシュとは異なるポイントがある。
| 項目 | 推奨 | 理由 |
|---|---|---|
| 第一セル高さ | 気泡離脱径 $d_w$ 以上 | RPIモデルの前提条件 |
| 壁面方向の成長率 | 1.1〜1.15 | ボイド分布の解像 |
| 軸方向要素サイズ | 管径の1/10〜1/5 | 沸騰開始点の捕捉 |
| 周方向分割 | 均等40分割以上 | 非対称ボイドパターンの解像 |
2D軸対称で計算してもいいですか?
低熱流束のサブクール沸騰なら2D軸対称でも合理的だ。ただしCHF近傍や膜沸騰に近づく場合は3D非対称な挙動が重要になるので、3D計算が必要になる。
物性値の設定
沸騰解析の物性値で注意すべき点はありますか?
最も重要なのは飽和温度 $T_{sat}$ と潜熱 $h_{fg}$ の圧力依存性だ。系圧が変化する場合、局所圧力に応じた飽和物性を使う必要がある。
水の場合、IAPWS-IF97(蒸気表の国際標準規格)の物性データを使うのが標準だ。Fluentでは組み込みの水蒸気物性テーブルがある。STAR-CCM+でもIAPWSベースの物性ライブラリが利用可能だ。
検証に使えるベンチマーク
沸騰解析の結果を検証するための実験データはありますか?
代表的なベンチマーク実験を挙げよう。
| 実験者 | 条件 | 計測量 |
|---|---|---|
| Bartolomej (1967) | 管内サブクール沸騰 | 断面ボイド率分布 |
| Lee & Mudawwar (1988) | 矩形流路沸騰 | 壁面温度、CHF |
| DEBORA実験 (CEA) | R-12冷媒沸騰 | 気泡径、ボイド率、速度 |
| PSBT (NRC) | PWR燃料集合体模擬 | サブチャンネルボイド率 |
DEBORA実験は冷媒を使っているんですね。
R-12は水より蒸気密度が高く密度比が小さいので、スケーリングパラメータで水の挙動を模擬できる。低圧・低温で実験できるので計測しやすいという利点があるんだ。
レイノルズの実験(1883年)——乱流発見の瞬間
オズボーン・レイノルズは、管内の水にインクを流す実験で「層流から乱流への遷移」を発見しました。流速を上げていくと、インクの線がある瞬間にグチャグチャに乱れる。この劇的な瞬間を、レイノルズは数学的に $Re = \rho uD/\mu$ という無次元数で表現した。100年以上経った今も、CFDエンジニアが最初に確認するのはこのレイノルズ数です。
実務者のための直感的理解
この解析分野のイメージ
CFDって、要は「デジタル風洞」です。自動車メーカーが巨大な風洞実験設備に何億円もかけるところを、PCの中で再現できる。でも1つ注意——風洞実験なら「風を当てれば結果が出る」けど、CFDでは「メッシュの品質」と「乱流モデルの選択」という見えない品質要因がある。ここを手抜きすると、きれいなコンター図が出ても中身はデタラメ…なんてことになりかねません。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め(計算領域)、水の入り口と出口を設計し(境界条件)、ポンプの強さを設定する(流量条件)。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか? 意外と直感的ではありませんか?
初心者が陥りやすい落とし穴
「y+って何ですか?」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない!」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た!」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか(低速)、全開にするか(高速)。でもCFDではもう一つ——「どのくらい暴れた水を出すか」(乱流強度)も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね? CFDでも入口条件のミスは下流全体に波及します。
CFDメッシュの品質管理や乱流モデルの選定に悩む時間を、もっと創造的な設計作業に使えたら。 — Project NovaSolverはそんな実務者の声から生まれました。
Project NovaSolver — CAE実務の課題に向き合う研究開発
「沸騰モデルをもっと効率的に解析できないか?」——私たちは実務者の声に耳を傾け、既存ワークフローの改善を目指す次世代CAEプロジェクトに取り組んでいます。具体的な機能はまだ公開前ですが、開発の進捗をお届けします。
進捗通知を受け取る →