The Boussinesq equation was first presented in the 1870s by Joseph Boussinesq to model the propagation of shallow water waves in multiple directions. It is a fourth-order nonlinear partial differential equation and belongs to KdV family, which has been widely used in math-physical field related to nonlinear wave phenomena, such as ion-sound in plasma, nonlinear lattice waves . The Boussinesq equation is also used to study the disturbance and repair of solitary waves in blood vessels with aneurysm  . A general form of the Boussinesq equation is
where is a sufficiently differentiable function, which for gives the good Boussinesq or well-posed equation, while for the bad or ill-posed Boussinesq equation.
In this study, we will take . If the term in Equation (1) is replaced with , it gives the so-called improved Boussinesq equation
Researchers do not observe any instability for the above improved Boussinesq equation . The improved Boussinesq equation (IBq) has been investigated theoretically and numerically by many people.
Yang  studied the existence, both locally and globally in time, the uniqueness of solutions and non-existence of global solutions to the initial boundary value problem of a generalized IBq equation. Iskandar and Jain  studied the dynamical behavior of IBq equation numerically. El-Zoheiry  designed a three-level iterative scheme based on the compact implicit method for solving the improved Boussinesq equation. Bratsos  proposed a second order implicit finite difference scheme by rational approximations and constructed a predictor-corrector scheme  to avoid the solution of a nonlinear system. Lin et al.  used the finite element method with linear B-spline basis functions to discretize the improved Boussinesq equation in space and derived a second-order system. Shokri and Dehghan  used collocation and approximated the solution by radial basis functions based on the third-order time discretization to solve the equation. Irk and Dăg  obtained numerical simulations of the improved Boussinesq equation by using two finite difference schemes and two finite element methods, based on the second- and third-order time discretization. Zhang and Lu  derived quadratic finite volume element schemes to solve the improved Boussinesq equation.
So far the main attention was paid to the improved Boussinesq equation in the absence of damping. However, the damping will change the amplitude of the waves and their shape. Yan et al.  studied the high-order modified Boussinesq equation with damping term. In , Arevalo et al. investigated the dynamics of a lattice soliton in the presence of damping. By using a multiple-scale perturbation expansion up to second-order, they derived a general expression for the first-order velocity correction. Chen et al.  studied the Cauchy problem for the generalized IBq equation with hydrodynamical damping.
The improved Boussinesq equation of interest is here
where is the positive constant. According to , the damping in the equation corresponds to the Stokes damping. For the improved Boussinesq equation with Stokes damping, we adopt the following periodic boundary conditions
As the analytic solutions for the Equation (3) are not available, it is necessary to design efficient numerical schemes to perform numerical simulation. To investigate the influence of damping on the solution, the proposed schemes should inherit characteristic properties of the system. The finite volume element method has been one of the most commonly used numerical methods for solving partial differential equations . One main attractive property of the method is that, the mass conservation law is maintained, which is fairly desirable for fluid and underground fluid computations. So it has been extensively used in computational fluid dynamics  - .
In this study, we will propose new finite volume element schemes for the improved Boussinesq equation with Stokes damping. The schemes can preserve mass conservation and energy dissipation of the improved Boussinesq equation. Furthermore, we will study numerically the rate of decay of solutions of the considered equation using our numerical schemes.
The paper is organized as follows: Our numerical method is described in Section 2, which will give good results for deterministic cases. In Section 3, we present some numerical examples to illustrate the accuracy and conservation property of the proposed numerical schemes. In Section 4, the numerical results for the improved Boussinesq equation with damping are given. The conclusions are summarized in the last section.
2. Numerical Method
In this section, new finite volume element schemes will be derived to solve the considered Equation (3).
2.1. Energy Dissipation Property
In order to see the energy dissipation property, it is more convenient to work with the following system representation
where is the Euler-Lagrange variational derivative   defined by
Then, the improved Boussinesq equation with Stokes damping has the following mass conservation property.
Proposition 1. Let u be the solution of Equation (5). Then the total mass is independent of t. Namely,
The conservation of mass can be easily proved with periodic boundary conditions.
Proposition 2. Let u be the solution of Equation (5). Then the total energy is independent of t. Namely,
Proof. According to the chain rule, we have
By the first equality of Equation (5), we can get
Using the integration by part formula, we have
At last, according to the second equality of Equation (5), we can complete the proof
2.2. Trial and Test Function Spaces
The region of interest can be decomposed into a grid with nodes
where and . Denote . Accordingly, we place a dual grid with nodes
where Denote , and .
Select the trial function space as the linear element space with respect to . The basis function of the node is
Then numerical solution U for Equation (5) can be uniquely written as , where . So in the element , we have
The test function space corresponding to is taken as the piecewise constant function space. The test function of the nodes is
In addition, we take as the time step and the numerical solution can be denoted by .
Recall the definition . According to the definition of Euler-Lagrange derivative, we have , and . So the Equation (5) can be rewritten as
Using any and to multiply both sides of Equation (14) and integrating on the interval , we can obtain the following semi-discrete scheme
where denotes the inner product of .
Now we discretize the above semi-discrete scheme in time. To simplify the notation, we will use the difference operator . Additionally, we define the discrete version of Euler-Lagrange derivative , and by
At last, we can get the following energy dissipation finite volume element schemes
Under the periodic boundary conditions, the Equation (19) has the following discrete mass conservation law and energy dissipation property.
Theorem 1. (Discrete mass conservation law). Under the discrete periodic boundary condition, then numerical solution by Equation (19) conserves mass
Theorem 2. (Discrete energy dissipation property). Under the discrete periodic boundary condition, then numerical solution by Equation (19) conserves energy
2.3. Finite Volume Element Schemes
For comparison, we employ the standard finite volume element scheme constructed based on the
Then we can obtain the following finite volume element schemes
Theorem 3. (Discrete mass conservation law). Under the discrete periodic boundary condition, then numerical solution by Equation (23) conserves mass
3. Validation of the Numerical Schemes
In this section, we present some numerical examples to illustrate the accuracy and conservation property of the proposed numerical schemes.
The initial conditions are taken as
where is amplitude of the solitary wave which is initially centered at travelling with velocity . In addition, we will take , , and the region in the following.
As the exact solution for the considered examples is not available, we will compute maximum absolute error for the considered examples using the half mesh principle. The maximum absolute error is defined as follows,
The maximum absolute errors Equation (3) and solitary wave with the above parameters have been computed at t= 1. They are displayed in Table 1 for the new finite volume element schemes (NFVEM). Examination of the table shows that the error measures of the finite volume element scheme diminish approximately quadratically as the space step size and time step size are simultaneously halved. The error table illustrates that the new numerical schemes have second-order accuracy in space and time direction.
Furthermore, we compute the problem using standard finite volume element schemes (FVEM). The obtained maximum absolute errors are displayed in Table 2. From the table, we can find that the standard finite volume element schemes also have second-order accuracy in space and time direction. The maximum absolute errors computed by the two numerical schemes are consistent. In other words, the proposed new schemes do not reduce the accuracy of the finite volume method.
3.2. Conservation Properties
Now we consider the conservation properties of the new finite volume element schemes and standard finite volume element schemes. The initial conditions and parameters are the same as the above except .
Figure 1 shows the evolutions of discrete mass of the proposed new finite volume element schemes and standard finite volume element schemes with .
Table 1. The maximum absolute error for the NFVEM at t = 1.
Table 2. The maximum absolute error for the FVEM at t = 1.
Figure 1. The evolutions of the discrete mass: (a) FVEM; (b) NFVEM. .
Figure 2. The evolutions of the discrete energy: (a) FVEM; (b) NFVEM.
From the figure, we can find that the discrete mass is well conserved. Actually it is conserved up to the machine accuracy, which well agrees with the discrete mass conservation law. In other words, the Stokes damping will not change the mass of system.
Figure 2 presents the evolutions of discrete energy of the proposed new finite volume element schemes and standard finite volume element schemes with . It is easy to find that the discrete energy obtained by standard finite volume element schemes decreases with oscillations. But the discrete energy obtained by new finite volume element schemes will decrease continuously. So our new finite volume element schemes can preserve the energy dissipation property of the improved Boussinesq equation with Stokes damping. From this section, we can conclude that the proposed numerical schemes are more suitable for solving improved Boussinesq equation with Stokes damping.
In this paper, we have proposed a new finite volume element method for solving the improved Boussinesq equation with Stokes damping. The proposed method is second-order accurate in time and space. Furthermore, the new method can preserve the mass conservation and energy dissipation property.
 Dehghan, M. and Shakeri, F. (2008) Use of He’s Homotopy Perturbation Method for Solving a Partial Differential Equation Arising in Modeling of Flow in Porous Media. Journal of Porous Media, 11, 765-778.
 Biswas, A., Milovic, D. and Ranasinghe, A. (2009) Solitary Waves of Boussinesq Equation in a Power Law Media. Communications in Nonlinear Science and Numerical Simulation, 14, 3738-3742.
 Christiansen, P.L., Muto, V. and Rionero, S. (1992) Solitary Wave Solutions to a System of Boussinesq-Like Equations. Chaos, Solitons and Fractals, 2, 45-50.
 Yang, Z. (1998) Existence and Non-Existence of Global Solutions to a Generalized Modification of the Improved Boussinesq Equation. Mathematical Methods in the Applied Sciences, 21, 1467-1477.
 Iskandar, L. and Jain, P.C. (1980) Numerical Solutions of the Improved Boussinesq Equation. Proceedings of the Indian Academy of Sciences-Mathematical Sciences, 89, 171-181.
 Lin, Q., Wu, Y., Loxton, R. and Lai, S. (2009) Linear B-Spline Finite Element Method for the Improved Boussinesq Equation. Journal of Computational and Applied Mathematics, 224, 658-667.
 Shokri, A. and Dehghan, M. (2010) A Not-a-Knot Meshless Method Using Radial Basis Functions and Predictor—Corrector Scheme to the Numerical Solution of Improved Boussinesq Equation. Computer Physics Communications, 181, 1990-2000.
 Zhang, Z. and Lu, F. (2012) Quadratic Finite Volume Element Method for the Improved Boussinesq Equation. Journal of Mathematical Physics, 53, Article ID: 013505.
 Yan, Z., Xie, F. and Zhang, H. (2001) Symmetry Reduction, Integrability and Solitary Wave Solutions to High-Order Modified Boussinesq Equation with Damping term. Communications in Theoretical Physics, 36, 1-6.
 Zhang, Z. (2008) Error Estimates of Finite Volume Element Method for the Pollution in Groundwater Flow. Numerical Methods for Partial Differential Equations, 25, 259-274.
 Li, J., Chen, Z. and He, Y. (2012) A Stabilized Multi-Level Method for Non-Singular Finite Volume Solutions of the Stationary 3D Navier-Stokes Equations. Numerische Mathematik, 122, 279-304.
 Luo, Z., Li, H., Sun, P., An, J. and Navon, I. (2013) A Reduced-Order Finite Volume Element Formulation Based on POD Method and Numerical Simulation for Two-Dimensional Solute Transport Problems. Mathematics and Computers in Simulation, 89, 50-68.
 Wang, Q., Zhang, Z. and Li, Z. (2013) A Fourier Finite Volume Element Method for Solving Two-Dimensional Quasi-Geostrophic Equations on a Sphere. Applied Numerical Mathematics, 71, 1-13.
 Wang, Q. and Zhang, Z. (2014) High-Order Upwind Finite Volume Element Schemes for Modeling of Neuronal Firing. International Journal of Computer Mathematics, 91, 625-640.
 Lin, Z., Ricardo, R. and Tian, C. (2014) Finite Volume Element Approximation of an Inhomogeneous Brusselator Model with Cross-Diffusion. Journal of Computational Physics, 256, 806-823.
 Matsuo, T. (2007) New Conservative Schemes with Discrete Variational Derivatives for Nonlinear Wave Equations. Journal of Computational and Applied Mathematics, 203, 32-56.