球体の発熱による温度分布-02
式(2)
熱拡散方程式
昨今,細胞内温度計測が盛んに行われています.
直径,L,の球体において,発熱量,p,が一様に発生している状況を考えます.
熱拡散方程式は,
\(\Large C_v \cdot \partial_t T(r,t) = \kappa \cdot \nabla^2 T(r,t) + p(r,t)\)
となります.
ここで,
Cv:単位体積熱容量 [J/(m3・K)] = C(比熱[J/(kg・K)])×ρ(密度[kg/m3])
κ:熱伝導率 [W/(m.K)]
p:単位体積,単位時間あたりの発熱量 [W/m3]
です.
ここでは定常状態を考えるので,第一項が0となり,さらに極座標で書き換えると,
\(\Large - \kappa \frac{1}{r^2} \frac{d}{dr} (r^2 \cdot \frac{dT}{dr}) = p(r,t) \)
となります.ただし,右辺,p,は細胞内であり,細胞外には発熱源はないので,
\(\Large - \kappa \frac{1}{r^2} \frac{d}{dr} (r^2 \cdot \frac{dT}{dr}) = \begin{pmatrix} p : (0<r<R) \\ 0 : (R<r) \end{pmatrix} \)
と場合分けして考えます.ここでRは半径です.
次元
ここで,次元を考えてみましょう.
一次元
・左辺
\(\Large C_v \cdot \partial_t T(r,t) : \left[ \frac{J}{m^3 \cdot K} \right] \left[ \frac{K}{s} \right] = \left[ \frac{W}{m^3} \right] \)
・右辺
\(\Large \kappa \cdot \nabla^2 T(r,t) : \left[ \frac{W}{m \cdot K} \right] \left[ \frac{K}{m^3} \right] = \left[ \frac{W}{m^3} \right]\)
となり,次元は,単位時間,単位体積あたりのエネルギーの変化 = 単位体積あたりの仕事率,となります.
従って,発熱量pの次元も,[W/m3]となる(トータルの仕事率でないことに注意).
極座標
\(\Large r^2 \frac{d}{dr} (r^2 \cdot \frac{dT}{dr}) : [m^2 \cdot \frac{K}{m} = m \cdot K] \)
\(\Large \frac{d}{dr} (r^2 \cdot \frac{dT}{dr}) :[ K] \)
\(\Large - \kappa : [ \frac{W}{m \cdot K}]\)
\(\Large - \kappa \frac{1}{r^2} \frac{d}{dr} (r^2 \cdot \frac{dT}{dr}) : [ \frac{W}{m \cdot K} \cdot \frac{1}{m^2} K] = [ \frac{J}{m^3 \cdot s}] \)
となり,単位時間,単位体積あたりのエネルギーの変化 = 単位体積あたりの仕事率,となります.
単位時間,単位体積あたりの発熱量と総発熱量
この二つの関係は,総発熱量をWとすると,
\(\Large W = p \cdot V\)
となります.ここでVは球の体積です.
また,球表面での熱の流出量,q,は,定常状態であるので総発熱量と一致します.従って,
\(\Large q = \frac{W}{S}\)
となります.ここでSは球の表面積です.従って,
\(\Large q = \frac{W}{S} = \frac{pV}{S} = p \frac{\frac{4}{3} \pi R^3}{4 \pi R^2} = \frac{p \ R}{3}\)
となります.
球の内外の環境のすりあわせ
方程式を解く上で,重要なのは,
界面での内外の温度は一致
界面での内外の熱の流出量は一致
という二つの条件が必要となります.
必要とならないのは,
温度の傾き
は連続的にならなくても問題ありません.界面を挟んで熱伝導率が異なる場合には不連続(温度は一致,温度がなめらかに変化しない)になります.
球の外と内側では温度は連続的につながっていると考えられるので,r=R (Rは半径)で値が一致する必要があります.
\(\Large T_{in} \vert_{r=R} = T_{out} \vert_{r=R} \)
の条件を満たす必要があります.
さらに,球表面での熱の流出量,q,は,温度の時間微分と熱伝導率との積とも記すことができるので,
\(\Large - \kappa_{in} \frac{dT_{in}}{dr} \vert_{r=R} = - \kappa_{out} \frac{dT_{out}}{dr} \vert_{r=R} = q\)
\(\Large \frac{dT}{dr} \vert_{r=R} = q = - \frac{p \ R}{3 \kappa} \)
と表すことができます.
次ページでは,球の内外での熱伝導方程式を考えていきましょう.