🙋
「半無限固体の非定常熱伝導」って、名前が難しすぎます…。これって結局どういう状況の話なんですか？
🎓
名前ほど難しくないよ。たとえば真冬の朝、分厚いコンクリートの床にお湯をこぼしたとする。表面はすぐ熱くなるけど、床のずっと奥はまだ冷たいまま。この「片側の表面だけ温度が変わって、反対側はうんと遠くてまだ影響が届いていない」状況を、理想化したのが半無限固体なんだ。「非定常」は時間とともに温度が変わっていく、という意味だよ。
🙋
なるほど。でも、内部の温度ってどうやって計算するんですか？場所も時間も変数なら、すごく複雑そうです。
🎓
ここがこの問題の美しいところでね。深さ x と時間 t を別々に追わなくても、η = x/(2√(αt)) という「相似変数」一つにまとめられるんだ。すると無次元の温度は θ = erf(η) という、誤差関数（erf）一発で書ける。左のスライダーで深さや時間を変えると、ツールがこの erf を計算して温度を出している。深さ50mm・600秒のデフォルトだと η≈0.30、温度は約141°Cになるよ。
🙋
θ や η が同じなら温度の状況も同じ、ということですか？じゃあ「浸透深さ δ」っていうのは何ですか？
🎓
そう、相似というのはまさにそれ。δ ≈ 4√(αt) は「熱がここまでは届いた」という目安の深さだ。δ より深いところは、まだほぼ初期温度のまま。おもしろいのは δ が時間の平方根で進むこと。4倍待っても深さは2倍にしかならない。だから「もう少し待てば芯まで温まる」と思っても、実際はなかなか進まないんだ。アルミの上で δ をスライダーで見ると、鋼よりずっと速く深くまで届くのが分かるよ。
🙋
材料を変えると速さが全然違うんですね。表面熱流束っていうのも気になります。なんで時間が経つと減っていくんですか？
🎓
qs = k(Ts−Ti)/√(παt) で、最初は理論上ほぼ無限大、その後 1/√t でどんどん減る。理由はシンプルで、加熱したての表面のすぐ内側にはまだ冷たい材料があって、温度の落差が急なんだ。だから熱がドッと流れ込む。時間が経つと温まった層が厚くなって落差がなだらかになり、流れ込む熱が減る。焼入れの冷却速度を見積もるときなんかに、この量がとても重要になるよ。
🙋
最後に一つ。「半無限」って現実にはありえないですよね。実際の有限な物体には使えないんですか？
🎓
いい質問だ。コツは「浸透深さ δ が物体の実際の厚さより小さいうちは半無限とみなせる」ということ。たとえば厚さ300mmのコンクリート壁でも、加熱から数分なら熱は表面付近にしか届かないから、半無限の式がそのまま使える。逆に δ が厚さを超え始めたら、もう反対側の境界が効いてくるので、有限厚さの解（フーリエ級数など）に切り替える必要があるんだ。
半無限固体の非定常熱伝導とは何ですか？
片側の表面だけが境界で、反対側が無限に遠い（厚さが十分大きい）とみなせる固体の、時間とともに変化する温度伝導です。表面温度が時刻ゼロで急に Ts に変わると、その熱が内部へ少しずつ拡散していきます。温度分布は誤差関数 erf を使って T(x,t) = Ts + (Ti − Ts)·erf(x/(2√(αt))) と表せます。x は表面からの深さ、t は経過時間、α は熱拡散率です。本ツールはこの解析解で内部温度をリアルタイムに計算します。
熱浸透深さ δ とは何を表していますか？
表面温度の変化が固体内部にどこまで届いたかを示す目安の深さで、δ ≈ 4√(αt) で見積もります。深さ δ より深い領域は、温度がまだほぼ初期温度 Ti のまま（約99%が未変化）と考えてよい範囲です。時間 t の平方根に比例して進むのが特徴で、4倍の時間が経つと浸透深さは2倍になります。半無限モデルは、この δ が実際の物体の厚さより小さい間だけ有効です。
熱拡散率 α と熱伝導率 k はどう違いますか？
熱伝導率 k は、定常状態でどれだけ熱を「通す」かを表す量です。一方、熱拡散率 α = k/(ρ·cp) は、温度変化がどれだけ「速く」伝わるかを表します。非定常問題では α が支配的で、α が大きいほど短時間で内部まで温度が届きます。例えばアルミは α が鋼の約8倍あり、同じ時間でずっと深くまで温まります。表面熱流束のように「熱量そのもの」を問うときは k も必要になります。
表面熱流束はなぜ時間とともに減るのですか？
表面熱流束は qs = k(Ts − Ti)/√(παt) で与えられ、時刻ゼロでは無限大、その後 1/√t で減衰します。これは、加熱の初期には表面のすぐ内側にまだ冷たい材料があり大きな温度勾配が立つのに対し、時間が経つと温まった層が厚くなって勾配がなだらかになるためです。焼入れの冷却速度や、瞬間的な接触の伝熱を評価するときに重要な量です。
金属の熱処理・焼入れ： 鋼の表面焼入れでは、加熱した部品を急冷したとき、表面から内部へ温度が下がっていく速さが組織と硬さを決めます。半無限固体モデルは、焼入れ初期に表面近傍がどれだけ速く冷えるか、所定の深さがいつ目標温度に到達するかの当たりづけに使われます。浸透深さ δ が部品の有効硬化層厚さの目安にもなります。
建築・地中の熱： 厚いコンクリート壁や地盤は、日射や気温変化に対してまさに半無限固体としてふるまいます。日中の表面加熱が壁内部のどこまで届くか、地表の季節変化が地中何メートルまで影響するかは、この erf 解で見積もれます。地中熱利用やヒートアイランド評価の基礎計算として広く使われます。
瞬間的な接触・火傷の評価： 熱い金属や熱湯に皮膚が短時間触れたときの温度上昇も、初期の数秒間は半無限固体の問題として扱えます。同じ表面温度でも、熱拡散率の大きい金属に触れたほうがアルミ・銅は速く熱が皮膚側へ流れ込むため、より危険になります。火傷リスクや保護具の設計の基礎になります。
CAE解析の検証・初期評価： 非定常熱伝導の有限要素解析を行う前に、この解析解で「表面付近の温度がいつ何度になるか」を概算しておくと、メッシュや時間刻みの妥当性を判断できます。FEM結果が erf 解と表面付近で大きくずれていれば、時間刻みが粗すぎる、または境界条件が誤っているサニティチェックとして機能します。
まず最大の落とし穴が、「半無限モデルを有限厚さの物体にいつまでも使い続ける」 ことです。このモデルは、熱浸透深さ δ ≈ 4√(αt) が物体の実際の厚さより小さい間だけ有効です。δ が厚さに達すると反対側の境界の影響が現れ、erf 解は実際より高い温度を予測してしまいます。経過時間を長く取るときは必ず δ をチェックし、それが物体厚さを超えたら有限厚さの解（フーリエ級数や数値解）に切り替えてください。本ツールは判定文でこの注意を表示します。
次に、「境界条件を取り違える」 こと。この式は表面温度が瞬間的に一定値 Ts に固定される（第一種境界条件）場合のものです。実際の加熱では、表面が一定の熱流束で加熱される場合や、対流で周囲流体と熱交換する（第三種境界条件、ビオ数で記述）場合があり、それぞれ温度分布の式がまったく異なります。表面温度が本当に一定とみなせるのか、まず物理状況を確認してください。対流が支配的なら、誤差関数ではなく集中熱容量法やビオ数による解が必要です。
最後に、「熱拡散率 α と熱伝導率 k を混同する」 という誤解。非定常問題で温度がどこまで速く届くかを決めるのは α であって k ではありません。コンクリートは木材より k が約10倍大きいですが、α の差はそこまで大きくなく、温度が浸透する速さの違いは k の比ほどではありません。一方、表面熱流束 qs のように「単位時間あたりの熱量」を問う量には k が直接効きます。どちらの量を求めたいのかで使う物性値が変わる点に注意してください。