Abstract

Heat transfer in fluid mechanisms has a stronghold in everyday activities. To this end, nanofluids take a leading position in the advent of the betterment of thermal conductivity. The present study examines numerical investigations of incompressible magnetohydrodynamic (MHD) flow of Carreau nanofluid by considering nonlinear thermal radiation, Joule heating, temperature-dependent heat source/sink, and chemical reactions with attached Brownian movement and thermophoresis above a stretching sheet that saturates the porous medium. Pertaining similarity assumptions are used to change the flow equations into tractable forms of higher order nonlinear ordinary differential equations (ODEs). The continuation technique is adopted in the MATLAB bvp4c package for the numerical outcomes. The velocity, temperature, and nanoparticle concentration distributions in contrast to the leading parameters are availed in graphical and tabular descriptions. Among the many outcomes, increasing the radiation parameter from 0.2 to 0.8 surged the heat transfer rate by at n = 1.5 and lifted it only by at n = 0.5. By boosting the magnetic parameter from 0 to 1.5, respective and rises in local drag forces are achieved in shear-thickening and thinning regions. On top of that, chemical reactions and Brownian motion parameters decay the concentration field. The distinctiveness of this method is that a solution is secured for the problem, which is highly sensitive to initial and boundary conditions. It will be worth mentioning that these fluid flow models will be applicable in various fields, such as engineering, petroleum, nuclear safety processes, and medical science.

1. Introduction

Despite their complex stress-strain relationships, in recent years, non-Newtonian fluids have attracted researchers’ attention on account of their multidisciplinary applications (see Ref. [14]). Non-Newtonian fluids are abundant in several engineering, industrial, and biological activities, such as in food processing, in arterial blood flow, in mud-drilling machines, in polymer suspensions, and in heavy oil lubrications. The notion of improving the thermal conductivity of nanofluids by dispersing nanoparticles in conventional fluids was first initiated by Choi and Eastmann [5]. Nanofluids stand for a mixture of nanosized particles (e.g., Cu, Ti, , and normally with size ) scattered stably in host fluids (e.g., water, blood, ethylene glycol, and vegetable oil). As compared to host fluids, nanofluids possess advanced thermophysical properties. Such fluids being enhanced heat transfer agents, their applications extend to cancer and hyperthermia treatment, solar energy collectors, safety management in nuclear reactors, and cooling processes of electronics, computer chips, and transformers. However, the study of nanofluid is still not at its end, and it seems there is still a need to have an outweighing theory on how to prepare and use nanoparticles to culminate the efficiency of thermal conductivity. A nonuniform heat source/sink in the form of surface temperature and wall heat flux was considered by Ramesh et al. [6] to investigate the stagnation point flow of MHD non-Newtonian dusty fluid flow. Growth in fluid-particle interaction parameter improves the velocity field of dusty fluid and diminishes this distribution for clean fluid as discussed by the authors. Bouslimi et al. [7] used the Runge–Kutta-based shooting method to investigate heat and mass transfer properties of mixed convection nanofluid under Soret diffusion and nonlinear heating caused by an extending surface. The thermal behavior of a hybrid nanofluid in a double-pipe heat exchanger was tested by Jalili et al. [8]. The result indicated that water/ nanofluid has superior convective heat transfer than water/ nanofluid and pure water. They also found that raising nanofluid concentration from to increased the heat transfer rate by . The use of nanoparticles is also an attractive and promising field in the development of diagnostics and treatments for cardiovascular disease. The authors in [9] examined the stenosed artery heat transfer behavior of hybrid nanofluid in a slanted orientation. Stenosis is a condition that results from the narrowing of blood vessels due to the buildup of cholesterol or any other plaque that blocks the smooth flow of blood in the artery. Enhancing the nanoparticle volume fraction decreases the value of the temperature profile, as elaborated in their work. Additional studies involving the importance of nanofluids can be taken care of through the investigations in [1013].

Magnetohydrodynamic (MHD) explains the study of dynamic interactions between magnetic fields and conducting fluids. An opposing electric current is induced owing to liquid motion, which alters mechanical properties of the fluid. MHD analysis regulates many natural phenomena and engineering and industrial problems on almost every scale, such as the formation of stars, proper mixing of alloying parts, electrolysis process, the MHD generator, and magnetic filtration and separation. Ishak et al. [14] studied hydromagnetic flow and the transfer of heat energy above a vertical and expanding surface with a varying magnetic field. They noticed that strengthening magnetic intensity brings a reduction of heat energy loss and local skin friction. The accounts of heat source/sink, chemical reaction, and nonlinear thermal radiation on the magnetohydrodynamic flow of Williamson nanofluid was explored by Bouslimi et al. [15]. Their result verified that the heat transfer rate of Williamson nanofluid is adversely affected by rising values of temperature ratios, thermal radiation, and heat generation/absorption parameters.

As a result of the diversity of flow in nature, investigators have recommended various non-Newtonian models based on their constitutive rheological variables, such as Casson fluid, power-law fluid, Williamson fluid, Oldroyd-B fluid, Maxwell fluid, Carreau fluid, and Eyring–Powell fluid. Among these models are Carreau fluids, which are realized as generalized Newtonian fluids having Newtonian , shear-thickening (or dilatant with ), and shear-thinning (or pseudoplastic with ) characteristics, where n stands for the power-law index. Examples of these fluids are detergents, fluid crystals, pulps, blood, and synovial fluids. Carreau fluids have gained awareness on account of their importance in the extrusion of polymers, tumor therapy, cosmetics, capillary electrophoresis, and petrochemical industries. The impact of melting and heat source/sink has been analyzed by Khan et al. [16], taking the Carreau nanofluid over a wedge. They recognized that the rise in the above parameters deteriorates the nanoparticle concentration curve in both shear-thinning and thickening regions. The Haar wavelet quasilinearization approach was utilized by Che Ghani and Siri [17], in contemplating the MHD Carreau nanofluid model with the impact of velocity slip boundary and suction/injection above an expanding surface. Babu et al. [18] depicted the nonlinear MHD convective Carreau nanofluid flow with the effect of viscous dissipation and Arrhenius activation energy above an exponentially expanding surface. Many other studies of Carreau nanofluids are detailed in [1922].

According to Buongiorno’s [23] inspection, there exist various slip mechanisms in nanoparticle-based fluid interactions. He identified in the laminar flow region that the prominent dissipative energy sources associated with a nanoparticle-based fluid slip are Brownian motion and thermophoretic diffusion. Contemplating the Rayleigh–Benard problem, the impact of thermophoresis and Brownian movement on CuO-water nanofluid heat transfer nature was studied by Haddad et al. [24]. Their study revealed the influence of thermophoresis and Brownian motion is more amplified at low-volume fractions of nanoparticles. The influence of thermophoretic diffusion and Brownian movement on the mixed convection flow of Carreau nanofluid was studied by Irfan [25]. The importance of thermophoretic and Brownian diffusion on MHD nanofluid flow over a stretching circular cylinder with the insertion of a variable magnetic field, free stream velocity, and multiple slips was examined by Majeed et al. [26]. Thermal analysis on the effect of thermophoresis, Brownian motion, and Hall currents in micropolar nanofluids was presented by Jalili et al. [27]. Lately, the authors in [28] considered the flow behavior of micropolar nanofluid subjected to electromagnetism, thermophoresis, and Brownian motion in a rotating realm. Further investigations on Brownian diffusion and thermophoresis in light of various non-Newtonian fluids can be addressed in [2931].

Many scientists are motivated to understand heat and mass transfer phenomena associated with chemical reactions and chemical transportation processes by considering their industrial applications in combustion, catalysis, and biochemical systems. These chemically reacting systems involve homogeneous (or bulk) and heterogeneous (or surface) reactions. Representing homogeneous reactions by cubic autocatalysis and heterogeneous reactions by a first-order process, Merkin [32] studied the boundary-layer flow of homogeneous-heterogeneous reactions. The author has indicated that the surface reaction is the most influential mechanism near the front edge of a uniform stream flow over a flat surface. In [33], Khan et al. picked the MHD flow of the Powell–Eyring fluid model to elaborate on the Newtonian heating impact on homogeneous and heterogeneous reactions. They observed that the skin friction coefficient advances for large values of magnetic and fluid parameters. In contrast, the mass transfer rate dampens for a homogeneous reaction parameter. Additional furtherance on chemical reactions for various non-Newtonian fluids is reported in [3437].

High-temperature gradient MHD convective flow problems are remarkably influenced by thermal radiation. Its significance can be seen in nuclear power plants, electrical power generators, solar power collectors, glassblowing, and manufacturing of plastic and rubber sheets. To analyze the energy production feature of flow between two circular plates, Jalili et al. [38] studied the thermal radiation influence of unsteady compressing flow of magnetohydrodynamic Casson fluid. They identified that non-Newtonian fluid temperature improves by , responding to a surge in the squeezing factor value. In light of the 2D cavity and 3D cavity, the impact of magnetic field and thermal radiation on the transport process of magnetohydrodynamic convection flow was exploited by Zhang et al. [39]. The study unfolded when strong thermal radiation and weak magnetic field are set forth, and the difference in flow and thermal radiation between 2D and 3D cavities is notable as compared to weak thermal radiation and strong magnetic field combination effects. Shaw et al. [40] investigated MHD hybrid nanofluid flow exposed to quadratic and nonlinear thermal radiations. By considering radiation properties and convective phenomena, the influence of Joule heating on magnetic Carreau nanofluid has been explored in [41]. Recently, employing hybrid analytical and numerical techniques, the authors in [42] presented the heat and mass transport features of axisymmetric micropolar fluid constricted to the magnetic field in a cylindrical polar system. Their work showed that increasing the radiation parameter from 0 to changes the shape of the temperature profile while keeping the maximum and minimum temperatures unaltered.

The transport process of non-Newtonian fluid via porous media is an attractive area of research owing to its vast applications, such as packed bed reactors, geothermal industries, enhanced oil recovery, drying of paper pulp, gel chromatography, and soil structures. A pioneering work regarding the flow of fluid in a porous medium was initiated by Darcy in 1856. A mixed convection power-law fluid flow drenching the porous space was discussed by authors in [43]. The nanofluid flow in a Darcy–Forchheimer porous medium was studied by Rasool and Zhang [44], taking into account the Cattaneo–Christov heat and mass flux model. They disclosed that resistive force resulting from the porosity factor increases the temperature field. Conversely, the concentration field of nanoparticles decays in response to an increase in the inertial force. The effect of using different hybrid nanofluids as solar energy absorbers in a Darcy–Forchheimer porous medium was studied by Alzahrani et al. [45]. The account of the non-Darcy–Forchheimer law on MHD Carreau fluid flow subjected to a heat source/sink and thermal radiation above a stretching sheet was anticipated by Siddiq et al. [46].

We have surveyed the above literature and many more not listed here and found that there does not exist a study that demonstrates the combined impacts of nonlinear thermal radiation, Ohmic heating, heat source/sink, chemical reaction, and Darcy–Forchheimer inertial effect on MHD flow of Carreau nanofluid embedding linearly stretching porous surface with Brownian motion and thermophoretic events.

2. Problem Formulation

The mathematical model of the present study is developed by assuming a steady, 2D, viscous, and incompressible nonlinear radiative flow of magneto Carreau nanofluid above a stretching sheet. The sheet is expanding with the velocity , where is a stretching constant. A constant magnetic intensity is administered orthogonal to the surface. We suppose a small magnetic Reynolds number in order to detach the influence of the induced magnetic field. Furthermore, the effects of Hall current and ion slip are ignored. The porous sheet is kept at a temperature , and the fluid’s free stream temperature is with . The x-axis is aligned in the direction of extension of the sheet, while the y-axis is normal to it. The physical sketch of the flow is represented in Figure 1.

The Cauchy stress tensor for Carreau model rheology and its corresponding shear rate are premeditated as follows (see [4749]):

Here, I, p, , , , n, and denote the identity tensor, pressure, infinite shear rate viscosity, zero-shear rate viscosity, first Rivlin–Ericksen tensor, power-law index, and time material parameter, respectively.

In two-dimensional steady flow, the velocity, temperature, and concentration fields take the following form:where and represent the respective x-and y-direction velocities.

Depending on the above assumptions, we execute the boundary-layer analysis that leads us to the set of equations governing the conservation of continuity, momentum, energy, and nanoparticle concentration expressed, respectively, as follows:

Here, denotes the effective heat capacity ratio of nanoparticles to base fluid. All the other notations are indicated in nomenclature. As in [17], the boundary conditions suitable for the current problem are as follows:where are the mass transfer velocity, heat transfer coefficient, and wall temperature of the fluid, respectively.

For an optically thick medium, the Rosseland approximation for the radiative heat flux provideswhere is the coefficient of the mean absorption and is the Stefan–Boltzmann constant. For nonlinear radiation,

We adopt the following similarity transformations:where is the similarity variable and denotes the stream function satisfying the equation of continuity , .

Using equations (1)–(3), (11), and (12), the transformed nondimensional form of governing equations (5)–(8) and corresponding boundary conditions (9) are expressed as follows:

With the boundary conditions,where is the Weissenberg number, is the Prandtl number, is the porosity parameter, is the local inertia coefficient parameter, is the Eckert number, is the magnetic parameter, is the thermal radiation parameter, is the ratio temperature parameter, is the Brownian motion parameter, is the thermophoresis parameter, is the heat source/sink parameter, is the Schmidt number, is the chemical reaction parameter, and is the Biot number.

3. Numerical Approaches

Bvp4c is a residual error-based mesh adaptive finite difference program that executes the three-stage Lobatto IIIa formula. The conventional method in Bvp4c provides a guess to the missed initial conditions. However, we found our problem very sensitive to the initial conditions and boundary values. Therefore, we adopted the continuation technique to resolve this problem. It is a tactic developed by Robert and Shipman [50]. They used the method to solve boundary value problems that cannot be addressed by conventional shooting methods. Introducing variable for , we reduce the order of equations (13)–(15) into a system of 7 ODEs:

For numerical computation, we have used the tolerance error, . Next, to implement the continuation method, the following steps are taken in the Bvp4c MATLAB program:(1)Using a dummy parameter , write the system as a sum of its linear and nonlinear parts(2)Using normal guess and transformed boundary conditions (19), approximate the solution for (18) by the linear part(3)Employing computed values of step 2 and taking a small positive fraction of as a coefficient of nonlinear part, approximate the solution(4)Repeat step 3 with small increment in until we come up with approximate solution to the problem with in the prescribed error tolerance

To examine the surface drag force, heat, and mass transfer rate, we have determined the skin friction, Nusselt number, and Sherwood number for Carreau fluid as follows:

In a dimensionless form, the local Reynolds number is ; the local skin friction, Nusselt, and Sherwood number, respectively, take the following form:

4. Results and Discussion

The transformed first-order ODE (18), along with the corresponding boundary conditions (19), has been solved using the continuation technique implemented on MATLAB package bvp4c. The influence of leading parameters on dimensionless fluid velocity, temperature, and nanoparticle concentration is discussed for both shearing cases in graphical and tabular forms. For the accuracy and verification of the method, a comparison in local shear stress has been made for the present and published results in [46] under some restricted conditions.

The Lorentz force effect resulting from the presence of applied magnetic field is noticed in Figures 24. This force regulates the way of fluid flow. Figure 2 indicates for both dilatant and pseudoplastic fluids, an increase in the magnetic parameter value deteriorates the velocity field. Here, the velocity curve downturn in the latter fluid type is higher than in the former one. Across the boundary layer, a sharp surge in the temperature field is achieved in response to increasing values of M. A similar behavior is noticed in the nanoparticle concentration at the free stream, and it follows the reverse direction near the wall (see Figures 3 and 4). These results coincide with the physical meaning of increasing magnetic field that intensifies the Lorentz force, which in turn has the capability of increasing the nanoparticle volume fraction in the motion of nanofluid and fluid temperature as well. Figures 57 illustrate the Weissenberg number effect. Practically, measures the time taken by the fluid to relax before regaining its original shape. An increase in amplifies the gap between shear-thinning and thickening regions of velocity, temperature, and concentration fields. Enhancing stimulates velocity in shear-thickening cases and declines it in shear-thinning regions. The opposite scenario is demonstrated in temperature and concentration fields.

The local inertia coefficient and porosity parameter influence on the momentum boundary layer are realized in Figures 8 and 9, respectively. For both regions, strengthening these parameters discouraged momentum thickness.

Temperature distribution responses to porosity, thermophoresis, and temperature ratio parameter values are demonstrated in Figures 1012. The increment in these parameters augmented the temperature field and its corresponding thermal thickness. The result in Figure 11 agrees with the fact that an increase in the thermophoresis parameter results in extra particles pushing away from the stretchable wall. As a result, the temperature distribution is facilitated. In Figure 12, the improvement of initiates the conductivity of the fluid flow. Consequently, uplifts. However, as noticed in the figures, shear-thinning fluids are more influenced than shear-thickening fluids for the mentioned parameters. The behavior of the temperature field corresponding to variations in the Eckert number is captured in Figure 13. The rise in Ec imparts excess viscous heating to the fluid. Following this, kinetic energy is transformed into internal energy, which in turn pushes up the temperature curve. Figures 1416 display the impacts of other pertinent parameters on temperature boundary layer thickness. Figure 14 describes that decreases with a rise in the value of the Prandtl number in the two fluid regions. Figure 15 shows us that a significant rise in temperature distribution can be achieved by a small increment in the Biot number, . The effect of thermal radiation is analyzed in Figure 16. Here, the temperature distribution and its corresponding thermal thickness are escalated by the rise of Rd. Substantially, the growth in radiation ray transports more heat energy to the running fluid, and in turn, a rise in thermal boundary layer thickness ensues. Here, we bear in mind that thermal boundary layer thickness is higher in the case of shear thinning than in shear-thickening behavior. The temperature field response for variations in the heat source/sink parameter is illustrated in Figure 17. The rise in values of causes temperature distribution growth.

Figures 1820 portray the twofold influence of parameter thermophoresis (Nt), Brownian motion (Nb), and the Schmidt number (Sc) in nanoparticle concentration. Near the wall and in the free stream, their variation brings reversed effects. Figure 18 displays a growing behavior of the concentration field for the rise in Nt. Amplifying the thermophoretic parameter causes the microscopic transfer of nanoparticles from warmer to cooler regions, which in turn grows the nanoparticle concentration. An increase in the Brownian motion parameter declines the nanoparticle concentration field, as witnessed in Figure 19. Boundary layer diffusion of nanoparticles is strengthened as the value in Brownian parameter surges. Consequently, nanoparticle concentration downturns. Figure 20 shows the expected response of nanoparticle concentration for an increase in the Schmidt number value. The coefficient of mass diffusion is inversely related to the Schmidt number. As a result, a decline in the distribution of nanoparticle concentration and its corresponding solutal thickness is realized. Figure 21 presents the nanoparticle concentration distribution to deliberate the influence of variations of the chemical reaction parameter. It reveals a decaying nature in the nanoparticle concentration field for the rise in .

Table 1 compares the shear stress values of the present result with the work of Siddiq et al. [46], under certain restricted conditions. It can be concluded that an excellent agreement is achieved. From Table 1, we notice that increases in parameters M, p, and inflated wall shear stress. On the contrary, an increase in n and deflated it.

Table 2 executes the local heat transfer rate against various values of listed parameters. Except for the Eckert number, the surge in the heat transfer rate occurred with the rise in parameters Pr, , , and for both shear-thinning and thickening behaviors. Here, we perceive that the heat transfer rate in the shear-thickening region is better than in the shear-thinning region.

Table 3 demonstrates the mass transfer rate response for various values of parameters Nb, Nt, , Sc, and pr. For both behaviors, the rise in Nt and pr encouraged the mass transfer rate, as opposed to the parameters Nb, , and Sc.

5. Conclusion and Future Directions

In the present work, we have numerically analyzed the boundary layer incompressible MHD Carreau nanofluid flow above an extending sheet under the effect of nonlinear thermal radiation, Joule heating, heat absorption/generation, chemical reaction, and Darcy–Forchheimer law with the incorporation of thermophoresis and Brownian motion, issuing shear-thinning and thickening behavior. The continuation technique in the MATLAB bvp4c algorithm has been utilized to obtain numerical results. The findings of this study are summarized as follows:(1)Increasing M strengthens thermal and concentration boundary layers, while in contrast, the momentum boundary layer is decayed.(2)As the Weissenberg number scales up, the velocity field for shear-thickening, the temperature and concentration fields for the shear-thinning case advance. On the other hand, the rise in this parameter results in a reverse effect on the velocity curve for the shear-thinning case and on the temperature and concentration profiles for the shear-thickening case.(3)The rise in porosity and local inertia coefficient parameter penalizes the velocity field.(4)An increment in from 0.1 to 0.4 lifted wall temperature by at and by at .(5)Increasing the values of p, , , and uplifts temperature distribution. On the contrary, an increase in Pr diminishes this profile.(6)The concentration field is decayed by the rising value on Nb, Sc, and , whereas a rise in Nt encouraged this profile.(7)It is witnessed that increasing Ec from 0 to 1.5, the average rise in temperature is for shear-thickening fluid and for shear-thinning fluid.(8)By increasing thermophoresis from 0.1 to 0.4, the average rise in boundary layer temperature is at and at .(9)When the values of M, p, and are uplifted, wall friction aggravates. The opposite result is achieved by surging the values of n and .

In our current work, we have employed the Fourier heat flux model and ignored the induced magnetic field. Nonetheless, in numerous empirical endeavors, such as high-temperature plasma, power generation, and purification of crude oil, the induced magnetic field effect plays a vital role. In future investigation, this problem can be extended for hybrid and trihybrid Carreau nanofluid convection and bioconvection by considering the induced magnetic field as well as the Cattaneo–Christov heat flux model.

Nomenclature

, :Velocity components
, :Space coordinates
a:Stretching constant
n:Power-law index
:Magnetic field
:Nanoparticle concentration
:Fluid temperature
:Wall temperature of the fluid
:Brownian diffusion coefficient
:Uniform volumetric heat source/sink
:Brownian motion parameter
:Free stream temperature
:Forchheimer coefficient
:Permeability of porous medium
:Thermophoresis parameter
:Local Reynolds number
:Radiation parameter
:Thermophoresis diffusion coefficient
:Temperature ratio
:Local Nusselt number
:Schmidt number
k:Thermal conductivity
Pr:Prandtl number
:Weissenberg number
:Eckert number
:Radiative heat flux
:Local Sherwood number.
Greek Letters
:Kinematic viscosity
:Local inertia coefficient
:Similarity variable
:Time constant
:Heat capacity ratio
:Fluid density
:Porosity of the medium
:Chemical reaction parameter
:Heat source/sink parameter
:Stream function
:Heat transfer Biot number
:Stefan–Boltzmann constant
:Dimensionless temperature
:Dimensionless concentration
:Electrical conductivity.
Abbreviations
ODE:Ordinary differential equation
MHD:Magnetohydrodynamics
Eq:Equation
Ref:Reference
2D, 3D:Two dimension, three dimension
Pr:Present
e.g.:Example.

Data Availability

The data used to support the findings of this investigation are incorporated in this article.

Conflicts of Interest

The authors declare that they have no conflicts of interest with regard to the publication of this article.

Acknowledgments

We express our heartfelt thankfulness to the authors of the literature used in this paper for the provided scientific aspects and ideas that support this work.