Laplace方程式

すでに学んだように,次の形の偏微分方程式をそれぞれ二次元,三次元のLaplace方程式といいます.

$\displaystyle \Delta_{2}u = u_{xx} + u_{yy} = 0,  \Delta_{3}u = u_{xx} + u_{yy} + u_{zz} = 0. $

これは非圧縮性流体の速度ポテンシャル,静電場のポテンシャル,熱伝導現象における定常温度分布などを表わす方程式です.$C^{2}$級の関数で,Laplace方程式をみたすものを調和関数(harmonic function)といいます.

二次元のLaplace方程式の解は次の形のいずれかを取ることを学びましたが,次の形の関数の有限個の一次結合でLaplace方程式の境界値問題の解を表わすのが適当でないことがよくあります.

$\displaystyle (A\cosh{kx} + B\sinh{kx})(C\cos{ky} + D\sin{ky}), $

$\displaystyle (A\cos{kx} + B\sin{kx})(C\cosh{ky} + D\sinh{ky}), $

$\displaystyle (Ax + B)(Cy + D). $

そんなときは,無限級数を用いて解を表わし,フーリエ級数の知識をもとに問題をときます.

例題 7..11  

2次元Laplace方程式の境界値問題

$\displaystyle u_{xx} + u_{yy} = 0  (0 < x < L, 0 < y < K), $

$\displaystyle u(x,0) = 0 ,  u(x,K) = x  (0 \leq x \leq L), $

$\displaystyle u_{x}(0,y) = 0 ,  u_{x}(L,y) = 0  (0 \leq y \leq K) $

を解け.

変数分離法と境界条件 $u_{x}(0,y) = 0, u(x,0) = 0$よりLaplace方程式の解は

$\displaystyle C\cosh{kx}\cos{ky},  C\cos{kx}\sinh{ky},  Cy $

のどれかである.3番目の形は境界条件 $u_{x}(L,y) = 0$をすべての$C$でみたし,また2番目の形も $\sin{kL} = 0$のとき,つまり $k = n\pi/L  (n = 1,2,\ldots)$のとき境界条件 $u_{x}(L,y) = 0$をすべての$C$でみたす.よって重ね合わせの原理より,

$\displaystyle u(x,y) = \frac{C_{0}}{2}y + \sum_{n=1}^{\infty}C_{n}\cos{\frac{n\pi x}{L}}\sinh{\frac{n\pi y}{L}} $

を境界値問題の解として考える.最後に境界条件 $u(x,K) = x$より

$\displaystyle u(x,K) = \underbrace{C_{0}K/2}_{a_{0}/2} + \sum_{n=1}^{\infty}\underbrace{C_{n}\sinh{\frac{n\pi K}{L}}}_{a_{n}}\cos{\frac{n\pi x}{L}} = x. $

この条件が満たされるためには,この級数が$x$に収束するように$a_{n}$を選ばなければならない.ところが,これは本質的にフーリエ余弦級数である.よって

$\displaystyle a_{0} = C_{0}K = \frac{2}{L}\int_{0}^{L}xdx = \frac{2}{L}\frac{x^{2}}{2}\mid_{0}^{L} = L, $


$\displaystyle a_{n} = C_{n}\sinh{\frac{n\pi K}{L}}$ $\displaystyle =$ $\displaystyle \frac{2}{L}\int_{0}^{L}x\cos{\frac{n\pi x}{L}}dx$  
  $\displaystyle =$ $\displaystyle \frac{2L}{n^{2}\pi^{2}}[(-1)^{n} - 1] = \left\{\begin{array}{cl}
...
...}{n^{2}\pi^{2}}, & n  \mbox{odd} \\
0, & n  \mbox{even}.
\end{array}\right .$  

よって
$\displaystyle u(x,y)$ $\displaystyle =$ $\displaystyle \frac{Ly}{2K}$  
  $\displaystyle -$ $\displaystyle \sum_{m=0}^{\infty}\frac{4L}{\pi^{2}}\frac{1}{(2m+1)^{2}}\frac{\cos{(2m+1)\pi x/L}\sinh{(2m+1)\pi y/L}}{\sinh{(2m+1)\pi K/L}}$  

となる. $ \blacksquare$