非线性薛定谔方程的呼吸子解及其怪波极限 下载: 1341次
1 引言
在现代科学与技术的发展过程中,系统的非线性特性一直受到广泛关注,非线性特性是理解很多自然现象的基础,在流体力学、化学、生物和非线性光学等领域均可以观察到,其重要性不容置疑。一般情况下,非线性现象利用线性和非线性效应的非线性演化方程进行建模,标准非线性薛定谔方程(NLSE)就是描述非线性现象的核心模型之一[1-2]。
自20世纪60年代以来,非线性薛定谔方程已成为广泛研究的偏微分方程,一维自聚焦的非线性薛定谔方程被认为是描述海洋波的动力学模型[3-10]。非线性薛定谔方程之所以被广泛研究,是因为其在海洋动力学之外的其他物理学分支中均有重要应用,如非线性光学[11-12]、Bose-Einstein凝聚[13]以及其他非线性物理学领域。NLSE的精确解包括孤子[1]、Peregrine怪波(RW)[14-15]、Akhmediev呼吸子(AB)[16-17]、Kuznetsov-Ma孤子(KMS)[18-19]以及这些解的不同组合[20]。
目前,非线性薛定谔方程的呼吸子解被认为是在海洋以及其他物理领域中怪波的潜在原型[21-23],当呼吸子周期趋于无穷大时,呼吸子解将成为有理分式的怪波解,因此,有关呼吸子的结果也将为怪波研究提供理论参考。1991年,呼吸子解首次被Kichenassamy[24]求得;随后几十年,呼吸子解分别在非线性薛定谔方程[25]、Korteweg-de Vries方程[26-27]、二维Fermi-Pasta-Ulam晶格[28]、向列相液晶盒[29]和掺杂光纤[14]等系统中获得;呼吸子由调制不稳定性(MI)[14,30-31]导致的周期性扰动所产生,这些扰动通常是混沌的,可能在频谱中包含多个频率成份,因此,如何构建高阶解成为研究者非常关心的问题,Akhmediev小组认为高阶解可以是几个低阶解的非线性叠加[32-33]。
在自然界中,怪波是一种突然出现但短时间内又突然消失并具有极强破坏力的波[34],然而这种破坏性至今无法完全解释,航海家们一致认为线性理论无法解释怪波[35],只有非线性理论才能解释这种能量集中且峰值幅度通常为有效波高两倍以上的现象。在最近的研究中,除NLSE方程之外,Hirota方程[36]、Fokas-Leneels方程[37]、Hirota Maxwell-Bloch(MB)方程[38]、Sasa-Satsuma方程[39]以及耗散系统[40]中均已证明怪波的存在。经过研究,科学家们一致认为调制不稳定性是形成怪波最根本的动力学原因[41],而高阶怪波是由具有相同或特定频率的多个呼吸子之间的相互作用形成[42],即当两个或多个呼吸子碰撞时,会产生时空均局域且幅值较大的怪波[43]。另外,对怪波形成机制的理解非常有用,在海洋中可以有效预防怪波的产生,从而避免造成危害;在光学领域的光纤系统中产生怪波,可以用于产生高能量脉冲。
本文基于标准非线性薛定谔方程,采用达布变换方法推导其一阶和高阶呼吸子的精确表达形式,并详细讨论怪波极限及其他动力学特性。
2 NLSE的一呼吸子解及其怪波极限
标准的无量纲NLSE的表达式为
式中
式中
当
式中
图 1. NLSE方程一呼吸子解。(a) a1=-0.5,b1=0.8;(b) a1=0.5,b1=1
Fig. 1. One-breather solution of NLSE. (a) a1=-0.5, b1=0.8; (b) a1=0.5, b1=1
为了在坐标平面内得到时间和空间上都局域的波——Peregrine怪波,NLSE方程的一阶呼吸子解的周期应趋于无穷,即当
图 2. Peregrine孤子解。(a) x1=0,t1=0;(b) x1=3,t1=3
Fig. 2. Peregrine soliton solution. (a) x1=0, t1=0; (b) x1=3, t1=3
3 NLSE的双呼吸子解及其怪波极限
为了得到NLSE方程的双呼吸子解,需要两个复特征值,即
当
式中下标
由
图 3. 双呼吸子的碰撞。(a) b1=0.8;(b) b1→1;(c)双呼吸子碰撞的避免
Fig. 3. Collision of double breather. (a) b1=0.8; (b) b1→1; (c) avoid the collision of double breather
当
图 4. 呼吸子的碰撞与分离。(a)双AB碰撞;(b)双AB分离;(c) AB与KM孤子碰撞
Fig. 4. Collision and separation of breathers. (a) Double Akhmediev breather collision; (b) double Akhmediev breather separation; (c) collision between Akhmediev breather and Kuznetsov-Ma soliton
图 5. 二阶怪波的产生与分裂。(a)二阶怪波;(b)二维图;(c)分裂
Fig. 5. Generation and separation of second-order rouge waves. (a) Second-order rouge waves; (b) two-dimensional diagram; (c) separation
当
由于
图 6. 双呼吸子碰撞。(a) κ1∶κ2=2κ∶κ; (b) κ1∶κ2=3κ∶2κ; (c) κ1∶κ2=4κ∶3κ
Fig. 6. Double breather collision. (a) κ1∶κ2=2κ∶κ; (b) κ1∶κ2=3κ∶2κ; (c) κ1∶κ2=4κ∶3κ
从
双呼吸子的碰撞特性具有极强的丰富性,本文只讨论几种具有代表性的碰撞关系。已经讨论了一般的双呼吸子的碰撞,接下来讨论双呼吸子碰撞的另一种特殊情况,即两个呼吸子并行传播时的碰撞特性[46]。在并行传输模式下,双呼吸子传播需要满足条件:
根据文献[ 47]中对各参数的定义可以得到:
式中
图 7. l1=-1.1705+1.4i,l2=-0.5+0.4i时双呼吸子的并行叠加。(a) x1=0,x2=0;(b)局部放大图;(c) x1=-2,x2=2
Fig. 7. Two breathers parallel superimposed with l1=-1.1705+1.4i and l2=-0.5+0.4i. (a) x1=0, x2=0; (b) partial enlargement; (c) x1=-2, x2=2
4 NLSE方程的三阶呼吸子解及其怪波极限
由(2)式可知,当
式中
当
同样地,当
图 8. 三阶呼吸子的碰撞。(a)双重碰撞;(b)三重碰撞;(c)简并态
Fig. 8. Collisions of third-order breathers. (a) Double collision; (b) triple collision; (c) degeneration
图 9. 三阶怪波的产生与分裂。(a)三阶怪波;(b)等高图;(c)分裂
Fig. 9. Generation and separation of third-order rouge waves. (a) Third-order rouge waves; (b) contours; (c) separation
图 10. 三阶呼吸子的碰撞。(a) κ1∶κ2∶κ3=2∶3∶5且κ1=0.8;(b) κ1∶κ2∶κ3=2∶3∶5且κ1=0.6;(c) κ1∶κ2∶κ3=3∶2∶1且κ1=1.8
Fig. 10. Third-order breather collision. (a) κ1∶κ2∶κ3=2∶3∶5 and κ1=0.8; (b) κ1∶κ2∶κ3=2∶3∶5 and κ1=0.6; (c) κ1∶κ2∶κ3=3∶2∶1 and κ1=1.8
由
由二阶呼吸子解及怪波极限的分析可知,三阶呼吸子并行传输需要满足条件:
假设
图 11. 三阶呼吸子的并行传输。(a) x1=0,x2=0,x3=0;(b)局部放大图;(c) x1=-4,x2=0,x3=4
Fig. 11. Third-order breather parallel transmission. (a) x1=0, x2=0, x3=0; (b) partial enlargement; (c) x1=-4, x2=0, x3=4
5 NLSE的四阶呼吸子解及其怪波极限
由(2)式可知,当
式中
图 12. 四阶呼吸子的相互作用。(a)双重碰撞;(b)四重碰撞
Fig. 12. Fourth-order breather interaction. (a) Double collision; (b) fourfold collision
特别地,当
当
当四阶怪波解具有较小的坐标偏移时,四阶怪波将发生分裂,为了方便分析,将特征频率以及坐标偏移量分别改写为
式中参数
图 13. 呼吸子的简并。(a) AB简并与KM孤子简并叠加;(b) AB的简并
Fig. 13. Degenerate breather. (a) Superposition of Akhmediev breather and Kuznetsov-Ma soliton degenerate state; (b) degenerate Akhmediev breather
图 14. 四阶怪波的产生。(a)四阶怪波;(b)二维图;(c)等高图
Fig. 14. Generation of fourth-order rouge wave. (a) Fourth-order rouge waves; (b) two-dimensional diagram; (c) contour map
图 15. 四阶怪波的分裂。(a)一次分裂;(b)两次分裂
Fig. 15. Separation of fourth-order rouge wave. (a) Once separation; (b) twice separation
同样地,由于
四阶呼吸子并行传输需要满足条件:
同样地,假设
由
图 16. 四阶呼吸子相互作用的等高图。(a) κ1∶κ2∶κ3∶κ4=1∶3∶5∶7;(b) κ1∶κ2∶κ3∶κ4=2∶4∶6∶8
Fig. 16. Contour map of fourth-order breather interaction. (a) κ1∶κ2∶κ3∶κ4=1∶3∶5∶7; (b) κ1∶κ2∶κ3∶κ4=2∶4∶6∶8
图 17. 四阶呼吸子的并行传输。(a) x1=0,x2=0,x3=0,x4=0;(b)局部放大图;(c) x1=-6,x2=-2,x3=2,x4=6
Fig. 17. Third-order breather parallel transmission. (a) x1=0, x2=0, x3=0, x4=0; (b) partial enlargement; (c) x1=-6, x2=-2, x3=2, x4=6
6 结论
基于标准无量纲的非线性薛定谔方程,采用达布变换方法得到了NLSE方程各阶呼吸子解的一般表达形式,并详细研究了怪波极限及其他动力学特性。当呼吸子在坐标平面内发生有效碰撞或周期趋于无穷时,会产生一个对应阶数的怪波。
然而,高能量脉冲在光纤光学中有用,海洋中高能量脉冲具有极强的破坏力,应当有效预防和避免。研究结果表明,合理设定参数既可以有效控制高阶怪波的产生,还可以选择怪波出现的位置。此外,还研究了高阶呼吸子的并行传输模式和频率比对传输特性的影响,以及高阶怪波的特性和一次分裂、二次分裂。本文只研究了标准的NLSE方程,后续工作将针对其他NLSE方程进行,期望详尽了解呼吸子以及怪波的动力学特性。
[2] AkhmedievN, AnkiewiczA. Solitons: nonlinear pulses and beams[M]. London: Chapman and Hall, 1997.
[10] PerrieW. Nonlinear ocean waves[M]. Southampton: Computational Mechanics Publications, 2006.
[12] Agrawal GP. Nonlinear fiber optics[M]. 4 th ed. Amsterdam: Elsevier , 2006.
[15] 武达, 王娟芬, 石佳, 等. 掺杂光纤中Peregrine孤子的产生和传输[J]. 光学学报, 2017, 37(4): 0406002.
[19] 杨光晔, 李禄, 田晋平. 基于谱过滤方法的Kuznetsov-Ma孤子向准基态孤子转化研究[J]. 光学学报, 2016, 36(6): 0619002.
[22] Shats M, Punzmann H, Xia H. Capillary rogue waves[J]. Physical Review Letters, 2010, 104(10): 104503.
[28] Sarkar R, Dey B. Exact compact breather-like solutions of two-dimensional Fermi-Pasta-Ulam lattice[J]. Journal of Physics A: Mathematical and General, 2006, 39(4): L99-L104.
[29] 朱叶青, 胡巍. 无外置偏压的向列相液晶盒中的呼吸子传输[J]. 光学学报, 2015, 35(9): 0919001.
[30] Benjamin T B, Feir J E. The disintegration of wave trains on deep water. Part 1. Theory[J]. Journal of Fluid Mechanics, 1967, 27(3): 417-430.
[31] 卢洵, 王东升. 色散管理孤子系统的调制不稳定性[J]. 激光技术, 2012, 36(4): 557-561.
[32] Akhmediev N N, Korneev V I, Mitskevich N V. Modulation instability in an optical fiber induced by cross-phase modulation[J]. Radiophysics and Quantum Electronics, 1991, 34(1): 73-77.
[33] Akhmediev N. Soto-Crespo J M, Ankiewicz A. Extreme waves that appear from nowhere: on the nature of rogue waves[J]. Physics Letters A, 2009, 373(25): 2137-2145.
[34] Akhmediev N, Ankiewicz A, Taki M. Waves that appear from nowhere and disappear without a trace[J]. Physics Letters A, 2009, 373(6): 675-678.
[35] Shukla P K, Kourakis I, Eliasson B, et al. Instability and evolution of nonlinearly interacting water waves[J]. Physical Review Letters, 2006, 97(9): 094501.
[36] Tao Y S, He J S. Multisolitons, breathers, and rogue waves for the Hirota equation generated by the Darboux transformation[J]. Physical Review E, 2012, 85(2): 026601.
[37] He J S, Xu S W, Porsezian K. Rogue waves of the Fokas-Lenells equation[J]. Journal of the Physical Society of Japan, 2012, 81(12): 124007.
[38] Li C, He J, Porsezian K, et al. Rogue waves of the Hirota and the Maxwell-Bloch equations[J]. Physical Review E, 2013, 87(1): 012913.
[39] Bandelow U, Akhmediev N. Sasa-Satsuma equation: soliton on a background and its limiting cases[J]. Physical Review E, 2012, 86(2): 026606.
[40] Soto-Crespo J M, Grelu P, Akhmediev N. Dissipative rogue waves: extreme pulses generated by passively mode-locked lasers[J]. Physical Review E, 2011, 84: 016604.
[41] Bespalov V I, Talanov V I. Filamentary structure of light beams in nonlinear liquids[J]. Zhetf Pisma Redaktsiiu, 1966, 3(12): 471.
[42] Voronovich V V, Shrira V I, Thomas G. Can bottom friction suppress 'freak wave' formation?[J]. Journal of Fluid Mechanics, 2008, 604: 263-296.
[43] He J S, Zhang H R, Wang L H, et al. Generating mechanism for higher-order rogue waves[J]. Physical Review E, 2013, 87(5): 052914.
[44] Kedziora D J, Ankiewicz A, Akhmediev N. Circular rogue wave clusters[J]. Physical Review E, 2011, 84(5): 056611.
[45] Kedziora D J, Ankiewicz A, Akhmediev N. Classifying the hierarchy of nonlinear-Schrödinger-equation rogue-wave solutions[J]. Physical Review E, 2013, 88: 013207.
[46] Chowdury A, Kedziora D J, Ankiewicz A, et al. Breather solutions of the integrable quintic nonlinear Schrödinger equation and their interactions[J]. Physical Review E, 2015, 91(2): 022919.
[47] Chowdury A, Krolikowski W. Breather-to-soliton transformation rules in the hierarchy of nonlinear Schrödinger equations[J]. Physical Review E, 2017, 95(6): 062226.
杜志峰, 宋丽军, 王艳. 非线性薛定谔方程的呼吸子解及其怪波极限[J]. 激光与光电子学进展, 2019, 56(5): 050601. Zhifeng Du, Lijun Song, Yan Wang. Breather Solutions and Their Rouge Wave Limits of NonlinearSchrödinger Equation[J]. Laser & Optoelectronics Progress, 2019, 56(5): 050601.