We study the initial-boundary problem of dissipative symmetric regularized long wave equations with damping term by finite difference method. A linear three-level implicit finite difference scheme is designed. Existence and uniqueness of numerical solutions are derived. It is proved that the finite difference scheme is of second-order convergence and unconditionally stable by the discrete energy method. Numerical simulations verify that the method is accurate and efficient.
A symmetric version of regularized long wave equation (SRLWE),
has been proposed to model the propagation of weakly nonlinear ion acoustic and space charge waves . The solitary wave solutions are
The four invariants and some numerical results have been obtained in , where is the velocity, . Obviously, eliminating from (1.1), we get a class of SRLWE:
Equation (1.3) is explicitly symmetric in the and derivatives and is very similar to the regularized long wave equation that describes shallow water waves and plasma drift waves [2, 3]. The SRLW equation also arises in many other areas of mathematical physics [4–6]. Numerical investigation indicates that interactions of solitary waves are inelastic ; thus, the solitary wave of the SRLWE is not a solution. Research on the wellposedness for its solution and numerical methods has aroused more and more interest. In , Guo studied the existence, uniqueness, and regularity of the numerical solutions for the periodic initial value problem of generalized SRLW by the spectral method. In , Zheng et al. presented a Fourier pseudospectral method with a restraint operator for the SRLWEs and proved its stability and obtained the optimum error estimates. There are other methods such as pseudospectral method, finite difference method for the initial-boundary value problem of SRLWEs (see [9–15]).
In applications, the viscous damping effect is inevitable, and it plays the same important role as the dispersive effect. Therefore, it is more significant to study the dissipative symmetric regularized long wave equations with the damping term
where are positive constants, is the dissipative coefficient, and is the damping coefficient. Equations (1.4)-(1.5) are a reasonable model to render essential phenomena of nonlinear ion acoustic wave motion when dissipation is considered. Existence, uniqueness, and wellposedness of global solutions to (1.4)-(1.5) are presented (see [16–20]). But it is difficult to find the analytical solution to (1.4)-(1.5), which makes numerical solution important.
To authors' knowledge, the finite difference method to dissipative SRLWEs with damping term (1.4)-(1.5) has not been studied till now. In this paper, we propose linear three level implicit finite difference scheme for (1.4)-(1.5) with
and the boundary conditions
We show that this difference scheme is uniquely solvable, convergent, and stable in both theoretical and numerical senses.
Suppose that , , the solution of (1.4)–(1.7) satisfies , , , and , where is a generic positive constant that varies in the context.
Multiplying (1.4) by and integrating over , we have
Then, multiplying (1.5) by and integrating over , we have
Adding (1.14) to (1.11), we obtain
So is decreasing with respect to , which implies that , . Then, it indicates that , , and . It is followed from Sobolev inequality that .
2. Finite Difference Scheme and Its Error Estimation
Let and be the uniform step size in the spatial and temporal direction, respectively. Denote , , , , , and . We define the difference operators as follows:
Then, the average three-implicit finite difference scheme for the solution of (1.4)–(1.7) is as follow:
There exist two constants and such that
Suppose that , are nonnegative functions and is nondecreasing. If and
If , , then the solution of (2.2)–(2.5) satisfies
Taking an inner product of (2.2) with (i.e., ) and considering the boundary condition (2.5) and Lemma 2.1, we obtain
where . Since
Taking an inner product of (2.3) with (i.e., ), we obtain
Adding (2.12) to (2.13), we have
Equation (2.14) can be changed to
Let , and (2.16) is changed to
If is sufficiently small which satisfies , then
Summing up (2.18) from 1 to , we have
From Lemma 2.3, we obtain , which implies that, , , and . By Lemma 2.2, we obtain .
Assume that , , the solution of difference scheme (2.2)–(2.5) satisfies:
Differentiating backward (2.2)–(2.5) with respect to , we obtain
Computing the inner product of (2.21) with (i.e., ) and considering (2.24) and Lemma 2.1, we obtain
where . It follows from Theorem 2.4 that
By the Schwarz inequality and Lemma 2.1, we get
it follows from (2.25) that
Computing the inner product of (2.22) with (i.e., ) and considering (2.24) and Lemma 2.1, we obtain
then (2.30) is changed to
Adding (2.29) to (2.32), we have
Leting , we obtain . Choosing suitable which is small enough to satisfy , we get
Summing up (2.34) from 1 to , we have
By Lemma 2.3, we get , which implies that , . It follows from Theorem 2.4 and Lemma 2.2 that , .
The solution of (2.2)–(2.5) is unique.
Using the mathematical induction, clearly, , are uniquely determined by initial conditions (2.4). then select appropriate second-order methods (such as the C-N Schemes) and calculate and (i.e. , , and , are uniquely determined). Assume that and are the only solution, now consider and in (2.2) and (2.3):
Taking an inner product of (3.1) with , we have
then it holds
Taking an inner product of (3.2) with and adding to (3.5), we have
which implies that (3.1)-(3.2) have only zero solution. So the solution and of (2.2)–(2.5) is unique.
4. Convergence and Stability
Let and be the solution of problem (1.4)–(1.7); that is, , , then the truncation of the difference scheme (2.2)–(2.5) is
Making use of Taylor expansion, it holds if .
Assume that , , then the solution and in the senses of norms and , respectively, to the difference scheme (2.2)–(2.5) converges to the solution of problem (1.4)–(1.7) and the order of convergence is .
Subtracting (2.2) from (4.1) subtracting (2.3) from (4.2), and letting , , we have
Computing the inner product of (4.3) with , we get
it follow from Lemma 1.1, Theorems 2.4, and 2.5 that
By the Schwarz inequality, we obtain
it follows from (4.9)–(4.10) and (4.6) that
Computing the inner product of (4.4) with , we obtain
Adding (4.12) to (4.11), we have
If is sufficiently small which satisfies , then
Summing up (4.16) from 1 to , we have
Select appropriate second-order methods (such as the C-N Schemes), and calculate and , which satisfies
we then have
By Lemma 2.3, we get
By Lemma 2.2, we have
Similarly to Theorem 4.1, we can prove the result as follows.
Under the conditions of Theorem 4.1, the solution and of (2.2)–(2.5) is stable in the senses of norm and , respectively.
5. Numerical Simulations
Since the three-implicit finite difference scheme can not start by itself, we need to select other two-level schemes (such as the C-N Scheme) to get , . Then, reusing initial value , , we can work out . Iterative numerical calculation is not required, for this scheme is linear, so it saves computing time.
When , the damping does not have an effect and the dissipative will not appear. So the initial conditions of (1.4)–(1.7) are same as those of (1.1):
Let , , , and . Since we do not know the exact solution of (1.4)-(1.5), an error estimates method in  is used: a comparison between the numerical solutions on a coarse mesh and those on a refine mesh is made. We consider the solution on mesh as the reference solution. In Table 1, we give the ratios in the sense of at various time steps.
Table 1. The error ratios in the sense of at various time steps.
Figure 1. When , the wave graph of at various times.
Figure 2. When , the wave graph of at various times.
From Table 1, it is easy to find that the difference scheme in this paper is second-order convergent. Figures 1 and 2 show that the height of wave crest is more and more low with time elapsing due to the effect of damping and dissipativeness. It simulates that the continue energy of problem (1.4)–(1.7) in Lemma 1.1 is digressive. Numerical experiments show that the finite difference scheme is efficient.
The work of Jinsong Hu was supported by the research fund of key disciplinary of application mathematics of Xihua University (Grant no. XZD0910-09-1). The work of Youcai Xu was supported by the Youth Research Foundation of Sichuan University (no. 2009SCU11113).
Seyler, CE, Fenstermacher, DL: A symmetric regularized-long-wave equation. Physics of Fluids. 27(1), 4–7 (1984). Publisher Full Text
Albert, J: On the decay of solutions of the generalized Benjamin-Bona-Mahony equations. Journal of Mathematical Analysis and Applications. 141(2), 527–537 (1989). Publisher Full Text
Amick, CJ, Bona, JL, Schonbek, ME: Decay of solutions of some nonlinear wave equations. Journal of Differential Equations. 81(1), 1–49 (1989). Publisher Full Text
Ogino, T, Takeda, S: Computer simulation and analysis for the spherical and cylindrical ion-acoustic solitons. Journal of the Physical Society of Japan. 41(1), 257–264 (1976). Publisher Full Text
Clarkson, PA: New similarity reductions and Painlevé analysis for the symmetric regularised long wave and modified Benjamin-Bona-Mahoney equations. Journal of Physics A. 22(18), 3821–3848 (1989). Publisher Full Text
Bogolubsky, IL: Some examples of inelastic soliton interaction. Computer Physics Communications. 13(3), 149–155 (1977). Publisher Full Text
Wang, T, Zhang, L, Chen, F: Conservative schemes for the symmetric regularized long wave equations. Applied Mathematics and Computation. 190(2), 1063–1080 (2007). Publisher Full Text
Wang, TC, Zhang, LM, Chen, FQ: Pseudo-compact conservative finite difference approximate solutions for symmetric regularized-long-wave equations. Chinese Journal of Engineering Mathematics. 25(1), 169–172 (2008)
Shang, YD, Guo, B: Global attractors for a periodic initial value problem for dissipative generalized symmetric regularized long wave equations. Acta Mathematica Scientia. Series A. 23(6), 745–757 (2003)
Guo, B, Shang, Y: Approximate inertial manifolds to the generalized symmetric regularized long wave equations with damping term. Acta Mathematicae Applicatae Sinica. 19(2), 191–204 (2003). Publisher Full Text
Shaomei, F, Boling, G, Hua, Q: The existence of global attractors for a system of multi-dimensional symmetric regularized wave equations. Communications in Nonlinear Science and Numerical Simulation. 14(1), 61–68 (2009). Publisher Full Text
Hu, B, Xu, Y, Hu, J: Crank-Nicolson finite difference scheme for the Rosenau-Burgers equation. Applied Mathematics and Computation. 204(1), 311–316 (2008). Publisher Full Text