TY - JOUR
AU1 - Zhang, Xiao-Hong
AU2 - Algehyne, Ebrahem A.
AU3 - Alshehri, Maryam G.
AU4 - Bilal, Muhammad
AU5 - Khan, Muhammad Altaf
AU6 - Muhammad, Taseer
AB - 1. Introduction Fluid flow over the substrate of a rotating disc has garnered considerable attention because of its consequences in real phenomena. Flow over the rotating disc surface is commonly used in electric power generation systems, co rotating machines, aerodynamics engineering, rotating equipment, geothermal sector, chemical reaction, and computer processing [1]. Accurate analysis and precise measurements of reaction are the prerequisite for better understanding the factors that overcome kinetics of reaction. Electrode system of rotating disk is broadly used to examine the reactions kinetics and that are suffering from mass transmission deficiency [2]. Tassaddiq et al. [3] examined the MHD, an incompressible iron ferrite and carbon nanotubes CNTs hybrid nanoliquid over an impermeable spinning disk. They perceived that the fluid temperature and velocity are considerably increased by increasing disc rotating velocity. Mustafa et al. [4] used the Bingham fluid model to analyze the heat transmission and viscous dissipation for swirling viscoelastic fluid flow crossing over a permeable spinning disc. Using a greater suction velocity at the disc allows for a faster thermal efficiency. However, as the fluid yield stress increases, the heat transfer is expected to deteriorate. Hafeez and Khan [5] scrutinized the mass and heat flux theories in Oldroyd-B fluid flow over an extending disk. It should be noted that the solutal relaxation time coefficient has a negative impact on mass transmission rate. Gul et al. [6] modeled the hybrid nanofluid flow between disk and cone with several cases. Based on the analysis of the results, it was concluded that a rotating disc that a rotating disc with a fixed cone would attain the intended cooling of disk-cone tools if the surface temperature remains steady. Fluid that is commonly used Ethanol, kerosene oil and water etc., which play a prominent role in energy transfer, heating and air conditioning processes, energy generation, and other small electronics mechanisms. However, since these liquids have a poor thermal energy transfer capacity, they are unable to accommodate the need for high heat exchange rates. Nanometer-sized particles (1-100nm), also recognized as nanomaterials, are applied to a natural fluid to substitute for this inadequacy and increase thermal conductivity. To quantify the thermophysical properties of various kinds of nanoparticles, the researcher used a variety of techniques [7]. Nanoparticles are used to clean surfaces in engineering because they have the ability to disperse and wet oil. It improves thermal efficiency, which is important in the fields of energy production, hypoxia, micro fabrication and metallurgy. Magnesium oxide MgO is a compound made up of Mg2+ and O2- ions that are bonded together through a powerful ionic bond. It is made by calcination of Mg (OH) 2 and MgCO3 at temperature range from 700°C-1500°C. Most of it is utilized in the hematologic and electrical industries. Similarly, Ag particles with antibacterial properties are used to control microorganism growth in wound, burns and surgical apparatus and infrastructure. Microorganisms are particularly toxic to compound based on silver and its ions. Hussanan et al. [8] investigated oxide nanomaterials in engine nanoliquids and water for energy up gradation. Using the finite element technique, Ghalambaz et al. [9] investigated convection flow of a hybrid nanofluid (Ag–MgO/water) within a cavity. Copper-Alumina micro materials were listed for the base fluid. Aluminum Oxide was used by Motlagh and Soltanpour [10] to elaborate the heat transition in a non-flat cavity. Acharya et al. [11] investigated the hybrid fluid flow over a spinning disc with the consequences of MHD and Hall current. They introduced the TiO2 and Cu nanomaterials to creative a class of nanoliquid. The magnetic interaction on flow with CNTs nano powders via a porous tube is studied by Akber et al. [12]. Shah et al. [13] reported that an angled magnetic field, performance of copper selenide and the efficacy of chemical catalytic reactors effects the energy and concentration distribution of nanofluids. Sowmya et al. [14] uncovered natural convective thermal transference in a rectangular chamber with two heated fins on the bottom wall for a hybrid nanofluid flow employing iron oxide and silver nanomaterials. Zhou et al. [15] Khan addressed Von Karman’s classic swirling flow for Maxwell fluid over a porous whirling disc with a constant suction/injection mechanism. In the engineering sectors, complex nonlinear boundary value problems that cannot be addressed usually are often observed. Convergence is sensitive to the relaxation parameters and initial approach for many problems that are commonly solved by other numerical algorithms. The PCM’s goal is to explore the method’s generalizability as a feasible alternative to nonlinear problems [16]. Shuaib et al. [17] emphasized the 3D unsteady fluid and heat propagation across the substrate of a non-flat stretchable rotating disc. The fluid has investigated under the effects of a magnetic strength from the outside. The feature of an ionic transition boundary layer flow over a gyrating disc was discovered by Shuaib et al. [18]. The Poisson’s, Nernst-Planck equation and Navier Stokes equations, were used to solve for the ionic compounds. Wang et al. [19] presented stability analysis of nonlinear problems for engineering applications through parametric continuation algorithm. They also looked at the static bifurcation that happens when solving nonlinear initial value problems with different characteristic roots and provided an algorithm for instantly calculating the bifurcation points. The aim of this study is to modify an idea of Ref. [20,21], by investigating the influence of two separate nanoparticles, Silver Ag/water and magnesium oxide MgO/water hybrid nanofluids, on a wavy spinning disc that moves upward and downward. This research is being considered in order to increase the thermal conductivity of the fluid flow. The modelled equations are solved numerically using the parametric continuation process (PCM), and the results have been validated and compared with Matlab’s package boundary value solver (bvp4c). All outcomes are in the closest possible agreement with one another. 2. Mathematical formulation The physical analysis of the problem, physicochemical characteristics and the equation of motion has been discussed in this section. 2.1 Physical description of the problem We have assumed a 3D flow of Ag-MgO/Water (Silver and magnesium-oxide) hybrid nanoliquid over a wavy up and down moving spinning disk. Initially the disk is at a(0) = h. At a vertical distance Z = a(t) the disk is moving with velocity ω = a(t). The consequences of buoyancy impacts are negligible. The disk is revolving about z-axis with angular velocity Ω(t), and it has been considered that the Ag-MgO nano-size particles are scattered consistent. The magnetic field of constant magnitude is uniformly applied with and respectively, where and are unit vectors. 2.2 Equation of motion The modeled equations, based on above presumption are expressed as [20,21]: (1) (2) (3) (4) (5) (6) (7) Here Fr and Fθ along x and z direction are the body forces, respectively, which can be stated as [21]: (8) (9) Here, Ha is , in which θ is the direction and B0 is the magnitude of magnetic field, where the velocity components are specified as u,v,w. 2.3 Boundary condition The initial and boundary conditions for wavy gyrating disk are: (10) 2.4 Thermophysical properties of nanoliquid The density and specific heat capacity of the hybrid nanofluid can be revealed as [21]: (11) (12) where (Cp)s1,(Cp)s2 are specific heat capacity, ρs1,ρs2 are the density and φ1,φ2 are the volume friction of the Ag-MgO hybrid nanoliquid respectively shown in Table 1. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 1. The numerical properties of water and hybrid nanofluid [20]. https://doi.org/10.1371/journal.pone.0254457.t001 The hybrid nanofluid viscosity μhnf is expressed as [22]: (13) Here, the thermal conductivity and Prandtl number is expressed as [21]: (14) 2.5 Karman’s approach We use the following similarity framework, to reduce Eqs (1)–(7) and (10) to the system of ODEs [20]: (15) As a result of Eq (15), we get the following: (16) (17) (18) (19) (20) (21) The transform conditions are: (22) Sign ω edify disk’s rotation, S controlled the up and down fluctuating of the disk and γ is the is the thermal energy parameter, which can be defined as [18]: (23) The dimensionless form of skin friction and Nusselt number is: (24) Where, τwr and τwϕ represent the radial and transverse stress respectively. 3. Numerical solution The following steps have been used, while solving the system of ODE (16–21) and their boundary conditions (22): Step 1: Reducing the system of modeled equations to the first order ODE we introduced the following similarity variables: (25) Using Eq (25) into the BVP (16–21) and (22), we get: (26) (27) (28) (29) (30) (31) the boundary conditions are: (32) Step 2: By introducing the embedding parameter p We introduced the parameter p in the above system of Eqs (26–31) as follow: (33) (34) (35) (36) (37) (38) Step 3: Differentiating w.r.t parameter ‘p’ Eqs (33–38) after differentiating w.r.t parameter p reform to the following structure with respect to parameter p (39) where, R is the remainder and A is the coefficient matrix: (40) where i = 1, 2, ………11. Step 4: Apply specify Cauchy problem and superposition principle for each component (41) where U, W are unknown vector functions and a is unknown blend coefficient. Solving the Cauchy problems for each component (42) (43) by putting Eq (41) into the original Eq (39), we get (44) Step 5: Solving the Cauchy problems The numerical implicit scheme is implemented in the proposed problem, presented as follow: From Eqs (42) and (43) (45) (46) here, we obtained the iterative form of the problem. (47)(48) 2.1 Physical description of the problem We have assumed a 3D flow of Ag-MgO/Water (Silver and magnesium-oxide) hybrid nanoliquid over a wavy up and down moving spinning disk. Initially the disk is at a(0) = h. At a vertical distance Z = a(t) the disk is moving with velocity ω = a(t). The consequences of buoyancy impacts are negligible. The disk is revolving about z-axis with angular velocity Ω(t), and it has been considered that the Ag-MgO nano-size particles are scattered consistent. The magnetic field of constant magnitude is uniformly applied with and respectively, where and are unit vectors. 2.2 Equation of motion The modeled equations, based on above presumption are expressed as [20,21]: (1) (2) (3) (4) (5) (6) (7) Here Fr and Fθ along x and z direction are the body forces, respectively, which can be stated as [21]: (8) (9) Here, Ha is , in which θ is the direction and B0 is the magnitude of magnetic field, where the velocity components are specified as u,v,w. 2.3 Boundary condition The initial and boundary conditions for wavy gyrating disk are: (10) 2.4 Thermophysical properties of nanoliquid The density and specific heat capacity of the hybrid nanofluid can be revealed as [21]: (11) (12) where (Cp)s1,(Cp)s2 are specific heat capacity, ρs1,ρs2 are the density and φ1,φ2 are the volume friction of the Ag-MgO hybrid nanoliquid respectively shown in Table 1. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 1. The numerical properties of water and hybrid nanofluid [20]. https://doi.org/10.1371/journal.pone.0254457.t001 The hybrid nanofluid viscosity μhnf is expressed as [22]: (13) Here, the thermal conductivity and Prandtl number is expressed as [21]: (14) 2.5 Karman’s approach We use the following similarity framework, to reduce Eqs (1)–(7) and (10) to the system of ODEs [20]: (15) As a result of Eq (15), we get the following: (16) (17) (18) (19) (20) (21) The transform conditions are: (22) Sign ω edify disk’s rotation, S controlled the up and down fluctuating of the disk and γ is the is the thermal energy parameter, which can be defined as [18]: (23) The dimensionless form of skin friction and Nusselt number is: (24) Where, τwr and τwϕ represent the radial and transverse stress respectively. 3. Numerical solution The following steps have been used, while solving the system of ODE (16–21) and their boundary conditions (22): Step 1: Reducing the system of modeled equations to the first order ODE we introduced the following similarity variables: (25) Using Eq (25) into the BVP (16–21) and (22), we get: (26) (27) (28) (29) (30) (31) the boundary conditions are: (32) Step 2: By introducing the embedding parameter p We introduced the parameter p in the above system of Eqs (26–31) as follow: (33) (34) (35) (36) (37) (38) Step 3: Differentiating w.r.t parameter ‘p’ Eqs (33–38) after differentiating w.r.t parameter p reform to the following structure with respect to parameter p (39) where, R is the remainder and A is the coefficient matrix: (40) where i = 1, 2, ………11. Step 4: Apply specify Cauchy problem and superposition principle for each component (41) where U, W are unknown vector functions and a is unknown blend coefficient. Solving the Cauchy problems for each component (42) (43) by putting Eq (41) into the original Eq (39), we get (44) Step 5: Solving the Cauchy problems The numerical implicit scheme is implemented in the proposed problem, presented as follow: From Eqs (42) and (43) (45) (46) here, we obtained the iterative form of the problem. (47)(48) 4. Analysis and discussion of results 4.1 Analysis The numerical outcomes of the system of DE (Differential Equations) are calculated via using PCM, and the bvp4c technique has been used to compare and validate the results. The results have been revealed through Figs (2–11). Fig 1 displays the mechanism of fluid flow over a curvy surface with up and down fluctuation under the consequences of magnetic field. In current paper, the intention of the analysis is to boost the thermal energy transmitting rates for a range of industrial implementations. When compared to a flat surface, energy transition is increased up to 15% due to the wavy swirling surface. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 1. https://doi.org/10.1371/journal.pone.0254457.g001 4.2 Discussion of results Figs (2–5) elaborate the upshot of volume friction parameters (ϕ1 or ϕAg and ϕ2 or ϕMgO) on axial velocity f(η) and energy profile θ(η) respectively. It can be perceived that the velocity and temperature transmission rate enhance with the rising values of ϕ1&ϕ2. Physically, the increasing number of nanoparticles (Ag-MgO) in base fluid reduces the specific heat capacity of water and improve the thermal diffusivity, because the specific heat capacity of silver and magnesium-oxide is less than water. Eventually, fluid losses its viscosity, which cause rises in temperature and velocity of fluid. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 2. https://doi.org/10.1371/journal.pone.0254457.g002 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 3. https://doi.org/10.1371/journal.pone.0254457.g003 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 4. https://doi.org/10.1371/journal.pone.0254457.g004 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 5. https://doi.org/10.1371/journal.pone.0254457.g005 The unsteadiness parameter S upshot on axial f(η) and azimuthal velocity h(η) profiles is demonstrated through Figs 6 & 7, respectively. It has been concluded that the axial velocity enhances with the positive increment of unsteadiness parameter S, while azimuthal velocity reduces. Because the kinetic viscosity of fluid deceases with its effect, that’s why such scenario has been observed in Fig 6. The effects of disk rotation ω parameter on radial velocity g(η) is shown via Fig 8. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 6. https://doi.org/10.1371/journal.pone.0254457.g006 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 7. https://doi.org/10.1371/journal.pone.0254457.g007 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 8. https://doi.org/10.1371/journal.pone.0254457.g008 It’s reasonable to assume that improvement rate of disk rotation significantly encourages the kinematic energy of the fluid ions, that leads to an increase in velocity g(η), which also generates heat. The disc surface eventually heats up and improves the fluid temperature θ(η). Fig 9 revealed that the effect of parameter γ reduces the temperature profile. Because the fluid particle releases thermal energy during upward/downward movement, which causes the declination of fluid temperature. The rotating disk’s upward and downward velocity is regulated by the parameter γ. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 9. https://doi.org/10.1371/journal.pone.0254457.g009 Fig 10 explored the consequences of Prandtl number Pr on energy profile. It can be noticed that the temperature declines by the action of Prandtl number. Physically, the action of Prandtl number directly affect the dynamic viscosity and specific heat capacity of fluid, as result the fluid thermal capability declines. Figs 11 & 12 demonstrate the influence of Batchlor number Bt on magnetic strength profile along axial m’(η) and radial direction n(η), respectively. The magnetic strength profile show declination versus the influence of Batchlor number, because the increment in Bt negatively affect the molecular diffusion rate of fluid. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 10. https://doi.org/10.1371/journal.pone.0254457.g010 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 11. https://doi.org/10.1371/journal.pone.0254457.g011 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 12. https://doi.org/10.1371/journal.pone.0254457.g012 Table 1 elaborate the thermophysical characteristics of silver, base fluid, and magnesium-oxide. Table 2 revealed the comparison between PCM and bvp4c method numerically for all velocities. It can be presumed that the parametric continuation method has strong convergence versus other numerical method. Tables 3 and 4 display the influence of nanoparticles quantity ϕ on velocity radial and azimuthal direction, skin fraction and Nusselt number. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 2. Comparative analysis between PCM and bvp4c method. https://doi.org/10.1371/journal.pone.0254457.t002 Download: PPT PowerPoint slide PNG larger image TIFF original image Table 3. The comparative analysis of silver and magnesium-oxide nanoparticles on velocity profile. https://doi.org/10.1371/journal.pone.0254457.t003 Download: PPT PowerPoint slide PNG larger image TIFF original image Table 4. Quantitively analysis for Nusselt number and Skin fraction. https://doi.org/10.1371/journal.pone.0254457.t004 4.1 Analysis The numerical outcomes of the system of DE (Differential Equations) are calculated via using PCM, and the bvp4c technique has been used to compare and validate the results. The results have been revealed through Figs (2–11). Fig 1 displays the mechanism of fluid flow over a curvy surface with up and down fluctuation under the consequences of magnetic field. In current paper, the intention of the analysis is to boost the thermal energy transmitting rates for a range of industrial implementations. When compared to a flat surface, energy transition is increased up to 15% due to the wavy swirling surface. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 1. https://doi.org/10.1371/journal.pone.0254457.g001 4.2 Discussion of results Figs (2–5) elaborate the upshot of volume friction parameters (ϕ1 or ϕAg and ϕ2 or ϕMgO) on axial velocity f(η) and energy profile θ(η) respectively. It can be perceived that the velocity and temperature transmission rate enhance with the rising values of ϕ1&ϕ2. Physically, the increasing number of nanoparticles (Ag-MgO) in base fluid reduces the specific heat capacity of water and improve the thermal diffusivity, because the specific heat capacity of silver and magnesium-oxide is less than water. Eventually, fluid losses its viscosity, which cause rises in temperature and velocity of fluid. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 2. https://doi.org/10.1371/journal.pone.0254457.g002 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 3. https://doi.org/10.1371/journal.pone.0254457.g003 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 4. https://doi.org/10.1371/journal.pone.0254457.g004 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 5. https://doi.org/10.1371/journal.pone.0254457.g005 The unsteadiness parameter S upshot on axial f(η) and azimuthal velocity h(η) profiles is demonstrated through Figs 6 & 7, respectively. It has been concluded that the axial velocity enhances with the positive increment of unsteadiness parameter S, while azimuthal velocity reduces. Because the kinetic viscosity of fluid deceases with its effect, that’s why such scenario has been observed in Fig 6. The effects of disk rotation ω parameter on radial velocity g(η) is shown via Fig 8. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 6. https://doi.org/10.1371/journal.pone.0254457.g006 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 7. https://doi.org/10.1371/journal.pone.0254457.g007 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 8. https://doi.org/10.1371/journal.pone.0254457.g008 It’s reasonable to assume that improvement rate of disk rotation significantly encourages the kinematic energy of the fluid ions, that leads to an increase in velocity g(η), which also generates heat. The disc surface eventually heats up and improves the fluid temperature θ(η). Fig 9 revealed that the effect of parameter γ reduces the temperature profile. Because the fluid particle releases thermal energy during upward/downward movement, which causes the declination of fluid temperature. The rotating disk’s upward and downward velocity is regulated by the parameter γ. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 9. https://doi.org/10.1371/journal.pone.0254457.g009 Fig 10 explored the consequences of Prandtl number Pr on energy profile. It can be noticed that the temperature declines by the action of Prandtl number. Physically, the action of Prandtl number directly affect the dynamic viscosity and specific heat capacity of fluid, as result the fluid thermal capability declines. Figs 11 & 12 demonstrate the influence of Batchlor number Bt on magnetic strength profile along axial m’(η) and radial direction n(η), respectively. The magnetic strength profile show declination versus the influence of Batchlor number, because the increment in Bt negatively affect the molecular diffusion rate of fluid. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 10. https://doi.org/10.1371/journal.pone.0254457.g010 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 11. https://doi.org/10.1371/journal.pone.0254457.g011 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 12. https://doi.org/10.1371/journal.pone.0254457.g012 Table 1 elaborate the thermophysical characteristics of silver, base fluid, and magnesium-oxide. Table 2 revealed the comparison between PCM and bvp4c method numerically for all velocities. It can be presumed that the parametric continuation method has strong convergence versus other numerical method. Tables 3 and 4 display the influence of nanoparticles quantity ϕ on velocity radial and azimuthal direction, skin fraction and Nusselt number. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 2. Comparative analysis between PCM and bvp4c method. https://doi.org/10.1371/journal.pone.0254457.t002 Download: PPT PowerPoint slide PNG larger image TIFF original image Table 3. The comparative analysis of silver and magnesium-oxide nanoparticles on velocity profile. https://doi.org/10.1371/journal.pone.0254457.t003 Download: PPT PowerPoint slide PNG larger image TIFF original image Table 4. Quantitively analysis for Nusselt number and Skin fraction. https://doi.org/10.1371/journal.pone.0254457.t004 5. Conclusion In the present study, we numerically scrutinized an unsteady 3D numerical model for Ag-MgO hybrid nanoliquid flow with mass and energy transition caused by the up and downward fluctuation of a wavy gyrating disc. Magnesium oxide and silver nanoparticles have been heavily reported to have broad-spectrum antibacterial operations among metal oxides and metals. However, in current paper, the study’s goal is to increase the rate of thermal energy propagation for a variety of industrial applications. For the purpose, the problem has been formulated as a system of partial differential equations (Navier Stokes and Maxwell). During the study, the following observation have been made: Silver nanoparticles are without a doubt the most used inorganic nanoparticles, with numerous innovations in biomaterial’s detection and antimicrobial operations. In comparison to a smooth substrate the wavy spinning surface raises the heat transfer rate by 15%. Magnesium oxide (MgO) is made up of Mg2+ and O2- ions that are held together by a strong ionic bond. Which can be synthesized at temperatures ranging from 700 0C to 1500 0C degrees Celsius and is primarily used in refractory and electrical applications. The up and downward fluctuation of the curly spinning disc affects fluid temperature and velocity in a constructive way. Water’s thermophysical properties are improved by the close bonds between water atoms (H++OH-) and silver ions Ag+. Hybrid nanoliquids are more effective at overcoming low energy transfer. For example, it improves carrier fluid thermal efficiency, which is essential in power generation, microfabrication, hyperthermia, metallurgical fields and air conditioning. The parametric continuation method is a strong numerical technique for highly nonlinear system of PDE than others numerical algorithm. Acknowledgments The second and third authors would like to acknowledge the support provided by the University of Tabuk. Furthermore, the second author would like to acknowledge the support provided by the Nanotechnology Research Unit (NRU), University of Tabuk, Saudi Arabia.
TI - The parametric study of hybrid nanofluid flow with heat transition characteristics over a fluctuating spinning disk
JF - PLoS ONE
DO - 10.1371/journal.pone.0254457
DA - 2021-08-16
UR - https://www.deepdyve.com/lp/public-library-of-science-plos-journal/the-parametric-study-of-hybrid-nanofluid-flow-with-heat-transition-Ont0nBfvoo
SP - e0254457
VL - 16
IS - 8
DP - DeepDyve
ER -