1 序論:様々な偏微分方程式

\[ \newcommand{\bm}[1]{\symbfit{#1}} \]

1.1 偏微分方程式の分類

一般に\(u = u(x, y)\)が独立変数\(x, y\)の関数であるとき,\(u\)\(x, y\)および偏導関数(\(\partial u/\partial x, \partial u/\partial y, \partial^2 u/\partial x \partial y, ...\))の間に成り立つ関係式を偏微分方程式と呼ぶ.物理的に興味のある系の多くは比較的低次の偏微分方程式で表される.よく現れる偏微分方程式として,例えば \[ a \frac{\partial^2 u}{\partial x^2} + b \frac{\partial^2 u}{\partial x \partial y} + c \frac{\partial^2 u}{\partial y^2} = F \left( u, x, y, \frac{\partial u}{\partial x}, \frac{\partial u}{\partial y} \right) \] の形で表される系を考えよう.このとき \(D = b^2 - 4 a c\)の値によって偏微分方程式は

  • \(D < 0\) : 楕円型
  • \(D = 0\) : 放物型
  • \(D > 0\) : 双曲型

と分類される.

この分類は2次曲線 \[ a x^2 + b xy + c y^2 = F(x, y) \]\(D\)の値によって楕円,放物線,双曲線に分類されることに由来するものであり,数学的にもこれらの分類によって解の性質は大きく異なる.

例1:Laplace方程式(楕円型方程式)

\(a = c = 1, b = 0, F = 0\)のとき, \[ \frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} = 0 \]

例2:拡散方程式(放物型方程式)

\(a = 1, b = c = 0, F = \partial u/\partial y\)\(y \rightarrow t\)と置き換えて \[ \frac{\partial u}{\partial t} = \frac{\partial^2 u}{\partial x^2} \]

例3:波動方程式(双曲型方程式)

\(a = 1, b = 0, c = -1, F = 0\)\(y \rightarrow t\)と置き換えて \[ \frac{\partial^2 u}{\partial t^2} = \frac{\partial^2 u}{\partial x^2} \]

1.2 初期条件と境界条件

  • 偏微分方程式の解は対象とする領域の境界で与えられる何らかの拘束条件(境界条件)によって決定される.変数の一つが時間(\(t\))であるときは,\(t = 0\)で与える境界条件を初期条件と呼ぶのが通例である.
  • 必要とされる境界条件の個数は微分の階数によって異なる.
  • 境界における解の値を与える境界条件をDirechlet型境界条件と呼ぶ.一方で境界における解の一階微分値を与える境界条件をNeumann型境界条件と呼ぶ.
  • 実用的には数値シミュレーションにおいて一番難しいのは適切な境界条件を設定することかもしれない.初心者は境界条件を適当に設定しがちであるが,物理的に意味のある境界条件を設定することは非常に重要である.

例1

一様重力場中の質点の運動方程式 \[ \frac{d^2 x}{d t^2} = -g \quad \left( \frac{dx}{d t} = v, \quad \frac{d v}{d t} = -g \right) \] の一般解は \[ x(t) = x_0 + v_0 t - \frac{1}{2} g t^2 \] である.この解を決定するには\(t = 0\)における条件(初期条件)\(x(0) = x_0\)および\(v(0) = v_0\)が必要.

この問題を境界値問題として見れば,\(t = 0\)および\(t = T\)の2点において一つずつ条件を与えることでも解を決定できる.

  • Direchlet型境界条件 \(x(0) = x_0\)\(x(T) = x_1\) を与える.

  • Neumann型境界条件 \(x'(0) = v_0\)\(x'(T) = v_1\) を与える.

偏微分方程式においても基本的に事情は同じであり,微分の階数に等しい条件を指定しなければならない.(境界条件を適当に与えると場合によっては解が無いということもあり得る!)

例2

波動方程式 \[ \frac{\partial^2 u}{\partial t^2} = \frac{\partial^2 u}{\partial x^2} \] は時間・空間ともに2階の偏微分方程式である. したがって,この方程式の解を決定するには初期条件として\(t=0\)における\(u\)および\(\partial u /\partial t\)の値,境界条件として両端での\(u\)の値を与える必要がある.

ここで,補助変数として\(\partial u/\partial t = \partial v/\partial x\)なる\(v\)を導入してみよう.このとき波動方程式は \[ \frac{\partial u}{\partial t} = \frac{\partial v}{\partial x}, \quad \frac{\partial v}{\partial t} = \frac{\partial u}{\partial x} \] のように,2変数\(u, v\)についての時間・空間ともに1階の偏微分方程式に書き換えられる. したがって,初期条件として\(t = 0\)における\(u, v\)の値が必要となる.一方で,境界条件としては,例えば両方の境界で\(u\)の値を与えてもよいし,一方の境界で\(u\)を,もう一方の境界で\(v\)を与えてもよい.
ただし,両端で\(u, v\)の両方の値を同時に指定することはできない(自由度が足りない)ことに注意しよう.物理的には右方向に伝播する波動については左側の境界が,左方向に伝播する波動については右側の境界がそれぞれ影響を与えるためと考えることができる.