[論文レビュー] Tensor calculus in spherical coordinates using Jacobi polynomials. Part-I: Mathematical analysis and derivations

Geoff Vasil, Daniel Lecoanet+3 more|arXiv (Cornell University)|Apr 27, 2018

Geophysics and Gravity Measurements参考文献 33被引用数 8

ひとこと要約

本稿では、スピン重み付き球面調和関数を用いた角方向の記述とスケーリングされたジャコビ多項式を用いた径方向変動に基づく、球座標系における偏微分方程式を解くスペクトル法を提示する。テンソル場を正則性空間に分解し、直交変換を活用することで、極および中心における特異的挙動を自然に取り入れた最小帯域幅のスパースな微分作用素を実現し、単位球内における複雑なPDEの高精度かつ効率的な解法を可能にする。

ABSTRACT

This paper presents a method for the accurate and efficient computations on scalar, vector and tensor fields in three-dimensional spherical polar coordinates. The methods uses spin-weighted spherical harmonics in the angular directions and rescaled Jacobi polynomials in the radial direction. For the 2-sphere, spin-weighted harmonics allow for automating calculations in a fashion as similar to Fourier series as possible. Derivative operators act as wavenumber multiplication on a set of spectral coefficients. After transforming the angular directions, a set of orthogonal tensor rotations put the radially dependent spectral coefficients into individual spaces each obeying a particular regularity condition at the origin. These regularity spaces have remarkably simple properties under standard vector-calculus operations, such as extit{grad} and extit{div}. We use a hierarchy of rescaled Jacobi polynomials for a basis on these regularity spaces. It is possible to select the Jacobi-polynomial parameters such that all relevant operators act in a minimally banded way. Altogether, the geometric structure allows for the accurate and efficient solution of general partial differential equations in the unit ball.

研究の動機と目的

• 3次元球座標系におけるテンソリアルPDEを数値的に安定かつ効率的に解くスペクトル法の開発。
• 高解像度シミュレーションにおいて北極および南極、および球の中心に生じる座標特異性に起因する数値的不安定性の解消。
• 特異点における正則性条件を内蔵した基底関数の構築により、滑らかな解に対してスペクトル収束を保証すること。
• 勾配や発散などの微分作用素（例：grad, div）を、最小帯域幅のスパース行列として導出し、正則性空間全体にわたって効率的に作用させる。
• 地球惑星および天体物理学的応用に適した一般PDEソルバーの構築を可能とし、問題タイプ間でのコード変更を最小限に抑える。

提案手法

• 角方向の依存性を表現するためにスピン重み付き球面調和関数を用い、波数乗算によるフーリエに類似したスペクトル演算を可能にする。
• テンソル場を異なる正則性空間に分解するための直交変換を適用し、各空間がr=0における特定の径方向挙動に対応する。
• 各正則性空間において、各テンソル成分の径方向正則性に一致するように調整されたパラメータを有するスケーリングされたジャコビ多項式を基底関数として採用する。
• ジャコビ多項式の代数的性質と再帰関係を活用して、微分作用素（例：grad, div）をスパース行列として導出する。
• 切断処理にτ法を用い、異なる残差選択が固有値問題における誤差の大きさに顕著な差をもたらすことを示す。
• すべてのテンソル成分にわたってスパarsityを保ち、計算効率を損なわずに最小帯域幅の方法で作用する作用素を構築する。

実験結果

リサーチクエスチョン

• RQ1 球座標系におけるテンソル場は、原点および極における境界条件を自然に満たす正則性空間にどのように分解可能か？
• RQ2 径方向および角方向に、スパースで最小帯域幅の微分作用素を構築するために必要な基底関数と変換は何か？
• RQ3 ジャコビ多項式の代数的性質は、テンソル場の径方向変動の効率的かつ高精度なスペクトル表現をどのように可能にするか？
• RQ4 τ法による異なる切断戦略が、この枠組みにおける固有値問題の精度に与える影響は何か？
• RQ5 この手法は、地球惑星および天体物理学的流体動力学に現れる一般PDEに対しても、高精度かつ高効率に適用可能か？

主な発見

• 適切なジャコビ多項式パラメータの選択により、最小帯域幅の微分作用素が得られ、計算効率が保証される。
• テンソル成分から正則性空間への直交変換は、球面調和関数の次数ℓにのみ依存し、単純な再帰式により計算可能であり、高ランクテンソルの取り扱いが可能になる。
• スピン重み付き球面調和関数により、角方向の演算をフーリエ級数と同様に実行でき、微分作用素は波数乗算として作用する。
• 基底関数は座標特異点（極および中心）における正則性条件を自動的に満たし、高解像度シミュレーションにおける数値的不安定性を排除する。
• PDEの種別（例：流体力学から磁気流体力学）を簡単に切り替えることが可能で、コード変更を最小限に抑えられるため、広範な適用可能性を示す。
• τ法による異なる切断戦略は、誤差の大きさに数個のオーダーの差をもたらすことが明らかとなり、スペクトル法における残差選択の重要性が強調される。