[論文レビュー] A structure-preserving algorithm for the fractional nonlinear Schrödinger equation based on the SAV approach

Yayun Fu, Wenjun Cai+1 more|arXiv (Cornell University)|Nov 18, 2019

Fractional Differential Equations Solutions参考文献 14被引用数 7

ひとこと要約

本稿では、空間分数非線形シュレーディンガー方程式に対する線形陰的エネルギー保存スキームを、スカラー補助変数（SAV）アプローチと空間方向のフーリエ擬スペクトル法、時間方向のクランク・ニコルソン時間積分法を組み合わせて提案する。この手法は、各時間ステップで結合されていない定数係数の連立一次方程式を解くことで、完全陰的スキームに比べて計算効率を著しく向上させつつ、離散的意味での正確なエネルギー保存を達成する。

ABSTRACT

The main objective of this paper is to present an efficient structure-preserving scheme, which is based on the idea of the scalar auxiliary variable approach, for solving the space fractional nonlinear Schrödinger equation. First, we reformulate the equation as a Hamiltonian system, and obtain a new equivalent system via introducing a scalar variable. Then, we construct a semi-discrete energy-preserving scheme by using the Fourier pseudo-spectral method to discretize the equivalent system in space direction. After that, applying the Crank-Nicolson method on the temporal direction gives a linear implicit scheme in the fully-discrete version. As expected, the proposed scheme can preserve the energy exactly and more efficient in the sense that only decoupled equations with constant coefficients need to be solved at each time step. Finally, numerical experiments are provided to demonstrate the effectiveness and conservation of the scheme.

研究の動機と目的

• 空間分数非線形シュレーディンガー方程式に対して、エネルギーを正確に保存する構造保存型数値スキームを開発すること。
• 従来の完全陰的スキームの高い計算コストを克服するため、線形陰的手法を構築すること。
• スカラー補助変数（SAV）アプローチを、これまで適用例のなかった分数階微分方程式に拡張すること。
• SAV法とフーリエ擬スペクトル法を組み合わせ、高い空間精度と計算効率を実現すること。
• 広範な数値実験を通じて、スキームの長時間安定性および保存性を実証すること。

提案手法

• 非線形項を分離するためにスカラー補助変数を導入し、分数NLS方程式をハミルトニアン系に再定式化する。
• 空間方向にフーリエ擬スペクトル法を用いて半離散化し、FFTを活用することでスペクトル精度と効率的な計算を達成する。
• 時間方向にクランク・ニコルソン法を適用して、非線形項の反復的解法を回避する線形陰的完全離散スキームを構築する。
• SAVアプローチによりエネルギー関数を二次形式に再定式化することで、構成上正確なエネルギー保存を保証する。
• 定数係数を持つ線形系を導出し、フーリエ空間における対角化により計算効率を向上させる。
• 周期的境界条件を実装し、さまざまなポテンシャル関数および分数量子数に対する1次元および2次元問題でスキームを検証する。

実験結果

リサーチクエスチョン

• RQ1 SAVアプローチは、分数階微分方程式に対しても、エネルギー保存スキームを構築するために成功裏に拡張可能か？
• RQ2 提案された線形陰的スキームは、分数NLS方程式の離散的エネルギーを正確に保存するか？
• RQ3 2次元ケースにおいて、分数量子数αがソリトン解の形状および減衰率に与える影響は何か？
• RQ4 異なる空間的・時間的離散化条件下での、スキームの収束速度および長期的保存性能は何か？
• RQ5 調和ポテンシャルや光格子ポテンシャルなどの異なるポテンシャル関数を用いた場合、スキームの性能はどのように変化するか？

主な発見

• 提案されたFSAVスキームは、離散的エネルギーを正確に保存し、十分に細かいグリッドでは相対的エネルギー誤差が1e-13未満に抑えられ、理論的エネルギー保存の妥当性が確認された。
• 空間方向にスペクトル収束を示し、1次元問題ではN=64で観測収束次数約10.56、N=128で最大20.30に達し、高い空間精度が裏付けられた。
• 時間方向収束次数は約2.0であり、クランク・ニコルソン法と整合的であることが、τ=0.1からτ=0.001に減少したL2誤差の低下から明らかになった。
• 分数量子数αが小さくなるほど、波動関数の減衰が速くなり、勾配が急になる傾向を示し、α→2の極限で古典的NLS方程式に近づくことが観察された。
• 質量保存は機械精度まで維持され、全テスト設定において相対的質量誤差は1e-14未満に保たれた。
• 計算効率が高く、各時間ステップで結合されていない定数係数の線形系を解くのみで、反復的非線形ソルバーを回避している。