円筒座標系の熱伝導 — トラブルシューティング
より充実した内容を cylindrical-conduction.html でご覧いただけます。
よくあるトラブルと対策
円筒の熱伝導解析で注意すべき点を教えてください。
頻出トラブルを整理しよう。
1. 2Dと3Dで結果が合わない
原因: 2Dモデルで軸対称設定を忘れている。平面応力/平面ひずみのデフォルトでは $r$ の効果が入らない。
対策: Ansys PLANE55ならKEYOPT(3)=1、Abaqusなら軸対称要素を選択。COMSOLならモデル作成時に「2D Axisymmetric」を選ぶ。
2. 多層構造で温度分布が不連続
各層の接続がうまくいかないことがありますか？
原因: 層間のノードが共有されていない、またはTied Contact/Bonded Contactの設定漏れ。
対策: 共通面でノードマージするか、接触条件を適切に設定。意図的に接触熱抵抗を入れる場合はGap Conductanceを設定する。
3. 内面対流係数の見積もりミス
配管内面の $h$ は流速、流体物性、管径で大きく変わる。
| 流体条件 | $h$ [W/(m$^2$ K)] |
|---|---|
| 空気・自然対流 | 5〜25 |
| 空気・強制対流 | 25〜250 |
| 水・強制対流 | 500〜10,000 |
| 蒸気凝縮 | 5,000〜50,000 |
| 沸騰 | 2,500〜75,000 |
水と空気で100倍以上違うんですね。
Dittus-Boelter式やGnielinski式で算出するのが基本だ。$\text{Nu} = 0.023 \text{Re}^{0.8} \text{Pr}^{0.4}$（Dittus-Boelter）。Re数を間違えると $h$ が桁で狂う。
4. 単位系の混乱
SI系（m）とmm系で $k$ の値が変わる。$k = 16.3$ W/(m K) はmm系では $k = 0.0163$ W/(mm K) だ。Ansys Mechanical のWorkbench環境はmm系がデフォルトなので要注意。
単位系は本当に要注意ですね。
理論解で検算する習慣をつければ、単位系のミスはすぐに発見できる。$q = 2\pi k L \Delta T / \ln(r_2/r_1)$ を手計算して結果と照合するだけで良い。
内径・外径の単位ミスが大事故に
配管の円筒熱伝導解析でriとroをmm入力すべきところをmで入力し、壁厚が1000倍に計算されてしまうヒューマンエラーは現場で繰り返される。1999年のNASA Mars Climate Orbiter墜落事故（ポンド力とニュートンの単位混在）を教訓に、CAEソルバーが入力単位を明示的に確認するダイアログを追加したメーカーが複数ある。
トラブル解決の考え方
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——円筒座標熱伝導の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告