複数の膜材料を同じ条件——同じ基板、同じスリットパターン——の上で走らせ、ストップウォッチ片手にどれが最も速く熱を逃がすかを見る。本ガイドではその原理、本ツールが前提としていること、そして熱工学の専門家でなくても結果を読める方法を説明します。
導電性膜の候補がいくつかある——Ginestium™、CVDグラフェン、CVDダイヤモンド——というとき、単純な問いが生まれます。この基板とこの加工パターンの上で、実際にどれが最も熱をうまく広げるのか？ 本ツールはプレートの片端を固定温度に保ち、各材料をそれぞれ単独で見たときに、熱がどれくらいの速さ・強さでもう一方の端へ進んでいくかを観察します。
材料は厳密に対等な条件で比較されます——同じ基板、同じスリットパターン、同じ熱条件——変わるのは膜だけです。
熱異方性シミュレーター（併設ツール）は「このパターンの実効伝導率はいくつか？」という定常状態の問いに答えます。本ツールはそれとは異なる、より具体的な問いに答えます——「ある時点で熱源が発生したとして、指定した距離にあるセンサーがその違いを感じ取るまでどれくらいかかるか？」。これは、現実的な伝播シナリオのもとで材料同士を比較するのに役立ちます。
本ツールは、非常に直接的な意味でデューデリジェンスのためのツールでもあります——すべての既定値には出典があり（§07参照）、記憶に頼るのではなく、技術パートナーの前で比較結果を根拠づけて説明できるようになっています。
最後に、3つの熱ツールの中で唯一時間軸を持つのが本ツールであり、したがって材料の比熱容量が関与する唯一のツールでもあります。この点はよく質問される部分であり、§10の注記で扱います。
シナリオは意図的にシンプルです——50 × 50 mmのプレート、固定温度（既定120 °C）に保たれた左端、そして熱が右へ浸透していくにつれて徐々に温まっていく残りの部分。端から選んだ距離に仮想的な「センサー」を置き、その温度が時間とともに上昇していく様子を観察します。
この進行過程は、閉形式の式によって直接計算できます（時間刻みでのシミュレーションは不要）。ただし材料ごとに1つの重要な数値——実効拡散率、αeffと表記——を知っている必要があります。これは伝導率と蓄熱能力の両方を考慮したうえで、熱がどれくらいの速さで伝わるかを示す値です。
比較対象の膜は、伝導率が驚くほど高いことがあります——しかしその厚みは数百ナノメートル程度にすぎず、下にある基板は数十から数千マイクロメートルにも及びます。システムの実効拡散率は、両者の厚み加重平均です。基板に対して膜が薄いほど、膜の性能にかかわらず、最終結果は基板自体の特性に近づいていきます。
だからこそ、本ツールは各膜のk·t（伝導率×厚み、x方向）も表示します——これは横方向に熱を広げる膜にとっての真の評価指標であり、伝導率単体よりも意味があります。伝導率が非常に高くても薄すぎる材料は、伝導率は劣っても厚みのある材料に負けることがあるのです。
スリットパターン（熱異方性シミュレーターと同じ幾何形状）は、同じ計算エンジンによって均質化されます——x方向を優先するデザインはここでも速いままですが、その効果はk_xx/k_yyの比としてではなく、伝播速度の差として現れます。本ツール固有のものは、時間経過に伴う伝播の式、それを支える比熱容量、そして最大4材料を同一の物理スケールの拡散マップ上で並べて比較できることです。
左端はその高さ全体にわたって均一な温度であると仮定されます——したがって、プレートのどの場所を見ても熱は同じように進みます（上下端の縁効果はありません）。これによって、時間刻みでのシミュレーションなしに閉形式の計算が可能になります。
プレートは右方向に無限に広がっているものとして扱われます——断熱された右端も、対流や放射による損失もありません。前線が反対側の端に到達するまでは妥当ですが、それを超えると、この式は実際の昇温を過小評価します。
膜と基板は、あらゆる位置xで同じ局所温度にあると仮定されます——本モデルは膜の昇温と基板の昇温を別々には追跡せず、一方が他方に対して遅れることも考慮しません。これは「等価平板」と呼ばれる前提であり、界面の扱い方（§10参照）と直接結びついています。
界面熱抵抗（膜／基板間）は計算され、面外方向の数値には表示されますが、横方向の伝播前線の計算には反映されません——これは意図的なモデル選択であり、本ツールが参照した元実装と同一です。詳細と正当化は§10にあります。
スリットによって除去された部分（空気、真空、エアロゲルなど）では、熱容量は無視できるものとして扱われます——実効拡散率の計算に入るのは、実際に中実な膜の割合のみです。
密度、熱容量、伝導率は室温での値として扱われ、過渡状態の間、温度に応じて変化することはありません——既定シナリオではプレートが20 °Cから120 °Cまで変化するにもかかわらずです。
カタログの材料については、厚みと面内伝導率は自由に編集できます。面外方向の伝導率、CTE、フォノン平均自由行程、密度、熱容量はカタログの値のままです。
カタログの各材料は、厚みと伝導率の入力欄の下に「文献値の範囲」という行を表示し、公開された測定値に基づく下限・上限の幅を示します。この上にカーソルを合わせると、完全な根拠（引用された測定値、前提条件、単位）がインターフェースの使用言語で表示されます。
3つの材料、3つのステータス：
ダイヤモンドとグラフェンについては、既定値はいまや実際の薄膜（ダイヤモンドはナノ結晶、グラフェンは基板上のfew-layer）のものであり、バルク結晶や実験室で測定された浮遊単層のものではありません——後者は、実際に成膜された膜の現実よりもはるかに有利なケースです。Ginestium™については、公開された測定値がないため、中心値は比較可能なグラファイト系炭素との類推によって設定されています。範囲の上限は社内で測定された最良ケースであり、既定として想定すべき平均値ではありません。
注意：「文献値」というステータスは、必ずしもその材料のすべての特性をカバーしているわけではありません。特にGinestium™については、密度と熱容量は代替値であり実測値ではありません——詳細は§10の注記にあります。
要約（CSV） — 直近の比較結果の要約表で、表計算ソフトですぐに使える形式です。
設定（JSON） — 本ツールのすべての設定（材料、パターン、基板、熱条件）をエクスポートします。結果は含まれません。そのまま再インポートできるため、基準となる構成を固定して共有するのに便利です。
このセクションは本ガイドの技術リファレンスです。サーバー側エンジンが実際に実行する計算式を示し、それぞれのモデル選択の理由を説明します。ツールの利用に必須ではありません。
時間の問題を扱う前に、スリットパターンは熱異方性シミュレーターと全く同じサーバーモジュール・同じ計算式で均質化されます——つまり、この材料でその シミュレーターを手動で操作した場合と同じ数値になります。その処理の流れは次の通りです。
これら6ステップの詳細と正当化は、熱異方性シミュレーターガイドの§09にあります。以下は本ツール固有の内容です。
体積熱容量は示量量です——積み重なった2つの材料は、それぞれが蓄えるものの和を蓄えます。したがって厚み加重平均が正しい組み合わせ方であり、これはK_xxで用いられているのと全く同じ構造です。分子と分母が同じ対象を扱っているからこそ、ステップ3での両者の商が意味を持つのです。
(1 − ψ)という係数は明示的な前提を表しています——スリットを満たす媒質（空気、真空、エアロゲル）は、この尺度では有意な蓄熱能力を持たないという前提です。残った中実な材料のみがカウントされます。
これは熱拡散率そのものの定義です——α = k/(ρ·c_p)。これは、材料が熱を運ぶ能力（分子）と、途中で熱を蓄える能力（分母）の間の競合を測る指標です。多くの熱を蓄える材料は、伝導が良くても前線の進行を遅らせます——これはまさに伝導率だけでは得られない情報であり、本ツールが存在する理由そのものです。
用いられるのは膜単体のk_xxではなくK_xxです——右方向に熱を伝えているのは膜だけでなく、基板を含めたプレート全体だからです。
熱源は左端の高さ全体にわたって均一であるため、問題は空間変数を1つしか持たなくなります。これは1次元の熱伝導方程式に帰着します。
この問題には、相似変数 η = x / (2√(α·t)) によって得られる厳密解があります。
これは過渡伝導問題の中でも、厳密な解析解を持つ数少ない問題の一つであり、均一熱源という前提こそがそれを可能にしています。その帰結は具体的です——時間刻みが不要なため、安定性条件を守る必要も、数値誤差が蓄積することもありません。3時間後の結果も3秒後の結果も、全く同じ計算コストで求められ、4材料すべてが数マイクロ秒で計算できます。erfc関数は標準ライブラリのものを使用しており、JavaScript実装であれば受け入れざるを得なかったであろう多項式近似の誤差も排除されています。
読み解く上でのポイント：この解はxとtに対して、x/√(α·t)というまとまりを通してのみ依存します。拡散率が4倍異なる2つの材料は、同じプロファイルを描きますが、一方が他方のちょうど半分の時間で到達するだけです。これがα_effを材料間で直接比較可能にしている理由です。
「10%到達遅延」と「50%到達遅延」は式(4)の逆問題です——到達すべき割合を固定し、時間を求めます。
標準ライブラリにはerfcの逆関数が存在しません。zは[0 ; 8]区間での二分法によって80回反復で求められます——erfcは狭義単調減少なので、二分法は80回目よりずっと前の段階で機械精度に達します。
遅延はx²に比例し、xには比例しません。センサーの距離を2倍にしても到達時間は2倍にはならず、4倍になります。これはあらゆる拡散現象に共通する特徴であり、優れた放熱膜であっても数ミリメートルを超えると効果が薄れてしまう理由を説明しています。
この量は過渡モデルには関与しません——指標として別途計算されるものです。これは、伝播シナリオとは無関係に、システムレベルで膜が何らかの重みを持つかどうかを決める値です——システムは膜と基板を厚みの比率で平均するため、k·tが基板よりはるかに小さい膜は、伝導率がいかに高くても見えなくなるまで希釈されてしまいます。
そして拡散マップは、式(4)をプレートの50 mmにわたって41点でサンプリングし、0から1の間に正規化した値です。マップのすべての行は同一です——これは表示上の簡略化ではなく、均一な熱源に対してこの解が実際に予測する結果そのものです。
成膜材料に関して繰り返し寄せられる2つの質問があります。比熱容量はどこに反映されているのか、そして界面熱抵抗はどう扱われているのか。3つの熱ツールで実際に採用されている前提とともに、率直にお答えします。
決め手となるのは体制（レジーム）です。定常状態では、比熱容量は一切効果を持ちません——温度場が落ち着いた後は、熱の分布を決めるのは伝導率のみです。逆に過渡状態では、伝導率とともに熱拡散率 α = k/(ρ·c_p) を決めるのは比熱容量であり、したがって温度が変化する速さを左右します。
3つのツールは次のように分かれます。
| ツール | 体制 | 膜のρ·c_p |
|---|---|---|
| 熱拡散コンパレーター | 過渡状態 | 明示的に関与する |
| 熱異方性シミュレーター | 定常状態 | 関与しない |
| パラメトリックスイープ | 定常状態 | 関与しない |
2つの定常状態ツールについては、これは見落としでも暗黙の組み込みでもありません。両者が解いている方程式——∇·(k∇T) = 0——にはρ·c_pが現れず、そもそもパラメータとしてこれらの項を持っていません。材料カタログは3つのツールで共有されており、ρとc_pをきちんと記録していますが、それを消費するのは本ツールだけです。
本ツールでは、比熱容量は実効拡散率を通じて関与します——これは§09の式(1)と式(2)であり、ここに合わせて再掲します。
成膜層についてここで明示的に置かれている前提は2つあります。まず、膜のアブレーションされた部分は何も蓄熱しません——スリット媒質（空気、真空、エアロゲル）はこの尺度では有意な蓄熱能力を持たないとみなされ、それが (1 − ψ) の係数です。次に、ρとc_pは300 Kの値として一定に扱われます——既定シナリオではプレートが20 °Cから120 °Cまで変化するにもかかわらず、温度依存性はモデル化されていません。
数値は共有サーバーカタログに由来します。
| 材料 | ρ (kg/m³) | c_p (J/kg·K) | データの位置づけ |
|---|---|---|---|
| Ginestium™ | 2,100 | 710 | 代替値 — 300Kにおける配向性合成熱分解炭素。Ginestium™自体の実測値ではない |
| グラフェンCVD few-layer | 2,200 | 700 | 文献値 |
| ダイヤモンドCVDナノ結晶 | 3,400 | 536 | NCD膜での実測値（302K） |
| 銅（対照） | 8,960 | 385 | 文献値、バルク銅 |
| 基板 | ユーザー入力（既定値 3,900 / 880） | — | |
「カスタム」材料にはρもc_pもありません。この場合、本ツールは値を推測する代わりに、「拡散率なし」と明示的に返し、マップも遅延値も表示しません。
厚みによる重み付けは膜の寄与を押し潰します。既定設定——Ginestium™ 400 nm、基板10 µm——では、膜が占めるのは(ρc)_sysのわずか約1.6%にすぎません。基板を2 mmにすると、その割合は0.01%未満まで下がります。実効拡散率は実質的に「システムの横方向伝導率を基板の熱容量で割ったもの」になります。
解釈における直接の帰結として、前線は膜のc_pに対してほとんど鈍感です。無視されているからではなく——式には確かに含まれています——積層構造の幾何学がその重みをほとんど与えないからです。膜が効くのは、式の分母にある比熱容量よりも、はるかに分子にあるk·tを通してです。その裏返しとして、基板について入力するρとc_pの値には、真に注意を払う価値があります——前線の速さを左右しているのは、まさにこれらの値だからです。
無視されているわけではありませんが、その扱いは部分的であり、面外方向に限られます。
実際に存在するもの。本ツールのR_th interfaceスライダー（0～100 ×10⁻⁶ m²·K/W、既定値10。熱異方性シミュレーターでも同じパラメータ）は、面外方向の抵抗に直列に加わります。
これは膜／基板間のただ1つのグローバルな界面です。幾何モデルは2層しか持たないため、界面ごとに抵抗を持つ多層積層モデルも、結晶粒間やパターン要素間の横方向抵抗も存在しません。
参考となるオーダー：グラフェン（k_z = 6 W/m·K、厚さ1.675 nm）の場合、膜の項は約3×10⁻¹⁰ m²·K/Wであり、典型的な界面抵抗より5桁小さい値です。R_zにおいては、常に界面と基板が支配的であり、膜が支配することはありません。
無視されているもの、そしてそれは意図的な前提です。R_th,ifは、k_xx/k_yy均質化にも、K_xxにも、したがってα_effにも一切関与しません——横方向拡散前線はこれを完全に無視します。本ツールは、この設定が実際にはマップに影響しないと誤解させないよう、該当するスライダーの下にこのことを明示的に表示しています。パラメトリックスイープに至っては、界面抵抗を一切持ちません——これは純粋に膜の面内のみを扱うツールです。
横方向輸送において、界面熱抵抗は熱の経路上に直列に入る抵抗ではありません——それは膜と基板という2つの並列経路の間の結合です。これを無視することは、z方向に完全な結合があると仮定すること——あらゆる位置xにおいて膜と基板が同じ局所温度を持つということ——と厳密に等価です。これは§03ですでに触れた等価平板の前提であり、モデルの形そのものと不可分です——式(4)が記述するのは1つの温度T(x, t)だけであり、膜用と基板用の2つの温度ではありません。
この前提が破綻する体制は特定できます——観測時間が短い場合、あるいは伝導性の高い膜が伝導性の低い基板の上にある場合です。その場合、膜は基板を引き連れる代わりに先行してしまい、センサーはモデルが予測するよりも早く熱の到達を検知することになります——このモデルが構造的に表現できない効果です。ここで対象としている厚みと観測時間については、この近似は妥当な範囲にとどまりますが、それでもこれは帰結ではなく前提であることに変わりはありません。