Effect of Couple Stress and Mass Transpiration on Ternary Hybrid Nanoliquid over a Stretching/Shrinking Sheet with Heat Transfer

The present article describes the unsteady flow of a couple stress via a ternary hybrid nanofluid on a stretching surface with porous media. The nanofluid exhibits important properties for increasing heat transmission, and it is widely used in manufacturing and industrial applications. The basic similarity equations have been discovered to accommodate both stretching/shrinking surfaces. Ternary hybrid nanofluid is a colloidal combination of three types of microspheres: Al2O3, single wall carbon nanotubes, and graphene. For investigating spherical, cylindrical, and platelet nanoparticles, the governing partial differential equations are converted into ordinary differential equations, expending appropriate transformations. The analytical solution can then be carried out using various forms of nanoparticles, such as spherical, cylindrical, and platelet, to obtain the solution domain. Heat transfer is used in an electrically conducting fluid and also including thermal radiation, as calculated with the Biot number. The focus of the present effort is the evaluation of the flow of ternary hybrid nanofluid over a porous media via thermal radiation, with couple stress, using an analytical process. For various physical parameters, the velocity and temperature conditions are shown graphically.


Introduction
The authors investigated the stretching sheet problem by examining its importance and applicability in industrial processes such as polymer extrusion, paper production, metal cooling, and glass blowing, among others [1]. The pioneering work on the surface stretching problem is by Sakiadis [2,3]. Crane [4] takes up the problem of the stretched surface, where the Newtonian flow is linearly changing from the slit. Later, several studies were conducted on the stretching sheet problem with several boundary conditions. In view of stream and water surface is investigated by Gu et al. [5]. The thermal conductivity and dynamic viscosity have calculated as these possibilities, the impact of mass transfer and energy properties concerning the air by Esfe et al. [6] using a hybrid nanofluid. Bhattacharyya et al. [7] studied the properties of hydrodynamic heat transport. Chamkha et al. [8] found that heat transfer across the dynamics of a hybrid nanofluid is higher than that of an alumina-water nanofluid, but lower than that of a copper-water nanofluid. Mahdavi et al. [9] investigated a fluid flow and heat transfer analysis of nanofluid jet cooling on a hot surface with various roughnesses. MHD nanofluid flow with suction, as well as magnetohydrodynamics a nanofluid flow in natural convection with a porous term, was studied by Benos et al. [10]. Besides these mixture, hybrid nanofluids are defined as fluids containing two types of nanoparticles. A hybrid nanofluid is a specific type of fluid that has superior thermal conductivity to either nanofluids or base fluids. Anusha et al. [11,12] investigated a hybrid nanofluid that included both transpiration and radiation. Sharifpur et al. [13] conducted an experimental investigation and developed a model for evaluating the thermal conductivity of α-Al 2 O 3 -glycerol nanofluids. In industrial settings, such as in the acceleration of sheets, where liquid cooling is crucial to achieving the desired result, the temperature distribution impacted by thermal radiation is highly significant. This type of flow problem can be seen in engineering applications such as extrusion operations, metal spinning, dye casting of metals, and polymer manufacture, where the maximum temperature and difference prediction should be controlled. This study has attracted the interest of many researchers due to the numerous exciting industrial and everyday applications of shear-thinning fluids. Wall paint, printing ink, nail polish, whipped cream, ketchup, and engine oil are a few examples of these. Shear thinning fluid, also known as pseudo-plastic fluid, is thought to exhibit characteristics of both plastic and Newtonian fluid. The more tension that is applied to a shear-thinning fluid, the more freely the fluid flows. This quality makes it a beneficial trait for products such as paint, oils, and cream [14][15][16].
According to a study by Zhang et al. [17] on thermal diffusivity and thermal conductivity, if a colloidal mixture of various nanofluids is generated by water conveying cylindrical and spherical nanoparticles, the normalized thermal conductivity of the cylindrical nanoparticles increases. Furthermore, water transporting CNT nanoparticles has a far higher thermal conductivity than nanofluids carrying spherical nanoparticles. Wakif et al. [18] explored the internal heating of spatially homogeneous and non-uniform nanofluids. Narayana et al. [19] studied an unstable MHD nanofluid subjected to non-uniform heat generation and absorption. Many equations have been proposed for to nanofluid fluid movement across porous media. The formal structure of issues regarding a porous medium was suggested by Darcy and Brinkmann [20,21]. Anuar et al. [22] studied the heat transfer of carbon nanotubes over an exponentially stretching/shrinking sheet using suction and the slip effect. Anwar et al. [23] investigated the unsteady MHD flow of Jeffery fluid flow using wall velocity and Newtonian fluid. The fluid flow and heat transfer of carbon nanotubes, along with a flat plate with a Navier slip boundary, were investigated by Khan et al. [24]. Shalini et al. [25] studied the unsteady MHD chemically reacting mixed convection nanofluid flow past an inclined pour stretching sheet, using the slip effect, as well as a variable thermal radiation and heat source The MHD flow of a Jeffery fluid with mixed convection on a porous medium was investigated by Yana et al. [26] using a radiative heat flux stretching/shrinking surface. Stokes is recognized for developing the concepts and equations for pair stress fluid flows, which is another important point to make. The governing equations are provided for an incompressible pair stress fluid flow with conservative body forces. Some insights regarding these ideas can be found in refs. [27,28].
The current work is an examination of ternary hybrid nanofluid flow under the major effects of a couple stress parameters, radiation, and transpiration, as seen in the above-mentioned studies. Constitutive law is one of the main ideas at work here; it is the relationship between the forces imposed on that substance and the resulting deformation at a microscopic level. It also has to do with how much strain a body is experiencing in relation to how much stress the material is under. A porous medium is used to stretch the fluid. Researchers are interested in conducting a stretching sheet experiment with a porous medium because fluid flow with a porous medium has several industrial applications. In addition, the current situation is an example of a couple stress with a porous media challenge. By using the appropriate similarity variables, the governing equations of the supplied PDEs are mapped into ODEs. The solution domain is then obtained by solving the ODEs analytically. The solution is achieved with the use of the incomplete gamma function regarding the heat transfer, which is induced by a stretching sheet with a linear velocity variation from the slit, with mass suction included. Throughout the equation, solid volume fractions of the three distinct nanoparticle structures are employed. In fact, there are few studies in the literature that cover so many aspects of this technique. These insights could be used in the development of industrial procedures for the creation of the best product. The novelty of the present work explains the flow of couple stress fluids, especially unsteady flow, with three types of different nanoparticles; thus, we extract the results on the basis of the analytical method. This work motivated researchers to conduct more studies on stretching sheet problems with the porous medium in the presence of ternary nanofluids, since there are many equations and derivations available to explain the porous medium and thermal radiation.

Physical Model and Formulations
The current study considers a couple stress fluid flow over a porous stretching sheet under the influence of radiation and mass transpiration. In the present work, we use non-Newtonian coupling. Moreover, it should be noted that the water containing spherical aluminum oxide, cylindrical single wall CNT, and platelet graphene nanoparticles, as shown in Figure 1, is transported. Table 1 shows the thermal characteristics of these three nanofluids. The following formula is used to calculate the unstable flow stream velocity far enough from the sheet U ∞ = − U 0 x 1−γt , suppose that U w and V w are the wall movement speeds, so the standard absorption speed is taken into account. The exterior free stream velocity at the wall is calculated as follows U w = λU 0 x 1−γt , where, λ represents the constant, and λ > 0 the sustaining flow, if λ < 0 be in conflict with the flow λ = 0.
Throughout the equation, solid volume fractions of the three distinct nanopartic structures are employed. In fact, there are few studies in the literature that cover so man aspects of this technique. These insights could be used in the development of industri procedures for the creation of the best product. The novelty of the present work explain the flow of couple stress fluids, especially unsteady flow, with three types of differe nanoparticles; thus, we extract the results on the basis of the analytical method. This wor motivated researchers to conduct more studies on stretching sheet problems with th porous medium in the presence of ternary nanofluids, since there are many equations an derivations available to explain the porous medium and thermal radiation.

Physical Model and Formulations
The current study considers a couple stress fluid flow over a porous stretching she under the influence of radiation and mass transpiration. In the present work, we use no Newtonian coupling. Moreover, it should be noted that the water containing spheric aluminum oxide, cylindrical single wall CNT, and platelet graphene nanoparticles, shown in Figure 1, is transported. Table 1 shows the thermal characteristics of these thre nanofluids. The following formula is used to calculate the unstable flow stream veloci far enough from the sheet Note that Equation (2) includes extra terms that provide a model beyond the standard Newtonian fluid, which are common for non-Newtonian fluids. In particular, the biharmonic term proportional to η 0 represents the connection between materials considering the couple stress. In addition, the imposed boundary conditions are as follows: The dimensionless variables are easier to locate via the similarity transformation.
The adopted similarity functions are as below: As a result, the velocity of the wall transpiration increases v w (x, 0, t) = U 0 ν f 1−γt f (0). The boundary layer hypothesis is used in the momentum and energy equations, which are two-dimensional equations. The thermal boundary conditions for two altered situations will be explained advanced. By conjecturing two similarity equations.
The quantities of ternary hybrid nanofluid are given below: Suganthi and Rajan et al. [31] discovered that the motion of solid particles in solid-liquid dispersion is dependent on particle form. Non-spherical particles have more motion and fluid flow resistance than spherical particles. The number of spherical nanoparticles is determined by: the quantities of cylindrical nanoparticles are given by: the quantities of platelet nanoparticles are given by To make further analysis easier, the following substitutions can be made, as follows:

Solution of Pressure Gradient
Its assets show that the pressure gradient ∂P ∂y is a consequence of t and y, as determined by the velocity elements, and Equation (4) is sovereign by x, namely, ∂P ∂y = F(t, y) and ∂x , which would be independent of y and remains equal inside the boundary layer and could be deduced in the free stream. Compute the x-momentum to the free stream with u = U ∞ and ∂p/∂x The velocity acceleration/deceleration feature is shown by β, where β = γ U 0 is the unsteadiness parameter. When β > 0, wall velocities accelerate, β < 0, the wall decelerates, β = 0 steady state, and p 0 is integration constant.

Solution of Momentum and Energy Equation
In order to determine the momentum equation, an additional similarity variable F(η) = f (η) + η is brought into the system, and the resultant governing equation is: In addition, the specified governing boundary conditions are transformed into: In the boundary conditions, we get solution p + q = V C , q = − (λ−1) α and: Apply the assumed solution to Equation (20) into Equation (18) to obtain the exponent α value as: Then we obtain the solution:

Solution of Heat Transfer
The heat transfer equation can be modified through the use of a significant nondimensional temperature as: We propose new variable ξ that satisfies the correlation coefficients to determine the exact solution of the temperature Equation (9) for the stretching/shrinking sheet.
Substituting F(η) and Equation (24) in Equation (9) yields the subsequent: The boundary conditions reduces to: The singular point of the transformation (25) is ξ = 0. As a result, the Frobenius method is used to analyze the power solutions of Equation (25): After derivative of the first and second order Equation (27) w. r. t. ξ and upon combining the above-mentioned θ and θ expressions, Equation (25) is defined as: (n + r) a r ξ n+r−1 = 0 . (28) Micromachines 2022, 13, 1694 7 of 21 here applying the recurrence relation and the following values of n = 0, n = 1 − k 1 k 2 respectively: where: The solution of Equation (28) is simplified as follows using this recurrence relation, as well as the variables of n.
Finding the values of C 1 and C 2 , we get: The conjecture in terms of η is given by substituting the C 1 and C 2 values in Equation (30) and the incomplete gamma function defined.

Results and Discussion
The current study examines a ternary hybrid nanofluid flow towards a porous stretching surface using water as the base fluid. To analyze the current situation analytically, the volume fraction of aluminum oxide, single wall CNT, and graphene nanoparticles are all included in the equation. The current problem's ODEs were derived from controlling PDEs with similar variables that were solved analytically to provide the solution domain. The solution domain can be used to investigate solution domain, axial velocity, and heat transfer. The results of the current effort can then be examined using various physical characteristics and graphical arrangements. In this unsteady case, the porous medium parameter and thermal radiation play a major role in the stretching sheet problems. The fluid flow due to the porous medium is applicable in many industrial applications, such as the geophysical and allied areas. The results of three different nanoparticles can be examined in the present analysis; these ternary nanoparticles exhibit better thermal efficiency than that of normal nanofluids. Moreover, the couple stress fluid parameter is used to characterize the non-Newtonian fluid behavior. Figure 2 depict the impact of solution domain α on mass transpiration V c with various values of couple stress parameters for the stretching sheet. In this graph, we can see the results of four roots α 1 , α 2 , α 3 and α 4 . α 1 , α 2 are were found on the x-axis, whereas α 3 , α 4 were on the positive y-axis. These four roots are completely contrary to one another. For the roots α 1 , and α 4 , it can be seen that the solution domain is greater for higher values of couple stress. However, when the couple stress parameter values increase in the roots α 2 , and α 3 , then the solution domain decreases. The results of a modeling of dispersion water with three different types of nanoparticles are water-Al 2 O 3 /single wall carbon nanotubes/graphene when Pr = 6.27 and φ 1 = φ 2 = φ 3 = 0.1 and the other parameters remain constant. Similar effects apply in Figure 3 for different parameters, such as the inverse Darcy number.    In this graph, we can see the results of the four roots α 1 , α 2 , α 3 and α 4 . α 1 , α 2 were found on the x-axis, whereas α 3 , α 4 were on the positive y-axis. These four roots are completely contrary to one another. For the roots α 1 , and α 4 it can be seen that solution domain is greater for higher values of couple stress. However, when the couple stress parameter values increase in the roots α 2 , and α 3 then the solution domain decrease. Figure 5 represents the effect of f (η) on η for different values of the couple stress parameter Λ for the stretching surface, keeping other parameters constant. When couple stress parameter Λ increases, the thermal boundary decreases. In fact, the force and impact of the boundary conditions override the couple stress parameter Λ in the case of stretching sheets. An increase in the value of the couple stress parameter Λ causes a comparable rise in the axial velocity profiles in the case of both porous stretching and shrinking. The boundary layer becomes thinner when non-Newtonian viscoelastic shear stress occurs. Additionally, the boundary layer is smaller when suction, as opposed to injection, is utilized. Figure 6 represents the effect of f (η) on η for different values of the inverse Darcy number for stretching surface, keeping the other parameters constant. When the Darcy number increases, the thermal boundary decreases. The inverse Darcy number, as previously stated, is an attractive body force whose projection on the x-axis is situated in the opposite direction of x. Therefore, a higher inverse Darcy number value will result in more resistance to the axial velocity, thus reducing it. On the other hand, for a stretching sheet, the magnetic field and the imposed boundary condition are both in the same direction and will reinforce one another. Figure 7 represents the effect of f η (η) on η for different values of couple stress parameter Λ for the stretching surface, keeping other parameters constant. When couple stress parameter Λ increases, the thermal boundary decreases. Figure 8 represents the effect of f η (η) on η for different values of mass transpiration V c for the stretching surface, keeping other parameters constant. When mass transpiration V c increases, the thermal boundary decreases. The noticed increase in influence is due to the dimensionless number's simple index of the ratio of thermal resistances within the fluid, and at its surface, the mass transpiration is only significant close to the wall. Figure 9 represents the effect of f η (η) on η for different values of λ for the stretching surface, keeping the other parameters constant. When λ increases, thermal boundary increases. Figure 10 represents the effect of f η (η) on η for different values of the Darcy number for the stretching surface, keeping the other parameters constant. When the Darcy number increases, the thermal boundary increases. In fact, the BC's power and influence are more significant while stretching a sheet, unlike the results for the fluids. Similar results are obtained by increasing the size of the Darcy number. Axial velocity profiles in the case of both mass transpiration parameters and stretching/shrinking increase. The dividing mass transpiration in the stretching/shrinking sheet causes the layer to narrow.                    Figure 11 represents the effect of f η (η) on η for different values of volume fraction φ 1 for the stretching surface, keeping the other parameters constant. When φ 1 increases, the thermal boundary increases. The ternary hybrid nanoparticles within the nanofluid increase as the volume fraction of the nanoparticles increases, allowing for an increased area for improved heat conduction. As can be seen in Figure 11, this raises the velocity of the nanofluid. Figure 12 shows the impact of θ(η) on η for various values of Λ for the stretching surface. θ(η) decreases with an increase the value of Λ, and it is observed that the boundary layer thickness also reduces in the stretching surface. Figure 13 indicates the impact of θ(η) on η for various values of the Biot number for the stretching case. Here, it seen that θ(η) decreases as the value of the Biot number increases. As a result, the Biot number enhances the thermal diffusivity of a nanofluid for rising quantities of Bi. The nanofluids are increase in the efficiency of the Bi as the temperature rises. The performance of nanofluids is affected by both the drop in particle volume fraction and the increase in temperature in the flow pattern. Figure 14 show the temperature distribution θ(η) versus the similarity variable η for various values of the mass transpiration V c for the stretching surface. This means that the thickness of the thermal boundary layer under suction is thinner than under blowing. Figure 15 show the temperature distribution θ(η) versus the similarity variable η for various values of the thermal radiation N R for the stretching surface. Thermal radiation, therefore, increases the thermal diffusivity of nanofluids; for emerging radiation parameter values, heat will be added to the regime, and temperatures will rise as a result. A fluid temperature greater than both the wall temperature and the surrounding ambient temperature is technically possible, as was mentioned for the heat transmission of flows over a stretching sheet. Here, we have discussed the effect of forced flow over a stretching sheet; we now look at heat transport in the presence of radiation. The effect of heat conductivity is amplified by the radiation. Radiation has the effect of dampening or enhancing heat transmission in a linear manner. Figure 16a,b show the streamline flow patterns for suction and injection situations, respectively. For suction and injection scenarios, we look at different line patterns. The streamline flow pattern for the stretched boundaries is shown in Figure 16a,b for varied levels of mass transpiration V c while other parameters remain constant. When wall suction occurs at specific locations, as predicted by the physical theory, the flow field becomes regularized. On the other hand, blowing (flow injection) eliminates the streamline in the circular shape. Figure 17 indicates the effect of the skin friction coefficient on the stretching/shrinking surface λ as a function of couple stress Λ, respectively. Here, it is noticed that the couple stress parameter is the reciprocal of thickness decreasing of the boundary layer occurs via the same process. While skin friction is linearly proportional to the thickness of the boundary layer, it is estimated to have an inverse proportionality λ, as depicted in Figure 17. In addition, upward pulling couple stress parameters will draw fluid toward the negative y direction and reduce skin friction. The thickness of the boundary layer will decrease when the suction parameter, which measures the strength of wall suction, is increased, and as a result, the gradient of the velocity on the wall will also increase. This impact will enhance skin friction, as shown in Figure 17. Figure 18 shows the effect of the skin friction coefficient on couple stress Λ as a function of the stretching/shrinking surface λ, respectively. Here, it should be noted that the process of the decrease of the boundary layer's thickness and the stretching/shrinking surface parameter are the same. Skin friction is thought to have an inverse relationship with boundary layer thickness, as shown in Figure 18, even though the two are linearly proportional λ. Additionally, fluid will be drawn toward the negative y direction, and the skin friction will be decreased by the upward pulling couple stress parameters. The suction parameter, which gauges the strength of the wall suction, causes a reduction in the thickness of the boundary layer, which causes an increase in the gradient of velocity on the wall. The increase in skin friction caused by this collision is depicted in Figure 18.

Conclusions
The present study examines a ternary hybrid nanofluid flow towards a porous stretching surface using water as the base fluid. The current problem's ODEs were derived by controlling PDEs with similar variables that were solved analytically to provide the solution domain. The solution domain can be used to investigate the solution domain, axial velocity, and heat transfer. Some comparisons of our results are given in Table 2. The results of the current effort can then be analyzed using various physical characteristics and graphical arrangements.
The values of 1 α and 2 α increase as the value of Λ increases, while the values of

Conclusions
The present study examines a ternary hybrid nanofluid flow towards a porous stretching surface using water as the base fluid. The current problem's ODEs were derived by controlling PDEs with similar variables that were solved analytically to provide the solution domain. The solution domain can be used to investigate the solution domain, axial velocity, and heat transfer. Some comparisons of our results are given in Table 2. The results of the