材料・パラメータ
初期条件
境界条件
格子数: 80
Δt: — s
長さ L: 1.0 m
中央点の温度応答 T_mid(t)
理論・主要公式
$$q = -k \frac{dT}{dx}$$
フーリエの熱伝導法則:熱流束 $q$(W/m²)は温度勾配に熱伝導率 $k$(W/mK)を掛けた値。
$$\frac{\partial T}{\partial t} = \alpha \frac{\partial^2 T}{\partial x^2}$$
1次元熱拡散方程式:$\alpha = k/(\rho c_p)$ は熱拡散率(m²/s)。
$$T(x,\infty) = T_L + \frac{T_R - T_L}{L}x$$
定常解析解:定常状態では温度は線形分布(境界条件が一定の場合)。
1次元非定常熱伝導とは
🙋
「非定常」熱伝導って何ですか? 普通の熱伝導とどう違うんですか?
🎓
大まかに言うと、温度が時間とともに変化する熱伝導のことだよ。例えば、冷蔵庫から出した鉄の棒を室温に置くと、最初は棒の真ん中が冷たいけど、時間が経つと全体が室温に近づくよね。この「時間とともに温度が変わっていく様子」を計算するのが非定常解析なんだ。このシミュレーターでは、上の「材料プリセット」を「銅」から「コンクリート」に変えると、温まるスピードが大きく異なるのが体感できるよ。
🙋
え、そうなんですか!「境界条件」って何を選べばいいんですか?「断熱」ってどういう状態?
🎓
端っこの状態を決める設定だね。例えば、左端を「固定温度」で高温に、右端を「断熱」にすると、魔法瓶のように右端から熱が逃げない状態を再現できる。逆に右端を「対流」にすると、空気や水で冷やされる様子が見れるよ。実際の製品設計では、エンジン部品(高温固定)と外気(対流)みたいな組み合わせが多いね。シミュレーターで「左端BC」と「右端BC」をいろいろ組み合わせて、温度分布がどう変わるか確認してみて。
🙋
「初期温度分布」で「中央高温」を選ぶと、途中で温度が下がっていくのはなぜですか?熱は端に逃げてるんですか?
🎓
その通り!真ん中だけ熱い状態は自然には安定しないんだ。熱は温度が高いところから低いところへ流れるから、真ん中の熱が両端に向かって「拡散」していく。これが「熱拡散」の現象だよ。材料を「グラスウール」にするとこの拡散がすごく遅くて、断熱材の効果が一目瞭然だ。シミュレーターのアニメーションを見ると、熱が波のように広がっていく様子がよくわかるよね。
よくある質問
時間刻み幅Δtが大きすぎると、有限差分法の数値的安定条件(CFL条件)を満たさず発散します。Δtを小さくするか、空間刻み幅Δxを大きくして、αΔt/(Δx)^2 ≤ 0.5程度になるよう調整してください。
初期条件は時刻0秒における棒全体の温度分布(例:一様20℃)です。境界条件は棒の両端での熱の出入りを指定します。固定温度(例:左端100℃)や断熱、対流などから選べます。まずは単純な固定温度から試すと挙動が理解しやすいです。
αが大きいほど熱が速く伝わります。例えば銅(α≈1.1×10⁻⁴ m²/s)は温度変化が急速に広がり、木材(α≈1×10⁻⁷ m²/s)は非常にゆっくりです。αを10倍にすると、同じ時間で温度が届く距離が約3倍になります。
熱伝達係数h [W/m²K]は、物体表面と周囲流体間の熱の伝わりやすさを表します。静止空気中では5〜25、強制対流(扇風機など)では10〜100、水中では100〜1000程度が目安です。値を大きくすると、表面温度が周囲温度に素早く近づきます。
実世界での応用
電子機器の熱設計:スマートフォンや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次元モデルはあくまで「最も熱が流れる経路」や「断面で代表的な挙動」を理解するための「第一近似」だと心得ておきましょう。