Approximate Solutions of the Model Describing Fluid Flow Using Generalized ρ-Laplace Transform Method and Heat Balance Integral Method
Approximate Solutions of the Model Describing Fluid Flow Using Generalized ρ-Laplace Transform...
Yavuz, Mehmet;Sene, Ndolane
axioms Article Approximate Solutions of the Model Describing Fluid Flow Using Generalized r-Laplace Transform Method and Heat Balance Integral Method 1,2, 3, Mehmet Yavuz * and Ndolane Sene * Department of Mathematics and Computer Sciences, Faculty of Science, Necmettin Erbakan University, 42090 Konya, Turkey Department of Mathematics, College of Engineering, Mathematics and Physical Sciences, University of Exeter, Cornwall TR10, UK Laboratoire Lmdan, Département de Mathématiques de la Décision, Faculté des Sciences Economiques et Gestion, Université Cheikh Anta Diop de Dakar, BP 5683 Dakar Fann, Senegal * Correspondence: email@example.com (M.Y.); firstname.lastname@example.org (N.S.) Received: 29 July 2020; Accepted: 21 October 2020; Published: 24 October 2020 Abstract: 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 r-Laplace homotopy transform method (r-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 r and j. Keywords: generalized fractional derivative; second-grade ﬂuid; r-Laplace homotopy perturbation transform method; heat integral balance method 1. Introduction Fractional calculus is a ﬁeld trying to understand the real-world phenomena modeled with non-integer-order derivatives [1,2], ﬁnance [3,4], biological processes and epidemic models [5–8], science and engineering [9–17], mechanics, etc. There exist at present numerous fractional operators which are deﬁned 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  and the Riemann–Liouville operator . Without singular kernels, we have two types: the fractional derivative with exponential kernel [19,20] and the fractional operator which has Mittag–Lefﬂer function as a kernel known as the Atangana–Baleanu operator [21–23]. We also summarized the existing FDs, which are the most popular in the literature. Some other types of derivatives as conformable , the Hilfer derivation , 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 brieﬂy enumerate some of them. Reference  proposed a discretization of the Rayleigh and Stokes equation for a heated generalized second-grade ﬂuid (HGSGF). Reference  presented the numerical Axioms 2020, 9, 123; doi:10.3390/axioms9040123 www.mdpi.com/journal/axioms Axioms 2020, 9, 123 2 of 18 discretization for Stokes’ ﬁrst problem for the HGSGF described by a fractional operator. Reference  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–33]. In addition, there are many other investigations for fractional heat equations; some numerical studies exist as well. Reference  investigated a numerical computation for the wave diffusion problem in fractional context. Reference  illustrated the numerical technique for the diffusion problem with fractional order derivative. Reference  proposed the numerical approximation with high order accuracy for the fractional reaction-sub-diffusion equation. In , authors considered a heat conduction equation with respect to the Caputo–Fabrizio fractional derivative (CFFD). In , 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 ﬂuids 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 ﬂuids. Our paper yields providing the application of the fractional analysis in mathematical physics and mechanics ﬂuids. One of the main novelties of this paper is to propose new procedures which give the semi-analytical solutions of the nanoﬂuid 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 nanoﬂuid equations). In our paper, we prove when we combine these two approaches, we obtain more accurate solutions for the fractional incompressible second-grade ﬂuids 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 r 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 r into the behavior of the solutions of the fractional the incompressible second-grade ﬂuids 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 ﬂuids 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 ﬂuids 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 r when we ﬁx the fractional order j. In addition, ﬁnally, we give the conclusions which we have pointed out in the paper and the corresponding perspectives in Section 8. 2. Basic Deﬁnitions on Fractional Derivations and Their Generalizations In this section, we call to mind the fundamental deﬁnitions and features of the fractional analysis and their generalized form recently proposed in the literature by Tomovski et al. , Sene and Srivastava , Jarad and Abdeljawad [43,44] and Jarad et al. . We ﬁrst deﬁne 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 j,< (j) > 0 is deﬁned by j 1 0 I f (t) = (g(t) g(l)) f (l)g (l)dl. (1) G(j) 0 Axioms 2020, 9, 123 3 of 18 It is clear that when g(t) = t, Equation (1) is the classical Riemann-Liouville fractional integral. Taking the function f : [0, ¥[ ! R, we deﬁne the integral of order j, of the function f starting at 0 as the following construction  j 1 j j,1 ( I f ) (t) = I f (t) = (t l) f (l)dl, (2) G(j) where G(.) represents the Euler Gamma function and t > 0, and 0 < j < 1. Its generalized form introduced in the literature by the authors in [41,43–45] is described as follows. From the function deﬁned by f : [0, ¥[ ! R, and taking g(t) = in Equation (1) we obtain a special case of the the generalized integral of order j, r > 0 of f starting at 0 as the following form Z j 1 t r r 1 t l dl j,r I f (t) = f (l) . (3) 1 r G(j) r l We use these deﬁnitions above to deﬁne 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 j,< (j) > 0 is deﬁned by 1 d dt j,r g (t) j 0 D f (t) = (g(t) g(l)) f (l)g (l)dl. (4) G(1 j) 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 deﬁne the Riemann-Liouville operator of order j, starting at 0 as the following form 1 d D f (t) = (t l) f (l)dl. (5) G(1 j) dt Its generalized form introduced in the literature by the authors in [41,43–45] is described as follows. From the function deﬁned by f : [0, ¥[ ! R, and taking g(t) = in Equation (4) we obtain a special case of the the generalized Riemann-Liouville derivative of order j, r > 0 of f starting at 0 as the following form. r r 1 d t l dl j,r 1 r D f (t) = t f (l) . (6) 1 r G(1 j) dt r 0 l In the same way, the Liouville-Caputo fractional derivative of f with respect to the function g of order j,< (j) > 0 is deﬁned by n 1 [k] + f (a ) j,r j,r D f (t) = D f (s) (g(s) g(0)) (t) . (7) g å k! k=0 It can be easily noticed that when g(t) = t, Equation (7) gives the classical Liouville-Caputo fractional derivative which is given j 0 D f (t) = (t l) f (l)dl. (8) G(1 j) From the function deﬁned by f : [0, ¥[ ! R, and taking g(t) = in Equation (7) we obtain a special case of the the generalized Liouville-Caputo derivative of order j, r > 0 of f starting at 0 as the following form: r r 1 t l dl j,r D f (t) = g f (l) . (9) 1 r G(1 j) r 0 Axioms 2020, 9, 123 4 of 18 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. j,1 j j,1 j In other words, D f (t) = D f (t) and D f (t) = D f (t). C C In what follows, we recall the r-Laplace transform (r-LT) of a real-valued function and r-LT of the Liouville-Caputo operator which we take into account throughout the paper in our investigations. The r-Laplace transform of the real-valued function f : [0, ¥) ! R is described by dt L f f (t)g (s) = e f (t) , r > 0, 1 r 0 t where for all s values the integral is valid. Theorem 1. Let the function f (t) be continuous and of exponential order e such that g f (t) is piecewise continuous over every ﬁnite interval [0, T] . Then r-Laplace transform of g f (t) exists for s > c and L fg f (t)g (s) = sL f f (t)g f (0), (10) r r 1 r where r > 0 and g = t . dt Proof. The proof can be found in . The r-LT of the Liouville-Caputo operator of order 0 < j < 1 is deﬁned by  with Corollary 3.3 as: " # n o n 1 j,r j k 1 k L D f (t) (s) = s L f f (t)g s g f (0) , s > c, r r 0 å k=0 where j > 0, f 2 AC [0, a] is the space of absolutely continuous functions on [0, a] for any a > 0 and g f , k = 0, 1, . . . , n is r-exponential order e . Then the following Remark may be given: Remark 2. r-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 r-LT and the classical Laplace transform is given by  with Theorem 3.2. n o L f f (t)g (s) = L f g (t + g(a)) (s). where f , g : [a, ¥) ! R are the real-valued functions such that g (t) is continuous and g (t) > 0 on [0, ¥) such that the generalized Laplace transform of f exists and and r > 0. We now give the following particular relationship which we will use in our calculations G 1 + L ft g (s) = r , p 2 R, s > 0. r p 1+ Remark 3. As the usual Laplace transform (1-Laplace) is a tool to solve classical fractional Riemann-Liouville and Liouville-Caputo derivatives, we use the r-Laplace transform to solve the incompressible second-grade ﬂuids model in the frame of Riemann-Liouville and Liouville-Caputo type fractional generalized operators. This conﬁrms 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 Axioms 2020, 9, 123 5 of 18 of arbitrary order and use them by a fractionalization process to produce new types of fractional derivatives of different kernels . The Mittag-Lefﬂer function with the parameters j and z is presented with the following sum E (c) = , j,z å G(jx + z) x=0 where j > 0, z 2 R and c 2 C. Note convergence of this series result from the assumptions j > 0, and z > 0. The Mittag-Lefﬂer approximation will be used for expressing the obtained approximate solution of the fractional incompressible second grade ﬂuid. 3. 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 ¶ u j,r j,r D u = 1 + h D + Grq, (11) C C 2 ¶x ¶ q j,r PrD q = , (12) ¶x where the conditions at t = 0 are u(x, 0) = q(x, 0) = 0, (13) and furthermore, the function u represents the velocity of the ﬂuid and q represents the temperature of the ﬂuid. Moreover, these functions satisfy the following relations u(0, t) = f H(t) sin wt, and q(0, t) = 1. (14) qC ugb (T T ) p w ¥ In the above equation, Gr = represents the Grashof number, Pr = represents the Prandtl number and H represents the Heaviside function. f indicates the constant which has the velocity’s dimension, C represents the heat capacity at constant pressure, k represents the thermal conductivity and q the constant density. u denotes the kinematic of the ﬂuid, g shows the gravitation acceleration, b represents the volumetric number of the thermal expansion, T denotes T w the plate temperature and T represents the ambient ﬂuid temperature of the plate . The main importance of our problem is two-equations constitute it. The ﬁrst equation is used in nanoﬂuid, and the second is a fractional diffusion equation. The fractional diffusion Equation (12) represents an exogenous input for the ﬁrst problem Equation (11). The governing Equations (11) and (12) model the heat transfer in a second grade ﬂuid over and oscillating vertical plate. For more pieces of information a f and the graphical description of the model can be found in . Furthermore, h = is a constant mu where a is the second grade parameter, and m 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 ﬁrst 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 r will be analyzed. Axioms 2020, 9, 123 6 of 18 4. 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 ﬁrst equation of the proposed model. We take into account the following function. 2 2 ¶ u ¶ u j,r Q(u, x, t, q) = + h D + Grq(x, t). (15) 2 2 ¶x ¶x Let us prove that Q is Lipschitz and continuous (L-C). Here for simpliﬁcation, we use the classical norm. By considering triangular inequality and taking the norm to both sides of Equation (15), we get 2 2 2 2 ¶ u ¶ v ¶ u ¶ v j,r kQ(u, x, t, q) Q(v, x, t, q)k + h D . 2 2
2 2 ¶x ¶x ¶x ¶x Let us assume that u is Lipschitz and continuous. We obtain the following expression after calculation when we consider the second term null 2 2 2 2 ¶ u ¶ v ¶ u ¶ v j,r kQ(u, x, t, q) Q(v, x, t, q)k + h D 2 2 2 2 ¶x ¶x ¶x ¶x j,r a u v + ah D u v (16) k k k k aku vk , where a is a constant. From the condition u(x, 0) = 0, we provide the Picard’s operator (PsO) as follows j,r Mu(x, t) = I Q(u, x, t, q), 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 j,r k Mu(x, t) u(x, 0)k = k I Q(u, x, t, q)k j,r I kQ(u, x, t, q)k j,r kQ(u, x, t, q)k I (1). Using the assumption that kQ (x, t, q)k is bounded by b and t T, we have the inequality 1 j r r T k Mu(x, t) u(x, 0)k b. (17) G(j) r 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 j,r k Mu(x, t) Mv(x, t)k = k I Q(u, x, t, q) Q(v, x, t, q)k j,r kQ(u, x, t, q) Q(v, x, t, q)k I (1). Axioms 2020, 9, 123 7 of 18 Using the fact that Q is L-C (Equation (16)), we obtain the following relationships 1 j r r T k Mu(x, t) Mv(x, t)k aku vk . G(j + 1) r From which the PsO is a contraction if the following relation is provided 1 j r r T 1 < . G(j + 1) r a 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 deﬁned. Theorem 3. The solution of the fractional differential equation described by Equation (12) is unique. Proof. For the heat Equation (12), we repeat the same procedures. We consider the following function ¶ q S(q, x, t) = . (18) ¶x Let us prove that S is Lipschitz and continuous with constant c. Here for simpliﬁcation, we use the classical norm. By considering triangular inequality and taking the norm of Equation (18), we get 2 2 0 ¶ q ¶ q S(q, x, t) S(q , x, t) . 2 2 ¶x ¶x Let us suppose that q is Lipschitz continuous, and c is a constant such that the following relationship holds 2 2 0 ¶ q ¶ q S(q, x, t) S(q , x, t) 2 2 ¶x ¶x c q q . From initial condition q(x, 0) = 0, we deﬁne the Picard’s operator as follows j,r Nq(x, t) = I S(q, x, t), where Nq : 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 j,r kNq(x, t) q(x, 0)k = k I S(q, x, t)k j,r I kS(q, x, t)k j,r kS(q, x, t)k I (1). By using the assumption that S(q, x, t) is bounded by d and t T, we get the k k following inequality 1 j r r T kNq(x, t) q(x, 0)k d. (19) G(j) r Axioms 2020, 9, 123 8 of 18 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 0 j,r 0 Nq(x, t) Nq (x, t) = I S(q, x, t) S(q , x, t) j,r 0 I S(q, x, t) S(q , x, t) 0 j,r S(q, x, t) S(q , x, t) I (1). Using S(q, x, t) is Lipschitz and continuous, we get the following connection 1 j r r T 0 0 Nq(x, t) Nq (x, t) c q q . G(j + 1) r From which the PsO is a contraction if the following relation is ensured 1 j r r T 1 < . G(j + 1) r c Finally, this shows that the solution to the incompressible second grade ﬂuid 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 deﬁned. 5. Description of Proposed Solution Methods 5.1. Heat Balance Integral Method (HBIM) In this subsection, we provide brieﬂy 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 q(x, t) = 1 . (20) where n represents the exponent and will be ﬁxed 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 . 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 d. That is Z Z d d ¶ q(x, t) j,r Pr D q(x, t)dx = dx. (21) ¶x 0 0 The next step is to replace Equation (20) into Equation (21), and solve the equation using r-Laplace of the obtained equation. The objective is to get the form of the depth d. The solution of the heat (12) is provided by replacing the expression of the d into Equation (20). 5.2. r-Homotopy Perturbation Laplace Transformation In this subsection, we deﬁne the r-Laplace homotopy transformation method (r-LHTM) to solve the problem mentioned in Section 3. This method is combined with the classical homotopy technique and r-Laplace transform and it was proposed ﬁrstly by Sene and Fall . Consider the following PDE in the generalized fractional sense: ¶ u(x, t) j,r j,r D u(x, t) = 1 + h D + Grq(x, t), (22) C C 2 ¶x Axioms 2020, 9, 123 9 of 18 with Dirichlet initial and boundary conditions given by u(x, 0) = 0, (23) and u(0, t) = f H(t) sin wt, (24) j,r where D shows the generalized-Caputo fractional operator (GCFO). Other variables stated in Equation (22) are the same with those which were deﬁned in Section 3. Using the r-Laplace transform, we assign the L fu(x, t)g = 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 1 1 1 j,r (25) U(x, s) = L f 1 + h D u(x, t) g + L fGrq(x, t)g + u(x, 0). r xx r j C j s s s We assume that U(x, s) = z U (x, s), (26) å m m=0 then substituting Equation (26) into Equation (25) and applying the homotopy steps, we have " # ¥ n ¥ o j,r m m z U (x, s) = z L z 1 + h D u(x, t) m r mxx å å j C m=0 m=0 1 1 + L fGrq(x, t)g + u(x, 0), (27) s s where u(x, t) is mth iteration of the expression u(x, t) . If we compare the values of the powers mxx xx of z, we generate the homotopies for the generalized Caputo fractional operator as follows: 1 1 z : U (x, s) = L fGrq(x, t)g + u(x, 0), 0 r s s j,r z : U (x, s) = L f 1 + h D u(x, t) g, 1 r 0xx j C j,r z : U (x, s) = L f 1 + h D u(x, t) g, 2 1xx j,r n+1 z : U (x, s) = L f 1 + h D u(x, t) g. n+1 r nxx j C Then the corresponding solution of Equation (22) is given by n o u(x, t) = L U (x, s) . (28) r å m m=0 6. Procedure Solutions In this Part, we give the solution of the incompressible ﬂuid described by the left GDD deﬁned 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 d. This technique is to integrate the heat diffusion equation from 0 to the depth d. We suppose the suggested solution of Equation (12) is expressed as q(x, t) = 1 . (29) d Axioms 2020, 9, 123 10 of 18 It follows the following identity Z Z d d ¶ q(x, t) j,r Pr D q(x, t)dx = dx, C 2 ¶x 0 0 j,r PrD q(x, t)dx = , (30) Pr n j,r D d = , n + 1 d n(n + 1) j,r PrD d = . Considering the integration in the interval (0, d), we arrive to the following relation Z Z d d n(n + 1) j,r Pr D ddx = dx, 0 0 j,r PrD d = 2n(n + 1). (31) For the rest of the paper we take Pr = 1 and f = 1, because in the present investigation its inﬂuence 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 d(0) = 0, 2n(n + 1) j 2 s d (s) = , 2n(n + 1) d (s) = . (32) 1+j Inverting Equation (32), we have the approximate equivalent 2n(n + 1) t d (t) = , G (1 + j) r j/2 2n(n + 1) t d(t) = . G (1 + j) r Thus, using Equation (29), the approximate solution of the heat Equation (12) can be formed in the following expression 0 1 B C B C q(x, t) = 1 . @ A j/2 2n(n+1) G(1+j) 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  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 proﬁle. In our study, the approximate solution of Equation (12) is clariﬁed as the expression of 0 1 B C q(x, t) = 1 . (33) @ A j/2 3 t G(1+j) Axioms 2020, 9, 123 11 of 18 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 . (34) j/2 3 t G 1+j r ( ) Another importance of our method, is it propose a classiﬁcation of the nature of the diffusion processes. For classiﬁcation, let’s the square penetration depth 2n(n + 1) t d (t) = . (35) G (1 + j) r Here we give a brief classiﬁcation. Some detail will be give in forthcoming paper. Note when 1 1 j = r = 1, we have a normal diffusive. When j < , we have a sub-diffusion process. When j > , r r the diffusion is super-diffusion. In addition, when j = , we have a ballistic diffusion. The solution of Equation (12) can be expressed in the following form j/2 j r r t t 1/2 1 2 q(x, t) = 1 2F x + F x , (36) j j r r where F = . Thus, the ﬁrst Equation (11), will have been solved using the r-LHTM. G(1+j) Now considering the solution obtained in Equation (36) we aim to solve by homotopy technique explained in Section 5.2 the following equation j/2 j 2 r r ¶ u t t j,r j,r 1/2 1 2 D u(x, t) = 1 + h D + 1 2F x + F x . (37) j j C C ¶x r r For this aim, we apply the homotopy steps to the last equation. Then we have j/2 j n o r r 1 t t 0 1/2 1 2 z : U (x, s) = L 1 2F x + F x , 0 j j s r r 1/r n o j,r 1 1/r r 1 z : U (x, s) = L 1 + h D 2.12 G(1 j) r G(j + 1) , 1 r j C j,r z : U (x, s) = L f 1 + h D u(x, t) g = 0, (38) 2 r 1xx j C j,r n+1 z : U (x, s) = L f 1 + h D u(x, t) gg = 0. n+1 r nxx j C s Axioms 2020, 9, 123 12 of 18 From the last homotopies and taking the inverse r-LT of each terms we can get the followings: j 1/r 1 r 1 2 j jr ( ) 2r 2 2 3 xr G 1 t 2 1 G(j+1) u (x, t) = G + 1 1/r j jr r t 1/r 2 r 1 +12 x G(1 j) r G(j + 1) + , G(j + 1) 1/r 1/r j r 1 jr 2 3 r G(1 j) r G(j + 1) t u (x, t) = , G(j + 1) j,r u (x, t) = L f 1 + h D u(x, t) g = 0, 2 r 1xx j C j,r u (x, t) = L f 1 + h D u(x, t) g = 0. n+1 r nxx j C Then the solution of the incompressible ﬂuid equation is given by 0 1 j 1/r r 1 2 j j 1 j/2 (jr) 2r @ A 2 u(x, t) = 3 xr G + 1 G 1 t (39) 2 2 G(j+1) 1/r j jr r t 1/r 2 r 1 +12 x G(1 j) r G(j + 1) + , G(j + 1) 2 1/r 1 j 1/r j r 1 jr 2 3 r G(1 j) r G(j + 1) t + . G(j + 1) 7. 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 ﬁrst 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 j = 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 r 1. We observe the solutions decrease considerably and converge the normal diffusion obtained with the values (j = 1 = r) when the order r increases to 1. We note the order r 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 r 1. In Figure 2, we present the solutions of the heat equation for varied values of the order r satisfying the condition r 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 satisﬁes the condition r 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 ﬁx j and r 1. Here we note the order r generates an acceleration impact. We ﬁx x = 0.3 and j = 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 inﬁnity. See the example in Figure 3. Please note that the effects of the order r are the same as in previous analysis: retardation and acceleration effects. Axioms 2020, 9, 123 13 of 18 0.9 0.8 =0.45 0.7 =0.5 0.6 0.5 =0.75 0.4 =0.85 0.3 =0.95 0.2 0.1 =1 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 Figure 1. Solution of Equation (12) for q versus x when j = 1 and r 1. 0.9 0.8 =1.45 0.7 =1.5 0.6 0.5 =1.75 0.4 =1.85 0.3 =1.95 0.2 =1 0.1 0 0.05 0.1 0.15 0.2 0.25 0.3 Figure 2. Solution of Equation (12) when j = 1 and r 1. 0.8 0.7 =0.45 0.6 =0.5 0.5 =0.75 0.4 =0.85 0.3 =0.95 0.2 =1 0.1 0 0.5 1 1.5 2 2.5 3 Figure 3. Solution of Equation (12) for q versus t when j = 1 and r 1. 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 deﬁned with the r-Laplace homotopy perturbation transform method and given in Equation (39). This mentioned method is very inﬂuential and accurate in ﬁnding the analytical solution of the ISGF equation of fractional order. Since, when looking at the solution, THETA Theta THETA Axioms 2020, 9, 123 14 of 18 one can understand that by using only the ﬁrst 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 . Figure 4. Super-diffusion case of incompressible second grade ﬂuid dynamics. Figure 5. Ballistic-diffusion case of incompressible second grade ﬂuid dynamics. Figure 6. Sub-diffusion case of incompressible second grade ﬂuid dynamics. In Figure 4, we obtain the super-diffusion-type process which is considered when j < . In this graphical representation we took j = 0.8 and r = 3. In Figure 5, we have the ballistic-type Axioms 2020, 9, 123 15 of 18 diffusion process which is obtained when j = and here we have j = 0.9 and r = 10/9. In Figure 6, we have regarded as sub-diffusion-type process which is taken into account when j > . For this representation we take j = 0.95 and r = 0.5. In Figure 7, we can see the effect of diffusion parameter r for different values of the space variable when j = 0.99 and t = 0.3. It is clear to recognize that when r values increase, the diffusion process approaches to the normal diffusion. Figure 7. Solution of Equation (11) when j = 0.99 and t = 0.3. 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. Figure 8. Solution of Equation (11) for various values of j and r. 8. 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 Axioms 2020, 9, 123 16 of 18 second-grade ﬂuid have been provided. The r-Laplace homotopy perturbation transform method and the HBIM have been combined for proposing a new procedure to obtain the solution of the second-grade ﬂuid model. In conclusion, our approach is useful and can be considered for the second grade ﬂuids 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 r 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 ﬂuid has an important effect on the behavior of the velocity of the ﬂuid. 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 ﬁnd 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. Author Contributions: Conceptualization, M.Y. and N.S.; methodology, N.S.; investigation, M.Y. and N.S.; writing–original draft preparation, M.Y. and N.S.; writing–review and editing, M.Y. and N.S.; visualization, M.Y. and N.S. All authors have read and agreed to the published version of the manuscript. Funding: This research received no external funding. Acknowledgments: M. Yavuz was supported by TUBITAK (The Scientiﬁc and Technological Research Council of Turkey). Conﬂicts of Interest: The authors declare no conﬂict of interest. References 1. Hristov, J. Approximate solutions to fractional subdiffusion equations. Eur. Phys. J. Spec. Top. 2011, 193, 229–243. 2. Darzi, R.; Agheli, B. An analytical approach for systems of fractional differential equations by means of the innovative homotopy perturbation method. Mathematica Moravica. 2018, 22, 93–105. 3. Yavuz, M.; Ozdemir, N. European vanilla option pricing model of fractional order without singular kernel. Fractal Fract. 2018, 2, 3. 4. Bas, E.; Acay, B.; Ozarslan, R. The price adjustment equation with different types of conformable derivatives in market equilibrium. AIMS Math. 2019, 47, 805. 5. Yavuz, M.; Bonyah, E. New approaches to the fractional dynamics of schistosomiasis disease model. Phys. A Stat. Mech. Appl. 2019, 525, 373–393. 6. Naik, P.A.; Owolabi, K.M.; Yavuz, M.; Zu, J. Chaotic dynamics of a fractional order HIV-1 model involving AIDS-related cancer cells. Chaos Solitons Fract. 2020, 140, 110272. 7. Yavuz, M.; Sene, N. Stability Analysis and Numerical Computation of the Fractional Predator–Prey Model with the Harvesting Rate. Fract. Fract. 2020, 4, 35. 8. Yavuz, M.; Ozdemir, N. Analysis of an epidemic spreading model with exponential decay law. Math. Sci. Appl. E-Notes 2020, 8, 142–154. 9. Santos, M.D. Fractional Prabhakar Derivative in Diffusion Equation with Non-Static Stochastic Resetting. Physics 2019, 1, 40–58. 10. Hristov, J. The heat radiation diffusion equation: Explicit analytical solutions by improved integral-balance method. Ther. Sci. 2018, 22, 777–788. 11. Khader, M.M. On the numerical solutions for the fractional diffusion equation. Commun. Nonlinear Sci. Numer. Simul. 2011, 16, 2535–2542. 12. Tasbozan, O.; Esen, A.; Yagmurlu, N.M.; Ucar, Y. A Numerical Solution to Fractional Diffusion Equation for Force-Free Case. Abstr. Appl. Anal. 2013, 2013, 6. 13. Yokus, A.; Bulut, H. On the numerical investigations to the Cahn-Allen equation by using ﬁnite difference method. Int. J. Optim. Control Theor. Appl. 2018, 9, 18–23. Axioms 2020, 9, 123 17 of 18 14. Sene, N. Second-grade ﬂuid model with Caputo–Liouville generalized fractional derivative. Chaos Solitons Fract. 2020, 133, 109631. 15. Pendra, S.K.; Abdeljawad, T.; Ravichandran, C.; Jarad, F. A complex valued approach to the solution of Riemann-Liouville integral, Atangana-Baleanu integral operator and non-linear Telegraph equation via ﬁxed point method. Chaos Solitons Fract. 2020, 130, 109439. 16. Evirgen, F.; Yavuz, M. An alternative approach for nonlinear optimization problem with Caputo-Fabrizio derivative. ITM Web Conf. 2018, 22, 01009. 17. Usta, F.; Budak, H.; Sarikaya, M.Z. Yang-Laplace Transform Method Volterra and Abels Integro-Differential Equations of Fractional Order. Int. J. Nonlinear Anal. Appl. 2017, 9, 203–214. 18. Caputo, M. Linear models of dissipation whose Q is almost frequency independent—II. Geophys. J. Int. 1967, 13, 529–539. 19. Caputo, M.; Fabrizio, M. A new deﬁnition of fractional derivative without singular kernel. Progr. Fract. Differ. Appl. 2015, 1, 1–15. 20. Hristov, J. Transient heat diffusion with a non-singular fading memory: From the Cattaneo constitutive equation with Jeffrey’s kernel to the Caputo-Fabrizio time-fractional derivative. Ther. Sci. 2016, 20, 757–762. 21. Atangana, A.; Baleanu, D. New fractional derivatives with nonlocal and non-singular kernel: Theory and application to heat transfer model. Ther. Sci. 2016, 20, 763–769. 22. Abro, K.A.; Gomez-Aguilar, J.F. A comparison of heat and mass transfer on a Walter ’s-B ﬂuid via Caputo-Fabrizio versus Atangana-Baleanu fractional derivatives using the Fox-H function. Eur. Phys. J. Plus 2019, 134, 101. 23. Yavuz, M.; Abdeljawad, T. Nonlinear regularized long-wave models with a new integral transformation applied to the fractional derivative with power and Mittag-Lefﬂer kernel. Adv. Differ. Equ. 2020, 2020, 367. 24. Khalil, R.; Al Horani, M.; Yousef, A.; Sababheh, M. A new deﬁnition of fractional derivative. J. Comput. Appl. Math. 2014, 264, 65–70. 25. Hilfer, R. Applications of Fractional Calculus in Physics; World Scientiﬁc: Singapore, 1999; pp. 87–130. 26. Abdeljawad, T.; Al-Mdallal, Q.M. Discrete Mittag-Lefﬂer kernel type fractional difference initial value problems and Gronwall’s inequality. J. Comput. Appl. Math. 2018, 339, 218–230. 27. Abdeljawad, T.; Baleanu, D. Discrete fractional differences with non-singular discrete Mittag-Lefﬂer kernels. Adv. Differ. Equ. 2016, 2016, 232. 28. Zhuang, P.H.; Liu, Q.X.; Numerical method of Rayleigh-Stokes problem for heated generalized second grade ﬂuid with fractional derivative. Appl. Math. Mech. 2009, 30, 1533–1546. 29. Wu, C. Numerical solution for Stokes’ ﬁrst problem for a heated generalized second grade ﬂuid with fractional derivative. Appl. Num. Math. 2009, 59, 2571–2583. 30. Ye, C.; Luo, X.; Wen, L. High-order numerical methods of fractional-order Stokes’ ﬁrst problem for heated generalized second grade ﬂuid. Appl. Math. Mech. Engl. 2012, 33, 65–80. 31. Hristov, J. Integral-balance solution to nonlinear subdiffusion equation. In Frontiers in Fractional Calculus, 1st ed.; Bhalekar, S., Ed.; Bentham Science Publishers: Sharjah, UAE, 2017; Volume 1, pp. 71–106. 32. Hashemi, M.S.; Baleanu, D.; Haghighi, M.P. Solving the time fractional diffusion equation using a lie group integrator. Ther. Sci. 2015, 19, 77–83. 33. Meerschaert, M.M.; Tadjeran, C. Finite difference approximations for fractional advection-dispersion ﬂow equations. J. Comput . Appl. Math. 2004, 172, 65–67. 34. Delic, A. Fractional in Time Diffusion-Wave Equation and its Numerical Approximation. Filomat 2016, 30, 1375–1385. 35. Sene, N. Analytical solutions and numerical schemes of certain generalized fractional diffusion models. Eur. Phys. J. Plus 2019, 134, 199. 36. Bhrawy, A.H.; Baleanu, D.; Mallawi, F. A new numerical technique for solving fractional sub-diffusion and reaction sub-diffusion equations with a nonlinear source term. Ther. Sci. 2015, 19, 25–34. 37. Avcı, D.; Yavuz, M.; Ozdemir, N. Fundamental Solutions to the Cauchy and Dirichlet Problems for a Heat Conduction Equation Equipped with the Caputo-Fabrizio Differentiation. In Heat Conduction: Methods, Applications and Research, 1st ed.; Hristov, J., Bennacer, R., Eds.; Nova Science Publishers: New York, NY, USA, 2019; Volume 1, pp. 95–107. Axioms 2020, 9, 123 18 of 18 38. Avcı, D.; Ozdemir, N.; Yavuz, M. Fractional Optimal Control of Diffusive Transport Acting on a Spherical Region. In Methods of Mathematical Modelling: Fractional Differential Equations, 1st ed.; Singh, H., Kumar, D., Baleanu, D., Eds.; CRC Press: Boca Raton, FL, USA, 2019; Volume 1, pp. 63–82. 39. Hristov, J. A transient ﬂow of a non-newtonian ﬂuid modelled by a mixed time-space derivative: An improved integral-balance approach. In Mathematical Methods in Engineering, 1st ed.; Tas, K., Baleanu, D., Tenreiro Machado, J.A., Eds.; Springer: Cham, Switzerland, 2019; Volume 1, pp. 153–174. 40. Ravichandran C.; Logeswari K.; Jarad, F. New results on existence in framework of Atangana-Baleanu derivative for fractional integro-differential equations. Chaos Solitons Fract. 2019, 125, 194–200. 41. Tomovski, Z.; Hilfer, R.; Srivastava, H.M. Fractional and operational calculus with generalized fractional derivative operators and Mittag–Lefﬂer type functions. Integral Transform. Spec. Funct. 2010, 21, 797–814. 42. Sene, N.; Srivastava, G. Generalized Mittag-Lefﬂer input stability of the fractional differential equations. Symmetry 2019, 11, 608. 43. Jarad, F.; Abdeljawad, T. A modiﬁed Laplace transform for certain generalized fractional operators. Res. Nonlinear Anal. 2018, 2, 88–98. 44. Jarad, F.; Abdeljawad, T. Generalized fractional derivatives and Laplace transform. Discret. Contin. Dyn. Syst.-S 2020, 13, 709–722. 45. Jarad, F.; Abdeljawad, T.; Baleanu, D. On the generalized fractional derivatives and their Caputo modiﬁcation. J. Nonlinear Sci. Appl. 2017, 10, 2607–2619. 46. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006; Volume 204. 47. Shah, N.A.; Khan, I. Heat transfer analysis in a second grade ﬂuid over and oscillating vertical plate using fractional Caputo–Fabrizio derivatives. Eur. Phys. J. C 2016, 76, 362. 48. Sene, N. Integral Balance Methods for Stokes’ First, Equation Described by the Left Generalized Fractional Derivative. Physics 2019, 1, 154–166. 49. Mitchell, S.L.; Myers, T.G. Improving the accuracy of heat balance integral methods applied to thermal problems with time dependent boundary conditions. Int. J. Heat Mass Transf. 2010, 53, 3540–3551. 50. Sene, N.; Fall, A.N. Homotopy Perturbation r-Laplace Transform Method and Its Application to the Fractional Diffusion Equation and the Fractional Diffusion-Reaction Equation. Fract. Fract. 2019, 3, 14. 51. Myers, T.G. Optimal exponent heat balance and reﬁned integral methods applied to Stefan problems. Int. J. Heat Mass Transf. 2010, 53, 1119–1127. 52. Yavuz, M.; Ozdemir, N. Numerical inverse Laplace homotopy technique for fractional heat equations. Ther. Sci. 2018, 22, 185–194. Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional afﬁliations. 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.pngAxiomsMultidisciplinary Digital Publishing Institutehttp://www.deepdyve.com/lp/multidisciplinary-digital-publishing-institute/approximate-solutions-of-the-model-describing-fluid-flow-using-liyV4tV4wk
Approximate Solutions of the Model Describing Fluid Flow Using Generalized ρ-Laplace Transform Method and Heat Balance Integral Method