基本中のキホン「熱伝導」をおさらいする：CAE解析とExcelを使いながら冷却系設計を自分でやってみる（2）（2/4 ページ）

CAE解析とExcelを使いながら冷却系の設計を“自分でやってみる／できるようになる”ことを目指す連載。連載第2回では、基本中のキホンとなる「熱伝導」についておさらいする。

LINE

Hatena

フーリエの法則

　次に、熱伝導の法則である「フーリエの法則（Fourier's law）」を説明します。図3は左右対称なので、これからは右半分だけを考えます。図4に右側の固体と固体内部の温度分布を示します。

図4　右側の固体と固体内部の温度分布［クリックで拡大］

　このとき、固体内部の小さな平面を考え、そこを短い時間Δt［s］の間に通過する熱量をΔQ［J］とします。小さな平面の面積をΔs［m²］、温度分布を座標xの関数T（x）とすると、これらには以下の関係があります。これがフーリエの法則です。

式6

　参考文献［1］では、式6の微分記号は偏微分記号ですが、1次元問題なので普通の微分記号としました。単位時間（1［s］です）に通過する熱量は、ΔQをΔtで割ったもので、また、単位面積、単位時間当たりに通過する熱量は、ΔQをΔsΔtで割ったもので、これが熱流束qです。式で表すと以下となります（式7）。

式7

　式7を式6に代入すると、式8となります。こちらの方が分かりやすいですね。

式8

　式8にマイナス（－）があるのは温度勾配dT／dxがマイナスのとき、－dT／dxはプラス、つまり熱流束はプラス値となります。温度勾配がマイナスということは温度が下がっていく方向で、その方向に熱が流れるということです。

　フーリエの法則の3次元版は熱流束qがベクトル量となり、次式のように表すことができます。ベクトル量は太字で表記します。

式9

式10

　grad Tは「スカラー場Tの勾配」を表しています。∇はハミルトン演算子でベクトル量です（参考文献［2］）。以降の連載の中で、∇²と表記した演算子が登場しますが、これは2つの∇同士の内積です。

　では、固体内部の温度分布を求めましょう。固体の断面積A_sは一定値なので、左側の流入熱流束q_sと右側の流出熱流束q_wは等しく、熱流束qをxの関数としたくても関数ではなく一定値となり、式5でした。では、式8を変形して積分しましょう。

式11

式12

　積分定数Cは、境界条件（x＝x_wでT＝T_w）から求めます。この条件を式12に代入すると、積分定数Cが求まります。以下です。

式13

式14

　式14を式12に代入します。ついでに、式5も代入します。

式15

　T_wは未知数なのですが、いずれ熱伝達を説明する際に流体温度から求まります。式15はxの1次式ですね。そういうわけで、図4の温度分布を直線状にしました。