固体・流体界面の取り扱い（CHT解析）
理論と物理
概要
先生、CHT解析で「界面の取り扱い」が大事だって聞いたんですけど、そもそもなぜ界面が問題になるんですか？
いい質問だね。CHT（Conjugate Heat Transfer）解析では、固体の熱伝導と流体の対流を同時に解くんだけど、この2つの領域が接する「界面」で温度と熱流束を正しく受け渡さないといけない。ここがCHT解析の精度を左右する最大のポイントなんだ。
なるほど。具体的にはどういう場面で問題になりますか？
例えばガスタービンブレードの冷却解析だと、ブレード内部（固体）は数mmの要素で十分なのに、外側の高温ガス（流体）は境界層を解像するために数十μmの要素が必要になる。こうなると固体と流体のメッシュが界面で合わなくなる。このとき「適合メッシュ」にするか「非適合メッシュ」にするかで、手間も精度も大きく変わるんだよ。
界面連続条件の数学的定式化
界面で守るべき数学的な条件って、具体的にはどうなっていますか？
界面 $\Gamma$ で満たすべき条件は2つある。まず温度連続条件（Dirichlet条件）：
次に熱流束連続条件（Neumann条件）：
ここで $k_s$, $k_f$ はそれぞれ固体・流体の熱伝導率、$n$ は界面の法線方向だ。この2つの条件を同時に満たすことで、界面でのエネルギー保存が保証される。
温度が等しいだけじゃダメなんですね。熱の出入りも合わせないと。
その通り。温度だけ合わせて熱流束が合わないと、界面からエネルギーが湧き出したり消えたりする「人工的な熱源」が生まれてしまう。実務では、温度分布は妥当に見えるのに全体の熱収支が2〜3%ずれるというトラブルの原因がこれだったりする。
適合メッシュと非適合メッシュ
適合メッシュと非適合メッシュの違いをもう少し詳しく教えてください。
まとめるとこんな感じだ。
| 項目 | 適合（Matching）メッシュ | 非適合（Non-matching）メッシュ |
|---|---|---|
| 界面の節点 | 完全に一致 | 一致しない |
| 補間誤差 | ゼロ（直接受渡し） | 補間手法に依存 |
| メッシュ生成 | 固体と流体を同時に考慮 | 各領域で独立に生成可能 |
| メッシュ効率 | 片方に無駄な細分化が発生 | 各領域で最適化可能 |
| 実務での採用 | 単純形状・高精度要求時 | 複雑形状で主流 |
実際、自動車のエンジンブロック冷却解析のように固体形状が複雑な場合は、固体と流体のメッシュを別々に作って非適合で繋ぐのが現場の常識だね。
非適合メッシュが実務では多いんですね。でも補間の精度が心配です…。
その心配はもっともだ。だからこそ「どの補間手法を使うか」が非適合メッシュCHTの命綱になる。
数値解法と実装
補間手法の比較
非適合メッシュで温度や熱流束を受け渡すとき、どんな補間手法があるんですか？
代表的な手法を4つ紹介しよう。
1. 最近接点補間（Nearest Neighbor）
最も近い節点の値をそのまま使う。実装は簡単だが、メッシュサイズ差が大きいと階段状のアーティファクトが出る。精度は $O(h)$ 程度。
2. 線形補間（Linear Interpolation）
周囲の節点から線形に内挿する。精度は $O(h^2)$ に向上するが、エネルギー保存が厳密には保証されない。
3. RBF補間（Radial Basis Function）
放射基底関数で滑らかな補間面を構築する。高精度だが計算コストが $O(N^2)$ で、大規模モデルではスケールしない場合がある。
4. 保存型補間（Conservative Interpolation）
GGI（General Grid Interface）やAMI（Arbitrary Mesh Interface）とも呼ばれる。面積重み付けで熱流束を分配し、エネルギー保存を厳密に満たす。CHT解析ではこれが第一選択だ。
保存型補間（GGI/AMI）の原理
保存型補間ってどういう仕組みですか？ 直感的に理解したいです。
考え方はシンプルだよ。界面の固体側メッシュ面 $A_s^i$ と流体側メッシュ面 $A_f^j$ の重なり面積 $A_{ij}$ を幾何学的に計算する。そして重み係数を作る：
この重みを使って流体側セル $j$ の界面温度を次のように算出する：
面積ベースの重みだから、全ての熱量が漏れなく保存される。「ピザを切り分けて、何枚受け取ったかで味の配分を決める」ようなものだと思えばいい。
面積で案分するから漏れがないんですね！ 各ソルバーでの名前が違うのがちょっと紛らわしいですけど…。
そうだね。OpenFOAMではAMI、STAR-CCM+ではIn-place Interface、Ansys FluentではCoupled Wall、COMSOLではIdentity Pair。名前は違うけど根底の考え方は同じ「面積ベースの保存型マッピング」だ。
実践ガイド
界面メッシュ設計の勘所
実務で界面メッシュを作るとき、教科書に載っていない「現場の知恵」を教えてください！
現場で蓄積された知見を5つ挙げよう。
- メッシュサイズ比は1:5以内に収める — 固体側と流体側の界面要素サイズ比が大きすぎると、1つの大きな固体セルに多数の流体セルが対応し、補間精度が低下する。経験的に1:5がボーダーライン。
- 流体側は必ず境界層メッシュを入れる — 界面近傍の温度勾配は流体側が支配する。y+管理と合わせて最低10層以上のインフレーション層を設ける。
- 界面形状はできるだけ単純化する — 複雑な曲面は重なり面積の計算精度が下がる。CADの段階でフィレットを除去する等の簡略化が有効。
- 界面の法線方向を統一する — 固体側と流体側で法線が逆向きになっていないか必ず確認。符号ミスで熱流束が反転する事故がかなり多い。
- 接触熱抵抗は明示的にモデル化する — 実際の界面にはグリースや空気層がある。無視すると温度差を過小評価し、実測と合わなくなる。
主要ソルバーでの設定
各ソルバーでの具体的な設定ポイントを教えてください。
| ソルバー | 界面設定 | 注意点 |
|---|---|---|
| Ansys Fluent | Coupled Wall境界条件 | Wall thickness=0を確認。Shadow faceの有無を必ずチェック |
| STAR-CCM+ | In-place Interface → CHT | Contact Resistance設定可。Mapped接続が推奨 |
| COMSOL | Heat Transferモジュール → Pair条件 | Non-conforming meshで自動補間。Identity Pair推奨 |
| OpenFOAM | chtMultiRegionFoam + AMI | createBaffles前にメッシュ品質確認。regionProperties設定が必須 |
OpenFOAMのchtMultiRegionFoamは設定ファイルが多くて大変そうですね…。
確かにハードルは高いけど、tutorials内のheatTransfer/chtMultiRegionFoamがよくまとまっているから、まずはそこから始めるのがおすすめだ。region間の境界条件定義（compressible::turbulentTemperatureCoupledBaffleMixed）が最大の肝だから、そこを理解すれば後はスムーズだよ。
ソフトウェア比較
結局、どのソルバーが界面処理に一番強いですか？
一長一短だから用途で選ぶのがベストだ。
| 観点 | Ansys Fluent | STAR-CCM+ | COMSOL | OpenFOAM |
|---|---|---|---|---|
| 非適合メッシュ対応 | 良好 | 非常に強い | 良好 | AMIで対応 |
| GUI操作の容易さ | 直感的 | 直感的 | 非常に簡単 | GUIなし |
| 接触熱抵抗 | 手動設定 | GUI内で設定可 | Pair条件で設定 | 辞書ファイルで設定 |
| 大規模並列 | 得意 | 非常に得意 | やや苦手 | 得意 |
どのソルバーでも「界面でのエネルギー収支を必ず確認する」という原則は共通だ。Post処理でFlux Reportを出して、固体入熱量と流体受熱量の差が0.1%以内に収まっているかチェックしよう。
先端技術
界面処理の研究で注目されている最近のトレンドってありますか？
3つの注目トレンドを紹介しよう。
1. オーバーセットメッシュ（Chimera）のCHT応用
回転機械で固体が動く場合、オーバーセットメッシュで界面処理する手法が実用化されつつある。従来のスライディングメッシュより大変位・大変形に強い。
2. 機械学習による補間高速化
ニューラルネットワークで界面温度場を予測する「データ駆動型補間」の研究が進んでいる。形状最適化ループの中で、メッシュ変更のたびにGGI計算をやり直す手間を省ける可能性がある。
3. Immersed Boundary法（IBM）
界面を陽にメッシュ化せず、背景メッシュ上でペナルティ項として処理する手法だ。メッシュ生成が激減するが、界面近傍の高精度化がまだ研究段階で、実務展開はこれからだね。
トラブルシューティング
CHT解析で界面まわりのトラブルって、どんなものが多いですか？
症状1: 界面で温度が不連続にジャンプする
原因：界面法線方向の反転、またはShadow faceの未定義。
対策：メッシュビューアで法線ベクトルを可視化し、両面の対応を確認する。
症状2: 全体の熱収支が数%以上ずれる
原因：非保存型補間の使用、またはメッシュサイズ比が1:10を超えている。
対策：保存型補間（GGI/AMI）に切り替え、サイズ比を1:5以内に調整。
症状3: 界面付近で振動的な温度分布が出る
原因：分割型カップリングの陽解法で界面データの更新が1ステップ遅れている。
対策：界面の緩和係数を0.5〜0.7に設定するか、陰的カップリングに変更する。
症状4: 計算が発散する
原因：固体と流体の熱伝導率差が極端（例：銅400 W/(m·K) vs 空気0.026 W/(m·K)で15,000倍）。
対策：初期値を解析解から設定し、カップリングの緩和係数を段階的に上げる。
界面の問題って見た目ではわかりにくいから怖いですね。まずは熱収支の確認を習慣にします！
それが正解だ。温度のコンター図がきれいに見えても、熱収支エラーが0.5%以上あったら必ず原因を追究すること。特に界面のフラックスレポートは毎回必ず確認する。これだけでトラブルの8割は未然に防げるよ。
関連トピック
なった
詳しく
報告