A numerical approach to the modeling of Thompson and Troian slip on magnetized flow of Al 2 O 3 – PAO nanolubricant over an inclined rotating disk

Abstract

In automotive fluids, hydraulic, gear, and bearing oils, as well as in applications operating in extremely high or cold temperatures, PAO is widely employed. In present work, we have made an attempt to develop a mathematical model to discuss the flow of magnetized $A l_{2} O_{3} - PAO$ nanolubricant over an inclined rotating disk in Darcy-Forchheimer porous medium with Thompson and Troian slip at the boundary. The effects of mixed convection, nonlinear heat radiation, viscous dissipation, Joule heating, and non-uniform heat source/sink are also included in the modeling. We have solved the proposed model numerically with the help of MATLAB built-in bvp-4c method after metamorphosing the PDEs to ODEs. The enhancing values of inertial parameter and velocity slip parameter decrease the tangential and radial velocities of the nanolubricant. The temperature of the nanolubricant $A l_{2} O_{3} - PAO$ enhance significantly by strengthening the magnetic field, whereas radial and tangential velocities get retarded. The non-uniform heat source/sink parameters play a vital role in controlling heat transmission phenomenon. The increasing values of Eckert number, radiation parameter, and non-uniform heat generation parameters tend to increase the value of Nusselt number. The value of Nusselt number drops with rising values of Biot number and non-uniform heat sink parameters.

Keywords

Darcy Forchheimer flow non-uniform heat source/sink Thompson and Troian slip convective conditions

Introduction

Lubricants usually include additives and might be based on mineral, synthetic, or vegetable oils. To lessen friction and wear, they are utilized as low-shear-resistance coatings between affecting surfaces. In addition to enhancing opposition to corrosion temperature and oxidation, the major goal of additives is to boost the performance and longevity of lubricants. Synthetic lubricants which have the ability to maintain their features at high temperatures without decomposing or combusting or experiencing chemical degradation have been developed in recent years in response to the need for high performance lubricants in engineering fields, primarily in space and the automotive industry.^1–6 Polyalphaolefins (PAOs) are among the synthetic lubricants which have gained more attention because of their intrinsic attributes such as wide operating temperature range, thermal stability and better oxidation, low viscosity, biodegradability, lower volatility, and high viscosity indices. PAOs have become new standard as a lubricant; according to the consumer reports, $70 %$ of the vehicles having either PAO or PAO based blended oils. The need for lubricants that efficiently reduce friction and wear under demanding circumstances while being less detrimental to the environment has also been supported by governmental agencies.⁷ Recent research indicates that nanotechnology can enhance lubricant enactment by lowering friction and system wear. The diffusion of nanoparticles in a base lubricant is known as a “nanolubricant,” which has caught the attention of researchers. Nanoparticles typically have diameters between 1 and 100 nm. The addition of nanoparticles to base oils has the potential to improve certain tribological qualities like friction and wear-resistance. The most significant aspect is that the majority of these particles are environmentally harmless. Despite being a relatively new field of research, nanotechnology has distinguished itself and demonstrated its potential in a number of fields, including tribology. The foremost purposes of using nanoparticles as lubricant additives are to give oil anti-wear qualities, modify friction, and increase resilience to extreme pressure.⁸ When employed in lubricants, nanoparticles can also facilitate heat transmission and bring down temperatures in the areas where they are used.^9,10 The nanolubricants have been attracting the attention of researchers and scientists due to their tremendous applications. Choi et al.¹¹ discussed heat dissipation experimentally and theoretically in nanolubricant evaporators. Gamaoun et al.¹² investigated conduction-convection heat transfer of a fin wetted with nanolubricant $ZnO - SAE 50$ . Riaz et al.¹³ discussed the flow of $ZnO - SAE 50$ nanolubricant over a Riga plate. Nayak et al.^14–16 studied transmission of heat in the flow of $ZnO - SAE 50$ nanolubricant over various surfaces. Khan¹⁷ analyzed the consequence of bioconvection in the flow of second grade nanofluid comprising nanoparticles and gyrotactic microorganisms. Zuhra et al.¹⁸ presented bioconvection’s simulation in the suspension of second grade nanofluid consisting of gyrotactic microorganisms and nanoparticles.

Magneto hydrodynamics (MHD) is the study of magnetic effects and the behavior of electrically conducting fluids. The fundamental theory underlying MHD holds that magnetic fields cause currents to pass through conductive fluids that are in motion. With the use of MHD, researchers have discovered a method for solving some of the issues encountered in engineering. Plasma confinement, nuclear reactor cooling with liquid metal, and electromagnetic casting are a few of these issues. MHD issues can arise in a wide range of circumstances, such as when forecasting space weather, controlling turbulent oscillations in semiconductor melt during crystal development, and measuring beverage flow rates in the food business. The studies^19–27 explain the flow of magnetized fluids over diverse geometries. Khan et al.²⁸ investigated second grade fluid’s thin film flow over a stretching sheet in porous medium along with heat transfer. The effects of heat source/sink on transfer of heat captured the attention of many researchers. Hsiao²⁹ analyzed MHD flow of viscoelastic fluid across a stretchable sheet with non-uniform heat source/sink. Ramesh et al.³⁰ have examined MHD flow across an inclined surface embedded in an incompressible fluid with the effect of non-uniform heat source/sink. Ramandevi et al.³¹ discussed combined effects of viscous dissipation and non-uniform heat sink/source on MHD flow of the fluid. Sulochana et al.³² examined the impact of non-uniform heat sink or source in $3 D$ flow of Casson fluid. Palwasha et al.³³ conducted a study to examine thin film fluid’s boundary layer flow having variable properties. Radiative heat, whether linear or nonlinear, has a substantial impact on a range of processes involving extremely high temperatures due to its many industrial and power plant applications. In business, the heat transfer procedure is crucial and necessary for final items to have the right qualities. Modern systems connected to spacecraft, power generator plasma, astrophysical flow, reactor cooling etc. are developed based on radiation applications. Khan et al.³⁴ discussed the impacts of thermal radiation and thermophoresis with mass and heat transfer in a magnetized thin film flow of second grade fluid having variable properties over a stretching sheet. Some recent articles discussing the effects of nonlinear thermal radiation can be seen in the references.^35–39

A thorough review of the literature reveals that no study is yet done to examine the flow characteristics and heat transference rate in the flow of magnetized $A l_{2} 0_{3} - PAO$ nanolubricant over a rotating inclined disk. The prime motive of this investigation is to fill this gap and introduce $A l_{2} 0_{3} - PAO$ nanolubricant over a rotating inclined disk subjected to Thompson and Troian slip and convective conditions. The flow takes place in Darcy-Forchheimer porous medium. The impacts of mixed convection, nonlinear heat radiation, viscous dissipation, Joule heating, and non-uniform heat source/sink are also included in the modeling. The proposed model is solved numerically via MATLAB built-in bvp-4c method after metamorphosing the PDEs to ODEs. The effects of imperative parameters on concerned profiles have been yakked via graphs.

Mathematical formulation

In current work, a mixed convective flow of magnetized $A l_{2} 0_{3} - PAO$ nanolubricant over a rotating inclined disk in a Darcy Forchheimer porous medium has been considered. At $z = 0$ , the disk swirls with uniform angular velocity $Ω$ . The components $(u, v, w)$ of velocity are taken along $(r, ϕ, z)$ directions respectively. In axial direction, a magnetic field having strength $B_{0}$ is implemented. The problem’s axial symmetry prevents us from including derivatives along $ϕ$ coordinate. The influences of non-uniform heat sink/source, nonlinear heat radiation, viscous dissipation, and Joule heating are considered in the modeling. The flow is subjected to the effects of Thompson and Troian slip and convective conditions. The depiction of problem is made in Figure 1.

Figure 1.

Flow geometry.

Governing equations

The flow phenomenon is expressed in this problem in the following ways^40,41:

\frac{\partial u}{\partial r} + \frac{u}{r} + \frac{\partial w}{\partial z} = 0,

(1)

\begin{matrix} u \frac{\partial u}{\partial r} - \frac{v^{2}}{r} + w \frac{\partial u}{\partial z} = - \frac{1}{ρ_{nf}} \frac{\partial p}{\partial r} + \frac{μ_{nf}}{ρ_{nf}} [\frac{\partial^{2} u}{\partial r^{2}} + \frac{1}{r} \frac{\partial u}{\partial r} - \frac{u}{r^{2}} + \frac{\partial^{2} u}{\partial z^{2}}] \\ + \frac{{(ρ β)}_{nf}}{ρ_{nf}} g (T - T_{\infty}) \cos ψ - \frac{σ_{nf} B_{0}^{2} u}{ρ_{nf}} - \frac{μ_{nf}}{ρ_{nf}} \frac{u}{K^{*}} - F^{*} u^{2}, \end{matrix}

(2)

\begin{matrix} u \frac{\partial v}{\partial r} + \frac{uv}{r} + w \frac{\partial v}{\partial z} = \frac{μ_{nf}}{ρ_{nf}} [\frac{\partial^{2} v}{\partial r^{2}} + \frac{1}{r} \frac{\partial v}{\partial r} - \frac{v}{r^{2}} + \frac{\partial^{2} v}{\partial z^{2}}] \\ - \frac{{(ρ β)}_{nf}}{ρ_{nf}} g (T - T_{\infty}) \sin ψ - \frac{σ_{nf} B_{0}^{2} v}{ρ_{nf}} - \frac{μ_{nf}}{ρ_{nf}} \frac{v}{K^{*}} - F^{*} v^{2}, \end{matrix}

(3)

\begin{matrix} u \frac{\partial w}{\partial r} + w \frac{\partial w}{\partial z} = - \frac{1}{ρ_{nf}} \frac{\partial p}{\partial z} + \frac{μ_{nf}}{ρ_{nf}} [\frac{\partial^{2} w}{\partial r^{2}} + \frac{1}{r} \frac{\partial w}{\partial r} + \frac{\partial^{2} w}{\partial z^{2}}] \\ - \frac{σ_{nf} B_{0}^{2} w}{ρ_{nf}} - \frac{μ_{nf}}{ρ_{nf}} \frac{w}{K^{*}} - F^{*} w^{2}, \end{matrix}

(4)

\begin{matrix} u \frac{\partial T}{\partial r} + w \frac{\partial T}{\partial z} = \frac{k_{nf}}{{(ρ C_{p})}_{nf}} [\frac{\partial^{2} T}{\partial r^{2}} + \frac{1}{r} \frac{\partial T}{\partial r} + \frac{\partial^{2} T}{\partial z^{2}}] \\ + 2 \frac{μ_{nf}}{{(ρ C_{p})}_{nf}} [{(\frac{\partial^{2} u}{\partial r^{2}})}^{2} + \frac{u^{2}}{r^{2}} + {(\frac{\partial^{2} w}{\partial z^{2}})}^{2}] \\ + \frac{σ_{nf} B_{0}^{2}}{{(ρ C_{p})}_{nf}} (u^{2} + v^{2}) + \frac{μ_{nf}}{{(ρ C_{p})}_{nf}} \\ [{(\frac{\partial^{2} v}{\partial z^{2}})}^{2} + {(\frac{\partial w}{\partial r} + \frac{\partial u}{\partial z})}^{2} + {r \frac{\partial}{\partial r} (\frac{v}{r})}^{2}] \\ - \frac{1}{{(ρ C_{p})}_{nf}} \frac{16 σ^{*} T_{\infty}^{3}}{3 κ^{*}} \frac{\partial^{2} T}{\partial z^{2}} + \frac{μ_{nf}}{{(ρ C_{p})}_{nf}} \frac{1}{K^{*}} (u^{2} + v^{2}) \\ + \frac{1}{{(ρ C_{p})}_{nf}} F^{*} (u^{3} + v^{3}) + \frac{1}{{(ρ C_{p})}_{nf}} q ″' . \end{matrix}

(5)

The conditions at the boundary for present problem are⁴²:

\begin{matrix} u = rS + γ^{*} {(1 - ξ^{*} \frac{\partial u}{\partial z})}^{- \frac{1}{2}} \frac{\partial u}{\partial z}, v = r Ω + {(1 - ξ^{*} \frac{\partial v}{\partial z})}^{- \frac{1}{2}} γ^{*} \frac{\partial v}{\partial z}, w = 0, \\ - k_{nf} (\frac{\partial T}{\partial z}) = h_{f} (T_{f} - T), at z = 0, \\ u \to 0, v \to 0, T \to T_{\infty}, p \to p_{\infty} as z \to \infty \end{matrix}

(6)

The following transformations are required for the current problem⁴³

\begin{matrix} u = r Ω F' (η), v = r Ω G (η), w = - \sqrt{2 Ω} υ_{SAE 50} F (η), \\ p = p_{\infty} - Ω μ_{f} P (η), \\ T = T_{\infty} + (T_{f} - T_{\infty}) θ (η), η = \sqrt{\frac{2 Ω}{υ_{f}} z} . \end{matrix}

(7)

Where $η$ is the non-dimensional distance along rotational axis.

$q ″'$ in non-dimensional for present situation can be expressed as

q ″' = \frac{k_{nf} Ω}{ν_{nf}} [A^{*} (T_{w} - T_{\infty}) F' + B^{*} (T - T_{\infty})],

(8)

where, $B^{*}$ and $A^{*}$ are temperature dependent and space-dependent coefficients. Additionally, $A^{*} > 0$ , $B^{*} > 0$ represents heat generation state and $A^{*} < 0$ , $B^{*} < 0$ , correspond to heat absorption of the system.

Equations (2)–(6) take the following form by using equations (7) and (8).

\begin{matrix} Ω_{2} (2 F ″' - KF') + 2 F F ″ - {F'}^{2} + G^{2} + (\frac{Gr}{R e^{2}}) Ω_{3} (\cos Ψ) θ \\ - (\frac{σ_{nf}}{σ_{f}}) Ω_{1} M F' - F_{r} {F'}^{2} = 0, \end{matrix}

(9)

\begin{matrix} Ω_{2} (2 G ″ - KG') + 2 F G' - F' G - (\frac{Gr}{R e^{2}}) Ω_{3} (\sin Ψ) θ \\ - (\frac{σ_{nf}}{σ_{f}}) Ω_{1} MG - F_{r} G^{2} = 0, \end{matrix}

(10)

\begin{matrix} (\frac{k_{nf}}{k_{f}} + \frac{4}{3} Rd) Ω_{4} + \Pr [F θ' + 6 (\frac{Ec}{Re}) Ω_{5} {F'}^{2} + Ec Ω_{5} ({F ″}^{2} + {G'}^{2}) \\ + Ec Ω_{5} K ({F'}^{2} + G^{2}) + Ec Ω_{4} F_{r} ({F'}^{3} + G^{3}) + \frac{1}{2} Ω_{4} \\ {(\frac{σ_{nf}}{σ_{f}}) MEc ({F'}^{2} + G^{2}) + \frac{1}{Ω_{2}} (\frac{k_{nf}}{k_{f}}) \frac{1}{\Pr} (A^{*} F' + B^{*} θ)}] = 0 . \end{matrix}

(11)

The boundary conditions are transformed as

\begin{matrix} F' = γ + γ_{1} {[1 - ξ F ″ (0)]}^{- \frac{1}{2}} F ″ (0), G = 1 \\ + γ_{1} {[1 - ξ G' (0)]}^{- \frac{1}{2}} G' (0), F = 0, \\ θ' = - (\frac{σ_{nf}}{σ_{f}}) B_{1} (1 - θ,) at η = 0 \\ F' \to 0, G \to 0, θ \to 0, as η \to 0, \end{matrix}

(12)

where

\begin{matrix} K = \frac{υ_{f}}{γ K^{*}}, Gr = \frac{g β_{f} (T_{f} - T_{\infty}) r^{3}}{υ_{f}^{2}}, Re = r (\frac{r Ω}{υ_{f}}), γ_{1} = γ^{*} \sqrt{\frac{2 Ω}{υ_{f}}}, \\ M = \frac{σ_{f} B_{0}^{2}}{Ω ρ_{f}}, F_{r} = \frac{C_{b}}{\sqrt{K^{*}}}, \Pr = \frac{υ_{f}}{α_{f}}, Rd = \frac{4 σ^{*} T_{\infty}^{3}}{k^{*} k_{f}}, Bi = \frac{h_{f}}{k_{f}} \sqrt{\frac{υ_{f}}{2 Ω}}, \\ Ec = \frac{r^{2} Ω^{2}}{{(C_{p})}_{f} (T_{f} - T_{\infty})}, γ = \frac{S}{Ω}, and ξ = ξ^{*} r Ω \sqrt{\frac{2 Ω}{υ_{f}}} \end{matrix}

We have made the following substitutions

\begin{matrix} Ω_{1} = \frac{ρ_{f}}{ρ_{nf}}, Ω_{2} = \frac{μ_{nf}}{μ_{f}} \frac{ρ_{f}}{ρ_{nf}}, Ω_{3} = \frac{{(ρ β)}_{nf}}{{(ρ β)}_{f}} \frac{ρ_{f}}{ρ_{nf}}, \\ Ω_{4} = \frac{{(ρ C_{p})}_{f}}{{(ρ C_{p})}_{nf}}, Ω_{5} = \frac{{(ρ C_{p})}_{f}}{{(ρ C_{p})}_{nf}} \frac{μ_{nf}}{μ_{f}} . \end{matrix}

Physical parameters

The engineering interest quantities such as circumferential wall shear stress $τ_{ϕ}$ and the radial wall stress $τ_{r}$ , can be expressed as follows with the help of Newtonian formula

At surface, the drag force is represented as

τ_{r} = μ_{f} {(\frac{\partial u}{\partial z} + \frac{\partial w}{\partial r})}_{z = 0} = μ_{f} (r Ω) \sqrt{\frac{2 Ω}{υ_{f}}} F ″ (0),

τ_{ϕ} = μ_{f} {(\frac{\partial v}{\partial z} + \frac{\partial w}{\partial r})}_{z = 0} = μ_{f} (r Ω) \sqrt{\frac{2 Ω}{υ_{f}}} G' (0) .

The skin friction coefficient is then given as

C_{f} = \frac{\sqrt{τ_{r}^{2} + τ_{ϕ}^{2}}}{ρ_{f} {(r Ω)}^{2}} .

The local skim friction in non-dimensional form is written as

\sqrt{R e_{r}} C_{f} = \sqrt{{[F ″ (0)]}^{2} + {[G' (0)]}^{2}} .

The required torque to rotate a disk of radius $r_{0}$ , can be determined as

T = - \int_{0}^{r_{0}} μ_{f} {(\frac{\partial u}{\partial z} + \frac{\partial w}{\partial r})}_{z = 0} 2 π r^{2} dr = - \frac{π ρ_{f} Ω r_{0}^{4}}{2} \sqrt{2 υ_{f} Ω} G' (0) .

The rate of heat transfer, can be determined as

N u_{r} = \frac{rq ″}{k_{f} (T_{w} - T_{\infty})},

where, $q ″ = - {(k_{f} + \frac{16 σ^{*} T_{\infty}^{3}}{3 k^{*}}) \frac{\partial T}{\partial z} |}_{z = 0}$ is wall heat flux. Equation (25) gives

N u_{r} = - 2 \sqrt{R e_{r}} (\frac{k_{nf}}{k_{f}} + \frac{4}{3} Rd) θ' (0) .

Solution methodology

The governing PDEs are transformed to ODEs with the introduction of appropriate similarity quantities. The resulting ODEs given in equations (9)−(11) along-with BCs given in equation (12), are solved numerically with MATLAB bvp-4c solver. We make use of following substitutions,

y_{1} = F, y_{2} = F', y_{3} = F ″, y_{3}' = F ″', y_{4} = G, y_{5} = G', y_{5}' = G ″,

y_{6} = θ, y_{7} = θ', y_{7}' = θ ″ .

The higher order equations are reduced to first order equations along with the initial conditions.

\begin{array}{l} y_{3^{'}} = \frac{K}{2} y_{2} - \frac{1}{2 Ω_{2}} \\ {2 y_{1} y_{3} - y_{2}^{2} + y_{4}^{2} + (\frac{G r}{R e^{2}}) Ω_{3} (c o s Ψ) y_{6} - Ω_{1} M (\frac{σ_{n f}}{σ_{f}}) y_{2} - F_{r} y_{2}^{2}} = 0, \end{array}

\begin{matrix} y_{5}' = \frac{K}{2} y_{4} - \frac{1}{2 Ω_{2}} \\ {2 y_{1} y_{5} - y_{2} y_{4} + (\frac{Gr}{R e^{2}}) Ω_{3} (\sin Ψ) y_{6} - Ω_{1} M (\frac{σ_{nf}}{σ_{f}}) y_{6} - F_{r} y_{4}^{2}} = 0, \end{matrix}

\begin{matrix} y_{7}' = - \frac{\Pr}{Ω_{4} (\frac{k_{nf}}{k_{f}} + \frac{4}{3} Rd)} [y_{1} y_{7} + 6 (\frac{Ec}{Re}) Ω_{5} y_{2}^{2} + Ec Ω_{5} (y_{3}^{2} + y_{5}^{2}) \\ + \frac{1}{2} Ω_{4} {MEc (\frac{σ_{nf}}{σ_{f}}) (y_{2}^{2} + y_{4}^{2}) + \frac{1}{Ω_{2} \Pr} (\frac{k_{nf}}{k_{f}}) (A^{*} y_{2} + B^{*} y_{6})} \\ + Ec Ω_{5} K (y_{2}^{2} + y_{4}^{2}) + Ec Ω_{4} F_{r} (y_{2}^{3} + y_{4}^{3})] . \end{matrix}

\begin{matrix} {\begin{matrix} y_{2} (0) = γ + γ_{1} {(1 - ξ y_{3} (0))}^{\frac{- 1}{2}} y_{3} (0), \\ y_{4} (0) = 1 + γ_{1} {(1 - ξ y_{5} (0))}^{\frac{- 1}{2}} y_{5} (0), \\ y_{1} (0) = 0, y_{7} (0) = - (\frac{k_{f}}{k_{nf}}) Bi (1 - y_{6} (0)), \\ y_{2} (\infty) = 0, y_{4} (\infty) = 0, y_{6} (\infty) = 0 . \end{matrix} \end{matrix}

Results and discussion

In order to make a strong understanding of the model, the solutions for velocity and temperature profiles are assembled numerically via MATLAB built-in bvp-4c solver and illustrated graphically. The relations for predicting thermo-physical properties of nanofluids are given in Table 1. The numerical values of thermo-physical properties of base lubricant (PAO) and nanoparticles $(A l_{2} O_{3})$ are given in Table 2. Figure 2 unveils the effects of inertial parameter $F_{r}$ on $F^{'} (η)$ , $G (η)$ , and $θ (η)$ . The increasing values of $F_{r}$ decline both the velocity profiles but favor the temperature of the nanolubricant $A l_{2} O_{3} - PAO$ . Physically, the boundary layer thickness depends upon the inertial parameter $F_{r}$ , which gets widen with an escalation in $F_{r}$ . Therefore, the movement of $A l_{2} O_{3} - PAO$ nanolubricant slows down along both the radial and tangential directions, whereas its temperature gets boosted. The behavior of radial velocity $F^{'} (η)$ and tangential velocity $G (η)$ of $A l_{2} O_{3} - PAO$ nanolubricant due to enhancing values of Grashof number $Gr$ is presented in Figure 3. The growing values of $Gr$ oppose the flow of nanolubricant along tangential direction but support the radial flow and hence reduction in tangential velocity $G (η)$ is resulted, while converse response occurs for radial velocity $F^{'} (η) .$ Figure 4 displays the impacts of velocity slip parameter $γ_{1}$ on tangential and radial velocities of $A l_{2} O_{3} - PAO$ nanolubricant. Both the velocity profiles of $A l_{2} O_{3} - PAO$ nanolubricant get decreased by increasing $γ_{1}$ . It is also perceived that when $γ_{1} \to \infty,$ the rotating disk stops rotating the nanolubricant. The flow becomes entirely potential for higher values of $γ_{1}$ due to which motion of the nanolubricant is resisted. The effects of magnetic parameter $M$ on $F^{'} (η)$ , $G (η)$ , and $θ (η)$ have been portrayed in Figure 5. The enhancing values of $M$ give rise to the Lorentz force which causes resistance to the flow of nanolubricant due to which the tangential and radial velocities of the nanolubricant decrease but its temperature goes up significantly. The imposed magnetic field plays a major role in adjusting and controlling movement of the nanolubricant over surface of the disk. Figures 6 and 7 unveil the influences of non-uniform parameters $A^{*}$ and $B^{*}$ respectively on temperature profile $θ (η)$ . When we enhance $A^{*}$ and $B^{*}$ positively they are known as heat sources and act as heat generators. They serve as releasing agents of heat energy into fluid flow and hence temperature of the nanolubricant gets improved. Additionally, when $A^{*}$ and $B^{*}$ are given negative values they are known to be heat sinks. They act as heat absorbers and cause a significant drop in temperature of the nanolubricant. Moreover, heat source or sink can be implemented as a physical tool to achieve desired cooling or heating rate in industrial sector and heat exchangers. The variation in $θ (η)$ corresponding to Radiation parameter $Rd$ is visualized in Figure 8. The thermal boundary gets enhanced by enhancing $Rd$ due to which temperature of the nanolubricant rises. Figure 9 exposes the effects of Biot number $Bi$ on temperature profile $θ (η) .$ The observation reveals that temperature of the nanolubricant falls by increasing $Bi$ . Physically, in present case the coefficient of heat transfer gets decreased by enhancing $Bi$ , which caused a decline in temperature of the nanolubricant. Table 3 presents the comparison of present computational results with the published work. It is pertinent to observe that present results excellently agree with the literature. Table 4 elaborates the effects of different parameters on Nusselt number. It is found that the growing values of radiation parameter and non-uniform heat generation parameters tend to increase value of Nusselt number. The value of Nusselt number drops with rising values of Biot number and non-uniform heat sink parameters.

Table 1.

Thermo-physical aspects of nanolubricant.

Properties	Nanofluid
Dynamic viscosity	$\frac{μ_{nf}}{μ_{f}} = {(1 - ϕ)}^{- 2.5}$
Density	$\frac{ρ_{nf}}{ρ_{f}} = {ϕ \frac{ρ_{s}}{ρ_{f}} + (1 - ϕ)}$
Heat capacitance	$\frac{{(ρ c_{p})}_{nf}}{{(ρ c_{p})}_{f}} = {ϕ \frac{{(ρ c_{p})}_{s}}{{(ρ c_{p})}_{f}} + (1 - ϕ)}$
Thermal expansion coefficient	$\frac{{(ρ β)}_{nf}}{{(ρ β)}_{f}} = {ϕ \frac{{(ρ β)}_{s}}{{(ρ β)}_{f}} + (1 - ϕ)}$
Electrical conductivity	$\frac{σ_{nf}}{σ_{f}} = (1 - ϕ) + ϕ \frac{σ_{s}}{σ_{f}}$
Thermal conductivity	$\frac{k_{nf}}{k_{f}} = \frac{((s - 1) k_{f} + k_{s}) - (s - 1) ϕ (k_{f} - k_{s})}{((s - 1) k_{f} + k_{s}) + ϕ (k_{f} - k_{s})}$

Table 2.

The features of $PAO$ (base fluid) and $A l_{2} O_{3}$ (nanoparticles).⁴⁴

S. no.		$A l_{2} O_{3}$	$PAO$
$1$	$ρ (Kg / m^{3})$	$3970$	$798$
$2$	$C_{p} (J / KgK)$	$765$	$2303$
$3$	$k (W / mK)$	$40$	$0.143$
$4$	$σ (S / m)$	$1.0 \times 10^{- 5}$	$1.0 \times 10^{- 11}$
5	$β (1 / K)$	$0.85 \times 10^{- 5}$	$3.5 \times 10^{- 4}$

Figure 2.

Effects of $F_{r}$ on $F' (η)$ , $G (η)$ , and $θ (η) .$

Figure 3.

Effects of $Gr$ on $F' (η)$ and $G (η) .$

Figure 4.

Effects of $γ_{1}$ on $F' (η)$ and $G (η) .$

Figure 5.

Effects of $M$ on $F' (η)$ , $G (η)$ , and $θ (η) .$

Figure 6.

Impact of $A^{*}$ on $θ (η) .$

Figure 7.

Impact of $B^{*}$ on $θ (η) .$

Figure 8.

Impact of $Rd$ on $θ (η) .$

Figure 9.

Impact of $Bi$ on $θ (η) .$

Table 3.

The comparison of present computational results of Skin friction with Nayak et al.,¹⁶ while additional assumptions are omitted.

$Γ$	Skin friction
	$M = 1.0$		$M = 3.0$		$M = 5.0$
	Nayak et al.¹⁶	Present work	Nayak et al.¹⁶	Present work	Nayak et al.¹⁶	Present work
$0.1$	$1.46870$	$1.468704$	$1.89639$	$1.896392$	$2.21627$	$2.216271$
$0.5$	$0.94252$	$0.94252$ 3	$1.12955$	$1.12955$ 4	$1.24770$	$1.24770$ 2
$0.8$	$0.74860$	$0.748601$	$0.87208$	$0.872083$	$0.94436$	$0.944360$

Table 4.

Effects of different parameters on Nusselt number.

$\Pr$	$Ec$	$Rd$	$A^{*}$	$B^{*}$	$Bi$	$N u_{r}$
$6.1$						$1.1318$
$6.3$	$0.4$	$2.2$	$0.4$	$0.3$	$0.8$	$1.0296$
$6.5$						$0.9441$
	$0.3$					$1.2691$
$6.3$	$0.5$	$2.2$	$0.4$	$0.3$	$0.8$	$2.0559$
	$0.7$					$3.3642$
		$2.1$				$0.9601$
$6.3$	$0.3$	$2.3$	$0.4$	$0.3$	$0.8$	$1.1633$
		$2.5$				$1.3788$
			$1.0$			$1.4271$
$6.3$	$0.3$	$2.2$	$2.0$	$0.3$	$0.8$	$1.5381$
			$3.0$			$1.6128$
				$- 0.1$		$0.8741$
$6.3$	$0.3$	$2.2$	$0.4$	$- 0.5$	$0.8$	$0.8221$
				$- 0.9$		$0.7743$
					$0.4$	$0.8616$
$6.3$	$0.3$	$2.2$	$0.4$	$0.3$	$0.6$	$0.8348$
					$0.8$	$0.8119$

Concluding remarks

In this work, an attempt is made to develop a mathematical model to discuss the flow of magnetized $A l_{2} O_{3} - PAO$ nanolubricant over an inclined rotating disk with Thompson and Troian slip at the boundary. The influences of mixed convection, nonlinear heat radiation, viscous dissipation, Joule heating, and non-uniform heat sink/source are also included in the modeling. The flow is supposed to take place in a Darcy-Forchheimer porous medium. We have solved the proposed model numerically with MATLAB built-in bvp-4c method after metamorphosing the PDEs to ODEs. The key findings of this work are as follows:

The enhancing values of velocity slip parameter decrease the tangential and radial velocities of the nanolubricant $A l_{2} O_{3} - PAO .$

The growing values of Grashof number oppose the flow of nanolubricant along tangential direction but support the radial flow.

With an escalation in inertial parameter, the movement of $A l_{2} O_{3} - PAO$ nanolubricant slows down along both the radial and tangential directions, whereas its temperature gets boosted.

The temperature of the nanolubricant $A l_{2} O_{3} - PAO$ enhances significantly by strengthening the magnetic field, whereas radial and tangential velocities get retarded.

The non-uniform heat source/sink parameters play a vital role in controlling heat transmission phenomenon.

The temperature of the nanolubricant $A l_{2} O_{3} - PAO$ falls by increasing Biot number and converse behavior is observed for radiation parameter.

The growing values of Eckert number, radiation parameter, and non-uniform heat source parameters tend to increase the value of Nusselt number.

The value of Nusselt number drops with rising values of Biot number and non-uniform heat sink parameters.

Footnotes

Appendix

Handling Editor: Chenhui Liang

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) received no financial support for the research, authorship, and/or publication of this article.

ORCID iD

Nargis Khan

References

Martin

Ohmae

. Nanolubricants. Chichester: John Wiley & Sons, 2008.

Rapoport

Moshkovich

Perfilyev

, et al. On the efficacy of IF–WS2 nanoparticles as solid lubricant: the effect of the loading scheme. Tribol Lett 2007; 28: 81–87.

Joly-Pottuz

Iwaki

. Superlubricity of tungsten disulfide coatings in ultra high vacuum. In: Erdemir

Martin

J-M

(eds) Superlubricity. Amsterdam: Elsevier Science BV, 2007, pp.227–236.

Bordignon

Salvaro

Binder

, et al. Tribological behaviour of plasma-functionalized graphene as low-viscosity oil additive. Tribol Lett 2018; 66: 1–14.

Wang

Liu

. Nanoparticle-based lubricant additives. In: Wang

Chung

(eds) Encyclopedia of tribology. New York, NY: Springer, 2013, pp.2369–2376.

Borda

FLG

de Oliveira

SJR

Lazaro

LMSM

, et al. Experimental investigation of the tribological behavior of lubricants with additive containing copper nanoparticles. Tribol Int 2018; 117: 52–58.

Zuin

Cousseau

Sinatora

, et al. Lipophilic magnetite nanoparticles coated with stearic acid: a potential agent for friction and wear reduction. Tribol Int 2017; 112: 10–19.

Ghalme

Koinkar

Bhalerao

. Effect of aluminium oxide (Al₂O₃) nanoparticles addition into lubricating oil on tribological performance. Tribol Ind 2020; 42: 494–502.

Guo

Liu

, et al. Performance of a domestic refrigerator using TiO₂-R600a nano-refrigerant as working fluid. Energy Convers Manag 2011; 52: 733–737.

10.

Pico

DFM

da Silva

LRR

Mendoza

OSH

, et al. Experimental study on thermal and tribological performance of diamond nanolubricants applied to a refrigeration system using R32. Int J Heat Mass Transf 2020; 152: 119493.

11.

Choi

Kim

Jang

, et al. Effect of polyolester oil-based multiwalled carbon-nanotube nanolubricant on the coefficient of performance of refrigeration systems. Appl Therm Eng 2021; 192: 116941.

12.

Gamaoun

Said

Makki

, et al. Energy transfer of a fin wetted with ZnO-SAE 50 nanolubricant: an application of α-parameterized differential transform method. Case Stud Therm Eng 2022; 40: 102501.

13.

Riaz

Khan

Shehzad

. Rheological behavior of magnetized ZnO – SAE50 nanolubricant over Riga plate: a theoretical study. Adv Mech Eng 2023; 15: 16878132231162305.

14.

Nayak

Mehmood

Muhammad

, et al. Entropy minimization in mixed convective Falkner-Skan flow of ZnO-SAE50 nanolubricant over stationary/moving Riga plate. Case Stud Therm Eng 2021; 26: 101176.

15.

Nayak

Prakash

Tripathi

, et al. 3D radiative convective flow of ZnO-SAE50 nano-lubricant in presence of varying magnetic field and heterogeneous reactions. Propulsion Power Res 2019; 8: 339–350.

16.

Nayak

Mehmood

Makinde

, et al. Magnetohydrodynamic flow and heat transfer impact on ZnO-SAE50 nanolubricant flow over an inclined rotating disk. J Cent South Univ 2019; 26: 1146–1160.

17.

Khan

. Bioconvection in second grade nanofluid flow containing nanoparticles and gyrotactic microorganisms. Braz J Phys 2018; 48: 227–241.

18.

Zuhra

Khan

Shah

, et al. Simulation of bioconvection in the suspension of second grade nanofluid containing nanoparticles and gyrotactic microorganisms. AIP Adv 2018; 8: 105210.

19.

Mehryan

SAM

Sheremet

Soltani

, et al. Natural convection of magnetic hybrid nanofluid inside a double-porous medium using two-eq28 energy model. J Mol Liq 2019; 277: 959–970.

20.

Yan

Izadi

Sheremet

, et al. Inclined Lorentz force impact on convective-radiative heat exchange of micropolar nanofluid inside a porous enclosure with tilted elliptical heater. Int Commun Heat Mass Transf 2020; 117: 104762.

21.

Abbasi

Shehzad

. Magnetized peristaltic transportation of boron-nitride and ethylene-glycol nanofluid through a curved channel. Chem Phys Lett 2022; 803: 139860.

22.

Lanjwani

Chandio

Anwar

, et al. MHD laminar boundary layer flow of radiative Fe-Casson nanofluid: stability analysis of dual solutions. Chin J Phys 2022; 76: 172–186.

23.

Sarada

Gowda

RJP

Sarris

, et al. Effect of magnetohydrodynamics on heat transfer behaviour of a non-Newtonian fluid flow over a stretching sheet under local thermal non-equilibrium condition. Fluids 2021; 6: 264.

24.

Punith Gowda

Sarris

Naveen Kumar

, et al. A three-dimensional non-Newtonian magnetic fluid flow induced due to stretching of the flat surface with chemical reaction. J Heat Transf 2022; 144: 113602.

25.

Kumar

Gamaoun

Abdulrahman

, et al. Heat transfer analysis in three-dimensional unsteady magnetic fluid flow of water-based ternary hybrid nanofluid conveying three various shaped nanoparticles: a comparative study. Int J Mod Phys B 2022; 36: 2250170.

26.

Khan

Gul

Islam

, et al. Brownian motion and thermophoresis effects on MHD mixed convective thin film second-grade nanofluid flow with Hall effect and heat transfer past a stretching sheet. J Nanofluids 2017; 6: 812–829.

27.

Zuhra

Khan

Islam

. Magnetohydrodynamic second-grade nanofluid flow containing nanoparticles and gyrotactic microorganisms. Comput Appl Math 2018; 37: 6332–6358.

28.

Khan

Islam

Gul

, et al. Thin film flow of a second grade fluid in a porous medium past a stretching sheet with heat transfer. Alex Eng J 2018; 57: 1019–1031.

29.

Hsiao

. Viscoelastic fluid over a stretching sheet with electromagnetic effects and nonuniform heat source/sink. Math Probl Eng 2010; 2010: 1–14.

30.

Ramesh

Chamkha

Gireesha

. MHD mixed convection flow of a viscoelastic fluid over an inclined surface with a nonuniform heat source/sink. Can J Phys 2013; 91: 1074–1080.

31.

Ramandevi

Reddy

JVR

Sugunamma

, et al. Combined influence of viscous dissipation and non-uniform heat source/sink on MHD non-Newtonian fluid flow with Cattaneo-Christov heat flux. Alex Eng J 2018; 57: 1009–1018.

32.

Sulochana

Payad

Sandeep

. Non-uniform heat source or sink effect on the flow of 3D Casson fluid in the presence of Soret and thermal radiation. Int J Eng Res Africa 2016; 20: 112–129.

33.

Palwasha

Khan

Shah

, et al. Study of two-dimensional boundary layer thin film fluid flow with variable thermo-physical properties in three dimensions space. AIP Adv 2018; 8: 105318.

34.

Khan

Gul

Islam

, et al. Thermophoresis and thermal radiation with heat and mass transfer in a magnetohydrodynamic thin-film second-grade fluid of variable properties past a stretching sheet. Eur Phys J Plus 2017; 132: 20.

35.

Khan

Hayat

Khan

, et al. Activation energy impact in nonlinear radiative stagnation point flow of Cross nanofluid. Int Commun Heat Mass Transf 2018; 91: 216–224.

36.

Hayat

Khan

Qayyum

, et al. New thermodynamics of entropy generation minimization with nonlinear thermal radiation and nanomaterials. Phys Lett A 2018; 382: 749–760.

37.

Ghadikolaei

Hosseinzadeh

Ganji

, et al. Nonlinear thermal radiation effect on magneto Casson nanofluid flow with Joule heating effect over an inclined porous stretching sheet. Case Stud Therm Eng 2018; 12: 176–187.

38.

Waqas

Jabeen

Hayat

, et al. Modeling and analysis for magnetic dipole impact in nonlinear thermally radiating Carreau nanofluid flow subject to heat generation. J Magn Magn Mater 2019; 485: 197–204.

39.

Farooq

Ijaz Khan

Waqas

, et al. Transport of hybrid type nanomaterials in peristaltic activity of viscous fluid considering nonlinear radiation, entropy optimization and slip effects. Comput Methods Programs Biomed 2020; 184: 105086.

40.

Mosayebidorcheh

Hatami

. Heat transfer analysis in carbon nanotube-water between rotating disks under thermal radiation conditions. J Mol Liq 2017; 240: 258–267.

41.

Usman

Hamid

Haq

, et al. Heat and fluid flow of water and ethylene-glycol based Cu-nanoparticles between two parallel squeezing porous disks: LSGM approach. Int J Heat Mass Transf 2018; 123: 888–895.

42.

Ramzan

Chung

Kadry

, et al. Nanofluid flow containing carbon nanotubes with quartic autocatalytic chemical reaction and Thompson and Troian slip at the boundary. Sci Rep 2020; 10: 18710.

43.

Mustafa

Khan

. Numerical study of partial slip effects on MHD flow of nanofluids near a convectively heated stretchable rotating disk. J Mol Liq 2017; 234: 287–295.

44.

Hatami

Farooq

ALİ

Alsabery

, et al. Mixed convection heat transfer of SiO₂-water and alumina-PAO nano-lubricants used in a mechanical ball bearing. J Therm Eng 2021; 7: 134–161.