Effect of inclined magnetic field on radiative heat and mass transfer in chemically reactive hybrid nanofluid flow due to dual stretching

This research analyzes the three-dimensional magneto hydrodynamic nanofluid flow through chemical reaction and thermal radiation above the dual stretching surface in the presence of an inclined magnetic field. Different rotational nanofluid and hybrid nanofluids with constant angular velocity \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\upomega }^{*}$$\end{document}ω∗ for this comparative study are considered. The constitutive relations are used to gain the equations of motion, energy, and concentration. This flow governing extremely non-linear equations cannot be handled by an analytical solution. So, these equations are transformed into ordinary differential equalities by using the similarity transformation and then handled in MATLAB by applying the boundary values problem practice. The outcomes for the considered problem are accessed through tables and graphs for different parameters. A maximum heat transfer amount is observed in the absence of thermal radiation and when the inclined magnetic field and axis of rotation are parallel.

Kinematic viscosity µ f Dynamic viscosity r, s Temperature and concentration profile 3D Three dimensional Greek symbols α 1 Temperature diffusivity Rotational velocity ρ Density σ Electrical conductivity τ w Shear stress ǫ Mixed convection parameter Nanofluid was first devised by Choi and Eastman 1 in the last decade of the 20th century. They revealed that when a nanoparticle having good thermal conductivity is dispersed into a base fluid, the thermal conductivity of the formed solution exceptionally increases. They remarked that pumping power for heat transport in a heat exchanger considerably declines when the nanofluid is utilized as an working fluid in a heat exchanger instead of conventional working fluids like water, oils, and ethylene glycol. Further, when two nanoparticles of nanometer are distributed into a host fluid, the formed mixture is known as a hybrid nanofluid 2 . Hayat and Nadeem 3 elaborated on the hybrid nanofluid for heat transfer improvement for rotating flow using silver and copper oxide nanoparticles. Shah and Ali 4 provided a comprehensive review of hybrid nanofluids and their applications which involve industrial applications, like wire drawing, coolant in engines of automobiles, nuclear reactors, micro-chips in computers, hot rolling, cancer therapeutics, glass fiber production, as a detergent, etc. Yasmin et al. 5 experimentally explored hybrid nanofluids for solar and thermal energy storage uses. Maddammasetty and Sireesha 6 used thermal systems to elaborate the heat transport applications of hybrid nanofluids. Numerous researchers have investigated the application of nano and hybrid nanoliquids in recent years. Some of the recently studied nano and hybrid nanofluids are given for knowledge gains [7][8][9][10][11][12] . There is vast room for research in this regard. Ghasemi et al. 13 studied nanofluid flow over the stretched surface with radiation and magnetic field effect using a novel spectral relaxation method. They remarked that high levels of magnetic field dramatically affect temperature and concentration profile. Sharma et al. 14 analyzed graphene Maxwell nanofluid past stretching surface. Viscous dissipation and unsteadiness reduce heat transfer rates. Hussain et al. 15 studied three-dimensional nanofluid flow with a magneto effect. Zinc nano-particle embedded nanofluid displayed a higher rate of heat transmission with an increment in magnetization force. Arshad and Hassan 16 explored hybrid nanofluids between rotating systems. They discovered that when two nanoparticles are present in a host fluid, the amount of heat transfer enhances. Hassan et al. 17 explored hybrid nanofluid for prescribed wall temperature cases with thermal radiation. Reduced drag and lift coefficients are obtained when the silver nanoparticle is used with a single-wall carbon nanotube. Hady et al. 18 elaborated on the radiation's influence on heat transfer in viscous movement across a stretching surface. TiO 2 -embedded nanofluids have higher cooling power as compared to other examined nanofluids. Ali et al. 19 studied electro-magneto-hydrodynamic nanofluid flow with variable heat fluxes. They remarked that the electric field directly influences the temperature profile. In recent years, Hassan 20 , Hussain 21 , Sheikholeslami 22 , Masood and Farooq 23 and Masood 24 have explored hybrid nanofluid in the presence of different body forces.
There are two types of chemical reactions specifically, irreversible, and reversible chemical reactions. Irreversible chemical reactions are chemical reactions that can not return to their initial stage. Whereas the reversible chemical reaction can return to its initial stage in the presence of a catalyst 25 . Recently, different-order chemical reactions are explored by researchers to investigate their impact on flow regimes. Anjum et al. 26 investigated binary chemical reaction effect fluid flow with double stratification impacts. They found out that skin friction enhanced with an augmentation in magnetizing force. Abbas et al. 27 explained the influence of chemical reactions on heat transmission in third-grade fluid over an exponential stretching surface. Temperature enhances with a rise in the level of chemical response they remarked. Elattar et al. 28 investigated the hybrid nanoliquid movement with hall current and chemical reaction effects. Species transportation enhanced with increment in chemical reaction they noted. Krisna et al. 29 explained the chemical consequence and radiation effect on convective stream with sucking and heat-making effects. The amount of heat transfer is highly influenced by the www.nature.com/scientificreports/ chemical reaction and magnetic force. Recently, second-order and higher-order chemical reactions have been investigated by numerous researchers 30-32 . In the 20th span, Prandtl 33 pioneered the idea of boundary layer flow in fluid dynamics. The boundary layer is the layer of the fluid that forms in the surrounding area of the surface bounding the fluid. Most common examples of boundary layer flows include, near the earth's surface, the interior of water pipes, and inside the blood vessels 34 . Khan and Pop 35 investigated the boundary layer flow past a stretching surface. They remarked that the Nusselt number is minimizing the phenomenon of thermophoresis, Brownian, Prandtl, and Lewis numbers. Crane et al. 36 analyzed flow past the stretching surface, Bongar 37 examined nanofluid flow past the stretching surface, Erickson 38 explored heat and mass transmission with suction effect on a continuous level plate, and Sakiadis 39 investigated axisymmetric flow in two-dimensional form for boundary layer flow behavior.
Chen and Stroble 40 explored the Buoyancy force on the boundary layer flow on an continuously moving plate. Takhar and Nath 41 discussed three-dimensional flow due to stretching surfaces. Wang et al. 42 examined viscous flow above a stretchable surface together with slip and suction effects. They stated that the width of the boundary layer enhanced with minimal mass suction effect Mehmood and Ali 43 explored analytical solutions of viscous flow with heat transfer. Shahzad et al. 44 examined heat transfer due to stretching surfaces using nanoparticles with the MHD effect. They remarked that rotation enhances the skin friction of rotating nanofluids. Hassan et al. 45 studied linear and nonlinear radiation effects on heat and mass transmission in a hybrid nanoliquid due to a stretchable surface. Hussain et al. 46 discussed the magnetic and nonlinear thermal radiative effects on three-dimensional movement due to stretching surfaces. Arshad et al. 47 examined magneto-hydrodynamic flow with chemical effect above the exponential stretching surface. Khan et al. 48 studied heat and mass transmission in Burger's nanofluid with magnetization and chemical reaction over an exponential stretchable surface. Numerous researchers have explored stretching surfaces for heat transfer and analysis of fluid flow [49][50][51][52][53] .
The above-conducted literature review suggests that numerous researchers have investigated the flow past the stretching surface. Arshad 53 explored the thermophoresis and Brownian motion with thermal radiative effect and uniform magnetic field. Arshad 47 investigated the chemical reaction effect over an exponential stretching surface. The novelty of this paper is to investigate the comparative dynamics of rotating water-based nano and hybrid nanofluids over dual stretching surfaces implanted in a permeable medium with radiative heat and mass transfer. Arshad 47,53 examined a uniform magnetic field, whereas in this study we have incorporated an inclined magnetic field with a chemical reaction. The prevailing equations are changed into the ordinary differential equation by employing a similarity transformation and attempted MATLAB by utilizing the boundary value problem method. The tolerance is set to be 10 −6 for obtained solutions. The flow behavior and characteristics of magneto hydrodynamic nanofluid and hybrid nanofluid are comparatively presented through tables and graphs for different parameters. By increasing parameters, the skin friction, Nusselt number, and Sherwood number are evaluated. This relative study helps to answer the following main research questions.
1. What is the influence of rotation parameters on velocity, temperature, and concentration profiles? 2. How does the increasing behavior of the magnetic force and porous medium parameter affect the primary and secondary velocity profile? 3. Do the increasing behavior of the radiation parameter and Prandtl number give the increased heat transmission rate and minimum skin friction? 4. How does chemical reaction, Lewis number change the temperature and concentration profiles? 5. What is the changing behavior on skin frictions along the x − axis , y − axis, Nusselt number as well as the Sherwood number of different parameters?

Basic equations
The basic flow governing equations for viscid incompressible liquid over the permeable surface in the existence of Rosseland radiation, and chemical reaction. (1) ∇.V = 0.

Mathematical formulation of problem
Consider an incompressible, steady, three-dimensional water-based nanofluid flowing above a porous stretchable sheet. A cartesian coordinate system x, y, z is considered to discuss the problem physically. The x-axis is taken in the horizontal direction, the z-axis is upward and the y-axis is perpendicular to both other axes. An inclined magnetic field B 0 in the z-axis direction with angle α is acting. The nanofluid is rotating at a constant speed with ω * along the z-axis, chemical reaction and, thermal radiation are considered. The flow is induced by a stretching sheet with speeds U w = ax and V w = by in the x-direction and the y-direction (see Fig. 1). Nanofluid and hybrid nanofluid are considered for this comparative study. The governing equation along with these considerations takes the following form 53 : Equation of continuity 53 : Momentum equations along the x and y axis 53 : Energy equation without q r relation 53 :

Concentration equation 53 :
Here u, v , and w are velocity components in the x, y , and z directions respectively. T and C are temperature and concentration.k 0 , q r , g * , B 0 , andK r is the porosity, radiation, gravitational acceleration, magnetic field, and chemical reaction respectively.ρ hnf , µ hnf , α 1 is the density, viscosity, and thermal diffusivity of hybrid nanofluid. By applying the Rosseland approximation, the radiative heat flux q r is defined by 53 : Here σ * is the Stefan-Boltzmann coefficient and k 1 is the mean absorption constant. Assuming that there is not much of a temperature differential within the flow, the expression of the term T 4 by using Taylor series expansion is expanded as follows: www.nature.com/scientificreports/ Consequently, by ignoring higher order terms above the first order in (T − T ∞ ) , we get.
By using Eqs. (11) and (12) Therefore, the energy equation takes the form.
The corresponding boundary conditions are: Similarity transformation. We define the following similarity transformation to transform the equation which are tackled numerically.
The non-dimensional quantities B 1 , B 2 , B 3 , B 4 and K 2 , are double nanoparticle relations (defined in Table 1), and , δ, Z, ǫ, M, Pr, π , Le, N t , τ , N b , and K c are defined as The modified boundary conditions are as follows: Here γ = b a is the dimensionless stretching ratio.
Physical quantities of interest. The most significant physical quantities of importance from an engineering perspective are the skin friction coefficients Cf x , Cf y , and Nusselt number Nu , which are defined as follows: www.nature.com/scientificreports/ Here τ zx and τ zy denote the shear stress along the stretched wall along the x-axis and y-axis and are defined as The dimensionless form of Eq. (18) is: By using the temperature field to define the thermal diffusion rate as a Nusselt number: Or By using the concentration field to define the mass transmission rate as Sherwood number: By applying the resemblance transformation Eq. (16), the non-dimensional form of the Sherwood number takes the form: The following Tables 1 and 2 show the thermophysical relations and values for the formulation of nanofluid and hybrid nanofluid. Table 1. Thermophysical relations of nanoparticles and base fluid 53 .

Properties Nanofluid
Density Heat capacity

Properties Hybrid nanofluid
Density www.nature.com/scientificreports/ Spherical-shaped nanoparticles are used having the shape effect 3 . Additionally, the volumetric concentration of the nanoparticle used is about 0.005% for each nanoparticle.
Numerical scheme and validation. The boundary value problem technique is used. Highly accurate and effective numeric outcomes are obtained (Table 3) when we use this technique. Equations (17)- (20) which are highly non-linear are changed to first-order ODEs by using a new set of variables defined as follows.
The equation y 3 ′ , y 6 ′ , y 8 ′ and y 10 ′ takes the following form: The transformed boundary conditions changed into the subsequent form: The skin friction, Nusselt, and Sherwood numbers changed into the following form:

Results and discussion
The interpretation of the problem introduced in the preceding section is hybrid nanofluid. The outputs of the present problem are obtained separately for each nanofluid i.e., copper-water nanofluid, copper/aluminum oxide-water based hybrid nanofluid. The obtained effects of different parameters by utilizing the boundary value problem technique at MATLAB are described in this segment.
y 1 = 0, y 2 = 1, y 4 = 0, y 5 = γ , y 7 = 1, atη = 0 y 2 → 0, y 5 → 0, y 7 → 0, asη → ∞    www.nature.com/scientificreports/ Effect of rotation and stretching ratio parameter. The effect of rotation constraint on velocity constituents p ′ (η), q ′ (η) , temperature r(η) , and concentration s(η) is presented in the following Fig. 2a-d respectively. The velocity profiles p ′ (η) and q ′ (η) are evident from the graph that it decays when rotation parameter increases for nanofluid and hybrid nanofluid. At the start when = 0 and = 1 , there is a minimum change in velocity profiles. When rotation increases by = 2 and = 3, these profiles decay rapidly. The specific reason for the occurrence of this event is that rotation is a direct function of the angular velocity. Additionally, the minimum resistive Lorentz force is present when the rotation is increased. This effect expands the momentum layer in the primary direction whereas it is reduced in the secondary direction. An opposite behavior is noted for the impact of rotation constraint on temperature profile r(η) and concentration profile s(η) respectively. Temperature and concentration profiles increase when the rotation parameter increases. Here least rise is noted for the concentration profile as compared to the temperature profile. It is worth mentioning here that when rotation is enhanced the associated thermal boundary layer contracts for both nanofluid and hybrid nanofluids. Additionally, the presence of thermal radiation has a direct impact on the temperature profile. Moreover, it declines when the rotation and magnetization force is high. It is noted that under the high influence of rotation, the concentration profile declines rapidly. Chemical reactions aid the smooth movement of the concentration of nano-particles when the rotational motion of the fluid is augmented. Figure 2e-h show the effect of stretching ratio constraints on velocities, temperature, and concentration profiles, respectively. When the stretching is enhanced velocity in the x-direction decline. Additionally, when stretching is enhanced the momentum layer associated with secondary velocity expands. Stretching has the opposite effect on the primary and secondary velocity profile. A slight difference has been observed in the dynamics of nano and hybrid nanoparticles. The temperature and concentration of both profiles show decreasing behavior for rising the stretching ratio parameter. It is fascinating to note that a high-temperature profile can be observed for a hybrid nanofluid, and a high-concentration profile can be observed for a single nanoparticle nanofluid. Additionally, the temperature boundary layer has contracted under the high increment of stretching ratio. Concentration profile decline with increment in the stretching ratio and the associated concentration layer contract due to the influence of chemical reaction.
Effect of magnetic force and porosity. The following Figs. 3a-c present the effect of magnetic force M on velocity profile p ′ (η), q ′ (η) , and temperature profile r(η) , respectively. In the non-existence of a magnetic force (M = 0), the fluid flow very smoothly, and when the magnetic field begins to work the fluid velocity profiles decays rapidly. This is owing to the Lorentz force acting on fluid which restricts the fluid to move and causes endurance and as a result, the velocity profiles decay under the rising behavior of the magnetic field constraint. Higher momentum border layer wideness is seen for hybrid nanoliquid in the non-existence of a magnetic force. An opposite behavior is presented for magnetic field constraint on temperature profile r(η). The temperature profile has a direct relation with magnetic field constraint. As the magnetic force parameter grows, the temperature profile increases in the same manner and higher thermal frontier layer wideness is detected for hybrid nanofluid. The reason is that the Lorentz force restricts the flow which allows the fluid to transfer a higher amount of heat transmission. The change of porous medium parameter Z on velocity profiles p ′ (η) , q ′ (η) , and temperature profile r(η) is represented in Figs. 3d-f. The momentum boundary layer thickness has an inverse relation with the porous medium parameter. Both the velocity profiles p ′ (η) and q ′ (η) decreases when the permeability of the medium increases because the velocity is dependent on the porosity of the surface. There is an inverse relation between the permeability parameter and temperature profile r(η) . When the porosity of the surface increases the thermal boundary layer expands as an outcome the fluids conduct more heat. In both cases, greater momentum and thermal boundary layer are noted for the hybrid nanofluid as associated with the nanofluid. This shows the high performance of hybrid nanoparticles nanofluid.
Effect of mixed convection and inclined magnetic field. The following Figs. 4a-d shows the change of mixed convection constraint ǫ on velocity, temperature, and concentration profiles p ′ (η), q ′ (η), r(η) and s(η) respectively. Both velocities profiles p ′ (η) and q ′ (η) has a linear relationship with mixed convection constraint. In the non-existence of mixed convection, the nanofluid and hybrid nanofluid move very smoothly. When the mixed convection parameter upsurges to a non-zero value the velocities profiles increase in the same way. Since mixed convection magnifies the buoyancy force and as a result, the velocity profiles increase. The higher momentum boundary layer is noted for nanofluid due to the presence of only a single nanoparticle in a base fluid. Due to the lower density of single nanoparticle nanofluid, it can move easily as compared to hybrid nanoparticle nanofluid. While inverse conduct is noted for temperature profile r(η) and concentration profile s(η) for the increasing estimates of mixed convection constraint. The temperature and concentration of both profiles decay by strengthening the mixed convection constraint. The buoyancy forces dominate the inertial forces, as a result, the temperature and concentration profile reduce. Moreover, it is interesting to note down in the temperature profile, a wider thermal boundary layer is detected for hybrid nanoparticle nanofluid as associated to the single nanoparticle nanofluid. Figure 4e shows the influence of the angle of inclination of the magnetic field acting on the rotational nanofluid and hybrid nanofluid. This provides a magnificent consequence of velocity profile p ′ (η) . There is a trigonometric "sin" function that is involved in magnetic field strength. As the angle of the inclined magnetic field increases from 0 0 → 90 0 , the inclined magnetic field strength terms get increasing values and as a result, it helps the nanofluid to reduce the speed because the magnitude of resistive force increases. Figure 5a shows the impression of chemical reaction on the concentration profile. In the absence of chemical response, the highest concentration border layer is noted for hybrid nanofluid. When the chemical reaction constraint begins to rise, www.nature.com/scientificreports/ the concentration profiles start to decay because when the chemical reaction rate increases, a higher volume fraction of nanoparticles undergoes the chemical reaction. A maximum decay is observed for the single-particle nanofluid. So, a consistent mass transfer rate is examined for the hybrid nanoliquid. Figure 5b indicates the change in Prandtl number on the temperature profile. As the significance of the Prandtl number boosts the temperature boundary layer declines since the thermal diffusivity decreases with a rise in the Prandtl number. A higher temperature boundary layer is observed for hybrid nanofluid as associated with nanofluid. Figure 5c shows the influence of Lewis number on the concentration profile. It is defined as the proportion of thermal diffusivity to mass diffusivity. So, the concentration profile decreases once the Lewis number rises. The most important factor of this research is the thermal radiation constraint whose influence on the temperature profile is indicated in Fig. 5d. In the absence of thermal radiation, maximum thermal boundary layer wideness is notable for hybrid nanoliquids over the permeable stretchable surface. When the value of the thermal radiation parameter π increases, a sudden fall in the thermal boundary layer is seen for a single nanoparticle nanofluid as compared to the hybrid nanofluid. This shows the efficiency of the hybrid nanofluid for a higher heat transfer amount and reduced skin friction. Figure 5e describes the influence of the thermophoresis parameter Nt on the    www.nature.com/scientificreports/

Conclusions
The current research investigates three-dimensional, rotating, incompressible, nanofluid, and hybrid nanofluid flow above the permeable stretchable surface for the heat and mass transmission rate. The governing equations are tackled at MATLAB through the bvp-4c algorithm after employing the similarity transformation. The major outcomes of the present study are stated below: 1. The velocity profiles decline with a rise in rotation, magnetic field, porosity, and increasing angle between the axis of rotation and the horizontal axis while increasing for mixed convection and stretching ratio parameter. 2. Prandtl number and thermal radiation, thermal slip, and stretching ratio parameters have decreasing effects on the temperature profile. 3. The concentration profile decreases under the increasing influence of thermophoresis and but has an increasing relation with the rotation parameter when it increases. 4. A maximum Nusselt and Sherwood number is noted in the absence of rotation and porosity of the medium. 5. The highest heat and mass transfer rate is noted when the inclined magnetic field and axis rotation are parallel to each other. 6. Thermophoresis and chemical reaction parameters increase the mass transfer rate when it increases. www.nature.com/scientificreports/ 7. Thermal radiation rapidly increases the heat and mass transmission when it increases. 8. Higher heat transfer rates and reduced skin friction are noted for hybrid nanofluid.  www.nature.com/scientificreports/

Data availability
All data generated or analyzed during this study are included in this published article.