混合対流

そのとおり。外部からの流れ（ファンやポンプ）と浮力による流れが同程度の強さで共存する場合を混合対流（mixed convection）と呼ぶ。たとえば垂直管内の上昇流で管壁を加熱する場合、ポンプによる強制流れと浮力による上昇流が重畳する。

Richardson数 $Ri$ で判断する。

非常に重要だ。加熱垂直管内の上昇流では浮力が流れを助ける方向に働く（助勢流 / aiding flow）。逆に下降流では浮力と流れが逆向きになる（対向流 / opposing flow）。

助勢流ではNu数が純粋な強制対流より増加する。対向流では流れの減速・逆流・再層流化が起こりうるため、Nu数が複雑に変化する。特に対向流の再層流化現象（laminarization）はCFDで正確に予測するのが難しいことで知られているよ。

水平管では浮力が二次流れ（縦渦）を生む。管断面の上部に高温流体が集まり、下部に低温流体が溜まる。この非対称な温度分布を正確に予測するには3D計算が必須で、2D軸対称近似は使えない。

浮力項を正しく扱うために、密度の温度依存性をどうモデル化するかが鍵だ。最も単純なのがBoussinesq近似で、密度を

温度差が大きくて $\beta \Delta T > 0.1$〜$0.2$ 程度になると精度が悪化する。空気の場合、$\Delta T > 30$度C程度で注意が必要だ。この場合は理想気体やpolynomial密度を使って密度の非線形温度依存性を直接扱うべきだ。FluentのIncompressible Ideal Gas設定が便利だよ。

浮力による乱流生成/減衰効果が重要になる。k-ε系モデルでは浮力生成項 $G_b = -g_i \frac{\mu_t}{\rho Pr_t} \frac{\partial \rho}{\partial x_i}$ が追加される。FluentではViscous ModelのOptionsで「Full Buoyancy Effects」をONにすることを強く推奨する。

また、対向流での再層流化現象を予測するにはTransition SSTモデルが有効だ。標準的な乱流モデルでは層流化を予測できず、Nu数を過大評価する。

壁面垂直方向のメッシュ要件は同等（$y^+ \approx 1$ 推奨）だが、浮力による二次流れを解像するために管断面方向のメッシュも十分に細かくする必要がある。水平管の場合、断面方向に少なくとも40〜60分割は必要だ。粗すぎると二次流れが解像できずNu数の非対称性を過小評価する。

非常に重要だ。オフィスの空調では天井からの給気（強制対流）と室内の発熱体（PC、人体、照明）による浮力が共存する。天井吹出し口からの冷気が十分な運動量を持たないと、暖気層（warm layer）が室上部に停滞して温度成層が発生する。

そのとおり。Fluentの場合、定常RANSではRealizable k-ε + Enhanced Wall Treatmentが室内環境CFDの標準選択肢だ。OpenFOAMではbuoyantSimpleFoamにkOmegaSSTを組合せる。人体の発熱は80〜120W程度のheat sourceとしてモデル化し、PMV-PPDによる快適性評価まで含めるのが実務的なワークフローだよ。

ファン付きの筐体では強制対流が主体だけど、ファンなし（自然空冷）やファン故障時は浮力駆動になる。設計上は両方のケースを評価する必要があり、これはまさに混合対流問題だ。

実務では3Dの筐体モデルを作成し、各部品をvolume heat sourceとしてモデル化する。基板はorthotropic（異方性）熱伝導体として面内/面外の熱伝導率を分けて入力する。STAR-CCM+やFluentのCHT機能で固体-流体を同時に解くのが標準的だよ。

鋭い指摘だ。密閉筐体内の自然対流は定常解が存在しない場合がある（Ra数が高いと非定常振動流になる）。定常計算で残差が振動する場合は、非定常計算に切り替えて時間平均を取るべきだ。FluentならTransient設定でadaptive time steppingを使うのが実用的だよ。

垂直管内の混合対流ではJacksonらの実験データが標準参照だ。矩形キャビティの混合対流ではHaidari et al.のデータがある。建築分野ではIEA Annex 20のベンチマークケース（室内環境CFD）が広く使われている。

重要なポイントを順に。(1) Operating Conditions: Operating Densityを適切に設定。Boussinesq近似なら参照温度での密度。(2) Gravity: 重力ベクトルを正しく設定（デフォルトは0）。(3) Material: DensityをBoussinesqまたはIncompressible Ideal Gasに設定。(4) Viscous Model: Full Buoyancy EffectsをON。

Boussinesq近似では圧力の計算から静水圧成分を除去するためにOperating Densityを使う。参照温度での密度に設定するのが標準。これを設定しないと圧力場に大きな静水圧成分が残り、数値精度が悪化する。Incompressible Ideal Gasを使う場合はOperating Densityは不要だ。

Physics ModelsでGravity ModelをONにし、重力ベクトルを指定する。密度はBoussinesq ModelまたはIdeal Gas（非圧縮性近似付き）を選択。乱流モデルでBuoyancy ProductionをONにするのを忘れないこと。Reference ValuesでReference Temperatureを適切に設定するのも重要だ。

buoyantSimpleFoam（定常）またはbuoyantPimpleFoam（非定常）を使う。constant/gファイルで重力加速度を設定。密度はconstant/thermophysicalPropertiesのequationOfStateでBoussinesq（Boussinesq近似）またはperfectGas（理想気体）を指定する。Boussinesqの場合、betaパラメータ（体膨張係数）の入力が必要だ。

浮力がある場合はSIMPLE系でも解けるが、収束が遅いことがある。Pseudo transient（FluentのPseudo Transient設定）やPIMPLE法（OpenFOAMのbuoyantPimpleFoamの定常モード活用）が安定する場合が多い。圧力のunder-relaxationは0.3程度から始めて、収束を見ながら調整しよう。

残差だけでなく、モニタリングポイントの温度と速度が定常に達していることを確認する。浮力流は残差が10のマイナス3乗程度で振動しても、積分量は収束している場合がある。逆に残差が下がっていても積分量が変動しているなら収束していない。

加熱された垂直管内の上昇流（助勢流）で、浮力による加速が壁面近傍の速度プロファイルを一様化させ、乱流エネルギーの生成が抑制されて流れが局所的に層流化する現象だ。Nu数が急激に低下するため、過剰な壁面温度上昇の原因になりうる。

再層流化のパラメータとして $Bo^ = Gr^ / (Re^{3.425} Pr^{0.8})$ が使われ、$Bo^* > 6 \times 10^{-7}$ 程度で層流化が始まるとされている。

通常のk-εやk-ωでは予測できない。Transition SSTモデルを使えばある程度予測可能だが、再層流化の開始位置が実験と10〜30%ずれることがある。正確な予測にはDNS（Direct Numerical Simulation）やWall-Resolved LESが必要だ。

混合対流では浮力プルームの非定常的な揺動や乱流との相互作用が重要で、RANSではこれを時間平均化してしまう。LESでは浮力による大規模渦構造や温度変動を直接解像できるので、非定常熱荷重の評価やThermal fatigue寿命予測に有用だ。

垂直管内の混合対流LESでは、$Re_{\tau}$ に基づく壁面解像のセル数が圧倒的に多くなる。$Re = 10000$ 程度の管流で数千万セル、$Re = 50000$ なら数億セルが必要。HPCクラスターでも1ケース数日〜数週間かかる。WMLESで壁面モデルを使えばセル数を1/10程度に削減できるよ。

原子炉の炉心では冷却材が下から上に流れ、燃料棒の発熱による浮力が流れに重畳する。通常運転時は強制対流支配だが、ポンプ喪失事故（LOCA）時には自然対流に遷移する。この遷移過程が混合対流であり、燃料被覆管の最高温度を正確に予測する安全評価に直結するんだ。NRC（米国原子力規制委員会）の安全解析ではCFDの検証・妥当性確認（V&V）が厳格に要求されるよ。

浮力流は初期条件との乖離が大きいと速度場が急激に変動して発散しやすい。以下の手順を試そう。(1) まずgravityをOFFにして等温流れ場を収束させる。(2) gravityをONにして温度差を段階的に増加させる（FluentならPatch機能で温度を設定）。(3) Pseudo Transientまたは非定常計算で安定化する。

密度のunder-relaxationを追加して0.8程度に設定する（Fluentでは通常1.0がデフォルトでリストに表示されない。Solution > Methods > Under-Relaxation FactorsでDensityを追加する）。Body Forceのunder-relaxationも0.8程度に下げるとよい。

浮力効果がモデルに反映されていない可能性がある。確認ポイント: (1) gravityベクトルの向きが正しいか。(2) DensityがConstantではなくBoussinesqまたはIdeal Gasになっているか。(3) Operating DensityまたはReference Densityが適切に設定されているか。(4) Full Buoyancy EffectsがONか。

非常に低いRi数（$Ri < 0.01$）なら浮力効果がほぼ無視できるので対称は正常。$Ri > 0.1$ で対称なら、初期条件が完全対称になっていて、浮力による摂動が数値的に発達しない可能性がある。わずかな非対称初期摂動（温度場に0.1度程度のランダムノイズ）を加えるか、非定常計算に切り替えるとよい。

浮力流では定常解が存在しない場合がある。まずRa数を確認しよう。$Ra > 10^8$ 程度（密閉キャビティの場合）では流れが本質的に非定常になる。buoyantPimpleFoamに切り替えて非定常計算し、時間平均を取るのが正しいアプローチだ。

浮力流の特性速度 $u_{buoy} = \sqrt{g \beta \Delta T L}$ とセルサイズ $\Delta x$ からCFL数が1以下になるように設定する。adaptive time steppingが使えるならMaxCo 0.5程度で自動調整させるのがよい。OpenFOAMのcontrolDictでadjustTimeStepをyesに設定し、maxCoを指定すればよい。