(1 + 1)维Benjiamin Omo方程的精确行波解研究
Study of the Exact Traveling Wave Solutions of (1+1) Dimensional Benjiamin Omo Equation

作者: 刘红霞 * , 韩青秀 , 伍 芸 :贵州师范大学数学科学学院,贵州 贵阳;

关键词: 动力系统Benjiamin Omo方程周期波孤立波Dynamical Systems Benjiamin Omo Equation Solitarywaves Periodic Waves

摘要: 通过行波变换,利用微分方程定性理论与动力系统分支方法,我们得到了(1 + 1)维Benjiamin Omo方程在不同参数形式下的线性波解。

Abstract: Through the traveling wave transformation, by exploiting the qualitative theory of differential equations and the bifurcation method of dynamical systems, we obtained some nonlinearwave solutions for Benjiamin Omo equation under different parameter forms after travelingwave transformation.

1. 研究背景和预备知识

我们考虑如下的Benjiamin Omo方程

u t t + β ( u 2 ) x x + γ u x x x x = 0. (1)

其中 β , γ 是任意常数。

(1 + 1) 维Benjiamin Omo方程 [1] [2] [3] 是物理学中一个十分重要的非线性发展方程,已经有许多人对这一模型进行了广泛的研究.如文献 [4] 中用基于改进的投影Riccati法求得了该模型的单变量行波解。文献 [5] 通过Jacobi椭圆函数展开法,求得了该方程的精确周期解。在文献 [6] [7] 中,利用双线性方法和扩展的测试函数求得了其周期孤波解和双周期孤波解。文献 [8] 中通过改进的辅助方程法,求得了该方程的精确解。而本文运用平面动力系统的分支方法定性理论 [9],讨论方程(1)的非线性波解。

为了研究方程(1)的非线性波,设 c > 0 是波速, u ( x , t ) = φ ( ξ ) ξ = x c t 。将 u = φ ( ξ ) 代入方程(1)可得

c 2 φ + 2 β ( φ ) 2 + 2 β φ φ + γ φ ' = 0. (2)

对(2)式关于 ξ 积分两次可得

c 2 φ + β φ 2 + γ φ = 0. (3)

φ = y ,得到下面的平面系统

{ d φ d ξ = y d y d ξ = 1 γ ( β φ 2 + c 2 φ ) . (4)

很明显,系统(4)是一个有着Hamiltonian函数的Hamiltonian系统

H ( φ , y ) = y 2 + 2 β 3 γ φ 3 + c 2 γ φ 2 . (5)

如果令 f ( φ ) = 1 γ ( β φ 2 + c 2 φ )

Δ = c 4 γ 2

则我们有下面的结果。

Δ > 0 时, f ( φ ) 有两个零点 φ 1 φ 2 ,它们的表达式为

φ 1 = 0 , φ 2 = c 2 β . (6)

Δ 0 时, f ( φ ) 没有零点。

如果 φ ˜ 是系统(4)的一个奇点,那么在 ( φ ˜ , 0 ) 的特征值为

λ = ± f ( φ ˜ ) .

利用微分方程动力系统的定性理论,我们有下面的结论。

(1) 如果 f ( φ ˜ ) > 0 ,则 ( φ ˜ , 0 ) 是一个鞍点。

(2) 如果 f ( φ ˜ ) < 0 ,则 ( φ ˜ , 0 ) 是一个中心。

(3) 如果 f ( φ ˜ ) = 0 ,则 ( φ ˜ , 0 ) 是一个退化鞍点。

2. Benjiamin Omo方程的行波解的参数表达式

命题1. 当 γ > 0 , β > 0 时,方程(1)有两个孤立波解 u 1 ± 和两个周期波解 u 2 ±

u 1 ± = p ( 1 3 tanh 2 ( 1 2 η 1 | ξ | ) ) (7)

u 2 ± = a 2 ( a 3 a 1 ) a 2 ( a 3 a 2 ) s n 2 ( β ( a 3 a 1 ) 6 γ | ξ | , a 3 a 2 a 3 a 1 ) a 3 a 1 ( a 3 a 2 ) s n 2 ( β ( a 3 a 1 ) 6 γ | ξ | , a 3 a 2 a 3 a 1 ) (8)

其中

p = c 2 2 β ,

η 1 = c 2 γ ,

a 1 + a 2 + a 3 = 3 c 2 2 β ,

a 1 a 2 + a 1 a 3 + a 2 a 3 = 0 ,

a 1 a 2 a 3 = 0.

图1所示

证明。在(5)式中,令 H ( φ 2 , 0 ) = h 1 ,则有

Γ 1 ± : y 2 = 2 β 3 γ ( φ 3 3 c 2 2 β φ 2 + 3 γ 2 β h 1 ) = 2 β 3 γ ( φ φ 2 ) 2 ( p φ ) (9)

由(4)式得

d φ d ξ = ± 2 β 3 γ ( φ φ 2 ) 2 ( p φ ) (10)

积分(10)式得

Figure 1. Track chart of the Γ 1 ± and Γ 2 ±

图1. Γ 1 ± Γ 2 ± 的轨道图

φ 2 φ d φ ( φ φ 2 ) p φ = ± 0 ξ 2 β 3 γ d ξ (11)

由(11)式得到方程的孤立波解的表达式为 u 1 ±

同理,在(5)式中,令 H ( φ 1 , 0 ) = 0 ,则有

Γ 2 ± : y 2 = 2 β 3 γ ( φ 3 3 c 2 2 β φ 2 ) = 2 β 3 γ ( φ a 1 ) ( φ a 2 ) ( a 3 φ ) (12)

由(4)式可得

d φ d ξ = ± 2 β 3 γ ( φ a 1 ) ( φ a 2 ) ( a 3 φ ) . (13)

积分(13)式得

a 2 φ d φ ( φ a 1 ) ( φ a 2 ) ( a 3 φ ) = ± 0 ξ 2 β 3 γ d ξ . (14)

由(14)式得到方程的周期波解的表达式为 u 2 ±

命题2. 当 γ < 0 , β > 0 时,方程(1)有两个孤立波解 u 3 ± 和两个周期波解 u 4 ±

u 3 ± = 3 c 2 2 β ( tanh 2 ( 1 2 η 2 | ξ | ) 1 ) (15)

u 4 ± = ( a 2 a 1 ) s n 2 ( β ( a 3 a 1 ) 6 γ | ξ | , a 2 a 1 a 3 a 1 ) + a 1 (16)

其中 η 2 = c 2 γ

图2所示

证明。在(5)式中,令 H ( φ 1 , 0 ) = 0 ,则有

Γ 3 ± : y 2 = 2 β 3 γ ( φ 3 + 3 c 2 2 β φ 2 ) = 2 β 3 γ φ 2 ( φ + 3 c 2 2 β ) (17)

Figure 2. Track chart of the Γ 3 ± and Γ 4 ±

图2. Γ 3 ± Γ 4 ± 的轨道图

由(4)式得

d φ d ξ = ± 2 β 3 γ φ 2 ( φ + 3 c 2 2 β ) (18)

积分(18)式得

3 c 2 2 β φ d φ φ φ + 3 c 2 2 β = ± 0 ξ 2 β 3 γ d ξ (19)

由(19)式得到方程的孤立波解的表达式为 u 3 ±

同理,在(5)式中,令 H ( φ 2 , 0 ) = h 1 ,则有

Γ 4 ± : y 2 = 2 β 3 γ ( φ 3 + 3 c 2 2 β φ 2 3 γ 2 β h 1 ) = 2 β 3 γ ( φ a 1 ) ( a 2 φ ) ( a 3 φ ) (20)

由(4)式可得

d φ d ξ = ± 2 β 3 γ ( φ a 1 ) ( a 2 φ ) ( a 3 φ ) (21)

积分(21)式得

a 1 φ d φ ( φ a 1 ) ( a 2 φ ) ( a 3 φ ) = ± 0 ξ 2 β 3 γ d ξ (22)

由(22)式可得到方程的周期波解的表达式为 u 4 ±

命题3。当 γ > 0 , β < 0 时,方程(1)有两个孤立波解(与 u 1 ± 有着相同形式)和两个周期波解(与 u 4 ± 有着相同形式)。

图3所示

证明。在(5)式中,令 H ( φ 2 , 0 ) = h 1 ,则有

Γ 5 ± : y 2 = 2 β 3 γ ( φ 3 + 3 c 2 2 β φ 2 3 γ 2 β h 1 ) = 2 β 3 γ ( φ 2 φ ) 2 ( φ p ) (23)

Figure 3. Track chart of the Γ 5 ± and Γ 6 ±

图3. Γ 5 ± Γ 6 ± 的轨道图

由(4)式得

d φ d ξ = ± 2 β 3 γ ( φ 2 φ ) 2 ( φ p ) (24)

由(24)式得到方程的孤立波解的参数表达式与 u 1 ± 式形式相同。

同理,在(5)式中,令 H ( φ 1 , 0 ) = 0 ,则有

Γ 6 ± : y 2 = 2 β 3 γ ( φ 3 + 3 c 2 2 β φ 2 ) = 2 β 3 γ ( φ a 1 ) ( a 2 φ ) ( a 3 φ ) (25)

由(4)式可得

d φ d ξ = ± 2 β 3 γ ( φ a 1 ) ( a 2 φ ) ( a 3 φ ) (26)

由(26)式可得到方程的周期波解的表达式为与 u 4 ± 形式相同 。

命题4。当 γ < 0 , β < 0 时,方程(1)有两个孤立波解(与 u 3 ± 有着相同形式)和两个周期波解(与 u 2 ± 有着相同形式)。

图4所示

证明。在(5)式中,令 H ( φ 1 , 0 ) = 0 ,则有

Γ 7 ± : y 2 = 2 β 3 γ ( φ 3 3 c 2 2 β φ 2 ) = 2 β 3 γ φ 2 ( 3 c 2 2 β φ ) (27)

由(4)式得

d φ d ξ = ± 2 β 3 γ φ 2 ( 3 c 2 2 β φ ) (28)

由(28)式得到方程的孤立波解的表达式与 u 3 ± 式形式相同。

同理,在(5)式中,令 H ( φ 2 , 0 ) = h 1 ,则有

Figure 4. Track chart of the Γ 7 ± and Γ 8 ±

图4. Γ 7 ± Γ 8 ± 的轨道图

Γ 8 ± : y 2 = 2 β 3 γ ( φ 3 3 c 2 2 β φ 2 + 3 γ 2 β h 1 ) = 2 β 3 γ ( φ a 1 ) ( φ a 2 ) ( a 3 φ ) (29)

由(4)式可得

d φ d ξ = ± 2 β 3 γ ( φ a 1 ) ( φ a 2 ) ( a 3 φ ) . (30)

由(30)式得到方程的周期波解的表达式与 u 2 ± 形式相同。

综上所述,在不同的参数条件下我们获得了Benjiamin Omo方程的多个行波解的精确解。

NOTES

*通讯作者。

文章引用: 刘红霞 , 韩青秀 , 伍 芸 (2021) (1 + 1)维Benjiamin Omo方程的精确行波解研究。 应用数学进展, 10, 3084-3090. doi: 10.12677/AAM.2021.109322

参考文献

[1] Benjiamin, T.B. (1967) Internal Waves of Permanent Formin Fluids of Great Depth. Journal of Fluid Mechanics, 29, 559-592.
https://doi.org/10.1017/S002211206700103X

[2] Davies, R.E. and Acrivos, A. (1967) Solitary Internal Waves Indeep Water. Journal of Fluid Mechanics, 29, 593-607.
https://doi.org/10.1017/S0022112067001041

[3] Ono, H. (1975) Algebraic Solitary Waves in Stratifiedfluids. Journal of the Physical Society of Japan, 39, 1082-1091.
https://doi.org/10.1143/JPSJ.39.1082

[4] Wang, Z., Li, D.-S., Lu H.-F., et al. (2005). A Method for Constructing Exact Solutions and Application to Benjamin Ono Equation. Chinese Physics, 14, 2158-2158.
https://doi.org/10.1088/1009-1963/14/11/003

[5] Fu, Z. and Liu, S. (2003) The JEFE Method and Periodic Solutions of Two Kinds of Nonlinear Wave Equations. Communications in Nonlinear Science and Numerical Simulation, 8, 67-75.
https://doi.org/10.1016/S1007-5704(02)00082-5

[6] Xu, Z.H., Xian, D.Q. and Chen, H.L. (2010) New Periodic Solitary-Wave Solutions for the Benjiamin Ono Equation. Applied Mathematics and Computation, 215, 4439-4442.
https://doi.org/10.1016/j.amc.2009.11.009

[7] 许瑞麟, 许晓革, 孟祥花. (1+1)维BenjiaminOno方程的精确显式解[J]. 北京信息科技大学学报(自然科学版), 2011, 26(3): 58-61.

[8] 杨娟, 冯庆江, 曾春花. 一类非线性数学物理方程的精确解[J]. 数学的实践与认识, 2018, 48(17): 309-313.

[9] Wu, Y., Liu, Z.R. and Bai, C.Z. (2013) New Types of Nonlinear Waves and Bifurcation Phenomena in Schamel-Korteweg-de Vries Equation. Abstract and Applied Analysis, 2013, 1056-1083.
https://doi.org/10.1155/2013/483492

分享
Top