Significance of non-uniform heat source/sink and cattaneo-christov model on hybrid nanofluid flow in a Darcy-forchheimer porous medium between two parallel rotating disks

Abstract

The suspension of nanoparticles in fluid influences several properties of the resulting fluid. Many production and manufacturing applications need knowledge of the heat transference mechanism in nanofluids. The current paper concerns the influence of non-uniform heat source/sink on (MoS₂-Go/water flow) hybrid nanofluid flow and (Go/water flow) nanofluid flow in a Darcy-Forchheimer porous medium between two parallel and infinite spinning disks in the occurrence of radiation. The Cattaneo-Christov model is utilized to analyze heat and mass transmission. The Cattaneo-Christov model introduces the time lag factors in the process of heat and mass transmission, known as the thermal relaxation parameter and solutal relaxation parameter, respectively. The governing equations are numerically solved employing the “bvp4c function in MATLAB.” The effect of the primary relevant parameters on the velocity, temperature, nanoparticle concentration, and is graphically depicted. Finally, a table is drawn to show the relationships of various critical factors on the Nusselt number, and Sherwood number. Results reveal that an increase in the thermal relaxation parameter reduces the heat transmission rate at both the upper and lower plate. Furthermore, an increase in the nanoparticle’s volume fraction causes enhancement in thermal conduction, which increases the heat transmission rate at the upper disk. The results of this study will be helpful to many transportation processes, architectural design systems, enhanced oil recovery systems, medical fields that utilize nanofluids, and so on.

1 Introduction

To address today’s escalating energy demands, the variety of industries has risen exponentially. Optimizing heating performance is an essential alternative for improving cooling and heating efficiency in nuclear and chemical reactors, electronic equipment, and so on. To extract the most out of these applications, scientists and researchers are very much concerned about boosting heating efficiency. Turkyilmazoglu (2022a) studied the influence of a horizontal uniform magnetic field in the flow instigated by a rotating disk. In another study, Turkyilmazoglu (2022b) investigated the influence of heat source/sink in the convective flow of fluids due to a vertical flat surface/cone immersed in porous media. Ali et al. (2019a) utilized the Bingham and Carreau models to investigate the complex rheology of slime. Asghar et al. (2020a) utilized the couple stress fluid model to study the movement of microorganisms in two-dimensional channel. Ali et al. (2019b) investigated the optimum speed of microorgamisms through Carreau fluid models with under magnetic and porous effects. Javid et al. (2019) studied and gave the numerical solution of magnetically induced flow of fluid confined inside two curled peristaltic walls. Asghar et al. (2020b) studied the non-Newtonian fluid flow confined within a complex wavy walls of a 2-dimensional channel using Taylor’s swimming sheet model. Asghar et al. (2022a) studied the non-Newtonian Couple stress fluid in 2-dimensional inclined channel with magnetic and electrical field. Some other remarkable studies on applications of fluid models on bio-medical fields can be refered from Refs. (Asghar et al., 2020c; Asghar et al., 2022d; Asghar et al., 2022b; Asghar et al., 2022c; Asghar et al., 2023; Wu et al., 2020; Shah et al., 2022). Two decades-long, exploration of nanofluids (Choi, 1995) has confirmed that they have superior heat transmission potential. The varied characteristics of nanoparticles (NPs) may be adjusted by modifying their size (diameter), substance (metallic oxides, non-metallic, metallic, etc.), and nanoparticle dispersion in the working fluid. Keeping in mind the properties, Rahman et al. (2022) studied the influence of suction and magnetic field on several water-based nanofluids over a decelerating rotating disk. They studied the different water-based nanofluids (NFs) with copper oxide, copper, alumina, silver, and titania. The primary issue with single nanoparticle NFs is that they have either better rheological characteristics or superior thermal networks. Mono nanoparticles lack all of the desirable characteristics needed for some applications. The features of nanofluids may be improved by adjusting the nanoparticle volume fraction; however, this has a limitation as the difficulty arises in the trade-off due rise in viscosity. This is a constraint, which may be overcome by combining more than one variety of NPs in the working fluid. Later, the researcher’s developed hybrid nanofluids (HNFs) that optimize the exclusive features of many varieties of nanoparticles.

Many studies have used hybrid nanofluid combinations of metallic and non-metal oxide nanoparticles to maximize thermal performance and nanoparticle stability. The development of nanoclusters increases the relative viscosity of HNFs, and HNFs have superior thermal conductivity than solitary nanoparticle nanofluids and base fluids (Ranga Babu et al., 2017). Ranga Babu et al. (2017) also marked out that metal nanoparticles make a nanolayer atop metallic oxide nanoparticles, and it produces a thermal interfacial layer between weak boundaries of the working fluid and hybrid NPs, resulting in significant thermal conductivity augmentation. Devi and Devi (2016) studied the significance of Newtonian heating in HNF flow over a three-dimensional stretched surface. They analyzed the comparative behavior of Cu–Al₂O₃/water and Cu/water. They concluded that HNF has a higher heat transmission rate (HTR) compared to NF. Waqas et al. (2021) also explained and gave the solution to a HNF flow problem over a rotating disk with non-linear radiation. They studied the model with water as a base fluid and SWCNT− TiO₂ and MWCNT− CoFe₂O₄ NPs. Khan et al. (2020c) expounded on the HNF flow over a thick moving surface instigated by mixed convection. They modeled the flow with water as a main fluid and SiO₂ and MoS₂ NPs. Masood et al. (2021) discussed the significance of a HNF flow over a stretched surface near a stagnation point. They studied the problem with heat generation/absorption and water as a working fluid with polystyrene and titanium oxide NPs. Alrabaiah et al. (2022) inspected the HNF flow inside a conical slit of a cone and a disk. They modeled the flow with water and magnesium oxide, and silver NPs. Yaseen et al. (2022c) expounded on the significance of the magnetic field on HNF flow and NF flow between two parallel plates. They studied the flow problem with MoS₂–SiO₂/H₂O–C₂H₆O₂ for HNF flow and MoS₂/H₂O for NF flow. Qureshi et al. (2021) investigated the HNF flow between two moving co-axial orthogonal disks. They considered the flow problem with water as a main fluid and titania and copper NPs and studied the significance of the magnetic field. Rashid et al. (2021) inspected the HNF flow over a cylinder in motion. They considered the flow problem with water as a main fluid and Titania and silver NPs. Khan et al. (2022) explicated the significance of suction in a HNF flow between two parallel plates. They considered the flow problem with water as a base fluid and carbon nanotubes and Fe₃O₄ NPs.

Technological advancements have considerably enhanced human interaction with fundamental tasks such as heating and cooling food and materials, transportation, and manufacturing. The heat altercation process in equipment or applications intended to suit social requirements is a key aspect typically seen in the aforementioned activities. This has been a significant issue for manufacturers over the decades, particularly recently, in improving the thermal control of various maneuvers in the electronic sector, power systems, thermal sector, and medicinal bids. For a long period, the Fourier law (Fourier, 1822) was used to analyze the heat transmission attributes. Years later, Cattaneo (Cattaneo, 1948) altered the Fourier law by including a parameter related to the time lag in the traditional Fourier law, which states that the heat transfer mechanism permits heat to be carried at a restricted pace through the transmission of heat waves. Later, Christov (Christov, 2009) used thermal relaxation time and the upper convective derivative of Oldroyd to overcome the limitation of the Cattaneo rule and arrive at a material invariant formulation. The developments by Cattaneo (Cattaneo, 1948) and Christov (Christov, 2009) have led the way for researchers to study the HTR with a time lag factor and it is called the Cattaneo-Christov heat flux model (CCM). Very recently, Turkyilmazoglu (2021) provided an analytical explanation of the application of CCM in cooling and its role in enhancing the HTR from surfaces. In recent times; Khan et al. (2020b) explicated the consequence of the thermal stratification and CCM in the water-based NF flow over a surface having an exponential stretching rate. They analyzed the behavior of NF with three diverse NPs, namely, copper, Titania, and alumina; Yaseen et al. (2022a) studied the difference in HTR of MoS₂/Kerosene oil and MoS₂-SiO₂/Kerosene oil flow between two spinning disks with CCM. They concluded that at the lower disk, the HTR of the hybrid nanofluid transcends the HTR of the nanofluid. (Ashraf and Ghehsareh (2021) explored the significance of CCM in the flow of Casson fluid past a heated surface; Kumar et al. (2021) studied the flow of water-based NF with carbon nanotubes (multi and single-walled) between two spinning disks in a porous medium with CCM. Bai et al. (2022) derived the analytical solution of an Oldroyd-B nanofluid flow problem through a vibrating tensile plate with CCM. Naz et al. (2022) also discussed the analytical solution of a Carreau nanoliquid flow problem over a rotating and stretchable disk with CCM. Some more interesting research on CCM can be found in Refs. (Irfan et al., 2018; Garia et al., 2021; Khan et al., 2021; Madhukesh et al., 2021).

According to the preceding discussion, it can be inferred that only a few studies are published to study the significance of HNF between two spinning disks. Authors have developed a model to analyze the HNF flow (MoS₂-Go/water flow) and NF flow (Go/water flow) between two parallel and infinite spinning disks. Furthermore, as a novelty, the significance of non-uniform heat source/sink, CCM, thermal radiation, magnetic field, and Darcy-Forchheimer porous medium are also studied. This research paper also aims to analyze the HTR and mass transmission rate (MTR) of HNF flow (MoS₂-Go/water flow) and NF flow (Go/water flow) at both disks. The numerical solution is sought via the bvp4c solver in MATLAB. As per the author’s information, such a comparative analysis of the HNF flow (MoS₂-Go/water) and NF flow (Go/water) between two parallel and infinite spinning disks has never been reported before and the results are new and novel.

2 Flow model and governing equations

2.1 Assumptions of the problem

Consider the three-dimensional incompressible, steady, and axisymmetric flow of the HNF (MoS₂-Go/water) and NF (Go/water) between the two spinning infinite disks placed along parallel lines (see Figure 1). The notation has been used for the cylindrical coordinate system (see Figure 1). The flow is exposed to a magnetic field B_o in the z-direction. As a novelty, to develop the model, the Darcy-Forchheimer law is used for porous medium. In addition, the non-uniform heat source/sink, Cattaneo-Christov double diffusion model, and thermal radiation is used to investigate HTR and MTR. Moreover, the last term in the energy equation (Eq. 5) denotes the non-uniform heat source/sink term and it is explained later. The lower disk and upper disk are placed at and (see Eq. 7). In the boundary conditions, a₁ and a₂ denote the shrinkage rates of the lower and upper disks in the radial direction, respectively. Furthermore, b₁ and b₂ are taken as the angular velocities of the lower and upper disks, respectively. The outer surfaces of the disks are assumed to be convectively heated by hot liquids having temperatures T_w (lower disk) and T_H (upper disk), with h₁ (lower disk) and h₂ (upper disk) as their heat transfer coefficients. In addition, C_w and C_H denote the NPs concentration at the lower and upper disks, respectively.

FIGURE 1

Based on the aforesaid points, the flow equations are (see Bhattacharyya et al., 2020; Mabood et al., 2021; Yaseen et al., 2022a):

Continuity equation:

Momentum equations:

Energy equation:

Concertation equation:

with Boundary conditions (BCs):

where “ are the velocity components along with directions,” respectively, “C is nanoparticles concentration,” “T is temperature,” “p is pressure,” “k^* is the mean absorption coefficient,” “σ^* is the Stefan-Boltzmann constant,” “F is the Forchheimer coefficient and k_fh is the porous medium permeability,” respectively, “D is diffusion coefficient,” “ is the thermal relaxation time,” and “ is the solutal relaxation time.” Furthermore, “k is thermal conductivity, is dynamic viscosity, is heat capacity, is density, is electrical conductivity.” Moreover, subscript hnf is for hybrid nanofluid, nf is for nanofluid, and f, bf is for fluid. In addition, the last term in the energy Eq. 5 i.e., denotes the significance of “non-uniform heat source/sink” and it is demarcated as follows (Yaseen et al., 2022b):where and the heat source/sink corresponding to space coefficients and corresponding to temperature dependence are signified by constants and , respectively. As a result, the heat source phenomena is characterized by and , whereas the heat sink phenomena is characterized by and .

2.2 Properties of hybrid nanofluid and nanofluid

This study describes the analysis of HNF and NF. The thermal characteristics of HNF and NF are influenced by the base fluid (water) and NPs; as well as by the volume fraction of MoS₂ and Go NPs, as shown in Table 1 (Khan et al., 2020d). The Shape of MoS₂ and Go NPs is taken as spherical. In this paper, the volume fraction of Go NPs and MoS₂ NPs is denoted by and , respectively. Furthermore, s₁ is used for Go NPs, and s₂ is used for MoS₂ NPs. The MoS₂ belongs to the transition metal class and Go belongs to the oxide class. The combination of MoS₂ and Go NPs is taken because metal nanoparticles make a nanolayer atop oxide nanoparticles, and it produces a thermal interfacial layer between weak boundaries of the working fluid and hybrid NPs, resulting in thermal conductivity augmentation.

2.2.1 Thermophysical correlations of Go/water nanofluid

The thermophysical correlations of the NF model used are as follows (Rawat and Kumar, 2020):

2.2.2 Thermophysical correlations of MoS2-Go/water hybrid nanofluid

The thermophysical correlations of the HNF model used are as follows (Khan et al., 2020c):

2.3 Similarity transformations

The transformations listed below are utilized to transform the governing equations (Yaseen et al., 2022a):

Eq. 1 is satisfied, and Eqs 2, 3, 4, 5, 6, 7 reduce to

The BCs are obtained as follows:

In the aforementioned equations, “ represents Reynold number, represents the inertial coefficient, is the Prandtl number, is the porosity parameter, represents the magnetic field parameter, and represents heat source/sink parameters, represents thermal radiation parameter, represents the thermal relaxation parameter, represents the solutal relaxation parameter, and represent shrinking parameters, is the rotation parameter, and are Biot numbers, represents Schmidt number and are constants”.

Eq. 16 is then differentiated w.r.t for the elimination of (pressure variable), and we get:

Eqs 16–21 define the (pressure variable) as follows:

To solve Eqn. (18) for P, use integration from 0 to and obtain the following form:

3 Engineering parameters

3.1 Nusselt numbers

The Nusselt numbers at the lower disk and upper disk are (Mabood et al., 2021) and (Yaseen et al., 2022a):

3.2 Sherwood numbers

The Sherwood numbers at the lower disk and upper disk are defined as (Khan et al., 2020a):

4 Methodology of numerical approach

This section focuses on the methodology used for deducing solutions as well as code validation. The equations are at first modeled as PDEs and later, converted into ODEs via similarity variables. The numerical solution of the Eqs 16, 17, 18, 19, 20 along with BCs (21) is deduced with the “bvp4c function” (a built-in package in MATLAB), the more specific details of the bvp4c function can be referred from the Shampine et al. (2003). The “bvp4c function” uses a finite difference scheme together with a precision of fourth order with the help of the “3-stage Lobatto IIIA formula”. To deduce the solution of the model, the ODEs obtained after similarity transformation are reduced into first-order ODEs by the following substitution:

Utilizing the new variables, the Eqs 16, 17, 18, 19, 20 are reduced to first-order ODEs and the following MATLAB syntax is used:

The above system Eq. 27 is subjected to the following BCs at the lower disk and upper disk :

From new variables in Eq. 26 and BCs Eq. 28, it is seen that the following conditions at the and are missing:

Furthermore, the values of missing conditions are guessed to initiate the process of finding the solution and other parameters present in the Eqs 27, 28 are set to find the desired solution. The process of iteration is repeated and the solution is accepted only when the conditions in Eq. 28 are satisfied. The process of finding the solution is shown via a flow chart in Figure 2. To validate the model and the code used to find the numerical solution, an assessment in Table 2 is outlined with the published computations of (Khan et al., 2018) study as a limiting case to validate the numerical code utilized to solve the current model. The comparative results are quite consistent, ensuring that the current conclusions are valid.

FIGURE 2

5 Results and discussion

This section focuses on the numerical outcomes and graphical and tabular results are interpreted physically. Authors have used the following general values in the numerical computations: , , , , , , , , , , , , , , , and . The analysis is performed for the porosity parameter , magnetic field parameter , nanoparticles volume fraction , thermal radiation parameter , thermal relaxation parameter , heat source/sink parameters , solutal relaxation parameter , and Schmidt number . The influence of the primary relevant factors on the velocity, temperature, and nanoparticle concentration is graphically depicted. Finally, a Table is drawn to show the relationships of various critical factors on the Nusselt and Sherwood number. The results are shown for both HNF (MoS₂-Go/water flow) and NF (Go/water flow) and in the figures, solid lines are drawn for HNF and dot lines are drawn for NF.

5.1 Velocity distribution

Figures 3, 4, 5, 6, 7, 8, 9, 10, 11 show the fluctuation of dimensionless velocity against the similarity variable while altering distinct flow-regulating parameters one at a time while leaving others unchanged. Figures 3, 4, 5 depict the influence of porosity parameter λ on axial velocity distribution , radial velocity distribution , and tangential velocity distribution , respectively. Within the core areas of the boundary layer region (BLR), there is a noticeable decrease in the axial velocity distribution and tangential velocity distribution with growing values of λ. The radial velocity shows transitioning nature in BLR and the transition point lies near . rises near the lower disk and behavior changes near the upper disk w.r.t to growing values of λ. Physically, the porosity parameter is linked to the friction force. Resistance increases with a rise in the porosity parameter, and therefore the velocity falls. Figures 6, 7, 8 elucidate the effectiveness of magnetic parameter M on , and . The increasing magnitude of parameter M is responsible for the decrease in velocity and . The radial velocity shows transitioning nature in BLR w.r.t parameter M. The decrease in the velocity and is related to the generation of Lorentz force. The presence of magnetic force generates the Lorentz force opposite to the flow direction, hence it opposes the motion, and hence velocity decreases.

FIGURE 3

FIGURE 4

FIGURE 5

FIGURE 6

FIGURE 7

FIGURE 8

FIGURE 9

FIGURE 10

FIGURE 11

Figures 9, 10, 11 elucidate the effectiveness of volume fraction on velocity distribution. On increasing the volume fraction of each nanoparticle Go, and MoS₂ in equal quantity , the axial velocity decreases and tangential velocity increases, whereas radial velocity shows transitioning nature in BLR, i.e., it first decreases and then increases. The increasing nanoparticles volume fraction in the base fluid causes the hindrance for the motion in the axial direction, hence the velocity decreases. On contrary, the tangential velocity increases, because the addition of nanoparticles causes the outward movement of fluid. In Figures 3, 4, 5, 6, 7, 8, 9, 10, 11, the authors have presented a comparison in velocity profiles of HNF (MoS₂-Go/water flow) and NF (Go/water flow). HNF is seen to have the higher axial velocity and NF has the higher tangential velocity . Near the lower disk, HNF has a higher radial velocity, but as fluid approaches the upper disk, NF has a higher radial velocity.

5.2 Temperature distribution

Figures 12, 13, 14, 15, 16, 17 show the fluctuation of dimensionless temperature against the similarity variable while altering distinct flow-regulating parameters one at a time while leaving others unchanged. Figures 12, 13, 14 depict the influence of porosity parameter λ, magnetic parameter M, and radiation parameter R_d on the temperature profile. It is observed that with an increment in the aforesaid parameters, the temperature shows transitioning nature in BLR. The temperature falls near the lower disk and rises near the upper disk. The effect of resistance forces due to the porosity parameter and the effect of generated Lorentz force due to the magnetic field is significant near the upper disk, which opposes the motion of the flow and increases friction. As a result, the temperature rises near the upper disk with increasing porosity parameter λ and magnetic parameter M. Figure depicts that augmenting parameter R_d causes the temperature to decrease first, then after a transition point, it increases. This means that near the upper disk, the heat carried by radiation supplies energy to the particles, increasing their kinetic energy and velocity, and causing an increase in the temperature.

FIGURE 12

FIGURE 13

FIGURE 14

FIGURE 15

FIGURE 16

FIGURE 17

Figure 15 depicts that an increase in the thermal relaxation parameter lead the temperature to rise near the lower disk but for the same, the temperature falls near the upper disk. The non-zero thermal relaxation parameter denotes time lag during heat transfer. For the zero value of , the heat transfer is governed by Fourier’s law. The results depict that near the lower disk the more time lag in heat transfer implies higher temperature but the same condition near the upper disk relates to the lower temperature. Figure 16 elucidates the effectiveness of volume fraction on temperature . On increasing the volume fraction of each nanoparticle Go , and MoS₂ in equal quantity , the temperature falls near the lower disk and surges near the upper one. The rise in the temperature near the upper disk can be attributed to the enhancement in the thermal conduction of the HNF or NF due to an increase in the nanoparticles in the base fluid. Figure 17 depicts that for the case of heat source , the temperature shows transitioning nature in BLR. Temperature first increases then decreases for increasing and . Whereas, for the case of the heat sink , the temperature increases for increasing magnitude of heat sink parameters. In Figures 12, 13, 14, 15, 16, 17, the authors have presented a comparison of temperature profiles of HNF (MoS₂-Go/water flow) and NF (Go/water flow). Near the lower disk, HNF has a higher thermal profile , but as fluid approaches the upper disk, NF has a higher thermal profile .

5.3 Nanoparticle concentration distribution

Figures 18, 19, 20 show the fluctuation of dimensionless concentration of nanoparticles against the similarity variable while altering distinct flow regulating parameters one at a time while leaving others unchanged. Figures 18, 19 depict the influence of the solutal relaxation parameter and Schmidt number Sc on the concentration . Figure 18 depicts that increasing values of the solutal relaxation parameter leads concentration to fall. The non-zero values of the solutal relaxation parameter denote the presence of time lag during mass transfer. The increasing time lag during mass transfer has an adverse effect on concentration. The increasing time lag during mass transfer leads to a fall in the concentration profile. Figure 19 depicts that increasing values of Sc cause concentration to fall. The concentration decreases as Sc increases. The Schmidt number is an important parameter to consider as it physically connects the depths of the hydrodynamic layer and the mass transfer boundary layer. The momentum diffusion of particles in fluid increases as the strength of the Schmidt number increases in fluid flow. As a result, concentration falls. Figure 20 depicts that on increasing the volume fraction of each nanoparticle Go , and MoS₂ in equal quantity , the concentration falls. The reason for this behavior can be attributed to the fact the mixing of more nanoparticles causes the HNF or NF to be more viscous and flow is resisted, hence the concentration of NPs falls. In Figures 18, 19, 20, the authors have presented a comparison in concentration profiles of HNF (MoS₂-Go/water flow) and NF (Go/water flow). It is seen that NF has a higher concentration profile as compared to HNF.

FIGURE 18

FIGURE 19

FIGURE 20

5.4 Streamlines and engineering parameters

Figure 21 visualizes the streamlines pattern of HNF (MoS₂-Go/water flow) and NF (Go/water flow). The streamlines represent the movement of the suspended particles in the stream and which are driven by it. The direction of velocity at each point in the streamline is given by the tangent at that location. The denser streamlines indicate that the fluid velocity is greater than when it is open out. The HNF is seen to have slightly less dense streamlines in the present case which means more nanoparticles cause fluid to be more viscous and flow is resisted.

FIGURE 21

Table 3 elucidates the effectiveness of pertinent parameters on the Nusselt and Sherwood number. Nusselt number corresponds to the heat transmission rate (HTR) and the Sherwood number corresponds to the mass transmission rate (MTR) in the BLR at both disks. At the higher magnitude of the heat sink parameter, the HTR rises for HNF and NF at the lower disk but contrary behavior is seen at the upper disk. Whereas, the HTR decreases for HNF and NF at both the upper and lower disk. At the higher values of radiation parameter R_d, HTR increases at the lower disk and decreases at the upper disk. An increase in the volume fraction of each nanoparticle Go, and MoS₂ in equal quantity enhances the thermal conductivity of the THNF due to an increase in the nanoparticles in the base fluid. The rising thermal conductivity enhances the MTR and HTR at the upper disk. It is also observed that at the higher values of the thermal relaxation parameter , HTR reduces at both disks. The non-zero values of the thermal relaxation parameter denote the time lag during heat transfer. Thus on augmenting , HTR at both disks reduces. Moreover, on increasing the magnitude of the Schmidt number Sc and solutal relaxation parameter , the MTR reduces at upper disks, whereas, on increasing the magnitude of the same parameters, the MTR rises at the lower disk.

6 Conclusion

6.1 Future scope of research

The present study discuss the flow of nanofluid and hybrid nanofluid without the nanoparticle aggregation effect. The study can de extended with proper modeling of fluid flow by utilization of validated thermophysical correlations for nanoparticle aggregation effect.

References

• 1

  AliN.AsgharZ.SajidM.Anwar BégO. (2019b). Biological interactions between Carreau fluid and microswimmers in a complex wavy canal with MHD effects. Journal of the Brazilian Society of Mechanical Sciences and Engineering. 41, 446–513. 10.1007/S40430-019-1953-Y

  • CrossRef
  • Google Scholar

• 2

  AliN.AsgharZ.SajidM.AbbasF. (2019a). A hybrid numerical study of bacteria gliding on a shear rate-dependent slime. Physica A: Statistical Mechanics and its Applications. 535, 122435. 10.1016/J.PHYSA.2019.122435

  • CrossRef
  • Google Scholar

• 3

  AlrabaiahH.BilalM.KhanM. A.MuhammadT.LegasE. Y. (2022). Parametric estimation of gyrotactic microorganism hybrid nanofluid flow between the conical gap of spinning disk-cone apparatus. Scientific Reports. 12, 59–14. 10.1038/s41598-021-03077-2

  • CrossRef
  • Google Scholar

• 4

  AsgharZ.AliN.JavidK.WaqasM.DogonchiA. S.KhanW. A. (2020a). Bio-inspired propulsion of micro-swimmers within a passive cervix filled with couple stress mucus. The Computer Methods and Programs in Biomedicine. 189, 105313. 10.1016/j.cmpb.2020.105313

  • CrossRef
  • Google Scholar

• 5

  AsgharZ.AliN.WaqasM.JavedM. A. (2020b). An implicit finite difference analysis of magnetic swimmers propelling through non-Newtonian liquid in a complex wavy channel. Computers and Mathematics with Applications. 79, 2189–2202. 10.1016/J.CAMWA.2019.10.025

  • CrossRef
  • Google Scholar

• 6

  AsgharZ.AliN.WaqasM.NazeerM.KhanW. A. (2020c). Locomotion of an efficient biomechanical sperm through viscoelastic medium. Biomechanics and Modeling in Mechanobiology. 19, 2271–2284. 10.1007/S10237-020-01338-Z

  • CrossRef
  • Google Scholar

• 7

  AsgharZ.Saeed KhanM. W.GondalM. A.GhaffariA. (2022a). Channel flow of non-Newtonian fluid due to peristalsis under external electric and magnetic field. Proceedings of the Institution of Mechanical Engineers, Part E. 236, 2670–2678. 10.1177/09544089221097693

  • CrossRef
  • Google Scholar

• 8

  AsgharZ.ShahR. A.AliN. (2022b). A computational approach to model gliding motion of an organism on a sticky slime layer over a solid substrate. Biomechanics and Modeling in Mechanobiology. 21, 1441–1455. 10.1007/S10237-022-01600-6

  • CrossRef
  • Google Scholar

• 9

  AsgharZ.ShahR. A.PashaA. A.RahmanM. M.KhanM. W. S. (2022c). Controlling kinetics of self-propelled rod-like swimmers near multi sinusoidal substrate. Computers in Biology and Medicine. 151, 106250. 10.1016/J.COMPBIOMED.2022.106250

  • CrossRef
  • Google Scholar

• 10

  AsgharZ.ShatanawiW.ShahR. A.GondalM. A. (2023). Impact of viscoelastic ooze slime on complex wavy gliders near a solid boundary. Chinese Journal of Physics. 81, 26–36. 10.1016/J.CJPH.2022.10.013

  • CrossRef
  • Google Scholar

• 11

  AsgharZ.WaqasM.Asif GondalM.Azeem KhanW. (2022d). Electro-osmotically driven generalized Newtonian blood flow in a divergent micro-channel. Alexandria Engineering Journal. 61, 4519–4528. 10.1016/J.AEJ.2021.10.012

  • CrossRef
  • Google Scholar

• 12

  BaiY.TangQ.ZhangY. (2022). Unsteady MHD oblique stagnation slip flow of Oldroyd-B nanofluids by coupling Cattaneo-Christov double diffusion and Buongiorno model. Chinese Journal of Physics. 79, 451–470. 10.1016/J.CJPH.2022.09.013

  • CrossRef
  • Google Scholar

• 13

  BhattacharyyaA.SethG. S.KumarR.ChamkhaA. J. (2020). Simulation of Cattaneo–Christov heat flux on the flow of single and multi-walled carbon nanotubes between two stretchable coaxial rotating disks. Journal of Thermal Analysis and Calorimetry. 139, 1655–1670. 10.1007/s10973-019-08644-4

  • CrossRef
  • Google Scholar

• 14

  Bilal AshrafM.Roohani GhehsarehH. (2021). Model of Casson fluid with Cattaneo–Chirstov heat flux and Hall effect. Indian Journal of Physics. 95, 1469–1477. 10.1007/s12648-020-01819-y

  • CrossRef
  • Google Scholar

• 15

  CattaneoC. (1948). Sulla Conduzione del Calore. Atti del Seminario Matematico e Fisico dell'Universita di Modena3, 83–101. 10.1007/978-3-642-11051-1_5

  • CrossRef
  • Google Scholar

• 16

  ChoiS. U. S. (1995). Enhancing thermal conductivity of fluids with nanoparticles. FED-American Society of Mechanical Engineers. Fluids Engineering Division Newsletter. 231, 99–105.

  • Google Scholar

• 17

  ChristovC. I. (2009). On frame indifferent formulation of the Maxwell-Cattaneo model of finite-speed heat conduction. Mechanics Research Communications. 36, 481–486. 10.1016/j.mechrescom.2008.11.003

  • CrossRef
  • Google Scholar

• 18

  DeviS. S. U. P. A. S. U.DeviS. S. U. P. A. S. U. (2016). Numerical investigation of three-dimensional hybrid Cu-Al2O3/water nanofluid flow over a stretching sheet with effecting Lorentz force subject to Newtonian heating. Canadian Journal of Physics. 94, 490–496. 10.1139/CJP-2015-0799

  • CrossRef
  • Google Scholar

• 19

  FourierJ. B. (1822). Théorie analytique de la chaleur. Wuhan, China; Scientific Research Publishing.

  • Google Scholar

• 20

  GariaR.RawatS. K.KumarM.YaseenM. (2021). Hybrid nanofluid flow over two different geometries with cattaneo–christov heat flux model and heat generation: A model with correlation coefficient and probable error. Chinese Journal of Physics. 74, 421–439. 10.1016/j.cjph.2021.10.030

  • CrossRef
  • Google Scholar

• 21

  IrfanM.KhanM.KhanW. A. (2018). Interaction between chemical species and generalized Fourier’s law on 3D flow of Carreau fluid with variable thermal conductivity and heat sink/source: A numerical approach. Results in Physics. 10, 107–117. 10.1016/j.rinp.2018.04.036

  • CrossRef
  • Google Scholar

• 22

  JavidK.AliN.AsgharZ. (2019). Rheological and magnetic effects on a fluid flow in a curved channel with different peristaltic wave profiles. Journal of the Brazilian Society of Mechanical Sciences and Engineering. 41, 483–514. 10.1007/S40430-019-1993-3

  • CrossRef
  • Google Scholar

• 23

  KhanM. I.QayyumS.ChuY. M.KhanN. B.KadryS.MackolilJ.et al (2021). Transportation of Marangoni convection and irregular heat source in entropy optimized dissipative flow. International Communications in Heat and Mass Transfer. 120, 105031–105110. 10.1016/j.icheatmasstransfer.2020.105031

  • CrossRef
  • Google Scholar

• 24

  KhanM. I.QayyumS.HayatT.AlsaediA. (2018). Entropy generation minimization and statistical declaration with probable error for skin friction coefficient and Nusselt number. Chinese Journal of Physics. 56, 1525–1546. 10.1016/j.cjph.2018.06.023

  • CrossRef
  • Google Scholar

• 25

  KhanM. S.MeiS.ShabnamAli ShahN.ChungJ. D.KhanA.et al (2022). Steady squeezing flow of magnetohydrodynamics hybrid nanofluid flow comprising carbon nanotube-ferrous oxide/water with suction/injection effect. Nanomaterials12, 660. 10.3390/nano12040660

  • CrossRef
  • Google Scholar

• 26

  KhanN. S.ShahQ.BhaumikA.KumamP.ThounthongP.AmiriI. (2020a). Entropy generation in bioconvection nanofluid flow between two stretchable rotating disks. Scientific Reports. 10, 4448–4526. 10.1038/s41598-020-61172-2

  • CrossRef
  • Google Scholar

• 27

  KhanU.AhmadS.HayyatA.KhanI.Sooppy NisarK.BaleanuD. (2020b). On the cattaneo-christov heat flux model and OHAM analysis for three different types of nanofluids. Applied Sciences. 10, 886. 10.3390/app10030886

  • CrossRef
  • Google Scholar

• 28

  KhanU.ShafiqA.ZaibA.BaleanuD. (2020c). Hybrid nanofluid on mixed convective radiative flow from an irregular variably thick moving surface with convex and concave effects. Case Studies in Thermal Engineering. 21, 100660. 10.1016/j.csite.2020.100660

  • CrossRef
  • Google Scholar

• 29

  KhanU.ZaibA.SheikholeslamiM.WakifA.BaleanuD. (2020d). Mixed convective radiative flow through a slender revolution bodies containing molybdenum-disulfide graphene oxide along with generalized hybrid nanoparticles in porous media. Crystal (Basel)10, 771–818. 10.3390/cryst10090771

  • CrossRef
  • Google Scholar

• 30

  KumarB.SethG. S.SinghM. K.ChamkhaA. J. (2021). Carbon nanotubes (CNTs)-based flow between two spinning discs with porous medium, Cattaneo–Christov (non-Fourier) model and convective thermal condition. Journal of Thermal Analysis and Calorimetry. 146, 241–252. 10.1007/s10973-020-09952-w

  • CrossRef
  • Google Scholar

• 31

  MaboodF.BerrehalH.YusufT. A.KhanW. A. (2021). Carbon nanotubes-water between stretchable rotating disks with convective boundary conditions: Darcy-forchheimer scheme. International Journal of Ambient Energy. 43, 3981–3994. 10.1080/01430750.2021.1874527

  • CrossRef
  • Google Scholar

• 32

  MadhukeshJ. K.Naveen KumarR.Punith GowdaR. J.PrasannakumaraB. C.RameshG. K.Ijaz KhanM.et al (2021). Numerical simulation of aa7072-aa7075/water-based hybrid nanofluid flow over a curved stretching sheet with Newtonian heating: A non-fourier heat flux model approach. Journal of Molecular Liquids. 335, 116103. 10.1016/j.molliq.2021.116103

  • CrossRef
  • Google Scholar

• 33

  MasoodS.FarooqM.AnjumA. (2021). Influence of heat generation/absorption and stagnation point on polystyrene-TiO 2/H 2 O hybrid nanofluid flow. Scientific Reports. 11, 22381. 10.1038/s41598-021-01747-9

  • CrossRef
  • Google Scholar

• 34

  NazR.RazaR.ZafranY.JavedM. (2022). Scrutiny of Carreau nanoliquid over rotating disk with radiation impacts and the Cattaneo–Christov flux theory. Waves Random Complex Media, 1–16. 10.1080/17455030.2022.2128231

  • CrossRef
  • Google Scholar

• 35

  QureshiM. Z. A.BilalS.MalikM. Y.RazaQ.SherifE. S. M.LiY. M. (2021). Dispersion of metallic/ceramic matrix nanocomposite material through porous surfaces in magnetized hybrid nanofluids flow with shape and size effects. Scientific Reports. 11, 12271–12319. 10.1038/s41598-021-91152-z

  • CrossRef
  • Google Scholar

• 36

  RahmanM.SharifF.TurkyilmazogluM.SiddiquiM. S. (2022). Unsteady three-dimensional magnetohydrodynamics flow of nanofluids over a decelerated rotating disk with uniform suction. Pramana96, 170–179. 10.1007/S12043-022-02404-0

  • CrossRef
  • Google Scholar

• 37

  Ranga BabuJ. A.KumarK. K.Srinivasa RaoS. (2017). State-of-art review on hybrid nanofluids. Pergamon.

  • Google Scholar

• 38

  RashidU.LiangH.AhmadH.AbbasM.IqbalA.HamedY. S. (2021). Study of (Ag and TiO2)/water nanoparticles shape effect on heat transfer and hybrid nanofluid flow toward stretching shrinking horizontal cylinder. Results in Physics. 21, 103812. 10.1016/j.rinp.2020.103812

  • CrossRef
  • Google Scholar

• 39

  RawatS. K.KumarM. (2020). Cattaneo–christov heat flux model in flow of copper water nanofluid through a stretching/shrinking sheet on stagnation point in presence of heat generation/absorption and activation energy. International Journal of Applied and Computational Mathematics. 6, 112. 10.1007/s40819-020-00865-8

  • CrossRef
  • Google Scholar

• 40

  ShahR. A.AsgharZ.AliN. (2022). Mathematical modeling related to bacterial gliding mechanism at low Reynolds number with Ellis Slime. The European Physical Journal Plus. 137, 600–612. 10.1140/EPJP/S13360-022-02796-3

  • CrossRef
  • Google Scholar

• 41

  ShampineL. F.GladwellI.ThompsonS. (2003). Solving ODEs with MATLAB. Cambridge, United States: Cambridge University Press.

  • Google Scholar

• 42

  TurkyilmazogluM. (2022a). Flow and heat over a rotating disk subject to a uniform horizontal magnetic field. Zeitschrift fur Naturforschung - Section A Journal of Physical Sciences. 77, 329–337. 10.1515/zna-2021-0350

  • CrossRef
  • Google Scholar

• 43

  TurkyilmazogluM. (2021). Heat transfer enhancement feature of the non-fourier cattaneo-christov heat flux model. Journal of Heat Transfer. 143, 094501. 10.1115/1.4051671

  • CrossRef
  • Google Scholar

• 44

  TurkyilmazogluM. (2022b). Multiple exact solutions of free convection flows in saturated porous media with variable heat flux. Journal of Porous Media. 25, 53–63. 10.1615/JPORMEDIA.2022041870

  • CrossRef
  • Google Scholar

• 45

  WaqasH.FarooqU.NaseemR.HussainS.AlghamdiM. (2021). Impact of MHD radiative flow of hybrid nanofluid over a rotating disk. Case Studies in Thermal Engineering. 26, 101015. 10.1016/j.csite.2021.101015

  • CrossRef
  • Google Scholar

• 46

  WuA.AbbasS. Z.AsgharZ.SunH.WaqasM.KhanW. A. (2020). A shear-rate-dependent flow generated via magnetically controlled metachronal motion of artificial cilia. Biomechanics and Modeling in Mechanobiology. 19, 1713–1724. 10.1007/S10237-020-01301-Y

  • CrossRef
  • Google Scholar

• 47

  YaseenM.RawatS. K.KumarM. (2022a). Cattaneo–Christov heat flux model in Darcy–Forchheimer radiative flow of MoS2–SiO2/kerosene oil between two parallel rotating disks. Journal of Thermal Analysis and Calorimetry. 147, 10865–10887. 10.1007/s10973-022-11248-0

  • CrossRef
  • Google Scholar

• 48

  YaseenM.RawatS. K.KumarM. (2022b). Hybrid nanofluid (MoS2–SiO2/water) flow with viscous dissipation and Ohmic heating on an irregular variably thick convex/concave-shaped sheet in a porous medium. Heat transfer. 51, 789–817. 10.1002/htj.22330

  • CrossRef
  • Google Scholar

• 49

  YaseenM.RawatS. K.ShafiqA.KumarM.NonlaoponK. (2022c). Analysis of heat transfer of mono and hybrid nanofluid flow between two parallel plates in a Darcy porous medium with thermal radiation and heat generation/absorption. Symmetry14, 1943. 10.3390/sym14091943

  • CrossRef
  • Google Scholar