[論文レビュー] How to overcome the Courant-Friedrichs-Lewy condition of explicit discretizations?

Denys Dutykh|arXiv (Cornell University)|Nov 24, 2016

Aquatic and Environmental Studies参考文献 6被引用数 7

ひとこと要約

本稿では、熱方程式の陽的差分離散化における厳格なコールラント＝フレイドリッヒス＝ルービン（CFL）条件を克服するためのハイブリッド陽的有限差分スキームを提案する。これにより、完全な陰的解法を必要とせずに、時間ステップを大きくとることができる。ハイブリッド化の手法として、制御された双曲的性質の導入、または複数時間レベルの更新（例：Dufort–Frankel、Saulyev、双曲化）を用いる。これらの手法により、無条件安定性または緩和されたCFL制限が達成され、陽的スキームの実装のしやすさや並列化の利便性を保ちながら、古典的な2次精度スキームに比べて精度が若干低下する代わりに、長時間シミュレーションに適した安定性が得られる。

ABSTRACT

This manuscript contains some thoughts on the discretization of the classical heat equation. Namely, we discuss the advantages and disadvantages of explicit and implicit schemes. Then, we show how to overcome some disadvantages while preserving some advantages. However, since there is no free lunch, there is a price to pay for any improvement in the numerical scheme. This price will be thoroughly discussed below. In particular, we like explicit discretizations for the ease of their implementation even for nonlinear problems. Unfortunately, when these schemes are applied to parabolic equations, severe stability limits appear for the time step magnitude making the explicit simulations prohibitively expensive. Implicit schemes remove the stability limit, but each time step requires now the solution of linear (at best) or even nonlinear systems of equations. However, there exists a number of tricks to overcome (or at least to relax) severe stability limitations of explicit schemes without going into the trouble of fully implicit ones. The purpose of this manuscript is just to inform the readers about these alternative techniques to extend the stability limits. It was not written for classical scientific publication purposes.

研究の動機と目的

• 長時間シミュレーションにおいて使用が制限される、熱方程式に対する陽的有限差分スキームにおける極めて厳しい時間ステップ制限を解消すること。
• 完全に陰的スキームが各時間ステップで大きな線形または非線形方程式系を解く必要があるのを避ける代替数値的手法を検討すること。
• 時間ステップを大きくとれるように離散演算子を変更することで、数値的安定性と精度のトレードオフを調査すること。
• 月や年単位の時間発展を必要とする、建築物理学やその他の長時間シミュレーション応用分野に実用的かつ実装可能な手法を提供すること。

提案手法

• Dufort–Frankel法を適用し、3層時間ステンシルを用い、現在時刻と直前の時刻における空間微分の平均をとることで、無条件安定性を達成する。
• Saulyev法を実装し、時間微分の近似を変更することで、時間ステップを拡大可能にしつつ、特定の離散化戦略によって安定性を維持する陽的スキームを採用する。
• 双曲化技術を用いて、放物型の熱方程式を双曲型系に再定式化し、有限速度での情報伝播を導入することで、CFL条件をΔt = O(Δx)に緩和する。
• 分散関係解析と誤差推定を実施し、とくに双曲化法の安定性と精度を評価する。
• 修正された方程式を導出し、導入された双曲的性質が、元の放物型演算子に比べて制御された精度の損失を引き起こすことを示す。
• 改善されたスキームを非線形問題およびソース項を含む多次元問題に一般化し、陽的構造と安定性を維持する。

実験結果

リサーチクエスチョン

• RQ1 熱方程式に対する陽的有限差分スキームを、各時間ステップでの非線形方程式系の解法を必要とせずに、無条件安定化できるか？
• RQ2 標準的な陽的スキームでは放物型性のため情報伝播が無限速となるが、その有限速度伝播を、陽的離散化に回復する方法は何か？
• RQ3 CFL条件を緩和するために離散演算子を変更する際の、安定性と精度のトレードオフはどのように評価できるか？
• RQ4 陽的スキームと複数時間レベル更新を組み合わせたハイブリッドスキームは、長時間シミュレーションにおいて安定性と計算効率の両立を達成できるか？
• RQ5 双曲化法は非線形および多次元問題にどの程度一般化可能であり、安定性と精度を維持できるか？

主な発見

• Dufort–Frankel法およびSaulyev法は無条件安定であり、時間ステップを任意に大きくとれるが、陽的かつ計算的に効率的である。
• 双曲化法は、放物型方程式を双曲型系に再定式化することで、有限速度での情報伝播を導入し、CFL条件を古典的なΔt = O(Δx²)からΔt = O(Δx)に緩和する。
• すべての改良スキームは、制御された双曲的性質または分散性を導入し、古典的な2次精度放物型スキームに比べて精度が低下するが、適切な条件下でO(Δx²)収束が確認される。
• 修正されたスキームは、陽的スキームの実装のしやすさや並列化の利便性を維持しており、建築物理学などの大規模シミュレーション応用に適している。
• 双曲化法は、圧縮性ナビエ＝ストークス方程式やMHD方程式のような複雑な系に対しても妥当性が確認されており、運動論的理論に基づく広範な適用可能性が示唆される。
• 安定性と精度のトレードオフは定量的に評価可能である。ユーザーは意図的に精度を低下させることで、著しい安定性の向上を達成でき、従来の陽的スキームでは不可能だった長時間シミュレーションが可能になる。