Home Physical Sciences Investigation of convective heat transport in a Carreau hybrid nanofluid between two stretchable rotatory disks
Article Open Access

Investigation of convective heat transport in a Carreau hybrid nanofluid between two stretchable rotatory disks

  • MD. Shamshuddin , Zahir Shah EMAIL logo , Usman , Anwar Saeed , Mansoor H. Alshehri , Narcisa Vrinceanu EMAIL logo and Elisabeta Antonescu
Published/Copyright: November 4, 2024

Abstract

Hybrid nanofluids (HNFs) have outstanding energy transfer capabilities that are comparable to mono-nanofluids. Materials had appliances in obvious fields such as heat generation, micropower generation, and solar collectors. The objective of this study is to investigate the new aspects of convective heat transfer in an electrically conducting Carreau HNF situated between two parallel discs. In addition to the presumed stretchability and rotation of the discs, physical phenomena like nonlinear radiation, viscous dissipation, Joule dissipation, and heat generation and absorption are considered. The Cu and TiO2 nanoparticles dispersed in engine oil to understand the intricate phenomenon of hybridization. The Tiwari and Das nanofluid model is employed to model the governing partial differential equations (PDEs) and then simplified using boundary layer approximation. The suitable transformations of similarity variables are defined and implemented to change the set of formulated PDEs into ordinary differential equations. The reduced system is solved semi-analytically by the homotopy analysis method. The influences of involving physical parameters on the velocity and temperature are plotted with the help of graphical figures. This study brings forth a significant contribution by uncovering novel flow features that have previously remained unexplored. By addressing a well-defined problem, our research provides valuable insights into the enhancement of thermal transport, with direct implications for diverse engineering devices such as solar collectors, heat exchangers, and microelectronics.

1 Introduction

The study of fluid flow and heat transmission in stretchy rotational discs has a broad spectrum of uses in industrial, mechanical, and thermal engineering. These applications include jet motors, spin coating, food processing, electric power generation systems, manufacturing, heat exchangers, gas turbines, medical equipment, and many others. As a result of its very valuable practical applications, researchers all over the world are actively exploring the movement of fluids and thermal exchange studies of rotating discs. The study of non-Newtonian fluid flow caused by rotating discs is an exciting and dynamic field that has piqued the interest of scientists from numerous disciplines. Non-Newtonian fluids find application across diverse industries, including paper manufacturing, petroleum extraction, and polymer sheet processing, among others. Because of the variety of non-Newtonian fluids, their behavior cannot be characterized by a single constitutive equation. Kármán [1] conducted substantial research on the fluid dynamics issue concerning rotating discs and provided a mathematical model in this context. He devised an interesting similarity transformation that converts the controlling partial differential equations (PDEs) to ordinary differential equations (ODEs). These ODEs are then resolved numerically using the momentum integral approach, allowing for more efficient system analysis. Von Karman’s work has established a solid foundation for understanding and investigating the intricate fluid dynamics and heat transport in spinning disc systems. Fang and Zhang [2] accomplished a great feat by getting the correct solution to the Navier–Stokes equations for a fluid flowing between two infinitely flexible discs undergoing stretching velocity acceleration. Imtiaz et al. [3] investigated the effect of single and multiwall carbon nanotubes on the flow and heat transmission between two spinning and stretching discs in great detail. Ahmed et al. [4] proposed a novel mathematical model to investigate the behavior of an electrically conductive Maxwell fluid. This fluid flows between two axially aligned rotating and stretching discs, keeping their spacing constant. Usman et al. [5] researched on increasing heat transport within a blood nanofluid. They used the power-law model to investigate the effect of a heat source/sink in the presence of two revolving stretchable discs to promote the heat transfer process. Usman et al. [6] explored the three-dimensional heat transfer and liquid flow processes that occur within a Carreau liquid contained between rotating stretchy discs. Important advances in the study of nanofluid flow over a variety of geometries were made by Das and Ali [7,8]. Yashkun et al. [9] scrutinized the heat transmission characteristics of magnetohydrodynamic (MHD) hybrid nanofluids (HNFs) on shrinking and stretching surfaces, taking into account the effects of suction and thermal radiation. Lund et al. [10] studied the effect of Joule heating and viscous dissipation on unsteady MHD 2D flow over a constantly parallel regular shrinking sheet.

In the past few years, scientists have been actively looking for new fluids to replace nanofluids and attain highly efficient heat transfer capacities. This pursuit has given rise to the phrase “hybrid nanofluids (HNFs).” By integrating two distinct nanoparticle variants into a foundational fluid, these novel fluids go beyond traditional nanofluids. When compared to ordinary nanofluids, the addition of numerous nanoparticles in HNFs results in superior thermal augmentation. Waini et al. [11], for example, investigated the impact of Cu–Al2O3–water HNF flow and heat transmission along a leaky moving sheet. In terms of increasing heat transfer qualities, their findings indisputably proved that HNFs outperform nanofluids having only a single type of nanoparticle. This growing interest in HNFs highlights its potential as a cost-effective solution for a wide array of applications involving heat transfer. Zainal et al. [12] explored mixed convection flow in an HNF across a vertical flat plate under the effect of MHD and convective boundary conditions. Their research found dual solutions in both the aiding and opposing flow areas. Waini et al. [13] investigated the consistent movement of HNF with mixed convection over a porous vertical plane. Within a given range of buoyancy factors, they discovered numerous dual solutions. Furthermore, the introduction of hybrid nanoparticles was discovered to slow down boundary layer separation. Rostami et al. [14] investigated sustained laminar MHD-mixed convection stagnation point flow in an HNF across a vertical permeable flat plate. They discovered two distinct solutions within a specified range of the buoyancy parameter and performed stability tests on these solutions. Shamshuddin et al. [15] examined the flow of HNF across a spinning surface and the effects of entropy formation. Farooq et al. [16] investigated the electromagnetic flow of HNF between two discs. The impacts of entropy formation on the flow of HNF across discs were investigated by Agrawal and Kaswan [17]. Waini et al. [18] used a Cu–Al2O3–water HNF to evaluate the unsteady flow and heat transport toward a shrinking–stretching surface. They discovered that as the suction parameter was increased, the velocity of the HNF increased. Furthermore, they revealed that the employment of the HNF, together with the stimulus of the suction force, resulted in an additional temperature reduction. Khashi’ie et al. [19] explored the outcome of convective circumstances and velocity slip on the 3D flow of HNF over a surface undergoing stretching and shrinking numerically. Increasing the Biot and slip parameters boosted heat transfer rates, whereas growing the Cu solid volume fraction had the reverse impact. Anuar et al. [20] investigated the stimulus of boundary layer flow and heat transmission over an inclined leaky surface that undergoes both stretching and shrinking, utilizing an HNF composed of Ag–MgO–water. They discovered that increasing the angle of inclination and suction parameters increased the local Nusselt number. Furthermore, increasing the volume percentage of Ag nanoparticles in the MgO–water nanofluid resulted in a decrease in the local Nusselt number.

Non-Newtonian fluids have gained a lot of interest in the industrial world because of their diverse applications. Much experimental and theoretical research has been conducted to better understand the transport behavior of these fluids, which are important in a variety of genetic and industrial processes, materials, and mechanical engineering. Non-Newtonian liquids can be found in a wide range of materials, including biomedical fluxes, splashes, bio-fluids in genetic material, and polymers. Despite extensive study in this field, many researchers continue to investigate non-Newtonian fluids from diverse angles, resulting in countless publications [21,22,23]. The well-designed Carreau fluid model [24] stands out among the constitutive models developed to represent the features of shear thinning/thickening non-Newtonian liquids. Sulochana et al. [25], for example, used the Runge–Kutta method in conjunction with a shooting operation to quantitatively analyze the effect of a stagnation point on the Carreau nanofluid. The magnetic parameter had a considerable influence on the flow field according to their findings. Khan et al. [26] recently proposed a new model for the three-dimensional flow of Carreau fluid that incorporates nonlinear thermal radiation. Irfan et al. [27] explored the 3D flow of Carreau fluid in more detail, taking into account numerous additional features. The thermal Biot number, Brownian motion, thermophoresis, and thermal radiation characteristics all had substantial effects on the temperature field, according to their research. Rooman et al. [28] investigated the heat transfer characteristics of Carreau nanofluid in the context of the renal tubule. There are many emerging streams that are associated with the different materials used in engineering science technology, which includes studies such as thermoelastic interaction in functional graded material [29], photo-thermal interactions in a semi-conductor material [30], and prothermal elastic materials [31].

This current analysis effort intends to analyze the importance of nonlinear radiation, heat, and Joule dissipation effects in the act of the phenomenon of surface heating by parallel rotatory disk in the presence of an HNF. By investigating the behavior of Cu and TiO2 nanoparticles dispersed in water, we aim to understand the intricate phenomenon of hybridization. This current novelty is to put effort to analyzing the importance of nonlinear radiation, heat, and Joule dissipation effects in the act phenomenon of surface heating by parallel rotatory disk in the presence of HNF.

2 Formulation of the problem

The following assumptions are made while formulating the problem:

  • Consider a 3D flow of the Carreau HNF flow between two parallel, rotating, and stretching discs, where the bottom disc is set at z = 0 and the upper disc is set at z = h .

  • The presumed axial-symmetric flow, which distributes radially on the discs and is electrically conducting, strikes them at a 90° angle.

  • The bottom and upper discs are rotating at a rate of r s 1 , r s 2 and stretching at a rate of r Ω 1 , r Ω 2 , respectively.

  • The constant magnetic field while preserving B 0 ’s strength is applied perpendicular to the rotating discs.

  • The axial-symmetric flow assumes that the pressure is constant and eliminates any fluctuation with respect to θ .

  • Under nonlinear radiation and heat generation/absorption concerning convective boundary conditions, the phenomena of heat transfer are examined.

  • In addition, T 1 and T 2 stand for the temperatures of the lower and top discs, respectively (Figure 1).

Figure 1 
               Flow configuration of magnetized Carreau HNF.
Figure 1

Flow configuration of magnetized Carreau HNF.

The leading PDEs of motion and energy for the time-independent 3D flow of a hybrid Carreau nanofluid between two stretchable spinning discs following Usman et al. [5], while taking into account viscous dissipation and joule heating due to convective heat transfer are expressed as follows:

(1) u r + u r + w z = 0 ,

(2) u u r v 2 r + w u z = υ hnf 2 u z 2 1 + Γ 2 u z 2 n 1 2 + υ hnf ( n 1 ) Γ 2 u z 2 2 u z 2 1 + Γ 2 u z 2 n 3 2 σ hnf B 0 2 ρ hnf u ,

(3) u v r + u v r + w v z = υ hnf 2 v z 2 1 + Γ 2 v z 2 n 1 2 + υ hnf ( n 1 ) Γ 2 u z 2 2 u z 2 1 + Γ 2 u z 2 n 3 2 σ hnf B 0 2 ρ hnf v ,

(4) u T r + w T z = α hnf 2 T z 2 + 16 σ 3 k ( ρ c p ) hnf z T 3 T z + μ hnf ( ρ c p ) hnf u z 2 + v z 2 + σ hnf B 0 2 ( ρ c p ) hnf ( u 2 + v 2 ) + Q 0 ( ρ c p ) hnf ( T T 2 ) .

The boundary conditions are as given as follows:

(5) u = r s 1 , v = r Ω 1 , w = 0 , k f T z = h 1 ( T 1 T ) , at z = 0 , u = r s 2 , v = r Ω 2 , w = 0 , k f T z = h 2 ( T T 2 ) , at z = h .

2.1 Hybrid nanomaterial properties

The thermophysical characteristics and relation of nano and hybrid nanomaterials are presented in Tables 1 and 2, respectively.

Table 1

Experimental values of engine oil, Cu, and TiO2 [32]

ρ ( kg / m 3 ) c p ( J / kg K ) k ( W / m K ) σ ( Ω m ) 1
Engine oil (EO) 884 1910 0.1440 5.5 × 10−6
Copper (Cu) 8,933 385.0 401.0 5.96 × 10−8
Titanium dioxide (TiO2) 4,250 686.20 8.9538 6.27 × 10−5
Table 2

Nano and hybrid nanofluids’ thermo-physical interactions [32]

Properties Nano and hybrid
Viscosity μ hnf μ bf = 1 ( 1 ϕ Cu ϕ TiO 2 ) 2.5
Density ρ hnf ρ bf = ϕ Cu ρ Cu ρ bf + ϕ TiO 2 ρ TiO 2 ρ bf + ( 1 ϕ Cu ϕ TiO 2 )
Thermal capacity ( ρ c p ) hnf ( ρ c p ) bf = ϕ Cu ( ρ c p ) Cu ( ρ c p ) bf + ϕ TiO 2 ( ρ c p ) Cu ( ρ c p ) bf + ( 1 ϕ Cu ϕ TiO 2 )
Thermal conductivity k hnf k bf = ( ϕ Cu k Cu + ϕ TiO 2 k TiO 2 ) ϕ Cu + ϕ TiO 2 + 2 k bf + 2 ( ϕ Cu k Cu + ϕ TiO 2 k TiO 2 ) 2 ( ϕ Cu + ϕ TiO 2 ) k bf ( ϕ Cu k Cu + ϕ TiO 2 k TiO 2 ) ϕ Cu + ϕ TiO 2 + 2 k bf ( ϕ Cu k Cu + ϕ TiO 2 k TiO 2 ) + ( ϕ Cu + ϕ TiO 2 ) k bf
Electrical conductivity σ hnf σ bf = ( ϕ Cu σ Cu + ϕ TiO 2 σ TiO 2 ) ϕ Cu + ϕ TiO 2 + 2 σ bf + 2 ( ϕ Cu σ Cu + ϕ TiO 2 σ TiO 2 ) 2 ( ϕ Cu + ϕ TiO 2 ) σ bf ( ϕ Cu k Cu + ϕ TiO 2 k TiO 2 ) ϕ Cu + ϕ TiO 2 + 2 σ b f ( ϕ Cu σ Cu + ϕ TiO 2 σ TiO 2 ) + ( ϕ Cu + ϕ TiO 2 ) σ bf

2.1.1 Similarity transformation

The similarity parameters for the current issue can be characterized as follows:

(6) η = z Ω 1 υ f , u = r Ω 1 F ( η ) , v = r Ω 1 G ( η ) , w = Ω 1 υ f H ( η ) , θ ( η ) = T T 2 T 1 T 2 .

By introducing (6) into (1)–(5), we achieved the ODEs:

(7) H + 2 F = 0 ,

(8) F 2 G 2 + F H = A 2 A 1 F { 1 + We 2 F 2 } n 1 2 + A 2 A 1 ( n 1 ) × We 2 F 2 F { 1 + We 2 F 2 } n 3 2 A 3 A 1 M F ,

(9) 2 F G + H G = A 2 A 1 G { 1 + We 2 G 2 } n 1 2 + A 2 A 1 ( n 1 ) × We 2 G 2 G { 1 + We 2 G 2 } n 3 2 A 3 A 1 M G ,

(10) Pr H θ = A 5 A 4 θ + 4 Rd 3 A 4 { ( 1 + ( θ w 1 ) θ ) 3 θ } + A 2 A 4 PrEc ( F 2 + G 2 ) + A 3 A 4 PrEc M ( F 2 + G 2 ) + Pr Q A 4 θ .

The transformed boundary constraints are presented by

(11) F ( 0 ) = A , G ( 0 ) = 1 , H ( 0 ) = 0 , θ ( 0 ) = C ( 1 θ ( 0 ) ) , F ( 1 ) = B , G ( 1 ) = Ω , H ( 1 ) = 0 , θ ( 1 ) = D θ ( 1 ) .

In this case, the dimensionless emergent parameters are specified as follows:

(12) A 1 = ρ hnf ρ f , A 2 = μ hnf μ f , A 3 = σ hnf σ f , A 4 = ( ρ c p ) hnf ( ρ c p ) f , A 5 = k hnf k f , M = B 0 2 σ f Ω 1 ρ f , A = s 1 Ω 1 , B = s 2 Ω 1 , Ω = Ω 2 Ω 1 , We = Ω 1 3 2 r Γ υ f , Ec = r 2 Ω 1 2 ( ρ c p ) f ( T 1 T 2 ) , Pr = υ f α f , Rd = 4 σ T 3 k f k , θ w = T 1 T 2 , Q = Q 0 ( ρ c p ) f Ω 1 , C = h 1 k f υ f Ω 1 , D = h 2 k f υ f Ω 1 ,

where A 1 , A 2 , A 3 , A 4 , A 5 are constants. M , A , B , Ω , We , Ec , Pr , Rd , θ w , Q , C and D are the magnetic field parameter, stretching parameters, Weissenberg number, rotation parameter, Eckert number, Prandtl number, radiation parameter, surface heating parameter, heat generation/absorption parameter, and lower and upper disk parameters.

2.2 Physically interesting quantities

2.2.1 Surface drag forces

The surface drag forces or skin friction coefficients at the heated surface are defined as follows in radial and azimuthal directions:

(13) C Fr = τ r z ρ f ( cr ) 2 , C G θ = τ G θ ρ f ( cr ) 2 ,

where τ w r and τ w θ are radial and azimuthal direction shear stresses, respectively, which are specified by:

(14) τ r z = μ hnf u z 1 + Γ 2 u z 2 n 1 2 z = 0 , τ θ z = μ hnf v z 1 + Γ 2 v z 2 n 1 2 z = 0 .

By applying (6), the dimensionless skin friction coefficient is reduced as follows:

(15) Re 1 2 C Fr = A 2 F ( 0 ) [ 1 + We 2 ( F ( 0 ) ) 2 ] n 1 2 , Re 1 2 C G θ = A 2 G ( 0 ) [ 1 + We 2 ( G ( 0 ) ) 2 ] n 1 2 .

2.2.2 Heat transfer rate

The heat transfer rate, or local Nusselt number, or temperature gradient is denoted as follows:

(16) Nu r = r q w k f ( T 1 T 2 ) .

Here, q w is the heat flux and is described as follows:

(17) q w = k hnf T z + ( q r ) w z = 0 .

The nondimensional will lead to

(18) Re 1 2 Nu r = A 5 [ { 1 + Rd ( 1 + ( θ w 1 ) θ ( η ) ) 3 } θ ( η ) ] η = 0 ,

where Re r = Ω 1 r 2 υ f is the local Reynolds number.

3 Homotopy analysis method (HAM) solution

The initial guesses are defined as follows:

(19) f 0 ( η ) = A ( 1 η ) + b η , f 0 ( η ) = 1 Ω η 2 , θ 0 ( η ) = C + C D C + D + C D ( 1 η ) .

The linear operators are defined as follows:

(20) L f = f ( η ) f ( η ) , L g = g ( η ) g ( η ) , L θ = θ ( η ) θ ( η ) ,

with properties:

(21) L f ( C 1 + C 2 η ) = 0 , L g ( C 3 + C 4 η ) = 0 , L θ ( C 5 + C 6 η ) = 0 ,

where C 1 C 6 are the constants of the general solution. Detailed procedures performed relying on the framework can be seen in the study by Lone et al. [33].

3.1 Convergence analysis

The HAM technique is the most useful technique for the solution of nonlinear higher-order boundary value problems. The auxiliary parameter f plays a vital role in calculating the convergence region of the series solutions of a boundary value problem. Here, the auxiliary parameters f , g , and θ are calculated for the modeled problem as shown in Figure 2. From this figure, we see that the convergence region for the axial velocity profile is 1.85 f 0.23 , radial velocity profile is 2.00 g 0.25 , and temperature profile is 1.00 θ 0.23 .

Figure 2 
                  
                     
                        
                           
                           
                              h
                           
                           h
                        
                      curves for the function 
                        
                           
                           
                              f
                              ,
                              g
                              ,
                              θ
                           
                           f,g,\theta 
                        
                     .
Figure 2

h curves for the function f , g , θ .

3.2 Validation

Table 3 reveals the comparison of computed results of the current investigation with previously reported cases in the study by Irfan et al. [34] keeping We = Q = 0 , n = 3 as particularized in table. The comparative results showed a strong level of consistency and were in good agreement, which is strong support for the validity and dependability of the present findings. We can thus say with confidence that the method developed for examining the flow of Carreau HNF is highly trustworthy and may be used with considerable confidence.

Table 3

Computed values of Re 1 2 Nu r compared to some limiting cases when We = Q = 0 , n = 3

Pr Irfan et al. [34] Present results
0.7 0.45445767 0.454458
1.0 0.58201280 0.582013
1.3 0.69302203 0.693022
1.6 0.79227119 0.792272
1.8 0.85344101 0.853441
2.0 0.91134742 0.911348
5.0 1.56806380 1.568063
7.0 1.89542450 1.895425
10.0 2.30801360 2.308014
20.0 3.35396200 3.353962
70.0 6.46238930 6.462390

4 Results and discussion

A time-independent 3D model for Carreau HNF between two parallel, rotating, and stretching discs is generated. It is important to mention that two types of nanoparticles are imposed in the suspension of engine oil. Mechanisms of nonlinear radiation, heat, and Joule dissipation are implemented in the heat energy mechanism. Such a complex model is solved using a semi-analytical scheme called HAM. Detailed discussions regarding thermal and velocity profiles are mentioned in the following sections.

4.1 Velocity profiles

The velocity profiles are graphed using the dimensionless similarity variable η to anticipate the responses of fluids exhibiting shear-thinning behavior with n = 0.5 and fluids showing shear-thickening behavior with n = 1.5.

Figures 3 and 4 show the effect of radial and azimuthal velocities for different values of the stretching parameter at the lower disc (denoted as A ). When the stretching rate parameter A of the lower disc increases in stretchable rotatory discs, it has significant effects on the velocity profiles along the radial and tangential directions. As the bottom disk’s stretching rate parameter A increases, the material in the radial direction suffers more dramatic stretching or expansion. This expansion makes more room for the fluid to flow outward in a radial direction. As a result of the increased volume available for fluid passage, the fluid’s velocity increases along the radial path. In contrast, increasing the stretching rate parameter A can result in tangential contraction or a reduction in the azimuthal dimension of the bottom disc. This azimuthal contraction influences fluid flow along the azimuthal direction, resulting in a drop in velocities along this path. Changes in the dimensions of the lower disc have an effect on the overall flow patterns in the system. Increased radial expansion and azimuthal contraction can cause flow patterns to shift, potentially resulting in lower azimuthal velocities and increased radial velocities.

Figure 3 
                  Influence of 
                        
                           
                           
                              A
                           
                           A
                        
                      and 
                        
                           
                           
                              n
                           
                           n
                        
                      on 
                        
                           
                           
                              F
                              (
                              η
                              )
                           
                           F(\eta )
                        
                     .
Figure 3

Influence of A and n on F ( η ) .

Figure 4 
                  Influence of 
                        
                           
                           
                              A
                           
                           A
                        
                      and 
                        
                           
                           
                              n
                           
                           n
                        
                      on 
                        
                           
                           
                              G
                              (
                              η
                              )
                           
                           G(\eta )
                        
                     .
Figure 4

Influence of A and n on G ( η ) .

Figures 5 and 6 show the effect of radial and azimuthal velocities for different values of the stretching parameter at the upper disc (denoted as B ). When the stretching rate parameter B of the upper disc increases in the setting of stretchable rotatory discs, it causes distinct changes in the velocity profiles along the radial and tangential directions. An increase in the stretching rate parameter B results in increased radial stretching or expansion of the material. This radial expansion allows more fluid to flow away from the center of rotation. As a result of the increased volume available for fluid movement, the fluid’s velocity increases in the radial direction. The increased stretching of the top disc along the azimuthal direction causes azimuthal contraction or dimension decrease. The fluid flow within the disc is affected by this azimuthal contraction, resulting in a drop in velocities along the azimuthal direction. The changes in dimensions caused by increasing the stretching rate parameter B have an effect on the overall flow patterns inside the system. The observed drop in velocities along the azimuthal direction and rise in velocities along the radial direction may be due to variations in flow patterns caused by radial expansion and azimuthal contraction.

Figure 5 
                  Influence of 
                        
                           
                           
                              B
                           
                           B
                        
                      and 
                        
                           
                           
                              n
                           
                           n
                        
                      on 
                        
                           
                           
                              F
                              (
                              η
                              )
                           
                           F(\eta )
                        
                     .
Figure 5

Influence of B and n on F ( η ) .

Figure 6 
                  Influence of 
                        
                           
                           
                              B
                           
                           B
                        
                      and 
                        
                           
                           
                              n
                           
                           n
                        
                      on 
                        
                           
                           
                              G
                              (
                              η
                              )
                           
                           G(\eta )
                        
                     .
Figure 6

Influence of B and n on G ( η ) .

Figures 7 and 8 depict the influence of radial and azimuthal velocities, respectively, for varying values of the rotation parameter (designated as Ω). The rotation parameter Ω has a considerable influence on the fluid flow behavior in the context of stretchy rotatory discs. When the rotation parameter increases, the rotational speed or angular velocity of the disc increases. This increase in rotational speed has a number of consequences on the fluid flow within the disc. As the rotational speed of the disc increases, so does the centrifugal force experienced by the fluid. This force pulls the fluid outward in the radial direction, resulting in larger radial velocities. The greater the centrifugal force, the faster the fluid moves away from the center of the revolution. The circular path around the center of rotation corresponds to the tangential direction. When the rotation parameter is increased, the disc spins faster, and this increased rotational speed transfers more kinetic energy to the fluid particles. As a result, the fluid particles in the tangential direction gain speed, increasing the tangential velocity. The impacts of n = 0.5 and n = 1.5 exhibit nearly comparable outcomes.

Figure 7 
                  Influence of 
                        
                           
                           
                              Ω
                           
                           \text{Ω}
                        
                      and 
                        
                           
                           
                              n
                           
                           n
                        
                      on 
                        
                           
                           
                              F
                              (
                              η
                              )
                           
                           F(\eta )
                        
                     .
Figure 7

Influence of Ω and n on F ( η ) .

Figure 8 
                  Influence of 
                        
                           
                           
                              Ω
                           
                           \text{Ω}
                        
                      and 
                        
                           
                           
                              n
                           
                           n
                        
                      on 
                        
                           
                           
                              G
                              (
                              η
                              )
                           
                           G(\eta )
                        
                     .
Figure 8

Influence of Ω and n on G ( η ) .

Figure 9 and 10 depict the influence of radial and azimuthal velocities, respectively, for varying values of the magnetic parameter (designated as M ). The velocity profiles in stretchy rotatory discs decrease noticeably when the magnetic parameter increases in both the radial and tangential directions. In stretchable rotatory disks, the influence of a magnetic field affects the flow behavior of the fluid. The magnetic forces applied to the fluid grow more prominent as the magnetic parameter is increased. These forces interact with the fluid’s motion, causing velocities to decrease in both the radial and tangential directions. Simply put, when the magnetic parameter increases, so do the magnetic forces, which impede fluid movement in the radial and tangential directions. As a result, velocities in these directions are slower than in situations with less magnetic impact. Furthermore, the enhancement of viscous and Joule dissipation, thickening of the boundary layer, modifications in flow structure, and Lorentz force acting as a magnetic drag are the main causes of the decrease in radial and azimuthal velocities with increasing magnetic parameter. These effects combined lead to the suppression of fluid motion in the presence of a strong magnetic field. The impacts of n = 0.5 and n = 1.5 exhibit nearly comparable outcomes.

Figure 9 
                  Influence of 
                        
                           
                           
                              M
                           
                           M
                        
                      and 
                        
                           
                           
                              n
                           
                           n
                        
                      on 
                        
                           
                           
                              F
                              (
                              η
                              )
                           
                           F(\eta )
                        
                     .
Figure 9

Influence of M and n on F ( η ) .

Figure 10 
                  Influence of 
                        
                           
                           
                              M
                           
                           M
                        
                      and 
                        
                           
                           
                              n
                           
                           n
                        
                      on 
                        
                           
                           
                              G
                              (
                              η
                              )
                           
                           G(\eta )
                        
                     .
Figure 10

Influence of M and n on G ( η ) .

4.2 Temperature profiles

Figure 11 depicts the effect of the Biot number at lower disc parameter C on the nondimensional temperature distribution. When the temperature ratio ( θ w ) is 1 (showing equal temperatures in both regions) and the Biot number at the lower disc parameter C increases, convective heat transfer at the surface is more effective than internal heat conduction within the solid. At θ w = 1 , where the temperatures are the same in both locations, an increase in the Biot number indicates that convective heat transfer is becoming more important relative to internal heat conduction. Improved convective heat transfer promotes improved heat exchange between the solid and the surrounding fluid. As a result, the temperature profiles in this situation rise faster because heat is transported more efficiently from the solid to the fluid. When the temperature ratio ( θ w ) is 1.5 (showing a temperature difference between the regions) and the Biot number at the lower disc parameter C rises, it indicates that convective heat transfer is becoming more important in comparison to internal heat conduction within the solid. At θ w = 1.5 , the higher temperature in one part compared to the other produces a temperature gradient, which causes heat transfer from the hotter to the colder region. A higher Biot number, indicating greater convective heat transfer, would result in more effective heat exchange between the solid and the surrounding fluid. Because of the improved convection, the heat is transferred away from the solid more effectively, resulting in a decrease in temperature profiles. It is also noted that when the temperature ratio parameter θ w increases, the temperature profile also increases because heat naturally flows from higher-temperature regions to lower-temperature regions. The pushing power for heat transmission grows as the temperature difference between two places increases. As a result, more heat moves from the hotter to the cooler zone.

Figure 11 
                  Influence of 
                        
                           
                           
                              C
                           
                           C
                        
                      and 
                        
                           
                           
                              
                                 
                                    θ
                                 
                                 
                                    w
                                 
                              
                           
                           {\theta }_{\text{w}}
                        
                      on 
                        
                           
                           
                              θ
                              (
                              η
                              )
                           
                           \theta (\eta )
                        
                     .
Figure 11

Influence of C and θ w on θ ( η ) .

Figure 12 depicts the effect of the Biot number at upper disc parameter D on the nondimensional temperature distribution. When the Biot number at the upper disc parameter D in stretchy rotatory discs increases at both θ w = 1 and θ w = 1.5 , the temperature profiles fall due to increased convective heat transfer. At θ w = 1 (identical temperatures in both regions) and a rising Biot number at the upper disc parameter D , convective heat transfer at the upper disc surface becomes more important than internal heat conduction within the solid. This improvement in convective heat transmission results in more effective heat removal from the upper disc. As a result, the upper disk’s temperature profiles decline. The convective heat transfer method effectively transports heat away from the upper disk’s surface, resulting in a decrease in temperature. Because the lower disk’s temperature is the same as the top disk’s initial temperature, it is unaffected. At θ w = 1.5 (temperature difference between sections) and a rising Biot number at the upper disc parameter D , convective heat transfer at the upper disk’s surface is intensified in comparison to internal heat conduction within the solid. Because of the increased convective heat transfer, the temperature profiles of the upper disc decline more quickly. Because the temperature of the upper disc is higher than that of the lower disc, the convective process successfully transports heat away from the upper disc, causing its temperature to fall. The temperature of the lower disc is virtually unaffected since it remains constant or is altered to a lesser extent. The fundamental component in both circumstances is the dominance of convective heat transfer over internal heat conduction. An increase in the Biot number at the upper disc parameter D results in stronger convective heat transport from the upper disk’s surface, resulting in a fall in its temperature profiles.

Figure 12 
                  Influence of 
                        
                           
                           
                              D
                           
                           D
                        
                      and 
                        
                           
                           
                              
                                 
                                    θ
                                 
                                 
                                    w
                                 
                              
                           
                           {\theta }_{\text{w}}
                        
                      on 
                        
                           
                           
                              θ
                              (
                              η
                              )
                           
                           \theta (\eta )
                        
                     .
Figure 12

Influence of D and θ w on θ ( η ) .

Figure 13 indicates the impact of heat sink/source parameter Q on the nondimensional temperature distribution. Temperature profiles often increase as the heat sink/source parameter Q increases in stretchy rotatory discs at both θ w = 1 and θ w = 1.5 . This is because increasing the heat sink/source parameter Q causes more heat to be added from the system, affecting the overall temperature distribution inside the discs. A rise in the heat source parameter Q indicates that more heat is being introduced into the system. This additional heat input increases the thermal energy within the discs, causing them to heat up. Furthermore, the direct input of thermal energy to the fluid, increased internal heat generation, and the positive source term in the energy balance equation all contribute to the fluid’s temperature increase as the heat source parameter rises. Higher overall temperatures are achieved inside the fluid domain as a result of this shifting of the thermal equilibrium. The heat source helps to raise the temperature profiles in the top and lower discs.

Figure 13 
                  Influence of 
                        
                           
                           
                              Q
                           
                           Q
                        
                      and 
                        
                           
                           
                              
                                 
                                    θ
                                 
                                 
                                    w
                                 
                              
                           
                           {\theta }_{\text{w}}
                        
                      on 
                        
                           
                           
                              θ
                              (
                              η
                              )
                           
                           \theta (\eta )
                        
                     .
Figure 13

Influence of Q and θ w on θ ( η ) .

Figure 14 illustrates the impact of magnetic parameter M on the nondimensional temperature distribution. Temperature profiles tend to grow in stretchy rotatory discs as the magnetic parameter M increases for both θ w = 1 and θ w = 1.5 . This behavior is caused by the magnetic field’s influence on heat transfer and fluid dynamics within the system. A greater magnetic field within the system correlates with an increase in the magnetic parameter M . This magnetic field has the potential to influence fluid flow patterns and heat transfer rates. Heat transport is frequently boosted in the presence of a magnetic field due to changes in the convective and conductive processes. Furthermore, the fluid’s temperature rises with a rise in the magnetic parameter because of increased viscous dissipation, Joule heating, and changes in the flow structure brought about by the Lorentz force.

Figure 14 
                  Influence of 
                        
                           
                           
                              M
                           
                           M
                        
                      and 
                        
                           
                           
                              
                                 
                                    θ
                                 
                                 
                                    w
                                 
                              
                           
                           {\theta }_{\text{w}}
                        
                      on 
                        
                           
                           
                              θ
                              (
                              η
                              )
                           
                           \theta (\eta )
                        
                     .
Figure 14

Influence of M and θ w on θ ( η ) .

Figure 15 indicates the influence of Eckert number Ec on the nondimensional temperature distribution. Temperature profiles tend to increase in stretchy rotatory discs as the Eckert number ( Ec ) increases for both θ w = 1 and θ w = 1.5 . The Eckert number is a dimensionless metric that characterizes the balance of kinetic and thermal energy in a fluid flow by representing the ratio of kinetic energy change to enthalpy change. A rise in the Eckert number suggests that there is more kinetic energy present than enthalpy change. This increased kinetic energy can improve heat mixing and transfer throughout the fluid. In other words, a higher Eckert number indicates that the fluid flow’s kinetic energy has a greater impact on the entire heat transfer process. Moreover, the temperature rises with the growing Eckert number due to increased convective heat transfer in the fluid flow; this means that there is more kinetic energy available for convective processes to convert into thermal energy, which raises the temperature of the fluid.

Figure 15 
                  Influence of Ec and 
                        
                           
                           
                              
                                 
                                    θ
                                 
                                 
                                    w
                                 
                              
                           
                           {\theta }_{\text{w}}
                        
                      on 
                        
                           
                           
                              θ
                              (
                              η
                              )
                           
                           \theta (\eta )
                        
                     .
Figure 15

Influence of Ec and θ w on θ ( η ) .

Figure 16 indicates the effect of the Radiation parameter Rd on the nondimensional temperature distribution. Temperature profiles tend to grow when the radiation parameter ( Rd ) in stretchy rotatory discs increases for both θ w = 1 and θ w = 1.5 . This is because of the effect of radiative heat transfer on the overall thermal energy distribution within the system. Radiative heat transfer becomes more prevalent in the overall heat transfer process as the radiation parameter Rd grows. Radiative heat transfer is less impacted by fluid characteristics and may transport heat over greater distances, making it especially important in systems with substantial temperature differences.

Figure 16 
                  Influence of 
                        
                           
                           
                              Rd
                           
                           \text{Rd}
                        
                      and 
                        
                           
                           
                              
                                 
                                    θ
                                 
                                 
                                    w
                                 
                              
                           
                           {\theta }_{\text{w}}
                        
                      on 
                        
                           
                           
                              θ
                              (
                              η
                              )
                           
                           \theta (\eta )
                        
                     .
Figure 16

Influence of Rd and θ w on θ ( η ) .

4.3 Skin friction and local Nusselt number

Table 4 indicates the influence of skin friction for numerous values of We , M , A , B , and Ω. From Table 4, we observe that an increase in the Weissenberg number ( We ) may result in a drop in skin friction coefficients for a shear-thinning fluid ( n = 0.5 ) that becomes less viscous with increasing shear rate. This is because a larger Weissenberg number shows that elastic effects have a greater influence on flow behavior. The elasticity of the fluid increases as the Weissenberg number increases. Higher elasticity in this instance could result in the fluid’s capacity to restore its shape more effectively following deformation. This has the ability to minimize shear stress and, as a result, skin friction on a solid surface. The influence of the Weissenberg number on skin friction coefficients may differ for a shear-thickening fluid (n = 1.5), which gets more viscous with the increasing shear rate. The elastic effects become more evident as the Weissenberg number increases, potentially resulting in increased shear stresses and skin friction. Increased Weissenberg number may result in increased deformation resistance and a tendency for the fluid to “stick” to the surface, resulting in increased friction along both the radial and azimuthal axes.

Table 4

Numerical outcomes of skin friction coefficients in radial and azimuthal directions for change values of We , M , A , B , Ω

We M A B Ω Re 1 2 C Fr Re 1 2 C G θ
n = 0.5 n = 1.5 n = 0.5 n = 1.5
0.1 0.1 0.1 0.1 0.1 0.463546 0.486463 0.323567 0.335425
0.4 0.444567 0.469866 0.318643 0.328646
0.7 0.430023 0.440754 0.309753 0.314678
0.1 0.2 0.696537 0.756864 0.213674 0.243567
0.6 0.705336 0.765357 0.235885 0.268975
1.2 0.713466 0.775319 0.249075 0.287523
0.1 0.1 0.697647 0.895458 0.298654 0.375367
0.2 0.680963 0.885436 0.287633 0.368775
0.3 0.678537 0.875421 0.275327 0.357543
0.1 0.1 0.457906 0.557547 0.254364 0.397635
0.2 0.435685 0.540862 0.236853 0.389863
0.3 0.424614 0.537995 0.219985 0.369949
0.1 0.2 0.544263 0.676467 0.356786 0.487554
0.3 0.557754 0.665374 0.368933 0.497527
0.4 0.564336 0.654356 0.375658 0.510096

Table 4 also indicates that when the magnetic parameter M increases then the skin friction along the redial and azimuthal direction grows because, as the magnetic parameter (M) grows, so does the magnetic field’s influence on fluid velocity. The magnetic field can cause Lorentz forces to change the flow pattern, thus increasing resistance to motion. In this instance, the magnetic field may obstruct fluid flow in both the radial and azimuthal directions, leading to higher skin friction coefficients.

Table 4 indicates that when the stretching parameters A and B are increased, the discs are stretched more. This increased stretching might result in improved fluid motion and higher shear rates. Higher shear rates can result in lower apparent viscosity in shear-thinning and shear-thickening fluids. As a result, the flow resistance and skin friction coefficients may drop.

In Table 4, we see that when the rotational parameter increases in the radial direction, it often leads to higher angular velocities, which might result in increased centrifugal forces. Higher shear rates can cause a drop in apparent viscosity in the case of a shear-thinning fluid ( n = 0.5 ) . This decrease in viscosity reduces flow resistance, resulting in reduced radial skin friction coefficients. Higher shear rates can cause an increase in apparent viscosity in the case of a shear-thickening fluid ( n = 1.5 ) . This increased viscosity increases flow resistance, resulting in larger radial skin friction coefficients. Increases in the rotational parameter in the azimuthal direction frequently result in more fluid motion and higher angular velocities. Increased mobility might result in more aggressive mixing and higher shear rates. Higher shear rates can cause variations in apparent viscosity in both shear-thinning ( n = 0.5 ) and shear-thickening ( n = 1.5 ) fluids. The viscosity of shear-thinning fluids may decrease, thus reducing resistance and leading to reduced skin friction coefficients. The viscosity of shear-thickening fluids may increase, potentially leading to increased resistance and skin friction coefficients.

Table 4 illustrates the influence of the Nusselt number for different numbers of Rd , Ec , M , C , and D . Increasing the radiation parameter, Eckert number, magnetic parameter, and Biot number at both the bottom and upper discs improve convective heat transmission and fluid mixing, according to Table 4. Higher Nusselt numbers result from enhanced convective heat transfer for both θ w = 1 and θ w = 1.5 . Increasing the radiation parameter ( Rd ) frequently improves radiative heat transfer. This new heat transfer mechanism boosts the convective heat transfer component, resulting in higher Nusselt numbers for both θ w = 1 and θ w = 1.5 . Increased radiative heat transfer provides a more efficient heat exchange process, which results in higher Nu. The Eckert number ( Ec ) connects the change in kinetic energy to the change in enthalpy in fluid flow. An increase in Ec can improve kinetic energy, resulting in more vigorous fluid mixing and better heat transmission. This increases convective heat transfer and, as a result, Nusselt values for both θ w = 1 and θ w = 1.5 . The magnetic parameter ( M ) increases the magnetic field’s influence on fluid velocity. This improves convective mixing and heat transport, resulting in greater Nusselt numbers for θ w = 1 and θ w = 1.5 . The interaction of the magnetic field with the fluid can change flow patterns and promote more effective heat exchange. The ratio of convective heat transfer at the surface to conductive heat transfer within the material is represented by the Biot number. The greater the Biot number, the greater the contribution of convective heat transfer. This improves heat exchange between the fluid and the discs, resulting in greater Nusselt numbers for θ w = 1 and θ w = 1.5 (Table 5).

Table 5

Numerical outcomes of Nusselt number coefficients for change in values of Rd , Ec , M , C , D

Rd Ec M C D Re 1 2 Nu r
θ w = 1.0 θ w = 1.5
0.1 0.1 0.1 0.1 0.1 0.555472 0.657757
0.2 0.586473 0.697853
0.3 0.609967 0.734578
0.1 0.2 0.546784 0.567853
0.3 0.554383 0.578653
1.4 0.569076 0.588595
0.1 0.1 0.497654 0.556848
0.2 0.507543 0.567896
0.3 0.514952 0.579087
0.1 0.1 0.690758 0.887646
0.2 0.747856 0.946785
0.3 0.896528 0.999963
0.1 0.2 0.577973 0.698658
0.3 0.608743 0.747996
0.4 0.636795 0.804474

5 Concluding remarks

In this work, we analyzed the consequences of convective heat transport in the MHD flow of Carreau HNF between two disks that are rotating as well as stretching. The Tiwari and Das nanofluid model is employed to model the fluid flow dynamics in terms of PDEs by taking into account the thermal radiation, viscous dissipation, joule heating, and heat generation/absorption. The resultant PDEs are then transmuted into ODEs by defining appropriate similarity variables. The converted complex ODEs are then solved by HAM passed through MATLAB. The significant outcomes of the current analysis are as follows:

  • When the stretching parameter at the bottom disc A and at the upper disc B of a stretchy rotatory disc rises, the fluid velocity profiles alter for n = 0.5 and n = 1.5 . Radial expansion causes greater speeds in the radial direction, whereas azimuthal contraction causes lower velocities in the azimuthal direction.

  • Increases in the rotation parameter of stretchy rotatory discs result in increased radial and azimuthal velocities for both n = 0.5 and n = 1.5 .

  • For both n = 0.5 and n = 1.5 , increases in the magnetic parameter produce decreases in radial and azimuthal velocities in stretchy rotatory discs due to higher magnetic forces restricting fluid motion.

  • The temperature profile increases at θ w = 1 and decreases at θ w = 1.5 when the Biot number at lower disc C increases.

  • The temperature profile decreases at both θ w = 1 and θ w = 1.5 when the Biot number at the upper disc D increases.

  • The temperature profile increases at both θ w = 1 and θ w = 1.5 when the heat source/sink parameter, magnetic parameter, Eckert number, and radiation parameter increase.

  • The skin friction coefficient decreases for both radial and azimuthal directions at n = 0.5 and n = 1.5 as the Weissenberg number and the stretching parameters at both the lower and upper discs increase.

  • The skin friction coefficient increases for both radial and azimuthal directions at n = 0.5 and n = 1.5 as the magnetic parameter increases.

  • For n = 0.5 , the skin friction coefficients increase in the radial direction as the rotational parameter increases, while for n = 1.5 , they decrease. Conversely, in the azimuthal direction, skin friction coefficients increase for both n = 0.5 and n = 1.5 as the rotational parameter increases.

  • Raising the radiation parameter, Eckert number, magnetic parameter, and Biot number at both the bottom and upper discs improves the Nusselt numbers for both θ w = 1 and θ w = 1.5 .

6 Future research direction

In the future, the current study can be extended by considering entropy generation in ternary hybrid nanoparticles over various plane and rough surfaces in the presence of various physical effects such as stagnation point flow, porous medium, bioconvection, and activation energy subject to slip and convective-Nield boundary conditions.

Acknowledgments

This project was supported by Researchers Supporting Project number (RSP2024R411), King Saud University, Riyadh, Saudi Arabia.

  1. Funding information: The project was financed by Lucian Blaga University of Sibiu through research grant LBUS – IRG – 2023 – 09.

  2. Author contributions: All authors have accepted responsibility for the entire content of this manuscript and approved its submission.

  3. Conflict of interest: The authors state no conflict of interest.

  4. Data availability statement: The datasets generated and/or analyzed during the current study are available from the corresponding author on reasonable request.

References

[1] Kármán TV. Über laminare und turbulente Reibung. ZAMM - J Appl Math Mech/Z Angew Math Mech. Jan. 1921;1(4):233–52. 10.1002/ZAMM.19210010401.Search in Google Scholar

[2] Fang T, Zhang J. Flow between two stretchable disks – An exact solution of the Navier–Stokes equations. Int Commun Heat Mass Transf. Oct. 2008;35(8):892–5. 10.1016/J.ICHEATMASSTRANSFER.2008.04.018.Search in Google Scholar

[3] Imtiaz M, Hayat T, Alsaedi A, Ahmad B. Convective flow of carbon nanotubes between rotating stretchable disks with thermal radiation effects. Int J Heat Mass Transf. Oct. 2016;101:948–57. 10.1016/J.IJHEATMASSTRANSFER.2016.05.114.Search in Google Scholar

[4] Ahmed J, Khan M, Ahmad L. MHD swirling flow and heat transfer in Maxwell fluid driven by two coaxially rotating disks with variable thermal conductivity. Chin J Phys. Aug. 2019;60:22–34. 10.1016/J.CJPH.2019.02.010.Search in Google Scholar

[5] Usman U, Ghaffari A, Kausar S. Numerical solution of the partial differential equations that model the steady three-dimensional flow and heat transfer of Carreau fluid between two stretchable rotatory disks. Numer Methods Partial Differ Equ. Sep. 2023;39(5):3532–60. 10.1002/NUM.22672.Search in Google Scholar

[6] Usman, Ghaffari A, Muhammad T, Mustafa I. Heat transfer enhancement in a power-law nanofluid flow between two rotating stretchable disks. Pramana - J Phys. Mar. 2022;96(1):1–11. 10.1007/S12043-021-02272-0/METRICS.Search in Google Scholar

[7] Das S, Ali A. Dynamical phenomena developed by a spiralling stretchable sheet in magnetized Casson-spinel ferrite nanofluid. Heliyon. Aug. 2023;9:e18376. 10.1016/j.heliyon.2023.e18376.Search in Google Scholar PubMed PubMed Central

[8] Das S, Das S. EDL Aspects in swirling ionic tribological fluid flow in a squeezed/split channel underlie a high-power magnetic field. Forces Mech. May. 2023;11:100196. 10.1016/j.finmec.2023.100196.Search in Google Scholar

[9] Yashkun U, Zaimi K, Abu Bakar NA, Ishak A, Pop I. MHD hybrid nanofluid flow over a permeable stretching/shrinking sheet with thermal radiation effect. Int J Numer Methods Heat Fluid Flow. 2021;31(3):1014–31.10.1108/HFF-02-2020-0083Search in Google Scholar

[10] Lund LA, Yashkun U, Shah NA. Multiple solutions of unsteady Darcy–Forchheimer porous medium flow of Cu–Al2O3/water based hybrid nanofluid with joule heating and viscous dissipation effect. J Therm Anal Calorim. 2024;149:2303–15.10.1007/s10973-023-12819-5Search in Google Scholar

[11] Waini I, Ishak A, Pop I. Flow and heat transfer of a hybrid nanofluid past a permeable moving surface. Chin J Phys. Aug. 2020;66:606–19. 10.1016/J.CJPH.2020.04.024.Search in Google Scholar

[12] Zainal NA, Nazar R, Naganthran K, Pop I. MHD mixed convection stagnation point flow of a hybrid nanofluid past a vertical flat plate with convective boundary condition. Chin J Phys. Aug. 2020;66:630–44. 10.1016/J.CJPH.2020.03.022.Search in Google Scholar

[13] Waini I, Ishak A, Groşan T, Pop I. Mixed convection of a hybrid nanofluid flow along a vertical surface embedded in a porous medium. Int Commun Heat Mass Transf. May 2020;114:104565. 10.1016/J.ICHEATMASSTRANSFER.2020.104565.Search in Google Scholar

[14] Rostami MN, Dinarvand S, Pop I. Dual solutions for mixed convective stagnation-point flow of an aqueous silica–alumina hybrid nanofluid. Chin J Phys. Oct. 2018;56(5):2465–78. 10.1016/J.CJPH.2018.06.013.Search in Google Scholar

[15] Shamshuddin MD, Saeed A, Asogwa KK, Jamshed W. A semi-analytical approach to investigate the entropy generation in a tangent hyperbolic magnetized hybrid nanofluid flow upon a stretchable rotating disk. J Magn Magn Mater. May 2023;574:170664. 10.1016/J.JMMM.2023.170664.Search in Google Scholar

[16] Farooq U, Imran M, Fatima N, Noreen S, Tahir M, Akgül A, et al. Cattaneo-Christov heat flux model in radiative flow of (Fe3O4– TiO2/Transformer oil) and (Cu–TiO2/Transformer oil) magnetized hybrid nanofluids past through double rotating disks. Case Stud Therm Eng. May 2023;45:102905. 10.1016/J.CSITE.2023.102905.Search in Google Scholar

[17] Agrawal R, Kaswan P. Entropy generation minimization of Ag-Fe3O4/water-ethylene glycol squeezed hybrid nanofluid flow between parallel disks . Int J Numer Methods Heat Fluid Flow. Jan. 2023;33(1):65–95. 10.1108/HFF-01-2022-0005/FULL/XML.Search in Google Scholar

[18] Waini I, Ishak A, Pop I. Unsteady flow and heat transfer past a stretching/shrinking sheet in a hybrid nanofluid. Int J Heat Mass Transf. Jun. 2019;136:288–97. 10.1016/J.IJHEATMASSTRANSFER.2019.02.101.Search in Google Scholar

[19] Khashi’ie NS, Arifin NM, Pop I, Nazar R, Hafidzuddin EH, Wahi N. Three-dimensional hybrid nanofluid flow and heat transfer past a permeable stretching/shrinking sheet with velocity slip and convective condition. Chin J Phys. Aug. 2020;66:157–71. 10.1016/J.CJPH.2020.03.032.Search in Google Scholar

[20] Anuar NS, Bachok N, Pop I. Influence of buoyancy force on Ag-MgO/water hybrid nanofluid flow in an inclined permeable stretching/shrinking sheet. Int Commun Heat Mass Transf. Apr. 2021;123:105236. 10.1016/J.ICHEATMASSTRANSFER.2021.105236.Search in Google Scholar

[21] Rooman M, Jan MA, Shah Z, Vrinceanu N. Entropy optimization on axisymmetric darcy-forchheimer powell-eyring nanofluid over a horizontally stretching cylinder with viscous dissipation effect. Coatings. May 2022;12(6):749. 10.3390/COATINGS12060749.Search in Google Scholar

[22] Rooman M, Jan MA, Shah Z, Alzahrani MR. Entropy generation and nonlinear thermal radiation analysis on axisymmetric MHD Ellis nanofluid over a horizontally permeable stretching cylinder. Waves Random Complex Media. 2022;1–15. 10.1080/17455030.2021.2020934.Search in Google Scholar

[23] Shutaywi M, Rooman M, Jan MA, Vrinceanu N, Shah Z, Deebani W. Entropy generation and thermal analysis on MHD second-grade fluid with variable thermophysical properties over a stratified permeable surface of paraboloid revolution. ACS Omega. Aug. 2022;7(31):27436–49. 10.1021/ACSOMEGA.2C02452.Search in Google Scholar PubMed PubMed Central

[24] Carreau PJ. Rheological equations from molecular network theories. Trans Soc Rheol. Mar. 1972;16(1):99–127. 10.1122/1.549276.Search in Google Scholar

[25] Sulochana C, Ashwinkumar GP, Sandeep N. Transpiration effect on stagnation-point flow of a Carreau nanofluid in the presence of thermophoresis and Brownian motion. Alex Eng J. Jun. 2016;55(2):1151–7. 10.1016/J.AEJ.2016.03.031.Search in Google Scholar

[26] Khan M, Irfan M, Khan WA, Alshomrani AS. A new modeling for 3D Carreau fluid flow considering nonlinear thermal radiation. Results Phys. Jan. 2017;7:2692–704. 10.1016/J.RINP.2017.07.024.Search in Google Scholar

[27] Irfan M, Khan M, Khan WA. Numerical analysis of unsteady 3D flow of Carreau nanofluid with variable thermal conductivity and heat source/sink. Results Phys. Jan. 2017;7:3315–24. 10.1016/J.RINP.2017.08.029.Search in Google Scholar

[28] Rooman M, Shah Z, Bonyah E, Jan MA, Deebani W. Mathematical modeling of Carreau fluid flow and heat transfer characteristics in the renal tubule. J Math. May. 2022;2022:2517933. 10.1155/2022/2517933.Search in Google Scholar

[29] Abbas IA. Generalized thermoelastic interaction in functional graded material with fractional order three-phase lag heat transfer. J Central South Univ. May. 2015;22(5):1606–13. 10.1007/s11771-015-2677-5.Search in Google Scholar

[30] Abbas I, Hobiny A, Marin M. Photo-thermal interactions in a semi-conductor material with cylindrical cavities and variable thermal conductivity. J Taibah Univ Sci. Sept. 2020;14(1):1369–76. 10.1080/16583655.2020.1824465.Search in Google Scholar

[31] Marin M, Hobiny A, Abbas I. The effects of fractional time derivatives in porothermoelastic materials using finite element method. Mathematics. June 2021;9(14):1606. 10.1080/16583655.2020.1824465.Search in Google Scholar

[32] Takabi B, Salehi S. Augmentation of the heat transfer performance of a sinusoidal corrugated enclosure by employing hybrid nanofluid. Adv Mech Eng. Jan. 2014;6:147059. 10.1155/2014/147059.Search in Google Scholar

[33] Lone SA, Shamshuddin MD, Shahab S, Iftikhar S, Saeed A, Galal AM. Computational analysis of MHD driven bioconvective flow of hybrid Casson nanofluid past a permeable exponential stretching sheet with thermophoresis and Brownian motion effects. J Magn Magn Mater. Aug. 2023;580:170959. 10.1016/J.JMMM.2023.170959.Search in Google Scholar

[34] Irfan M, Khan M, Khan WA. 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 Phys. Sep. 2018;10:107–17. 10.1016/J.RINP.2018.04.036.Search in Google Scholar

Received: 2024-01-27
Revised: 2024-07-18
Accepted: 2024-07-28
Published Online: 2024-11-04

© 2024 the author(s), published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

  1. Regular Articles
  2. Numerical study of flow and heat transfer in the channel of panel-type radiator with semi-detached inclined trapezoidal wing vortex generators
  3. Homogeneous–heterogeneous reactions in the colloidal investigation of Casson fluid
  4. High-speed mid-infrared Mach–Zehnder electro-optical modulators in lithium niobate thin film on sapphire
  5. Numerical analysis of dengue transmission model using Caputo–Fabrizio fractional derivative
  6. Mononuclear nanofluids undergoing convective heating across a stretching sheet and undergoing MHD flow in three dimensions: Potential industrial applications
  7. Heat transfer characteristics of cobalt ferrite nanoparticles scattered in sodium alginate-based non-Newtonian nanofluid over a stretching/shrinking horizontal plane surface
  8. The electrically conducting water-based nanofluid flow containing titanium and aluminum alloys over a rotating disk surface with nonlinear thermal radiation: A numerical analysis
  9. Growth, characterization, and anti-bacterial activity of l-methionine supplemented with sulphamic acid single crystals
  10. A numerical analysis of the blood-based Casson hybrid nanofluid flow past a convectively heated surface embedded in a porous medium
  11. Optoelectronic–thermomagnetic effect of a microelongated non-local rotating semiconductor heated by pulsed laser with varying thermal conductivity
  12. Thermal proficiency of magnetized and radiative cross-ternary hybrid nanofluid flow induced by a vertical cylinder
  13. Enhanced heat transfer and fluid motion in 3D nanofluid with anisotropic slip and magnetic field
  14. Numerical analysis of thermophoretic particle deposition on 3D Casson nanofluid: Artificial neural networks-based Levenberg–Marquardt algorithm
  15. Analyzing fuzzy fractional Degasperis–Procesi and Camassa–Holm equations with the Atangana–Baleanu operator
  16. Bayesian estimation of equipment reliability with normal-type life distribution based on multiple batch tests
  17. Chaotic control problem of BEC system based on Hartree–Fock mean field theory
  18. Optimized framework numerical solution for swirling hybrid nanofluid flow with silver/gold nanoparticles on a stretching cylinder with heat source/sink and reactive agents
  19. Stability analysis and numerical results for some schemes discretising 2D nonconstant coefficient advection–diffusion equations
  20. Convective flow of a magnetohydrodynamic second-grade fluid past a stretching surface with Cattaneo–Christov heat and mass flux model
  21. Analysis of the heat transfer enhancement in water-based micropolar hybrid nanofluid flow over a vertical flat surface
  22. Microscopic seepage simulation of gas and water in shale pores and slits based on VOF
  23. Model of conversion of flow from confined to unconfined aquifers with stochastic approach
  24. Study of fractional variable-order lymphatic filariasis infection model
  25. Soliton, quasi-soliton, and their interaction solutions of a nonlinear (2 + 1)-dimensional ZK–mZK–BBM equation for gravity waves
  26. Application of conserved quantities using the formal Lagrangian of a nonlinear integro partial differential equation through optimal system of one-dimensional subalgebras in physics and engineering
  27. Nonlinear fractional-order differential equations: New closed-form traveling-wave solutions
  28. Sixth-kind Chebyshev polynomials technique to numerically treat the dissipative viscoelastic fluid flow in the rheology of Cattaneo–Christov model
  29. Some transforms, Riemann–Liouville fractional operators, and applications of newly extended M–L (p, s, k) function
  30. Magnetohydrodynamic water-based hybrid nanofluid flow comprising diamond and copper nanoparticles on a stretching sheet with slips constraints
  31. Super-resolution reconstruction method of the optical synthetic aperture image using generative adversarial network
  32. A two-stage framework for predicting the remaining useful life of bearings
  33. Influence of variable fluid properties on mixed convective Darcy–Forchheimer flow relation over a surface with Soret and Dufour spectacle
  34. Inclined surface mixed convection flow of viscous fluid with porous medium and Soret effects
  35. Exact solutions to vorticity of the fractional nonuniform Poiseuille flows
  36. In silico modified UV spectrophotometric approaches to resolve overlapped spectra for quality control of rosuvastatin and teneligliptin formulation
  37. Numerical simulations for fractional Hirota–Satsuma coupled Korteweg–de Vries systems
  38. Substituent effect on the electronic and optical properties of newly designed pyrrole derivatives using density functional theory
  39. A comparative analysis of shielding effectiveness in glass and concrete containers
  40. Numerical analysis of the MHD Williamson nanofluid flow over a nonlinear stretching sheet through a Darcy porous medium: Modeling and simulation
  41. Analytical and numerical investigation for viscoelastic fluid with heat transfer analysis during rollover-web coating phenomena
  42. Influence of variable viscosity on existing sheet thickness in the calendering of non-isothermal viscoelastic materials
  43. Analysis of nonlinear fractional-order Fisher equation using two reliable techniques
  44. Comparison of plan quality and robustness using VMAT and IMRT for breast cancer
  45. Radiative nanofluid flow over a slender stretching Riga plate under the impact of exponential heat source/sink
  46. Numerical investigation of acoustic streaming vortices in cylindrical tube arrays
  47. Numerical study of blood-based MHD tangent hyperbolic hybrid nanofluid flow over a permeable stretching sheet with variable thermal conductivity and cross-diffusion
  48. Fractional view analytical analysis of generalized regularized long wave equation
  49. Dynamic simulation of non-Newtonian boundary layer flow: An enhanced exponential time integrator approach with spatially and temporally variable heat sources
  50. Inclined magnetized infinite shear rate viscosity of non-Newtonian tetra hybrid nanofluid in stenosed artery with non-uniform heat sink/source
  51. Estimation of monotone α-quantile of past lifetime function with application
  52. Numerical simulation for the slip impacts on the radiative nanofluid flow over a stretched surface with nonuniform heat generation and viscous dissipation
  53. Study of fractional telegraph equation via Shehu homotopy perturbation method
  54. An investigation into the impact of thermal radiation and chemical reactions on the flow through porous media of a Casson hybrid nanofluid including unstable mixed convection with stretched sheet in the presence of thermophoresis and Brownian motion
  55. Establishing breather and N-soliton solutions for conformable Klein–Gordon equation
  56. An electro-optic half subtractor from a silicon-based hybrid surface plasmon polariton waveguide
  57. CFD analysis of particle shape and Reynolds number on heat transfer characteristics of nanofluid in heated tube
  58. Abundant exact traveling wave solutions and modulation instability analysis to the generalized Hirota–Satsuma–Ito equation
  59. A short report on a probability-based interpretation of quantum mechanics
  60. Study on cavitation and pulsation characteristics of a novel rotor-radial groove hydrodynamic cavitation reactor
  61. Optimizing heat transport in a permeable cavity with an isothermal solid block: Influence of nanoparticles volume fraction and wall velocity ratio
  62. Linear instability of the vertical throughflow in a porous layer saturated by a power-law fluid with variable gravity effect
  63. Thermal analysis of generalized Cattaneo–Christov theories in Burgers nanofluid in the presence of thermo-diffusion effects and variable thermal conductivity
  64. A new benchmark for camouflaged object detection: RGB-D camouflaged object detection dataset
  65. Effect of electron temperature and concentration on production of hydroxyl radical and nitric oxide in atmospheric pressure low-temperature helium plasma jet: Swarm analysis and global model investigation
  66. Double diffusion convection of Maxwell–Cattaneo fluids in a vertical slot
  67. Thermal analysis of extended surfaces using deep neural networks
  68. Steady-state thermodynamic process in multilayered heterogeneous cylinder
  69. Multiresponse optimisation and process capability analysis of chemical vapour jet machining for the acrylonitrile butadiene styrene polymer: Unveiling the morphology
  70. Modeling monkeypox virus transmission: Stability analysis and comparison of analytical techniques
  71. Fourier spectral method for the fractional-in-space coupled Whitham–Broer–Kaup equations on unbounded domain
  72. The chaotic behavior and traveling wave solutions of the conformable extended Korteweg–de-Vries model
  73. Research on optimization of combustor liner structure based on arc-shaped slot hole
  74. Construction of M-shaped solitons for a modified regularized long-wave equation via Hirota's bilinear method
  75. Effectiveness of microwave ablation using two simultaneous antennas for liver malignancy treatment
  76. Discussion on optical solitons, sensitivity and qualitative analysis to a fractional model of ion sound and Langmuir waves with Atangana Baleanu derivatives
  77. Reliability of two-dimensional steady magnetized Jeffery fluid over shrinking sheet with chemical effect
  78. Generalized model of thermoelasticity associated with fractional time-derivative operators and its applications to non-simple elastic materials
  79. Migration of two rigid spheres translating within an infinite couple stress fluid under the impact of magnetic field
  80. A comparative investigation of neutron and gamma radiation interaction properties of zircaloy-2 and zircaloy-4 with consideration of mechanical properties
  81. New optical stochastic solutions for the Schrödinger equation with multiplicative Wiener process/random variable coefficients using two different methods
  82. Physical aspects of quantile residual lifetime sequence
  83. Synthesis, structure, IV characteristics, and optical properties of chromium oxide thin films for optoelectronic applications
  84. Smart mathematically filtered UV spectroscopic methods for quality assurance of rosuvastatin and valsartan from formulation
  85. A novel investigation into time-fractional multi-dimensional Navier–Stokes equations within Aboodh transform
  86. Homotopic dynamic solution of hydrodynamic nonlinear natural convection containing superhydrophobicity and isothermally heated parallel plate with hybrid nanoparticles
  87. A novel tetra hybrid bio-nanofluid model with stenosed artery
  88. Propagation of traveling wave solution of the strain wave equation in microcrystalline materials
  89. Innovative analysis to the time-fractional q-deformed tanh-Gordon equation via modified double Laplace transform method
  90. A new investigation of the extended Sakovich equation for abundant soliton solution in industrial engineering via two efficient techniques
  91. New soliton solutions of the conformable time fractional Drinfel'd–Sokolov–Wilson equation based on the complete discriminant system method
  92. Irradiation of hydrophilic acrylic intraocular lenses by a 365 nm UV lamp
  93. Inflation and the principle of equivalence
  94. The use of a supercontinuum light source for the characterization of passive fiber optic components
  95. Optical solitons to the fractional Kundu–Mukherjee–Naskar equation with time-dependent coefficients
  96. A promising photocathode for green hydrogen generation from sanitation water without external sacrificing agent: silver-silver oxide/poly(1H-pyrrole) dendritic nanocomposite seeded on poly-1H pyrrole film
  97. Photon balance in the fiber laser model
  98. Propagation of optical spatial solitons in nematic liquid crystals with quadruple power law of nonlinearity appears in fluid mechanics
  99. Theoretical investigation and sensitivity analysis of non-Newtonian fluid during roll coating process by response surface methodology
  100. Utilizing slip conditions on transport phenomena of heat energy with dust and tiny nanoparticles over a wedge
  101. Bismuthyl chloride/poly(m-toluidine) nanocomposite seeded on poly-1H pyrrole: Photocathode for green hydrogen generation
  102. Infrared thermography based fault diagnosis of diesel engines using convolutional neural network and image enhancement
  103. On some solitary wave solutions of the Estevez--Mansfield--Clarkson equation with conformable fractional derivatives in time
  104. Impact of permeability and fluid parameters in couple stress media on rotating eccentric spheres
  105. Review Article
  106. Transformer-based intelligent fault diagnosis methods of mechanical equipment: A survey
  107. Special Issue on Predicting pattern alterations in nature - Part II
  108. A comparative study of Bagley–Torvik equation under nonsingular kernel derivatives using Weeks method
  109. On the existence and numerical simulation of Cholera epidemic model
  110. Numerical solutions of generalized Atangana–Baleanu time-fractional FitzHugh–Nagumo equation using cubic B-spline functions
  111. Dynamic properties of the multimalware attacks in wireless sensor networks: Fractional derivative analysis of wireless sensor networks
  112. Prediction of COVID-19 spread with models in different patterns: A case study of Russia
  113. Study of chronic myeloid leukemia with T-cell under fractal-fractional order model
  114. Accumulation process in the environment for a generalized mass transport system
  115. Analysis of a generalized proportional fractional stochastic differential equation incorporating Carathéodory's approximation and applications
  116. Special Issue on Nanomaterial utilization and structural optimization - Part II
  117. Numerical study on flow and heat transfer performance of a spiral-wound heat exchanger for natural gas
  118. Study of ultrasonic influence on heat transfer and resistance performance of round tube with twisted belt
  119. Numerical study on bionic airfoil fins used in printed circuit plate heat exchanger
  120. Improving heat transfer efficiency via optimization and sensitivity assessment in hybrid nanofluid flow with variable magnetism using the Yamada–Ota model
  121. Special Issue on Nanofluids: Synthesis, Characterization, and Applications
  122. Exact solutions of a class of generalized nanofluidic models
  123. Stability enhancement of Al2O3, ZnO, and TiO2 binary nanofluids for heat transfer applications
  124. Thermal transport energy performance on tangent hyperbolic hybrid nanofluids and their implementation in concentrated solar aircraft wings
  125. Studying nonlinear vibration analysis of nanoelectro-mechanical resonators via analytical computational method
  126. Numerical analysis of non-linear radiative Casson fluids containing CNTs having length and radius over permeable moving plate
  127. Two-phase numerical simulation of thermal and solutal transport exploration of a non-Newtonian nanomaterial flow past a stretching surface with chemical reaction
  128. Natural convection and flow patterns of Cu–water nanofluids in hexagonal cavity: A novel thermal case study
  129. Solitonic solutions and study of nonlinear wave dynamics in a Murnaghan hyperelastic circular pipe
  130. Comparative study of couple stress fluid flow using OHAM and NIM
  131. Utilization of OHAM to investigate entropy generation with a temperature-dependent thermal conductivity model in hybrid nanofluid using the radiation phenomenon
  132. Slip effects on magnetized radiatively hybridized ferrofluid flow with acute magnetic force over shrinking/stretching surface
  133. Significance of 3D rectangular closed domain filled with charged particles and nanoparticles engaging finite element methodology
  134. Robustness and dynamical features of fractional difference spacecraft model with Mittag–Leffler stability
  135. Characterizing magnetohydrodynamic effects on developed nanofluid flow in an obstructed vertical duct under constant pressure gradient
  136. Study on dynamic and static tensile and puncture-resistant mechanical properties of impregnated STF multi-dimensional structure Kevlar fiber reinforced composites
  137. Thermosolutal Marangoni convective flow of MHD tangent hyperbolic hybrid nanofluids with elastic deformation and heat source
  138. Investigation of convective heat transport in a Carreau hybrid nanofluid between two stretchable rotatory disks
  139. Single-channel cooling system design by using perforated porous insert and modeling with POD for double conductive panel
  140. Special Issue on Fundamental Physics from Atoms to Cosmos - Part I
  141. Pulsed excitation of a quantum oscillator: A model accounting for damping
  142. Review of recent analytical advances in the spectroscopy of hydrogenic lines in plasmas
  143. Heavy mesons mass spectroscopy under a spin-dependent Cornell potential within the framework of the spinless Salpeter equation
  144. Coherent manipulation of bright and dark solitons of reflection and transmission pulses through sodium atomic medium
  145. Effect of the gravitational field strength on the rate of chemical reactions
  146. The kinetic relativity theory – hiding in plain sight
  147. Special Issue on Advanced Energy Materials - Part III
  148. Eco-friendly graphitic carbon nitride–poly(1H pyrrole) nanocomposite: A photocathode for green hydrogen production, paving the way for commercial applications
Downloaded on 4.2.2026 from https://www.degruyterbrill.com/document/doi/10.1515/phys-2024-0078/html?lang=en
Scroll to top button