Thermal transport analysis of squeezing hybrid nanofluid flow between two parallel plates

Hybrid nanofluids outperform mono nanofluids in terms of heat transmission. They can be found in heat exchangers, the automobile industry, transformer cooling, and electronic cooling, in addition to solar collectors and military equipment. The primary goal of this study is to scrutinize the magnetohydrodynamic hybrids nanofluid (Copper-oxides and Titanium dioxide/water) ( CuO + Ti O 2 ) and ( Ti O 2 / H 2 O ) nanofluid flow in two parallel plates using walls suction/injection under the influence of a magnetic field and thermal radiation. Implementing the appropriate transformation, the governing partial differential equations are converted into equivalent ordinary differential equations. In MATLAB software, the built-in numerical technique bvp4c is used to evaluate the final system to generate the graphical results against velocity and temperature gradient for various parameters. The increasing evaluation of Reynolds number and magnetic parameter influenced velocity concentration. The response surface method is used to generate the sensitivity and contour graphs and tables. Skin friction and the Nusselt number have an outcome on flow characteristics as well. This type of research is critical in the coolant process, aviation engineering, and industrial cleaning processes, among other fields. In several cases, comparing current and historical findings reveals good agreement.


Introduction
Nanofluids are utilized to evolve the thermal conduction and heat transport rate of the main fluid. Base fluids include tri ethylene's glycols, water, coolants, ethylene's, lubricating oils and oil, thermoplastic remedies, and biological-fluids. Nanofluids have distinct properties such as homogeneity, highly thermal conduction at low nanoparticle accordion, long-term stability, and minimum obstructing flow paths. As a result, above mention liquids are found in a vast variety of electrical devices, like as micro-electro-mechanical circuits, microchip coolants, micro reactors, and electronic fluid flow projection systems. Nanofluid technologies include structure heating and cooling, heat transfer, inspected the impact of MHD and thermal source or sink on the movement of hybrid fluid. Bilal et al. 3 are researching the use of thermal radiation to squeeze and dilate permeable walls, allowing a hybrid magnetohydrodynamic nanofluid (Carbon nanotubes and ferrous oxide/water) to flow into a vertical channel. Gul et al. 4 studied the movement of a hybrids nanofluids made up of copper & aluminum oxide molecules. Waini et al. 5 described the hybrid nanofluid, turbulent movement and heat transferal past a stretching or shrinking sheet. Huminic and Huminic 6 investigate the effects of nanofluids and hybrid nanofluids on the creation of volatility in various energy processes with valid boundary constraints and physical settings. Salman et al. 7 address recent breakthroughs in microscale-faced producing heat transfer enhancement in this study. In Nadeem et al., 8 we investigated the hybrid nanofluid flow via a permeable exponentially stretching tube. The goal of this study, according to Yashkun et al.,9 was to analyze the heat transport parameters of a magnetohydrodynamic hybrid nanofluid across a longitudinal stretching and shrinking surface in the absence of suction and radiation heat impacts. Jalili et al., [10][11][12][13] Talarposhti et al., 14 Sheikholeslami et al., 15 and Jalili et al. 16 list other studies in this field. The study focuses on the boundary layer movement of a hybrid nanofluid across a revolving cylinder in the appearance of oblique magnetism, according to Abbas et al. 17 and Aladdin et al. 18 described the properties of suction and magnetization on a revolving plate of hybrid nanofluid flow consisting of water as the main fluid, Aluminum Oxide as nanomaterials, and Copper (Cu) as nanoparticles. Alsaedi et al. 19 quantitatively estimated the flow of a hybrid nanofluid via two coaxially constructed cylinders. According to Acharya,20 fins of various forms are employed in a range of industrial and technological applications, including semiconductor devices, steam turbines, electrical converters, vehicle radiators, heat transfer, and hydrogen fuel cells, air-cooled automobiles, and so on. Heat generation in high-temperature technological processes under the effect of radiation is a significant event in heat transferable. This phenomenon is momentous in heat transmission and thermal system design. Ashwinkumar et al. 21 investigated the upshot of nonlinear heat radiation on the twodimensional magnetohydrodynamic flow of a hybrid nanofluid in two different configurations. The linear form is determined to be incorrect for the hightemperature disparity. In the modeled Rosseland estimate, a non-dimensional variable is utilized, and has no effect on the Prandtl number, according to Magyari and Pantokratoras. 22 Recently, nonlinear thermally radiative has been examined, and numerous hybrid nanofluid movements with varied geometries have been reported. [23][24][25][26][27] The magneto-hydrodynamic flow with fluctuating fluid properties has piqued the interest of the scientific community. Technology, physical science, and chemical engineering all benefit from research on magnetization consequences. Contact between an electrical conductor fluid and a magnetization has an impact on a variety of technical systems, including magnetohydrodynamics (MHD) turbines, motors, exchanges, and flow separation processes. A great quantity of research on the unique fluid characteristics across numerous geometries with varying flow conditions has been published due to several applications in geophysics, magnetohydrodynamic power production, and so on. Safdar et al. 28 investigated theoretical and computational models for stable MHD Maxwell nanofluid flow along a porosity expanded sheet including gyrotactic organisms. Hussain et al. 29 detect entropy in a hybrid nanofluid flow governed by magnetohydrodynamics, variable viscosity, and coupled convection. The radiative and MHD micro-polar fluid movement on stretching or shrinking sheets of CuO and TiO 2 nanomaterials is defined by Waini et al. 30 Thanesh Kumar et al. 31 reported on the consequence of titanium oxide nanoparticle morphologies on the magnetohydrodynamic steady movement of Maxwell fluid after passing through a wedge with nonlinear thermal radiation influences. Yahaya et al. 32 investigated heat transport on a stretched plate using a Cu À Al 2 O 3 =H 2 O hybrid nanofluid. According to Waini et al., 33 this study investigates the compressed dual nanofluids movement on a porous sensing area using magnetohydrodynamic (MHD) and radiative influence. The sensitivity and contour graphs and tables are created using the response surface method. In the presence of quadratic velocity, Zainal et al. 34 studied the MHD flow of a hybrid Al 2 O 3 À Cu=H 2 O through a permeability stretching or shrinking sheet. Refs. [35][36][37][38][39][40][41][42][43][44][45][46] offer other works on this topic.
By examining the laminar and steady flow of a hybrid nanofluid between two parallel permeable plates with a fluctuating magnetization, the continuing research closes a gap in the investigation. The subject at hand will be investigated for the initial period and has never been studied before. The most current study is a first for the field, to the best of the researcher's knowledge. The impacts of suction or injection through the channel, as well as the effects of other significant elements on the velocity distribution, temperature distribution, Nusselt number, and skin friction, were also evaluated by the researcher.

Mathematical modeling
Hybrid nanofluid flow is considered through two parallel horizontally plates diverged by a distance of H, where H denotes the space across plates at a time t = 0.
H defines as H t ð Þ = H ffiffiffiffiffiffiffiffiffiffiffiffiffi 1 À at p . For t = 1 a and a.0, the plates are crammed together till these collide to each other but both plates are disjoint as a\0, as drawn in Figure 1. Moreover, the higher plate with velocity v = dH=dt is rapidly approaching the bottom plate. The inspired magnetic field is created using a per- To make a hybrid nanofluid CuO + TiO 2 =H 2 O ð Þ , single and multiwall copper oxides are mixed with major fluid water. The xaxis is parallel to the horizontal axis, while the z-axis is orthogonal to the plate in the geometrical coordinates system. T 0 and T H are the thermal effect of the bottom and top plate, correspondingly. Because the flow parameters of the hybrid nanofluid CuO + TiO 2 =water ð Þ under study are not time-dependent, the fluctuating magnetic field has a greater consequence.
Flow is considered through two horizontal parallel plates. The effects of thermal radiation and suction/ injection are taken under consideration.

Governing equations
The following are the continuity, momentum, magnetic field, and energy conservation calculations for the abovementioned hybrid nanofluids. [47][48][49] Maxwell equations. 49 Where b 1 &b 2 are the parameters of the magnetic field, the pressure of the fluid is p, the temperature is T and s hnf is the electrical conductivity of hybrid nanofluid.
The radiative heat flux q r is defined as

Thermophysical properties
The hybrid nanofluid heat capacity is represented by rC p À Á hnf and the kinematic viscosity of hybrid nanofluid is given as m hnf . 46 Thermos-physicals properties of base fluid and nano particles are given in Table 1. and f 1 , f 2 is the volume fraction of CuO, TiO 2 respectively.

Boundary conditions
The boundary constrains are delineated as follows 48 : The transformations of similarity listed below 46 Transformed set of equation The resulting group of ordinary differential equations after the operation of factors of resemblance (9) into (2-6) As well as the evaluate boundary conditions Where the magnetic variable is signified by Reynolds's Magnetized factor, the temperature ratio parameter is

Physical components
The Nusselt numbers and skin frictions factors are described below, Numerical scheme The arrangement of ordinary differential equations is not acutely sensitive to proper clarification to accomplish an appropriate explanation. The MATLAB computational software bvp4c function and shooting method are used to code the above-mentioned combined ODEs (10)(11)(12) and the corresponding boundary constraints (13) Bvp4c is essentially a Labotto-III collocation formula that is used to generate numerical results. Higher-order derivatives of the controlling equations have deviated into a first-order set of ODEs Applying these translations into equations (10)-(12), we get As well as the modified boundary constraints Mesh size and error tolerance are both 0.001 and 10 À6 , respectively, assuring precise numerical outputs. The velocity and temperature outlines that satisfy the boundary constraints at infinite are determined based on this estimate result. We can assume a simple starting estimate for the first response for the bvp4c procedure will relate to the first solution even with bad predictions. We produced ODEs in a q-component group by inserting factor q. We employ q-variable in equations (17)- (19) to generate ODEs in a q-variable group, and therefore Differentiate by q outcomes in the following method in terms of factor-q values. Differentiate equations (21)-(23) with respect to q Where coefficient of the matrix is denoted by h 1 , e 1 is the reminder, and d 1 = dP i dt , 1 ł i ł 8: Cauchy Problem is Vector functions are denoted by y 1 , v 1 . We can automatically fulfill the ODEs by determining the two Cauchy problems for each characteristic.
With the help of numerical solution: For the solution to the problem, a complete scheme was adopted.
Taking the appropriate factors, Demarcated restrictions are commonly implemented for P i , where 1 ł i ł 8, for the solution of ordinary differential equations, we are instructed to apply d 2 = 0, this appears to be a matrix as Where a 1 = Àl 1 :y 1 l 1 :v 1 :

Response surface method (RSM)
The RSM is used to establish empirical relationships between a variety of input parameters and a variety of outputs. Given that it provides information on the least and most dominant input factors influencing the answers, it is a particularly excellent tool for evaluating multi-response and multi-variable processes. RSM was founded in 1951 by Box and Wilson, 54 and its main objective is response optimization. A collection of mathematical and statistical tools known as the Response Surface Methodology (RSM) can be used to model and analyze problems when a variety of variables affect the interest response. The goal is to make this answer better. The link between the responses and the independent factor in RSM problems is frequently ambiguous. Finding a good estimation for the best feasible connection between y and the group of independent variables is thus the first stage in RSM. In several regions of the independent factors, a moderate polynomial is widely used. If the answer can be precisely characterized by a linear function of the independent variables, then the first-order model is used as the approximation function. The response surface methodology was used to finally establish an empirical relationship between the convective heat transfer coefficient and optimization. To plot contour, sensitivity outcomes and tables, the response surface method is employed. 55 Coding for factors is used.

Sum of square is types III -Partial
The methodology is proposed to be meaningful by the simulation F-value of 3543.59. These values can be seen in Table 2. An F-value this strong might happen owing to interference only 0.01% of the time.
Simulation aspects are considered meaningful when the p-value is lower than 0.0500. A, B, C, AC, BC, A2, B, and C2 are essential structural components in this instance. Regression coefficients are not meaningful if the value is higher than 0.1000. Simulation minimization may enhance your concept if it has a lot of unnecessary phrases (except those needed to maintain hierarchy). Table 3 is Boxes-Behnken's Concept against the real and the coded values, represent the combination of the parameter of volume fraction of nano particles f, thermal radiation parameter Rd and the Biot number Bi response of local skin friction and local Nusselt number.
For the goodness of fit persistence, relay on several indicators is considered. For instance, the first one is the lack of ft which has a p-value \0:001 for all three fitted models. Secondly adjusted R 2 which presents how much the models explain the variation in response is used. For Cf the adjusted R 2 is 96.2%, and for Nu adjusted R 2 is 99.9%. As a result, all models account for a large amount of the variation in the corresponding responses. Thirdly, a typical residual quantilequantile plot is used to evaluate the goodness of ft. In an appropriate model that successfully demonstrates the functional relationship between the input parameters and response, the relationship between theoretical and observed quantiles is proved to be one-to-one.
For three responses, three response surface equations are considered. For each model parameters A, B, and C are Nanoparticle volume fraction f ð Þ, thermal radiation Rd, and Biot number Bi. For each of these parameters, three levels are low, medium and high levels coded as (21, 0, 1) are chosen. The input parameters together with respective notations are presented in Table 3.

Results and discussions
The influence of precise emerging variables in the governing flow framework on the velocity concentration, temperature panel and magnetic distribution is evaluated through graphically and tabular form. The influence of different included components such as Prandtl number, magnetic Reynolds number Re m , magnetic parameter, Eckert number squeezing component and nanomaterials volume fractions f have analyzed the behavior of heat and mass transfer. The bvp4c methodology, which is a significant numerical method, is used to decode the difficulties at hand. The response surface method is used to plot the sensitivity and contour graphs and tables. The ranges of pertinent parameters are given as 0:1łSł1:2,0:1łM ł1:5,1:0łRe m ł4:0,0:1 łEcł1:0, 0:1łRdł1:5, and 1:5łu w ł1:8: The consequence of squeezing components S via velocity panel f is demonstrated in Figure 2. As the movement on the surface, the fluids are squeezed into the channel. From this figure, the velocity of fluids at the boundary region is increased due to increasing outcomes of the squeeze variable. Figure    formed as a result of the magnetic field. The Lorentz force had an effect on the channel's boundary area. This downward tendency is noticed because the introduction of a transversal magnetization causes the energy of the resistance nature known as Lorentz force, which tends to hinder fluid flow and so decelerates velocity. Figure 6 has drawn the implication of the Reynolds number variable Re m on the axial flow panel f 0 . By enhancing outcomes of the Reynolds number factor Re m the velocity f 0 of fluids is increased. The influence of the squeezing component S on the axial velocity f 0 is shown in Figure 7. By increasing the outcomes of the squeezing variable, S the velocity of fluids is raised. The kinematic viscosity of the fluid divided by the top plate's flow rate is known as the squeezing factor. It's worth noting that large or small S values denote the sluggish or quick flow. Greater S outcomes suggest that the lower plate is moving away from the top plate, or that the distance between the two plates is growing, whereas diminishing S values show that the top plates is revolving outside of the bottom plate, or that the gap among the two plates is decreasing.    The impact of the squeezing factor S on the thermal gradient u is drawn in Figure 8. By enhancing the consequences of the squeezing component, the thermal implication of fluids is downward. Practically, the increasing outcome of S the fluid flow causes friction, which induces heat, as well as enhances the fluid's kinematic energy, which also generates heat. Figure 9 plotted the impact of nanomaterials volume fraction f on the temperature panel u. By rising the outcomes of nanomaterials volume fraction f the fluid's temperature is upward. The temperature gradient u of nanofluids becomes parabolic as the nanomaterials volume fraction f increases, approaching the bottom in the center of the channel. The consequence of the Eckert number Ec on the temperature panel u is sketched in Figure 10. The thermal impact u of nano-fluids is enhanced by the increment in consequences of the Eckert number Ec. Flow friction in nanoparticles is created with heightened efficiency as the Ec increases. Kinetic energy is transferred into thermic energy in this physical process, which aids in the increment in a temperature gradient, and for greater results of Eckert number, the temperature panel turn into parabolic and reaches its extreme value in the center of the channel. Furthermore, boosting the volume proportion of nanomaterials causes the density of fluid particles to increase. Fluid flow is valued in this physical process,     Figure 11. By enhancing outcomes of radiation parameters Rd the temperature panel u of fluids is upsurge. Temperature and associated boundary layer thickness decrease as heat radiation increases. Modifications of the value of thermal radiation are due to the existence of conduction over radiation. Figure 12 exemplified the behavior of u w on the thermal distribution u. The intensifying values of the temperature ratio factor have a positive effect on the temperature profile as observed from the graph. Figure 13 illustrated the features of the squeezing component S on the K. The magnetic profile K of liquids is decline by increasing impact of the squeezing variable S. It shows that the K is mostly linear at S = 0:1, but becomes parabolic as S rises, approaching the maximal outcomes in the channel's middle. Figure 14 described the residual plots for squeezing flow. This figure describes the percentage of residual on normal probability, residual via fitted value on versus fits, frequency of residual on the histogram, and residual via observation order on versus order. Figure 15 Discussed the contour plot of response via thermal radiation parameter and volume fraction of nano particles. Figure 16 defines the Contour plot of response volume fraction and Biot number. Figure 17 depicted the contour plot of response via radiation parameter and Biot number. Figure 18 mentioned the 3-D surface plot between thermal radiation parameter and volume fraction parameter. Figure 19 represented the 3-D surface plot between radiation and volume fraction. Figure 20 highlighted the 3-D surface plot between Biot number and volume fraction.

Conclusion
A hybrid liquid flow between two plates is investigated in the existence of a magnetized factor. Utilizing resemblance modification, the suggested flow model's controlling equations are turned into highly nonlinear structures of ordinary differential equations. For the flow model's formulation, the simulation model bvp4c was operated. On the flow panel, thermal profile, and magnetic panel, the effects of several related factors are presented using various graphs and tables. The response   surface method is employed to plot the sensitivity and contour graphs and tables. The key values of the current research are given following: v Due to the improving value of nanoparticle volume fraction, the boundary layer's consistency is diminished, and the velocity of fluids is raised. v Due to rises outcomes of the Eckert number the temperature of fluids is boosted.
v The magnetic distribution of fluids is reduced by increment in the squeezing factor. v Residual plots for squeezing flow, third surface plot via different parameters, and contour plots via various factors are also discussed. v The insertion of nanoparticles into the base fluid significantly develops the velocity of the hybrid nanofluid and the rate of energy transfer.   Applications for hybrid nanofluid flow include petroleum engineering, geothermal engineering, automobile manufacturing, nuclear waste storage, thermal extrusion systems, heat exchangers, energy resources, ventilation, grain storage, and other industries. Because of their high thermal conductivity and low viscosity, hybrid nanofluids may have an impact on heat transfer systems. 54,55 Declaration of conflicting interests The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding
The author(s) received no financial support for the research, authorship, and/or publication of this article.