[論文レビュー] Efficient Scaling and Moving Techniques for Spectral Methods in Unbounded Domains

Mingtao Xia, Sihong Shao+1 more|arXiv (Cornell University)|Sep 28, 2020

Ocean Waves and Remote Sensing参考文献 23被引用数 6

ひとこと要約

本稿では、周波数および外部誤差指標を用いて、非有界領域におけるスペクトル法の動的スケーリングおよび移動手法を提案する。これらの手法により、拡散および移動を伴う時間依存問題において高速なスペクトル収束が達成され、細胞増殖モデルにおける爆発の追跡に成功し、固定スケーリング手法を上回る性能を示した。

ABSTRACT

When using Laguerre and Hermite spectral methods to numerically solve PDEs in unbounded domains, the number of collocation points assigned inside the region of interest is often insufficient, particularly when the region is expanded or translated to safely capture the unknown solution. Simply increasing the number of collocation points cannot ensure a fast convergence to spectral accuracy. In this paper, we propose a scaling technique and a moving technique to adaptively cluster enough collocation points in a region of interest in order to achieve a fast spectral convergence. Our scaling algorithm employs an indicator in the frequency domain that is used to determine when scaling is needed and informs the tuning of a scaling factor to redistribute collocation points to adapt to the diffusive behavior of the solution. Our moving technique adopts an exterior-error indicator and moves the collocation points to capture the translation. Both frequency and exterior-error indicators are defined using only the numerical solutions. We apply our methods to a number of different models, including diffusive and moving Fermi-Dirac distributions and nonlinear Dirac solitary waves, and demonstrate recovery of spectral convergence for time-dependent simulations. Performance comparison in solving a linear parabolic problem shows that our frequency scaling algorithm outperforms the existing scaling approaches. We also show our frequency scaling technique is able to track the blowup of average cell sizes in a model for cell proliferation.

研究の動機と目的

• スペクトル法を用いて非有界領域の偏微分方程式を解く際、注目領域におけるコロケーション点の不足という課題に対処すること。
• 解の挙動の変化に適応できない固定スケーリングおよび移動戦略の限界を克服すること。
• 解の導関数や解析的知識に依存せずに、数値解のみに基づく指標を提案し、事前知識なしにコロケーション点を自動的に再配置できること。
• スケーリングと移動の手法を組み合わせることで、拡散および移動を含む時間依存問題においてスペクトル収束を確保すること。
• 非線形および爆発を示す問題、特に細胞増殖モデルに対して、提案手法の有効性を実証すること。

提案手法

• 中間的な数値解のフーリエ変換から導出される周波数指標を導入し、高周波成分を検出し、近似誤差を評価する。
• 周波数指標を用いて時間依存的な座標系のスケーリングを制御し、誤差を低く保ちながらスペクトル収束を維持するため、スケーリング係数β(t)を動的に調整する。
• 外領域における誤差を推定するための外部誤差指標を提案し、移動する解に追随できるようにコロケーション点の再配置を適応的に行う。
• 「まず移動してからスケーリングする」戦略を採用し、翻訳と拡散を併存する問題に対応する。
• 空間離散化に一般化ラゲルおよびエルミート関数を用いた擬スペクトルフレームワークにこれらの手法を適用する。
• 指標の計算に解析的解やその導関数の知識を一切使用せず、数値解のみに依存する。

実験結果

リサーチクエスチョン

• RQ1 非有界領域におけるスペクトル法を、拡散や移動といった解の特徴の変化に適応可能にするにはどうすればよいか？
• RQ2 数値解から導出される周波数ドメイン指標を用いて、動的スケーリングを制御し、スペクトル収束を向上させることは可能か？
• RQ3 非有界領域における外領域の誤差をどのように推定し、それをもとにメッシュの移動を誘発するか？
• RQ4 適応的スケーリングおよび移動手法は、時間依存偏微分方程式において固定スケーリング手法をどの程度上回るか？
• RQ5 提案手法は、細胞増殖モデルのような非有界領域における爆発挙動を正確に追跡できるか？

主な発見

• 周波数依存スケーリングアルゴリズムは、線形放物型問題を解く際、固定スケーリング手法よりも高速なスペクトル収束を達成し、N=20のときt=10で誤差が1.143e-02から8.662e-06未満に低下した。
• 時間発展中に周波数指標が約10^-6を維持するように、スケーリング係数βxを0.9から0.2766に動的に調整することで、高い精度が確保された。
• 細胞増殖モデルにおける解の平均サイズは、t=10で解析解⟨x(t)⟩=5+tと15.001の誤差内で正確に追跡された。
• 構造的細胞集団モデルにおける爆発挙動を効果的に捉えた。平均サイズは時間とともに線形に増加し、適応的スケーリングが解の解像度を維持するために不可欠であった。
• 移動とスケーリングの組み合わせにより、翻訳と拡散を併存する問題を安定して解像でき、解のダイナミクスに関する事前知識なしにスペクトル収束を維持した。
• 外部誤差指標は、解の移動が必要な場合に適切に検出され、メッシュの移動を誘発することができ、コロケーション点が注目領域に適応的に移動した。