Variational Iteration Method for Solving Boussinesq Equations Using Maple

Abstract

In this study, we applied the variational iteration method to solve the Boussinesq time equation. Bossiness’s article from 1872 introduced the equations that are now known as the Boussinesq equations. Numerical methods are commonly utilized to solve nonlinear equation systems. Several research papers have documented the values of the variational iteration method and its applications for various categories of differential equations. A comparison of the exact and numerical solutions was obtained using the variational iteration method. The variational iteration method shows that the proposed method is very effective and convenient. The results are shown for different specific cases of the problem. The variational iteration method is useful in numerical simulations and approximate analytical solutions, and it is used to resolve nonlinear differential equations in various situations using Maple. For example, the linear Boussinesq equation was resolved using the variational iteration method. By comparing the numerical results, we found that the variable repetition method produced accurate results and was close to the exact solution, allowing it to be widely applied to the Boussinesq equation. This proves the effectiveness of the method and the capability to quickly and effectively obtain the numerical number solution related to the exact solution using the Maple 18 program. Additionally, the outcomes are extremely precise.

Keywords

Share and Cite:

Aljuhani, A. , Maturi, D. and Alshehri, H. (2022) Variational Iteration Method for Solving Boussinesq Equations Using Maple. Applied Mathematics, 13, 960-967. doi: 10.4236/am.2022.1312060.

1. Introduction

The variational iteration method (VIM) was proposed by Ji-Huan He in 1999   . This method has been used in many studies to solve different types of equations. Abdul-Majid Wazwaz used it to determine rational solutions for the KdV, the K(2, 2), the Burgers, and the cubic Boussinesq equations  . M. Javidi constructed solitary wave solutions for the Boussinesq equation using VIM . In addition, the VIM has been proven by numerous authors to be an effective mathematical tool for resolving various nonlinear problems. Tamer Abassy applied the modified VIM to solve the Boussinesq equation . S. A. El-Wakil also used the improved VIM to solve Boussinesq-like B(m, n) equations . It was successfully applied to solve complex ordinary differential equation problems , partial differential equations , and fractional nonlinear differential equations . E. Rama presented a study for solving various types of problems . Elcin deals with the blow-up solutions of the generalized Boussinesq equation  and uses effective programs (Maple) to construct a numerical simulation of the Boussinesq equations . Many applications have been introduced in the literature    . Numerical computer programs are developed for these problems, and corresponding examples are provided to verify the accuracy and efficiency of the simplified VIM in solving the Boussinesq equation.

VIM effectively analyzes close-form analytical or approximate solutions of nonlinear evolution equations. If there is an exact solution, this method provides fast, convergent, and successive approximations of it; otherwise, approximations can be used for numerical purposes. Another significant benefit of VIM is its ability to drastically reduce calculation size while retaining the high accuracy of numerical results.

Scientific fields as diverse as fluid dynamics, solid-state physics, plasma physics, and chemical kinetics use nonlinear phenomena. The examination of these scientific models utilized various analytical and numerical techniques because of the rising interest in the theory of solitary waves. The classical Boussinesq equation,

${u}_{tt}={\left(u+{u}^{2}+{u}_{xx}\right)}_{xx}$ (1)

which was derived in 1872 to describe shallow water waves, has the flaw that the Cauchy problem is improperly posed. Therefore, it cannot be used to analyze numerical wave propagation issues.

The equation also occurs in various physical applications, such as vibrations in a nonlinear string, iron sound waves in plasma, and nonlinear lattice waves. Additionally, it was used to address issues with water percolation in porous subsurface strata. Recently, certain novel approaches to solving nonlinear equations, such as the VIM, have garnered considerable interest.

Boussinesq put forth a well-known model of nonlinear dispersive waves in a generalized form:

${u}_{tt}={\left[f\left(u\right)\right]}_{xx}+{u}_{xxxx}+h\left(x,t\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}-\infty 0,$ (2)

where u, f, and h are sufficiently differentiable functions and $f\left(0\right)=0$. The initial conditions associated with the Boussinesq Equation (1) have the following form:

$u\left(x,0\right)=a\left(x\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}{u}_{t}\left(x,0\right)=b\left(x\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}-\infty (3)

With $a\left(x\right)$ and $b\left(x\right)$ given ${C}^{\infty }$.

This paper aims to solve the Boussinesq equation and compare the exact and numerical solutions obtained using the method of variable iteration by the Maple 18 program.

To present a clear overview of method, we have chosen two examples, to illustrate the variational iteration method and the obtained solutions are compared with the exact solutions.

2. He’s Variational Iteration Method

The variational iteration method is easy to compute the successive approximations. In a fewer iterations we get exact solutions of problems if solution of those problems exist.

To illustrate the methodology of the proposed method, using the VIM, we consider

$Lu+Nu=g\left(x\right),$ (4)

where L is a linear operator, N is a nonlinear operator, and $g\left(x\right)$ is a known analytic function. According to the VIM, we can construct a correction functional as follows:

${u}_{n+1}\left(x\right)={u}_{n}\left(x\right)+{\int }_{0}^{x}\lambda \left[L{u}_{n}\left(\xi \right)+{\left(N{u}_{n}\left(\xi \right)\right)}^{-}-g\left(\xi \right)\right]\text{d}\left(\xi \right)$ (5)

where $\lambda$ is a general Lagrangian multiplier, which can be identified optimally via the variational theory, the subscript n denotes the nth order approximation, and ${u}_{n}^{-}$ is considered a restricted variation, i.e., $\delta {\left[{u}_{n}\right]}^{-}=0$.

Equation (5) is called a correction functional. He’s VIM has been shown to solve a large class of nonlinear problems effectively, easily, and accurately, with approximations converging rapidly to accurate solutions.

For linear problems, its exact solution can be obtained by only one iteration step as the Lagrange multiplier can be exactly identified.

3. Application

Two examples are given in this section to illustrate the effects of the proposed method.

3.1. Example 1

We consider the cubic Boussinesq equation:

${u}_{tt}-{u}_{xx}+2{\left({u}^{3}\right)}_{xx}-{u}_{xxxx}=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}u\left(x,0\right)=\frac{1}{x},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{u}_{t}\left(x,0\right)=-\frac{1}{{x}^{2}}$ (6)

The exact solution is:

$u\left(x,t\right)=\frac{1}{x+t}$

The correction functional for (6) reads

$\begin{array}{c}{u}_{n+1}\left(x,t\right)={u}_{n}\left(x,t\right)+{\int }_{0}^{t}\lambda \left(\xi \right)\left(\frac{{\partial }^{2}{u}_{n}\left(x,\xi \right)}{\partial {\xi }^{2}}-\frac{{\partial }^{2}\left({\stackrel{˜}{u}}_{n}\right)\left(x,\xi \right)}{\partial {x}^{2}}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+2\frac{{\partial }^{2}{\left({u}_{n}\left(x,\xi \right)\right)}^{3}}{\partial {x}^{2}}+\frac{{\partial }^{4}{u}_{n}\left(x,\xi \right)}{\partial {x}^{4}}\right)\text{d}\xi \text{.}\end{array}$ (7)

This yields the stationary conditions

$\begin{array}{l}1-{\lambda }^{\prime }=0,\\ \lambda \left(\xi =t\right)=0,\\ {\lambda }^{″}=0.\end{array}$ (8)

This in turn gives

$\lambda =\xi -t\text{.}$ (9)

Substituting this value of the Lagrangian multiplier into functional (7) gives the iteration formula

$\begin{array}{c}{u}_{n+1}\left(x,t\right)={u}_{n}\left(x,t\right)+{\int }_{0}^{t}\left(\xi -t\right)\left(\frac{{\partial }^{2}{u}_{n}\left(x,\xi \right)}{\partial {\xi }^{2}}-\frac{{\partial }^{2}\left({\stackrel{˜}{u}}_{n}\right)\left(x,\xi \right)}{\partial {x}^{2}}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+2\frac{{\partial }^{2}{\left({u}_{n}\left(x,\xi \right)\right)}^{3}}{\partial {x}^{2}}+\frac{{\partial }^{4}{u}_{n}\left(x,\xi \right)}{\partial {x}^{4}}\right)\text{d}\xi \text{.}\end{array}$ (10)

The given initial values admit the use of

${u}_{0}\left(x,t\right)=\frac{1}{x}-\frac{t}{{x}^{2}}\text{.}$ (11)

Using (10) we obtain the following successive approximations:

$\begin{array}{l}{u}_{0}\left(x,t\right)=\frac{1}{x}-\frac{t}{{x}^{2}},\\ {u}_{1}\left(x,t\right)=\frac{1}{x}-\frac{t}{{x}^{2}}+\frac{{t}^{2}}{{x}^{3}}-\frac{{t}^{3}}{{x}^{4}}+\text{smallterms},\\ {u}_{2}\left(x,t\right)=\frac{1}{x}-\frac{t}{{x}^{2}}+\frac{{t}^{2}}{{x}^{3}}-\frac{{t}^{3}}{{x}^{4}}+\frac{{t}^{4}}{{x}^{5}}-\frac{{t}^{6}}{{x}^{7}}+\text{smallterms},\\ {u}_{3}\left(x,t\right)=\frac{1}{x}-\frac{t}{{x}^{2}}+\frac{{t}^{2}}{{x}^{3}}-\frac{{t}^{3}}{{x}^{4}}+\frac{{t}^{4}}{{x}^{5}}-\frac{{t}^{6}}{{x}^{7}}+\frac{{t}^{7}}{{x}^{8}}-\frac{{t}^{8}}{{x}^{9}}+\text{smallterms},\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}⋮\\ {u}_{n}\left(x,t\right)=\frac{1}{x}-\frac{t}{{x}^{2}}+\frac{{t}^{2}}{{x}^{3}}-\frac{{t}^{3}}{{x}^{4}}+\frac{{t}^{4}}{{x}^{5}}-\frac{{t}^{6}}{{x}^{7}}+\frac{{t}^{7}}{{x}^{8}}-\frac{{t}^{8}}{{x}^{9}}+\cdots \end{array}$ (12)

And in a closed form by

$u\left(x,t\right)=\frac{1}{x+t}\text{.}$ (13)

3.2. Example 2

We consider the cubic Boussinesq equation:

${u}_{tt}={u}_{xx}+3{\left({u}^{2}\right)}_{xx}+{u}_{xxxx},$ (14)

$u\left(x,0\right)=2\frac{a{k}^{2}{\text{e}}^{kx}}{{\left(1+a{\text{e}}^{kx}\right)}^{2}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{u}_{t}\left(x,0\right)=-2\frac{a{k}^{3}\sqrt{1+{k}^{2}}{\text{e}}^{kx}\left(a{\text{e}}^{kx}-1\right)}{{\left(1+a{\text{e}}^{kx}\right)}^{3}}$

The exact solution is:

$u\left(x\right)=2\frac{{\text{e}}^{kx+\sqrt{1+{k}^{2}}t}}{{\left(1+{\text{e}}^{kx+k\sqrt{1+{k}^{2}t}}\right)}^{2}}$

By same way in example 1 we solved example 2.

Figure 1 and Figure 2 show the exact and approximate solutions. This problem was solved by VIM and their results are shown in Table 1 and Table 2 using maple.

Figure 1. Graph showing the correspondence between exact solution and numerical solution in Example 1.

Figure 2. Graph showing the correspondence between exact solution and numerical solution in Example 2.

Table 1. Comparison of VIM with exact solution for Example 1.

Table 2. Comparison of VIM with exact solution for Example 2.

4. Conclusion

The VIM was used to solve the Boussinesq equations using Maple 18 software. The results were created using tables and figures. Table 1 shows the numerical solution and the right solution. We can see that the numerical solution is generally relevant to the precise answer by comparing the numerical results, proving the method’s efficacy and the ability to swiftly and easily obtain the numerical solution relating to the exact solution using Maple 18 software. Furthermore, the results obtained are quite precise.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

  He, J.-H. and Wu, X.-H. (2007) Variational Iteration Method: New Development and Applications. Computers & Mathematics with Applications, 54, 881-894.https://doi.org/10.1016/j.camwa.2006.12.083  He, J.H. (1999) Variational Iteration Method—A Kind of Nonlinear Analytical Technique: Some Examples. International Journal of Nonlinear Mechanics, 34, 699-708.https://doi.org/10.1016/S0020-7462(98)00048-1  He, J.H. (2007) Variational Iteration Method—Some Recent Results and New Interpretations. Journal of Computational and Applied Mathematics, 207, 3-17.https://doi.org/10.1016/j.cam.2006.07.009  Wazwaz, A.M. (2011) Linear and Nonlinear Integral Equations. Springer, Berlin, 35-36.https://doi.org/10.1007/978-3-642-21449-3  Wazwaz, A.M. (2007) The Variational Iteration Method for Rational Solutions for KdV, K(2, 2), Burgers, and Cubic Boussinesq Equations. Journal of Computational and Applied Mathematics, 207, 18-23.https://doi.org/10.1016/j.cam.2006.07.010  Javidi, M. and Jalilian, Y. (2008) Exact Solitary Wave Solution of Boussinesq Equation by VIM. Chaos, Solitons & Fractals, 36, 1256-1260.https://doi.org/10.1016/j.chaos.2006.07.046  Abassy, T.A., El-Tawil, M.A. and El-Zoheiry, H. (2007) Modified Variational Iteration Method for Boussinesq Equation. Computers & Mathematics with Applications, 54, 955-965.https://doi.org/10.1016/j.camwa.2006.12.040  El-Wakil, S.A., Abulwafa, E.M. and Abdou, M.A. (2009) An Improved Variational Iteration Method for Solving Coupled KdV and Boussinesq-Like B(m, n) Equations. Chaos, Solitons & Fractals, 39, 1324-1334.https://doi.org/10.1016/j.chaos.2007.05.020  Pan, C., Zhang, R., Luo, H. and Shen, H. (2016) Simplified Variational Iteration Method for Solving Ordinary Differential Equations and Eigenvalue Problems. Advances in Mechanical Engineering, 8, 1-13.https://doi.org/10.1177/1687814016681469  Elham, S. and Hossein, J. (2011) A Modified Variational Iteration Method for Solving Generalized Boussinesq Equation and Lieâ€™ Nard Equation. International Journal of Physical Sciences, 6, 5406-5411.  Wu, G.C. (2011) A Fractional Variational Iteration Method for Solving Fractional Nonlinear Differential Equations. Computers & Mathematics with Applications, 61, 2186-2190.https://doi.org/10.1016/j.camwa.2010.09.010  Rama, E., Somaiah, K. and Sambaiah, K. (2021) A Study of Variational Iteration Method for Solving Various Types of Problems. Malaya Journal of Matematik, 9, 701-708.https://doi.org/10.26637/MJM0901/0123  Yusufoğlu, E. (2008) Blow-Up Solutions of the Generalized Boussinesq Equation Obtained by Variational Iteration Method. Nonlinear Dynamics, 52, 395-402.https://doi.org/10.1007/s11071-007-9287-1  Yusufoglu, E. and Bekir, A. (2009) Numerical Simulations of the Boussinesq Equation by He’s Variational Iteration Method. International Journal of Computer Mathematics, 86, 676-683.https://doi.org/10.1080/00207160701689542  Abdou, M.A. and Soliman, A.A. (2005) New Applications of Variational Iteration Method. Physica D: Nonlinear Phenomena, 211, 1-8.https://doi.org/10.1016/j.physd.2005.08.002  Alavi, A. and Heidari, A. (2012) Parametric Iteration Method for Solving Linear Optimal Control Problems.Applied Mathematics, 3, 1059-1064.https://doi.org/10.4236/am.2012.39155  Othman, M.I.A. and Mahdy, A.M.S. (2010) Differential Transformation Method and Variation Iteration Method for Cauchy Reaction-Diffusion Problems. Journal of Mathematics and Computer Science, 1, 61-75.https://doi.org/10.22436/jmcs.001.02.01  Elsheikh, A.M. and Elzaki, T.M. (2015) Variation Iteration Method for Solving Porous Medium Equation. International Journal of Development Research, 5, 4677-4680.     customer@scirp.org +86 18163351462(WhatsApp) 1655362766  Paper Publishing WeChat 