白 雪，李 宏，方志朝

内蒙古大学数学科学学院，内蒙古 呼和浩特

收稿日期：2018年3月1日；录用日期：2018年3月20日；发布日期：2018年3月27日

本文研究了正则长波Burgers方程的混合有限体积元方法。引入迁移算子把试探函数空间映射为检验函数空间，构造了半离散和线性向后Euler全离散混合有限体积元格式，证明了两种离散格式解的存在唯一性，并得到了最优阶误差估计。最后，给出数值算例来验证理论分析结果和数值格式的有效性。

关键词 :混合有限体积元方法，正则长波Burgers方程，最优阶误差估计

Copyright © 2018 by authors and beplay安卓登录

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

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

本文考虑如下正则长波Burgers方程

{ ( a )     u t + σ u x + δ u u x − α u x x − β u x x t = 0 ,   ( x , t ) ∈ Ω × Γ , ( b )     u ( a , t ) = u ( b , t ) = 0 ,   t ∈ Γ ¯ , ( c )     u ( x , 0 ) = u 0 ( x ) ,   x ∈ Ω , (1)

其中 Ω = ( a , b ) ， Γ = ( 0 , T ] ， 0 < T < ∞ ，系数 α , δ , β 和s是给定的正常数。a和b分别为耗散系数和色散系数。 u 0 ( x ) 为给定的初始函数。

RLW-Burgers方程(1)又称为Benjamin-Bona-Mahony-Burgers方程，是带有耗散项 − α u x x 的正则长波方程，这类数学物理方程在非线性波问题中有着十分重要的地位。因此，求解此类问题的方法受到广泛的重视，文献 [ 1 ] 利用未知函数变换的思想，得到了RLW-Burgers方程的解析解，文献 [ 2 ] 运用一类含参数的G展开法得到了多种函数形式的显式行波解。这些解析解是在限定了初边值条件下给出的，因此研究该方程的数值计算方法是很有意义的。文献 [ 3 ] 研究了RLW-Burgers方程的有限差分法，给出了收敛性和稳定性分析。文献 [ 4 ] 构造了RLW-Burgers方程的H¹-Galerkin混合有限元格式，给出了半离散和全离散格式的最优阶误差估计，并给出数值实验算例。文献 [ 5 ] 运用径向基函数法数值求解了高维广义Benjamin-Bona-Mahony-Burgers方程。

混合有限体积元方法又称为混合控制体积方法，最早是由Russel在文献 [ 6 ] 中求解线性椭圆问题时提出，随后文献 [ 7 ] [ 8 ] 通过数值实验验证该方法的可行性。Chou等人在文献 [ 9 ] [ 10 ] 中分别在三角形网格和矩形网格剖分下，使用最低阶Raviart-Thomas有限元空间作为试探函数空间，构造了椭圆方程的混合有限体积元格式并给出了最优阶误差估计。该方法不仅可以利用有限体积元方法的简便性同时求解多个不同的物理量，而且还可以保持某些物理量的局部守恒性，因此该方法自提出以来得到了很大的关注和发展，已经被用来求解二阶拟线性椭圆方程 [ 11 ] 、正则长波方程 [ 12 ] 、抛物方程 [ 13 ] [ 14 ] 、阻尼sine-Gordon方程 [ 15 ] 等数学物理方程。本文研究了RLW-Burgers方程(1)的混合有限体积元方法，构造了半离散和线性向后Euler全离散混合有限体积元格式，给出了理论分析并通过数值算例验证了方法的有效性。

本文余下部分安排如下：第2节，给出混合有限元形式，并引入迁移算子构造了半离散混合有限体积元格式；第3节，给出了迁移算子的一些性质，并给出半离散格式解的存在唯一性和稳定性以及误差估计；第4节，构造了线性向后Euler全离散格式，给出了全离散格式解的存在唯一性和最优阶误差估计；第5节，给出了数值算例验证该方法的有效性。本文理论分析时所用的C都是与空间网格参数h和时间离散步长 Δ t 无关的正常数。

引入流量函数 f ( u ) = σ u + δ 2 f ( u ) − α u x − β u x t ，其中 f ( u ) = u 2 ，于是可以把原方程(1)改写为下面的形式

{ ( a )     u t + p x = 0 ,   ( x , t ) ∈ Ω × Γ , ( b )     p = σ u + δ 2 f ( u ) − α u x − β u x t ,   ( x , t ) ∈ Ω × Γ , ( c )     u ( a , t ) = u ( b , t ) = 0 ,   t ∈ Γ ¯ , ( d )     u ( x , 0 ) = u 0 ( x ) ,   x ∈ Ω . (2)

则相应的混合变分形式为：求 { u ( t ) , p ( t ) } ∈ H 0 1 ( Ω ) × L 2 ( Ω ) ，使得

{ ( a )     ( u t , v ) − ( v x , p ) = 0 ,   ∀ v ∈ H 0 1 ( Ω ) , ( b )     ( p , q ) − σ ( u , q ) + α ( u x , q ) + β ( u x t , q ) = δ 2 ( f ( u ) , q ) ,   ∀ q ∈ L 2 ( Ω ) , ( c )     u ( x , 0 ) = u 0 ( x ) ,   x ∈ Ω . (3)

现在对区域 Ω ¯ = [ a , b ] 使用以下节点做原始网格剖分： a = x 0 < x 1 < x 2 < ⋯ < x N = b 。对应的对偶网格剖分为： a = x 0 < x 1 / 2 < x 3 / 2 < ⋯ < x N − 1 / 2 < x N = b ，其中 x i + 1 / 2 = ( x i + x i + 1 ) / 2 , ( 0 ≤ i ≤ N − 1 ) 。则原始网格剖分 T h = { A i = [ x i , x i + 1 ] ; i = 0 , 1 , ⋯ , N − 1 } 。记 h i = x i + 1 − x i ， h = max 0 ≤ i ≤ N − 1 h i 。假设 T h 为拟一致剖分，即存在正常数a使得 h i ≥ α h , ( i = 0 , 1 , ⋯ , N − 1 ) 。记 A 0 * = [ x 0 , x 1 / 2 ] ， A N * = [ x N − 1 / 2 , x N ] ， A i * = [ x i − 1 / 2 , x i + 1 / 2 ] , ( i = 1 , 2 , ⋯ , N − 1 ) ，则对偶剖分为 T h * = { A i * : i = 0 , 1 , ⋯ , N } 。基于剖分 T h ，选取有限元空间 H 0 h × L h 作为试探函数空间，其中

H 0 h = { ω h ∈ H 0 1 ( Ω ) : ω h ∈ P 1 ( A ) , ∀ A ∈ T h } ,   L h = { v h ∈ L 2 ( Ω ) : v h | A ∈ P 0 ( A ) , ∀ A ∈ T h } .

记 { φ i : i = 1 , 2 , ⋯ , N − 1 } 和 { χ A i : i = 0 , 1 , ⋯ , N − 1 } 分别为有限元空间 H 0 h 和 L h 的基，其中 ϕ i 是分段线性多项式(见文献 [ 12 ] )， χ A i 是集合 A i 的特征函数。

定义迁移算子 γ h : H 0 h → L 2 ( Ω ) 如下

γ h ω h = ∑ i = 1 N − 1 ω h ( x i ) χ A i * ,     ∀ ω h ∈ H 0 h .

则有 γ h φ i = χ A i * ( i = 1 , 2 , ⋯ , N − 1 ) 。选取检验函数空间 Y h = R ( γ h ) 为 γ h 的值域。系统(1)中(a)和(b)方程分别在对偶剖分和原始剖分上进行积分，得到

{ ( a )     ∫ A i * u t d x + p ( x i + 1 / 2 , t ) − p ( x i − 1 / 2 , t ) = 0 ,     t ∈ Γ ,     ( 1 ≤ i ≤ N − 1 ) , ( b )     ∫ A i ( p − σ u + α u x + β u x t ) d t = δ 2 ∫ A i f ( u ) d x ,     t ∈ Γ ,     ( 0 ≤ i ≤ N − 1 ) . (4)

利用迁移算子 γ h ，可以把(4)改写为下列的形式：

{ ( a )     ( u t + p x , γ h v h ) = 0 ,     ∀ v h ∈ H 0 h , ( b )     ( p − σ u + α u x + β u x t , q h ) = δ 2 ( f ( u ) , q h ) ,     ∀ q h ∈ L h . (5)

注意到

( p x , γ h v h ) = ∑ i = 1 N − 1 ∫ A i * ( p x v h ( x i ) ) d x = ∑ i = 1 N − 1 v h ( x i ) ( p ( x i + 1 / 2 ) − p ( x i − 1 / 2 ) ) d x = b ( γ h v h , p ) . (6)

通过计算可知 b ( γ h v h , q h ) = − ( v h x , q h ) , ∀ v h ∈ H 0 h , ∀ q h ∈ L h ，利用形式(5)可以得到半离散混合有限体积元格式：求 { u h ( t ) , p h ( t ) } ∈ H 0 h × L h ，使得

{ ( a )     ( u h t , σ h v h ) − ( v h x , p h ) = 0 ,     ∀ v h ∈ H 0 h , ( b )     ( p h , q h ) − σ ( u h , q h ) + α ( u h x , q h ) + β ( u h x t , q h ) = δ 2 ( f ( u h ) , q h ) ,     ∀ q h ∈ L h , ( c )     u h ( 0 ) = π h u 0 , (7)

其中 π h 是椭圆投影算子，具体定义见第3节。

为了进行误差估计，首先给出迁移算子 γ h 的性质。

引理3.1 [ 12 ] ：迁移算子 γ h 是有界的，且满足

‖ γ h ω h ‖ ≤ 3 ‖ ω h ‖ ,     ∀ ω h ∈ H 0 h . (8)

引理3.2 [ 12 ] ：迁移算子 γ h 满足对称性，即有

( ω h , γ h v h ) = ( γ h ω h , v h ) ,     ∀ v h , ω h ∈ H 0 h . (9)

引理3.3 [ 12 ] ：迁移算子 γ h 满足下述正定性

( γ h ω h , ω h ) ≥ 3 4 ‖ ω h ‖ 2 ,     ∀ ω h ∈ H 0 h . (10)

引理3.4 [ 12 ] ：设I是恒等算子，则迁移算子 γ h 满足下述性质

‖ ( I − γ h ) ω h ‖ ≤ 12 12 h | ω h | 1 ,     ∀ ω h ∈ H 0 h . (11)

| ( v h , ( I − γ h ) ω h ) | ≤ 12 12 h | p h | 1 ‖ ω h ‖ ,     ∀ p h , ω h ∈ H 0 h . (12)

下面引入椭圆投影算子 π h : H 0 1 ( Ω ) → H 0 h ，满足

( ( ω − π h ω ) x , v h ) = 0 ,   ∀ v h ∈ L h , 其中 ω ∈ H 0 1 ( Ω ) . (13)

同时引入 L 2 正交投影算子 R h : L 2 ( Ω ) → L h ，满足

( φ − R h φ , v h ) = 0 ,     ∀ v h ∈ L h , 其中 φ ∈ L 2 ( Ω ) . (14)

由文献 [ 16 ] [ 17 ] 可知，对 i = 0 , 1 ，下列估计式成立

‖ ∂ i ω ∂ t i − ∂ i π h ω ∂ t i ‖ ≤ C h ‖ ∂ i ω ∂ t i ‖ 1 , 其中 ∂ i ω ∂ t i ∈ H 0 1 ( Ω ) , (15)

‖ ∂ i ω ∂ t i − ∂ i π h ω ∂ t i ‖ 1 ≤ C h ‖ ∂ i ω ∂ t i ‖ 2 , 其中 ∂ i ω ∂ t i ∈ H 2 ( Ω ) ∩ H 0 1 ( Ω ) , (16)

‖ ∂ i φ ∂ t i − ∂ i R h φ ∂ t i ‖ ≤ C h ‖ ∂ i φ ∂ t i ‖ 1 , 其中 ∂ i φ ∂ t i ∈ H 1 ( Ω ) . (17)

利用算子 γ h ，可以得到下述引理。

引理3.5 [ 12 ] ：存在与h无关的正常数C使得

‖ v − γ h π h v ‖ ≤ C h ‖ v ‖ 1 ,     ∀ v ∈ H 0 1 ( Ω ) , (18)

| ( v , ( I − γ h ) ω h ) | ≤ C h ‖ v ‖ 1 ‖ ω h ‖ ,     ∀ v ∈ H 0 1 ( Ω ) ,     ∀ ω h ∈ H 0 h . (19)

为了对半离散格式进行误差估计，将误差写为

u − u h = u − π h u + π h u − u h = ξ + η ,

p − p h = p − R h p + R h p − p h = ρ + ε .

格式(3)和格式(7)作差，利用投影算子可得误差方程为

{ ( a )     ( η t , γ h v h ) − ( v h x , ε ) = − ( ξ t , γ h v h ) − ( u t , ( I − γ h ) v h ) ,     ∀ v h ∈ H 0 h , ( b )     α ( η x , q h ) + β ( η x t , q h ) + ( ε , q h ) − σ ( η , q h ) = ( ξ , q h ) − δ 2 ( f ( u h ) − f ( u ) , q h ) ,   ∀ q h ∈ L h , ( c )     η ( 0 ) = π h u ( 0 ) − u h ( 0 ) = 0. (20)

类似于文献 [ 12 ] 中的证明，可以得到下述半离散格式解的存在唯一性和稳定性。

定理3.1：半离散格式(7)存在唯一解 { u h , p h } ，且存在与h和t无关的正常数C，使得

‖ u h ( t ) ‖ 1 ≤ C ‖ u 0 ‖ 1 ,     ‖ u h ( t ) ‖ ∞ ≤ C ‖ u 0 ‖ 1 ,

‖ u h t ( t ) ‖ 1 + ‖ p h ( t ) ‖ ≤ C ( ‖ u 0 ‖ 1 + ‖ u 0 ‖ 1 2 + ‖ u 0 ‖ 1 3 ) .

下面给出半离散格式(7)的误差估计。

定理3.2：设 { u , p } 是混合变分形式(3)的解， { u h , p h } 是半离散格式(7)的解，则存在与h和t无关的正常数C，使得

‖ u ( t ) − u h ( t ) ‖ 1 ≤ C h ( ‖ u ( t ) ‖ 2 + ( ∫ 0 t ( ‖ u ‖ 1 2 + ‖ u t ‖ 1 2 ) d s ) 1 2 ) ,

‖ u t ( t ) − u h t ( t ) ‖ 1 ≤ C h ( ‖ u t ( t ) ‖ 2 + ‖ u ( t ) ‖ 1 + ( ∫ 0 t ( ‖ u ‖ 1 2 + ‖ u t ‖ 1 2 ) d s ) 1 2 ) ,

‖ p ( t ) − p h ( t ) ‖ ≤ C h ( ‖ p ( t ) ‖ 1 + ‖ u ( t ) ‖ 1 + ‖ u t ‖ 1 + ( ∫ 0 t ( ‖ u ‖ 1 2 + ‖ u t ‖ 1 2 ) d s ) 1 2 ) .

证明：在误差方程(20)中取 v h = η , q h = η x ，则有

α ( η x , η x ) + β ( η x t , η x ) + ( η t , γ h η ) − σ ( η , η x ) = − ( ξ t , γ h η ) − ( u t , ( I − γ h ) η ) + ( ξ , η x ) − δ 2 ( f ( u h ) − f ( u ) , η x ) . (21)

则由引理3.5及引理3.1，有

α ‖ η x ‖ 2 + 1 2 d d t ( β ‖ η x ‖ 2 + ( η , γ h η ) ) − σ ( η , η x ) = − ( ξ t , γ h η ) − ( u t , ( I − γ h ) η ) + ( ξ , η x ) − δ 2 ( f ( u h ) − f ( u ) , η x ) . (22)

易知 ‖ f ( u ) − f ( u h ) ‖ ≤ C ( ‖ ξ ‖ + ‖ η ‖ ) ,   ( η , η x ) = 0 。则(22)式可变为

α ‖ η x ‖ 2 + 1 2 d d t ( β ‖ η x ‖ 2 + ( η , γ h η ) ) ≤ α 2 ‖ η x ‖ 2 + C ( ‖ ξ ‖ 2 + ‖ ξ t ‖ 2 + h 2 ‖ u t ‖ 1 2 + ‖ η ‖ 2 ) . (23)

对(23)两边关于t积分，可得

α ∫ 0 t ‖ η x ‖ 2 d t + β ‖ η x ‖ 2 + 3 4 ‖ η ‖ 2 ≤ C ∫ 0 t ( ‖ ξ ‖ 2 + ‖ ξ t ‖ 2 + h 2 ‖ u t ‖ 1 2 + ‖ η ‖ 2 ) d s   . (24)

由Gronwall引理，可得

‖ η ‖ 1 2 + ∫ 0 t ‖ η x ‖ 2 d s ≤ C ∫ 0 t ( ‖ ξ ‖ 2 + ‖ ξ t ‖ 2 + h 2 ‖ u t ‖ 1 2 ) d s   . (25)

为估计 ‖ u t − u h t ‖ 1 ，在(20)中取 v h = η t , q h = η x t ，则有

( η t , γ h η t ) + α ( η x , η x t ) + β ( η x t , η x t ) − σ ( η , η x t ) = − ( ξ t , γ h η t ) − ( u t , ( I − γ h ) η t ) + ( ξ , η x t ) − δ 2 ( f ( u h ) − f ( u ) , η x t ) . (26)

利用迁移算子的性质，并应用Cauchy-Schwarz不等式和Young不等式，可得

3 4 ‖ η t ‖ 2 + β ‖ η x t ‖ 2 ≤ C ( ‖ ξ t ‖ 2 + ‖ η ‖ 2 + h 2 ‖ u t ‖ 1 2 + ‖ ξ ‖ 2 + ‖ η x ‖ 2 ) + 3 8 ‖ η t ‖ 2 + β 2 ‖ η x t ‖ 2 . (27)

整理可得

3 8 ‖ η t ‖ 2 + β 2 ‖ η x t ‖ 2 ≤ C ( ‖ η ‖ 1 2 + ‖ ξ ‖ 2 + ‖ ξ t ‖ 2 + h 2 ‖ u t ‖ 1 2 ) . (28)

将(25)式代入到(28)式，得到

‖ η t ‖ 1 2 ≤ C ( ‖ ξ ‖ 2 + ‖ ξ t ‖ 2 + h 2 ‖ u t ‖ 1 2 ) + C ∫ 0 t ( ‖ ξ ‖ 2 + ‖ ξ t ‖ 2 + h 2 ‖ u t ‖ 1 2 ) d s   . (29)

为估计 ‖ p − p h ‖ ，在(20)式中取 q h = ε ，则有

‖ ε ‖ 2 ≤ C ( ‖ ξ ‖ 2 + ‖ η ‖ 2 + ‖ η x ‖ 2 + ‖ η x t ‖ 2 ) + 1 2 ‖ ε ‖ 2 . (30)

将(25)式和(29)式代入到(30)式，得到

‖ ε ‖ 2 ≤ C ( ‖ ξ ‖ 2 + ‖ ξ t ‖ 2 + h 2 ‖ u t ‖ 1 2 ) + C ∫ 0 t ( ‖ ξ ‖ 2 + ‖ ξ t ‖ 2 + h 2 ‖ u t ‖ 1 2 ) d s   . (31)

最后利用投影估计和三角不等式完成定理的证明。

令 0 = t 0 < t 1 < ⋯ < t M = T 是时间区域 [ 0 , T ] 的剖分， t n = n Δ t ，其中时间步长 Δ t = T / M ，M是某一正整数。为了构造全离散格式，对于 [ 0 , T ] 上函数j，记 φ n = φ ( t n ) ， ∂ t φ n = ( φ n − φ n − 1 ) / Δ t 。并记 u h n 和 p h n 分别是u和p在 t = t n 处的全离散解，现在给出线性向后Euler全离散混合有限体积元格式：求 { u h n , p h n } ∈ H 0 h × L h , ( n = 1 , 2 , ⋯ , M ) ，满足

{ ( a )     ( ∂ t u h n , γ h v h ) − ( v h x , p h n ) = 0 ,   ∀ v h ∈ H 0 h , ( b )     ( p h n , q h ) − σ ( u h n , q h ) + α ( u h x n , q h ) + β ( ∂ t u h x n , q h ) = δ 2 ( f ( u h n − 1 ) , q h ) ,   ∀ q h ∈ L h , ( c )     u h n = π h u 0 . (32)

首先给出全离散格式的解的存在唯一性证明。

定理4.1：设 u h 0 和 { u h i , p h i } i = 1 n − 1 是已知的，则线性向后Euler混合有限体积元格式(32)存在唯一解 { u h n , p h n } 。

证明：全离散格式(32)可以改写为，对 ∀ v h ∈ H 0 h , ∀ q h ∈ L h ，有

Δ t ( p h n , q h ) − Δ t σ ( u h n , q h ) + α Δ t ( u h x n , q h ) + β ( u h x n , q h ) = β ( u h x n − 1 , q h ) + δ Δ t 2 ( f ( u h n − 1 ) , q h ) , (33)

( u h n , γ h v h ) − Δ t ( v h x , p h n ) = ( u h n − 1 , γ h v h ) . (34)

由于 { φ j } j = 1 N − 1 和 { χ A j } j = 1 N − 1 分别是有限元空间 H 0 h 和 L h 的基函数，则 u h n ∈ H O h , p h n ∈ L h 可以表示为

u h n ( x ) = ∑ j = 1 N − 1 u h j n φ j ( x ) ,     p h n ( x ) = ∑ j = 1 N − 1 p h j n χ A j .

将上述表达式应用在(33)式和(34)式上，并取 v h = ϕ i ( i = 1 , 2 , ⋯ , N − 1 ) , v h = χ A i ( i = 0 , 1 , 2 , ⋯ , N − 1 ) ，则全离散格式可以写为下述矩阵形式：求 { u ˜ h n , p ˜ h n } , ( n ≥ 1 ) ，使得

( Δ t A ( β + α Δ t ) B − σ Δ t C − Δ t B T D ) ( p ˜ h n u ˜ h n ) = ( β B u ˜ h n − 1 + 0.5 δ Δ t F ( u ˜ h n − 1 ) D u ˜ h n − 1 ) , (35)

其中

u ˜ h n = ( u h 1 n , u h 2 n , ⋯ , u h N − 1 n ) T ,     P ˜ h n = ( P h 1 n , P h 2 n , ⋯ , P h N − 1 n ) T ,

A = ( χ A i , χ A j ) i = 0 , ⋯ , N − 1 ; j = 0 , ⋯ , N − 1 ,     B = ( ϕ j x , χ A i ) i = 0 , ⋯ , N − 1 ; j = 0 , ⋯ , N − 1 ,

C = ( ϕ j , χ A i ) i = 0 , ⋯ , N − 1 ; j = 0 , ⋯ , N − 1 ,     D = ( ϕ j , γ h ϕ i ) i = 1 , ⋯ , N − 1 ; j = 1 , ⋯ , N − 1 ,

F ( u ˜ h n − 1 ) = ( f ( u h n − 1 ) , χ A i ) i = 0 , ⋯ , N − 1 T

注意到矩阵A和D是对称正定矩阵，则(35)可变形为

其中 H = D u ˜ h n − 1 + β B T A − 1 B u ˜ h n − 1 + δ Δ t 2 B T A − 1 F ( u ˜ h n − 1 ) 。当 Δ t 充分小的时候，系统(36)的系数矩阵是可逆的，从而有唯一解，于是全离散格式(32)有唯一解。

下面给出全离散格式(32)的误差估计。将误差写为

u ( t n ) − u h n = u ( t n ) − π h u ( t n ) + π h u ( t n ) − u h n = ξ n + η n ,

p ( t n ) − p h n = p ( t n ) − R h p ( t n ) + R h p ( t n ) − p h n = ρ n + ε n .

将格式(7)与格式(32)作差，利用投影算子可得误差方程为

其中 τ n = u t n − ∂ t u n , λ n = u x t n − ∂ t u x n 。

定理4.2：设 { u h n , p h n } 是系统(32)的解，系统(3)的解 { u , p } 满足正则性质：

u , u t ∈ L ∞ ( H 2 ( Ω ) ∩ H 0 1 ( Ω ) ) ,     P ∈ L ∞ ( H 1 ( Ω ) ) ,     u t t ∈ L ∞ ( H 0 1 ( Ω ) ) ,

则存在与h和 Δ t 无关的正常数C，使得

max 1 ≤ n ≤ M ( ‖ u ( t n ) − u h n ‖ 1 + ‖ u t ( t n ) − ∂ t u h n ‖ 1 + ‖ p ( t n ) − p h n ‖ ) ≤ C ( h + Δ t ) . (38)

证明：在(37)中取 v h = η n , q h = η x n ，得到

( ∂ t η n , γ h η n ) − σ ( η n , η x n ) + α ( η x n , η x n ) + β ( ∂ t η x n , η x n ) = − ( u t n , ( I − γ h ) η n ) − ( τ n , γ h η n ) − β ( λ n , η x n ) − ( ∂ t ξ n , γ h η n )       + σ ( ξ n , η x n ) − δ 2 ( f ( u h n − 1 ) − f ( u n ) , η x n ) . (39)

通过简单的计算，可知

( ∂ t η x n , η x n ) ≥ 1 2 Δ t [ ‖ η x n ‖ 2 − ‖ η x n − 1 ‖ 2 ] , (40)

( ∂ t η n , γ h η n ) ≥ 1 2 Δ t [ ( η n , γ h η n ) − ( η n − 1 , γ h η n − 1 ) ] . (41)

应用上述不等式及Cauchy-Schwarz不等式和Young不等式，可得

通过计算可得

‖ f ( u h n − 1 ) − f ( u n ) ‖ 2 ≤ C ( ‖ ξ n − 1 ‖ 2 + ‖ η n − 1 ‖ 2 ) + C Δ t ∫ t n − 1 t n ‖ u t ‖ 2 d t (43)

在(42)式两端同乘 2 Δ t ，并令n由1到m求和，可得

( η m , γ h η m ) + ‖ η x m ‖ 2 + α Δ t ∑ n = 1 m ‖ η x n ‖ 2 ≤ C ( ‖ η 0 ‖ 1 2 + ‖ ξ 0 ‖ 2 ) + C Δ t ∑ n = 1 m ( ‖ η n ‖ 1 2 + h 2 ‖ u t n ‖ 1 2 + ‖ τ n ‖ 2 + ‖ λ n ‖ 2 + ‖ ξ n ‖ + ‖ ∂ t ξ n ‖ 2 + Δ t ∫ t n − 1 t n ‖ u t ‖ d t )     . (44)

注意到 η 0 = 0 ，选取 p h 0 满足 ξ 0 = 0 ，应用引理3.3及离散的Gronwall引理，可知

‖ η m ‖ 1 2 ≤ C Δ t ∑ n = 1 m ( h 2 ‖ u t n ‖ 1 2 + ‖ τ n ‖ 2 + ‖ λ n ‖ 2 + ‖ ξ n ‖ + ‖ ∂ t ξ n ‖ 2 + Δ t ∫ t n − 1 t n ‖ u t ‖ d t )     . (45)

下面来估计 ‖ u t ( t n ) − ∂ t u h n ‖ 1 ，在(37)式中取 v h = ∂ t η n , q h = ∂ t η x n ，可得

( ∂ t η n , γ h ∂ t η n ) − σ ( η n , ∂ t η x n ) + α ( η x n , ∂ t η x n ) + β ( ∂ t η x n , ∂ t η x n ) = − ( u t n , ( I − γ h ) ∂ t η n ) − ( τ n , γ h ∂ t η n ) − ( ∂ t ξ n , γ h ∂ t η n )     + σ ( ξ n , ∂ t η x n ) − δ 2 ( f ( u h n − 1 ) − f ( u n ) , ∂ t η x n ) − ( λ n , ∂ t η x n ) . (46)

根据(43)式及引理3.1，可得

3 4 ‖ ∂ t η n ‖ 2 + β ‖ ∂ t η x n ‖ 2 ≤ C ( ‖ η n ‖ 1 2 + h 2 ‖ u t n ‖ 1 2 + ‖ ξ n ‖ + ‖ ∂ t ξ n ‖ 2 + ‖ τ n ‖ 2 + ‖ λ n ‖ 2   + Δ t ∫ t n − 1 t n ‖ u t ‖ 2 d t + ‖ η n − 1 ‖ 2 + ‖ ξ n − 1 ‖ 2 ) + 3 8 ‖ ∂ t η n ‖ 2 + β ‖ ∂ t η x n ‖ 2 . (47)

将式(45)代入式(47)，可得

‖ ∂ t η n ‖ 1 2 ≤ C ( ‖ ξ n ‖ + ‖ ∂ t ξ n ‖ 2 + ‖ τ n ‖ 2 + ‖ λ n ‖ 2 + h 2 ‖ u t n ‖ 1 2 + Δ t ∫ t n − 1 t n ‖ u t ‖ 2 d t + ‖ ξ n − 1 ‖ 2 )                             + C Δ t ∑ k = 1 n ( ‖ ξ k ‖ + ‖ ∂ t ξ k ‖ 2 + ‖ τ k ‖ 2 + ‖ λ k ‖ 2 + h 2 ‖ u t k ‖ 1 2 + Δ t ∫ t k − 1 t k ‖ u t ‖ 2 d t )     . (48)

下面估计 ‖ p n − p h n ‖ 。在(37)式中取 q h = ε n ，利用引理3.1，可得

把(45)式和(48)式代入(49)式中，有

‖ ε n ‖ 2 ≤ C ( ‖ ξ n ‖ 2 + ‖ ∂ t ξ n ‖ 2 + ‖ τ n ‖ 2 + ‖ λ n ‖ 2 + h 2 ‖ u t n ‖ 1 2 + Δ t ∫ t n − 1 t n ‖ u t ‖ 2 d t + ‖ ξ n − 1 ‖ 2 )                     + C Δ t ∑ k = 1 n ( ‖ ξ k ‖ 2 + ‖ ∂ t ξ k ‖ 2 + ‖ τ k ‖ 2 + ‖ λ k ‖ 2 + h 2 ‖ u t k ‖ 1 2 + Δ t ∫ t k − 1 t k ‖ u t ‖ 2 d t )     . (50)

为了(45)式和(50)式的右端项，利用Taylor展开的积分余项公式，可以得到

‖ ∂ t ξ n ‖ 2 ≤ C Δ t ∫ t n − 1 t n ‖ ξ t ‖ 2 d s ,   n ≥ 1 ,

‖ τ n ‖ 2 ≤ C Δ t ∫ t n − 1 t n ‖ u x t t ‖ 2 d s ,   ‖ λ n ‖ 2 ≤ C Δ t ∫ t n − 1 t n ‖ u x t t ‖ 2 d s ,   n ≥ 1.

最后使用上述不等式以及三角不等式完成定理的证明.

例1：本例考虑方程(1)中取 σ = δ = α = β = 1 ，计算区域 Ω = ( − 10 , 10 ) ， Γ = ( 0 , 10 ] ，初始函数 u 0 ( x ) = e − x 2 。空间和时间步长为 h = Δ t = 0.1 。图1为数值解 u h ( x , t ) 在时间 t = 0 , 1 , 2 , ⋯ , 10 时的图像，图2为数值解 p h ( x , t ) 在时间 t = 0 , 1 , 2 , ⋯ , 10 时的图像。

例2：文中的格式和理论分析的结果也适合于下述广义正则长波Burgers方程

{ ( a )     u t + σ u x + δ u u x − α u x x − β u x x t = g ( x , t ) ,     ( x , t ) ∈ Ω × Γ , ( b )     u ( a , t ) = u ( b , t ) = 0 ,     t ∈ Γ ¯ , ( c )     u ( x , 0 ) = u 0 ( x ) ,     x ∈ Ω , (51)

在(51)中取 σ = δ = α = β = 1 ，计算区域 Ω = ( 0 , 1 ) ， Γ = ( 0 , 10 ] ，初始函数 u 0 ( x ) = sin ( π x ) ，源项函数 g ( x , t ) = e − t ( π cos ( π x ) − sin ( π x ) + π 2 e − t sin ( 2 π x ) ) ，则精确解为

u ( x , t ) = e − t sin ( π x ) ,   p ( x , t ) = e − t sin ( π x ) + 1 2 e − 2 t ( sin ( π x ) ) 2 .

图1. 数值解u_h(x,t)

图2. 数值解p_h(x,t)

表1. 误差估计和收敛阶

表1给出了上述广义正则长波Burgers方程(51)的混合有限体积元格式在不同的空间和时间等距剖分下的误差估计和收敛精度，可以看出收敛精度与理论分析结果是一致的，这表明混合有限体积元方法求解广义正则长波Burgers方程和正则长波Burgers方程是有效的。

国家自然科学基金(11701299, 11761053)，内蒙古自治区自然科学基金(2016BS0105, 2017MS0107)，内蒙古自治区高等学校青年科技英才支持计划(NJYT-17-A07)。

白 雪,李 宏,方志朝. 正则长波Burgers方程的混合有限体积元方法 Mixed Finite Volume Element Method for Regularized Long Wave Burgers Equation[J]. 应用数学进展, 2018, 07(03): 257-268. https://doi.org/10.12677/AAM.2018.73032