[论文解读] 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.

研究动机与目标

• 开发一种在三维球坐标系中求解张量型PDE的数值稳定且高效的谱方法。
• 解决在高分辨率模拟中由北极和南极坐标奇点以及球心处奇点引起的数值不稳定性问题。
• 构建一种能自然满足奇点处正则性条件的基函数，确保光滑解的谱收敛性。
• 推导出作用于不同正则性空间的稀疏、最小带宽的微分算子（例如梯度、散度），实现高效计算。
• 为地球物理与天体物理应用中的通用PDE求解器提供支持，且在不同问题类型间仅需极少代码修改。

提出的方法

• 使用自旋加权球谐函数表示角向依赖，通过波数乘法实现类似傅里叶的谱运算。
• 应用正交变换将张量场分解到不同的正则性空间，每个空间对应于在 r=0 处特定的径向行为。
• 在每个正则性空间上采用一组参数可调的缩放雅可比多项式作为基函数，使其与每个张量分量的径向正则性相匹配。
• 通过利用雅可比多项式的代数性质及其递推关系，将导数算子（如梯度、散度）表示为稀疏矩阵。
• 采用τ方法进行截断，表明不同的残差选择在特征值问题中会导致截然不同的误差量级。
• 构建在所有张量分量间以最小带宽方式作用的算子，保持稀疏性与计算效率。

实验结果

研究问题

• RQ1 如何将球坐标系中的张量场分解为正则性空间，使其自然满足原点与极点处的边界条件？
• RQ2 何种基函数与变换可实现径向与角向方向上稀疏、最小带宽的微分算子构造？
• RQ3 雅可比多项式的代数性质如何实现张量场中径向变化的高效且精确的谱表示？
• RQ4 通过τ方法实现的不同截断策略对本框架中特征值问题精度有何影响？
• RQ5 该方法能否在球几何时实现通用PDE的高精度与高效率求解，包括地球物理与天体物理流体动力学中的问题？

主要发现

• 通过选择适当的雅可比多项式参数，该方法实现了最小带宽的导数算子，确保了计算效率。
• 从张量分量到正则性空间的正交变换仅依赖于球谐函数的阶数 ℓ，且可通过简单递推计算，从而支持高阶张量处理。
• 自旋加权球谐函数使得角向运算可类比傅里叶级数，导数作用等价于波数乘法。
• 该基函数能自动满足坐标奇点（极点与中心）处的正则性条件，消除了高分辨率模拟中的数值不稳定性。
• 该框架可实现PDE类型间的无缝切换（如流体动力学与磁流体动力学），仅需极少代码修改，展现出广泛适用性。
• 通过τ方法实现的不同截断策略导致误差相差多个数量级，凸显了残差选择在谱方法中的关键作用。