Journal of Function Spaces
Volume 2021 (2021), Article ID 5588601, 13 pages
https://doi.org/10.1155/2021/5588601

Research Article

The New Semianalytical Technique for the Solution of Fractional-Order Navier-Stokes Equation

Nehad Ali Shah,¹ Mounirah Areshi,² Jae Dong Chung,¹ and Kamsing Nonlaopon³

¹Department of Mechanical Engineering, Sejong University, Seoul 05006, Republic of Korea
²Mathematics Department, College of Science, University of Tabuk, Tabuk, Saudi Arabia
³Department of Mathematics, Faculty of Science, Khon Kaen University, Khon Kaen 40002, Thailand

Correspondence should be addressed to Kamsing Nonlaopon; nkamsi@kku.ac.th

Received 1 March 2021; Revised 29 March 2021; Accepted 3 April 2021; Published 20 April 2021

Academic Editor: Qifeng Zhang

Copyright © 2021 Nehad Ali Shah 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.

Abstract.

In this paper, we introduce a modified method which is constructed by mixing the residual power series method and the Elzaki transformation. Precisely, we provide the details of implementing the suggested technique to investigate the fractional-order nonlinear models. Second, we test the efficiency and the validity of the technique on the fractional-order Navier-Stokes models. Then, we apply this new method to analyze the fractional-order nonlinear system of Navier-Stokes models. Finally, we provide 3-D graphical plots to support the impact of the fractional derivative acting on the behavior of the obtained profile solutions to the suggested models.

1. Introduction

The fractional-order Navier-Stokes equation (NSE) has been extensively analyzed. These equations model the fluid motion defined by several physical processes, such as the movement of blood, the ocean’s current, the flow of liquid in vessels, and the airflow around an aircraft’s arms [1–3]. The classical NSEs were generalized by El-Shahed and Salem [4] by replacing the first time derivative with a Caputo fractional derivative of order , where . Using Hankel transform, Fourier sine transform, and Laplace transform, the researchers achieved the exact solution for three different equations. In 2006, Momani and Odibat [5] solve fractional-order NSEs using the Adomian decomposition method. Ganji et al. [6] applied an analytical technique, the homotopy perturbation method, for solving the fractional-order NSEs in polar coordinates, and the results achieved were expressed in a closed form. Singh and Kumar [7] solved the fractional-order reduced differential transformation method (FRDTM) to achieve an approximated analytical result of fractional-order multidimensional NSE. Oliveira and Oliveira [8] analyzed the residual power series method (RPSM) to find the result of the nonlinear fractional-order two-dimensional NSEs. Zhang and Wang [9] suggested numerical analysis for a class of NSEs with fractional-order derivatives; Ravindran, the exact boundary controllability of Galerkin approximations of a Navier-Stokes system for soret convection [10]; and Cibik and Yilmaz, the Brezzi-Pitkaranta stabilization and a priori error analysis for the Stokes control [11].

Some researchers mix two powerful techniques to achieve another result technique to solve equations and systems of fractional-order NSEs. Below, we define some of these combinations: a combination of the Laplace transformation and Adomian decomposition method; Kumar et al. [12] introduced the homotopy perturbation transform method (HPTM), combined Laplace transformation with the homotopy perturbation method, and solved fractional-order NSEs in a tube. Jena and Chakraverty [13] implemented the homotopy perturbation transformation method (HPETM), and this technique consists in the mixture of Elzaki transformation technique and homotopy perturbation technique; Prakash et. al [1] suggested -homotopy analysis transformation technique to achieve a result of coupled fractional-order NSEs. This technique mixture of the Laplace transformation and residual power series method is defined:

where is the Caputo derivative of order , is the velocity vector, is the time, is kinematics viscosity, is the pressure, and is the density.

In this work, we consider two special cases. First, we consider unsteady, one-dimensional motion of a viscous fluid in a tube. The fractional-order Navier-Stokes equations in cylindrical coordinates that governs the flow field in the tube are given by

with initial condition

where and is a function depending only on .

Consider that the fractional-order two-dimensional Navier-Stokes equations is defined as

with initial conditions

where denote as constant density, time, and pressure, respectively. are the spatial components, and and are two functions depending only on and .

The residual power series method (RPSM) is a simple and efficient technique for constructing a power series result for extremely linear and nonlinear equations without perturbation, linearization, and discretization. Unlike the classical power series technique, the RPS approach does not need to compare the coefficients of the corresponding terms and a recursion relation is not required. This approach calculates the power series coefficients by a series of algebraic equations of one or more variables, and its reliance on derivation, which is much simpler and more precise than integration, which is the basis of most other solution approaches, is the main advantage of this methodology. This method is, in effect, an alternative strategy for obtaining theoretical results for the fractional-order partial differential equations [14].

The RPSM was introduced as an essential tool for assessing the power series solution’s values for the first and second-order fuzzy DEs [15]. It has been successfully implemented in the approximate result of the generalized Lane-Emden equation [16], which is a highly nonlinear singular DE, in the inaccurate work of higher-order regular DEs [17], in the solution of composite and noncomposite fractional-order DEs [18], in predicting and showing the diversity of results to the fractional-order boundary value equations [19], and in the numerical development of the nonlinear fractional-order KdV and Burgers equation [20], in addition to some other implementations [21–23], and recently, it has been applied to investigate the approximate result of a fractional-order two-component evolutionary scheme [24].

This paper introduces the modified analytical technique: the residual power series transform method (RPSTM) is implemented to investigate the fractional-order NS equations. The result of certain illustrative cases is discussed to explain the feasibility of the suggested method. The results of fractional-order models and integral-order models are defined by using the current techniques. The new approach has lower computing costs and higher rate convergence. The suggested method is also constructive for addressing other fractional orders of linear and nonlinear PDEs.

2. Preliminaries

Definition 1. The Abel-Riemann of fractional operator of order is given as [25–27]

where and

Definition 2. The fractional-order Abel-Riemann integration operator is defined as [25–27]

The operator of basic properties

Definition 3. The Caputo fractional operator of is defined as [25–27]

Definition 4. The fractional-order Caputo operator of Elzaki transform is given as [25–27]

Definition 5. A power series definition of the form [14]

where and is called fractional power series (FPS) about , where are the constants called the coefficients of the series. If , then the fractional power series will be reduced to the fractional Maclaurin series.

Theorem 6. Assume that has a fractional power series representation at of the form [14]

For , if are continuous on , then the coefficients can be written as

where , and is the radius of convergence.

Definition 7. The expansion of power series of the form [14]

is said to be multifractional power series at , where are the coefficients of multifractional power series.

Theorem 8. Let us assume that which has the multifractional power representation at can be written as [14]

For , if are the continuous on , then the coefficients are given by

where and .
So, we can write the fractional power expansion of of the form

which is the generalized Taylor expansion. If we consider , then the generalized Taylor formula will be converted to classical Taylor series.

Corollary 9. Let us assume that has a multifractional power series representation about as [14]

For if are continuous on , then

where .

3. The Procedure of RPSTM

In this section, we explain the steps of RPSTM for solving the fractional-order partial differential equation

with initial condition

First, we use the Elzaki transform to (21); we get

By the fact that and using the initial condition (22), we can write (23) as

where

Second, we define the transform function as the following formula:

We write the th truncated series of (25) as

As stated in [25], the definition of Elzaki residual function to (25) is

and the th Elzaki residual function of (27) is

Third, we expand a few of the properties arising in the basic RPSM to find out certain facts:

(i) and for each

(ii), ,

Furthermore, to evaluate the coefficient functions we can recursively solve the following scheme

Finally, we implemented the Elzaki inverse to to achieve the th approximate supportive solution .

4. Numerical Results

Example 1. Consider the time-fractional-order one-dimensional NS equation of the form

Subject to the initial condition

Applying Elzaki transform to (31) and using the initial condition given in (32), we get

The th truncated term series of (33) is

and the th Elzaki residual function is

Now, to determine , we substitute the th truncated series (34) into the th Elzaki residual function (35), multiply the resulting equation by , and then solve recursively the relation for . The following are the first several components of the series :

Putting the values of in (34), we get

Using inverse Elzaki transform to (38), we get

Putting , we have

In Figure 1, the RPSTM and the exact results of Example 1 at are shown by plots (a) and (b), respectively. From the given figures, it can be seen that both the exact and the EDM results are in close contact with each other. Also, in the Figure 2 subgraph, the RPSTM results of Example 1 are calculated at different fractional-order and . It is investigated that fractional-order problem results are convergent to an integer-order effect as fractional-order analysis to integer-order. The same phenomenon of convergence of fractional-order solutions towards integral-order solutions is observed.

Example 2. Consider the fractional-order one-dimensional NS equation of the form

Figure 1: Graph of exact and analytical results of Problem 1.

Figure 2: The fractional order of and of Problem 1.

Subject to the initial condition,

Applying Elzaki transform to (41) and using the initial condition given in (42), we get

The th truncated term series of (43) is

and the th Elzaki residual function is

Now, to determine we substitute the th truncated series (44) into the th Elzaki residual function (45), multiply the resulting equation by , and then solve recursively the relation , for . The following are the first several components of the series :

Putting the values of in (44), we get

Using inverse Elzaki transform to (48), we get

Putting , we have

In Figure 3, the RPSTM and the exact results of Example 2 at are shown by graphs, respectively. From the given figures, it can be seen that both the exact and the EDM results are in close contact with each other. Also, in the Figure 4 subgraph, the RPSTM results of Example 2 are calculated at different fractional-order and . It is investigated that fractional-order problem results are convergent to an integer-order effect as fractional-order analysis to integer-order. The same phenomenon of convergence of fractional-order solutions towards integral-order solutions is observed.

Example 3. Consider the fractional-order two-dimensional NS equation of the form

with initial condition

Figure 3: Graph of exact and analytical results of Problem 2.

Figure 4: The fractional order of and of Problem 2.

Applying Elzaki transform to (51) and using (52), we get

The th truncated term series of (53) is

and the th Elzaki residual function is

Now, to determine and we substitute the th truncated series (54) into the th Elzaki residual function (55), multiply the resulting equation by , and then solve recursively the relation for and . The following are the first several components of the series and :

Putting the values of and in (54), we have

Using inverse Elzaki transform, we get

Putting , we get the solution in closed form

In Figures 5 and 6, the RPSTM and the exact results of Example 3 at are shown by graphs, respectively. From the given figures, it can be seen that both the exact and the RPSTM results are in close contact with each other. Also, in the Figure 7 and 8 subgraph, the RPSTM results of Example 3 are calculated at different fractional-order and . It is investigated that fractional-order problem results are convergent to an integer-order effect as fractional-order analysis to integer-order. The same phenomenon of convergence of fractional-order solutions towards integral-order solutions is observed.

Example 4. Consider the fractional-order two-dimensional NS equation as

with initial condition

Figure 5: Graph of exact and analytical results of Problem 3.

Figure 6: Graph of exact and analytical results of Problem 3.

Figure 7: The fractional order of and of Problem 3.

Figure 8: The fractional order of and of Problem 3.

Applying Elzaki transform to (60) and using (61), we get

The th truncated term series of (62) is

and the th Elzaki residual function is

Now, to determine and , , we substitute the th truncated series (63) into the th Elzaki residual function (64), multiply the resulting equation by , and then solve recursively the relation , for and . The following are the first several components of the series and :

Putting the values of and in (63), we have

Using inverse Elzaki transform, we get

Putting , we get the solution in closed form

In Figures 9 and 10, the RPSTM and the exact results of Example 4 at are shown by graphs, respectively. From the given figures, it can be seen that both the exact and the RPSTM results are in close contact with each other. Also, in the Figure 11 and 12 subgraph, the RPSTM results of Example 4 are calculated at different fractional-order and . It is investigated that fractional-order problem results are convergent to an integer-order effect as fractional-order analysis to integer-order. The same phenomenon of convergence of fractional-order solutions towards integral-order solutions is observed.

Figure 9: Graph of exact and analytical results of Problem 4.

Figure 10: Graph of exact and analytical results of Problem 4.

Figure 11: The fractional order of and of Problem 4.

Figure 12: The fractional order of and of Problem 4.

5. Conclusions

In this article, a modified method constructed by a mixture of the residual power series and Elzaki transformation operator is presented to solve fractional-order Navier-Stokes models. The merit of the modified technique is to reduce the size of computational work needed to find the result in a power series form whose coefficient to be calculated is in successive algebraic steps. The technique gives a series form of results that converges very fast in physical models. It is predicted that this article achieved results which will be useful for further analysis of the complicated nonlinear physical problems. The calculations of this technique are very straightforward and simple. Thus, we deduce that this technique can be implemented to solve several schemes of nonlinear fractional-order partial differential equations.

Data Availability

The numerical data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this article.

Acknowledgments

This work was supported by the Korea Institute of Energy Technology Evaluation and Planning (KETEP) grant funded by the Korea government (MOTIE) (20202020900060, The Development and Application of Operational Technology in Smart Farm Utilizing Waste Heat from Particulates Reduced Smokestack).

References

A. Prakash, P. Veeresha, D. G. Prakasha, and M. Goyal, “A new efficient technique for solving fractional coupled Navier-Stokes equations using -homotopy analysis transform method,” Pramana, vol. 93, no. 1, pp. 1–10, 2019.

A. Rauf, Y. Mahsud, I. A. Mirza, and Q. Rubbab, “Multi-layer flows of immiscible fractional Maxwell fluids with generalized thermal flux,” Chinese Journal of Physics, vol. 62, pp. 313–334, 2019.

A. Shah, N. A. Nehad, T. Elnaqeeb, and M. M. Rashidi, “Magnetohydrodynamic free convection flows with thermal memory over a moving vertical plate in porous medium,” Journal of Applied and Computational Mechanics, vol. 5, no. 1, pp. 150–161, 2019.

M. El-Shahed and A. Salem, “On the generalized Navier-Stokes equations,” Applied Mathematics and Computation, vol. 156, no. 1, pp. 287–293, 2004.

S. Momani and Z. Odibat, “Analytical solution of a time-fractional Navier-Stokes equation by Adomian decomposition method,” Applied Mathematics and Computation, vol. 177, no. 2, pp. 488–494, 2006.

Z. Z. Ganji, D. D. Ganji, A. D. Ganji, and M. Rostamian, “Analytical solution of time-fractional Navier-Stokes equation in polar coordinate by homotopy perturbation method,” Numerical Methods for Partial Differential Equations, vol. 26, no. 1, pp. 117–124, 2010.

B. K. Singh and P. Kumar, “FRDTM for numerical simulation of multi-dimensional, time-fractional model of Navier-Stokes equation,” Ain Shams Engineering Journal, vol. 9, no. 4, pp. 827–834, 2018.

D. S. Oliveira and E. C. de Oliveira, “Analytical solutions for Navier-Stokes equations with Caputo fractional derivative,” SeMA Journal, vol. 78, no. 1, pp. 137–154, 2021.

J. Zhang and J. R. Wang, “Numerical analysis for Navier-Stokes equations with time fractional derivatives,” Applied Mathematics and Computation, vol. 336, pp. 481–489, 2018.

S. S. Ravindran, “Exact boundary controllability of Galerkin approximations of Navier-Stokes system for soret convection,” An International Journal of Optimization and Control: Theories & Applications (IJOCTA), vol. 5, no. 2, pp. 41–49, 2015.

A. Cibik and F. Yilmaz, “Brezzi-Pitkaranta stabilization and a priori error analysis for the Stokes control,” An International Journal of Optimization and Control: Theories & Applications (IJOCTA), vol. 7, no. 1, pp. 75–82, 2017.

D. Kumar, J. Singh, and S. Kumar, “A fractional model of Navier-Stokes equation arising in unsteady flow of a viscous fluid,” Journal of the Association of Arab Universities for Basic and Applied Sciences, vol. 17, no. 1, pp. 14–19, 2015.

R. M. Jena and S. Chakraverty, “Solving time-fractional Navier-Stokes equations using homotopy perturbation Elzaki transform,” SN Applied Sciences, vol. 1, no. 1, pp. 1–13, 2019.

M. Alquran, M. Ali, M. Alsukhour, and I. Jaradat, “Promoted residual power series technique with Laplace transform to solve some time-fractional problems arising in physics,” Results in Physics, vol. 19, p. 103667, 2020.

A. Çibik, “The effect of a sparse grad-div stabilization on control of stationary Navier-Stokes equations,” Journal of Mathematical Analysis and Applications, vol. 437, no. 1, pp. 613–628, 2016.

O. A. Arqub, A. El-Ajou, A. S. Bataineh, and I. Hashim, “A representation of the exact solution of generalized Lane-Emden equations using a new analytical method,” Abstract and Applied Analysis, vol. 2013, Article ID 378593, 10 pages, 2013.

A. Arqub, Z. A.-H. Omar, R. Al-Badarneh, and S. Momani, “A reliable analytical method for solving higher-order initial value problems,” Discrete Dynamics in Nature and Society, vol. 2013, Article ID 673829, 12 pages, 2013.

A. El-Ajou, O. A. Arqub, Z. Al Zhour, and S. Momani, “New results on fractional power series: theories and applications,” Entropy, vol. 15, no. 12, pp. 5305–5323, 2013.

O. A. Arqub, A. El-Ajou, Z. Al Zhour, and S. Momani, “Multiple solutions of nonlinear boundary value problems of fractional order: a new analytic iterative technique,” Entropy, vol. 16, no. 1, pp. 471–493, 2014.

A. El-Ajou, O. A. Arqub, and S. Momani, “Approximate analytical solution of the nonlinear fractional KdV-Burgers equation: a new iterative algorithm,” Journal of Computational Physics, vol. 293, pp. 81–95, 2015.

O. A. Arqub, A. El-Ajou, and S. Momani, “Constructing and predicting solitary pattern solutions for nonlinear time-fractional dispersive partial differential equations,” Journal of Computational Physics, vol. 293, pp. 385–399, 2015.

M. Alquran, “Analytical solutions of fractional foam drainage equation by residual power series method,” Mathematical Sciences, vol. 8, no. 4, pp. 153–160, 2014.

S.-D. Lin and L. Chia-Hung, “Laplace transform for solving some families of fractional differential equations and its applications,” Advances in Difference Equations, vol. 2013, no. 1, p. 9, 2013.

M. Alquran, “Analytical solution of time-fractional two-component evolutionary system of order 2 by residual power series method,” Journal of Applied Analysis & Computation, vol. 5, no. 4, pp. 589–599, 2015.

T. M. Elzaki, “The new integral transform Elzaki transform,” Global Journal of Pure and Applied Mathematics, vol. 7, no. 1, pp. 57–64, 2011.

T. M. Elzaki, “On the connections between Laplace and Elzaki transforms,” Advances in Theoretical and Applied Mathematics, vol. 6, no. 1, pp. 1–11, 2011.

T. M. Elzaki, “On the new integral transform"Elzaki transform"fundamental properties investigations and applications,” Global Journal of Mathematical Sciences: Theory and Practical, vol. 4, no. 1, pp. 1–13, 2012.