[論文レビュー] A spectral method in space and time to solve the advection-diffusion and wave equations in a bounded domain

J. Pétri|arXiv (Cornell University)|Jan 7, 2014

Differential Equations and Numerical Methods参考文献 2被引用数 3

ひとこと要約

本稿では、有界領域における線形移流拡散方程式および波動方程式を解くために、空間および時間の両方でチェビシェフ多項式を用いた完全スペクトル法を提示する。この手法は、空間と時間の対称的取り扱いによりスペクトル精度を達成するが、行列逆行列計算の制限により、特に3次元または高解像度の場合のスケーラビリティに制限がある。

ABSTRACT

The advection-diffusion and wave equations are the fundamental equations governing any physical law and therefore arise in many areas of physics and astrophysics. For complex problems and geometries, only numerical simulations can give insight into quantitative and accurate behavior of the seeked solutions. The standard numerical algorithm to solve partial differential equations is to split the space and time discretisation separately into different uncorrelated methods. Time is usually advanced by explicit schemes, or, for too restrictive time steps, by implicit or semi-implicit algorithms. This separate time and space slicing is artificial and sometimes unpractical. Indeed, treating space and time directions symmetrically and simultaneously without splitting is highly recommended in some problems like diffusion. It is the purpose of this work to present a simple numerical algorithm to solve the standard linear scalar advection-diffusion and wave equations using a fully spectral method in a two-dimensional Cartesian $(x,t)$ bounded space-time domain. Generalization in three-dimensions $(x,y,t)$ is shown for the pure diffusion problem. The basic idea is to expand the unknown function in Chebyshev polynomials for the spatial variables $(x,y)$ as well as for the time variable $t$. We show typical examples and demonstrate the spectral accuracy of the method. The great advantage of fully spectral methods resides in their high-accuracy for a relatively small number of grid points (for sufficiently smooth solutions) compared to standard time-stepping techniques.

研究の動機と目的

• 有界領域における線形偏微分方程式を解くために、空間と時間の両方を対称的に取り扱う完全スペクトル数値法の開発。
• チェビシェフ展開を用いた空間変数および時間変数の両方におけるスペクトル法の実現可能性と精度の検証。
• 空間と時間の離散化に分離演算子を用いないことで、従来の時間積分スキームの限界を克服すること。
• さまざまな境界条件の下で、標準的なテスト問題（移流、拡散、波動方程式）における本手法の性能評価。
• 特に行列逆行列計算に伴う高メモリおよび高計算コストという計算上のボトルネックを特定し、反復解法による未来の改善策を提案すること。

提案手法

• 2次元または3次元の直角座標系領域において、空間変数（x, y）および時間変数（t）の両方をチェビシェフ多項式で展開する。
• チェビシェフ＝ガウス＝ロバットノードにおけるコロケーション点を用いて、偏微分方程式をスペクトル係数に関する代数方程式系に変換する。
• ガラーキン法またはコロケーション法にタウ法またはペナルティ法を適用して、境界条件および初期条件を強制する。
• 空間および時間の微分をスペクトル空間で表すための微分行列（乗算行列）を構築する。
• 得られた大規模でスパースな線形方程式系を直接行列逆行列計算により解き、解係数を求める。
• スペクトル基底および行列演算を（x, y, t）領域に拡張することで、2次元空間への手法の一般化を行う。

実験結果

リサーチクエスチョン

• RQ1 空間および時間の両方で完全スペクトル法を用いることで、有界領域における線形移流拡散方程式および波動方程式に対して高い精度が達成可能か？
• RQ2 スペクトル法における空間と時間の対称的取り扱いは、従来の時間積分アプローチと比較して、精度および安定性においてどのように異なるか？
• RQ3 多次元問題におけるスペクトル法において、直接行列逆行列計算の計算上の制限は何か？
• RQ4 2次元および3次元設定において、滑らかな解に対して、本手法がどの程度スペクトル収束性を保っているか？
• RQ5 非一様な拡散係数および出射波や反射波を含む複雑な境界条件に対しても、本手法を拡張可能か？

主な発見

• 本手法はスペクトル収束性を達成しており、(x,t)空間で129×129グリッドポイントにおいて波動方程式の相対誤差が約10⁻¹¹まで低下した。
• 2次元拡散方程式において、17×17×17解像度で誤差はたったの1.065×10⁻¹³にまで低下し、スペクトル精度が確認された。
• 行列逆行列計算に伴うメモリおよび計算コストが解像度の向上を著しく制限しており、3次元では33×33×33が最大で実行可能な解像度であった。
• 本手法は、わずかなグリッドサイズでも高い精度を維持しており、標準的な有限差分スキームと比較してスペクトル法の優位性を示している。
• 空間と時間の対称的取り扱いにより、境界条件の一貫した適用が可能であり、時間と空間の演算子を人工的に分離する必要がなくなった。
• アルゴリズムは非一様な拡散係数に対しても容易に拡張可能であり、さらなる開発により円柱座標系や球座標系への応用も可能となる。