Next Article in Journal
Special Issue Editorial: “Discrete and Continuous Memristive Nonlinear Systems and Symmetry”
Next Article in Special Issue
Experimental Study on the Influence of Interfacial Energy Instability on the Flow Pattern Spatiotemporal Evolution of Thermal- Buoyant Capillary Convection
Previous Article in Journal
Soliton Waves with the (3+1)-Dimensional Kadomtsev–Petviashvili–Boussinesq Equation in Water Wave Dynamics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analytical Study of the Energy Loss Reduction during Three-Dimensional Engine Oil-Based Hybrid Nanofluid Flow by Using Cattaneo–Christov Model

by
Ramadan A. ZeinEldin
1,
Asad Ullah
2,*,
Hamiden Abd El-Wahed Khalifa
3,4 and
Muhammad Ayaz
5
1
Deanship of Scientific Research, King Abdulaziz University, Jeddah 21589, Saudi Arabia
2
Department of Mathematical Sciences, University of Lakki Marwat, Lakki Marwat 28420, Khyber Pakhtunkhwa, Pakistan
3
Department of Operations and Management Research, Faculty of Graduate Studies for Statistical Research, Cairo University, Giza 12613, Egypt
4
Department of Mathematics, College of Science and Arts, Qassim University, Al-Badaya 51951, Saudi Arabia
5
Department of Mathematics, Abdul Wali Khan University Mardan, Mardan 23200, Khyber Pakhtunkhwa, Pakistan
*
Author to whom correspondence should be addressed.
Symmetry 2023, 15(1), 166; https://doi.org/10.3390/sym15010166
Submission received: 8 December 2022 / Revised: 27 December 2022 / Accepted: 4 January 2023 / Published: 6 January 2023
(This article belongs to the Special Issue Fluid Flow and Heat Transfer, Symmetry and Asymmetry)

Abstract

:
In this work, we analyzed the hybrid nanofluid (Ag+CuO+kerosene oil) flow past a bidirectionally extendable surface in the presence of a variable magnetic field. The hybrid nanofluid flow considered is electrically conductive and steady. For the simulation of the problem, the Cattaneo–Christov double-diffusion (CCDD) model was considered, which generalizes Fourier’s and Fick’s laws. The impact of the Hall current produced was taken into account. The physical problem was transformed into a mathematical form with the help of suitable transformations to reduce the complexity of the problem. The transformed system of coupled ordinary differential equations (ODEs) was solved with the semi-analytical method. The results are plotted in comparison with the ordinary nanofluid (CuO+kerosene oil) and hybrid nanofluid (Ag+CuO+kerosene oil). The impact of various parameters ( P r , S c , γ 0 , m , M , N b , N t , ϵ 1 , ϵ 2 ) on the state variables is described. The velocity gradient under the impact of the mass flux and magnetic parameter shows a decreasing behavior, while the Hall parameter and the stretching ratio show an increasing behavior. Moreover, the skin friction, rate of heat, and mass transfer are numerically displayed through tables. In this work, we found that the thermal and concentration relaxation coefficients showed a decreasing behavior for their increasing trends. For the validation of the implemented technique, the squared residuals are computed in Table 2, which shows that the increasing number of iterations decreases the squared residual error. The results show that Ag+CuO+kerosene oil has good performance in the reduction of the heat transfer rate.

1. Introduction

In the last few years, the applications of heat transfer in different fields (engineering and technology) have attracted many researchers. One such important application is the nature of the electronic devices used, where the thermal conductivity is very low, due to the ordinary mixtures (ethylene glycol, water, oil, etc.). However, now, almost every researcher has access to increase the thermal conductivity of these fluids (materials). This is thanks to Choi and Eastman [1,2], who introduced the idea by the mixing of small-sized particles (having a size below 100 nm) with the available fluid. This available fluid is known as the base fluid, while the mixture formed by a single-type mixing procedure is known as the nanofluid. This new mixture formed has many physical advantages over the ordinary fluid. This procedure opened the door to a new era in the field of fluid mechanics. Initially, this work looked to be not very interesting, but as time passed, this work has put the interest of the researchers on a single track. Currently, two different techniques are used for the preparation of these nanofluids. These fluids can be prepared at the laboratory, as well as at the industry level. Maiga et al. [3] conducted a numerical simulation of the nanofluid’s behavior. They estimated the heat transfer increase of the nanofluid in an evenly heated tube under turbulent and laminar flow conditions and employed approximated correlations for the experimental data. They reported that the increasing value of the Reynolds number in the turbulent case increases the transfer of heat. This rapid increase was found to be due to the dimension of the nanoparticles used in the formation of the nanofluid. A new theoretical model for nanoliquids’ behavior through an expanding surface was analyzed by Khan and Pop [4]. In this work, they reported the stabilization of the stream and thermal behavior. The magnetohydrodynamics (MHD) mixed convective stagnation point flow of a micro-polar nanoliquid was performed by Seth et al. [5]. They used a new empirical relation for this study and implemented the successive linearization method for solution purposes. The study of nanofluids subject to different conditions and geometries can be found in [6,7,8].
Researchers are always trying to investigate problems and find a feasible and suitable solution. In light of the above literature, researchers have taken keen interest in advancing the existing work to a next level. In the preparation of nanofluids instead of single-element nanoparticles, two or more than two nanoparticles of different elements are used. Therefore, a single-type nanoparticle is replaced with a hybrid one. These nanofluids have various applications in agriculture production, applied and biological sciences, and material sciences. These are accessible with a very low cost and have many applications in the transfer of heat and other coolant projects. On the other hand, the advantages of these fluids include increased thermal conductivity, good transfer of heat, and other physical significance. Due to these properties, they are used for coolant purposes to increase efficiency and reduce energy loss. In 2007, for the first time, Jana et al. [9] investigated the thermal conductivity of hybrid nanofluids. They performed a comparative survey with a single-phase common nanofluid, where they found a remarkable increase in the transfer of heat for computing the thermal conductivity of the hybrid nanofluid. Madhesh et al. [10] studied a copper-titanium hybrid nanofluid to analyze the rheological characteristic, as well as the transfer of heat. In this work, in the preparation of the hybrid nanofluid, 55nm particles with concentrations varying from 0.1 vol.% to 2.0 vol.% were used. In the experiment, a 52% increase in the coefficient of the convective heat, a 49% increase in the Nusselt number, an average increase of 68% in the coefficient of heat transfer, and a 1.0 % increase in the volume concentration were observed. Later on, in 2010, Ho et al. [11] performed an investigation of the increase of the thermal conductivity by utilizing a water-based hybrid nanofluid having a microencapsulated phase change material (MEPCM) and Al2O3 nanoparticles. In this work, higher thermal conductivity and concentration as compared to the phase change materials (PCMs) were reported. Suresh et al. [12] studied the nanoparticles’ concentration effect over the thermal conductivity for the hybrid nanofluid Al2O3-Cu/water. In this study, for synthesis purposes, the two-step method was implemented. Labib et al. [13] considered the Al2O3–CNT/water nanofluid and studied the impact of the Reynolds number and the nanoparticles’ concentration for the transfer of heat. They reported that the increase of Al2O3 and carbon nanotubes (CNTs) inside the medium increased the thermal conductivity. In this work, the concentration of the nanoparticles increased the transfer of heat coefficient and work, in addition to the friction factor. The Al2O3-Cu/water hybrid nanofluid flow by considering the forced convection was utilized by Moghadassi et al. [14]. Both single- and two-phase analysis were performed for synthesis purposes, and both phases were analyzed. They reported a 4.73 % increase in the single-phase model, while 13.46 % in the two-phase model for the hybrid nanofluid in comparison with water. Furthermore, a 4% increase in the Nusselt number for the two-phase model was observed.
In recent years, various uses for electrically conducting fluid and MHD boundary flow have been reported in engineering and industry [15]. Similar applications of the MoS2 nanoparticles for the MHD slip flow were reported by Khan [16]. His results showed the heat transfer increase for the MoS2 nanoparticles. Chamkaa et al. [17] analyzed the Cu-water nanofluid filled square cavity and studied the mixed convection with the impact of the magnetic field. The mathematical problem was tackled with the finite-difference method. The results were presented for the inclination angle, the Hartmann and Richardson parameters, and the nanofluid volume fraction under the impact of the applied magnetic field. Alsabery et al. [18] considered the two-phase Buongiorno model for the square enclosure in the mixed convection of Al2O3-water nanofluid flow with the impact of the magnetic field. Various obstacles (circular, triangular, and square) were assumed inside the cavity. The Galerkin method was employed for the solution of the problem, and the results were presented in the form isotherms and streamlines. They found that the heat produced inside the cavity was higher for the triangular obstacles compared to the others. Ellahi [19] showed that the velocity profile decreased with the variation in the MHD parameter for the non-Newtonian MHD fluid flow inside a pipe. The pipe was considered at a higher temperature compared to the fluid flowing through it. Explicit expressions for the state variables were derived with the semi-analytical method, the HAM. Sheikholeslami et al. [20] reported the MHD fluid flow inside a permeable channel. They found that, with the increasing values of the Hartmann number, the velocity boundary layer declined, while on the other hand, the boundary layer thickness decreased with the increasing values of the Reynolds number and volume fraction. Seyyedi et al. [21] used the ANSYS Fluent software and the control volume finite-element method for the natural convection MHD fluid flow through a semi-annulus enclosure. In their study, they reported the entropy generation due to the MHD flow. In this work, they showed that the entropy generation decreased for various pertinent parameters. Rashidi et al. [22] studied the entropy generation and the second law of thermodynamics for the MHD nanofluid flow inside a porous rotating disk. In this study, the applications of the rotating disk with entropy generation were presented. In the investigation of the MHD (three-dimensional) 3D flow for the Lorentz force impact over the temperature profile, Sheikholeslami and Ellahi used the lattice Boltzmann method (LBM) [23]. The viscosity and thermal conductivity were calculated with the correlation model, Koo–Kleinstreuer–Li (KKL). The generation of the entropy and heat transfer both were reported to decline with the increasing values of the Rayleigh parameter, while the opposite behavior was reported for the Hartmann parameter. Recently, Biswas et al. [24] studied the application of Cu-water nanofluid inside a square cavity, where the left wall was considered hot. The problem was solved with the finite-volume method, and various parameters’ impacts were studied. The wavy wall and convection impacts were studied in detail. The bio-convection impacts over the heat and mass transfer were reported for cleaner energy and bio-cell production. Manna et al. [25] reported the entropy generation and Lorentz forces’ impact inside a quarterly squared cavity. The work was carried out for various dimensionless parameters, where the results described with the practical applications of energy storage, the cooling of electronic devices, and other food and technological purposes. Biswas et al. [26] also reported partial magnetic impacts on the nanoliquid flow past a porous medium numerically. The magnetic field impact showed a 17.15 % or 9.71 % decline in the heat transfer rate. More recent studies of the flow inside a square cavity can be found in [27,28,29,30].
The boundary layer flow over different geometries and, especially, in the case of a stretching sheet has been analyzed by various researchers. The role of this type of flow is encountered in wire drawing, copper wire tinning, metallurgy, glass fiber, rubber and plastic sheets, glass blowing, etc. The basic idea of the boundary layer flow was presented by Sakiadis [31,32,33] in the case of a solid surface motion at a fixed speed. This idea was further utilized by Crane [34] for a two-dimensional (2D) flow past a stretching sheet. The scientist continued their work on this new idea and extended the work further to other physical constraints [35,36]. The boundary layer flow past an exponentially stretchable sheet by considering the suction and injection of the sheet for the heat and mass transfer analysis was utilized by Magyari et al. [37]. In the boundary layer flow analysis, the shrinking and stretching sheet caused the flow to move back and forth towards the origin. This unusual behavior was reported by Wang [38] for the first time by analyzing the liquid film’s unsteady flow past a stretching sheet. The exact solution of the Navier–Stokes equations for a stretching sheet was presented by Miklavcic and Wang [39]. They also reported that the stability of the flow over the sheet requires the suction of mass. Bhattacharyya [40] investigated fluid flow and heat transfer across an exponentially shrinking sheet. Their findings showed that the feasibility of the steady flow depends on the mass suction parameter when exceeding a specific critical value. In another study, Bhattacharyya et al. [41] examined the stagnation point flow and heat transfer across a sheet subject to exponential shrinkage. Bachok et al. [42] examined the steady two-dimensional stagnation point flow of an aqueous nanofluid over an exponentially stretching/shrinkage surface. A substantial amount of research has been dedicated to the analysis of fluid flow across the exponentially stretching/shrinking plate. Some recent examples are mentioned here [43,44,45,46].
The Buongiorno nanoliquid model is modified with the help of the Oldroyd convection time derivative to obtain the Cattaneo–Christov model [47]. This modification is mainly concerned with the Fourier–Fick laws [48]. Zheng and Zhang [49] investigated the Maxwell nanofluid flow past a stretching sheet by considering the heat and mass transfer together with the slip effects. They used the Cattaneo–Christov double-diffusion (CCDD) model in their analysis. Hayat et al. [50] studied the steady, 3D, and bidirectional Newtonian fluid flow by considering the CCDD model. In the very next year, Hayat et al. [51] analyzed analytically the viscoelastic nanofluid flow by considering the CCDD model. Recently, Hayat et al. [52,53] investigated the CCDD model for the Prandtl and Walters-B fluids by considering the heat source/sink and chemical reaction effects past an expanding surface. The Eyring–Powell nanomaterial was analyzed by Ibrahim [54] considering the CCDD model for the flow over a bidirectionally expanding sheet. More recent surveys of the Cattaneo–Christov model can be found in [55,56,57,58]. From the literature, it is clear that the impact of the Hall current on the CCDD model with a variable magnetic field has not been reported for the hybrid nanofluid Ag+CuO+kerosene oil.
This work aimed to analytically investigate the hybrid nanofluid flow based on kerosene oil past a bidirectionally expanding surface in the presence of a variable magnetic field. The nanofluid flow considered is electrically conductive and steady. For the exploration and modeling of the physical problem, the CCDD model was used for the optimal solution. Moreover, the current produced due to the magnetic parameter was also counted in the physical problem. In this article, Section 2 is devoted to the problem description and physical geometry with a detailed mathematical description. In Section 3, the analytical method used is described in detail with residual error analysis, and in Section 4, the results obtained are described and discussed in detail. Finally, the discussion of the tables and conclusions are presented in Section 5 and Section 6, respectively.

2. Mathematical Model of the Problem

Consider a steady, bidirectional and incompressible flow of a hybrid nanofluid (Ag+CuO+kerosene oil) over an extendable surface. The geometry of the problem was chosen in the Cartesian coordinate system, centered at ( x , y , z ) = ( 0 , 0 , 0 ) , as can be seen in Figure 1[50]. The geometry of the problem was chosen in such away that the u w = c exp x + y L and v w = d exp x + y L velocities act along the x- and y-axis, respectively. The wall temperature and concentration of the hybrid nanoparticles were assumed in the exponential forms T w T = T 0 exp a ( x + y ) 2 L and C w C = C 0 exp a ( x + y ) 2 L , respectively. A variable magnetic field B ( x , y ) = B 0 exp x + y 2 L was applied perpendicular to the surface, where B 0 shows the strength of the applied magnetic field. The fluid is electrically conductive, in which the induced current produced was considered in the nanofluid fluid flow. Moreover, the heat and mass transfer were investigated through the double-diffusion Cattaneo–Christov theory, which acts as a generalization of the famous Fourier and Fick laws [47,50]:
q ˜ + Λ e q ˜ t + V ˜ · q ˜ q ˜ · V ˜ + ( · V ˜ ) q ˜ + k T = 0 ,
J ˜ + Λ c J ˜ t + V ˜ · J ˜ J ˜ · V ˜ + ( · V ˜ ) J ˜ + D B C = 0 .
In Equations (1) and (2), q ˜ and J ˜ show the heat and mass fluxes, T and C denote the temperature and concentration, respectively, k is the thermal conductivity, V ˜ is the velocity field, D B is the Brownian motion parameter, and Λ e and Λ c describe the energy and concentration relaxation factors, respectively. By letting Λ e = 0 = Λ c , the famous Fourier and Fick laws can be obtained. Considering the steadiness condition of the nanofluid flow, we have from Equations (1) and (2) the theory, which acts as a generalization of the famous Fourier and Fick laws [47,50].
q ˜ + Λ e V ˜ · q ˜ q ˜ · V ˜ + ( · V ˜ ) q ˜ + k T = 0 ,
J ˜ + Λ c V ˜ · J ˜ J ˜ · V ˜ + ( · V ˜ ) J ˜ + D B C = 0 .
Since we considered the nanofluid flow as electrically conductive under the impact of the applied variable magnetic parameter, after some time, the magnetic field produces a current known as the Hall current. This whole process is governed by the generalized Ohm’s law given below [59,60,61]:
j ˜ + ω e B 0 t e × ( j ˜ × B ˜ ) + σ P e e n e = σ ( E ˜ + V ˜ × B ˜ ) ,
where n e , e, σ , j ˜ , V ˜ , and B 0 describe the number density, charge of electrons, electrical conductivity and electrical current density, velocity of the fluid, and strength of the magnetic field, respectively. In light of the above assumptions, the modified Buongiorno nanofluid model for the boundary layer nanofluid flow takes the following form [52,55,62]:
u x + v y + w z = 0 ,
u u x + v u y + w u z = ν h n f u z z σ h n f ρ h n f B 2 1 + m 2 ( u m v ) ,
u v x + v v y + w v z = ν h n f v z z σ h n f ρ h n f B 2 1 + m 2 ( v + m u ) ,
u T x + v T y + w T z = α m T z z Φ e Λ e + τ T z C z + D T T T z 2 ,
u C x + v C y + w C z = D B C z z Φ c Λ c + D T T C z z .
Here, in Equations (9) and (10), the terms Φ e and Φ c are the steady Cattaneo–Christov expressions defined as [59]:
Φ e = u 2 T x x + v 2 T y y + w 2 T z z + 2 ( u v T x y + v w T y z + u w T x z ) + T x ( u u x + v u y + w u z ) + T y ( u v x + v v y + w v z ) + T z ( u w x + v w y + w w z ) ,
Φ c = u 2 C x x + v 2 C y y + w 2 C z z + 2 ( u v C x y + v w C y z + u w C x z ) + C x ( u u x + v u y + w u z ) + C y ( u v x + v v y + w v z ) + C z ( u w x + v w y + w w z ) .
Here, u, v, and w are the velocity components along the x, y, and z directions, respectively, ρ is the density, ν is the kinematic viscosity, Λ e and Λ c represent the relaxation time of the energy and concentration fluxes, τ is the adiabatic index, D T represents the thermophoretic parameter, and D B indicates the diffusion parameter of the nanoparticles. The boundary restrictions from the physical problem take the following form.
u = u w ( x , y ) , v = v w ( x , y ) , w = w w ( x , y ) = w 0 exp x + y 2 L , T = T w , C = C w , as z 0 , u 0 , v 0 , C C , T T as z .
where u w = c exp x + y L and v w = d exp x + y L are the bidirectional velocities along the x- and y a x i s , respectively. Here, c and d are the references of the velocities and L is the length (characteristic) of the sheet. Introduce the following similarity variable [62]:
u = c exp x + y 2 L f , v = c exp x + y 2 L g , w = c ν f 2 L 0.5 exp x + y 2 L ( f + g + ζ f + ζ g ) , T = T + T 0 exp a ( x + y ) 2 L , C = C + C 0 exp a ( x + y ) 2 L , ζ = c 2 ν f L 0.5 exp x + y 2 L , T T = ( T f T ) θ ( ζ ) .
Here, the prime shows the derivative with respect to ζ and w 0 = c ν f 2 L 0.5 S 0 , where S 0 represents the dimensionless flux of the mass, which further explains the suction/injection.
Equations (6)–(10) together with the boundary conditions (13), by implementing the conditions presented in Equation (14), can be written in the following form [50,62]:
μ h n f μ f ρ h n f ρ f f + f ( g + f ) 2 f 2 2 f g σ h n f σ f ρ h n f ρ f M 1 + m 2 ( f m g ) = 0 ,
μ h n f μ f ρ h n f ρ f g + g ( g + f ) 2 g 2 2 f g σ h n f σ f σ h n f σ f M 1 + m 2 ( g + m f ) = 0 ,
k h n f k f ( ρ C p ) h n f ( ρ C p ) f θ P r + θ ( f + g ) a ( f + g ) θ + ϵ 1 [ ( ζ ( f + g ) + ( 2 a + 1 ) ( f + g ) ) ( f + g ) θ a ( ( a + 2 ) ( f + g ) 2 ( f + g ) ( f + g ) ) θ ( f + g ) 2 θ ] + N b θ ϕ + N t θ 2 = 0 ,
ϕ S c + N t N b θ S c + ϕ ( f + g ) a ( f + g ) ϕ + ϵ 2 [ ( ζ ( f + g ) + ( 2 a + 1 ) ( f + g ) ) ( f + g ) ϕ a ( ( a + 2 ) ( f + g ) 2 ( f + g ) ( f + g ) ) ϕ ( f + g ) 2 ϕ ] = 0 ,
With the B.Cs.:
f = 0 = g , f = 1 , g = γ 0 , θ = 1 , ϕ = 1 at ζ = 0 , f 0 , g , θ 0 , ϕ 0 as ζ .
Here, P r = ν f α f is the Prandtl number, N t = τ D T ( T w T ) ν T and N b = τ D B ( C w C ) ν are the thermophoretic and Brownian motion parameters, γ 0 = d c is the stretching ratio, S c = ν D B is the Schmidt number, ϵ 1 = Λ e a and ϵ 2 = Λ e a are the thermal relaxation parameters for the energy and mass fluxes, respectively, M = 2 L σ f B 0 2 ρ C p is the magnetic parameter, m is the Hall parameter, and μ h n f , ρ h n f , and ( ρ C p ) h n f show the hybrid nanofluid dynamic viscosity, density, and specific heat, respectively. The hybrid nanofluid models for the thermal and electrical conductivity, specific heat, density, and dynamic viscosity are defined as [63,64]:
k h n f k b f = ( 1 φ 2 ) + 2 φ 2 k m 2 k m 2 k b f l n k m 2 + k b f 2 k b f , k b f k f = ( 1 φ 1 ) + 2 φ 1 k m 1 k m 1 k f l n k m 1 + k f 2 k f ,
σ h n f σ b f = 1 + 3 σ m 2 σ b f 1 φ 2 σ m 2 σ b f + 2 σ m 2 σ b f 1 φ 2 , σ b f σ f = 1 + 3 σ m 1 σ f 1 φ 1 σ m 1 σ f + 2 σ m 1 σ f 1 φ 1 ,
( ρ C p ) h n f ( ρ C p ) f = ( 1 φ 2 ) 1 1 ( ρ C p ) m 1 ( ρ C p ) f φ 1 + φ 2 ( ρ C p ) m 2 ( ρ C p ) f ,
μ h n f μ f = 1 ( 1 φ 1 ) 2.5 ( 1 φ 2 ) 2.5 ,
ρ h n f ρ f = ( 1 φ 2 ) 1 1 ρ m 1 ρ f φ 1 + φ 2 ρ m 2 ρ f .
The basic engineering quantities of interest take the following forms:
S h x ( R e x ) 1 / 2 = ϕ ( 0 ) ,
N u x ( R e x ) 1 / 2 = k h n f k f θ ( 0 ) ,
C f x ( 2 R e x ) 1 / 2 = μ h n f μ f f ( 0 ) ,
C f y ( 2 R e y ) 1 / 2 = μ h n f μ f d c 1.5 g ( 0 ) .
Here, R e x = u w x ν and R e y = V w y ν demonstrate the local Reynolds numbers. Table 1 was used in simulation of this work, where the basic physical and dynamical properties of the hybrid nanofluid used are presented.

3. Problem Solution by Homotopy Analysis Method

In this section, we solve the reduced system of ODEs together with the boundary conditions defined in Equations (15)–(19). These problems are normally very highly nonlinear, and that is one of the reasons that researchers use the semi-analytical techniques for solution purposes. One of the techniques we use here is the homotopy analysis method (HAM) introduced by Liao [67]. Liao used the concept of topology by constructing a continuous mapping between two functions ξ 1 ( x ¯ ) and ξ 2 ( x ¯ ) , which are themselves continuous, defined over the two topological spaces X ¯ and Y ¯ , on which this mapping is constructed. This mapping is mathematically defined as [68]:
Φ ˜ : X ¯ × [ 0 , 1 ] Y ¯ ,
where Φ ˜ [ X ¯ , 0 ] = ξ 1 ( X ¯ ) and Φ ˜ [ X ¯ , 1 ] = ξ 2 ( X ¯ ) hold x ¯ X ¯ .
Now, to solve Equations (15)–(18) together with the boundary conditions (19), we chose f 0 = exp ( ζ ) + 1 , g 0 = γ 0 ( exp ( ζ ) + 1 ) , θ 0 = exp ( ζ ) , and ϕ 0 = exp ( ζ ) as the initial guesses with the corresponding operators, defined as:
L f ¯ ( f ¯ ) = f ¯ f ¯ , L g ¯ ( g ¯ ) = g ¯ g ¯ , L θ ¯ ( θ ¯ ) = θ ¯ θ ¯ , and L ϕ ¯ ( ϕ ¯ ) = ϕ ¯ ϕ ¯ ,
which agree with
L f ¯ ( a 1 + a 2 exp ( ζ ) + a 3 exp ( ζ ) ) = 0 , L g ¯ ( a 4 + a 5 exp ( ζ ) + a 6 exp ( ζ ) ) = 0 , L θ ¯ ( a 7 exp ( ζ ) + a 8 exp ( ζ ) ) = 0 , L ϕ ¯ ( a 9 exp ( ζ ) + a 10 exp ( ζ ) ) = 0 ,
where a i , i = 1 10 are arbitrary constants.

3.1. Convergence Analysis

In this section, we briefly explain the convergence of the method used in Section 3. For this purpose, we present the residual errors computed for the state variables in the form of Table 2. The residuals were computed subject to variations in the convergence control parameters used in the homotopy. These parameters are f , g , θ , and ϕ and were chosen in such away that the square residual is minimized. Using the concept of Liao [69], the mathematical relations for the square residuals take the following forms:
ϵ f = 1 p + 1 l = 0 p N f r = 0 n f ¯ ( ζ ) r = 0 n g ¯ ( ζ ) l Δ ζ 2 ,
ϵ g = 1 p + 1 l = 0 p N g r = 0 n f ¯ ( ζ ) r = 0 n g ¯ ( ζ ) l Δ ζ 2 ,
ϵ θ = 1 p + 1 l = 0 p N θ r = 0 n f ¯ ( ζ ) r = 0 n g ¯ ( ζ ) r = 0 n θ ¯ ( ζ ) l Δ ζ 2 ,
ϵ ϕ = 1 p + 1 l = 0 p N ϕ r = 0 n f ¯ ( ζ ) r = 0 n g ¯ ( ζ ) r = 0 n θ ¯ ( ζ ) r = 0 n ϕ ¯ ( ζ ) l Δ ζ 2 .
The total residual can be obtained from the following empirical relation by combining Equations (32)–(35), and we have [69]:
ϵ s = ϵ f + ϵ g + ϵ θ + ϵ ϕ .
In Equation (36), ϵ s shows the sum of the residuals obtained. For the second-order deformation, we chose Δ ζ = 0.45 and p = 20 for the convergence control parameters f = 2.21453 , g = 2.56672 , θ = 1.56345 , and ϕ = 1.56345 . In Table 2, the various residuals computed are presented for the corresponding state variables. From the table, it is clear that, when the number of iterations increases, the residual error decreases.

4. Results and Discussion

The results obtained by implementing the HAM in the above equations are presented in the form of graphs and tables in this section. The graphical description in Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11, Figure 12, Figure 13, Figure 14, Figure 15, Figure 16 and Figure 17 expressed here shows the impact of various parameters on the state variables. Furthermore, the engineering quantities of interest are expressed in Table 3 and Table 4.
The impact of the mass flux S 0 in the velocity f ( ζ ) of the considered nanofluid (CuO+kerosene oil) and the hybrid nanofluid (Ag+CuO+kerosene oil) is depicted in Figure 2. The values of S 0 were taken to be 0.2 , 0.4 , 0.6 . The graph displays a decreasing dependence of the fluid velocity with increasing S 0 . Both fluids exhibit almost the same dependence. Thus, the increasing S 0 causes a lower distribution in f ( ζ ) .
The dependence of both nanofluids’ dynamics on the varying magnetic field strength M is displayed in Figure 3. The graph shows that the increasing magnetic field strength M results in a reduction in the f ( ζ ) profile. The values of M were taken to be 0.2 , 0.4 , 0.6 . The drop in f ( ζ ) is more prominent at the intermediate values of ζ . The spacing between different curves decreases with the increasing values of M. The increasing M values causes f ( ζ ) to reduce at a greater rate in the hybrid nanofluid compared to the ordinary nanofluid. This shows that the increasing M values result in a reduction of the fluid velocity distribution.
The dependence of f ( ζ ) on the Hall effect parameter m is depicted in Figure 4. The chosen values of m were 0.2 , 0.6 , 0.9 . The graph shows that the increasing values of m increase the profile of f ( ζ ) for both types of nanofluid. The increase is more prominent at the middle values of ζ . The increase in the hybrid nanofluid is the largest one for the highest of m, that is at m = 0.9 . This shows that the increasing Hall effect strength results an increase in the velocity distribution of both types of nanofluid.
The impact of the stretching ratio γ 0 on the velocity distribution f ( ζ ) is plotted in Figure 5. The chosen values of γ 0 were 0.2 , 0.6 , 0.8 . It is apparent that the changing of the γ 0 strength results in an increase in f ( ζ ) . The increasing γ increases the velocity distribution of the simple nanofluid at a higher rate compared to the chosen hybrid nanofluid. The increase is more obvious at the intermediate values of γ 0 .
The impacts of the mass flux parameter S 0 and stretching ratio γ 0 on the other component of the velocity g ( ζ ) are respectively exhibited in Figure 6 and Figure 7. The chosen values of S 0 were 0.1 , 0.3 , 0.5 , and those of γ 0 were 0.2 , 0.6 , 0.8 . It is clear from Figure 6 that the increasing mass flux causes the velocity distribution to reduce along the vertical direction. The drop is more drastic for the hybrid nanofluid compared to the simple nanofluid. The increasing stretching ratio strength γ 0 causes the velocity distribution to increase along the vertical direction, as displayed in Figure 7. The increase in the velocity distribution is more prominent for the hybrid nanofluid compared to the ordinary nanofluid. The spacing between different curves changes with the higher values of γ 0 at the intermediate values of ζ .
The dependence of the chosen simple and hybrid nanofluids temperature θ ( ζ ) on Prandtl number P r and thermophoresis parameter N t is respectively plotted in Figure 8 and Figure 9. The values of both parameters were chosen to be 0.5 , 0.9 , 1.3 . The temperature of both fluids drops with the increasing values of P r . This is due to the thermal diffusion, because it has an inverse relation with P r . On the other hand, when the viscous diffusion increases, the Prandtl number jumps up. As a result, the temperature profile decreases. The temperature drop shows a complex dependence on the increasing P r . At a smaller value of P r , the simple nanofluid temperature drops at a higher rate compared to the hybrid nanofluid temperature drop. At higher values, the temperature drop decreases compared to the smaller separation between the different curves. Thus, the smaller thermal diffusivity associated with a higher P r results in a decrease in the temperature distribution of both types of fluid. Figure 9 shows that the changing thermophoresiscauses the temperature of the fluids to increase. The increasing nature of the fluids’ temperature with the increasing N t is higher for the hybrid nanofluid compared to the single-component nanofluid.
The variation of the fluids’ temperature with the changing Brownian motion parameter N b , stretching ratio γ 0 , and thermal energy relaxation parameter ϵ 1 is displayed in Figure 10, Figure 11 and Figure 12, respectively. Figure 10 shows that the increasing randomness associated with the higher values of N b causes the fluids’ temperature to change. The rate of increase of the hybrid nanofluid is higher with the increasing N b compared to the ordinary nanofluid. The increasing stretching ratio strength associated with the higher γ 0 decreases the fluid temperature, as exhibited in Figure 11. The dropping rate for both fluids almost shows the same trend. The reduction in the fluid temperature is more drastic at the intermediate values of ζ . The thermal energy relaxation exhibited by the increasing values of ϵ 1 also results in the reduction of the fluid temperature distribution θ ( ζ ) . The temperature dropping rate is higher for the hybrid nanofluid compared to the ordinary nanofluid. This is more obvious at the middle values of ζ . Thus, the increasing thermal energy relaxation causes the temperature distribution to reduce.
The variation of the fluid concentration with the changing values of the mass flux relaxation parameter ϵ 2 and stretching ratio γ 0 are plotted respectively in Figure 13 and Figure 14. The values of ϵ 2 were chosen to be 0.0 , 0.3 , 0.6 , whereas those of γ 0 were taken to be 0.2 , 0.4 , 0.6 . The changing mass flux relaxation drops the fluids’ temperature. Initially, the dropping rate for the ordinary nanofluid is higher. As the values of ϵ 2 increase, the dropping rate decreases, as is clear from the reducing separation between the different curves. At the highest ϵ 2 , the concentration profiles for both fluids almost coincide. The drop in the fluids’ concentration is more visible at the middle values of ζ . The increasing stretching ratio γ 0 also results in a reduction in the fluids’ concentration, as shown in Figure 14. The drop in the hybrid nanofluid concentration is more obvious compared to the ordinary nanofluid.
The dependence of the fluid concentration ϕ on the increasing strength of N b , N t , and the Schmidt number S c is respectively plotted in Figure 15, Figure 16 and Figure 17. Figure 15 displays a strange dependence of the fluid concentration on N b . The chosen values of N b were 0.2 , 0.6 , 1.2 . At smaller ζ , the fluid concentration first drops, then rises with increasing ζ and reaches a maximum, then drops again with higher values of ζ . The rate of the concentration drop decreases with the increasing strength of the randomness associated with the higher N b values, as exhibited by the decreasing separation between the different curves. The changing thermophoresis associated with the higher values of N t causes the fluid concentration to change. The increasing rate of the fluid concentration is more prominent up to about ζ = 3.0 . Thus, the increasing thermophoresis affects the hybrid nanofluid concentration more drastically compared to the ordinary fluid. The impact of the increasing Schmidt number S c on the fluid concentration ϕ ( ζ ) is displayed in Figure 17. The chosen values of S c were 0.6 , 1.2 , 1.8 . The graph exhibits a drop in the concentration profiles of both fluids. The rate of decrease of the ordinary nanofluid is more drastic with the increasing S c compared to the hybrid nanofluid. Thus, the increasing kinematic viscosity associated with higher values of S c affects the concentration profile of the ordinary nanofluid more drastically compared to the hybrid nanofluid. The variation is more prominent at the intermediate values of ζ .

5. Tables’ Discussion

In this section, various coefficients of engineering interest are displayed in Table 3 and Table 4 numerically. The variations of these coefficients were analyzed under the variations of the ratio parameter γ , thermal relaxation parameter ϵ 1 , and concentration relaxation parameter ϵ 2 . Table 3 shows the opposite behavior with the increasing values of the ratio parameter. In Table 4, the relaxation parameter of heat with increasing values shows that the heat transfer is increasing. A quite similar trend for the mass transfer rate is observed in Table 4 for the mass relaxation parameter.

6. Conclusions

This section is concerned with the results obtained in this work. In this work, the CCDD model was used to simulate the Ag+CuO+kerosene oil hybrid nanofluid passing a bidirectionally stretching sheet, under the influence of the variable magnetic parameter. The fluid was chosen to be electrically conductive, where the current produced (Hall current) was taken into account. The results obtained through graphs and tables are summarized below:
  • When the Hall current increases, the velocity gradient increases, while the dimensionless mass flux and the magnetic parameter decrease this profile.
  • Both the velocity and temperature profiles fall with the increasing values of the ratio parameter.
  • The higher values of the Prandtl number decrease the layer thickness of the thermal boundary.
  • Both the temperature and concentration boundary layer show an opposite trend for the larger values of the Brownian motion parameter.
  • The thermal and concentration parameters decrease the temperature and concentration boundary layers with their increasing values.
  • At the surface, both the temperature and concentration profiles show an increasing trend for various values of the thermal and concentration boundary layers.
  • The skin friction shows the opposite trend for the larger values of the ratio parameter γ 0 .
  • The graphs show that the hybrid nanofluid shows a fast increasing/decreasing trend, which implies further applications in cooling.
  • The efficiency of the implemented technique is proven in Section 3.1 through the squared residual errors presented in Table 2.

Author Contributions

Conceptualization, H.A.E.-W.K., and A.U.; software, M.A., R.A.Z., and A.U.; writing—original draft preparation, A.U., R.A.Z., and H.A.E.-W.K.; writing—review and editing, H.A.E.-W.K., M.A., and A.U.; formal analysis, M.A., H.A.E.-W.K., and R.A.Z.; validation, H.A.E.-W.K., A.U., M.A., and R.A.Z.; methodology, A.U., M.A., and H.A.E.-W.K.; investigation, H.A.E.-W.K., M.A., and A.U.; resources, A.U., and R.A.Z.; project administration, A.U., and H.A.E.-W.K.; funding acquisition, R.A.Z., and H.A.E.-W.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research work was funded by Institutional Fund projects under Grant No. (IFPIP: 1588-305-1443). The authors gratefully acknowledge the technical and financial support provided by the Ministry of Education and King Abdulaziz University, DSR, Jeddah, Saudi Arabia.

Data Availability Statement

The data used to support the findings of this study are available within the article.

Acknowledgments

This research work was funded by Institutional Fund projects under Grant No. (IFPIP: 1588-305-1443). The authors gratefully acknowledge the technical and financial support provided by the Ministry of Education and King Abdulaziz University, DSR, Jeddah, Saudi Arabia.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

σ Electrical conductivity
CConcentration
BUniform magnetic parameter strength T
N u x Local Nusselt number
R e x Local Reynolds number
C f x , C f y Local skin frictions
S h x Sherwood number along the x-axis
P r Prandtl number
u w Stretching velocity along the x-axis
v w Stretching velocity along the y-axis
TFluid temperature (K)
θ Dimensionless temperature
fDimensionless velocity along the x-axis
gDimensionless velocity along the y-axis
θ Dimensionless temperature
ϕ Dimensionless concentration
ρ Density
μ Dynamic viscosity
C p Specific heat
( C p ) h n f Hybrid nanofluid specific heat
μ h n f Hybrid nanofluid dynamic viscosity
ρ h n f Hybrid nanofluid density
J w Mass flux
Condition at infinity
0Reference condition
wCondition at the wall
x, y, and zCoordinates
ζ Similarity variable
γ 0 Stretching ratio parameter
S c Schmidt number
tTime (s)
k h n f Hybrid nanofluid thermal conductivity
k b f Base fluid thermal conductivity
σ b f Base fluid electrical conductivity
σ h n f Hybrid nanofluid electrical conductivity
MMagnetic parameter
mHall parameter
V ˜ Velocity field
q ˜ Heat flux
j ˜ Mass flux
Λ e Energy relaxation parameter
c, dReferences of velocity
LCharacteristic length
B ( x , y ) Variable magnetic field
B 0 Strength of the magnetic field
Λ c Concentration relaxation parameter
D B Brownian motion parameter
τ Adiabatic index
S 0 Suction/injection parameter
N t Thermo-migration parameter
ϵ 1 Thermal relaxation parameter
ϵ 2 Mass relaxation parameter
D T Thermophoretic parameter

References

  1. Choi, S.U.; Eastman, J.A. Enhancing Thermal Conductivity of Fluids with Nanoparticles; Technical report; Argonne National Lab.(ANL): Argonne, IL, USA, 1995.
  2. Eastman, J.A.; Choi, S.; Li, S.; Yu, W.; Thompson, L. Anomalously increased effective thermal conductivities of ethylene glycol-based nanofluids containing copper nanoparticles. Appl. Phys. Lett. 2001, 78, 718–720. [Google Scholar] [CrossRef]
  3. Maïga, S.E.B.; Nguyen, C.T.; Galanis, N.; Roy, G. Heat transfer behaviors of nanofluids in a uniformly heated tube. Superlattices Microstruct. 2004, 35, 543–557. [Google Scholar] [CrossRef]
  4. Khan, W.; Pop, I. Boundary-layer flow of a nanofluid past a stretching sheet. Int. J. Heat Mass Transf. 2010, 53, 2477–2483. [Google Scholar] [CrossRef]
  5. Seth, G.; Kumar, B.; Nandkeolyar, R. MHD mixed convection stagnation point flow of a micropolar nanofluid adjacent to stretching sheet: A revised model with successive linearization method. J. Nanofluids 2019, 8, 620–630. [Google Scholar] [CrossRef]
  6. Rashidi, M.M.; Ghahremanian, S.; Toghraie, D.; Roy, P. Effect of solid surface structure on the condensation flow of Argon in rough nanochannels with different roughness geometries using molecular dynamics simulation. Int. Commun. Heat Mass Transf. 2020, 117, 104741. [Google Scholar] [CrossRef]
  7. Mansoury, D.; Ilami Doshmanziari, F.; Rezaie, S.; Rashidi, M.M. Effect of Al2O3/water nanofluid on performance of parallel flow heat exchangers. J. Therm. Anal. Calorim. 2019, 135, 625–643. [Google Scholar] [CrossRef]
  8. Hatami, M.; Zhou, J.; Geng, J.; Song, D.; Jing, D. Optimization of a lid-driven T-shaped porous cavity to improve the nanofluids mixed convection heat transfer. J. Mol. Liq. 2017, 231, 620–631. [Google Scholar] [CrossRef]
  9. Jana, S.; Salehi-Khojin, A.; Zhong, W.H. Enhancement of fluid thermal conductivity by the addition of single and hybrid nano-additives. Thermochim. Acta 2007, 462, 45–55. [Google Scholar] [CrossRef]
  10. Madhesh, D.; Parameshwaran, R.; Kalaiselvam, S. Experimental investigation on convective heat transfer and rheological characteristics of Cu–TiO2 hybrid nanofluids. Exp. Therm. Fluid Sci. 2014, 52, 104–115. [Google Scholar] [CrossRef]
  11. Ho, C.; Huang, J.; Tsai, P.; Yang, Y. Preparation and properties of hybrid water-based suspension of Al2O3 nanoparticles and MEPCM particles as functional forced convection fluid. Int. Commun. Heat Mass Transf. 2010, 37, 490–494. [Google Scholar] [CrossRef]
  12. Suresh, S.; Venkitaraj, K.; Selvakumar, P.; Chandrasekar, M. Synthesis of Al2O3–Cu/water hybrid nanofluids using two step method and its thermo physical properties. Colloids Surfaces A: Physicochem. Eng. Asp. 2011, 388, 41–48. [Google Scholar] [CrossRef]
  13. Labib, M.N.; Nine, M.J.; Afrianto, H.; Chung, H.; Jeong, H. Numerical investigation on effect of base fluids and hybrid nanofluid in forced convective heat transfer. Int. J. Therm. Sci. 2013, 71, 163–171. [Google Scholar] [CrossRef]
  14. Moghadassi, A.; Ghomi, E.; Parvizian, F. A numerical study of water based Al2O3 and Al2O3–Cu hybrid nanofluid effect on forced convective heat transfer. Int. J. Therm. Sci. 2015, 92, 50–57. [Google Scholar] [CrossRef]
  15. Daniel, Y.S.; Aziz, Z.A.; Ismail, Z.; Salah, F. Thermal radiation on unsteady electrical MHD flow of nanofluid over stretching sheet with chemical reaction. J. King Saud Univ.-Sci. 2019, 31, 804–812. [Google Scholar] [CrossRef]
  16. Khan, I. Shape effects of MoS2 nanoparticles on MHD slip flow of molybdenum disulphide nanofluid in a porous medium. J. Mol. Liq. 2017, 233, 442–451. [Google Scholar] [CrossRef]
  17. Chamkha, A.J.; Ismael, M.A. Magnetic field effect on mixed convection in lid-driven trapezoidal cavities filled with a Cu–water nanofluid with an aiding or opposing side wall. J. Therm. Sci. Eng. Appl. 2016, 8, 031009. [Google Scholar] [CrossRef] [Green Version]
  18. Alsabery, A.I.; Ismael, M.A.; Chamkha, A.J.; Hashim, I. Mixed convection of Al2O3-water nanofluid in a double lid-driven square cavity with a solid inner insert using Buongiorno’s two-phase model. Int. J. Heat Mass Transf. 2018, 119, 939–961. [Google Scholar] [CrossRef]
  19. Ellahi, R. The effects of MHD and temperature dependent viscosity on the flow of non-Newtonian nanofluid in a pipe: Analytical solutions. Appl. Math. Model. 2013, 37, 1451–1467. [Google Scholar] [CrossRef]
  20. Sheikholeslami, M.; D Ganji, D. Magnetohydrodynamic flow in a permeable channel filled with nanofluid. Sci. Iran. 2014, 21, 203–212. [Google Scholar]
  21. Seyyedi, S.M.; Dogonchi, A.; Hashemi-Tilehnoee, M.; Asghar, Z.; Waqas, M.; Ganji, D. A computational framework for natural convective hydromagnetic flow via inclined cavity: An analysis subjected to entropy generation. J. Mol. Liq. 2019, 287, 110863. [Google Scholar] [CrossRef]
  22. Rashidi, M.; Abelman, S.; Mehr, N.F. Entropy generation in steady MHD flow due to a rotating porous disk in a nanofluid. Int. J. Heat Mass Transf. 2013, 62, 515–525. [Google Scholar] [CrossRef]
  23. Sheikholeslami, M.; Ellahi, R. Three dimensional mesoscopic simulation of magnetic field effect on natural convection of nanofluid. Int. J. Heat Mass Transf. 2015, 89, 799–808. [Google Scholar] [CrossRef]
  24. Biswas, N.; Mandal, D.K.; Manna, N.K.; Benim, A.C. Enhanced energy and mass transport dynamics in a thermo-magneto-bioconvective porous system containing oxytactic bacteria and nanoparticles: Cleaner energy application. Energy 2023, 263, 125775. [Google Scholar] [CrossRef]
  25. Manna, N.K.; Biswas, N.; Mandal, D.K.; Sarkar, U.; Öztop, H.F.; Abu-Hamdeh, N. Impacts of heater-cooler position and Lorentz force on heat transfer and entropy generation of hybrid nanofluid convection in quarter-circular cavity. Int. J. Numer. Methods Heat Fluid Flow, 2022; ahead-of-print. [Google Scholar] [CrossRef]
  26. Biswas, N.; Mondal, M.K.; Manna, N.K.; Mandal, D.K.; Chamkha, A.J. Implementation of partial magnetic fields to magneto-thermal convective systems operated using hybrid-nanoliquid and porous media. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2022, 236, 5687–5704. [Google Scholar] [CrossRef]
  27. Mandal, D.K.; Biswas, N.; Manna, N.K.; Gorla, R.S.R.; Chamkha, A.J. Hybrid nanofluid magnetohydrodynamic mixed convection in a novel W-shaped porous system. Int. J. Numer. Methods Heat Fluid Flow, 2022; ahead-of-print. [Google Scholar] [CrossRef]
  28. Biswas, N.; Mandal, D.K.; Manna, N.K.; Benim, A.C. Magneto-hydrothermal triple-convection in a W-shaped porous cavity containing oxytactic bacteria. Sci. Rep. 2022, 12, 1–30. [Google Scholar] [CrossRef] [PubMed]
  29. Mandal, D.K.; Biswas, N.; Manna, N.K.; Gayen, D.K.; Gorla, R.S.R.; Chamkha, A.J. Thermo-fluidic transport process in a novel M-shaped cavity packed with non-Darcian porous medium and hybrid nanofluid: Application of artificial neural network (ANN). Phys. Fluids 2022, 34, 033608. [Google Scholar] [CrossRef]
  30. Mandal, D.; Biswas, N.; Manna, N.; Gorla, R. Magneto-thermal convection of hybrid nanofluid in a non-Darcian porous complex wavy enclosure. Eur. Phys. J. Spec. Top. 2022, 213, 2695–2712. [Google Scholar] [CrossRef]
  31. Sakiadis, B.C. Boundary-layer behavior on continuous solid surfaces: I. Boundary-layer equations for two-dimensional and axisymmetric flow. AIChE J. 1961, 7, 26–28. [Google Scholar] [CrossRef]
  32. Sakiadis, B. Boundary-layer behavior on continuous solid surfaces: II. The boundary layer on a continuous flat surface. AiChE J. 1961, 7, 221–225. [Google Scholar] [CrossRef]
  33. Sakiadis, B. Boundary-layer behavior on continuous solid surfaces: III. The boundary layer on a continuous cylindrical surface. AiChE J. 1961, 7, 467–472. [Google Scholar] [CrossRef]
  34. Crane, L.J. Flow past a stretching plate. Z. Angew. Math. Phys. ZAMP 1970, 21, 645–647. [Google Scholar] [CrossRef]
  35. Gupta, P.; Gupta, A. Heat and mass transfer on a stretching sheet with suction or blowing. Can. J. Chem. Eng. 1977, 55, 744–746. [Google Scholar] [CrossRef]
  36. Devi, R.; Neeraja, A.; Reddy, N.B. Radiation effect on MHD slip flow past a stretching sheet with variable viscosity and heat source/sink. Int. J. Sci. Innov. Math. Res. 2015, 3, 8–17. [Google Scholar]
  37. Magyari, E.; Keller, B. Heat and mass transfer in the boundary layers on an exponentially stretching continuous surface. J. Phys. D: Appl. Phys. 1999, 32, 577. [Google Scholar] [CrossRef]
  38. Wang, C. Liquid film on an unsteady stretching surface. Q. Appl. Math. 1990, 48, 601–610. [Google Scholar] [CrossRef] [Green Version]
  39. Miklavčič, M.; Wang, C. Viscous flow due to a shrinking sheet. Q. Appl. Math. 2006, 64, 283–290. [Google Scholar] [CrossRef] [Green Version]
  40. Bhattacharyya, K. Boundary layer flow and heat transfer over an exponentially shrinking sheet. Chin. Phys. Lett. 2011, 28, 074701. [Google Scholar] [CrossRef]
  41. Bhattacharyya, K.; Vajravelu, K. Stagnation-point flow and heat transfer over an exponentially shrinking sheet. Commun. Nonlinear Sci. Numer. Simul. 2012, 17, 2728–2734. [Google Scholar] [CrossRef]
  42. Bachok, N.; Ishak, A.; Pop, I. Boundary layer stagnation-point flow and heat transfer over an exponentially stretching/shrinking sheet in a nanofluid. Int. J. Heat Mass Transf. 2012, 55, 8122–8128. [Google Scholar] [CrossRef]
  43. Sharma, R.; Ishak, A.; Nazar, R.; Pop, I. Boundary layer flow and heat transfer over a permeable exponentially shrinking sheet in the presence of thermal radiation and partial slip. J. Appl. Fluid Mech. 2014, 7, 125–134. [Google Scholar]
  44. Jusoh, R.; Nazar, R.; Pop, I. Magnetohydrodynamic boundary layer flow and heat transfer of nanofluids past a bidirectional exponential permeable stretching/shrinking sheet with viscous dissipation effect. J. Heat Transf. 2019, 141, 012406. [Google Scholar] [CrossRef]
  45. Dero, S.; Rohni, A.M.; Saaban, A. MHD micropolar nanofluid flow over an exponentially stretching/shrinking surface: Triple solutions. J. Adv. Res. Fluid Mech. Therm. Sci. 2019, 56, 165–174. [Google Scholar]
  46. Lund, L.A.; Omar, Z.; Khan, I. Quadruple solutions of mixed convection flow of magnetohydrodynamic nanofluid over exponentially vertical shrinking and stretching surfaces: Stability analysis. Comput. Methods Programs Biomed. 2019, 182, 105044. [Google Scholar] [CrossRef] [PubMed]
  47. Christov, C. On frame indifferent formulation of the Maxwell–Cattaneo model of finite-speed heat conduction. Mech. Res. Commun. 2009, 36, 481–486. [Google Scholar] [CrossRef]
  48. Rauf, A.; Abbas, Z.; Shehzad, S. Utilization of Maxwell-Cattaneo law for MHD swirling flow through oscillatory disk subject to porous medium. Appl. Math. Mech. 2019, 40, 837–850. [Google Scholar] [CrossRef]
  49. Sui, J.; Zheng, L.; Zhang, X. Boundary layer heat and mass transfer with Cattaneo–Christov double-diffusion in upper-convected Maxwell nanofluid past a stretching sheet with slip velocity. Int. J. Therm. Sci. 2016, 104, 461–468. [Google Scholar] [CrossRef]
  50. Hayat, T.; Muhammad, T.; Alsaedi, A.; Ahmad, B. Three-dimensional flow of nanofluid with Cattaneo–Christov double diffusion. Results Phys. 2016, 6, 897–903. [Google Scholar] [CrossRef] [Green Version]
  51. Hayat, T.; Aziz, A.; Muhammad, T.; Alsaedi, A. Model and comparative study for flow of viscoelastic nanofluids with Cattaneo–Christov double diffusion. PLoS ONE 2017, 12, e0168824. [Google Scholar] [CrossRef] [Green Version]
  52. Hayat, T.; Qayyum, S.; Shehzad, S.A.; Alsaedi, A. Chemical reaction and heat generation/absorption aspects in flow of Walters-B nanofluid with Cattaneo–Christov double-diffusion. Results Phys. 2017, 7, 4145–4152. [Google Scholar] [CrossRef]
  53. Hayat, T.; Aziz, A.; Muhammad, T.; Alsaedi, A. Three-dimensional flow of Prandtl fluid with Cattaneo–Christov double diffusion. Results Phys. 2018, 9, 290–296. [Google Scholar] [CrossRef]
  54. Ibrahim, W. Three dimensional rotating flow of Powell-Eyring nanofluid with non-Fourier’s heat flux and non-Fick’s mass flux theory. Results Phys. 2018, 8, 569–577. [Google Scholar] [CrossRef]
  55. Irfan, M.; Khan, M.; Khan, W. On model for three-dimensional Carreau fluid flow with Cattaneo–Christov double diffusion and variable conductivity: A numerical approach. J. Braz. Soc. Mech. Sci. Eng. 2018, 40, 1–10. [Google Scholar] [CrossRef]
  56. Iqbal, Z.; Khan, M.; Ahmed, A.; Ahmed, J.; Hafeez, A. Thermal energy transport in Burgers nanofluid flow featuring the Cattaneo–Christov double diffusion theory. Appl. Nanosci. 2020, 10, 5331–5342. [Google Scholar] [CrossRef]
  57. Ahmad, I.; Faisal, M.; Javed, T. Bi-directional stretched nanofluid flow with Cattaneo–Christov double diffusion. Results Phys. 2019, 15, 102581. [Google Scholar] [CrossRef]
  58. Rauf, A.; Shehzad, S.; Abbas, Z.; Hayat, T. Unsteady three-dimensional MHD flow of the micropolar fluid over an oscillatory disk with Cattaneo–Christov double diffusion. Appl. Math. Mech. 2019, 40, 1471–1486. [Google Scholar] [CrossRef]
  59. Fiza, M.; Alsubie, A.; Ullah, H.; Hamadneh, N.N.; Islam, S.; Khan, I. Three-dimensional rotating flow of MHD Jeffrey fluid flow between two parallel plates with impact of hall current. Math. Probl. Eng. 2021, 2021, 6626411. [Google Scholar] [CrossRef]
  60. Shamshuddin, M.; Salawu, S.; Ogunseye, H.; Mabood, F. Dissipative power-law fluid flow using spectral quasi linearization method over an exponentially stretchable surface with Hall current and power-law slip velocity. Int. Commun. Heat Mass Transf. 2020, 119, 104933. [Google Scholar] [CrossRef]
  61. Ramzan, M.; Gul, H.; Chung, J.D.; Kadry, S.; Chu, Y.M. Significance of Hall effect and Ion slip in a three-dimensional bioconvective Tangent hyperbolic nanofluid flow subject to Arrhenius activation energy. Sci. Rep. 2020, 10, 1–15. [Google Scholar] [CrossRef]
  62. Zangooee, M.R.; Hosseinzadeh, K.; Ganj, D.D. Investigation of three-dimensional hybrid nanofluid flow affected by nonuniform MHD over exponential stretching/shrinking plate. Nonlinear Eng. 2022, 11, 143–155. [Google Scholar] [CrossRef]
  63. Turkyilmazoglu, M.; Pop, I. Heat and mass transfer of unsteady natural convection flow of some nanofluids past a vertical infinite flat plate with radiation effect. Int. J. Heat Mass Transf. 2013, 59, 167–171. [Google Scholar] [CrossRef]
  64. Turkyilmazoglu, M. On the transparent effects of Buongiorno nanofluid model on heat and mass transfer. Eur. Phys. J. Plus 2021, 136, 1–15. [Google Scholar] [CrossRef]
  65. Algehyne, E.A.; Alharbi, A.F.; Saeed, A.; Dawar, A.; Ramzan, M.; Kumam, P. Analysis of the MHD partially ionized GO-Ag/water and GO-Ag/kerosene oil hybrid nanofluids flow over a stretching surface with Cattaneo–Christov double diffusion model: A comparative study. Int. Commun. Heat Mass Transf. 2022, 136, 106205. [Google Scholar] [CrossRef]
  66. Ahmad, S.; Ali, K.; Ashraf, M.; Jamshed, W. Thermal Characteristics of Kerosene Oil based Hybrid Nanofluids (Ag-MnZnFe2O4): A Comprehensive Study. Front. Energy Res. 2022, 10, 978819. [Google Scholar] [CrossRef]
  67. Liao, S.J. The Proposed Homotopy Analysis Technique for the Solution of Nonlinear Problems. Ph.D. Thesis, Shanghai Jiao Tong University, Shanghai, China, 1992. [Google Scholar]
  68. Liao, S. Homotopy analysis method: A new analytical technique for nonlinear problems. Commun. Nonlinear Sci. Numer. Simul. 1997, 2, 95–100. [Google Scholar] [CrossRef]
  69. Liao, S. An optimal homotopy-analysis approach for strongly nonlinear differential equations. Commun. Nonlinear Sci. Numer. Simul. 2010, 15, 2003–2016. [Google Scholar] [CrossRef]
Figure 1. Geometrical description of the study.
Figure 1. Geometrical description of the study.
Symmetry 15 00166 g001
Figure 2. Influence of S 0 on f ( ζ ) .
Figure 2. Influence of S 0 on f ( ζ ) .
Symmetry 15 00166 g002
Figure 3. Influence of M on f ( ζ ) .
Figure 3. Influence of M on f ( ζ ) .
Symmetry 15 00166 g003
Figure 4. Dependence of f ( η ) on m.
Figure 4. Dependence of f ( η ) on m.
Symmetry 15 00166 g004
Figure 5. Variation of f ( η ) with γ 0 .
Figure 5. Variation of f ( η ) with γ 0 .
Symmetry 15 00166 g005
Figure 6. Impact of S 0 on g ( ζ ) .
Figure 6. Impact of S 0 on g ( ζ ) .
Symmetry 15 00166 g006
Figure 7. Dependence of g ( ζ ) on γ 0 .
Figure 7. Dependence of g ( ζ ) on γ 0 .
Symmetry 15 00166 g007
Figure 8. Variations in θ ( ζ ) with P r .
Figure 8. Variations in θ ( ζ ) with P r .
Symmetry 15 00166 g008
Figure 9. Impact of N t on θ ( ζ ) .
Figure 9. Impact of N t on θ ( ζ ) .
Symmetry 15 00166 g009
Figure 10. θ ( ζ ) versus N b .
Figure 10. θ ( ζ ) versus N b .
Symmetry 15 00166 g010
Figure 11. Dependence of θ ( ζ ) on γ 0 .
Figure 11. Dependence of θ ( ζ ) on γ 0 .
Symmetry 15 00166 g011
Figure 12. Influence of ϵ 1 on θ ( ζ ) .
Figure 12. Influence of ϵ 1 on θ ( ζ ) .
Symmetry 15 00166 g012
Figure 13. Influence of ϵ 2 on ϕ ( ζ ) .
Figure 13. Influence of ϵ 2 on ϕ ( ζ ) .
Symmetry 15 00166 g013
Figure 14. Influence of γ 0 on ϕ ( ζ ) .
Figure 14. Influence of γ 0 on ϕ ( ζ ) .
Symmetry 15 00166 g014
Figure 15. Influence of N b on ϕ ( ζ ) .
Figure 15. Influence of N b on ϕ ( ζ ) .
Symmetry 15 00166 g015
Figure 16. Influence of N t on ϕ ( ζ ) .
Figure 16. Influence of N t on ϕ ( ζ ) .
Symmetry 15 00166 g016
Figure 17. Influence of S c on ϕ ( ζ ) .
Figure 17. Influence of S c on ϕ ( ζ ) .
Symmetry 15 00166 g017
Table 1. Some basic physical and dynamical properties of the nanofluid used in this work [65,66].
Table 1. Some basic physical and dynamical properties of the nanofluid used in this work [65,66].
PropertiesAgCuOKerosene Oil
k W m K 429180.145
C p J k g K 2355402090
ρ k g m 3 10,5006500783
Table 2. Squared residual error values subject to the auxiliary variables’ values.
Table 2. Squared residual error values subject to the auxiliary variables’ values.
n ϵ f ϵ g ϵ θ ϵ ϕ
2 5.26 × 10 4 4.34 × 10 5 5.78 × 10 4 5.78 × 10 4
5 6.31 × 10 5 2.49 × 10 6 4.87 × 10 5 4.87 × 10 5
10 5.47 × 10 5 6.91 × 10 6 3.14 × 10 5 3.14 × 10 5
15 6.69 × 10 6 8.76 × 10 7 5.43 × 10 6 5.43 × 10 6
20 7.98 × 10 6 4.28 × 10 7 4.90 × 10 6 4.90 × 10 6
25 6.27 × 10 7 3.76 × 10 8 5.67 × 10 7 5.67 × 10 7
30 4.14 × 10 7 2.93 × 10 8 3.49 × 10 7 3.49 × 10 7
Table 3. Computed values of C f x and C f y , under the variation of the ratio parameter ( γ 0 ) .
Table 3. Computed values of C f x and C f y , under the variation of the ratio parameter ( γ 0 ) .
γ 0 0.20.40.60.81.0
C f x −1.04025−1.07483−1.09753−1.11643−1.16989
C f y −2.35546−1.99025−1.50623−1.31678−1.17889
Table 4. Computed values of the heat and mass transfer rate for various values of ϵ 1 and ϵ 2 .
Table 4. Computed values of the heat and mass transfer rate for various values of ϵ 1 and ϵ 2 .
ϵ 1 , ϵ 2 0.00.20.40.60.8
θ ( 0 ) −0.52003−0.52098−0.52941−0.53452−0.53742
ϕ ( 0 ) −0.51341−0.51661−0.51946−0.52367−0.53854
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

ZeinEldin, R.A.; Ullah, A.; Khalifa, H.A.E.-W.; Ayaz, M. Analytical Study of the Energy Loss Reduction during Three-Dimensional Engine Oil-Based Hybrid Nanofluid Flow by Using Cattaneo–Christov Model. Symmetry 2023, 15, 166. https://doi.org/10.3390/sym15010166

AMA Style

ZeinEldin RA, Ullah A, Khalifa HAE-W, Ayaz M. Analytical Study of the Energy Loss Reduction during Three-Dimensional Engine Oil-Based Hybrid Nanofluid Flow by Using Cattaneo–Christov Model. Symmetry. 2023; 15(1):166. https://doi.org/10.3390/sym15010166

Chicago/Turabian Style

ZeinEldin, Ramadan A., Asad Ullah, Hamiden Abd El-Wahed Khalifa, and Muhammad Ayaz. 2023. "Analytical Study of the Energy Loss Reduction during Three-Dimensional Engine Oil-Based Hybrid Nanofluid Flow by Using Cattaneo–Christov Model" Symmetry 15, no. 1: 166. https://doi.org/10.3390/sym15010166

APA Style

ZeinEldin, R. A., Ullah, A., Khalifa, H. A. E. -W., & Ayaz, M. (2023). Analytical Study of the Energy Loss Reduction during Three-Dimensional Engine Oil-Based Hybrid Nanofluid Flow by Using Cattaneo–Christov Model. Symmetry, 15(1), 166. https://doi.org/10.3390/sym15010166

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop