熱伝導方程式の導出は？考え方や手順も！（フーリエの法則：微小要素：エネルギー保存則：境界条件など）

熱伝導方程式の導出は？考え方や手順も！（フーリエの法則：微小要素：エネルギー保存則：境界条件など）

熱が固体の中をどのように移動するのかを数式で表したものが、熱伝導方程式です。

一見すると偏微分や二階微分が登場して難しく見えますが、出発点は身近な熱の流れとエネルギーの出入りにあります。

本記事では、フーリエの法則、微小要素、エネルギー保存則を順番に結び付けながら、熱伝導方程式が導かれる考え方をわかりやすく解説します。

一次元の場合だけでなく、境界条件や初期条件の役割、実務で式を扱う際の注意点にも触れていきます。

熱伝導方程式の全体像

それではまず熱伝導方程式の全体像について解説していきます。

温度変化を表す基本式

熱伝導方程式は、物体の各場所における温度が、時間の経過とともにどう変わるかを示す方程式です。

熱い部分から冷たい部分へ熱が流れ、その結果として温度差が小さくなっていく現象を扱います。

一次元の棒を考える場合、温度を位置xと時間tの関数として表し、T(x,t)と置きます。

ρは密度、cは比熱、kは熱伝導率、qは単位体積当たりの発熱量です。

発熱がない条件ではqがゼロとなり、式はよりシンプルになります。

左辺は物体内部に熱エネルギーが蓄えられる速さを表し、右辺は熱伝導と内部発熱の影響を示しています。

導出で押さえる三つの要素

熱伝導方程式を導くために必要な中心要素は、フーリエの法則、微小要素に対する熱収支、エネルギー保存則の三つです。

フーリエの法則によって、温度の傾きがどれだけ熱流束を生むかを表します。

次に、非常に小さな直方体や薄片を取り出し、その中へ流れ込む熱と外へ流れ出る熱を比べます。

最後に、差し引きされた熱量が微小要素の内部エネルギーの変化に等しいと考えれば、偏微分方程式へ到達します。

この順番を意識すると、記号だけを暗記するよりも式の意味を理解しやすくなるでしょう。

結論としての物理的意味

熱伝導方程式の要点は、温度が周囲との関係によって変化するという点です。

ある位置の温度が高くても、近くの温度分布が一様なら、そこへ向かう熱の出入りは釣り合う場合があります。

重要なのは温度そのものだけではなく、温度勾配と、その勾配が場所ごとにどのように変化しているかです。

この見方を持つと、二階微分がなぜ式に含まれるのかを自然に理解できます。

フーリエの法則と熱流束

続いてはフーリエの法則と熱流束を確認していきます。

温度勾配と熱の移動方向

フーリエの法則は、熱伝導を考える際の出発点です。

一次元では、単位面積当たりに流れる熱量である熱流束をQxとすると、温度勾配に比例する関係で表せます。

マイナス記号は、熱が温度の高い方向から低い方向へ流れることを示します。

右に行くほど温度が高くなるなら、熱は左向きに流れるため、符号が反転する仕組みです。

熱流束は温度差ではなく温度勾配に対応する量である点を押さえてください。

熱伝導率の役割

熱伝導率kは、材料が熱を伝えやすいかどうかを示す物性値です。

銅やアルミニウムのような金属は熱伝導率が大きく、熱が速やかに広がります。

一方で、木材、空気、発泡材などは熱伝導率が小さく、断熱材として利用されます。

同じ温度勾配であっても、kが大きければ熱流束は大きくなります。

項目	意味	温度変化への影響
熱伝導率 k	熱の伝わりやすさ	大きいほど熱が移動しやすい
密度 ρ	単位体積当たりの質量	大きいほど蓄熱量に関係する
比熱 c	温度を上げるために必要な熱	大きいほど温度が変わりにくい
熱拡散率 α	温度変化の広がりやすさ	大きいほど温度分布が均されやすい

材料選定では熱伝導率だけでなく、密度と比熱も含めて確認する必要があります。

熱拡散率による見通し

熱伝導率を密度と比熱で割った値は、熱拡散率と呼ばれます。

熱拡散率が大きい材料では、局所的な温度差が比較的短い時間で周囲へ広がります。

この値は、熱がどれだけ流れるかではなく、温度分布が変化する速さを考えるために便利です。

発熱がない場合の熱伝導方程式は、∂T／∂t ＝ α∂²T／∂x² と表せます。

式を簡潔にしたい計算問題では、熱拡散率を使った形がよく採用されます。

微小要素による熱収支

続いては微小要素による熱収支を確認していきます。

一次元の微小区間

導出では、断面積A、長さΔxの小さな棒状要素を考えます。

左側の面から入る熱量と、右側の面から出る熱量には、温度分布の違いによってわずかな差が生じます。

左面における熱流束をQx、右面における熱流束をQx＋Δxと表すと、微小要素に正味で入る熱は両者の差から求められます。

単位時間当たりに流入する熱量は、A Qx － A Qx＋Δx と書けます。

ここでは、熱流束の向きと符号を一貫して扱うことが大切です。

流入熱量と流出熱量

右側の熱流束は、左側の値をテイラー展開によって近似できます。

Δxが十分に小さいとき、Qx＋ΔxはQxに∂Qx／∂xとΔxを掛けた項を加えた形になります。

このため、正味の流入熱量は、熱流束の空間微分に比例する量へ整理されます。

熱流束が右へ行くほど小さくなる場所では、入ってくる熱のほうが多くなり、要素の温度は上がりやすくなります。

反対に、流出熱量が大きい場所では、内部の熱エネルギーが減少します。

内部発熱の追加

ヒーター、電気抵抗、化学反応、摩擦などがある場合、物体内部で熱が生まれます。

単位体積、単位時間当たりの発熱量をqとすれば、微小要素の発熱量はqAΔxです。

この項を熱収支に加えることで、温度変化をより現実的に扱えます。

発熱項は、熱源がどこにどれほど分布しているかによって位置や時間の関数になる場合もあります。

内部発熱は境界から入る熱とは別に加えるという整理をしておくと、式の組み立てで混乱しにくくなります。

エネルギー保存則による導出

続いてはエネルギー保存則による導出を確認していきます。

内部エネルギーの増加量

微小要素の体積はAΔxであり、その質量はρAΔxです。

温度が微小時間Δtの間にΔTだけ変化するとき、内部エネルギーの増加量はρcAΔxΔTで表せます。

単位時間当たりに直すと、ρcAΔx∂T／∂tとなります。

比熱cは、物質の温度を変えるために必要な熱量に関係するため、この項は蓄熱の大きさを示しています。

相変化や大きな温度範囲で物性値が変わる場合には、ρやcを一定としてよいか別途検討が必要です。

収支式から偏微分方程式まで

流入熱量、流出熱量、内部発熱量の合計は、微小要素の内部エネルギー増加量に等しくなります。

フーリエの法則を熱流束の式へ代入し、面積Aで割った後、Δxを限りなく小さくします。

すると、熱流束の変化が温度の二階微分として表現されます。

熱伝導率が位置や温度によって変化する場合は、kを二階微分の外へ出せません。

この違いは、複合材料や高温環境を扱う際に特に重要です。

三次元への拡張

実際の部品はx方向だけで熱が移動するとは限りません。

幅、奥行き、高さの各方向に熱が伝わる場合は、x、y、z方向の熱収支をすべて加えます。

その結果、温度の二階微分はラプラシアンと呼ばれる形へ拡張されます。

等方性材料でkが一定なら、ρc∂T／∂tはk∇²T＋qとなります。

ここで∇²Tは、各方向の二階微分を足し合わせた量です。

三次元化しても基本原理は微小要素のエネルギー保存であり、考え方そのものは変わりません。

初期条件と境界条件

続いては初期条件と境界条件を確認していきます。

初期条件の設定

熱伝導方程式は時間を含むため、計算を始める時点の温度分布が必要です。

これが初期条件です。

たとえば、長さLの棒全体が20度で始まるなら、tがゼロのときにT(x,0)＝20と設定します。

加熱直後の部品を冷却する問題では、場所ごとに異なる初期温度を与えることもあります。

初期条件が曖昧だと、同じ方程式でも温度履歴を一意に決められません。

温度指定境界条件

端部の温度を直接指定する条件は、第一種境界条件またはディリクレ境界条件と呼ばれます。

棒の左端を常に100度、右端を常に20度に保つ場合、T(0,t)とT(L,t)をそれぞれ指定します。

恒温槽に接している面や、十分に大きな熱容量を持つ熱源に密着する面では、この条件が近似として使われます。

境界条件の種類	代表的な表現	想定される状況
第一種境界条件	温度を指定	恒温面、温度制御面
第二種境界条件	熱流束を指定	断熱面、一定加熱面
第三種境界条件	対流熱伝達を指定	空気や冷却水に接する面

境界の物理状況を正しく読み替えることが、解析精度を左右します。

断熱条件と対流条件

熱が境界を通過しない面では、熱流束がゼロとなります。

これは第二種境界条件の一例であり、温度勾配がゼロという形で表せます。

一方、空気や水に接する表面では、物体内部の伝導だけでなく、周囲流体との対流も考慮します。

境界条件は数式を閉じるための条件ではなく、現実の熱の逃げ道や供給源を表す設定です。

導出と計算における注意点

続いては導出と計算における注意点を確認していきます。

符号と単位の確認

熱流束の符号は、どの方向を正と定めたかによって変わります。

途中で基準方向を変えると、流入と流出の関係が逆転し、式全体の符号を誤る原因になります。

熱伝導率の単位はワット毎メートル毎ケルビン、密度はキログラム毎立方メートル、比熱はジュール毎キログラム毎ケルビンです。

これらを式に入れると、左辺と右辺はいずれも単位体積当たり、単位時間当たりのエネルギーにそろいます。

次元が一致するかを確認する作業は、導出ミスを早期に発見する有効な方法です。

定常問題と非定常問題

温度が時間によって変わらない状態を定常状態と呼びます。

定常問題では∂T／∂tがゼロとなるため、熱伝導方程式は空間方向だけの方程式になります。

壁を通る熱損失の概算や、十分な時間が経過した装置内部の温度分布では、定常解析が役立ちます。

加熱開始直後、冷却途中、周期的に温度が変わる設備では、時間変化を含む非定常解析が必要です。

問題文や実測状況から、どちらを扱うべきかを最初に見極めましょう。

数値計算での扱い

複雑な形状や材料ごとに異なる物性値を持つ場合、解析解を求めることは簡単ではありません。

そのようなときは、有限差分法や有限要素法を用いて温度場を離散化します。

時間刻みや要素サイズを粗くしすぎると、温度変化を正しく追えなかったり、数値計算が不安定になったりします。

境界条件、接触熱抵抗、発熱量、物性値の温度依存性を適切に設定することが重要です。

熱伝導方程式そのものが正しくても、入力条件が現実と違えば解析結果の信頼性は下がります。

熱伝導方程式の導出に関するまとめ

熱伝導方程式は、フーリエの法則と微小要素に対するエネルギー保存則を組み合わせて導きます。

まず温度勾配から熱流束を表し、微小要素に流入する熱、流出する熱、内部発熱を整理します。

その差が内部エネルギーの増加量に等しいと置くことで、温度の時間微分と空間二階微分を含む式が得られます。

特に重要なのは、温度勾配が熱の流れを決め、熱流束の変化が温度変化を生むという関係です。

さらに、初期条件と境界条件を与えることで、実際の加熱、冷却、断熱、対流といった現象に対応できます。

式を暗記するだけでなく、微小要素へ入る熱と出る熱を図でイメージしながら追うと、導出の理解は大きく深まるでしょう。