拡散方程式の意味や基本概念を理解したら、次に気になるのは「実際にどのような解を持つのか」という点ではないでしょうか。
拡散方程式の一般解は、フーリエ変換やグリーン関数といった強力な数学ツールを使って表現されます。
本記事では、拡散方程式の基本解・一般解・解の特徴を、数式の意味を丁寧に追いながら詳しく解説していきます。
目次
拡散方程式の一般解とは？基本解とグリーン関数による表現
それではまず、拡散方程式の一般解の意味と、基本解・グリーン関数を用いた表現について解説していきます。
拡散方程式の一般解は、初期条件をグリーン関数（基本解）で畳み込み積分することで表されます。
1次元の無限領域における拡散方程式の基本解（グリーン関数）は次のガウス関数です。
G(x, t) = (1/√(4πDt)) exp(-x²/(4Dt))
これはt=0においてデルタ関数δ(x)を初期条件とした場合の解です。
この基本解は時間とともに広がるガウス分布（正規分布）の形をしており、拡散現象の本質を端的に表しています。
基本解の特徴と時間発展
基本解G(x,t)は時間が経つにつれて幅が広がり、高さが低くなります。
広がりの特徴量である標準偏差σは次のように表されます。
σ = √(2Dt)
これは「拡散距離がtの平方根に比例する」ことを意味します。
これは拡散の重要な特徴であり、波動（距離∝時間）とは異なる性質です。
たとえば拡散係数D = 1 m²/s のとき、1秒後の広がりは約1.4 m、4秒後は約2.8 mとなるでしょう。
一般解の畳み込み積分表現
任意の初期条件 u(x, 0) = u₀(x) に対する一般解は、基本解との畳み込みで表されます。
u(x, t) = ∫_{-∞}^{∞} G(x-y, t) u₀(y) dy
= (1/√(4πDt)) ∫_{-∞}^{∞} exp(-(x-y)²/(4Dt)) u₀(y) dy
これは「初期条件の各点からの寄与をガウス関数で重み付けして足し合わせる」という物理的に直感的な表現です。
解の一意性と初期条件の役割
適切な成長条件の下で、拡散方程式の解は初期条件によって一意に決まります。
初期条件の形の違いが解の時間発展にどう影響するかを理解することは、拡散現象の物理的直感を養ううえで重要です。
特殊な初期条件（矩形波・正弦波・デルタ関数など）に対して解析的に解が求まる場合は、それぞれの解の特徴を比較するとよいでしょう。
フーリエ変換を用いた拡散方程式の一般解の導出
続いては、フーリエ変換を使った一般解の導出法を確認していきます。
フーリエ変換は、空間変数xに関する微分をk（波数）の代数演算に変換し、偏微分方程式を常微分方程式に簡略化します。
フーリエ変換の適用手順
u(x,t)のフーリエ変換を û(k,t) = ∫ u(x,t) e^(-ikx) dx とすると、∂²u/∂x² は -k² û に対応します。
拡散方程式はk空間では次の常微分方程式に変換されます。
∂û/∂t = -Dk² û
解：û(k,t) = û(k,0) exp(-Dk²t)
k空間では各フーリエ成分が指数関数的に減衰することがわかります。
高波数成分（細かい空間変動）ほど速く減衰し、低波数成分（なだらかな変動）は長く残ります。
逆フーリエ変換で実空間に戻す
û(k,t)を逆フーリエ変換すると実空間での解u(x,t)が求まります。
exp(-Dk²t)の逆フーリエ変換がガウス関数 (1/√(4πDt))exp(-x²/(4Dt)) になることを利用すると、先ほどの畳み込み積分の一般解が導かれます。
これによりフーリエ変換法と基本解法の整合性が確認できます。
フーリエ級数展開（有限領域の場合）
有限区間 [0, L] 上の拡散方程式では、フーリエ変換の代わりにフーリエ級数展開を用います。
u(x,t) = Σ_{n=1}^{∞} Bₙ sin(nπx/L) exp(-D(nπ/L)²t)
Bₙ：初期条件から決まるフーリエ係数
各項は時間とともに指数関数的に減衰し、長時間では最低次（n=1）の項が支配的になります。
解の特徴と物理的解釈
続いては、拡散方程式の解が持つ重要な特徴と、その物理的な意味を確認していきます。
情報の瞬時伝播と実際の限界
拡散方程式の数学的な解は、初期条件が有限サポートを持つ場合でも、t>0 では空間全体で非ゼロの値を取ります。
これは情報が瞬時に全空間に伝わることを意味し、特殊相対性理論の光速不変の原理と矛盾するように見えます。
実際には、伝播速度が非常に速いだけで、遠方への影響は指数的に小さく実用上は問題になりません。
より精密な記述には、相対論的拡散方程式（電信方程式）が使われることがあります。
エントロピーと不可逆性
拡散方程式は時間反転に対して非対称であり、物理的な不可逆過程を表します。
時間を反転させた方程式（∂u/∂t = -D∂²u/∂x²）は不安定であり、拡散の逆過程（「集中」）は自然には起こらないことを示しています。
これは熱力学の第二法則（エントロピー増大の原則）と整合しています。
保存則と全量の保存
拡散方程式では、空間全体での全量（∫u dx）が保存されます。
d/dt ∫_{-∞}^{∞} u dx = D ∫_{-∞}^{∞} ∂²u/∂x² dx = 0
（境界での消滅がない場合）
これは質量保存・エネルギー保存に対応する物理的に重要な性質です。
物質拡散では総物質量が保たれ、熱伝導では断熱系での総熱エネルギーが保存されます。
まとめ
本記事では、拡散方程式の基本解・一般解（畳み込み積分とフーリエ変換による導出）・解の特徴について解説しました。
基本解はガウス分布の形を持ち、拡散距離が√(Dt)に比例するという本質的な性質を表します。
フーリエ変換を使うと各波数成分の減衰率が明確になり、解の構造を視覚的に理解しやすくなるでしょう。
保存則・不可逆性・最大値原理など、解の持つ物理的特徴も合わせて理解することで、拡散現象への洞察が深まるはずです。
