[論文レビュー] A Modulus-Squared Dirichlet Boundary Condition for Time-Dependent Complex Partial Differential Equations and its Application to the Nonlinear Schrödinger Equation

Ronald M. Caplan, R. Carretero-González|arXiv (Cornell University)|Oct 4, 2011

Numerical methods for differential equations参考文献 12被引用数 6

ひとこと要約

本稿では、時間依存複素PDE、特に非線形シュレーディンガー方程式（NLSE）に対して、単純で高次精度に適合する、モジュラス二乗ディリクレ（MSD）境界条件を導入する。この境界条件は、内部点における解の時間微分を用いて境界における時間微分を計算することで、境界でモジュラス二乗値を一定に保つ。これにより、渦やダークソリトンのようなコherent構造の安定的かつ高精度なシミュレーションが、従来のディリクレ条件やゼロラプラシアン条件に比べて小さな計算領域で可能になる。

ABSTRACT

An easy to implement modulus-squared Dirichlet (MSD) boundary condition is formulated for numerical simulations of time-dependent complex partial differential equations in multidimensional settings. The MSD boundary condition approximates a constant modulus-square value of the solution at the boundaries. Application of the MSD boundary condition to the nonlinear Schrödinger equation is shown, and numerical simulations are performed to demonstrate its usefulness and advantages over other simple boundary conditions.

研究の動機と目的

• 時間依存複素PDEに対して、位相の振動により失敗する標準ディリクレ条件の実装における課題に対処すること。
• ボーズ＝アインシュタイン凝縮（BEC）や非線形光学のようないずれの系でも一定密度背景を持つ系に適した、境界で解のモジュラス二乗値を一定に保つ境界条件を開発すること。
• 高次精度数値スキームの正確性を保ちつつ、容易に実装可能な単純な境界条件を提供すること。
• MSD境界条件がNLSEにおけるコherent構造（例えばヴォルテックスやダークソリトン）のシミュレーションにおいて有効であることを示すこと。
• MSD条件が複雑な人工境界処理を要せず、不自然な反射を避けつつ安定性を維持することを示すこと。

提案手法

• MSD境界条件は、境界における時間微分を隣接する内部点における解とその時間微分の関数として表現することにより定式化される：$ \left.\frac{d\Psi}{dt}\right|_b \approx i\,\text{Im}\left[\frac{1}{\Psi_{b-1}}\left.\frac{d\Psi}{dt}\right|_{b-1}\right]\Psi_b $。
• この手法により、内部解の位相の進化に基づいて動的に時間微分を調整することで、境界における解のモジュラス二乗値が一定に保たれる。
• 多次元空間において、4次精度ルンゲ＝クッタ有限差分法を用いてNLSEにこの手法を適用する。
• 境界条件は、元の数値スキームの形式的精度を保つように実装され、モジュラス二乗値が境界で一定であるという仮定が成り立つ限り有効である。
• 線形化を用いた安定性解析により、時間ステップの上限が導かれる。この上限は境界値と解の周波数に依存し、実用的な時間ステップは理論的上限の80%に設定される。
• MATLABによる実装例が提供されており、時間ステッピングスキームにおいて境界ごとにコードを2行のみで適用可能であることが示されている。

実験結果

リサーチクエスチョン

• RQ1 時間依存複素PDEに対して、境界で解のモジュラス二乗値を一定に保つ単純な境界条件を考案可能か？
• RQ2 MSD境界条件は、NLSEにおけるヴォルテックスやダークソリトンのようなコherent構造のダイナミクスをどの程度正確に保持できるか？
• RQ3 MSD条件は、人工的な反射や不安定性を導入せずに、高次精度数値スキームの正確性を維持できるか？
• RQ4 標準ディリクレ条件やゼロラプラシアン条件に比べ、MSD条件は必要なグリッドサイズを小さくできるか？
• RQ5 MSD条件の安定性挙動は何か？理論的上限と比較してどうか？

主な発見

• MSD境界条件により、境界でモジュラス二乗値が一定である解（例えばダークソリトンやヴォルテックス）のシミュレーションが、標準ディリクレ条件やゼロラプラシアン条件に比べて顕著に小さな計算領域で正確に可能になる。
• 解のモジュラス二乗値が境界で概ね一定である限り、MSD手法は元の数値スキームの形式的精度を保つ。
• 数値シミュレーションでは、MSD条件が不自然な反射を効果的に抑制し、内部ダイナミクスに干渉しないことが示された。これは、渦度や共動する背景流れが存在する場合でも同様に成り立つ。
• 線形化から導かれた安定性限界は実用的であり、理論的上限の80%に時間ステップを設定することで、大多数のケースで安定性が保証された。
• MSD条件は、境界値が一定の問題では標準ディリクレ条件と同等に有効であるが、実装がはるかに簡単で、非ゼロの背景密度を持つ系に適している。
• この手法は頑健でスケーラブルであり、コードの断片が最小限に抑えられる。時間ステッピングスキームに境界条件を適用するには、境界ごとにたった2行のコードで十分である。