Ramiﬁcation of Hall and Mixed Convective Radiative Flow towards a Stagnation Point into the Motion of Water Conveying Alumina Nanoparticles Past a Flat Vertical Plate with a Convective Boundary Condition: The Case of Non-Newtonian Williamson Fluid

: Heat transfer technologies are experiencing rapid expansion as a result of the demand for efﬁcient heating and cooling systems in the automotive, chemical, and aerospace industries. Therefore, the current study peruses an inspection of mixed convective radiative Williamson ﬂow close to a stagnation point aggravated by a single nanoparticle (alumina) from a vertical ﬂat plate with the impact of Hall. The convective heating of water conveying alumina (Al 2 O 3 ) nanoparticles, as appropriate in engineering or industry, is investigated. Using pertinent similarity variables, the dominating equations are non-dimensionalized, and after that, via the bvp4c solver, they are numerically solved. We extensively explore the effects of many relevant parameters on axial velocity, transverse velocity, temperature proﬁle, heat transfer, and drag force. In the opposing ﬂow, there are two solutions seen; in the aiding ﬂow, just one solution is found. In addition, the results designate that, due to nanoﬂuid, the thickness of the velocity boundary layer decreases, and the thermal boundary layer width upsurges. The gradients for the branch of stable outcome escalate due to a higher Weissenberg parameter, while they decline for the branch of lower outcomes. Moreover, a magnetic ﬁeld can be used to inﬂuence the ﬂow and the properties of heat transfer.


Introduction
As research of subject, the mechanics of shear-thinning (pseudoplastic) fluids has grown in popularity because of its enormous significance in biosciences, engineering, and industrial sectors, such as in extruded polymer sheets, high melts of natural polymers, and sheets coated with emulsion-like photographic films. To examine the rheological characteristics of all varieties of pseudoplastic fluids, numerous models have been developed (paints, blood, honey, etc.). The Cross model, Power-law model, Carreau model, Ellis model, and Williamson model are a few of them. For further research, we have taken into account the model of Williamson fluid in this analysis. The advantage of this model is that it takes into account both minimum µ 0 (shear stress goes to zero) and maximum µ ∞ (shear stress goes to infinity) viscosities. Williamson [1] considered a model for analyzing the behavior of pseudoplastic flows in 1929 and provided experimental support for it. He discussed the significance of shear-thinning flows and pointed out how they differ greatly from viscous flows. In their study of the Williamson fluid flow across a stretchy sheet, Nadeem et al. [2] found that the velocity distribution and friction factor degrade by increasing the Williamson parameter. Hayat et al. [3] explored the Williamson fluid passing via a stretchy sheet and looked at the degradation of wall shear stress for a greater amount of Williamson fluid when heat radiation and a magnetic field were present. The features of mass and heat transfer comprising Williamson fluid across a permeable stretchy sheet subject to a magnetic field were inspected by Shawky et al. [4]. The rate at which heat is transferred to the stretched surface increases as the Williamson and porosity parameters are increased. Significant studies regarding distinct aspects of Williamson fluid are found in references [5][6][7][8].
Nanofluids have been widely used in industry. Nanofluid, characterized by a substantial increase in molecular mobility compared to unadventurous engineered fluid [9], is found to serve in many practical applications, such as porous materials [10,11], and has drawn particular attention for encapsulation in the medical industry in recent years. The addition of nanoparticles or metal oxide is a renowned tactic for modifying the heat transfer, as well as the thermal conductivity, of fluids. Metals are the most valuable nanoparticles due to their greater thermal conductivity correlation than traditional fluids [12]. By suspending a specific kind of nanoparticle in working fluids, a persistent fluid, known as a nanofluid, is created. When it comes to heat transfer, nanofluids usually perform better than more common fluids such as water, oil, and ethylene glycol. Due to their superior heat transmission qualities, nanofluids have attracted a lot of attention. The inherent difficulty of heat transmission is the most well-known concern with heating systems. A better solvent that is less expensive and more widely available is water. As a result, water is frequently utilized in industry for the dispensation of nano-liquids due to its weak thermal conductivity (TCN), which requires perfection if employed as a heat exchanger by dispersion of nanomaterials [13]. Masuda et al. [14] demonstrated that the viscosity and TCN of common fluids can be altered by the addition of nanoparticles. Choi and Eastman [15] explored that the TCN of fluid substantially increases when the nanoparticles are detached in regular fluids. Beg et al. [16] utilized Soret and Dufour impacts to achieve a non-identical solution of buoyancy flow over an inclined surface. They observed that the fluid flow is increased due to the negative inclination of the plate, whilst the reverse trend is observed for the positive inclination. Anbuchezhian et al. [17] investigated the impact of temperature stratification on the flow of a nanofluid across a stretchy sheet when powered by solar radiation. It has been demonstrated that base liquids' thermal conductivity and convective heat transfer abilities are improved by nanofluids. Nasrin and Alim [18] scrutinized free convective flow and heat transport-incorporated nanofluids in a wavy solar collector. According to Kandasamy et al. [19], solar radiation submerged in a non-Darcy porous medium caused an erratic magnetic flow of a nanofluid across a permeable wedge. They have shown that the copper nanoparticle plays a vital role in absorbing solar radiation. By using solar radiation, Khan et al. [20] calculated the numerical outcomes of a 3D flow via a bi-directional exponentially stretchy surface that was filled with nanofluid. The model for solar energy was created by Shehzad et al. [21] by taking into account a non-Newtonian nanofluid from a stretchy sheet. Zeeshan and Majeed [22] investigated the buoyancy flow in a non-Darcy porous media using a stretchy sheet and magnetic ferroparticles. They concluded that the frictional reluctance of a Lorentz force diminishes with a higher electric field, whereas the opposite behavior is seen for temperature. Madhukesh et al. [23] examined the flow with heat transfer stimulated through water-based AA7072/AA7075 hybrid nanofluids past a stretchable curved surface. The heat transfer is improved due to thermal relaxation and Newtonian heating parameters. The flow of hybrid nanofluid scattered ferrite particles through a contraction/expansion rotating disc was inspected by Gowda et al. [24]. More about the significance of nanofluids can be observed in references [25][26][27]. Recently, Animasaun et al. [28] considered the impact of the magnetic field inspired by three kinds of nanofluids through a heated convective surface with a source/sink. A rise in the wall's convective heating is a factor that can improve the temperature distribution.
Magnetohydrodynamic (MHD) electrically conducting fluid is widely used in many engineering and industrial processes. Flow meters, MHD power plants, nuclear reactors, stirring, and MHD pumps are a few examples. This is also relevant to how the solar cycle and sunspots grow in terms of solar physics. The idea of the magnetic flow of an electrically conducting liquid from a flexible sheet was first discussed by Pavlov [29]. Many researchers have indeed discussed the electrically conducting fluid transient from a stretchy surface with a magnetic field [30][31][32]. On the other hand, in the rarefied medium, with a low density or a strong magnetic field, the fluid conductivity is anisotropic and the Hall current cannot be ignored. The Hall current plays a vital role in the Hall accelerator problem and flight MHD. Due to the low density of plasma in the magnetosphere and the solar wind, the Hall influence may be significant during the restoration of the magnetic field near the Earth's magnetopause [33]. Su and Zheng [34] looked at how velocity slip affected the magneto flow of a nanofluid while transporting heat through a stretchy wedge using Joule heating. Zaib and Shafie [35] examined the mutual influences of MHD and Hall current on viscous, dissipative, time-dependent free convective flow via a stretchy sheet. Hall current and radiation absorption were used by Sreedevi et al. [36] to evaluate the magneto free convective flow from a stretchy surface. The magneto flow with heat transport of nanoparticles over a nonlinear stretchable/shrinking sheet was investigated by Pal and Mandal [37] using Hall current (HC). The HC impact on the radiative flow included by nanofluid across a revolving disc was prompted by Khan et al. [38]. In addition to using sensitivity analysis, Rana et al. [39] investigated the importance of Hall current, activation energy, and nanoparticles in enhancing the phenomenon of heat transfer through a nanofluid. A water-based hybrid nanofluid's quadratic buoyancy convection flow was recently studied by Rana and Gupta [40] using the surface of a rotating cone with Hall current.
Numerous engineering and industrial processes use the thermal radiation impact with convective boundary conditions in die forging, storage of thermal energy, chemical reactions, and nuclear turbines. Aziz [41] inspected the viscous convective flow across a plate under convectively heated conditions. Makinde and Aziz [42] prompted the magnetic impact on buoyancy flow from a vertically heated plate saturated in a porous medium. Ishak [43] used numerical analysis to investigate the flow and heat transfer toward a stretchy, highly permeable, convectively heated sheet. The specific issue of viscous flow through a penetrable, convectively heated, elongating wall was resolved by Yao et al. [44]. Rahman et al. [45] provoked the convective heat transport phenomenon through a vertical plate. Maxwell fluid in a nanofluid was assessed by Mustafa et al. [46] using an exponentially stretchy sheet under convective conditions. Ibrahim and Haq [47] examined the magnetic flow of a nanofluid past a convectively heated stretched sheet as it travelled toward an SP. Recently, Makinde et al. [48] examined the effects of slip and irregular radiation on magneto flow-integrated nanoparticles near an SP toward a connective heated stretchy sheet.
The novelty of this work is to propose the mutual impacts of MHD and Hall current on the mixed convection Williamson flow induced by nanofluid over a heated vertical flat plate with a convective boundary condition and irregular radiation. Nanofluid is basically the combination of water-based and alumina nanoparticles which are also investigated in the current problem. The numerical, transmogrified, nonlinear ODEs are solved using the bvp4c solver. One of the goals of this challenge is to find multiple solutions for the scenario of opposing flow. According to the available information, no one has yet considered this particular problem.

Problem Formulation
The magneto-mixed convective Williamson fluid flow via a stagnation point stimulated by nanoparticles via a vertical flat plate with thermal radiation and Hall current effects is considered. The convective boundary condition is also examined. The problem of nanofluid flow is schematically shown in Figure 1, where the positive x− and y−coordinates deliberate along and normal to the flat plate in the trend of motion and toward the fluid, respectively, while the z−axis coordinate is transverse to the xy−plane. There is a significant external magnetic impact B 0 in the corresponding y−direction. Due to the exceedingly low persuasive Reynolds number, the induced magnetic influence is also disregarded. The Hall impact cannot be disregarded because the strength of electron-atom collisions is assumed to be quite strong [49]. The flow becomes three-dimensional as the Hall current becomes stronger in the respective z−direction, causing a crossflow in that direction. To get around this issue, we presupposed that flow amounts in the z−direction had not changed. The flat plate must have infinite lengths in that direction for this to be true. It is also believed that the plate is electrically non-conducting and that J y = 0 is applied uniformly across the flow according to the generalized Ohms' law [50]. The governing equations for controlling the flow, along with Boussinesq approximations, are provided as per these presumptions [35,37].
∂u ∂x subject to the boundary conditions are Lubricants 2022, 10, x FOR PEER REVIEW 4 of 21 flat plate with a convective boundary condition and irregular radiation. Nanofluid is basically the combination of water-based and alumina nanoparticles which are also investigated in the current problem. The numerical, transmogrified, nonlinear ODEs are solved using the bvp4c solver. One of the goals of this challenge is to find multiple solutions for the scenario of opposing flow. According to the available information, no one has yet considered this particular problem.

Problem Formulation
The magneto-mixed convective Williamson fluid flow via a stagnation point stimulated by nanoparticles via a vertical flat plate with thermal radiation and Hall current effects is considered. The convective boundary condition is also examined. The problem of nanofluid flow is schematically shown in Figure 1, where the positive x − and y − coordinates deliberate along and normal to the flat plate in the trend of motion and toward the fluid, respectively, while the z − axis coordinate is transverse to the xy − plane. There is a significant external magnetic impact 0 B in the corresponding y − direction. Due to the exceedingly low persuasive Reynolds number, the induced magnetic influence is also disregarded. The Hall impact cannot be disregarded because the strength of electron-atom collisions is assumed to be quite strong [49]. The flow becomes three-dimensional as the Hall current becomes stronger in the respective z − direction, causing a crossflow in that direction. To get around this issue, we presupposed that flow amounts in the z − direction had not changed. The flat plate must have infinite lengths in that direction for this to be true. It is also believed that the plate is electrically non-conducting and that 0 y J = is applied uniformly across the flow according to the generalized Ohms' law [50].The governing equations for controlling the flow, along with Boussinesq approximations, are provided as per these presumptions [35,37].   Fourier's law of heat conduction, which is expressed in differential form by Equation (5), asserts that the rate of heat transmission through a medium is inversely proportional to the negative gradient in the temperature per area. Further, in these above equations, u, w, v are the components of velocity in the x−, z−, and y−directions, respectively, m is the Hall current parameter, Γ is the time constant, µ n f is the viscosity of the nanofluid, ρ n f is the density of the nanofluid (NF), T and T ∞ signify the temperature and the free-stream temperature, q r the radiative heat-flux, h f the coefficient of heat transfer, σ n f the nanofluid EC, k n f is the NF thermal conductivity, g the gravitational acceleration, and (ρc p ) n f the NF heat capacity. These are demonstrated as [23,24,30,51]: where k s1 , k f signpost the thermal conductivity of the solid nanoparticles and the regularbased fluid (water), respectively. Additionally, ϕ signifies the solid nanoparticles volume fraction, which is here the sum of the single nanoparticles and the regular fluid (water). In addition, β s1 , β f designate the coefficient of thermal expansion of the nanoparticles and base fluid respectively, σ s1 , σ f signpost the electrical conductivity of the solid nanoparticles and the basic fluid, and ρ s1 and ρ f signify the density of the nanoparticle and base fluid, respectively. Moreover, the remaining subscripts in the above equations defined as f and n f are called the based fluid and the nanofluid. The thermophysical data of the base fluid (water) and the (Al 2 O 3 /water) nanofluid are written in Table 1. The radiative heat-flux q r is articulated as [26]: where k * and σ * indicate the absorption coefficient and the Stefan-Boltzmann constant, respectively. Utilizing Equation (12), Equation (4) can be written as To ease the analysis of the given nanofluid problem, the following self-similarity variables are introduced [35,39]: where ψ presents the stream function and it is distinguished as u = ∂ψ ∂y and v = − ∂ψ ∂x . With the aid of Equation (14), the nonlinear Equations (2), (3), and (13) yield the following nonlinear ODEs: 1 and the transformed subjected boundary conditions are: where ( ) denote differentiation with respect to the pseudo-similarity variable η. In addition, the aforementioned similarity equations comprised different influential controlling parameters which are symbolically and namely written as follows: We a = √ 2aΓ √ Re x , and Bi N = h f υ f /a/k f imply the magnetic parameter, the mixed convective parameter (which is the ratio of the Grash of number and the square of the local Reynolds number), the Prandtl number, the radiation parameter, the Weissenberg parameter, and the convective parameter, respectively.
The important physical quantities of interest are the skin frictions in axial and transverse directions and the local Nusselt Number, which is stated as [35] where τ w and m w are the shear stresses for the Williamson nanofluid flow in the axial and transverse directions, respectively. Meanwhile, q w is the heat flux, which is stated as: , and Using (14), we get where Re x = u e (x)x/v f is called the local Reynolds number.

Results and Discussion
This section illustrates the analysis of the solutions both graphically and in the form of several numerical tables for the branches of unstable and stable solution. Thus, to understand the behavior of the dual results for a limiting case, we first discuss in detail the numerical scheme, as well as the code reliability and validations for the phenomenon of available results. The distorted system of nonlinear Equations (15)-(17) and (18) were eased numerically via an efficient bvp4c solver for dissimilar values of pertaining parameters that were apprehensive in the ensuing given the (Al 2 O 3 /water) nanofluid flow problem. According to Oztop and Abu-Nada [53], the value of the Prandtl number Pr = 6.2 is taken into account throughout the investigation and the volume fraction of nanoparticles ranges from 0 to 0.2, which ϕ = 0 correlates to the base fluid. Table 2 depicts the assessment of the friction factor and heat transfer for the varied values Σ a with the available reported work (Ishak et al. [54]) for the limiting cases when m = 0, ϕ = 0, u e (x) = 0, We a = 0, N r = 0, Pr = 1.0, γ = 1.0, and θ(0) = 1. The results show that the upper branch's current fallouts are in good arrangement with the published research of Ishak et al. [54]. Moreover, Table 3 is also constructed to compare the heat transfer numerical upper branch solution values with the published work of Grubka and Bobba [55], Ali [56], and Yih [57] for the several values of Pr when m = 0, ϕ = 0, u e (x) = 0, We a = 0, N r = 0, Pr = 1.0, γ = 0.0, G (0) = 1 and Bi N → ∞ . We can be confident that the supplied code is adequately developed to identify the unknown results of the given problem because the solutions are in exceptional alignment in both tables. To observe the impressions of different prominent constraints on the shear stress in the axial and transverse directions, the heat transfer and the streamline profiles of the alumina (Al 2 O 3 /water) nanoparticles for the unstable and stable branches of outcome are graphically presented in Figures 2-14, and also in Tables 4-6, respectively. The given problem simulations were performed by varying one parameter at a time and fixing the rest of the comprised parameters. The default values of the fixed embraced parameters are the following: ϕ = 0.025, We a = 0.05, γ a = −3.5, Σ a = 0.15, m = 0.5, N r = 2.0 and Bi N = 1.3. In the whole study, the stable and unstable solutions are signified by the solid blue and dashed lines, respectively. Meanwhile, the small, solid black, red, pink, and blue balls symbolize the critical or bifurcation point, where the two solutions converge to one location.      stresses and heat transfer rate for both branch solutions. These solutions describe that the shear stresses in both directions augment, along with the heat transfer, in both branch solutions. Furthermore, the critical values indicate that the Willaimson parameter delays the separation of boundary layers. In addition, it can be seen from these portraits that the heat transfer values are positive, suggesting that heat is transferring from a heated plate to a cool liquid.   The thermal radiation r N and convective parameter N Bi influences on the temperature distribution for both solution branches are illustrated in Figures 8 and 9, re-  The thermal radiation r N and convective parameter N Bi influences on the temperature distribution for both solution branches are illustrated in Figures 8 and 9, respectively. Figure 8 discloses that the temperature distribution uplifts with growing   In addition, the temperature profiles enlarged for both solution branches, with larger impacts of the solid nanoparticles' volume fraction. Physically, it is worth mentioning that, due to the increasing impacts of the nanoparticles' volume fraction, the thermal conductivity rises as a result, which produces a high improvement in the thermal boundary layer, as well as the temperature profiles.  In addition, the temperature profiles enlarged for both solution branches, with larger impacts of the solid nanoparticles' volume fraction. Physically, it is worth mentioning that, due to the increasing impacts of the nanoparticles' volume fraction, the thermal conductivity rises as a result, which produces a high improvement in the thermal boundary layer, as well as the temperature profiles.      Figures 13a,b and 14a,b present the patterns of streamlines with and without the Hall parameter. It is clear from these graphs that the first branch solution's patterns seem to be simple and straight, while the second branch solution's patterns are complicated. Tables 4 and 5 display the quantitative outcome of the shear stress in both the axial and transverse axes directions for the stable and unstable solutions, with variations in the solid nanoparticle volume fractions ϕ, the magnetic parameter Σ a , the Weissenberg parameter We a , and the Hall parameter m. Since the output of the products signifies that the shear stress in both directions (axial and transverse) escalates for the branch of stable outcomes, with sophisticated values of ϕ and Σ a , the impact of the same constraints can reverse the alumina nanomaterials flow for the unstable branch solutions. However, the shear stress augments and declines in the corresponding transverse and axial directions with higher influences of the Hall parameter. Physically, the useful conductivity decreases due to the Hall parameter, which lessens the magnetic damping force on transverse and axial velocities. Consequently, the shear stress increases in the transverse direction. Conversely, the heat transfer values for the stable and unstable solution branches are quantitatively bounded in Table 6. Here, the heat transfer behaves magnificently for the stable branch solution, compared quantitatively to the unstable solution branch, with superior impacts of ϕ and N r . The Nusselt number is quantitatively augmented in both distinct solutions as we increase the implementation of the convective parameter. In addition, the heat transfer coefficient is directly proportional to the dimensionless convective parameter. According to basic physics, if we increase the values of the convective parameter, as a result, the heat transfer also increases; therefore, the Nusselt number is significantly boosted for both branches of the solution.  Figures 13a,b and 14a,b present the patterns of streamlines with and without the Hall parameter. It is clear from these graphs that the first branch solution's patterns seem to be simple and straight, while the second branch solution's patterns are complicated.

Conclusions
This enquiry has looked at the cumulative effects of MHD and Hall current on the flow of a Williamson fluid incorporated by nanofluid through a vertical flat plate that has been convectively heated while being exposed to thermal radiation. The converted system of ODEs has been numerically solved using the bvp4c solver for various values of the pertinent parameters. From this study, the following findings can be made:

•
For the stable branch solutions, the magnetic parameter causes an increase in the axial shear stress and the rate of heat transfer, while it has an entirely different effect on the branch of unstable solutions. • Due to the higher values of the magnetic parameter, the transverse shear stress increases for upper and lower branches.

•
Due to the bigger values of the magnetic parameter, the magnitude of the bifurcation values was increased.

•
The radiation parameter is more significantly impacted by the escalating temperature distribution in both solution branches.

•
Due to the greater Hall parameter, the axial shear stress falls, while the transverse shear stress increases.

•
The heat transfer rate is increased by the nanoparticles for stable branch solutions, while it is decreased for unstable branch solutions.

•
The Williamson constraint augments the shear stresses, as well as the heat transfer. For the branch of unstable solutions, however, the behavior of the results is the same in the directions of transverse shear stress and inverted in the axial shear stress. Physically, there is a drag-type force known as a Lorentz force that tends to impede the velocity of the flow, and, consequently, the thicknesses of the velocity boundary layer become depressed. Consequently, the velocity profiles behave inversely to the coefficient of shear stresses. Hence, the shear stress escalates due to the larger consequences of Σ a . As illustrated in Figure 4, the heat transmission decreases with Σ a for the branch of stable solutions, but it increases with Σ a for the branch of unstable solutions. Since there is a large resistance to the motion of fluid due to the Lorentz force, which in turn produces the heat, as an artefact, the thickness or width of the thermal boundary layer, as well as the Nusselt number profiles, increases. Further, the temperature of fluid disappears from the flat plate at some large distance. This is not surprising, since there is a magnetic field that exerts retarding force on the free convection flow. Several physical properties of such fluids can be adjusted by varying the magnetic field. The obtained solutions undoubtedly demonstrate that the flow and the characteristics of heat transfer can be controlled using the magnetic field. Moreover, the changing critical values are obtained for each selected default values of Σ a (see Figures 2-4). With growing values of Σ a , the critical or bifurcation value elevates magnitude-wise. As a result, this pattern suggests that the flow separation decreases as the magnetic parameter's effects increase. Moreover, for negative values of γ a , there exists a critical value γ a C1, with dual solutions for γ a > γ a C1: a saddle node bifurcation at γ a = γ a C1, and no outcome for γ a < γ a C1. The separation of the boundary layer occurs at γ a = γ a C1; therefore, it is not possible to obtain the solution beyond this amount. Table 4. Numerical values of the skin friction coefficient along the x− axis direction for the several values of ϕ, Σ a , We a and m when N r = 2.0, Bi N = 1.3, and γ a = −3.5.  Table 5. Numerical values of the skin friction coefficient along the z− axis direction for the several values of ϕ, Σ a , We a and m when N r = 2.0, Bi N = 1.3, and γ a = −3.5. Figures 5-7 demonstrate the impact of the Williamson parameter on the shear stresses and heat transfer rate for both branch solutions. These solutions describe that the shear stresses in both directions augment, along with the heat transfer, in both branch solutions. Furthermore, the critical values indicate that the Willaimson parameter delays the separation of boundary layers. In addition, it can be seen from these portraits that the heat transfer values are positive, suggesting that heat is transferring from a heated plate to a cool liquid.

Parameters
The thermal radiation N r and convective parameter Bi N influences on the temperature distribution for both solution branches are illustrated in Figures 8 and 9, respectively. Figure 8 discloses that the temperature distribution uplifts with growing values of N r in both solutions. This is because higher impressions of radiation indicate the supremacy of conduction and improve the thermal conductivity of the alumina nanoparticles, resulting in a rise in the thickness of the TBL, as well as an increase in the profile of temperature. Figure 9 discloses that the augmentation in the amount of Bi N comes with the stronger convective-type heating at the wall surface of the plate, which increases the temperature slope at the flat plate for the upper branch, and a slight transformation behavior is seen for the branch of the unstable solution. The thermal efficiency can now dive deeper into the static fluid. As a result, as Bi N increases, so does the temperature and thickness of the TBL. It is worth noting that, by using reasonably large influences of the convective restriction, the uniform wall temperature θ(0) = 1 can be recovered. Furthermore, Bi N = 0 pertains to the insulated sheet case.
Figures 10-12 display the impression of the solid nanoparticle volume fraction ϕ on the axial and transverse velocity profiles and the dimensionless temperature distribution profiles for the stable and unstable solution branches, respectively. Outcomes indicate that the axial velocity portrays an intensifying tendency for the upper branch solutions initially and then starts to decline continuously for both branches due to the augmentation in the values of the solid nanoparticle volume fraction. On the other hand, the transverse velocity profile curves are amplified for the branch of upper solutions due to the larger impacts of the nanoparticle volume fraction, while it is reduced for the branch of lower solutions owing to the superior impact of nanoparticles. In addition, the temperature profiles enlarged for both solution branches, with larger impacts of the solid nanoparticles' volume fraction. Physically, it is worth mentioning that, due to the increasing impacts of the nanoparticles' volume fraction, the thermal conductivity rises as a result, which produces a high improvement in the thermal boundary layer, as well as the temperature profiles. Figure 13a,b and Figure 14a,b present the patterns of streamlines with and without the Hall parameter. It is clear from these graphs that the first branch solution's patterns seem to be simple and straight, while the second branch solution's patterns are complicated.

Conclusions
This enquiry has looked at the cumulative effects of MHD and Hall current on the flow of a Williamson fluid incorporated by nanofluid through a vertical flat plate that has been convectively heated while being exposed to thermal radiation. The converted system of ODEs has been numerically solved using the bvp4c solver for various values of the pertinent parameters. From this study, the following findings can be made:

•
For the stable branch solutions, the magnetic parameter causes an increase in the axial shear stress and the rate of heat transfer, while it has an entirely different effect on the branch of unstable solutions. • Due to the higher values of the magnetic parameter, the transverse shear stress increases for upper and lower branches. • Due to the bigger values of the magnetic parameter, the magnitude of the bifurcation values was increased.

•
The radiation parameter is more significantly impacted by the escalating temperature distribution in both solution branches.

•
Due to the greater Hall parameter, the axial shear stress falls, while the transverse shear stress increases.

•
The heat transfer rate is increased by the nanoparticles for stable branch solutions, while it is decreased for unstable branch solutions. Data Availability Statement: Did not report any data.