自然対流のCFD解析
理論と物理
自然対流の基礎
先生、自然対流って外部動力なしで流れが起きるんですよね。何が駆動力なんですか？
浮力だよ。流体の密度が温度によって変化し、重力場の中で密度差が圧力不均衡を生むことで流れが発生する。高温部の流体は低密度になって上昇し、低温部の流体が下降する。この循環が自然対流だ。
自然対流の強さを表す指標は？
Rayleigh数 $Ra$ が基本指標だ。
ここで $g$ は重力加速度、$\beta$ は体膨張係数、$\Delta T$ は温度差、$L$ は代表長さ、$\nu$ は動粘度、$\alpha$ は温度拡散率。Grashof数 $Gr$ とPrandtl数 $Pr$ の積だ。
$Ra$ によって流れのレジームが分かれる。垂直平板の場合、$Ra < 10^9$ が層流、$Ra > 10^9$ が乱流域。密閉キャビティでは $Ra > 10^8$ 程度で非定常流れに遷移する。
Boussinesq近似
Boussinesq近似って具体的にどういう仮定ですか？
密度変動が参照密度に比べて十分小さい場合（$\beta \Delta T \ll 1$）に使える近似で、運動方程式の浮力項にのみ密度変動を考慮し、連続の式では密度を一定と扱う。
$$ \rho \approx \rho_0 [1 - \beta(T - T_0)] $$
Boussinesq近似って具体的にどういう仮定ですか？
密度変動が参照密度に比べて十分小さい場合（$\beta \Delta T \ll 1$）に使える近似で、運動方程式の浮力項にのみ密度変動を考慮し、連続の式では密度を一定と扱う。
空気なら $\Delta T < 30$度C程度が安全な適用範囲。水なら体膨張係数が温度依存性を持つので、広い温度範囲ではpolynomial密度を使うべきだよ。
Rayleigh数10⁹を超えると「乱流自然対流」の世界
自然対流はRayleigh数（Ra = Gr × Pr）によって流れのレジームが劇的に変わる。Ra ＜ 10⁸では安定した層流対流セルが形成され、Nusselt数の相関式も比較的シンプルだ。しかしRa ＞ 10⁹になると乱流化が始まり、Nu数の予測精度が急落する。これは「太陽光で熱せられた建物外壁」や「大型変圧器の油冷却」に相当するスケールで、現場では普通に出てくる条件だ。乱流自然対流では壁面近傍の羽状構造（プルーム）がランダムに発生・消滅するため、定常解析が収束しない場合も多い。そういう場合は非定常計算で時間平均を取る戦略に切り替えると良い。
各項の物理的意味
- 時間項 $\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$。時間刻みの安定性に直結 |
数値解法と実装
自然対流の乱流モデリング
自然対流では乱流モデルの選択が難しいと聞きました。
自然対流は強制対流と比べて乱流強度が低く、遷移領域が広い。標準k-εモデルは過剰な乱流拡散を生じやすく、Nu数を過大評価する傾向がある。推奨は以下の順だ。
| Rayleigh数範囲 | 推奨モデル | 備考 |
|---|---|---|
| $Ra < 10^9$ | Laminar（層流モデル） | 乱流モデル不要 |
| $10^9 < Ra < 10^{12}$ | SST k-ω + Low-Re壁面処理 | 壁面第一層 $y^+ < 1$ 必須 |
| $Ra > 10^{12}$ | SST k-ω or Realizable k-ε | Enhanced Wall Treatment |
自然対流で壁面関数は使えないんですか？
基本的に避けるべきだ。自然対流の壁面近傍の速度・温度プロファイルは強制対流の対数則とは異なるため、標準的な壁面関数の適用性が低い。$y^+ \approx 1$ のLow-Re approachが推奨される。Fluentの Enhanced Wall Treatment は $y^+$ に応じて自動切り替えするが、自然対流では $y^+ < 1$ を確保するのがベストプラクティスだ。
メッシュ設計
自然対流のメッシュ設計で気をつけることは？
温度境界層と速度境界層の厚さを見積もって、それぞれに十分な数のセルを配置する。自然対流の境界層厚さは
で概算できる。たとえば $L = 0.1$m、$Ra = 10^9$ の垂直平板なら $\delta_T \approx 0.6$mm。この薄い境界層内に少なくとも10〜15層のセルを配置する必要がある。
境界層の外側は粗くていいんですか？
コア領域は比較的粗くて構わない。ただし急激なセルサイズ変化は避けるべきで、隣接セルの体積比を1.2以下に抑えるのが目安だ。STAR-CCM+のPrism Layer MesherやFluentのInflation Layerで壁面近傍を自動的に細分化できるよ。
自然対流CFDのBoussinesq近似——いつ使えて、いつ破綻するか
自然対流解析の定番「Boussinesq近似」は、密度変化を浮力項だけに線形近似（ρ ≈ ρ₀(1-βΔT)）し、他の項は一定密度として扱う手法だ。計算が安定し収束しやすいが、温度差ΔTが大きくなると誤差が無視できなくなる。目安として「βΔT < 0.1（約10〜20℃差まで）」が適用限界の経験則とされる。炉内燃焼や太陽熱集熱器など数百℃の温度差が生じる系では、密度を温度の完全関数として扱う「非Boussinesq（variable density）」モデルが必須。PressureベースのNavier-Stokesソルバで低Mach数圧縮性を考慮するか、密度依存の状態方程式を直接組み込む必要がある。
風上差分（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次だが、流れの方向を正しく捕捉するため安定性が高い。
実践ガイド
密閉キャビティのベンチマーク
自然対流CFDの検証に使えるベンチマークはありますか？
最も有名なのがde Vahl Davis（1983）の差分加熱矩形キャビティ問題だ。左壁が高温、右壁が低温、上下壁が断熱で、$Ra = 10^3$〜$10^6$ の範囲でNu数と流れ場の参照解が与えられている。CFDコードの検証では、壁面平均Nu数がde Vahl Davisの値と1%以内で一致することを確認するのが標準だ。
具体的にどんな値になりますか？
| $Ra$ | 平均 $Nu$ (de Vahl Davis) |
|---|---|
| $10^3$ | 1.118 |
| $10^4$ | 2.243 |
| $10^5$ | 4.519 |
| $10^6$ | 8.800 |
CFDの結果がこれらと2%以上ずれていたら設定やメッシュに問題がある可能性が高い。
電子機器の自然空冷設計
ファンなしの自然空冷ではどんなCFD設計をするんですか？
密閉筐体内の自然対流CFDでは、(1) 各部品の発熱量を体積ソースとして設定、(2) 筐体壁面の外側に自然対流+輻射の境界条件を設定、(3) 筐体内部の空気を流体として解く、という手順になる。多くの場合、輻射の寄与が全放熱量の30〜50%に達するので、Surface-to-Surface輻射モデル（FluentのS2S、STAR-CCM+のSurface-to-Surface Radiation）の併用が必須だ。
基板の熱伝導をモデル化するのは大変ですよね。
PCB（プリント基板）は銅層とガラスエポキシ層の積層構造で、面内方向と厚さ方向で熱伝導率が10倍以上異なる。Ansys Icepakなど電子機器熱設計専用のツールではPCBの自動異方性モデル化機能がある。汎用CFDソルバーではorthotropic thermal conductivityを手動設定する必要があるよ。
通風孔がある筐体では？
通風孔（vent）がある場合は外気との空気交換が起こるので、筐体外部の空間もモデル化する必要がある。通風孔は多数の小さな開口の集合なので、直接メッシュを切ると膨大なセル数になる。Porous jump（多孔板モデル）で開口率と圧損係数を指定して模擬するのが実用的だ。
電子機器の「密閉筐体CFD」が難しいわけ
電子機器の密閉筐体内の自然対流解析は、実務で最も計算コストがかかる部類に入る。理由は「細長い空間に弱い浮力流」という最悪の組み合わせだからだ。アスペクト比の大きい密閉空間では、解の収束が非常に遅く、定常解析でも数千ステップかかることがある。実際の設計現場では、まず簡略化した2Dモデルで流れのパターンを確認し、発熱源の配置が決まった段階で3Dの詳細計算に移行するという「2段階アプローチ」が採用されることが多い。また、密閉空間では初期条件の設定（基板温度の初期値など）が収束性に大きく影響するため、まず静止流体状態から加熱して流れを起こすシーケンス的な初期化が有効だ。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ CFDでも入口条件のミスは下流全体に波及します。
ソフトウェア比較
Nu数の相関式
CFD結果を検証するためのNu数相関式を教えてください。
代表的なものを挙げよう。
垂直平板（層流）:
垂直平板（全域） (Churchill & Chu):
水平加熱板（上面）: $Nu_L = 0.54 Ra_L^{1/4}$ ($10^4 < Ra < 10^7$)
これらの相関式とCFDの乖離はどの程度が普通ですか？
適切に設定されたCFDなら5〜10%以内。20%以上ずれているなら設定やメッシュに問題がある。
ソルバー別の設定ポイント
各ソルバーで特に気をつけるべき点は？
Fluent: Pressure-Based Coupled solver推奨。Body Force WeightedのPressure Interpolationが浮力流で安定性を改善する。Pseudo Transientも有効。
STAR-CCM+: Segregated Flow + Gravity Model ON。Reference ValuesのReference Temperatureを流体の平均温度に設定。Boussinesq Modelの場合、Thermal Expansion Coefficientの入力が必要。
OpenFOAM: buoyantBoussinesqSimpleFoam（Boussinesq近似、定常）が最も簡単。非定常の場合はbuoyantBoussinesqPimpleFoam。constant/transportPropertiesでbeta（体膨張係数）とTRef（参照温度）を設定する。
OpenFOAMのbuoyantBoussinesqSimpleFoamは初心者にも使いやすいですか？
Galerkin FEM系のソルバーより設定項目が明示的なので、何をやっているかは分かりやすい。ただしメッシュ生成（blockMeshやsnappyHexMesh）のハードルが高い。まずはtutorials/heatTransfer/buoyantBoussinesqSimpleFoam/hotRoomを動かしてみるのがおすすめだ。
自然対流ソルバー——「収束性」がツール選定の決め手になる
自然対流の解析ではソルバーの収束性がツール選定に直結する。浮力と圧力のカップリングは強制対流より不安定になりやすく、SIMPLE系アルゴリズムでは緩和係数の調整が欠かせない。Fluent は圧力緩和係数0.3、運動量0.7が初期値だが、自然対流では圧力0.1〜0.2に下げないとflopする案件が多い。STAR-CCM+のCOUPLEDソルバーは自然対流との相性が良く、緩和係数の手動調整なしで安定するケースが多い。OpenFOAMはbuoyantBoussinesqSimpleとbuoyantBoussinesqPisoの使い分けが重要で、定常解析でも収束しない場合はPisoFoam（非定常）に切り替えるのが鉄則だ。「どのソルバーが自然対流で安定するか」は実務ナレッジとして社内に蓄積しておく価値がある。
選定で最も重要な3つの問い
- 「何を解くか」：自然対流のCFD解析に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
高Ra数の乱流自然対流
Rayleigh数が非常に高い場合（$Ra > 10^{12}$）のCFDはどうなりますか？
大規模建築空間（アトリウム、倉庫）、原子炉格納容器内の自然対流、マントル対流などがこの領域だ。流れは完全に乱流で、大規模な渦構造（plume、thermal）が支配する。RANSでは定性的な傾向はつかめるが、定量的な予測にはLESが望ましい。
Rayleigh-Benard対流もCFDで解けますか？
解ける。水平面を下から加熱する古典的なRayleigh-Benard対流は、$Ra > 1708$ でセル構造が現れ、$Ra > 10^7$ 程度で乱流化する。DNSデータが豊富に存在するので（Verziccoらのデータベースなど）、CFDコードの検証に最適だ。Nu数のスケーリング則 $Nu \sim Ra^{0.31}$ が $10^7 < Ra < 10^{14}$ で成り立つことが知られている。
DNS/LESデータの活用
DNSデータをRANSモデルの改良に使う研究もあるんですか？
ある。自然対流の壁面近傍の乱流統計量（レイノルズ応力、乱流熱流束）のDNSデータを使って、RANSモデルの浮力項の係数を最適化する研究が進んでいる。PINNs（Physics-Informed Neural Networks）を使って乱流モデルの補正項を学習させるアプローチもある。
輻射との連成
自然対流では輻射の影響が大きいんですよね？
常温の室内環境でも壁面間の輻射が全熱伝達の30〜50%を占める。高温環境（炉、高温配管）ではさらに支配的になる。CFDではSurface-to-Surface（S2S）輻射モデルが標準で、面間の形態係数（view factor）を計算して輻射熱交換を求める。
透明媒体（空気など）ならS2Sで十分だが、煙やガスを含む場合はParticipating media（参加性媒体）モデルが必要。Discrete Ordinates（DO）モデルやP1モデルが使われる。FluentではRadiation ModelsからDOやS2Sを選択できるよ。
Coffee Break よもやま話
高Ra数自然対流——原子炉冷却の最終防衛線
原子炉の炉心溶融事故シナリオで、最後の頼みの綱になるのが「自然対流による受動冷却」だ。外部電源が全て失われても、溶融した燃料から発生する崩壊熱を自然対流だけで除去するという設計思想で、Ra数は10¹²〜10¹⁴という超高Ra数の領域になる。この領域ではLESでも計算コストが天文学的になるため、スケールモデル実験とRANSモデルの組み合わせ、さらにBEACH（Boussinesq-based Entropy-Augmented CHT）などの特殊なモデルが研究されている。福島第一原発事故後、この分野の研究投資が世界規模で急増し、高Ra数CFDの精度向上が安全解析の最重要課題のひとつになっている。
DNSデータをRANSモデルの改良に使う研究もあるんですか？
ある。自然対流の壁面近傍の乱流統計量（レイノルズ応力、乱流熱流束）のDNSデータを使って、RANSモデルの浮力項の係数を最適化する研究が進んでいる。PINNs（Physics-Informed Neural Networks）を使って乱流モデルの補正項を学習させるアプローチもある。
自然対流では輻射の影響が大きいんですよね？
常温の室内環境でも壁面間の輻射が全熱伝達の30〜50%を占める。高温環境（炉、高温配管）ではさらに支配的になる。CFDではSurface-to-Surface（S2S）輻射モデルが標準で、面間の形態係数（view factor）を計算して輻射熱交換を求める。
透明媒体（空気など）ならS2Sで十分だが、煙やガスを含む場合はParticipating media（参加性媒体）モデルが必要。Discrete Ordinates（DO）モデルやP1モデルが使われる。FluentではRadiation ModelsからDOやS2Sを選択できるよ。
高Ra数自然対流——原子炉冷却の最終防衛線
原子炉の炉心溶融事故シナリオで、最後の頼みの綱になるのが「自然対流による受動冷却」だ。外部電源が全て失われても、溶融した燃料から発生する崩壊熱を自然対流だけで除去するという設計思想で、Ra数は10¹²〜10¹⁴という超高Ra数の領域になる。この領域ではLESでも計算コストが天文学的になるため、スケールモデル実験とRANSモデルの組み合わせ、さらにBEACH（Boussinesq-based Entropy-Augmented CHT）などの特殊なモデルが研究されている。福島第一原発事故後、この分野の研究投資が世界規模で急増し、高Ra数CFDの精度向上が安全解析の最重要課題のひとつになっている。
トラブルシューティング
定常計算が収束しない
自然対流の定常計算で残差が振動して収束しません。
自然対流では物理的に非定常な流れが発生することが多い。残差振動は「定常解が存在しない」ことを意味している場合がある。$Ra > 10^8$（密閉キャビティ）や $Ra > 10^9$（垂直平板）では非定常計算に切り替えるべきだ。
定常計算のまま収束させるテクニックはありますか？
(1) FluentのPseudo Transient設定を有効にする。(2) under-relaxation factorsを大幅に下げる（圧力0.2、運動量0.3程度）。(3) First Order Upwindで初期収束させてからSecond Orderに切り替える。ただしこれらは人工的に定常化しているだけなので、物理的妥当性の確認が必要だよ。
Nu数がde Vahl Davisのベンチマーク値と合わない
差分加熱キャビティで $Ra = 10^6$ のNu数がベンチマークの8.8ではなく6.5になります。
メッシュが粗すぎる可能性が高い。$Ra = 10^6$ では壁面近傍の境界層が非常に薄く、不均一メッシュでの壁面第一層が厚すぎると熱流束を過小評価する。壁面法線方向に最低20層、壁面第一層厚さが境界層の1/10以下になるようにメッシュを調整しよう。
メッシュを細かくしたら8.6になりました。まだ2%ずれています。
Second Order Upwindを使っているか確認。First Orderだと数値拡散でNu数が低く出る。また、Boussinesqモデルの参照温度 $T_0$ が高温壁と低温壁の平均温度に設定されているか確認。参照温度の不適切な設定はNu数に影響する。
輻射モデルを入れると結果が大きく変わる
輻射モデルを追加したらNu数が50%増加しました。妥当ですか？
壁面放射率（emissivity）の設定を確認しよう。金属光沢面なら $\varepsilon \approx 0.1$、塗装面なら $\varepsilon \approx 0.9$。デフォルトの $\varepsilon = 1.0$（完全黒体）のままだと輻射を過大評価する。実際の壁面材料に合わせて設定すべきだ。
実験と比較する場合、輻射は入れるべきですか？
実験データが対流成分のみのNu数を報告しているのか、輻射込みの全熱伝達を報告しているのかを確認する必要がある。多くの実験では輻射を分離できていないので、CFDでも輻射を含めて比較するのが妥当だ。
自然対流CFDが収束しない——浮力主導流れの典型的な罠
自然対流CFDは「収束が遅い」「残差が振動する」という悩みが絶えない。主因は浮力と流れの強い非線形結合にある。典型的な失敗例：①初期条件を静止状態（ゼロ速度）から始めると浮力で発生した流れが非現実的に大きくなり発散する——対策は緩和係数を0.1〜0.3程度に絞った段階的起動。②粗いメッシュで定常解析を試みる——自然対流は本質的に非定常な振動解の場合があり（Ra>10⁸で乱流移行）、定常仮定自体が成立しないことがある。③重力方向の設定ミス——ANSYSではデフォルト重力がゼロで、浮力が全く働かない設定になりがち。起動時のgベクトル設定は必ず確認しよう。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——自然対流のCFD解析の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告