Analytical and Gegenbauer spectral solutions of damped wave equations with singular potentials in ℝn

技术摘要：具有奇异势能的阻尼波方程在 $\mathbb{R}^n$ 中的解析与 Gegenbauer 谱解法

问题表述
本研究探讨了定义在矩形区域 $\Omega = [0, \pi]^n$ 上的多维阻尼波方程的存在性、唯一性及其数值近似问题。所控制的方程是一个包含形式为 $|x|^{-\alpha}$ 的奇异源项的混合双曲-抛物型算子：
$$\frac{\partial^2 u}{\partial t^2} + \gamma \frac{\partial u}{\partial t} - \Delta u + \frac{u}{|x|^\alpha} = f(x, t)$$
该问题在 Cauchy 初值条件（$u(x,0)=\phi(x)$，$u_t(x,0)=\psi(x)$）和非齐次 Neumann 边界条件（$\partial u / \partial \nu = g_j$）下进行表述。分析重点在于奇异参数满足 $0 < \alpha < n$ 且阻尼系数 $\gamma > 0$ 的情形。奇异势能的存在为可积性和正则性带来了显著的分析挑战，而该混合算子的性质则要求同时处理波传播与扩散耗散。

方法论
本文采用了双重解析与数值方法：

1. 解析框架（特征函数展开）：

   • 将解分解为瞬态齐次部分 ($u_h$) 和稳态特解部分 ($u_p$)。
   • 利用时间傅里叶级数构建特解，将随时间变化的 PDE 简化为每个时间谐波 $k$ 对应的耦合椭圆方程组。
   • 通过引入复值空间函数 $U_k(x) = A_k(x) - iB_k(x)$，将系统简化为复型 Helmholtz 型方程：$(L - \lambda_k)U_k = F_k$，其中 $L = -\Delta + |x|^{-\alpha}$ 且 $\lambda_k = k^2 - i\gamma k$。
   • 解通过使用 Neumann 拉普拉斯算子的特征函数 $\Phi_m(x) = \prod \cos(m_j x_j)$ 表示为双级数。系数通过向正交基进行投影来确定，并考虑了内部源项和边界通量。

2. 数值框架（Gegenbauer 谱配置法）：

   • 开发了一种使用参数 $\lambda=1$ 的 Gegenbauer 多项式（对于最优节点分布，这会退化为切比雪夫多项式）的高阶谱方法。
   • 空间离散利用在张量积网格上的切比雪夫-高斯-洛巴托 (CGL) 节点来处理定义域几何结构。
   • 使用 $V_\epsilon(x) = (|x|^2 + \epsilon^2)^{-\alpha/2}$ 对原点附近的奇异势 $|x|^{-\alpha}$ 进行正则化，以防止数值溢出同时保持渐近行为。
   • 时间积分采用二阶隐式-显式 (IMX) 方案，对刚性扩散项进行隐式处理，对双曲项进行显式处理，以确保在不受到严格 CFL 条件限制的情况下保持稳定性。

主要贡献与理论结果
本文确立了该解的严谨数学性质：

• 系数收敛性： 证明了展开系数渐近衰减为 $O(|m|^{-n-\alpha})$，从而确保尽管存在奇异性，级数仍然是良定义的，前提是 $0 < \alpha < n$。
• 稳定性与非共振性： 证明了阻尼系数 $\gamma > 0$ 至关重要。它确保了复谱分母 $|\mu_m - \lambda_k| \geq \gamma k > 0$，从而防止了除以零的情况并消除了纯共振现象。
• 正则性与 Sobolev 空间：
  • 解 $u$ 被证明在 $L^2(\Omega)$ 中是有界的且适定的。
  • 空间梯度 $\nabla u$ 仅在更严格的条件 $0 < \alpha < n-2$ 下保持平方可积（即 $u \in H^1(\Omega)$）。对于 $n-2 \leq \alpha < n$，梯度在原点处表现出奇异性。
• 全局稳定性： 推导出了一个显式的范数估计，表明 $L^2$ 范数随时间呈指数级衰减 $e^{-\gamma t/2}$，证实了全局稳定性和能量耗散。
• 误差分析： 建立了 Gegenbauer 配置法的理论误差界。空间离散误差以 $O(N^{-(n-\alpha-1/2)})$ 的速率代数收敛，反映了由奇异性引起的正则性损失。对于弱奇异性，可以观察到近谱收敛。

数值结果
在 一、二、三维空间中进行了全面的数值实验。

• 验证： 模拟验证了理论收敛率。对于弱奇异性（$\alpha = 0.3$），该方法实现了高阶收敛（速率 $> 3.5$）。对于接近临界阈值的强奇异性（$\alpha \approx n - 0.5$），收敛率平滑地降至理论代数极限。
• 能量行为： 数值图表确认了能量随时间的指数衰减，这与阻尼参数 $\gamma$ 一致。
• 维度缩放： 结果表明，该方法在 $n=1, 2, 3$ 维下均保持鲁棒且稳定，计算成本随 $O(N^n)$ 缩放。

意义
本文声称填补了文献中的空白，为具有奇异势能的 $n$ 维混合双曲-抛物型算子提供了显式的解析表示，而标准存在性定理通常无法提供此类构造性解。研究表明，阻尼系数在存在奇异性的情况下起到了根本的稳定作用，能够抵御共振。此外，这种集成的解析-数值方法为表征非均匀介质中具有局部奇异性的波传播与扩散提供了一个鲁棒的模型。作者将这些发现定位为可直接应用于热弹性学、流固耦合以及耗散介质中的电磁波传播等物理领域。