ポアソン方程式は電磁気学・熱伝導・流体力学など多くの分野で登場する重要な偏微分方程式ですが、「実際にどうやって解くのか」が気になる方も多いのではないでしょうか。
解法には大きく分けて解析的手法と数値的手法があり、問題の形状や条件によって使い分けが必要です。
本記事では、変数分離法・グリーン関数法・有限差分法を中心に、ポアソン方程式の求解手順をわかりやすく解説していきます。
目次
ポアソン方程式の解き方：基本方針と解法の選択
それではまず、ポアソン方程式の解法の全体像と、どのような場合にどの解法を選ぶべきかについて解説していきます。
ポアソン方程式 Δu = f を解くには、問題の領域形状・境界条件・ソース項の形によって最適な解法が異なります。
| 解法 | 適した条件 | 特徴 |
|---|---|---|
| 変数分離法 | 単純な領域・特殊な境界条件 | 解析的・厳密解が得られる |
| グリーン関数法 | 自由空間・特殊領域 | 任意のソース項に対応 |
| フーリエ変換法 | 無限領域・周期的な問題 | 解析的・変換域での計算 |
| 有限差分法 | 一般的な領域・複雑なソース | 数値的・汎用性が高い |
| 有限要素法 | 複雑な形状の領域 | 数値的・形状の自由度が高い |
境界条件の確認と問題の定式化
どの解法を選ぶにせよ、まず行うべきは境界条件の明確な設定です。
ディリクレ境界条件（境界での値の指定）、ノイマン境界条件（法線方向微分の指定）、またはその混合を確認します。
境界条件が不適切だと、解が存在しなかったり一意に定まらなかったりするため、定式化の段階が非常に重要でしょう。
また、ソース項f(x)の形（定数・関数・デルタ関数など）も解法の選択に大きく影響します。
解の重ね合わせと特解・斉次解
ポアソン方程式の一般解は、特解（ポアソン方程式を満たす一つの解）と斉次解（ラプラス方程式の一般解）の和として表せます。
u = u_p + u_h
u_p：ポアソン方程式Δu = fの特解
u_h：ラプラス方程式Δu = 0の解（斉次解）
u_hは境界条件を満たすように決定されます。
この分解は解法を整理するうえで非常に有効な考え方です。
解の存在・一意性の確認
実際に解を求める前に、解が存在し一意であることを確認することが望ましいでしょう。
ディリクレ境界条件の場合、ポアソン方程式の解は一意に存在することが数学的に保証されています。
ノイマン境界条件では、整合条件（ガウスの発散定理から導かれる条件）が必要であり、解は定数の任意性を残します。
変数分離法によるポアソン方程式の解き方
続いては、解析的手法の代表である変数分離法を用いたポアソン方程式の解き方を確認していきます。
変数分離法は、解を複数の変数の積の形に仮定し、偏微分方程式を常微分方程式の組に分解する方法です。
変数分離法の基本手順
2次元矩形領域 [0,a]×[0,b] 上のポアソン方程式を例に手順を説明します。
∂²u/∂x² + ∂²u/∂y² = f(x,y)
境界条件：u = 0（全境界上）
まず固有関数展開を利用します。ラプラス演算子の固有関数（ディリクレ境界条件の場合はsin関数）を用いて、uとfをそれぞれフーリエ級数展開します。
展開係数を比較することで、各モードの係数が代数的に求まります。
この方法は領域が矩形・円・球などの単純な形状のときに有効です。
円形領域での変数分離（極座標系）
円形領域では極座標系（r, θ）を用いた変数分離が有効です。
ラプラシアンを極座標で表すと次のようになります。
Δu = ∂²u/∂r² + (1/r)∂u/∂r + (1/r²)∂²u/∂θ²
u(r,θ) = R(r)Θ(θ) と仮定して変数分離すると、Θ(θ)に対する常微分方程式とR(r)に対するオイラー型方程式が得られます。
フーリエ級数の形で展開し、境界条件を用いて係数を決定することで解が求まるでしょう。
固有値問題との関係
変数分離法は、ラプラシアンの固有値問題と密接に関係しています。
固有値λと固有関数φを使うと、ポアソン方程式の解はそれらの線形結合で表されます。
固有値問題を解くことがポアソン方程式の求解の鍵となる場合が多く、スペクトル解法の基礎にもなっています。
グリーン関数法と有限差分法による数値解法
続いては、より汎用性の高い解法であるグリーン関数法と有限差分法について解説していきます。
グリーン関数法の原理と手順
グリーン関数G(x; y)は、点yに置かれた点源（デルタ関数 δ(x-y)）に対するポアソン方程式の解です。
ΔG(x;y) = δ(x-y)
グリーン関数が求まれば、任意のソース項f(x)に対する解は次の積分で表せます。
u(x) = ∫ G(x;y)f(y)dy
自由空間（無限領域）における3次元のグリーン関数は G(x;y) = -1/(4π|x-y|) で与えられます。
これはクーロンポテンシャルの形と一致しており、電磁気学との深い関係を示しています。
半空間や球などの特殊な領域では、鏡像法（Method of Images）を使ってグリーン関数を構成することができます。
有限差分法の基本手順
複雑な形状の領域や一般的なソース項に対しては、数値解法である有限差分法が有効です。
有限差分法では、連続な空間を格子点に離散化し、微分を差分商で近似します。
∂²u/∂x² ≈ (u_{i+1,j} – 2u_{i,j} + u_{i-1,j}) / h²
同様にy方向も差分化し、全格子点でポアソン方程式の離散版を立てる
これにより、ポアソン方程式は連立一次方程式の形に変換されます。
この連立方程式を解くことで、各格子点での近似解が得られます。
格子幅hを小さくするほど精度が上がりますが、計算コストも増大するため、問題の精度要求と計算資源のバランスが重要でしょう。
数値解法の精度と安定性
有限差分法の精度は格子幅hに依存し、標準的な中心差分を用いると2次精度（誤差がO(h²)のオーダー）が得られます。
大規模な連立方程式を効率よく解くために、共役勾配法・マルチグリッド法・FFT（高速フーリエ変換）などのアルゴリズムが利用されます。
特に周期境界条件の場合、FFTを使ったスペクトル法が非常に効率的です。
まとめ
本記事では、ポアソン方程式の解き方として変数分離法・グリーン関数法・有限差分法を中心に解説しました。
解析的手法は特殊な領域や単純なソース項のときに厳密解を与え、数値的手法は複雑な問題に広く対応できます。
問題の条件に応じて適切な解法を選択することが、ポアソン方程式を使いこなす鍵です。
境界条件の設定・特解と斉次解の分離・解の一意性の確認という基本の流れを押さえたうえで、各解法の手順を習得していきましょう。
