非定常熱伝導ヒーズラー線図計算機 戻る
熱工学

非定常熱伝導ヒーズラー線図計算機

平板・無限円筒・球の非定常熱伝導を1次近似(ヒーズラー線図)で解析。ビオ数・フーリエ数・中心温度・表面温度・総熱移動量をリアルタイム計算します。

形状・材料設定
α = 8.44×10⁻⁵ m²/s
熱的条件

1次近似解 (Fo ≥ 0.2)

平板: \(\theta_0 = C_1 e^{-\zeta_1^2 Fo}\)
円筒: \(\theta_0 = C_1 e^{-\zeta_1^2 Fo}\)
球: \(\theta_0 = C_1 e^{-\zeta_1^2 Fo}\)

ここで \(\zeta_1\) はBi数に依存する超越方程式の第1根
ビオ数 Bi
フーリエ数 Fo
中心温度 T₀ (°C)
表面温度 Ts (°C)
Q/Qmax (%)
第1根 ζ₁
温度分布(t = 300 s)

非定常熱伝導ヒーズラー線図とは

🧑‍🎓
ヒーズラー線図って何ですか?教科書に出てくるグラフですが、どうやって使うんですか?
🎓
ざっくり言うと、物が温まったり冷えたりするときの「内部の温度」をサッと調べるための便利なチャートだよ。例えば、焼き入れした金属部品がどのくらい時間で冷めるか、とか。このシミュレーターでは、上の「形状」と「材料プリセット」を選んで、ビオ数やフーリエ数をスライダーで動かすと、リアルタイムで中心温度が計算されるんだ。
🧑‍🎓
え、ビオ数とフーリエ数って何が違うんですか?両方ともスライダーで変えられますよね。
🎓
良い質問だね。ビオ数は「表面での熱の逃げやすさ」と「内部の熱の伝わりやすさ」の比だ。例えば、水冷(熱伝達率hが大きい)するとビオ数は大きくなる。フーリエ数は「経過時間」を表す無次元数だ。シミュレーターでビオ数を大きくしてフーリエ数を動かすと、表面と中心の温度差がどう変わるか、すぐに体感できるよ。
🧑‍🎓
なるほど!で、この「1次近似解」って何ですか?式が3つとも同じ形に見えますけど。
🎓
鋭い!厳密解は無限級数なんだけど、時間がある程度経つ(Fo ≥ 0.2)と、第1項だけでほぼ正確な値が出せるんだ。これが1次近似。式の形は同じでも、係数$C_1$や$\zeta_1$の値が、形状(平板・円筒・球)とビオ数によって変わる。このシミュレーターは、その面倒な係数の計算を全部裏でやって、パッと答えを出してくれる便利ツールなんだ。

物理モデルと主要な数式

無次元中心温度を求める1次近似式です。フーリエ数$Fo$がある値以上になると、この近似が非常に有効です。

$$\theta_0 = C_1 e^{-\zeta_1^2 Fo}$$

ここで、
$\theta_0 = (T_0 - T_\infty)/(T_i - T_\infty)$:無次元中心温度
$C_1, \zeta_1$:形状(平板/円筒/球)とビオ数$Bi$で決まる係数(超越方程式の根)
$Fo = \alpha t / L_c^2$:フーリエ数(無次元時間)
$T_0$:中心温度、$T_\infty$:周囲流体温度、$T_i$:初期温度
$\alpha$:熱拡散率、$t$:経過時間、$L_c$:特性長さ(板厚さの半分など)

ビオ数は、固体内部の熱抵抗と表面の熱抵抗の比を表す重要な無次元数です。

$$Bi = \frac{h L_c}{k}$$

ここで、
$h$:熱伝達率 [W/m²K](シミュレーターの「表面の熱伝達」に相当)
$L_c$:特性長さ [m]
$k$:熱伝導率 [W/mK](「材料プリセット」で選択)
$Bi \ll 1$:内部の温度がほぼ一様(集中定数系)
$Bi \gg 1$:表面と中心で大きな温度差が生じる

実世界での応用

金属の焼入れ・熱処理:高温から急冷する過程で、部品の中心が目標温度まで冷えるのに要する時間を推定します。材料(鋼、アルミなど)と冷却媒体(水、油、空気)の組み合わせでビオ数が大きく変わり、冷却速度が決まります。

食品の加熱殺菌・冷却:レトルトパウチ食品の中心まで加熱する時間、またはチルド食品を冷却する時間の設計に利用されます。円筒形状(缶詰)や平板形状(パウチ)に対する解析が可能です。

電子部品・バッテリーの熱設計:発熱するチップやバッテリーセルが、使用開始後どれくらいで定常状態に近づくかの初期評価に使われます。放熱条件(自然対流か強制冷却か)がビオ数に直接影響します。

建築建材の断熱性能評価:外気温の変動(日射や夜間冷却)に対して、壁体内部の温度がどのように応答するかを解析します。コンクリート壁(熱容量大)と断熱材の組み合わせ効果を評価する際の基礎計算として活用されます。

よくある誤解と注意点

このツールを使いこなす上で、特に初心者がハマりがちなポイントをいくつか挙げておくよ。まず「特性長さ Lc の設定ミス」だ。平板なら厚さの半分、無限円筒や球なら半径そのものがLcになる。ここを間違えるとビオ数もフーリエ数も全て狂う。例えば、厚さ20mmの板を冷却する場合、Lcは10mm(0.01m)だ。板厚全体を入力しないよう注意しよう。

次に「1次近似が使える範囲の見落とし」。ツールはFo≥0.2あたりから精度が高まるけど、冷却初期(Foが非常に小さい時)の結果はあくまで目安だ。例えば、金属を水にドブ漬けした直後の数秒間の温度は、この計算よりも実際は表面近くがもっと急激に冷える。初期過渡現象を詳細に知りたいなら、別の手法が必要になる。

最後に「材料プリセットと現実の乖離」。同じ「鋼」でも組成や処理で熱伝導率kは変わる。重要な設計判断に使うなら、必ず自社材料の実測値でkを確認し、ツールの「カスタム」設定で入力する癖をつけよう。プリセットは便利だが、あくまで初期検討用だ。

関連する工学分野

このヒーズラー線図の考え方は、熱伝導に留まらず、拡散現象全般に応用できるのが面白いところだ。例えば、物質拡散(質量伝達)。半導体製造でのドーパント拡散や、食品の塩分浸透も、熱伝導方程式と全く同じ形(フィックの第二法則)で記述される。この場合、熱拡散率αの代わりに拡散係数D、温度Tの代わりに濃度Cを考えれば、同じツールのロジックで「中心部の濃度がどれだけ上がったか」を評価できる。

また、地下水流動や地盤工学における間隙水圧の消散過程も、同様の非定常拡散方程式に従う。コンソリデーションテスト(圧密試験)のデータ解析では、時間係数(フーリエ数に相当)が重要なパラメータになる。さらに電気工学では、大電流が流れた後のケーブル導体の温度上昇(ジュール熱による過渡加熱)を評価する際にも、同じ「集中定数系」や「分布定数系」の考え方が生きてくるんだ。

発展的な学習のために

このツールの背後にある理論をもっと深掘りしたいなら、まずは「偏微分方程式の変数分離法」を学ぶのが第一歩だ。ヒーズラー線図の元になる厳密解が、なぜ無限級数になるのかを理解できる。教科書では「非定常熱伝導の厳密解」という章を探せばいい。

次に、ツールが内部で解いている超越方程式 $\zeta \tan \zeta = Bi$(平板の場合)に注目してみよう。この方程式の第1根が$\zeta_1$で、これとBiから係数$C_1$が決まる。この関係をグラフ化したものが、まさに教科書にあるヒーズラー線図そのものだ。数値計算ソフト(例えばExcelのソルバー機能でも可)でこの超越方程式を自分で解いて係数を求めてみると、理解が一気に深まる。

さらに実務的な次のステップとしては、「2次元・3次元の非定常熱伝導」「熱弾性応力の連成解析」を学ぶことを勧める。現実の部品は複雑な形状が多く、単純な平板や円筒だけでは評価できない。また、急激な加熱冷却では温度勾配に伴う熱応力が発生し、破損の原因になる。ヒーズラー線図で温度分布がわかれば、それを入力としておおよその熱応力も見積もることができるんだ。まずは一次元の原理をしっかり押さえた君なら、その次の複雑な世界への入り口もきっとスムーズに開けるはずだ。