[论文解读] Numerical study of blow-up in solutions to generalized Kadomtsev-Petviashvili equations

Christian Klein, Ralf Peter|arXiv (Cornell University)|Oct 19, 2013

Nonlinear Waves and Solitons参考文献 13被引用 4

一句话总结

本文通过动态再缩放方法对广义 Kadomtsev-Petviashvili (gKP) 方程在临界和超临界非线性情况下的解的爆破现象进行了数值研究，以分类爆破类型。研究识别出不同的爆破情形：对于 gKP I 方程在 $ n = 4/3 $ 时，爆破遵循 $ L^2 $-临界孤立子轮廓，且 $ \|u\|_{\infty} \propto (t^* - t)^{-3/4} $；而对于 $ n > 4/3 $，爆破由指数型缩放律描述，且 $ \|u\|_{\infty} \propto (t^* - t)^{-2/(3n)} $，与理论预期一致。

ABSTRACT

We present a numerical study of solutions to the generalized Kadomtsev-Petviashvili equations with critical and supercritical nonlinearity for localized initial data with a single minimum and single maximum. In the cases with blow-up, we use a dynamic rescaling to identify the type of the singularity. We present a discussion of the observed blow-up scenarios.

研究动机与目标

• 研究广义 Kadomtsev-Petviashvili (gKP) 方程解发生有限时间爆破的条件。
• 利用动态再缩放技术对爆破情形中的奇点类型进行分类。
• 确定在 $ n > 2 $ 时 gKP II 方程是否发生爆破，若发生，则表征其渐近行为。
• 研究初始数据的质量与能量在触发爆破中的作用，特别是与孤立子解的关系。
• 通过数值方法验证 gKP I 和 gKP II 方程中爆破的猜想渐近缩放律。

提出的方法

• 使用谱方法在空间上求解 gKP 方程，并采用分裂步 Fourier 方法结合自适应时间步长进行时间推进。
• 应用动态再缩放技术，将爆破问题转化为稳态问题，从而识别爆破类型。
• 使用再缩放变换 $ u(x,y,t) \to L(t)^{-2/n} u(x/L(t), y/L^2(t), t^*) $，其中 $ L(t) \propto (t^* - t)^{\alpha} $，以提取缩放指数。
• 将 $ \ln \|u\|_{\infty} $ 和 $ \ln \|u_y\|_2^2 $ 拟合为 $ c \ln(t^* - t) + C $ 形式，以确定爆破速率。
• 构造初始数据为 $ u_0(x,y) = \beta \partial_{xx} \exp(-x^2 - y^2) $，确保 $ \partial_x^{-1} u \in \mathcal{S}(\mathbb{R}^2) $，并满足约束条件 $ \int \partial_{yy} u \, dx = 0 $。
• 分析 $ u_y $ 的 $ L^2 $-范数以确认爆破，因为 $ \|u_y\|_2 \to \infty $ 在 $ t^* $ 处成立，与理论预期一致。

实验结果

研究问题

• RQ1 在初始数据满足什么条件下，gKP I 和 gKP II 方程在 $ n > 4/3 $ 和 $ n > 2 $ 时会发生爆破？
• RQ2 在爆破时间 $ t^* $ 附近，$ \|u\|_{\infty} $ 和 $ \|u_y\|_2 $ 的渐近缩放行为如何？
• RQ3 在 $ n = 4/3 $ 时，gKP I 的爆破轮廓是否与再缩放后的孤立子解一致？
• RQ4 对于 $ n > 4/3 $ 的 gKP I 和 $ n > 2 $ 的 gKP II，爆破是否遵循如猜想所示的指数型缩放？
• RQ5 初始数据的质量与能量如何影响爆破的发生，特别是与孤立子解相比？

主要发现

• 对于 gKP I 在 $ n = 4/3 $ 时，当初始数据具有足够小的能量和较大的质量时发生爆破，且 $ \|u\|_{\infty} \propto (t^* - t)^{-3/4} $，与 $ L^2 $-临界孤立子轮廓一致。
• $ u_y $ 的 $ L^2 $-范数以 $ \|u_y\|_2 \propto (t^* - t)^{-1} $ 的方式爆破，与 $ L^2 $-临界爆破的理论预测一致。
• 对于 $ n > 4/3 $ 的 gKP I，爆破速率为 $ \|u\|_{\infty} \propto (t^* - t)^{-2/(3n)} $，表明其为非孤立子型、自相似的轮廓，且在再缩放变量中呈现指数型缩放。
• 对 $ \ln \|u\|_{\infty} $ 拟合为 $ c_1 \ln(t^* - t) + C_1 $ 得到 $ c_1 = -0.1623 $，非常接近预期的 $ -1/6 $，从而证实了 $ n = 4/3 $ 时的 $ -2/(3n} $ 缩放律。
• 对于 $ n > 2 $ 的 gKP II，当质量较大且能量较小时观察到爆破，且 $ \|u\|_{\infty} \propto (t^* - t)^{-2/(3n)} $，表明两种 gKP 类型中存在统一的爆破机制。
• 对于 $ n > 4/3 $ 的 $ \|u_y\|_2 $，其缩放行为为 $ \|u_y\|_2 \propto (t^* - t)^{-(1+4/n)/6} $，拟合结果得到 $ c_2 = -0.858 $，接近预期的 $ -2/3 $，支持了猜想的缩放律。