有限差分法で1次元非定常熱伝導をリアルタイムアニメーション。初期条件・境界条件を変えて温度プロファイルの変化を体験。
このシミュレーターの根幹をなすのは、1次元非定常熱伝導方程式（拡散方程式）です。温度$T$が時間$t$と位置$x$によってどう変化するかを記述します。
$$\frac{\partial T}{\partial t}= \alpha \frac{\partial^2 T}{\partial x^2}$$ここで、$T(x, t)$は温度[°C]、$t$は時間[s]、$x$は位置[m]です。$\alpha$は熱拡散率 [m²/s]で、材料によって決まる値です。この値が大きい（銅など）ほど温度変化が速く伝わります。右辺の「温度の2階微分」が、温度勾配の勾配、つまり熱の“広がりやすさ”を表しています。
シミュレーターで選択できる「対流」境界条件は、ニュートンの冷却則でモデル化されます。例えば右端($x=L$)では次の式が適用されます。
$$-\lambda \frac{\partial T}{\partial x}\bigg|_{x=L}= h \, (T(L, t) - T_{\infty})$$ここで、$\lambda$は熱伝導率[W/(m·K)]、$h$は熱伝達係数[W/(m²·K)]、$T_{\infty}$は周囲流体の温度[°C]です。左辺は固体内部から端へ流れ込む熱流束、右辺は端から流体へ逃げる熱流束です。$h$が大きい（強制対流など）ほど、端の温度は速く$T_{\infty}$に近づきます。
電子機器の熱設計：スマートフォンやPCのCPUは発熱体です。基板上での熱の広がり（1次元的な近似も可能）と、筐体表面からの放熱（対流境界条件）をシミュレーションし、オーバーヒートを防ぐ設計に役立ちます。
建築断熱の評価：壁の内部の温度分布の時間変化を解析できます。外気温の変動（非定常）に対して、室内側の温度がどれだけ安定しているか、断熱材の効果を「断熱」境界条件などを用いて評価します。
鋳造・熱処理プロセス：金属を鋳型に流し込む際や、焼入れなどの熱処理時に、製品内部の温度が時間とともにどう変化するかを予測します。材料の熱拡散率$\alpha$が冷却速度を決定する鍵となります。
地中熱の利用：地中はある深度でほぼ一定温度を保ちます。地中に埋設したパイプの周りの温度分布の時間変化を解析し、効率的な地中熱ヒートポンプの設計に応用されます。
まず、「熱拡散率」と「熱伝導率」を混同しないでください。シミュレーターで直接調整する「熱拡散率α」は、温度変化の「速さ」を決めます。一方、対流条件で裏側に登場する「熱伝導率λ」は、熱の「流れやすさ」そのものです。例えば、断熱材のグラスウールは熱伝導率が非常に小さい（熱を通しにくい）ので、結果として熱拡散率も小さくなり、温度変化が遅くなります。この違いを理解しておかないと、材料データシートを見てパラメータを設定する際に混乱します。
次に、「対流」条件の熱伝達係数hの値はシビアです。デフォルト値で遊ぶのはいいですが、実務で使うには適切な値の調査が必須です。例えば、自然対流（空気）なら5〜25 W/(m²·K)、強制対流（ファンによる空気）なら25〜250 W/(m²·K)、水冷なら500〜10000 W/(m²·K)と、桁が全然違います。ここを「とりあえず100」で統一すると、現実とかけ離れた結果になります。
最後に、1次元モデルの限界を常に意識しましょう。このツールは棒や厚板の「厚さ方向」の解析には最適ですが、現実の熱は3次元的に広がります。例えば、スマホの基板でチップから発生した熱は、基板面内（2次元）にも広がります。1次元モデルはあくまで「最も熱が流れる経路」や「断面で代表的な挙動」を理解するための「第一近似」だと心得ておきましょう。
この1次元非定常熱伝導の計算は、「拡散」という普遍的な物理過程の基礎です。そのため、熱伝導以外の多くの工学分野と数式の形がそっくりで、考え方を応用できます。
まず、材料工学や化学工学での「物質拡散」です。熱伝導方程式の温度Tを濃度Cに、熱拡散率αを拡散係数Dに置き換えると、それが「フィックの第二法則」になります。例えば、金属表面への炭素の浸透（浸炭処理）や、半導体製造での不純物の拡散プロセスをシミュレートする際の基礎となります。
次に、土木・地盤工学における「地中熱伝達」や「コンクリートの水和熱解析」です。地中は巨大な熱容量を持つため、地中熱ヒートポンプの設計では、季節変動に伴う非定常な地中温度分布の予測が重要です。また、ダムなどの大規模コンクリート構造物では、固化時の発熱（水和熱）とその後の冷却によるひび割れを防ぐため、内部温度分布の時間変化を詳細に追う必要があります。
さらに、電気工学における「伝送線路の信号劣化」にも関連します。ケーブルや基板配線を長くすると信号が遅れてぼやける現象がありますが、これを記述する方程式も拡散型の式になります。熱の代わりに電圧や電流が「拡散」していくと考えれば、直感的に理解の助けになるでしょう。
このシミュレーターに慣れたら、次のステップとして「離散化」と「数値解法」の概念に触れてみましょう。コンピュータは連続的な微分方程式をそのまま解けません。どうしているかというと、棒をたくさんの小さなセル（メッシュ）に分割し、時間も細かいステップに区切って、近似的に計算しています。この方法を「有限差分法」と呼びます。例えば、熱伝導方程式の時間微分を $$\frac{\partial T}{\partial t} \approx \frac{T_{\text{new}} - T_{\text{old}}}{\Delta t}$$ と近似するのが第一歩です。この考え方がわかると、シミュレーション結果の精度や計算時間について、より深く議論できるようになります。
数学的には、偏微分方程式の分類を学ぶと視野が広がります。今回扱ったのは「放物型」偏微分方程式です。これに対して、振動を表す「双曲型」や、静的な場を表す「楕円型」があります。それぞれ解法や性質が異なるので、違いを知っておくことは、より複雑なCAEソフトウェアを使いこなす上で強力な基礎体力になります。
実践的な次のトピックとしては、「2次元版の熱伝導シミュレーション」に進むことをお勧めします。1次元から2次元になると、熱が平面状に広がる様子や、異なる材料が接する「界面条件」など、現実に即したより豊富な現象を扱えるようになります。それが理解できれば、市販の本格的なCAEソフトの出力結果を読み解く力も、ぐっと身につくはずです。