拡散方程式を理解したら、次は実際に解いてみたいと思う方も多いのではないでしょうか。
拡散方程式の解き方にはいくつかのアプローチがあり、問題の設定（初期条件・境界条件・領域の形）によって最適な解法が異なります。
本記事では、変数分離法・フーリエ級数展開・フーリエ変換・数値解法という代表的な解法を、具体的な手順とともにわかりやすく解説していきます。
目次
拡散方程式の解き方の全体像：解法の選択基準
それではまず、拡散方程式のさまざまな解法と、どの場合にどの解法を使うべきかの基準について解説していきます。
| 解法 | 適した条件 | 得られる解の形 |
|---|---|---|
| 変数分離法 | 有限領域・単純な境界条件 | フーリエ級数（固有関数展開） |
| フーリエ変換法 | 無限領域・任意の初期条件 | ガウス関数との畳み込み積分 |
| ラプラス変換法 | 半無限領域・特殊な境界条件 | 補誤差関数を含む表式 |
| 有限差分法 | 一般的な領域・数値解が必要 | 格子点での数値解 |
問題の定式化：初期条件と境界条件の確認
拡散方程式を解く際の最初のステップは、問題を正確に定式化することです。
初期条件（t=0での空間分布）と境界条件（空間端での条件）を明確にすることが求解の大前提です。
たとえば「[0,L]上のディリクレ問題」は次のように定式化されます。
∂u/∂t = D ∂²u/∂x²（方程式）
u(0,t) = u(L,t) = 0（ディリクレ境界条件）
u(x,0) = f(x)（初期条件）
斉次化と特解の分離
境界条件が非斉次（ゼロでない）場合は、まず斉次化（境界をゼロにする変換）を行います。
u = v + w として、wは非斉次境界条件を満たす定常解、vは斉次境界条件を持つ問題の解と分離します。
この分離によって問題が扱いやすくなります。
対称性の活用
問題が空間的対称性を持つ場合（偶関数・奇関数の初期条件など）、それを活用すると計算量を大幅に削減できます。
余弦展開・正弦展開の選択は、境界条件の種類（ノイマン・ディリクレ）に対応しています。
変数分離法による解き方の詳細手順
続いては、変数分離法の具体的な手順を一ステップずつ確認していきます。
ステップ1：解を積の形に仮定する
変数分離法では、u(x,t) = X(x)T(t) という積の形の解が存在すると仮定します。
拡散方程式に代入すると次のようになります。
X(x)T'(t) = D X”(x)T(t)
両辺をX(x)T(t)で割ると：T'(t)/(DT(t)) = X”(x)/X(x) = -λ（分離定数）
ステップ2：常微分方程式を解く
分離定数λを用いると二つの常微分方程式が得られます。
T'(t) = -λD T(t) → T(t) = C exp(-λDt)
X”(x) = -λ X(x) → X(x) = A sin(√λ x) + B cos(√λ x)
ディリクレ境界条件 X(0) = X(L) = 0 を適用すると、固有値 λₙ = (nπ/L)² が離散的に定まります。
ステップ3：フーリエ係数の決定
一般解は固有関数の重ね合わせです。
u(x,t) = Σ_{n=1}^{∞} Bₙ sin(nπx/L) exp(-D(nπ/L)²t)
係数Bₙは初期条件 u(x,0) = f(x) からフーリエ正弦展開で決まります。
Bₙ = (2/L) ∫_0^L f(x) sin(nπx/L) dx
数値解法：有限差分法の手順
続いては、複雑な問題に対して強力な数値解法である有限差分法の具体的手順を確認していきます。
離散化と差分スキームの構成
有限差分法では空間をΔxの格子、時間をΔtの格子で離散化します。
最も単純な陽解法（前進オイラー法）では次の差分スキームを使います。
u_i^{n+1} = u_i^n + r(u_{i+1}^n – 2u_i^n + u_{i-1}^n)
r = DΔt/Δx²（クーラント数）
安定性条件（CFL条件）
陽解法では、安定に計算するための条件があります。
CFL条件（安定性条件）：r = DΔt/Δx² ≤ 1/2
この条件を満たさないと数値解が発散してしまいます。安定な計算のためにΔtとΔxを適切に選ぶことが重要です。
より安定な陰解法（クランク-ニコルソン法）は無条件安定ですが、各時間ステップで連立方程式を解く必要があります。
数値解の精度評価
有限差分法の精度は格子幅Δxと時間刻みΔtに依存します。
前進時間中心空間差分（FTCS）スキームは時間1次・空間2次精度です。
格子幅を半分にすると空間誤差は1/4になるため、精度要求に応じて格子を細かくすることが必要でしょう。
まとめ
本記事では、拡散方程式の解き方として変数分離法・フーリエ変換法・有限差分法の手順を解説しました。
解析的解法は問題の構造理解に、数値解法は実用的な複雑問題に、それぞれ強みがあります。
問題の領域・境界条件・精度要求に応じて最適な解法を選択することが重要で、複数の手法を組み合わせることも実際の工学問題では一般的です。
まずは変数分離法で基本を固め、フーリエ変換法・数値解法へと順番に習得していくことをおすすめします。
