# Dufour Effect on Transient MHD Double Convection Flow of Fractionalized Second-Grade Fluid with Caputo–Fabrizio Derivative

Authors: Imran Siddique; Sehrish Ayaz; Fahd Jarad
Journal: Complexity
 (2021)

---

## Abstract

This article presents the problem, in which we study the unsteady double convection flow of a magnetohydrodynamics (MHD) differential-type fluid flow in the presence of heat source, Newtonian heating, and Dufour effect over an infinite vertical plate with fractional mass diffusion and thermal transports. The constitutive equations for the mass flux and thermal flux are modeled for noninteger-order derivative Caputo–Fabrizio (CF) with nonsingular kernel, respectively. The Laplace transform and Laplace inversion numerical algorithms are used to derive the analytical and semianalytical solutions for the dimensionless concentration, temperature, and velocity fields. Expressions for the skin friction and rates of heat and mass transfer from the plate to fluid with noninteger and integer orders, respectively, are also determined. Furthermore, the influence of flow parameters and fractional parametersα and β on the concentration, temperature, and velocity fields are tabularly and graphically underlined and discussed. Furthermore, a comparison between second-grade and viscous fluids for noninteger and integer is also depicted. It is observed that integer-order fluids have greater velocities than noninteger-order fluids. This shows how the fractional parameters affect the fluid flow.


---

## Body

## 1. Introduction

The interest in fluid mechanics is truly significant within the sight of transport phenomena, which is a critical element in thermal, chemical, and mechanical engineering science. A few actual systems exist which can be utilized to move thermal energy and compound species through a phase and across limits of the phase. The three mechanisms for heat transfer are diffusion, convection, and radiation. The classification of convection of heat transfer into three consequent branches are natural (free), forced, and mixed convection, which is essential for the physical system that takes up the motion of the fluid. Free convection flows ensuing from the heat and mass transfer directed by the combined buoyancy effects because of temperature and concentration variations have been widely studied due to their applications in geotechnical engineering and chemical and bioengineering and in industrial activities [1]. Usually, the mass transfer due to the concentration disparity influences the rate of heat transfer. The driving force for the free convection is buoyancy, so its effects cannot be neglected whether the velocity of the fluid is small and change in temperature between the ambient fluid and surface is large enough [2–4].Electrically conducting fluids also have accepted enough consideration from the researchers due to their extensive applications in industrial appliances. The MHD has its own practical implication, such as the tumor treating fields and power generation and earthquake assumption [5]. Parvin and Nasrin [6] have presented the analysis of the flow and heat transfer characteristics for MHD-free convection in inclusion with heated difficulties. They showed that the influence of the magnetic parameter on streamlines and isotherms is significant.The energy flux which is due to a composition gradient is said to be the Dufour or diffusion-thermo effect. Such influences are significant when density differences occur in the flow regime. Such as when species are introduced at a surface in the fluid domain with different (lower) density than the surrounding fluid, Dufour effects can be beneficial. Also, when heat and mass transfers take place simultaneously in a moving fluid, it has a relationship between the fluxes and the driving potentials are of more twisting nature. It has been analyzed that an energy flux can be generated not only by temperature gradients but also by composition gradients as well. The diffusion-thermo effect was found to be of a considerable magnitude such that it cannot be negligible [7]. Dufour effects are essential in geothermal energy, hydrology, and nuclear waste disposal. In view of the importance of the diffusion-thermo effect, Kafoussias and Williams [8] studied the effect of thermal diffusion and diffusion-thermo on the mixed free forced convective and mass transfer boundary layer flow with temperature-dependent viscosity. Babu et al. [9] studied the diffusion-thermo and radiation effects on MHD-free convective heat and mass transfer flow past an infinite vertical plate in the presence of a chemical reaction of the first order. The dimensionless governing equations were solved with the Laplace transform technique. Rajput and Gupta [10] investigated the diffusion-thermo effect on unsteady free convection MHD flow past an exponentially accelerated plate through porous media with variable temperature and constant mass diffusion in an inclined magnetic field. Sharma and Buragohain [11] examined the Soret and Dufour effects on unsteady flow past an oscillating vertical plate with the help of numerical technique. Postelnicu [12] studied simultaneous heat and mass transfer by natural convection from a vertical plate embedded in an electrically conducting fluid saturated porous medium in the presence of Soret and Dufour effects using the Darcy–Boussinesq model. Gaikwad et al. [13] investigated the onset of double diffusive convection in two component couple of the stress fluid layer with Soret and Dufour effects using both linear and nonlinear stability analysis. Prakash et al. [14] considered the Dufour effects on unsteady MHD natural convection flow past a spontaneously started infinite vertical plate with variable temperature and constant mass diffusion through a permeable medium, and the dimensionless governing equations were solved in a closed form by utilizing the Laplace transform technique. They tracked down that the Dufour impact has critical effect on the velocity and temperature fields.Numerous fluids in practical developments show non-Newtonian behavior because the consistent Newtonian fluids do not explicitly clarify the attributes of real fluids. Among non-Newtonian fluids, second-grade fluid is one of the viscoelastic fluids which were introduced by Rivlin [15] and Rivlin and Erickson [16]. Beard and Walters [17] are considered the pioneer of viscoelastic fluids. They developed the boundary layer theory for the second-grade fluids. This boundary layer theory for the second-grade fluids has motivated many researchers to really explore this kind of fluids with various conditions. Ariel [18] attained an interpretive solution for an incompressible laminar second-grade fluid between the plates. Kecebas and Yurusoy [19] analyzed an unsteady two-dimensional power law fluid of the second grade and used a finite difference approach to solve reduced governing equations. Raftari [20] considered the MHD steady flow and heat transfer of a second-grade fluid and obtained an analytical solution. Aman et al. [21] analyzed the unsteady heat and mass transfer in second-grade fluid over a flat plate with wall suction and injection.Recently, it has progressively been seen as a dynamic tool through which a beneficial generalization of physical ideas can be obtained. Most fractional derivatives used are the Riemann–Liouville (RL) fractional derivative and the Caputo fractional derivative [22, 23]. It is observed that these operators exhibit obstacle in applications, such as the RL derivative of a constant is not zero, and the Laplace transform of the RL derivative involves terms which have no physical signification. The Caputo fractional derivative has excluded these difficulties, but the kernel of the definition is a singular function. Caputo and Fabrizio have introduced recently a new definition of the fractional derivatives with an exponential kernel without singularities [24]. The results that are been analyzed using these operators are expressed in complicated forms involving some generalized functions [25–34].The innovation of the present paper is to examine the double convection flow of an incompressible differential-type fluid near a vertical plate with heat source, Newtonian heating, and diffusion-thermo effect. Fractional derivative CF with nonsingular kernel is used in the constitutive equations of the mass flux and thermal flux to describe the diffusion and thermal processes, respectively. Semianalytical solutions of the dimensionless problems are established by virtue of the Laplace inversion numerical algorithm Gaver–Stehfest [35, 36]. Expressions of skin friction, Sherwood and Nusselt numbers with fractional, and ordinary cases, respectively, are also determined. The results which we attained here are new and can be applied to other viscoelastic fluids. Applications of this research would be helpful in magnetic material processing and chemical engineering systems. At the end, the influence of flow parameters and the fractional parameter on the temperature and concentration field as well as on the velocity field are tabularly and graphically analyzed.


## 2. Mathematical Model

Let us consider the double convection flow of an electrically conducting incompressible differential-type fluid lying over an infinite vertical plate occupying in thexξ︹1-plane with Newtonian heating as shown in Figure 1. Initially, the fluid and the plate are at rest and its temperature is M (ambient fluid temperature) and the concentration level on the plate is L−11/s2+b1s+b2=2/b12−4b2sinhb12−4b2/2te−b1/2t (ambient fluid concentration). After time t︹1=0+, the heat transfer from the plate to the fluid is proportional to the local surface temperature T︹1, and the concentration level on the plate is C︹w (wall concentration) which is thereafter kept constant. Presume that the influence of double convection and viscous dissipation in momentum and energy equations are insignificant, respectively. Also, the direction of flow has no pressure gradient. We assume that the velocity, temperature, and concentration are functions of ξ︹1 and t︹1 only. For such a flow, the constraint of incompressibility is identically satisfied. Taking the consistent Boussinesq approximation, the convection flow is governed by the following set of partial differential equations [32, 33, 37]:Figure 1
Flow geometry.(1)∂u︹1ξ︹1,t︹1∂t︹1=ν∂2u︹1ξ︹1,t︹1∂ξ︹12+α1ρ∂3u︹1ξ︹1,t︹1∂ξ︹12∂t︹1+gβTT︹1ξ︹1,t︹1−T︹∞+gβCC︹1ξ︹1,t︹1−C︹∞−σB02u︹1ξ︹1,t︹1ρ,ξ︹1,t︹1>0,(2)ρCP∂T︹1ξ︹1,t︹1∂t︹1=−∂q︹∂ξ︹1−QT︹1ξ︹1,t︹1−T︹∞−DmKTρCs∂j︹∂ξ︹1,ξ︹1,t︹1>0,(3)q︹=−k1∂T︹1∂ξ︹1,(4)∂C︹1ξ︹1,t︹1∂t︹1=−∂j︹∂ξ︹1,ξ︹1,t︹1>0,(5)j︹=−Dm∂C︹1∂ξ︹1.The appropriate initial and boundary conditions are as follows:(6)u︹1ξ︹1,0=0,T︹1ξ︹1,0=T︹∞,C︹1ξ︹1,0=C︹∞,ξ︹1≥0,(7)u︹10,t︹1=0,∂T︹1∂ξ︹1|ξ︹1=0=−h1k1T︹10,t︹1,C︹10,t︹1=C︹w,t︹1>0,(8)u︹1ξ︹1,t︹1⟶0,T︹1ξ︹1,t︹1⟶T︹∞,C︹1ξ︹1,t︹1⟶C︹∞,ξ︹1⟶0,t︹1>0.On introducing the nondimensional quantities from Appendix A into equations (1)–(8), we get the following nondimensional partial differential equations:(9)∂u∂t=∂2uy,t∂y2+γ∂3uy,t∂y2∂t+Grθy,t+GmCy,t−Muy,t,y,t>0,(10)Pr∂θ∂t=−∂q∂y−Sθy,t−PrDu∂j∂y,y,t>0,(11)q=−∂θ∂y,(12)∂C∂t=−∂j∂y,y,t>0,(13)j=−1Sc∂C∂y,with the initial and boundary conditions(14)uy,0=0,θy,0=0,Cy,0=0,y≥0,(15)u0,t=0,∂θ∂y|y=0=−θ0,t+1,C0,t=1,(16)uy,t⟶0,θy,t⟶0,Cy,t⟶0,asy⟶∞,t>0.To establish a model with time-fractional derivatives, we assume a thermal process with memory illustrated by the next generalized fractional constitutive equation for thermal flux and mass diffusion, respectively [38, 39]:(17)jy,t=−1ScDtαCF∂Cy,t∂y,α∈0,1,(18)qy,t=−DtβCF∂θy,t∂y,β∈0,1,where the CF time-fractional derivative DtpCF⋅ of order p is defined by [24](19)DtpCFuy,t=11−p∫0texp−pt−τ1−p∂uy,τ∂τdτ,p∈0,1.The Laplace transform of the CF time derivative is as follows:(20)LDtpCFuy,t=sLuy,t−uy,01−ps+p.Remark 1.
Ifuy,0=0 and p⟶0, equation (20) becomes LDtpCFuy,t=Luy,t. In this case, the generalized mass flux and thermal flux equations (17) and (18) reduce to the classical Fourier’s law and Fick’s law equations (11) and (13), respectively. For p⟶1,(21)limp⟶1LDCFtpuy,t=sLuy,t−uy,0=L∂uy,t∂t.
For the correspondence, it will be as follows:(22)∂uy,t∂t=limp⟶1DtpCFuy,t.


## 3. Solution of the Problem

### 3.1. Concentration Field

Applying the Laplace transform to equations (12), (13), third equation in (15), third equation in (16), and (17), keep in mind the initial condition (third equation in (14)), and after simplification, we obtain the transformed problem:(23)∂2C¯y,s∂y2−Sc1−αs+αC¯y,s=0,y>0,(24)C¯0,s=1s,C¯y,s⟶0,asy⟶∞.The differential equation (24) gives the solution with respect to condition (24):(25)C¯y,s=1se−ya0s+a,where a0=Sc1−αanda=a/1−α.The inverse Laplace transform of the above equation is perceived using equation (A.3) from Appendix C.(26)Cy,t=12eya0aerfcya02t+at+e−ya0aerfcya02t−at.The local mass transfer coefficient from the plate to the fluid, that is, Sherwood number, is taken by the subsequent relation:(27)Shg=j0,t=−11−αScL−1ss+a∂C¯y,s∂y|y=0=a01−αSce−atπt.The obtained results in equations (26) and (27) identical results exist in [29] (equations (3.9) and (3.22)).


### 3.2. Concentration Field for an Ordinary Caseα⟶0

In special case whenα⟶0, we obtain the ordinary concentration field by means of the equation (A.4) from Appendix C.(28)Cy,t=erfcya12t.Similarly, we obtain the expression for the Sherwood number as follows:(29)Sho=Scπt.


## 3.1. Concentration Field

Applying the Laplace transform to equations (12), (13), third equation in (15), third equation in (16), and (17), keep in mind the initial condition (third equation in (14)), and after simplification, we obtain the transformed problem:(23)∂2C¯y,s∂y2−Sc1−αs+αC¯y,s=0,y>0,(24)C¯0,s=1s,C¯y,s⟶0,asy⟶∞.The differential equation (24) gives the solution with respect to condition (24):(25)C¯y,s=1se−ya0s+a,where a0=Sc1−αanda=a/1−α.The inverse Laplace transform of the above equation is perceived using equation (A.3) from Appendix C.(26)Cy,t=12eya0aerfcya02t+at+e−ya0aerfcya02t−at.The local mass transfer coefficient from the plate to the fluid, that is, Sherwood number, is taken by the subsequent relation:(27)Shg=j0,t=−11−αScL−1ss+a∂C¯y,s∂y|y=0=a01−αSce−atπt.The obtained results in equations (26) and (27) identical results exist in [29] (equations (3.9) and (3.22)).


## 3.2. Concentration Field for an Ordinary Caseα⟶0

In special case whenα⟶0, we obtain the ordinary concentration field by means of the equation (A.4) from Appendix C.(28)Cy,t=erfcya12t.Similarly, we obtain the expression for the Sherwood number as follows:(29)Sho=Scπt.


## 4. Temperature Field

Applying the Laplace transform to equations (10), (11), second equation in (15), second equation in (16), and (18), using the initial condition (second equation in (14)), after simplification, we get the transformed problem:(30)∂2θ¯∂y2−b0s2+b1s+b2sθ¯y,s=b3s+bse−ya0s+a,(31)∂θ¯y,s∂y|y=0=−θ¯0,s+1s,θ¯y,s⟶0,asy⟶∞.Differential equation (30) gives the following solution with the boundary condition (31):(32)θ¯y,s=1s1b01w1s+ξe−yb0w1s+sb6w2se−ya0s+as+1s1b0e−yb0w1sw1s+ξ1−a0b6s+a,where bi,i=0,…,6 and ξ are constants given in Appendix B(33)b=β1−β,w1s=b0s2+b1s+b2s,w2s=s2+b4s+b5s+b.To find the inverse Laplace transform, equation (34) can be written as follows:(34)θ¯y,s=θ¯1s⋅θ¯2y,s+θ¯3s⋅θ¯4y,s+θ¯6s⋅θ¯5y,s,where(35)θ¯1s=1s1b01w1s+ξ,θ¯2y,s=e−yb0w1s,θ¯3y,s=sb6w2s,θ¯4y,s=e−ya0s+as,θ¯5s=1s1b0e−yb0w1sw1s+ξ,θ¯6y,s=1−a0b6s+a.The inverse Laplace transform of equation (35) is obtained using (A.3), (A.5)–(A.11) from Appendix C, and by taking convolution theorem, we will get the following equation:(36)θy,t=θ1t∗θ2y,t+θ3t∗θ4y,t+θ6t∗θ5y,t,where ∗ denotes the convolution product.(37)θ1t=−1b0−1b0∫0tG2τdτ,θ2y,t=e−b1te−b0y2/4tπt−∫0∞e−b1u−b0y2/4uπub2ut−uJ12b2ut−udu,θ3t=δt+2b−b4b42−4b5∫0∞δ′t−τe−b4/2τsinhb42−4b52τdτ−2b5b42−4b5sinhb42−4b52te−b4/2τ,θ4y,t=equation27,θ5y,t=θ1t∗θ2y,t=∫0tθ1t−τθ2y,τdτ,θ6t=δt+a0b6e−at2tπt,here(38)G2t=e−b1t1πt−ξeξ2terfcξt−∫0∞e−b1u1πu−ξeξ2uerfcξub2ut−uJ12b2ut−udu.The local coefficient of heat transfer from the plate to the fluid, in terms of the Nusselt number, is as follows:(39)Nug=q0,t=−11−βL−1∂θ¯y,s∂y|y=0=11−βL−1w1ss+bw1s+ξ1+b6w2s1+a0b6s+a+a0b6w2ss+a.The inverse Laplace transform of equation (39) is established numerically and is described in Section 6 in the tabular form.


### 4.1. Temperature Field for an Ordinary Caseα⟶0,β⟶0

In special case, that is,α⟶0, β⟶0, to obtain the ordinary temperature field by means of equations (A.10) and (A.11) from Appendix C,(40)θy,t=1Pr∫0t1πτe−a1τ−Pry2/4τ−ςePrςy+ς2−a1τerfcPry2τ+ςτdτ+a22∫0te−a3τ−iySca3erfcySc2τ−ia3τ+e−a3τ+iySca3erfcySc2τ+ia3τdτ,where ai,i=1,…,3 and ζ are constants given in Appendix B.The above equation would be expressed as the Nusselt number using (A.12) from Appendix C as follows:(41)Nuo=1−ςa1−ς2a1erfa1t−e−a1−ξ2terfcςt−1+a2Sc∫0∞1πuδt−u−a3e−a3t−udu.


## 4.1. Temperature Field for an Ordinary Caseα⟶0,β⟶0

In special case, that is,α⟶0, β⟶0, to obtain the ordinary temperature field by means of equations (A.10) and (A.11) from Appendix C,(40)θy,t=1Pr∫0t1πτe−a1τ−Pry2/4τ−ςePrςy+ς2−a1τerfcPry2τ+ςτdτ+a22∫0te−a3τ−iySca3erfcySc2τ−ia3τ+e−a3τ+iySca3erfcySc2τ+ia3τdτ,where ai,i=1,…,3 and ζ are constants given in Appendix B.The above equation would be expressed as the Nusselt number using (A.12) from Appendix C as follows:(41)Nuo=1−ςa1−ς2a1erfa1t−e−a1−ξ2terfcςt−1+a2Sc∫0∞1πuδt−u−a3e−a3t−udu.


## 5. Velocity Field

Applying the Laplace transform to equation (9), with the initial condition (first equation in (14)) and using the expressions of equations (25) and (32), we obtain the following Laplace transform velocity:(42)∂2u¯∂y2−s+M1+γsu¯y,s=−Gr1+γs1s1b01w1s+ξe−yb0w1s+sb6w2se−ya0s+as+1s1b0e−yb0w1sw1s+ξ1−a0b6s+a−Gm1+γs1se−ya0s+a.Here,u¯y,s represents the function of the Laplace transform and uy,t that has to satisfy the conditions.(43)u¯0,s=0,u¯y,s⟶0,asy⟶∞.The ordinary differential equation (42) gives the solution with subject to conditions (43):(44)u¯y,s=1s3+d0s2+d1s+d2.1w1s+ξA1+A2sw2s1se−yc0s+M/s+d3−e−yb0w1s+A2s3+d0s2+d1s+d2.ss+aw1s+ξ.a0w2s1se−yc0s+M/s+d3+e−yb0w1s+1s2+a4s+a5A3+A41w2s1se−yc0s+M/s+d3−e−ya0s+a,where c0=1/γai,i=4,5,Aj,j=0,…,4 and dk,k=0,…,3 are all constants given in Appendix B.The skin friction coefficient corresponding to this motion is as follows:(45)τg=∂u0,t∂y=L−1∂u¯0,s∂y=L−11s3+d0s2+d1s+d2.1w1s+ξA1+A2.sw2sc0ss+Ms+d3−b0w1s+A2s3+d0s2+d1s+d2ss+aw1s+ξ.a0w2sc0ss+Ms+d3+b0w1s+1s2+a4s+a5.A3w2s+A4c0ss+Ms+d3−a0s+a.The inverse Laplace transform of equations (44) and (45) will be found numerically in Section 6 by applying the Stehfest’s algorithm [35].


### 5.1. Velocity Field for Fractional Viscous Fluidγ⟶0

In this special case, that is,γ⟶0, we will obtain the velocity field for fractional viscous fluid from equation (44) as follows:(46)u¯2y,s=1s2+p1s+p2.1w1s+ξA5+A6sw2s1se−ys+M−e−yb0w1s+A6s2+p1s+p2.ss+aw1s+ξ.a0w2s1se−−ys+M+e−yb0w1s+1s+p3A7+A81w2s1se−ys+M−1se−ya0s+a,where Ai,i=5,…,8 and pj,j=0,…,3 are all constants given in Appendix B.The Laplace transform of equation (46) is established numerically and described in Section 6.


### 5.2. Velocity Field for Ordinary Second-Grade Fluidα⟶0,β⟶0

In this ordinary case, whereα⟶0,β⟶0, the expression of flow of velocity for the second-grade fluid given in equation (44) would be illustrated as follows:(47)u¯y,s=A9s2+p4s+p5.1s+a1+ζ.1se−yc0s+M/s+d3−1se−yPrs+a1+1s2+p4s+p5A10s+a3+A11.1se−yc0s+M/s+d3−1se−yScs,where Ai,i=9,…,11 and pj,j=4,5 are constants given in Appendix B.Equation (47) can be expressed as follows:(48)u¯y,s=u¯1s⋅u¯2y,s−u¯3y,s+u¯4s⋅u¯5y,s−u¯6y,s.The inverse Laplace transform of equation (48) using equations (A.3), (A.4), (A.7), (A.9), and (A.13) from Appendix C as well as the convolution theorem is given as follows:(49)uy,t=u1t∗u2y,t−u3y,t+u4t∗u5y,t−u6y,t,where ∗ denotes the convolution product.(50)u1t=2A9p42−4p5e−p4/2t∫0tsinhp42−4p52t−τ1πτ−ζeζ2τerfcζτep4/2−a1τdτ,u2y,t=e−yc0−yc0M−d32π∫0∞∫0t1te−d3t−yc02/4τ−τI12M−d3τtdtdτ,u3y,s=12eyPra1erfcPry2t+a1t+e−yPra1erfcPry2t−a1t,u4y,t=2A10p42−4p5e−p4/2t∫0tsinhp42−4p52t−τep4/2−a3τdτ+A112p42−4p5sinhp42−4p52te−p4/2t,u5y,t=u2y,t,u6y,t=erfcySc2t.


### 5.3. Velocity Field for Ordinary Viscous Fluidα⟶0,β⟶0,γ⟶0

In this special case whereα⟶0,β⟶0, and γ⟶0, the velocity expression for ordinary viscous given in equation (44), takes a form as follows:(51)u4y,s=A12s+p61s+a1+ζ1se−ys+M−1se−yPrs+a1+1s+p7A13s+a3+A141se−ys+M−1se−yScs,where Ai,i=12,…,14 and pj,j=6,7 are all constants given in Appendix B.Equation (51) can be expressed as follows:(52)u¯4y,s=u¯1s⋅u¯2y,s−u¯3y,s+u¯4s⋅u¯5y,s−u¯6y,s.The inverse Laplace transform of equation (52), using (A.3), (A.4), and (A.9) from Appendix C, then by the convolution theorem, is given as follows:(53)uy,t=u1t∗u2y,t−u3y,t+u4t∗u5y,t−u6y,t,where(54)u1t=A12e−p6t∫0tep6−a1τ1πτ−ζeζ2τerfcζτdτ,u2y,t=12eyp0Merfcyp02t+Mt+e−yp0Merfcyp02t−Mt,u3y,t=12eyPra1erfcyPr2t+a1t+e−yPra1erfcya22t−a1t,u4t=A13e−p7t−e−a5ta3−p7+A14e−a3t,u5y,t=u2y,t,u6y,t=erfcySc2t.


## 5.1. Velocity Field for Fractional Viscous Fluidγ⟶0

In this special case, that is,γ⟶0, we will obtain the velocity field for fractional viscous fluid from equation (44) as follows:(46)u¯2y,s=1s2+p1s+p2.1w1s+ξA5+A6sw2s1se−ys+M−e−yb0w1s+A6s2+p1s+p2.ss+aw1s+ξ.a0w2s1se−−ys+M+e−yb0w1s+1s+p3A7+A81w2s1se−ys+M−1se−ya0s+a,where Ai,i=5,…,8 and pj,j=0,…,3 are all constants given in Appendix B.The Laplace transform of equation (46) is established numerically and described in Section 6.


## 5.2. Velocity Field for Ordinary Second-Grade Fluidα⟶0,β⟶0

In this ordinary case, whereα⟶0,β⟶0, the expression of flow of velocity for the second-grade fluid given in equation (44) would be illustrated as follows:(47)u¯y,s=A9s2+p4s+p5.1s+a1+ζ.1se−yc0s+M/s+d3−1se−yPrs+a1+1s2+p4s+p5A10s+a3+A11.1se−yc0s+M/s+d3−1se−yScs,where Ai,i=9,…,11 and pj,j=4,5 are constants given in Appendix B.Equation (47) can be expressed as follows:(48)u¯y,s=u¯1s⋅u¯2y,s−u¯3y,s+u¯4s⋅u¯5y,s−u¯6y,s.The inverse Laplace transform of equation (48) using equations (A.3), (A.4), (A.7), (A.9), and (A.13) from Appendix C as well as the convolution theorem is given as follows:(49)uy,t=u1t∗u2y,t−u3y,t+u4t∗u5y,t−u6y,t,where ∗ denotes the convolution product.(50)u1t=2A9p42−4p5e−p4/2t∫0tsinhp42−4p52t−τ1πτ−ζeζ2τerfcζτep4/2−a1τdτ,u2y,t=e−yc0−yc0M−d32π∫0∞∫0t1te−d3t−yc02/4τ−τI12M−d3τtdtdτ,u3y,s=12eyPra1erfcPry2t+a1t+e−yPra1erfcPry2t−a1t,u4y,t=2A10p42−4p5e−p4/2t∫0tsinhp42−4p52t−τep4/2−a3τdτ+A112p42−4p5sinhp42−4p52te−p4/2t,u5y,t=u2y,t,u6y,t=erfcySc2t.


## 5.3. Velocity Field for Ordinary Viscous Fluidα⟶0,β⟶0,γ⟶0

In this special case whereα⟶0,β⟶0, and γ⟶0, the velocity expression for ordinary viscous given in equation (44), takes a form as follows:(51)u4y,s=A12s+p61s+a1+ζ1se−ys+M−1se−yPrs+a1+1s+p7A13s+a3+A141se−ys+M−1se−yScs,where Ai,i=12,…,14 and pj,j=6,7 are all constants given in Appendix B.Equation (51) can be expressed as follows:(52)u¯4y,s=u¯1s⋅u¯2y,s−u¯3y,s+u¯4s⋅u¯5y,s−u¯6y,s.The inverse Laplace transform of equation (52), using (A.3), (A.4), and (A.9) from Appendix C, then by the convolution theorem, is given as follows:(53)uy,t=u1t∗u2y,t−u3y,t+u4t∗u5y,t−u6y,t,where(54)u1t=A12e−p6t∫0tep6−a1τ1πτ−ζeζ2τerfcζτdτ,u2y,t=12eyp0Merfcyp02t+Mt+e−yp0Merfcyp02t−Mt,u3y,t=12eyPra1erfcyPr2t+a1t+e−yPra1erfcya22t−a1t,u4t=A13e−p7t−e−a5ta3−p7+A14e−a3t,u5y,t=u2y,t,u6y,t=erfcySc2t.


## 6. Numerical Results and Discussion

The MHD second-grade fluid on an infinite vertical plate is considered with Newtonian heating, heat source, and diffusion-thermo effects. Time-fractional derivative CF with a nonsingular kernel is used in the constitutive equations of the mass flux and thermal flux to describe the diffusion and thermal processes, respectively. The magnetic field is introduced in the fluid flow which acts as opposing force to fluid motion. The expressions for dimensionless concentration, temperature, velocity fields, skin friction, and Sherwood and Nusselt numbers are obtained by means of the Laplace transform technique. Solutions for the classical model corresponding to the integer-order derivative are also obtained as limiting cases. All the parameters and profiles which are used here are dimensionless. It is observed initially in a consequence of the fractionalize parametersαandβ on the concentration, temperature, and velocity of fluid flow. Along with the effect of the absorption parameter S, the Prandtl number Pr, Dufour number Du, and time t on the temperature as well as the impact of Schmidt number Sc, and Reynolds number Re on the concentration and fluid velocity are studied. The consequence of the thermal and mass Grashof number Gr and Gm second-grade parameter γ and magnetic parameter Mfor velocity is also presented.Figures2(a) and 3(a) present the dimensionless temperature and concentration profiles for distinct values of the fractional parameters α and β. As probable, the fluid temperature and concentration are decreasing functions with respect to their fractional parameters. Their values are maximum near the plate and smoothly decrease to zero for increasing y. Figure 4(a) was drawn to interpret the effect of the fractional parameters α and β on the fluid velocity. If we give the same values to fractional parameters, the fluid flow velocity raises by increasing the values of α and β. The influence of the Schmidt number Sc on the fluid concentration is presented in Figure 2(b). It can be clearly seen from the figure that the concentration level of the fluid decreases whenever Sc is increasing. By the increase in time t, we observed from Figure 2(c) that the concentration profile is increasing.Figure 2
Concentration profiles vs. fory at α=0.3,Sc=4,t=3.
(a)(b)(c)Figure 3
Temperature profiles vs. for y atα=0.3, β=0.3, Pr=0.6, S=3,Sc=4, Du=0.05, t=2.
(a)(b)(c)(d)(e)(f)Figure 4
Velocity profilesvs for y at α=0.3,β=0.3,Pr=0.6,S=3,Sc=4,Du=0.05Du=0.05,Gm=1,Gr=2, γ=5,M=1,t=2.
(a)(b)(c)(d)(e)(f)(g)(h)(i)(j)The influence of the Prandtl numberPrand absorption parameter S is shown in Figures 3(b) and 3(c). It is noticed that by increasing the values of PrandS, the temperature profile is decreasing while increasing corresponding to the values of the Schmidt number Sc as presented in Figure 3(d). Figures 3(e) and 4(e) illustrate the temperature and velocity profiles for different values of the Dufour number Du. It is observed that the thermal diffusivity as well as the velocity profile both are increasing and the boundary layer thickness gets maximized. In addition, by increasing in time t, the temperature profile and velocity profile both are increasing as shown in Figures 3(f) and 4(j).It is seen in Figures4(b)–4(d) that by assigning the higher values to Pr,S, Sc, and γ, the fluid velocity profile decreases. The MHD principle is used for controlling the flow field in the essential direction by changing the making of the boundary layer. The variation of the velocity profile with different values of the magnetic parameter M is shown in Figure 4(i). Increase in the values of M shows the reduction in velocity. We agree with this result as expected that the magnetic field exerts a retarding effect on the mixed convection flow. Figure 4(g) is plotted to see the impact of thermal Grashof Gr. Gr is the ratio of buoyancy forces to viscous forces on the motion of the fluid, which stimulates free or inner convection. It is found that the fluid flow velocity is increasing by increasing the values of Gr. Figure 4(f) is plotted to allocated the influence of the mass Grashof Gm. It generates due to the change in concentration by a change in the density of a fluid, and it is the ratio of buoyancy forces to viscous forces; we found that by increasing Gm, the velocity profile increases.The impact of fractional parametersα and β is described in Table 1. It is attained that the concentration, rate of mass, temperature, rate of heat transfer, velocity, and skin friction are decreasing for the large variation of fractional parameters α and β.Table 1
Validation of the model for various values of fractional parameters.
Fractional parametersα and βConcentration (α)Sherwood numberTemperature (α and β)Nusselt numberVelocity (α and β)Skin friction0.01.0003.9504.4442.45945.5831.1960.10.7531.1423.4682.414155.1881.0100.20.5670.7572.7062.307215.8240.8240.30.4260.5702.1122.153245.0010.6460.40.3210.4531.6481.967254.2740.4820.50.2410.3711.2851.761251.2460.3370.60.1810.3101.0021.548240.9030.2160.70.1350.2630.7821.341226.4900.1210.80.1010.2250.6091.001210.1000.0540.90.0750.1940.4750.829193.0640.015Furthermore, we have drawn a comparison between fractional second-grade and fractional viscous fluids ordinary fluid models in Figure5. It is investigated that ordinary fluids have higher velocities as compare to the fractional fluids. It reveals in what way noninteger-order fractional parameters influence the flow of fluid.Figure 5
Comparison between ordinary and fractional fluids.Furthermore, to see the validity of our results for concentration, temperature, and velocity profiles graphically, we plotted Figures6(a)–6(c). It can be seen from these figures that by ignoring the effects of Du andM, our results are identical to those obtained by Vieru et al. [40], Imran et al. [25], and Siddique et al. [33] for fractional fluids.Figure 6
Comparison between the present results with the concentration, temperature, and velocity profiles of Vieru et al. [40], Imran et al. [25], and Siddique et al. [33], respectively.
(a)(b)(c)


## 7. Conclusion

In this paper, we analyzed the double convective flow of an incompressible differential-type fluid near a vertical plate with heat absorption, Newtonian heating, and diffusion-thermo effect. Time-fractional derivative CF is used in the constitutive equations of the mass flux and thermal flux to describe the diffusion and thermal processes, respectively. Semianalytical solutions of the dimensionless problems are obtained by virtue of the Laplace inversion numerical algorithm Stehfest’s. The computations and discussion graphically and numerically have formed to distinguish the effect of CF time-fractional parameters and the second-grade parameterγ. From numerical simulation and graphical interpretation, the findings are summarized as follows:(i)
For greater values of fractional parameterα and flow parameter Sc, the concentration profile decreases, whereas it increases due to the increasing values of time t(ii)
For larger values of fractional parametersαandβ and flow parameters Pr and S, the temperature profile decreases, whereas it increases due to increasing values of the Dufour number Du, Sc, and time t(iii)
Fluid velocity increases with the increasing values of fractional and flow parametersα,β,Du,Gm,Gr, and time t, and it is observed that the boundary layer thickness increases and velocity is maximum near the plate(iv)
Velocity field as well as the boundary layer thickness decreases near the plate as we increase the values of flow parametersPr,S,Sc,M, and γ(v)
Ordinary fluids (Newtonian and second grade) have greater velocities than fractional fluids(vi)
Skin friction, Nusselt numbers, and Sherwood numbers decrease by increasing the fractional parametersαandβ(vii)
The solutions obtained by Vieru et al. [40], Imran et al. [25], and Siddique et al. [33] for fractional fluids are the particular case of our general results for fractional second-grade fluid when Du=M=0, and they are in good agreement graphically



---

*Source: 7625031-2021-09-11.xml*