Advances in Applied Mathematics
Vol. 09  No. 03 ( 2020 ), Article ID: 34460 , 8 pages
10.12677/AAM.2020.93033

Exact Solutions of the General Variable Coefficient Kdv-Burgers Equation

Ping Zhang1, Feng Luo2, Yuhuai Sun3

1The Engineering & Technical College of Chengdu University of Technology, Leshan Sichuan

2Leshan Vocational & Technical College, Leshan Sichuan

3Sichuan Normal University, Chengdu Sichuan

Received: Feb. 21st, 2020; accepted: Mar. 4th, 2020; published: Mar. 11th, 2020

ABSTRACT

In this paper, taking the generalized variable coefficient Kdv-Burgers equation as an example, the specific process of solving the exact solution of nonlinear partial differential equation by the expansion method of three auxiliary differential equations is introduced, and some three soliton solutions of the equation are obtained. These new soliton solutions include: hyperbolic function form solution, trigonometric function form solution, mixed action solution of hyperbolic function and exponential function, mixed action solution of trigonometric function and hyperbolic function, mixed action solution of trigonometric function and rational function, etc. Because of the arbitrariness of coefficients, the three auxiliary equations expansion method can construct more exact solutions of partial differential equations with variable coefficients.

Keywords:Expansion Method of Three Auxiliary Differential Equations, Generalized Variable Coefficient Kdv-Burgers Equation, Homogeneous Equilibrium Method, Exact Solution

广义变系数Kdv-Burgers方程的精确解

张萍1,罗缝2,孙峪怀3

1成都理工大学工程技术学院,四川 乐山

2乐山职业技术学院,四川 乐山

3四川师范大学,四川 成都

收稿日期:2020年2月21日;录用日期:2020年3月4日;发布日期:2020年3月11日

摘 要

本文以广义变系数Kdv-Burgers方程为例,介绍了三辅助微分方程展开法求非线性偏微分方程精确解的具体过程,并且由此得到了该方程的一些三孤子解,这些新的孤子解包含了:双曲函数形式解、三角函数形式解、双曲函数与指数函数混合作用解、三角函数与双曲函数混合作用解,三角函数与有理函数混合作用解等等。由于系数的任意性,使得三辅助方程展开法能够构造更多的变系数偏微分方程的精确解。

关键词 :三辅助微分方程展开法,广义变系数Kdv-Burgers方程,齐次平衡法,精确解

Copyright © 2020 by author(s) and Hans Publishers Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

1. 引言

随着科学技术的飞快发展,非线性偏微分方程的精确解广泛地应用于物理、数学、非线性光学、生物工程、通信、流体力学等领域,如非线性光学中对于光孤子的研究,使得暗孤子和亮孤子在光纤光学中被广泛应用;流体力学中,对于木星的红斑现象的解释等等。

人们在研究液体内含有气泡的流动、弹性管内液体的流动等问题的过程中,常常将这些问题的控制方程归结于变系数Kdv-Burgers方程 [1] [2],形如:

u t + a u u x + b u x x + u x x x = 0 (1)

其中 a b 0 ,a表示耗散项系数,b表示色散项系数。如果考虑更多的实际物理条件时,耗散项和色散项都与时间有关。本文主要研究系数与时间有关的广义变系数Kdv-Burgers方程 [3] [4] [5]:

u t + α ( t ) u u x + β ( t ) u x x + γ ( t ) u x x x = 0 (2)

这里的 α ( t ) β ( t ) γ ( t ) 为关于t的函数,其中 α ( t ) 表示耗散项系数, β ( t ) 表示色散项系数, γ ( t ) 表示对流项系数。

徐 [6] 使用新的三辅助微分方程展开法得到了变系数(2 + 1)维Nizhnik-Novikov-Vesselov方程的一些新的三孤子解,本文将以广义变系数Kdv-Burgers方程为例,来介绍新的三辅助微分方程展开法求非线性偏微分方程精确解的过程,从而获得它的一些新的精确解。

2. 三辅助方程展开法

步骤1:下列非线性偏微分方程

H ( u , u t , u x , u x t , u t t , u x x , ) = 0 (3)

假设方程(3)有如下形式的解:

u ( x , t ) = a 0 ( x , t ) + k = 1 n i + j + l = k a i j l ( x , t ) ϕ i ( ξ ) ( G ( θ ) G ( θ ) ) j ( F ( η ) F ( η ) ) l (4)

其中 a 0 ( x , t ) , a i j l ( x , t ) 是关于t的待定函数,通过齐次平衡法 [7] 平衡方程(3)可以确定正整数n的取值。

步骤2:辅助微分方程中的 ϕ ( ξ ) G ( θ ) F ( η ) 则满足下列三个常微分辅助方程:

ϕ ( ξ ) = A + B ϕ ( ξ ) + C ϕ 2 ( ξ ) (5)

G ( θ ) + λ G ( θ ) + μ G ( θ ) = 0 (6)

F ( η ) F ( η ) = Q F 2 ( η ) + R F ( η ) F ( η ) + P ( F ( η ) ) 2 (7)

其中 A , B , P , Q , R , λ , μ 为任意常数, C 0 , ξ = k 1 x + λ 1 ( t ) , θ = k 2 x + λ 2 ( t ) , η = k 3 x + λ 3 ( t ) k 1 , k 2 , k 3 为待定系数, λ 1 ( t ) , λ 2 ( t ) , λ 3 ( t ) 为关于t的待定函数。

步骤3:借助方程(5)、(6)和(7),将方程(4)代入到(3)中,此时,通过计算便可以得到一个关于 ϕ i ( ξ ) ( G ( θ ) G ( θ ) ) j ( F ( η ) F ( η ) ) l ( i , j , l = 0 , 1 , 2 , ) 的多项式,令这个多项式的每一项幂系数为零,得到一个关于 k 1 k 2 k 3 λ 1 ( t ) λ 2 ( t ) λ 3 ( t ) a 0 ( x , t ) a i j l ( x , t ) ( i , j , l = 0 , 1 , 2 , ) 的线性方程组,借助计算机软件Maple求解线性方程组,可以得到 k 1 k 2 k 3 λ 1 ( t ) λ 2 ( t ) λ 3 ( t ) a 0 ( x , t ) a i j l ( x , t ) ( i , j , l = 0 , 1 , 2 , ) 的取值。

步骤4:将所求得的 k 1 k 2 k 3 λ 1 ( t ) λ 2 ( t ) λ 3 ( t ) a 0 ( x , t ) a i j l ( x , t ) ( i , j , l = 0 , 1 , 2 , ) 的值和方程(5)、(6)、(7)的解代入到(4)中,从而得到方程(3)的解。

方程(5)、(6)、(7)的解,见表1~3。

Table 1. The solutions of Equation (5) with Δ = B 2 − 4 A C

表1. 方程(5)的解,其中 Δ = B 2 4 A C

Table 2. The G ′ ( θ ) G ( θ ) -solutions of Equation (6) with Δ 1 = λ 2 − 4 μ

表2. 方程(6)的形如 G ( θ ) G ( θ ) 解,其中 Δ 1 = λ 2 4 μ

Table 3. The F ′ ( η ) F ( η ) -solutions of Equation (7) with Δ 2 = Q ( P − 1 ) and Δ 3 = R 2 + 4 Q − 4 P Q

表3. 方程(7)的形如 F ( η ) F ( η ) 解,其中 Δ 2 = Q ( P 1 ) , Δ 3 = R 2 + 4 Q 4 P Q

3. 广义变系数Kdv-Burgers方程的解

下面用此方法来求解广义变系数Kdv-Burgers方程的三孤子解,通过齐次平衡法[7]平衡方程(2)中的 u u x u x x x ,有 n = 2 ,为了简化计算,假设方程(2)有如下形式的解:

u ( x , t ) = a 0 ( t ) + a 1 ( t ) ϕ ( ξ ) + a 2 ( t ) G ( θ ) G ( θ ) + a 3 ( t ) F ( η ) F ( η ) (8)

借助辅助微分方程(5)、(6)和(7),将方程(8)代入到方程(2)中,此时,通过计算可得到一组关于 ϕ i ( ξ ) ( G ( θ ) G ( θ ) ) j ( F ( η ) F ( η ) ) j ( i , j , l = 0 , 1 , 2 , 3 ) 的多项式,令每一项幂系数为零,得到一个关于 a 0 ( t ) , a 1 ( t ) , a 2 ( t ) , a 3 ( t ) , λ 1 ( t ) , λ 2 ( t ) , λ 3 ( t ) , k 1 , k 2 k 3 的线性方程组,借助Maple求解此方程组,可得到:

k 1 = k 1 , k 2 = k 2 , k 3 = k 3

α ( t ) = α ( t ) , β ( t ) = β ( t ) , γ ( t ) = γ (t)

a 0 ( t ) = a 0 ( t ) , a 1 ( t ) = a 3 ( t ) k 3 Q k 1 A , a 2 ( t ) = a 3 ( t ) k 3 k 2 , a 3 ( t ) = a 3 ( t ) (9)

λ 1 ( t ) = C 1 , λ 2 ( t ) = C 2 , λ 3 ( t ) = C 3

将式(9)代入到式(8)中,从而得到广义变系数Kdv-Burgers方程的如下形式的解:

情形1:当 Δ > 0 Δ 1 > 0 时,方程(2)有如下形式的解:

(1.1)当 R 0 , Δ 3 0 时:

u 1.1 ( x , t ) = a 0 ( t ) a 3 ( t ) k 3 Q k 1 A ( B + Δ 2 C tanh ( Δ 2 ξ ) ) + a 3 ( t ) k 3 k 2 ( λ 2 + Δ 1 2 c 1 sinh ( Δ 1 2 θ ) + c 2 cosh ( Δ 1 2 θ ) c 2 sinh ( Δ 1 2 θ ) + c 1 cosh ( Δ 1 2 θ ) ) + a 3 ( t ) ( R 2 ( 1 P ) + R Δ 3 2 ( 1 P ) c 3 exp ( Δ 3 2 η ) + c 4 exp ( Δ 3 2 η ) c 3 exp ( Δ 3 2 η ) c 4 exp ( Δ 3 2 η ) ) (10)

此时 u 1.1 ( x , t ) 为变系数Kdv-Burgers方程的双曲函数与指数函数混合形式解。

注:当 c 1 , c 2 为非零常数且满足如下关系时,可得到 u 1.1 ( x , t ) 有如下形式的解:

(1.1.1)当 c 2 2 > c 1 2 cosh ( β ) = c 2 c 2 2 c 1 2 sinh ( β ) = c 1 c 2 2 c 1 2 时:

u 1.1.1 ( x , t ) = a 0 ( t ) a 3 ( t ) k 3 Q k 1 A ( B + Δ 2 C tanh ( Δ 2 ξ ) ) + a 3 ( t ) k 3 k 2 ( λ 2 + Δ 1 2 coth ( Δ 1 2 θ + β ) ) + a 3 ( t ) ( R 2 ( 1 P ) + R Δ 3 2 ( 1 P ) c 3 exp ( Δ 3 2 η ) + c 4 exp ( Δ 3 2 η ) c 3 exp ( Δ 3 2 η ) c 4 exp ( Δ 3 2 η ) ) (11)

(1.1.2)当 c 2 2 < c 1 2 cosh ( β ) = c 1 c 1 2 c 2 2 sinh ( β ) = c 2 c 1 2 c 2 2 时:

u 1.1.2 ( x , t ) = a 0 ( t ) a 3 ( t ) k 3 Q k 1 A ( B + Δ 2 C tanh ( Δ 2 ξ ) ) + a 3 ( t ) k 3 k 2 ( λ 2 + Δ 1 2 tanh ( Δ 1 2 θ + β ) ) + a 3 ( t ) ( R 2 ( 1 P ) + R Δ 3 2 ( 1 P ) c 3 exp ( Δ 3 2 η ) + c 4 exp ( Δ 3 2 η ) c 3 exp ( Δ 3 2 η ) c 4 exp ( Δ 3 2 η ) ) (12)

(1.1.3)当 c 2 = 0 c 1 0 λ > 0 μ = 0 时:

u 1.1.3 ( x , t ) = a 0 ( t ) a 3 ( t ) k 3 Q k 1 A ( B + Δ 2 C tanh ( Δ 2 ξ ) ) + a 3 ( t ) k 3 k 2 ( λ 2 + λ 2 tanh ( λ 2 θ ) ) + a 3 ( t ) ( R 2 ( 1 P ) + R Δ 3 2 ( 1 P ) c 3 exp ( Δ 3 2 η ) + c 4 exp ( Δ 3 2 η ) c 3 exp ( Δ 3 2 η ) c 4 exp ( Δ 3 2 η ) ) (13)

(1.1.4)当 c 1 = 0 c 2 0 λ > 0 μ = 0 时:

u 1.1.4 ( x , t ) = a 0 ( t ) a 3 ( t ) k 3 Q k 1 A ( B + Δ 2 C tanh ( Δ 2 ξ ) ) + a 3 ( t ) k 3 k 2 ( λ 2 + λ 2 coth ( λ 2 θ ) ) + a 3 ( t ) ( R 2 ( 1 P ) + R Δ 3 2 ( 1 P ) c 3 exp ( Δ 3 2 η ) + c 4 exp ( Δ 3 2 η ) c 3 exp ( Δ 3 2 η ) c 4 exp ( Δ 3 2 η ) ) (14)

(1.2)当 R 0 , Δ 3 < 0 时:

u 1.2 ( x , t ) = a 0 ( t ) a 3 ( t ) k 3 Q k 1 A ( B + Δ 2 C tanh ( Δ 2 ξ ) ) + a 3 ( t ) k 3 k 2 ( λ 2 + Δ 1 2 c 1 sinh ( Δ 1 2 θ ) + c 2 cosh ( Δ 1 2 θ ) c 2 sinh ( Δ 1 2 θ ) + c 1 cosh ( Δ 1 2 θ ) ) + a 3 ( t ) ( R 2 ( 1 P ) + R Δ 3 2 ( 1 P ) i c 3 cos ( Δ 3 2 η ) c 4 sin ( Δ 3 2 η ) i c 3 sin ( Δ 3 2 η ) + c 4 cos ( Δ 3 2 η ) ) (15)

(1.3)当 R = 0 , Δ 2 0 时:

u 1.3 ( x , t ) = a 0 ( t ) a 3 ( t ) k 3 Q k 1 A ( B + Δ 2 C tanh ( Δ 2 ξ ) ) + a 3 ( t ) k 3 k 2 ( λ 2 + Δ 1 2 c 1 sinh ( Δ 1 2 θ ) + c 2 cosh ( Δ 1 2 θ ) c 2 sinh ( Δ 1 2 θ ) + c 1 cosh ( Δ 1 2 θ ) ) + a 3 ( t ) ( Δ 2 1 P c 3 cos ( Δ 2 η ) + c 4 sin ( Δ 2 η ) c 3 sin ( Δ 2 η ) c 4 cos ( Δ 2 η ) ) (16)

(1.4)当 R = 0 , Δ 2 < 0 时:

u 1.4 ( x , t ) = a 0 ( t ) a 3 ( t ) k 3 Q k 1 A ( B + Δ 2 C tanh ( Δ 2 ξ ) ) + a 3 ( t ) k 3 k 2 ( λ 2 + Δ 1 2 c 1 sinh ( Δ 1 2 θ ) + c 2 cosh ( Δ 1 2 θ ) c 2 sinh ( Δ 1 2 θ ) + c 1 cosh ( Δ 1 2 θ ) ) + a 3 ( t ) ( Δ 2 1 P i c 3 cosh ( Δ 2 η ) c 4 sinh ( Δ 2 η ) i c 3 sinh ( Δ 2 η ) c 4 cosh ( Δ 2 η ) ) (17)

情形2:当 Δ < 0 R = 0 , Δ 2 0 时,方程(2)有如下形式的解:

(2.1)当 Δ 1 < 0 时:

u 2.1 ( x , t ) = a 0 ( t ) a 3 ( t ) k 3 Q k 1 A ( B Δ 2 C cot ( Δ 2 ξ ) ) + a 3 ( t ) k 3 k 2 ( λ 2 + Δ 1 2 c 1 sin ( Δ 1 2 θ ) + c 2 cos ( Δ 1 2 θ ) c 2 sin ( Δ 1 2 θ ) + c 1 cos ( Δ 1 2 θ ) ) + a 3 ( t ) ( Δ 2 1 P c 3 cos ( Δ 2 η ) + c 4 sin ( Δ 2 η ) c 3 sin ( Δ 2 η ) c 4 cos ( Δ 2 η ) ) (18)

(2.2)当 Δ 1 > 0 时:

u 2.2 ( x , t ) = a 0 ( t ) a 3 ( t ) k 3 Q k 1 A ( B Δ 2 C cot ( Δ 2 ξ ) ) + a 3 ( t ) k 3 k 2 ( λ 2 + Δ 1 2 c 1 sinh ( Δ 1 2 θ ) + c 2 cosh ( Δ 1 2 θ ) c 2 sinh ( Δ 1 2 θ ) + c 1 cosh ( Δ 1 2 θ ) ) + a 3 ( t ) ( Δ 2 1 P c 3 cos ( Δ 2 η ) + c 4 sin ( Δ 2 η ) c 3 sin ( Δ 2 η ) c 4 cos ( Δ 2 η ) ) (19)

(2.3)当 Δ 1 = 0 时:

u 2.3 ( x , t ) = a 0 ( t ) a 3 ( t ) k 3 Q k 1 A ( B Δ 2 C cot ( Δ 2 ξ ) ) + a 3 ( t ) k 3 k 2 ( λ 2 + c 2 c 1 + c 2 θ ) + a 3 ( t ) ( Δ 2 1 P c 3 cos ( Δ 2 η ) + c 4 sin ( Δ 2 η ) c 3 sin ( Δ 2 η ) c 4 cos ( Δ 2 η ) ) (20)

表1~3还可以得到更多广义变系数Kdv-Burgers方程的精确解,这里就不再一一列出。

其中 ξ = k 1 x + C 1 θ = k 2 x + C 2 η = k 3 x + C 3 C 1 C 2 C 3 为任意常数。

4. 结论

本文以广义变系数Kdv-Burgers方程为例,介绍了三辅助方程展开法求解非线性偏微分方程的过程,并且得到了该方程的一些新的三孤子解。这些新解不仅比以前所得的解更加丰富,而且含有更多的函数类型,使得这些孤子新解在物理中有更重要的实际意义。在这些新解中,由于系数的任意性,使得三辅助方程展开法能够构造更多的变系数微分方程的精确解。

文章引用

张 萍,罗 缝,孙峪怀. 广义变系数Kdv-Burgers方程的精确解
Exact Solutions of the General Variable Coefficient Kdv-Burgers Equation[J]. 应用数学进展, 2020, 09(03): 277-284. https://doi.org/10.12677/AAM.2020.93033

参考文献

  1. 1. 王岗伟, 张颖元. 变系数Kdv-Burgers方程的精确解[J]. 聊城大学学报, 2011, 24(2): 9-12.

  2. 2. 王燕, 孙福伟. 变系数Kdv-Burgers方程的精确解及其Bäcklund变换[J]. 北方工业大学学报, 2009, 21(3):45-50.

  3. 3. 刘式达, 刘式适. 湍流的Kdv-Burgers模型[J]. 中国科学(A辑), 1991(9): 938-946.

  4. 4. 史秀珍, 斯仁道尔吉. 变系数Burgers方程与Kdv-Burgers方程的试探函数解[J]. 内蒙古大学学报(自然科学版), 2012, 43(1): 23-26.

  5. 5. Lü, H.L., Ming, Q.H. and Yu, J.Q. (2013) Exact Solutions of the Generalized Kdv-Burgers Equation. Journal of Qufu Normal University, 39, 33-37.

  6. 6. 徐兰兰, 陈怀堂. 变系数(2 + 1)维Nizhnik-Novikov-Vesselov方程的三孤子解[J]. 物理学报, 2013(9): 27-32.

  7. 7. 张国栋, 秦清锋. 齐次平衡法在微分方程中的应用[J]. 中国新技术新产品, 2010(21): 243-244.

期刊菜单