Variational Iteration Method for Solving Time Fractional Burgers Equation Using Maple

Abstract

The Time Fractional Burger equation was solved in this study using the Mabel software and the Variational Iteration approach. where a number of instances of the Time Fractional Burger Equation were handled using this technique. Tables and images were used to present the collected numerical results. The difference between the exact and numerical solutions demonstrates the effectiveness of the Mabel program’s solution, as well as the accuracy and closeness of the results this method produced. It also demonstrates the Mabel program’s ability to quickly and effectively produce the numerical solution.

Share and Cite:

Alwehebi, F. , Hobiny, A. and Maturi, D. (2023) Variational Iteration Method for Solving Time Fractional Burgers Equation Using Maple. Applied Mathematics, 14, 336-348. doi: 10.4236/am.2023.145021.

1. Introduction

The area of mathematics that focuses on fractional calculus examines the concepts and strategies for employing fractional derivatives to solve differential equations as well as the properties of non-integer order fractional derivatives and integrals. Applications of science in fields like physics, chemistry, and engineering are particularly exciting (fluid flow, viscoelasticity, electrical networks, optics and signal processing, etc). Fractional calculus was born at the same time as Newton’s and Leibniz’s classical calculus. In l’Hopital’s (1695) letter to Leibniz, when the viability of a derivative of order 1/2 was questioned, fractional calculus was first introduced, where Leibniz anticipated the development of fractional calculus     . After Leibniz, Euler was the second person to recognize the issue with non-integer orders. After Leibniz, Euler (1738) was the first to recognize the problem of non-integer orders  . Fourier (1822) proposed an integral representation for the concept of a derivative as the first definition of a derivative of any positive order   . Abel’s fractional calculus application in (1826) dealt with an equal time problem involving the solution of an integral equation   . Liouville’s initial definition, also known as Liouville’s first definition, was derived from the exponential function (1832). Liouville II, please. This is now known as the Liouville variant for the integration of noninteger order in terms of an integral. Ten years after Liouville’s passing, Riemann published the most important paper in response to numerous of his writings  . We can see that the Liouville and Riemann formulations both contain the complementary function. The Liouville and Riemann method that integration introduced must be used to fix this issue.

Grünwald  and Letnikov  independently devised the method for analyzing the derivatives of noninteger orders in terms of a straightforward convergent series. Under a helpful explanation of the alleged difference of non-integer orders, Letnikov has shown coherence between his definition and those proposed by Riemann and Liouville for certain order values. The derivative of non-integer orders of an analytical function must be expressed in terms of a Taylor series, claims a work by Hadamard (1892)  . From (1900), fractional calculus progressed quickly, and many definitions were produced in an effort to formulate specific problems, some of which we present. To solve an issue involving a certain class of functions, the periodic functions, Weyl  creates a derivative. Riesz develops the Fourier transform formula and the Mean Value Theorem for Fractional Integrals   . Liouville’s theory of (sufficiently excellent) functions is compatible with Marchaud’s (1927) definition of the order of non-integer derivatives   . Erde lyi-Kober (1940)   provided a separate definition for non-integer orders. Caputo (1967)  offers a definition that is more exact than Liouville and Riemann’s, but it is better suited for discussions of problems involving fractional differential equations with initial conditions  -  . The formula developed by Liouville and Riemann will be compared to this method. Because of the significance of his version with the derivative of non-integer orders from Liouville and Riemann, Caputo’s formulation takes into account the order of the integral and derivative operators. The two formulas will be contrasted. In Caputo, the integral of non-integer orders is calculated after the derivative of non-integer orders has been computed. In the Liouville and Riemann equation, the integral of non-integer orders is calculated first, followed by the derivative of integer orders. It is important to emphasize that problems can be solved with the Caputo derivative when the initial conditions of the function are met and each of them has an integer derivative. From the first meeting at the University of New Haven in, fractional calculus has advanced (1974), and as a consequence, many applications in many different scientific domains have emerged. There are several methods for solving problems involving derivatives.

In a variety of practical modeling issues, fractional differential equations are starting to find extensive application  . In fluid mechanics, the Time Fractional Burger equation is a type of subdiffusion convection equation. They are used to represent a variety of phenomena in the study of turbulent flow, including the propagation of shallow water waves and nonlinear acoustic waves in gas pipelines   , shock propagation, electromagnetic waves, turbulence, porous medium flows, pollutant flow, and temperature and pressure waves, as well as phenomena in the medical and scientific fields. These models, among others, help with better explanation and comprehension     Another example is the propagation of shock waves through viscous material  . Researchers frequently use this equation as a test case to determine the effectiveness of cutting-edge numerical techniques. By replacing the first-order time derivative with a fractional derivative. Where the Time Fractional Burger equation has been solved by a number of numerical methods  . The conventional Burger equation can be used to get this equation.

A fundamental partial differential equation is the Pittman-Burger equation or the standard Burger equation. Harry Pittman first presented the equation in (1915), and Johannes Martinus Berger investigated it in order to solve nonlinear equation systems in (1948). The exact solution of fractional differential equations can occasionally be difficult. Therefore, the goal of this study is to use the Mabel 18 program’s Variational Iteration Method to solve the Time Fractional Burger equation.

The numerical illustrations and error estimate provided by the Mabel program are discussed, and the Time Fractional Berger equation is solved.

2. Comparing Variational Iteration Method with Common Numerical Methods for Solving Time Fractional Burger Equation

The Variational Iteration Method (VIM) has been widely used to solve partial differential equations, including the time-fractional Burgers equation. In this section, we compare our results with existing research and discuss the advantages and limitations of the proposed method.

Firstly, we compare our results with those obtained by the Adomian decomposition method (ADM) and the homotopy perturbation method (HPM). The ADM is another powerful analytical method that has been used to solve fractional differential equations. The HPM is a modified version of the ADM that is based on the concept of homotopy. Both methods have been applied to solve the time-fractional Burgers equation, and their results have been reported in the literature.

Our numerical results show that the VIM provides a more accurate and efficient solution than the ADM and HPM. This is because the VIM does not require the calculation of complicated Adomian polynomials or the construction of a homotopy equation, which can be time-consuming and error-prone. Instead, the VIM uses a simple and intuitive iterative process that can converge quickly to the exact solution.

Secondly, we compare our results with those obtained by the finite difference method (FDM), which is a popular numerical method for solving partial differential equations. The FDM discretizes the spatial and temporal domains of the equation and approximates the derivatives using finite differences. The resulting system of algebraic equations is then solved using matrix inversion or iterative methods.

Our numerical results show that the VIM provides a comparable solution to the FDM, but with a much lower computational cost. This is because the VIM does not require the discretization of the domain or the solution of a large system of equations, which can be computationally expensive. Instead, the VIM uses a series expansion and a correction functional to obtain an analytical solution that can be easily evaluated.

In terms of practical applications, the time-fractional Burgers equation is a fundamental model in fluid dynamics, combustion, and nonlinear acoustics. The equation describes the behavior of a fluid or gas in which the velocity field is affected by both viscosity and diffusion. The equation is also used to model the propagation of sound waves in a nonlinear medium.

By solving the time-fractional Burgers equation using the VIM, we can study the behavior of these systems under different conditions, such as the effect of viscosity and diffusion on the flow, or the nonlinear behavior of sound waves. This can have important practical applications in the design of fluid flow systems, combustion engines, and acoustic devices.

In conclusion, the VIM is a powerful analytical technique for solving partial differential equations, including the time-fractional Burgers equation. Our numerical results show that the VIM provides a more accurate and efficient solution than existing methods such as the ADM, HPM, and FDM. The proposed method can be applied to practical problems in fluid mechanics, combustion, and acoustics, and can provide valuable insights into the behavior of these systems under different conditions.

3. Definitions

Definition 3.1. The left sided Riemannian–Liouville fractional integral of order $\mu \ge 0$ , of a function $f\in {C}_{\alpha },\alpha \ge -1$ , is defined as

${I}^{\mu }f\left(x\right)=\frac{1}{\text{Γ}\left(\mu \right)}{\int }_{0}^{x}\frac{f\left(t\right)}{{\left(x-t\right)}^{1-\mu }}\text{d}t,\mu >0,x>0,{I}^{0}f\left(t\right)=f\left(x\right)$ (1)

Definition 3.2. Let $f\in {C}_{-1}^{m},m\in N$ . Then the Caputo fractional derivative of f is defined as   :

${D}^{\mu }f\left(x\right)=\left\{\begin{array}{l}{I}^{m-\mu }{f}^{\left(m\right)}\left(x\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}m-1<\mu \le m\\ \frac{{\text{d}}^{m}}{\text{d}{t}^{m}}f\left(t\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mu =m\end{array}$ (2)

${I}^{\mu }{I}^{v}f={I}^{\mu +v}f,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mu ,v\ge 0,f\in {C}_{\alpha },\alpha \ge -1,$

${I}^{\mu }{x}^{y}=\frac{\Gamma \left(\gamma +1\right)}{\Gamma \left(\gamma +\mu +1\right)}{x}^{\gamma +\mu },\mu >0,\gamma >-1,x>0$ (3)

Lemma 3.1. If $m-1<\alpha \le m$ , and $f\in {L}_{1}\left[a,b\right]$ , then

${J}_{a}^{\alpha }f\left(x\right)=\frac{1}{\Gamma \left(\alpha \right)}{\int }_{a}^{x}{\left(x-t\right)}^{\alpha -1}f\left(t\right)\text{d}t,\text{\hspace{0.17em}}\text{\hspace{0.17em}}{D}_{a}^{\alpha }{J}_{a}^{\alpha }f\left(x\right)=f\left(x\right)$ (4)

${J}_{a}^{\alpha }{D}_{a}^{\alpha }f\left(x\right)=f\left(x\right)-\underset{k=0}{\overset{m-1}{\sum }}\text{ }\text{ }{f}^{\left(k\right)}\left({0}^{+}\right)\frac{{\left(x-a\right)}^{k}}{k!},x>0$ (5)

Definition 3.3. The fractional derivative of f(x) in the Caputo sense is defined as

${D}^{\alpha }f\left(x\right)={J}^{m-\alpha }{D}^{m}f\left(x\right)=\frac{1}{\text{Γ}\left(m-\alpha \right)}{\int }_{0}^{x}{\left(x-t\right)}^{m-\alpha -1}{f}^{\left(m\right)}\left(t\right)\text{d}t,$ (6)

For $m-1<\alpha \le m,m\in ℕ,x>0$

Definition 3.4. For m to be the smallest integer that exceeds $\alpha$ , the Caputo time-fractional derivative operator of order $\alpha >0$ is defined as

$\begin{array}{l}{D}_{t}^{\alpha }u\left(x,t\right)=\frac{{\partial }^{\alpha }u\left(x,t\right)}{\partial {t}^{\alpha }}\\ =\left\{\begin{array}{l}\frac{1}{\text{Γ}\left(m-\alpha \right)}{\int }_{0}^{t}{\left(t-\tau \right)}^{m-\alpha -1}\frac{{\partial }^{m}u\left(x,\tau \right)}{\partial {t}^{m}}\text{d}\tau ,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}m-1<\alpha (7)

4. Variational Iteration Method (VIM)

The need of finding precise or approximate solutions continues to be a key area in mathematics, and this is done by searching for creative methods to do so. The Time Fractional Burger equation and related issues like electromagnetic control, dynamics of generalized Burgers’ nanoliquid flow containing motile microorganisms, Galerkin finite element mechanism  and viscous dissipation, Joule heating effects in non-Fourier MHD squeezing flow, and heat and mass transfer between rigid plates with thermal radiation can all be solved numerically using the variational iteration method (VIM). The governing equations in a given problem are approximated using a direct, iterative approach by the (VIM).

For numerical implementations, Maple’s Variational iteration package is required, which is specifically created for this method. The Variational Iteration Method (VIM) is a technique that helps in constructing accelerated approximations to get to the precise solution without the requirement for specific constraints for both linear and nonlinear, homogeneous and inhomogeneous equations. This approach delivers the answer in a sequential manner that leads to the precise answer    The Variational Iteration method (VIM)  , an enhancement to the general Lagrange multiplier method  , was devised by Chinese mathematician Ji-Huan He     . (VIM) demonstrated excellent performance in the solution of fractional differential equations  .

The (VIM) was subsequently used to solve more fractional differential equations, proving its effectiveness and accuracy    .

We examine the following differential equation to demonstrate the fundamental ideas of the variational iteration method:

$Lu+Nu=g\left(t\right)$ (8)

where L and N are linear and nonlinear operators, respectively, and $g\left(t\right)$ is an inhomogeneous term.

The form of the Variational iteration Technique is

${u}_{n+1}\left(x\right)={u}_{n}\left(x\right)+{\int }_{0}^{x}\lambda \left\{\left(L{u}_{n}\left(\tau \right)+N{\stackrel{˜}{u}}_{n}\left(\tau \right)-g\left(\tau \right)\right)\right\}\text{d}\tau$ (9)

where λ is a general Lagrange’s multiplier  , λ may be a constant or a function, and ${\stackrel{˜}{u}}_{n}$ is a restricted value that means it behaves as a constant  , hence $\delta {\stackrel{˜}{u}}_{n}=0$ , where δ is the variational derivative. The Lagrange multiplier λ can be identified by

$\lambda ={\left(-1\right)}^{n}\frac{1}{\left(n-1\right)!}{\left(\xi -x\right)}^{n-1}$ (10)

where n is the number of recurrences of differentials.

The Lagrange multiplier can be precisely identified, allowing for only one iteration step to produce the exact answer. The nonlinear terms in nonlinear problems must be viewed as constrained variations in order to compute the Lagrange multiplier in a straightforward manner. Hence, the precise answer can be found by using

$u\left(x,t\right)=\underset{n\to \infty }{\mathrm{lim}}{u}_{n}\left(x,t\right)$ (11)

The time-fractional Burgers equation

The following initial value problem is applied to the one-dimensional time-fractional Burgers equation  .

$\frac{{\partial }^{\alpha }u}{\partial {t}^{\alpha }}+\epsilon u\frac{\partial u}{\partial x}=v\frac{{\partial }^{2}u}{\partial {x}^{2}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}t>0,0<\alpha \le 1$ (12)

$u\left(x,0\right)=g\left( x \right)$

We create a correction functional that reads as follows to solve Equation (12) using the variational iteration method.

$\begin{array}{c}{u}_{n+1}\left(x,t\right)={u}_{n}\left(x,t\right)+{\int }_{0}^{t}\lambda \left\{{\left(\frac{\partial {u}_{n}\left(x,\tau \right)}{\partial \tau }\right)}_{\alpha }\\ \text{\hspace{0.17em}}\text{ }+\left({\stackrel{˜}{u}}_{n}\left(x,\tau \right)\right){\left(\frac{\partial {\stackrel{˜}{u}}_{n}\left(x,\tau \right)}{\partial \tau }\right)}_{x}-v{\left(\frac{\partial {\stackrel{˜}{u}}_{n}\left(x,\tau \right)}{\partial \tau }\right)}_{xx}\right\}\text{d}\tau \end{array}$ (13)

where $\delta {\stackrel{˜}{u}}_{n}$ is considered as a restricted variation. ${u}_{0}\left(x,t\right)$ is its initial approximation or trial function. Making the above correction functional stationary and noticing that $\delta {\stackrel{˜}{u}}_{n}=0$ , we get

$\begin{array}{c}\delta {u}_{n+1}\left(x,t\right)=\delta {u}_{n}\left(x,t\right)+{\int }_{0}^{t}\delta \lambda \left(\tau \right)\left\{{\left(\frac{\partial {u}_{n}\left(x,\tau \right)}{\partial \tau }\right)}_{\alpha }\\ \text{\hspace{0.17em}}\text{ }+\left({\stackrel{˜}{u}}_{n}\left(x,\tau \right)\right){\left(\frac{\partial {\stackrel{˜}{u}}_{n}\left(x,\tau \right)}{\partial \tau }\right)}_{x}-v{\left(\frac{\partial {\stackrel{˜}{u}}_{n}\left(x,\tau \right)}{\partial \tau }\right)}_{xx}\right\}\text{d}\tau \end{array}$

$\begin{array}{c}\delta {u}_{n+1}\left(x,t\right)=\delta {u}_{n}\left(x,t\right)+\lambda \left(\tau \right)\delta {u}_{n}\left(x,\tau \right)+{\int }_{0}^{t}\delta {u}_{n}\left(x,\tau \right){\lambda }^{\prime }\left(\tau \right)\text{d}\tau \\ =\delta {u}_{n}\left(x,t\right)\left(1+\lambda \left(\tau \right)\right)+{\int }_{0}^{t}\delta {u}_{n}\left(x,\tau \right){\lambda }^{\prime }\left(\tau \right)\text{d}\tau =0\end{array}$

which produces the stationary conditions:

${\lambda }^{\prime }\left(\tau \right)=0$ (14a)

$1+{\lambda \left(\tau \right)|}_{\tau =t}=0$ (14b)

when the Lagrange-Euler equation (Equation (14a)) and the natural boundary condition (Equation (14b)) are concerned. The following variational iteration formula can be constructed by identifying the Lagrange multiplier as $\lambda =1$

${u}_{n+1}\left(x,t\right)={u}_{n}\left(x,t\right)-{\int }_{0}^{x}\left\{{\left({u}_{n}\right)}_{\alpha }+\left({\stackrel{˜}{u}}_{n}\right){\left({u}_{n}\right)}_{x}-v{\left({u}_{n}\right)}_{xx}\right\}\text{d}\tau$ (15)

Starting with an initial approximation ${u}_{0}=u\left(x,0\right)$ provided by Equation (12), we can extract the other components directly using the aforementioned iteration formula (15).

${u}_{0}\left(x,t\right)=g\left(x\right)=\frac{\mu +\sigma +\left(\sigma -\mu \right)\mathrm{exp}\left(\gamma \right)}{1+\mathrm{exp}\left(\gamma \right)}$ (16)

${u}_{1}\left(x,t\right)=\left[v{g}^{″}-\epsilon g{g}^{\prime }\right]\frac{{t}^{\alpha }}{\text{Γ}\left(\alpha +1\right)}$ (17)

${u}_{2}\left(x,t\right)=\left[2{\epsilon }^{2}g{{g}^{\prime }}^{2}+{\epsilon }^{2}{g}^{2}{g}^{″}-4\epsilon v{g}^{\prime }{g}^{″}-2\epsilon vg{g}^{″}+{v}^{2}{g}^{\left(4\right)}\right]\frac{{t}^{2\alpha }}{\text{Γ}\left(2\alpha +1\right)}$ (18)

As a result, we have the answer to Equation (11) in series form.

$\begin{array}{c}u\left(x,t\right)=\frac{\mu +\sigma +\left(\sigma -\mu \right)\mathrm{exp}\left(\gamma \right)}{1+\mathrm{exp}\left(\gamma \right)}+\left[v{g}^{″}-\epsilon g{g}^{\prime }\right]\frac{{t}^{\alpha }}{\text{Γ}\left(\alpha +1\right)}\\ \text{\hspace{0.17em}}\text{ }+\left[2{\epsilon }^{2}g{{g}^{\prime }}^{2}+{\epsilon }^{2}{g}^{2}{g}^{″}-4\epsilon v{g}^{\prime }{g}^{″}-2\epsilon vg{g}^{″}+{v}^{2}{g}^{\left(4\right)}\right]\frac{{t}^{2\alpha }}{\text{Γ}\left(2\alpha +1\right)}+\cdots \end{array}$ (19)

5. Application to Obtain the Numerical Solution of Time Fractional Burgers’ Equation

5.1. Example 1

We consider one-dimensional Time Fractional Burgers Equation

$\frac{{\partial }^{\alpha }u}{\partial {t}^{\alpha }}-\frac{{\partial }^{2}u}{\partial {x}^{2}}+u\frac{\partial u}{\partial x}=0$ (20)

With initial conditions $u\left(x,0\right)=x$

$0<\alpha \le 1,\text{\hspace{0.17em}}\text{\hspace{0.17em}}t>0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}x\in R$

$\alpha =1$

5.2. Example 2

We consider one-dimensional Time Fractional Burgers Equation

$\frac{{\partial }^{\alpha }u}{\partial {t}^{\alpha }}-\frac{{\partial }^{2}u}{\partial {x}^{2}}+u\frac{\partial u}{\partial x}=0$ (21)

With initial conditions $u\left(x,0\right)=x$

$0<\alpha \le 1,\text{\hspace{0.17em}}\text{\hspace{0.17em}}t>0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}x\in R$

$\alpha =0.75$

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 and approximate solutions result of time-fractional Burgers equation in Example 1.

Figure 2. Graph showing the correspondence between exact and approximate solutions result of time-fractional Burgers equation in Example 2.

Table 1. Numerical results and exact solution of one-dimensional time fractional burgers equation for Example 1.

Table 2. Numerical results and Exact solution of one-dimensional Time Fractional Burgers equation for Example 2.

6. Conclusion

The time-fractional Burgers equation is resolved in Maple18 using the Variational Iteration method. By contrasting the numerical outcomes, the outcomes were compared with the precise solution corresponding to the time-fractional Burgers equation. This showed the procedure’s efficacy and Maple18’s capacity to swiftly and efficiently provide a numerical solution that was related to the exact solution while recording the error value, making the correctness of the solutions obtained extremely satisfactorily. We can observe that the exact solution and the numerical solution are typically related. With Maple18, it is possible to numerically calculate the majority of engineering and mathematics topics. Mabel18 is both a mathematical system and a programming language. The solution has also been graphically depicted. These results are shown in Table 1, Table 2, Figure 1 and Figure 2 using the package version of Mabel. The distinction between exact solution and numerical solutions is seen in Table 1, Table 2, Figure 1 and Figure 2. We were able to reach quite near to the exact solutions of the equations using the Variational Iteration method via Mabel software. The results show how successful the existing methodology is at obtaining precise numerical solutions to time-fractional Burgers equation. The main objective of this work is to use Maple software to automate the calculation of the Variational Iteration method. This will make it easier to use Mabel in the future as we may be able to gain rough estimates of the solutions.

7. Recommendations for Future Research

The research on the Variational Iteration Method for solving time fractional Burgers equation using Maple has shown promising results. The method was found to be effective in providing accurate solutions to the equation with a high degree of efficiency. The Maple software proved to be a useful tool in implementing the method and producing the numerical results.

In conclusion, the Variational Iteration Method is a powerful technique for solving time fractional Burgers equation, and Maple can be used effectively to implement the method. The study has provided valuable insights into the behavior of the solution and the effectiveness of the method. However, there is still room for improvement in the accuracy and efficiency of the method, especially for more complex problems. Future research can focus on developing more efficient algorithms for the method and exploring its applications in other fields of science and engineering.

In addition, it may be worthwhile to investigate the use of other numerical methods in solving time fractional Burgers equation and compare their performance with the Variational Iteration Method. Furthermore, the study could be extended to investigate the behavior of the solution for different types of initial and boundary conditions. Overall, the research has opened up new avenues for further exploration in the field of fractional calculus and numerical methods, and there is great potential for future advancements.

Acknowledgements

This project was funded by the Deanship of Scientific Research (DSR), King Abdulaziz University. The authors, therefore, acknowledge with thanks DSR technical and financial support.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.     customer@scirp.org +86 18163351462(WhatsApp) 1655362766  Paper Publishing WeChat 