Abstract

In this exploration, we decided to investigate the significance of prescribed thermal conditions on unsteady 3D dynamics of water-based radiative hybrid nanofluid with the impact of cylindrical-shaped nanosized particles (alumina () and titania ()). For physical relevancy, the impact of the Lorentz force is also included. The combination of suitable variables has been used to transform the transport equations into the system of ordinary differential equations and then numerically solved via the Keller-Box approach. Graphical illustrations have been used to predict the impact of the involved parameters on the thermal setup. Convergence analysis is presented via the grid independence approach. Skin frictions and local Nusselt numbers against various choices of involved parameters are plotted and arranged in tabular forms. It is observed through the present investigation that temperature distribution is increased with the higher choices of radiation parameter and decreased with the improvement in the choices of temperature maintaining indices (i.e., ). Moreover, the thermophysical properties except specific heat for hybrid nanofluid are improved with the involvement of cylindrical-shaped nanoparticles. The temperature of the hybrid nanofluid is observed to be higher for variable thermal conditions as compared to uniform thermal conditions. Outfalls for a limited version of the report have been compared with a previous published paper.

1. Introduction

Nanofluids have been widely used in many technological and industrial processes like nuclear reactors, automobile radiators, and solar aircrafts because their superior thermal conductivity as compared to conventional fluids as the rate of heating/cooling is extremely reliant on the performance of thermal conductivity of the nanosized particles. Various mathematical relations [1, 2] have been adopted to investigate the thermophysical properties of nanosized particles but the most appropriate mathematical relations regarding these properties are developed/discussed by Masoumi et al. [3]. The boundary regime flow of nanofluid due to an expanding device is firstly numerically deliberated by Khan and Pop [4] with random motion and thermodiffusion effects of tiny-sized particles. Sheikholeslami and Rokni et al. [5] elaborated the heat transfer mechanism by using single-phase and double-phase estimations of magneto-nanofluid and proved that higher estimation of tiny particle concentration augments the temperature gradient. Ganvir et al. [6] summarized the combined performance of convective heat transference, particle size, thermophysical properties, inlet velocity, volume concentration, and liquid temperature, both analytically and numerically.

An improved approach to overcome the thermal needs of industrial and engineering processes is famous with the name of hybrid nanofluid. The investigation about the importance of hybrid nanofluid has gained a tremendous height due to their wide applications in domestic refrigerators, engine cooling, microelectronics, heat exchanger devices, pharmaceutical processes, fuel cells, grinding processes, ultrasonic radiations, thermal diffusion processes, and many more. Hybrid nanofluid can be formed by submerging two or more nanoparticles into the host liquid. These nanoparticles may include metals, dielectrics, liquid materials, polymeric, lipids, and semiconductors. Some common examples of nanoparticles are copper, zinc-oxide, carbon nanotubes, phosphates, zinc sulfide, cadmium telluride, etc. Alumina and titania nanoparticles are widely used in engineering applications, and these are prepared from metal precursors. These tiny-sized particles can be synthesized by electrochemical, chemical, or photochemical methods. These nanoparticles have commercial applications due their structure, high strength, electron affinity, and electrical conductivity. The most widely used oxide ceramic material is alumina, and it has applications in cutting tools, tap washers, spark plugs, etc. The most abundantly used nanomaterial for human life is titania, and it is used in sunscreen, biomedical applications, photovoltaic devices, pharmaceutical drugs, and waste water treatment and as a food additive. Sarkar et al. [7] reviewed the advantages and disadvantages of nano as well as hybrid nanofluids and recommended that hybrid nanofluids have various advantages as compared to conventional nanofluids because of their auspicious heat transference enhancement, pressure drop ability, improved thermal network, and favorable aspect ratio. Sidik et al. [8] disclosed the recent progress related to the field of hybrid nanofluids by discussing the factors affecting their thermal properties and stability. Unsteady dynamics of hybrid nanofluid with heat transference characteristic due to the expanding/contracting device is numerically explored by Waini et al. [9] and presented the stability analysis regarding the dual solutions. Numerical evaluation regarding pure water-based hybrid nanofluids with the combination of nanoparticles (alumina, titania, and silica) is explored by Minea [10]. Investigation concerning water-based hybrid nanofluids with the optical properties of titania and alumina nanoparticles is made by Leong et al. [11] along the stability of the obtained solution. Moldoveanu et al. [12] glimpsed the experimental evaluation regarding the hybrid mixture of alumina and titania nanoparticles with viscosity as the main focus of the investigation. Moldoveanu et al. [13] also evaluated the hybrid mixture of titania and alumina nanoparticles with thermal conductivity as the foremost point of the exploration. Shirazi et al. [14] experimentally discussed the mixture of titania nanoparticles and water in order to observe the level of the oil recovery enhancement process. Khan et al. [15] elaborated the mixed convective stagnation dynamics of the radioactive mixture of titania, copper, and water in the magnetic environment towards an expandable device. An experimental evaluation regarding the stability, characterization, and dynamic viscosity of the hybrid mixture of titania and cupric oxide with water as working fluid is completed by Asadi et al. [16]. Ahmad et al. [17] explained the heat/mass transference attributes in the hybrid interpretation of alumina and copper nanoparticles through porous media. A comparative depiction regarding unsteady transport of magnetically influenced water-based fluid with the hybrid mixture of nanoparticle combinations (copper-alumina and alumina-titania) towards an expanding device using finite element approach is elaborated by Ali et al. [18]. Some recent scientific contributions about hybrid nanofluids have been addressed by some scholars/researchers (refs. [1921]) and their applications (refs. [2226]).

In the last few years, researchers are interested to discuss the novel impact of the shape of nanoparticles in the improvement of heating/cooling processes, solar aircrafts, effective thermal conductivity, and other thermophysical characteristics. Zhang et al. [27] measured the values of effective thermal conductivity and thermal diffusivity of nanofluids by considering the spherical as well as cylindrical-shaped nanoparticles using transient technique and predicted that Hamilton and Crosser models provide the best approximation for these thermophysical characteristics. Later on, Timofeeva et al. [28] theoretically explained the effect of the nanoparticle shape on thermophysical behaviour of alumina with the help of experimental data. Yang and Ma [29] provided the computer simulation for the understanding of translocation processes of nanoparticles with the usages of different shapes (ellipsoids, discs, rods, and spheres) of nanoparticles across a lipid bilayer. Maheshwary et al. [30] experimentally discussed the significance of the particle shape, particle size, and concentration on thermal conductivity of water-conveying titania nanofluid and predicted that thermal conductivity of the mixture is improved by intensifying the choices of the particle shape, size, and concentration. Sheikholeslami [31] discussed the effect of various shapes (platelet, brick, cylinder, and sphere) of nanoparticles on the forced convective flow of water-conveying cupric oxide nanofluid within a permeable lid-driven enclosure in the magnetic environment by using CVFEM and indicated that the Nusselt number declines with the augments of the Lorentz force. Rashid and Liang [32] numerically investigated the implication of the nanoparticle shape (sphere and lamina) on the dynamics of magnetized nanofluid with heat transference phenomenon, thermal radiation effect, and joule heating process towards an expanding disk through a porous space in a rotating frame. Dinarvand and Rostami [33] analyzed the shape factor influence of nanoparticles (graphene oxide and magnetite) for bidirectional unsteady dynamics of water-conveying hybrid nanofluid squeezed between two parallel surfaces using the Tiwari-Das model and predicted that the shape factor effect of nanoparticles has a crucial role in food processing, polymer processing, injection modeling, lubrication, etc. Bhattad and Sarkar [34] theoretically examined the significances of the nanoparticle size and shapes (brick, platelet, sphere, and cylinder) on the thermohydraulic enactment of a sheet evaporator by using hybrid nanofluids having various combinations of mixtures.

The combined significances of thermal radiation and Lorentz force space contribute a vital role in the development of combustion processes, nuclear weapons, electron ramifications, polarization process, stellar evolution, heat conduction process, petroleum reservoirs, etc. Turkyilmazoglu and Pop [35] numerically addressed the natural convection dynamics of radiative water-conveying nanofluid containing copper, alumina, titania, cupric oxide, and silver nanoparticles across a flat device. Devi and Devi [36] numerically examined the effects of Newtonian heating and Lorentz force on bidirectional dynamics of water-conveying hybrid nanofluid with the mixture of copper and alumina. From this exploration, it is conveyed that hybrid nanofluid provides the better rate of heat transference than conventional nanofluid. Sheikholeslami and Sadoughi [37] numerically disclosed the MHD effect on the flow of nanofluid inside a porous enclosure with four square heat sources by considering the importance of nanoparticle shapes. Hayat et al. [38] discussed the radiative and heat transfer characteristics for hybrid mixture of silver and cupric oxide nanoparticles with water as base liquid in a rotating frame. It is deduced that rotation and radiation phenomena boost the thermal environment of the hybrid mixture. Some more recent exploration regarding the implications of the Lorentz force and porous media is found in the refs. [3943].

The variation in the temperature fluctuation at the geometric surfaces is beneficial for several industrial and engineering applications. Liu and Andersson [44] implemented the variable thermal conditions to investigate the heat transference characteristics for 3D dynamics of a liquid towards a bidirectional expanding device. Both the heating processes, namely, PST (prescribed surface temperature) and PHF (prescribed heat flux), have been discussed by Liu and Andersson, and it was predicted that variable thermal conditions provide an improved rate of heat transference than arbitrary thermal conditions at the geometric surface. Oliveira et al. [45] discussed the practical applications of variable thermal conditions in the engineering processes like power converters, motor controllers, passive thermosyphons, and air conditioning process. Waini et al. [46] discussed the heat transference characteristic for the steady dynamics of hybrid nanofluid past a vertical thin needle by considering the variable heat flux at the geometric surface. Waini et al. [47] also investigated the heat transference process along with variable heat flux for the stagnated dynamics of hybrid nanofluid (alumina and copper) with water as working liquid on a contracting cylindrical geometry. Some more recent contributions related to the variable temperature of the geometrical surfaces are found in the refs. [4852].

In the view of abovementioned comprehensive literature survey, it is noticed that much attention has not been given to the dynamics of hybrid nanofluid towards bidirectional elongating geometry. The main theme of the present contribution is to predict the effects of cylindrical-shaped nanoparticles (alumina and titania) for radiative water-conveying hybrid nanofluid flow towards an unsteady bidirectional elongating device with prescribed thermal conditions, and this type of contribution is not found in literature to the best of the author’s knowledge. Additionally, influence of the Lorentz force is also incorporated in the mathematical model. Suitable mathematical relations have been used to transform the transport equations into dimensionless forms, and then, computer simulation is made via the Keller-Box method [5359]. Finally, the foremost outcomes obtained through present numerical investigation are presented through various plots and tables.

2. Mathematical Formulation

In order to frame the unsteady mathematical model for bidirectional dynamics of water-conveying hybrid nanofluid with cylindrical-shaped nanoparticles (alumina () and titania ()), the Cartesian configuration is adopted. The mathematical relation of the Lorentz force is used to inspect the MHD (magnetohydrodynamics) effects with strength ( represents the initial strength of the magnetic field). The mathematical relation of Rosseland approximation is then followed to examine the effects of thermal radiation. The tiny particles are considered in thermal equilibrium. The no-slip phenomenon is considered at the surface to keep the flow incompressible as well as laminar. Expansion velocity is opted along the -axis, and the expansion velocity is decided along the -axis, whereas is the region covered by the hybrid nanofluid (as sketched in Figure 1). In order to provide the variable temperature mechanism at the surface, two types of thermal conditions, namely, PST and PHF, are applied. Table 1 is constructed to summarize the thermophysical characteristics of water , alumina, and titania.

In the continuation of the abovementioned assumptions with a boundary layer theory, the transport equations are manifested as follows (refs. [38, 52]):

The velocity and thermal conditions for equation (1) are conveyed as follows (ref. [41]):

Here, designates the velocity components along the -, -, and -directions, respectively, shows the temperature at the surface, time factor is symbolized by , are the indices that are used to control the temperature at the surface, and are dimensional constants, is opted to describe the effective viscosity of the hybrid mixture, is chosen to label the density of the hybrid mixture, is taken to mark the thermal conductivity of the hybrid mixture, is picked to state the thermal diffusivity of the hybrid mixture, is selected to represent the specific heat capacity, and is typified to express the influence of electrical conductivity of the hybrid nanofluid.

The mathematical relations to introduce the cylindrical shaped nanoparticles for present evolution of hybrid nanomaterial are composed as follows (refs. [2, 27, 31]):

Here, quantities of volume fractions for alumina and titania nanoparticles are expressed through and , respectively. The case of conventional fluid can be recovered by considering .

The equation for radiative heat transfer is defined as follows (refs. [35, 38]):

Here, describes the radiative heat transference, illustrates the Stefan Boltzmann factor, and explains the effect of the mean absorption factor.

The set of relations used to nondimensionalize the present mathematical model is conveyed as follows (refs. [41, 52]):

With the involvement of equations (6) and (7), the transport equations become with boundary restrictions

Here, the Hartmann number is recognized by , the unsteady factor is stated by , the elongation ratio is expressed by , the Prandtl factor is indicated by , is the radiation factor, and are the relations for the present hybrid mixture and these relations are elaborated as follows:

The most fascinating quantities for thermal processes and most used in the improvement of heat exchanger devices are termed as skin-friction coefficients (i.e., and ) and the local Nusselt number (i.e., ). The mathematical formulations of these quantities are communicated as follows (refs. [36, 44, 51]):

The dimensionless formulations of the abovementioned quantities using Reynold’s numbers are transported as follows:

3. Keller-Box Simulation

The final form of the system of equations obtained through the aforementioned modeling along with boundary conditions is simulated via the Keller-Box approach. This numerical approach has accuracy of up to second order and has rapid convergence ability than other routine work numerical approaches (shooting method, RK-method, BVP4c, etc.). This approach also provides the flexibility about adoption of the step size for the computational domain and is more appropriate for the solution of boundary layer flow problems. The major advantage of this method over other numerical methods is its unique conversion procedure of differential equations into algebraic equations using central difference approximations. Foremost steps to implement this numerical approach are stated below and summarized via the flow chart (Figure 2): (i)The first step is to transform the higher-order differential equations into first-order differential equations(ii)The next step is to transmute the obtained differential system into the difference equation system via central difference numerical approximations(iii)Linearization of the system of difference equations is completed with the courtesy of Newton’s linearization standard method(iv)The linearized equations are then arranged into matrix-vector forms(v)The LU decomposition technique is opted to solve the obtained matrix-vector problem(vi)Finally, the value of the unknown vector provides the numerical solution of the aforementioned mathematical problem

During the implementation of the abovementioned steps, the computational domain is truncated into the finite domain . In order to obtain the first approximation of the numerical solution, we selected and then, the desired accuracy, i.e., , is achieved by varying the value of (the numbers of grid points) with the reduction in the value of (step size).

Table 2 is designed to estimate the rate of convergence of the Keller-Box simulation as well as to find the best choice of for the simulation of the local Nusselt number and skin-friction coefficients. It is deduced through Table 2 that one thousand grid points are enough for the convergent solution of and five hundred grid points are sufficient to achieve the convergence criteria for , whereas five thousand grid points are necessary to attain the convergent approximation for both as well as . In order to check the stability of the Keller-Box solution, the value of is increased up to ten thousand and the solution is found consistent. The convergent solution obtained through Table 2 is used for further manipulations in order to find the impact of involved parameters on the thermal setup, local Nusselt number, and skin-friction coefficients.

4. Code Validation

In order to validate the numeric code for the solution of the considered problem, the outfalls reported in the present contribution have been compared with the outfalls discussed in the article (Liu and Andersson [44]) for and in the absence of nanoparticles. A convincing scientific connection has been found between the present scrutiny and the published activity. In this regard, Table 3 is arranged in the analysis.

5. Results and Discussion

The present heading describes the importance of the involved important parameters like radiation factor and temperature maintaining indices on thermal setups via various plots, whereas the physical quantities like the local Nusselt number and skin friction coefficients are discussed for various estimations of involved constraints via various graphs and tables. Moreover, the thermophysical properties for the present hybrid mixture are also computed at the end of this section and discussed deeply. Figure 3 describes the impact of radiation factor on and on . It is detected through Figures 3(a) and 3(b) that the escalating choice of enhances the worth of thermal setups. The thickness of the thermal layer is observed larger for smaller choices of as compared to larger estimations of . The maximum temperature for the PST mechanism is observed as one, whereas the maximum temperature for the PHF case is noticed as 0.5 when is increased from 0 to 2.0. Overall, the strength of the thermal setup for the PHF case is observed to be more prominent than that for the PST case. Physically, the abovementioned changes are produced because is the mathematical ratio of Stefan Boltzman number to mean absorption factor . The value of is increased whereas the value of is diminished with the positive tendency of . Figure 4 explains the influence of temperature maintaining index corresponding to the -direction on and on while other parameters are being fixed. The temperature of hybrid nanofluid is diminished for both the PST and PHF cases with the escalating amount of . The value of temperature is detected higher for the PHF case than the PST case for and is deduced higher for the PST case than the PHF case for other choices of . Physically, the temperature distribution is well dominant for variable thermal conditions than arbitrary thermal conditions. The wideness of the thermal layer is larger in the PHF case than the PST case. Figure 5 is plotted to explain the worth of temperature maintaining index corresponding to the -direction on and on by keeping the other involved factors fixed. The temperature with the higher estimations of is reduced for both the PHF and PST mechanisms. The value of temperature fluctuation is attained higher for PST conditions than PHF conditions for all the choices of , but the thermal thickness is achieved quite dominant for the PHF case as associated with the PST case.

Figure 6(a) manifests the combined influence of and on for fixed values of other involved parameters. Skin-friction coefficient is reduced with increasing amounts of both and . Physically, more electric conduction is produced in the flow of the hybrid mixture with the positive growth in from 0.5 to 3.5. Skin-friction coefficient is reduced with the phenomenon of electric conduction. Moreover, the electrical conductivity of alumina is much higher than the electrical conductivity of host liquid and this phenomenon produces the reduction in the value of . The smaller selection of produces the stream of moderate thickness, whereas the thickness of the stream is observed to be double when reached 0.10. Figure 6(b) explains the combined influence of and on for fixed selections of other involved constraints. Skin-friction coefficient is reduced with the higher estimation of . Physically, expansion rate is reduced with the improvement in the value of , and therefore, reduction in the computation of is attained. Figure 7(a) reveals the collective effect of and on with other parameters being kept fixed. Skin-friction coefficient is abridged with growing amounts of both and . Actually, a less elongation rate is produced with the development in from 0.2 to 1.4. Skin-friction coefficient is diminished with the productions of a less elongation rate. Moreover, the electrical conductivity of titania is much higher than the electrical conductivity of host liquid and this phenomenon produces the reduction in the computation of . As electrical conductivity of titania is higher than the electrical conductivity of alumina, therefore, faster reduction in is achieved than . The smaller assortment of yields the stream of minute thickness, whereas the thickness of the stream is observed triple when reached 0.10. Figure 7(b) elucidates the joined impact of and on for stationary values of other mathematical constraints. Skin-friction coefficient is enhanced with the advancement in . Precisely, expansion rate is enhanced and expansion rate is diminished with the development in , and therefore, enhancement in the computation of is attained.

Figure 8(a) explores the graphical assessment of Nusselt number against the variations in and with other parameters being retained as fixed. Nusselt number is improved with the enhancement in from −2.0 to 2.0, and also, it is increased with the escalation in from 0.0 to 2.0. Energy is transformed in the form of electromagnetic waves with the involvement of thermal radiation, and it is more efficient in porous space or vacuum because it is the ideal situation for full transmission of the radiation energy. The value of is augmented with the escalation in from 0.0 to 2.0. The rate of heat transference is observed slower for negative values of than positive values of because the surface temperature is mentioned higher for the negative value of than the positive value of , and therefore, the thermal flux across the surface will be higher for the positive value of than its negative value. As a result, the Nusselt number is tremendously improved.

Figure 8(b) describes the graphical evaluation of Nusselt number against the variations in and with other parameters being taken as fixed. Nusselt number is improved with the enhancement in from −2.0 to 2.0, and also, it is increased with the escalation in from 0.0 to 2.0. The transmission of heat is perceived slower for negative selections of than positive selections of because the stretching device is maintained at higher temperature for the negative value of than the positive value of , and therefore, thermal flux across the stretching device will be higher for the positive value of than its negative value. As a result, the Nusselt number is extremely enriched. Physically, unsteady expansion parameter is involved in equation (9) with the sum of dimensionless temperature and its derivative with respect to . As a result, the rate of heat transference across the plane is increased.

Roles of solid volume fractions for alumina and for titania on the local Nusselt number as well as on skin friction coefficients are discussed in Table 4. It is attained in Table 4 that escalating selections of and improve the skin-friction coefficients for the present model. It is also observed that the impact of and is more dominant for the flow along the -direction as compared to the -direction. Mathematically, the involvement of is absent in equation (14), whereas it is involved in the reciprocal form as described in equation (15). As in our study, the value of is selected to be 0.5, so the numerical value of is observed to be higher than the numerical value of . Moreover, the roles of and on the local Nusselt number are also discussed in Table 3 and it is deduced that higher estimations of and in the range improve the rate of heat transference which is beneficial for many industrial and engineering applications. Furthermore, the rate of heat transference is observed to be equal quantitatively for both the PST and PHF mechanisms. The thermophysical features for the present hybrid mixture are discussed in Table 5 for various choices of (the weightage of titania) by adjusting . The escalating estimation of enhances the density, thermal conductivity, and electrical conductivity of the hybrid mixture whereas the heat capacity is reduced with the present development in . Physically, the mixture of alumina and titania with the working fluid, i.e., water, provides the improvement in the thermophysical features as compared to that in Table 1. In Table 6, thermophysical characteristics for present hybrid combination are computed for different selections of (the solid volume fraction of alumina) by taking . With the increase of , the density, the electrical conductivity, and the thermal conductivity are improved but the heat capacity is condensed with this variation. The base fluid, i.e., water, delivers enlargement in thermophysical properties of hybrid nanofluid associated to Table 1 with the chemically mixture of alumina and titania.

6. Conclusions

This study provides the mathematical analysis for unsteady bidirectional dynamics of radiative water-conveying hybrid nanofluid (i.e., the combination of titania and alumina) with the significance of variable thermal conditions (PST and PHF) and impact of the cylindrical shape of nanoparticles. The influence of the Lorentz force is also incorporated to make the investigation more impactful. Numerical simulation is made via the Keller-Box approach, and key observations are listed as follows: (i)The thermophysical features of the hybrid nanofluid except specific heat are improved with the involvement of the cylindrical shape of nanoparticles, i.e., titania and alumina(ii)The rate of heat transference is observed identical for both the PST and PHF cases(iii)The magnitude of the stress applied on the -direction is observed to be higher than the magnitude of the stress applied along the -direction with the positive estimations of volume fractions of nanoparticles(iv)The temperature of the hybrid nanofluid is increased with the development in the choice of and decreased with the improvement in the choices of temperature maintaining indices (v)Skin-friction coefficients are reduced with progressions in the values of the Hartmann number, unsteady parameter, and volume percentages of nanoparticles(vi)The local Nusselt number is improved with developments in the amounts of radiation and unsteady parameters

This scientific contribution has many mechanical, biomedical, and commercial applications. These are coating a sheet with hybrid nanomaterials, manufacturing of printing ink, degrading organic contaminants, manufacturing of sodium vapour lamps, etc. This study is also really helpful for the researchers working in the field of nanomaterials and can be protracted in the future by considering different geometries.

Nomenclature
:Strength of magnetic field
:Stretching rates
:Time
:Time coefficient
:Space coordinates
:Stretching velocities
PST:Prescribed surface temperature
PHF:Prescribed heat flux
:Surface temperature
:Thermal indices
:Dimensional constants
:Ambient temperature
:Velocity field components
:Temperature
:Hybrid nanofluid
:Dynamic viscosity of
:Density of
:Electrical conductivity of
:Thermal conductivity of
:Thermal diffusivity of
:Specific heat capacity of
:Volume fraction of alumina
:Volume fraction of titania
:Fluid
:Base fluid
:Alumina nanoparticles
:Titania nanoparticles
:Radiative heat transfer
:Stefan Boltzmann constant
:Mean absorption coefficient
:Similarity variable
:Dimensionless velocities
:Dimensionless temperatures
:Prandtl number
:Radiation parameter
:Unsteady parameter
:Hartmann number
:Dimensionless quantities
:Reynolds numbers
:Nusselt number
:Skin friction coefficients
:Stretching ratio parameter
:Thermophysical properties
:Step size
:Convergence criterion
:Numbers of grid points.

Data Availability

The raw data supporting the conclusion of this report will be made available by the corresponding author without undue reservation.

Conflicts of Interest

The authors declare that they have no competing interests.

Authors’ Contributions

All authors contributed equally to this work, and all the authors have read and approved the final version of the report.