Influence of buoyancy forces in MHD non-Newtonian convective nanofluid utilizing Buongiornos Model i
* Author to whom correspondence should be addressed.
Journal of Thermal Engineering 2024, Vol. 10, Issue 5, pp. 1107-1119; doi.org/10.14744/thermal.0000854
Abstract
Keywords: Buongiorno's Model; Buoyancy forces; Casson Fluid; Chemical Reaction; Nanofluid; Shooting Method
Introduction
Stream conduct over extending surfaces has drawn attention of many researchers because of its wide space of modern and fabricating applications like fake strands, petrol
ventures, metal turning, polymer handling and so forth. Crane [1] concentrated on the stream towards an extending sheet. Sarma and Rao [2] explored the stream conduct in viscoelastic liquid over an extending sheet. Additionally, Radiation impact with hotness source was considered by
*Corresponding author. *E-mail address: saloni.20.jindal@gmail.com This paper was recommended for publication in revised form by Editor-in-Chief Ahmet Selim Dalkılıç Published by Yıldız Technical University Press, İstanbul, Turkey Copyright 2021, Yıldız Technical University. This is an open access article under the CC BY-NC license (http://creativecommons.org/licenses/by-nc/4.0/).
Manjunatha et al. [3] over permeable medium. In ongoing many years, heat move angles are substantially more significant for issues identified with extending sheet because of its warming and cooling essential variables for making superior grade of eventual outcome. Fourier [4] proposed law of heat conduction very first and revealed that this law gives premise to discover the hotness move conduct under various conditions. Yet, this model has significant disadvantage that it offers numerical expression of energy in illustrative structure. From that point onward, Cattaneo [5] has acquired new statures of Fourier’s law with warm unwinding time because of the conditions that are in illustrative structure which moved into exaggerated structure. Very recently, heat dissipation impact under different physical conditions has been analyzed by researchers [6-9]. Researchers are paying more attention to nanofluid because of its physical properties and uses in many areas, mainly production, industry, and medicine. These tiny fluids have lots of important properties for research, like low resistivity and unique thermophysical aspects. In the last few decades, the word “nanofluid” has been a topic of interest in a number of scientific groups because it is able to transfer heat quickly. Choi and Eastman [10] identified a brand-new class of fluids called nanofluids. Because they have a high rate of heat transmission, nanofluids are of interest to many scientists and academics. Compared to solids, fluids are less effective in transferring heat. The creation of ultrafine particles enhances the thermal conductivity of nanofluids, according to Xuan and Li [11]. Through Brownian motion and thermophoresis diffusion, Buongiorno [12] created a methodical strategy that boosted the heat exchange phenomenon. Moreover, influence of activation energy over porous surface in presence of Carreau nanofluid has been presented by Shahid et al. [13]. Din et al. [14] studied the significance of bio-convection and slip conditions over wedge by utilizing tangential hyperbolic and Carreau nanofluids. Three-dimensional second-class nanofluid flow across a stretched surface has been studied by Ahmad et al. [15] using the Optimal Homotopy Analysis approach. The investigation of non-Newtonian fluids has numerous real-world uses in technology and other fields of study. Non-Newtonian fluid mechanics has a lot of interesting and important uses in the fields of engineering, healthcare, and the applied sciences. Since both the surrounding air and the ocean are fluids, the studies of fluids are essential to metrology, oceanography, and hydrology. Identifying the flow of various biological fluids is important for successful medical treatment. When hybrid nanoliquids were subjected to nonlinear solar radiation, Acharya et al. [16] investigated how various solar thermal devices affected the flow fluctuations and heat exchange properties of the materials. Very recently, Rana et al. [17] uses ANN to predict the stable solutions and critical points in viscoelastic fluids induced by horizontal sheet. In recent years, nanofluid flow in presence of Magnetohydrodynamics under different physical aspects has been explored by authors [18-21].
The field of MHD has important applications in the medical and industrial sciences. MHD has a wide range of potential uses, including endoscopy, cell separation, cancer tumor treatment, drug targeting and surgical blood flow control. MHD utilized in medication i.e., growth or disease treatment. A stagnation point MHD nanofluid stream was mathematically examined by Anwar et al. [22] induced by stretching surface. Then Shawky et al. [23] discussed Williamson nanofluid towards an extending plot. Vajravelu and Cannon [24] inspected the stream conduct of liquid on non-direct extending sheet. Additionally, a blended MHD nanofluid stream alongside entropy and convective investigation on a non-straight extending sheet was addressed by Matin et al. [25]. Jain and Choudhary [26] examined the Soret and Dufour impacts with MHD stream utilizing compound responses. Siddheshwar and Mahabaleshwar [27] concentrated on the issue of hotness and stream transportation on non-straight extending sheet with pull/infusion. By using hybrid nanofluids flow in presence of Stefan blowing bio-convection has been analyzed by Rana et al. [28]. The opinions of a few different specialists [29-36], considered a different aspect of those kinds of problems under diverse circumstances. Motivated by the already existing literature uncovers that no such study has been made to the till now. Current research focuses on the effect of transport phenomenon by creating outcome of MHD nanofluid flow and free stream in presence of compound response towards an extending surface. The philosophy embraced addresses whole framework through Runge Kutta Fehlberg strategy by following shooting procedure utilizing ODE45 solver. The numerical troubles showing up in the nanofluid conditions drove us to utilize the mathematical methodology. Utilizing the concept of buoyancy forces and nanoparticles, current research is relevant to novel microbial fuel cell technologies in exponential surface. The current research’s results assist the biological, technological, and manufacturing domains in producing desired products that make this study to novel one.
Materials And Methods
3D non-Newtonian nanofluid model with magnetohydrodynamic flow induced by stretching surface is considered in the current research (Figure 1). In this model, thermophoresis and Brownian motion impact of nanofluid has been incorporated among stretched velocities ,
direction serially. Buoyancy forces along with free stream velocity are incorporated in the current study. Modeled equations are shown as ( [14], [37-39]): (1)
(10) (11) Reduced B.C’s are: (12) Cfx, Cfy, Nux and Shx are defined as:
Now, a defined new set of variables is provided below to convert the equations previously defined into the first-order O.D.E’s:
Where Rex shows Reynolds’s number. Numerical Solution The system of D.E’s (8)-(11) and B.C’s (12) has been solved with Runge-Kutta-Fehlberg approach by following shooting procedure, as indicated in a flowchart in Figure
2. To get a numerical solution, the step size is set to 0.01
and the maximum value to 10. Main advantage of this technique is that it has 5th order truncation error over other numerical techniques. Additionally, as compared to other numerical techniques, the computation of the solution is simpler and easier. The well-known programme MATLAB is used in the current analysis to do computations utilizing the ODE45 solver. 1st order D.E’s have been created from the coupled D.E’s (8)-(11) in this case. The controlling nonlinear O.D.E’s are reconfigured as follows for this purpose: (15)
Equations (15)-(18) are changed into the following system of D.E’s by using Equation (19) to them: (20) (21)
parameter, L2 = 0.3; solutal buoyancy parameter, M = 0.1; magnetic parameter, c = 0.1; stretching ratio parameter, β = 0.1; Casson Fluid parameter, Pr = 0.72; Prandtl number, Nb = 2.5; Brownian motion, Nt = 0.5; thermophoresis, Le = 2; Lewis number, Cr = 0.5; chemical reaction parameter are shown in Table 2-4. Influence of free stream velocity λ over velocity profile f '(ξ) has been presented in Figure 3(a) and it shows that velocity profile enhances with enhancement in free stream velocity ranging from 0.00 ≤ λ ≤ 0.20. Figure 3(b) shows variety in liquid speed against attractive boundary M. Presence of attractive boundary M opposes the liquid molecule to move uninhibitedly and fundamental purpose for the obstruction is that attractive boundary M produces Lorentz power and this attraction conduct can be taken on for controlling the smooth motion. Accordingly, improvement in the worth of attractive boundary M causes the declination of speed dissemination. Figure 4(a) inspects temperature appropriation variety against the liquid boundary Brownian movement boundary Nb. Liquid temperature upgrades for higher Nb and in result neighborhood Nusselt number declines. Figure 4(b) ponders the effect of liquid temperature under the result of Nt. Temperature slope tumbles down for higher upsides of
Results And Discussion
For adjusting fluid parameter Pr in the absence of magnetohydrodynamic flow, the estimated findings for Nusselt number are compared with those from Abolbashari et al. [40] and Khan and Pop [41] in the current analysis (Table 1). When compared to the conclusions of the current literature, it is found that the results are very accurate. The results (f "(0), g"(0)), -θ’(0), -ϕ’(0), respectively, for fixed values of governing fluid parameters as e = 0.1; thermal Biot number, h = 0.1; concentration Biot number, λ = 0.1; free stream velocity parameter, L1 = 0.3; thermal buoyancy
Table 1. Comparison of Nusselt number along with residual error against Prandtl number Pr
Table 2. Values of f '' (0) and g'' (0) for variables β, λ and M β
thermophoresis boundary Nt that decrease of conduction of nanoparticles. Thus, temperature ascends for higher thermophoresis boundary Nt as illustrated in Figure 4(b). Figure 4(c) describes the variation in temperature profile
as a function of e; (0.15 ≤ e ≤ 0.75) and it can be seen that temperature and convective heating increase as e increases. Additionally, plots diverge for around 0.0 to 3.0 and convergent with ξ → ∞. Additionally, when e changes between
0.3. and 0.7, the Nusselt number increases in the range of
0.125323233938653 to 0.366860190093477. Further, Pr is the “dimensionless” amount which is the proportion of energy diffusivity to the warm diffusivity. Assuming that the liquid is taken to be moderately more thick then the worth of Pr increments and there is less convection in heat move rate. As the worth of Pr expands it will decrease the warm diffusivity. This low warm diffusivity declines temperature field as shown via Figure 4(d). Influence of concentration distribution ϕ(ξ) against fluid parameters Nt, h, Cr, Pr has been displayed via Figure 5(a)-5(d) in a serial manner. Figure 4(a) shows the influence of concentration distribution ϕ(ξ) against thermophoresis Nt within the range 1 ≤ Nt ≤9. This diagram demonstrates how the nanoparticle volume fraction increases with increasing thermophoresis Nt. Fundamentally, when a molecule applies thermophoresis to another molecule, the result is the development of particle movement from a hotter to a cooler region, which shows strengthening of the
nanoparticle volume fraction as observed in Figure 5(a). Additionally, nanoparticle concentration is directly influenced by concentration Biot number. Figure 5(b) shows how the concentration profile against h has changed over time. When h changes from 0.1 to 0.5, enhancement in h increases concentration profile and decreases Sherwood number from 0.093328137143792 to 0.373362399233673. Figure 5(c) illustrates the variation in nanoparticle concentration in relation to Cr (a chemical reaction parameter). This graph demonstrates how concentration distributions gradually decrease as Cr increases (0.0 ≤ Cr ≤ 2.0) and become stable distant from the surface to satisfy boundary conditions. Additionally, when Cr increases from 0.0 to 2.0, a stronger chemical reaction cause a Sherwood number augmented in the range of 0.091974413822926 to 0.095246246753709. Variation in concentration against Pr (1.0 ≤ Pr ≤ 9.0). As shown in Figure 5(d), concentration distribution decreases as Pr increases. Figure 6(a) manipulates effect of skin
friction coefficient against Pr in the range 1.0 to 5.0 for different β (0.08, 0.09, 0.10, 0.11, 0.12) and it is found that skin friction declines with inclination in Pr while augmented for higher β. Figure 6(b) shows graph of Nusselt number against Pr (1 ≤ Pr ≤5) for different thermal Biot number e (0.1 ≤ e ≤ 0.5) and this plot shows that Nusselt number enhances for greater Pr and e. Influence of Sherwood number against fluid parameter Pr for different concentration Biot number h (0.1 ≤ h ≤ 0.5) has been displayed via Figure 6(c). Sherwood number augmented with augmentation in both Pr and h as illustrated in Figure 6(c). Figure 6(d) illustrates impact of Sherwood number against Brownian motion Nb (1.0 ≤ Nb ≤ 3.0) for different Concentration Biot number h (0.1 ≤ h ≤ 0.5). Sherwood number increases very slowly and thus very minor effect of Sherwood number is noticed out for Pr and simultaneously it is found that higher h helps to increase in Sherwood number as illustrated in Figure 6(d).
Figure 7(a)-7(d) shows influence of skin friction coefficient against Le & Cr, Pr & Le , Nb & Nt and Pr & Nt respectively. Figure 7(a) manifests skin friction coefficient under the influence of Lewis number Le (1 ≤ Le ≤ 3) & chemical reaction Cr (1 ≤ Le ≤ 2) and this contour plot shows that skin friction coefficient declines with rise in values of Le (1 ≤ Le ≤ 3) and Pr (1 ≤ Pr ≤ 4). Influence of skin friction against Pr & Le has been illustrated in Figure 7(b) via contour plot. In this plot, Pr plays a significant role to falls down skin friction coefficient in comparison to Lewis number. Figure 7(c) manifests impact of skin friction across Nb (1 ≤ Nb ≤ 3) & Nt (1 ≤ Nt ≤ 2) and this contour shows that skin friction coefficient falls down with rise in both nanofluid parameters Nb & Nt. Moreover, combined impact of Pr (1≤ Pr ≤ 4) & Nt (1 ≤ Nt ≤ 3) over skin friction coefficient has been displayed via Figure 7(d). It is noticed that skin friction coefficient declines due to more dominancy of nanofluid parameter Nt in comparison to Prandtl number Pr.
Figure 6. (a) Skin friction coefficient against β (0.08 ≤ β ≤ 0.12) & Pr (1.0 ≤ Pr ≤ 5.0), (b) Nusselt number against e (0.1 ≤ e ≤ 0.5) & Pr (1.0 ≤ Pr ≤ 5.0), (c) Sherwood number against h (0.1 ≤ e ≤ 0.5) & Pr (1.0 ≤ Pr ≤ 5.0) and (d) Sherwood number against h (0.1 ≤ h ≤ 0.5) & Nb (1.0 ≤ Nb ≤ 3.0).
Conclusion
Current investigation deals with the study of 3D magnetohydrodynamic flow with buoyancy forces induced by stretching surface. Major outcomes of the study are: • Augmentation in concentration is for larger thermophoresis parameter Nt. • Velocity lessens down for larger β and magnetic parameter M due to produced Lorentz drag force that ultimately enhances the quality of final product. • Temperature and concentration rises with rise in thermal and concentration Biot numbers serially.
Larger Prandtl number Pr declines temperature and concentration distribution. • Skin friction coefficient declines with rise in values of Le and Pr and falls down with rise in both the values of nanofluid parameters Nb &Nt. • Comparative study has been performed under [40, 41] to get higher accuracy and validation of the current research. Results of the current survey have important implications for the regulation of transport phenomena and fluid velocity in a wide variety of formulation processes, which in turn benefits production, implements of technology, and •
Figure 7. Skin friction coefficient for (a) Lewis number Le & Chemical reaction Cr, (b) Prandtl number Pr & Lewis number Le, (c) Brownian motion Nb & Thermophoresis parameter Nt and (d) Prandtl number Pr & Thermophoresis parameter Nt.
biological networks in their pursuit of the ideal product standard.
Nomenclature
Cartesian coordinates Casson fluid parameter Concentration Biot number Free stream velocity parameter Thermophoresis parameter Chemical reaction parameter Magnetic parameter Thermal Biot number
Thermal buoyancy parameter Temperature profile Gravitational acceleration Skin friction coefficient along x-axis Solutal buoyancy parameter Prandtl number Horizontal velocity profile Concentration profile Skin friction coefficient along y-axis Sherwood number Wall shear along x- direction Nusselt number Convective mass transfer coefficient
T τyz C∞ Tw T∞ Cw ρ Uw, Vw u, v, w hf C a L U0, V0 ν ξ k DT (ρCp)nf Nb DB (ρCp)p σ B0 βC ' U∞ βT Le
Temperature Wall shear along y- direction Ambient concentration Wall temperature Ambient temperature Wall concentration density Stretching velocity in x- and y- direction Velocity components Convective heat transfer coefficient Concentration Constant Characteristic Length Constant Kinematic viscosity Similarity variable Thermal conductivity Thermophoretic diffusion coefficient Specific heat capacity of nanofluid Brownian motion parameter Brownian diffusion coefficient Specific heat capacity of particle Electrical conductivity Magnetic field intensity Volumetric coefficient of solutal expansion Prime denotes derivative w.r.t. ξ Stream velocity Volumetric coefficient of thermal expansion Lewis number
Abbreviations
B.C’s Boundary Conditions D.E’s Differential Equations O.D.E’s Ordinary Differential Equations P.D.E’s Partial Differential Equations MHD Magnetohydrodynamics ANN Artificial Neural Network 3D Three Dimensional
Ethics
There are no ethical issues with the publication of this manuscript.
References
- Kumar SK, Muniamuthu S, Mohan A, Amirthalingam βC Volumetric coefficient of solutal expansion P, Anbu Muthuraja M. Effect of charging and dis- ' Prime denotes derivative w.r.t. ξ charging process of PCM with paraffin and Al2O3 U∞ Stream velocity additive subjected to three point temperature loca- βT Volumetric coefficient of thermal expansion tions. J Ecol Engineer 2022;23:34-42. [CrossRef] Le Lewis number
- Xuan Y, Li Q. Heat transfer enhancement of nano- Authors equally contributed to this work. fluids. Int J Heat Fluid Flow 2000;21:58-64. [CrossRef]
- Buongiorno J. Convective transport in nanofluids. J DATA AVAILABILITY STATEMENT Heat Transf 2006;128:240-250. [CrossRef]
- Shahid A, Bhatti MM, Ellahi R, Mekheimer KS. The authors confirm that the data that supports the Numerical experiment to examine activation findings of this study are available within the article. Raw energy and bi-convection Carreau nanofluid flow data that support the finding of this study are available from on an upper paraboloid porous surface: Application the corresponding author, upon reasonable request. in solar energy. Sustain Energy Technol Assess 2022;52:102029. [CrossRef]
- Din ISU, Siddique I, Ali R, Jarad F, Abdal S, Hussain CONFLICT OF INTEREST S. On heat and flow characteristics of Carreau nano- The author declared no potential conflicts of interest fluid and tangent hyperbolic nanofluid across a with respect to the research, authorship, and/or publication wedge with slip effects and bioconvection. Case Stud of this article. Ther Engineer 2022;39:102390. [CrossRef] 1118 J Ther Eng, Vol. 10, No. 5, pp. 1107−1119, September, 2024
- Ahmad M, Muhammad T, Ahmad I, Aly S. Time- dependent 3D flow of viscoelastic nanofluid over an heat transfer to a newtonian fluid over non-linear unsteady stretching surface. Phys A Stat Mech Appl extrusion stretching sheet. Int J Appl Comput Math 2020;551:124004. [CrossRef] 2018;4:35. [CrossRef]
- Acharya N. On the flow patterns and thermal behaviour of hybrid nanofluid flow inside a micro- of bio-convective Stefan blowing Ag-MgO/water channel in presence of radiative solar energy. J hybrid nanofluid induced by stretching cylin- Therm Anal Calorim 2020;141:1425-1442. [CrossRef] der utilizing non-Fourier and non-Fick's laws.
- Rana P, Bhardwaj A, Makkar V, Pop I, Gupta G. Nanomaterials 2021;11:1735. [CrossRef] Critical points and stability analysis in MHD [29] Kandasamy R, Muhaimin I, Ram NS, Prabhu KKS. radiative non‐Newtonian nanoliquid transport Thermal stratification effects on hiemenz flow of phenomena with artificial neural network predic- nanofluid over a porous wedge sheet in the presence tion. Math Method Appl Sci 2023;46:11726-11746. of suction/injection due to solar energy: Lie group [CrossRef] transformation. Transp Porous Media 2012;94:399-
- Lu D, Ramzan M, Mohammad M, Howari F, Chung 416. [CrossRef] JD. A thin film flow of nanofluid comprising carbon [30] Rana P. MHD convective heat transfer in the annu- nanotubes influenced by Cattaneo-Christov heat lus between concentric cylinders utilizing nanopar- flux and entropy generation. Coatings 2019;9:296. ticles and non-uniform heating. AIP Conf Proc [CrossRef] 2020;2214:020013. [CrossRef]
- Alamri SZ, Khan AA, Azeez M, Ellahi, R. Effects of mass transfer on MHD second grade fluid towards bioconvective Casson nanofluid fluid flow over a stretching cylinder: A novel perspective of Cattaneo- nonlinear convective stretching sheet in suspen- Christov heat flux model. Phys Lett A 2019;383:276- sion of gyrotactic microorganism. Heat Transf
- [CrossRef] 2020;49:2449-2471. [CrossRef]
- Rana P, Shukla N, Bég OA, Bhardwaj A. Lie group analysis of nanofluid slip flow with Stefan blow- A, Waqas M, Ganji D. Investigation of cross-fluid ing effect via modified Buongiorno's Model: flow containing motile gyrotactic microorganisms Entropy generation analysis. Differ Equ Dyn Syst and nanoparticles over a three-dimensional cylin- 2021;29:193-210. [CrossRef] der. Alex Engineer J 2020;59:3297-3307. [CrossRef]
- Rana P, Sharma PK, Kumar S, Makkar V, Mahanthesh B. Multiple solutions and stability analysis in MHD KS, Ismail AI, Saxena H. The role of quadratic-lin- non‐Newtonian nanofluid slip flow with convec- early radiating heat source with Carreau nano- tive and passive boundary condition: Heat transfer fluid and exponential space-dependent past a cone optimization using RSM‐CCD. J Appl Math Mech and a wedge: A medical engineering application 2023;104:e202200145. [CrossRef] and renewable energy. J Comput Biophys Chem
- Anwar MI, Shafie S, Hayat T, Shehzad SA, Salleh 2023;22:997-1011. [CrossRef] MZ. Numerical study for MHD stagnation-point [34] Rana P. Heat transfer optimization and rheological flow of a micropolar nanofluid towards a stretching features of Buongiorno nanofluid in a convectively sheet. J Braz Soc Mech Sci Engineer 2017;39:89-100. heated inclined annulus with nonlinear thermal [CrossRef] radiation. Propuls Power Res 2023;12:539-555.
- Shawky HM, Eldabe NTM, Kamel KA, Abd-Aziz [CrossRef] EA. MHD flow with heat and mass transfer of [35] Bahrami HR, Ghaedi M, Attarzadeh A. Employing Williamson nanofluid over stretching sheet through nonuniform magnetic fields to improve energy porous medium. Microsyst Technol 2019;25:1155- transfer of flow after a sudden expansion inside a 1169. [CrossRef] miniature channel: A hydrothermal study. Engineer
- Vajravelu K, Cannon J. Fluid flow over a nonlinearly Rep 2024;e12847. [CrossRef] stretching sheet. Appl Math Comput 2006;181:609- [36] Mebarek-Oudina F, Chabani I, Vaidya H, Ismail
- Matin MH, Nobari MRH, Jahangiri P. Entropy and porous media contribution to entropy gen- analysis in mixed convection MHD flow of nano- eration. Int J Numer Method Heat Fluid Flow fluid over a non-linear stretching sheet. J Ther Sci 2024;34:809-836. [CrossRef] Technol 2012;7:104-119. [CrossRef] [37] Chamkha J, Rashad A Unsteady heat and mass trans-
- Jain S, Choudhary R. Soret and dufour effects on fer by mhd mixed convection flow from a rotating thermophoretic MHD flow and heat transfer over a vertical cone with chemical reaction and Soret and non-linear stretching sheet with chemical reaction. Dufour effects. Can J Chem Engineer 2014;92:758- Int J Appl Comput Math 2018;4:50. [CrossRef] 767. [CrossRef] J Ther Eng, Vol. 10, No. 5, pp. 1107−1119, September, 2024 1119
- RamReddy CH, Murthy PVSN, Chamkha AJ, Int J Numer Method Heat Fluid Flow 2015;25:762- Rashad AM. Soret effect on mixed convection flow 781. [CrossRef] in a nanofluid under convective boundary con- [40] Abolbashari MH, Freidoonimehr N, Nazari F, Rashidi dition. Int J Heat Mass Transf 2013;64:384-392. MM. Analytical modeling of entropy generation for [CrossRef] Casson nano-fluid flow induced by a stretching sur-
- Hayat T, Muhammad T, Shehzad SA, Alsaedi A. face. Adv Powder Technol 2015;26:542-552. [CrossRef] Soret and dufour effects in three dimensional flow [41] Khan W, Pop I. Boundary-layer flow of a nano- over an exponentially stretching surface with porous fluid past a stretching sheet. Int J Heat Mass Transf medium, chemical reaction and heat source/sink. 2010;53:2477-2483. [CrossRef]
Share and Cite
GUPTA, S.; SHARMA, P.K.; KUMAR, S.; TIWARI, C.M. Influence of buoyancy forces in MHD non-Newtonian convective nanofluid utilizing Buongiornos Model i. Journal of Thermal Engineering 2024, Vol. 10, pp. 1107-1119. https://doi.org/10.14744/thermal.0000854

