本文在降阶法的基础上运用变量代换法研究了一类广义非线性KdV方程的精确解,获得了这个方程具有不同物理结构包括紧孤子、孤立子、孤立波相似解和周期解在内的行波解。对于这类具有正或负n指数的广义KdV方程,方程中各项的系数连同波速一起决定着解的物理结构。 The aim of this paper is devoted to the study of physical structures of solutions for generalized nonlinear KdV-type equations in two-dimensional space. The variable replacement method is used to get compactons, solitons, solitary patterns and periodic solutions. Furthermore, we point out that the different exponents and coefficients lead to different results of physical structures for this kind of equations with positive or negative exponents.
陈静
西南科技大学理学院,四川 绵阳
收稿日期:2018年11月23日;录用日期:2018年12月17日;发布日期:2018年12月24日
本文在降阶法的基础上运用变量代换法研究了一类广义非线性KdV方程的精确解,获得了这个方程具有不同物理结构包括紧孤子、孤立子、孤立波相似解和周期解在内的行波解。对于这类具有正或负n指数的广义KdV方程,方程中各项的系数连同波速一起决定着解的物理结构。
关键词 :广义非线性KdV方程,紧孤子,孤立子
Copyright © 2018 by author and beplay安卓登录
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
近几十年来,对刻画自然界中非线性波现象的偏微分方程的研究已引起众多学者的高度关注,特别是为平衡弱非线性和色散而产生的孤立子解更是被广泛地研究。为了获得非线性方程的精确解,大量方法已经被用于处理相应的偏微分方程(参见文 [
通过对非线性色散方程的研究,Rosenau和Hyman在文 [
u t + ( u m ) x + ( u m ) x x x = 0 , m > 1 , 1 < n ≤ 3. (1)
在文 [
Wazwaz在文 [
u t t − u x x − u y y − a ( u 2 n ) x x − b [ u n ( u n ) x x ] x x = 0 , n > 1 , (2)
u t t − u x x − u y y − a ( u − 2 n ) x x − b [ u − n ( u − n ) x x ] x x = 0 , n > 1. (3)
Wazwaz总结了正余弦拟设法的主要步骤,并用它获得了方程(2) (3)的紧孤子、孤立子、孤立波相似解和周期解。他进一步指出,方程中函数的指数和系数比a/b、以及正负号的变化都会导致解的物理结构产生质的变化。在文 [
本文的主要工作是寻求另外的方法去扩展Wazwaz在文 [
u t + a u x + b ( u n ) x + b ( u n ) y + ( u n ) x x x + ( u n ) y y y = 0 , n > 0 , (4)
u t + a u x + b ( u − n ) x + b ( u − n ) y + ( u − n ) x x x + ( u − n ) y y y = 0 , n > 1 , (5)
其中a,b为非零常数。方程(4) (5)可看作是K(n,n)方程添加了新的一项 a u x 并将它们推广到(2 + 1)维。运用变量 ξ = μ x + η y − c t 代换,我们把非线性偏微分方程(4) (5)转化成更容易求解的关于 ξ 的ODES,并因此得到方程(4) (5)行波解的解析表达式。
本文中我们引入了如下q型双曲函数
sinh q ( x ) = exp ( x ) − q exp ( − x ) 2 , cosh q ( x ) = exp ( x ) + q exp ( − x ) 2 ,
csch q ( x ) = 1 sinh q ( x ) , sech q ( x ) = 1 cosh q ( x ) ,
这里 q > 0 ,且为常数。
我们首先考虑具有正指数的方程(4)。假设方程(4)具有形如 u ( x , y , t ) = u ( ξ ) 的解,其中
ξ = μ x + η y − c t , μ , η , c 是不为0的常数,且 μ + η ≠ 0 。易知
∂ ∂ t = − c d d ξ , ∂ ∂ x = μ d d ξ , ∂ ∂ x 3 = μ 3 d 3 d ξ 3 , ∂ ∂ y = σ d d ξ , ∂ ∂ y 3 = σ 3 d 3 d ξ 3 . (6)
将(6)式代入方程(4)中得到下面这个非线性ODE
( a μ − c ) d u d ξ + b ( μ + η ) d u n d ξ + ( μ 3 + η 3 ) d 3 u n d ξ 3 = 0 (7)
为方便起见,在本文中我们令 ρ = μ 2 − μ η + η 2 ,且对任意的 μ ≠ 0 , η ≠ 0 都有 ρ > 0 ,下面不再赘述。
1) 当 n = 1 时,方程(7)变形为
d 3 u d ξ 3 + 1 ρ ( b + a μ − c μ + η ) d u d ξ = 0 , (8)
若 b + a μ − c μ + η ≥ 0 ,解得紧孤子解
u ( x , y , t ) = { ( λ 1 2 + λ 2 2 ) − 1 2 sin ( 1 ρ ( b + a μ − c μ + η ) ξ + θ ) + A , | φ | ≤ π 0 , 其 它
这里 φ = 1 ρ ( b + a μ − c μ + η ) ξ + θ , λ 1 , λ 2 , A 为任意常数,且当 λ 1 = λ 2 = 0 时, θ = 0 ;否则 θ 由式子 sin θ = λ 1 ( λ 1 2 + λ 2 2 ) − 1 2 , cos θ = λ 2 ( λ 1 2 + λ 2 2 ) − 1 2 确定。
若 b + a μ − c μ + η > 0 ,我们得到下列形式的孤立波相似解(下式中 λ 3 , λ 4 都是常数)
u ( x , y , t ) = λ 3 cosh ( λ 4 λ 3 ) ( − 1 ρ ( b + a μ − c μ + η ) ξ + A ) , 当 k 4 k 3 > 0 ; 或 u ( x , y , t ) = λ 3 sinh ( − λ 4 λ 3 ) ( − 1 ρ ( b + a μ − c μ + η ) ξ + A ) , 当 k 4 k 3 < 0.
2) 当 n ≥ 2 时,我们考虑下面两种情况:
情况I、 μ ≠ c a 时,对方程(7)积分一次并令积分常数为0得
( a μ − c ) u + b ( μ + η ) u n + ( μ 3 + η 3 ) d 2 u n d ξ 2 = 0. (9)
令 d u n d ξ = Z , d 2 u n d ξ 2 = Z d Z d u n ,代入式(9)中分离变量得
( μ 3 + η 3 ) Z d Z = n [ ( c − a μ ) u n − b ( μ + η ) u 2 n − 1 ] d u . (10)
进一步对方程 积分并令积分常数为0有
μ 3 + η 3 2 Z 2 = n [ c − a μ n + 1 u n + 1 − b ( μ + η ) 2 n u 2 n ] . (11)
由此有
n u n − 3 2 d u 2 n ( c − a μ ) ( n + 1 ) ( μ + η ) − b u n − 1 = ± d ξ ρ . (12)
令 V = 2 n ( c − a μ ) ( n + 1 ) ( μ + η ) − b u n − 1 ,
则当 b > 0 时,方程(12)变形为
d V 2 n ( c − a μ ) ( n + 1 ) ( μ + η ) − V 2 = ± n − 1 2 n b ρ d ξ . (13)
当 a μ − c μ + η ≠ 0 时,对式(13)积分解得如下紧孤子解
u ( x , y , t ) = { [ 2 n ( c − a μ ) b ( n + 1 ) ( μ + η ) cos 2 ( n − 1 2 n b ρ ξ ± r ) ] 1 n − 1 , | ϕ | ≤ π 2 0 , 其 它
这里 ϕ = n − 1 2 n b ρ ξ ± r 。在本文中,为方便起见我们用r表示任意积分常数。
类似地,在 b < 0 的情况下,我们解得孤立波相似解
u ( x , y , t ) = [ 2 n ( c − a μ ) b ( n + 1 ) ( μ + η ) cosh 2 ( n − 1 2 n − b ρ ξ ± i r ) ] 1 n − 1 , i = − 1 .
情况II、 μ = c a 。
方程(7)变形为线性方程
b d u n d ξ + ( μ 2 − μ η + η 2 ) d 3 u n d ξ 3 = 0. (14)
当 b > 0 时,解方程(14)得 u n = ( k 1 2 + k 2 2 ) 1 2 sin ( b ρ ξ + α ) + A ,其中 k 1 , k 2 为任意常数,且当 k 1 = k 2 = 0 时, α = 0 ,否则 α 满足 sin α = k 1 ( k 1 2 + k 2 2 ) − 1 2 , cos α = k 2 ( k 1 2 + k 2 2 ) − 1 2 。
由此得到紧孤子解
u ( x , y , t ) = { [ ( k 1 2 + k 2 2 ) 1 2 sin ( b ρ ξ + α ) + A ] 1 n , | ψ | ≤ π 0 , 其 它
这里 ψ = b ρ ξ + α 。
当 b < 0 ,我们有解
u ( x , y , t ) = [ k 3 exp ( − b ρ ξ ) + k 4 exp ( − − b ρ ξ ) + A ] 1 n , k 3 , k 4 为任意常数。
综上,我们整理得到下列形式孤立波相似解
u ( x , y , t ) = [ 2 k 3 cosh ( k 4 k 3 ) ( − b ρ ξ ) + A ] 1 n , 当 k 4 k 3 > 0 ; 或 u ( x , y , t ) = [ 2 k 3 cosh ( − k 4 k 3 ) ( − b ρ ξ ) + A ] 1 n , 当 k 4 k 3 < 0.
考虑具有负指数的方程(5)。按照前面所述方法,我们将行波变量 ξ = μ x + η y − c t 代入方程(5),使其变为下列ODE
( a μ − c ) d u d ξ + b ( μ + η ) d u − n d ξ + ( μ 3 + η 3 ) d 3 u − n d ξ 3 = 0 , (15)
其中 μ , η , c 为不为0的常数,且 μ , η , c 。
情况1、 μ ≠ c a 。
对方程(15)积分并令积分常数为0可得
( a μ − c ) u + b ( μ + η ) u − n + ( μ 3 + η 3 ) d 2 u − n d ξ 2 = 0. (16)
进一步变形有
− n u − n − 3 2 d u 2 n ( c − a μ ) ( n − 1 ) ( μ + η ) − b u − n − 1 = ± d ξ ρ , (17)
当 b > 0 时,解得周期解
u ( x , y , t ) = [ b ( n − 1 ) ( μ + η ) 2 n ( c − a μ ) sec 2 ( n + 1 2 n b ρ ξ ± r ) ] 1 n + 1 .
当 b < 0 时,解得孤立子解
u ( x , y , t ) = [ b ( n − 1 ) ( μ + η ) 2 n ( c − a μ ) sech 2 ( n + 1 2 n − b ρ ξ ± i r ) ] 1 n + 1 .
情况2、 μ = c a 。
此时我们考虑如下方程
b d u − n d ξ + ( μ 2 − μ η + η 2 ) d 3 u − n d ξ 3 = 0. (18)
当 b > 0 时,解之得周期解
u ( x , y , t ) = [ ( l 1 2 + l 2 2 ) 1 2 sin ( b ρ ξ + β ) + A ] − 1 n ,
这里 l 1 , l 2 , A 为常数, l 1 2 + l 2 2 ≠ 0 或 A ≠ 0 ,且 β 满足等式 sin β = l 1 ( l 1 2 + l 2 2 ) − 1 2 , cos β = l 2 ( l 1 2 + l 2 2 ) − 1 2 。
当 b < 0 时,解之得一般解
u ( x , y , t ) = [ l 3 exp ( − b ρ ξ ) + l 4 exp ( − − b ρ ξ ) + A ] − 1 n .
当 A = 0 时则有两个特殊的孤立子解
u ( x , y , t ) = [ 1 2 l 3 sech ( l 4 l 3 ) ( − b ρ ξ ) ] 1 n , 当 l 4 l 3 > 0 ; 或 u ( x , y , t ) = [ 1 2 l 3 csc h ( − l 4 l 3 ) ( − b ρ ξ ) ] 1 n , 当 l 4 l 3 < 0.
本文的主要目的是研究广义非线性具有正负n指数的KDV型方程(4)和(5)的精确解。我们运用变量代换法获得了包括紧孤子、孤立子、孤立波相似解和周期解在内的行波解。在下表1中,我们容易观察
指数 | n | μ | 系数 | 解的结构 |
---|---|---|---|---|
正指数 | n = 1 | b + a μ − c μ + η ≥ 0 | 紧孤立子解 | |
b + a μ − c μ + η < 0 | 孤立波相似解 | |||
n ≥ 2 | μ ≠ c a | b > 0 | 紧孤立子解 | |
b < 0 | 孤立波相似解 | |||
μ = c a | b > 0 | 紧孤立子解 | ||
b < 0 | 孤立波相似解 | |||
负指数 | n ≥ 2 | μ ≠ c a | b > 0 | 周期解 |
b < 0 | 孤立子解 | |||
μ = c a | b > 0 | 周期解 | ||
b < 0 | 孤立子解 |
表1. 广义色散方程解的物理结构
到不同的指数和系数决定了具有不同物理结构的解。
研究偏微分方程以获得其精确解进而研究非线性波现象的运动规律和特点是众多科学工作者一直以来的重要工作,值得指出的是,文中所使用的方法对于获得这些非线性偏微分方程的一些特殊解有一定的优势。
陈 静. 一类广义色散方程的行波解A Study on Exact Travelling Wave Solutions of Generalized Dispersive Equations[J]. 应用数学进展, 2018, 07(12): 1593-1599. https://doi.org/10.12677/AAM.2018.712186
https://doi.org/10.1016/S0096-3003(02)00610-0
https://doi.org/10.1103/PhysRevLett.70.564
https://doi.org/10.1016/S0378-4754(01)00291-9
https://doi.org/10.1016/j.camwa.2004.06.029
https://doi.org/10.1016/S0096-3003(02)00005-X
https://doi.org/10.1016/S0096-3003(03)00238-8
https://doi.org/10.1007/978-1-4899-1343-2_38
https://doi.org/10.1103/PhysRevE.60.6218
https://doi.org/10.1016/S0375-9601(00)00577-6
https://doi.org/10.1007/978-94-015-8289-6