Next Article in Journal
Generalized Galbraith’s Test: Characterization and Applications to Anonymous IBE Schemes
Next Article in Special Issue
Defect Detection in Atomic Resolution Transmission Electron Microscopy Images Using Machine Learning
Previous Article in Journal
Investigation of the Effect of a Flapping Propeller on Its Aerodynamic Performance
Previous Article in Special Issue
Numerical Investigation of MHD Pulsatile Flow of Micropolar Fluid in a Channel with Symmetrically Constricted Walls
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Conservative and Implicit Second-Order Nonlinear Numerical Scheme for the Rosenau-KdV Equation

Harbin Engineering University, Harbin 150001, China
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(11), 1183; https://doi.org/10.3390/math9111183
Submission received: 2 March 2021 / Revised: 28 April 2021 / Accepted: 4 May 2021 / Published: 24 May 2021
(This article belongs to the Special Issue Applied Mathematics and Computational Physics)

Abstract

:
In this paper, for solving the nonlinear Rosenau-KdV equation, a conservative implicit two-level nonlinear scheme is proposed by a new numerical method named the multiple integral finite volume method. According to the order of the original differential equation’s highest derivative, we can confirm the number of integration steps, which is just called multiple integration. By multiple integration, a partial differential equation can be converted into a pure integral equation. This is very important because we can effectively avoid the large errors caused by directly approximating the derivative of the original differential equation using the finite difference method. We use the multiple integral finite volume method in the spatial direction and use finite difference in the time direction to construct the numerical scheme. The precision of this scheme is O ( τ 2 + h 3 ) . In addition, we verify that the scheme possesses the conservative property on the original equation. The solvability, uniqueness, convergence, and unconditional stability of this scheme are also demonstrated. The numerical results show that this method can obtain highly accurate solutions. Further, the tendency of the numerical results is consistent with the tendency of the analytical results. This shows that the discrete scheme is effective.

1. Introduction

Proposed by Korteweg and de Vries, the Korteweg–de Vries (KdV) equation,
u t + u u x + u x x x = 0
has been widely studied. It can describe ion–phonon waves, magnetic fluid waves in cold plasma, unidirectional shallow water waves with small amplitude and long waves, and other wave processes in some physical and biological systems.
It has a wide range of physical applications, so there is great interest in this equation. A great many numerical methods have been proposed to obtain the numerical solutions of KdV equations [1,2,3,4,5]. In addition, [6] developed a new integral equation using the negative-order KdV equation and derived multiple soliton solutions, while [7] created various negative-order KdV equations in (3 + 1) dimensions and discussed the solutions for each derived model.
Given the shortcomings of the KdV equation in describing wave–wave and wave–wall interactions, Rosenau [8,9] proposed the Rosenau equation to cope with the compact discrete dynamic system.
u t + u x x x x t + u x + u u x = 0
The existence, uniqueness, and regularity of solutions were derived by Park [10]. Since then, several numerical methods have been studied for the Rosenau equation. For example, ref [11] used the Petviashvili iteration method to construct numerical solitary wave solutions; ref [12] applied Galerkin mixed finite element methods to (2) by employing a splitting technique; ref [13] discussed new methods to expand solutions for wave equations like Rosenau-type equations with damping terms; ref [14] constructed an implicit Crank–Nicolson formula of the mixed finite element method for nonlinear fourth-order Rosenau equations; and [15] proposed a meshfree method based on the radial basis function for the Rosenau equation and other higher-order partial differential equations(PDEs). The long-time behavior of solutions was investigated in [16].
To better study nonlinear waves, the viscous term   u x x x needs to be included.
u t + u x x x x t + u x + u u x + u x x x = 0
Equation (3) is usually called the Rosenau-KdV equation. The authors of [17,18] proposed conservative schemes for the Rosenau-KdV equation based on the finite difference method. The authors of [19] proposed a Crank–Nicolson meshless spectral radial point interpolation (CN-MSRPI) method for the nonlinear Rosenau-KdV equation. The authors of [20] solved the equation by the first-order Lie–Trotter and second-order Strang time-splitting techniques combined with quintic B-spline collocation, while [21] studied numerical solutions by using the subdomain method based on sextic B-spline basis functions. Although various methods have been proposed, we wonder whether there might be a new method with higher accuracy and efficiency that can keep some properties of the original partial differential equation. Further, research on the Rosenau-KdV equation under certain conditions is relatively lacking.
In this paper, we consider the Rosenau-KdV Equation (4) with initial condition
u ( x , 0 ) = u 0 ( x ) ,     x [ x l , x r ]
and boundary conditions
u ( x l , t ) = u ( x r , t ) = 0 ,     u x ( x l , t ) = u x ( x r , t ) = 0 ,   t [ 0 , T ]
Here,   u 0 ( x )   is a known smooth function, and x l and x r are, respectively, the left border and the right border of   x .
Theorem 1.
The system (3)–(5) satisfies the following conservative property:
E ( t ) = u L 2 2 + u x x L 2 2 = E ( 0 ) = Const .
Here, u   L 2 2 = x l x r u 2   d x .
Proof. 
Integrate both sides of Equation (3) from x l to x r and apply (5); we thus obtain
x l x r ( u t + u x x x x t + u x x + u x + u u x ) u   d x = 1 2 t ( u   L 2 2 + u x x L 2 2 ) = 0
Let E ( t ) = u L 2 2 + u x x L 2 2 . Then we get
E ( t ) = u L 2 2 + u x x L 2 2 = E ( 0 ) = Const .
Hence, the system (3)–(5) meets the conservative property. □
In this paper, we present a two-level implicit nonlinear discrete scheme for the Rosenau-KdV Equations (3)–(5) by using a new method named the multiple integral finite volume method (MIFVM). The remaining contents of this paper are arranged as follows: In Section 2, we introduce MIFVM in detail and propose a numerical scheme. The conservative property of this scheme is also discussed. In Section 3, the solvability of this numerical scheme is derived. Then, in Section 4, we show some prior estimates. According to the prior estimates, we demonstrate the convergence with order   O ( τ 2 + h 3 ) and unconditional stability of this numerical scheme in Section 5. In Section 6, the uniqueness of this numerical solution is verified with the classic theorem. Finally, we verify the effectiveness of the numerical scheme via some numerical experiments in Section 7.

2. A Two-Level Implicit Nonlinear Discrete Scheme and Its Conservative Law

2.1. Notation

Let   h   and   τ   be uniform step sizes in the spatial and temporal directions, respectively. Let   x j = x l + j h ( j = 0 , 1 , , J ) , t n = n τ ( n = 0 , 1 , , N ) , where h = ( x r x l ) / J , τ = T / N . Further, let u j = u j ( t ) = u ( x l + j h , t ) , u j n = u ( x l + j h , t n ) , Z h 0 = { u j | u 0 = u J = 0 , j = 0 , 1 , , J } , and Ω h = { x j | j = 0 , 1 , , J } . In this paper, we let   C   denote a generic positive constant independent of   h   and   τ   . The difference operators, inner product, and norms we defined are shown below.
u j n + 1 2 = u j n + 1 + u j n 2 ( u j n ) x = u j + 1 n u j n h ,   ( u j n ) x ¯ = u j n u j 1 n h ,   ( u j n ) x ^ = u j + 1 n u j 1 n 2 h , ( u j n ) x x ¯ = ( u j n ) x ¯ x = u j + 1 n 2 u j n + u j 1 n h 2 ,   ( u j n ) x x ¯ x ^ = u j + 2 n 2 u j + 1 n + 2 u j 1 n u j 2 n 2 h 3 , ( u j n ) x x x ¯ x ¯ = u j + 2 n 4 u j + 1 n + 6 u j n 4 u j 1 n + u j 2 n h 4 ,   ( u j n + 1 2 ) t ^ = u j n + 1 u j n τ , u n = ( u n , u n ) ,   u n = max x j Ω h | u j n | ,   ( u n , v n ) = j = 0 J u j n v j n h
However, we should note that if the inner product operates on different functions, there will be different ranges of values of j, for example, ( u x n , v x n ) = j = 0 J 1 ( u j n ) x ( v j n ) x h and ( u x ¯ n , v x ¯ n ) = j = 1 J ( u j n ) x ¯ ( v j n ) x ¯ h .
Lemma 1.
For any two mesh functions u ,   v Z h 0 , the following equations hold.
( u x , v ) = ( u , v x ¯ ) , ( u x x ¯ , v ) = ( u x , v x )   , ( u x ^ , v ) = ( u , v x ^ )
Furthermore, if ( u 0 n ) x x ¯ = ( u J n ) x x ¯ = 0 , then ( ( u n ) x x x ¯   x ¯ , u n ) = u x x n 2 .
Lemma 2.
For any mesh function u Z h 0 , the following equation holds.
u x ^ 2 u x 2
Lemma 3.
For any discrete function u Z h 0 , we have
( φ ( u ) , u ) = j = 1 J 1 1 3 ( u j ) x ^ ( u j 1 + u j + u j + 1 ) u j h = 0 ,
where φ ( u j ) = 1 3 ( u j ) x ^ ( u j 1 + u j + u j + 1 ) .
Proof. 
Because u Z h 0 , we have
( φ ( u ) , u ) = 1 6 h j = 1 J 1 [ u j + 1 u j u j u j 1 + ( u j + 1 ) 2 ( u j 1 ) 2 ] u j = 1 6 h [ j = 1 J 2 ( u j + u j + 1 ) u j u j + 1 j = 2 J 1 ( u j 1 + u j ) u j 1 u j ] = 0

2.2. The Multiple Integral Finite Volume Method(MIFVM)

In this paper, we use a method named MIFVM to construct a two-level implicit nonlinear scheme for the Rosenau-KdV Equations (3)–(5).The method uses multiple integrals and combines the finite difference method with the finite volume method. We thus discretize the original PDE into separate spatial and temporal directions.
In the spatial x   direction, firstly, by multiple integrals, we turn the original differential Equation (3), with unknown function u and its derivative, into an integral equation with only the unknown function. This is very important because we can effectively avoid the large errors caused by directly approximating the derivative of the original differential equation using the finite difference method. We use the multiple integral finite volume method in the spatial direction and use finite difference in the time direction to construct the numerical scheme. Firstly, in the spatial direction, the number of integration steps     m   depends on the order of the highest derivative in the   x   direction of the original PDE. Considering the original Equation (3), the order of the highest derivative in the   x   direction is four, so m = 2 4 1 = 15 . Now, we define a 15-time integral,
x x x x u = x x x x u ( x , t )   d e f _ _ x j + ε 7 x j + ε 8 d x f 2 x j + ε 6 x j + ε 7 d x f 1 x j + ε 5 x j + ε 6 d x e 2 x j x j + ε 5 d x e 1 x j ε 4 x j d x d 2 x j ε 3 x j ε 4 d x d 1 x j ε 2 x j ε 3 d x c 2 x j ε 1 x j ε 2 d x c 1 x f 1 x f 2 d x f   x e 1 x e 2 d x e x d 1 x d 2 d x d x c 1 x c 2 d x c x e x f d x b x c x d d x a x a x b u ( x , t )   d x
and we treat original Equation (2) using integral (12). Then, we can get
x j + ε 7 x j + ε 8 d x f 2 x j + ε 6 x j + ε 7 d x f 1 x j + ε 5 x j + ε 6 d x e 2 x j x j + ε 5 d x e 1 x j ε 4 x j d x d 2 x j ε 3 x j ε 4 d x d 1 x j ε 2 x j ε 3 d x c 2   x j ε 1 x j ε 2 d x c 1 x f 1 x f 2 d x f x e 1 x e 2 d x e x d 1 x d 2 d x d x c 1 x c 2 d x c x e x f d x b x c x d d x a x a x b u t   d x + x j + ε 7 x j + ε 8 d x f 2 x j + ε 6 x j + ε 7 d x f 1 x j + ε 5 x j + ε 6 d x e 2 x j x j + ε 5 d x e 1 x j ε 4 x j d x d 2 x j ε 3 x j ε 4 d x d 1 x j ε 2 x j ε 3 d x c 2   x j ε 1 x j ε 2 d x c 1 x f 1 x f 2 d x f x e 1 x e 2 d x e x d 1 x d 2 d x d x c 1 x c 2 d x c x e x f d x b x c x d d x a x a x b u x x x x t   d x + x j + ε 7 x j + ε 8 d x f 2 x j + ε 6 x j + ε 7 d x f 1 x j + ε 5 x j + ε 6 d x e 2 x j x j + ε 5 d x e 1 x j ε 4 x j d x d 2 x j ε 3 x j ε 4 d x d 1 x j ε 2 x j ε 3 d x c 2   x j ε 1 x j ε 2 d x c 1 x f 1 x f 2 d x f x e 1 x e 2 d x e x d 1 x d 2 d x d x c 1 x c 2 d x c x e x f d x b x c x d d x a x a x b u x   d x + x j + ε 7 x j + ε 8 d x f 2 x j + ε 6 x j + ε 7 d x f 1 x j + ε 5 x j + ε 6 d x e 2 x j x j + ε 5 d x e 1 x j ε 4 x j d x d 2 x j ε 3 x j ε 4 d x d 1 x j ε 2 x j ε 3 d x c 2   x j ε 1 x j ε 2 d x c 1 x f 1 x f 2 d x f x e 1 x e 2 d x e x d 1 x d 2 d x d x c 1 x c 2 d x c x e x f d x b x c x d d x a x a x b u u x   d x + x j + ε 7 x j + ε 8 d x f 2 x j + ε 6 x j + ε 7 d x f 1 x j + ε 5 x j + ε 6 d x e 2 x j x j + ε 5 d x e 1 x j ε 4 x j d x d 2 x j ε 3 x j ε 4 d x d 1 x j ε 2 x j ε 3 d x c 2   x j ε 1 x j ε 2 d x c 1 x f 1 x f 2 d x f x e 1 x e 2 d x e x d 1 x d 2 d x d x c 1 x c 2 d x c x e x f d x b x c x d d x a x a x b u x x x   d x = 0  
We then use Lagrange interpolation to approximate u ( x j ± ε i , t ) ( i = 1 , 2 , , 8 ) , because they aren’t defined on grid nodes. In addition, to obtain a high-precision numerical scheme, the following Lagrange interpolation polynomials are used.
u ( x , t ) = ( x x j ) ( x x j + 1 ) 2 h 2 u j 1 ( t ) ( x x j 1 ) ( x x j + 1 ) h 2 u j ( t ) + ( x x j 1 ) ( x x j ) 2 h 2 u j + 1 ( t ) + O ( h 3 )
u ( x , t ) = ( x x j 1 ) ( x x j + 1 ) ( x x j + 2 ) 12 h 3 u j 2 ( t ) ( x x j 2 ) ( x x j + 1 ) ( x x j + 2 ) 6 h 4 u j 1 ( t ) + ( x x j 2 ) ( x x j 1 ) ( x x j + 2 ) 6 h 4 u j + 1 ( t ) ( x x j 2 ) ( x x j 1 ) ( x x j + 1 ) 12 h 4 u j + 2 ( t ) + O ( h 4 )
and
u ( x , t ) = ( x x j 1 ) ( x x j ) ( x x j + 1 ) ( x x j + 2 ) 24 h 4 u j 2 ( t ) ( x x j 2 ) ( x x j ) ( x x j + 1 ) ( x x j + 2 ) 6 h 4 u j 1 ( t ) + ( x x j 2 ) ( x x j 1 ) ( x x j + 1 ) ( x x j + 2 ) 4 h 4 u j ( t ) ( x x j 2 ) ( x x j 1 ) ( x x j ) ( x x j + 2 ) 6 h 4 u j + 1 ( t ) + ( x x j 2 ) ( x x j 1 ) ( x x j ) ( x x j + 1 ) 24 h 4 u j + 2 ( t ) + O ( h 5 )
Secondly, in the temporal direction, we use center difference,
( u j n + 1 2 ) t = u j n + 1 u j n τ + O ( τ 2 )
to approximate the one-order derivative. Then, the numerical scheme will possess two-order accuracy in the temporal direction.
With the 15-time integral, Lagrange interpolation, and center difference, we obtain a series of numerical schemes with eight parameters,   ε i ( i = 1 , 2 , , 8 ) . As soon as we identify the eight parameters, we obtain a specific scheme. In fact, finally, we want to obtaina specific scheme that can keep some properties of the original PDE, such as the conservative property.

2.3. A Two-Level Implicit Nonlinear Discrete Scheme

According to the specific steps introduced above, to retain theenergy conservative property of problem (3)–(5), we choose ε 1 = ε 4 = ε 5 = ε 8 = 3 h and ε 2 = ε 3 = ε 6 = ε 7 = 3 h / 3 . Now, let us substitute the eight parameters and (17) into (13). After simplifying, we obtain a two-level implicit nonlinear discrete scheme for (3)–(5). This is presented below.
        1 9 ( ( u j 1 n + 1 2 ) t ^ + 7 ( u j n + 1 2 ) t ^ + ( u j + 1 n + 1 2 ) t ^ ) + ( u j n + 1 2 ) x x x ¯ x ¯ t ^ + ( u j n + 1 2 ) x ^ + 1 3 ( u j n + 1 2 ) x ^ ( u j 1 n + 1 2 + u j n + 1 2 + u j + 1 n + 1 2 ) + ( u j n + 1 2 ) x x ¯ x ^ = 0 ,                                                                                                                                                               1 j J 1 ,       0 n N 1
u j 0 = u 0 ( x j ) ,   1 j J 1
u 0 n = u J n = 0 ,   ( u 0 n ) x = ( u J n ) x = 0 ,     0 n N 1

2.4. Conservative Law of the Discrete Scheme

Theorem 2.
The two-level implicit nonlinear numerical scheme (18) possesses the following property:
E n = 7 9 u n 2 + 2 h 9 j = 0 J 1 u j n u j + 1 n + u x x n 2 = E n 1 = = E 0
Proof. 
Computing the inner product of (18) with   2 u n + 1 2 ( i . e . u n + 1 + u n ) , we have
    7 9 τ ( u n + 1 2 u n 2 ) + 2 h 9 τ ( j = 0 J 1 u j n + 1 u j + 1 n + 1 j = 0 J 1 u j n u j + 1 n ) + 1 τ ( u x x n + 1 2 u x x n 2 ) + ( φ ( u n + 1 2 ) , 2 u n + 1 2 ) = 0
Let   E n = 7 9 u n 2 + 2 h 9 j = 0 J 1 u j n u j + 1 n + u x x n 2 . Applying Lemma 3, we have
E n + 1 = E n .
Thus, we obtain     E n = = E 0 , which proves Theorem 2. It shows that this numerical scheme can retain the conservation property of the original PDE. □

3. Solvability

The following lemmas will be very helpful for proving the solvability of the discrete scheme (17)–(19).
Lemma 4.
Ref [22] Let   H   be a finite-dimensional inner product space; suppose that   g : H H , is continuous and there exists an α > 0 such that ( g ( x ) , x ) > 0 for all x H with x = α . Then there is x H such that g ( x ) = 0 and x α .
It is a classic theory and comes from the paper Existence and uniqueness theorems for solutions of nonlinear boundary value problems. This article was published in the Proceedings of Symposia in Applied Mathematics in 1965.
Lemma 5.
  2 M E is a positive definite matrix, where   E   is an identity matrix and
M = [ 1 0 0 0 0 0 0 7 1 0 0 0 0 1 7 0 0 0 0 0 0 7 1 0 0 0 0 1 7 0 0 0 0 0 0 1 ] ( J + 1 ) × ( J + 1 )
Proof. 
We know that
2 M E = [ 1 0 0 0 0 0 0 13 2 0 0 0 0 2 13 0 0 0 0 0 0 13 2 0 0 0 0 2 13 0 0 0 0 0 0 1 ] ( J + 1 ) × ( J + 1 )
Let P i ( 1 i J + 1 ) be ordered principal minor determinants of 2 M E . Obviously, we have P 1 = 1 , P 2 = 13 , P 3 = | 1 0 0 0 13 2 0 2 13 | = 165 , and P J = P J + 1 . In addition, from 2 M E , we have
P i = 13 P i 1 4 P i 2 ,       3 i J .
So, when i = 4 , we have P 4 = 13 P 3 4 P 2 > P 3 . Similarly, when 5 i J , we have
P J > P J 1 > > P 5 > P 4
Then, we have
P J + 1 = P J > P J 1 > > P 4 > P 3 > P 2 > P 1 > 0
Hence, 2 M E is a positive definite matrix. □
Theorem 3.
There is a u n + 1 Z h 0 that satisfies the discrete scheme (18)–(20).
Proof. 
Suppose that u 0 , u 1 , , u n 1 and u n satisfy (18)–(20) for n N 1 . Next, we prove that there is a u n + 1 that satisfies the discrete scheme (18)–(20).
Let g be an operator on Z h 0 defined by
g ( v ) = 2 9 A ( v u n ) + 2 ( v u n ) x x x ¯   x ¯ + τ v x ^ + τ 3 v x ^ ( v j 1 + v j + v j + 1 ) + τ v x x ¯ x ^
where
A = [ 1 0 0 0 0 0 1 7 1 0 0 0 0 1 7 0 0 0 0 0 0 7 1 0 0 0 0 1 7 1 0 0 0 0 0 1 ] ( J + 1 ) × ( J + 1 ) , N = [ 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 ] ( J + 1 ) × ( J + 1 ) , v = ( v 0 v 1 v 2 v J ) Z h 0
Obviously, g is continuous, A = M + N , and ( N   v , v ) = v 0 v 1 + v J 1 v J = 0 . Let λ 0 , λ 1 , , λ J be the eigenvalues of M and let λ min = { λ 0 , λ 1 , , λ J } . Take the inner product of (24) with v . By Lemma 1 and Lemma 3, we have
( g ( v ) , v ) = 2 9 ( A v , v ) 2 9 ( A u n , v ) + 22 v x x 2 22 ( v x x , u x x n ) 2 9 ( M v , v ) 2 9 A u n v + 2 v x x 2 2 v x x u x x n 2 9 ( λ 0 v 0 2 + + λ J v J 2 ) 1 9 A u n 2 1 9 v 2 + v x x 2 u x x n 2 ( 2 λ min 1 ) 9 v 2 1 9 A u n 2 u x x n 2
From Lemma 2, we can guarantee that 2 λ min 1 > 0 . Therefore, let
v 2 = A u n 2 + 9 u x x n 2 + 1 2 λ min 1
For all v Z h 0 , we have ( g ( v ) ,   v ) > 0 . From Lemma 4, there is a v * = u n + u n + 1 2 Z h 0 such that g (   v * ) = 0 . So, there is a u n + 1 = 2 v * u n that satisfies the scheme (18)–(20). □

4. Some Prior Estimates for the Discrete Scheme

Lemma 6.
Suppose that u 0 H 0 2 [ x l , x r ] ; then the solution of (3)–(5) satisfies
u C ,     u x C ,     u C ,     u x C
Proof. 
From (16), we have
u C ,     u x x C
Then, by the Holder inequality and the Schwarz inequality, we obtain
u x 2 = x l x r u x u x d x = u u x | x l x r x l x r u u x x d x = x l x r u u x x d x u u x x 1 2 ( u 2 + u x x 2 )
Thus, u x C . By the Sobolev inequality we have u C ,     u x C . □
Lemma 7.
[Discrete Sobolev Inequality] [22]. There are two constants C 1 and C 2 such that
u n C 1 u n + C 2 u x n
Lemma 8.
Assume that u Z h 0 ; then the solution of the discrete scheme (18)–(20) satisfies
u x x n C ,     u n C ,     u x n C , u n C , u x n C .
Proof. 
From (21) we have
u x x n C , u n C
By Lemma 1 and the Cauchy–Schwarz inequality, we obtain
u x 2 u n     u x x n 1 2 ( u x x n 2 + u n 2 ) C
Applying Lemma 7, we also obtain
u n C , u x n C .

5. Convergence and Stability of the Discrete Scheme

Let v j n + 1 2 = v ( x j , t n + 1 2 ) be the solution of (3)–(5). By substituting this into (18), we obtain the truncation error of scheme (17)–(19)
E r j n + 1 2 = 1 9 ( ( v j 1 n + 1 2 ) t ^ + 7 ( v j n + 1 2 ) t ^ + ( v j + 1 n + 1 2 ) t ^ ) + ( v j n + 1 2 ) x x x ¯ x ¯ t ^ + ( v j n + 1 2 ) x ^ + 1 3 ( v j n + 1 2 ) x ^ ( v j 1 n + 1 2 + v j n + 1 2 + v j + 1 n + 1 2 ) + ( v j n + 1 2 ) x x ¯ x ^
By Taylor expansion and Lagrange interpolation, we know that E r j n + 1 2 = O ( τ 2 + h 3 ) .
Theorem 4.
Suppose u 0 H 0 2 [ x l , x r ] and u ( x , t ) C 5 , 3 ; then the numerical solution u j n of scheme (17)–(19) converges to the solution v j n + 1 2 of the initia lboundary value problem (3)–(5) with order O ( τ 2 + h 3 ) by the norm     .
Proof. 
Let e j n + 1 2 = v j n + 1 2 u j n + 1 2 and subtract (18) from (33); we then have
  E r j n + 1 2 = 1 9 ( ( e j 1 n + 1 2 ) t ^ + 7 ( e j n + 1 2 ) t ^ + ( e j + 1 n + 1 2 ) t ^ ) + ( e j n + 1 2 ) x x x ¯ x ¯ t ^       + ( e j n + 1 2 ) x ^ + ( e j n + 1 2 ) x x ¯ x ^ + ( φ ( v j n + 1 2 ) φ ( u j n + 1 2 ) ) .
Computing the inner product of (34) with 2 e n + 1 2 ( i . e . e n + 1 + e n ) , we have
( E r j n + 1 2 , 2 e n + 1 2 ) = 7 9 τ e n + 1 2 + 2 h 9 τ j = 1 J 1 e j n + 1 e j + 1 n + 1 7 9 τ e n 2 2 h 9 τ j = 1 J 1 e j n e j + 1 n + 1 τ ( e x x n + 1 2 e x x n 2 ) + ( φ ( v n + 1 2 ) φ ( u n + 1 2 ) , 2 e n + 1 2 )
From Lemmas 6 and 7 and the Cauchy–Schwarz inequality, we obtain
  ( φ ( v n + 1 2 ) φ ( u n + 1 2 ) , 2 e n + 1 2 ) = 2 3 h j = 1 J 1 ( v j 1 n + 1 2 + v j n + 1 2 + v j + 1 n + 1 2 ) ( v j n + 1 2 ) x ^ e j n + 1 2 + 2 3 h j = 1 J 1 ( u j 1 n + 1 2 + u j n + 1 2 + u j + 1 n + 1 2 ) ( u j n + 1 2 ) x ^ e j n + 1 2 = 2 3 h j = 1 J 1 ( e j 1 n + 1 2 + e j n + 1 2 + e j + 1 n + 1 2 ) ( v j n + 1 2 ) x ^ e j n + 1 2 + 2 3 h j = 1 J 1 ( u j 1 n + 1 2 + u j n + 1 2 + u j + 1 n + 1 2 ) ( u j n + 1 2 ) x ^ e j n + 1 2       2 3 h j = 1 J 1 ( u j 1 n + 1 2 + u j n + 1 2 + u j + 1 n + 1 2 ) ( v j n + 1 2 ) x ^ e j n + 1 2 = 2 3 h j = 1 J 1 ( e j 1 n + 1 2 + e j n + 1 2 + e j + 1 n + 1 2 ) ( v j n + 1 2 ) x ^ e j n + 1 2   2 3 h j = 1 J 1 ( u j 1 n + 1 2 + u j n + 1 2 + u j + 1 n + 1 2 ) ( e j n + 1 2 ) x ^ e j n + 1 2 2 3 C h j = 1 J 1 ( | e j 1 n + 1 2 | + | e j n + 1 2 | + | e j + 1 n + 1 2 | ) | e j n + 1 2 | + 2 3 C h j = 1 J 1 | ( e j n + 1 2 ) x ^ | | e j n + 1 2 | C ( e n + 1 2 + e n 2 + e x ^ n + 1 2 + e x ^ n 2 )
In addition, we have
( E r j n + 1 2 , 2 e n + 1 2 ) = ( E r j n + 1 2 , e n + 1 + e n ) E r n + 1 2 2 + e n + 1 2 + e n 2 2
Substituting (36) and (37) into (35), with Lemma 3, we have
7 9 e n + 1 2 + 2 h 9 j = 1 J 1 e j n + 1 e j + 1 n + 1 7 9 e n 2 2 h 9 j = 1 J 1 e j n e j + 1 n + e x x n + 1 2 e x x n 2 τ E r n + 1 2 2 + C τ ( e n + 1 2 + e n 2 + e x n + 1 2 + e x n 2 + e x x n 2 e x x n + 1 2 )
Let B n = 7 9 e n 2 + 2 h 9 i = 1 J 1 e i n e i + 1 n + e x x n 2 + e x n 2 . Obviously, B 0 = 0 . Then, (38) can be rewritten as
B n + 1 B n τ E r n + 1 2 2 + C τ ( B n + 1 + B n )
When τ is sufficiently small that 1 C τ > 0 , we have
B n + 1 1 + C τ 1 C τ B n + τ 1 C τ E r n + 1 2 2 τ 1 C τ k = 0 n ( 1 + C τ 1 C τ ) n k E r k + 1 2 2 O 2 ( τ 2 + h 3 ) k = 1 n + 1 ( 1 + C τ 1 C τ ) k
Then we have
B n O 2 ( τ 2 + h 3 ) k = 1 n ( 1 + C τ 1 C τ ) k O 2 ( τ 2 + h 3 ) k = 1 n ( 1 + 2 C τ 1 C τ ) k O 2 ( τ 2 + h 3 )
That is, e n O ( τ 2 + h 3 ) , e x n O ( τ 2 + h 3 ) . Using Lemma 8, we have
e n O ( τ 2 + h 3 )
Similarly, we can prove the following theorem.
Theorem 5.
Under the conditions of Theorem 4, the solution u j n of discrete scheme (18)–(20) is unconditionally stable by the norm     .

6. Uniqueness of the Numerical Solution

Theorem 6.
The solution of the discrete scheme (18)–(20) is unique.
Proof. 
We assume that   u n   and w n are two different solutions of (18)–(20). Let S j n + 1 2 = w j n + 1 2 u j n + 1 2 . Then, we have
    1 9 ( ( S j 1 n + 1 2 ) t ^ + 7 ( S j n + 1 2 ) t ^ + ( S j + 1 n + 1 2 ) t ^ ) + ( S j n + 1 2 ) x x x ¯ x ¯ t ^ + ( S j n + 1 2 ) x ^ + ( S j n + 1 2 ) x x ¯ x ^ + 1 3 ( w j n + 1 2 ) x ^ ( w j 1 n + 1 2 + w j n + 1 2 + w j + 1 n + 1 2 ) 1 3 ( u j n + 1 2 ) x ^ ( u j 1 n + 1 2 + u j n + 1 2 + u j + 1 n + 1 2 ) = 0
By computing the inner product of (42) with   2 S n + 1 2 ( i . e . S n + 1 + S n ) , we obtain
    7 9 τ S n + 1 2 + 2 h 9 τ j = 1 J 1 S j n + 1 S j + 1 n + 1 7 9 τ S n 2 2 h 9 τ j = 1 J 1 S j n S j + 1 n + 1 τ ( S x x n + 1 2 S x x n 2 ) + ( φ ( w ) φ ( u ) , 2 e n + 1 2 ) = 0
Let   Z n = 7 9 S n 2 + 2 h 9 j = 1 J 1 S j n S j + 1 n + S x x n 2 + S x n 2 ; we know that   Z 0 = 0 . From (43) we obtain
Z n + 1 Z n C τ ( S n + 1 2 + S n 2 + S x x n + 1 2 + S x x n 2 + S x ^ n + 1 2 + S x ^ n 2 )
Similarly, while   1 2 C τ > 0 , we have
Z n + 1 ( 1 + β τ ) Z n ( 1 + β τ ) n + 1 Z 0 = 0
Hence, we have   S n 2 = 0 , where   β = 4 C 1 2 C τ . This implies that   u n = w n . The discrete scheme (18)–(20) is thus uniquely solvable. □

7. Results

7.1. Example

We consider the Rosenau-KdV equation
u t + u x x x x t + u x + u u x + u x x x = 0 ,       ( x , t ) [ 40 , 40 ] × [ 0 , 10 ]
with initial condition
u ( x , 0 ) = ( 35 312 313 35 24 ) sech 4 [ 1 24 2 313 26 x ] , x [ 40 , 40 ]
and boundary conditions
u ( 40 ,   t ) = u ( 40 ,   t ) = 0 ,     u x ( 40 ,   t ) = u x ( 40 ,   t ) = 0 ,   t [ 0 , 10 ]
The exact solution is given by
u ( x ,   t ) = ( 35 312 313 35 24 ) sech 4 { 2 313 26 24 [ x ( 1 2 + 1 26 313 ) t ] }

7.2. Figures, Tables, and Schemes

We discretize the problem (46)–(48) using the numerical scheme (18)–(20).
From Figure 1, Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6, we can see that the numerical solution is consistent with the exact solution.
In Table 1, the errors with various     h   and   τ   are given. It is obvious that the errors are reducing with decreasing     h   and   τ . Hence, our discrete scheme is reasonable. e n ( h ,   τ ) / e n ( h / 2 ,   τ / 2 ) and e n ( h ,   τ ) / e n ( h / 2 ,   τ / 2 ) are given in Table 2, which interprets the convergence rates of the numerical scheme with various h and τ and various norms. From Table 3, we can see that the discrete E n is conservative. This property is consistent with the original equation. The numerical experiment shows that our discrete scheme is efficient.

8. Conclusions

In this paper, a second-order implicit nonlinear discrete scheme for the Rosenau-KdV equation is proposed via the multiple integral finite volume method (MIFVM). The discrete scheme possesses the conservative property of the original equation. The solvability, uniqueness, convergence, and unconditional stability of the scheme were demonstrated in detail. Numerical experiments verified that the discrete scheme given by MIFVM is effective.

Author Contributions

Y.L. analyzed and interpreted the new numerical method, the multiple integral finite volume method. C.G. and Y.W. obtained this conservative nonlinear implicit numerical scheme and demonstrated the existence, uniqueness, convergence, and stability of this numerical scheme. Y.W. proved that the numerical scheme maintained the energy property of the original equation and verified the feasibility of this numerical scheme with a numerical experiment. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (No. 11526064) and the Fundamental Research Fund for the Central Universities (No. 3072020CF2408).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors are grateful to Sun Yue for her work and help. I thank the referees for their valuable work in their paper.

Conflicts of Interest

The authors’ interpretation of data or presentation of information wasn’t influenced by any personal or other organizations. The interpretation and presentation will not influence the interests of any personal or other organizations, either.

References

  1. Yokus, A.; Bulut, H. Numerical simulation of KdV equation by finite difference method. Indian J. Phys. 2018, 92, 1571–1575. [Google Scholar] [CrossRef]
  2. Başhan, A. A novel approach via mixed Crank–Nicolson scheme and differential quadrature method for numerical solutions of solitons of mKdV equation. Pramana-J. Phys. 2019, 92, 84. [Google Scholar] [CrossRef]
  3. Başhan, A.; Yağmurlu, N.M.; Uçar, Y.; Esen, A. A new perspective for the numerical solutions of the cmKdV equation via modified cubic B-spline differential quadrature method. Int. J. Mod. Phys. C 2018, 29, 1850043. [Google Scholar] [CrossRef]
  4. Kong, D.; Xu, Y.; Zheng, Z. Numerical method for generalized time fractional KdV-type equation. Numer. Methods Partial. Differ. Equ. 2019, 36, 906–936. [Google Scholar] [CrossRef]
  5. Cairone, F.; Anandan, P.; Bucolo, M. Nonlinear systems synchronization for modeling two-phase microfluidics flows. Nonlinear Dyn. 2017, 92, 75–84. [Google Scholar] [CrossRef]
  6. Wazwaz, A.-M. A new integrable equation that combines the KdV equation with the negative-order KdV equation. Math. Methods Appl. Sci. 2018, 41, 80–87. [Google Scholar] [CrossRef]
  7. Wazwaz, A.-M. Negative-order KdV equations in (3 + 1) dimensions by using the KdV recursion operator. Waves Random Complex Media 2017, 27, 768–778. [Google Scholar] [CrossRef]
  8. Rosenau, P. A Quasi-Continuous Description of a Nonlinear Transmission Line. Phys. Scr. 1986, 34, 827–829. [Google Scholar] [CrossRef]
  9. Rosenau, P. Dynamics of Dense Discrete Systems: High Order Effects. Prog. Theor. Phys. 1988, 79, 1028–1042. [Google Scholar] [CrossRef] [Green Version]
  10. Park, M.A. On the Rosenau equation. Matemática Aplic. Comput. 1990, 9, 145–152. [Google Scholar]
  11. Erbay, H.; Erbay, S.; Erkip, A. Numerical computation of solitary wave solutions of the Rosenau equation. Wave Motion 2020, 98, 102618. [Google Scholar] [CrossRef]
  12. Atouani, N.; Ouali, Y.; Omrani, K. Mixed finite element methods for the Rosenau equation. J. Appl. Math. Comput. 2018, 57, 393–420. [Google Scholar] [CrossRef]
  13. Michihisa, H. New asymptotic estimates of solutions for generalized Rosenau equations. Math. Methods Appl. Sci. 2019, 42, 4516–4542. [Google Scholar] [CrossRef]
  14. Shi, D.; Jia, X. Superconvergence analysis of the mixed finite element method for the Rosenau equation. J. Math. Anal. Appl. 2020, 481, 123485. [Google Scholar] [CrossRef]
  15. Safdari-Vaighani, A.; Larsson, E.; Heryudono, A. Radial Basis Function Methods for the Rosenau Equation and Other Higher Order PDEs. J. Sci. Comput. 2018, 75, 1555–1580. [Google Scholar] [CrossRef] [Green Version]
  16. Wang, Y.; Feng, G. Large-time behavior of solutions to the Rosenau equation with damped term. Math. Methods Appl. Sci. 2017, 40, 1986–2004. [Google Scholar] [CrossRef]
  17. Wang, X.; Dai, W. A conservative fourth-order stable finite difference scheme for the generalized Rosenau–KdV equation in both 1D and 2D. J. Comput. Appl. Math. 2019, 355, 310–331. [Google Scholar] [CrossRef]
  18. Luo, Y.; Xu, Y.; Feng, M. Conservative Difference Scheme for Generalized Rosenau-KdV Equation. Adv. Math. Phys. 2014, 2014, 986098. [Google Scholar] [CrossRef] [Green Version]
  19. Hussain, M.; Haq, S. Numerical simulation of solitary waves of Rosenau–KdV equation by Crank–Nicolson meshless spectral interpolation method. Eur. Phys. J. Plus 2020, 135, 98. [Google Scholar] [CrossRef]
  20. Kutluay, S.; Karta, M.; Yağmurlu, N.M. Operator time-splitting techniques combined with quintic B-spline collocation method for the generalized Rosenau–KdV equation. Numer. Methods Partial. Differ. Equ. 2019, 35, 2221–2235. [Google Scholar] [CrossRef]
  21. Browder, F.E. Existence and uniqueness theorems for solutions of nonlinear boundary value problems. In Proceedings of the Sum of Squares: Theory and Applications; American Mathematical Society: Providence, RI, USA, 1965; Volume 17, pp. 24–49. [Google Scholar]
  22. Zhou, Y.L. Applications of Discrete Functional Analysis to the Finite Difference Method, 1st ed.; International Academic Publishers: Beijing, China, 1991; pp. 3–20. [Google Scholar]
Figure 1. Numerical solution and exact solution with h = τ = 1 / 4 , t = 0 .
Figure 1. Numerical solution and exact solution with h = τ = 1 / 4 , t = 0 .
Mathematics 09 01183 g001
Figure 2. Numerical solution and exact solution with h = τ = 1 / 4 , t = 5 .
Figure 2. Numerical solution and exact solution with h = τ = 1 / 4 , t = 5 .
Mathematics 09 01183 g002
Figure 3. Numerical solution and exact solution with h = τ = 1 / 4 , t = 10 .
Figure 3. Numerical solution and exact solution with h = τ = 1 / 4 , t = 10 .
Mathematics 09 01183 g003
Figure 4. Numerical solution and exact solution with   h = τ = 1 / 8   , t = 0 .
Figure 4. Numerical solution and exact solution with   h = τ = 1 / 8   , t = 0 .
Mathematics 09 01183 g004
Figure 5. Numerical solution and exact solution with   h = τ = 1 / 8   , t = 5 .
Figure 5. Numerical solution and exact solution with   h = τ = 1 / 8   , t = 5 .
Mathematics 09 01183 g005
Figure 6. Numerical solution and exact solution with   h = τ = 1 / 8   , t = 10 .
Figure 6. Numerical solution and exact solution with   h = τ = 1 / 8   , t = 10 .
Mathematics 09 01183 g006
Table 1. The errors at different times with various h and τ .
Table 1. The errors at different times with various h and τ .
h = τ = 1 / 4 h = τ = 1 / 8
e n e n e n e n
t = 2 2.31861708 × 10 4 1.23912354 × 10 3 5.80241511 × 10 5 4.38690718 × 10 4
t = 4 4.68968494 × 10 4 2.46559087 × 10 3 1.17406720 × 10 4 8.73135975 × 10 4
t = 6 7.06502002 × 10 4 3.66896388 × 10 3 1.77078291 × 10 4 1.30146495 × 10 3
t = 8 9.33077015 × 10 4 4.84157574 × 10 3 2.34850462 × 10 4 1.72786004 × 10 3
t = 10 1.14711186 × 10 3 5.97898965 × 10 3 2.91676007 × 10 4 2.16772647 × 10 3
Table 2. The convergence rateswith various h and τ and various norms.
Table 2. The convergence rateswith various h and τ and various norms.
e n ( h ,   τ ) / e n ( h / 2 ,   τ / 2 ) e n ( h ,   τ ) / e n ( h / 2 ,   τ / 2 )
τ = h = 1 / 2 τ = h = 1 / 4 τ = h = 1 / 8 τ = h = 1 / 2 τ = h = 1 / 4 τ = h = 1 / 8
t = 2 2.813341342.824594843.984081693.99595175
t = 4 2.812871472.823833793.968571533.99439224
t = 6 2.812101182.819103103.980330463.98977197
t = 8 2.811121832.802064763.963164633.97306867
t = 10 2.810011482.758184543.974377013.93282902
Table 3. Discrete   E n   values at different times with various h and τ .
Table 3. Discrete   E n   values at different times with various h and τ .
h = τ = 1 / 2 h = τ = 1 / 4 h = τ = 1 / 8
t = 2 3.086750126.1734919912.34697937
t = 4 3.086766516.1735009512.34698364
t = 6 3.086790876.1735143212.34698593
t = 8 3.086819186.1735299612.34696462
t = 10 3.086848446.1735462212.34685552
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Guo, C.; Wang, Y.; Luo, Y. A Conservative and Implicit Second-Order Nonlinear Numerical Scheme for the Rosenau-KdV Equation. Mathematics 2021, 9, 1183. https://doi.org/10.3390/math9111183

AMA Style

Guo C, Wang Y, Luo Y. A Conservative and Implicit Second-Order Nonlinear Numerical Scheme for the Rosenau-KdV Equation. Mathematics. 2021; 9(11):1183. https://doi.org/10.3390/math9111183

Chicago/Turabian Style

Guo, Cui, Yinglin Wang, and Yuesheng Luo. 2021. "A Conservative and Implicit Second-Order Nonlinear Numerical Scheme for the Rosenau-KdV Equation" Mathematics 9, no. 11: 1183. https://doi.org/10.3390/math9111183

APA Style

Guo, C., Wang, Y., & Luo, Y. (2021). A Conservative and Implicit Second-Order Nonlinear Numerical Scheme for the Rosenau-KdV Equation. Mathematics, 9(11), 1183. https://doi.org/10.3390/math9111183

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop