广义BBM-KdV方程的一个守恒C-N差分格式 A Conservative C-N Difference Scheme for the Generalized BBM-KdV Equation

doi:10.12677/PM.2021.114055

Pure Mathematics
Vol. 11 No. 04 ( 2021 ), Article ID: 41507 , 8 pages
10.12677/PM.2021.114055

广义BBM-KdV方程的一个守恒C-N差分格式

何丽，王希，胡劲松^*

●How to Cite this Article

西华大学理学院，四川成都

收稿日期：2021年3月2日；录用日期：2021年4月2日；发布日期：2021年4月12日

摘要

在进行非线性扩散波的研究时，BBM-KdV方程因能描述大量的物理现象如浅水波和离子波等而占有重要的地位，其数值研究少有涉及。本文研究了一类带有齐次边界条件的广义BBM-KdV方程的初边值问题，提出了一个具有二阶理论精度的两层非线性有限差分格式，合理模拟了问题本身的一个守恒量，并给出差分格式的先验估计，讨论其差分解的存在唯一性，并用离散泛函分析方法证明该格式的收敛性和无条件稳定性，最后通过数值模拟验证了该数值方法的可靠性。

关键词

广义BBM-KdV方程，差分格式，守恒，收敛性，稳定性

A Conservative C-N Difference Scheme for the Generalized BBM-KdV Equation

Li He, Xi Wang, Jinsong Hu^*

School of Science, Xihua University, Chengdu Sichuan

Received: Mar. 2^nd, 2021; accepted: Apr. 2^nd, 2021; published: Apr. 12^th, 2021

ABSTRACT

In the study of nonlinear diffusion waves, the BBM-KdV equation occupies an important position because it can describe a large number of physical phenomena such as shallow water waves and ion waves, and its numerical research is rarely involved. This paper studies the initial-boundary value problem of a generalized BBM-KdV equation with homogeneous boundary conditions, and proposes a two-level nonlinear finite difference scheme with second-order theoretical accuracy, which reasonably simulates a conserved quantity of the problem itself. A priori estimation of the difference scheme is given, and the existence and uniqueness of the difference decomposition is discussed. Discrete functional analysis is used to prove the convergence and unconditional stability of the scheme. Finally, the reliability of the numerical method is verified by numerical simulation.

Keywords:Generalized BBM-KdV Equation, Difference Scheme, Conservation, Convergence, Stability

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

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

1. 引言

本文考虑广义BBM-KdV方程

$u_{t} - u_{x x t} + β u_{x x x} + u_{x} + γ u^{p} u_{x} = 0$ ， $(x, t) \in (x_{L}, x_{R}) \times (0, T]$ ， (1)

的如下初边值问题：

$u (x, 0) = u_{0} (x)$ ， $x \in [x_{L}, x_{R}]$ ， (2)

$u (x_{L}, t) = u (x_{R}, t) = 0$ ， $u_{x} (x_{L}, t) = u_{x} (x_{R}, t) = 0$ ， $t \in [0, T]$ ， (3)

其中， $u_{0} (x)$ 是一个已知的初值函数， $β, γ$ 为正实数， $p \geq 1$ 为整数。问题(1)~(3)具有如下守恒律 [1]：

$E (t) = {‖ u ‖}_{L_{2}} + {‖ u_{x} ‖}_{L_{2}} = E (0)$ ， (4)

其中 $E (0)$ 为与初始条件有关的常数。

当 $p = 1$ 时，方程(1)即为通常的BBM-KdV方程 [1]

$u_{t} - u_{x x t} + β u_{x x x} + u_{x} + γ u u_{x} = 0$ 。 (5)

作为BBM方程 [2] [3] [4] 和KdV方程 [5] 推广，BBM-KdV方程(1)或(5)在进行非线性扩散波的研究时非常重要。BBM方程和KdV方程已引起了广泛地研究 [6] - [17]，而对于BBM-KdV方程的研究甚少，仅有文献 [18] [19] 通过数值模拟方法证实了BBM-KdV方程的解存在性，并讨论了其边界条件的物理意义，文献 [1] 进一步对BBM-KdV方程(5)给出了两个二阶精度的数值求解算法。

本文考虑更一般的情形，对广义BBM-KdV方程的初边值问题(1)~(3)进行数值求解研究，并进行数值模拟验证。

2. 数值格式和守恒性

剖分区域 $[x_{L}, x_{R}] \times [0, T]$ ，设 $τ$ 为时间步长， $t_{n} = n τ$ ， $0 \leq n \leq N$ ， $N = [\frac{T}{τ}]$ ； $x_{j} = x_{L} + j h$ ， $0 \leq j \leq J$ ， $h = \frac{x_{R} - x_{L}}{J}$ 为空间步长；记 $u_{j}^{n} = u (x_{j}, t_{n})$ ， $U_{j}^{n} \approx u (x_{j}, t_{n})$ ， $Z_{h}^{0} = {U = (U_{j}) | U_{- 1} = U_{0} = U_{J} = U_{J + 1} = 0, j = - 1, 0, 1, \dots, J, J + 1}$ 。规定C为与时间步长和空间步长均无关的常数，且 $C > 0$ 。并定义以下符号：

${(U_{j}^{n})}_{x} = \frac{U_{j + 1}^{n} - U_{j}^{n}}{h}$ ， ${(U_{j}^{n})}_{\bar{x}} = \frac{U_{j}^{n} - U_{j - 1}^{n}}{h}$ ， ${(U_{j}^{n})}_{\hat{x}} = \frac{U_{j + 1}^{n} - U_{j - 1}^{n}}{2 h}$ ， ${(U_{j}^{n})}_{t} = \frac{U_{j}^{n + 1} - U_{j}^{n}}{τ}$ ，

$U_{j}^{n + \frac{1}{2}} = \frac{U_{j}^{n + 1} + U_{j}^{n}}{2}$ ， $〈 U^{n}, V^{n} 〉 = h \sum_{j = 1}^{J - 1} U_{j}^{n} V_{j}^{n}$ ， ${‖ U^{n} ‖}^{2} = 〈 U^{n}, U^{n} 〉$ ， ${‖ U^{n} ‖}_{\infty} = \max_{1 \leq j \leq J - 1} | U_{j}^{n} |$ 。

由于 $u^{p} u_{x} = \frac{1}{p + 2} [u^{p} u_{x} + {(u^{p + 1})}_{x}]$ ，于是在数值离散时，对问题(1)~(3)提出如下两层C-N有限差分数值求解格式：

${(U_{j}^{n})}_{t} - {(U_{j}^{n})}_{x \bar{x} t} + β {(U_{j}^{n + \frac{1}{2}})}_{x \bar{x} \hat{x}} + {(U_{j}^{n + \frac{1}{2}})}_{\hat{x}} + φ (U_{j}^{n + \frac{1}{2}}) = 0$ ， $j = 1, \dots, J - 1$ ， $n = 1, \dots, N - 1$ ， (6)

$U_{j}^{0} = u_{0} (x_{j})$ ， $j = 0, \dots, J$ ， (7)

$U^{n} \in Z_{h}^{0}$ ， ${(U_{0}^{n})}_{\hat{x}} = {(U_{j}^{n})}_{\hat{x}} = 0$ ， $n = 0, \dots, N$ (8)

其中： $φ (U_{j}^{n + \frac{1}{2}}) = \frac{γ}{p + 2} {{(U_{j}^{n + \frac{1}{2}})}^{p} {(U_{j}^{n + \frac{1}{2}})}_{\hat{x}} + {[{(U_{j}^{n + \frac{1}{2}})}^{p + 1}]}_{\hat{x}}}$ 。

C-N差分格式(6)~(8)对不变量(4)有如下的数值模拟结果：

定理1 若定义离散能量 $E^{n} = {‖ U^{n} ‖}^{2} + {‖ U_{x}^{n} ‖}^{2}$ ，则C-N差分格式(6)~(8)关于 $E^{n}$ 是守恒，即

$E^{n} = E^{n - 1} = \dots = E^{0}$ ， (9)

其中， $n = 1, 2, \dots, N - 1$ 。

证明：以向量 $2 U^{n + \frac{1}{2}}$ 对(6)式取内积，由(8)式和离散分部求和公式 [20]，可以得到

${‖ U^{n} ‖}_{t}^{2} + {‖ U_{x}^{n} ‖}_{t}^{2} + 2 β 〈 U_{x \bar{x} \hat{x}}^{n + \frac{1}{2}}, U^{n + \frac{1}{2}} 〉 + 2 〈 U_{\hat{x}}^{n + \frac{1}{2}}, U^{n + \frac{1}{2}} 〉 + 2 〈 φ (U^{n + \frac{1}{2}}), U^{n + \frac{1}{2}} 〉 = 0$ ， (10)

又

$〈 U_{x \bar{x} \hat{x}}^{n + \frac{1}{2}}, U^{n + \frac{1}{2}} 〉 = 0$ ， $〈 U_{\hat{x}}^{n + \frac{1}{2}}, U^{n + \frac{1}{2}} 〉 = 0$ ， (11)

$\begin{matrix} 〈 φ (U^{n + \frac{1}{2}}), U^{n + \frac{1}{2}} 〉 = \frac{γ h}{p + 2} \sum_{j = 1}^{J - 1} {{(U_{j}^{n + \frac{1}{2}})}^{p} {(U_{j}^{n + \frac{1}{2}})}_{\hat{x}} + {[{(U_{j}^{n + \frac{1}{2}})}^{p + 1}]}_{\hat{x}}} U_{j}^{n + \frac{1}{2}} \\ = \frac{γ h}{p + 2} \sum_{j = 1}^{J - 1} {(U_{j}^{n + \frac{1}{2}})}^{p + 1} {(U_{j}^{n + \frac{1}{2}})}_{\hat{x}} + \frac{γ h}{p + 2} \sum_{j = 1}^{J - 1} {[{(U_{j}^{n + \frac{1}{2}})}^{p + 1}]}_{\hat{x}} U_{j}^{n + \frac{1}{2}} \\ = \frac{γ h}{p + 2} \sum_{j = 1}^{J - 1} {(U_{j}^{n + \frac{1}{2}})}^{p + 1} {(U_{j}^{n + \frac{1}{2}})}_{\hat{x}} - \frac{γ h}{p + 2} \sum_{j = 1}^{J - 1} {(U_{j}^{n + \frac{1}{2}})}^{p + 1} {(U_{j}^{n + \frac{1}{2}})}_{\hat{x}} \\ = 0 \end{matrix}$ ， (12)

于是将(11)式和(12)式带入(10)式，有

${‖ U^{n} ‖}_{t}^{2} + {‖ U_{x}^{n} ‖}_{t}^{2} = 0$

然后将上式两端乘以 $τ$ ，再对时间层n递推即可得(9)式成立。

3. 差分格式的可解性

以下我们用Brouwer不动点定理 [21] 和数学归纳法来证明C-N差分格式(6)~(8)的数值解是存在的，即

定理2 存在 $U^{n}$ 满足数值差分格式(6)~(8)，其中 $U^{n} \in Z_{h}^{0}$ ，且 $1 \leq n \leq N$ 。

证明：由(7)式知，显然 $U^{0}$ 是C-N差分格式(6)~(8)的解；现假设 $U^{n}$ 是，其中 $n \leq N - 1$ ；下面归纳证明存在 $U^{n + 1}$ 满足C-N差分格式(6)~(8)。

定义g是 $Z_{h}^{0}$ 上一个算子，且满足如下条件：

$g (V) = 2 (V - U^{n}) - 2 {(V - U^{n})}_{x \bar{x}} + τ V_{\hat{x}} + β τ V_{x \bar{x} \hat{x}} + τ φ (V)$ ， (13)

由(11)式和(12)式可知

$〈 V_{x \bar{x} \hat{x}}, V 〉 = 0$ ， $〈 V_{\hat{x}}, V 〉 = 0$ ， $〈 φ (V), V 〉 = 0$ ；

于是，以向量V对(13)式取内积，得

$\begin{matrix} 〈 g (V), V 〉 = 2 {‖ V ‖}^{2} - 2 〈 U^{n}, V 〉 + 2 {‖ V_{x} ‖}^{2} - 2 〈 U_{x}^{n}, V_{x} 〉 \\ \geq 2 {‖ V ‖}^{2} - 2 ‖ U^{n} ‖ \cdot ‖ V ‖ + 2 {‖ V_{x} ‖}^{2} - 2 ‖ U_{x}^{n} ‖ \cdot ‖ V_{x} ‖ \\ \geq 2 {‖ V ‖}^{2} - ({‖ U^{n} ‖}^{2} + {‖ V ‖}^{2}) + 2 {‖ V_{x} ‖}^{2} - ({‖ U_{x}^{n} ‖}^{2} + {‖ V_{x} ‖}^{2}) \\ \geq {‖ V ‖}^{2} - ({‖ U^{n} ‖}^{2} + {‖ U_{x}^{n} ‖}^{2}) + {‖ V_{x} ‖}^{2} \\ \geq {‖ V ‖}^{2} - ({‖ U^{n} ‖}^{2} + {‖ U_{x}^{n} ‖}^{2}) \end{matrix}$ ，

显然，如果取 ${‖ V ‖}^{2} = {‖ U^{n} ‖}^{2} + {‖ U_{x}^{n} ‖}^{2} + 1$ ，则有 $〈 g (V), V 〉 > 0$ 。从而由Brouwer不动点定理 [21] 可知，必存在 $V^{*} \in Z_{h}^{0}$ ，满足 $g (V^{*}) = 0$ ，即 $U^{n + 1} = 2 V^{*} - U^{n}$ 满足C-N差分格式(6)~(8)。从而由归纳假设知，C-N差分格式(6)~(8)的数值解是存在的。

4. 收敛性、稳定性与数值解的唯一性

我们将C-N差分格式(6)~(8)的截断误差定义如下：

$r_{j}^{n} = {(u_{j}^{n})}_{t} - {(u_{j}^{n})}_{x \bar{x} t} + β {(u_{j}^{n + \frac{1}{2}})}_{x \bar{x} \hat{x}} + {(u_{j}^{n + \frac{1}{2}})}_{\hat{x}} + φ (u_{j}^{n + \frac{1}{2}})$ ， (14)

且由Taylor公式可知， $r_{j}^{n} = O (τ^{2} + h^{2})$ 。

引理1初边值问题(1)~(3)的连续解满足如下估计：

${‖ u ‖}_{L_{2}} \leq C$ ， ${‖ u_{x} ‖}_{L_{2}} \leq C$ ， ${‖ u ‖}_{L_{\infty}} \leq C$ 。

证明：由不变量(4)式可得：

${‖ u ‖}_{L_{2}} \leq C$ ， ${‖ u_{x} ‖}_{L_{2}} \leq C$ ，

再根据Sobolev不等式即有， ${‖ u ‖}_{L_{\infty}} \leq C$ 。

定理3 设 $u_{0} \in H_{}^{2} [x_{L}, x_{R}]$ ，则C-N差分格式(6)~(8)的数值解满足以下估计：

$‖ U^{n} ‖ \leq C$ ， $‖ U_{x}^{n} ‖ \leq C$ ， ${‖ U^{n} ‖}_{\infty} \leq C$ ( $n = 1, 2, \dots, N$ )。

证明：由(9)式(即定理1)可得

$‖ U^{n} ‖ \leq C$ ， $‖ U_{x}^{n} ‖ \leq C$ ，

再根据离散Sobolev不等式 [20] 即有： ${‖ U^{n} ‖}_{\infty} \leq C$ 。

由定理3可知，C-N差分格式(6)~(8)的数值解 $U^{n}$ 以范数 ${‖ \cdot ‖}_{\infty}$ 关于初始值绝对稳定。

以下我们用能量方法来证明C-N差分格式(6)~(8)的数值解的收敛性，即

定理4 C-N差分格式(6)-(8)的数值解以范数 ${‖ \cdot ‖}_{\infty}$ 收敛到原初边值问题(1)~(3)的连续解，并且收敛阶为 $O (τ^{2} + h^{2})$ 。

证明：记 $e_{j}^{n} = u_{j}^{n} - U_{j}^{n}$ ，用(6)式去减(14)式，可得

$r_{j}^{n} = {(e_{j}^{n})}_{t} - {(e_{j}^{n})}_{x \bar{x} t} + β {(e_{j}^{n + \frac{1}{2}})}_{x \bar{x} \hat{x}} + {(e_{j}^{n + \frac{1}{2}})}_{\hat{x}} + φ (u_{j}^{n + \frac{1}{2}}) - φ (U_{j}^{n + \frac{1}{2}})$ ， (15)

以向量 $2 e^{n + \frac{1}{2}}$ 对(15)式两端取内积，并由离散分部求和公式 [20]，有

${‖ e^{n} ‖}_{t}^{2} + {‖ e_{x}^{n} ‖}_{t}^{2} = - 2 β 〈 e_{x \bar{x} \hat{x}}^{n + \frac{1}{2}}, e^{n + \frac{1}{2}} 〉 - 2 〈 e_{\hat{x}}^{n + \frac{1}{2}}, e^{n + \frac{1}{2}} 〉 - 2 〈 φ (u^{n + \frac{1}{2}}) - φ (U^{n + \frac{1}{2}}), e^{n + \frac{1}{2}} 〉 + 〈 r^{n}, 2 e^{n + \frac{1}{2}} 〉$ ， (16)

由(11)式，有

$〈 e_{\hat{x}}^{n + \frac{1}{2}}, e^{n + \frac{1}{2}} 〉 = 0$ ， $〈 e_{x \bar{x} \hat{x}}^{n + \frac{1}{2}}, e^{n + \frac{1}{2}} 〉 = 0$ ， (17)

根据引理1和定理3，利用离散Cauchy-Schwarz不等式，有

$\begin{array}{l} - 〈 φ (u^{n + \frac{1}{2}}) - φ (U^{n + \frac{1}{2}}), e^{n + \frac{1}{2}} 〉 \\ = \frac{- γ h}{p + 2} \sum_{j = 1}^{J - 1} [{(u_{j}^{n + \frac{1}{2}})}^{p} {(u_{j}^{n + \frac{1}{2}})}_{\hat{x}} - {(U_{j}^{n + \frac{1}{2}})}^{p} {(U_{j}^{n + \frac{1}{2}})}_{\hat{x}}] e_{j}^{n + \frac{1}{2}} \\ - \frac{γ h}{p + 2} \sum_{j = 1}^{J - 1} {{[{(u_{j}^{n + \frac{1}{2}})}^{p + 1}]}_{\hat{x}} - {[{(U_{j}^{n + \frac{1}{2}})}^{p + 1}]}_{\hat{x}}} e_{j}^{n + \frac{1}{2}} \\ = \frac{- γ h}{p + 2} \sum_{j = 1}^{J - 1} [{(u_{j}^{n + \frac{1}{2}})}^{p} {(e_{j}^{n + \frac{1}{2}})}_{\hat{x}} + e_{j}^{n + \frac{1}{2}} {(U_{j}^{n + \frac{1}{2}})}_{\hat{x}} \sum_{k = 0}^{p - 1} {(u_{j}^{n + \frac{1}{2}})}^{p - 1 - k} {(U_{j}^{n + \frac{1}{2}})}^{k}] e_{j}^{n + \frac{1}{2}} \\ + \frac{γ h}{p + 2} \sum_{j = 1}^{J - 1} [e_{j}^{n + \frac{1}{2}} \sum_{k = 0}^{p} {(u_{j}^{n + \frac{1}{2}})}^{p - k} {(U_{j}^{n + \frac{1}{2}})}^{k}] {(e_{j}^{n + \frac{1}{2}})}_{\hat{x}} \\ \leq C ({‖ e_{^{x}}^{n + 1} ‖}^{2} + {‖ e_{^{x}}^{n} ‖}^{2} + {‖ e^{n + 1} ‖}^{2} + {‖ e^{n} ‖}^{2}) \end{array}$ ， (18)

$〈 r^{n}, 2 e^{n + \frac{1}{2}} 〉 \leq {‖ r^{n} ‖}^{2} + \frac{1}{2} ({‖ e^{n + 1} ‖}^{2} + {‖ e^{n} ‖}^{2})$ ， (19)

将(17)~(19)三式一起代入(16)式，令 $B^{n} = {‖ e^{n} ‖}^{2} + {‖ e_{x}^{n} ‖}^{2}$ ，然后整理得

$B^{n + 1} - B^{n} \leq τ {‖ r^{n} ‖}^{2} + C τ (B^{n + 1} + B^{n})$ ， (20)

将(20)式从0到 $N - 1$ 递推并求和，即得

$B^{N} \leq B^{0} + C τ \sum_{n = 0}^{N - 1} {‖ r^{n} ‖}^{2} + C τ \sum_{n = 0}^{N} B^{n}$ ，

又 $τ \sum_{n = 0}^{N - 1} {‖ r^{n} ‖}^{2} \leq N τ \max_{0 \leq n \leq N - 1} {‖ r^{n} ‖}^{2} \leq T \cdot O {(τ^{2} + h^{2})}^{2}$ ， $B^{0} = O {(τ^{2} + h^{2})}^{2}$ ，于是

$B^{N} \leq O {(τ^{2} + h^{2})}^{2} + C τ \sum_{n = 0}^{N - 1} B^{n}$ ，

从而由离散Gronwall不等式 [20]，有

$‖ e^{N} ‖ \leq O (τ^{2} + h^{2})$ ， $‖ e_{x}^{N} ‖ \leq O (τ^{2} + h^{2})$ ，

最后再根据由离散Sobolev不等式 [20]，即得

${‖ e^{N} ‖}_{\infty} \leq O (τ^{2} + h^{2})$ 。

用类似定理的证明方法，有：

定理5 C-N差分格式(6)~(8)的数值解唯一。

5. 算例验证

为了验证本文数值方法的可靠性，考虑 $p = 3$ 、 $p = 5$ 两种情形分别进行数值实验。固定参数 $β = 1$ ， $γ = 0.25$ ，其中，当 $p = 3$ 时，广义BBM-KdV方程(1)的孤波解为：

$u (x, t) = 10^{\frac{1}{3}} {sech}^{\frac{2}{3}} (\frac{1}{2} x - \frac{5}{8} t)$ ，

当 $p = 5$ 时，广义BBM-KdV方程(1)的孤波解为：

$u (x, t) = 7^{\frac{1}{5}} {sech}^{\frac{2}{5}} (\frac{1}{2} x - \frac{13}{24} t)$ 。

取初值函数 $u_{0} (x) = u (x, 0)$ 进行数值求解计算，固定 $x_{L} = - 40$ ， $x_{R} = 60$ ， $T = 20$ 。对时间步长 $τ$ 和空间步长h的取不同值时，将数值解和孤波解在不同时刻误差的 $l_{\infty}$ 模列于表1；将C-N差分格式对不变量(4)的数值模拟能量 $E^{n}$ 在不同时刻的数据列于表2。

Table 1. l ∞ -errors of C-N differential scheme at different moments

表1. C-N差分格式在不同时刻的 $l_{\infty}$ -误差

Table 2. Partial data of C-N difference scheme for analog E n for invariant (4)

表2. C-N差分格式对不变量(4)的模拟 $E^{n}$ 的部分数据

6. 结论

本文对一类带有齐次边界条件的广义BBM-KdV方程的初边值问题(1)~(3)进行了数值方法研究，提出了一个两层非线性差分格式(6)~(8)，该格式是无条件稳定的。从表1可以看出，本文的格式明显具有二阶精度；从表2可以看出，数值式格式对原问题的守恒性质(4)也进行了合理有效地模拟。所以本文数值求解方法是可靠的。

基金项目

国家自然科学基金项目(11701481)；四川应用基础研究项目(2019JY0387)。

文章引用

何丽,王希,胡劲松. 广义BBM-KdV方程的一个守恒C-N差分格式
A Conservative C-N Difference Scheme for the Generalized BBM-KdV Equation[J]. 理论数学, 2021, 11(04): 428-435. https://doi.org/10.12677/PM.2021.114055

参考文献

1. Rouatbi, A. and Omrani, K. (2017) Two Conservative Difference Schemes for a Model of Nonlinear Dispersive Equa-tions. Chaos, Solitons & Fractals, 104, 516-530. https://doi.org/10.1016/j.chaos.2017.09.006

2. Peregrine, D.H. (1966) Calculations of the Development of an Undular Bore. Journal of Fluid Mechanics, 25, 321-330.

3. Peregrine, D.H. (1967) Long Waves on Beach. Journal of Fluid Mechanics, 27, 815-827. https://doi.org/10.1017/S0022112067002605

4. Benjamin, T.B., Bona, J.L. and Mahony, J.J. (1972) Model Equations for Long Waves in Nonlinear Dispersive Systems. Philosophical Transactions of the Royal Society London, Series A, 272, 47-48. https://doi.org/10.1098/rsta.1972.0032

5. Kortewag, D.J. and De Vries, G. (1985) On the Change of Form of Long Waves Advancing in a Rectangular Canal, and on a New Type of Long Stationary Waves. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, 39, 422-443. https://doi.org/10.1080/14786449508620739

6. Achouri, T., Khiari, N. and Omrani, K. (2006) On the Conver-gence of Difference Schemes for the Benjamin-Bona-Mahony (BBM) Equation. Applied Mathematics and Computation, 182, 999-1005. https://doi.org/10.1016/j.amc.2006.04.069

7. Omrani, K. and Ayadi, M. (2008) Finite Difference Discretization of the Benjamin-Bona-Mahony-Burgers Equation. Numerical Methods for Partial Differential Equations, 24, 239-248. https://doi.org/10.1002/num.20256

8. Khalifa, A.K., Raslan, K.R. and Alzubaidi, H.M. (2007) A Finite Differ-ence Scheme for the MRLW and Solitary Waves Interactions. Applied Mathematics and Computation, 189, 346-354. https://doi.org/10.1016/j.amc.2006.11.104

9. Omrani, K. (2006) The Convergence of Fully Discrete Galerkin Approximations for the Benjamin-Bona-Mahony (BBM) Equation. Applied Mathematics and Computation, 180, 614-621. https://doi.org/10.1016/j.amc.2005.12.046

10. Achouri, T., Ayadi, M. and Omrani, K. (2009) A Fully Galerkin Method for the Damped Generalized Regularized Long-Wave (DGRLW) Equation. Numerical Methods for Partial Differential Equations, 25, 668-684. https://doi.org/10.1002/num.20367

11. Kadri, T., Khiari, N., Abidi, F. and Omrani, K. (2008) Methods for the Numerical Solution of the Benjamin-Bona-Mahony-Burgers Equation. Numerical Methods for Partial Differential Equations, 24, 1501-1516. https://doi.org/10.1002/num.20330

12. Achouri, T. and Omrani, K. (2010) Application of the Homotopy Per-turbation Method to the Modified Regularized Long-Wave Equation. Numerical Methods for Partial Differential Equations, 26, 399-411. https://doi.org/10.1002/num.20441

13. Abbasbandy, S. and Shirzadi, A. (2010) The First Integral Method for Modified Benjamin-Bona-Mahony Equation. Communications in Nonlinear Science and Numerical Simulation, 15, 1759-1764. https://doi.org/10.1016/j.cnsns.2009.08.003

14. Dehghan, M.., Abbaszadeh, M. and Mohebbi, A. (2014) The Numerical Solution of Nonlinear High Dimensional Generalized Benjamin-Bona-Mahony-Burgers Equation via the Meshless Method of Radial Basis Functions. Computers & Mathematics with Applications, 68, 212-237. https://doi.org/10.1016/j.camwa.2014.05.019

15. Saka, B., Şahin, A. and Dağ, İ. (2011) B-Spline Collocation Algorithms for Numerical Solution of the RLW Equation. Numerical Methods for Partial Differential Equations, 27, 581-607. https://doi.org/10.1002/num.20540

16. Mei, L. and Chen, Y. (2012) Numerical Solutions of RLW Equation Using Galerkin Method with Extrapolation Techniques. Computer Physics Communications, 183, 1609-1616. https://doi.org/10.1016/j.cpc.2012.02.029

17. Mohammadi, M. and Mokhtari, R. (2011) Solving the Generalized Regularized Long Wave Equation on the Basis of a Reproducing Kernel Space. Journal of Computational and Applied Mathematics, 235, 4003-4014. https://doi.org/10.1016/j.cam.2011.02.012

18. Dutykh, D. and Pelinovsky, E. (2014) Numerical Simulation of a Solitonic Gas in KDV and KDV-BBM Equations. Physics Letters A, 378, 3102-3110. https://doi.org/10.1016/j.physleta.2014.09.008

19. Asokan, R. and Vinodh, D. (2018) Soliton and Exact Solutions for the KdV-BBM Type Equations by Tanh-Coth and Transformed Rational Function Methods. International Journal of Applied and Computational Mathematics, 4, Article No. 100. https://doi.org/10.1007/s40819-018-0533-7

20. Chou, Y.L. (1991) Application of Discrete Functional Analysis to the Finite Difference Methods. Pergamon Press, Elmsford, 260 p.

21. Browder, F.E. (1965) Existence and Uniqueness Theorems for Solutions of Nonlinear Boundary Value Problems. Proceedings of Symposia in Applied Mathematics, 17, 24-49. https://doi.org/10.1090/psapm/017/0197933

期刊菜单