[论文解读] Decay of the Kolmogorov $N$-width for wave problems

Constantin Greif, Karsten Urban|arXiv (Cornell University)|Mar 20, 2019

Numerical methods in inverse problems参考文献 9被引用 4

一句话总结

本文证明了具有不连续初值的双曲波方程解流形的柯尔莫哥洛夫 $N$-宽度的衰减速率至多为 $\mathcal{O}(N^{-1/2})$，表明基于投影的模型降阶方法具有较慢的逼近速率。该结果通过在解空间中构造正交函数并建立有限维子空间中最佳逼近误差的界推导得出。

ABSTRACT

The Kolmogorov $N$-width $d_N(\mathcal{M})$ describes the rate of the worst-case error (w.r.t.\ a subset $\mathcal{M}\subset H$ of a normed space $H$) arising from a projection onto the best-possible linear subspace of $H$ of dimension $N\in\mathbb{N}$. Thus, $d_N(\mathcal{M})$ sets a limit to any projection-based approximation such as determined by the reduced basis method. While it is known that $d_N(\mathcal{M})$ decays exponentially fast for many linear coercive parametrized partial differential equations, i.e., $d_N(\mathcal{M})=\mathcal{O}(e^{-βN})$, we show in this note, that only $d_N(\mathcal{M}) =\mathcal{O}(N^{-1/2})$ for initial-boundary-value problems of the hyperbolic wave equation with discontinuous initial conditions. This is aligned with the known slow decay of $d_N(\mathcal{M})$ for the linear transport problem.

研究动机与目标

• 分析具有不连续初值的参数化双曲波问题的柯尔莫哥洛夫 $N$-宽度的衰减速率。
• 确定基于投影的模型降阶方法（如降基法）是否能为这类问题实现快速（指数级）收敛。
• 将已知的线性输运问题中缓慢的 $N$-宽度衰减结果扩展到二阶双曲波方程。
• 为解流形的任意 $N$-维线性子空间建立最坏情况逼近误差的下界。
• 证明 $\mathcal{O}(N^{-1/2})$ 的缓慢衰减速率是具有不连续数据的波问题的固有特性，与输运方程中的行为一致。

提出的方法

• 在 $\mathbb{R} \times \mathbb{R}^+$ 上建立二阶波方程的初值-边值问题，初值不连续。
• 使用分布解的概念，将解流形 $\mathcal{M}$ 定义为所有满足 $\mu \geq 0$ 的解 $u_\mu$ 的集合。
• 在 $L^2(\Omega_I)$ 中构造一组有限正交函数 $\tilde{\psi}_i$，以从下方界定 $N$-宽度。
• 应用引理 4.1，将有限集的 $N$-宽度与标准基 $\{e_1, \dots, e_{2N}\}$ 的 $N$-宽度联系起来，得到 $d_N(\tilde{\Psi}_{2N}) = \frac{1}{\sqrt{2}}$。
• 通过不等式链 $d_N(\mathcal{M}) \geq d_N(\Phi_{2N}) \geq \frac{1}{2} d_N(\Psi_{2N})$ 推导出最终的下界。
• 结合正交函数构造与尺度分析，得到 $d_N(\mathcal{M}) \geq \frac{1}{4} N^{-1/2}$。

实验结果

研究问题

• RQ1 具有不连续初值的双曲波方程解流形的柯尔莫哥洛夫 $N$-宽度的衰减速率是什么？
• RQ2 在线性输运问题中观察到的缓慢 $N^{-1/2}$ 衰减是否也出现在二阶双曲PDE中？
• RQ3 基于投影的逼近方法能否为这类波问题实现指数级收敛？
• RQ4 $N^{-1/2}$ 衰减速率是否是双曲问题模型降阶的根本限制？
• RQ5 初值中存在不连续性如何影响解流形的逼近性质？

主要发现

• 具有不连续初值的波方程解流形的柯尔莫哥洛夫 $N$-宽度至多以 $\mathcal{O}(N^{-1/2})$ 速率衰减。
• 对于具有不连续初值的波方程，严格建立了下界 $d_N(\mathcal{M}) \geq \frac{1}{4} N^{-1/2}$。
• 衰减速率 $\mathcal{O}(N^{-1/2})$ 与线性输运问题一致，表明其在逼近速度上存在共同的根本限制。
• 该结果意味着此类波问题的降基法无法实现指数级收敛，限制了其效率。
• 这种缓慢的衰减是解流形的内在特性，即使在具有平滑跃变的连续初值下依然存在。
• 证明依赖于在 $L^2$ 中构造一组有限正交函数，并利用其性质从下方界定 $N$-宽度。