Significant Production of Thermal Energy in Partially Ionized Hyperbolic Tangent Material Based on Ternary Hybrid Nanomaterials

: Nanoparticles are frequently used to enhance the thermal performance of numerous materials. This study has many practical applications for activities that have to minimize losses of energy due to several impacts. This study investigates the inclusion of ternary hybrid nanoparticles in a partially ionized hyperbolic tangent liquid passed over a stretched melting surface. The fluid motion equation is presented by considering the rotation effect. The thermal energy expression is derived by the contribution of Joule heat and viscous dissipation. Flow equations were modeled by using the concept of boundary layer theory, which occurs in the form of a coupled system of partial differential equations (PDEs). To reduce the complexity, the derived PDEs (partial differential equations) were transformed into a set of ordinary differential equations (ODEs) by engaging in similarity transformations. Afterwards, the converted ODEs were handled via a finite element procedure. The utilization and effectiveness of the methodology are demonstrated by listing the mesh-free survey and comparative analysis. Several important graphs were prepared to show the contribution of emerging parameters on fluid velocity and temperature profile. The findings show that the finite element method is a powerful tool for handling the complex coupled ordinary differential equation system, arising in fluid mechanics and other related dissipation applications in applied science. Furthermore, enhancements in the Forchheimer parameter and the Weissenberg number are necessary to control the fluid velocity.


Introduction
The modeling of real-world phenomena is a charming field of research for the researchers working on modeling and simulations, engineers, physicists and mathematicians because of their utilization to predict and monitor several mechanisms. Several relations have been proposed to study the characteristics and nature of different materials. The hyperbolic tangent liquid is an important non-Newtonian material with a constitutive relation of for ∞ * = 0, Γ * ̇< 1, we obtain So far, extensive studies have been conducted on this model to monitor the aspects of thermal and mass transport. For instance, Khan et al. [1] studied the involvement of double stratification on the thermally radiated flow of chemically reactive species of hyperbolic tangent liquid in the presence of heat absorption/generation and first-order chemical reaction embedded in a porous stretching sheet. They solved the derived boundary value problem via an analytical approach package in the Mathematica symbolic computational software by using the concept of residual error minimization. They presented the effectiveness of the BVPh2.0 tool by tabulating the error analysis. They found that a decline in fluid velocity against the Weissenberg number and an enhancement in temperature were achieved against the radiation parameter. Rehman et al. [2] modeled the hyperbolic tangent fluid with an incompressibility assumption flowing over a stretched cylinder. They utilized the shooting approach to solve nonlinear-modeled equations. They found that augmenting the values of the Prandtl number and the curvature parameter increased the heat transfer rate but reversed the bearing conducted for the stratification parameter. Khan et al. [3] implemented an analytical technique proposed by Liao for the modeled problem of a thermally stratified, hyperbolic tangent liquid. They plotted the obtained solution against different involved parameters. They recorded the direct relationship between the temperature and radiation parameters, whereas an escalation in the chemical reaction parameter lessened the concentration field. Gharami et al. [4] studied heat and mass transport on a magnetized hyperbolic tangent fluid model passed over a stretched cylinder. They solved the modeled nonlinear flow equations numerically via an explicit finite difference procedure coded in a FORTAN package. They presented the comparative study for the effectiveness of the approach and for the authenticity of the obtained solution by computing the dimensionless stress and heat transfer coefficient. They observed an increase in the fluid velocity for the Grashof number. The magnetohydrodynamic (MHD) slip flow of a hyperbolic tangent model in an inclined channel, obeying peristaltic transport, was explored by Hayat et al. [5]. They used a lubrication approach to formulate the model problem. The perturbation scheme was used to show the bearing of involved parameters on velocity and temperature. They noticed a reduction in fluid velocity against the Froude number. Additionally, they recorded an enhancement in the temperature against the thermal slip parameter but a decline in the heat transfer rate.
Mechanisms of heat transportation and the contribution of nanoparticles in fluid have many applications that cannot be overlooked. The dispersion of nanoparticles is essential to monitor their thermal performance. For instance, Nawaz et al. [6] examined dust particles in a magnetized hyperbolic tangent model passed over a moveable heated sheet. They engaged in a finite element procedure to handle the nonlinear equations. They monitored the decline in the temperature profile against the Prandtl number. Wang and Mujumdar [7] presented a review on the material characteristics of nanofluid theoretically and numerically. They discussed the merits and demerits of different models proposed by several researchers, which deal with thermophysical features. Couette flow and Poiseuille flow of hyperbolic tangent models with Navier slip conditions were studied by Tlau [8]. They handled the modeled equations analytically. They recorded the decline in dimensionless stress against the Weissenberg parameter. Kebede et al. [9] studied unsteady, buoyancy-driven nanofluid flow passed over a moving wedge, obeying dissipation phenomena. They solved the resulting equations analytically. They presented the convergence region for a solution through h-curves. Additionally, they presented the error analysis, along with an approximation order, for the convergent series solution. They noticed an increase in the concentration field against the Biot and Darcy numbers. The MHD nanofluid flow over an inclined sheet with radiation and heat source/sink consideration was numerically reported by Saidulu et al. [10]. They presented the comparative study and found that their results were consistent with those reported in the open literature. They monitored the decline in the heat transfer rate against thermophoresis and a radiation parameter. The involvement of dust particles in a mixed convective MHD nanofluid passed over a horizontal elongating sheet was examined by Mahdy and Hashoudy [11]. They recorded an increase in the skin friction coefficient and a decline in the heat transfer rate against the power law index. The bio-convection phenomenon for the MHD nanofluid model with mixed boundary conditions passed over an exponentially stretched sheet was analyzed by Shafiq et al. [12]. They established a theoretical study by considering the contribution of Joule heating. Hayat et al. [13] analyzed the contribution of a chemical reaction, thermal radiation and heat generation on a non-Newtonian stretched flow. They modeled the equation via the boundary layer approach and solved it by using HAM. They found a decrease in the temperature against the Prandtl number, and an escalation is recorded for the radiation parameter. The involvement of Newtonian heating in an inclined magnetohydrodynamic (MHD), non-Newtonian material was analytically studied by Shafiq et al. [14]. They found a decline in the velocity for the angle of inclination and the Weissenberg number. Swain et al. [15] discussed the role of hybrid nanoparticles in the concentration of heat energy using slip conditions passed over a porous surface. Shafiq et al. [16] studied features of entropy generation in Walters' B nanofluid in a Riga plate under the action of thermal radiation.
Mebarek-Oudina et al. [17] discussed the influences of heat convection in a trapezoidal cavity by inserting hybrid nanoparticles. Marzougui et al. [18] studied transfer phenomena in nanoparticles, considering the magnetic field inside a lid-driven cavity. Pushpa et al. [19] discussed the role of heat dissipation through convective heat energy, including nanoparticles and Bouncy force. Upreti et al. [20] formulated a 3D model related to variable heat sources and Ohmic dissipation, including aspects of the Darcy-Forchheimer flow and nanoparticles passed through a moveable surface. In another contribution, Upreti et al. [21] discussed the features of suction/injection and thermophoresis under the action of a magnetohydrodynamic. Three-dimensional heat transfer flow and entropy generation in hybrid nanoparticles over a convective surface were studied by [22]. Several important works have been reported by discussing the transport mechanism maintained in [23][24][25][26][27][28][29][30][31][32][33].
So far, extensive contributions to nanofluid research have been made, considering several effects. This report is prepared to show the contribution of ternary hybrid nanoparticles on a hyperbolic tangent model. Literature has not investigated this yet. This research will serve as a basis for researchers to work and explore further on nanoparticles. This report covers the introduction in Section 1. Modeling with thermos-physical features and dimensionless is contained in Section 2. Section 3 comprises the methodology. Sections 4 and 5 describe the solution and key findings. The scheme of hybrid nanostructures is prepared in Figure 1.

Statement of Developing Problem
Thermal features of a hyperbolic tangent martial towards a rotating frame are considered. The motion of tri-hybrid nanoparticles is induced by the movement of the wall and the rotation of the melting surface. The law related to Hall and ion-slip forces is incorporated to analyze the flow behavior in the presence of a Forchheimer porous medium. Heat characterizations occur due to Joule heating and viscous dissipation. A new theoretical scheme is implemented based on tri-hybrid nanoparticles. The composite of 2 3 , 2 and 2 in 2 6 2 is called tri-hybrid nanoparticles, while the composite relation among 2 3 and 2 is known as hybrid nanoparticles, and 2 is nanofluid. A magnetic field is considered in the normal direction of surface. A continuity equation (mentioned in Equation (1)) for steady and incompressible flow is formulated in -, -and -directions using boundary layer approximations and the law of conservation of mass. Momentum equations (mentioned in Equations (2) and (3)) are modeled, considering the Darcy-Forchheimer law and generalized Ohm's law. An energy equation is derived using the Darcy-Forchheimer law, generalized Ohm's law, Joule heating and viscous dissipation. The geometry of this model is shown in Figure 2. A set of nonlinear partial differential equations (PDEs) is modeled under boundary layer approximations (BLAs) and assumptions, while the simplified PDEs [30,31] are Here, space coordinates ( , , ), velocities ( ̅, ̅ , ̅ ), kinematic viscosity ( ), Hall number ( ), ion-slip number ( ), magnetic number ( 0 ), time constant number (Γ), electrical conductivity ( ), fluid density ( ), power law number ( ), thermal conductivity ( ), specific heat capacitance ( ), fluid dynamic viscosity ( ), permeability related to porous medium ( * ), inertia coefficient in view of porous medium ( ) and tri-hybrid nanoparticles ( ℎ ) are captured above. A set of PDEs is formulated according to required boundary conditions, and boundary conditions (BCs) are modeled as ̅ = , ̅ = 0, ̅ = 0, ̅ = , when = 0, Transformations of the current model are Using transformations, a set of non-linear ordinary differential equations (ODEs) are obtained as Dimensionless (boundary equations) BCs via variables and desired (boundary equations) BCs are defined as Correlations among tri-hybrid nanoparticles, nanofluid and hybrid nanoparticles in view of the base fluid, ethylene glycol, are defined [17] as follows and the Thermal properties are listed in Table 1. Here, volume fractions ( 1 , 3 , 2 ), wall temperature ( ), ambient temperature ( ∞ ), naofluid ( ), base fluid ( ), ethylene glycol ( 2 6 2 ), silicon dioxide ( 2 ), aluminum oxide ( 2 3 ) and titanium dioxide ( Skin friction coefficients in the presence of a hyperbolic tangent martial is modeled as The Nusselt number is a dimensionless parameter. It is noticed that the Nusselt number is the ratio of convective heat transfer and conductive heat transfer. From a mathematical point of view, the Nusselt number is defined as The rate of heat transfer, called the temperature gradient, is

Numerical Approach
It is challenging for mathematicians and scientists to derive the exact solution of differential equations (DEs). Solutions related to some of the DEs are handled easily, but some higher order (nonlinear) problems are not as easily tackled due to the complexities of DEs. Various numerical methods are used to find solutions to such problems. The Keller box method, shooting method and RK method are numerical methods with numerous advantages. However, some limitations exist, such as stability and convergence of problems. There is agreement that the finite element approach [26,29] is the best technique to solve problems. The finite element method has a great ability to handle complex geometries, along with various types of boundary conditions (BCs). A strong approach, called the finite element approach, is imposed to obtain a numerical solution of formulated ODEs associated with BCs. The steps related to the discussion of the finite element approach (FEA) are discussed below.
Step-I: Residuals are made by adjusting all terms in one place within BCs. Equations (7)-(9) are known as strong form, while the conversion from strong form into weak form is conducted via the weighted residual method. The weighted residuals of the present problem are: where 1 , 2 , 3 and 4 are weight functions.
Step-II: An approach related to the Galerkin finite element scheme is implemented to make weak forms in terms of the shape functions.
Step-III: An assembly approach is utilized for the development of the stiffness element, and the assembly approach is performed via the assembly procedure of the finite element approach (FEA). Stiffness matrices are derived as Step-IV: A Picard linearization approach provides a transformed algebraic system (linear equations).
Step-VII: 300 elements are required to obtain convergence analysis.

Graphical Outcomes and Discussion
Ternary hybrid nanoparticles are inserted into the motion of the fluid particles in the presence of ion-slip and Hall currents towards a 3D-stretching heated surface. The theory of Darcy-Forchheimer is carried out in the present analysis. Viscous dissipation and Joule heating are considered. The numerical results and graphical investigations are simulated by FEM. The details and outcomes of the present analysis are discussed below.

Outcomes of Velocity Profiles versus Parameters
The distribution of motion into fluid particles is analyzed versus the variation of the Weissenebrg number, ion-slip number, Hall number and Forchheimer number. Figures  3-12 are plotted to show the distribution of motion (for vertical and horizontal velocities) into fluid particles. Figures 3 and 4 show the influence of on the motion of fluid particles (for vertical and horizontal velocities). The motion into nanoparticles is reduced when is increased. It is deduced that the existence of is based on a hyperbolic tangent material in dimensionless momentum equations. Moreover, a decline in motion happened due to an inverse relation of versus the viscous force. Therefore, more viscosity is produced in the motion of fluid particles when is increased. Motion for the appearance of is higher than motion for the disappearance of . Vertical and horizontal velocities are boosted when is enhanced. Physically, is the ratio among elastic force and viscous force, while fluid becomes more viscous when is increased. Hence, the viscosity of fluid particles becomes thicker due to the retardation of force. Thickness associated to the boundary layers is reduced using higher values of . The influence of is visualized in the motion of fluid particles (for vertical and horizontal motions), and this impact is addressed in Figures 5 and 6. Velocity curves are reduced versus variation in . Higher values of increase the thickness of the motion of the fluid particles. It is revealed that the parameter related to is modeled due to the Forchheimer theory, and the Darcy-Forchheimer flow has a nonlinear relationship against the flow of fluid. A retardation force is visualized when the Forchheimer flow is established. Moreover, is associated with the porosity of surfaces. The flow slows down due to the porosity of the surface in the -and -directions. The role of the Hall number ( ) is shown in Figures 7 and 8 (for vertical and horizontal velocities). The motion into tri-hybrid nano-structures is boosted when the Hall number is enhanced. The Hall number exists because of the generalized Ohm's theory, and the Hall number occurs in momentum and energy equations. Moreover, a direct relationship is tabulated among the Hall number and flow distribution. The Lorentz force is upgraded when the Hall number is boosted. Therefore, a higher Lorentz force introduces the argument of the motion of fluid particles. Physically, this occurs to manufacture more speed of fluid particles. Figures 9 and 10 (for vertical and horizontal velocities) plot the relationship between fluid particles and the ion-slip number ( ). The law related to the generalized Ohm's law is used in the motion of particles, which establishes the ion-slip number ( ). The thinning layers are produced when the ion-slip number ( ) is increased due to an inverse relationship between the Lorentz force and flow distribution. Moreover, the flow distribution for = 0 is higher, as compared to ≠ 0. Both and appear in the denominator (in momentum equations). Therefore, an inverse proportional relationship is captured among velocities and ( and ). The thickness of the momentum boundary layers declined but increased in and . Figures 11 and 12 reveal the role of velocity curves in -and -directions against the variation in Ω. It is estimated that the velocity curves are reduced when Ω is increased, including the influence of tri-hybrid nanoparticles. The parameter associated with Ω is generated due to the appearance of the rotating surface. Therefore, the role of the rotating frame, brings declination into the motion of fluid particles in both horizontal and vertical directions.

Outcomes of Temperature Profile versus Parameters
Heat production is measured, considering the impacts of the Eckert number, Hall number, ion-slip number, magnetic number and Forchheimer number, as shown in Figures 13-18. Figure 16 visualizes the impact of nanoparticles, fluid, hybrid nanoparticles and tri-hybrid nanoparticles on the temperature profile. Figure 13 manufactures the production of heat energy in the motion of particles versus the variation in the Eckert number. The existence of is formulated due to the influence of viscous dissipation because occurred in a coefficient form of the viscous dissipation term. Viscous dissipation is known as work conducted, while the rate of work conducted is increased via large values . Further, work conducted is utilized in the motion of fluid particles, and the internal energy of fluid particles is boosted with the help of work conducted. Meanwhile, the temperature of particles increased. is used to determine the kinetic energy during flows in nanoparticles and hybrid nanoparticles. Further, this parameter is used to conduct enthalpy differences across boundary layers related to thermal energy, while is significantly used to measure heat dissipation in the motion of fluid particles. Figures 14 and 15 are plotted to show heat production against the impacts of the ion-slip and Hall numbers, inserting tri-hybrid nanoparticles. The occurrence of the ion-slip and Hall numbers are modeled due to the existence of a Joule heating term. Heat energy decays when the ion-slip and Hall numbers are boosted. Mathematically, the ion-slip and Hall numbers are modeled in the Joule heating term (appearing as the denominator). Therefore, an inverse relationship is captured among heat energy, and ion-slip and Hall numbers. Consequently, the production of heat energy declines when the ion-slip and Hall numbers are boosted. Thickness regarding thermal layers also decreases via higher values of ion-slip and Hall numbers. Physically, this declining trend in producing heat energy occurs due to the phenomenon of Ohmic dissipation. The thickness associated within thermal boundary layers is controlled using partially ionized fluid under the action of a magnetic field. The thermal boundary layer (TBL) decreases against the higher impact of ion-slip and Hall numbers. Figure 16 reveals the vital role of tri-hybrid nanoparticles in the heat energy of particles. The comparative consequences among fluid, hybrid nanoparticles, nanoparticles and tri-hybrid nanoparticles are analyzed by the thermal performance of heat energy. The heat production for the case of tri-hybrid nanoparticles is investigated more efficiently than the heat production for fluid, nanoparticles and hybrid nanoparticles. Finally, it is concluded that the appearance of tri-hybrid nanoparticles plays an important role for the maximum production of heat energy. In view of thermal layers, thermal layers for a case of fluid are less than thermal layers for a case of nanoparticles, hybrid nanoparticles and tri-hybrid nanoparticles. Argumentation in the thermal energy of tri-hybrid nanoparticles is measured against variations in the magnetic number, as shown in Figure 17. The appearance of 2 happens due to the effect of Joule heating in the energy equation. 2 brings argumentation into the thermal energy when Joule heating is significant. The impact of is visualized on the temperature profile, while these effects are carried out by Figure 18. It is estimated that argumentation is investigated in heat energy against the distribution in . The direct relation is formulated among and thermal energy.

Variation of Temperature Gradient and Surface Forces versus Parameters
The prediction of a temperature gradient and drag forces is visualized against variations in , and . This prediction is highlighted in Table 4. Numerical results are compared with the studies reported in [33], which are shown with the help of Table 3. Drag forces declined versus the variation in ion-slip and Hall numbers, while drag forces inclined when the impact of was enhanced in the existence of tri-hybrid nanoparticles. Moreover, the temperature gradient enhanced using variation in and . However, the production of the temperature gradient is boosted via variation in , which is monitored in Table 4.

Remarks and Key Judgments
In the present analysis, tri-hybrid nano-structures are inserted into the motion and heat energy of fluid particles. The generalized Ohm's theory is considered under the action of a constant magnetic field. Viscous dissipation and Joule heating terms are modeled. The numerical scheme, known as the finite element approach, is implemented to show the numerical and graphical consequences. The main points of this analysis are listed below:  Ternary hybrid nanoparticles play a vital role in enhancing the performance of heat energy. Moreover, heat energy is produced by hybrid nanoparticles, nanoparticles and nanofluid;  Argumentation in the motion of the fluid is predicted using the higher values of ionslip and Hall numbers, while the opposite treatment is found in the motion of fluid particles via variation in the Weissenberg and Forchheimer numbers;  The maximum production of thermal energy occurs when the Eckert number is increased, and declination in heat energy occurs via higher values of ion-slip and Hall numbers;  The maximum production of a temperature gradient can be achieved for the case of ternary hybrid nanoparticles, rather than the production of a temperature gradient for the cases of nanoparticles and hybrid nanoparticles;  Ternary hybrid nanoparticles boost the surface force, as compared to nanoparticles and hybrid nanoparticles;  Hall and ion-slip currents are significantly useful to enhance the temperature gradient, but the surface force is boosted when Hall and ion-slip numbers are enhanced.