Next Article in Journal
A Parallel Computational–Statistical Framework for Simulation of Turbulence: Applications to Data-Driven Fractional Modeling
Previous Article in Journal
Multicorn Sets of z¯k+cm via S-Iteration with h-Convexity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Time Two-Mesh Finite Difference Numerical Scheme for the Symmetric Regularized Long Wave Equation

School of Mathematics and Big Data, Hohhot Minzu College, Hohhot 010051, China
*
Authors to whom correspondence should be addressed.
Fractal Fract. 2023, 7(6), 487; https://doi.org/10.3390/fractalfract7060487
Submission received: 23 May 2023 / Revised: 16 June 2023 / Accepted: 16 June 2023 / Published: 19 June 2023

Abstract

:
The symmetric regularized long wave (SRLW) equation is a mathematical model used in many areas of physics; the solution of the SRLW equation can accurately describe the behavior of long waves in shallow water. To approximate the solutions of the equation, a time two-mesh (TT-M) decoupled finite difference numerical scheme is proposed in this paper to improve the efficiency of solving the SRLW equation. Based on the time two-mesh technique and two time-level finite difference method, the proposed scheme can calculate the velocity u ( x , t ) and density ρ ( x , t ) in the SRLW equation simultaneously. The linearization process involves a modification similar to the Gauss-Seidel method used for linear systems to improve the accuracy of the calculation to obtain solutions. By using the discrete energy and mathematical induction methods, the convergence results with O ( τ C 2 + τ F + h 2 ) in the discrete L -norm for u ( x , t ) and in the discrete L 2 -norm for ρ ( x , t ) are proved, respectively. The stability of the scheme was also analyzed. Finally, some numerical examples, including error estimate, computational time and preservation of conservation laws, are given to verify the efficiency of the scheme. The numerical results show that the new method preserves conservation laws of four quantities successfully. Furthermore, by comparing with the original two-level nonlinear finite difference scheme, the proposed scheme can save the CPU time.

1. Introduction

The regularized long wave (RLW) equation [1,2] is a nonlinear partial differential equation that mainly describes the evolution of waves in shallow water channels and ion acoustics, etc. It is a simplified version of the more complex Korteweg-de Vries (KdV) equation [3], which includes higher-order nonlinearities and dispersion effects. The symmetric regularized long wave (SRLW) equation [4] is a modified version of the RLW equation that includes a symmetry-breaking term. This term allows for the formation of asymmetric solutions, making the SRLW equation a more realistic model for waves in shallow water channels.
In this paper, the following initial boundary value problem of the SRLW equation is considered:
u t + ρ x + u u x u x x t = 0 , x L x x R , 0 < t T , ρ t + u x = 0 , x L x x R , 0 < t T , u x L , t = u x R , t = 0 , ρ x L , t = ρ x R , t = 0 , 0 < t T , u ( x , 0 ) = u 0 ( x ) , ρ ( x , 0 ) = ρ 0 ( x ) , x L x x R ,
where u ( x , t ) and ρ ( x , t ) are the fluid velocity and the density, respectively.
The SRLW equation has attracted significant attention and has been extensively studied in the literature. The existence of global attractors of the SRLW equation was studied in [5]. The travelling and solitary wave solutions of the SRLW equation were investigated by several methods including the exp-function method [6], the ( G / G )-expansion method [7,8], the tanh-function method [9,10], the Lie symmetry approach [11], the extended simple equation method [12], the analytical method [13], the tan ( ( ϕ ( ξ ) / 2 ) )-expansion method [14] and the sine–cosine method [15]. In [16,17], the spectral method and the Fourier pseudo-spectral method with a constraint operator were developed as approximations for the nonlinear term of the SRLW equation. Numerical solutions to the SRLW equation have also been studied by numerous methods, ranging from conservative finite difference schemes to mixed finite element methods. In their study, Wang et al. [18] introduced three conservative finite difference schemes that achieve second-order accuracy in both spatial and temporal domains. Their results indicate that each scheme is able to preserve energy conservation, but only the first scheme is capable of preserving mass conservation. Yimnet et al. [19] presented a novel finite difference method for the SRLW equation that utilizes a four-level average difference technique for solving the fluid velocity independently from the density. A coupled conservative three-level implicit scheme achieving fourth-order convergence rate was developed by Hu et al. [20]. Li [21] considered a weighted and compact conservative difference scheme that is decoupled and linearized in practical computation, thus requiring only the solution of two tridiagonal systems of linear algebraic equations at each time step. Bai et al. [22] investigated a two-layer conservative finite difference scheme for the SRLW equation with homogeneous boundary conditions and analyzed the scheme’s convergence and stability using a discrete functional analysis method. Xu et al. [23] applied a mixed finite element method to solve the dissipative SRLW equations with damping term. He et al. [24] developed a fourth-order accurate compact difference scheme for the SRLW equation for a single nonlinear velocity form and conducted theoretical analysis using the discrete energy method.
In terms of numerical computation, the time two-mesh (TT-M) method, when combined with either the finite element method or the finite difference method, offers better computational efficiency in solving a broad range of nonlinear partial differential equations. For instance, Liu et al. [25] proposed the fast TT-M finite element method for solving the fractional water wave model and applied it successfully to other fractional models. Yin et al. [26] developed a TT-M finite element algorithm for solving a space fractional Allen-Cahn model and analyzed the problem of parameter selection in detail. The TT-M finite element method was also leveraged by Liu et al. [27] to numerically solve the 2D Gray-Scott model with space fractional derivatives. A nonlinear distributed order diffusion model was efficiently solved using the TT-M algorithm in conjunction with the H 1 -Galerkin mixed finite element method by Wen et al. [28], both smooth and non-smooth solutions were considered. Additionally, Tian et al. [29] developed a finite element method equipped with the TT-M technique to solve the coupled Schrödinger-Boussinesq equations. Moreover, some studies have investigated combining the TT-M and finite difference methods to solve nonlinear fractional partial differential equations, such as the works of Qiu and Xu et al. [30,31], who proposed and analyzed a TT-M algorithm based on finite difference methods. The TT-M technique was also employed by Niu et al. [32] to develop a fast high-order compact difference scheme for the nonlinear distributed order fractional Sobolev model in porous media. Furthermore, He et al. [33] extended the application of the TT-M method by studying a primary scheme of second-order convergence in time and fourth-order in space for solving the nonlinear Schrödinger equation with a time two-mesh high-order compact difference scheme. Despite the extensive research on the TT-M method in various fields, to the best of our knowledge, no study on the application of the TT-M method combined with finite difference to the SRLW equation has been discovered. Hence, investigations on the TT-M finite difference method’s performance when applied to the SRLW equation are still required.
The SRLW Equation (1) is a coupled system and most existing schemes found in the literature are also coupled and require nonlinear implementation, which results in prolonged CPU processing time. The objective of this paper is to develop a novel difference scheme that offers the following three main advantages for solving the SRLW equation: (i) Based on the time two-mesh technique, we proposed a scheme that achieves decoupling and the nonlinear term of the system is linearized by using Taylor’s formula for a function with three variables, which is different from the literature [32,33]. In [32,33], the time two-mesh scheme is formulated by using Taylor’s formula for a function with one or two variables. As a result, our scheme becomes a linearized system in approximate numerical solution and can reduce the computational time. (ii) The scheme’s convergence and stability have been verified through detailed proof, and the theoretical analysis process is more complex than that of existing methods, since the linearized scheme contains a function of three variables. (iii) The linearization process involves a modification similar to the Gauss-Seidel method used for linear systems to improve the accuracy of the calculation for solutions. The conservation law of four quantities is preserved in the new scheme at the discrete level.
The remaining part of this article is organized as follows. In Section 2, some notations and useful lemmas are given. In Section 3, the TT-M finite difference numerical scheme is presented. In Section 4, the convergence and stability of the scheme is analyzed. In Section 5, some numerical results are provided to test the theoretical results and computational efficiency of the scheme. Finally, in Section 6, we provide the conclusions of the paper.

2. Notations and Some Lemmas

As usual, the time interval ( 0 , T ] and spatial interval x L , x R are divided into N and J uniform partitions. The following notations will be used in this paper:
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 t = u j n + 1 u j n τ , u j n + 1 2 = 1 2 ( u j n + 1 + u j n ) , u j n x x ¯ = u j + 1 n 2 u j n + u j 1 n h 2 ,
where τ , h denote the uniform time and spatial step length, respectively, x j = x L + j h ,   j = 0 , 1 , 2 , , J t n = n τ ,   n = 1 , 2 , , [ T / τ ] = N , superscript n denotes a quantity associated with the time level t n , subscript j denotes a quantity associated with space mesh point x j . In this paper, M denotes general constant, which may have a different value in a different place.
Since u 0 for x + or x , we may assume u 1 n = u J + 1 n = 0 , 1 n N for simplicity, where j = 1 and J + 1 are ghost points. Let H h , 0 denote the set of mesh functions u n defined on I h with boundary conditions u 1 n = u 0 n = u J n = u J + 1 n = 0 . For any two mesh functions u n , w n H h , 0 , we define the discrete inner product and norms as follows:
( u n , w n ) = h j = 1 J 1 u j n w j n , u n = ( u n , u n ) , u n = max 1 j J 1 | u j n | .
Next, we present some useful lemmas.
Lemma 1
(See [24]). For any mesh functions u n , w n H h , 0 , we have
( a ) u x n , w n = u n , w x ¯ n = u n , w x n , ( b ) u x x ¯ n , w n = u x n , w x n , ( c ) u x ^ n , w n = u n , w x ^ n .
Lemma 2
(See [33,34]). Assume that a sequence of non-negative real numbers a j j = 0 satisfying
a n + 1 α + β j = 0 n a j τ , n 0 ,
then there has the inequality a n + 1 α + τ β a 0 e β ( n + 1 ) τ , where α 0 , β and τ are positive constants.
Lemma 3
(See [22,34]). For any discrete mesh function u n H h , 0 , there exists constants C 1 and C 2 , such that
u n C 1 u n + C 2 u x n .

3. The TT-M Finite Difference Scheme

In this paper, we studied a TT-M finite difference fast numerical method for the SRLW Equation (1). In order to give the TT-M finite difference scheme, firstly, the time interval ( 0 , T ] is partitioned uniformly into P coarse time intervals and then each of them is divided into s 2 s Z + fine time intervals. The coarse time mesh with the nodes t k s = k τ C ( k = 1 , , P ) satisfying 0 = t 0 < t s < t 2 s < < t P s = T and the fine time mesh with the nodes t n = n τ F ( n = 1 , 2 , , P s = N ) satisfying 0 = t 0 < t 1 < t 2 < < t P s = T , where τ C = s τ F and τ F are the coarse time and the fine time step size, respectively.
Secondly, the truncation errors of the problem (1) are considered, let v j n = u x j , t n ,   φ j n = ρ x j , t n be the exact solutions of u x , t and ρ x , t in term of the point x j , t n , then we have
E r j n = ( v j n ) t + ( φ j n + 1 ) x ^ ( v j n ) x x ¯ t + 1 3 v j n + 1 2 ( v j n + 1 2 ) x ^ + ( v j n + 1 2 ) 2 x ^ ,
E s j n = ( φ j n ) t + ( v j n ) x ^ ,
v 0 n = v J n = 0 , φ 0 n = φ J n = 0 , v j 0 = v 0 x L + j h , φ j 0 = φ 0 x L + j h .
By Taylor series expansion, we have
E r j n = ( u t + ρ x u x x t + u u x ) ( x j , t n ) = O ( h 2 + τ ) , E s j n = ( ρ t + u x ) ( x j , t n ) = O ( h 2 + τ ) .
Next, based on Equations (2) and (3), a TT-M finite difference scheme for problem (1) is constructed with three steps.
Step 1: on the coarse time mesh, let u C , j k s = u x j , t k s , ρ C , j k s = ρ x j , t k s be the numerical solutions of of u x , t and ρ x , t in term of the point x j , t k s , then the coarse time nonlinear finite difference scheme is given as
( u C , j k s ) t + ( ρ C , j ( k + 1 ) s ) x ^ ( u C , j k s ) x x ¯ t + 1 3 u C , j k s + 1 2 ( u C , j k s + 1 2 ) x ^ + ( u C , j k s + 1 2 ) 2 x ^ = 0 ,
( ρ C , j k s ) t + ( u C , j k s ) x ^ = 0 ,
u C , 0 k s = u C , J k s = 0 , ρ C , 0 k s = ρ C , J k s = 0 , k = 0 , 1 , , P , u C , j 0 = u 0 x L + j h , ρ C , j 0 = ρ 0 x L + j h , j = 1 , 2 , , J 1 ,
where u C , j k s + 1 2 = 1 2 ( u C , j ( k + 1 ) s + u C , j k s ) .
Step 2: based on the solutions u C k s , ρ C k s at time levels t k s obtained from step 1, we apply the Lagrange’s linear interpolation formula to compute u C k s l , ρ C k s l at time levels t k s l ( l = 1 , 2 , , s 1 and k = 1 , 2 , , P , k s l = n ) , we have
u C k s l = t k s l t k s t ( k 1 ) s t k s u C ( k 1 ) s + t k s l t ( k 1 ) s t k s t ( k 1 ) s u C k s = l s u C ( k 1 ) s + ( 1 l s ) u C k s ,
ρ C k s l = t k s 1 t k s t ( k 1 ) s t k s ρ C ( k 1 ) s + t k s 1 t ( k 1 ) s t k s t ( k 1 ) s ρ C k s = l s ρ C ( k 1 ) s + ( 1 l s ) ρ C k s .
Remark 1.
The Equation (7) is only employed for theoretical analysis of the scheme. In numerical simulation, the coarse numerical solutions ρ C k s l are not needed for computation since they are not used in step 3.
Step 3: based on all the coarse numerical solutions u C , j n ( n = 0 , 1 , 2 , , P s = N , j = 1 , 2 , , J 1 ) obtained in the first two steps, Taylor’s formula is used to construct a linearized system on the fine time mesh, which is expressed as follows. Let u F , j n = u x j , t n , ρ F , j n = ρ x j , t n be the numerical solutions of u x , t and ρ x , t in term of the point x j , t n on the fine time mesh, then
( u F , j n ) t + ( ρ F , j n + 1 ) x ^ ( u F , j n ) x x ¯ t + 1 6 h f ( u C , j 1 n + 1 2 , u C , j n + 1 2 , u C , j + 1 n + 1 2 ) + f x ( u C , j 1 n + 1 2 , u C , j n + 1 2 , u C , j + 1 n + 1 2 ) ( u F , j 1 n + 1 2 u C , j 1 n + 1 2 ) + f y ( u C , j + 1 n + 1 2 , u C , j n + 1 2 , u C , j + 1 n + 1 2 ) ( u F , j n + 1 2 u C , j n + 1 2 ) + f z ( u C , j 1 n + 1 2 , u C , j n + 1 2 , u C , j + 1 n + 1 2 ) ( u F , j + 1 n + 1 2 u C , j + 1 n + 1 2 ) = 0 ,
( ρ F , j n ) t + ( u F , j n ) x ^ = 0 ,
u F , 0 n = u F , J n = 0 , ρ F , 0 n = ρ F , J n = 0 , u F , j 0 = u 0 ( x L + j h ) , ρ F , j 0 = ρ 0 ( x L + j h ) , j = 1 , , J 1 , n = 0 , 1 , 2 , , N ,
where f ( x , y , z ) = ( z x ) y + z 2 x 2 and
f x ( x , y , z ) = y 2 x , f y ( x , y , z ) = z x , f z ( x , y , z ) = y + 2 z
are the three partial derivatives of f ( x , y , z ) with respect to x , y , z .
Remark 2.
Our method, similarly to the Gauss-Seidel method applied to linear systems, has been modified in order to enhance the accuracy of fine mesh solutions u F n + 1 by using u F n in calculation.

4. Convergence Analysis and Stability of the TT-M Finite Difference Scheme

The focus of this section is on performing convergence analysis of the nonlinear system specifically on the coarse time mesh.
Theorem 1.
Suppose that the exact solutions v n , φ n to the initial boundary value problem Equation (1) is sufficiently smooth and let u C n , ρ C n be the numerical solutions on the coarse time mesh. Then,
v n u C n O ( h 2 + τ C ) , φ n ρ C n O ( h 2 + τ C ) .
Proof. 
Denote e C , j k s = v j k s u C , j k s , η C , j k s = φ j k s ρ C , j k s , 1 j J 1 , 0 k P . Subtracting Equation (4) from Equation (2) and Equation (5) from Equation (3), we obtain
E r C , j k s = ( e C , j k s ) t + ( η C , j ( k + 1 ) s ) x ^ ( e C , j k s ) x x ¯ t + 1 3 v j k s + 1 2 ( v j k s + 1 2 ) x ^ + ( v j k s + 1 2 ) 2 x ^ 1 3 u C , j k s + 1 2 ( u C , j k s + 1 2 ) x ^ + ( u C , j k s + 1 2 ) 2 x ^ ,
E s C , j k s = ( η C , j k s ) t + ( e C , j k s ) x ^ .
The proof contains two cases. Firstly, we consider the case of n = k s ( k = 0 , 1 , 2 , , P ) , then n + 1 = ( k + 1 ) s . The initial and boundary condition satisfies
e C , j 0 = 0 , η C , j 0 = 0 , u C , 0 n = u C , J n = 0 , ρ C , 0 n = ρ C , J n = 0 .
Taking the inner product (·, ·) on both sides of Equation (10) with e C n + 1 + e C n , we have
( E r C n , e C n + 1 + e C n ) = ( ( e C n + 1 ) t , e C n + e C n ) + ( ( η C n + 1 ) x ^ , e C n + 1 + e C n ) ( ( e C n ) x x ¯ t , e C k s + 1 + e C n ) + h j = 1 J 1 ( I + II ) ( e C , j n + 1 + e C , j n ) ,
where
I = 1 3 v j n + 1 2 ( v j n + 1 2 ) x ^ u C , j n + 1 2 ( u C , j n + 1 2 ) x ^ , II = 1 3 ( v j n + 1 2 ) 2 x ^ ( u C , j n + 1 2 ) 2 x ^ .
Notice that
( ( e C n ) t , e C n + 1 + e C n ) = 1 τ C ( e C n + 1 2 e C n 2 ) ,
( ( η C n + 1 ) x ^ , e C n + 1 + e C n ) = h j = 1 J 1 η C , j n + 1 ( e C , j n + 1 + e C , j n ) x ^ ,
( ( e C n ) x x ¯ t , e C n + 1 + e C n ) = 1 τ C ( e C , x n + 1 2 e C , x n 2 ) ,
h j = 1 J 1 I · ( e C , j n + 1 + e C , j n ) = 2 3 h j = 1 J 1 ( e C , j n + 1 2 · v j n + 1 2 ) x ^ e C , j n + 1 2 2 3 h j = 1 J 1 ( e C , j n + 1 2 ) x ^ ( e C , j + 1 n + 1 2 + e C , j 1 n + 1 2 ) u C , j n + 1 2 ,
h j = 1 J 1 II · ( e C , j n + 1 + e C , j n ) = 2 3 h j = 1 J 1 ( v j n + 1 2 · e C , j n + 1 2 ) x ^ e C , j n + 1 2 2 3 h j = 1 J 1 u C , j n + 1 2 e C , j n + 1 2 ( e C , j n + 1 2 ) x ^ ,
then substituting Equations (13)–(17) into (12), we have
e C n + 1 2 +   e C , x n + 1 2 =   e C n 2 +   e C , x n 2 + τ C h j = 1 J 1 η C , j n + 1 ( e C , j n + 1 + e C , j n ) x ^ + 2 3 τ C h j = 1 J 1 ( e C , j n + 1 2 ) x ^ ( e C , j + 1 n + 1 2 + e C , j 1 n + 1 2 ) u C , j n + 1 2 + u C , j n + 1 2 e C , j n + 1 2 ( e C , j n + 1 2 ) x ^ + τ C ( E r C n , e C n + 1 + e C n ) .
From Cauchy–Schwarz inequality, we obtain
( E r C n , e C n + 1 +   e C n )   E r C n 2 +   1 2 ( e C n + 1 2 +   e C n 2 ) .
Using Lemma 1, the Equation (18) can be rewritten as
e C n + 1 2 + e C , x n + 1 2 e C n 2 + e C , x n 2 + M τ C ( η C n + 1 2 + e C n + 1 2 + e C n 2 + e C , x n + 1 2 + e C , x n 2 ) + τ C E r C n 2 .
Similarly, taking the inner product (·, ·) on both sides of Equation (11) with η C n + 1 + η C n , we obtain
η C n + 1 2 = η C n 2 τ C h j = 1 J 1 ( e C , j n ) x ^ ( η C , j n + 1 + η C , j n ) + τ C ( E s C n , η C n + 1 + η C n ) .
From the Cauchy–Schwarz inequality, we have
( E s C n , η C n + 1 + η C n )   E s C n 2 +   1 2 ( η C n + 1 2 +   η C n 2 ) .
Using Lemma 1, the Equation (20) can be rewritten as
η C n + 1 2   η C n 2 +   M τ C ( e C , x n 2 +   η C n + 1 2 +   η C n 2 ) + τ C E s C n 2 .
Add Equations (19) and (21), we obtain
e C n + 1 2 + e C , x n + 1 2 + η C n + 1 2 e C n 2 + e C , x n 2 + η C n 2 + M τ C ( η C n + 1 2 + η C n 2 + e C n + 1 2 + e C n 2 + e C , x n + 1 2 + e C , x n 2 ) + τ C E r C n 2 + τ C E s C n 2 .
Let B C n = e C n 2 + e C , x n 2 + η C n 2 , then Equation (22) becomes
B C n + 1 B C n M τ C ( B C n + 1 + B C n ) + M τ C ( h 2 + τ C ) 2 ,
and obtain
( 1 M τ C ) ( B C n + 1 B C n ) 2 M τ C B C n + M τ C ( h 2 + τ C ) 2 .
By taking τ C small enough so that ( 1 M τ C ) > λ > 0 , then
B C n + 1 B C n M τ C B C n + M τ C ( h 2 + τ C ) 2 .
Summing from 0 to P 1 inequalities in Equation (23), we have
B C P B C 0 M τ C n = 1 P 1 B C n + M ( h 2 + τ C ) 2 ,
and using Lemma 2, we obtain
B C P [ B C 0 + M ( h 2 + τ C ) 2 ] e M P τ C .
From Equation (24) and the initial and boundary condition, we have
e C n < O ( h 2 + τ C ) , e C , x n < O ( h 2 + τ C ) , η C n < O ( h 2 + τ C ) .
Then using Lemma 3, we obtain
e C n < O ( h 2 + τ C ) .
Secondly, we consider the case of n = k s l , ( l = 1 , 2 , , s 1 and k = 1 , 2 , , P , k s l = n ) . Based on the Lagrange’s interpolation formula, we obtain
v k s l = t k s l t k s t ( k 1 ) s t k s v ( k 1 ) s + t k s l t ( k 1 ) s t k s t ( k 1 ) s v k s = l s v ( k 1 ) s + ( 1 l s ) v k s + v ( θ 1 ) 2 ( t t ( k 1 ) s ) ( t t k s ) , θ 1 ( t ( k 1 ) s , t k s ) ,
φ k s l = t k s l t k s t ( k 1 ) s t k s φ ( k 1 ) s + t k s l t ( k 1 ) s t k s t ( k 1 ) s φ k s = l s φ ( k 1 ) s + ( 1 l s ) φ k s + φ ( θ 2 ) 2 ( t t ( k 1 ) s ) ( t t k s ) , θ 2 ( t ( k 1 ) s , t k s ) .
Subtracting Equation (27) from (6), we have
v k s l u C k s l = l s ( v ( k 1 ) s u C ( k 1 ) s ) + ( 1 l s ) ( v k s u C k s ) + v ( θ 1 ) 2 ( t t ( k 1 ) s ) ( t t k s ) .
Subtracting Equation (28) from (7), we obtain
φ k s l ρ C k s l = l s ( φ ( k 1 ) s ρ C ( k 1 ) s ) + ( 1 l s ) ( φ k s ρ C k s ) + φ ( θ 2 ) 2 ( t t ( k 1 ) s ) ( t t k s ) .
Using (25), (26) and triangle inequality, we conclude
e C k s l O ( h 2 + τ C ) , η C k s l O ( h 2 + τ C ) .
We obtain the result of Theorem 1 by synthesizing the aforementioned two cases. □
Next, we give the convergence analysis of the scheme on the fine time mesh.
Theorem 2.
Suppose that the exact solutions v n , φ n to the initial boundary value problem Equation (1) is sufficiently smooth and let u F n , ρ F n be the numerical solutions on the fine time mesh. Then,
v n u F n O ( h 2 + τ C 2 + τ F ) , φ n ρ F n O ( h 2 + τ C 2 + τ F ) .
Proof. 
Assume e F , j n = v j n u F , j n , η F , j n = φ j n ρ F , j n , 1 j J 1 , 0 n N , Subtracting Equation (8) from Equation (2) and Equation (9) from Equation (3), we obtain
E r F , j n = ( e F , j n ) t + ( η F , j n + 1 ) x ^ ( e F , j n ) x x ¯ t + 1 6 h [ ( f x e F , j 1 n + 1 2 + f y e F , j n + 1 2 + f z e F , j + 1 n + 1 2 ) + Q ] ,
E s F , j n = ( η F , j n ) t + ( e F , j n ) x ^ ,
where
Q = 1 2 [ ( f x x ( e C , j 1 n + 1 2 ) 2 + f y y ( e C , j n + 1 2 ) 2 + f z z ( e C , j + 1 n + 1 2 ) 2 ] + [ f x y e C , j 1 n + 1 2 e C , j n + 1 2 + f x z e C , j 1 n + 1 2 e C , j + 1 n + 1 2 + f y z e C , j n + 1 2 e C , j + 1 n + 1 2 ] ,
and f x x = f x x ( ξ , ε , δ ) , f y y = f y y ( ξ , ε , δ ) , f z z = f z z ( ξ , ε , δ ) , f x y = f x y ( ξ , ε , δ ) , f x z = f x z ( ξ , ε , δ ) , f y z = f y z ( ξ , ε , δ ) are the second order partial derivatives of f ( x , y , z ) , ξ ( v j 1 n , u C , j 1 n ) , ε ( v j n , u C , j n ) , δ ( v j + 1 n , u C , j + 1 n ) .
Taking the inner product (·, ·) on both sides of Equation (29) with e F n + 1 + e F n , we have
e F n + 1 2 + e F , x n + 1 2 = e F n 2 + e F , x n 2 + τ F h j = 1 J 1 η F , j n + 1 ( e F , j n + 1 ) x ^ + ( e F , j n ) x ^ τ F 3 j = 1 J 1 ( f x e F , j 1 n + 1 2 + f y e F , j n + 1 2 + f z e F , j + 1 n + 1 2 ) e F , j n + 1 2 τ F 3 j = 1 J 1 Q e F , j n + 1 2 + 2 τ F ( E r F n , e F n + 1 2 ) .
Using f y = 1 2 ( f x + f z ) and f x x = 2 , f y y = 0 , f z z = 2 , f x y = 1 , f x z = 0 , f y z = 1 , we obtain
j = 1 J 1 ( f x e F , j 1 n + 1 2 + f y e F , j n + 1 2 + f z e F , j + 1 n + 1 2 ) e F , j n + 1 2 = ( f x e F , x ¯ n + 1 2 , e F n + 1 2 ) + 3 h ( f y e F n + 1 2 , e F n + 1 2 ) + ( f z e F , x n + 1 2 , e F n + 1 2 ) ,
j = 1 J 1 Q e F , j n + 1 2 = 1 2 j = 1 J 1 [ f x x ( e C , j 1 n + 1 2 ) 2 + f y y ( e C , j n + 1 2 ) 2 + f z z ( e C , j + 1 n + 1 2 ) 2 ] e F , j n + 1 2 + j = 1 J 1 [ f x y e C , j 1 n + 1 2 e C , j n + 1 2 + f x z e C , j 1 n + 1 2 e C , j + 1 n + 1 2 + f y z e C , j n + 1 2 e C , j + 1 n + 1 2 ] e F , j n + 1 2 = ( ( e C n + 1 2 ) x 2 , e F n + 1 2 ) ( ( e C n + 1 2 ) 2 , e F , x n + 1 2 ) + ( e C , x n + 1 2 , e C n + 1 2 e F n + 1 2 ) ( e C n + 1 2 , ( e C n + 1 2 e F n + 1 2 ) x ) .
Using Lemma 1 and the Cauchy–Schwarz inequality, we have
τ F 3 ( f x e F , x ¯ n + 1 2 , e F n + 1 2 ) τ F h ( f y e F n + 1 2 , e F n + 1 2 ) τ F 3 ( f z e F , x n + 1 2 , e F n + 1 2 ) M τ F ( e F , x n + 1 2 2 + e F n + 1 2 2 ) ,
τ F 3 ( ( e C n + 1 2 ) x 2 , e F n + 1 2 ) + τ F 3 ( ( e C n + 1 2 ) 2 , e F , x n + 1 2 ) τ F 3 ( e C , x n + 1 2 , e C n + 1 2 e F n + 1 2 ) + τ F 3 ( e C n + 1 2 , ( e C n + 1 2 e F n + 1 2 ) x ) M τ F ( e C n + 1 2 2 e C n + 1 2 2 + e C n + 1 2 2 e C , x n + 1 2 2 + e F n + 1 2 2 + e F , x n + 1 2 2 ) ,
2 τ F ( E r F n , e F n + 1 2 ) τ F E r F n 2 + M τ F e F n + 1 2 2 .
Substituting Equations (34) and (35) into (31), then
e F n + 1 2 + e F , x n + 1 2 e F n 2 + e F , x n 2 + M τ F ( η F n + 1 2 + e F n + 1 2 + e F n 2 + e F , x n + 1 2 + e F , x n 2 ) + M τ F ( e C n + 1 2 2 e C n + 1 2 2 + e C n + 1 2 2 e C , x n + 1 2 2 ) + τ F E r F n 2 .
Taking the inner product (·, ·) on both sides of Equation (30) with η F n + 1 + η F n , we obtain
η F n + 1 2 η F n 2 + M τ F ( e F , x n 2 + η F n + 1 2 + η F n 2 ) + τ F E s F n 2 .
Add Equations (37) and (38), we have
e F n + 1 2 + e F , x n + 1 2 + η F n + 1 2 e F n 2 + e F , x n 2 + η F n 2 + M τ F ( η F n + 1 2 + η F n 2 + e F n + 1 2 + e F n 2 + e F , x n + 1 2 + e F , x n 2 ) + M τ F ( e C n + 1 2 2 e C n + 1 2 2 + e C n + 1 2 2 e C , x n + 1 2 2 ) + τ F E r F n 2 + τ F E s F n 2 .
Let B F n = e F n 2 + e F , x n 2 + η F n 2 , then
B F n + 1 B F n M τ F ( B F n + 1 + B F n ) + M τ F ( e C n + 1 2 2 e C n + 1 2 2 + e C n + 1 2 2 e C , x n + 1 2 2 ) + τ F E r F n 2 + τ F E s F n 2 ,
and obtain
( 1 M τ F ) ( B F n + 1 B F n ) 2 M τ F B F n + M τ F ( h 4 + τ C 4 + τ F 2 ) .
By taking τ F small enough so that ( 1 M τ F ) > λ > 0 , then
B F n + 1 B F n M τ F ( h 4 + τ C 4 + τ F 2 ) + M τ F B F n .
Summing from 0 to N 1 inequalities in Equation (40), we obtain
B F N B F 0 + M ( h 4 + τ C 4 + τ F 2 ) + M τ F n = 0 N 1 B F n .
Using Lemma 2, we obtain
B F N [ B F 0 + M ( h 4 + τ C 4 + τ F 2 ) ] e M N τ F .
From Equation (42) and the initial and boundary condition, we have
e F n O ( h 2 + τ C 2 + τ F ) , e F , x n O ( h 2 + τ C 2 + τ F ) , η F n < O ( h 2 + τ C 2 + τ F ) .
Using Lemma 3, it led to
e F n O ( h 2 + τ C 2 + τ F ) .
This completes the proof of Theorem 2. □
Next, we prove stability of the scheme on the coarse time mesh.
Theorem 3.
Suppose that u 0 H h , 0 x L , x R , if τ C and h are sufficiently small, then the proposed scheme (4)–(5) is stable with respect to the initial conditions.
Proof. 
The proof contains two cases. Firstly, we consider the case of n = k s ( k = 0 , 1 , 2 , , P ) . Taking the inner product (·, ·) on both sides of Equation (4) with u C n + 1 + u C n , we have
u C n + 1 2 u C n 2 + u C , x n + 1 2 u C , x n 2 + τ C h j = 1 J 1 [ ρ C , j n + 1 ( u C , j n + 1 ) x ^ ρ C , j n + 1 ( u C , j n ) x ^ ] = 0 .
Taking the inner product (·, ·) on both sides of Equation (5) with ρ C n + 1 + ρ C n , we obtain
ρ C n + 1 2 | ρ C n 2 + τ C h j = 1 J 1 [ ( u C , j n ) x ^ ρ C , j n + 1 + ( u C , j n ) x ^ ρ C , j n ] = 0 .
Add Equations (45) and (46), we obtain
u C n + 1 2 + u C , x n + 1 2 + ρ C n + 1 2 = u C n 2 + u C , x n 2 + ρ C n 2 + τ C h j = 1 J 1 ρ C , j n + 1 ( u C , j n + 1 ) x ^ τ C h j = 1 J 1 ( u C , j n ) x ^ ρ C , j n u C n 2 + u C , x n 2 + | ρ C n 2 + M τ C ( ρ C n + 1 2 + ρ C n 2 + u C n + 1 2 + u C n 2 + u C , x n + 1 2 + u C , x n 2 ) .
Let D C n = u C n 2 + u C , x n 2 + ρ C n 2 , then Equation (47) becomes
D C n + 1 D C n M τ C ( D C n + 1 + D C n ) ,
and obtain
( 1 M τ C ) ( D C n + 1 D C n ) 2 M τ C D C n .
By taking τ C small enough so that 1 M τ C > λ > 0 , then
D C n + 1 D C n M τ C D C n .
Summing from 0 to P 1 inequalities in Equation (48), we have
D C N D C 0 + M τ C n = 0 N 1 D C n ,
and using Lemma 2, we obtain
D C N D C 0 e M T ,
that is
u C n 2 + u C , x n 2 + ρ C n 2 M ( u C 0 2 + u C , x 0 2 + ρ C 0 2 ) .
Secondly, we consider the case of n = k s l , ( l = 1 , 2 , , s 1 and k = 1 , 2 , , P , k s l = n ) . From Equations (6) and (7), we obtain
u C k s l l s u C ( k 1 ) s + ( 1 l s ) u C k s M ( u C 0 + u C , x 0 + ρ C 0 ) , ρ C k s l l s ρ C ( k 1 ) s + ( 1 l s ) ρ C k s M ( u C 0 + u C , x 0 + ρ C 0 ) .
The results obtained in the aforementioned two cases mean that the solution of scheme (4)–(5) is stable on the coarse time mesh. This completes the proof of Theorem 3. □
Next, we present stability analysis of the scheme on the fine time mesh.
Theorem 4.
Suppose that u 0 H h , 0 x L , x R , if τ F and h are sufficiently small, then the proposed scheme (8)–(9) is stable with respect to the initial conditions.
Proof. 
Taking the inner product (·, ·) on both sides of Equation (8) with u F n + 1 + u F n , we have
1 τ F ( u F n + 1 2 + u F , x n + 1 2 ) = 1 τ F ( u F n 2 + u F , x n 2 ) h j = 1 J 1 [ ρ F , j n + 1 ( u F , j n + 1 ) x ^ ρ F , j n + 1 ( u F , j n ) x ^ ] 1 3 j = 1 J 1 [ ( u C , j + 1 n + 1 2 u C , j 1 n + 1 2 ) u C , j n + 1 2 + ( u C , j + 1 n + 1 2 ) 2 ( u C , j 1 n + 1 2 ) 2 ] u F , j n + 1 2 1 3 j = 1 J 1 ( f x u F , j 1 n + 1 2 + f y u F , j n + 1 2 + f z u F , j + 1 n + 1 2 ) u F , j n + 1 2 + 1 3 j = 1 J 1 ( f x u C , j 1 n + 1 2 + f y u C , j n + 1 2 + f z u C , j + 1 n + 1 2 ) u F , j n + 1 2 .
From the Cauchy–Schwarz inequality, we obtain
1 3 j = 1 J 1 [ ( u C , j + 1 n + 1 2 u C , j 1 n + 1 2 ) u C , j n + 1 2 + ( u C , j + 1 n + 1 2 ) 2 ( u C , j 1 n + 1 2 ) 2 ] u F , j n + 1 2 M ( u C n + 1 2 2 u C n + 1 2 2 + u C n + 1 2 2 u C , x n + 1 2 2 + u F n + 1 2 2 + u F , x n + 1 2 2 ) ,
1 3 j = 1 J 1 ( f x u F , j 1 n + 1 2 + f y u F , j n + 1 2 + f z u F , j + 1 n + 1 2 ) u F , j n + 1 2 M ( u F n + 1 2 2 + u F , x n + 1 2 2 ) ,
1 3 j = 1 J 1 ( f x u C , j 1 n + 1 2 + f y u C , j n + 1 2 + f z u C , j + 1 n + 1 2 ) u F , j n + 1 2 M ( u C n + 1 2 2 + u C , x n + 1 2 2 + u F n + 1 2 2 + u F , x n + 1 2 2 ) ,
then substituting Equations (51)–(53) into (50), we have
u F n + 1 2 + u F , x n + 1 2 u F n 2 + u F , x n 2 + τ F h j = 1 J 1 [ ρ F , j n + 1 ( u F , j n + 1 ) x ^ + ρ F , j n + 1 ( u F , j n ) x ^ ] + M τ F ( u C n + 1 2 2 u C n + 1 2 2 + u C n + 1 2 2 u C , x n + 1 2 2 + u C n + 1 2 2 + u C , x n + 1 2 2 ) + M τ F ( u F n + 1 2 2 + u F , x n + 1 2 2 ) .
Taking the inner product (·, ·) on both sides of Equation (9) with ρ F n + 1 + ρ F n , we obtain
ρ F n + 1 2 = ρ F n 2 τ F h j = 1 J 1 [ ( u F , j n ) x ^ ρ F , j n + 1 + ( u F , j n ) x ^ ρ F , j n ] .
Add Equations (54) and (55), we obtain
u F n + 1 2 + u F , x n + 1 2 + ρ F n + 1 2 u F n 2 + u F , x n 2 + ρ F n 2 + M τ F ( u F n + 1 2 + u F , x n 2 + u F , x n + 1 2 + u F , x n 2 + ρ F n + 1 2 + ρ F n 2 ) + M τ F ( u C n + 1 2 2 u C n + 1 2 2 + u C n + 1 2 2 u C , x n + 1 2 2 + u C n + 1 2 2 + u C , x n + 1 2 2 ) .
Let D F n = u F n 2 + u F , x n 2 + ρ F n 2 , then combine with the result of Theorem 3, the Equation (56) becomes
D F n + 1 D F n M τ F ( D F n + 1 + D F n ) + M τ F ( u C 0 2 + u C , x 0 2 ) ,
and obtain
1 M τ F ( D F n + 1 D F n ) 2 M τ F D F n + M τ F ( u C 0 2 + u C , x 0 2 ) .
By taking τ F small enough so that 1 M τ F > λ > 0 , then
D F n + 1 D F n M τ F D F n + M τ F ( u C 0 2 + u C , x 0 2 ) .
Summing from 0 to N 1 inequalities in Equation (57), we obtain
n = 0 N 1 ( D F n + 1 D F n ) M τ F n = 0 N 1 D F n + M ( u C 0 2 + u C , x 0 2 ) ,
and then
D F N D F 0 + M τ F n = 0 N 1 D F n + M ( u C 0 2 + u C , x 0 2 ) .
Using Lemma 2, we obtain
D F N [ D F 0 + M u C 0 2 + u C , x 0 2 ) ] e M T ,
that is
u F n 2 + u F , x n 2 + ρ F n 2 M ( u F 0 2 + u F , x 0 2 + ρ F 0 2 + u C 0 2 + u C , x 0 2 ) .
Notice that u C 0 2 = u F 0 2 , u C , x 0 2 = u F , x 0 2 , then
u F n 2 + u F , x n 2 + ρ F n 2 M ( u F 0 2 + u F , x 0 2 + ρ F 0 2 ) ,
which indicates that the solution of scheme (8)–(9) is stable on the fine time mesh. This completes the proof of Theorem 4. □

5. Numerical Results

This section provides some numerical examples aimed at demonstrating the accuracy and computational time of the TT-M finite difference scheme that was discussed in Section 3. All simulations were implemented using a Intel(R) i7-10710U 1.61 GHz CPU and 16 GB memory personal computer running Windows 10 and Matlab R2019b. The SRLW equation is represented in the following formation:
u t + ρ x + u u x u x x t = 0 , 40 x 40 , 0 < t 4 , ρ t + u x = 0 , 40 x 40 , 0 < t 4 , u x L , t = u x R , t = 0 , ρ x L , t = ρ x R , t = 0 , 0 < t 4 , u ( x , 0 ) = u 0 ( x ) , ρ ( x , 0 ) = ρ 0 ( x ) , 40 x 40 .
and consider the following initial conditions
u 0 ( x ) = 5 2 sech 2 5 6 x , ρ 0 ( x ) = 5 3 sech 2 5 6 x .
The exact solitary wave solution [4] of the SRLW Equation (1) has the following form
u ( x , t ) = 3 v 2 1 v sech 2 v 2 1 4 v 2 ( x v t ) , ρ ( x , t ) = 3 v 2 1 v 2 sech 2 v 2 1 4 v 2 ( x v t ) ,
which takes the form of
u ( x , t ) = 5 2 sech 2 5 6 x 3 2 t , ρ ( x , t ) = 5 3 sech 2 5 6 x 3 2 t ,
when we set v = 1.5 .
Next, we presented some values of error, convergence rate, conservation laws and long time simulation of the proposed scheme.

5.1. Error and Convergence Rate

We define the error and convergence rate by the following formula:
e m ( h , τ ) = v n u m n , η m ( h , τ ) = φ n ρ m n , u R a t e m x = log 2 e m ( 2 h , 4 τ ) e m ( h , τ ) , ρ R a t e m x = log 2 η m ( 2 h , 4 τ ) η m ( h , τ ) , u R a t e m t = log 2 e m ( 2 h , 2 τ ) e m ( h , τ ) , ρ R a t e m t = log 2 η m ( 2 h , 2 τ ) η m ( h , τ ) ,
where m represents the TT-M finite difference scheme or the standard nonlinear finite difference (SNFD) scheme in [22]. We set τ C = 4 τ F in the whole numerical illustration process.
In order to illustrate the advantage of this approach, we calculated the error, convergence rate and computational time of the presented scheme and compared these numerical results with those obtained from scheme in [22]. Table 1 and Table 2 present discrete norm errors, the corresponding convergence rates and the time cost for both the TT-M finite difference scheme and the SNFD scheme in [22] under τ F = h 2 and h = τ F , respectively. We recorded the error of both schemes at the final time t = 4 for various mesh steps. The results from the new scheme show nearly identical important values as those obtained by the SNFD scheme in [22]. In term of the convergence rate, the results indicate that both the TT-M finite difference scheme and the SNFD scheme achieve approximately second-order convergence in space when τ F = h 2 and first-order in time when h = τ F , which confirming the theoretical results in Theorem 2. Especially, the proposed scheme achieves a remarkable reduction in computation time of over 30 percent as compared to the SNFD scheme under τ F = h 2 and nearly half the computational time of the SNFD scheme under h = τ F . The obtained results demonstrate a clear improvement achieved by our scheme over the previous method reported by [22].
The 3D images of numerical solutions of u ( x , t ) and ρ ( x , t ) computed by the present scheme were shown in Figure 1. From the graph, one can see the propagation state of waves over a period of time. Figure 2 illustrates the exact and numerical solutions of u ( x , t ) and ρ ( x , t ) at t = 4 obtained by the present scheme. It is evident that our numerical solutions exhibit an excellent correspondence with the exact solution. Furthermore, the CPU times of the two schemes are plotted in Figure 3 under τ F = h 2 and h = τ F , respectively. From the figure, it is also can be seen that the present scheme can significantly decrease computation time. To sum up, compared with the method in [22], the scheme proposed in this paper not only ensures error and convergence order, but also greatly reduces computational time.

5.2. Conservative Approximations

To further verify the accuracy of the new scheme, we calculate four conservation laws of the SRLW Equation (1), such as:
Q 1 = 1 2 ρ d x , Q 2 = 1 2 u d x , Q 3 = 1 2 ( u 2 + u x 2 + ρ 2 ) d x , Q 4 = 1 2 ( ρ u + 1 6 u 3 ) d x .
Afterwards, employing discrete forms, we are able to compute four approximate conservative quantities which can be represented as
Q 1 = h 2 j = 0 J 1 ρ j n , Q 2 = h 2 j = 0 J 1 u j n , Q 3 = h 2 j = 0 J 1 ( u j n ) 2 + 1 2 h j = 0 J 1 ( u j + 1 n u j n ) 2 + h 2 j = 0 J 1 ( ρ j n ) 2 , Q 4 = h 2 j = 0 J 1 ρ j n u j n + h 12 j = 0 J 1 ( u j n ) 3 .
The values of four quantities are recorded in Table 3, Table 4, Table 5 and Table 6. In Table 3 and Table 4, regardless of the time step and grid spacing, the quantities Q 1 and Q 2 remain well-preserved at various times. In Table 5, for the case h = 1 / 2 and τ F = 1 / 4 , one can see that the quantity Q 3 experiences a slight increase as time increases, however, as the spatial and temporal step sizes decrease, the variation of Q 3 becomes extremely small. In Table 6, it has been found that for quantity Q 4 , there was a minor decline under different mesh steps, but it gradually rebounded over time. Meanwhile, as the spatial and temporal step sizes decrease, the Q 4 increases slightly. Figure 4 plots the variation curves of four quantities for the case h = 1 / 8 and τ F = 1 / 64 , which visually demonstrate that our scheme preserves the four conservation laws.

5.3. Long-Time Simulation

Using the parameters h = 0.1 , τ F = 0.01 , x L = 20 , x R = 100 and T = 40 , the long time waveforms of u ( x , t ) and ρ ( x , t ) generated by the present scheme are illustrated in Figure 5. The waveforms at t = 20 and 40 demonstrate a significant level of agreement with the waveforms at t = 0, which also indicates the accuracy of the scheme. Figure 6 displays the computational times of both the SNFD scheme and proposed scheme at long time. From the graph, it can be seen that the computational time difference between the two schemes grows increasingly as the simulation time increases, which indicates that the longer the simulation time, the more computational time can be saved by the scheme proposed in this paper. This once again proves the advantage of our scheme compared to the SNFD scheme in [22].

6. Conclusions

This paper presents a new time two-mesh finite difference scheme for solving the symmetric regularized long wave equation, which contains a nonlinear derivative term in its formulation. Based on the previous work, our aim is to obtain a scheme that is not only error-preserving but also can save more computational time compared to the scheme in [22]. To achieve this goal, we have constructed a TT-M finite difference scheme that mainly consists of three steps. Firstly, the time interval is divided into coarse and fine time meshes, then the nonlinear system is solved on the coarse time mesh; secondly, coarse numerical solutions on the fine time mesh are computed using an interpolation formula based on the solutions derived in the step one; lastly, the nonlinear term of the SRLW equation is linearized using Taylor’s formula for a function with three variables and constructed the TT-M finite difference linear numerical scheme on the fine time mesh. In terms of theoretical results, we investigated the convergence and stability of numerical solutions. We observed that the solutions obtained from our proposed scheme converged to the exact solutions of the SRLW equation, and the rate of convergence is O ( τ C 2 + τ F + h 2 ) . Additionally, τ F and h are sufficiently small, then the proposed scheme (8)–(9) is stable with respect to the initial conditions on the fine time mesh. The numerical results demonstrated that the proposed method effectively supported the analysis of the convergence rate. Meanwhile, one can observe that the errors of u ( x , t ) and ρ ( x , t ) in both the SNFD scheme in [22] and the present scheme are nearly identical, but the proposed scheme can reduce the computational time compared with the SNFD scheme. In addition, we calculated four conservation laws of the SRLW equation, regardless of the time step and grid spacing; the four quantities remain well-preserved at various times. The long time simulation example of the proposed scheme is also conducted in this part, the result shows that the waveform always maintains its original state in a long time. The error increases gently as simulation time is extended, which is a limitation of the proposed scheme and it will be improved by other methods in our future work. These findings demonstrate that the proposed method is more effective and yields a better improvement in solving the SRLW equation.

Author Contributions

Conceptualization, J.G.; methodology, J.G.; software, J.G. and S.H.; validation, S.H., Q.B. and J.L.; formal analysis, J.G. and S.H.; writing—original draft preparation, J.G.; writing—review and editing, J.G. and S.H.; funding acquisition, S.H., Q.B and J.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by National Natural Science Foundation of China (Nos. 12161034 and 11961022), Natural Science Foundation of Inner Mongolia (No. 2021MS01017), Natural Scientific Research Innovation Team of Hohhot Minzu College (No. HMTD202005).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All the data were computed using our method.

Acknowledgments

We are grateful to the anonymous reviewers for their valuable suggestions and comments.

Conflicts of Interest

The author declare no conflict of interest.

References

  1. Peregrine, D.H. Calculations of the development of an undular bore. J. Fluid Mech. 1966, 25, 321–330. [Google Scholar] [CrossRef]
  2. Peregrine, D.H. Long waves on beach. J. Fluid Mech. 1967, 27, 815–827. [Google Scholar] [CrossRef]
  3. Korteweg, D.J.; de Vries, G. On the change of form of long waves advancing in a rectangular canal, and on a new type of long stationary waves. Philos. Mag. 1895, 39, 422–443. [Google Scholar] [CrossRef]
  4. Seyler, C.E.; Fenstermacher, D.L. A symmetric regularized-long-wave equation. Phys. Fluids 1984, 27, 4–7. [Google Scholar] [CrossRef]
  5. Fang, S.; Guo, B.; Qiu, H. The existence of global attractors for a system of multi-dimensional symmetric regularized wave equations. Commun. Nonlinear Sci. Numer. Simul. 2009, 14, 61–68. [Google Scholar]
  6. Xu, F. Application of Exp-function method to symmetric regularized long wave (SRLW) equation. Phys. Lett. A 2008, 372, 252–257. [Google Scholar] [CrossRef]
  7. Abazari, R. Application of (G/G)-expansion method to travelling wave solutions of three nonlinear evolution equation. Comput. Fluids 2010, 39, 1957–1963. [Google Scholar] [CrossRef]
  8. Chand, F.; Malik, A.K. Exact traveling wave solutions of some nonlinear equations using(G/G)-expansion method. Int. J. Nonlinear Sci. 2012, 14, 416–424. [Google Scholar]
  9. Yalçınkaya, İ.; Ahmad, H.; Tasbozan, O.; Kurt, A. Soliton solutions for time fractional ocean engineering models with Beta derivative. J. Ocean Eng. Sci. 2022, 7, 444–448. [Google Scholar] [CrossRef]
  10. Mamun, A.A.; Ananna, S.N.; An, T.; Shahen, N.H.M.; Asaduzzaman, M. Dynamical behaviour of travelling wave solutions to the conformable time-fractional modified Liouville and mRLW equations in water wave mechanics. Heliyon 2021, 7, e07704. [Google Scholar] [CrossRef]
  11. Hussain, A.; Jhangeer, A.; Abbas, N.; Khan, I.; Nisar, K.S. Solitary wave patterns and conservation laws of fourth-order nonlinear symmetric regularized long-wave equation arising in plasma. Ain Shams Eng. J. 2021, 12, 3919–3930. [Google Scholar] [CrossRef]
  12. Lu, D.; Seadawy, A.R.; Ali, A. Applications of exact traveling wave solutions of modified Liouville and the symmetric regularized long wave equations via two new techniques. Results Phys. 2018, 9, 1403–1410. [Google Scholar] [CrossRef]
  13. Manafian, J.; Zamanpour, I. Exact travelling wave solutions of the symmetric regularized long wave (SRLW) using analytical methods. Stat. Optim. Inf. Comput. 2014, 2014, 47–55. [Google Scholar]
  14. Sendi, C.T.; Manafian, J.; Mobasseri, H.; Mirzazadeh, M.; Zhou, Q.; Bekir, A. Application of the ITEM for solving three nonlinear evolution equations arising in fluid mechanics. Nonlinear Dyn. 2019, 95, 669–684. [Google Scholar] [CrossRef]
  15. Bekir, A. New solitons and periodic wave solutions for some nonlinear physical models by using the sine–cosine method. Phys. Scr. 2008, 77, 045008. [Google Scholar] [CrossRef]
  16. Guo, B. The spectral method for symmetric regularized wave equations. J. Comput. Math. 1987, 5, 297–306. [Google Scholar]
  17. Zheng, J.; Zhang, R.; Guo, B. The Fourier pseudo-spectral method for the SRLW equation. Appl. Math. Mech. 1989, 10, 801–810. [Google Scholar]
  18. Wang, T.C.; Zhang, L.M.; Chen, F.Q. Conservative schemes for the symmetric regularized long wave equations. Appl. Math. Comput. 2007, 190, 1063–1080. [Google Scholar] [CrossRef]
  19. Yimnet, S.; Wongsaijai, B.; Rojsiraphisal, T.; Poochinapan, K. Numerical implementation for solving the symmetric regularized long wave equation. Appl. Math. Comput. 2016, 273, 809–825. [Google Scholar] [CrossRef]
  20. Hu, J.S.; Zheng, K.L.; Zheng, M.B. Numerical simulation and convergence analysis of high-order conservative difference scheme for SRLW equation. Appl. Math. Model. 2014, 38, 5573–5581. [Google Scholar] [CrossRef]
  21. Li, S.G. Numerical study of a conservative weighted compact difference scheme for the symmetric regularized long wave equations. Numer. Methods Partial. Differ. Equ. 2019, 35, 60–83. [Google Scholar] [CrossRef] [Green Version]
  22. Bai, Y.; Zhang, L.M. A conservative finite difference scheme for symmetric regularized long wave equations. Acta Math. Appl. Sin. 2007, 30, 248–255. [Google Scholar]
  23. Xu, Y.C.; Hu, B.; Xie, X.P.; Hu, J.S. Mixed finite element analysis for dissipative SRLW equations with damping term. Phys. Fluids 2012, 218, 4788–4797. [Google Scholar] [CrossRef]
  24. He, Y.Y.; Wang, X.F.; Cheng, H.; Deng, Y.Q. Numerical analysis of a high-order accurate compact finite difference scheme for the SRLW equation. Appl. Math. Comput. 2022, 418, 126837. [Google Scholar] [CrossRef]
  25. Liu, Y.; Yu, Z.D.; Li, H.; Liu, F.W.; Wang, J.F. Time two-mesh algorithm combined with finite element method for time fractional water wave model. Int. J. Heat Mass Transf. 2018, 120, 1132–1145. [Google Scholar] [CrossRef]
  26. Yin, B.L.; Liu, Y.; Li, H.; He, S. Fast algorithm based on TT-M FE system for space fractional Allen-Cahn equations with smooth and non-smooth solutions. J. Comput. Phys. 2019, 379, 351–372. [Google Scholar] [CrossRef] [Green Version]
  27. Liu, Y.; Fan, E.Y.; Yin, B.L.; Li, H.; Wang, J.F. TT-M finite element algorithm for a two-dimensional space fractional Gray-Scott model. Comput. Math. Appl. 2020, 80, 1793–1809. [Google Scholar] [CrossRef]
  28. Wen, C.; Liu, Y.; Yin, B.L.; Li, H.; Wang, J.F. Fast second-order time two-mesh mixed finite element method for a nonlinear distributed-order sub-diffusion model. Numer. Algorithms 2021, 88, 523–553. [Google Scholar] [CrossRef]
  29. Tian, J.L.; Sun, Z.Y.; Liu, Y.; Li, H. TT-M Finite Element Algorithm for the Coupled Schrödinger–Boussinesq Equations. Axioms 2022, 11, 314. [Google Scholar] [CrossRef]
  30. Qiu, W.L.; Xu, D.; Guo, J.; Zhou, J. A time two-grid algorithm based on finite difference method for the two-dimensional nonlinear time-fractional mobile/immobile transport model. Numer. Algorithms 2020, 85, 39–58. [Google Scholar] [CrossRef]
  31. Xu, D.; Guo, J.; Qiu, W.L. Time two-grid algorithm based on finite difference method for two-dimensional nonlinear fractional evolution equations. Appl. Numer. Math. 2019, 152, 169–184. [Google Scholar] [CrossRef]
  32. Niu, Y.X.; Liu, Y.; Li, H.; Liu, F.W. Fast high-order compact difference scheme for the nonlinear distributed-order fractional Sobolev model appearing in porous media. Math. Comput. Simul. 2023, 203, 387–407. [Google Scholar] [CrossRef]
  33. He, S.; Liu, Y.; Li, H. A time two-mesh compact difference method for the one-dimensional nonlinear schrödinger equation. Entropy 2022, 24, 806. [Google Scholar] [CrossRef] [PubMed]
  34. Zhou, Y.L. Application of Discrete Functional Analysis to the Finite Difference Method; International Academic Publishers: Beijing, China, 1990. [Google Scholar]
Figure 1. Three−dimensional images of u ( x , t ) (a) and ρ ( x , t ) (b) with h = 1 / 8 , τ F = 1 / 64 .
Figure 1. Three−dimensional images of u ( x , t ) (a) and ρ ( x , t ) (b) with h = 1 / 8 , τ F = 1 / 64 .
Fractalfract 07 00487 g001
Figure 2. Exact and numerical solution of u ( x , t ) (a) and ρ ( x , t ) (b) at t = 4 with h = 1 / 8 , τ F = 1 / 64 .
Figure 2. Exact and numerical solution of u ( x , t ) (a) and ρ ( x , t ) (b) at t = 4 with h = 1 / 8 , τ F = 1 / 64 .
Fractalfract 07 00487 g002
Figure 3. Comparison of CPU times with τ F = h 2 (a) and h = τ F (b) based on the data in Table 1 and Table 2.
Figure 3. Comparison of CPU times with τ F = h 2 (a) and h = τ F (b) based on the data in Table 1 and Table 2.
Fractalfract 07 00487 g003
Figure 4. Quantities Q 1 (a), Q 2 (b), Q 3 (c) and Q 4 (d) under mesh steps h = 1 / 8 , τ F = 1 / 64 .
Figure 4. Quantities Q 1 (a), Q 2 (b), Q 3 (c) and Q 4 (d) under mesh steps h = 1 / 8 , τ F = 1 / 64 .
Fractalfract 07 00487 g004
Figure 5. Numerical solutions of u ( x , t ) (a) and ρ ( x , t ) (b) at long time with h = 1 / 10 , τ F = 1 / 100 .
Figure 5. Numerical solutions of u ( x , t ) (a) and ρ ( x , t ) (b) at long time with h = 1 / 10 , τ F = 1 / 100 .
Fractalfract 07 00487 g005
Figure 6. Comparison of computation times of the SNFD scheme in [22] with present scheme at long time under h = 1 / 10 , τ F = 1 / 100 .
Figure 6. Comparison of computation times of the SNFD scheme in [22] with present scheme at long time under h = 1 / 10 , τ F = 1 / 100 .
Fractalfract 07 00487 g006
Table 1. The errors and convergence rates with τ F = h 2 .
Table 1. The errors and convergence rates with τ F = h 2 .
SNFD Scheme in [22]
h , τ F e SNFD ( h , τ F ) uRate S N F D x η SNFD ( h , τ F ) ρ Rate S N F D x CPU(s)
1 2 , 1 4 8.2798 × 10 2 2.9797 × 10 1 0.14
1 4 , 1 16 2.0361 × 10 2 2.02 7.3503 × 10 2 2.020.51
1 8 , 1 64 5.0629 × 10 3 2.01 1.8306 × 10 2 2.0112.30
1 16 , 1 256 1.2641 × 10 3 2.00 4.5719 × 10 3 2.00136.31
1 32 , 1 1024 3.1597 × 10 4 2.00 1.1427 × 10 3 2.001943.47
Present scheme
h , τ F e TT M ( h , τ F ) uRate TT M x η TT M ( h , τ F ) ρ Rate TT M x CPU(s)
1 2 , 1 4 8.3934 × 10 2 3.0015 × 10 1 0.10
1 4 , 1 16 2.0360 × 10 2 2.04 7.3555 × 10 2 2.030.30
1 8 , 1 64 5.0616 × 10 3 2.01 1.8308 × 10 2 2.017.53
1 16 , 1 256 1.2639 × 10 3 2.00 4.5721 × 10 3 2.0076.99
1 32 , 1 1024 3.1596 × 10 4 2.00 1.1427 × 10 3 2.001297.88
Table 2. The errors and convergence rates with h = τ F .
Table 2. The errors and convergence rates with h = τ F .
SNFD Scheme in [22]
h , τ F e SNFD ( h , τ F ) uRate S N F D t η SNFD ( h , τ F ) ρ Rate S N F D t CPU(s)
1 8 , 1 8 3.1920 × 10 2 1.1381 × 10 1 2.47
1 16 , 1 16 1.5552 × 10 2 1.04 5.5256 × 10 2 1.0414.09
1 32 , 1 32 7.6748 × 10 3 1.02 2.7221 × 10 2 1.02109.20
1 64 , 1 64 3.8121 × 10 3 1.01 1.3509 × 10 2 1.01770.01
1 128 , 1 128 1.8997 × 10 3 1.00 6.7298 × 10 3 1.005792.53
Present scheme
h , τ F e TT M ( h , τ F ) uRate TT M t η TT M ( h , τ F ) ρ Rate TT M t CPU(s)
1 8 , 1 8 3.1962 × 10 2 1.1414 × 10 1 1.49
1 16 , 1 16 1.5540 × 10 2 1.04 5.5323 × 10 2 1.047.10
1 32 , 1 32 7.6694 × 10 3 1.02 2.7237 × 10 2 1.0251.40
1 64 , 1 64 3.8102 × 10 3 1.01 1.3513 × 10 2 1.01414.33
1 128 , 1 128 1.8991 × 10 3 1.00 6.7308 × 10 3 1.013289.71
Table 3. Quantity Q 1 under different mesh steps h and τ F at various times.
Table 3. Quantity Q 1 under different mesh steps h and τ F at various times.
TT-M Finite Difference Scheme
( 1 2 , 1 4 ) ( 1 4 , 1 16 ) ( 1 8 , 1 64 ) ( 1 16 , 1 256 )
t = 0.0 4.47213595494.47213595494.47213595494.4721359549
t = 0.5 4.47213595494.47213595494.47213595494.4721359549
t = 1.0 4.47213595494.47213595494.47213595494.4721359549
t = 1.5 4.47213595494.47213595494.47213595494.4721359549
t = 2.0 4.47213595494.47213595494.47213595494.4721359549
t = 2.5 4.47213595494.47213595494.47213595494.4721359549
t = 3.0 4.47213595494.47213595494.47213595494.4721359549
t = 3.5 4.47213595494.47213595494.47213595494.4721359549
t = 4.0 4.47213595494.47213595494.47213595494.4721359549
Table 4. Quantity Q 2 under different mesh steps h and τ F at various times.
Table 4. Quantity Q 2 under different mesh steps h and τ F at various times.
TT-M Finite Difference Scheme
( 1 2 , 1 4 ) ( 1 4 , 1 16 ) ( 1 8 , 1 64 ) ( 1 16 , 1 256 )
t = 0.0 6.70820393246.70820393246.70820393246.7082039324
t = 0.5 6.70820393246.70820393246.70820393246.7082039324
t = 1.0 6.70820393246.70820393246.70820393246.7082039324
t = 1.5 6.70820393246.70820393246.70820393246.7082039324
t = 2.0 6.70820393246.70820393246.70820393246.7082039324
t = 2.5 6.70820393246.70820393246.70820393246.7082039324
t = 3.0 6.70820393246.70820393246.70820393246.7082039324
t = 3.5 6.70820393246.70820393246.70820393246.7082039324
t = 4.0 6.70820393236.70820393236.70820393236.7082039323
Table 5. Quantity Q 3 under different mesh steps h and τ F at various times.
Table 5. Quantity Q 3 under different mesh steps h and τ F at various times.
TT-M Finite Difference Scheme
( 1 2 , 1 4 ) ( 1 4 , 1 16 ) ( 1 8 , 1 64 ) ( 1 16 , 1 256 )
t = 0.0 17.381436010017.389076477817.390998209517.3914793723
t = 0.5 17.387936445917.389492542717.391024752417.3914810404
t = 1.0 17.402157526017.390443445217.391085054417.3914848239
t = 1.5 17.417894819517.391408984317.391145366017.3914885933
t = 2.0 17.427957514917.391992404517.391180361817.3914907570
t = 2.5 17.432062473617.392108994017.391185444517.3914910387
t = 3.0 17.429885221317.391907241917.391171002417.3914901056
t = 3.5 17.425724228917.391591623117.391150403617.3914888028
t = 4.0 17.420800446017.391304870417.391132470117.3914876800
Table 6. Quantity Q 4 under different mesh steps h and τ F at various times.
Table 6. Quantity Q 4 under different mesh steps h and τ F at various times.
TT-M Finite Difference Scheme
( 1 2 , 1 4 ) ( 1 4 , 1 16 ) ( 1 8 , 1 64 ) ( 1 16 , 1 256 )
t = 0.0 29.864568509529.827080011129.817468377729.8150480627
t = 0.5 29.519499827929.737796480529.794967795429.8094115019
t = 1.0 29.295309435229.680430505629.780509942529.8057894879
t = 1.5 29.237315733929.666565115029.777099727829.8049405034
t = 2.0 29.305203410129.686437337129.782246897329.8062391542
t = 2.5 29.433300260329.720401894929.790926448429.8084218678
t = 3.0 29.552089751129.751920910329.798937602529.8104338587
t = 3.5 29.635532719629.773318952129.804361420029.8117950932
t = 4.0 29.676564410829.784346385929.807161305429.8124980105
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gao, J.; He, S.; Bai, Q.; Liu, J. A Time Two-Mesh Finite Difference Numerical Scheme for the Symmetric Regularized Long Wave Equation. Fractal Fract. 2023, 7, 487. https://doi.org/10.3390/fractalfract7060487

AMA Style

Gao J, He S, Bai Q, Liu J. A Time Two-Mesh Finite Difference Numerical Scheme for the Symmetric Regularized Long Wave Equation. Fractal and Fractional. 2023; 7(6):487. https://doi.org/10.3390/fractalfract7060487

Chicago/Turabian Style

Gao, Jingying, Siriguleng He, Qingmei Bai, and Jie Liu. 2023. "A Time Two-Mesh Finite Difference Numerical Scheme for the Symmetric Regularized Long Wave Equation" Fractal and Fractional 7, no. 6: 487. https://doi.org/10.3390/fractalfract7060487

APA Style

Gao, J., He, S., Bai, Q., & Liu, J. (2023). A Time Two-Mesh Finite Difference Numerical Scheme for the Symmetric Regularized Long Wave Equation. Fractal and Fractional, 7(6), 487. https://doi.org/10.3390/fractalfract7060487

Article Metrics

Back to TopTop