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

*is 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. *e constitutive equations for the mass flux and thermal flux are modeled for noninteger-order derivative Caputo–Fabrizio (CF) with nonsingular kernel, respectively. *e 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. *is shows how the fractional parameters affect the fluid flow.


Introduction
e 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. e three mechanisms for heat transfer are diffusion, convection, and radiation. e 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. e 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][3][4].
Electrically conducting fluids also have accepted enough consideration from the researchers due to their extensive applications in industrial appliances. e 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. ey showed that the influence of the magnetic parameter on streamlines and isotherms is significant. e 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.
e 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. e 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. ey 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. ey developed the boundary layer theory for the second-grade fluids. is 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 secondgrade 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 secondgrade 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. e 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]. e results that are been analyzed using these operators are expressed in complicated forms involving some generalized functions [25][26][27][28][29][30][31][32][33][34]. e innovation of the present paper is to examine the double convection flow of an incompressible differentialtype 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  Figure 1: Flow geometry. 2 Complexity 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. e 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.
On introducing the nondimensional quantities from Appendix A into equations (1)-(8), we get the following nondimensional partial differential equations: Pr with the initial and boundary conditions 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]: where the CF time-fractional derivative CF D p t (·) of order p is defined by [24] e Laplace transform of the CF time derivative is as follows: 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, For the correspondence, it will be as follows:
e local mass transfer coefficient from the plate to the fluid, that is, Sherwood number, is taken by the subsequent relation: e obtained results in equations (26) and (27) identical results exist in [29] (equations (3.9) and (3.22)).

Concentration Field for an Ordinary Case
In special case when α ⟶ 0, we obtain the ordinary concentration field by means of the equation (A.4) from Appendix C.
Similarly, we obtain the expression for the Sherwood number as follows:

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: Differential equation (30) gives the following solution with the boundary condition (31): where b i , i � 0, . . . , 6 and ξ are constants given in Appendix B To find the inverse Laplace transform, equation (34) can be written as follows: where e 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: where * denotes the convolution product. (27), here 6 Complexity e local coefficient of heat transfer from the plate to the fluid, in terms of the Nusselt number, is as follows: e inverse Laplace transform of equation (39) is established numerically and is described in Section 6 in the tabular form.
e above equation would be expressed as the Nusselt number using (A.12) from Appendix C as follows:

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: Here, u(y, s) represents the function of the Laplace transform and u(y, t) that has to satisfy the conditions. u(0, s) � 0, u(y, s) ⟶ 0, as y ⟶ ∞.
Complexity 7 e ordinary differential equation (42) gives the solution with subject to conditions (43): e skin friction coefficient corresponding to this motion is as follows: e inverse Laplace transform of equations (44) and (45) will be found numerically in Section 6 by applying the Stehfest's algorithm [35].

Velocity Field for Fractional Viscous Fluid (c ⟶ 0).
In this special case, that is, c ⟶ 0, we will obtain the velocity field for fractional viscous fluid from equation (44) as follows: where A i , i � 5, . . . , 8 and p j , j � 0, . . . , 3 are all constants given in Appendix B.
e Laplace transform of equation (46) is established numerically and described in Section 6.

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: 8 Complexity where A i , i � 9, . . . , 11 and p j , j � 4, 5 are constants given in Appendix B.

Numerical Results and Discussion
e 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. e magnetic field is introduced in the fluid flow which acts as opposing force to fluid motion. e 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. e consequence of the thermal and mass Grashof number Gr and Gm second-grade parameter c and magnetic parameter Mfor velocity is also presented. Figures 2(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. eir 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 β. e 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. 10 Complexity e influence of the Prandtl number Prand absorption parameter S is shown in Figures 3(b) and 3(c). It is noticed that by increasing the values of Pr and S, 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 Figures 4(b)-4(d) that by assigning the higher values to Pr, S, Sc, and c, the fluid velocity profile decreases. e MHD principle is used for controlling the flow field in the essential direction by changing the making of the boundary layer. e 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. Complexity 13 e 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 β.
Furthermore, we have drawn a comparison between fractional second-grade and fractional viscous fluids ordinary fluid models in Figure 5. 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.
Furthermore, to see the validity of our results for concentration, temperature, and velocity profiles graphically, we plotted Figures 6(a)-6(c). It can be seen from these figures that by ignoring the effects of Du and M, our results are identical to those obtained by Vieru et al. [40], Imran et al. [25], and Siddique et al. [33] for fractional fluids.

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 diffusionthermo 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. e computations and discussion graphically and numerically have formed to distinguish the effect of CF time-fractional parameters and the second-grade parameter c. From numerical simulation and graphical interpretation, the findings are summarized as follows:  (c) 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.

Data Availability
No data were used to support this study.

Conflicts of Interest
e authors declare that they have no conflicts of interest.