共役熱伝達（CHT）

大ありだよ。たとえばタービン翼の冷却設計では、翼内部の冷却通路を流れる冷却空気と、翼外面を通過する高温燃焼ガスが、翼金属を介して熱的に結合している。固体側の温度分布を正確に得るには、流体側の熱伝達係数 $h$ を仮定するのではなく、流体の温度場と固体の熱伝導を同時に解く必要があるんだ。

もちろん。ヒートシンクとファン冷却の組合せ、パワー半導体モジュールの冷却板設計、LEDパッケージの放熱設計など、固体の熱伝導と流体の対流伝熱が密に結合する問題はすべてCHTの出番だ。

流体領域ではNavier-Stokes方程式とエネルギー方程式を解く。固体領域では熱伝導方程式を解く。そしてこの2つの領域の界面で温度と熱流束の連続条件を課すんだ。

界面条件を数式で書くとこうなる。

そのとおり。固体側の熱伝導方程式は定常の場合

理想的な界面なら上記でよい。実際のTIM（Thermal Interface Material）やボルト締結面では接触熱抵抗を界面条件に追加する必要がある。Ansys FluentやSTAR-CCM+ではthin wallやcontact resistanceとして設定できるよ。

大きく2つある。一体型（monolithic）と分離型（partitioned）だ。一体型は固体と流体を1つのソルバーで同時に解く。Ansys Fluent、STAR-CCM+、OpenFOAMのchtMultiRegionFoamがこの方式だ。

分離型は固体ソルバーと流体ソルバーを別々に走らせて、界面データを交換する。たとえばAnsys MechanicalとFluent間、あるいはAbaqusとSTAR-CCM+間のco-simulationがこれにあたる。FSI（流体構造連成）でも使われる手法だね。

一体型は界面の整合性が高く収束が早い。分離型は既存ソルバーを組合せられる柔軟性がある反面、界面補間の精度や収束の安定性に注意が必要だ。

大いになる。固体側は比較的粗くて済むことが多いけど、流体側は壁面境界層を解像する必要がある。Ansys Fluentの場合、壁面第一層の $y^+$ を1程度に設定して壁面関数を使わない（低Re乱流モデル）か、$y^+ \approx 30$ で壁面関数を使うかを選択する。

そのとおり。特に局所的な熱伝達係数やNu数分布を議論するならprism layer（インフレーション層）を十分な層数配置すべきだ。STAR-CCM+なら壁面からの総厚さと層数を指定して自動生成できる。OpenFOAMではsnappyHexMeshのaddLayersControlで設定する。

残差だけでなく、界面での温度と熱流束のモニタリングが重要だ。界面の平均温度や最高温度が反復間で変動しなくなったことを確認する。Fluentならsurface monitorで界面の面積加重平均温度を追跡するのが定石だよ。

一体型なら通常のCFD計算と同程度。分離型co-simulationでは各ステップ内でのサブイテレーション回数を3〜10程度に設定し、界面値の変動が十分小さくなることを確認する。緩和係数（under-relaxation）の調整もポイントだ。

典型的な手順はこうだ。(1) CADから固体領域と流体領域を定義する。(2) 流体領域のメッシュを生成し、壁面境界層をprism layerで解像する。(3) 固体領域のメッシュを生成する。(4) 界面をconformal（節点一致）またはnon-conformal（非一致）で接続する。(5) 材料物性値を設定する。(6) 境界条件と初期条件を設定して計算を実行する。

conformalは界面の節点を完全に一致させる方式で、補間誤差がない。ただしメッシュ生成の制約が大きい。non-conformalは界面の節点が一致しなくてもよく、メッシュ生成が容易だが、界面での補間による微小なエネルギー不整合が生じうる。実務ではnon-conformalが圧倒的に多いよ。

古典的なベンチマークとしては、平板上のCHTで解析解と比較する方法がある。流体側のNu数相関式（たとえば $Nu = 0.332 Re^{1/2} Pr^{1/3}$ ）と組合せて固体温度分布を理論的に予測し、CFD結果と照合する。

また、ECSのようなラウンドロビンテストの実験データも利用できる。具体的には、矩形チャネル内のピン付き壁面CHT問題は複数の研究グループが実験結果を公開しているので、検証に適している。

3水準以上のメッシュ密度（粗・中・細）でRichardson外挿を行い、Grid Convergence Index（GCI）を算出する。界面温度と界面熱流束の両方でGCIを評価するのが理想的だ。

一番多いのは固体と流体の材料物性値の単位系不整合だ。特にOpenFOAMでは全物性をSI基本単位で入力する必要がある。次に多いのが界面設定の不備で、界面としてマッピングされていないfaceが存在すると、そこが断熱壁になってしまう。STAR-CCM+ではpart contactの設定漏れが典型例だよ。

固体と流体の熱容量比が極端に大きい場合（たとえば銅とガスの組合せ）、エネルギー方程式の緩和係数を下げる必要がある。Fluentならunder-relaxation factorを0.8〜0.9程度に調整する。それでもダメなら擬似非定常法（pseudo transient）を使うと安定することが多い。

かなり違う。主要なソルバーを比較してみよう。

たしかにディレクトリ構造が複雑で、各regionごとにsystem/、constant/、0/を用意する必要がある。ただし一度テンプレートを作ってしまえば使い回しが利く。OpenFOAM v2306以降はチュートリアルが充実しているので、まずはheatTransfer/chtMultiRegionFoam/multiRegionHeaterを動かすのが早道だ。

ある。代表的なのがAnsys System Couplingで、Fluent（流体）とMechanical（固体）を接続する。Siemens側ではSTAR-CCM+とAbaqusのco-simulationが可能だ。この方式は固体側で熱応力まで一貫して解けるのが利点だよ。

一体型CHTは流体のみの計算の1.1〜1.5倍程度。固体領域は自由度が追加されるけど、計算量の大半は流体側で決まるからね。Co-simulation型は通信オーバーヘッドがあるぶん1.5〜3倍程度になることが多い。

Ansys Fluent 2024R1以降はGPUソルバーがCHTに対応している。STAR-CCM+もGPU accelerationでCHTを扱える。OpenFOAMはAmgXやPETScバックエンドでGPU解法を使えるけど、マルチリージョンでの安定性は発展途上だ。

あるよ。固体と流体で時間スケールが大きく異なる場合、非定常CHTが必要になる。たとえばディーゼルエンジンのシリンダヘッドでは、燃焼サイクル（ミリ秒オーダー）と金属の温度変動（秒〜分オーダー）が共存する。

このようなマルチスケール問題では、流体側と固体側で異なる時間刻みを使うdual time steppingが有効だ。STAR-CCM+のImplicit Unsteady CHT機能やFluentのDual Time Stepping機能がこれに対応している。

流体側はCFL条件やターンオーバー時間で決まるけど、固体側は固体のFourier数 $Fo = \alpha \Delta t / L^2$ で制約される。典型的には固体側の $\Delta t$ は流体側の10〜1000倍大きくできる。超周期（super cycling）テクニックとして、流体を何ステップか進めてから固体を1ステップ進める手法もある。

代表的な応用を挙げよう。

そのとおり。翼表面温度を10度下げるだけで翼寿命が2倍になるとも言われる。CHTの精度が設計に直結する典型例だよ。

鋭いね。1000K以上の環境ではStefan-Boltzmann則 $q_{rad} = \varepsilon \sigma (T^4 - T_{surr}^4)$ に基づく表面間輻射やガス輻射も重要になる。FluentのS2SモデルやDOモデル、STAR-CCM+のSurface-to-Surface Radiation、OpenFOAMのviewFactorモデルなどをCHTと組合せることで、対流・伝導・輻射の三位一体の解析が可能だ。

それは界面マッピングの設定ミスが最も疑わしい。まず界面がちゃんとcoupled wallとして認識されているか確認しよう。Fluentなら境界条件でwallタイプがcoupledになっているか、STAR-CCM+ならContact Interfaceが正しく生成されているかをチェック。

non-conformal interfaceで、マッピング先が見つからないfaceが存在する場合がある。Fluentのmesh interfaceでorphan faceが報告されていないか確認。STAR-CCM+ではContact設定のTolerance値が小さすぎて界面が部分的にしか接続されていないケースがある。

CHTでよくある症状だ。原因と対策を整理しよう。

FluentのPseudo Transient設定は定常CHTの収束改善にかなり効果的だ。Time Scale FactorをAutoにしておけば、ソルバーが自動的に適切な擬似時間刻みを選んでくれる。特に発熱量が大きい問題で有効だよ。

まず乱流モデルの妥当性を疑おう。標準k-εは衝突噴流領域や再付着点付近で熱伝達を過小評価する傾向がある。SST k-ωモデルやRealizable k-εモデルに切り替えてみよう。$y^+$が十分小さい（1以下）ことも再確認すること。

固体側の熱伝導率が正しいか、界面の接触熱抵抗が考慮されているか、流体の物性値（特にPrandtl数に効く粘度と熱伝導率）が正確かを順にチェックする。また、2D軸対称で解いている場合に実は3D効果が重要なケースもある。最終手段としてLES（Large Eddy Simulation）を検討する価値があるよ。