「拡散」という現象は、インクが水の中で広がったり、熱が伝わったりする日常的な現象に密接に関わっています。
この「拡散」という現象を数学的に記述するのが拡散方程式です。
拡散方程式は物理・化学・生物・金融工学など非常に幅広い分野で活躍する偏微分方程式であり、その基本概念を押さえることは理工系の学習において欠かせないでしょう。
本記事では、拡散方程式の意味・導出のアイデア・熱伝導との関係・物質拡散への応用までをわかりやすく解説していきます。
目次
拡散方程式とは何か？基本的な定義と意味
それではまず、拡散方程式の定義と基本的な意味について解説していきます。
拡散方程式とは、ある物理量（温度・濃度など）が時間とともに空間的に広がっていく現象を記述する偏微分方程式です。
最も基本的な1次元の拡散方程式は次のように表されます。
∂u/∂t = D ∂²u/∂x²
u：温度または濃度などの物理量
t：時間、x：空間座標
D：拡散係数（正の定数）
この方程式は「物理量の時間変化率が空間の二階微分に比例する」ことを表しています。
数学的には放物型偏微分方程式に分類され、ポアソン方程式（楕円型）や波動方程式（双曲型）と並ぶ代表的な偏微分方程式のひとつです。
拡散係数Dの物理的意味
拡散方程式における拡散係数Dは、拡散の速さを表すパラメータです。
Dが大きいほど物理量が素早く広がり、Dが小さいほどゆっくりと広がります。
熱伝導の場合、Dは熱拡散率（熱伝導率・密度・比熱から決まる量）であり、物質の種類によって異なる値を取ります。
D = κ/(ρ c_p)
κ：熱伝導率、ρ：密度、c_p：定圧比熱
Dの単位はm²/sであり、長さの二乗を時間で割った次元を持ちます。
放物型方程式としての性質
拡散方程式が放物型に分類される意味は、解の性質に深く関わっています。
放物型方程式の解は、初期条件が与えられると時間が進むにつれて滑らかになっていく性質があります。
これは波動方程式（不連続が保たれる）とは対照的な性質でしょう。
また、情報が瞬時に（無限大の速度で）空間全体に伝わるという特徴もあります（熱は瞬時にすべての場所に影響を及ぼすという意味で、実際には近似ですが）。
熱方程式としての表現
拡散方程式は熱伝導の文脈では熱方程式（heat equation）とも呼ばれます。
3次元での表現は次のようになります。
∂T/∂t = D ΔT
T：温度、Δ：ラプラシアン（∂²/∂x² + ∂²/∂y² + ∂²/∂z²）
時間が十分に経過して定常状態（∂T/∂t = 0）に達すると、ΔT = 0 すなわちラプラス方程式に帰着します。
これはポアソン方程式・ラプラス方程式との重要なつながりです。
拡散方程式の物理的背景と応用分野
続いては、拡散方程式がどのような物理的原理に基づき、どの分野で活用されているかを確認していきます。
フィックの法則と物質拡散
物質拡散の分野では、フィックの第一法則と第二法則が拡散方程式の基礎となっています。
フィックの第一法則は「拡散フラックスは濃度勾配に比例する」という法則です。
J = -D ∇c
J：拡散フラックス（単位面積・単位時間当たりの物質量）
c：濃度、D：拡散係数
フィックの第二法則は、連続の式（質量保存則）とフィックの第一法則を組み合わせることで導かれる拡散方程式です。
∂c/∂t = D Δc
これは熱方程式と全く同じ数学的構造を持っています。
多様な応用分野
拡散方程式はさまざまな分野で活用されています。
| 分野 | 拡散方程式が記述するもの | 具体例 |
|---|---|---|
| 熱工学 | 熱の伝わり方 | 電子部品の冷却設計 |
| 化学・化学工学 | 物質の拡散・混合 | 膜透過・反応拡散系 |
| 生物学 | 細胞内物質の輸送 | 酸素の組織内拡散 |
| 金融工学 | オプション価格の変動 | ブラック-ショールズ方程式 |
| 材料科学 | 半導体中の不純物拡散 | ドーピングプロセス |
ブラウン運動との関係
拡散方程式は確率論とも深く関わっています。
ランダムウォーク（ブラウン運動）する粒子の確率密度関数が、拡散方程式を満たすことが数学的に示されます。
この関係はアインシュタインが1905年に示したもので、拡散方程式の確率論的解釈として非常に重要です。
金融工学でのブラック-ショールズ方程式もこの延長線上にあります。
拡散方程式の初期条件・境界条件と解の性質
続いては、拡散方程式を解くために必要な初期条件・境界条件の設定と、解の特徴的な性質を確認していきます。
初期条件の設定
拡散方程式は時間を含む方程式であるため、初期条件（t=0 における物理量の空間分布）が必要です。
u(x, 0) = u₀(x)
u₀(x)：初期分布（既知の関数）
初期条件が与えられると、その後の時間発展が方程式によって決定されます。
典型的な初期条件として、デルタ関数（点源）・矩形分布・正弦波などがよく使われます。
境界条件の種類と設定
有限の領域で拡散方程式を解く場合、空間的な境界条件も必要です。
ディリクレ境界条件（境界での値を固定）・ノイマン境界条件（境界でのフラックスを指定）・周期境界条件などが代表的です。
断熱壁では法線方向の熱流束がゼロ（ノイマン条件）、一定温度の壁ではディリクレ条件が使われるでしょう。
解の滑らか化と最大値原理
拡散方程式の解は、時間が経過するにつれて滑らかになっていく性質を持ちます。
初期条件が不連続であっても、t>0 では解は無限回微分可能な滑らかな関数になります。
また、最大値原理が成立し、任意の時刻での最大値は初期条件の最大値を超えることはありません。
これは「熱は源なしに自然に高まることはない」という物理的直感と一致しています。
まとめ
本記事では、拡散方程式の意味・拡散係数の物理的意義・フィックの法則との関係・応用分野・解の性質について解説しました。
拡散方程式は熱伝導・物質拡散・金融工学・生物学など幅広い分野で活躍する放物型偏微分方程式です。
初期条件と境界条件を正しく設定し、方程式の数学的性質を理解することが、実際の問題を解くうえでの土台となるでしょう。
次のステップとして、一般解の導出や具体的な解法（変数分離法・フーリエ変換）に取り組むと、理解がより深まるはずです。
