🙋
このシミュレーターで熱源をクリックして置いたら、色がじわーっと広がっていくのが面白いです！あの広がり方って、何で決まるんですか？
🎓
大まかに言うと、「熱拡散率 α」という数値で決まるんだ。このシミュレーターの左パネルで「材料」を切り替えてみて。銅（α≈112）を選ぶと一瞬で色が広がるけど、木材（α≈0.13）にするとなかなか広がらないだろ？これは銅が木の約860倍も熱を伝えやすいからなんだ。実際にフライパンの取っ手が木やプラスチックなのは、この熱拡散率の差を利用しているんだよ。
🙋
え、そんなに違うんですか！ところで、「境界条件」っていうボタンがあるけど、ノイマンとディリクレって何が違うんですか？
🎓
いい質問だね。ディリクレ条件は「壁の温度を固定する」ということ。例えば、壁を0°Cに保つ冷蔵庫の壁みたいなイメージだ。一方、ノイマン条件は「壁を通る熱の流れをゼロにする」断熱の条件。シミュレーターで試すと、ディリクレだと端から熱が逃げて温度が下がるけど、ノイマンだと熱が逃げないから全体がゆっくり均一化するのがわかるよ。
🙋
なるほど！熱源と冷源を両方置いて、しばらく待つと温度が直線的に変化するような模様になりました。これって何か意味があるんですか？
🎓
それは「定常状態」って呼ばれる状態で、時間が経っても温度分布が変化しなくなった状態だ。1次元だと温度は直線分布になる。例えば、左壁が100°C、右壁が0°Cなら、定常状態では温度が50°Cの線が真ん中にくる。これはフーリエの法則から導かれる基本中の基本で、建築の断熱設計や電子機器の放熱設計でもこの定常解を基に計算するんだよ。
厳密には異なります。熱伝導はフーリエの法則 $q = -k \nabla T$ で記述される熱の移動現象そのものを指し、熱拡散は温度場の時間変化 $\partial T/\partial t = \alpha \nabla^2 T$ を表します。熱拡散率 $\alpha = k/(\rho c_p)$ は熱伝導率 $k$ だけでなく、密度 $\rho$ と比熱 $c_p$ も含むため、「熱を伝えやすいか」と「温度が変わりやすいか」は別の概念です。
このツールは陽解法（FTCS法：前方差分時間・中心差分空間）を使用しています。安定性条件 $\alpha \Delta t / (\Delta x)^2 \leq 0.25$（2次元の場合）を満たすようにタイムステップが設定されています。空間分解能は100×100グリッドで、定性的な現象理解には十分ですが、工学的な精密計算にはより細かいメッシュと陰解法が必要です。
赤いブラシで描いた「熱源」は、その位置の温度を常に100°Cに固定するディリクレ境界条件として機能します。物理的には、一定温度に保たれたヒーターや恒温槽と接触している状態に相当します。青い「冷源」は0°C固定で、冷却装置や氷水との接触を模擬しています。消しゴムを使うと、その点は自由に温度変化できるようになります。
このツールは等間隔格子の有限差分法（FDM）を使っていますが、実際のCAEソフト（Ansys、COMSOL等）では有限要素法（FEM）や有限体積法（FVM）が主流です。FEM/FVMは複雑な形状や非一様メッシュに対応でき、対流や輻射も含めた連成解析が可能です。ただし、熱拡散の基本原理（フーリエの法則）は同じなので、このツールで培った直感はCAE解析の結果解釈にも役立ちます。
電子機器の熱設計：CPUやパワー半導体の発熱を基板やヒートシンクに効率よく拡散させる設計では、まさに2次元（3次元）の熱拡散解析が必須です。このシミュレーターで「熱源を1点に置いて広がりを観察する」操作は、チップからの放熱パスの設計と同じ感覚を体験できます。
建築の断熱・結露対策：壁や窓を通じた熱の流出入は、熱伝導方程式で記述されます。冬場に窓の室内側が結露するのは、ガラスを通じた熱拡散により表面温度が露点以下に下がるためです。このツールで冷源と熱源を壁の両側に置くと、その温度勾配を直感的に確認できます。
金属加工・溶接の温度管理：溶接時の熱影響部（HAZ）の広がりは、母材中の熱拡散で決まります。材料ごとの熱拡散率の違い（鋼 vs アルミ vs 銅）は、溶接条件の設定に直結します。シミュレーターで材料を切り替えて拡散速度を比較すると、なぜアルミの溶接には大電流が必要なのかが実感できます。
食品工学・殺菌プロセス：缶詰やレトルト食品の加熱殺菌では、容器の中心部まで十分な温度が到達する時間を計算する必要があります。食品の熱拡散率は水に近い値で、金属よりはるかに小さいため、厚い食品ほど加熱に時間がかかります。
まず一つ目の注意点は、「熱伝導率が高い＝温度が上がりやすい」は間違いということ。銅は熱伝導率が高いけれど、比熱容量も大きいため、「温度の変わりやすさ」を表す熱拡散率は熱伝導率ほど極端には大きくならない。シミュレーターで鋼と銅を比べると、銅の方が熱が広がるのは速いけど、「温度が上がりにくい」のは比熱の影響なんだ。
次に、「シミュレーションの速度＝実際の時間」と思い込まないで。このツールの「速度」スライダーは1フレームあたりの計算ステップ数を変えているだけで、実際の物理時間のスケールは材料のαに依存する。鋼で1秒かかる拡散が、木材では何十秒もかかる。画面上では同じ速度で動いて見えるけど、物理的な時間スケールは「シミュレーション時間 t」の表示を必ず確認しよう。
最後に、2次元シミュレーションの限界を理解しておこう。現実の熱拡散は3次元で起こるため、薄い板の面内拡散には良い近似だけど、厚みのある物体（例えば丸太の断面）の温度分布を正確に予測するには3次元解析が必要。また、このツールは対流や輻射を無視した純粋な熱伝導のみを扱っている。実際の放熱設計では、空気の対流（自然対流・強制対流）が支配的な場合も多いことを忘れずに。