[论文解读] Projected dynamics of kinetic equations with energy diffusion in spaces of orthogonal polynomials

Jon Wilkening, Antoine Cerfon+1 more|arXiv (Cornell University)|Feb 12, 2014

Model Reduction and Neural Networks被引用 4

一句话总结

本文提出使用非经典正交多项式——特别是对非经典权函数正交的多项式——来对速度空间中具有能量扩散的动能方程进行投影。结果表明，对于轻微奇异的初始条件，该方法在相同计算成本下相比经典Hermite多项式可实现数量级更高的精度，且精确解能以舍入误差范围内的精度返回投影子空间，而Hermite多项式则会出现显著偏差。

ABSTRACT

We examine a partial differential equation describing energy diffusion in velocity space, analyzing its projected dynamics in finite-dimensional subspaces of orthogonal polynomials. Comparison with an expensive but highly accurate spectral transform approach reveals that for a large class of (mildly singular) initial conditions, using polynomials that are orthogonal with respect to a non-classical weight function can be many orders of magnitude more accurate than classical Hermite polynomials for the same computational work. For the new polynomials, the exact solution is found to leave the subspace for a time, but return (up to roundoff accuracy) to the same point evolved to by the projected dynamics in that time. By contrast, using classical polynomials, the exact solution differs significantly from the projected dynamics solution when it arrives in the subspace. We also explore the connection between eigenfunctions of the projected evolution operator and (non-normalizable) eigenfunctions of the full evolution operator, as well as the effect of truncating the computational domain. Techniques for suppressing floating-point arithmetic errors are discussed, including use of the singular value decompositon to solve a better-conditioned version of the stiffness-mass generalized eigenvalue problem.

研究动机与目标

• 通过使用非经典正交多项式，提高速度空间中具有能量扩散的动能方程数值解的精度。
• 解决经典Hermite多项式在长时间投影解时出现显著误差的局限性。
• 分析投影系统的动力学行为，并与精确解进行比较，特别关注子空间返回与稳定性。
• 研究投影算子的本征函数与完整演化算子的本征函数之间的关系。
• 通过奇异值分解缓解刚度-质量广义特征值问题中的浮点数算术误差。

提出的方法

• 将动能方程投影到由非经典权函数导出的正交多项式张成的有限维子空间上。
• 将投影动力学与高精度谱变换方法进行比较，以评估解的精度。
• 使用奇异值分解将刚度-质量广义特征值问题重新表述为更良态的形式，以减少数值误差。
• 分析精确解随时间相对于投影子空间的行为，特别是其在离开子空间后返回子空间的情况。
• 研究投影演化算子的谱性质，并将其与完整算子的非可归一化本征函数相关联。
• 截断计算域并评估其对解精度与稳定性的影响。

实验结果

研究问题

• RQ1 对于相同计算量，非经典正交多项式能否在求解能量扩散方程时显著优于经典Hermite多项式？
• RQ2 为何使用非经典多项式时，精确解在短暂离开子空间后能以舍入误差范围内的精度返回子空间，而使用Hermite多项式时却不能？
• RQ3 投影演化算子的本征函数与完整演化算子的非可归一化本征函数之间有何关系？
• RQ4 域截断对投影解的精度与稳定性有何影响？
• RQ5 如何有效抑制刚度-质量广义特征值问题中的浮点数算术误差？

主要发现

• 对于一大类轻微奇异的初始条件，非经典正交多项式在相同计算成本下产生的解比经典Hermite多项式高出多个数量级的精度。
• 当使用非经典多项式进行投影时，精确解在短暂离开子空间后能以舍入误差范围内的精度返回同一子空间点，表明投影具有高度保真性。
• 相比之下，使用Hermite多项式时，精确解在重新进入子空间后与投影动力学解出现显著偏离，表明近似质量较差。
• 投影算子的本征函数被发现与完整演化算子的非可归一化本征函数密切相关，表明该方法具有谱一致性。
• 通过奇异值分解预处理刚度-质量广义特征值问题，能有效抑制浮点数算术误差，提升数值稳定性。
• 域截断会引入误差，但非经典多项式方法相比经典方法仍保持鲁棒性与高精度。