A numerical treatment of MHD radiative flow of Micropolar nanofluid with homogeneous-heterogeneous reactions past a nonlinear stretched surface.

The impact of nonlinear thermal radiation in the flow of micropolar nanofluid past a nonlinear vertically stretching surface is investigated. The electrically conducting fluid is under the influence of magnetohydrodynamics, heat generation/absorption and mixed convection in the presence of convective boundary condition. The system of differential equations is solved numerically using the bvp4c function of MATLAB. To authenticate our results, two comparisons with already studied problems are also conducted and an excellent concurrence is found; hence reliable results are being presented. Complete deliberation for magnetite nanofluid with Ferric Oxide (Fe3O4) nanoparticles in the water-based micropolar nanofluid is also given to depict some stimulating phenomena. The effect of assorted parameters on velocity, homogeneous-heterogeneous reactions, temperature and micropolar velocity profiles are discussed and examined graphically. Moreover, graphical illustrations for the Nusselt number and Skin friction are given for sundry flow parameters. It is examined that temperature distribution and its associated boundary layer thickness increase for mounting values of the magnetic parameter. Additionally, it is detected that the Nusselt number decays when we increase the values of the Biot number.

Introduction

Flows over stretched surfaces have various engineering and industrial applications like the extrusion of plastic sheets, extraction of polymer, glass blowing, drawing of wires, paper production and rubber sheets[1]. The pioneering work of flow of Newtonian fluid over a moving rigid surface by Sakiadis[2] has motivated the follower scientists and researchers to reveal more thought-provoking ideas in this core area. Sakiadis’s idea was improved by Crane[3] who considered a stretching sheet instead of a rigid surface. This was followed by an excellent research article by Gupta and Gupta[4] who studied the effect of heat transfer on a viscid fluid past a stretching sheet. Later, Chakrabarti and Gupta[5] improved the work of Gupta and Gupta[4] by analyzing the flow of viscid fluid past a stretching sheet in the presence of magnetohydrodynamics. Devi and Ganga[6] discussed dissipation impacts on nonlinear flow with magnetohydrodynamics past a permeable medium. Considering the importance of flows over stretched surfaces, scientists and researchers are encouraged even today to look for new ideas. Some salient contributions in this regard include a study by Hsiao[7] who discussed the impacts of Ohmic dissipation and MHD on the flow of viscoelastic fluid through a stretched sheet. Ibrahim et al.[8] deliberated the magnetohydrodynamics nanofluid flow over a stretching surface. Lin et al.[9] inspected the Pseudoplastic flow of nanofluid past a time-dependent stretching surface under the influence of variable thermal conductivity, viscous dissipation, and Joule heating. The nanofluid flow over the stretching sheet, with the effects of stagnation point, heat source/sink and electro magnetohydrodynamic with slip boundary condition, was analyzed by Hsiao[10]. Hayat et al.[11] found the series solution for the third-grade fluid flow with an impact of chemical reaction past an unsteady stretching surface. Malvandi et al.[12] obtained the numerical solution of the flow of variable nanofluid in the neighborhood of a stagnation point with Navier’s slip condition. Numerical solution of flow of MHD nanofluid past a shrinking/stretching surface in a spongy medium near a stagnation point was discussed by Khalili et al.[13]. The flow of Jeffrey fluid over a linearly stretching sheet was examined by Hayat et al.[14]. Hayat et al.[15] also analyzed the Powell Eyring fluid flow past a stretched surface mounted at an angle with effects radiation and non-uniform sink/source. Ahmadi et al.[16] inspected a detailed study of nanofluid flow past an unsteady extended surface. Chen et al.[17] found the numerical solution of Maxwell fluid flow past a time-dependent stretched surface. Hua and Su[18] studied the time-dependent flow over a stretched surface in a moving fluid with impacts of frictional and Ohmic heating. Hayat et al.[19] explained the flow of viscous fluid past an exponential stretched surface with slip condition and magnetohydrodynamics. Waqas et al.[20] found numerical solution of the micropolar fluid flow with effects of mixed convection, magnetohydrodynamics and convective boundary conditions past a nonlinear stretched surface. Lu et al.[21] debated the flow of magneto-hydrodynamic Carreau nanofluid flow numerically with impacts of nonlinear thermal radiation with zero mass flux condition past a radially stretching sheet.

Flows involving effects of the chemical reactions have gained the attraction of the researchers and scientists. Applications of these flows in various processes may include the production of ceramics and polymers, chemical processing, fibrous insulations, water and air pollutions and molecular diffusion of species. Chemical reactions are labeled as homogeneous- heterogeneous (h & h) reactions. The domain in the case of heterogeneous reactions is restricted however it covers the entire phase uniformly. Combustion and catalysis cover both h & h reactions. Moreover, the heterogeneous catalyst is found in the solid phase while homogeneous catalyst exists in gaseous and liquid phases. Recent investigations involving the impacts of h & h reactions include a study by Imran et al.[22] who discussed the flow of Casson fluid numerically with effects of h & h reactions with magnetohydrodynamics and viscous dissipation. Hayat et al.[23] examined analytical solution of 2D micropolar fluid with h & h reactions and magnetohydrodynamics. The flow of 3D Oldroyd-B fluid flow with nonlinear thermal radiation and homogeneous-heterogeneous (h-h) reactions is discussed by Lu et al.[24]. Abbas and Sheikh[25] numerically found the solution of Ferrofluid with nonlinear slip condition and h & h reactions. Nadeem et al.[26] explored the influence of magnetic dipole numerically past a stretched cylinder with h & h reactions. Raju et al.[27] found dual solutions of Jeffrey fluid flow numerically with the impact of h & h reactions. Imad et al.[28] examined the Prandtl fluid flow with h & h reactions past a linearly stretched surface.

During the last few years, the heat transport phenomenon in fluid flows got broad consideration owing to its varied applications in engineering and technology. Nanofluids are comprised of nanometer-sized particles, named as nanoparticles, and some usual base fluids. The nanoparticles used in nanofluids are characteristically composed of oxides, metals, carbon nanotubes, or carbides. The concept of nanofluids was first time given by Choi[29] Later, Keblinski et al.[30] experimentally and numerically proved that the thermal conductivity of nanofluids is prominent in comparison to the ordinary base fluids. Consequently, usage of nanofluids in engineering applications like in domestic refrigerators, nuclear reactors, chiller, heat exchangers, engine cooling and vehicle thermal management has immensely increased. More about applications of nanofluids may be found in[31]. Our intention here is to report different latest studies highlighting varied aspects of nanofluids. These may include an investigation by Nadeem and Noor[32] who investigated the impact of Cattaneo-Christov mass and heat flux in a nanofluid implanted in a spongy medium. Nadeem and Lee[33] used an exponentially stretching surface to deliberate the flow of nanofluid. Sheikholeslami[34] described the ethylene glycol nanofluid flow in the presence of thermal radiation past a porous enclosure. Nadeem et al.[35] considered the influence of mass and heat transfer on peristaltic flow of a nanofluid between eccentric cylinders. Sheikholeslami[36] used complex geometry to discuss the impact of the electric field on the hydrothermal conduct of nanofluid. In the presence of porous medium Raju and Sandeep[37] investigated radiative magnetohydrodynamic Jeffrey nanofluid flow past a cone with chemical reaction.

A literature survey indicates that abundant studies are available discussing the flow of micropolar fluid flow with different geometries. Comparatively less work is done with micropolar nanofluid flow and as far as our knowledge is concerned no study so far is conducted for the flow of micropolar nanofluid with h-h reactions over a nonlinear vertically stretched sheet. Additional features of nonlinear thermal radiation and heat absorption/generation coefficient are also added characteristics towards its novelty. Numerical solution of the transformed equations is gained by utilizing the built-in function bvp4c. A complete parametric investigation is conducted using the graphical depiction of different controlling flow parameters on axial velocity, h-h reactions, temperature, microrotation velocity, and concentration profile. Additionally, the local Nusselt and friction factor numbers are graphed and dissected for increasing dynamical parameters. A comparison to a previously done study is also conducted to validate our presented results. An excellent alignment is obtained in this regard which proves that our results are indubitable.

Theory and Flow Field Analysis

We assume an incompressible, steady two-dimensional boundary layer flow of micropolar nanofluid past a nonlinear vertically stretching surface with stretching velocity u(x) = U(x) = cx along x-axis, where c is a constant. The flow is under the impacts of the nonlinear magnetic field, nonlinear thermal radiation, heat absorption/generation coefficient, and h & h reactions. The strength of the magnetic field is taken as B(x) = B0x( acting in a normal direction to y-axis. We also assume that in qusient micropolar nanofluid both base fluid (water) and nanoparticles (Fe3O4) are in thermal equilibrium and there is no chance of slip occurrence. The temperature at the surface is T where the ambient temperature is taken as T∞ as y→∞. Schematic diagram of the flow problem is depicted in Fig. 1.

Figure 1

Schematic diagram of the model.

The model for homogeneous (isothermal cubic autocatalytic) and heterogeneous reactions with two chemical species A* and B* was proposed by Merkin[38] and Chaudhary and Merkin[39,40] is given bywhere the concentration of chemical species B and A are given by b and a and k,(i = c, s) are the rate quantities. Both reaction forms are considered as isothermal. Applying the boundary layer approximation, the continuity, micropolar, energy and concentration equations can be stated as follows:These equations are supported by the following boundary conditionswhere shows the spin gradient viscosity. The symbols g1, β, κ, B(x), N, j, Q0, (D, D), (k, k), m0, k*, h, T, n are the gravitational acceleration, thermal expansion coefficient, magnetic parameter, coefficient of volumetric thermal expansion, microrotation velocity, micro inertia, the variable volumetric rate of heat source, diffusion coefficients, rate constants, boundary parameter, mean absorption coefficient, convective heat transfer coefficient, convective fluid temperature, and power law index respectively. It is pertinent to mention that the value of n directly relates to the shape of the surface. For n = 1, the surface is flat, for n > 1, the surface is inner concave and for n < 1, the surface is outer convex. Similarly, the values n can also predict the motion of the fluid. The values n = 1, n < 1, and n > 1, represent the linear motion of the fluid with constant velocity decelerated motion and accelerated motion respectively. In Eq. (6), nonlinear thermal radiation heat flux q via Rosseland’s approximation is given by

The values of the specific heat C, the density ρ, the thermal conductivity K1 and the electric conductivity σ of the base fluid (water) and the magnetite nanoparticles (Fe3O4) (water) are given in Table 1.

Table 1

Thermo-physical characteristics of the base fluid (water) and nanoparticles (Fe3O4).

Physical properties	Base fluid (water)	Fe3O4
Cp (J/kg K)	4179.00	670
ρ (kg/m3)	997.100	5180
K1 (W/mK)	0.61300	9.7
σ (Ωm)	1.4700	1163.1

The mathematical form for the thermophysical properties are

Similarity transformation

Using the following transformations

The requirement of continuity equation is inevitably satisfied, and Eqs (4) to (8) becomeand the boundary conditions in Eq. (9) yield the following formwhere prime designate the derivative with respect to η. The parameters k2, M, K, R, γ, k1, S, δ, B, λ represent the strength of heterogeneous reaction, magnetic parameter, micropolar parameter, radiation parameter, heat generation parameter, the strength of homogeneous reaction, Schmidt number, the ratio of diffusion coefficient, Biot number and mixed convection parameter respectively are defined as follows

In general, the chemical species A and B won’t be identical, nevertheless, we can expect that these will be equivalent in size[38]. Here, we suppose that the diffusion species coefficients D and D are identical, i.e., δ = 1, thus we have

Now applying the above property, Eqs (16) and (17) and their corresponding boundary condition take the form

Friction factor and local Nusselt number

The dimensional forms of the C (skin friction coefficient) and Nu (the local Nusselt number) are categorized as

Now the surface heat flux (q) and surface shear stress (τ) are assumed as:

Using Eqs (12) and (24) in Eq. (23), we acquirewhere is the local Reynolds number.

Results and Discussion

The numerical solution, of Eqs (13, 14, 15, 21) with associated boundary conditions Eqs (18) and (22), is obtained using MATLAB built-in function bvp4c. For this purpose, initially, we have converted the third and second order differential equations to the first-order ordinary differential equations using new parameters. The built-in function bvp4c requires an initial supposition and the tolerance for the current problem is given as 10−7. Considering the values of the involved parameters, we have used suitable finite values of η→∞, as η = η∞ = 3 to 6. For the computation of the numerical solution, we have selected the following fixed values of involved parameters M = 0.1, Pr = 6.2, K = 0.5, ϕ = 0.1, m0 = 0.5, B = 0.5, S = 3.

Tables 2 and 3 are erected numerically to validate our presented results by fixing values of different parameters. It is found that all obtained results are in good alignment to those given in Reddy[41] and Cortell[42].

Table 2

Comparison of −f ″(0) with[41] for numerous values of n and M when ϕ = 0.0.

n	M	−f ″(0)
		Reddy[41]	Present result
1	0.5	0.865956	0.86595
	01	1.097058	1.09705
	02	1.753714	1.75371
	03	2.744580	2.74458
2	0.5	0.950290	0.95029
	01	1.101523	1.038900
	02	1.572680	1.57268
	03	2.429416	2.42941

Table 3

Comparison of −f ″(0) with[42] for various values of n in the absence of the micropolar nanofluid and h-h reactions when M = 0 = ϕ = λ.

n	−f ″(0)
	Cortell[42]	Present result
0.0	0.627547	0.627547
0.2	0.766758	0.766758
0.5	0.889477	0.889477
0.75	0.953786	0.953786
1.0	1.0	1.0
1.5	1.061587	1.061587
3.0	1.148588	1.148588

The impact of solid volume fraction ϕ on axial velocity, microrotation and temperature distributions is depicted in Figs 2–4. Figure 2 shows that the axial velocity profile decline versus increasing values of volume fraction. The effect of ϕ on microrotation profile is given in Fig. 3. It is perceived that the microrotation profile decays near and grows far away from the boundary layer and satisfies it asymptotically. While the temperature distribution diminishes, and its associated thermal boundary layer thickness augments when values of ϕ are increased as given in Fig. 4.

Figure 2

Effect of ϕ on f′(η).

Figure 4

Effect of ϕ on θ(η).

Figure 3

Effect of ϕ on g′(η).

Figures 5–7 examine the effect of micropolar material parameter K on axial velocity, microrotation, and temperature profiles. Figure 5 demonstrates the effect of K on velocity field. It is witnessed that the velocity field and momentum boundary layer augment versus growing values of the K. This is because of the reason that viscosity declines with growing the values of K. Due to this low viscosity, the temperature distribution also decreases which is illustrated in Fig. 6. Moreover, the value of microrotation velocity at the wall diminishes continuously from maximum to zero far away from the wall as illustrated in Fig. 7.

Figure 5

Effect of K on f′(η).

Figure 7

Effect of K on g′(η).

Figure 6

Effect of K on θ(η).

Figures 8–10 outline the behavior of axial velocity, microrotation velocity, and temperature profile for different values of the magnetic parameter M respectively. Figures 8 and 9 display the impact of M on axial velocity and temperature distribution. It is witnessed that for higher the values of M, the velocity profile declines whereas temperature field and thermal boundary layer thickness increase. Actually, the increase in values of M means stronger Lorentz force which offers more resistance to the fluid’s motion. That is why decline in fluid’s velocity is witnessed. Also, because of this fact more heat is produced and as a result, the temperature profile and thermal boundary layer thickness are increased. Figure 10 shows the impact of magnetic parameter M on microrotation velocity. For gradually mounting values of magnetic parameter M, the microrotation velocity profile diminishes near the boundary and an opposite behavior is observed far-off from the boundary.

Figure 8

Effect of M on f′(η).

Figure 10

Effect of M on g′(η).

Figure 9

Effect of M on θ(η).

Figures 11–14 analyze the axial velocity, microrotation velocity, temperature and concentration fields for numerous values of power-law index parameter n respectively. Figure 11 depicts the upshot of n on velocity field. It is witnessed that the velocity distribution rises when values of n are increased. This is due to the fact that enhanced values of n are the main cause for the decrease in the thickness of the wall and as a result, more stretching of the sheet is observed which is responsible for enhanced velocity. However, the microrotation velocity and temperature profiles diminish by enhancing the values of n as given in Figs 12 and 13 respectively. Figure 14 illustrates the concentration profile for numerous values of n. It is seen that for the higher values of power-law index parameter n, the concentration field boosts.

Figure 11

Effect of n on f′(η).

Figure 14

Effect of n on h(η).

Figure 12

Effect of n on g′(η).

Figure 13

Effect of n on θ(η).

Figure 15 specifies the microrotation velocity field for distinct values of boundary parameter m0. It is noticed that the microrotation velocity profile is enhanced for growing values of m0. Figure 16 is graphed to perceive the variation in temperature profile for different values of Biot number B. It is seen that temperature profile is increased for growing estimations of Biot number. It is due to the fact that increasing values of Biot number B rises the heat transfer coefficient which ultimately enhances the temperature. The action of γ (heat generation parameter) on the temperature profile is portrayed in Fig. 17. For increasing values of heat generation parameter γ means T > T∞ i.e., more heat transfer from the surface to the fluid and as a result, the fluid temperature is enhanced. Figures 18 and 19 portrayed to depict the influence of radiation parameter R and temperature ratio parameter θ on temperature field respectively. It is apparent from the figure that the temperature field is higher for values of both the radiation parameter and temperature ratio parameter. Physically this is due to the fact that with the increase in radiation parameter, the mean absorption coefficient decreases. Hence the rate of radiative heat transfers to the fluid increases.

Figure 15

Effect of m on g′(η).

Figure 16

Effect of B on θ(η).

Figure 17

Effect of γ on θ(η).

Figure 18

Effect of R on θ(η).

Figure 19

Effect θ on θ(η).

Figure 20 specifies the concentration profile for varying values of Schmidt number S. As we know that the Schmidt number is inversely proportional to mass diffusivity. That is why concentration distribution enhances and its corresponding boundary layer thickness diminishes versus gradual increasing values of the S. In Fig. 21 impact of the strength of homogeneous reaction k1 on concentration field is described. As the reactants are depleted in a chemical reaction. Because of this fact concentration field shows decreasing the tendency for improving values of k1. From Fig. 22, it is noticed that concentration distribution is the decreasing function of the strength of heterogeneous reaction k2. Higher values of k2 weaken the diffusion coefficient and as a result, less diffused particles abate the concentration field.

Figure 20

Effect of S on h(η).

Figure 21

Effect of K1 on h(η).

Figure 22

Effect of of K2 on h(η).

Finally, the effect of physical quantities of interest such as the Skin friction coefficient and the Nusselt number for various values of dimensionless governing flow field parameters are graphed in Figs 23–26 respectively. Figure 23 portrays the impact on Skin friction coefficient for various values of Biot number B and solid volume fraction ϕ. It can be seen that the Skin friction coefficient decreases versus augmented values of the Biot number B and an opposite behavior are observed in case of solid volume fraction ϕ. Figure 24 displays effect on the Skin friction coefficient for varied estimated values of the magnetic parameter M and solid volume fraction ϕ. It is observed that the skin friction coefficient enhances for escalating values of the magnetic parameter and solid volume fraction ϕ. Figure 25 demonstrates the effect of the magnetic parameter and the volume fraction on the Nusselt number. It is found that the Nusselt number diminishes with enhancing the value of magnetic parameter however the contradictory impact is observed in the case of volume fraction parameter. Further, the Nusselt number increases with increasing the Biot number and the volume fraction parameter. This effect is displayed in Fig. 26.

Figure 23

Effects of ϕ and B on .

Figure 26

Effects of ϕ and B on .

Figure 24

Effects of ϕ and M on .

Figure 25

Effects of ϕ and M on .

Table 4 is erected to display the numerical values of Nusselt number for the varied values of the different parameter. It is seen that Nusselt number increases and Skin friction coefficient decrease for growing values of solid volume fraction ϕ and radiation parameter R. By mounting the values of heat generation parameter (γ), the Nusselt number diminishes while it enhances for increasing values of Biot number B. Moreover, the skin friction coefficient diminishes with enhancing the value of Biot number and it has no effect on skin coefficient with increasing the value of heat generation parameter. Table 5 depicts the behavior of h′(0) and φ′(0) when δ = 1. It is witnessed that for increasing values of k2, both h′(0) and φ′(0) show increasing behavior.

Table 4

Numerical values of Nusselt number , skin friction and g′(0) when .

ϕ	R d	γ	B i	Re_x^−1/2Nu_x	Re_x^1/2C_f	g′(0)
0.1	03	0.1	0.5	1.70870	1.14570	0.62173
0.2				1.73360	1.11390	0.51651
0.3				1.76080	1.06140	0.42920
0.1	0.3	0.1	0.5	0.53758	1.14640	0.62174
	01			0.84084	1.14600	0.62174
	02			1.27460	1.14580	0.62173
0.1	03	0.1	0.5	1.70870	1.14570	0.62173
		0.3		1.61730	1.14570	0.62164
		0.5		1.38010	1.14570	0.62133
0.1	03	0.1	0.5	1.70870	1.14570	0.62173
			1.0	2.89670	1.13450	0.62164
			1.5	3.74190	1.12700	0.62157

Table 5

Numerical values of h′(0) and −φ′(0) for different values of k2 when δ = 1.

k 2	h′(0)	−φ′(0)
0.2	0.16977	0.16977
0.3	0.23660	0.23660
0.4	0.29451	0.29451
0.5	0.34515	034515

Concluding Remarks

In this article, a numerical solution for the effects of h-h reactions on the flow of micropolar nanofluid with water as the base fluid and Ferric oxide as nanoparticles past a nonlinear stretched surface is obtained and discussed. Additional impacts like nonlinear thermal radiation, magnetohydrodynamics, and heat generation/absorption with convective boundary condition are added features towards the novelty of the problem. The transformed system of ordinary differential equations is solved by MATLAB built-in function bvp4 for varied values of arising parameters. The salient points of the present study are given as follows:

The velocity distribution with its associated boundary layer thickness increase for enhanced values of micropolar parameter K while the temperature profile and microrotation velocity diminish.

For higher values of k1 and k2 (strength of homogeneous and heterogeneous reactions), concentration profile reduces.

For gradually mounting values of magnetic parameter M, the microrotation velocity profile contracts near the boundary and opposite trend are witnessed far away from the boundary.

The Nusselt number and Skin friction coefficient increase and decrease for growing values of solid volume fraction ϕ and radiation parameter R respectively.