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.

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 [24].

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 [2534].

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 the x-plane with Newtonian heating as shown in Figure 1. Initially, the fluid and the plate are at rest and its temperature is (ambient fluid temperature) and the concentration level on the plate is (ambient fluid concentration). After time , the heat transfer from the plate to the fluid is proportional to the local surface temperature , and the concentration level on the plate is (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 and 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]:

The appropriate initial and boundary conditions are as follows:

On introducing the nondimensional quantities from Appendix A into equations (1)–(8), we get the following nondimensional partial differential equations: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 of order is defined by [24]

The Laplace transform of the CF time derivative is as follows:

Remark 1. If and , equation (20) becomes . 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 ,For the correspondence, it will be as follows:

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:

The differential equation (24) gives the solution with respect to condition (24):where .

The inverse Laplace transform of the above equation is perceived using equation (A.3) from Appendix C.

The local mass transfer coefficient from the plate to the fluid, that is, Sherwood number, is taken by the subsequent relation:

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

In special case when , 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:

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:

Differential equation (30) gives the following solution with the boundary condition (31):where and are constants given in Appendix B

To find the inverse Laplace transform, equation (34) can be written as follows:where

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:where denotes the convolution product.here

The local coefficient of heat transfer from the plate to the fluid, in terms of the Nusselt number, is as follows:

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

In special case, that is, , , to obtain the ordinary temperature field by means of equations (A.10) and (A.11) from Appendix C,where 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:

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:

Here, represents the function of the Laplace transform and that has to satisfy the conditions.

The ordinary differential equation (42) gives the solution with subject to conditions (43):where and are all constants given in Appendix B.

The skin friction coefficient corresponding to this motion is as follows:

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

In this special case, that is, , we will obtain the velocity field for fractional viscous fluid from equation (44) as follows:where and 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

In this ordinary case, where , the expression of flow of velocity for the second-grade fluid given in equation (44) would be illustrated as follows:where and are constants given in Appendix B.

Equation (47) can be expressed as follows:

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:where denotes the convolution product.

5.3. Velocity Field for Ordinary Viscous Fluid

In this special case where , and , the velocity expression for ordinary viscous given in equation (44), takes a form as follows:where and are all constants given in Appendix B.

Equation (51) can be expressed as follows:

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:where

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 on the concentration, temperature, and velocity of fluid flow. Along with the effect of the absorption parameter , the Prandtl number , Dufour number , and time on the temperature as well as the impact of Schmidt number , and Reynolds number on the concentration and fluid velocity are studied. The consequence of the thermal and mass Grashof number and second-grade parameter and magnetic parameter for 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. 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 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 is increasing. By the increase in time t, we observed from Figure 2(c) that the concentration profile is increasing.

The influence of the Prandtl number and absorption parameter is shown in Figures 3(b) and 3(c). It is noticed that by increasing the values of , the temperature profile is decreasing while increasing corresponding to the values of the Schmidt number 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 . 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 , 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 , , 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 is shown in Figure 4(i). Increase in the values of 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 . 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 . Figure 4(f) is plotted to allocated the influence of the mass Grashof . 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 , 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

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 , our results are identical to those obtained by Vieru et al. [40], Imran et al. [25], and Siddique et al. [33] for fractional fluids.

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 , the concentration profile decreases, whereas it increases due to the increasing values of time t(ii)For larger values of fractional parameters and flow parameters and , the temperature profile decreases, whereas it increases due to increasing values of the Dufour number , , and time (iii)Fluid velocity increases with the increasing values of fractional and flow parameters , and time , 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 parameters , 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 (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 , and they are in good agreement graphically

Appendix

A. Nomenclature

: velocity field of fluid in the x direction: coordinate axis normal to the plate: dimensionless coordinate axis normal to the plate: time: concentration of the fluid far away from the plate: temperature of the fluid near the plate: temperature of the fluid far away from the plate: specific heat at a constant pressure: mass flux diffusion: dimensionless heat absorption parameter: thermal conductivity of the fluid: coefficient of mass diffusivity: thermal Grashof number: magnetic field parameter: Schmidt number: dimensionless velocity: species concentration: concentration of the plate: dimensionless time: dimensionless concentration: dimensionless temperature: magnetic field parameter: heat flux: heat absorption parameter: heat transfer coefficient: concentration susceptibility: thermal diffusion ratio: mass Grashof number: Prandtl number: Dufour parameterGreek symbols: kinematic viscosity: acceleration due to gravity: dimensionless temperature: volumetric coefficient of thermal expansion: volumetric coefficient of expansion with concentration: dynamic viscosity: fluid density: electrical conductivity: second-grade fluid parameter: second-grade coefficient

B. Nondimensional Quantities

C. Some Constants Involved in the Text

D. Some Inverse Laplace Formulas

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.