key: cord-0019559-3djyykt7
authors: Yu, Bai; Ramzan, Muhammad; Riasat, Saima; Kadry, Seifedine; Chu, Yu-Ming; Malik, M. Y.
title: Impact of autocatalytic chemical reaction in an Ostwald-de-Waele nanofluid flow past a rotating disk with heterogeneous catalysis
date: 2021-07-30
journal: Sci Rep
DOI: 10.1038/s41598-021-94918-7
sha: 17d4cc5e690b824f7e884d77505336a32ed783a1
doc_id: 19559
cord_uid: 3djyykt7

The nanofluids owing to their alluring attributes like enhanced thermal conductivity and better heat transfer characteristics have a vast variety of applications ranging from space technology to nuclear reactors etc. The present study highlights the Ostwald-de-Waele nanofluid flow past a rotating disk of variable thickness in a porous medium with a melting heat transfer phenomenon. The surface catalyzed reaction is added to the homogeneous-heterogeneous reaction that triggers the rate of the chemical reaction. The added feature of the variable thermal conductivity and the viscosity instead of their constant values also boosts the novelty of the undertaken problem. The modeled problem is erected in the form of a system of partial differential equations. Engaging similarity transformation, the set of ordinary differential equations are obtained. The coupled equations are numerically solved by using the bvp4c built-in MATLAB function. The drag coefficient and Nusselt number are plotted for arising parameters. The results revealed that increasing surface catalyzed parameter causes a decline in thermal profile more efficiently. Further, the power-law index is more influential than the variable thickness disk index. The numerical results show that variations in dimensionless thickness coefficient do not make any effect. However, increasing power-law index causing an upsurge in radial, axial, tangential, velocities, and thermal profile.

The customary fluids like oil and water are poor heat transfer liquids as they hold low thermal conductivities. As the role of thermal conductivity is vital in heat transfer processes, therefore, a variety of methods are devised to enhance the thermal conductivity. The efficient technique for the enhancement of thermal conductivity is to insert the nanoparticles into the base liquids. The nanoparticles may be from metals, oxides, nitrides, or carbides. The nanofluids with unique characteristics have wide-ranging applications including fuel cells, hybrid-powered engines, and pharmaceuticals, etc. The novel notion of nanofluids is introduced by Choi 1 in 1995. Later, many investigations have been reported [2] [3] [4] [5] [6] . Lately, Ram and Kumar 7 analyzed the fluid flow past a rotating disk with viscosity as a function of temperature. Rashidi et al. 8 investigated the fluid flow past a rotating spongy disk by using numerical method. They focused on the temperature-dependent viscosity, density, and thermal conductivity. Sheikholeslami et al. 9 examined the three-dimensional condensation nanofluid thin-film flow by employing the effects of thermophoresis and Brownian motion with normalized thickness over an inclined rotating disk. Bachok et al. 10 considered the Maxwell-Garnett model and the Patel model to study the impact of effective thermal conductivity of nanofluid flow past a rotating porous disk. Kendoush 11 obtained the similarity solution to visualize the heat transfer rate for rotational and flow Reynolds number. Turkyilmazoglu 12 performed the boundary layer flow analysis by considering copper and Aluminum oxide. The main conclusion of their finding was that copper has much heat transfer rate than Aluminum oxide. The study of heat augmentation by employing nanofluid has www.nature.com/scientificreports/ been the main concern of researchers 13-17 . Hayat et al. 18 analyzed the increasing effect of heat transfer rate for elevating nanoparticle volume fraction for MHD nanofluid flow due to rotating disk. Various researchers have put much effort into identifying the beneficial properties of nanofluid [19] [20] [21] [22] [23] [24] .

Presently the researchers are more interested in studying non-Newtonian fluids for their abundant industrial and engineering applications in comparison to the Newtonian fluid. Therefore, the class of Non-Newtonian fluid has gained considerable attention from many researchers. Mitschka et al. 25, 26 analyzed the power-law Non-Newtonian fluids to derive the relation for the frictional resistance of the disk. Andreson et al. 27 investigated that boundary layer thickness decreases for diminishing values of power-law index. They eliminated the ambiguity in the previous results obtained for shear-thinning fluids by obtaining the results for shear thickening fluids. Attia 28 studied the heat transfer process of unsteady Reiner-Rivlin flow past a rotating disk. Sahoo 29 studied the Von Karman flow of non-Newtonian fluid. Ahmadpour and Sadeghy 30 discussed the swirling flow of Bingham fluid. They studied the impact of yield stress on boundary layer thickness and volumetric flow rate. Griffiths 31 studied the power-law and Bingham fluid due to rotating disks. Griffiths et al. 32 studied the neutral curve for power-law fluids on a rotating disk. Lin et al. 33 considered the convective heat transfer phenomenon in powerlaw fluids along with an inclined plate. They concluded that the heat transfer process is highly dependent on the power-law exponent. Following their models, more thermal conductivity models were addressed 34, 35 . Ming et al. 36 evaluated the steady heat transfer phenomenon of power-law fluid past a rotating disk. They studied the significant impacts of the power-law index on thermal and radial velocity profiles.

Due to abundant civil, aeronautical, mechanical applications, the deformable or elements of variable thickness have been taken into consideration. Shufrin 37 analyzed the stability of deformable plates. To improve the utilization ability and to reduce the weight of structural elements, the plate elements have been considered. Fang et al. 38 explored the sheet of variable thickness with power-law surface velocity. Stretching sheets of variable thickness have been reexamined by obtaining the dual solution in a thermal diffusive flow 39 . Hayat et al. 40 investigated the same problem by considering the Cattaneo-Christov heat flux model. Wahed et al. 41 deliberated the research on nanofluid flow with non-linear velocity over a moving surface of varying thickness. Li et al. 42 studied the nanofluid flow past a rotating disk of variable thickness for power-law fluid. Xun et al. 43 studied the Ostwald-de Waele fluid past a rotating disk of variable thickness. Nanofluid flow past a rotating disk with variable thickness influenced by melting heat transfer is examined by Hayat et al. 44 .

The homogeneous-heterogeneous reaction is executed in the presence of a catalyst. The rate of the chemical reaction is enhanced in the attendance of a catalyst. However, in the presence of surface catalyzed chemical reaction the rate of reaction is more accelerated. Chaudhary and Merkin 45 were the pioneers who developed the model for homogeneous-heterogeneous reactions. Recently homogeneous heterogeneous in the disk problem has been studied copiously 46, 47 . Liu et al. 48 studied the impact of the surface catalyzed parameter by considering porous media. Hayat et al. 49 studied the Darcy-Forchheimer nanofluid flow with homogeneous heterogeneous reactions in the presence of carbon nanotubes.

Given the foregoing, it is witnessed that abundant studies are accessible in the literature focusing on non-Newtonian nanofluid flows over rotating disks. This geometry even becomes narrower if we talk about the flow of non-Newtonian nanofluid flow over a rotating disk with variable thickness in a permeable media with homogeneous-heterogeneous chemical reactions. But here in this study, the novelty is manifold including the surface catalyzed chemical reaction in addition to the homogeneous-heterogeneous chemical reactions. Secondly, the consideration of variable thermal conductivity and the viscosity instead of the constant and the non-Newtonian fluid Ostwald-de-Wale blended with melting heat boundary condition. The problem is solved numerically employing the bvp4c MATLAB built-in function. The salient outcomes of the model are discussed via graphs and tables. A comparison with a previously published paper is also added and an excellent concurrence is achieved in this regard.

The goal of this research is to answer the subsequent pertinent questions:

i. How surface catalyzed parameter influence the concentration profile? ii. What is the association of the power-law index and variable viscosity? iii. Which is more influential on the thickness of the boundary layer either the power-law index or variable thickness index? iv. How do velocity and thermal profiles are influenced by disk thickness index? v. How does the Schmidt number affect the concentration profile of chemical species for pseudoplastic fluid? vi. How does the heat transfer process is affected by varying Prandtl number?

The flow under consideration assumes the axisymmetric, laminar, steady flow of Ostwald-de Waele fluid in porous media driven by the rotation of disk of variable thickness with angular velocity along z-direction with homogeneous-heterogeneous reactions. The heterogeneous reaction occurring on the surface of porous media is termed a surface catalyzed chemical reaction. Melting surface temperature is maintained at T m on the disk, while T ∞ is representing the ambient fluid's temperature (Fig. 1) .

The homogeneous chemical reaction governed by the isothermal cubic autocatalytic reaction is written by the following chemical equation 45, 46 While the heterogeneous reaction proceeds on the catalyst's surface and fluid-solid interface and its order are one. The chemical equation for this reaction is given by www.nature.com/scientificreports/ here A * and B * are representing chemical species involved in a chemical reaction. If we consider the concentration of the chemical species A * and B * is a and b respectively, then the rate of homogeneous reaction can be determined by using the following equation 40 The rate at which surface catalyzed chemical reaction proceeds can be estimated by the following equation 48 here A * and B * have the diffusion coefficient given by D A * and D B * respectively. n is the normal vector of unit magnitude in the fluid direction.

The rate of reaction occurring in porous media is given by 48 where the porous media has the interfacial surface area S. The flow geometry is analyzed by considering the cylindrical coordinates (r, ϕ, z) under the assumption of With boundary conditions are 44 consistency coefficient respectively. n is the power-law index. For n = 1 , we have Newtonian fluid. 0 < n < 1 corresponds to pseudo-plastic fluid while for dilatant fluid we have n > 1 . We shall define the dimensionless radius by r * = r/R o . Assume the following similarity transformation:

Consider deformations as:

The transformed partial differential equation including continuity and momentum equations are

Assuming the diffusion coefficient of chemical species same we take ψ + φ = 1 and δ = 1 , we obtain, For simplicity we denote A = r n−1 (1 + r) m(n+1) , where here α is the dimensionless thickness coefficient of the disk, Re is the Reynolds number, δ is the quotient of diffusion coefficients,K 1 and K 2 are the measurement of the strength of homogeneous and heterogeneous reaction respectively. S ν parameter of interfacial area, K, k vs and Ma are the porosity, surface catalyzed parameter, and melting heat parameter. Pr is the Prandtl number. Moreover, radial and tangential shear stress can be estimated by the following equation

Heat flux is defined as

The system of Eqs. (18)-(24) is transformed to a system of first-order differential equations and solved using the MATLAB software function bvp4c. To get a numerical solution, a tolerance of 10 -6 is set for the initial approximations. This presumed prior guess must satisfy Eq. (25) . The transformed coupled non-linear ordinary differential equations are computed using the bvp4c technique. To get the system of first-order equations, new variables are first introduced:

On applying the transformation we obtained the following system of differential equations as. 

This section studies the axial, radial tangential, thermal, and concentration profile for arising pertinent parameters with disk thickness index and power-law index. Figure 2 represents the radial velocity profile for increasing power-law index by keeping the thickness index of the disk constant. Non-Newtonian rheology modestly affects the radial velocity profile for small values of η . For pseudo-plastic fluids, the peak values are attained for the radial velocity profile. Increasing pseudo-plasticity indicates that the accuracy of boundary layer approximations deteriorates. In fact, increasing the power-law index leads to monotonic thickening of boundary layers. Rotating disks exert the centrifugal force in the outward direction causing the outward radial flow and inward axial flow. Moreover, peak values tend to reduce as we check for dilatant fluids. Figure 3 is the depiction of decreasing shear-driven motion by fixing the thickness index of the disk m = 0.8 and increasing the power-law index. The non-Newtonian rheology is exhibited by varying power-law index which causes the fluids parallel to the disk to perforate otherwise stagnant fluid. Figure 4 represents the graphical behavior axial velocity profile for increasing power-law index to study the non-Newtonian rheology. The decline in axial velocity is due to the centrifugal force causing the axial inflow. The viscosity function is dependent upon viscous consistency coefficient. Increasing power-law index results in the increase of variable viscosity, which in turn responsible for the thickening of the boundary layer. And hence the significant enhancement of the flow occurs in the axial direction to compensate the radial outflow. Figure 5 exhibits the decrease of the thermal profile from shear-thinning www.nature.com/scientificreports/ fluid to shear-thickening fluid. Increasing the power-law index by fixing disk thickness m = 0.8 will cause the reduction in thermal boundary layer thickness. Which causes the upsurge in heat transfer efficiency. Figure 6 is drawn for radial velocity profile of dilatant and pseudoplastic fluid by varying thickness of the disk. The graphical analysis of non-Newtonian rheology reveals that the power-law index is more influential than the thickness index of the disk. In the regime of the disk, the radial velocity upshoot occurs increasing the thickness index of disk. In the boundary layer region, the shear stress is firstly positive and increasing and reaches to its maximum values and then it falls, causing the radial velocity profile to go to zero. Increasing the thickness index of disk, the boundary layer gets thicker. For pseudo-plasticity, the peak value is small as compared to the dilatant fluids. Figure 7 is drawn for tangential velocity profile for pseudoplastic and dilatant fluids by varying index thickness of the disk. The tangential velocity profile escalates by increasing the thickness index of disk. Figure 8 represents the axial inflow by varying wall thickness parameter for dilatant and pseudoplastic fluids. The power-law index is more influential as compared to the wall thickness index. Figure 9 is the depiction of the thermal profile for varying thickness index of the disk. The trend for increasing disk thickness index is again increasing. The reason for up shooting in the tangential, axial, and thermal profile is that the velocity and thermal boundary layer thickness enhances for increasing disk index thickness. Figure 10 exhibits the thermal profile for escalating Pr . The increase in thermal profile causes the increase in thermal conductivity which leads to more surface heat www.nature.com/scientificreports/ transfer rate. Resultantly thermal boundary layer thickness gets thicker. Figure 11 depicts the concentration profile for homogeneous reaction parameter and surface catalyzed parameter for shear thickening fluids. The enhancement in surface catalyzed and homogeneous reaction parameter causes the thinning of the concentration boundary layer for non-Newtonian rheology. The reactants consume during homogeneous reaction proceeds. By increasing the homogeneous reaction parameter, the reaction proceeds more efficiently in the presence of porous media with surface catalyzed reaction. Hence, a decline in concentration profile is seen. Figure 12 displays the impact of the melting heat parameter on the axial velocity profile. The escalating axial velocity profile indicates that upon increasing melting heat parameter, molecular motion is enhanced. Figure 13 is sketched for a thermal profile by varying Brownian motion parameter for pseudoplastic fluid. The physical significance of enlargement in the thermal profile is the heating up of the thermal boundary layer by the Brownian motion of nanoparticles. Figure 14 is drawn for altered values of thermophoresis parameter to visualize thermal profile. The upsurge thermal profile for increasing thermophoresis parameter reduces heat exchange in the thermal boundary layer of pseudoplastic fluid. Similarly, Fig. 15 is sketched for concentration profile by taking the variation in thermophoresis parameter. The amount of mass exchange has been reduced causing the escalating thermal profile. Figure 16 depicts the variation in concentration profile by an increase in Schmidt number along with www.nature.com/scientificreports/ variation in k vs = 0.5 to k vs = 1 . The reaction rate boosts up because of the wider absorption interfacial syrface area on porous media. Additionally, the nanoparticles of reactants become more accelerated and collide much faster than before, causing the thickness of the concentration boundary layer. Therefore, a decline in concentration profile is witnessed. Figure 17 is the depiction of variation in drag force coefficient in the radial direction for non-Newtonian rheology by taking the thickness index of the disk m = 0.8, 1.5, 2.5 . Figure 18 represents the drag force coefficient in tangential direction for non-Newtonian rheology by taking the thickness index of disk m = 0.8, 1.5, 2.5 . Both figures show the enhancement in drag force coefficient in radial and tangential direction respectively. Heat transfer rate for non-Newtonian rheology is presented in Fig. 19 by escalating thickness index of disk m = 0.8, 1.5, 2.5 . The decline in rate of heat transfer is witnessed for the power-law index n throughout 0.3 to 1.2. Figure 20 is the sketch of heat transfer rate for non-Newtonian rheology for mounting values of Prandtl number. Heat transfer rate is seen to increase by increasing Prandtl number. Figure 21 presents the heat transfer rate for increasing value of the power-law index. Table 1 is developed for numerical values f ′ (0), g ′ (0), θ ′ (0) and −h(∞) for Newtonian fluid with thickness index of disk m = 0 and Ma, K, k 1 , k 2 , k vs = 0, and Pr = 2, r = 0.6 . The results obtained are found in excellent agreement with the previous literature. Table 2 i. Increasing surface catalyzed parameter causes the decline in concentration profile more efficiently, as it causes a boost in the pace of reaction rate. The reaction rate enhances owing to the broader absorption interfacial area on permeable media. Furthermore, the nanoparticles of reactants become more acceler- www.nature.com/scientificreports/ ated and collide much faster than earlier, instigating the thickness of the concentration boundary layer. Therefore, a decrease in concentration profile is observed. ii. The viscosity function is reliant on the viscous consistency coefficient. Improving power-law index results in the rise in variable viscosity, which in turn liable for the thickening of the boundary layer. Thus, significant augmentation in the flow occurs in the axial direction to compensate for the radial outflow. iii. The power-law index is more influential than the variable thickness disk index. iv. Velocity and temperature boundary layer gets thicker for increasing disk thickness index. v. The concentration profile for pseudoplastic fluid of a chemical species is a decreasing function of Schmidt number. vi. Larger values of the Prandtl number make the heat transfer process more efficient. vii. Increasing power-law index leads to monotonic thickening of boundary layers. Table 1 .

Numerical results of f ′ (0), g ′ (0), θ ′ (0) and −h(∞) for Newtonian fluid with thickness index of disk m = 0 and Ma, K, K 1 , K 2 , k vs = 0 and Pr = 2, r = 0.6. Table 2 .

Numerical results of f ′ (0), g ′ (0), θ ′ (0) and −h(∞) for increasing power-law index n by taking Pr = 2, r = 0.6 and disk thickness index m = 0.8. 

Enhancing thermal conductivity of fluid with nanoparticles

Numerical simulation for homogeneous-heterogeneous reactions and Newtonian heating in the silver-water nanofluid flow past a nonlinear stretched cylinder

Onset of gyrotactic microorganisms in MHD Micropolar nanofluid flow with partial slip and double stratification

Impact of hall and ion slip in a thermally stratified nanofluid flow comprising Cu and Al 2 O 3 nanoparticles with nonuniform source/sink

Comparative analysis of magnetized partially ionized copper, copper oxide-water and kerosene oil nanofluid flow with Cattaneo-Christov heat flux

Nonlinear radiation effect on MHD Carreau nanofluid flow over a radially stretching surface with zero mass flux at the surface

Heat transfer in FHD boundary layer flow with temperature dependent viscosity over a rotating disk

Investigation of entropy generation in MHD and slip flow over a rotating porous disk with variable properties

Numerical investigation of nanofluid spraying on an inclined rotating disk for cooling process

Flow and heat transfer over a rotating porous disk in a nanofluid

Similarity solution for heat convection from a porous rotating disk in a flow field

Nanofluid flow and heat transfer due to a rotating disk

Reliability of thermal desalination (solar stills) for water/wastewater treatment in light of COVID-19 (novel coronavirus "SARS-CoV-2") pandemic: What should consider

A renewable energy-driven thermoelectric-utilized solar still with external condenser loaded by silver/nanofluid for simultaneously water disinfection and desalination

Effect of Ag, Au, TiO 2 metallic/metal oxide nanoparticles in double-slope solar stills via thermodynamic and environmental analysis

Experimental investigation at a summit above 13,000 ft on active solar still water purification powered by photovoltaic: A comparative study

Experimental assessment on passive solar distillation system on Mount Tochal at the height of 3964 m: Study at high altitude

MHD) flow of Cu-water nanofluid due to a rotating disk with partial slip

Numerical investigation of natural convection and irreversibilities between two inclined concentric cylinders in presence of uniform magnetic field and radiation

An updated review on application of nanofluids in heat exchangers for saving energy

Natural convective heat transfer and entropy generation of alumina/water nanofluid in a tilted enclosure with an elliptic constant temperature: Applying magnetic field and radiation effects

Influence of a membrane on nanofluid heat transfer and irreversibilities inside a cavity with two constant-temperature semicircular sources on the lower wall: applicable to solar collectors

Effect of replacing nanofluid instead of water on heat transfer in a channel with extended surfaces under a magnetic field

First approach on nanofluid-based solar still in high altitude for water desalination and solar water disinfection (SODIS)

Non-Newtonian fluids II. Rotational currents Ostwald-de Waelescher non-Newtonian fluids

Non-Newtonian fluids v frictional resistance of discs and cones rotating in power-law non-Newtonian fluids

Flow of a power-law fluid over a rotating disk revisited

Unsteady flow of a non-Newtonian fluid above a rotating disk with heat transfer

Effects of partial slip, viscous dissipation and Joule heating on Von Kármán flow and heat transfer of an electrically conducting non-Newtonian fluid

Swirling flow of Bingham fluids above a rotating disk: an exact solution

Flow of a generalised Newtonian fluid due to a rotating disk

The neutral curve for stationary disturbances in rotating disk flow for power-law fluids

Magnetohydrodynamics thermocapillary Marangoni convection heat transfer of power-law fluids driven by temperature gradient

Mixed convection heat transfer in power law fluids over a moving conveyor along an inclined plate

Marangoni convection of power law fluids driven by power-law temperature gradient

Steady flow and heat transfer of the power-law fluid over a rotating disk

Stability of variable thickness shear deformable plates first order and high order analyses

Boundary layer flow over a stretching sheet with variable thickness

Dual solutions in a thermal diffusive flow over a stretching sheet with variable thickness

Impact of Cattaneo-Christov heat flux in the flow over a stretching sheet with variable thickness

Flow and heat transfer over a moving surface with non-linear velocity and variable thickness in a nanofluids in the presence of Brownian motion

Precipitation phenomenon of nanoparticles in power-law fluids over a rotating disk

Flow and heat transfer of Ostwald-de Waele fluid over a variable thickness rotating disk with index decreasing

Melting heat transfer and induced magnetic field effects on flow of water based nanofluid over a rotating disk with variable thickness

A simple isothermal model for homogeneous-heterogeneous reactions in boundary-layer flow. I. Equal diffusivities

Upshot of heterogeneous catalysis in a nanofluid flow over a rotating disk with slip effects and Entropy optimization analysis

Homogeneous and heterogeneous reactions in a nanofluid flow due to a rotating disk of variable thickness using HAM

Effects of heterogeneous catalysis in porous media on nanofluid-based reactions

Forchheimer flow of carbon nanotubes due to a convectively heated rotating disk with homogeneous-heterogeneous reactions

The authors extend their appreciation to the Deanship of Scientific Research at King Khalid University, Abha 61413, Saudi Arabia for funding this work through research groups program under Grant No. RGP-1-75-42.

M.R. supervised and conceived the idea; S.R. wrote the manuscript; S.K. did the software work; Y.M.C. helped in graphical depiction and funding arrangements; M.Y.M. did software work; B.Y. helped in revising the manuscript.

The research was supported by the National Natural Science Foundation of China (Grant Nos. 11971142, 11871202, 61673169, 11701176, 11626101, 11601485).

The authors declare no competing interests.

Correspondence and requests for materials should be addressed to Y.-M.C.Reprints and permissions information is available at www.nature.com/reprints.Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.