Numerical Analysis of an Unsteady, Electroviscous, Ternary Hybrid Nanofluid Flow with Chemical Reaction and Activation Energy across Parallel Plates

Despite the recycling challenges in ionic fluids, they have a significant advantage over traditional solvents. Ionic liquids make it easier to separate the end product and recycle old catalysts, particularly when the reaction media is a two-phase system. In the current analysis, the properties of transient, electroviscous, ternary hybrid nanofluid flow through squeezing parallel infinite plates is reported. The ternary hybrid nanofluid is synthesized by dissolving the titanium dioxide (TiO2), aluminum oxide (Al2O3), and silicon dioxide (SiO2) nanoparticles in the carrier fluid glycol/water. The purpose of the current study is to maximize the energy and mass transfer rate for industrial and engineering applications. The phenomena of fluid flow is studied, with the additional effects of the magnetic field, heat absorption/generation, chemical reaction, and activation energy. The ternary hybrid nanofluid flow is modeled in the form of a system of partial differential equations, which are subsequently simplified to a set of ordinary differential equations through resemblance substitution. The obtained nonlinear set of dimensionless ordinary differential equations is further solved, via the parametric continuation method. For validity purposes, the outcomes are statistically compared to an existing study. The results are physically illustrated through figures and tables. It is noticed that the mass transfer rate accelerates with the rising values of Lewis number, activation energy, and chemical reaction. The velocity and energy transfer rate boost the addition of ternary NPs to the base fluid.


Introduction
The squeezing flow is important in lubrication equipment, polymer processing, molding compaction, and injection, all of which use hydrodynamical technologies generated from moving surfaces. Jackson [1] highlights the connection between loaded bearings and compressing flow operation in engines, which includes the occurrence of adhesion. Muhammad et al. [2] describe the unsteady squeezed flow of a hybrid nanofluid (HNF), made up of CNTs and CuO, using a numerical technique. The fluid velocity improves as the squeezing intensity and volume fraction of nanomaterials increases. Ramesh et al. [3] explore the unsteady squeezing flow of hydromagnetic and Casson NFs using enclosed parallel disks. Selimefendigil et al. [4] conduct a numerical analysis of nanofluid forced circulation inside a branching channel, under the effect of a constant magnetic flux. Xu et al. [5] study stable power law NF flow, including gyrotactic microorganisms that transmit energy between two parallel plates. Shuaib et al. [6] use 3D flow characteristics to display a 3D NF flow across two contemporaneous circular plates. Their purpose is to see how magnetic induction affects NF flow with heat transmission qualities in the long run. To examine the MHD NF flow through extendable spinning discs, Ahmadian et al. [7] employ the numerical approach bvp4c. The disc stretching process, according to the findings, opposes the flow tendency. Bilal et al. [8] consider the effects of MHD and EHD parameters on the flow of water-based hybrid NFs across two circular plates. With the growing Reynolds number, magnetic, and electric effects, heat transmission is estimated to rise [9,10]. Khan et al. [11] use Fourier's and Fick's laws to explore the channel flow of a second-grade viscoelastic fluid between two plates, formed by a vibrating wall with mass and energy transport characteristics. Alsallami et al. [12] develop an Maxwell nanofluid flow with arrhenius activation energy over a rotating disk. Dawar et al. [13] deploy freshwater as a conventional fluid across two surfaces in order to study copper oxide and copper nanoparticles. The results demonstrate that the particulate concentration factor has a dual influence on velocity distribution. Some recent studies are found in [14][15][16][17].
A hybrid nanofluid is a new type of fluid that outperforms when compared to regular fluids, such as ethyl alcohol, water, nanofluids, and ethylene, during energy transitions. HNFs have a huge spectrum of thermal properties, including the ability to freeze at high temperatures [18][19][20]. Hybrid NFs are used in energy generation, heat transfers, heat pumps, air conditioners, the automotive industry, electrical appliances, turbines, nuclear reactors, broadcasting, spacecraft, and biotechnology [21]. Coolant and freezing solutions containing ethylene glycol and water provide significant advantages, such as reducing corrosion and acid deterioration, as well as suppressing the growth of most bacteria and fungi. In the industry, ethylene glycol and water mixes are referred to as glycol concentrates, mixtures, solutions, and compounds. We employed the TiO 2 , SiO 2 , and Al 2 O 3 NPs in the working fluid in this study. TiO 2 is a white inorganic material that has been used for over a century in a variety of foodstuffs [22]. It is the brightest and whitest pigment known, with reflective qualities, and the ability to absorb and scatter UV radiation [23,24]. Silica is one of the most abundant class of compounds, and because of its inexpensive cost of manufacture, high surface area, and wettability, it has the potential to be an outstanding choice for commercial usage [25,26]. Zhang et al. [27] inspect the entropy maximization in a hydromagnetic HNF flow of SiO2 and (MoS2 (molybdenum disulfide) NPs flowing toward a stretchy surface. Ahmed et al. [28] inspect the Ag-MgO HNF flow with heat propagation generated by a curved spinning disc that rotates in three dimensions, both vertically and horizontally. Chu et al. [29] evaluate flow kinematics and heat transfer from the perspective of TiO 2 and Al 2 O 3 NPs used to rise the thermal characteristics of the base fluid. The efficiency of an HNF containing TiO 2 and MWCNTs is analyzed by Chu et al. [30]. Long [33] provide a concise overview of most of TiO 2 s notable properties, as well as a summary of innovative thermal applications, including its high refractive index, overwhelmingly high boiling and melting points, high stiffness and hardness, ability to absorb or reflect ultraviolet radiation, and photocatalytic nature. Recently, many researchers reported significant contributions to ternary hybrid nanofluid flow [34][35][36][37].
Magnetism is one of the key significant features of engineering and drug distribution due to its wide range of applications; heat exchangers, clutches, and compressors, to name a few major commodities, are all modified by the combination of flowing fluid under a magnetic field [38][39][40]. Magnetic fields have the potential to regulate and make the working temperature of a number of industrial devices more convenient. Magnetic fields are used in interplanetary and extremely high magneto applications, as well as in aerodynamics and chemical chemistry. The strength and scattering of the applied magnetics have an impact on the flow behavior. To describe the flow properties under the upshot of magnetic flux, many researchers committed to fluid mechanics. Hayat et al. [41] look at how specific heat and a produced magnetosphere affected the sinusoidal flow of an HNF flow, via a lateral tube. Raza et al. [42] examine the effect of a molybdenum disulfide nanofluid exchanger and MHD on free convective flow through a channel. Dezfulizadeh et al. [43] explore the performance of MHD ternary HNF flow in a thermal exchanger, using a unique compoundperverted turbulator and spiral rotors. Per the PEC indices, the twisted spinning bar with ellipsoidal surface obstacles has the maximum exergy effectiveness, increasing by 7% in Re, to 12,000. References [44][45][46] contains some recent literature on MHD HNF.
The present study aimed to numerically examine the cumulative influence of the electromagnetic force, chemical reaction, suction/injection, inertia force, activation energy, ionized fluid, and magnetic field on the squeezing flow of ternary hybrid nanofluids across parallel plates. We supposed that the lower plate is permeable and stretching with a uniform velocity. For this purpose, the phenomena were modeled and formulated in the form of a system of PDEs, which are solved through the parametric continuation method. The results are shown through figures and tables. In the above-described, ionized, ternary nanofluid model, the effects of activation energy, heat source, and chemical reactions in the uses of ternary hybrid nanofluid are the main novelty of the proposed model. Furthermore, in the next section, the problem was articulated, resolved, and discoursed.

Governing Equations
The ternary hybrid nanofluid flow across two parallel infinite plates, consisting of titanium dioxide, silicon dioxide, and aluminum oxide is reported. The flow mechanism is graphically depicted in Figure 1. The upper plate is located at . The lower plate is permeable, which allows suction/injection effect, signified as V w = −V 0 (1−αt) . Both plates are assumed at constant temperatures T 1 and T 2 . The lower plate is expanding with the linear velocity u w = −bx (1−αt) . Furthermore, the time-dependent magnetic field is characterized as (1−αt) . The basic flow equations are communicated as [47][48][49]: ∂u ∂x Micromachines 2022, 13, 874 5 of 17 The similarity variables are: Therefore, the transformed set of ODEs is: , 2

SqPr
Pr Equation (1) is the continuity equation, Equation (2) is the momentum equation with electroviscous and uniform magnetic effect, Equation (3) is the Poisson equation, Equations (4) and (5) are Nernst-Planck equations, while Equations (6) and (7) are the energy and mass distribution equations, respectively.
In Equations (1)- (7), n − and n + are the negative and positive charged ions; U = ∂v ∂x − ∂u ∂y are the associate's condition for the upper and lower plate; T and C are the temperature and the concentration, respectively; Fr = C * b √ K * is the porous media non-inertial coefficient, where K* and C * b are the permeability factor and drag force constant, respectively; Q 0 is the heat generation term; K 2 = 2z 2 e 2 n 0 ε 0 εk b T is the inverse Debye factor; W is the electric potential of ions; and B = ρk 2 T 2 ε 0 ε 2z 2 e 2 µ 2 is fixed at a constant temperature. E a is the activation energy, and k T is the chemical reaction rate. Furthermore, ρC p Thn f , µ Thn f , σ Thn f , κ Thn f and ρ Thn f are the heat capacity, dynamic viscosity, electrical conductivity, thermal conductivity, and density of ternary HNF, respectively.
The boundary conditions are: The similarity variables are: Therefore, the transformed set of ODEs is: Here, . M is the magnetic term, Pr is the Prandtl number, Sq is the squeezing constraint, K * 1 is the local porosity term, λ > 0 is the stretching parameter for the lower plate (λ = 0) for the fixed plate), Fr is the Forchheimer number, Sc is the Schmidt number, Nt is the thermophoresis constant, S is the suction/injection term, Le is the Lewis number, Q is the heat source/sink term, Nb is the Brownian motion constant, E is the Arrhenius activation energy coefficient, and σ is the chemical reaction term, defined as: The transform boundary conditions are: The Nusselt number and the skin friction are characterized as: where

Numerical Solution
This section shows how to use the algorithm of the numerical scheme to solve the numerical solutions of the suggested mathematical model. The main steps for dealing with the parametric continuation method scheme and future direction are as follows [50][51][52][53][54][55][56][57]: Step 1: Simplifying the modeled equations to 1st order: By putting Equation (19) in Equations (10)- (15) and (17), we achieve: with the corresponding boundary conditions: Step 2: Introducing the embedding parameter p: Step 3: Differentiating by parameter 'p': where ∆ is the coefficient matrix.
Step 4: Apply the Cauchy principal: where W and U are the indefinite vector functions.
By putting the approximate solution Equation (26) into the original Equation (24), we obtain: Step 5: Solving the Cauchy problems: Finally, we obtain:

Results and Discussion
This section reveals the physical description of the obtained results in form of figures and tables for velocity, energy, and mass transfer profiles versus several physical constraints. Velocity Profile (f (η)): Figure 2a-e displays the tendency of velocity profile(f (η)) versus suction parameter S > 0, Darcy-Forchheimer Fr, magnetic field M, injection S < 0, and volume friction of ternary nanoparticles Ψ = (φ 1 = φ 2 = φ 3 ), respectively. Figure 2a-c shows that the velocity field is lessened with the upshot of suction parameter, Darcy-Forchheimer, and magnetic field. Physically, the rising effect of the suction factor diminishes the motion of fluid particles, which causes a reduction in the velocity profile, as shown in Figure 2a. The permeability of the plate surface enhances with the variation of the Darcy effect, which also encourages more suction from the plate surface and, as a result, fluid velocity (f (η)) declines, as depicted in Figure 2b. The resistive force, which is created due to the magnetic effect M, opposes the fluid motion, similarly also deducing the velocity boundary layer, as shown in Figure 2c. Figure 2d,e illustrate that the velocity distribution accelerates with the flourishing values of injection and volume friction of ternary nanoparticles. Physically, due to the injection effect of fluid particles, the fluid moves fast; as a consequence, the velocity of fluid flow elevates, as elaborated in Figure 2d. The addition of ternary nanoparticles (TiO 2 , SiO 2 , Al 2 O 3 ) to the base fluid magnifies its thermal conduction, which also causes the inclination in the velocity field, as revealed in Figure 2e.
Electric Field (g(η), h(η)): Figure 3a-d report the presentation of the electric field (g(η), h(η)) profile versus the Schmidt number Sc, and squeezing term Sq, respectively. The kinetic viscosity of fluid enhances with the effect of the Schmidt number, which diminishes the molecular dissemination, and causes the lessening of the electric field (g(η), h(η)), as shown in Figure 3a,b. The influence of the squeezing variable fluctuates the fluid particles, which enhances its velocity, and as a result, the electric profile is also boosted, as elaborated in Figure 3c,d.
Energy Profile θ(η): Figure 4a-c represents the tendency of energy profile θ(η) versus the heat source Q, volume friction of ternary nanoparticles Ψ, and injection S < 0 constraints, respectively. As a consequence of the heat generation term, thermal energy is generated inside the fluid flow, which causes the elevation of the energy profile θ(η), as seen in Figure 4a. Figure 4b illustrates that the inclusion of nanomaterials (TiO 2 , SiO 2 , Al 2 O 3 ) in the base fluid augments the thermal conduction of the base fluid, as well as reduces the average heat capacity, because the specific heat capacity of ethylene glycol/water is much higher than ternary nanoparticles. This is why the energy propagation rate of ternary nanofluid magnifies with the rising quantity in the concentration of nanoparticles Ψ. The energy transfer rate of ternary NFs declines with the upshot of the injection term, as shown in Figure 4c. volume friction of ternary nanoparticles ψ , and injection S < 0 constraints, respectively.
As a consequence of the heat generation term, thermal energy is generated inside the fluid flow, which causes the elevation of the energy profile ( ) θ η , as seen in Figure 4a. Figure   4b illustrates that the inclusion of nanomaterials (TiO2, SiO2, Al2O3) in the base fluid augments the thermal conduction of the base fluid, as well as reduces the average heat capacity, because the specific heat capacity of ethylene glycol/water is much higher than ternary nanoparticles. This is why the energy propagation rate of ternary nanofluid magnifies with the rising quantity in the concentration of nanoparticles ψ . The energy transfer rate of ternary NFs declines with the upshot of the injection term, as shown in Figure 4c.   Mass Profile φ(η): Figure 5a-e displays the trend of mass profile φ(η) versus Le, Nt, Nb, activation energy E, and chemical reaction rate σ, respectively. Figure 5a-c elaborates that the transfer rate boosts with the rising values of Lewis number, which remarkably declines with the upshot of Nt and Nb. Physically, the molecular diffusion rate reduces with the variation of Le, which results in the reduction in the concentration boundary layer, as seen in Figure 5a. Furthermore, we are interested in investigating the influence of Nb and Nt on the flow mechanism, as these are two important factors that govern nanofluid movement. Brownian motion is a haphazard motion occurring as a result of nanomaterials in a fluid flow. Brownian motion is more powerful in fluids with low viscosity and elevated heat, as well as in fluids with tiny particles. However, their effect reduces the mass proportion ratio, as manifested in Figure 5b,c. The increment in activation energy constraint E and chemical reaction term dramatically elevate the mass transmission ratio, as publicized in Figure 5d,e. The effect of both factors accelerates the kinetic energy inside the fluid, which encourages fluid particles to move fast; as a result, the concentration profile of ternary nanofluid enhances.  (TiO2, SiO2, Al2O3), and the physical model of ternary hybrid nanofluid, respectively. Table 3 reveals the statistical assessment of current results with the available work for validity purposes. It is perceived that both results show Tables 1 and 2 expose the experimental values of the base fluid (ethylene glycol/water) and ternary nanoparticles (TiO 2 , SiO 2 , Al 2 O 3 ), and the physical model of ternary hybrid nanofluid, respectively. Table 3 reveals the statistical assessment of current results with the available work for validity purposes. It is perceived that both results show the best settlement. Furthermore, the influence of magnetic term and suction constraints enhances the skin friction of both the upper and lower plate. Table 4 exhibits the comparative valuation of present outcomes with the published literature for the Nusselt number −(θ (1)).   [23].
Viscosity µ Thn f

Conclusions
We studied the properties of transient, electroviscous, ternary hybrid nanofluid flow through squeezing parallel infinite plates. The ternary HNF was manufactured by adding the TiO 2 , SiO 2, and Al 2 O 3 to the conventional fluid glycol/water. The ternary hybrid nanofluid flow was modeled in the form of the system of partial differential equations, which were subsequently simplified to a set of ODEs through resemblance substitution. The obtained nonlinear set of dimensionless ODEs is further solved via the parametric continuation method. The key findings are: • The velocity field f (η) is reduced with the effect of the suction parameter, Darcy-Forchheimer, and magnetic field.

•
The flourishing values of injection and volume friction of ternary nanoparticles (TiO 2 , SiO 2 , Al 2 O 3 ) accelerate the velocity distribution.

•
The electric field (g(η), h(η)) declines with the upshot of Schmidt number Sc, while enhancing with the increment of squeezing term Sq.

•
The thermal energy field θ(η) is elevated versus the variation of heat source and the inclusion of nanomaterials to the base fluid, while reducing with injection effect.

•
The mass allocation rate boosts with the rising values of Lewis number, activation energy constraint E, and chemical reaction, while declines with the upshot of thermophoresis and Brownian motion.