Abstract
This article investigates the pulsatile flow of viscous incompressible MHD nanofluid in a rectangular channel. At the upper and lower walls, the channel has symmetrical constrictions. The goal is to analyze the heat transfer features of the nanofluid flow under the effect of the magnetic field and thermal radiation. Five different nanofluids, formed with nanoparticles of copper , magnetite , silver , titanium oxide , and single wall carbon nanotube in the base fluid of water, are considered in the study. The unsteady governing equations are transformed using the vorticity-stream function approach. The solution is obtained using the finite difference technique (FDM). The effect of various flow controlling parameters on velocity, temperature, Nusselt, and skin-friction profiles is inspected by using graphs. Across the channel, graphs of vorticity, streamlines, and temperature distribution are also displayed. The thickness of the thermal boundary layer upsurges with escalating values of the magnetic field, radiation parameters, and solid volume fraction, whereas it declines with escalating values of the Strouhal and Prandtl numbers. The profiles are usually found to have a more regular pattern upstream of the constriction than that downstream of the constriction. At the throat of the constriction, the carbon nanotube-based nanofluid attains higher temperatures than the other nanofluids downstream of the constriction. However, in the lee of the constriction, silver-based nanofluid attains higher temperatures than the other nanofluids downstream of the constriction. The behavior, in most cases, is opposite upstream of the constriction. The findings of the study can be utilized to cure stenosis in blood vessels, design biomechanical devices, and employ flow pulsation to control industrial operations.
1 Introduction
Conventional heat transfer liquid sources are incapable of satisfying current cooling requirements largely due to their poor convective heat transfer coefficients. Researchers have demonstrated that nanoparticles (NPs) usually made of metals or oxides enhance the coefficients of heat convection and conduction in fluids, allowing for higher heat transfer rates for the coolants. Thus, the heat transfer and thermal system effectiveness can be improved significantly by mixing NPs in pure fluid, forming nanofluid (NF). NFs have promising applications in many areas of industry and biomedicine due to their enhanced thermal conductivity. The characteristic of improved thermal conductivity might serve as a major factor for performance improvement. As NFs can enhance heat transfer, heat exchangers can be designed to be both energy efficient and small.
Wang and Mujumdar (2008) investigated the convective heat transfer of NF flow (NFF). Saidur et al. (2011a), Saidur et al. (2011b), Mahian et al. (2013), and Kasaeian et al. (2015) worked to enhance the thermophysical properties as well as heat transfer capacity of the fluids using NPs. The study of these flows through stretching surfaces got motivation from the perspective of its application, particularly in plastic film drawing. Hence, several researchers paid a lot of attention to this issue and studied the movement of boundary layers over different forms of stretching surfaces. Akbar et al. (2014) used a homogeneous model to study NFF at stagnation point above a stretching plate with slip boundary conditions. Analysis of pulsatile flow in a constricted channel under the impact of magnetic field was presented by Bandyopadhyay and Layek (2011) and Bandyopadhyay and Layek (2012). Nasir et al. (2019) explained the Darcy–Forchheimer 2D thin-film fluid of NF. Mustafa et al. (2011) considered the NFF at the stagnation point above a stretching sheet. Wong and Leon (2010) reported that the thermophysical properties of fluids are improved significantly, even at moderate NP concentrations. Haq et al. (2016) investigated the fully developed flow of water-functionalized magnetite NPs among two parallel disks by taking water as the base fluid. The peristaltic flow of incompressible viscous fluid having metallic NPs was examined by Akbar (2014)via an irregular duct. Aly (2020) explained a non-homogeneous two-phase model for –water NF-filled annulus used for the simulation of an incompressible smoothed particle hydrodynamic system (ISPH) between a wavelength rectangle and the square cavity. Two-dimensional –water NF under the combined impact of Lorentz and Kelvin forces was inspected by Sheikholeslami et al. (2017a). Said et al. (2015) analyzed the influence of short suspended SWCNT thermophysical properties in water and enhanced the thermal productivity of a flat plate solar collector. Sheikholeslami and Ganji (2013) explained the - NFF between parallel plates. Yang et al. (2020) explained air purification with total heat recovery using NFs for the first time, demonstrating the synergistic application of NFs for heat, mass transfer, photocatalysis, and sterilization.
In the magnetohydrodynamic (MHD) flow, the magnetic area induces an electric current in a moving conductive fluid. The induced current causes force on conductive fluid ions. Natural convection under the impact of the magnetic field remained a topic of great interest because of its comprehensive applications in the design of liquid-metal cooling systems, accelerators, and pumps, as well as MHD generators and flow meters (Cha et al., 2002). In a single constricted channel, Ali et al. (2020) investigated the steady and pulsatile modes of non-Newtonian MHD Casson fluid. Sheikholeslami (2017a) researched the hydrothermal study of MHD nanofluid in an open porous cavity by using the Lattice–Boltzmann method with the Brownian motion impact on NF properties. Shah et al. (2019) addressed the idea of the electrical MHD rotational flow of and for engine oils. Sheikholeslami et al. (2017b) researched the forced convection of MHD NFF. They considered the power of Brownian motion for modeling NF. Haq et al. (2017) supposed the heat transfer efficiency of engine oil between two dispersed cylinders with MHD effects in the presence of both and , to monitor the spontaneous motion of the NPs. There are several studies, as well, regarding MHD NF flow problems past a plate/sheet under various physical conditions and objectives (Narayana and Venkateswarlu, 2016; Babu et al., 2018; Tarakaramu et al., 2019; Devaki et al., 2020; Tarakaramu et al., 2020; Venkateswarlu and Narayana, 2021; Alzahrani et al., 2022; Mahmood et al., 2022; Puneeth et al., 2022; Qi et al., 2022; Ramadan et al., 2022; Tlili and Alharbi, 2022; Tlili et al., 2022).
Kakarantzas et al. (2009) investigated the natural convection of MHD at the sinusoidal upper wall temperature in a vertical cylindrical cavity. Rashidi et al. (2016) examined the convective heat transfer of MHD NFF in a vertical channel with sinusoidal walls. Sheikholeslami (2017b) tested the Buongiorno model to investigate the melting heat transfer of NFF due to the magnetic field. The impacts of heat and mass transfer of NFF flow over a vertical infinite flat plate were investigated by Turkyilmazoglu and Pop (2013). They derived the exact analytical solutions for various water-based NFs containing , , , , and . Pakdaman et al. (2012) explored the thermophysical properties and overall efficiency of NFs flowing within vertical helically coiled tubes. In the occurrence of a non-uniform magnetic field, Sheikholeslami et al. (2015) explored the transport of forced convection heat within a lid-driven semi-annulus enclosure filled with –water NF. Using a continuously stretching porous layer, Gopal et al. (2021) investigated the thermo-physical properties characteristics of complex order chemical processing and viscous dissipation on NF. Two space coordinates are used to model the porous medium, laminar, time-invariant, and MHD incompressible Newtonian NF. The role of heat source and Soret impacts on MHD convective Ferro-NF (-water) flow across an inclined channel with a porous media were examined by Sabu et al. (2021). The heat transfer of NFF from an inclined channel in the occurrence of a magnetic field, heat flux, and mass diffusion was explored using the FDM by Kumar et al. (2021).
Krishna et al. (2021a) and Krishna et al., 2021b) investigated the radiative MHD flow of an incompressible viscous electrically conducting non-Newtonian Casson hybrid NF over vertical moving porous surface under the influence of slip velocity in a rotating frame. Kavya et al. (2022) explore the varying fluid momentum and thermal energy characteristics of the laminar, steady, incompressible, 2D, non-Newtonian pseudo-plastic Williamson hybrid NF over a stretching cylinder with MHD effects.
Shah et al. (2022) determined the heat transfer properties of a MHD Prandtl hybrid NF over a stretched surface in the presence of bioconvection and chemical reaction effects. Zhang et al. (2022) presented the solar source, although it can reduce energy consumption (EC) for buildings on cold days; in the summer, its presence on the envelopes intensifies EC. Gao et al. (2022) examined the hydrodynamic and thermal performances of liquid water in the attendance of two carbon structures by molecular dynamics simulation (MDS). Nayak et al. (2022) conducted a numerical treatment on flow and heat transfer of radiative hybrid NF past an isothermal stretched cylinder set in a porous medium.
The goal of the present work is to numerically investigate the heat transfer characteristics of pulsatile viscous MHD NFF in a constricted channel. The study is conducted for five different nano-impurities (, , , , and ), which have been chosen for the preparation of NFs. The mathematical analysis is done using the transformation of the vorticity-stream function, and numerical computations are made using FDM. The impact of the flow controlling parameters on the flow profiles of the five NFs, as well as on the skin-friction and Nusselt number profiles, are also observed. The objective is to examine the cumulative impact of the applied magnetic field and thermal radiation on the wall shear stress (WSS), velocity, and temperature profiles in 2D. The Newtonian pulsatile NF flows help in understanding the influence of various metallic NPs homogeneously suspended in the blood, which is driven by drug trafficking (pharmacology) applications.
The following is how the article’s next section is organized. Section 2 defines the mathematical model and transformation. The outcomes and relevant discussions are presented in Section 3. The final remarks are eventually presented in Section 4.
2 Mathematical model
We are working at an incompressible electrically conducting NF flowing via a rectangular channel that should be laminar at Reynolds number 700. The channel walls have a pair of symmetrical constrictions. The resulting electric field is normal, and the direction of is perpendicular to the plane of flow. We take a Cartesian coordinate system (x,y) in which the flow direction and the direction of are, respectively, along the -axis and -axis. Because the magnetic Reynolds number for the flow is so small, the induced electric can be deemed insignificant. In the transformed coordinate system (to be discussed later on), as illustrated in Figure 1, the constrictions span from to , with the center at where and symbolize the lower and upper walls, respectively, and and are the heights of the constrictions at the lower and upper walls, respectively (Bandyopadhyay and Layek, 2012; Ali et al., 2020).
FIGURE 1
The flow phenomenon is represented by the unsteady incompressible viscous flow equations as follows.
The continuity equation:The momentum equation:The energy equation:where is the radiative heat flux. We get by expanding about and omitting higher-order terms
Then, and .
Eq. 5 becomeswhere (along -axis) and (along -axis) show the velocity components and the subscript “” represents the nanofluid, symbolizes the pressure, symbolizes the density, symbolizes the characteristic flow velocity, symbolizes the kinematic viscosity, and symbolizes the temperature. symbolizes the thermal conductivity, symbolizes the specific heat, symbolizes current density, symbolizes the magnetic field, symbolizes the strength of the uniform magnetic field, symbolizes electric conductivity, and symbolizes the magnetic permeability of the medium. As the electric current direction is normal to the channel plane, therefore . Ohm’s law gives
Maxwell’s equation (constant) for the steady flow. For the present study, is assumed to be zero. Then, Eq. 7 gives . Therefore, . Hence, Eq. 3 becomes
The dimensionless version of the governing model is created by considering the following quantities:
Here, symbolizes the maximum width of the channel, symbolizes the characteristic flow velocity, symbolizes the period of the pulse, symbolizes the Reynolds number, symbolizes the Strouhal number, symbolizes the Hartmann number, symbolizes the Prandtl number, and symbolizes the radiation parameter.
The effective dynamic viscosity , effective density , heat capacitance , and effective thermal conductivity of a NF are given as:where symbolizes the volume fraction of NF, and the subscripts correspond to NF, corresponds to the pure fluid, and corresponds to the solid state. For the present work, the following five different NFs are considered:-, -, -, -, and -. The thermophysical properties of the five different NPs and water are shown in Table 1 (Sheikholeslami and Ganji, 2013; Sheikholeslami et al., 2017a). Using the quantities from Eqs 9, 10 in Eqs 2–4, 6 giveswhere , , , , .
TABLE 1
| Physical properties | Base Fluid | |||||
|---|---|---|---|---|---|---|
| (J/kg K) | 4076.4 | 385 | 235 | 670 | 686.2 | 600 |
| (Kg/) | 997.8 | 8933 | 10,500 | 5180 | 4250 | 2100 |
| K (W/m K) | 0.60475 | 401 | 429 | 9.7 | 8.9538 | 3500 |
Thermophysical properties of different NPs and the base fluid.
2.1 Boundary conditions
For the steady case of the flow problem under consideration, the boundary conditions, obtained by solving the dimensionless form of Eq. 8 and performing some manipulations using Eq. 10, are given as follows:When , the -velocity at the inlet takes the form:The flow is classified as sinusoidal for the pulsatile flow:Furthermore, and (i.e., no-slip conditions) are assumed on the walls. For the outlet boundary, the fully developed flow conditions are incorporated. In the dimensionless form, the temperature condition at the lower wall, and at the upper wall, .
2.2 Vorticity-stream function formulation
We differentiate Eqs 12, 13 with respect to and , respectively. Then, from their subtraction and using Eq. 18, we get transformed equations from the variables to vorticity-stream functions as follows:And the stream function equation (Poisson equation) is given as
2.3 Transformation of coordinates
For treating the channel walls as the straighten ones, the coordinates are transformed asThus, the lower wall symbolizes by and the upper wall symbolizes by . In the updated coordinate system , Eqs 14, 19, 20 arewhereThe velocity components and in terms of take the formsFurthermore, the wall boundary conditions in system for and are given bywhere symbolizes the pulsating amplitude. If , we get steady flow
The temperature’s converted boundary conditions are , at ; , at .
The other concerning non-dimensional physical quantities comprise the skin-friction coefficient and Nusselt number, defined bywhere and are defined as
By using dimensionless variables from Eq. 9 and the coordinate transformation from Eq. 21, we getwhere and .
3 Results and discussion
The problem Eqs 22–24 is computed using a numerical scheme based on FDM, subject to the relevant boundary conditions in Sections 2.1, 2.3. The numerical scheme adopts a standard approach, as used by Bandyopadhyay and Layek (2011), Bandyopadhyay and Layek (2012), and Ali et al. (2020). The computational domain is assumed as . The resulting linear systems are solved using TDMA during each of the ADI process’s two half-steps. A Cartesian grid of elements is formed for and , with and . For time integration, a step size of is considered. The constriction length is assumed as The constrictions on the lower and upper walls are assumed to have heights ; hence, at the constricted part of the channel, the minimum width is still 30% of the maximum channel width . The computations for the current study are performed in a sequential fashion. The results can be found by parallel computing for time-efficient solutions (Ali and Syed, 2013).
The pulsatile motion is modeled by adding in the inflow boundary condition the sinusoidal time-dependent function . For comparison of the five NFs, the effects of the physical parameters such as , , , and on the existence of the dimensionless streamwise velocity and temperature are analyzed. We perform simulations for a long enough time but display the results graphically only for , the point at which the flow rate is at its maximum and taking the location (throat of the constriction) as well as (in the lee of the constriction) where the fluid has entered in the low-pressure zone from the high-pressure zone. For a pulse cycle, is the acceleration phase, and is the deceleration phase.
For validation, the present results for the pulsatile flow condition are compared with those obtained by Bandyopadhyay and Layek (2012) without the heat effect. Figure 2 shows a good agreement of the present results, specifically for the wall shear stress (WSS), with (Bandyopadhyay and Layek, 2012) for at .
FIGURE 2
The effects of varying the volume fraction on the wall shear stresses (WSS) at the upper wall are computed for with , , , and of the five NFs. The results are shown in Figure 3 at and . The WSS is maximum at as the flow rate is maximum. The WSS decreases for -, -, and - on increasing , whereas the WSS slightly increases for - and - on increasing .
FIGURE 3
The velocity and temperature profiles for each of the five NFs for are shown in Figures 4, 5 and Figures 6, 7, respectively, by setting the other parameters as , , and . The slope/gradient of the velocity profile rises with rising values of , i.e., the viscous boundary layer declines with rising values of . In Figure 4, profiles at and are shown. It has been noted that the peak value of rises with and exhibits a parabolic profile at . The velocity reaches its maximum at the middle of the constriction. -based NF attains the peak value of the velocity higher than that of the other NFs. In Figure 5, profiles at and are shown. The profiles are not parabolic as some backflow in the vicinity of the walls is observed. The backflow reduces with an increase in . -based NF attains the peak value of velocity higher than that of the other NFs. In Figures 6, 7, the temperature profiles for are shown at and , respectively. It is seen that the thickness of the thermal boundary layer increases, resulting in the rise of the temperature as is increased. The temperature distribution across the channel at for distinct values of is shown in Figure 8.
FIGURE 4
FIGURE 5
FIGURE 6
FIGURE 7
FIGURE 8
The temperature profiles for each of the five NFs for are computed by setting the other parameters as , , and . In Figures 9, 10, the temperature profiles for are shown at and , respectively. It is seen that the thickness of the thermal boundary layer declines, resulting in a decline of the temperature profile as the rises.
FIGURE 9
FIGURE 10
The temperature profiles for each of the five NFs for are computed by setting the other parameters as , , and . In Figures 11, 12, the temperature profiles for are shown at and , respectively. It is seen that the thickness of the thermal boundary layer declines, resulting in a decline of the temperature profile as the increases. The temperature distribution across the channel at for distinct values of is shown in Figure 13.
FIGURE 11
FIGURE 12
FIGURE 13
The temperature profiles for each of the five NFs for are computed by setting the other parameters as , , and . In Figures 14, 15, the temperature profiles for are shown at and , respectively. It is seen that the thickness of the thermal boundary layer rises, resulting in the rise of the temperature as rises. The temperature distribution across the channel at for distinct values of is shown in Figure 16.
FIGURE 14
FIGURE 15
FIGURE 16
Nearly for all of the cases discussed above, for the variation of any of , , , and , a common behavior can be noticed. At the throat of the constriction, -based NF attains higher temperatures than the other NFs downstream of the constriction. In the lee of the constriction, -based NF attains higher temperatures than the other NFs downstream of the constriction. The behavior, in most cases, is opposite upstream of the constriction.
The velocity and temperature profiles for each of the five NFs for the four selected time levels: , are shown in Figures 17, 18, respectively, by setting the parameters as , , , and . The selected time levels are related to the start of pulsatile motion, the maximum flow rate, the minimum flow rate, and the instant zero net flow, respectively. In Figure 17, it is observed that the streamwise velocity achieves its peak value when the flow rate is maximum during the pulse cycle. The profiles are parabolic at . It is found that the profile is symmetric for all . -based NF attains the peak value of the velocity higher than that of the other NFs. In Figure 18, it is detected that during the acceleration () and deceleration ) phases of the pulsation cycle, - and -based NFs attain higher temperature values than that of the other NFs downstream of the constriction. Whereas, for the same phases, -based NF attains higher temperature values than that of the other NFs upstream of the constriction.
FIGURE 17
FIGURE 18
The effects of the variation of the solid volume fraction for and on the temperature profiles are presented in Figure 19 for -based NF. It is evident that the heat transfer rate is higher for the NF with a higher concentration of NPs. Moreover, the temperature gradients are higher near the lower wall of the channel.
FIGURE 19
The effects of , , and on the Nusselt number, as well as skin-friction profiles, are shown in Figures 20, 21, respectively. It is observed that the Nusselt number rises for higher values of , , and . The Nusselt number grows as the flow becomes more turbulent due to the rising number of collisions among the fluid particles. Also, the skin friction coefficient as the values of upsurges. However, it remains unchanged for and .
FIGURE 20
FIGURE 21
The streamlines, vorticity, and temperature distribution plots at distinct instants of a pulse cycle with , , , , and are shown in Figure 22 for -based NF. In Figure 22, the formation of vertical eddies in the vicinity of the walls can be observed. Over time, the eddies grow for a specific value of and slowly occupy a greater part of the channel downstream. At , the presence of the backflow is observed, and even before the constriction bump, the vertical eddies are formed. The formation of these vertical eddies and flow separation is essential. However, these can be handled by increasing the strength of the magnetic field.
FIGURE 22
4 Concluding remarks
In this research, the numerical analysis of the pulsatile flow of five different NFs in a channel impacted by the magnetic field and thermal radiation through a rectangular channel having constricted walls is presented. The NPs include
,
,
,
, and
that are mixed in the base fluid, water
, to form the five NFs. The impacts of each of
,
,
, and
on the flow profiles are studied. In this analysis, the following major outcomes were observed:
• The streamwise velocity escalates as the values of upsurge. The slope/gradient of the velocity profile is higher for higher values of , i.e., the viscous boundary layer declines with rising values of . At the throat of the constriction, - attains the velocity’s peak value higher than that of the other NFs. However, in the lee of the constriction, - attains the velocity’s peak value higher than that of the other NFs.
• The thickness of the thermal boundary layer rises as any of , , and is increased. This results in the rise of temperature.
• The thickness of the thermal boundary layer declines as any of , and is increased. This results in a decline in the temperature profile.
• At the throat of the constriction, - attains higher temperatures than the other NFs downstream of the constriction.
• In the lee of the constriction, - attains higher temperatures than the other NFs downstream of the constriction. The behavior, in most cases, is opposite upstream of the constriction.
• The streamwise velocity attains its maximum at the middle of the constriction at during the pulse cycle. - attains the velocity’s peak value higher than that of the other NFs.
• During the acceleration and deceleration phases of the pulsation cycle, - and -based NFs attain higher temperature values than that of the other NFs downstream of the constriction. Whereas, for the same phases, -based NF attains higher temperature values than that of the other NFs upstream of the constriction.
• The heat transfer rate is higher for the NF with a higher concentration of the NPs.
• The skin-friction coefficient escalates with escalating values of . However, no difference is found in the case of variation in and .
Statements
Data availability statement
The original contributions presented in the study are included in the article/supplementary material; further inquiries can be directed to the corresponding author.
Author contributions
AA: conceptualization; ZB: methodology; MA: validation; S.A.: Coding; ESM Tag El Din: Analysis; SM Hussain: Data Curation.
Funding
The authors are grateful to the Deanship of Scientific Research, Islamic University of Madinah, Ministry of Education, KSA, for supporting this research work through research project grant under Research Group Program/1/804.
Acknowledgments
The authors are thankful to the reviewer for the comments and suggestions. The authors highly appreciate the reviewer for the high standard review.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
AkbarN. S.KhanZ. H.NadeemS. (2014). The combined effects of slip and convective boundary conditions on stagnation-point flow of CNT suspended nanofluid over a stretching sheet. J. Mol. Liq.196, 21–25. 10.1016/j.molliq.2014.03.006
2
AkbarN. S. (2014). Metallic nanoparticles analysis for the peristaltic flow in an asymmetric channel with MHD. IEEE Trans. Nanotechnol.13 (2), 357–361. 10.1109/tnano.2014.2304362
3
AliA.FarooqH.AbbasZ.BukhariZ.FatimaA. (2020). Impact of Lorentz force on the pulsatile flow of a non-Newtonian Casson fluid in a Constricted channel using Darcy’s law: A numerical study. Sci. Rep.10, 10629. 10.1038/s41598-020-67685-0
4
AliA.SyedK. S. (2013). Adv. Comp.91, 87–118.
5
AlyA. M. (2020). ZAMM J. App. Math. Mech.100 (9).
6
AlzahraniJ.VaidyaH.PrasadK. V.RajashekharC.MahendraD. L.TliliI. (2022). Micro-polar fluid flow over a unique form of vertical stretching sheet: Special emphasis to temperature-dependent properties. Case Stud. Therm. Eng.34, 102037. 10.1016/j.csite.2022.102037
7
BabuH.AjmathK. A.VenkateswarluB.NarayanaP. V. S. (2018). Thermal radiation and heat source effects on MHD non-Newtonian nanofluid flow over a stretching sheet. J. nanofluids8, 1085–1092. 10.1166/jon.2019.1666
8
BandyopadhyayS.LayekG. C. (2011). Numerical computation of pulsatile flow through a locally constricted channel. Commun. Nonlinear Sci. Numer. Simul.16, 252–265. 10.1016/j.cnsns.2010.03.017
9
BandyopadhyayS.LayekG. C. (2012). Study of magnetohydrodynamic pulsatile flow in a constricted channel. Commun. Nonlinear Sci. Numer. Simul.17, 2434–2446. 10.1016/j.cnsns.2011.09.040
10
ChaJ. E.AhnY. C.KimM. H. (2002). Flow measurement with an electromagnetic flowmeter in two-phase bubbly and slug flow regimes. Flow. Meas. Instrum.12 (5-6), 329–339. 10.1016/s0955-5986(02)00007-9
11
DevakiP.VenkateswarluB.SrinivasS.SreenadhS. (2020). Nonlin. Eng.9, 51–59.
12
GaoJ.LiuJ.YueH.ZhaoY.TliliI.KarimipourA. (2022). Effects of various temperature and pressure initial conditions to predict the thermal conductivity and phase alteration duration of water based carbon hybrid nanofluids via MD approach. J. Mol. Liq.351, 118654. 10.1016/j.molliq.2022.118654
13
GopalD.SaleemS.JagadhaS.AhmadF.AlmatroudA. O.KishanN. (2021). Numerical analysis of higher order chemical reaction on electrically MHD nanofluid under influence of viscous dissipation. Alexandria Eng. J.60, 1861–1871. 10.1016/j.aej.2020.11.034
14
HaqR. U.NoorN. F. M.KhanZ. H. (2016). Adv. Powder Technol.27 (4), 1568–1575.
15
HaqR. U.ShahzadF.AlmdallalQ. (2017). MHD pulsatile flow of engine oil based carbon nanotubes between two concentric cylinders. Results Phys.7, 57–68. 10.1016/j.rinp.2016.11.057
16
KakarantzasS. C.SarrisI. E.GrecosA. P.VlachosN. S. (2009). Magnetohydrodynamic natural convection in a vertical cylindrical cavity with sinusoidal upper wall temperature. Int. J. Heat. Mass Transf.52, 250–259. 10.1016/j.ijheatmasstransfer.2008.06.035
17
KasaeianA.EshghiA. T.SametiM. (2015). A review on the applications of nanofluids in solar energy systems. Renew. Sustain. Energy Rev.43, 584–598. 10.1016/j.rser.2014.11.020
18
KavyaS.NagendrammaV.AhammadN. A.AhmadS.RajuC. S. K.ShahN. A. (2022). Magnetic-hybrid nanoparticles with stretching/shrinking cylinder in a suspension of MoS4 and copper nanoparticles. Int. Commun. Heat Mass Transf.136, 106150. 10.1016/j.icheatmasstransfer.2022.106150
19
KrishnaM. V.AhammadN. A.ChamkhaA. J. (2021). Ain. Sh. Eng. J.12 (3), 3043–3056.
20
KrishnaM. V.AhammadN. A.ChamkhaA. J. (2021). Radiative MHD flow of Casson hybrid nanofluid over an infinite exponentially accelerated vertical porous surface. Case Stud. Therm. Eng.27, 101229. 10.1016/j.csite.2021.101229
21
KumarM. A.ReddyY. D.RaoV. S.GoudB. S. (2021). Thermal radiation impact on MHD heat transfer natural convective nano fluid flow over an impulsively started vertical plate. Case Stud. Therm. Eng.24, 100826. 10.1016/j.csite.2020.100826
22
MahianO.KianifarA.KalogirouS. A.PopI.WongwisesS. (2013). A review of the applications of nanofluids in solar energy. Int. J. Heat. Mass Transf.57, 582–594. 10.1016/j.ijheatmasstransfer.2012.10.037
23
MahmoodR. T.AsadM. J.HadriS. H.El-ShorbagyM. A.MousaA. A. A.et al (2022). Hum. Ecol. risk assess.
24
MustafaM.HayatT.PopI.AsgharS.ObaidatS. (2011). Stagnation-point flow of a nanofluid towards a stretching sheet. Int. J. Heat. Mass Transf.54, 5588–5594. 10.1016/j.ijheatmasstransfer.2011.07.021
25
NarayanaP. V. S.VenkateswarluB. (2016). Front. Heat. Mass Trans.7.
26
NasirS.ShahZ.IslamS.BonyahE.GulT. (2019). Darcy Forchheimer nanofluid thin film flow of SWCNTs and heat transfer analysis over an unsteady stretching sheet. AIP Adv.9, 015223–15310. 10.1063/1.5083972
27
NayakM. K.MaboodF.DogonchiA. S.RamadanK. M.TliliI.KhanW. A. (2022). Wave random complex.
28
PakdamanM. F.Akhavan-BehabadiM. A.RaziP. (2012). Exp. Therm. Fluid Sci.40, 103–111.
29
PuneethV.AliF.KhanM. R.AnwarM. S.AhammadN. A. (2022). Bio. Conv. Ref.21.
30
QiX.SidiM. O.TliliI.IbrahimT. K.ElkotbM. A.El-ShorbagyM. A.et al (2022). Optimization and sensitivity analysis of extended surfaces during melting and freezing of phase changing materials in cylindrical Lithium-ion battery cooling. J. Energy Storage51, 104545. 10.1016/j.est.2022.104545
31
RamadanK. M.QisiehO.TliliI. (2022). Thermal creep effects on fluid flow and heat transfer in a microchannel gas cooling. Proc. Institution Mech. Eng. Part C J. Mech. Eng. Sci.236, 5033–5047. 10.1177/09544062211057039
32
RashidiM. M.NasiriM.KhezerlooM.LaraqiN. (2016). Numerical investigation of magnetic field effect on mixed convection heat transfer of nanofluid in a channel with sinusoidal walls. J. Magn. Magn. Mat.401, 159–168. 10.1016/j.jmmm.2015.10.034
33
SabuA. S.MathewA.NeethuT. S.GeorgeK. A. (2021). Statistical analysis of MHD convective ferro-nanofluid flow through an inclined channel with hall current, heat source and soret effect. Therm. Sci. Eng. Prog.22, 100816. 10.1016/j.tsep.2020.100816
34
SaidZ.SaidurR.SabihaM. A.RahimN. A.AnisurM. R. (2015). Thermophysical properties of Single Wall Carbon Nanotubes and its effect on exergy efficiency of a flat plate solar collector. Sol. Energy115, 757–769. 10.1016/j.solener.2015.02.037
35
SaidurR.KaziS. N.HusainM. S.RahmanM. M.MohammedH. A. (2011). Renew. Sust. Energy Rev.15, 310–323.
36
SaidurR.LeongK. Y.MohammedH. A. (2011). A review on applications and challenges of nanofluids. Renew. Sustain. Energy Rev.15 (3), 1646–1668. 10.1016/j.rser.2010.11.035
37
ShahS. A. A.AhammadN. A.DinE. M. T. E.GamaounF.AwanA. U.AliB., Nanomater. 12, 2174, (2022).
38
ShahZ.BonyahE.IslamS.GulT. (2019). Impact of thermal radiation on electrical MHD rotating flow of Carbon nanotubes over a stretching sheet. AIP Adv.9 (1), 015115–15213. 10.1063/1.5048078
39
SheikholeslamiM.GanjiD. D. (2013). Heat transfer of Cu-water nanofluid flow between parallel plates. Powder Technol.235, 873–879. 10.1016/j.powtec.2012.11.030
40
SheikholeslamiM.GanjiD. D.MoradiR. (2017). Forced convection in existence of Lorentz forces in a porous cavity with hot circular obstacle using nanofluid via Lattice Boltzmann method. J. Mol. Liq.246, 103–111. 10.1016/j.molliq.2017.09.053
41
SheikholeslamiM.HayatT.AlsaediA. (2017). Numerical study for external magnetic source influence on water based nanofluid convective heat transfer. Int. J. Heat. Mass Transf.106, 745–755. 10.1016/j.ijheatmasstransfer.2016.09.077
42
SheikholeslamiM. (2017). Influence of magnetic field on nanofluid free convection in an open porous cavity by means of Lattice Boltzmann method. J. Mol. Liq.234, 364–374. 10.1016/j.molliq.2017.03.104
43
SheikholeslamiM. (2017). Numerical simulation of magnetic nanofluid natural convection in porous media. Phys. Lett. A381, 494–503. 10.1016/j.physleta.2016.11.042
44
SheikholeslamiM.RashidiM. M.GanjiD. D. (2015). Effect of non-uniform magnetic field on forced convection heat transfer ofFe3O4–water nanofluid. Comput. Methods Appl. Mech. Eng.294, 299–312. 10.1016/j.cma.2015.06.010
45
TarakaramuN.NarayanaP. V. S.VenkateswarluB. (2019). App. Math. Sci. Comp.87-97.
46
TarakaramuN.NarayanaP. V. S.VenkateswarluB. (2020). Numerical simulation of variable thermal conductivity on 3D flow of nanofluid over a stretching sheet. Nonlinear Eng.9, 233–243. 10.1515/nleng-2020-0011
47
TliliI.AlharbiT. (2022). Investigation into the effect of changing the size of the air quality and stream to the trombe wall for two different arrangements of rectangular blocks of phase change material in this wall. J. Build. Eng.52, 104328. 10.1016/j.jobe.2022.104328
48
TliliI.SajadiS. M.BaleanuD.GhaemiF. (2022). Sustain. Energy Technol. Assess.52, 102100.
49
TurkyilmazogluM.PopI. (2013). Heat and mass transfer of unsteady natural convection flow of some nanofluids past a vertical infinite flat plate with radiation effect. Int. J. Heat. Mass Transf.59, 167–171. 10.1016/j.ijheatmasstransfer.2012.12.009
50
VenkateswarluB.NarayanaP. V. S. (2021). Heat. Trans.50, 432–449.
51
WangX. Q.MujumdarA. (2008). A review on nanofluids - part II: Experiments and applications. Braz. J. Chem. Eng.25 (4), 631–648. 10.1590/s0104-66322008000400002
52
WongK. V.LeonO. D. (2010). Adv. Mech. Eng.2, 1–11.
53
YangL.HuangJ. N.JiW.MaoM. (2020). Investigations of a new combined application of nanofluids in heat recovery and air purification. Powder Technol.360, 956–966. 10.1016/j.powtec.2019.10.053
54
ZhangJ.SajadiS. M.TliliI.YangChen, I.Tlili and M. A. Fagiry (2022). Effects of Al2O3 and TiO2 nanoparticles in order to reduce the energy demand in the conventional buildings by integrating the solar collectors and phase change materials. Sustain. Energy Technol. Assessments52, 102114. 10.1016/j.seta.2022.102114
Summary
Keywords
pulsatile flow, constricted channel, heat transfer analysis, nanofluid, carbon nanotube
Citation
Ali A, Bukhari Z, Amjad M, Ahmad S, Tag El. Din ESM and Hussain SM (2022) Newtonian heating effect in pulsating magnetohydrodynamic nanofluid flow through a constricted channel: A numerical study. Front. Energy Res. 10:1002672. doi: 10.3389/fenrg.2022.1002672
Received
25 July 2022
Accepted
22 August 2022
Published
29 September 2022
Volume
10 - 2022
Edited by
Adnan, Mohi-ud-Din Islamic University, Pakistan
Reviewed by
Iskander Tlili, National School of Engineers of Monastir, Tunisia
N. Ameer Ahammad, University of Tabuk, Saudi Arabia
Updates
Copyright
© 2022 Ali, Bukhari, Amjad, Ahmad, Tag El. Din and Hussain.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Muhammad Amjad, muhammadamjad@cuivehari.edu.pk
This article was submitted to Process and Energy Systems Engineering, a section of the journal Frontiers in Energy Research
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.