[论文解读] Local structure of singular profiles for a Derivative Nonlinear Schrödinger Equation

Yuri Cher, Gideon Simpson+1 more|arXiv (Cornell University)|Feb 7, 2016

Advanced Mathematical Physics Problems参考文献 13被引用 9

一句话总结

本文通过渐近分析与数值延拓，研究了 $L^2$-超临界导数非线性薛定谔方程（gDNLS）中奇异轮廓的局部结构。推导出临界附近的通用爆破轮廓 $Q$ 及其参数 $a$、$b$ 和 $\text{Re}(Q)$，表明当 $\sigma \to 1^+$ 时 $a \to 0$，并利用牛顿求解器与理查德森外推法提高精度，数值计算了轮廓。

ABSTRACT

The Derivative Nonlinear Schrödinger equation is an $L^2$-critical nonlinear dispersive equation model for Alfvén waves in space plasmas. Recent numerical studies on an $L^2$-supercritical extension of this equation provide evidence of finite time singularities. Near the singular point, the solution is described by a universal profile that solves a nonlinear elliptic eigenvalue problem depending only on the strength of the nonlinearity. In the present work, we describe the deformation of the profile and its parameters near criticality, combining asymptotic analysis and numerical simulations.

研究动机与目标

• 理解 $L^2$-超临界导数非线性薛定谔方程（gDNLS）中爆破解在奇点附近的局部结构。
• 确定当非线性强度 $\sigma$ 趋近临界值 $\sigma = 1$ 时，爆破轮廓 $Q$ 及其参数 ($a$、$b$) 如何变形。
• 通过数值求解非线性特征值问题，建立轮廓 $Q$ 在不同初值条件下的普适性。
• 利用理查德森外推法提高 $a$ 和 $b$ 的数值精度，并与现有求解器进行结果验证。

提出的方法

• 使用牛顿求解器与稀疏直接线性求解器，在均匀网格上求解爆破轮廓 $Q$ 的非线性特征值问题。
• 通过 $W(x) = Q(-x)$ 的变换处理原点处的对称性与连续性，实现在 $[0, x_{\max}]$ 上的求解。
• 采用二阶中心差分法处理空间导数，并通过渐近衰减行为施加远场边界条件。
• 在原点处施加对称与反对称条件，并施加零相位条件 $v_0 = 0$，以确保 $|Q|$ 的最大值位于 $x=0$。
• 在 $\sigma$ 上采用延拓法，从 $\sigma = 2$ 开始逐步减小，将 $\sigma_{j-1}$ 处的解作为 $\sigma_j$ 的初始猜测。
• 应用理查德森外推法提升收敛阶数，使用公式 $a^{\text{Rich.}} = \frac{4}{3}a^{\Delta x/2} - \frac{1}{3}a^{\Delta x}$，并在 $\sigma$ 值之间使用三次样条插值。

实验结果

研究问题

• RQ1 在 $L^2$-超临界 gDNLS 方程中，当 $\sigma \to 1^+$ 时，爆破轮廓 $Q$ 及其参数 $a$、$b$ 如何演化？
• RQ2 如数值模拟所暗示，爆破轮廓 $Q$ 是否在不同初值条件下具有普适性？
• RQ3 在 $\sigma = 1$ 附近解析 $Q$ 时会遇到哪些数值挑战，如何缓解？
• RQ4 特别是原点附近与远场区域，$Q$ 的振幅与相位行为如何？
• RQ5 理查德森外推法在多大程度上能提升计算参数 $a$ 和 $b$ 的精度？

主要发现

• 爆破轮廓 $Q$ 具有普适性，与初值条件无关，当适当选择 $b$ 时，其最大振幅位于原点。
• 当 $\sigma \to 1^+$ 时，参数 $a \to 0$，表明爆破轮廓中的特征长度尺度缩小，与临界性一致。
• 数值解已计算至 $\sigma = 1.044$，受限于 $\Delta x \ll a$ 的要求，当 $a \to 0$ 时需更高分辨率。
• 远场区域中，轮廓 $Q$ 衰减至低于 $10^{-300}$，达到双精度浮点数精度的极限。
• 理查德森外推法显著提升了 $a$ 和 $b$ 的收敛阶数，结果经与 BVP_SOLVER-2 对比验证。
• 延拓法成功追踪了从 $\sigma = 2$ 到 $\sigma = 1.044$ 的轮廓，近临界区域步长自适应减小至 $\Delta\sigma = 0.00078125$。