On linear waves with memory in a Bessel-like medium

基本信息

• 论文ID: 2510.11493
• 标题: On linear waves with memory in a Bessel-like medium
• 作者: A. Giusti, I. Colombaro, A. Mentrelli
• 分类: math-ph math.MP
• 发表时间: October 14, 2025
• 论文链接: https://arxiv.org/abs/2510.11493

摘要

本文讨论了由含记忆的波动方程控制的系统中调和波和瞬态波的传播，该方程的积分核在拉普拉斯域中涉及第一类修正贝塞尔函数的比值。特别地，瞬态波的研究通过基于Talbot方法的全数值方法进行，该方法用于拉普拉斯变换的数值反演。

研究背景与动机

问题描述

本研究要解决的核心问题是在具有记忆效应的贝塞尔类介质中波的传播行为。具体而言，研究对象是形如以下的含记忆波动方程：

$\frac{\partial^2 Y}{\partial t^2} - c^2\left[1 - \Phi_\nu(t) * \right]\frac{\partial^2 Y}{\partial x^2} = 0$

其中记忆核在拉普拉斯域中表示为： $\tilde{\Phi}_\nu(s) = \frac{2(\nu+1)}{\sqrt{s\tau}}\frac{I_{\nu+1}(\sqrt{s\tau})}{I_\nu(\sqrt{s\tau})}$

重要性与应用

这类方程在多个物理领域都有重要应用：

1. 血液动力学：描述充满液体的弹性管道中的波传播
2. 粘弹性力学：材料的记忆效应建模
3. 电气系统：具有记忆特性的电路分析

现有方法的局限性

之前的研究主要采用Buchen-Mainardi半解析算法，但该方法仅在波前附近可靠，无法提供完整的时空演化图像。

研究动机

本文旨在通过数值方法克服现有方法的局限性，提供对整个时空域内波传播行为的完整分析。

核心贡献

1. 建立了完整的色散分析框架：针对ν=0情况，详细分析了空间衰减时间周期(SATP)调和波的色散关系
2. 推导了精确的色散律和衰减因子：通过Kelvin函数给出了显式表达式
3. 实现了全数值瞬态波分析：采用Talbot方法实现拉普拉斯反变换的数值计算
4. 揭示了异常色散特性：证明了系统表现出相速度小于群速度的异常色散现象

方法详解

任务定义

研究在贝塞尔类记忆介质中波的传播特性，包括：

• 调和波的色散和衰减特性分析
• 瞬态波的阶跃响应计算
• 相速度和群速度的确定

模型架构

1. 基本波动方程

考虑ν=0的特殊情况： $\frac{\partial^2 Y}{\partial t^2} - c^2\left[1 - \Phi(t) * \right]\frac{\partial^2 Y}{\partial x^2} = 0$

其中： $\tilde{\Phi}(s) = \frac{2}{\sqrt{s\tau}}\frac{I_1(\sqrt{s\tau})}{I_0(\sqrt{s\tau})}$

2. 色散分析方法

采用标准调和波解形式： $Y(t,x) = A e^{i(\omega t - kx)}$

得到色散关系： $(i\omega)^2 + k^2 c^2[1 - \hat{\Phi}(\omega)] = 0$

3. SATP波分析

对于空间衰减时间周期波，设定：

• $\omega \in \mathbb{R}$（实频率）
• $k = \kappa - i\delta_{att}$（复波数）

通过求解非线性方程组： $\begin{cases} \kappa^2 - \delta_{att}^2 = A(\omega) \\ -2\kappa\delta_{att} = B(\omega) \end{cases}$

其中A(ω)和B(ω)通过Kelvin函数表示。

技术创新点

1. Kelvin函数表示法

将复数贝塞尔函数比值转换为实Kelvin函数的组合，避免了复数运算的数值不稳定性。

2. Talbot方法应用

采用Talbot方法进行拉普拉斯反变换： $Y(t,x) = \frac{1}{2\pi i}\int_C e^{st}\tilde{Y}(s,x)ds$

通过变形积分路径，将Bromwich积分转换为实轴上的积分，提高数值稳定性。

3. 全域数值解

克服了传统方法仅在波前有效的局限，实现了整个时空域的精确计算。

实验设置

初边值问题

考虑半无限域上的阶跃响应问题： $\begin{cases} Y_{tt} - c^2[1 - \Phi(t)*]Y_{xx} = 0, & x > 0, t > 0 \\ Y(t,0) = H(t), & t > 0 \\ Y(0,x) = 0, Y_t(0,x) = 0, & x > 0 \end{cases}$

数值实现

• 采用Wolfram Mathematica®进行符号计算
• 有限差分法计算群速度的数值导数
• Talbot方法的高精度数值积分

参数设置

使用无量纲参数：

• 频率参数：$\omega\tau$
• 空间参数：$x/(c\tau)$
• 时间参数：$(ct-x)/(c\tau)$

实验结果

主要结果

1. 色散律和衰减因子

图1显示了无量纲色散律$c\tau\kappa$和衰减因子$c\tau\delta_{att}$随$\omega\tau$的变化：

• 色散律呈非线性增长
• 衰减因子在中等频率达到峰值后下降

2. 相速度和群速度

图2展示了相速度和群速度的比较：

• 相速度：$v_p(\omega) = \frac{\omega}{\kappa(\omega)}$
• 群速度：$v_g(\omega) = \left[\frac{d\kappa(\omega)}{d\omega}\right]^{-1}$
• 数值结果显示：$v_p(\omega) < v_g(\omega) < c$，证实了异常色散特性

3. 瞬态波传播

图3显示了不同位置的阶跃响应：

• 波前以速度c传播
• 波形在传播过程中发生显著色散
• 远离波前的区域呈现复杂的振荡衰减行为

重要发现

1. 异常色散确认

系统表现出典型的异常色散特征：

• $v_p < v_g < c$
• 两个速度在高频极限下都趋近于c

2. 记忆效应的影响

记忆项导致：

• 波的衰减和色散
• 复杂的瞬态行为
• 非指数型的波形演化

相关工作

历史发展

1. Mainardi等(1982-1983)：最早在血液动力学中引入此类模型
2. Giusti & Mainardi(2016)：建立了动态粘弹性类比
3. Colombaro等(2017)：基于贝塞尔函数的线性粘弹性模型

方法比较

• Buchen-Mainardi算法：半解析，仅波前有效
• 本文方法：全数值，整个域有效
• 渐近方法：高频近似，局限性强

结论与讨论

主要结论

1. 色散特性：系统表现出异常色散，相速度始终小于群速度且都不超过波前速度c
2. 数值方法有效性：Talbot方法成功克服了传统方法的局限性
3. 记忆效应：贝塞尔型记忆核导致复杂的波传播行为

局限性

1. 计算复杂度：瞬态阶段需要大量节点，计算成本高
2. 参数限制：仅考虑了ν=0的特殊情况
3. 理论分析不足：群速度的严格性质难以解析证明

未来方向

1. 混合方法：结合半解析方法处理波前，数值方法处理远场
2. 一般化：扩展到任意ν值的情况
3. 应用拓展：在具体物理系统中的应用验证

深度评价

优点

1. 方法创新：首次将Talbot方法应用于贝塞尔记忆介质的波传播问题
2. 理论完备：提供了从色散分析到瞬态计算的完整框架
3. 数值精度：克服了现有方法的局限性，实现全域高精度计算
4. 物理洞察：揭示了异常色散等重要物理现象

不足

1. 计算效率：瞬态计算的高复杂度限制了实际应用
2. 参数范围：仅限于特定的记忆核形式
3. 实验验证：缺乏与实际物理实验的对比

影响力

1. 理论贡献：为记忆介质中的波传播提供了新的数值分析工具
2. 应用价值：在生物医学、材料科学等领域具有潜在应用
3. 方法论意义：Talbot方法的应用为类似问题提供了参考

适用场景

1. 学术研究：波传播理论、记忆材料建模
2. 工程应用：管道流动分析、振动控制
3. 数值计算：拉普拉斯反变换的高精度计算

参考文献

论文引用了23篇重要文献，涵盖：

• 贝塞尔函数理论（Abramowitz & Stegun, 1972）
• 波传播理论（Mainardi, 2022; Whitham, 2011）
• 数值方法（Talbot, 1979; Garrappa等, 2021）
• 应用背景（血液动力学、粘弹性等相关文献）

本论文在理论分析和数值方法两方面都有重要贡献，为记忆介质中的波传播研究提供了新的工具和洞察。尽管在计算效率和参数范围方面存在一定局限，但其方法论价值和物理洞察使其成为该领域的重要进展。