Approximate Solutions of the Model Describing Fluid Flow Using Generalized ρ -Laplace Transform Method and Heat Balance Integral Method

: This paper addresses the solution of the incompressible second-grade ﬂuid models. Fundamental qualitative properties of the solution are primarily studied for proving the adequacy of the physical interpretations of the proposed model. We use the Liouville-Caputo fractional derivative with its generalized version that gives more comprehensive physical results in the analysis and investigations. In this work, both the ρ -Laplace homotopy transform method ( ρ -LHTM) and the heat balance integral method (HBIM) are successfully combined to solve the fractional incompressible second-grade ﬂuid differential equations. Numerical simulations and their physical interpretations of the mentioned incompressible second-grade ﬂuid model are ensured to illustrate the main ﬁndings. It is also proposed that one can recognize the differences in physical analysis of diffusions such as ballistic diffusion, super diffusion, and subdiffusion cases by considering the impact of the orders ρ and ϕ .


Introduction
Fractional calculus is a field trying to understand the real-world phenomena modeled with non-integer-order derivatives [1,2], finance [3,4], biological processes and epidemic models [5][6][7][8], science and engineering [9][10][11][12][13][14][15][16][17], mechanics, etc. There exist at present numerous fractional operators which are defined by virtue of singular and nonsingular kernels. Since fractional terminology began with Leibniz's question in 1695, the list of the existing fractional operators is naturally long. With singular kernels, we have the Caputo operator [18] and the Riemann-Liouville operator [18]. Without singular kernels, we have two types: the fractional derivative with exponential kernel [19,20] and the fractional operator which has Mittag-Leffler function as a kernel known as the Atangana-Baleanu operator [21][22][23]. We also summarized the existing FDs, which are the most popular in the literature. Some other types of derivatives as conformable [24], the Hilfer derivation [25], the Erdélyi-Kober derivation, and others. Note the discrete forms of the fractional derivatives previously cited exist in the literature; for more information, see the following investigations [26,27].
In the literature, there exist several investigations related to these fractional differential equations. We briefly enumerate some of them. Reference [28] proposed a discretization of the Rayleigh and Stokes equation for a heated generalized second-grade fluid (HGSGF). Reference [29] presented the numerical 2 of 18 discretization for Stokes' first problem for the HGSGF described by a fractional operator. Reference [30] provided a discretization for Stokes' equation for the HGSGF that contains the Riemann-Liouville derivative. For further information about the HGSGF, one could see works in [1,[31][32][33]. In addition, there are many other investigations for fractional heat equations; some numerical studies exist as well. Reference [34] investigated a numerical computation for the wave diffusion problem in fractional context. Reference [35] illustrated the numerical technique for the diffusion problem with fractional order derivative. Reference [36] proposed the numerical approximation with high order accuracy for the fractional reaction-sub-diffusion equation. In [37], authors considered a heat conduction equation with respect to the Caputo-Fabrizio fractional derivative (CFFD). In [38], authors formulated a fractional optimal control problem for an anomalous diffusion process. We conclude the literature review by giving the next remark. Many other works in terms of analytical solution and numerical discretization for the fractional diffusion equation may also be found in the literature [35,39,40].
In this paper, we address the diffusion processes. We mainly investigate the solution of the incompressible second-grade fluids model which is constructed with the generalization of a fractional operator. The primary importance of this paper is the fractional aspect of the presented model. We apply the generalized fractional operator in mechanics fluids. Our paper yields providing the application of the fractional analysis in mathematical physics and mechanics fluids.
One of the main novelties of this paper is to propose new procedures which give the semi-analytical solutions of the nanofluid equations as well. The process of getting the solutions combines the heat integral balance method and the homotopy method (which is used classically for getting the approximate solutions for nanofluid equations). In our paper, we prove when we combine these two approaches, we obtain more accurate solutions for the fractional incompressible second-grade fluids models. In addition, the structure of the obtained solutions is more useful, and the graphical representations are most straightforward. As another contribution of the paper, we analyze the behavior of the obtained solutions when the order ρ is into (0, 1) or outside (0, 1). In other words, we consider the subdiffusion process and the superdiffusion process. In other words, what is the impact of the order ρ into the behavior of the solutions of the fractional the incompressible second-grade fluids models?
In Section 2, we offer the fractional-order derivatives that are used in our paper. In Section 3, we provide and give the model presentations which are constructed with a fractional operator. In Section 4, we present the fundamental qualitative properties of the fractional incompressible second-grade fluids model and its solutions as well. In Section 5, we describe the solution methods which are combination of two different methods that get the approximate solutions of the fractional incompressible second-grade fluids model. In Section 6, we provide the solutions of the mentioned problem. In Section 7, we illustrate our main results by the graphical representations and discuss the impact of the order ρ when we fix the fractional order ϕ. In addition, finally, we give the conclusions which we have pointed out in the paper and the corresponding perspectives in Section 8.

Basic Definitions on Fractional Derivations and Their Generalizations
In this section, we call to mind the fundamental definitions and features of the fractional analysis and their generalized form recently proposed in the literature by Tomovski et al. [41], Sene and Srivastava [42], Jarad and Abdeljawad [43,44] and Jarad et al. [45]. We first define the so-called fractional integral proposed by Riemann and Liouville. Let g(t) be a strictly increasing function with continuous derivative g on the interval (a, b). The left Riemann-Liouville fractional integral of f with respect to the function g of order ϕ, (ϕ) > 0 is defined by It is clear that when g(t) = t, Equation (1) is the classical Riemann-Liouville fractional integral. Taking the function f : [0, ∞[−→ R, we define the integral of order ϕ, of the function f starting at 0 as the following construction [46] ( where Γ(.) represents the Euler Gamma function and t > 0, and 0 < ϕ < 1. Its generalized form introduced in the literature by the authors in [41,[43][44][45] is described as follows. From the function defined by f : [0, ∞[−→ R, and taking g(t) = t ρ ρ in Equation (1) we obtain a special case of the the generalized integral of order ϕ, ρ > 0 of f starting at 0 as the following form We use these definitions above to define the Riemann-Liouville fractional operator in its classical form and its generalized form. The left Riemann-Liouville fractional derivative of f with respect to the function g of order ϕ, (ϕ) > 0 is defined by It can be easily noticed that when g(t) = t, Equation (4) is the classical Riemann-Liouville fractional derivative. By using f : [0, ∞[−→ R, we define the Riemann-Liouville operator of order ϕ, starting at 0 as the following form Its generalized form introduced in the literature by the authors in [41,[43][44][45] is described as follows. From the function defined by f : [0, ∞[−→ R, and taking g(t) = t ρ ρ in Equation (4) we obtain a special case of the the generalized Riemann-Liouville derivative of order ϕ, ρ > 0 of f starting at 0 as the following form.
In the same way, the Liouville-Caputo fractional derivative of f with respect to the function g of order ϕ, (ϕ) > 0 is defined by It can be easily noticed that when g(t) = t, Equation (7) gives the classical Liouville-Caputo fractional derivative which is given From the function defined by f : [0, ∞[−→ R, and taking g(t) = t ρ ρ in Equation (7) we obtain a special case of the the generalized Liouville-Caputo derivative of order ϕ, ρ > 0 of f starting at 0 as the following form: Remark 1. Each of these generalized forms given in Equation (6) and Equation (9) is derivable from the classical Riemann-Liouville derivative given by Equation (5) and the Liouville-Caputo fractional derivative given by Equation (8), respectively, by changing the variable, the parameter and the functional notation appropriately. In other words, In what follows, we recall the ρ-Laplace transform (ρ-LT) of a real-valued function and ρ-LT of the Liouville-Caputo operator which we take into account throughout the paper in our investigations. The ρ-Laplace transform of the real-valued function f : [0, ∞) → R is described by where for all s values the integral is valid.
Proof. The proof can be found in [43].
The ρ-LT of the Liouville-Caputo operator of order 0 < ϕ < 1 is defined by [43] with Corollary 3.3 as: is the space of absolutely continuous functions on [0, a] for any a > 0 and γ k f , k = 0, 1, . . . , n is ρ-exponential order e c t ρ ρ . Then the following Remark may be given: Remark 2. ρ-generalized variation of the Laplace transform is derivable from the classical Laplace transform itself by suitably changing the variable and the index and the functional notation.
The main relationship between the ρ-LT and the classical Laplace transform is given by [44] with such that the generalized Laplace transform of f exists and and ρ > 0. We now give the following particular relationship which we will use in our calculations This confirms that these fractional generalized operators that can be used to produce more general types of fractional derivatives with memory effect. Therefore, it is always of interest to introduce new local derivatives of arbitrary order and use them by a fractionalization process to produce new types of fractional derivatives of different kernels [42].
The Mittag-Leffler function with the parameters ϕ and ζ is presented with the following sum where ϕ > 0, ζ ∈ R and χ ∈ C. Note convergence of this series result from the assumptions ϕ > 0, and ζ > 0. The Mittag-Leffler approximation will be used for expressing the obtained approximate solution of the fractional incompressible second grade fluid.

Model Presentation and Work Project
In this part, we combine two useful methods of obtaining the solutions of FDEs. We propose in our investigations the approximate solution of the ISGF model which is covered with the left generalized fractional operator expressed by the equations PrD ϕ,ρ where the conditions at and furthermore, the function u represents the velocity of the fluid and θ represents the temperature of the fluid. Moreover, these functions satisfy the following relations In the above equation, represents the Grashof number, Pr = qC p k represents the Prandtl number and H represents the Heaviside function. f indicates the constant which has the velocity's dimension, C p represents the heat capacity at constant pressure, k represents the thermal conductivity and q the constant density. υ denotes the kinematic of the fluid, g shows the gravitation acceleration, β T represents the volumetric number of the thermal expansion, T w denotes the plate temperature and T ∞ represents the ambient fluid temperature of the plate [47]. The main importance of our problem is two-equations constitute it. The first equation is used in nanofluid, and the second is a fractional diffusion equation. The fractional diffusion Equation (12) represents an exogenous input for the first problem Equation (11). The governing Equations (11) and (12) model the heat transfer in a second grade fluid over and oscillating vertical plate. For more pieces of information and the graphical description of the model can be found in [47]. Furthermore, η = α 1 f 2 µυ is a constant where α 1 is the second grade parameter, and µ is the diffusion term. The solution of the model presented by Equations (11) and (12) can be approximated by many methods: the homotopy perturbed method and the homotopy analysis method. Equation (12) is a heat equation and many methods can be considered to determine its analytical or approximate solutions. We can use the Fourier sine transform, the Fourier transform, the integral balance method, and many others. Our motivation in this paper is to bring a more precise solutions of the model presented by Equations (11) and (12) by first solving Equation (11) using Homotopy perturbation method and the HBIM to give a more precise approximate solution for Equation (12). The impact of the order ρ will be analyzed.

Fundamental Qualitative Properties of the Solutions
In this section, we consider the Banach sketch (BS) to provide the conditions to obtain the qualitative properties to the incompressible SGF constructed with the left generalized fractional derivative.

Theorem 2.
The solution of the fractional differential equation described by Equation (11) exists.
Proof. We begin with the first equation of the proposed model. We take into account the following function.
Let us prove that Θ is Lipschitz and continuous (L-C). Here for simplification, we use the classical norm. By considering triangular inequality and taking the norm to both sides of Equation (15), we get Let us assume that u is Lipschitz and continuous. We obtain the following expression after calculation when we consider the second term null where a is a constant. From the condition u(x, 0) = 0, we provide the Picard's operator (PsO) as follows where Mu : H → H and H is a compact set. The operator M should be bounded before continuing our reasoning. We recall the Euclidean norm, we have the expression Using the assumption that Θ (x, t, θ) is bounded by b and t ≤ T, we have the inequality Equation (17) proves that the PsO is bounded. Let's prove this operator is now a contraction or we will provide a condition under which the PsO is a contraction. Applying the classical norm, we obtain the following inequality Using the fact that Θ is L-C (Equation (16)), we obtain the following relationships From which the PsO is a contraction if the following relation is provided We conclude the proof by taking into account the Banach Theorem. Thus, we prove the qualitative properties of the solution to Equation (11). In conclusion, we can now give a semi-analytical solution to Equation (11) because the problem of getting the exact solution is well defined. Theorem 3. The solution of the fractional differential equation described by Equation (12) is unique. (12), we repeat the same procedures. We consider the following function

Proof. For the heat Equation
Let us prove that Σ is Lipschitz and continuous with constant c. Here for simplification, we use the classical norm. By considering triangular inequality and taking the norm of Equation (18), we get Let us suppose that θ is Lipschitz continuous, and c is a constant such that the following relationship holds From initial condition θ(x, 0) = 0, we define the Picard's operator as follows where Nθ : H → H and H is a Banach space. The operator N should be bounded before continuing our reasoning. We recall the classical norm; we have the following expression By using the assumption that Σ(θ, x, t) is bounded by d and t ≤ T, we get the following inequality Equation (19) proves that the PsO is bounded. We prove that this operator is now a contraction or we will provide a condition under which the PsO is a contraction. If we apply the norm, we obtain the following inequality Using Σ(θ, x, t) is Lipschitz and continuous, we get the following connection From which the PsO is a contraction if the following relation is ensured Finally, this shows that the solution to the incompressible second grade fluid model which is given in Equations (11) and (12) is unique.
In conclusion, we can now give an approximate solution of Equation (12) because the problem of getting the exact solution is well defined.

Heat Balance Integral Method (HBIM)
In this subsection, we provide briefly the method that is used for obtaining the solution of the second equations. The method exists in the literature and provided in different studies as in [1,48]. According to the HBIM method, we suppose the solution is regarded as the following form where n represents the exponent and will be fixed to n = 2. However, as we will discuss later, there exist techniques to get the approximate value of the exponent n proposed by Myers and Mitchell [49]. The technique which is used in the study provides a physical concept. The main argument is to integrate the heat Equation (12) between 0 to the depth δ. That is The next step is to replace Equation (20) into Equation (21), and solve the equation using ρ-Laplace of the obtained equation. The objective is to get the form of the depth δ. The solution of the heat (12) is provided by replacing the expression of the δ into Equation (20).

ρ-Homotopy Perturbation Laplace Transformation
In this subsection, we define the ρ-Laplace homotopy transformation method (ρ-LHTM) to solve the problem mentioned in Section 3. This method is combined with the classical homotopy technique and ρ-Laplace transform and it was proposed firstly by Sene and Fall [50]. Consider the following PDE in the generalized fractional sense: with Dirichlet initial and boundary conditions given by and where D ϕ,ρ C shows the generalized-Caputo fractional operator (GCFO). Other variables stated in Equation (22) are the same with those which were defined in Section 3. Using the ρ-Laplace transform, we assign the L ρ {u(x, t)} = U(x, s). Then applying the homotopy to Equation (22) we derive the homotopies for the GCFO. Taking the LT of both sides of Equation (22), yields We assume that then substituting Equation (26) into Equation (25) and applying the homotopy steps, we have where u(x, t) mxx is mth iteration of the expression u(x, t) xx . If we compare the values of the powers of z, we generate the homotopies for the generalized Caputo fractional operator as follows: Then the corresponding solution of Equation (22) is given by

Procedure Solutions
In this Part, we give the solution of the incompressible fluid described by the left GDD defined in Equation (12). We begin the resolution by solving the fractional heat equation described by Equation (12). We use the HBIM which consists of expressing the similarity variable of the heat equation using the penetration depth denoted in this paper by δ. This technique is to integrate the heat diffusion equation from 0 to the depth δ. We suppose the suggested solution of Equation (12) is expressed as It follows the following identity Considering the integration in the interval (0, δ), we arrive to the following relation For the rest of the paper we take Pr = 1 and f = 1, because in the present investigation its influence in the dynamics of the considered model is not under consideretion. The LT to both sides of Equation (31), we get the relationships with the condition δ(0) = 0, Inverting Equation (32), we have the approximate equivalent δ 2 (t) = 2n(n + 1) Thus, using Equation (29), the approximate solution of the heat Equation (12) can be formed in the following expression For the heat diffusion equation, we have a parabolic equation; thus, in this study, we stipulate the exponent n = 2. There exist in the literature many discussions related to the good choose of the exponent for the diffusion and the fractional diffusion equations. In this context, Myers [51] provided an excellent method for getting the exponent n by minimizing the residual term of the diffusion equation. However, in many cases, we choose n = 2, to obtain a more accurate profile. In our study, the approximate solution of Equation (12) is clarified as the expression of The importance of our method is the physical propose. Our method gives the estimate of the similarity variable of the diffusion equation of fractional order. That is Another importance of our method, is it propose a classification of the nature of the diffusion processes. For classification, let's the square penetration depth Here we give a brief classification. Some detail will be give in forthcoming paper. Note when ϕ = ρ = 1, we have a normal diffusive. When ϕ < 1 ρ , we have a sub-diffusion process. When ϕ > 1 ρ , the diffusion is super-diffusion. In addition, when ϕ = 1 ρ , we have a ballistic diffusion. The solution of Equation (12) can be expressed in the following form where F ϕ = 12 Γ(1+ϕ) . Thus, the first Equation (11), will have been solved using the ρ-LHTM. Now considering the solution obtained in Equation (36) we aim to solve by homotopy technique explained in Section 5.2 the following equation For this aim, we apply the homotopy steps to the last equation. Then we have . . .
From the last homotopies and taking the inverse ρ-LT of each terms we can get the followings: Then the solution of the incompressible fluid equation is given by

Graphics and Discussions
In this part, we have the graphical representations of super-diffusion, ballistic-diffusion, and sub-diffusion cases. Also, we discuss the results of the solutions obtained in this study.
We first depict and analyze the solution of fractional diffusion equation (FDE) (12), which has been represented in Equation (33). We make the following assumptions t = 0.3 and ϕ = 1, and we depict the solution according to the state variable x of the diffusion process. In Figure 1, we demonstrate the solutions of FDE for different values of the order ρ ≤ 1. We observe the solutions decrease considerably and converge the normal diffusion obtained with the values (ϕ = 1 = ρ) when the order ρ increases to 1. We note the order ρ has a prominent effect on the diffusion processes. In general, it generates a retardation impact on the diffusion processes. Physically, the behaviors are explained by the sub-diffusion process generated by the fractional heat equation when ρ ≤ 1.
In Figure 2, we present the solutions of the heat equation for varied values of the order ρ satisfying the condition ρ ≥ 1. We note the solutions decrease according to the state variable x. Analyzing the behaviors of the solutions according to one another, we observe when the order increases and satisfies the condition ρ ≥ 1; all the curves increase as well and converge to the normal diffusion. The arrow in Figure 2 indicates these observations. The super-diffusion process explains the behaviors of the solutions when we fix ϕ and ρ ≥ 1. Here we note the order ρ generates an acceleration impact.
We fix x = 0.3 and ϕ = 1, and we depict the solutions according to time t of the diffusion process. We obtain the opposite behaviors, contrary to the previous cases. We note after certain times, all the curves increase and converge to infinity. See the example in Figure 3. Please note that the effects of the order ρ are the same as in previous analysis: retardation and acceleration effects.  In Figures 4-6, we point out the solutions with respect to the Super-diffusion, Ballistic-diffusion and Sub-diffusion cases, respectively. These mentioned graphical representations have been obtained by using the exact solution which is defined with the ρ-Laplace homotopy perturbation transform method and given in Equation (39). This mentioned method is very influential and accurate in finding the analytical solution of the ISGF equation of fractional order. Since, when looking at the solution, one can understand that by using only the first two iterations the solution has been provided. The truncation error does not need to have occurred for this problem. This situation may be possible for the solution of some kind of linear problems as well as it is may differ according to the nature of the problem [52].   In Figure 4, we obtain the super-diffusion-type process which is considered when ϕ < 1 ρ . In this graphical representation we took ϕ = 0.8 and ρ = 3. In Figure 5, we have the ballistic-type diffusion process which is obtained when ϕ = 1 ρ and here we have ϕ = 0.9 and ρ = 10/9. In Figure 6, we have regarded as sub-diffusion-type process which is taken into account when ϕ > 1 ρ . For this representation we take ϕ = 0.95 and ρ = 0.5.
In Figure 7, we can see the effect of diffusion parameter ρ for different values of the space variable when ϕ = 0.99 and t = 0.3. It is clear to recognize that when ρ values increase, the diffusion process approaches to the normal diffusion. In Figure 8, we represent the solution of the problem in Equation (11) for different values of the fractional parameter and diffusion term. It is concluded that both of the fractional parameter and diffusion term approach to the unit, the process designates to the normal diffusion case.

Concluding Remarks
In this study, the existence and uniqueness have been investigated for the FDD described by a fractional operator with its generalized form. The approximate solutions of the incompressible second-grade fluid have been provided. The ρ-Laplace homotopy perturbation transform method and the HBIM have been combined for proposing a new procedure to obtain the solution of the second-grade fluid model. In conclusion, our approach is useful and can be considered for the second grade fluids models, because the obtained solutions converge as well. They approach the exact solution as well. We represented the solutions of the equations composing the second-grade model studied in this paper. We pointed out the physical aspect of the considered model. We note the order ρ has retardation or acceleration impact in the diffusion processes. The numerical simulations and interpretations of the main results were presented. According to the numerical computations we have pointed out that the behavior of the temperature of the fluid has an important effect on the behavior of the velocity of the fluid. In other words, the type of the diffusion process in the fractional heat equation generates the same diffusion process in the fractional velocity equation. Another result in this paper concerns the exponent n = 2 for the heat balance integral method. The exponent n = 2 considered in this paper can be revised according to Myers method related to the exponent in future works. Future direction of investigation is to find a the best value of the exponent n for the semi-analytical solution of the second grade model considered in this paper using Myers method.