耦合kdv方程的相互作用解
Interaction Solutions of the Coupled KdV Equations
DOI: 10.12677/AAM.2023.124185, PDF, HTML, XML, 下载: 155  浏览: 203  科研立项经费支持
作者: 董 露, 陈怀堂*:临沂大学数学与统计学院,山东 临沂
关键词: 耦合kdv方程相互作用解辅助方程法Coupled KdV Equations Interaction Solution Auxiliary Equation Method
摘要: 随着人们对非线性方程的深入研究,经常通过探求耦合kdv方程的相互作用解来描述水波等自然运动规律。本文提出了一种求耦合kdv方程相互作用解的辅助方程新方法。该新方法可以很容易得到三角函数、指数函数、双曲函数和其他函数的混合函数解。利用该方法,我们成功地得到了耦合kdv方程的相互作用解。这些解在帮助物理学家准确分析相关领域中的特殊现象方面具有十分重要的理论意义和应用价值。
Abstract: In this paper, based on the auxiliary equation method, we obtain new interaction solutions of the coupled KdV equations, these solutions are degenerated to the solitary wave solutions, the triangle function solutions and other function solutions. It is significant to help physicists to analyze special phenomena in their relevant fields accurately.
文章引用:董露, 陈怀堂. 耦合kdv方程的相互作用解[J]. 应用数学进展, 2023, 12(4): 1785-1791. https://doi.org/10.12677/AAM.2023.124185

1. 引言

非线性偏微分方程的相互作用解在非线性光学、理论物理、等离子体物理、流体动力学、半导体等领域有着重要的应用。非线性偏微分方程的相互作用解包括了三角函数、指数函数、双曲函数和其他函数的混合函数解对于一些物理现象的研究有重要意义,而对这些解析解的研究有助于我们更好地理解复杂的物理现象。近几十年来,人们提出了许多求非线性偏微分方程精确解的方法,如:1) 逆散射理论 [1] ,Hirotas双线性方法 [2] ;相对于逆散射理论方法而言,Hirotas双线性方法也被称为直接方法,这种方法的优点在于它是一种代数而不是解析的方法,这种方法已经从求KdV方程、MKDV方程等的孤立子解而发展成一种求解一大批非线性偏微分方程孤子解的相当普遍的方法,而这种方法的关键在于寻求相关变量变换。2) 截断Painlev展开法 [3] 、Darboux变换法 [4] ;Darboux变换法相对于截断Painlev展开法更加实用,Darboux变换是求解非线性微分方程的一种非常实用的方法,比多尺度法、双线性变换法等其他方法更实用、更简单。近年来,人们研究了大量求解非线性偏微分方程的有效方法,它们是齐次平衡法 [5] 、正弦余弦法 [6] 、sech函数法 [7] 、双曲正切函数法 [8] [9] 、多重exp函数法 [10] [11] 、雅可比椭圆函数展开法 [12] [13] 。而辅助方程法也是一种重要的方法,它在非线性微分方程的研究中扮演着很重要的角色,它不仅能求解一大批非线性偏微分方程的孤子解,也可以求几种类型的特殊解。它以简洁易懂的特质而备受关注。

众所周知,复杂的自然现象常常用非线性偏微分方程来描述。其中最具代表性的非线性方程是耦合KdV方程方程 [14] [15] 。本文利用一种新的辅助方程方法并结合齐次平衡法研究了耦合KdV方程:

u t + α 1 u x x x α 2 u u x + α 3 v v x = 0 (1.1)

v t + β 1 v x x x β 2 u v x β 3 v u x = 0 (1.2)

的精确解。

方程(1.1)和(1.2)描述了具有不同色散关系的两长波的相互作用。其中 α 1 , α 2 , α 3 β 1 , β 2 , β 3 都是常数。

2. 新辅助方程的新解

假设新的辅助方程为:

F ( ξ ) = A 1 + B 1 F ( ξ ) + C 1 F 3 ( ξ ) (2.1)

如果 A 1 = B 1 = 0 , C 1 0 ,则我们由maple软件得到方程(2.1)的新的解为:

F 1 ( ξ ) = 2 C 1 3 { ln [ ( tanh ( ξ ) + tan ( ξ ) ) e ξ 1 + tan ( ξ ) coth ( ξ ) + tanh ( ξ ) + tan ( ξ ) ] } 1 (2.2)

F 2 ( ξ ) = 2 C 1 3 { ln [ ( tanh ( ξ ) + coth ( ξ ) ) e ξ 1 + tanh ( ξ ) tan ( ξ ) + tan ( ξ ) + coth ( ξ ) ] } 1 (2.3)

F 3 ( ξ ) = 2 C 1 3 { ln [ ( tanh ( ξ ) + cot ( ξ ) ) e ξ 1 + cot ( ξ ) coth ( ξ ) + tanh ( ξ ) + cot ( ξ ) ] } 1 (2.4)

F 4 ( ξ ) = 2 C 1 3 { ln [ ( cot ( ξ ) + coth ( ξ ) ) e ξ 1 + tanh ( ξ ) cot ( ξ ) + cot ( ξ ) + coth ( ξ ) ] } (2.5)

F 5 ( ξ ) = 2 C 1 3 { ln [ ( 1 + tan ( ξ ) coth ( ξ ) ) e ξ 1 + tanh ( ξ ) + tan ( ξ ) + tan ( ξ ) coth ( ξ ) ] } 1 (2.6)

F 6 ( ξ ) = 2 C 1 3 { ln [ ( 1 + cot ( ξ ) coth ( ξ ) ) e ξ 1 + tanh ( ξ ) + cot ( ξ ) + cot ( ξ ) coth ( ξ ) ] } 1 (2.7)

F 7 ( ξ ) = 2 C 1 3 { ln [ ( 1 + tanh ( ξ ) cot ( ξ ) ) e ξ 1 + tan ( ξ ) + coth ( ξ ) + tanh ( ξ ) tan ( ξ ) ] } 1 (2.8)

F 8 ( ξ ) = 2 C 1 3 { ln [ ( 1 + tanh ( ξ ) cot ( ξ ) ) e ξ 1 + cot ( ξ ) + coth ( ξ ) + tanh ( ξ ) cot ( ξ ) ] } 1 (2.9)

在求解的过程中我们忽略了一些虚解,因此将该辅助方程应用到非线性偏微分方程中时,我们就仅得到了非线性偏微分方程的实相互作用解。

3. 新的辅助方程法

在上一节中,我们已经求出了辅助方程的新解,为了更好地求解耦合KdV方程的相互作用解,本文介绍了一种新的辅助方程法。这种辅助方程法的步骤如下所示:

第一步:对于给定的具有自变量 x , y , t 等的非线性偏微分方程:

{ P ( u t , u x x , u y y , v , ) = 0 , Q ( v t , v x x , v y y , u , ) = 0. (3.1)

我们作如下变换:

{ u = a 0 + a 1 F 2 ( ξ ) + a 2 G 2 ( ξ ) + a 3 F ( ξ ) G ( ξ ) , v = a 10 + a 11 F 2 ( ξ ) + a 12 G 2 ( ξ ) + a 13 F ( ξ ) G ( ξ ) . ξ = k 1 ( x c 1 t ) η = k 2 ( x c 2 t ) (3.2)

其中 a i 为待定系数, F ( ξ ) 满足(2.1), G ( η ) 满足下式。

G ( η ) = A 2 + B 2 G 2 ( η ) (3.3)

第二步:将变换(3.2)代入方程(3.1),再利用(2.1)和(3.1),则方程左端化为关于 F i ( ξ ) G j ( η ) ( F ( ξ ) ) k 的多项式,令 F i ( ξ ) G j ( η ) ( F ( ξ ) ) k 的所有系数都等于零,则得到关于待定系数的代数方程组。因此, a i 将通过解代数方程组来确定。然后我们将这个方法应用到耦合KdV方程中去,就能得到耦合KdV方程的相互作用解。

4. 耦合KdV方程的新相互作用解

对于求解耦合KdV方程的新相互作用解,

我们利用上述方法求解含有待定系数的代数方程组,

得到以下结果:

1)

a 0 = 1 β 1 ( α 2 β 2 ) ( 4 α 1 B 1 k 1 2 β 1 4 B 1 β 1 2 k 1 2 a 10 α 1 β 2 α 3 β 1 α 1 β 2 α 1 β 3 + α 2 β 1 a 10 α 1 β 3 α 3 β 1 α 1 β 2 α 1 β 3 + α 2 β 1 + a 10 α 2 β 1 α 3 β 1 α 1 β 2 α 1 β 3 + α 2 β 1 + a 10 β 1 β 3 α 3 β 1 α 1 β 2 α 1 β 3 + α 2 β 1 )

a 1 = a 11 α 3 β 1 α 1 β 2 α 1 β 3 + α 2 β 1

C 1 = α 11 ( β 2 + β 3 ) α 3 β 1 α 1 β 2 α 1 β 3 + α 2 β 1 6 β 1 k 1 2

A 1 = a 12 = a 13 = a 2 = a 3 = 0

c 1 = 1 β 1 ( α 2 β 2 ) ( 4 α 1 B 1 k 1 2 β 1 β 2 4 α 2 B 1 β 1 2 k 1 2 a 10 α 1 β 2 2 α 3 β 1 α 1 β 2 α 1 β 3 + α 2 β 1 a 10 α 1 β 2 β 3 α 3 β 1 α 1 β 2 α 1 β 3 + α 2 β 1 + a 10 α 2 β 1 β 2 α 3 β 1 α 1 β 2 α 1 β 3 + α 2 β 1 + a 10 α 2 β 1 β 3 α 3 β 1 α 1 β 2 α 1 β 3 + α 2 β 1 ) (4.1)

其中, B 1 , k 1 , a 0 , a 10 a 11 都是任意常数。

2)

A 2 = N 1 96 k 2 β 1 M

B 2 = a 12 2 α 3 β 2 2 + 2 a 12 2 α 3 β 2 β 3 + a 12 2 α 3 β 3 2 144 α 1 β 1 β 2 144 α 1 β 1 β 3 + 144 β 1 2 α 2 k 2

a 2 = 12 β 1 a 12 2 α 3 β 2 2 + 2 a 12 2 α 3 β 2 β 3 + a 12 2 α 3 β 3 2 144 α 1 β 1 β 2 144 α 1 β 1 β 3 + 144 β 1 2 α 2 β 2 + β 3

c 2 = N 2 12 M

A 1 = a 1 = a 11 = a 13 = a 3 = 0

其中,

M = ( α 1 β 2 + α 1 β 3 α 2 β 1 ) a 12 2 α 3 β 2 2 + 2 a 12 2 α 3 β 2 β 3 + a 12 2 α 3 β 3 2 144 α 1 β 1 β 2 144 α 1 β 1 β 3 + 144 β 1 2 α 2 ( α 1 β 1 )

N 1 = 12 a 0 α 1 α 2 β 1 β 2 a 12 2 α 3 β 2 2 + 2 a 12 2 α 3 β 2 β 3 + a 12 2 α 3 β 3 2 144 α 1 β 1 β 2 144 α 1 β 1 β 3 + 144 β 1 2 α 2 + 12 a 0 α 1 α 2 β 1 β 3 a 12 2 α 3 β 2 2 + 2 a 12 2 α 3 β 2 β 3 + a 12 2 α 3 β 3 2 144 α 1 β 1 β 2 144 α 1 β 1 β 3 + 144 β 1 2 α 2 12 a 0 α 1 β 1 β 2 2 a 12 2 α 3 β 2 2 + 2 a 12 2 α 3 β 2 β 3 + a 12 2 α 3 β 3 2 144 α 1 β 1 β 2 144 α 1 β 1 β 3 + 144 β 1 2 α 2 12 a 0 α 1 β 1 β 2 β 3 a 12 2 α 3 β 2 2 + 2 a 12 2 α 3 β 2 β 3 + a 12 2 α 3 β 3 2 144 α 1 β 1 β 2 144 α 1 β 1 β 3 + 144 β 1 2 α 2

12 a 0 α 2 2 β 1 2 a 12 2 α 3 β 2 2 + 2 a 12 2 α 3 β 2 β 3 + a 12 2 α 3 β 3 2 144 α 1 β 1 β 2 144 α 1 β 1 β 3 + 144 β 1 2 α 2 + 12 a 0 α 2 β 1 2 β 2 a 12 2 α 3 β 2 2 + 2 a 12 2 α 3 β 2 β 3 + a 12 2 α 3 β 3 2 144 α 1 β 1 β 2 144 α 1 β 1 β 3 + 144 β 1 2 α 2 a 10 a 12 α 3 α 1 β 2 2 2 a 10 a 12 α 3 α 1 β 2 β 3 a 10 a 12 α 3 α 1 β 3 2 + a 10 a 12 α 3 α 2 β 2 β 1 + a 10 a 12 α 3 α 2 β 1 β 3 + a 10 a 12 α 3 β 2 β 1 β 3 + a 10 a 12 α 3 β 1 β 3 2 ,

N 2 = 12 a 0 α 1 2 β 2 2 a 12 2 α 3 β 2 2 + 2 a 12 2 α 3 β 2 β 3 + a 12 2 α 3 β 3 2 144 α 1 β 1 β 2 144 α 1 β 1 β 3 + 144 β 1 2 α 2 + 12 a 0 α 1 2 β 2 β 3 a 12 2 α 3 β 2 2 + 2 a 12 2 α 3 β 2 β 3 + a 12 2 α 3 β 3 2 144 α 1 β 1 β 2 144 α 1 β 1 β 3 + 144 β 1 2 α 2 24 a 0 α 1 α 2 β 1 β 2 a 12 2 α 3 β 2 2 + 2 a 12 2 α 3 β 2 β 3 + a 12 2 α 3 β 3 2 144 α 1 β 1 β 2 144 α 1 β 1 β 3 + 144 β 1 2 α 2 12 a 0 α 1 β 1 α 2 β 3 a 12 2 α 3 β 2 2 + 2 a 12 2 α 3 β 2 β 3 + a 12 2 α 3 β 3 2 144 α 1 β 1 β 2 144 α 1 β 1 β 3 + 144 β 1 2 α 2 + 12 a 0 α 2 2 β 1 2 a 12 2 α 3 β 2 2 + 2 a 12 2 α 3 β 2 β 3 + a 12 2 α 3 β 3 2 144 α 1 β 1 β 2 144 α 1 β 1 β 3 + 144 β 1 2 α 2 + a 10 a 12 α 3 α 1 β 2 2 + a 10 a 12 α 3 α 1 β 2 β 3 a 10 a 12 α 3 α 2 β 2 β 1 a 10 a 12 α 3 α 2 β 3 β 1 , (4.2)

其中, a 0 , a 10 , a 12 , k 1 , k 2 , c 1 , B 1 , C 1 都是任意常数。

利用解1和(2.1)的解,我们得到耦合KdV方程的相互作用解如下所示:

{ u 1 = a 0 + a 1 { 2 C 1 3 { ln [ ( tanh ( ξ ) + tan ( ξ ) ) e ξ 1 + tan ( ξ ) coth ( ξ ) + tanh ( ξ ) + tan ( ξ ) ] } 1 } 2 , v 1 = a 10 + a 11 { 2 C 1 3 { ln [ ( tanh ( ξ ) + tan ( ξ ) ) e ξ 1 + tan ( ξ ) coth ( ξ ) + tanh ( ξ ) + tan ( ξ ) ] } 1 } 2 . (4.3)

{ u 2 = a 0 + a 1 { 2 C 1 3 { ln [ ( tanh ( ξ ) + coth ( ξ ) ) e ξ 1 + tanh ( ξ ) tan ( ξ ) + tan ( ξ ) + coth ( ξ ) ] } 1 } 2 , v 2 = a 10 + a 11 { 2 C 1 3 { ln [ ( tanh ( ξ ) + coth ( ξ ) ) e ξ 1 + tanh ( ξ ) tan ( ξ ) + tan ( ξ ) + coth ( ξ ) ] } 1 } 2 . (4.4)

{ u 3 = a 0 + a 1 { 2 C 1 3 { ln [ ( tanh ( ξ ) + cot ( ξ ) ) e ξ 1 + cot ( ξ ) coth ( ξ ) + tanh ( ( ξ ) + cot ( ξ ) ] } 1 } 2 , v 3 = a 10 + a 11 { 2 C 1 3 { ln [ ( tanh ( ξ ) + cot ( ξ ) ) e ξ 1 + cot ( ξ ) coth ( ξ ) + tanh ( ( ξ ) + cot ( ξ ) ] } 1 } 2 . (4.5)

{ u 4 = a 0 + a 1 { 2 C 1 3 { ln [ ( cot ( ξ ) + coth ( ξ ) ) e ξ 1 + tanh ( ξ ) cot ( ξ ) + cot ( ξ ) + coth ( ξ ) ] } 1 } 2 , v 4 = a 10 + a 11 { 2 C 1 3 { ln [ ( cot ( ξ ) + coth ( ξ ) ) e ξ 1 + tanh ( ξ ) cot ( ξ ) + cot ( ξ ) + coth ( ξ ) ] } 1 } 2 . (4.6)

{ u 5 = a 0 + a 1 { 2 C 1 3 { ln [ ( 1 + tan ( ξ ) coth ( ξ ) ) e ξ 1 + tanh ( ξ ) + tan ( ξ ) + tan ( ξ ) coth ( ξ ) ] } 1 } 2 , v 5 = a 10 + a 11 { 2 C 1 3 { ln [ ( 1 + tan ( ξ ) coth ( ξ ) ) e ξ 1 + tanh ( ξ ) + tan ( ξ ) + tan ( ξ ) coth ( ξ ) ] } 1 } 2 . (4.7)

其中, ξ = k 1 ( x c 1 t )

至此,我们已经求得了耦合KdV方程相互作用解。如(4.3)和(4.7)所示。若用代数方程组的解2,再利用(3.3)的解,还可得到耦合KdV方程的周期解和更多的孤立波解。

5. 总结

本文给出了一种新的辅助方程法,这种辅助方程法在求解耦合KdV方程的相互作用解以及其他非线性偏微分方程方面都有重要的应用。

致谢

本文受到临沂大学创新训练项目(X202210452171)资助。

参考文献

NOTES

*通讯作者。

参考文献

[1] Gardner, C.S., Greene, J.M., Kruskal, M.D. and Miura, R.M. (1967) Method for Solving the Korteweg-deVries Equa-tion. Physical Review Letters, 19, 1095-1097.
https://doi.org/10.1103/PhysRevLett.19.1095
[2] Hirota, R. (1971) Exact Solution of the Korteweg-deVries Equation for Multiple Collisions of Solutions. Physics Review Letters, 27, 1192-1194. 7
https://doi.org/10.1103/PhysRevLett.27.1192
[3] Weiss, J., Tabor, M. and Garnevale, G. (1983) The Painlevé Property for Partial Differential Equations. Journal of Mathematical Physics, 24, 522-526.
https://doi.org/10.1063/1.525721
[4] Li, Y.S. and Zhang, J.E. (2001) Darboux Transformations of Classical Boussinesq System and Its Multi-Soliton Solutions. Physics Letters A, 284, 253-258.
https://doi.org/10.1016/S0375-9601(01)00331-0
[5] Wang, M.L., Zhou, Y.B. and Li, Z.B. (1996) Application of a Homogeneous Balance Method to Exact Solutions of Nonlinear Equations in Mathematical Physics. Physics Letters A, 216, 67-75.
https://doi.org/10.1016/0375-9601(96)00283-6
[6] 徐炳振, 李悦科, 阎循领. 一类五阶非线性演化方程的新孤波解[J]. 物理学报, 1998, 47(12): 1946-1951.
[7] 陈德芳, 楼森岳. KdV方程与高阶KdV方程行波解之间的形变理论[J]. 物理学报, 1991, 40(4): 513-521.
[8] Ma, W.X. and Fuchssteiner, B. (1996) Explicit and Exact Solutions to a Kolmogorov-Petrovskii-Piskunov Equation. International Journal of Non-Linear Mechanics, 31, 329-338.
https://doi.org/10.1016/0020-7462(95)00064-X
[9] Blickle, T., Lakatos, B.G, Mihálykó, C. and Ulbert, Z. (1998) The Hyperbolic Tangent Distribution Family. Powder Technology, 97, 100-108.
https://doi.org/10.1016/S0032-5910(97)03393-7
[10] Ma, W.X. and Zhu, Z.N. (2012) Solving the (3 + 1)-Dimensional Generalized KP and BKP Equations by the Multiple Exp-Function Algorithm. Applied Mathematics and Computation, 24, 11871-11879.
https://doi.org/10.1016/j.amc.2012.05.049
[11] Ma, W.-X., Huang, T. and Zhang, Y. (2010) A Multiple Exp-Function Method for Nonlinear Differential Equations and Its Application. Physica Scripta, 82, 5468-5478.
https://doi.org/10.1088/0031-8949/82/06/065003
[12] Fu, Z.T., Liu, S.K., Liu, S.D. and Zhao, Q. (2001) New Ja-cobi Elliptic Function Expansion and New Periodic Solutions of Nonlinear Wave Equations. Physics Letters A, 290, 72-76.
https://doi.org/10.1016/S0375-9601(01)00644-2
[13] Chen, H.T. and Zhang, H.Q. (2004) New Double Pe-riodic and Multiple Soliton Solutions of the Generalized (2 + 1)-Dimensional Boussinesq Equation. Chaos, Solitons & Fractals, 20, 765-769.
https://doi.org/10.1016/j.chaos.2003.08.006
[14] 许斌, 刘希强, 刘玉堂. 耦合KdV方程组的对称,精确解和守恒律[J]. 应用数学学报, 2010, 33(1): 118-123.
[15] Chen, H.T. and Zhang, H.Q. (2003) Improved Jacobin Elliptic Function Method and Its Applications. Chaos, Solitons & Fractals, 15, 585-591.
https://doi.org/10.1016/S0960-0779(02)00147-9