[論文レビュー] A new mixed potential representation for the equations of unsteady, incompressible flow

Leslie Greengard, Shidong Jiang|arXiv (Cornell University)|Sep 22, 2018

Numerical methods in engineering参考文献 34被引用数 5

ひとこと要約

本稿では、熱ポテンシャルと調和ポテンシャルの組み合わせを用いて、非定常かつ非圧縮性のストークスおよびナビエ＝ストークス流れの新しい混合ポテンシャル定式化を提案する。これにより、境界積分方程式の結合系が得られ、その系は結合源積分方程式（CSIE）と呼ばれる良好に条件付けられた系となる。この手法により、スペクトル遅延補正を用いた予測子・修正子スキームを介して高次精度が達成され、大きな時間刻みでも、高周波数の圧力場に対しても安定した性能を示す。

ABSTRACT

We present a new integral representation for the unsteady, incompressible Stokes or Navier-Stokes equations, based on a linear combination of heat and harmonic potentials. For velocity boundary conditions, this leads to a coupled system of integral equations: one for the normal component of velocity and one for the tangential components. Each individual equation is well-condtioned, and we show that using them in predictor-corrector fashion, combined with spectral deferred correction, leads to high-order accuracy solvers. The fundamental unknowns in the mixed potential representation are densities supported on the boundary of the domain. We refer to one as the vortex source, the other as the pressure source and the coupled system as the combined source integral equation.

研究の動機と目的

• 従来の投影法や渦度に基づく手法が抱える数値的困難を回避する、非定常非圧縮性ストークス方程式およびナビエ＝ストークス方程式の新しい数学的表現の開発。
• 境界条件の下で拡散する速度場を進める際の非圧縮性の強制の難しさ、特に圧力勾配が存在する場合の取り扱いの改善。
• 調和ポテンシャルと熱ポテンシャルに基づく境界積分定式化を構築し、圧力と渦度源の物理的役割を明確に分離する。
• 良好に条件付けられた個々の式を保証する予測子・修正子フレームワーク上で、スペクトル遅延補正（SDC）を用いて高次時間精度を達成する。
• 高速多重極法やその他の高速アルゴリズムと互換性を持つ、効率的で高次数の数値ソルバーの基盤を提供する。

提案手法

• 粘性拡散のための熱ポテンシャルと、圧力および非圧縮性のための調和ポテンシャルを組み合わせた混合ポテンシャル表現を用いて、非定常ストークス方程式を定式化する。
• 境界上に定義された、法線速度成分（渦源）と接線成分（圧力源）の両方を含む、結合した境界積分方程式系を導出する。
• 境界上の密度（渦源および圧力源を表す）を未知数とする結合源積分方程式（CSIE）を導入する。
• CSIEに予測子・修正子戦略を適用し、2次時間離散化を予測子として用い、スペクトル遅延補正（SDC）により高次精度を達成する。
• 時間方向の積分を積分積分法（product integration）を用いて実装し、ガウスカーネルを含む主要な積分を解析的に評価することで、時間方向に4次精度を達成する。
• 結合系内の個々の式が良好に条件付けられていることを保証し、大きな時間刻みでも安定かつ正確な解が得られる。

実験結果

リサーチクエスチョン

• RQ1 熱ポテンシャルと調和ポテンシャルに基づく混合ポテンシャル定式化は、非定常非圧縮性流れに対して良好に条件付けられた高次境界積分方程式系を提供できるか？
• RQ2 分数ステップ法や圧力ポisson方程式の解法を用いずに、境界積分表現を通じて非圧縮性制約を自然に強制できるか？
• RQ3 スペクトル遅延補正を用いた予測子・修正子スキームは、大きな時間刻みに対しても安定性と精度を保ちながら高次時間精度を達成できるか？
• RQ4 境界表現における渦源と圧力源の分離が、特に周期的圧力場において解の条件数と精度に与える影響は何か？
• RQ5 この新しい定式化は、非定常形状や複雑な境界条件（多重連結領域を含む）へと拡張可能か？

主な発見

• 提案された混合ポテンシャル表現により、各成分（渦源および圧力源）が独立に適切に定式化された良好に条件付けられた境界積分方程式系が得られた。
• 数値実験の結果、SDC⁰₅（1回の補正スイープ）では2次収束が確認され、SDC¹₅（2回の補正スイープ）では著しく高い収束率が得られ、最も細かいグリッドでは観測された精度の次数が約2.6に達した。
• 圧力場の振動を十分に解像しない大きな時間刻みでも、高精度を維持しており、急速な圧力変動に対して強いロバストネスを示した。
• 積分積分法と十分な補正スイープ回数を用いたSDCを適用することで、時間方向に4次精度が達成され、収束表から一貫した誤差低減が確認された。
• 渦源と圧力源の分離により、境界密度を介して直接的に速度境界条件を強制できるため、実装が簡素化され、安定性が向上した。
• 標準的な熱ポテンシャルと調和ポテンシャルに依存しているため、高速多重極法などの高速アルゴリズムと互換性がある。