[论文解读] Range-separated tensor representation of the discretized multidimensional Dirac delta and elliptic operator inverse

Boris N. Khoromskij|arXiv (Cornell University)|Dec 6, 2018

Tensor decomposition and applications参考文献 19被引用 4

一句话总结

本文提出了一种针对离散化多维狄拉克δ函数和椭圆算子逆的、与算子相关的范围分离（RS）张量近似方法，利用格林核的低秩张量分解。通过将格林核分解为短程和长程分量，并施加椭圆算子，该方法实现了$O(n)$复杂度量级，并可高效、准确地求解含多个奇异电荷的势方程——这对通过泊松-玻尔兹曼方程建模生物分子静电势至关重要。

ABSTRACT

In this paper, we introduce the operator dependent range-separated tensor approximation of the discretized Dirac delta in $\mathbb{R}^d$. It is constructed by application of the discrete elliptic operator to the range-separated decomposition of the associated Green kernel discretized on the Cartesian grid in $\mathbb{R}^d$. The presented operator dependent local-global splitting of the Dirac delta can be applied for solving the potential equations in non-homogeneous media when the density in the right-hand side is given by the large sum of pointwise singular charges. We show how the idea of the operator dependent RS splitting of the Dirac delta can be extended to the closely related problem on the range separated tensor representation of the elliptic resolvent. The numerical tests confirm the expected localization properties of the obtained operator dependent approximation of the Dirac delta represented on a tensor grid. As an example of application, we consider the regularization scheme for solving the Poisson-Boltzmann equation for modeling the electrostatics in bio-molecules.

研究动机与目标

• 为解决传统基于网格的方法在求解具有奇异狄拉克δ源项的椭圆PDE时出现的数值不稳定性和高计算成本问题。
• 开发一种基于张量的、与算子相关的狄拉克δ函数分解方法，将其划分为短程和长程分量，以改善条件数和精度。
• 将该方法扩展至椭圆算子的逆以及非均匀介质中含多个点电荷的势方程求解。
• 利用低秩张量格式，实现基于泊松-玻尔兹曼方程的大规模生物分子系统可扩展模拟。
• 证明集体狄拉克δ函数的长程部分随粒子数呈对数增长，从而在大规模N粒子系统中实现线性-对数复杂度。

提出的方法

• 通过范围分离（RS）张量格式，将椭圆算子的离散格林核分解为短程和长程分量。
• 对每个分量应用椭圆算子$\mathcal{L}$，以生成狄拉克δ函数的相应RS分解：$\delta = \mathcal{L}p_s + \mathcal{L}p_l = \delta_s + \delta_l$。
• 以低秩张量格式表示$\delta_s$和$\delta_l$，确保总近似误差受$\varepsilon$-截断秩和网格尺寸$n$的控制。
• 以$d$重$n \times \cdots \times n$笛卡尔网格上牛顿核（$1/\|x\|$）的规范RS张量表示作为分解基础。
• 利用径向函数如$1/\|x\|$和$1/\|x\|^3$可实现低秩RS张量分解的事实，从而高效计算流体动力学和弹性势能。
• 通过使用长程平滑分量对奇异电荷分布进行正则化，将该方法应用于泊松-玻尔兹曼方程，减轻了由点电荷奇异性引起的条件数问题。

实验结果

研究问题

• RQ1 在$\mathbb{R}^d$中，离散狄拉克δ函数能否通过与算子相关的范围分离张量格式实现高效近似？
• RQ2 格林核的RS分解是否能产生一种稳定且低秩的狄拉克δ表示，从而将奇异分量与光滑分量分离？
• RQ3 对于$N$个点电荷的集体狄拉克δ函数的长程部分，其张量秩是否仅随$N$的对数增长？
• RQ4 该张量方法如何改善含多个奇异源项的势方程求解的条件数和可扩展性？
• RQ5 该方法能否扩展至具有常系数的椭圆算子逆，如双调和算子或斯托克斯算子？

主要发现

• 通过将椭圆算子$\mathcal{L}$作用于格林核的RS分解，构建了与算子相关的狄拉克δ函数的RS分解，得到短程和长程张量分量之和。
• 长程分量$\delta_l$是光滑的，适合低秩张量近似，而短程分量$\delta_s$通过局部校正捕捉奇异性。
• 由于张量格式的特性，该方法在网格尺寸$n$上实现了$O(n)$复杂度量级，从而可在大规模3D网格上实现高效计算。
• 集体狄拉克δ函数长程部分的秩仅随粒子数$N$的对数增长，确保了对大规模$N$-粒子系统的可扩展处理。
• 数值测试验证了近似的局部化和精度，误差受$\varepsilon$-截断和网格尺寸的控制。
• 该方法成功应用于生物分子静电势的泊松-玻尔兹曼方程，提供了一种正则化方案，避免了点电荷奇异性导致的数值不稳定性。