Finite Element Method of BBM-Burgers Equation with Dissipative Term Based on Adaptive Moving Mesh

Lu, Changna; Gao, Qianqian; Fu, Chen; Yang, Hongwei

doi:https://doi.org/10.1155/2017/3427376

Discrete Dynamics in Nature and Society

On this page

Abstract Introduction Conclusion Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2017 | Article ID 3427376 | https://doi.org/10.1155/2017/3427376

Finite Element Method of BBM-Burgers Equation with Dissipative Term Based on Adaptive Moving Mesh

Changna Lu,¹Qianqian Gao,¹Chen Fu,¹and Hongwei Yang^2,3

Academic Editor: Zuonong Zhu

Received17 Jul 2017

Revised24 Oct 2017

Accepted06 Nov 2017

Published29 Nov 2017

Abstract

A finite element model is proposed for the Benjamin-Bona-Mahony-Burgers (BBM-Burgers) equation with a high-order dissipative term; the scheme is based on adaptive moving meshes. The model can be applied to the equations with spatial-time mixed derivatives and high-order derivative terms. In this scheme, new variables are needed to make the equation become a coupled system, and then the linear finite element method is used to discretize the spatial derivative and the fifth-order Radau IIA method is used to discretize the time derivative. The simulations of 1D and 2D BBM-Burgers equations with high-order dissipative terms are presented in numerical examples. The numerical results show that the method keeps a second-order convergence in space and provides a smaller error than that based on the fixed mesh, which demonstrates the effectiveness and feasibility of the finite element method based on the moving mesh. We also study the effect of the dissipative terms with different coefficients in the equation; by numerical simulations, we find that the dissipative term plays a more important role than in dissipation.

1. Introduction

In mathematics, physics, and engineering, nonlinear partial differential equations play an important role due to their rich mathematical structures and features [1–4] as well as important applications in fluid dynamics, plasma physics, fiber optics, condensed matter physics, and chemistry [5–11]. It is important to search for solutions. These nonlinear partially differential equations can be solved exactly or numerically to understand the properties of solutions.

In this paper, we consider the nonlinear multidimensional Benjamin-Bona-Mahony-Burgers (BBM-Burgers) equation with a dissipative term and use the finite element method based on the moving mesh to solve its numerical solution. The equation has the following form:where , stands for the dimension, stands for the corresponding spatial vector, the coefficients , , , is a Laplace operator, and is one constant to define the order of the dissipative term. In the equation, is a sufficiently smooth function and stands for the source term, and and are given functions.

Equation (1) is related to the well-known Benjamin-Bona-Mahony (BBM) equation, which was advocated by Benjamin et al. [12] in 1972 as a refinement of the Korteweg-de-Vries (KdV) equation. There are several different transformations of this equation; when the coefficients , , the equation becomes , which is the Burgers equation with a nonlinear term. When parameters , (1) becomes , which is called Regularized Long Wave (RLW) equation, also called BBM equation. The RLW equation is usually used to analyze small amplitude long waves on the surface of water in fluid dynamics. When , (1) becomes a generalized BBM-Burgers equation, which is a mathematical model of propagation of small-amplitude long wave in nonlinear dispersive media and better than Korteweg-de-Vries-Burgers (KdV-B) equation in fluid mechanics; when , , and , the BBM-Burgers equation satisfies two conservation laws:If , , (1) becomes In physical phenomena, the effect of dispersion in the BBM-Burgers equation with a high-order dissipative term is similar to the RLW equation; the effect of dissipation is similar to the Burgers equation.

To get the approximate solution of the BBM-Burgers equation with a high-order dissipative term, some numerical and analytical methods have been proposed in recent years. In the aspect of mathematical theory, Zhao et al. consider the existence of the generalized BBM-Burgers equation with a dissipative term and verify the corresponding convergence in H-measures [13, 14]. Kondo and Webler apply kinetic decomposition for the generalized BBM-Burgers equations with a dissipative term to obtain the approximate transport equation, and then take advantage of the averaging lemma to get the convergence [15]. Seok states the existence and conservation laws of the generalized BBM-Burgers equation with a dissipative term [16]. For more theoretical analyses, refer to the references in [13–16]. In the aspect of numerical solutions, many methods have been used, like Painlev’e test, Darboux Transformation, bilinear method, symmetry method [17], and so on [18–21]. Hong et al. apply the homotopy perturbation method and He’s variational iteration method to solve the BBM-Burgers equation. Ganji et al. put forward the notion that the method they used in [22] is better than the homotopy analysis method proposed in [23]. Finite element Galerkin methods have been discussed by Kadri et al. [24], Lee [25], Yin and Piao [26], and Dehghan et al. [27] to solve the BBM-Burgers equations with a dissipative term. Finite difference methods have also been widely used to solve the equations in [28–30]; as for other numerical methods, refer to [31–34].

The references introducing the numerical simulations above are mostly based on fixed meshes; in this paper, we consider the finite element method based on the adaptive moving mesh to solve the BBM-Burgers equation with a high-order dissipative term. The moving mesh method is a kind of adaptive mesh method which adds elements near the place where the numerical solution changes rapidly and decreases elements in the solution changing slowly; the total element number still remains unchanged. When we compare the finite element method based on adaptive moving mesh with methods based on fixed mesh or proposed in [26], we find that the moving mesh method takes fewer elements to get the same distinguishability and simulates steep waves and the transition of gap distinctly. In recent years, the moving mesh method has witnessed further development and extensive application [35]. This method was applied to approximate the anisotropic diffusion-convection-reaction problems to satisfy a discrete maximum principle, and it was adopted to study the porous medium equation and approximate the solution [36]. For more information about the moving mesh method, please refer to [37–41].

The remainder of this paper is divided into three sections. In Section 2, we introduce the method we used, the discretization of the BBM-Burgers equation with a fourth-order dissipative term in one dimension, and the implementation of the moving mesh method. Section 3 presents 1D and 2D numerical examples of BBM-Burgers equations with a dissipative term. Finally, Section 4 shows the conclusion we draw from the numerical examples.

2. Finite Element Method Based on Moving Mesh

In this section, we mainly introduce the linear finite element method based on moving mesh and make use of the method to solve the BBM-Burgers equation with a high-order dissipative term. For convenience, in Section 2.2, we introduce the finite element method of the one-dimensional equation with a fourth-order dissipative term ; the other parts of this section will focus on two dimensions.

2.1. Adaptive Moving Mesh Method

In this subsection, we mainly want to show the generation of physical mesh at a new time level by an adaptive moving mesh strategy based on the Moving Mesh Partial Differential Equation (MMPDE) [42]; at the current time level, the physical mesh and the numerical solution are known (the mesh and value will be introduced clearly in the following parts). During the strategy, one of the crucially important rules is to view the nonuniform mesh (i.e., the physical mesh in this paper) as a uniform mesh under the metric tensor ; the tensor is a symmetric and positive definite matrix. Under the metric tensor , the physical mesh and the computational mesh can be transformed into each other if and only if it satisfies (the following notations and concepts about the mesh are the same as those in [36])where is an element of the physical mesh , is an element of the computational mesh corresponding to , and and , respectively, show the volume of the corresponding element. Meanwhile, , denotes the determinant of a matrix, and stands for the trace of a matrix. represents the total error of the numerical solution, is the Jacobian matrix of the affine mapping : , and is the average of the metric tensor over the elements . In this paper, we take , where is an identity matrix and stands for a Hessian matrix. Equations (4) and (5) also satisfy the equidistribution principle and alignment principle on moving mesh, which means that, under the metric tensor , each element has the same numerical error and is similar to a reference element ; the reference element is uniform. From (4) and (5), we find that the metric tensor plays an important role in the adaptive moving mesh and controls the adaptation.

If we want to get the mesh which satisfies the two principles, we need to minimize the energy function that contains the metric tensor : where and are dimensionless parameters. When we consider the equation to be just one-dimensional, (6) will have a relationship with (4) and ; when the equation is two-dimensional, parameters are usually defined as and , and stands for the dimension in space. The energy function is also a function about the physical mesh coordinate and the computational mesh coordinate , where stands for the number of mesh vertices. The coordinate is an inverse transformation of ; when or , the coordinate is the same as . The energy function is too long to write; we write it in a simple form:Assume that the current physical mesh and are known; by minimizing the energy function, we can get the corresponding computational mesh of the next time level, which can be expressed by In order to solve the minimum of energy function in (6), we need to take the derivative of the coordinate with respect to the energy function; that is,It is not difficult to find that the equations are nonlinear and difficult to solve directly; in order to solve this problem, many scholars put forward a large number of methods to improve it. In this paper, we use the MMPDE [42] method to get the minimum value of the energy function (6). For convenience and efficiency, - MMPDE is made use of rather than - MMPDE; during the calculation, - only involves the differentiation of with respect to its first two arguments in the energy function, and there is no need for the metric tensor to update. The minimization process can be replaced as the derivative of the computational mesh coordinate with respect to , which is an ordinary differential equation and has the following form:where the derivative is seen as a row vector, is a relaxation factor which is used for adjusting the time scale on moving mesh, and is a positive function used to make (10) have some desired properties above; in this paper, , .

At the current physical mesh and computational solution , we define a reference computational mesh , which has uniform elements; the mesh also has the same relationship, grid connection, and number of mesh vertices as the physical mesh and computational mesh . The three kinds of meshes just differ in the location of the mesh coordinate. During the mesh generation, by integrating (10) from to , we get the new computational mesh ; there is no need for the metric tensor to update. The relationship between the physical mesh and the computational mesh can be expressed by an affine mapping: ; the physical mesh in the next time level can take the reference computational mesh by using a form of linear interpolation to get the approximate result: .

2.2. Linear Finite Element Discretization of the 1D BBM-Burgers Equation with a Dissipative Term

In this subsection, we mainly introduce the linear finite element discretization of the 1D BBM-Burgers equation with a dissipative term. When , , and the discretization of the equation are known, we will get the value at the new time level. The 1D BBM-Burgers equation with a fourth-order dissipative term has the formulationBefore solving (11) by the finite element method based on moving mesh, we have to introduce a new variable rather than solving the equation directly; because a mixed derivative of the finite element approximation based on the moving mesh method is unable to be defined, then (11) becomes a coupled system. Before the discretization of the 2D equation with a fourth-order dissipative term, we need to introduce two new different variables to replace and ; the equation then becomes three coupled equations; the other steps in the 2D case are similar to the 1D case. The new formulation of (11) is as follows: We divide the domain into uniform elements and mesh vertices; the initial physical mesh is . We define the linear basis function of the finite element method as . On a fixed mesh, the basis function is just related to ; on a moving mesh, it is also related to and . The basis function spans the space , . We need to find whether and with the Dirichlet boundary condition equal to zero can satisfy the following system:Take , , where the derivative of the basis function [43] with respect to is , and has the following formulation: (in the two-dimensional equation, ), so can be expressed asTake the basis function and substitute (14) into (13); we can getWe rewrite (15) into a matrix form:where and stand for and and the vector can be expressed as = . On a moving mesh, with the boundary condition, the solution of (11) can be obtained by solving an ordinary differential equation firstly:Only if the ODE (17) is solved can the solution of the BBM-Burgers equation with a dissipative term be expressed, where .

All the process of discretization is based on the homogeneous boundary condition ; when the boundary condition is nonhomogeneous, we just need to subtract the integration at the boundary.

2.3. The Implementation of the Method Combined with Finite Element and Moving Mesh

In this subsection, we mainly introduce the implementation of the two-dimensional method. The ways to solve the equation usually can be divided into two kinds: a simultaneous solution procedure and an alternate solution procedure. The simultaneous solution procedure combines the finite element discretization with the mesh equation (10) into a large system and solves them together; the alternate solution procedure solves one by one. In this paper, we introduce a formulation of the second method which is called MP (the method solves the moving equation firstly and solves the physical equation later), which is more flexible and efficient. The method contains four important parts: initialize the mesh, move the mesh, discretize the equation, and integrate time from to .

Step 1 (initialize the mesh). In the whole computational domain, we first divide the domain into uniform elements; the corresponding initial physical mesh is , and the initial coordinate, initial value of the BBM- Burgers equation with a dissipative term, initial time step, and initial time are known; they are , , , and .

Step 2 (move the mesh). (a) Suppose the current time , the numerical solution of the equation, the coordinate , of the mesh, and the time step are known; (b) we make use of the energy function (9) and ODE function (10) to obtain the coordinate , and then we derive the relationship between the new computational mesh and the physical mesh , where can be expressed by . When the reference computational mesh and are known, we get the new physical mesh approximately; the new coordinate is ; (c) when , using linear interpolation = , the velocity of the corresponding grid node is = . The fifth-order Radau IIA method is used in time discretization; if you are interested in more details, please refer to [44]. From the method, we can determine the time step ; if the real time step , the next time level, the mesh coordinate, and the time step of next time level update: , , and .

Step 3. On the new physical mesh , use the discrete version (16) based on the linear finite element method to solve (12) from time to ; we get the new numerical solution at the new time level; the solution is .

Step 4. When , , go to Step 2; else, stop the computation and save the value of solution.

3. Numerical Examples

Example 1 (a single solitary wave in 1D BBM-Burgers equation). In the one-dimensional model, where the coefficients , , the final time , the exact solution of the equation is , and the source term .

We solve the numerical solution with the method we introduced; with the different element numbers, the corresponding error and the order of convergence for both moving mesh and fixed mesh are reported in Table 1. When the element number is the same, the error based on moving mesh is always smaller than that on fixed mesh, and the convergence order of the two kinds of meshes can lead to the second order in the spatial direction. Comparison between the exact solution and the numerical solution based on moving mesh at the final time is shown in Figure 1. From the figure, we can see that the grid points mainly focus on the domain, where the solution fluctuates greatly, and match the exact solution well. With time going on, the moving mesh trajectory is plotted in Figure 1; the feature of the moving process indicates that most of the points focus on nearby solitary waves and move with the propagation of the solitary wave, which also demonstrates the superiority of the moving mesh method.

(a) Computed solution and exact solution

(b) Mesh trajectories

Example 2 (the elastic collision of three solitary waves in 1D BBM-Burgers equation). In this example, we study the elastic collision of three solitary waves and simulate numerical solutions based on moving mesh and fixed mesh; the BBM-Burgers equation has the following form: with the initial condition , where , , , , , , , and the final time .
In this example, the exact solution of the equation is unknown; we mainly want to show the superiority of the moving mesh. In Figure 2, the conservation laws (2) are plotted; we can see that is close to zero and is a little bigger than zero. The two conservation laws are not conserved by finite element approximation on the moving mesh, because, in this case, and both matrices and are time-dependent (if you want to know more, read [45]). Since the exact solution is unknown, we take the element number on the fixed mesh as a reference solution and plot the numerical solutions of moving mesh and fixed mesh with element number at final time in Figure 3. It can be seen that the numerical solution of fixed mesh is on the left of the reference solution and the wave peaks are shorter than those in the reference solution; there are several small waves in the left of the first wave on the fixed mesh. The result also shows that the numerical solution of moving mesh can match the reference better than that based on the fixed mesh. If we consider the computer memory and computational efficiency, the moving mesh is better for that it moves its mesh points quickly to the place where the numerical solution changes rapidly to get higher distinguishability with fewer elements. The propagation process of the three solitary waves is plotted in Figure 4.

Example 3 (different coefficients in 1D BBM-Burgers equation with a high-order dissipative term). In this case, we mainly consider the following BBM-Burgers equation with a fourth-order dissipative term; the equation is The initial condition , where ; when the coefficients change, the maximum value of the numerical solution changes as well.
In this example, we mainly want to talk about the dissipation effect of the dissipative terms and in the BBM equation. Numerical solutions obtained with different coefficients and one of the coefficients fixed at are shown in Figure 5. When two BBM equations with one dissipative term or have the same dissipation effect or similar wave peak at the final time, the corresponding coefficient of dissipative term is plotted in Figure 6. From these figures, we find that when the coefficient of the second-order dissipative term increases and the coefficient of the fourth-order dissipative term is fixed, the numerical solution decreases a lot; however, when the coefficient increases and the coefficient is fixed, the numerical solution changes a little. The numerical simulation provides three sets of data to show the result in Figure 5: ; when , the corresponding maximum values of the numerical solutions are , ; when , the corresponding values are . A lot of numerical tests have been carried out to adjust and confirm the corresponding coefficient of different dissipative terms to keep the equation with dissipation effect or similar wave peak at . The value of the wave peak is at , and the corresponding data are shown in Figure 6. When , , with time going on, drops steadily; however, when , , the value drops quickly before time , and then the value drops slowly and steadily. We simulate a series of numerical examples with different coefficients; all of these data show us that the dissipative term plays a more important role in dissipation than the dissipative term .

(a) , changes

(b) , changes

(a) ,

(b) ,

Example 4 (two solitary waves in 2D BBM-Burgers equation). This example has two solitary waves and we want to verify the convergence order for a 2D BBM-Burgers equation; the equation has the following formulation: where , , and the exact solution of the equation is , where , , , . To maintain the balance of the equation, = = − − + .

Numerical solutions are shown in Table 2 and Figure 7. From Table 2, we can find that the finite element method leads to the same second order of convergence for both fixed and moving meshes, and the convergence result is similar to Example 1. The error of the numerical solution on the moving mesh is smaller than that on the fixed mesh with the same element number. What is more, with fewer mesh elements, the moving mesh reaches second-order convergence more quickly than the fixed mesh. The trajectory and the numerical solution of both the moving mesh and the fixed mesh at final time can be seen in Figure 7. From the figure, we can see that most points focus on the place where the numerical solution changes a lot.

(a) Moving mesh

(b) The corresponding numerical solution

(c) Fixed mesh

(d) The corresponding numerical solution

4. Conclusion

In this paper, we study the numerical model based on adaptive moving mesh finite element method to simulate 1D and 2D BBM-Burgers equations with a high-order dissipative term. The method can be used in equations with mixed derivatives and high-order derivative terms; on the moving mesh, new variables need to be introduced to transform the equation into coupled equations. A series of numerical examples demonstrate that the finite element method based on moving mesh keeps a second-order convergence in space as that based on fixed mesh. With the same element number, the finite element moving mesh method improves the distinguishability and provides a smaller error than the fixed mesh. Meanwhile, numerical tests indicate that the dissipative term plays a more important role than in dissipation. We compare two BBM equations with different dissipative terms and ; to keep the maximum value of the numerical solutions the same at the final time, the coefficient of is around 50 times that of in the same BBM equation with different dissipative terms in the dissipation aspect.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was partially supported by NSFC (41375115), Open Fund of the Key Laboratory of Meteorological Disaster of the Ministry of Education (Nanjing University of Information Science and Technology) (KLME1507), and Open Fund of the Key Laboratory of Ocean Circulation and Waves, Chinese Academy of Sciences (KLOCAW1401).

References

X. Xu X, “An integrable coupling hierarchy of the Mkdvintegrable systems, its Hamiltonian structure and corresponding nonisospectral integrable hierarchy,” Applied Mathematics and Computation, vol. 216, no. 1, pp. 344–353, 2010.
View at: Google Scholar
X. Li Y, Q. Zhao L, Y. Li X, and etal., “Binary Bargmann symmetry constraint associated with 3×3 discrete matrix spectral problem,” Journal of Nonlinear Sciences and Applications, vol. 8, no. 5, pp. 496–506, 2015.
View at: Google Scholar
Y. Zhang, H. Dong, X. Zhang, and H. Yang, “Rational solutions and lump solutions to the generalized (3+1)-dimensional Shallow Water-like equation,” Computers & Mathematics with Applications, vol. 73, no. 2, pp. 246–252, 2017.
View at: Publisher Site | Google Scholar | MathSciNet
M. P. Barnett, J. F. Capitani, J. von zur Gathen, and J. Gerhard, “Symbolic calculation in chemistry: selected examples,” International Journal of Quantum Chemistry, vol. 100, no. 2, pp. 80–104, 2004.
View at: Publisher Site | Google Scholar
M. F. Saleh, W. Chang, J. C. Travers, P. S. J. Russell, and F. Biancalana, “Plasma-induced asymmetric self-phase modulation and modulational instability in gas-filled hollow-core photonic crystal fibers,” Physical Review Letters, vol. 109, no. 11, Article ID 113902, pp. 1–5, 2012.
View at: Publisher Site | Google Scholar
H. W. Yang, D. Z. Yang, Y. L. Shi, S. S. Jin, and B. S. Yin, “Interaction of algebraic Rossby solitary waves with topography and atmospheric blocking,” Dynamics of Atmospheres and Oceans, vol. 71, pp. 21–34, 2015.
View at: Publisher Site | Google Scholar
H. W. Yang, Z. H. Xu, D. Z. Yang, X. R. Feng, B. S. Yin, and H. H. Dong, “ZK-Burgers equation for three-dimensional Rossby solitary waves and its solutions as well as chirp effect,” Advances in Difference Equations, vol. 2016, article no. 167, 2016.
View at: Publisher Site | Google Scholar
S.-P. Gorza, B. Deconinck, P. Emplit, T. Trogdon, and M. Haelterman, “Experimental demonstration of the oscillatory snake instability of the bright soliton of the (2+1)D hyperbolic nonlinear Schrödinger equation,” Physical Review Letters, vol. 106, no. 9, Article ID 094101, 2011.
View at: Publisher Site | Google Scholar
Y. Shi, B. Yin, H. Yang, D. Yang, and Z. Xu, “Dissipative nonlinear Schrödinger equation for envelope solitary Rossby waves with dissipation effect in stratified fluids and its solution,” Abstract and Applied Analysis, vol. 2014, Article ID 643652, 9 pages, 2014.
View at: Publisher Site | Google Scholar
H. Yang, Q. Zhao, B. Yin, and H. Dong, “A new integro-differential equation for rossby solitary waves with topography effect in deep rotational fluids,” Abstract and Applied Analysis, vol. 2013, Article ID 597807, 8 pages, 2013.
View at: Publisher Site | Google Scholar
H. Yang, B. Yin, Y. Shi, and Q. Wang, “Forced ILW-Burgers equation as a model for Rossby solitary waves generated by topography in finite depth fluids,” Journal of Applied Mathematics, vol. 2012, Article ID 491343, 17 pages, 2012.
View at: Publisher Site | Google Scholar
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, vol. 272, no. 1220, pp. 47–78, 1972.
View at: Publisher Site | Google Scholar | MathSciNet
H. Zhao and B. Xuan, “Existence and convergence of solutions for the generalized BBM-Burgers equations with dissipative term,” Nonlinear Analysis Theory Methods and Applications, vol. 75, no. 1-2, pp. 107–135, 2000.
View at: Publisher Site | Google Scholar | MathSciNet
C. I. Kondo and C. M. Webler, “The generalized BBM-Burgers equations: convergence results for conservation law with discontinuous flux function,” Applicable Analysis: An International Journal, vol. 95, no. 3, pp. 1–21, 2016.
View at: Publisher Site | Google Scholar | MathSciNet
C. I. Kondo and C. M. Webler, “Higher-order for the multidimensional generalized bbm-burgers equation: existence and convergence results,” Acta Applicandae Mathematicae, vol. 88, no. 7, pp. 977–995, 2010.
View at: Publisher Site | Google Scholar
H. Seok, “Kinetic decomposition for the generalized BBM-Burgers equations with dissipative term,” in Proceedings of the Royal Society of Edinburgh, vol. 134, pp. 1149–1162, 2004.
View at: Google Scholar
H. Dong, T. Chen, L. Chen, and Y. Zhang, “A new integrable symplectic map and the Lie point symmetry associated with nonlinear lattice equations,” Journal of Nonlinear Sciences and Applications, vol. 9, no. 7, pp. 5107–5118, 2016.
View at: Google Scholar | MathSciNet
A. Bekir, “Painlevé test for some (2+1)-dimensional nonlinear equations,” Chaos Solitons and Fractals, vol. 32, no. 2, pp. 449–455, 2007.
View at: Google Scholar
X.-X. Xu, “A deformed reduced semi-discrete Kaup-Newell equation, the related integrable family and Darboux transformation,” Applied Mathematics and Computation, vol. 251, pp. 275–283, 2015.
View at: Publisher Site | Google Scholar | MathSciNet
X. Guo, “On bilinear representations and infinite conservation laws of a nonlinear variable-coefficient equation,” Applied Mathematics and Computation, vol. 248, pp. 531–535, 2014.
View at: Publisher Site | Google Scholar | MathSciNet
Q.-L. Zhao, X.-Y. Li, and F.-S. Liu, “Two integrable lattice hierarchies and their respective Darboux transformations,” Applied Mathematics and Computation, vol. 219, no. 10, pp. 5693–5705, 2013.
View at: Publisher Site | Google Scholar | MathSciNet
D. D. Ganji, H. Babazadeh, M. H. Jalaei, and H. Tashakkorian, “Application of He's variational iteration method for solving nonlinear BBMB equations and free vibration of systems,” Acta Applicandae Mathematicae, vol. 106, no. 3, pp. 359–367, 2009.
View at: Publisher Site | Google Scholar | MathSciNet
B. Hong and D. Lu, “Homotopic approximate solutions for the general perturbed burgers-BBM equation,” Journal of Information and Computational Science, vol. 11, no. 11, pp. 4003–4011, 2014.
View at: Publisher Site | Google Scholar
T. Kadri, N. Khiari, F. Abidi, and K. Omrani, “Methods for the numerical solution of the Benjamin-BONa-Mahony-BURgers equation,” Numerical Methods for Partial Differential Equations, vol. 24, no. 6, pp. 1501–1516, 2008.
View at: Publisher Site | Google Scholar | MathSciNet
H. C. Lee, “Distributed feedback control of the Benjamin-Bona-Mahony-Burgers equation by a reduced-order model,” East Asian Journal on Applied Mathematics, vol. 5, no. 1, pp. 61–74, 2015.
View at: Google Scholar
Y. Yin and G. R. Piao, “Quadratic B-spline finite element method for the generalized Benjamin-Bona-Mahony-Burgers equation,” Journal of the Korean Mathematical Society, vol. 29, no. 5, 2013.
View at: Google Scholar
M. Dehghan, M. Abbaszadeh, and A. Mohebbi, “The use of interpolating element-free Galerkin technique for solving 2D generalized Benjamin-BONa-Mahony-BURgers and regularized long-wave equations on non-rectangular domains with error estimate,” Journal of Computational and Applied Mathematics, vol. 286, Article ID 10067, pp. 211–231, 2015.
View at: Publisher Site | Google Scholar | MathSciNet
K. Omrani and M. Ayadi, “Finite difference discretization of the Benjamin-BONa-Mahony-BURgers equation,” Numerical Methods for Partial Differential Equations, vol. 24, no. 1, pp. 239–248, 2008.
View at: Publisher Site | Google Scholar | MathSciNet
A. Mohebbi and Z. Faraz, “Solitary wave solution of nonlinear Benjamin-BONa-Mahony-BURgers equation using a high-order difference scheme,” Computational & Applied Mathematics, pp. 1–13, 2015.
View at: Google Scholar | MathSciNet
H.-T. Che, X.-T. Pan, L.-M. Zhang, and Y.-J. Wang, “Numerical analysis of a linear-implicit average scheme for generalized Benjamin-Bona-Mahony-Burgers equation,” Journal of Applied Mathematics, vol. 2012, Article ID 308410, pp. 295–305, 2012.
View at: Publisher Site | Google Scholar
G. Arora, R. C. Mittal, and B. K. Singh, “Numerical solution of BBM-Burger equation with quartic B-spline collocation method,” Journal of Engineering Science and Technology, vol. 9, pp. 104–116, 2014.
View at: Google Scholar
M. Alquran and K. Al-Khaled, “Sinc and solitary wave solutions to the generalized Benjamin-Bona-Mahony- Burgers equations,” Physica Scripta, vol. 83, no. 6, Article ID 065010, 2011.
View at: Publisher Site | Google Scholar
Ö. Oruç, F. Bulut, and A. Esen, “Numerical solutions of regularized long wave equation by haar wavelet method,” Mediterranean Journal of Mathematics, vol. 13, no. 5, pp. 1–19, 2016.
View at: Google Scholar
Y. Fang, H. Dong, Y. Hou, and Y. Kong, “Frobenius integrable decompositions of nonlinear evolution equations with modified term,” Applied Mathematics and Computation, vol. 226, no. 1, pp. 435–440, 2014.
View at: Publisher Site | Google Scholar | MathSciNet
C. Lu, W. Huang, and J. Qiu, “Maximum principle in linear finite element approximations of anisotropic diffusion-convection-reaction problems,” Numerische Mathematik, vol. 127, no. 3, pp. 515–537, 2014.
View at: Publisher Site | Google Scholar | MathSciNet
C. Ngo and W. Huang, “A study on moving mesh finite element solution of the porous medium equation,” Journal of Computational Physics, vol. 331, pp. 357–380, 2017.
View at: Publisher Site | Google Scholar | MathSciNet
W. Huang and R. D. Russell, Adaptive Moving Mesh Methods, vol. 174 of Applied Mathematical Sciences, Springer Science & Business Media, New York, NY, USA, 2011.
View at: Publisher Site | MathSciNet
R. Li, T. Tang, and P. Zhang, “Moving mesh methods in multiple dimensions based on harmonic maps,” Journal of Computational Physics, vol. 170, no. 2, pp. 562–588, 2001.
View at: Publisher Site | Google Scholar | MathSciNet
C. J. Budd, W. Huang, and R. D. Russell, “Adaptivity with moving grids,” Acta Numerica, vol. 18, pp. 111–241, 2009.
View at: Publisher Site | Google Scholar | MathSciNet
T. Tang, “Moving mesh methods for computational fluid dynamics,” in Recent Advances in Adaptive Computation, vol. 383 of Contemporary Mathematics, pp. 141–174, 2005.
View at: Publisher Site | Google Scholar | MathSciNet
W. Huang and L. Kamenski, “A geometric discretization and a simple implementation for variational mesh generation and adaptation,” Journal of Computational Physics, vol. 301, pp. 322–337, 2015.
View at: Publisher Site | Google Scholar | MathSciNet
W. Huang, Y. Ren, and R. D. Russell, “Moving mesh partial differential equations (MMPDES) based on the equidistribution principle,” SIAM Journal on Numerical Analysis, vol. 31, no. 3, pp. 709–730, 1994.
View at: Publisher Site | Google Scholar | MathSciNet
P. K. Jimack and A. J. Wathen, “Temporal derivatives in the finite-element method on continuously deforming grids,” SIAM Journal on Numerical Analysis, vol. 28, no. 4, pp. 990–1003, 1991.
View at: Publisher Site | Google Scholar | MathSciNet
S. Gonzalez-Pinto, J. Montijano I, and S. Perez-Rodriguez, “Two-step error estimators for implicit Runge-Kutta methods applied to stiff systems,” Acm Transactions on Mathematical Software, vol. 30, pp. 1–18, 2004.
View at: Google Scholar
C. Lu, W. Huang, and J. Qiu, “An adaptive moving mesh finite element solution of the regularized long wave equation,” Journal of Scientific Computing, pp. 1–23, 2017.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2017 Changna Lu 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.

PDF Download Citation

Download other formats

Order printed copies

Views

1886

Downloads

1513

Citations