Thermally stratified flow of hybrid nanofluids with radiative heat transport and slip mechanism: multiple solutions

Rai Sajjad Saif, HashimMaria Zaman and Muhammad Ayaz

1 Department of Humanities and Sciences, School of Electrical Engineering and Computer Science(SEECS), National University of Sciences and Technology (NUST), Islamabad, Pakistan

2 Department of Pure & Applied Mathematics, University of Haripur, 22620, Haripur, Pakistan

3 Department of Mathematics & Statistics, Riphah International University Islamabad, Pakistan

4 Department of Mathematics, Abdul Wali Khan University, Mardan, 23200, Pakistan

Abstract Research on flow and heat transfer of hybrid nanofluids has gained great significance due to their efficient heat transfer capabilities.In fact,hybrid nanofluids are a novel type of fluid designed to enhance heat transfer rate and have a wide range of engineering and industrial applications.Motivated by this evolution, a theoretical analysis is performed to explore the flow and heat transport characteristics of Cu/Al2O3 hybrid nanofluids driven by a stretching/shrinking geometry.Further, this work focuses on the physical impacts of thermal stratification as well as thermal radiation during hybrid nanofluid flow in the presence of a velocity slip mechanism.The mathematical modelling incorporates the basic conservation laws and Boussinesq approximations.This formulation gives a system of governing partial differential equations which are later reduced into ordinary differential equations via dimensionless variables.An efficient numerical solver, known as bvp4c in MATLAB, is utilized to acquire multiple (upper and lower)numerical solutions in the case of shrinking flow.The computed results are presented in the form of flow and temperature fields.The most significant findings acquired from the current study suggest that multiple solutions exist only in the case of a shrinking surface until a critical/turning point.Moreover, solutions are unavailable beyond this turning point, indicating flow separation.It is found that the fluid temperature has been impressively enhanced by a higher nanoparticle volume fraction for both solutions.On the other hand, the outcomes disclose that the wall shear stress is reduced with higher magnetic field in the case of the second solution.The simulation outcomes are in excellent agreement with earlier research,with a relative error of less than 1%.

Keywords: multiple numerical solutions, hybrid nanofluids, thermal radiation, thermal stratification, slip conditions, heat transfer

1.Introduction

The thermal performance of fluids employed for heat transport analysis is of main concern in various thermal processes used in industry.In recent years,several researchers put their efforts into improving the thermal performance of these processes by employing different techniques.Correspondingly,it was shown that the thermal features of ordinary fluids were substantially changed with the addition of nanosized particles of single type, known as nanoparticles.Such nanoparticles could be made up of oxides, carbides, carbon nanotubes,and metals.Nanofluids are a stable mixture of base fluid and suspended nanoparticles(up to 100 nm).Engine oil,water, and ethylene glycol may be taken as base fluids.As expected, nanofluids displayed higher heat dissipation performance.In addition, they have shown greater increases in their thermal physical characteristics, like thermal conductivity, etc.Treatment of cancer, medicine, transformer cooling, solar collectors, heat exchangers, nuclear reactor cooling, freezers, electronics cooling, and automobiles are just a few examples of nanofluid applications.Due to the successful performance of nanofluids in the heat transport mechanism, numerous studies have reported on the heat transfer of nanofluids.Choi [1] was the first one to spotlight the notion and implications of using nanoparticles in a conventional fluid to improve heat transfer.After that, Buongiorno [2] and Tiwari and Das [3] proposed two common hydrodynamic models for nanofluids.Buongiorno’s model considered the impacts of Brownian motion and thermophoresis.On the other hand, the Tiwari and Das model explored the nature of nanofluids by volume fraction considering the nanoparticles.Several authors, including Kuznetsov and Nield[4],Makinde and Aziz[5],Mustafa et al[6],Khan and Aziz [7], and Hashim and Khan [8], investigated the various physical characteristics of nanofluids using Buongiorno’s model.On the other hand, the single-phase model proposed by Tiwari and Das [3] has been applied by various authors.Rashad et al [9] utilized the numerical technique to explore the mixed convection flow of copper–water nanofluid in a rectangular cavity within a porous medium.The natural convection flow of an Al2O3/water nanofluid along with heat transport in an L-shaped enclosure has been presented by Mohebbi and Rashidi [10].Later,Bhatti et al [11] numerically investigated the 3D unsteady flow of a nanofluid with gyrotactic microorganisms driven by cylindrical geometry.Many researchers have explored the flow and heat transport properties in various geometries using nanofluids in recent years due to the effectiveness of nanofluids, including Hashim et al [12], Hamid et al [13], Hafeez et al [14], etc.

These days,the research community focuses on radiative heat transfer and flow processes in energy conversion systems that run at high temperatures, owing to their remarkable performance in a variety of fields of science, including satellites, missiles, various aircraft propulsion devices, gas turbines, and nuclear power plants.For a device operating at an above-average temperature where radiation from heated walls and working fluid is different, the effect of thermal radiation is highly important in the transfer of heat flow.Heat transfer features under the influence of thermal radiation have been investigated by several researchers.The influences of buoyancy force and thermal radiation on stagnation point flow past a stretching sheet were examined by Pal [15].Furthermore, Bidin and Nazar [16] numerically investigated the effects of thermal radiation on flow and heat transport analysis due to an exponentially extending surface.Later,Aziz [17] conducted a numerical study on the flow and heat transport mechanism for viscous fluid flow over an unsteady stretching surface.Pal and Mondal [18] analysed non-Darcy flow over a stretching plate incorporating thermal radiation effects.Dogonchi and Ganji[19]addressed the heat transport characteristics of a nanofluid past a stretching surface in the presence of thermal radiation.Lin et al[20]explored the flow and heat transport characteristics of copper–water nanofluid subject to thermal radiation and nanoparticle shape factor effects.Waqas et al [21] presented a numerical study to address the flow of non-Newtonian Carreau fluid subject to thermal radiation using revised nanofluid model.Later,Sheikholeslami et al[22]employed the control volume-based finite element numerical method to study the flow of a nanofluid through a wavy chamber with thermal radiation impacts.The closed-form solutions for two-dimensional flow and heat transfer analysis in the presence of thermal radiation past a vertical plate have been examined by Turkyilmazoglu [23].

Stratification occurs because of temperature gradients or the mixture of multiple fluids having different densities in engineering and industrial mechanisms.Exploring the mechanism of thermal stratification during flow and heat transport in nanofluids is of tremendous physical interest.This process is essential in the disciplines of lake thermohydraulics, salinity and thermal stratification mechanisms in oceans, heat rejection into the surrounding environment via rivers, agriculture fields, volcanic flows, and industries such as reservoirs, industrial food and salinity, atmosphere involving heterogeneous mixtures,and groundwater reservoirs.As a pioneer, Yang et al [24] investigated the free convective flow of a thermally stratified fluid due to a non-isothermal plate.Ishak et al [25] examined the time-independent mixed convection flow through a stable stratified medium near a vertical flat plate.Cheng [26] employed the cubic spline collocation method to explore the impacts of double stratification on natural convection flow of a non-Newtonian fluid near a vertical wavy surface.Mukhopadhyay and Ishak [27]presented a numerical study to discuss the mixed convection axisymmetric flow of thermally stratified viscous fluid over a stretching cylinder.Mishra et al [28] conducted a numerical study to probe a steady flow due to a vertical surface subject to double stratified micropolar fluid.An identical study of nanofluid flow and heat transport with thermal stratification has been presented by Abbasi et al [29], Hayat et al [30],Eswaramoorthi et al [31], Jabeen et al [32].

According to the literature review, there have been several investigations on the flow and transfer of hybrid nanofluids using various mechanisms.In most of these works, the authors have studied the hybrid nanofluid heat transport features by considering constant wall and free stream temperatures and computed single solutions by using analytical and numerical techniques.However, in a variety of real-world circumstances, these temperatures do not remain constant,and we must treat them as a function of space and time variables.Moreover, the energy transport phenomenon of a hybrid nanofluid driven by a stretching/shrinking geometry with variable temperatures has various realistic industrial,engineering and biomechanical applications, like polymer technology, blood flow, treatment of several diseases,metallurgical processes, and annealing and thinning of wires.As a result of these applications in numerous disciplines of science and technology, it is worthwhile to discuss and explore thermally stratified flow of hybrid nanofluids with thermal radiation.

As per the authors’ knowledge and based on the open literature review,it is noticed that multiple numerical solutions for thermally stratified flow of hybrid nanofluid driven by a shrinking surface have not been reported yet.The core novelty of the current study is to perform a numerical simulation to predict the multiple solutions for thermally stratified flow of Cu − Al2O3− hybri d nanomaterials along with heat transport analysis in the presence of thermal radiation and slip mechanism.In this research, the authors formulated the problem of two-dimensional time-dependent magnetohydrodynamic(MHD)hybrid nanofluid flow over a flat sheet with the help of conservation laws and Boussinesq approximations in the form of partial differential equations.Moreover, the current investigation employs an efficient and versatile numerical method,the bvp4c routine in MATLAB, to acquire the problem solutions.The outcomes are compared to those obtained without hybrid nanofluids as well as those obtained with nanofluids.

2.Mathematical model

2.1.Formulation and basic equations

As shown in figure 1, the considered physical situation involves unsteady,two-dimensional,and incompressible flow of hybrid nanofluids past a flat surface with stretching and shrinking characteristics.During the mathematical formations, the following assumptions have been kept in mind:

Figure 1.Flow geometry and coordinate system.

Therefore,in view of the Boussinesq approximations and above stated restrictions,the basic conservation equations for MHD unsteady flow of hybrid nanofluids subject to the Tiwari and Das[3]model can be acknowledged as(see Fang et al [33], Rohini et al [34], Devi and Devi [35], Ismail et al [36]):

Continuity equation:

Momentum equation:

Energy equation:

The following boundary conditions are put on the surface walls as well as on the free stream:

Table 1.Thermo-physical characteristics of hybrid nanofluids and base fluid (see Khanafer et al [37], Oztop and Abu-Nada [38]).

In the above equation, the subscript"hnf" represents the hybrid nanofluids, while"f" denotes the base fluid.

The radiative heat fluxqrin equation (3) is depicted as:

where(k*)signifies the Stefan–Boltzmann constant and coefficient of mean absorption.Now, pluggingqrfrom equation (7) into equation (3) , we get

The dimensionless form of the modelled problem is obtained by utilizing the following subsequent transformations:

Substituting equation (9) into equations (2) and (8), the following converted system of ordinary differential equations is derived:

subject to the conditions

The wall mass transfer velocity (or fluid suction velocity)becomes=S,whereSis a constant that specifies the wall mass transfer parameter, withS> 0 indicating suction,S< 0 indicating injection andS= 0 indicating impermeability.

The other dimensionless parameters are written as

In this analysis,the variables of physical curiosity are the drag force and the heat transfer rate,which are written in their dimensionless form as follows:

3.Numerical method

Figure 2.Variation of f″ (0) with magnetic parameter M againstλ.

The governing set of ordinary differential equations (10) and(11) along with boundary conditions (12) is numerically integrated by employing the built-in MATLAB solver bvp4c.The flow and heat transport characteristics are observed by the virtue of non-dimensional velocity, temperature, skinfriction, and heat transfer rate computed by the above-mentioned numerical method.The main purpose of this analysis is to predict multiple solutions for flow fields in the case of a shrinking surface.These solutions can be achieved by two different sets of initial guesses which fulfil the far-field boundary conditions.In this method, higher-order non-linear ODEs are first converted into first-order differential equations by choosing suitable variables.For this,let us assumef=y1,f′ =y2,f″ =y3,θ=y4andθ′ =y5..The new system of first-order differential equations is given by

and the boundary conditions become

An important step in this methodology is to give a suitable finite value to(η= ∞), so that the field boundary conditions(17)are asymptotically satisfied.The tolerance criterion is set as 10−6to gain the correct numerical outcomes.

4.Computed results

4.1.Code validation

To begin, tables 2 and 3 highlight the validity of the current results.In these table,the computed results of the skin-friction coefficient (upper and lower branch solutions) are compared with the published works of Wang et al[39],Mahapatra et al[40], Ismail et al [36] and Bachok et al [41].In this regard,remarkable consistency has been accomplished.This proves that the proposed model and present findings are accurate.

?

?

Figure 3.Variation of −θ′ (0) with radiation parameter R againstλ.

4.2.Discussion

Various numerical and graphical outcomes of the problem have been prepared and presented in detail.The simulated outcomes are displayed through velocity and temperature distributions together with the skin-friction coefficient and heat transfer rate for different flow parameters to better understand the flow and heat transport features, for instance,for magnetic parameterM, stretching/shrinking parameterλ,suction parameterS, velocity slip parameterα1, nanoparticle volume fractionsφ1andφ,2unsteadiness parameterβ, radiation parameterR, thermal stratification parameterst,heat generation/absorption parameterδand Prandtl numberPr.

Figure 4.Variation of f″ (0) with suction parameterS againstλ.

Figure 5.Variation of −θ′ (0) with stratification parameter st againstλ.

Figure 6.Velocity fieldsf ′(η) for distinct nanoparticle volume fractionφ.1

Figure 2 is plotted to witness the impact of magnetic parameterMon wall shear stressf″ (0) against shrinking parameterλ.It is shown that a dual solution exists for shrinking parameter(λ< −1) and suction parameter (S>2)when keeping all other parameters fixed.Further, we observed that these dual solutions occur up to a critical valueλcof the shrinking parameterλ.The critical valueλcplays an important role in predicting the nature of computed solutions.It is seen that a unique solution exists whenλ=λ cand dual solutions are possible forλ<λ cwhile no solution exists whenλ>λ.cThe respective critical valuesλcfor distinct magnetic parametersM=(0 , 0.2, 0.6) areλc= −3.3556,−3.3840 and−3.4019, as displayed in figure 2.A significant decline in the critical value is noted with increasing value of the magnetic parameter which enhances the existence domain of dual solutions.From a physical point of view,this happens because greater values of the magnetic parameter create resistance in the flow, which results in smothering of the velocity,as shown in the first solution.Moreover,the value off″(0) with respect to the first solution is always increasing with increasing magnetic parameter.On the other hand, an opposite behaviour is depicted by the second solutions.The impact of radiation parameterRon heat transfer rate−θ′ (0)against shrinking parameterλis displayed in figure 3.Again,we depict the existence of dual solutions in the case of shrinking flow.The magnitude of the critical value∣λc∣increases with increasing radiation parameterR.The wall shear stressf″ (0) is plotted in figure 4 for distinct values of suction parameterSas a function of shrinking parameterλ.It is worth mentioning that dual solutions exists for the shrinking parameter(λ< 0) within a specific range of the suction parameter,i.e.S= 2.1, 2.3, 2.5.One can clearly see that the local wall shear stressf″ (0) increases with increasing suction parameter for the first solution, while it decreases in the case of the second solution.Figure 5 portrays the variation of local heat transfer rate−θ′ (0) with varying values of thermal stratification parameters.tWe observe that as the stratification parameter increases the rate of heat transfer decreases for both solutions.

Figure 7.Temperature fieldsθ (η) for distinct nanoparticle volume fractionφ.1

The impacts of the volume fraction of copper nanoparticlesφ1on the velocity and temperature distributions for fixed values of other parameter are displayed in figures 6 and 7.The outcomes uncovered that the variation in nanoparticle volume fraction has less effect on velocity fields.However,figure 6 suggests that the fluid velocity shows a decreasing behaviour with higher nanoparticle volume fraction for both solutions.On the other hand, it is observed through figure 7 that a substantial rise in temperature distribution is noted for greater volume fraction.The velocity distributions off′(η)for various values of suction parameterSare plotted in figure 8.The plotted graphs show the existence of dual solutions for shrinking flow inside the boundary layer region.As expected, an enhancement in velocity profiles is seen for greaterSin the upper solutions, while the inverse is noted in the lower solutions.Figure 9 portrays the relationship between the dimensionless temperature profilesθ(η) and the nanoparticle volume fractionφ2within the boundary layer.It should be noted that both the temperature profiles and associated boundary layer thickness increase with growing values ofφ2.The velocity distributionsf′(η) for varyingMare illustrated in figure 10, which shows that as the magnetic parameterMincreases the velocity profiles show a decreasing trend for the second solution and the opposite is noted for the first solutions.The dual temperature profilesθ(η) for precise entries of the radiation parameterRare presented in figure 11.With an increment in radiation parameter, both solutions showed that the temperature distribution increases.On the other hand, the corresponding thermal boundary layer thickness is higher for the second solution.Figure 12 depicts the effect of velocity-slip parameterα1on the dual velocity profilesf′(η) inside the boundary layer region.Once again,an opposite behaviour is displayed by the first and second branch solutions within the boundary layer region.The larger values ofα1cause a reduction in fluid velocity for the first branch and enhance the fluid velocity in the second branch.Finally, figure 13 is outlined for the behaviour temperature distributionθ(η) under the influence of stratification parameters.tThe hybrid nanofluid temperature is seen to rise for increasing values of the thermal slip parameter.

Figure 8.Velocity fieldsf ′(η) for distinct suction parameterS.

Figure 9.Temperature fieldsθ (η) for distinct nanoparticle volume fractionφ.2

Figure 10.Velocity fieldsf ′(η) for distinct magnetic parameter M.

Figure 11.Temperature fieldsθ (η) for distinct radiation parameter R.

Figure 12.Velocity fieldsf ′(η) for distinct velocity-slip parameter α.1

Figure 13.Temperature fieldsθ (η) for distinct stratification parameter st.

5.Main findings

Numerical simulations for an unsteady thermally stratified flow of MHD hybrid nanofluid across a stretching/shrinking surface with thermal radiation and slip mechanism were carried out in this study.The main feature of the current study was the depiction of multiple branches,namely the upper and lower branch,for flow and temperature fields in the case of a shrinking surface.The following major outcomes can be summarized.

1.As the suction parameter was increased, the existence domain of the dual solution was increased with higher critical values of the shrinking parameter.

2.A higher skin-friction coefficient was noted for larger values of the magnetic parameter in the upper branch.

3.At higher values of the stratification parameter, a substantial rise in fluid temperature was observed for both solutions.

4.Higher values of the radiation and thermal stratification parameters decreased the Nusselt number for both the upper and lower branch solutions.

5.A decreasing tendency was observed for velocity curves with increased values of the velocity-slip parameter in the case of the second solution.

6.The hybrid nanofluid temperature was significantly increased by a greater thermal radiation parameter in both solutions.

ORCID iDs