MHD Mixed Convection Hybrid Nanoﬂuids Flow over a Permeable Moving Inclined Flat Plate in the Presence of Thermophoretic and Radiative Heat Flux Effects

: Recent nanotech advancements have created a tremendous platform for the development of a superior ultrahigh performance coolant referred to as nanoﬂuid for several industrial and engineering technologies. In this research, the impact of thermophoretic and viscous dissipation on the radiative mixed convective ﬂow comprising hybrid nanoﬂuid through an inclined permeable moving ﬂat plate with a magnetic ﬁeld is examined numerically. A model of non-linear differential equations is derived based on some realistic assumptions and tackled numerically using the bvp4c technique. The impact of the speciﬁc set of distinguished parameters on the velocity proﬁles, shear stress, temperature distribution proﬁles, heat transfer, concentration distribution proﬁle, and mass transfer for the two dissimilar branch solutions are discussed in detail. In addition, it has been discovered that double solutions exist in the case of an opposing ﬂow, while a single solution is observed in the case of an assisting ﬂow. The temperature distribution proﬁle escalates with the radiation parameter, while decelerating the velocity and concentration proﬁles.


Introduction
Nanotechnology is a rapidly expanding field in this modern age of advancement, where things are becoming smaller in size and getting better in the terminology of appropriate features. An area of particular importance is the preparation of products at the atomic and molecular levels for specific industrial purposes. One of the main ingredients of nanotechnology is nanofluid, which is primarily utilized to handle heat transfer challenges effectively. Thermal management of extremely sensitive systems is a major challenge in today's engineering and science operations, such as thermal power plants, nuclear reactions, and technologically advanced procedures. A nanofluid is formed by combining some regular fluids (kerosene oil, water, alcohol, etc.) with solid nanoparticles ranging in size from 1 to 100 nm. The credit for this large-scale advancement is given to the initial research of Choi and Eastman [1], where they developed the nanofluid concept. According to Eastman et al. [2], the deferment of Cu nanoparticles increases the conduction capacity of engine oil by 40%. The primary reason for this finding is that Cu has 3000 times the thermal conductivity of engine oil. Vajravelu et al. [3] elaborated the concept of the boundary layer (BL) flow of nanofluids using Cu and Ag as nanoparticles. Makinde and Aziz [4] investigated the behavior of nanofluids flowing past a stretchable sheet while taking the convective boundary condition into account. The time-dependent flow of a nanofluid near a stagnation point was inspected by Bachok et al. [5]. The unsteady/steady flow of a nanofluid through a movable surface in a consistent exterior free stream was investigated by Roşca and Pop [6]. Reddy and Chamkha [7] examined the Dufour and Soret effects on magneto flow through a stretchable sheet induced by water-based Al 2 O 3 and TiO 2 nanomaterials in porous media. The features of heat transfer for forced convection flow conveying nanomaterials through a moving plate immersed in a porous media with a heat sink/source were inspected by Ghosh and Mukhopadhyay [8]. Hussain et al. [9] presented a statistical analysis of the magnetized flow of a nanofluid through a stretchable sheet with a dissipation effect and also obtained the numerical solution. Waini et al. [10] utilized the Tiwari and Das model to scrutinize the Soret and Dufour effects on the flow of a nanofluid past a moving slim needle and presented double solutions. Hussain et al. [11] scrutinized the impact of solar radiation on the fluid flow of a non-Newtonian nanofluid immersed in a porous medium. Recently, Alazwari et al. [12] inspected the entropy generation of first-grade viscoelastic nanofluid past a porous plate. They observed that the entropy uplifts due to the volume fraction and Deborah number.
Several kinds of research have been carried out on nanofluids up to this point, but a binary hybrid nanofluid (HBN), one of the novel sorts of nanoliquid has recently piqued the interest of researchers. HBNs are made in two ways: (a) two or more dissimilar kinds of nanoparticles scattered in a regular or working base fluid, or (b) composite nanomaterials deferred in the regular fluid (conjoined). In fact, the application of nanofluid can improve heat transfer rates and lower production costs, which is why researchers are interested in this topic. The influences of convective heat transport through the tube employing a water-based Al 2 O 3 /Cu hybrid nanofluid were scrutinized by Suresh et al. [13] in which he discovered that hybrid nanofluids have a higher Nusselt number than regular fluids. Generally, the necessary cooling important for the electronic ingredients to avoid overheating is accomplished through the use of air-cooled heat sinks (ACHS) or liquid-CHS, with liquid-CHS being more suitable in terms of thermal efficiency. Selvakumar and Suresh [14] measured the pressure difference and rate of heat transport in an EHS using a water-based Al 2 O 3 /Cu hybrid nanofluid. According to their experimental outcomes, the coefficient of convective heat transport of a heat sink is augmented when hybrid nanoparticles are utilized as the working liquid. The features of heat loss through a riser pipe past a solar collector of the flat plate filled with a water-based hybrid nanofluid were inspected by Nasrin and Alim [15]. Devi and Devi [16] utilized the water-based hybrid nanofluid as well as a normal nanofluid to inspect the parametric characteristics of fluid flow past a porous stretchable sheet with a magnetic field. Ghadikolaei et al. [17] examined the transitive magnetic effect on the flow of water-based HBN past a stretching sheet with an ambient velocity and different kinds of shape factors. Huminic and Huminic [18] developed and intensely inspected the hybrid nanofluids which led to an enhancement of the thermal conductivity and ultimately an enhancement of heat transfer in the heat exchanger. The impact of the convective condition on the steady flow induced by HBN with heat transfer over a porous shrinking/stretching sheet was investigated by Waini et al. [19]. They observed that the rate of heat transfer is higher for HBN relative to a normal nanoliquid. Khan et al. [20] examined the mixed convective flow of HBN through a moving wedge with a magnetic effect and found double solutions. Roy and Pop [21] obtained an analytic solution of Stoke's second problem induced by a water-based Cu/Al 2 O 3 HBN with a heat source. The impact of a heat source/sink on magnetized flow through an exponentially stretchable sheet induced by a hybrid nanofluid was examined by Khan et al. [22]. They scattered alumina and copper nanoparticles in the blood base fluid. Recently, Khan et al. [23] discussed the axisymmetric flow and heat transport impinging on a permeable shrinkable/stretchable rotating disk conveying hybrid nanofluid and discovered a double branch outcome. Hussain [24] and Hussain et al. [25] examined the significance of a hybrid nanofluid with different aspects.
Thermophoresis refers to the phenomenon in which tiny micron-sized particles scattered in a non-isothermal gas accumulate velocity in the diminishing temperature direction. The velocity of gas molecules emerging from the warm end of the particles is higher than that of those emerging from the cold end. Collision of the molecules with the particles occurs quite vigorously, as they move more quickly. Because of the discrepancy in momentum, the particles develop a velocity in the cooler temperature direction. At first, Hales et al. [26] examined the impact of a thermophoretic assumption in the interest of engineering geometry. They inspected laminar equations through steam and aerosol transfer towards a vertical isothermal sheet. Chamkha and Pop [27] explored the role of thermophoresis particle deposition in a natural convection flow layer generated by a flat vertical plate in a porous media. Chamkha et al. [28] investigated the effect of aerosol particle thermophoresis in the laminar boundary layer through a vertical plate. Kandasamy et al. [29] analyzed the impact of irregular viscosity on magneto flow past a permeable wedge with a chemical reaction. The time-dependent flow of a micropolar fluid near a stagnation-point past a flat plate with thermophoresis and suction effects was examined by Zaib and Shafie [30]. Animasaun [31] inspected the Dufour and thermophoresis effects on magneto flow of a Casson fluid through a porous vertical plate with dissipation and temperature-dependent conductivity. Recently, Chamkha and Issa [32] scrutinized the 2D flow past a permeable semi-infinite surface subjected to a heat sink/source and thermophoresis effects.
Thermal radiation is critical in several engineering and industrial procedures such as nondestructive testing, electric power, panels of solar cells, and others. As a result, understanding the characteristic of thermal radiation is critical to achieving the preferred class of goods in industrial processes. The impact of radiation through the features of a non-Newtonian fluid with thermal diffusion and mixed convection was inspected by Mahmoud and Megahed [33]. Shehzad et al. [34] studied the thermal radiative impact on the 3D flow and heat transfer of a Jeffrey nanofluid with magnetism. Sheikholeslami and Rokni [35] investigated the roles of Coulomb force and thermal radiation in the trajectories of nanofluids in an enclosure of porous space. They discovered that the process of radiation aids in increasing the thickness of the thermal layer. Wakif et al. [36] explored the significance of the radiative effect on magneto flow of a hybrid nanofluid with effects of surface roughness. Recently, Khan et al. [37] discussed the impressions of porous media on radiative mixed convective flow induced by a hybrid nanofluid with an erratic heat sink/source from a vertical porous cylinder and presented double solutions.
As a result of the above reviews, the novelty of the current work lies in determining the numerical dual solutions of the MHD stagnation point flow induced by an HBN towards a porous moving inclined flat plate with mixed convection and thermal radiation. Moreover, the impact of viscous dissipation and thermophoretic effects are also incorporated which have not been investigated previously. The Tiwari-Das model is utilized to calculate the fine point of fluid flow. The hybrid nanofluid is created by suspending two distinct nanoparticles in the working fluids, namely silver (Ag) and titanium dioxide (TiO 2 ) (water). Using the similarity variables, the leadings equations and boundary conditions are later transformed into a system of ordinary differential equations, and then analysis is carried out in MATLAB using the bvp4c solver. The impacts of numerous values of the pertinent parameters on the physical quantities are demonstrated graphically, as well as in tabular form.

Mathematical Formulation
Consider a steady two-dimensional (2D) laminar flow of an incompressible electrically conducting hybrid nanofluid over a continuously moving semi-infinite inclined permeable flat plate with an acute angle Γ to the vertical, as shown in Figure 1, where, (x r , y r ) are the Cartesian coordinates with x r − and y r − being measured along the flat plate and normal to it, respectively. We assume that the velocities of the moving plate and the free stream are U a and U b , respectively (in other words, the plate moves away or toward the origin with constant velocities U a and ambient velocity U b ), while v w (x r ) is the mass flux velocity with v w (x r ) > 0 for suction and v w (x r ) < 0 for injection, respectively, T w is the constant temperature and C w is the constant concentration of the plate, while T ∞ and C ∞ are the temperature and concentration of the free stream, Q ra is the radiative heat flux (RHF) in the y r -direction. In addition, T w (x) < T ∞ represents the cooled plate, whereas T w (x) > T ∞ denotes the heated plate. For thermal enhancement, two different nanoparticles are considered, namely Ag (silver) and TiO 2 (titanium dioxide), diluted in the base working fluid (water) to form a novel class of hybrid nanofluid. A magnetic field (MF) B(x r ) is applied in the y r -direction. The following requisite assumptions are made: • The particles of concentration flux are adequately undersized such that neither the temperature field nor the main velocity stream are exaggerated via the thermo-physical processes practiced by the moderately small number of particles.

•
Owing to the behavior of the BL, the gradient of temperature is greater in the y rdirection than the x r -direction. As a result, only the component of thermophoretic velocity perpendicular to the surface is significant.

•
The particle diffusivity is assumed to be constant, and the particle concentration is low enough that particle clotting in the BL is assumed to be negligible.

•
The Bossiness approximation may be adopted for steady laminar flow.

•
The magnetic Reynolds number (MRN) is taken to be small enough that the IMF is insignificant in comparison to the AMF. r r r r flat plate and normal to it, respectively. We assume that the velocities of the moving plate and the free stream are a U and b U , respectively (in other words, the plate moves away or toward the origin with constant velocities a U and ambient velocity b U ), while   The particles of concentration flux are adequately undersized such that neither the temperature field nor the main velocity stream are exaggerated via the thermophysical processes practiced by the moderately small number of particles.  Owing to the behavior of the BL, the gradient of temperature is greater in the  r y direction than the  r x direction. As a result, only the component of thermophoretic velocity perpendicular to the surface is significant.  The particle diffusivity is assumed to be constant, and the particle concentration is low enough that particle clotting in the BL is assumed to be negligible.  The Bossiness approximation may be adopted for steady laminar flow.  The magnetic Reynolds number (MRN) is taken to be small enough that the IMF is insignificant in comparison to the AMF. The considered fluid is measured to be gray, representing an absorbing-emitting but non-scattering medium, and the Rosseland approximation (RAN) is exercised to describe a ra Q in the  r x direction that is measured to be negligible relative to the  r y direction. It is also assumed that     The considered fluid is measured to be gray, representing an absorbing-emitting but non-scattering medium, and the Rosseland approximation (RAN) is exercised to describe a Q ra in the x r -direction that is measured to be negligible relative to the y r -direction. It is also assumed that B(x r ) = B 0 / 2(x r /l r ) where B 0 is the constant strength of the magnetic field and l r is the characteristic length of the plate. Under the above stated conditions, the steady continuity, the momentum, the energy, and the concentration boundary layer equations of the hybrid nanofluid, in Cartesian coordinates (x r , y r ), can be written as (see Noor et al. [38]): Mathematics 2022, 10, 1164 subject to the boundary conditions where u r and v r are the velocity components along (x r , y r )-axes, T r the temperature, and C r the concentration of the hybrid nanofluids, g r is the acceleration due to gravity, D B the molecular diffusivity of the species concentration, and V TM the thermophoretic velocity.
Further, k hbn f is the thermal conductivity, ρ hbn f is the density, σ hbn f is the electrical conductivity, µ hbn f is the absolute viscosity, ρc p hbn f is the heat capacity, and (ρβ) hbn f is the coefficient of thermal expansion, which are given by (see Ho et al. [41]; Sheremet et al. [42]).
ρc p hbn f (ρβ) hbn f In Equations (6)-(10), the symbol φ is called the solid nanoparticles volume fraction (where φ = 0 corresponds to a regular (viscous) fluid) and D A represents the shape factor parameter, and here the value of the shape factor is considered as D A = 3 (spherical). Moreover, the other notations such as are the respective electrical conductivity, thermal conductivity, viscosity, specific heat capacity, density, and thermal expansion coefficient of the base working fluid while φ a signifies the Ag (silver) nanoparticles and φ b denotes the TiO 2 (titanium dioxide) nanoparticles. The physical data of the base working fluid and the two distinct hybrid nanoparticles (Ag (silver) and TiO 2 (titanium dioxide)) are written in Table 1. Now, employing the RAN, the Q ra (RHF) can simply take the following form (see Bataller [43]; Ishak [44]; Magyari and Pantokratoras [45]); where k A and σ A signify the mean absorption coefficient and the corresponding Stefan-Boltzmann constant, respectively. In addition, by applying a Taylor series to the fourth power of T 4 r about the point T ∞ , and ignoring the higher power order terms, a simplified form can be obtained as T 4 Using this, Equation (3) can be re-written as: The thermophoretic velocity V TM which appears in Equation (4) can be written as where T b is a reference temperature and k is the thermophoretic coefficient with a range of values from 0.2 to 1.2 as designated by Batchelor and Shen [46]. Therefore, Equation (1) is fulfilled for the choice of stream function ψ(x r , y r ) defined as u r = ∂ψ/∂y r and v r = −∂ψ/∂x r .
Guided by the boundary conditions (5), we introduce the following similarity variables: Thus, the velocity components (u r , v r ) are given by Here, the prime symbolizes the relative derivative with respect to ξ, f w is the constant mass flux velocity. Therefore, the cases f w > 0 for suction and f w < 0 for injection are observed, respectively.
Substituting Equation (14) into Equations (2), (4) and (12) generates the following system of ordinary (similarity) coupled equations as follows: Mathematics 2022, 10, 1164 with associated boundary conditions: To obtain the similarity solution of Equation (3), we presume the required case is given In addition, Pr is the Prandtl number, Sc a is the Schmidt number, Ec a is the Eckert number, Λ a is the magnetic parameter, R ad is the radiation parameter, δ a is the thermophoretic parameter, and Σ a is the mixed convection parameter, which are defined as with as the local Grashof number and Re as the local Reynolds number. It should be noticed that Σ a > 0 corresponds to assisting flow, Σ a < 0 corresponds to opposing flow, and Σ a = 0 corresponds to forced convection flow. In addition, ε r is the constant moving plate parameter, with ε r > 0 for the moving plate in the up direction, ε r < 0 for the plate moving in the down direction, and ε r = 0 for the static plate.

The Engineering Interest Quantities
The gradients of practical interest are the local shear stress C f , local heat transfer Nu x r , and the local mass transfer Sh x r , which are defined as: Substituting Equation (14) into the aforementioned Equation (22), the following obtained dimensionless form of the gradients can be written as: where Re is the local Reynolds number.

Numerical Solution Procedure and Validation of the Scheme
This section discusses mainly the complete solution procedure of the scheme as well as the confirmation of the code. The problem is initially bounded in the form of PDEs and then changed into ODEs via executing the similarity variables. Consequently, the given model accomplished the similarity dual solutions (first solution and second solution) of Equations (17)- (19) along with BCs (20) via executing a bvp4c built-in package working in the MATLAB program, and the specifics are reported in Shampine et al. [48]. However, this package further uses a finite difference scheme with fourth-order precision that incorporates the three-stage Lobatto IIIA formula. For the working system of the technique, the set of similarity ODEs from higher second and third-order are altered into first-order ODEs via introducing new variables. Let the new variables be as follows: By making use of the new variables in the requisite posited similarity equations, the required set of first-order ODEs can be written as: The following known and unknown ICs are used to solve the above set of Equation (25) as a one-point boundary value problem (IVP): Furthermore, the arbitrary constants E 1 , E 2 , and E 3 are guessed during the process of bvp4c simulations and the comprised parameters are also set accordingly to execute the acceptable solution. For the acceptable outcome, we can estimate the value of the aforementioned constants in such a way that the working process of iteration is repeated until, once satisfied, the subsequent BCs are as follows: To put it another way, the computational solution is deemed reasonable when no warnings or errors are generated during execution and the far-field BCs (27) are met. In addition, the present problem consists of more than one (first and second) solution. Therefore, the code needed two dissimilar choices or guesses for the output. The guess or estimate for the first solution is very simple while for the second solution the guess is quite complex to appropriately select. Meanwhile, the guess needed more time and effort for the second branch solution which satisfies the far-field BCs asymptotically. First and foremost, the validation technique is carried out to ensure that the current model is reliable. In this respect, a comparative analysis was conducted towards the solo numerical value of the shear stress (|F (0)|) with the available published works of Cortell [49], Yazdi et al. [50], and Javed et al. [51] for the FBS due to the limiting cases of the stretching parameter ε r = 1.0 and without the influence of the stagnation point (u r → 0 as y r → ∞), f w = 0, Γ = 0, Λ a = 0, Σ a = 0, φ a = 0, and φ b = 0. The comparison is shown computationally in Table 2, where the solutions are exceptionally well-matched and give assurance on the considered scheme applied in this problem. Hence, this gives further confidence and power that the unavailable results of the existing problem are new and original. Table 2. Comparison of the shear stress (|F (0)|) value for the first branch solution (FBS) due to the limiting case of the stretching parameter ε r = 1.0 when f w = 0, Γ = 0, Λ a = 0, Σ a = 0, φ a = 0, φ b = 0, and without the influences of stagnation point (u r → 0 as y r → ∞).

Graphical Results and Discussion
In this portion of work, we need to discuss the physical interpretation of the graphical and tabular outcomes of the (Ag-TiO 2 /water) hybrid nanoparticles for the two dissimilar branch solutions due to the influence of the various influential embedded flow parameters. Tables 1 and 2 signify the thermo-physical data of the hybrid nanoparticles and the comparison of the considered scheme. Furthermore, Table 3 explicates the values of the quantity C f (i.e., the shear stress) towards the variation in the suitable flow parameters for the first branch solution (FBS) as well as the second branch solution (SBS) when ε r = −2.5 and Γ = π/6. The outcomes reveal that the shear stress for the FBS enriches owing to the larger values of f w , Λ a , and Σ a , while it is reduced for the SBS. Meanwhile, the impacts of the hybrid nanoparticles can boost the behavior of the quantity C f (i.e., the shear stress) for both dissimilar branch solutions. Table 4 clarifies the computational values of the quantity Nu x r (i.e., the heat transfer rate) owing to the variation in several flow parameters for the first and second branch solutions when ε r = −2.5, f w = 2.0, Λ a = 0.05, and Γ = π/6. It is observed that the magnitude of the heat transfer upsurges for the FBS and for the SBS due to the superior values of Σ a . Meanwhile, the escalating values of the hybrid nanoparticles, the radiation parameter, and the Eckert number shrink the heat transfer behavior for both solution branches. Similarly, Table 5 gives the quantitative values of the quantity Sh x r (i.e., the Sherwood number) due to the varying values of the different flow parameters for the FBS as well as for the SBS when ε r = −2.5, f w = 2.0, Λ a = 0.05, R ad = 1.0, Ec a = 0.10, Σ a = −1.0, and Γ = π/6. From the results, it can be observed that the mass transfer enhances the FBS and shrinks for the SBS due to the mounting values of φ a , φ b . Alternatively, the upsurging role of the parameters Sc a and δ a monotonically escalates the mass transfer rate for both solutions (FBS and SBS). However, the shear stress is not directly affected by Sc a and δ a because these parameters appeared in the concentration equation and play a key role in the energy equation due to the presence of thermophoretic values.   Figures 2-4 elucidate the variation in the C f , Nu x r , and Sh x r due to the larger values of Λ a against Σ a , respectively. Multiple branch (FBS and SBS) solutions are accomplished in all the physical quantities for the only case of the buoyancy opposing flow (Σ a < 0), while in the case of the buoyancy assisting flow (Σ a > 0), only a single solution is possible. Further, these graphs illustrate the values of C f and Sh x r increase for the FBS and decline for the SBS due to the larger impact of Λ a , while the Nu x r abruptly shrinks for both solution branches. In the physical scenario, the larger influence of the magnetic parameter Λ a can create a strong drag force which shrinks the motion of the fluid flow of the hybrid nanoparticles. Therefore, the motion of the fluid flow is inversely related to the friction of the fluid. As a result, the shear stress upsurges because of the larger impacts of the magnetic parameter. Since multiple solutions are possible here only for the case of the buoyancy opposing flow (Σ a < 0) where the FBS represents the physically stable or trustworthy solution and the SBS corresponds to the unstable (not physically reliable) solution. Moreover, the unique solution is observed in these graphs for the case when Σ a = Σ ac , no solution is seen for the case when Σ a < Σ ac , and the possible outcomes are detected for the case when Σ a > Σ ac . According to the graphs, it is observable that for each distinct value of Λ a , the following change critical values are found (see  and are written in each window of the graphs. These critical points are highlighted by the small solid balls that the readers can easily see. Meanwhile, the magnitude of the critical values is also enriched monotonically with growing values of Λ a . This tendency further demonstrates that the boundary layer separation diminishes due to the influences of the larger magnetic parameter. a cally stable or trustworthy solution and the SBS corresponds to the unstable (not physically reliable) solution. Moreover, the unique solution is observed in these graphs for the case when    a ac , no solution is seen for the case when    a ac , and the possible outcomes are detected for the case when    a ac . According to the graphs, it is observable that for each distinct value of  a , the following change critical values are found (see  and are written in each window of the graphs. These critical points are highlighted by the small solid balls that the readers can easily see. Meanwhile, the magnitude of the critical values is also enriched monotonically with growing values of  a . This tendency further demonstrates that the boundary layer separation diminishes due to the influences of the larger magnetic parameter.      The impacts of the radiation parameter R ad on the physical quantities such as C f , Nu x r , and Sh x r of the (Ag-TiO 2 /water) hybrid nanoparticles for the FBS as well as the SBS towards the buoyancy parameter Σ a are graphically portrayed in Figures 5-7. It is observed that the C f declines in the FBS but inclines in the SBS due to the augmentation in the values of R ad , while Nu x r diminishes in the outcomes of the FBS as well as in the SBS due to the larger influence of the radiation parameter. In general, the larger impact of R ad gives a lower quantity of heating to the adjustable (Ag-TiO 2 /water) hybrid nanofluid owing to the presence of the Eckert number in the similarity Equation (18). As a result, the heat transfer decelerates with the superior influences of R ad . On the other hand,Sh x r initially improves for the FBS and then significantly starts decaying for the same branch as R ad increases. Meanwhile, the behavior of the solution continuously develop higher and higher for the SBS owing to the larger impressions of R ad . For a particular domain of, unique solution is observed when, no solution when, and two solutions when. Moreover, the obtained critical values are Σ ac = −1.7383, −1.3687, and −1.1164 for the corresponding dissimilar values of R ad (=0, 0.5 and 1.0), respectively. From these outputs, it can be seen that the behavior of the magnitude of the critical values decreases if we escalate the values of R ad . Hence, the current pattern suggests that the separation of the boundary layer augments with larger values of R ad . Figure 8 expounds the variations in the thermophoretic parameter δ a on the local Sherwood number Sh x r of the (Ag-TiO 2 /water) hybrid nanoparticles for both solution branches versus Σ a . For rising values of δ a , the mass transfer rate enriches for the FBS as well as for the SBS, while the gap between the outcomes looks similar in both branches with δ a . Generally, this behavior is due to the fact that higher consequences of the parameter δ a can improve the thermophoretic coefficient k; as a result, the mass transfer rate would be significantly developed. Moreover, higher values of the thermophoretic parameter δ a can give us a single bifurcation value as shown in the graph (see Figure 8).        rate would be significantly developed. Moreover, higher values of the thermophoretic parameter  a can give us a single bifurcation value as shown in the graph (see Figure 8). of the (Ag-TiO2/water) hybrid nanoparticles for the two dissimilar branch solutions due to the greater values of  a against the pseudo-similarity parameter  , respectively. Further, to deal with these graphical outcomes, it is observed that the velocity and concentration distribution profiles increase for the FBS and falloffs the curves for the SBS due to the larger values of  a . Moreover, this dynamic growth of the (Ag-TiO2/water) hybrid nanoparticles' fluid flow for the FBS indicates that the velocity and concentration boundary layer become progressively thinner specifying that the magnetic parameter has a sensational influence and supports the motion of the hybrid nanoparticles near the surface of Figures 9-11 display the variation in the F (ξ), G(ξ), and S(ξ) of the (Ag-TiO 2 /water) hybrid nanoparticles for the two dissimilar branch solutions due to the greater values of Λ a against the pseudo-similarity parameter ξ, respectively. Further, to deal with these graphical outcomes, it is observed that the velocity and concentration distribution profiles increase for the FBS and falloffs the curves for the SBS due to the larger values of Λ a . Moreover, this dynamic growth of the (Ag-TiO 2 /water) hybrid nanoparticles' fluid flow for the FBS indicates that the velocity and concentration boundary layer become progressively thinner specifying that the magnetic parameter has a sensational influence and supports the motion of the hybrid nanoparticles near the surface of the thermophoretic flat plate. On the other hand, the temperature distribution behaves differently in both solution branches as compared to the outcomes of F (ξ) and S(ξ) for the larger impact of the magnetic parameter. As a general explanation, the larger magnetic parameter produces the conception of the Lorentz forces which creates a noticeable boost in the motion of the fluid and ultimately affects or decreases the temperature distribution.     The impacts of the radiation parameter ad R on the profiles of the (Ag-TiO2/water) hybrid nanoparticles for the FBS as well as for the SBS are graphically established in Figures 12-14, respectively. The radiation parameter has a significant influence on the velocity, concentration, and temperature distributions. The rise in the radiation parameter accelerates the temperature but declines the fluid velocity and concentration for the FBS, while the tendency for the SBS is the opposite. Moreover, the thickness of the thermal boundary layer enriches with higher radiation parameters and with lower velocity or momentum and concentration boundary layers. Physically, the maximum impacts of radiation cause a considerable quantity of heating to the (Ag-TiO2/water) hybrid nanoparticles that augments the temperature distribution as well as the thickness of the thermal boundary layer. The impacts of the radiation parameter R ad on the profiles F (ξ), G(ξ), and S(ξ) of the (Ag-TiO 2 /water) hybrid nanoparticles for the FBS as well as for the SBS are graphically established in Figures 12-14, respectively. The radiation parameter has a significant influence on the velocity, concentration, and temperature distributions. The rise in the radiation parameter accelerates the temperature but declines the fluid velocity and concentration for the FBS, while the tendency for the SBS is the opposite. Moreover, the thickness of the thermal boundary layer enriches with higher radiation parameters and with lower velocity or momentum and concentration boundary layers. Physically, the maximum impacts of radiation cause a considerable quantity of heating to the (Ag-TiO 2 /water) hybrid nanoparticles that augments the temperature distribution as well as the thickness of the thermal boundary layer.

Conclusions
The impact of thermophoretic and radiative heat flux on MHD mixed convection water-based hybrid Ag-TiO2 nanofluid flow over a permeable moving inclined flat plate with uniform free-stream velocity has been studied numerically. It has been found that for a specific set of embedded influential flow parameters, the shear stress, velocity profiles, heat transfer, temperature distribution profiles, concentration distribution, and the mass transfer reveal multiple (FBS and SBS) solutions over a broad range of the mixed or buoyancy convection parameter. The main findings are concluded as follows:


The mass transfer and the shear stress enlarge owing to the higher magnetic parameter while the heat transfer reduces.  Increasing the solid volume fraction of silver and titanium dioxide nanoparticles enhances the shear stress and mass transfer but reduces the heat transfer.

Conclusions
The impact of thermophoretic and radiative heat flux on MHD mixed convection water-based hybrid Ag-TiO 2 nanofluid flow over a permeable moving inclined flat plate with uniform free-stream velocity has been studied numerically. It has been found that for a specific set of embedded influential flow parameters, the shear stress, velocity profiles, heat transfer, temperature distribution profiles, concentration distribution, and the mass transfer reveal multiple (FBS and SBS) solutions over a broad range of the mixed or buoyancy convection parameter. The main findings are concluded as follows:

•
The mass transfer and the shear stress enlarge owing to the higher magnetic parameter while the heat transfer reduces.

•
Increasing the solid volume fraction of silver and titanium dioxide nanoparticles enhances the shear stress and mass transfer but reduces the heat transfer.

•
The augmented value of the thermophoretic parameter develops the mass transfer while the heat transfer rate shrinks with Eckert number.

•
The magnitude of the critical values surges due to the higher influences of the magnetic parameter but reduces with the radiation parameter.

•
The higher influence of the thermophoretic parameter enriches the mass transfer rate.

•
The boundary layer separation enhances with radiation, but shrinks owing to the Lorentz forces.

•
The velocity and concentration decelerate with the radiation parameter, while the temperature significantly increases.

Acknowledgments:
The financial support received from the Universiti Kebangsaan Malaysia (Project Code: DIP-2020-001) and the Universiti Teknikal Malaysia Melaka is gratefully acknowledged.

Conflicts of Interest:
The authors declare no conflict of interest. Titanium dioxide a, b

Nomenclature
Solid nanoparticles b f Base fluid hbn f Hybrid nanofluid w Wall boundary condition ∞ Far-field condition ' Derivative with respect to ξ