AMP Advances in Mathematical Physics 1687-9139 1687-9120 Hindawi Publishing Corporation 423718 10.1155/2013/423718 423718 Research Article Conservative Linear Difference Scheme for Rosenau-KdV Equation Hu Jinsong 1 Xu Youcai 2 Hu Bing 2 Neidhardt Hagen 1 School of Mathematics and Computer Engineering Xihua University Chengdu 610039 China xhu.edu.cn 2 School of Mathematics Sichuan University Chengdu 610064 China scu.edu.cn 2013 27 3 2013 2013 05 02 2013 22 03 2013 2013 Copyright © 2013 Jinsong Hu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

A conservative three-level linear finite difference scheme for the numerical solution of the initial-boundary value problem of Rosenau-KdV equation is proposed. The difference scheme simulates two conservative quantities of the problem well. The existence and uniqueness of the difference solution are proved. It is shown that the finite difference scheme is of second-order convergence and unconditionally stable. Numerical experiments verify the theoretical results.

1. Introduction

KdV equation has been used in very wide applications and undergone research which can be used to describe wave propagation and spread interaction as follows : (1)ut+uux+uxxx=0. In the study of the dynamics of dense discrete systems, the case of wave-wave and wave-wall interactions cannot be described using the well-known KdV equation. To overcome this shortcoming of the KdV equation, Rosenau [5, 6] proposed the so-called Rosenau equation: (2)ut+uxxxxt+ux+uux=0. The existence and the uniqueness of the solution for (2) were proved by Park , but it is difficult to find the analytical solution for (2). Since then, much work has been done on the numerical method for (2) ( and also the references therein). On the other hand, for the further consideration of the nonlinear wave, the viscous term +uxxx needs to be included  (3)ut+uxxxxt+ux+uux+uxxx=0. This equation is usually called the Rosenau-KdV equation. Zuo  discussed the solitary wave solutions and periodic solutions for (2). Recently,  discussed the solitary solutions for the generalized Rosenau-KdV equation with usual power law nonlinearity. In [15, 16], the authors also gave the two invariants for the generalized Rosenau-KdV equation. In particular,  not only derived the singular 1-solition solution by the ansatz method but also used perturbation theory to obtain the adiabatic parameter dynamics of the water waves. In , The ansatz method is applied to obtain the topological soliton solution of the generalized Rosenau-KdV equation. The G/G method as well as the exp-function method are also applied to extract a few more solutions to this equation. But the numerical method to the initial-boundary value problem of Rosenau-KdV equation has not been studied till now. In this paper, we propose a conservative three-level finite difference scheme for the Rosenau-KdV equation (3) with the boundary conditions (4)u(xL,t)=u(xR,t)=0,ux(xL,t)=ux(xR,t)=0,uxx(xL,t)=uxx(xR,t)=0,t[0,T], and an initial condition (5)u(x,0)=u0(x),x[xL,xR]. The initial boundary value problem (3)–(5) possesses the following conservative properties : (6)Q(t)=xLxRu(x,t)dx=xLxRu0(x)dx=Q(0),(7)E(t)=uL22+uxxL22=E(0).

The solitary wave solution for (3) is [14, 15] (8)u(x,t)=(-3524+35312313)×sech4[124-26+2313×sech4×(x-(12+126313)t)124-26+2313]. When -xL0,xR0, the initial-boundary value problem (3)–(5) and the Cauchy problem (3) are consistent, so the boundary condition (4) is reasonable.

It is known the conservative scheme is better than the nonconservative ones. The nonconservative scheme may easily show nonlinear blow up. A lot of numerical experiments show that the conservative scheme can possesses some invariant properties of the original differential equation . The conservative scheme is more suitable for long-time calculations. In , Li and Vu-Quoc said “… in some areas, the ability to preserve some invariant properties of the original differential equation is a criterion to judge the success of a numerical simulation.” In this paper, we propose a three-level linear finite difference scheme for the Rosenau-KdV equation (3)–(5). The difference scheme is conservative which simulates conservative properties (6) and (7) at the same time.

The rest of this paper is organized as follows. In Section 2, we propose a three-level linear finite difference scheme for the Rosenau-KdV equation and discuss the discrete conservative properties. In Section 3, we show that the scheme is uniquely solvable. Then, in Section 4, we prove that the finite difference scheme is of second-order convergence, unconditionally stable. Finally, some numerical tests are given in Section 5 to verify our theoretical analysis.

2. Finite Difference Scheme and Conservation Properties

Let h=(xR-xL)/J and τ be the uniform step size in the spatial and temporal direction, respectively. Denote xj=xL+jh(j=-1,0,1,2,,J,J+1), tn=nτ(n=0,1,2,,N, N=[T/τ]), ujnu(xj,tn) and Zh0={u=(uj)u-1=u0=uJ=uJ+1=0,j=-1, 0,1,2,,J,J+1}. Throughout this paper, we denote C as a generic positive constant independent of h and τ, which may have different values in different occurrences. We introduce the following notations: (9)(ujn)x=uj+1n-ujnh,(ujn)x¯=ujn-uj-1nh,(ujn)x^=uj+1n-uj-1n2h,(ujn)t^=ujn+1-ujn-12τ,u¯jn=ujn+1+ujn-12,un,vn=hj=1J-1ujnvjn,un2=un,un,un=max1jJ-1ujn.

We propose a three-level linear finite difference scheme for the solution of (3)–(5) as follows: (10)(ujn)t^+(ujn)xxx¯x¯t^+(u¯jn)x^+(u¯jn)xx¯x^+13[ujn(u¯jn)x^+(ujnu¯jn)x^]=0,(11)j=1,2,3,,J-1,n=1,2,3,,N-1,(12)uj0=u0(xj),j=0,1,2,3,,J,(13)unZh0,(u0n)x^=(uJn)x^=0,(u0n)xx¯=(uJn)xx¯=0,n=1,2,3,,N. From the boundary conditions (4), we know that (13) is reasonable.

Lemma 1.

It follows from summation by parts that for any two mesh functions u,vZh0, (14)ux,v=-u,vx,uxx¯,v=-ux,vx. Then one has (15)ux,u=-u,ux¯,uxx¯,u=-ux,ux=-ux2. Furthermore, if (u0)xx¯=(uJ)xx¯=0, then (16)uxxx¯x¯,u=  uxx2.

The difference scheme (10)–(13) simulates two conservative properties of the problems (6) and (7) as follows.

Theorem 2.

Suppose that u0H02[xL,xR],u(x,t)C5,3, then the difference scheme (10)–(13) is conservative: (17)Qn=h2j=1J-1(ujn+1+ujn)+h6τj=1J-1ujn(ujn+1)x^=Qn-1==Q0,(18)En=12(un+12+un2)+12(uxxn+12+uxxn2)=En-1==E0.

Proof.

Multiplying (10) with h, summing up for j from 1 to J-1, and considering the boundary condition (13) and Lemma 1, we get (19)hj=1J-1ujn+1-ujn-12τ+h6j=1J-1[ujn(ujn+1)x^-ujn-1(ujn)x^]=0. Then, (17) is gotten from (19).

Taking an inner product of (10) with 2u¯n (i.e., un+1+un-1), considering the boundary condition (13) and Lemma 1, we obtain (20)12τ(un+12-un-12)+12τ(uxxn+12-uxxn-12)+2u¯x^n,u¯n+2u¯xx¯x^n,u¯n+2P,u¯n=0, where Pj=(1/3)[ujn(u¯jn)x^+(ujnu¯jn)x^]. According to (21)u¯x^n,u¯n=0,u¯xx¯x^n,u¯n=0,(22)P,u¯n=13hj=1J-1[ujn(u¯jn)x^+(ujnu¯jn)x^]u¯jn=112j=1J-1[ujn(uj+1n+1+uj+1n-1-uj-1n+1-uj-1n-1)=112j=1J-1+uj+1n(uj+1n+1+uj+1n-1)-uj-1n(uj-1n+1+uj-1n-1)]=112j=1J-1×(ujn+1+ujn-1)=112j=1J-1(ujn+uj+1n)(uj+1n+1+uj+1n-1)(ujn+1+ujn-1)-112j=1J-1(ujn+uj-1n)(ujn+1+ujn-1)(uj-1n+1+uj-1n-1)=0, we have (23)(un+12-un-12)+(uxxn+12-uxxn-12)=0. Then, (18) is gotten from (23).

3. Solvability Theorem 3.

There exists unZh0 which satisfies the difference scheme (10)–(13), (1nN).

Proof.

Use mathematical induction to prove it. It is obvious that u0 is uniquely determined by the initial condition (12). We also can get u1 in order O(h2+τ2) by two-level C-N scheme (i.e., u0 and u1 are uniquely determined). Now suppose u0,u1,,un  (1nN-1) is solved uniquely. Consider the equation of (10) for un+1:(24)12τujn+1+12τ(ujn+1)xxx¯x¯+12(ujn+1)x^+12(ujn+1)xx¯x^+16[ujn(ujn+1)x^+(ujnujn+1)x^]=0. Taking an inner product of (24) with un+1, we get (25)12τun+12+12τuxxn+12+12ux^n+1,un+1+12uxx¯x^n+1,un+1+h6j=1J-1[ujn(ujn+1)x^+(ujnujn+1)x^]ujn+1=0. Similar to the proof of (21), we have (26)ux^n+1,un+1=0,uxx¯x^n+1,un+1=0. By (27)h6j=1J-1[ujn(ujn+1)x^+(ujnujn+1)x^]ujn+1=112j=1J-1[ujn(uj+1n+1-uj-1n+1)+(uj+1nuj+1n+1-uj-1nuj-1n+1)]ujn+1=112j=1J-1[ujnujn+1uj+1n+1+uj+1nujn+1uj+1n+1]  -112j=1J-1[uj-1nuj-1n+1ujn+1+ujnuj-1n+1ujn+1]=0, and from (25)–(27), we have (28)un+12+uxxn+12=0. That is, (24) has only a trivial solution. Therefore, (10) determines ujn+1 uniquely. This completes the proof.

4. Convergence and Stability

Let v(x,t) be the solution of problem (3)–(5), vjn=v(xj,tn), then the truncation error of the difference scheme (10)–(13) is as follows: (29)rjn=(vjn)t^+(vjn)xxx¯x¯t^+(v¯jn)x^+(v¯jn)xx¯x^+13[vjn(v¯jn)x^+(vjnv¯jn)x^]. Making use of Taylor expansion, we know that rjn=O(τ2+h2) holds if h,τ0.

Lemma 4.

Suppose that u0H02[xL,xR], then the solution un of (3)–(5) satisfies (30)uL2C,uxL2C,uLC,uxLC.

Proof.

It is follows from (7) that (31)uL2C,uxxL2C. By Holder inequality and Schwarz inequality, we get (32)uxL22=xLxRuxuxdx=uux|xLxR-xLxRuuxxdx=-xLxRuuxxdxuL2·uxxL212(uL22+uxxL22), which implies that (33)uxL2C. Using Sobolev inequality, we get that uLC,  uxLC.

Lemma 5 (discrete Sobolev's inequality [<xref ref-type="bibr" rid="B27">27</xref>]).

There exist two constants C1 and C2 such that (34)unC1un+C2uxn.

Lemma 6 (discrete Gronwall inequality [<xref ref-type="bibr" rid="B27">27</xref>]).

Suppose that w(k) and  ρ(k) are nonnegative function and ρ(k) is nondecreasing. If C>0, and (35)w(k)ρ(k)+Cτl=0k-1w(l),k, then (36)w(k)ρ(k)eCτk,k.

Theorem 7.

Suppose u0H02[xL,xR], then the solution of (10)–(13) satisfies: unC,  uxnC,  uxxnC, which yield unC,  uxnC  (n=1,2,,N).

Proof.

It is follows from (18) that (37)unC,uxxnC. According to (15) and Schwarz inequality, we get (38)uxn2un·uxxn12(un2+uxxn2)C. Using Lemma 5, we have unC,  uxnC.

Theorem 8.

Suppose u0H02[xL,xR], u(x,t)C5,3, then the solution un of the difference scheme (10)–(13) converges to the solution v(x,t) of the problem (3)–(5) with order O(τ2+h2) in norm ·.

Proof.

Subtracting (10) from (29) and letting ejn=vjn-ujn, we have (39)rjn=(ejn)t^+(e¯jn)xxx¯x¯t^+(e¯jn)x^+(e¯jn)xx¯x^+R1,j+R2,j, where R1,j=(1/3)[vjn(v¯jn)x^-ujn(u¯jn)x^],  R2,j=(1/3)[(vjnv¯jn)x^-(ujnu¯jn)x^]. Computing the inner product of (39) with 2e¯n, we obtain (40)rn,2e¯n=12τ(en+12-en-12)+12τ(exxn+12-exxn-12)+e¯x^n,2e¯n+e¯xx¯x^n,2e¯n+R1,2e¯n+R2,2e¯n. Similar to the proof of (21), we have (41)e¯x^n,2e¯n=0,e¯xx¯x^n,2e¯n=0. Then, (40) can be rewritten as follows: (42)(en+12-en+12)+(exxn+12-exxn-12)=2τrn,2e¯n+2τ-R1,2e¯n+2τ-R2,2e¯n. Using Lemma 4 and Theorem 7, we get (43)|vjn|C,|ujn|C,|(ujn)x^|  C,(j=0,1,2,,J;  n=1,2,,N). According to the Schwarz inequality, we obtain (44)-R1,2e¯n=-23hj=1J-1[vjn(v¯jn)x^-ujn(u¯jn)x^]e¯jn=-23hj=1J-1[vjn(e¯jn)x^+ejn(u¯jn)x^]e¯jn23Chj=1J-1(|(e¯jn)x^|+|ejn|)|e¯jn|C[e¯xn2+en2+e¯n2]C[en+12+en-12+en2+exn+12+exn-12],-R2,2e¯n=-23hj=1J-1[(vjnv¯jn)x^-(ujnu¯jn)x^]e¯jn=23hj=1J-1[vjnv¯jn-ujnu¯jn](e¯jn)x^=23hj=1J-1[vjne¯jn+ejnu¯jn](e¯jn)x^23Chj=1J-1(|(e¯jn)|+|ejn|)|(e¯jn)x^|C[e¯xn2+en2+e¯n2]C[en+12+en-12+en2+exn+12+exn-12]. Noting that (45)rn,2e¯n=rn,en+1+en-1rn2+12[en+12+en-12], and from (42)–(45), we have (46)(en+12-en-12)+(exxn+12-exxn-12)Cτ[en+12+en2+en-12+exn+12+exn2+exn-12]+2τrn2. Similar to the proof of (38), we have (47)exn+1212(en+12+exxn+12),exn212(en2+exxn2),exn-1212(en-12+exxn-12). Then, (46) can be rewritten as (48)(en+12-en-12)+(exxn+12-exxn-12)Cτ[en+12+en2+en-12Cτ+exxn+12+exxn2+exxn-12]+2τrn2. Let Bn=en+12+exxn+12+en2+exxn2. Then, (48) can be rewritten as follows: (49)Bn-Bn-1Cτ(Bn+Bn-1)+2τrn2, which yields (50)(1-Cτ)(Bn-Bn-1)2CτBn-1+2τrn2. If τ is sufficiently small, which satisfies 1-Cτ>0, then (51)Bn-Bn-1CτBn-1+Cτrn2. Summing up (51) from 1 to n, we have (52)BnB0+Cτl=1nrl2+Cτl=0n-1Bl. First, we can get u1 in order O(τ2+h2) that satisfies B0=O(τ2+h2)2 by two-level C-N scheme. Since (53)τl=1nrl2nτmax1lnrl2T·O(τ2+h2)2, then we obtain (54)BnO(τ2+h2)2+Cτl=0n-1Bl. From Lemma 6 we get (55)BnO(τ2+h2)2, which implies that (56)enO(τ2+h2),exxnO(τ2+h2). From (47) we have (57)exnO(τ2+h2). By Lemma 5 we obtain (58)enO(τ2+h2).

Finally, we can similarly prove results as follows.

Theorem 9.

Under the conditions of Theorem 8, the solution un of (10)–(13) is stable in norm ·.

5. Numerical Simulations

Since the three-level implicit finite difference scheme cannot start by itself, we need to select other two-level schemes (such as the C-N Scheme) to get u1. Then, be reusing initial value u0, we can work out u2,u3,. Iterative numerical calculation is not required, for this scheme is linear, so it saves computing time. Let xL=-70,xR=100,andT=40, (59)u0(x)=(-3524+35312313)sech4(124-26+2313x).

In Table 1, we give the error at various time steps. Using the method in [30, 31], we verified the second convergence of the difference scheme in Table 2. Numerical simulations on two conservation invariants Qn and En are given in Table 3.

The error at various time steps.

τ = h = 0.1 τ = h = 0.05 τ = h = 0.025
e n e n e n e n e n e n
t = 10 1.641934 e - 3 6.314193 e - 4 4.113510 e - 4 1.582641 e - 4 1.028173 e - 4 3.965867 e - 5
t = 20 3.045414 e - 3 1.131442 e - 3 7.631169 e - 4 2.835874 e - 4 1.905450 e - 4 7.097948 e - 5
t = 30 4.241827 e - 3 1.533771 e - 3 1.062971 e - 3 3.843906 e - 4 2.650990 e - 4 9.610332 e - 5
t = 40 5.297873 e - 3 1.878952 e - 3 1.327645 e - 3 4.709118 e - 4 3.306738 e - 4 1.176011 e - 4

The verification of the second convergence.

e n ( h , τ ) / e 2 n ( h / 2 , τ / 2 ) e n ( h , τ ) / e 2 n ( h / 2 , τ / 2 )
τ = h = 0.1 τ = h = 0.05 τ = h = 0.025 τ = h = 0.1 τ = h = 0.05 τ = h = 0.025
t = 10 3.991564 4.000797 3.989657 3.990655
t = 20 3.990757 4.004916 3.989749 3.995343
t = 30 3.990539 4.009713 3.990136 3.999764
t = 40 3.990427 4.014970 3.990030 4.004314

Numerical simulations on the two conservation invariants Qn and En.

τ = h = 0.1 τ = h = 0.05 τ = h = 0.025
Q n E n Q n E n Q n E n
t = 0 5.497722548019 1.984553365290 5.498060684522 1.984390175264 5.498145418391 1.984349335263
t = 10 5.497724936513 1.984595075859 5.498060837192 1.984401029470 5.498145479109 1.984352109750
t = 20 5.497728744900 1.984645964099 5.498061080542 1.984414367496 5.498145545374 1.984355520610
t = 30 5.497731963790 1.984679827211 5.498061287046 1.984423270337 5.498145609535 1.984357811266
t = 40 5.497734235191 1.984701501262 5.498061398506 1.984428974030 5.498145659050 1.984359292230

The wave graph comparison of u(x,t) between τ=h=0.1 and τ=h=0.025 at various times is given in Figures 1 and 2.

When τ=h=0.1, the wave graph of u(x,t) at various times.

When τ=h=0.025, the wave graph of u(x,t) at various times.

Numerical simulations show that the finite difference scheme is efficient.

Acknowledgments

The work was supported by Scientific Research Fund of Sichuan Provincial Education Department (11ZB009) and the fund of Key Disciplinary of Computer Software and Theory, Sichuan, Grant no. SZD0802-09-1.

Cui Y. Mao D.-k. Numerical method satisfying the first two conservation laws for the Korteweg-de Vries equation Journal of Computational Physics 2007 227 1 376 399 10.1016/j.jcp.2007.07.031 MR2361527 ZBL1131.65073 Zhu S. Zhao J. The alternating segment explicit-implicit scheme for the dispersive equation Applied Mathematics Letters 2001 14 6 657 662 10.1016/S0893-9659(01)80022-7 MR1836066 ZBL0996.65083 Bahadır A. R. Exponential finite-difference method applied to Korteweg-de Vries equation for small times Applied Mathematics and Computation 2005 160 3 675 682 10.1016/j.amc.2003.11.025 MR2111671 ZBL1062.65087 Özer S. Kutluay S. An analytical-numerical method for solving the Korteweg-de Vries equation Applied Mathematics and Computation 2005 164 3 789 797 10.1016/j.amc.2004.06.011 MR2135019 ZBL1070.65077 Rosenau P. A quasi-continuous description of a nonlinear transmission line Physica Scripta 1986 34 827 829 10.1088/0031-8949/34/6B/020 Rosenau P. Dynamics of dense discrete systems Progress of Theoretical Physics 1988 79 1028 1042 10.1143/PTP.79.1028 Park M. A. On the Rosenau equation Matemática Aplicada e Computacional 1990 9 2 145 152 MR1084695 ZBL0723.35071 Chung S. K. Ha S. N. Finite element Galerkin solutions for the Rosenau equation Applicable Analysis 1994 54 1-2 39 56 10.1080/00036819408840267 MR1382206 ZBL0830.65097 Omrani K. Abidi F. Achouri T. Khiari N. A new conservative finite difference scheme for the Rosenau equation Applied Mathematics and Computation 2008 201 1-2 35 43 10.1016/j.amc.2007.11.039 MR2432579 ZBL1156.65078 Chung S. K. Finite difference approximate solutions for the Rosenau equation Applicable Analysis 1998 69 1-2 149 156 10.1080/00036819808840652 MR1708193 ZBL0904.65093 Chung S. K. Pani A. K. Numerical methods for the Rosenau equation Applicable Analysis 2001 77 3-4 351 369 10.1080/00036810108840914 MR1975741 ZBL1021.65048 Manickam S. A. V. Pani A. K. Chung S. K. A second-order splitting combined with orthogonal cubic spline collocation method for the Rosenau equation Numerical Methods for Partial Differential Equations 1998 14 6 695 716 MR1653330 ZBL0930.65111 Kim Y. D. Lee H. Y. The convergence of finite element Galerkin solution for the Roseneau equation The Korean Journal of Computational & Applied Mathematics 1998 5 1 171 180 MR1600314 ZBL0977.65080 Zuo J.-M. Solitons and periodic solutions for the Rosenau-KdV and Rosenau-Kawahara equations Applied Mathematics and Computation 2009 215 2 835 840 10.1016/j.amc.2009.06.011 MR2561541 ZBL1175.65124 Esfahani A. Solitary wave solutions for generalized Rosenau-KdV equation Communications in Theoretical Physics 2011 55 3 396 398 10.1088/0253-6102/55/3/04 MR2884092 Razborova P. Triki H. Biswas A. Perturbation of dispersive shallow water waves Ocean Engineering 2013 63 1 7 10.1016/j.oceaneng.2013.01.014 Ebadi G. Mojaver A. Triki H. Yildirim A. Biswas A. Topological solitons and other solutions of the Rosenau-KdV equation with power law nonlinearity Romanian Journal of Physics 2013 58 1-2 1 10 Wang T. Zhang L. Chen F. Conservative schemes for the symmetric regularized long wave equations Applied Mathematics and Computation 2007 190 2 1063 1080 10.1016/j.amc.2007.01.105 MR2339701 ZBL1124.65080 Li S. Vu-Quoc L. Finite difference calculus invariant structure of a class of algorithms for the nonlinear Klein-Gordon equation SIAM Journal on Numerical Analysis 1995 32 6 1839 1875 10.1137/0732083 MR1360462 ZBL0847.65062 Chang Q. Jia E. Sun W. Difference schemes for solving the generalized nonlinear Schrödinger equation Journal of Computational Physics 1999 148 2 397 415 10.1006/jcph.1998.6120 MR1669707 ZBL0923.65059 Wang T.-C. Zhang L.-M. Analysis of some new conservative schemes for nonlinear Schrödinger equation with wave operator Applied Mathematics and Computation 2006 182 2 1780 1794 10.1016/j.amc.2006.06.015 MR2282620 ZBL1161.65349 Wang T. Guo B. Zhang L. New conservative difference schemes for a coupled nonlinear Schrödinger system Applied Mathematics and Computation 2010 217 4 1604 1619 10.1016/j.amc.2009.07.040 MR2727403 ZBL1205.65242 Zhang L. A finite difference scheme for generalized regularized long-wave equation Applied Mathematics and Computation 2005 168 2 962 972 10.1016/j.amc.2004.09.027 MR2171754 ZBL1080.65079 Fei Z. Vázquez L. Two energy conserving numerical schemes for the sine-Gordon equation Applied Mathematics and Computation 1991 45 1 17 30 10.1016/0096-3003(91)90087-4 MR1122000 ZBL0732.65107 Wong Y. S. Chang Q. Gong L. An initial-boundary value problem of a nonlinear Klein-Gordon equation Applied Mathematics and Computation 1997 84 1 77 93 10.1016/S0096-3003(96)00065-3 MR1444445 ZBL0884.65091 Chang Q. S. Guo B. L. Jiang H. Finite difference method for generalized Zakharov equations Mathematics of Computation 1995 64 210 537 553 10.2307/2153438 MR1284664 ZBL0827.65138 Hu J. Zheng K. Two conservative difference schemes for the generalized Rosenau equation Boundary Value Problems 2010 18 543503 MR2600713 ZBL1187.65090 Pan X. Zhang L. On the convergence of a conservative numerical scheme for the usual Rosenau-RLW equation Applied Mathematical Modelling 2012 36 8 3371 3378 10.1016/j.apm.2011.08.022 MR2920897 ZBL1252.65144 Pan X. Zhang L. Numerical simulation for general Rosenau-RLW equation: an average linearized conservative scheme Mathematical Problems in Engineering 2012 2012 15 517818 MR2922125 Wang T. Guo B. A robust semi-explicit difference scheme for the Kuramoto-Tsuzuki equation Journal of Computational and Applied Mathematics 2009 233 4 878 888 10.1016/j.cam.2009.07.058 MR2557281 ZBL1181.65117 Hu B. Hu J. Xu Y. C-N difference schemes for dissipative symmetric regularized long wave equations with damping term Mathematical Problems in Engineering 2011 2011 16 2-s2.0-79959272870 10.1155/2011/651642 651642 ZBL1213.76125