Abstract

In this investigation, thermal radiation effect on MHD nonlinear convective micropolar couple stress nanofluid flow by non-Fourier’s-law heat flux model past a stretching sheet with the effects of diffusion-thermo, thermal-diffusion, and first-order chemical reaction rate is reported. The robust numerical method called the Galerkin finite element method is applied to solve the proposed fluid model. We applied grid-invariance test to approve the convergence of the series solution. The effect of the various pertinent variables on velocity, angular velocity, temperature, concentration, local skin friction, local wall couple stress, local Nusselt number, and local Sherwood number is analyzed in both graphical and tabular forms. The range of the major relevant parameters used for analysis of the present study was adopted from different existing literatures by taking into consideration the history of the parameters and is given by The result obtained in this study is compared with that in the available literatures to confirm the validity of the present numerical method. Our result shows that the heat and mass transfer in the flow region of micropolar couple stress fluid can be enhanced by boosting the radiation parameters.

1. Introduction

Numerous mathematical models were proposed to study the rheological properties of non-Newtonian fluids. The fluid model pioneered by Eringen [1] in 1996 revealed the existence of microscopic effects resulting from the local structure and micromotion of the fluid constituents. Moreover, they can sustain couple stresses and comprise the Newtonian models as a special case. The importance of heat and mass transfer and micropolar fluid flow is particularly evident in new and emerging areas of materials processing. Materials such as polymers, alloys, ceramics, composites, semiconductors, and optical materials need thermal energy for fabrication. For instance, temperature control helps to ensure the product quality and consistent production capacity in the polymer extrusion process. It is due to this fact that the heat and mass transfer in the boundary layer flow of non-Newtonian fluid with diverse effects of parameters attracted scholars from all corners [26]. These relevant parameters are used to control the heat and mass transfer during the extrusion process. In some studies, non-Fourier’s-law heat flux model is applied to govern the heat and mass transfer in the boundary layer flow region [7, 8].

Thermally radiative fluid flows are usually encountered when the difference between the temperature at the surface of the sheet and the ambient temperature is high. In numerous industrial processes, the thermal boundary layer thickness can be altered by the use of the thermal radiation. Examples of such industrial processes include missile technology, nuclear reactors, satellites, power plants, and gas turbines. The effect of radiation parameters in different boundary layer flow regions has been introduced by numerous scholars [914]. Mixed convection can be considered as a combination of free and forced convection which occur due to a significant difference in temperature between the surface (wall) and the ambient fluid. Mixed convection has an essential role when the buoyancy force considerably disrupts the flow and thermal fields. Ramzan et al. [15] studied a mixed convection viscoelastic nanofluids past porous media considering Soret-Dufour effect. They employed a homotopy analysis method to solve the proposed problem. Some authors [1618] analyzed the influence of mixed convection parameters on the boundary layer flow of Oldroyd-B fluid. Most recently, Ibrahim and Gadisa [19] reported the nonlinear convective flow of a couple stress-micropolar nanofluid with the effects of slip and convective boundary conditions.

On the other hand, in chemical process engineering, Dufour and Soret have a vital application. Moorthy and Senthilvadivu [20] suggested that when heat and mass transfer processes take place at the same time between the fluxes, the driving potential is of more complex nature, as energy flux can be generated not only by temperature gradients but also by composition gradients. According to the assumption of Fick’s law or Fourier’s, Soret and Dufour effects are typically dilapidated in heat and mass transfer processes. There are, however, exceptions in certain circumstances. For instance, the Soret effect can be utilized for isotope separation and in mixtures between gases with extremely light molecular weight like H2 and He. For average molecular weights like N2 and air, the Dufour effect was found to be of a significant magnitude such that it cannot be neglected [21]. Ahammad and Mollah [22] introduced the concepts of MHD free convection flow and mass transfer over a stretching sheet with Dufour and Soret effects. They solved numerically by applying Runge-Kuta with the shooting technique. Soret and Dufour effect on MHD Casson fluid past a stretching sheet was studied by Hayat et al. [23]. Later on, Ali and Shah [24] reported free-convection MHD micropolar fluid considering Soret and Dufour effects. Following this, different scholars [2529] scrutinized the effects of Soret and Dufour on micropolar fluid. Most recently, Bhubaneswar et al. [30] forwarded the concept of cross diffusion effects on MHD convection of Casson-Williamson fluid past a stretching surface with radiation and chemical reaction. Ibrahim and Gadisa [19] also analyzed the nonlinear convective flow of a couple stress-micropolar nanofluids with non-Fourier heat flux model past the stretching surface in the presence of slip and convective boundary conditions.

The present scrutiny is motivated by the results of the last two papers. As far as we revised and mentioned above, the problem of nonlinear convective flow of micropolar couple stress nanofluid using the Cattaneo-Christov model past the stretching sheet with the effects of thermal radiation and Soret, Dufour, and chemical reaction is still unnoticed. Thus, the main objective of the present study is to fill this gap. We employed the robust numerical technique called GFEM explained in equations (21)–(31). We performed the grid-independence test or grid convergence test to confirm the convergence of the series solution. The effect of the relevant parameters on linear velocity, angular velocity (microrotation), temperature, concentration, local skin friction, local wall couple stress, local Nusselt number, and local Sherwood number is elaborated in both graphical and tabular forms.

2. Problem Formulation

In this study, we consider the two-dimensional steady incompressible laminar MHD boundary layer flow of a nonlinear convective micropolar couple stress nanofluid using the Cattaneo-Christov heat flux model in the presence of thermal radiation with the effects of Soret, Dufour, and chemical reactions past the stretching sheet as plotted in Figure 1 below. The sheet is stretched linearly with velocity , where is constant. Applying these suppositions, the governing boundary layer equations with Boussinesq approximations are as follows (Wubshet and Gosa [19]):with the following boundary conditions at the surface (i.e., ). , and boundary conditions at the far field (i.e., ):where are the velocity components in the x and y directions, is the kinematic viscosity, is the couple stress viscosity, is the couple stress viscosity parameter, is the dynamic viscosity, is the vortex viscosity, is the thermal conductivity of the fluid, is the microinertia density, is the mass diffusivity, and is the spin-gradient viscosity and given by , where , is the stretching velocity, is the gravitational acceleration, and are the relaxation time and the retardation time, respectively, are the Deborah numbers with respect to relaxation time of heat flux and nanoparticles concentration, respectively, are the linear and nonlinear thermal expansion coefficients due to temperature, are the linear and nonlinear thermal expansion coefficients due to concentration, is the density of base liquid, is the concentration susceptibility, is the specific heat capacity of the base fluid, is the temperature, is the concentration, is the electric conductivity, is the magnetic parameter, is the reaction rate, is the thermal-diffusion ratio, is the mean fluid temperature, is the radiative heat flux, where is the Stefan-Boltzmann constant, and is the mean absorption coefficient.

The coupled nonlinear partial differential equations (PDEs) with boundary conditions (6) above can be reduced to the appropriate coupled nonlinear ordinary differential equations (ODEs) by using the following similarity and dimensionless variables.

is the similarity variable and the dimensionless variables are defined as follows:

The continuity equation can be satisfied if we define the stream functions as follows:

Using the introduced similarity transformation above (7), we reduced equations (1)–(5) with the boundary condition (6) to the following coupled highly nonlinear ODEs:with the corresponding boundary conditions:

Here, we affirm the dimensionless parameters as follows:where is the Prandtl number, is the Deborah number with respect to the relaxation time of the heat flux, is the magnetic field parameter, is the mixed convection parameter/local buoyancy parameter, is the nonlinear convection parameter due to temperature, is the nonlinear convection parameter due to concentration, is the ratio of concentration to thermal buoyancy forces, is the Grashof number in terms of temperature, is the Grashof number in terms of concentration, is the Schmidt number, is the Soret, is the Dufour number, is the chemical reaction rate, is the material parameter, is the couple stress parameter, is the Brownian motion parameter, and is the thermophoresis parameter.

The engineering physical quantities of interest in this paper are skin friction coefficient, wall couple stress, Nusselt number, and Sherwood number defined as follows:where the wall shear stress , wall heat flux , and wall mass flux are defined as follows:

After substituting the values introduced in the similarity transformation above, we get

3. Numerical Simulations

3.1. Finite Element Method

We applied the robust numerical technique called the Galerkin finite element method to obtain the solution of coupled nonlinear partial differential equations governing the boundary layer flow. The fundamental steps needed to apply GFEM are dividing the domain into elements/discretization of the domain, the element formulation/derivation of the element equation, the assemblage of the element equation into its global form, and the imposition of boundary conditions and solving a system of linear equations, respectively, (see [19]). The following steps are crucial to apply the finite element method.

3.1.1. Discretization of the Domain

The fundamental concept of the FEM is to divide the domain or region of the problem into small connected parts, called finite elements. The collection of elements is called the finite element mesh. These finite elements are connected in a nonoverlapping manner, such that they completely cover the entire space of the problem.

3.1.2. Generation of the Element Equations
(i)A typical element is isolated from the mesh and the variational formulation of the given problem is constructed over the typical element.(ii)Over an element, an approximate solution of the variational problem is supposed, and by substituting this in the system, the element equations are generated.(iii)The element matrix, which is also known as the stiffness matrix, is constructed by using the element interpolation functions.
3.1.3. Assembly of the Element Equations

The algebraic equations obtained are assembled by imposing the interelement continuity conditions. This yields a large number of algebraic equations known as the global finite element model, which govern the whole domain.

3.1.4. Imposition of Boundary Conditions

On the assembled equations, both Dirichlet and Neumann boundary conditions are imposed.

3.1.5. Solution of the Assembled Equations

The assembled equations so obtained can be solved by any of the numerical techniques, namely, LU decomposition method, Gauss elimination method, and so forth.

Assuming,

The system of differential equations above will be reduced to the following equations:with the following corresponding boundary conditions:

3.2. Variational Formulation

The variational formulations used in solving differential equations by the finite element method are considered in detail by different scholars [19, 3134].

The variational formulation related to the equations (15)–(19) over a typical element is given bysubjected to the boundary condition (20), where , , and are arbitrary weight functions may be regarded as variations in , respectively.

3.3. Finite Element Formulation

At the third step, we look for the approximation solution of the following form:with , the quadratic shape functions are defined aswhere

At the fourth step, substituting the approximate solution of equation (26) into equations (21)–(25), we acquire the finite element model equation which is given bywhere denotes the elemental stiffness matrix, is the vector of elemental nodal variables (unknowns), and is the force vector expressed as follow:where each is of the order and is of order . These matrices are defined aswhere

3.4. Assembly of the System of Equations

At this step, we assemble the elemental systems to obtain the global system, given as follows:

This results in a large number of algebraic equations that govern the entire domain. And then, the global matrix will be modified by imposing the boundary conditions defined in equation (20). The last step is solving the assembled system of equations by the standard techniques like the Gauss elimination method, LU decomposition, Gauss Jordan method, or any iterative scheme.

4. Results and Discussion

The main target of the present study is to analyze the effects of thermal radiation, diffusion-thermo (Dufour), thermal-diffusion (Soret), chemical reaction, and Cattaneo-Christov model on nonlinear convective MHD micropolar couple stress nanofluids past a linearly stretching surface. The robust numerical method called the Galerkin finite element method (GFEM) is applied to solve the proposed model. We performed grid-invariance test or grid convergence test to confirm the convergence of the series solution. The impact of these pertinent parameters on velocity, angular velocity, temperature, concentration, local skin friction, local wall couple stress, local Nusselt number, and local Sherwood number was analyzed in both graphical and tabular forms. The default values of the present parameters used to plot the graphs are chosen based on the existing literature and parameter history and given below (Wubshet and Gosa [19]):

Figures 2 and 3 are plotted to anticipate the control of radiation parameter on linear velocity and temperature distributions. Enlarging radiation parameter in the boundary layer flow of micropolar couple stress nanofluid is to enhance the velocity of the fluid flow as revealed in Figure 2. Radiation in the boundary layer flow region rises as increases; this in turn result enlargement in thermal boundary layer thickness as indicated in Figure 3. That is, to increase the radiation parameter is to initiate the temperature rise of the fluid flow. The Dufour parameter also produces similar effects on the velocity and temperature profiles of the laminar flow as illustrated in Figures 4 and 5. It is shown that the larger values of influenced the fluid to flow with faster speed and higher temperature. The enhancement of the caused increases in the concentration gradient which caused mass diffusion taking place more rapidly. In this circumstance, the rate of energy transfer associated with the particles became higher. That is why the temperature profile was boosted in the boundary layer flow region as plotted in Figure 5. The Dufour number has influenced insignificantly the concentration species to be lower in the laminar flow as indicated in Figure 6.

Figures 7 and 8 notice the impacts of Soret number on the concentration profile and temperature profile, respectively. For larger values of Soret parameter, the concentration contour increases significantly whereas the temperature contour decreases insignificantly as noted in Figures 7 and 8, respectively. Figures 911 illustrated the impacts of the material parameter on linear velocity, angular velocity, and temperature of the micropolar-couple stress fluid, respectively. It is revealed that increasing the material parameter is to increase the linear velocity of the fluid and lower temperature in the flow region. In this circumstance, the angular velocity has not shown consistence as plotted in Figure 10. This result is in good agreement with the study reported by Wubshet Ibrahim and Gosa Gadisa [19]. As plotted in Figures 12 and 13, the couple stress parameter has quite opposite effect on velocity and temperature of the boundary layer flow. The larger value in the couple stress parameter has the tendency to resist the fluid flow as noted in Figure 12 and caused higher temperature as shown in Figure 13. The control of mixed convection parameter on velocity, angular velocity, temperature, and concentration distributions is plotted in Figures 1417. It is revealed that, in Figure 14, the enhancement in the mixed convection parameter initiates the fluid to flow more rapidly. This is due to the fact that the higher mixed convection parameter associates with the larger thermal buoyancy force which is responsible for the improvement of the linear velocity distribution and decline in temperature and concentration distribution contours (Figures 16 and 17). However, for the increasing values of mixed convection parameter, the angular velocity profile contour has not shown consistency as indicated in Figure 15. It is observed that the microrotation distribution very close to the sheet declines and at some distant from the sheet varies quite opposite with the larger values of mixed convection parameter. Figure 18 demonstrates the influence of first-order chemical reaction rate on concentration contour. It shows that concentration diminishes on elevated values of chemical reaction parameters. Figure 19 inspects that the temperature in the boundary layer region is higher in Fourier’s model than in the Cattaneo-Christov heat flux model.

Figures 2023 are plotted to analyse the effects of radiation parameter versus mixed convection parameter on local skin friction, local wall stress, local Nusselt number, and local Sherwood number, respectively. It is concluded that, for higher values of radiation parameter in the boundary layer region of the flow, both local skin friction and Local wall couple stress decrease insignificantly, while both heat and mass transfer significantly increase. Figure 24 illustrates the effects of Dufour number versus on heat transfer rate in the laminar flow. It is observed that, with the rise of Dufour number (diffusion-thermo), the heat transfer decrease for micropolar couple stress nanofluid and Soret number (thermal-diffusion) has a decreasing impact on the mass transfer rate as indicated in Figure 25.

The grid-invariance test is performed to maintain the four-decimal point accuracy. It is also called the grid-invariance test or grid convergence test. We used this test to improve the results using a successively smaller step size for the calculations. We started by choosing a coarser mesh with 100 number of elements having step size of . Then, enhancing the number of elements ten times, we obtained a medium mesh with 1000 elements having a step size of . Finally, we have a fine mesh of 1500 elements with step size of and get four-decimal point-accuracy in velocity, angular velocity, temperature, and nanoparticle concentration values. After increasing the number of elements more than 1500, the accuracy is not affected but only to enlarge the compilation time. This is shown in Tables 14. Figures 2629 are plotted to show the course, medium, and fine meshes for every fifth element of mesh.

The assembled system of equations above is nonlinear in nature, and therefore, the iterative scheme is used to find the numerical solution. The system is linearized after incorporating the functions which are expected to be known at the beginning of the iteration. The iterative process is completed or terminated when the following condition (convergence formula) is satisfied:

Excellent convergence has been realized for all results, and it has been confirmed through the tabular and graphical forms for grid-invariance test mentioned above.

Table 5 indicates that our result is in good agreement with that in the existing literature, whereas Table 6 elaborates the effects with different parameters on local skin friction, local wall couple stress, heat, and mass transfer in the boundary layer flow region of micropolar couple stress fluid.

5. Conclusion

In this study, micropolar couple stress nanofluid flow past the stretching surface with the impact of relevant parameters is analyzed. The heat transfer in the boundary layer flow is modeled by the Cattaneo-Christov heat flux model. The robust numerical method called the Galerkin finite element method (GFEM) is applied to solve the proposed model. We performed grid-invariance test or grid convergence test to confirm the convergence of the series solution. The effect of numerous pertinent variables on velocity, angular velocity, temperature, concentration, local skin friction, local wall couple stress, local Nusselt number, and local Sherwood number is analyzed in both graphical and tabular forms, and the following remarks are forwarded:(1)Both velocity and temperature distributions are increasing functions of radiation parameter and Dufour number.(2)Chemical reaction and mixed convection parameters have a tendency to retard the concentrations of the species while Soret number revealed quite opposite effect.(3)Material parameter and couple stress parameter effects are reversed on the velocity and temperature profiles.(4)Heat and mass transfer in the flow region can be enhanced by boosting the radiation parameters.

Nomenclature

:Constant
:Dimensionless stream function
:Concentration
:Galerkin finite element method
:Specific heat at constant pressure
:Dimensionless microrotation function
:Skin friction coefficient
:Temperature
:Mass diffusivity
:Gravitational acceleration
:Mean fluid temperature
:Velocity components
:Kinematic viscosity
:Stretching velocity
:Linear and nonlinear thermal expansion coefficients due to temperature
:Linear and nonlinear thermal expansion coefficients due to concentration
:Spin gradient
:Microinertia density
:Grashof number in terms of temperature
:Grashof number in terms of concentration
:Reynolds number
:Surface heat flux
:Couple stress viscosity
:Thermal conductivity
:Density of base liquid
:Dufour number
:Thermal diffusivity of the base fluid
:Material parameter
 Sc:Soret number
:Schmidt number
:Couple stress parameter
:Nusselt number
:Vortex viscosity
:Sherwood number
:Chemical reaction term
:Condition at the free stream
:Electric conductivity
:Condition at the surface
:Nonlinear convection parameter due to concentration
:Ratio of concentration to thermal buoyancy forces
:Wall shear stress
:Radiation parameter
:Wall heat flux
:Stream function (nonlinear convection parameter due to temperature)
:Deborah number with respect to the relaxation time of the heat flux
:Concentration susceptibility
:Prandtl number
:Magnetic field parameter
:Stefan-Boltzmann constant
:Wall mass flux
:Wall couple stress
:Mixed convection parameter
:Dimensionless similarity variable.

Data Availability

The data used in this article are freely available for the user.

Conflicts of Interest

The authors declare that they have no conflicts of interest.