- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Recently Accepted Articles ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
Abstract and Applied Analysis
Volume 2013 (2013), Article ID 768976, 12 pages
Explicit Multistep Mixed Finite Element Method for RLW Equation
1School of Mathematical Sciences, Inner Mongolia University, Hohhot 010021, China
2School of Statistics and Mathematics, Inner Mongolia University of Finance and Economics, Hohhot 010070, China
Received 15 February 2013; Revised 22 April 2013; Accepted 30 April 2013
Academic Editor: Marco Donatelli
Copyright © 2013 Yang Liu 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.
An explicit multistep mixed finite element method is proposed and discussed for regularized long wave (RLW) equation. The spatial direction is approximated by the mixed Galerkin method using mixed linear space finite elements, and the time direction is discretized by the explicit multistep method. The optimal error estimates in and norms for the scalar unknown and its flux based on time explicit multistep method are derived. Some numerical results are given to verify our theoretical analysis and illustrate the efficiency of our method.
In this paper, we consider the following initial boundary problem of RLW equation: where is a bounded open interval, with . The initial value is given function, and the coefficients , , are all positive constants.
In recent years, large nonlinear phenomena are found in many research fields, for example, physics, biology, fluid dynamics, and so forth. These phenomena can be described by the mathematical model of some nonlinear evolution equations. In particular, some attention has also been paid to nonlinear RLW equations [1, 2] which play a very important role in the study of nonlinear dispersive waves. Solitary waves are wave packet or pulses, which propagate in nonlinear dispersive media. Due to dynamical balance between the nonlinear and dispersive effects, these waves retain an unchanged waveform. A soliton is a very special type of solitary wave, which also keeps its waveform after collision with other solitons. The regularized long wave (RLW) equation is an alternative description of nonlinear dispersive waves to the more usual Korteweg-de Vries (KdV) equation. Mathematical theories and numerical methods for (1) were considered in [1–24]. The existence and uniqueness of the solution of RLW equation are discussed in . Their analytical solutions were found under restricted initial and boundary conditions, and therefore they got interest from a numerical point of view. Several numerical methods for the solution of the RLW equation have been introduced in the literature. These include a variety of difference methods [5–8, 17], finite element methods based on Galerkin and collocation principles [9–13], mixed finite element methods [14–16], meshfree method , adomian decomposition method , and so on.
In , Chatzipantelidis studied the explicit multistep methods for some nonlinear partial differential equations and discussed some mathematical theories. Akrivis et al.  studied the multistep method for some nonlinear evolution equations. Mei and Chen  presented the explicit multistep method based on Galerkin method for regularized long wave (RLW) equation. In this paper, our purpose is to propose and study an explicit multistep mixed method, which combines a mixed Galerkin method in the spatial direction and the explicit multistep method in the time direction, for RLW equation. We derive optimal error estimates in and norms for the scalar unknown and its flux for the fully discrete explicit multistep mixed scheme and compare our method’s accuracy with some other numerical schemes. Compared to the numerical methods in [20, 25, 26], we not only obtain the approximation solution for , but also get the approximation solution for .
The layout of the paper is as follows. In Section 2, an explicit multistep mixed scheme and numerical process are given. The optimal error estimates in and norms for the scalar unknown and its flux for the fully discrete explicit multistep mixed scheme are proved in Section 3. In Section 4, some numerical results are shown to confirm our theoretical analysis. Finally, some concluding remarks are given in Section 5. Throughout this paper, will denote a generic positive constant which does not depend on the spatial mesh parameter or time discretization parameter .
2. The Mixed Numerical Scheme
With the auxiliary variable , we reformulate (1) as the following first-order coupled system: We consider the following mixed weak formulation of (2). Find satisfying: Noting the Dirichlet boundary conditions and , we can get easily and then get the scheme (4).
Let and be finite dimensional subspaces of and , respectively, defined by where is a partition of into subintervals , , , , and denotes the polynomials of degree less than or equal to in .
In the following discussion, we will give an explicit multistep mixed scheme. We take linear finite element spaces and , and then and can be expressed as the following formulation: Substitute (8) into (6) and (7), and take and in (6) and (7), respectively, to obtain where .
We subdivide the space variable domain into uniform subintervals with grid points , such that , . Using the local coordinate transformation , , we transform a subinterval into a standard interval . Furthemore, we have We take linear basis functions defined as follows: and then the variables and over the element are written as Then, we get the following equations: Then, the system (13) has the following matrix form: with the following element matrices: Assembling contributions from all elements, we obtain the following coupled system of nonlinear matrix equations: To formulate a fully discrete scheme, we consider a uniform partition of with time step length , , and time levels , . We now discuss a fully discrete scheme based on a linear explicit multistep method. We now define and as approximations to and , respectively, and formulate the following fully discrete linear explicit multistep mixed scheme: with given the initial approximations and . In the explicit multistep mixed system (17), the parameter variable and is described by the coefficients of the term , for the following polynomials and , respectively:
Remark 1. There have been many numerical schemes for the RLW equation, but we have not seen the related research on explicit multistep mixed element method for RLW equation in the literature. From the viewpoint of numerical theory, we propose a mixed element scheme (6) and (7), which is different from some other mixed finite element methods in [14–16], for the RLW equation and derive some a priori error estimates based on the explicit multistep mixed element method. From the perspective of numerical calculation, our method is efficient for RLW equation.
3. Two-Step Mixed Scheme and Optimal Error Estimates
3.1. Two-Step Mixed Scheme and Some Lemmas
In this section, we will discuss some a priori error estimates based on explicit 2-step mixed finite element method for the RLW equation. For the fully discrete procedure, let be a given partition of the time interval with step length , for some positive integer . For a smooth function on , define .
For the theoretical analysis of a priori error estimates, we define the following projections.
Lemma 3 (see [15, 27, 28]). Furthermore, one also defines a Ritz projection of as the solution of
where , and is taken appropriately so that
where is a positive constant. Moreover, it is easy to verify that is bounded.
With , the following estimates hold:
Lemma 4. For , , , and , the following estimates hold:
Proof. Using the Taylor expansion, we have
Combining (37) and (38) and noting that , we obtain
From (39), we have
Using a similar estimate as the one for , we have
where , .
From (41), we have Using the Taylor expansion and noting that , we have Using (43), we obtain By (44), we have Using the similar method to the estimate for and (31), we obtain
3.2. Optimal Error Estimates
In this subsection, we derive the fully discrete optimal error estimates and obtain the following theorem.
Theorem 5. Assuming that , and , are given, then for , , one has
Proof. Taking in (33) and using Cauchy-Schwarz’s inequality and Poincaré’s inequality, we get Set in (34) to obtain Use (48) as well as the Cauchy-Schwarz and Young’s inequalities to obtain Note that We now estimate , , , and as Substituting (53) into (52), we obtain Substituting (36), (51), and (54) into (50), using (48), and summing from , the resulting inequality becomes Choose in such a way that for , . Then, as an application of Gronwall’s lemma, we obtain Combine (28), (31), and (56) with the triangle inequality to complete the and error estimates for . Furthermore, use (48) and the triangle inequality to complete the optimal error estimates for and .
Remark 7. Although some convergence proofs of multistep methods for RLW/BBM are provided in [20, 25, 29], our convergence results of multistep methods are proved based on a mixed finite element scheme. Based on the current discussion, we have to provide the detailed proofs for multistep mixed finite element methods in this paper.
4. Numerical Results
In order to test the viability of the proposed method, we consider a test problem. We write conservation laws as [3, 4] where , , and are usually called mass, momentum, and energy, respectively, which are observed to check the conservation of the numerical scheme.
We consider the motion of a single solitary wave and take as initial condition, with and , The corresponding exact solution with initial condition (60) is
In this procedure, we take space-time domain as and .
In Table 1, we take spatial mesh parameter and time discretization parameter and list the three invariants , , and and the optimal error estimate in and norms for at different times , , , , , and . At the same time, we show some numerical results at time obtained by other numerical methods in Table 1. From Table 1, we find that our method is more accurate than the numerical methods in [17, 18, 28] but is less than the numerical methods in [23, 24]. From the shown results in Table 1, we can see that , and keep almost constants, so the conservation for our method is very well.
In Tables 2 and 3, we take spatial mesh parameter and obtain the optimal error results in and norms for at different times , , , , , and with different time discretization parameters , and . From Tables 2 and 3, we see easily that the convergence rate for time is close to order 2. Similarly, the results for are shown in Tables 4 and 5.
In Tables 6 and 7, the optimal error results in and norms for at different times , , , , , and with different spatial mesh parameters , , and and time discretization parameter are shown. It is easy to see that the convergence rate for space is close to order 2. The similar results for are listed in Tables 8 and 9.
Figure 1 shows the surface for the exact solution in space-time domain (), and the corresponding surface for the numerical solution with and is described in Figure 2. From Figures 1 and 2, we see easily that the exact solution is approximated very well by the numerical solution . In Figures 3 and 4, we show the surface for the exact solution and the numerical solution , respectively, with and and obtain a good approximation solution for the exact solution .
The comparison between the exact solution and the numerical solution is described at different times , and with and in Figure 5. The similar comparison for the exact solution and the numerical solution is shown in Figure 6. Figures 5 and 6 show that the solitary wave for and moves to the right with unchanged form and velocity, respectively. Furthermore, the exact solutions and are approximated well by the numerical solutions and , respectively.
5. Concluding Remarks
In this paper, we propose and analyze an explicit multistep mixed finite element method, which combines spatial mixed finite element method and time explicit multistep method, for RLW equation. We discuss the numerical process for our method, prove the theoretical results for the fully discrete explicit multistep mixed scheme, obtain the optimal convergence order, and compare our method’s accuracy with some other numerical schemes. Compared with the numerical method in [20, 25, 26], our method can obtain the optimal error estimates in and norms for the scalar unknown and its flux simultaneously. From our numerical results, we can see that our method is efficient for RLW equation.
This work was supported by the National Natural Science Fund of China (11061021), the Scientific Research Projection of Higher Schools of Inner Mongolia (NJZZ12011 and NJZY13199), the Natural Science Fund of Inner Mongolia Province (2012MS0108 and 2012MS0106), and the Program of Higher-Level Talents of Inner Mongolia University (125119).
- D. H. Peregrine, “Calculations of the development of an undular bore,” Journal of Fluid Mechanics, vol. 25, no. 2, pp. 321–330, 1966.
- T. B. Benjamin, J. L. Bona, and J. J. Mahony, “Model equations for long waves in nonlinear dispersive systems,” Philosophical Transactions of the Royal Society of London A, vol. 272, no. 1220, pp. 47–78, 1972.
- P. J. Olver, “Euler operators and conservation laws of the BBM equation,” Mathematical Proceedings of the Cambridge Philosophical Society, vol. 85, no. 1, pp. 143–160, 1979.
- B. Saka and I. Dağ, “Quartic B-spline collocation algorithms for numerical solution of the RLW equation,” Numerical Methods for Partial Differential Equations, vol. 23, no. 3, pp. 731–751, 2007.
- J. L. Bona and P. J. Bryant, “A mathematical model for long waves generated by wavemakers in non-linear dispersive systems,” Mathematical Proceedings of the Cambridge Philosophical Society, vol. 73, no. 5, pp. 391–405, 1973.
- D. Bhardwaj and R. Shankar, “A computational method for regularized long wave equation,” Computers & Mathematics with Applications, vol. 40, no. 12, pp. 1397–1404, 2000.
- X. H. Zhao, D. S. Li, and D. M. Shi, “A finite difference scheme for RLW-Burgers equation,” Journal of Applied Mathematics & Informatics, vol. 26, no. 3-4, pp. 573–581, 2008.
- L. Zhang, “A finite difference scheme for generalized regularized long-wave equation,” Applied Mathematics and Computation, vol. 168, no. 2, pp. 962–972, 2005.
- I. Dağ, A. Dogan, and B. Saka, “B-spline collocation methods for numerical solutions of the RLW equation,” International Journal of Computer Mathematics, vol. 80, no. 6, pp. 743–757, 2003.
- A. Esen and S. Kutluay, “Application of a lumped Galerkin method to the regularized long wave equation,” Applied Mathematics and Computation, vol. 174, no. 2, pp. 833–845, 2006.
- I. Dağ, B. Saka, and D. Irk, “Galerkin method for the numerical solution of the RLW equation using quintic B-splines,” Journal of Computational and Applied Mathematics, vol. 190, no. 1-2, pp. 532–547, 2006.
- K. R. Raslan, “A computational method for the regularized long wave (RLW) equation,” Applied Mathematics and Computation, vol. 167, no. 2, pp. 1101–1118, 2005.
- L. Wahlbin, “Error estimates for a Galerkin method for a class of model equations for long waves,” Numerische Mathematik, vol. 23, pp. 289–303, 1975.
- Z. D. Luo and R. X. Liu, “Mixed finite element analysis and numerical solitary solution for the RLW equation,” SIAM Journal on Numerical Analysis, vol. 36, no. 1, pp. 89–104, 1999.
- L. Guo and H. Chen, “-Galerkin mixed finite element method for the regularized long wave equation,” Computing, vol. 77, no. 2, pp. 205–221, 2006.
- H. M. Gu and N. Chen, “Least-squares mixed finite element methods for the RLW equations,” Numerical Methods for Partial Differential Equations, vol. 24, no. 3, pp. 749–758, 2008.
- J. L. Bona, W. G. Pritchard, and L. R. Scott, “Numerical schemes for a model for nonlinear dispersive waves,” Journal of Computational Physics, vol. 60, no. 2, pp. 167–186, 1985.
- S. U. Islam, S. Haq, and A. Ali, “A meshfree method for the numerical solution of the RLW equation,” Journal of Computational and Applied Mathematics, vol. 223, no. 2, pp. 997–1012, 2009.
- D. Kaya, “A numerical simulation of solitary-wave solutions of the generalized regularized long-wave equation,” Applied Mathematics and Computation, vol. 149, no. 3, pp. 833–841, 2004.
- L. Q. Mei and Y. P. Chen, “Explicit multistep method for the numerical solution of RLW equation,” Applied Mathematics and Computation, vol. 218, no. 18, pp. 9547–9554, 2012.
- L. R. T. Gardner, G. A. Gardner, and A. Dogan, “A least-squares finite element scheme for the RLW equation,” Communications in Numerical Methods in Engineering, vol. 12, no. 11, pp. 795–804, 1996.
- I. Dağ, “Least squares quadratic B-spline finite element method for the regularized long wave equation,” Computer Methods in Applied Mechanics and Engineering, vol. 182, no. 1-2, pp. 205–215, 2000.
- I. Dağ and M. N. Özer, “Approximation of the RLW equation by the least square cubic B-spline finite element method,” Applied Mathematical Modelling, vol. 25, no. 3, pp. 221–231, 2001.
- A. Dogan, “Numerical solution of RLW equation using linear finite elements within Galerkin’s method,” Applied Mathematical Modelling, vol. 26, no. 7, pp. 771–783, 2002.
- P. Chatzipantelidis, “Explicit multistep methods for nonstiff partial differential equations,” Applied Numerical Mathematics, vol. 27, no. 1, pp. 13–31, 1998.
- G. Akrivis, O. Karakashian, and F. Karakatsani, “Linearly implicit methods for nonlinear evolution equations,” Numerische Mathematik, vol. 94, no. 3, pp. 403–418, 2003.
- A. K. Pani, R. K. Sinha, and A. K. Otta, “An -Galerkin mixed method for second order hyperbolic equations,” International Journal of Numerical Analysis and Modeling, vol. 1, no. 2, pp. 111–129, 2004.
- M. F. Wheeler, “A priori error estimates for Galerkin approximations to parabolic partial differential equations,” SIAM Journal on Numerical Analysis, vol. 10, pp. 723–759, 1973.
- J. L. Bona, W. G. Pritchard, and L. R. Scott, “An evaluation of a model equation for water waves,” Philosophical Transactions of the Royal Society of London A, vol. 302, no. 1471, pp. 457–510, 1981.