[论文解读] Solitary waves in one-dimensional pre-stressed lattice and its continual analog

Vsevolod Vladimirov, Sergii Skurativskyi|arXiv (Cornell University)|Dec 18, 2015

Advanced Mathematical Physics Problems参考文献 8被引用 3

一句话总结

本文研究了在一维预压颗粒链的连续类比中，高度局域化且稳定的行进波——紧凑子解。通过多尺度渐近展开推导出非线性PDE，作者识别出亮型与暗型紧凑子解，利用分析与数值方法证明其稳定性，并表明这些紧凑子在碰撞后仍能保持形状，与离散颗粒链中观测到的动力学行为高度一致。

ABSTRACT

One of the most interesting phenomena occuring in nonlinear media models is the existence of wave patterns, such as kinks, solitons, compactons, peakons and many others. There are known numerous nonlinear evolutionary PDEs, supporting soliton (multi-soliton) and compacton traveling wave (TW) solutions. Unfortunately, the vast majority of the models, with the exception of completely integrable ones, do not enable to analyze the properties of solitary waves interaction using only qualitative methods. Therefore it is instructive, when dealing with the non-integrable PDEs, to combine the qualitative treatment with numerical simulations. In this report we are going to present the results of studying compacton solutions in the continual models for granular pre-stressed chains. The model is shown to possess a pair of compacton TW solutions which are the bright and dark compactons. First we consider the stability properties of the compacton solutions and show that both the bright and the dark compactons pass the stability test. Next we analyze the dynamics of the compactons, simulating numerically the temporal evolution of a single compacton, a well as the interection of pairs of compactons, including bright-bright, dark-dark and bright-dark pairs. To be able to simulate the evolutin of interacting compactons, we have modified the numerical scheme built by J. de Frutos, M. A. Lopez-Marcos, and J. M. Sanz-Serna. Results of simulations are compared with that of evolution of corresponding impulse in the granular pre-stressed chain.

研究动机与目标

• 分析从一维预压颗粒链导出的连续模型中紧凑子解的存在性与稳定性。
• 研究非可积PDE中的紧凑子是否能表现出孤子样行为，如碰撞后形状保持不变。
• 比较连续模型中紧凑子的动力学与底层离散颗粒链的数值模拟结果。
• 评估连续近似的有效性与准确性，以捕捉离散预压系统中局域脉冲的行为。

提出的方法

• 使用多尺度渐近展开与连续极限，从预压颗粒链的离散ODE系统推导出非线性PDE。
• 应用形式化的多尺度分解，获得支持紧凑子解的修正PDE，其时间尺度较长而空间尺度较短。
• 构建哈密顿表示并利用变分原理，将紧凑子解识别为拉格朗日泛函的极值点。
• 使用[2, 7, 6]中的方法进行稳定性分析，检测随时间增长的扰动。
• 基于de Frutos、López-Marcos与Sanz-Serna的方法，实现一种改进的数值格式，用于模拟紧凑子的演化与碰撞。
• 将连续PDE的数值结果与相同初始条件下离散颗粒链的直接模拟结果进行比较。

实验结果

研究问题

• RQ1 在预压颗粒链的连续类比中是否存在紧凑子解？其解析形式为何？
• RQ2 所导出的紧凑子解在小扰动下是否稳定？在传播过程中是否保持形状？
• RQ3 相互作用的紧凑子——亮-亮、暗-暗及亮-暗对——的动力学行为如何，特别是碰撞后如何？
• RQ4 连续PDE模型在多大程度上能准确再现离散颗粒链中观测到的局域脉冲动力学？
• RQ5 为何标准Nesterenko方程中的紧凑子解不通过稳定性检验，而多尺度分析导出的PDE中的解却保持稳定？

主要发现

• 多尺度推导出的PDE支持一对紧凑子解：亮型与暗型紧凑子，二者均通过稳定性检验。
• 数值模拟表明，单个紧凑子及相互作用对（亮-亮、暗-暗、亮-暗）在碰撞后均能完全恢复其形状。
• 在物理空间中，紧凑子解的宽度为Δx ≈ 4.96a（当n = 3/2时，对应Hertzian接触），与实验和数值结果一致。
• 通过直接比较数值模拟结果，证实连续模型能准确再现离散颗粒链中局域脉冲的动力学。
• 标准Nesterenko方程中的紧凑子解不稳定，但通过多尺度分析得到的修正PDE却能产生稳定的紧凑子。
• 数值格式成功适配以模拟紧凑子相互作用，从而能够详细研究其碰撞动力学。