Heat and mass flux through a Reiner–Rivlin nanofluid flow past a spinning stretching disc: Cattaneo–Christov model

The current work scrutinizes a non-Newtonian nanofluid free convective flow induced by a rotating stretchable disc. The examination surveys the Stefan blowing and Cattaneo–Christov mass and heat fluxes, as a precise illustrative model. The innovative aspects of the ongoing project include the analysis of the border sheet nanofluid flow near a revolving disc through thermophoresis, Reiner–Rivlin prototype features, and random nanoparticle motion. The Reiner–Rivlin non-Newtonian model is considered together with the effect of an unvarying axial magnetic strength. The constitutive formulae of a Reiner–Rivlin liquid have been reproduced in the cylindrical coordinates. Through implementing the applicable relationship transformations, the controlling partial differential equations are transferred to ordinary differential equations (ODE). This procedure yields a group of coupled nonlinear ordinary differential equations in relation to speed, heat, and nanoparticle concentration profiles that are impacted by several physical characteristics. These equations are analyzed by using the homotopy perturbation method (HPM). Due to the analytical solution given by HPM, the current work enables us to take the infinity of the layer as a parameter of the problem and discuss its variation in the obtained distributions. Consequently, a physical significant graphical visualization of the data is emphasized. The rates of mass and temperature transmission are examined to understand if any of the relevant parameters may improve these rates. Additionally, the Stefan blowing causes extra particles diffusion, which enhances heat transfer and raises the nanoparticles concentration and could be useful in some medical therapies. Furthermore, the stretching of the rotating disc is concluded, which improves the fluid heat transfer.


Heat flux relaxation time (ρc) f
Fluid heat capacity (ρc) p Nanoparticles heat capacity β 1 , β 2 Coefficients of heat and mass slip θ Non-dimensional heat ϕ Non-dimensional nanoparticle volume fraction γ Thermal slip factor 1 Energy diffusion Deborah number 2 Mass diffusion Deborah number α Rate of stretching parameter δ Slip factor of nanoparticle volume fraction The innovative potential in several biological sciences, biomedical (medicine) technological aspects, and non-Newtonian are the (Reiner-Rivlin) nanofluid mechanisms. Additionally, a variety of biological liquids including blood, synovial fluid, and saliva, exhibits a non-Newtonian behavior, and significant viscoelastic properties are demonstrated. Hayat et al. 1 investigated the hydromagnetic movement of a Reiner-Rivlin nano liquid over a revolving disc. In their work, the physical aspect of entropy optimization was addressed. Numerical analysis was done on the thermodynamic characteristics of the Reiner-Rivlin nano liquid movement caused by a rotating disc 2 . To examine how nanoparticles affect the thermodynamics of the Reiner-Rivlin nanomaterial, the nonhomogeneous two-phase nano-liquid model was studied. Numerous limiting instances were reported to be in strong accord with the bench marking investigations 3 . The drag force and the normalized force of a spherical bubble comprising a Reiner-Rivlin flow of a sheet of a surfactant film flowing in a Newtonian fluid were estimated 4 . The results of this framework can be employed in a wide range of engineering settings, in addition to chemical, biochemical, and environmental disciplines. Using the cell model technique, Jaiswal and Gupta 5 provided an examination of viscous incompressible liquid Stokes flowing through a swarm of immiscible Reiner-Rivlin fluid droplets in a cell. A Reiner-Rivlin viscoelastic flowing incorporating Cattaneo-Christov (CC) temperature distribution in a permeable material across a roughly rotating disc was examined 3 . Concentration and temperature profiles were improved as a result of Reiner-Rivlin material parameters. He et al. 6 looked at a planar boundary separating two bounded horizontal magnetic liquids, with the act of a uniform longitudinal magnetic strength.
The insinuation of the linear as well as nonlinear curves showed that the linear curve governs the stability zones. A well-known approximation for a viscous flow between two flat plates, one of which rotates, and the other is stationary, was taken into consideration 7 . The topic has drawn considerable interest, particularly for Navier-Stokes flow. Nevertheless, the description of a relatively non-Newtonian fluid as elucidated by Reiner-Rivlin theory seems to be of more significance. Practically speaking, organized convective cooling technologies involve a coolant imposed on a flat plate and can use these fluxes. This phenomenon also happens frequently in a wide variety of technical and manufacturing uses, including treating wastewater, revolving machinery, medical instruments, spinning blades, computer storage systems, etc. Additionally, all rounded spinning bodies moving in a fluid experience a forced flow caused by a rotating disc, hence the current analysis satisfies the indicated requirements to assess the precise analytical solution to real-world issues. Therefore, the goal of the current investigation is to broaden this research to non-Newtonian cases, while taking the Reiner-Rivlin fluid model into consideration. Nanofluids are fluids, including nanoparticles, which demonstrate considerable thermal improvement at low volume fractions. Many of the works on nanofluids aim to gain a better knowledge of their behavior in such a way as to employ them in such fields where heat transfer enhancement is critical, such as nuclear reactors, transport, technology, pharmacology, and foodstuffs, to mention but a few. It has also been suggested that nanofluids can be used as a smart liquid, permitting heat transmission to be controlled or increase suddenly. Wong  www.nature.com/scientificreports/ such requirements. In a few words, nanofluids are liquid concentrations of nanomaterials having a fundamental dimension of fewer than 100 nm. An investigation of the stagnation point flow of a Cu-water nanofluid across an extending and contracting layer was conducted 9 . The case when there is a delay before feeling the beginning of disruptions across the material is the subject of this study. It was claimed that meteorological and environmental research, biochemical engineering, the generation of power and mobility, the translation of solar energy, electronic properties, sensor microfluidics, tumbler in polymer manufacture, and other sectors might benefit from this phenomenon 10 . A mathematical framework of an MHD mixed convective Cu-water nanofluid border layer movement across a surface stretched surface was examined 11 . The innovative feature of this project is the use of a modified Buongiorno type to account for Brownian motion, thermophoresis, and volume fraction for nanofluids. An attempt was undertaken to theoretically explore the influence of volume fraction on a mixed convective Cu-water nanofluid flow across a stretching surface with activation energy and thermal radiation, due to the growing demand for regulated cooling systems 12 . The influences of mixed convection, thermal radiation, and chemical reaction were numerically investigated for the steady Cu-water nanofluid flow in the presence of a magnetic field 13 . Heat transmission in a tube with non-parallel sides was analyzed 14 . Cu and Ag nanoparticles were used as the base fluid, and water was utilized as a solvent. The thermocapillary flow qualities and heat exchange efficiency of nanofluid droplets on a range of surfaces were examined in order to identify the effects of the underlying structure of the nanofluid evaporating droplets 15 . Kumar et al. 16 used hybrid nanofluids to investigate the stability of heat transfer caused by large characteristics as well as the variation of fluids under the influence of thermal radiation. According to the graphical data, large amounts of the radiation factor result in more heat transfer, meanwhile, greater values of the Prandtl number result in heat transfer suppression. An overview of several nanofluid instability concerns as a function of different hydrodynamic and hydromagnetic characteristics was provided 17 . Nanofluid investigations are gaining more attention from scientists in recent days due to their multiple applications that make it possible for many industrial heat transfer processes. The field of medicine also uses nanofluids, some of which contain nanodrugs. Additionally, they are used in hard drives, heating and cooling equipment, jet engines, turbine systems, optical sensors, and hard discs. As a result, the current work is compatible with these fluids. The Cattaneo-Christov (CC) model successfully modified the Fourier model to include the crucial component of thermal and mass relaxation times. This demonstrates how a hyperbolic energy equation is created for a temperature field, allowing for the conduction of heat at a constrained velocity through the use of heat waves. Interesting real-world uses of this kind of heat transmission include nanofluidic flows and skin burns models. The CC with a varying entropy production and thermal relaxation time was studied 18 . The micropolar fluid that absorbs heat in the presence of, mixed convection and partial slip was examined. To examine the feature of heat and mass transfer, two different nanoparticles, a single-wall carbon nanotube and a multi-wall carbon nanotube, were taken into consideration. Using the effect of instigation energy with CC, an angled cylinder was inspected to investigate the stagnation place of the flowing of an organic molecule Powell-Eyring nanofluid 19 . The flow data was designed using the Buongiorno nanofluid framework, which took into consideration such nonlinear effects as infrared energy, binary chemical processes, and non-Fourier heat transfer. The movement, temperature distribution, and entropy formation of a non-Newtonian hybrid nanofluid through a lubricated surface with the concept of CC heat flux have not been yet studied, according to a thorough survey of the bibliography. The Casson fluid framework of the hybrid nanofluid through the lubricated surface with thermal radiation and viscous dissipation was therefore evaluated in a previous study 20 . The impact of the chemical reaction of the upperconnected Maxwell liquid was the focus of Khan et al. 21 . In such transmission, entropy minimizing in a hybrid nanofluid stagnation point flow across a nonlinearly extending layer in the presence of Thomson and Troian border requirements was considered 22 . The Darcy-Forchheimer relationship was established as a result of the porous medium. Using the CC heat flux framework, the nonlinear thermal propagation, heat production, and viscous dissipation were addressed. Heat and mass transfer processes were investigated using the CC methodology. Using mathematical formalism, Raju et al. 23 explored the heat movement transmission parameters of a Maxwell liquid over a stretched surface with CC, heat source or sink, and suction/injection were surveyed. Subbarayudu et al. 24 used initiating energies with a twofold chemical process across a stretched sheet to construct a conceptual model to examine the hydro-magnetic movement of carbon nanotubes immersed in a Maxwell nanofluid. Heat transfer processes in energy expression were investigated using both the dynamic absorption and nonlinear heat flux models through an original approach, as well as the CC prototype of temperature diffusion, which is a more advanced version than Fourier's heat flux formula. Shankar et al. 25 expounded the magnetic effects of CC double diffusion patterns on heat and mass transmission behavior of a viscous incompressible, time-dependent, two-dimensional Casson nanofluid movement across the channel with Joule heating and viscous dissipation impacts. The validity of the findings and the approaches employed to accomplish the goal of the present work was clearly demonstrated by a comparison between the results provided and the findings of earlier research. In light of the advantages of employing the CC model, the current problem is performed along with this approach.
The occurrence of ODEs/PDEs plays an essential role in several zones of science, physics, engineering, and ecology. The physical significance of the dynamical behaviors requires understanding the ODEs/PDEs. Many scientists have investigated nonlinear phenomena significantly in recent decades. In engineering implementations, the nonlinear behavior is more frequent than the linear one. Although linear systems can yield analytical solutions simply, only limited particular nonlinear schemes have analytical solutions. The perturbation theory is the first approach to yield approximate analytical nonlinear system solutions. As a new methodology of the perturbation theories, He 26 is accredited as being the first to solve an ODE by introducing a synthetic incorporated parameter q ∈ [0, 1] into the differential equation. Along with this parameter, the problem may be divided into two parts. Consequently, the artificial parameter divides the previous two parts. Therefore, He [26][27][28][29] proposed a prospective and effective analytical method in this approach, which is abbreviated as (HPM). This technique is used to provide an analytical approximate solution for a variety of ODEs encountered in science www.nature.com/scientificreports/ and mathematics. When the HPM is compared to other perturbation methods, it is clear that it produces more accurate results. Therefore, the HPM has all of the advantages of the classical perturbation approaches, with no need for small parameter assumptions in the solution procedure. The approach has a simpler structure since it needs less calculation time and has rapid computational accuracy than the preceding classical techniques.
To solve any problem, this method just needs primary circumstances in such a way as to produce an analytical approximate solution as an infinite power series. HPM is often used to solve NODEs by a number of researchers whose findings are easily presented. One of them is referred to as the linear part which must have an exact solution, whereas the other is known as the nonlinear part. The nonlinear EHD stability of a cylinder-shaped interface between two Walters' B-type viscoelastic liquids was examined in a permeable media 30 . The extended frequency theory was merged with the HPM to produce a distribution of the surface displacement in an analytical periodical approximation. A motile microbe behaving in an MHD movement of an incompressible nano liquid that followed the non-Newtonian Jeffrey prototype was investigated 31 . Using HPM, the fundamental equations of motion are analytically resolved. The non-Newtonian flow is followed by an incompressible nanofluid flow. The Casson prototype was employed for its non-Newtonian fluid behavior 32 . The flow fills up the conical space among the cone stationary, the revolving surfaces, and the horizontal disc. The HPM was used to analytically solve these equations. Owing to the highly nonlinear fundamental equations in the present study, the effective HPM is used to analyze the problem at hand. In a conclusion, the HPM can be used to investigate a wide range of additional nonlinear issues of scientific and technological significance, especially in linear and nonlinear wave mechanical and vibrational issues. The current study investigates the Reiner-Rivlin fluid model induced by a rotating stretchy disc with Stefan blowing and the CC modeling, in addition to the above-mentioned characteristics. It is considered that there is a constant magnetic field. The controlling equations of motion are converted into a new group of ODEs using appropriate similarity modifications. Food manufacturing, heat transfer absorption throughout paper drying operations, and separating processes are affected by the main challenge of the present paper.
The answers to the following relevant queries are anticipated at the conclusion of this investigation: • What are the mechanisms of heat and mass transmission through nanofluids under the effects of Stefan blowing and the CC model? • Does the Reiner-Rivlin model cause major deviations in the thermal boundary of a rotating disc?
• Does the expanding nature of the revolving disc enhance the flow, heat, and nanoparticle transfer rates?
• What are the effects of thermophoresis and Brownian motion features on the flow?
The remainder of the paper is structured as follows: The approach taken to this issue is outlined in section "Problem structure". The essential boundary conditions are included in this section together with the governing equations of motion. In section "Straightforward changes of correlation", the boundary-value problem is covered utilizing simple similarity transformations. How to use the HPM to develop the analytical solution is described in section "Methodology of solution". The conclusions and comments are introduced in section "Results and discussion". The main findings are subsequently reported as closing remarks in section "Conclusions".

Problem structure
The configuration addresses a two-dimensional MHD non-Newtonian nanofluid obeying a Reiner-Rivlin representation, flowing past a rotating stretching disc. The model utilizes cylindrical coordinates (r, θ , z) for ease of use. The disc is located at z = 0 and spins with a uniform angular velocity across the axial direction, while it stretches in the radial path with a stable stretched parameter a. It is planned to provide a magnetic field B 0 parallel to the z-axis that is strong and consistent. The revolving surface is excited by a warmer fluid of temperature T r and meditation C r , whilst in free stream; the surrounding fluid is kept at constant values of temperature T ∞ and nanoparticles C ∞ . Additionally, along with the normal velocity condition, the flat surface is disposed to border circumstances, having mass concentration with the Stefan blowing effect. In accordance with the abovementioned construction, the configuration model is clarified and shown in Fig. 1a. www.nature.com/scientificreports/ The spinning disc reactor (SDR), shown in Fig. 1b, is one of the most significant applications of our model. Producing nanoparticles with regulated qualities, such as particle size and distribution, is easy with this equipment. Over the past few decades, this technique has been successfully established. The velocity of the disc, its surface type (smooth or grooved), the feed liquid flow rates, the number of feeding points and their positions, and the presumptive stabilization manager are the main SDR control strictures. The liquid combinations are concentrated in the disc center and rotate at high speeds (between 300 and 3000 rpm), forming a thin fluid layer (1-200 nm). Large heat and mass transmission is induced by the fluid layer thickness and the wide contact area it has with the disc surface. Finally, highly quick and efficient mixing (micro-mixing) of the components of the liquid flowing over the disc surface is made possible by the drag forces between the fluid flowing and the disc surface 33 .
Controlling equations of motion. The Cauchy stress tensor τ is used for the Reiner-Rivlin fluid 34 and is defined as follows: where A = ∇V + ∇V T is the first order Rivlin-Ericksen tensor, V represents the fluid velocity and I is the identity tensor.
The conservation of mass and momentum equations of an incompressible non-Newtonian liquid may be described as follows 2,35 : and where ρ f is the fluid density, J is the current density and B is the magnetic strength.
The temperature and nanoparticles dissipation equations 25,36 of fluid may be written as: and where T is the fluid heat, Q is the thermal flux, C is nanoparticles concentration and j is the mass flux. The CC diffusion theory 37,38 is employed to analyse the heat and mass flow performance of nano-fluids instead of the conventional transport models. To investigate this impact, the upper-convective material time derivative of a vector F has the form: In terms of CC representation, the Fourier's law for both the thermal flux Q and the mass change j have been generalized as follows: where E is the heat flux relaxation time, α m is the thermal conductivity, C is the mass flux relaxation time, and D B is the Brownian diffusivity factor. On the other hand, the standard Fourier and Fick laws may be required in return by putting the coefficients E = C = 0 in the above equations, correspondingly.
Utilizing the situation of incompressibility ∇ · V = 0 and the situation of stationary flow that is ∂Q ∂t = 0 and ∂j ∂t = 0 , the Eqs. (7) and (8) change over to By reducing Q and j in Eqs. (4) and (5), and using Eqs. (9) and (10), the heat and nanoparticles equations with the Brownian and thermophoresis properties become www.nature.com/scientificreports/ and such that E and C are the modified expressions for CC theory. In light of the above-mentioned description, these expressions 35 are defined as follows: Simultaneously, the governing equations of motion concerning continuity and momentum of the Reiner-Rivlin fluid may be formulated as follows: where the mathematical expression of the components τ ij of the Reiner-Rivlin fluid may be listed as follows: and The temperature equation then becomes Additionally, the nanoparticles concentration equation gives www.nature.com/scientificreports/ The aforementioned main equations of motion should be valid to the suitable border situations. The suitable border situations were given by 35 . These circumstances may be listed as follows: where β 1 and β 2 are the slip coefficients of heat and mass, respectively. Now, the boundary-value problem of this manuscript becomes well-defined. The method of solution will be outlined in section "Methodology of solution".
Supplementary amounts of attention. The skin friction factor C f , the Nusselt numeral Nu and the Sherwood numeral Sh are among the most important physical amounts that are related to the topic. These amounts are mathematically well-described as: and From a physical standpoint: • The skin friction factor C f is produced by the viscosity of fluids. It has been recognized as the laminar drag in the turbulent drag. Skin friction drag is usually connected to families of the Reynolds numeral, which is the relationship between the inertial force and the viscous force. • The Nusselt number Nu is the proportion between convection and thermal conduction over a boundary.
• Finally, the Sherwood number Sh represents the proportion between mass transmission by convection and mass transfer by dispersion.

Straightforward changes of correlation
The foremost nonlinear partial equations are simply converted into ODEs through an appropriate similarity conversion. In view of Refs. 2,35 , the appropriate similarity transformations may be expressed as: where the dash signifies the differentiation with reference to the independent parameter η. Inserting these transformations into the aforementioned governing equations of motion, one finds and (26) www.nature.com/scientificreports/ where Re = √ �/υ r is the Reynolds numeral. Moreover, the border restrictions as given in Eq. (27) have been transformed to: where, η ∞ represents the infinity value.

Methodology of solution
As known, the main aim of all perturbation methods is transforming the NDEs into normal equations. Because the HPM 39,40 is significant, effective, and promising, it will be adopted to analyze the fundamental organization of NDEs given by Eqs. (32)(33)(34)(35) along with the border restrictions (39). The differential equation will be divided into linear and nonlinear portions. The two branches are disconnected by the embedded limitation, specifically q ∈ [0, 1] to construct what is recognized by the Homotopy equation. Therefore, the above-mentioned equations may be formulated as follows and L 4 ≡ d 2 dη 2 are the linear operators. According to the aforementioned technique, the dependent variable may be expanded as follows: where χ(η, q) stands for any of the functions f ′ , g, θ and ϕ.
Replacing Eq. (44) into Eqs. (40)-(43), after explanation and comparison of the similar powers of q-terms, one gets: Zero-order system. and with the suitable border circumstances given as:  The constants c 10 , c 11 , . . . , c 36 will be moved to the Appendix to follow the manuscript easily. Lastly, in view of the HPM, the estimated solution of the non-dimensional velocity components, temperature and nanoparticles concentration may be written as follows: Consequently, the profiles of these functions can be represented as: (59) (63) www.nature.com/scientificreports/ and The following section seeks to explore the impacts of the parameter settings in the problem at hand on the various distributions of nanoparticles velocity, heat, and concentrations, with some diagrams to further demonstrate these influences.

Results and discussion
Under the influence of Stefan blowing, the stationary two-dimensional MHD border layer flow of Reiner-Rivlin nanoparticles over a revolving extending warmed disc is explored. The HPM is utilized to analyze the nondimensional ODEs (32)- (35) with border restrictions (39).
In order to elucidate the problem physically, the results are discussed to demonstrate the effects of numerous restrictions. These factors include the Reiner-Rivlin fluid factor K , the Brownian movement limit N B , the thermophoresis limitation N T , the rate of the stretching constraint α , the nanoparticle volume fraction slip parameter δ, the Schmidt numeral S c and the Prandtl number Pr . The study at hand concentrates on the impacts of these limitations on the velocity, temperature, nanoparticles profiles, as well as skin-friction C f , Nusselt Nu and Sherwood Sh quantities. These distributions are plotted in accordance with the data mentioned in Figs. 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33. For more practicality, the above-mentioned distributions are graphed along with the non-dimensional parameter η for some standard values of the above factors, which change according to the factor illustrated in every diagram, as follows: Speed profile. The speed profile is of a considerable interest when examining fluid motion, where it gives many details about the flow field. Therefore, Figs. 2, 3, 4, 5 are plotted to display the changes of the radial as well as azimuthal velocities versus the axis of the disk for various standards of the substantial limitation of the Reiner-Rivlin fluid K and the magnetic factor M . It is found that the radial and azimuthal velocities increase near the disk and attain maximum values at a short distance away, and then they decrease as they move far away from it to vanish at the free stream. In details, it is understood from Figs. 2 and 3 that the radial and azimuthal velocities decrease with the rise in the value of K . Physically, the non-Newtonian factor K , as an amount of an additional     www.nature.com/scientificreports/ Figures 6 and 7 display the profiles of the radial and azimuthal velocities for diverse standards of the rate of the stretching factor α . Figure 6 shows that there is a dual nature in the radial velocity f ′ (η) with the rise of α , where f ′ (η) is enhanced at small values of η ; whereas it decreases at large values of η . Naturally, nearest the disc, the rate of the stretching parameter improves the radial velocity due to the dependence of stretching of the surface on the radius r . After that, this effect is reversed when moving away from the disc. This result corresponds   www.nature.com/scientificreports/ to that previously described 35,42,43 . By contrast, since the azimuthal velocity takes the tangential direction, this dual effect does not occur. Therefore, the growth in the rate of the extending parameter yields a decrease in the azimuthal speed as seen from Fig. 7. This is similar with the preceding depicted findings 35,43 . Figures 8 and 9 demonstrate the influence of the Stefan blowing limitation S w on the radial and tangential velocities. It is concluded that as the Stefan blowing limitation S w grows, both the radial velocity as well as the tangential velocity increase. Physically, this is due to the fact that the blowing effect causes more strong movement of the fluid, and eventually a great amount of the fluid moves faster in the radial and tangential directions sufficiently   www.nature.com/scientificreports/ away from the disc. This finding is consistent with that previously illustrated 35,42,44 . Lattif et al. 42 contended that the mutual impact of the rising blowing factor and the stretching stricture cause the radial speed distribution to rise. As a conclusion, the Stefan blowing effect outweighs that of the rate of the stretching parameter. Numerous works related to the border sheet of non-Newtonian liquid flows, as given in Refs. 2,35,41,42,45 , numerically investigated all the physical distributions. These works couldn't take infinity greater than 10 in the best   www.nature.com/scientificreports/ cases. Fortunately, in the analytical solutions of the present study, we could take it as larger values η ∞ , where η ∞ is an arbitrary number. Consequently, the current study scrutinizes the effect of η ∞ on all the previous physical distributions. Furthermore, the improvement in η ∞ supports the accidental movement of the nanoparticles, and henceforth it improves the radial motion and reduces the tangential motion. Accordingly, the growth of f ′ (η)   Temperature distribution. As shown in the methodology of the problem, the current work investigates the MHD Reiner-Rivlin nanofluid. Accordingly, this mechanism has enhanced the heat transfer rate. Therefore,   Just as before, the magnetic strength produces a Lorentz force that contradicts movement, which in turn improves the thermal diffusion as found in Fig. 12. This conclusion is consistent with Mabood et al. 35 and Alebraheem and Ramzan 41 . In addition, the Brownian movement factor N B , as a measure of the random motion of the nanoparticles, raises the temperature   www.nature.com/scientificreports/ in the neighborhood of the boundary layer as observed by Fig. 13. Physically, the random or accidental motion of particles is a logic reason for the rise of heat of the fluid because of the growth of the friction force between particles. This finding is in accord with the earlier result 46,47 . Figures 14 and 15 show the enhancement in the temperature due to the growth in the Stefan blowing factor S w and the thermophoresis factor N T . Actually, the growth of the blowing parameter causes more collision between      www.nature.com/scientificreports/ Figures 16 and 17 show the relation between heat and both the stretching factor α and the thermal slip factor γ , correspondingly. It is revealed that with the rise of α and γ a reduced effect on the thermal dispersal is found. Really, the increase of α leads to a wide dispersion of heat since the extending superficial of the ambient fluid develops, which brings about a reduction in the liquid heat as depicted in Fig. 16. Figure 17 shows the same performance of γ , where the thermal slip parameter is a measure of thermal change at the edge. Therefore, it   www.nature.com/scientificreports/ leads to a reduction in the liquid heat from the disc temperature owing to the existence of thermal opposition. These conclusions are parallel to those obtained by Sabu et al. 2 and Mabood et al. 35 . Figures 18 and 19 show the reduction of heat transfer in light of the increase of both the Deborah numeral owing to heat dispersion 1 and the Schmidt numeral S c . Physically, materials that have low relaxation times flow easily, which means that the thermal transfer in the fluid with the increase of the Deborah number 1 needs additional time to transmit temperature in the vicinity of liquid elements. Therefore, heat reduces as 1 increases, as shown in Fig. 18. This result agrees with that earlier obtained 3,35,43 . Finally, the increase of the Schmidt numeral S c , as the ratio between momentum and mass diffusivity, means a reduction in the molecular diffusivity and the random motion of molecules, which yields a decrease in the liquid temperature, as indicated by Fig. 19. The combined effects of those two parameters cause a reduction in the temperature profile as stated by Lattif et al. 42 .

Nanoparticles volume fraction distribution.
In what follows, an interesting phenomenon concerning the nanoparticles volume fraction ϕ is displayed throughout Figs. 20, 21, 22 and 23 to illustrate the effects of the limitations: the rate stretching parameter α , the nanoparticles volume fraction slip parameter δ , the Deborah numeral owing to the mass dispersion 2 and the Schmidt numeral S c , correspondingly. As shown in Fig. 20, the distribution ϕ decreases with the increase of the value of the parameter α . Physically, as α rises, an extensive mass diffuses from the stretching sheet into the surrounding liquid. Henceforth, it produces nanoparticles scattering away from the boundary surface, and brings about a reduction of the nanoparticles concentration. The growth in the factor δ provides a objective justification to the decrease in ϕ , where the nanoparticles accelerate in their random flow with the increase of δ , as demonstrated in Fig. 21. Figure 22 shows the impact of the factor 2 on the ϕ profile. It is found that the profile of ϕ decreases through the rise of 2 . Actually, the increase of 2 means more time to diffuse mass to the vicinity fluid; therefore, the nanoparticles volume fraction ϕ minimizes through the flow. The increase of Schmidt numeral S c reduces the nanoparticles volume fraction, as displayed in Fig. 23, owing to the fact that the Schmidt numeral is the proportion of momentum to mass diffusivity. Consequently, mass diffusivity diminishes for large values of S c , which requires a reduction of ϕ . All of these findings are consistent with those of 2,35 .
Finally, Figs. 24 and 25 show the impact of both parameters M and S w on the nanoparticles volume fraction ϕ . Figure 24 indicates that ϕ rises with the rise of the blowing factor S w . This occurs owing to the growing of particles diffusion by the blowing factor, which means an enhancement of ϕ . This outcome is in accord with Mabood et al. 35 and Beg et al. 46 . In addition, Fig. 25 indicates that the ϕ distribution has more increasing rate with the induced magnetic parameter M . The speed profile decreases with the intensity magnetic force, and this slow motion leads to more accumulation of nanoparticles and hence an increase in ϕ . The same behavior was observed by Mabood et al. 35 . Figures 26,27,28,29,30,31,32 and 33 discuss the influences of diverse factors on the skin friction constant C f , Nusselt Nu and Sherwood Sh numbers. Figures 26 and 27 display the variation of C f versus M . It is noticed that the skin friction grows through the rise of M . Additionally, Fig. 26 illustrates that the skin friction grows with the growth of the Schmidt numeral S c and drops slightly with the Stefan blowing coefficient. Figure 27 indicates that the skin friction grows through the rise of both the stretching factor α and the material factor K . These results are physically expected, owing to the appearances of the intersecting velocity gradients of the Reiner-Rivlin nanofluid sheets, which are produced by the enhanced cross-viscosity. These results agree with those obtained earlier 35,42,43 . Figures 28, 29, 30 and 31 illustrate the behavior of the Nusselt number Nu versus the Prandtl number Pr . It is concluded that the influence of Pr fluctuates in accordance with the values of the other parameters. As observed from all these graphs, the Nusselt number drops by the rise of the Stefan blowing coefficient S w , the Brownian movement factor N B , the thermophoresis factor N T and M ; whereas it rises with the rise of the Schmidt numeral S c , the stretching parameter α , the slip parameter δ and Deborah number due to energy diffusion 1 . www.nature.com/scientificreports/ Finally, the behavior of the Sherwood number Sh versus the magnetic parameter M is illustrated throughout Figs. 32 and 33. It is found that the Sherwood numeral Sh decreases with the rise of the magnetic factor M , the Stefan numeral S w and the Schmidt numeral S c ; whereas it improves with the growth of the stretching factor α and the Deborah numeral owing to the mass dispersion 2 . It is worth mentioning that most of these results are consistent with those observed by Mabood et al. 35 and Lattif et al. 42 .

Conclusions
The foremost purpose of this work is to scrutinize the straightforward MHD movement of a non-Newtonian laminar hydrodynamic two-dimensional flow that follows the Reiner-Rivlin nanofluid. The cylindrical coordinates are adequate for describing the problem configuration. The border sheet movement lies down on a revolving disc through a stretched parameter in its radial direction. The results of this model may be valuable in many industrial and engineering applications like SDR as mentioned in section "Problem structure". The physical phenomenon follows the standard approximations of the boundary layer flow. The plate surface does not promote mass concentration and the Stefan blowing effect due to boundary circumstances. The CC heat and mass diffusion concept has also been exploited. The fundamental partial differential equations are transformed into conventional standard ones via correspondence transformations. Several non-dimensional physical numbers are gained. The analysis is performed in light of the HPM to create an estimated distribution of all typical relevant functions concerning velocity, temperature, and concentration, which simplifies the mathematical manipulation. The main findings of the work may be summarized in the following points: • In accordance with the profiles of the velocity components, it is observed that the rise of the factor S w raises the distribution of these components. Simultaneously, the increase of the parameters M, α and K reduces the velocity distributions. • Heat transmission is examined to specify the impact of the revealed parameters on this distribution. It is found that the increase of the parameters α, 1 , and S c decreases the energy distribution. On the other hand, the increase of the parameters M, N B , S w and N T improves heat transfer. This means that the usage of magnetized nanofluids has an effective role in the applications that need more improvement in heat exchange. • The nanoparticles volume fraction is investigated regarding the various physical parameters. It is shown that the increase of the parameters S w and M increases this distribution. Meanwhile, the increase of the parameters α, δ, 2 and S c decreases it. • The skin friction, Nusselt and Sherwood numerals are illustrated graphically with the variation of the relevant parameters.