Analyzing the impact of induced magnetic flux and Fourier's and Fick's theories on the Carreau-Yasuda nanofluid flow.

The current study analyzes the effects of modified Fourier and Fick's theories on the Carreau-Yasuda nanofluid flow over a stretched surface accompanying activation energy with binary chemical reaction. Mechanism of heat transfer is observed in the occurrence of heat source/sink and Newtonian heating. The induced magnetic field is incorporated to boost the electric conductivity of nanofluid. The formulation of the model consists of nonlinear coupled partial differential equations that are transmuted into coupled ordinary differential equations with high nonlinearity by applying boundary layer approximation. The numerical solution of this coupled system is carried out by implementing the MATLAB solver bvp4c package. Also, to verify the accuracy of the numerical scheme grid-free analysis for the Nusselt number is presented. The influence of different parameters, for example, reciprocal magnetic Prandtl number, stretching ratio parameter, Brownian motion, thermophoresis, and Schmidt number on the physical quantities like velocity, temperature distribution, and concentration distribution are addressed with graphs. The Skin friction coefficient and local Nusselt number for different parameters are estimated through Tables. The analysis shows that the concentration of nanoparticles increases on increasing the chemical reaction with activation energy and also Brownian motion efficiency and thermophoresis parameter increases the nanoparticle concentration. Opposite behavior of velocity profile and the Skin friction coefficient is observed for increasing the stretching ratio parameter. In order to validate the present results, a comparison with previously published results is presented. Also, Factors of thermal and solutal relaxation time effectively contribute to optimizing the process of stretchable surface chilling, which is important in many industrial applications.

Introduction

Most commonly working liquids and materials, in many engineering disciplines, such as material and chemical processing, possess multifaceted rheological properties, whose viscosity and viscoelasticity can be continuously deformed and reshaped by imposing some forces and external conditions, such as temperature, timescale, stress, and strain. These fluids have the property to be used as heat exchangers and coolant to reduce pumping power. Other than the Newtonian model, these fluids also exhibit a relationship between shear-stress—a strain that makes them completely different. Nowadays, non-Newtonian fluid dynamics are involved in abundant researches due to their practical application. Such liquids have a shear dependent viscosity. The Carreau-Yasuda model is one of the non-Newtonian liquid models which predicts the shear-thinning/thickening behavior and also shows both elevated and low shear levels relationship. Because of this fact, the Carreau-Yasuda model has great industrial and technological applications which include drilling muds, molten polymers, oils, foodstuff, volcanic lava, liquid suspensions, certain paints, polycrystal melts, cosmetic products, and many more. Tanveer et al.[1]] discussed the Carreau-Yasuda nanofluid model in a curved channel with porous space and mixed convection. Carreau-Yasuda nanofluid flow over a nonlinearly stretching sheet in the presence of Joule heating is addressed by Shahid et al.[2]. Khechiba et al.[3] studied the effect of Carreau-Yasuda rheological parameters in a horizontal porous cavity on subcritical Lapwood convection. Abbasi et al.[4] explored entropy generation for Carreau-Yasuda nanofluid flow. Hayat et al.[5] discussed in detail the Carreau-Yasuda fluid flow with Soret and Dufour effects theoretically and numerically. Ahmed et al.[6] investigated entropy generation to analyze the peristaltic motion of Carreau-Yasuda fluid. Different other features of Carreau-Yasuda nanofluid are elaborated by different researchers[7-9].

A Magnetic field not only exhibits significant features to regulate the cooling rate but also gives high-quality industrial outputs. Magnetohydrodynamics (MHD) is particularly concerned with the dynamics of electrically conducting fluids. MHD has gained a wide range of attention from many scientists[10-15] in the field of medicine and science. Ali et al.[16] highlighted the numerical investigation of the induced magnetic field in mixed convection in the stagnation flow point. They explained when an induced magnetic field is enhanced the transfer rate of heat as well as Skin friction is enhanced. Ramzan et al.[17] illustrated MHD Carreau flow with nonlinear radiation effects with zero mass surface flux. The significance of MHD and heat sink/source has been observed in many engineering and physical processes. Salahuddin et al.[18] used a sensor surface in MHD squeezed flow to understand Carreau-Yasuda fluid flow. Hayat et al.[19] described the radial magnetic field effects on Carreau-Yasuda nanofluid flow using a curved geometry. Tanveer et al.[20] addressed the influence of the radial magnetic field with velocity and thermal slip effects on the peristaltic flow of the Carreau-Yasuda fluid. They observed the opposite impacts of velocity slip on temperature and velocity distributions. Hayat et al.[21] explained joule heating and MHD effects on an inclined channel. Malik et al.[22] elaborated on the Darcy–Forchheimer Carreau-Yasuda nanofluid flow with magnetohydrodynamics. Abdul Hakeem et al.[23] used a porous medium to explained MHD effects with the heat source and sink term. Hayat et al.[24] addressed how MHD affects Carreau-Yasuda fluid peristaltic transport in the presence of curved channel and Hall effects. The heat sink and heat source effects considering MHD flow over the vertical stretched sheet was studied by Ibrahim et al.[25]. They examined that the velocity profile and the transfer rate of heat are increased by enhancing the Hartmann number and parameter of mixed convection. Bilal et al.[26] discussed the exclusive impact of a magnetized viscous fluid in the presence of Hall current over a variably thicked sheet and non-Fourier flux theory. Khan et al.[27] elaborated the bioconvection MHD Carreau nanofluid flow for the parabolic surface. He used generalized Fourier’s and Fick’s laws to discover the mass and heat flux phenomena.

In fluid mechanics, different researchers[28-32] put remarkable contributions to understand the mass transfer process by adding Arrhenius activation energy and chemical reaction. It has a large variety of applications in recovery of thermal oil, chemical engineering, in geothermal reservoirs, and cooling of nuclear reactors. Alghamdi[33] discussed the significance of binary chemical reaction on a nanofluid flow in the presence of Arrhenius Activation Energy with a rotating disk geometry with mixed convection. Dhlamini et al.[34] studied nanofluid flow with Arrhenius Activation Energy and binary chemical reaction in presence of mixed convection and convective boundary conditions. The meaning of Arrhenius Activation Energy and binary chemical reaction in the presence of heat source/sink was elaborated by Hayat et al.[35]. Analysis of Arrhenius Activation Energy and binary chemical reaction in Couette-Poiseuille nanofluid flow is reported by Ellahi et al.[36]. Hayat et al.[37] explained binary chemical reaction and Arrhenius activation energy in MHD nanofluid flow with entropy generation minimization. Some recent investigations featuring aspects of activation energy may be found in[38-45].

In light of the above literature review, this current study is focused to scrutinize the detailed aspects of induced magnetic flux and modified Fourier and Fick's theories on Carreau-Yasuda nanofluid flow with binary chemical reaction and activation energy. Besides, the Newtonian heating condition is incorporated on the boundary to analyze the behavior of the flow. The impact of different physical parameters is also analyzed. Furthermore, the effects of the local Nusselt number and the Skin friction coefficient are addressed. The transport phenomenon is explained by governing equations that are developed including effects of heat generation/absorption in the energy equation. The complicated nonlinear equations are solved using MATLAB solver bvp4c. Graphical results are displayed to discuss the behavior of involved parameters.

Mathematical formulation

The proposed mathematical model is considered under following assumptions:

Here, steady Carreau-Yasuda nanofluid's incompressible flow to a stretching sheet is considered.

This current investigation is obtained by assuming the flow with chemical reaction and the Arrhenius activation energy process.

To analyze the heat mechanism heat generation/absorption is introduced.

The stretched surface is located on - axis and Carreau-Yasuda nano liquid is occupied in the region .

Let and . The “” and “” correspond to stretching and free stream velocities respectively.

An induced magnetic field with its parallel component and normal component is considered here. In Fig. 1 configuration of the flow is plotted.

Figure 1

Schematic diagram of the model.

The Carreau-Yasuda nanofluid model[2,18,19] is explained by the following equation:where and are the parameters used for Carreau-Yasuda nanofluid, shows first Rivlin-Ericksen tensor, is represented by here Assuming transform the above equation into the following form:

Under the above assumptions the model equations are written as[2,18,19,46]:with the following set of conditions on the boundary

Here Eqs. (3) and (4) are the continuity and its corresponding Induced magnetic field equations. Similaraly, Eqs. (5) and (6) represent the momentum and its associated Induced magnetic field equation. However, the heat and concentration equations are numbered as Eqs. (7) and (8), respectively.

Applying transformation

The model equations are transfigured under the aforementioned transformation into the following form:

With

The non-dimensional variables are defined as:

Physical quantities

Physically interesting quantities are to be determined for practical applications in engineering. For instance, the rate of heat transfer and the coefficient of Skin friction is mathematically expressed as:in the above equation, the shear stress is and the heat flux is defined as:

The final expression for the and is written as:in which local Reynold number is

Numerical scheme

To test the translated coupled non-linear ordinary differential equations, MATLAB program bvp4c is implemented. Using bvp4c, we calculate residuals and carry out the computations for different step sizes As in this process, the absolute convergence requirements were taken . It is most critical that the necessary finite values of being chosen. For this computational purpose, the asymptotic boundary conditions at for a given case is constrained to that is enough to better illustrate the solution's asymptotic behavior, of the governed equations. The assumed initial approximation must encounter the BCs without disturbing the solution process. For this first and foremost, fresh variables are added as:

Using the above expressions in MATLAB bvp4c, a new form of first-order equations is:with the transformed BCs

Graphical discussion

A graphical interpretation is captured in this section for different embedded physical parameters like and [2,18,19,46] on the physical quantities like velocity, temperature, and concentration. The desired ranges of the these parameters .

Velocity and Induced magnetic field profile

The velocity profile is analyzed in Fig. 2 for different values of . Increasing trends in the velocity are examined for the mounting values of . It is because of an increase in the value of exhibit more pressure initially, which ultimately increases the velocity and momentum boundary layer. is examined for different values of in Fig. 3. Enhanced behavior of momentum is noticed on increasing induced magnetic field parameters. Generally, an electric current is developed on an increasing magnetic field which causes the electric force to increase, which, ultimately responsible for the increase in the thickness of both thermal as well as momentum boundary layers. Figure 4 elucidates that the increasing value of increases . Since Weissenberg is measured as a ratio of the liquid's relaxation time to a given process time. So, its increasing value produces enhancement in relaxation time of the fluid which offers more pressure in flow direction and causes enhancement in the velocity. Figure 5 increases for mounting which causes more resistance in fluid flow. As a result, axial velocity flow declines upon increasing . The note represents Newtonian flow behavior. Also, the effects are highlighted on in this section. Clearly, Fig. 6 shows that for increasing the first increases near the surface from (0 to 1.9) and then decrease from (2 to 5). Induced magnetic boundary layer thickness also increases but decreases far from the boundary as the magnetic field is responsible for producing Lorentz force in the flow direction which increases for higher .

Figure 2

versus .

Figure 3

versus .

Figure 4

versus .

Figure 5

versus .

Figure 6

versus .

Temperature profile

The effects of on are explained in this section. The influence of on for heat sink and the heat source is plotted in Figs. 7 and 8. An increase in the temperature of nanofluid is measured for enhancing the value of the heat sink . Whereas, enhancing the value of the heat source reduces the fluid temperature. Sink term puts more energy into layers of a thermal boundary, which leads to a rise in . Figure 9 shows when increases the thermal distribution also increases. Enhancing thermal relaxation time heat transfer from one particle to another is fast which becomes the reason for the increase in temperature. Figure 10 shows the effects of . Temperature distribution boosts up due to the random motion of particles for growing values of . The response of temperature distribution against the thermophoresis parameter is seen in Fig. 11. Physically when rises thermophoretic force increases due to which hot particles move toward cold particles hence temperature rises. Increasing trends of temperature distribution are observed for increasing values of in Fig. 12. Heat transfer and thermal boundary layer have a direct relation with . Therefore, when is increased heat transfer coefficient is increased, and ultimately heat is transferred from the heated source to the cooler surface, which transfers extra heat from the surface to the nanofluid. Hence, of nanofluid increases on increasing . Similarly, Fig. 13 highlights the effects of . The temperature profile is enhanced for boosting the in the flow region. This temperature profile variation is because of large values of considerably reduces the thermal diffusivity and as a result, the thermal boundary layer thickness reduces which produces an increment in thermal profile.

Figure 7

versus .

Figure 8

versus .

Figure 9

versus .

Figure 10

versus .

Figure 11

versus .

Figure 12

versus .

Figure 13

versus .

Concentration profile

The importance of activation energy is displayed in Fig. 14. The figure shows that the enhancing , excite thickness of the boundary layer increases the concentration . Figure 15 highlights the impact of on . Observation shows that increases both the solutal layer and the concentration field. A higher estimation of is a reason for solutal layer thickness. Therefore, increases. The effect of on is graphed in Fig. 16. The figure indicates, increasing the value of reduces . Figure 17 highlights the effect of the solutal relaxation parameter on . A dropped is observed for large values of . Both and have opposite effects on the (see Figs. 18, 19). A considerable increase in is observed when increases. Moreover, the impact of on is graphed in Fig. 20. The figure indicates that a gradual increase in estimates of the Schmidt number results in a thicker boundary layer of concentration. For a large , concentration diffusivity of fluid is increased which is suppressed for the increasing values of the Schmidt number.

Figure 14

versus .

Figure 15

versus .

Figure 16

versus .

Figure 17

versus .

Figure 18

versus .

Figure 19

versus .

Figure 20

versus .

Skin friction coefficient and heat transfer rate

The impact of different emerging thermophysical parameters on the heat transfer rate as well as on the skin friction coefficient is presented in Tables 1 and 2. Friction factor coefficient increases on escalating Weissenberg number and induced magnetic field parameter whereas reduces for the growing value of stretching ratio parameter . Similarly, the heat transfer rate increases on boosting Thermal relaxation factor , , and the Prandtl number and declines for increasing heat source/sink parameter .

Table 1

Numerical estimation of for different parameters.

β	n	A	W_e	f^″(0)1+(f^″)^dn-1dWed
0.1	1.1	0.3	1.5	− 0.4172815
0.2				− 0.49195911
0.3				− 0.57265426
0.4				− 0.66085834
	1.2			− 0.53327392
	1.3			− 0.53897413
	1.4			− 0.54478907
		0.4		− 0.35475913
		0.5		− 0.35037484
		0.6		− 0.27708379
			1.6	− 0.42489466
			1.7	− 0.43325409
			1.8	− 0.44241708

Table 2

Numerical estimation of for different parameters.

λ1∗	N_t	N_b	Q	Pr	Bi	-θ^″(0)
0.2	0.1	0.1	0.3	1.1	0.9	0.13792991
0.3						0.13804022
0.4						0.13818509
0.5						0.13838548
	0.2					0.27877992
	0.3					0.13790234
	0.4					0.13788854
		0.2				0.13795968
		0.3				0.13798945
		0.4				0.13801907
			0.4			0.13784733
			0.5			0.13776472
			0.6			0.13768207
				1.2		0.13794798
				1.3		0.13796686
				1.4		0.13798662
					1.0	0.20425922
					1.1	0.20969592
					1.2	0.2145712

Also, for the accuracy and correctness of the numerical scheme, grid-free analysis for the Nusselt number is presented in Table 3. The current findings of Skin friction coefficient are compared for different values of by Mahapatra and Gupta[47], Ishak et al.[48], Nazar et al.[49], Ali et al.[50] and Gireesha et al.[51] in Table 4 by ignoring induced magnetic field effects and Careau Yasuada nanofluid. For each considered value, comparison table displays good agreement, which provides the validity of the correctness and reliability of the latest results.

Table 3

Grid free analysis for the Nusselt number.

S. no.	Grid size	(Nu_x)_ave,θ
1	10×10	0.17827252
2	50×50	0.17821146
3	100×100	0.17826705
4	300×300	0.12915111
5	500×500	0.12915111
6	600×600	0.12915111
7	1000×1000	0.12915111

Table 4

Comparison of outcomes with Mahapatra and Gupta[47], Ishak et al.[48], Nazar et al.[49], Ali et al.[50], and Gireesha et al.[51]for co-efficient Skin friction considering by varying .

A=a/c	[47]	[48]	[49]	[50]	[51]	Present study
0.1	− 0.9694	− 0.9694	− 0.9694	− 0.9694	− 0.96938	− 0.96933
0.2	− 0.9181	− 0.9181	− 0.9181	− 0.9181	− 0.91810	− 0.91811
0.5	− 0.6673	− 0.6673	− 0.6673	− 0.6673	− 0.66723	− 0.66724
1.0	–	–	–	–	0.90852	0.90853
2.0	2.0175	2.0175	2.0176	2.0175	2.01750	2.01752
3.0	4.7293	4.7294	4.7296	4.7293	4.72928	4.72930
4.0	–	–	–	–	8.00043	8.00046

Conclusions

In this exploration, we have studied the influence of inclined magnetic flux and modified Fourier and Fick's theories are examined on Carreau-Yasuda nanofluid flow induced by a stretching sheet using the Buongiorno model. Additionally, activation energy with binary chemical reaction is introduced to examine the concentration field. Furthermore, heat source/sink effects are considered along with Newtonian heating on the boundary to analyze the heating mechanism. The final remarks drawn from this study are as follows:

The growing values of decline the thickness of the momentum boundary layer and . To boost the rotation parameter, an increase in the axial velocity profile is seen.

boosts on increasing the heat source/sink parameter.

For broad parameters of and , the temperature of nanofluids is increased.

A large thermal relaxation factor rises and also the thermal boundary layer, whereas larger solutal relaxation factor drops the concentration profile.

The opposite performance of the concentration boundary layer is observed on increasing and .

Increasing , boosts the concentration profile.

The friction factor coefficient increases on escalating Weissenberg . Whereas, its value reduces for increasing stretching ratio parameter .

Increasing values of the Thermal relaxation factor and the Prandtl number boosts the heat transfer rate.