[论文解读] Numerical study of blow-up in solutions to generalized Korteweg-de Vries equations

Christian Klein, Ralf Peter|arXiv (Cornell University)|Jul 2, 2013

Advanced Mathematical Physics Problems参考文献 4被引用 19

一句话总结

本论文对具有临界（$n=4$）和超临界（$n=5$）非线性项的广义Korteweg-de Vries（gKdV）方程中的爆破现象进行了数值研究。通过动态重标度和谱方法，证实了解析解表现出与Martel、Merle和Raphaël所预测一致的自相似爆破形态，爆破时间 $t^*$ 对色散参数 $\epsilon$ 呈指数依赖关系，且在 $\epsilon \to 0$ 的极限下，最小爆破时间 $t_0^* \approx 1.0536$ 大于Hopf临界时间。即使初始数据的质量远超孤立子质量，研究也证实了爆破动力学的普遍性。

ABSTRACT

We present a detailed numerical study of solutions to general Korteweg-de Vries equations with critical and supercritical nonlinearity. We study the stability of solitons and show that they are unstable against being radiated away and blow-up. In the $L_{2}$ critical case, the blow-up mechanism by Martel, Merle and Raphaël can be numerically identified. In the limit of small dispersion, it is shown that a dispersive shock always appears before an eventual blow-up. In the latter case, always the first soliton to appear will blow up. It is shown that the same type of blow-up as for the perturbations of the soliton can be observed which indicates that the theory by Martel, Merle and Raphaël is also applicable to initial data with a mass much larger than the soliton mass. We study the scaling of the blow-up time $t^{*}$ in dependence of the small dispersion parameter $ε$ and find an exponential dependence $t^{*}(ε)$ and that there is a minimal blow-up time $t^{*}_{0}$ greater than the critical time of the corresponding Hopf solution for $ε o0$. To study the cases with blow-up in detail, we apply the first dynamic rescaling for generalized Korteweg-de Vries equations. This allows to identify the type of the singularity.

研究动机与目标

• 研究广义Korteweg-de Vries（gKdV）方程在 $n \geq 4$ 条件下的解的稳定性和爆破动力学，尤其关注 $L^2$-临界（$n=4$）和超临界（$n=5$）情形。
• 通过数值方法验证Martel、Merle和Raphaël针对 $L^2$-临界gKdV方程提出的理论爆破机制。
• 分析爆破时间 $t^*$ 对小色散参数 $\epsilon$ 的依赖关系，特别是在 $\epsilon$ 较小时的行为。
• 确定在初始质量远大于孤立子质量的条件下，是否也能观察到与孤立子扰动中相同的爆破形态，从而验证爆破动力学的普遍性。

提出的方法

• 应用动态重标度方法，将gKdV方程转化为可追踪爆破奇点的形式。
• 采用带去混叠处理的伪谱傅里叶谱方法求解空间变量，以高精度解析振荡和爆破动力学。
• 实施自适应时间步长和分辨率控制，以在解趋近爆破时保持精度，特别关注傅里叶系数模长以检测分辨率损失。
• 将对数范数（$\ln||u_x||_2^2$，$\ln||u||_\infty$）拟合至 $\alpha \ln(t^* - t) + \beta$，以估计 $t^*$ 并验证自相似爆破行为。
• 对 $\ln(t^*/t_c)$ 进行最小二乘拟合至 $\gamma \epsilon + \delta$，以确定 $t^*$ 在 $\epsilon \to 0$ 时的渐近行为，揭示其指数缩放特性。
• 对 $n=4$ 和 $n=5$ 情况下，初始数据 $u_0 = \text{sech}^2 x$ 的数值模拟，通过改变 $\epsilon$ 观察色散冲击波的形成和爆破序列。

实验结果

研究问题

• RQ1 Martel、Merle和Raphaël针对 $L^2$-临界gKdV方程（$n=4$）所预测的爆破机制，是否在初始数据不必然接近孤立子的解中通过数值方法得以显现？
• RQ2 在 $n=4$ 和 $n=5$ 情况下，爆破时间 $t^*$ 对小色散参数 $\epsilon$ 的依赖关系如何？当 $\epsilon \to 0$ 时，是否表现出指数缩放？
• RQ3 在初始质量远大于孤立子质量的解中，是否也能观察到与孤立子扰动中相同的自相似爆破形态，从而表明爆破动力学具有普遍性？
• RQ4 在小-$\epsilon$ 范围内，色散冲击波是否总是先于爆破出现？与相应Hopf解的临界时间 $t_c$ 相比，这如何影响爆破时间？
• RQ5 在 $\epsilon \to 0$ 的极限下，最小爆破时间 $t_0^*$ 是否严格大于相应Hopf解的临界时间 $t_c$？

主要发现

• 对于 $n=4$，数值结果证实了Martel、Merle和Raphaël所预测的自相似爆破形态，满足 $||u_x||_2 \sim ||Q'||_2 / (l_0 (t^* - t))$，当 $t \to t^*$ 时，解收敛至缩放后的孤立子形态。
• 对于 $n=5$，数值结果表明爆破动力学遵循与 $n=4$ 情况相同的普遍形态，表明Martel、Merle和Raphaël的理论在初始质量远大于孤立子质量的情况下依然适用。
• 对于 $n=4$ 和 $n=5$，爆破时间 $t^*$ 对 $\epsilon$ 呈指数依赖关系，$\ln(t^*/t_c)$ 拟合为 $\gamma \epsilon + \delta$ 的线性关系，相关系数 $r=0.9999$，表明 $\lim_{\epsilon \to 0} t^* = t_0^* \approx 1.0536 > t_c \approx 0.5341$。
• 在小-$\epsilon$ 极限下，色散冲击波总是先于爆破形成，且冲击结构中首先出现的孤立子正是最终发生爆破的那个。
• 对 $\ln||u_x||_2^2$ 和 $\ln||u||_\infty$ 拟合至 $\alpha \ln(t^* - t) + \beta$ 所得的 $t^*$ 估计值在不同范数下具有一致性，且由于重标度时间中的指数缩放，$n=5$ 情况下的收敛性优于 $n=4$ 情况。
• 当 $\epsilon = 0.01$ 时，$||u_x||_2^2$ 的爆破时间估计为 $t^* \approx 1.1475$，$||u||_\infty$ 的估计为 $t^* \approx 1.1434$，证实了结果的一致性及渐近区的开始。