拡散方程式の一般解とは？求め方と解の性質を解説！（特解・同次方程式・非同次方程式・変数分離・固有関数など）というテーマでは、まず拡散方程式の解が何を意味するのかを押さえることが大切です。
拡散方程式は、温度、濃度、密度などの分布が時間とともに広がり、なめらかになっていく現象を表す偏微分方程式です。
その一般解は、初期条件や境界条件に応じて変わり、単純な一つの式だけで完結するとは限りません。
特に有限区間の問題では、変数分離法によって固有関数を求め、それらを重ね合わせることで一般解を表します。
非同次方程式では、同次方程式の一般解に特解を加える考え方が重要になります。
拡散方程式の一般解は、基本的に境界条件に合う固有関数と時間減衰する指数関数の重ね合わせとして表されます。
本記事では、拡散方程式の一般解の意味、変数分離法による求め方、同次方程式と非同次方程式の違い、特解の考え方、解の性質までわかりやすく解説していきます。
拡散方程式の一般解は境界条件に合う固有関数の重ね合わせで表されます
それではまず、拡散方程式の一般解の結論に相当する内容について解説していきます。
一般解はすべての解を含む形です
一般解とは、与えられた方程式を満たす解をできるだけ広く含む形のことです。
ただし、拡散方程式のような偏微分方程式では、一般解という言葉の意味が常微分方程式より少し複雑になります。
なぜなら、拡散方程式の解は初期条件や境界条件によって大きく変わるからです。
たとえば、同じ一次元の拡散方程式でも、両端の温度をゼロに固定する場合と、両端を断熱する場合では、使う固有関数が変わります。
また、無限領域で考える場合には、ガウス関数を用いた基本解や畳み込み積分が登場します。
有限区間では、正弦関数や余弦関数の級数として解を表すことが多いです。
つまり、拡散方程式の一般解は、問題設定とセットで考える必要があります。
拡散方程式では、方程式だけでなく、初期条件と境界条件を含めて解の形が決まります。
有限区間では固有関数展開を使います
有限区間の拡散方程式では、変数分離法によって解を求めるのが代表的です。
一次元の拡散方程式 ∂u/∂t = D∂²u/∂x² を考えます。
ここで、解を u(x,t) = X(x)T(t) と仮定します。
この形を方程式に代入すると、空間だけの式と時間だけの式に分けることができます。
空間方向の式では、境界条件に合う関数を求めます。
この関数が固有関数です。
たとえば、区間 0
時間方向の式では、各固有関数に対応して指数関数的に減衰する項が得られます。
両端がゼロの一次元拡散方程式では、u(x,t) = Σ Bn sin(nπx/L) exp(-D(nπ/L)²t) の形で解を表せます。
係数 Bn は初期条件から決まります。
このように、固有関数と時間減衰項を重ね合わせることで、初期分布に合う解を作ります。
各モードは時間とともに減衰します
拡散方程式の解では、固有関数ごとの成分をモードと呼ぶことがあります。
各モードは、時間とともに指数関数的に減衰します。
高いモードほど空間的な変化が細かく、減衰が速くなります。
低いモードほどゆっくり減衰し、長時間後まで残りやすくなります。
これは拡散方程式の重要な性質です。
初期分布に急な変化や細かな凹凸が含まれていても、それらは時間とともにすばやく消えていきます。
その結果、解はだんだんなめらかな形になります。
長時間後には、最も減衰の遅い低次モードが支配的になることがあります。
この性質により、拡散方程式の解は時間とともに平滑化されるといえます。
拡散方程式の一般解では、細かい空間変動ほど早く消えます。
そのため、時間が進むほど解はなめらかになり、初期分布の細かな形は失われていきます。
変数分離法では空間の式と時間の式に分けて解きます
続いては、拡散方程式の一般解を求める代表的な方法である変数分離法を確認していきます。
解を位置の関数と時間の関数の積と仮定します
変数分離法では、解 u(x,t) を位置だけに依存する関数 X(x) と、時間だけに依存する関数 T(t) の積として仮定します。
つまり、u(x,t) = X(x)T(t) と置きます。
この仮定は一見すると限定的に見えますが、得られた解を重ね合わせることで広い初期条件に対応できます。
線形の拡散方程式では、解の和もまた解になります。
この重ね合わせの原理があるため、個別の分離解を組み合わせて一般的な解を作ることができます。
u = XT を方程式に代入すると、X と T を含む式になります。
両辺を XT で割ると、時間だけの式と空間だけの式が等しい形になります。
時間だけに依存する量と空間だけに依存する量が常に等しいため、それらは定数に等しいと考えます。
この定数を分離定数と呼びます。
空間方向の問題から固有値と固有関数が得られます
変数分離法で重要なのが、空間方向の固有値問題です。
境界条件を満たす X(x) を求めることで、許される固有値と固有関数が決まります。
たとえば、両端で u = 0 となる境界条件では、X(0) = 0、X(L) = 0 を満たす必要があります。
この条件を満たす代表的な関数が sin(nπx/L) です。
このとき、n は正の整数であり、それぞれの n に対応する固有値が決まります。
固有値は、時間方向の減衰速度にも関係します。
大きな固有値に対応するモードほど速く減衰します。
境界条件が変われば、固有関数も変わります。
たとえば、断熱条件では余弦関数が現れることがあります。
時間方向の式は指数関数になります
空間方向の固有値が決まると、時間方向の式も解けます。
拡散方程式の場合、時間方向の解は指数関数的に減衰する形になります。
たとえば、固有値が λ の場合、時間項は exp(-Dλt) のような形になります。
D は拡散係数です。
拡散係数が大きいほど、減衰は速くなります。
固有値が大きいほど、やはり減衰は速くなります。
このため、細かな空間変化を表す高次モードは、短時間で小さくなります。
低次モードはゆっくり減衰し、長時間後の解に影響しやすくなります。
変数分離法の結果、拡散方程式の解は、空間的な形と時間的な減衰を組み合わせたものとして理解できます。
| 手順 | 内容 | ポイント |
|---|---|---|
| 解の仮定 | u(x,t) = X(x)T(t) と置きます | 位置と時間を分けます |
| 代入 | 拡散方程式へ代入します | 空間式と時間式に分離します |
| 固有値問題 | 境界条件に合う X を求めます | 固有関数が決まります |
| 時間項 | T(t) を指数関数で求めます | 各モードは減衰します |
| 重ね合わせ | 各モードを足し合わせます | 初期条件に合わせます |
同次方程式と非同次方程式では一般解の作り方が変わります
続いては、同次方程式と非同次方程式における一般解の違いを確認していきます。
同次方程式は外部からの発生項がない場合です
同次の拡散方程式とは、外部からの発生項や強制項がない方程式です。
代表的には、∂u/∂t = D∂²u/∂x² のような形です。
この場合、分布は時間とともに拡散し、境界条件に応じた定常状態やゼロ状態へ近づきます。
両端をゼロに固定した同次問題では、すべてのモードが時間とともに減衰し、最終的にゼロへ近づきます。
断熱境界の場合には、全体の平均値が保存され、最終的に一様な分布へ近づくことがあります。
同次方程式の一般解は、固有関数と指数減衰項の重ね合わせで表されることが多いです。
初期条件をフーリエ級数のように展開し、その係数によって各モードの大きさを決めます。
この考え方は、熱伝導や濃度拡散の基礎問題でよく使われます。
非同次方程式では特解が必要です
非同次の拡散方程式とは、外部からの発生項や強制項を含む方程式です。
たとえば、内部発熱がある熱伝導問題では、右辺に発熱項が加わります。
濃度拡散でも、化学反応や物質の供給がある場合には、非同次項が現れることがあります。
非同次方程式では、一般解は同次方程式の一般解と特解の和で表されます。
非同次方程式の一般解は、一般解 = 同次解 + 特解 という形で考えます。
同次解は境界条件に合う自由な減衰成分、特解は外部からの強制項に対応する成分です。
特解は、非同次項を満たすために必要な特別な解です。
特解の求め方は、非同次項の形によって変わります。
時間に依存しない発熱であれば、まず定常解を特解として求めることがあります。
時間に依存する強制項であれば、畳み込み積分やグリーン関数を使う場合もあります。
定常解と過渡解に分けると理解しやすいです
非同次問題では、解を定常解と過渡解に分けると見通しがよくなります。
定常解は、時間が十分に経った後に残る成分です。
過渡解は、初期条件とのずれを調整する成分で、時間とともに減衰します。
たとえば、一定の内部発熱がある棒では、長時間後に一定の温度分布が形成されます。
この温度分布が定常解です。
最初の温度分布が定常解と違う場合、その差が過渡解として時間とともに消えていきます。
この考え方を使うと、非同次方程式の解を物理的に理解しやすくなります。
非同次拡散方程式では、外部から与えられる影響を特解が表し、初期条件による一時的な変化を同次解が表します。
拡散方程式の解には平滑化と最大値原理という重要な性質があります
続いては、拡散方程式の解が持つ代表的な性質について確認していきます。
解は時間とともになめらかになります
拡散方程式の解は、時間とともになめらかになる性質を持っています。
初期分布に急な変化や不連続に近い形があっても、時間が少し進むと分布は平滑化されます。
これは、高次モードが速く減衰することと関係しています。
細かな凹凸ほど早く消えるため、全体として滑らかな形になるのです。
物理的には、濃度や温度の急な差が強い流れを生み、その差をすばやく小さくするという意味です。
この平滑化性は、拡散方程式の大きな特徴です。
波動方程式のように形が伝わるのではなく、拡散方程式では形がぼやけて広がります。
そのため、長時間後には初期分布の細かい情報が失われやすくなります。
最大値原理により極端な値は自然に抑えられます
拡散方程式には最大値原理と呼ばれる性質があります。
簡単にいうと、内部で突然初期値や境界値を超えるような極端な最大値が生まれにくいという性質です。
熱伝導で考えると、外部から熱源がないのに、物体内部の温度が最初や境界よりも急に高くなることは自然ではありません。
拡散は高いところから低いところへ広がるため、極端な山は削られ、谷は埋められます。
この性質により、解は安定した振る舞いを示します。
最大値原理は、拡散方程式の理論的な解析でも重要です。
解の一意性や安定性を示すときにも使われます。
物理的にも、拡散がむらを増やすのではなく、むらを減らす現象であることを示しています。
保存量がある場合とない場合があります
拡散方程式では、境界条件によって全体量が保存される場合と保存されない場合があります。
たとえば、断熱境界や物質が外へ出ない境界では、全体の熱量や物質量が保存されることがあります。
この場合、分布は広がりますが、総量は変わりません。
最終的には、保存された総量が全体に均一に分布する形へ近づきます。
一方、境界から熱や物質が出入りする場合には、全体量は保存されません。
両端の温度をゼロに固定する熱伝導問題では、内部の熱は境界から逃げ、最終的に温度がゼロへ近づきます。
このように、解の長時間挙動は境界条件に大きく左右されます。
拡散方程式の一般解を見るときは、保存量があるかどうかも重要な視点です。
拡散方程式の解は、時間とともに細かな変化を失い、なめらかで安定した分布へ向かいます。
境界条件によって、総量が保存される場合と外へ逃げる場合があるため、長時間後の姿も変わります。
まとめ
拡散方程式の一般解は、問題設定、初期条件、境界条件によって形が変わります。
有限区間の同次問題では、変数分離法によって固有関数を求め、それらを時間減衰する指数関数と組み合わせて解を表します。
代表的には、正弦関数や余弦関数の級数展開として一般解が得られます。
係数は初期条件から決まり、各モードは時間とともに指数関数的に減衰します。
高次モードほど速く減衰するため、解は時間とともになめらかになります。
非同次方程式では、同次解に特解を加えることで一般解を作ります。
特解は外部からの発熱、供給、強制項などに対応し、同次解は初期条件とのずれを調整する過渡成分になります。
拡散方程式の解には、平滑化、最大値原理、長時間での定常状態への接近といった重要な性質があります。
拡散方程式の一般解を理解するポイントは、固有関数の重ね合わせ、指数減衰、特解と同次解の分解を押さえることです。
この考え方を身につけると、熱伝導、濃度拡散、物質移動、確率分布の広がりなど、さまざまな現象を共通した数学の視点で理解しやすくなるでしょう。
