NAFEMS T4: 対流境界条件付き熱伝導

ベンチマークの仕様と理論的背景

NAFEMSベンチマークとは

NAFEMS（英国発の解析技術者協会）の標準ベンチマーク集は、参照解付きの小さな検証問題を提供する、ソルバー検証の世界共通言語です。T4はその熱伝導編の代表問題で、「対流境界条件（第3種・Robin境界）を正しく実装・設定できているか」を1問で確認できるように設計されています。新しいソルバーの導入時・バージョンアップ時・熱解析の教育時に、最初に解くべき定番です。

T4の問題設定

定常2次元熱伝導で、形状は幅0.6 m×高さ1.0 mの長方形板です。

支配方程式はラプラス方程式 \( \nabla^2 T = 0 \)、対流辺の境界条件は \( -k\,\partial T/\partial n = h\,(T - T_\infty) \) です。

なぜ「良い検証問題」なのか

実施手順と収束確認

標準的な実施手順

1. 要素選択 — 2次の四辺形要素（8節点）を推奨。1次要素でも解けるが収束が遅い
2. メッシュ系列 — 例えば6×10、12×20、24×40の3水準（細分比2）。評価点Eが節点上に乗る分割を選ぶと補間誤差の議論が消える
3. 境界設定 — 下辺に100°C、右辺・上辺に対流（h=750、T∞=0°C）、左辺は何も設定しない（自然に断熱）
4. 評価 — 各水準の \( T_E \) を表にし、参照値18.3°Cへの収束を確認

合否判定の目安

2次要素の細かいメッシュで \( T_E \) は18.2〜18.3°Cに収束します。実務の合格基準は参照値±0.5%（±0.1°C）以内への単調収束です。1水準だけ解いて「18.6°Cだったので概ね一致」と締めるのは検証としては不十分で、メッシュ収束の傾向まで見せて初めてT4を解いたことになります。収束次数の評価（GCI）まで行えば、そのまま社内の検証テンプレートになります。

対流境界の離散化を意識する

対流境界条件は、有限要素法では境界要素上の積分 \( \int_\Gamma h\,N_i N_j\,dS \) が剛性行列に、\( \int_\Gamma h\,T_\infty N_i\,dS \) が荷重ベクトルに入る形で実装されます。つまりhは「面に付く」量であり、誤って節点や体積に与える・面の選択を間違える・要素の表裏を取り違える、といったミスがそのまま温度場に現れます。T4で参照値が再現できれば、この境界実装チェーンが一通り正しいことの証明になります。

実務での活用法

受け入れ試験・回帰テストとしての使い方

T4の実務価値は「一度解いて終わり」ではなく繰り返し使う基準器にあります。①新ソルバー・新バージョンの受け入れ試験（アップデートで境界処理が変わっていないか）、②解析環境の移行時の同等性確認、③新人教育の第1課題（境界3種の設定を体で覚える）、④入力テンプレートの妥当性確認。入力ファイルと期待値（18.3°C±0.1）をスクリプト化しておけば、環境更新のたびに数分で回帰確認できます。

単位系・温度基準の落とし穴

T4は摂氏で定義されていますが、放射を含まない線形問題なので絶対温度でも同じ答えになります（T∞=0°C=273.15 Kと100°C=373.15 Kを使えば \( T_E = 291.45 \) K = 18.3°C）。ただしソルバーの既定温度単位と混ぜると即死します——「T∞に0を入れたが単位はKだった（=−273°C）」は定番ミスで、結果が大きく負側にずれたらまずこれを疑います。SI以外の単位系（mm系）ではhの単位換算（W/m²K → mW/mm²K）も要注意です。

拡張のさせ方——T4を起点にした段階検証

T4が通ったら、同じモデルを使って検証範囲を段階的に広げられます。①温度依存物性を入れて非線形熱伝導の反復収束を確認、②非定常化して時間積分を検証（初期温度0°Cからの過渡）、③対流hを温度依存にして境界の非線形性を検証、④固体を2材料に分けて界面連続性を検証。「参照値のある問題からの1ステップ拡張」は、複雑モデルでバグと物理を切り分ける最短経路です。

ツール別の設定方法

対流境界条件の設定名対比

モデル化の細部で揃えるべきこと

ツール間で結果を突き合わせる際は、①評価点の値の取り方（節点値か要素内補間か）、②2次要素の中間節点の扱い、③2D定式化の種類（純2D・単位厚さの平面・軸対称の誤選択がないか）を揃えます。T4程度の問題でツール間に0.1°C以上の差が残る場合、物理ではなくこれらの後処理・定式化の差がほぼ全てです。

ベンチマーク文化の現在

標準ベンチマークの体系とT4の位置

NAFEMSベンチマーク集は、熱（Tシリーズ）・線形弾性（LEシリーズ）・自由振動（FVシリーズ）などの体系で、いずれも「参照解があり、特定の機能を狙い撃ちで検証する」思想で設計されています。T4は対流境界、LE5はそり拘束ねじり、LE11は熱応力というように、1問1機能の割り切りが特徴です。ASME V&V規格の言葉では、これらはコード検証（code verification）の一部として位置づけられ、実験比較による妥当性確認（validation)の前段に置かれます。

自動化された検証スイートへ

現代のソルバー開発・運用では、NAFEMS型ベンチマークをCI（継続的インテグレーション）に組み込み、コード変更のたびに自動実行する運用が標準化しています。ユーザー側でも、解析環境（ソルバー・OS・ライブラリ）の更新のたびに社内ベンチマーク集を自動実行し、参照値との差分をレポートする仕組みを持つ組織が増えています。T4のような軽量問題はこの用途に最適で、「環境が変わっても答えが変わらないこと」を機械的に保証する安全網になります。

参照解の高精度化

古典ベンチマークの参照値は当時の計算資源で作られたもので、現代の高精度計算（hp-FEM・級数解・外挿）でより多くの桁が確定しています。T4の18.3°Cも実用上は十分ですが、より厳密な収束次数の議論をしたい場合は、自前の超細密メッシュ解をRichardson外挿した値を「作業参照値」として併用する方法が実務的です。

トラブル対応

「18.3°Cにならない」ときの診断表

ベンチマークの真の使い方

関連記事：NAFEMSベンチマーク一覧、NAFEMS LE5、GCIによるメッシュ収束性検証、平板の強制対流（hの物理的背景）。