Numerical analysis of the chemically reactive EMHD flow of a nanofluid past a bi-directional Riga plate influenced by velocity slips and convective boundary conditions.

This report presents the three-dimensional electromagnetohydrodynamic flow of a zinc-oxide-water nanofluid past a bidirectional Riga plate with velocity slips and thermal and mass convection conditions. The Cattaneo-Christov heat and mas flux model, thermal radiation, chemical reaction and activation energy are considered to analyze the flow problem. The volume fraction of the ZnO nanoparticles is taken 6% in this analysis. An appropriate set of similarity variables is used to transform the partial differential equations into ordinary differential equations. During this process, some parameters are found and influences of these factors on the flow profiles are shown and discussed in detail. A numerical technique called NDSolve is considered for the solution of the nanofluid flow problem. The results showed that higher solid volume fraction and slip parameter have reduced velocities profiles and the increasing solid volume fraction and thermal Biot number have increased the temperature profile. Additionally, the concentration Biot number has increased the concentration profile. The modified Hartmann number has significantly increased the velocity profile. Dual impacts in velocity profiles along primary and secondary direction has been observed due to stretching ratio parameter. A comparison of current results has been carried with a fine agreement amongst current and established results.

Introduction

Nowadays researchers have created numerous easier and less expensive methods to create nanoparticles of economically significant materials as a result of the special features of nanotechnology. A number of metal oxide nanoparticles are created with potential future uses. One of them is zinc oxide which is one of the best utilized at the nanoscale. ZnO is an inorganic substance having numerous applications at the industrial level such as capacitors, cerium, porcelain, glass, concretes, polymers, polyester, adhesives, and colors etc. ZnO is a semiconductor in nature and has various distinctive qualities, including strong room-temperature fluorescence, high mobility of electrons and superb transparency etc. Zinc oxide is crucial for both scientific and industrial applications due to its wide band gap and high excitonic binding energy.

In past few decades the rheological and thermal characteristics of nanoparticles base fluids have been captured more consideration due to their important applications at industrial level. Such fluids are termed as nanofluid and can be manufactured by suspending a single type of nanoparticles in pure fluids. It has been shown that such fluids have tremendous thermal conductivities due to which these fluids are using for coolant purposes such as, coolant of heat exchangers, nuclear reactors, electrical devices and coolant of auto engines etc. Choi and Eastman[1] has presented the nano particulates to pure fluid for the enhancement its thermal conductivities. Ayub et al.[2] have discussed electrically conducted nanoparticles-based fluid over a Riga plate and have revealed that flow of fluid has declined while thermal flow has supported by the augmenting values of nanoparticles volumetric fractions. Ramesh et al.[3] have assessed the bioconvective Maxwell nanofluid flow past a Riga sheet by incorporating nonlinear thermal radiations and activation energy effects. It has highlighted by the authors that greater activated energy factor has enhanced the concentration of fluid. Ali et al.[4] have introduced rotary transient nanofluid flow upon Riga plate by considering the effects of microbes, chemical reaction and non-Fourier thermal flux and noticed that the thermal curve rises with increment of thermophoresis and Brownian effects and have retarded with higher values of magnetic effects. Singh and Ghosh et al.[5] provided a thorough numerical and experimental analysis, in order to assess the thermal functionality of the 60° and 30° Chevron plate geometries employing (MWCNT)/distilled water nanoliquid as refrigerant. Pal and Mandal[6] have conducted melted heat transportation phenomenon for Sisko nanofluid flow upon stretching surface. Waqas et al.[7] inspected numerically the bioconvective flow of nanofluid flow with impact of activated energy past a Riga surface and noticed that energy communication rate decayed by augmentation in Prandtl number. Singh et al.[8,9] reported the computational and experimental study employing graphene oxide nanoliquids at various compositions and flow rates. The research revealed that utilising nanofluid caused a pressure losses and increased pouring energy, but improved thermal conduction, thermal efficiency, performance, and ultimate energy transfer index by 13%, 14%, 9%, and 10%, respectively. The dynamic viscosity and density of Al2O3/distilled water, (MWCNT)/DW, and graphene nano particulate nanoliquids was described by Singh and Ghosh[10,11].

The Riga surface is an eminent actuator that comprises of ever-holding constant electrodes and magnet producing Lorentz forces, where such forces weaken swiftly at some distance from the surface known as Riga Plate. For its important applications at industrial level, researchers have been convinced to conduct numerous investigations for fluid flow upon Riga plates. Bhatti and Michaelides[12] discussed the influences of activated energy on thermal bioconvective nanoparticles flow over a Riga plate and have revealed that higher values of Peclet and bioconvective Schmidt numbers have declined the profiles of microbes. Shafiq et al.[13] have studied nanoliquid flow past a radiated Riga plate and have revealed that the augmenting values of radiated factor have enhanced the thermal profiles. Rasool et al.[14] have investigated Marangoni based forced convection nanofluid flow past a Riga plate influenced by magnetic effects. It has shown in this study that effective Lorentz force has controlled the fluid motion and thermal flow. Mburu et al.[15] have inspected numerically the generation of irreversibility and heat transmission for nanoparticles flow upon Riga plate. Pal and Mandal[16] have inspected mixed convection electrically conducted nanofluid flow past a shrinking/stretching surface with heat source and viscous dissipation.

The collective deliberation of mass as well as thermal transportation has a considerable importance in industrial and engineering applications for instance, nuclear reactor coolant, thermal conduction in tissues, electronic devices coolant, food processing and air conditioning etc. The classical theories of Fourier[17] and Fick[18] have provided a base for heat and mass transportation problems. But since relaxation time are changing that has to affect thermal as well as concentration fields, so these two concepts were not more suitable for thermal and mass transportation problems. Cattaneo[19] primarily has modified the theory of Fourier’s for thermal diffusions. This idea of Cattaneo was further modifies by Christov[20] by exchanging time-based partial derivative with Oldroyd upper convection derivatives. Afterwards, a number of studies have been published in literature by considering the influence of heat and mass flux model of Cattaneo–Christov in the fluid flow models. Naseem et al.[21] have inspected analytically the nanoliquid flow past a Riga plate using Cattaneo–Christov (CC) approach and have revealed that the thermal as well as concentration characteristics have been augmented with upsurge in relaxation parameters. Pal and Mandal et al.[22] examined the magnetic flux in a permeable media under the impact of thermal radiation and varying thermal conductivity, to investigate the thermal performance of nanofluid over a stretched surface. Rasool and Wakif[23] have inspected EMHD mixed convection flow of nanoparticles-based fluid towards vertically placed Riga plate using CC model. Their outcomes have been indicated that the escalating values of Hartman number have augmented the horizontal motion of fluid while thermal flow and concentration of nanoparticles have been reduced with higher values of relaxation time factor.

The smallest quantity of energy that is required for different compounds to commence a chemical reaction is termed as activation energy introduced first by Svante Arrhenius for instance, the energy required by a car to start its engine etc. Due to its important industrial and engineering applications many investigations have been conducted by various researchers to discuss the mass diffusions for fluid flow problems. Ali et al.[24] have discussed bioconvective rotary motion of Maxwell nanoparticles over a Riga plate with influence of activation energy. It has established by the authors that the thermal flow rate has been augmented with upsurge in the values of Brownian and thermophoresis factors while concentration of fluid has been reduced by higher values of activation factor of energy. Shahid et al.[25] have discussed experimentally the influence of activation energy nanofluid flow upon a porous surface and have determined the numerical outcomes for considered flow problem. Xia et al.[26] discoursed the energy transmissions for nonlinear mixed convection nanofluid flow with influence of activation energy subject to multiple slip conditions at boundaries. Their outcomes have depicted that the thermal flow rate have been upsurge while fluid motion have been declined with higher values of Hartmann number and nanoparticles. Riaz et al.[27] have studied the Cattaneo–Christov model for nonlinear convection nanofluid flow over an extending surface subject to activation energy and have revealed that concentration of fluid flow has been upsurge with growth in activation energy factor.

Thermal radiations play a pivotal role in transportation of heat. The effects of these radiations are pretty imperative in high thermal processes. Rehman et al.[28] have deliberated the influence of thermal radiations upon fluid flow pattern over flat cylindrical surface and have concluded that thermal flow rate has been amplified with upsurge in the values of thermal radiation factor. Ahmad et al.[29] have inspected the impacts of thermal radiations upon mass and heat transfers for fluid flowing in an infinite circular pipe. It has concluded by the authors of the work that fluid motion has retarded whereas the thermal flow rate has been enlarged with augmenting values of thermal radiation factors. Rehman et al.[30] have conducted a comparative study for thermal transmissions of MHD Jeffery fluid flow with thermal radiations over cylinder shaped and plane sheets. Ashraf et al.[31] have discussed thermally radiative nanofluid flow with influence of Darcy-Forchheimer model. Eid and Nafe[32] inspected the variations in thermal conductivities and production of heat upon MHD nanofluid flow through porous medium and have determined that augmentation in thermal radiations and nanoparticles have augmented Nusselt number and while skin friction has retarded in this phenomenon. The readers can further study literature upon flow control with influence of thermal radiations in Refs.[33-38].

After a careful analysis of the above-mentioned literature survey, we have observed that very less work based on the three-dimensional electromagnetohydrodynamic nanofluid flow past a bi-directional Riga plate with velocity slip and thermal and mass convective conditions is considered. Therefore, the authors have proposed this model. The flow contains the ZnO nanoparticles, which are suspended into water and analyzed through a bi-directional Riga plate with velocity slips and convection conditions. The consequences of Cattaneo–Christov heat and mas flux model, thermal radiation, chemical reaction and Arrhenius activation energy are studied. An appropriate set of variables is considered to convert the equations which governed the flow problem into dimension free form. A numerical technique called NDSolve is considered for the solution of the nanofluid flow problem. This work is composed of several sections which compiles the flow problem. “Problem formulation” section  represents the problem formulation. “Numerical solution” section  displays the numerical investigation of the present problem. “Results and discussion” section  shows the results and discussion of the flow problem. “Conclusion” section  is composed of concluding remarks.

Problem formulation

Consider the three dimensional EMHD flow of nanofluid through a bi-directional Riga plate as shown in Fig. 1b. The stretching velocities along x- and y- axes are respectively described as and , where and are constants such that both quantities are positive. The temperature and concentration of nanofluid at the surface of the plate are represented by and , respectively while at free stream these values are depicted by and , respectively such that , . Furthermore, the following assumptions are taken into consideration.

Figure 1

(a) Flow chart. (b) Flow geometry.

Riga plate

Non-Fick’s and non-Fourier model

Arrhenius activation energy and thermal radiation

Slip and convection conditions

With above stated suppositions the leading equations are[39-41]:where

The boundary conditions are:

The radiative heat flux is defined as:

Expending by Taylor series, we have:

By ignoring higher terms, Eq. (10) can be reduced as:

Replacing Eq. (11) into (4), we have:

The thermophysical properties are defined as:

The thermophysical properties of pure fluid and nanoparticles are described numerically in Table 1.

Table 1

Numerical values of thermophysical properties of and [42–45].

Base fluid and nanoparticle	ρKgm^-3	C_pJKg^-1K^-1	kWm^-1K^-1
H_2O	997.1	4179	0.613
ZnO	5600	495.2	13

The set of appropriate variables are described as:

Making use of Eq. (14) the leading equations are reduced as:

With conditions at boundaries as:

Above the dimension-free parameters are mathematically described as:

Quantities Interest

The skin friction coefficients, local Nusselt and Sherwood numbers are described mathematically as:where

Using Eq. (14), the quantities of the interest are reduced as:

Above and are local Reynolds numbers.

Also, and .

Numerical solution

When a physical phenomenon is modelled mathematically, it gave rise to nonlinear equation. Sometimes the nonlinearity of resultant equations is much higher. So, such equations are very difficult to solve analytically. In current problem the resultant Eqs. (15–18) are nonlinear for which NDsolve[46,47] has used which is a numerical technique in Mathematica and is used to solve such problems efficiently.

Results and discussion

This work addresses the 3D electromagnetohydrodynamic nanofluid flow past a bi-directional Riga plate with slip and convective boundary conditions. Nanofluid is composed of zinc-oxide nanoparticles which are suspended in water. The CC model for heat and mass flux is used in a flow problem. Moreover, thermal radiation is incorporated in energy equation and the effects of E and K are employed in concentration equation. The numerical investigation of the flow is incorporated with the help of NDSolve technique. Figure 1a shows the Flow chart. Figure 2 presents the influence of volumetric fraction of nanoparticles on velocity profile . Since with a growing , the density of fluid augments due to which the fluid becomes denser. In this phenomenon fluid experienced more friction due to resistive force in the direction of motion. As a result, velocity profile reduces. A similar impact of is depicted against velocity profile as shown in Fig. 3. Figure 4 shows the impact of on . The increasing escalates . Physically, the thermal conductivity of enhances with the increasing which results the enhancement in the thermal boundary layer thickness and temperature profile. Therefore, the increasing increases . Figure 5 shows the effect of modified Hartmann number on velocity profile . The increasing increases . The momentum boundary layer thickness also increases with the increasing . Because the modified Hartmann number has bigger values, the external electric field is increased, maximizing the velocity distribution . Similar impacts were found in[39,48,49]. The influences of stretching ratio factor on and are depicted in Figs. 6 and 7. It is observed that the increasing reduces the velocity profile , while increases the velocity profile . The reason is that the stretching parameter has direct relation with the stretching velocity constant along y-direction and inverse relation with the stretching velocity constant along x-direction. Therefore, a decreasing impact is found for , while an increasing impact found for . Similar impacts were found in[48,49]. Figures 8 and 9 depict the influence of slip parameter on and . It has observed that and declines with higher values of . Actually, when the values of become higher, then some of the stretching velocities transformed to the fluid particles that results in retardation of fluid motion in all directions. Due to the decline in momentum boundary layer thicknesses the velocities and are declined. Figure 10 depicts that with augmentation in thermal Biot number there is a minimization in thermal resistance amongst nanoparticles at the surface of a Riga plate. Hence the thermal distribution at the Riga plate upsurge due to augmentation in , that results in maximizing the thermal flow profile as portrayed in Fig. 10. Similarly the growth in the values of concentration Biot number leads to augmentation in concentration distribution as shown in Fig. 11. Figure 12 presents the influence of thermal relaxation factor on thermal flow profile. It is noticed from this figure that the upsurge in retards thermal characteristics. Actually, this parameter is employed as an indicator for estimating the quantity of time that is required for heat transmission from one zone to another one of the same materials. When is enhanced then more time is required for thermal flow from one zone to another, hence there is a retarding behavior in heat transportation. Therefore, the increasing decreases the thermal profile. The influence of mass relaxation factor on concentration profile is depicted in Fig. 13. The factor is actually an indicator used for estimation of amount of time needed for mass diffusions from higher to lower concentration zone. The augmentation in indicates that more time is required for transmission of mass from one region to another that ultimately declines the concentration of nanofluid. Figure 14 depicts the influence of radiation factor on thermal profile. It is obvious from the figure that heat transmission jumps up with higher values of . Physically this can be interpreted as the growth in augments the heat energy amongst nanoparticles closed to the Riga plate. In this process the charge on these nanoparticles augments, that ultimately grows up the thermal profile. Figure 15 portrays the influence of activation energy factor on concentration profile. For higher values of , a bulk of molecules that entails small quantity of energy supports diffusions of mass. Therefore, higher values of results a growth in the concentration profile. The influence of chemical reaction parameter on concentration profile is depicted in Fig. 16. Intensifying values of results the retardation of mass diffusivity that weakens the thickness of concentration boundary layer thickness. Hence, higher values of causes a decline in concentration profile. Using variations in different parameters the results of current investigation have been compared with established results in the literature. In Tables 2 and 3, a comparison is carried out for and with the results of Iqbal et al.[50], Khan et al.[51] and Makinde and Aziz[52] for different values of stretching ratio parameter, while keeping other parameters zero. A fine agreement of current results with those established in literature has been found. In Table 4, it has noticed that the values of and have been retarded with augmentation in the values of solid volume fraction, modified Hartmann number, stretching ratio and slip factors. The numerical impacts of radiation parameter, solid volume fraction and thermal Biot number on Nusselt number have been described. In Table 5, it is noticed that the resistive force to fluid motion has been increased for augmentation of these stated parameters. Therefore, the values of Nusselt number has been upsurge due to augmenting values of radiation parameter, solid volume fraction and thermal Biot number.

Figure 2

Effect of on .

Figure 3

Effect of on .

Figure 4

Effect of on .

Figure 5

Effect of on .

Figure 6

Effect of on .

Figure 7

Effect of on .

Figure 8

Effect of slip parameter on .

Figure 9

Effect of slip parameter on .

Figure 10

Influence of thermal Biot number upon .

Figure 11

Influence of concentration Biot number upon .

Figure 12

Effect of on .

Figure 13

Effect of on .

Figure 14

Effect of on .

Figure 15

Effect of on .

Figure 16

Effect of on .

Table 2

Comparison of current results of with previously reported results, when .

λ	Iqbal et al.[50]	Khan et al.[51]	Makinde and Aziz[52]	Present results
0.00	1.0000	1.0000	1.0000	1.00000
0.25	1.048813	1.048813	1.048812	1.05428
0.50	1.093097	1.093095	1.093095	1.09731
0.75	1.134485	1.134485	1.134485	1.13782
1.00	1.173720	1.173721	1.173721	1.17641

Table 3

Comparison of current results of with previously reported results, when .

λ	Iqbal et al.[50]	Khan et al.[51]	Makinde and Aziz[52]	Present results
0.00	0.000000	0.000000	0.000000	0.000000
0.25	0.194564	0.194564	0.194564	0.197344
0.50	0.465205	0.465205	0.465205	0.468343
0.75	0.794622	0.794622	0.794622	0.797595
1.00	1.173720	1.173720	1.173720	1.176410

Table 4

Impacts of the embedded parameters on and .

Π	α	H	λ	C_sx	C_sy
0.01	0.1	0.9	0.5	0.817115	0.409982
0.03				0.795033	0.394869
0.04				0.782609	0.386907
0.06				0.755505	0.370360
	0.2			0.663296	0.330800
	0.3			0.592753	0.299782
	0.4			0.536768	0.274672
		0.2		0.589796	0.272797
		0.5		0.566982	0.273606
		0.7		0.551849	0.274141
			0.2	0.518715	0.095370
			0.4	0.531060	0.211124
			0.6	0.542225	0.341353

Table 5

Impacts of the embedded parameters on .

Π	Rd	Bi_T	Nu
0.01	0.1	0.1	0.104389
0.03			0.105025
0.04			0.106852
0.06			0.108529
	0.2		0.111480
	0.3		0.119405
	0.4		0.127304
		0.2	0.197853
		0.3	0.284146
		0.4	0.363343

Conclusion

This work addresses the three-dimensional electromagnetohydrodynamic flow of a water-based zinc-oxide nanofluid over a bi-directional Riga plate with velocity slips and thermal convective boundary conditions. The Cattaneo–Christov heat and mass flux model is taken into consideration in order to analyze the thermal and concentration profiles. Moreover, thermal radiation effect is consideration in energy equation and the effects of Arrhenius activation energy and chemical reaction are taken in concentration equation. A suitable set of similarity variables is used to transform the equations that governed the flow problem, into dimension free form. During this process some parameters are created and influences of these factors on flow profiles are discussed. A numerical technique called NDSolve is considered for the solution of the nanofluid flow problem. The results showed that higher solid volume fraction and slip parameter have reduced velocities profiles and the increasing solid volume fraction and thermal Biot number have increased the temperature profile. Additionally, the concentration Biot number has increased the concentration profile. The modified Hartmann number has significantly increased the velocity profile. Dual impacts in velocity profiles along primary and secondary direction has been observed due to stretching ratio parameter. The increasing thermal and mass relaxation factors have augmented the temperature and concentration profiles, respectively. The concentration profile is increased due to chemical reaction parameter while increased due to activation energy parameter. A comparison of current results has been carried with a fine agreement amongst current and established results. The Nusselt number is upsurge due to augmenting values of radiation parameter, solid volume fraction and thermal Biot number.