スペクトル法によるPDE解法

なぜ微分を「掛け算」に置き換えられるのか

有限差分法や有限要素法が解を局所的な多項式でつなぐ「局所近似」なのに対し、スペクトル法は解を領域全体にわたる大域的な基底関数の和で表現します。周期境界条件の問題であれば、解 u(x) を複素指数関数（フーリエ基底）の重ね合わせで書きます。

横にスクロール

フーリエ級数展開（周期 L、波数 k）
  u(x) = Σ_k  û_k * exp( i k x )

  û_k : 波数 k のフーリエ係数（複素数）
  k   : 2π/L の整数倍（離散波数）

ここで x について微分すると、各項は exp(i k x) の形なので、微分は指数関数を保ったまま前に係数 i k を落とすだけになります。

微分のフーリエ空間での作用
  d/dx  exp(i k x) = (i k) * exp(i k x)

  → u'(x)  のフーリエ係数は  (i k)   * û_k
  → u''(x) のフーリエ係数は  (i k)^2 * û_k = -k^2 * û_k

つまり物理空間では微分（差分近似）だった操作が、波数空間では係数ベクトルへの対角な掛け算に化けます。これがスペクトル法の核心で、微分演算子がフーリエ基底に対して対角化されるという事実に尽きます。実際の計算はFFTで係数 {û_k} を得て、波数 i k を掛け、逆FFTで物理空間へ戻すという流れになり、微分1回あたり O(N log N) で済みます。差分のような打ち切り誤差は原理的に入らず、離散化された波数の範囲内では微分が厳密です。

スペクトル収束 ── 少ない点で高精度が出る理由

スペクトル法が有限差分・有限要素と決定的に違うのは収束の速さです。解が滑らか（無限回微分可能で周期的）なとき、フーリエ係数 û_k は k が大きくなるにつれて指数的に速く0へ減衰します。したがって有限個の波数で打ち切ったときの誤差も指数的に小さくなり、これをスペクトル収束（指数収束）と呼びます。

例えば差分法で誤差を4分の1にするには格子点をおよそ倍にする必要がありますが、スペクトル法では数点足すだけで誤差が桁で落ちることがあります。この「解像度あたりの情報密度の高さ」が、同じ精度をはるかに少ないメモリと演算で達成できる理由です。乱流の直接数値シミュレーション（DNS）でスペクトル法が長く標準であり続けてきたのも、限られた格子点で最大限の有効波数域を確保できるからにほかなりません。

非周期問題とチェビシェフ基底

フーリエ基底が使えるのは周期境界条件のときだけです。壁のある流路や有限区間の問題では、両端で周期性が成り立たないため、フーリエ級数をそのまま当てると端点でギブス現象が起きます。そこで非周期の有界区間ではチェビシェフ多項式を基底に用います。

チェビシェフ多項式 T_n(x) は区間 [-1, 1] 上で定義され、変数変換 x = cos(θ) を通すと余弦級数、すなわちフーリエ級数に帰着します。この関係のおかげで、チェビシェフ係数の計算も余弦変換（FFTの一種）で O(N log N) にできます。さらに重要なのは格子点の取り方で、等間隔ではなく端点近くに密集するチェビシェフ・ガウス・ロバット点（両端で密、中央で粗）を使います。

チェビシェフ基底では微分はフーリエほど単純な対角行列にはならず、係数間の漸化式（または微分行列との積）で表されますが、波数空間的な扱いで微分を高精度化できる点は共通です。周期方向はフーリエ、壁に垂直な方向はチェビシェフ、というように基底を軸ごとに使い分ける構成が、チャネル乱流などの実務で広く使われています。

擬スペクトル法とエイリアシング（2/3則）

線形項は波数空間で対角な掛け算で済みますが、実務のPDEには非線形項（例えば移流項 u・∂u/∂x）が付きものです。非線形項を波数空間で直接畳み込みとして計算すると O(N^2) になり、スペクトル法の速度優位が消えてしまいます。そこで擬スペクトル法（pseudospectral method）を使います。考え方は単純で、掛け算だけを物理空間で行うというものです。

擬スペクトル法での非線形項 u * v の計算
  1. û_k, v̂_k を逆FFTで物理空間の u(x), v(x) に戻す
  2. 物理空間で点ごとに積 w(x) = u(x) * v(x) を計算（O(N)）
  3. w(x) をFFTで波数空間 ŵ_k に戻す

  畳み込みを O(N^2) で解く代わりに、
  FFT2回＋逆FFT1回（各 O(N log N)）で済む

ただしこの手続きには落とし穴があります。物理空間での積は波数空間では畳み込みに相当し、2つの波数 p と q の積は波数 p+q の成分を生みます。p+q が離散格子で表現できる最大波数を超えると、その高波数成分は折り返して低波数の偽の成分として混入します。これがエイリアシングで、放置すると解に非物理的なエネルギーが蓄積し、時間発展が不安定化することがあります。

対策の定番が2/3則（two-thirds rule）です。2次の非線形項（積が2つの関数の掛け算）であれば、保持する波数を全体の下側3分の2に切り詰め、上側3分の1をゼロで埋めておけば、積で生じる折り返し成分が有効波数域に入り込まなくなります。

2/3則によるデアリアシング（2次非線形の場合）
  総モード数 N のうち
    |k| が N/3 まで        → 保持
    |k| が N/3 を超える帯域 → 0 でパディング（切り捨て）

  → 積 p+q による折り返しが有効域を汚染しなくなる
  → 実効的な解像度は 2/3 に減るが、解は汚染から守られる

代償として実効的な解像度が3分の2に落ちるため、必要な有効波数を確保するには格子点をその分多めに取る必要があります。それでも、エイリアシング除去なしに安定な長時間積分は困難なため、DNSなどでは2/3則、あるいは位相をずらして2回計算し平均するフェーズシフト法が実務の標準対策になっています。

実務での位置づけ

スペクトル法は「滑らかな解・単純形状・高精度要求」という条件が揃う場面で圧倒的に強く、乱流DNS、地球流体、気象・気候モデルの一部などで使われ続けています。一方で複雑形状や不連続解には向かず、そこでは有限差分法と有限要素法や、両者の中間として要素内をスペクトル的に高次化するスペクトル要素法が選ばれます。並列化の観点では、時間発展の各ステップで複数回の多次元FFTを要するため、実行時間は領域分割法（並列PDE計算）で扱う分割戦略とFFT転置の通信コストに強く支配されます。

演算精度そのものも見落とせません。スペクトル法は微分誤差が極めて小さいぶん、丸め誤差やFFTの累積誤差が相対的に無視できなくなる領域まで踏み込むため、実装では数値安定性と浮動小数点誤差（HPC文脈）の考慮が精度の頭打ちを左右します。

まとめ

• スペクトル法は解を大域的な基底で展開し、微分を波数空間での対角な掛け算に置き換える手法で、微分を厳密かつ O(N log N) で計算できる。
• 解が滑らかならスペクトル収束（指数収束）が働き、有限差分・有限要素の代数的収束よりはるかに少ない格子点で高精度を達成できるが、不連続があればギブス現象で崩れる。
• 周期境界にはフーリエ基底、非周期・有界区間には端点で密なチェビシェフ基底を使い、軸ごとに基底を使い分ける構成が実務で一般的である。
• 非線形項は物理空間で積を取る擬スペクトル法が実用的だが、積が生むエイリアシングを2/3則などで除去しないと偽エネルギーが溜まり不安定化するため、デアリアシングは必須である。