List of symbols:
1.
Introduction
Nanofluid is an engineered colloidal mixture regarding base fluid with single nanoparticles. Such behavior of fluid is applicable in cooling process, thermal process, thermal enhancement, thermal reduction and microelectronics. Deep exploration regarding case of mixture of base fluid with two kind's nanoparticles is termed as hybrid nanofluid. Moreover, Casson fluid behaves like a non-Newtonian liquid and shear thinning. Such type of non-Newtonian rheology is applicable in petroleum products, syrup drugs, production of plastic materials and engineering field etc. Rheology related to Casson fluid was investigated by Casson [1]. He has investigated that Casson model displays role like a shear thinning, plastic liquid model and high shear viscosity. Chabani et al. [2] discussed features of hybrid nanoparticles involving Lorentz force in triangular enclosure. They have studied convective flow considering the Darcy Forchheirmer model and included that flow is reduced when Lorentz force is implemented. Chu et al. [3] investigated a study related to hybrid nanofluid in the presence of magnetohydrodynamics for unsteady flow in two parallel plates including shape effects. They considered engine oil as base fluid and model is numerical solved by numerical scheme. Selimefendigil et al. [4] studied the role of various shapes regarding nanoparticles in ventilated cavities considering artificial neural networks. They implemented a finite element approach to find numerical results. It was included that maximum temperature was obtained using L-shaped cavity. Saleem et al. [5] estimated performance of several shapes of nanoparticles in horizontal surface considering base fluid (water). They have used a shooting approach to find parametric study on temperature and flow profiles. It was established that mass transport rate is enhanced versus erratic motion. Saleem et al. [6] established modeling related to Eyring-Powell material considering magnetized effect in mass species and heat transfer involving chemical reaction using finite element approach. Algehyne et al. [7] studied thermal features of Maxwell fluid in hybrid particles over a stretching frame. They have utilized finite element approach to find computational study and parametric study. It was investigated that highest thermal performance can be achieved by hybrid nanofluid rather than nanofluid. Imran et al. [8] performed thermal aspects of the solar collector using approach of hybrid nanoparticles via several shape effects implementing a numerical approach. Several studies regarding hybrid nanofluid are mentioned in [9]. Khan et al. [10] discussed features of Lorentz force in Casson fluid inserting dust particles along with hybrid nanoparticles past a stretchable sheet using an analytic approach.
Ternary hybrid nanofluid is a mixture of three kinds of nanoparticles in base fluid. It is an efficient approach to achieve the highest temperature rather than nanofluid and hybrid nanofluid. The process regarding tri-hybrid nanoparticles is most efficient for development of thermal enhancement and thermal reduction. Sarada et a. [11] estimated enhancement process in energy efficient using applications of ternary hybrid nanoparticles in curved stretching sheet implementing shooting approach. Nazir et al. [12] visualized comparison enhancement into different fluids using Hall currents inserting mixture of multiple nanoparticles in base fluid via finite element methodology. Dezfulizadeh et al. [13] discussed impacts of Lorentz force in MHD flow using exergy an efficiency approach in the presence of ternary hybrid nanofluid. Oke [14] determined performance of tri-hybrid nanofluid in mass dissuasion and thermal energy considering base fluid (EG) past a rotating surface. Xiu et al. [15] studied dynamic behavior of ternary-hybrid nanoparticles in mass diffusion and thermal field considering small and large volume of nanoparticles with platelet, spherical and cylindrical nanoparticles over wedge using shooting approach. They have found that tri-hybrid nanofluid was observed very significantly for thermal performance rather than hybrid nanofluid. Animasaun and Asogwa [16] discussed comparative influences among alumina and water-based nanoparticles including based cupric nanoparticles (water-based) over stretching surface. Rasool et al. [17] experienced multiple aspects based on radiation and viscous dissipation using theory of non-Fourier's in second grade rheology using analytic approach on porous surface. Ashraf et al. [18] captured thermal features of magnetohydrodynamic flow in peristaltic fluid inserting magnetite nanoparticles in blood using the D-quadrature algorithm.
Nowadays, several industrial applications in pigments, electronics, cosmetics, food processing and utilizing engineering process metal oxide nanomaterials such as silicon dioxide, aluminum oxide and titanium dioxide. The retention and transport of nanoparticles (Al2O3,SiO2,TiO2) ware utilized to resolve issues of petroleum engineering and environmental. Nanoparticles (Al2O3,SiO2,TiO2) were introduced for enhancement of oil recovery by Hendraningrat et al. [19]. Al2O3 nanoparticles have higher boiling and melting points which are useful for enhancement of thermal inertia (see Song et al. [20]). In addition, Al2O3 nanoparticles are applicable to enhance the strengthens of smoothness, ceramics, creep resistance and crack hardness. Dynamic features of thermal transport in the presence of partially ionized liquid subjected to magnetic field is significantly exaggerated by magnetic field. Animasaun et al. [21] investigated that the Hall effect occurs when a magnetic field is applied perpendicular to the current flow through a thin sheet. As a result, an electric field is created parallel to the magnetic and current fields and directly proportional to the product of the magnetic induction and current density. Animasaun et al. [22] discussed comparison impacts among 47 nm and 36 nm alumina nanoparticles in base fluid (water) under action of Hall currents. Farooq et al. [23] developed model of entropy production rate in mass diffusion and thermal field under effects of Hall currents, viscous dissipation, magnetic field and Joule heating in cylindrical tubules.
The purpose of the present flow problem is to study characterizations of Casson fluid involving Lorentz force and thermal radiation (non-linear) over a 3D vertical surface accompanied with buoyancy force, Forchheirmer porous model and Joule heating. In additionally, mixture of three kinds of nanoparticles (Al2O3,SiO2,TiO2) in ethylene glycol is inserted to visualize thermal performance among cylindrical and platelet nanoparticles. The current problem fills the gap while not adequately studied with the following points.
● Analysis of heat transfer in Casson liquid under effects of Joule heating and viscous dissipation, thermal radiation (non-linear), Hall and ion slip forces using Darcy's Forchheirmer theory past a 3D surface has not been investigated yet.
● Previous to this investigation, no comparative analysis on Casson fluid among cylindrical and platelet nanoparticles over 3D surface using Forchheirmer porous model has been existed in literature works computed by finite element method.
● In existence literature, no investigations have been reported including ternary hybrid nanoparticles and shape effects rather than [24] to visualize comparative analysis among cylindrical and platelet nanoparticles.
Existence works demonstrate that development of Casson liquid inserting ternary hybrid nanofluid using thermal properties of various shapes of nanoparticles over 3-dimensional vertical using Forchheirmer porous model frame is not studied yet. Moreover, ion slip and Hall currents are considered with viscous dissipation and variable thermal conductivity. Present complex development is numerically handled by a finite element approach. In addition, section Ⅰ, section Ⅱ, section Ⅳ and section Ⅴ are based on literature review, problem modeling, numerical approach, results and conclusions, respectively.
2.
Modeling of developed model
A 3D model of Casson fluid is developed over a stretching surface using theories of generalized Ohm's law and Darcy's Forchheirmer model. Three types of nanoparticles are inserted into EG (ethylene glycol) while various shapes of nanoparticles (cylindrical and platelet) are addressed. Motion of nanoparticles is produced using movement of wall velocity. Buoyancy force is considered whereas viscous dissipation and thermal radiation (non-linear) are added in the energy equation. Suspension of nanoparticles regarding aluminum oxide, titanium oxide and silicon dioxide in working fluid is considered. Reduced form regarding Navier-Stoke equations is obtained and physical model is considered by Figure 1. Thermal radiation in terms of non-linear is utilized to determine thermal fields. Surface is assumed as non-conducting in terms of electrically and non-isothermal. Figure 1 is 3D vertical surface in y-, x- and z-components while non-uniform magnetic field is considered along z-component and gravitational force acts along downward direction. Partially flow over surface is induced by bidirectional movement of wall velocities Uw and Vw at z=0. A 3D surface lies in plane at z=0 and flow occupies region in field at z>0. Wall temperature and ambient temperature are considered as Tw and T∞. Further, it is studied that ethylene glycol is shear thinning [25], non-Newtonian liquid [26] and plasma [26]. Casson fluid model is also known as non-Newtonian liquid and shear thinning [1]. Table 1 reveals thermal properties of various nanoparticles. So, its rheology is prescribed by the Casson fluid stress tensor relation.
Conservation laws for energy equation, momentum equation for incompressible and steady flow [24] using Casson liquid are
Boundary conditions [24] of developing model are defined as
Similarity variables [24] are delivered as
Equations (1)–(4) into dimensionless form [24] using variables similarity are
BCs in term of dimensionless [24] are
Correlations regarding ternary hybrid nanofluid [28] for density, thermal conductivity and viscosity are defined as
Thermal conductivity and viscosity models for cylindrical and platelet (see Figure 2) [28] are defined as
Wall shear stresses due to ternary hybrid nanofluid [24] are derived as
Temperature gradient due to ternary hybrid nanofluid [24] is
3.
Computational approach
Finite element method [6,7] is employed to construct numerical solution of desired problem. It is described under.
Step 1: Firstly, the problem domain is transformed into a finite number of elements and derived linear type polynomial over each element. Weak form can be achieved from strong form using weighted residual function. Weighted residuals [6,7] are derived as
Here, w3,w1,w2,w4 and w5 are weighted functions. Unknown variables and shape functions [6] are defined as
Step 2: Boundary integral vectors, stiffness matrices and force vectors are derived. After it, the global stiffness matrix is derived. In additionally, Picard linearization approach is utilized for achieving linear form of algebraic equations. Derivation of stiffness elements are
Step 3: System of linear equations is solved in form of iteratively under numerically tolerance (10−5). Convergence and error analysis is mentioned below
Step 4: MAPLE 18 is used to design programming of finite element approach. Mesh free study is recorded in Table 2. Flow chart of numerical methodology is addressed by Figure 3. Validation results of present model is verified with already published work [24] considering absence of ternary hybrid nanofluid, power law index number, Darcy's Forchheirmer model and Casson fluid. Validation of numerical results with published work is mentioned in Table 3.
4.
Results and discussion consequences
The parametric study is captured to determine flow behavior and transport of thermal energy into a mixture of ternary hybrid nanoparticles in base and working fluids. The graphical outcomes are displayed in Figures 4–12. It is noticed that ternary hybrid nanofluid is a mixture of aluminum dioxide, titanium dioxide, silicon dioxide in base fluid and working fluid. Hybrid nanofluid is considered a mixture of silicon dioxide, titanium dioxide while nanofluid is a mixture of titanium dioxide in Casson fluid. In additionally, thermal and velocity fields are plotted versus different parameters implementing ranges of 0.5≤β≤4, 0.1≤n≤3, 0.0≤M≤2.0, ,0.0≤βe≤3.0, ,0.0≤βi≤4.0, ,204≤Pr≤208, 0.0≤Ec≤3.0, 0.0≤α1≤1.5, 0.0≤ϕ1≤0.04, 0.0≤λ≤0.5, 0.0≤ϕ2≤0.07, 0.0≤ϕ2≤0.07, 0.0≤ϕ3≤0.06, 0.0≤α1≤2.0, 0.0≤ϵ≤1.5, 0.0≤Fr≤3.0, 0.0≤Nr≤3.0 and 0.0≤θω≤2.0.
Figures 4–7 are plotted to determine the effect of several parameters on flow involving dispersion of ternary hybrid nanoparticles. It is visualized that platelet nanoparticles are represented by solid curves on graphs and dot curves are associated with cylindrical nanoparticles on graphs. Figure 4 clearly reveal that there is reduction into flow (y- and x-directions) when β is changed. The effect of the Casson parameter for cylindrical nanoparticles is less than the effect of the Casson number for platelet nanoparticles. This effect is created using concept of Casson fluid in momentum equations whereas existence of β is formulated considering tensor of Casson fluid in momentum equations. Flow becomes reduced in y- and x-directions because of β is appeared in form of 1β. Therefore, an inverse proportional relation can be investigated versus velocity fields. Consequently, flow slows down in terms of y- and x-directions. Casson fluid behaves as a shear thinning and non-Newtonian liquid. Flow becomes significant for higher values of β rather flow for minimize values of β. Thickness for momentum boundary layers is declined when β is gradually enhanced. It is studied that current analysis is transformed into Newtonian model for β→∞ while other values of β is case of non-Newtonian model. In momentum equation, β occurs in the denominator. Hence, an inclination in β results velocity field is enhanced. Role βe (Hall parameter) is predicted on flow behavior in term of y- and x-directions. An enhancement is studied into flow (x-and y-direction) versus variation in βe. It is dimensionless parameter which is formulated using the theory of generalized Ohm's law in current analysis. Further, flow for platelet nanoparticles is produced significantly rather than flow for cylindrical particles. Figure 5 visualize the role of ion slip number (βi) on velocity fields in terms of y- and x-directions. Similar impact is conducted when βi is increased. It is observed that ethylene glycol was predicted as plasma and its behavior of flow is exposed as a non-uniform magnetic field. An interaction among plasma and magnetic field produces ion slip and Hall currents which is responsible for ion slip and Hall forces. Applied magnetic field is opposite against action of ion slip and Hall forces. Hence, Lorentz force reduces when βe and βi are increased. Physically, magnetic field was placed perpendicular along the z-direction of the surface and against flow. Therefore, flow was reduced by applying higher values of Lorentz force because Lorentz force is addressed as opposite against applied magnetic field. Thickness based on flow curves for the case of platelet nanoparticles is higher than thickness generated by cylindrical nanoparticles. Role of Fr on flow in y- and x-directions is addressed by Figure 7. It was found that velocity distribution in y- and x-directions is decreased in presence of higher numerical values of Fr. It is noticed that Fr is dimensionless parameter which is produced considering Forchheirmer-porous theory. Mathematically, Fr is inversely proportional versus velocity field. Consequently, an inclination in Fr results flow increases in y- and x-directions. Thickness for the momentum layers is declined when Fr is changed. Flow for case of Fr=0 is greater than flow for case Fr≠0.
Figures 8–11 reveal impact of temperature profile carrying three kind nanoparticles having cylindrical and platelet shapes against change in Nr,Fr,βe and βi. Reduction into temperature versus implication of Nr has been estimated by Figure 8. Highest performance of thermal process for case of platelet nanoparticles has been obtained rather than for cylindrical particles. It is addressed that Nr is termed as dimensionless which is produced using the concept of non-uniform thermal radiation in energy equation. It was noticed that thermally energy reduces when Nr is distributed. Physically, heat energy moves away from z=0 along with thermal radiations. Therefore, temperature decreases at walls of plate. Thickness of thermally layers for case of platelet nanoparticles is greater than thermally thickness for case of cylindrical nanoparticles. Observations of Fr on thermally field was investigated by Figure 9 including ternary hybrid nanoparticles. Maximum production of thermal field was estimated versus distribution in Fr. Thermal enhancement for the case Fr=0 is less than thermal enhancement for case of Fr≠0. Thermal production for cylindrical particles in presence of variable thermal conductivity is less than thermal production for platelet particles. Figures 10 and 11 are plotted to estimate distribution into thermal energy using parameters regarding ion slip (βi) and Hall parameters (βe). Reduction into thermal distribution was visualized using higher values of ion slip (βi) and Hall parameters (βe). Ion slip (βi) and Hall (βe) are opposite against Lorentz force. Further, collision rate into nanoparticles is inversely directly proportional to ion slip and Hall forces result in reduction into magnetic intensity. This reduction in magnetic intensity produces declination on thermal profile. Platelet nanoparticles experience significant thermal energy rather than experience of thermal energy for cylindrical nanoparticles. Physically, βi and βe are associated in the denominator of the Joule heating term in the energy equation. Consequent, heat energy generates because of the Joule heating process. Figure 12 reveals comparison analysis among alumina-silica-titania nanoparticles/EG, silica-titania nanoparticles/EG, silica-titania nanoparticles/EG, titania nanoparticles/EG and working fluid. It was observed that alumina-silica-titania nanoparticles/EG experiences higher thermal energy rather than for the cases of silica-titania nanoparticles/EG, silica-titania nanoparticles/EG, titania nanoparticles/EG and working fluid. Therefore, ternary hybrid nanofluid is observed as most significant for producing maximum production of thermal energy as well as reduction in thermal energy.
Table 4 is prepared to determine production of wall shear stresses and thermal rate using change in Hall parameter, magnetic parameter, ion slip parameter, Casson parameter and Fr. It was numerically investigated that production of thermal rate ingresses when Hall parameter, magnetic parameter, ion slip parameter are changed. However, the same behavior is investigated versus β and thermal rate decreases versus change in Fr. Wall shear stresses are inversely proportional against distribution in ion slip parameter, Casson number and Hall parameter. But inverse behavior for case of wall shear stresses is estimated when Lorentz force is changed.
5.
Conclusions
Thermal aspects of Casson fluid containing ternary hybrid nanofluid over stretching sheet utilizing Darcy's Forchheirmer theory are addressed. Various shapes regarding nanoparticles (platelet, cylindrical and spherical) are implemented in EG. Thermal radiation (non-linear) is utilized in the energy equation considering effects of viscous dissipation and variable wall velocity. Performance of several parameters on flow and temperature profile is addressed and numerical study is tackled by FEM. Key findings of the problem are listed below.
Hall and ion slip numbers play a remarkable role in controlling momentum and thermal boundary layers. However, declination into motion has been investigated versus β and Fr.
Temperature decreases against variation in βe and βi but opposite trend is noticed in thermal layers against distribution in Fr.
Ethylene glycol is studied as partially ionized and shear thinning and its interaction with magnetic field produces ion slip and Hall currents. Due to this process, EG experiences ion slip and Hall forces which are against Lorentz force. Therefore, flow of EG was accelerated versus change in βe and βi.
Convergence of current analysis is ensured using finite element methodology.
Maximum production for thermal energy can be achieved rather than nanofluid and hybrid nanofluid. This process is applicable in coolants (in automobiles), fuel dynamics, solar system, microelectronics, temperature reduction, cooling process and pharmaceutical process.
It was included that platelet nanoparticles experience more heat energy rather than thermal energy for cylindrical nanoparticles. Similar behavior was analyzed in velocity fields.
Nusselt number increases versus change in ion slip, Hall and M parameters but thermal rate decreases versus distribution in Fr.
Acknowledgements
This research was supported by the Postdoctoral Research Fellowship Training Program from Khon Kaen University (PD2565-02-05).
Conflicts of interest
The authors declare that they have no conflicts of interest.