1次元非定常熱伝導シミュレーター 戻る
熱解析

1次元 非定常 熱伝導シミュレーター — 初期境界問題と陽解法

有限差分法で1次元非定常熱伝導をリアルタイムアニメーション。初期条件・境界条件を変えて温度プロファイルの変化を体験。

材料・パラメータ

初期条件

境界条件

計算結果
中央温度 T_mid
最高温度 T_max
時刻 t
0 s
フーリエ数 Fo
熱拡散率 α
左端温度
右端温度
対流 h
格子数: 80 Δt: s 長さ L: 1.0 m
温度プロファイル T(x, t)
中央点の温度応答 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次元モデルはあくまで「最も熱が流れる経路」や「断面で代表的な挙動」を理解するための「第一近似」だと心得ておきましょう。

使い方ガイド

  1. 材料プリセットまたは熱拡散率α(m²/s)を設定:鋼0.12×10⁻⁴、アルミ0.97×10⁻⁴、コンクリート0.3×10⁻⁶などから選択
  2. 左端(x=0)と右端(x=L)の境界条件を指定:固定温度TL・TR、または対流係数h(W/m²K)による自然対流条件
  3. 初期温度分布と時間ステップを入力後、シミュレーション開始ボタンをクリック;陽解法(FTCS差分スキーム)で温度場u(x,t)の時間発展をアニメーション表示

具体的な計算例

厚さL=0.1mのステンレス鋼板(α=4.2×10⁻⁶m²/s)の初期温度100℃。左端を0℃(TL=0)、右端を100℃(TR=100)で固定し、Δt=1秒、Δx=0.01mで計算すると、中央部x=0.05mの温度は初期値50℃から約120秒後に62℃に上昇。対流条件h=25W/m²Kを加えると周辺環境(Ta=20℃)への放熱により温度上昇が25%抑制される。

実務での注意点