Effect of Thermal Radiation on Three-Dimensional Magnetized Rotating Flow of a Hybrid Nanofluid

The effect of thermal radiation on the three-dimensional magnetized rotating flow of a hybrid nanofluid has been numerically investigated. Enhancing heat transmission is a contemporary engineering challenge in a range of sectors, including heat exchangers, electronics, chemical and biological reactors, and medical detectors. The main goal of the current study is to investigate the effect of magnetic parameter, solid volume fraction of copper, Eckert number, and radiation parameter on velocity and temperature distributions, and the consequence of solid volume fraction on declined skin friction and heat transfer against suction and a stretching/shrinking surface. A hybrid nanofluid is a contemporary type of nanofluid that is used to increase heat transfer performance. A linear similarity variable is–applied to convert the governing partial differential equations (PDEs) into corresponding ordinary differential equations (ODEs). Using the three-stage Labatto III-A method included in the MATLAB software’s bvp4c solver, the ODE system is solved numerically. In certain ranges of involved parameters, two solutions are received. The temperature profile θη upsurges in both solutions with growing values of EC and Rd. Moreover, the conclusion is that solution duality exists when the suction parameter S≥Sci, while no flow of fluid is possible when S


Introduction
Fluid dynamics study has attracted the interest of experts, scholars, and researchers from numerous fields in current times owing to its multiple applications in engineering, science, and technology, as well as biopolymer-based detectors for medical sensing applications and diagnosis. Sakiadis [1] initially suggested the notion of a boundary layer steady flow on a stretching two-dimensional surface. Later, Crane [2] modernized the thought of Sakiadis. Then he applied it on steady-flow two-dimensional linear exponentially stretching surfaces. He claimed that the speed with which a sheet is stretched from a slit is related to the distance between them. Due to the fact that nanoparticles can pass through tissues and cells, chemical researchers state that nanotechnology can offer innovative solutions for treating stenosis, which is harmful and can cause death. It is noteworthy to mention the advanced evolvement of nanoparticles in drugs [3][4][5][6][7]. Choi [8] initiated research with a focus on nanoparticles by revealing their dynamic and abnormal attributes. Furthermore, suction and the Newtonian heating effect in a hybrid nanofluid. Furthermore, Devi and Devi [28], with Newtonian heating, expanded work to three-dimensional time-independent. They reported that a hybrid nanofluid transfers heat more quickly than a normal nanofluid. After that, Devi and Devi [29] examined the two-dimensional steady flow of a hybrid nanofluid over a stretched sheet. In a hybrid nanofluid, Aly and Pop [30] studied the two-dimensional steady flow over a stretched and shrunk plate using suction, the biot number, and MHD effects. Wainiet al. [31] explored the influence of hybrid nanoparticles on a fluid's steady flow in an exponentially stretching/shrinking two-dimension allayer. It was found that suction effects and shrinking effects had dual solutions. Lund et al. [32] analyzed two-dimensional unsteady flow of a hybrid nanofluid on stretching and shrinking sheets under the influence of radiation and MHD. Teh and Asghar [33] studied the MHD hybrid nanofluid flow under the influence of Joule heating on a three-dimensional rotating stretching/shrinking surface. The act of creating heat by passing current through a conductor is known Joule heating, also known as ohmic heating, resistance heating, or resistive heating. Incandescent lightbulbs, resistance ovens, electric stoves, soldering irons, and cartridge heaters are all examples of Joule heating. Additionally, Khashi et al. [34] assessed the two-dimensional steady flow behavior of a hybrid Cu − Al 2 O 3 /water nanofluid associated with a radially stretched/shrunk sheet under the impact of Joule heating, MHD, and suction. Moreover, Yan et al. [35] explored the impact of Joule heating on electrically conducting hybrid nanofluid flow moving over an exponential surface.
Thermal radiation has a significant influence at high operating temperatures that cannot be overlooked. Radiation sensitivity is crucial in the design of a suitable method since several industrial processes take place at extremely high temperatures. It also plays a substantial role in a variety of industrial applications, such as glass processing, furnace construction, internal combustion engines, plasma physics, and spacecraft [36].Sreedevi et al. [37] investigated the two-dimensional unsteady flow of mass and heat transfer of a hybrid nanofluid on a stretched surface by radiation. Recently, Dero et al. [38] studied the numerical analysis of the Cu + Al 2 O 3 /water hybrid nanofluid in terms of cross-flow and stream wise under the influence of thermal radiation. Some study has been done on flow and thermal radiation. Some current publications can be found in these studies [39][40][41].
In the two-dimensional case, research on hybrid nanofluids under MHD, radiation, and Joule heating has been rigorous (see Waini et al. [18], Khashi'ie et al. [34], and Yan et al. [35]). According to the literature survey on previous studies, it is critical to address such a problem as three-dimensional flow since this setting enables more extensive and realistic real-world applications. The present study expanded on the work of Anuar et al. [22] by incorporating the influence of radiation, MHD, and Joule heating using the research of Tiwari and Das [42]. Therefore, for this study, a novel physical model of the influence of thermal radiation on the three-dimensional magnetized rotating flow of a hybrid nanofluid with Joule heating is developed. To reach high convective heat transfer efficiency, in the hybrid nanofluid measured in this study, alumina (Al 2 O 3 ) and copper (Cu) nanoparticles are used. By suspending two nanoparticles, a hybrid nanofluid (Al 2 O 3 − Cu/water) is designed. Profiles of temperature and velocity for several values of volume fraction copper, magnetic, Eckert number, and radiation parameters of solid nanoparticles are examined in this study. Furthermore, the effect of solid nanoparticle copper against stretched/shrunk and suction parameters in terms of declined heat transfer and change in skin friction are also included in this study. The present numerical findings are compared to the results of prior investigations for comparison purposes. To the best of our knowledge, this model is different and new and no related article has been found in the literature.

Mathematical Formulation
Three-dimensional hybrid nanofluid steady flow with heat transfer in stretched and shrunk sheets is shown in Figure 1, which has Cartesian coordinates with the x-axis and the y-axis measured in the z = 0 plane. The fluid that occupies half of the area is at z ≥ 0. In the x-direction, the surface is deemed to be stretched/shrunk only when velocity U w = ax and temperature is T w = T ∞ + T 0 x 2 [43]. The entire system rotates at a constant velocity Ω along a perpendicular to the surface z-axis. In the z-axis, there is a variable magnetic field B and radiation q r . The surface of the hybrid nanofluid is constant and has ambient temperature T ∞ .T w is the wall temperature, and T 0 is the characteristic temperature. The mathematical governing equations help understand the physical model by indicating that the flow is incompressible, Newtonian, and laminar. Figure 1 shows that fluid stages and nanoparticles are considered in a thermal equilibrium state. They are distinctly tiny and uniform in shape, allowing the slip velocity among phases to be discounted.

Mathematical Formulation
Three-dimensional hybrid nanofluid steady flow with heat transfer in stretched and shrunk sheets is shown in Figure 1, which has Cartesian coordinates with the x-axis and the y-axis measured in the = 0 plane. The fluid that occupies half of the area is at ≥ 0. In the x-direction, the surface is deemed to be stretched/shrunk only when velocity = and temperature is = + [43]. The entire system rotates at a constant velocity Ω along a perpendicular to the surface z-axis. In the z-axis, there is a variable magnetic field and radiation . The surface of the hybrid nanofluid is constant and has ambient temperature T .T is the wall temperature, and is the characteristic temperature. The mathematical governing equations help understand the physical model by indicating that the flow is incompressible, Newtonian, and laminar. Figure 1 shows that fluid stages and nanoparticles are considered in a thermal equilibrium state. They are distinctly tiny and uniform in shape, allowing the slip velocity among phases to be discounted. The governing equation for a hybrid nanofluid is as follows [22]: The boundary conditions are = , = 0 , = , = , = 0, → 0, → 0, → , → ∞ The hybrid nanofluid velocity components , , and are along the x-axis, the yaxis, and the z-axis, respectively. is the stretched/shrunk parameter such that < 0 shows the shrinking surface, = 0 represents the static surface, and > 0 indicates the The governing equation for a hybrid nanofluid is as follows [22]: The boundary conditions are The hybrid nanofluid velocity components u, v, and w are along the x-axis, the y-axis, and the z-axis, respectively. λ is the stretched/shrunk parameter such that λ < 0 shows the shrinking surface, λ = 0 represents the static surface, and λ > 0 indicates the stretching sheet. T represents temperature. Hybrid nanofluid thermophysical properties were used to determine the above equation, as mentioned in Tables 1 and 2. Alumina (Al 2 O 3 ) nanoparticles are denoted by φ 1 , and copper (Cu) nanoparticles are indicated by φ 2 . It should be noted that the current model is used for only spherical nanoparticles because it gives better heat transfer performance. Moreover, c p , ρ, µ, σ, and k correspond with specific heat capacity, density, dynamic viscosity, electrical conductivity, and thermal conductivity, respectively. Fluid, nanofluid, hybrid nanofluid, solid nanoparticles 1 (Al 2 O 3 ), and solid nanoparticles 2 (Cu) are denoted by subscripts f , n f , hn f , s1, and s2, respectively.
The similarity factors indicated below are employed in this example [44].
Here, prime shows the differentiation for η, a is a stretching constant, and ν f is the kinematic viscosity. Although w w = − √ aν f S, S stands for the injection/suction parameter.
When S < 0, the flow is injection, and when the S > 0, the flow is suction.
Subject to conditions: where ω is the rotation, M the magnetic parameters, Ec is the Eckert number, q r is the radiative heat flux, and Pr is Prandtl number, which is described by: The skin friction coefficients C f x and C f y show the x-axis and the y-axis, respectively. The local Nusselt number Nu x is specified as: By using Equations (6) and (12), we get: Rd θ (0) (13) Here, Re = U w x ν f is the local Reynolds number. Table 1. Thermophysical properties of a hybrid nanofluid [33].

Names Properties
Dynamic viscosity Density Thermal conductivity Electrical conductivity The physical properties of the solid nanoparticles and the base fluid(water) are presented in Table 2 [27,33].

Stability Analysis
The system of Equations (7)- (9) indicates that more than one numerical solution exists in a specific array of stretched/shrunk and suction parameters. Therefore, an analysis of stability is needed to determine which of the solutions is stable [45,46]. That being stated, the transformation requires the introduction of new similarity variables, which are defined as where τ represents the dimensionless time (t). When Equation (14) is substituted into the unsteady Equations (2)-(4), the following equations can be acquired.
subject to boundary conditions The boundary conditions stated in Equation (23) must be changed prior to resolving the problem of eigenvalue. Harris et al. [47] proposed relaxing and replacing one of the far-field conditions. So, in this research, relaxing F 0 (∞) → 0 is impossible since it will be substituted by F 0 (0) = 1, which can influence the other numerical solutions and boundary conditions. As a result, we selected F 0 (∞) → 0 to be loosened and substituted with F 0 (0) = 1. This allows the stability of the solution to be determined, where ε 1 > 0 indicates that the solution is stable, whereas ε 1 < 0 shows that the solution is unstable owing to the development of the disturbance of the solution.

Results and Discussion
The solution's duality has been achieved in the figures by using various initial guesses for f (0), g (0), and θ (−0), with the outcome that both velocities and temperature profiles satisfied the boundary condition η → ∞ asymptotically. By using the software MAT-LAB bv4pc solver, we numerically solved nonlinear ODEs (Equations (7)-(9)) with the boundary equation (Equation (10)). The bvp4c solver was created by Jacek Kierzenka and Lawrence F. Shampine of Southern Methodist University, Texas [48]. The bvp4c solver is a finite difference method that uses the three-stage Lobatto IIIA algorithm to provide numerical solutions with fourth-order precision. Before starting to discuss the results of the current study, we have compared the coding of a numerical method to make sure that our computer code is working properly. First, to validate the coding of a numerical scheme in this study, the reduced skin friction f (0) and g (0) for pure waterwhen φ 1 = φ 2 = 0, S = 0, M = Rd = Ec = 0, and Pr = 6.2 is compared with values issued by Nazar et al. [49] and Anuaret al. [22], as in Table 3, for stretching sheet λ = 1. The findings are consistent with those of previous studies. To yield the wanted hybrid nanofluid, i.e., Al 2 O 3 − Cu/water, first, we show that alumina is isolated in the base fluid (water). After that, copper (Cu) is dispersed into the nanofluid Al 2 O 3 /water. Here, φ 1 represents alumina (Al 2 O 3 ) and φ 2 represents copper (Cu). Table 3. The comparison results of rotation parameter ω at different values when φ 1 = φ 2 = 0, λ = 1, S = 0, M = Rd = Ec = 0, and Pr = 6.2.

Nazar et al. [49]
Anuar et al. [22] Current Results The suction parameter ( S) of various values with the variation in f (0), g (0), and -θ (0) with three values of solid nanoparticle volume fraction copper φ 2 ispresented in  Figure 2b, the first solution declines and the second solution of g (0) is improved when the values of φ 2 increase. As can be seen in Figure 2c, when the value of φ 2 rises, the magnitude of −θ (0) is also decreased in both first and second solutions.
Due to the existence of multiple solutions, we investigated which of the governing parameters can play a role in the formation of dual solutions. Against varying values of solid nanoparticles, Figure 3a-c demonstrates volume fraction φ 2 = 0.0, 0.03, and 0.06 along with reduced skin friction f (0) and g (0) and reduced heat transfer rate −θ (0) through the presence of parameters φ 1 = 0.01, Pr = 6.2, M = 0.01, Rd = 0.1, Ec = 0.01, Ω = 0.01, and S = 2.2. For shrinking sheet λ ≤ −1, we can understand that the presence of nonunique solutions is notable. When i = 1, 2, 3,λ ci represents the critical point at which both solutions meet at the shrinking sheet parameter. It can be shown that no solution exists when λ < λ ci but a unique solution is observed when λ > −1. Beyond these critical values, the boundary layer estimate is no longer defensible. Furthermore, the values of λ c1 = −1.2019, λ c2 = −1.4011, and λ c3 = −1.5601 are the corresponding critical points of φ 2 = 0.0, 0.03, and 0.06, respectively. This is proof that increasing the value of φ 2 can cause the boundary layer separation to be delayed. Figure 3a shows that in the first solution, the value of f (0) increases when λ ≤ 0 and decreases when λ > 0, while the second solution declines after the increment of φ 2 . In Figure 3b, when the influence of φ 2 improves, the value of g (0) decreases in the first and increases in the second solution. In Figure 3c, −θ (0) declines in both solutions when φ 2 increases. Similar findings can be seen in Anuar et al. [22,26]. Figure 4a-c depicts different values of M= 0, 1, and 3 when the velocity and temperature profiles are ( f (η), g(η)) and θ(η), respectively, under different parameters, such as φ 1 = 0.01, φ 2 = 0.03, Pr = 6.2, S = 1.9, Rd = Ec = 0.1, Ω = 0.01, and λ = −1. Figure 4a illustrates that f (η) profiles decrease with an upsurged value of M in both (first and second) solutions due to the high drag force created by the magnetic number. The effect of M on g(η) profiles is shown in Figure 4b. It is observed that g(η) profiles reduce with an upsurged value of M in both solutions. Figure 4c shows that θ(η) increases as the value of M increases in both solutions. As a consequence, the influence of heat transfer is increased as M increases. The Lorentz force, which is caused by the magnetic field, makes the rate of transfer less resistant [49]. can cause the boundary layer separation to be delayed. Figure 3a shows that in the first solution, the value of ′′(0) increases when λ ≤ 0 and decreases when λ > 0, while the second solution declines after the increment of . In Figure 3b, when the influence of improves, the value of (0) decreases in the first and increases in the second solution. In Figure 3c, − ′(0) declines in both solutions when increases. Similar findings can be seen in Anuar et al. [22,26].   Figure 4b. It is observed that ( ) profiles reduce with an upsurged value of in both solutions. Figure 4c shows that ( ) in-  only exists in Equation (9). As a result, the values of have no influence on the velocity profile since they are uncoupled from the momentum Equations (7) and (8).The thermal boundary layer thickness constantly rises in the first and second solutions with a growing value of , which means that higher values result in a reduced temperature gradient at the surface.Due to the existence of high radiation, a big quantity of heat energy is produced in the system, which causes the temperature to increase and implies that the temperature of the fluid ( ) rises.  The plots of the velocity profile ( f (η), g(η)) and the temperature profile θ(η) against solid nanoparticle volume fractions φ 2 = 0.02, 0.04, and 0.06 of copper are portrayed in Figure 5a-c, respectively. Figure 5a shows that in thefirst and second solutions, the fluid velocity f (η) decreases when the value of φ 2 is enhanced. From Figure 5b, the g(η) profile in both solutions declines with a growing value of φ 2 . Figure 5c show that in the first and second solutions, the temperature profile θ(η) increases when the value of φ 2 is enhanced. Figure 6 illustrates the value of Ec at a different point for the temperature profile θ (η). In the first and second solutions, the values of Ec increase, causing θ (η) to rise as well. The intensity of heat transfer rises as the value of Ec grows due to the increasing heat created by Joule heating. It is worth noting that the influence of the increasing values of Ec on ( f (η), g(η)) is not important because, as seen in Table 4, the values remain the same for the increasing values of Ec. Some variations in the value of Ec are provided in Table 4.  Figure 7 shows the various values of radiation parameter Rd = 0.0, 0.2, and 0.4 for the temperature profile θ (η). The radiation parameter Rd only exists in Equation (9). As a result, the values of Rd have no influence on the velocity profile since they are uncoupled from the momentum Equations (7) and (8).The thermal boundary layer thickness constantly rises in the first and second solutions with a growing value of Rd, which means that higher Rd values result in a reduced temperature gradient at the surface.Due to the existence of high radiation, a big quantity of heat energy is produced in the system, which causes the temperature to increase and implies that the temperature of the fluid θ (η) rises.      Because dual solutions are obvious in the data, we ran a stability analysis to assess the solution's practicability. Based on our findings, we determined that only the first option is stable and viable for use as guidance in real-world applications. This claim is supported by tabulation (see Table 5), which shows that the solution is stable when the executed minimum eigenvalues are positive and non-stable otherwise.  Because dual solutions are obvious in the data, we ran a stability analysis to assess the solution's practicability. Based on our findings, we determined that only the first option is stable and viable for use as guidance in real-world applications. This claim is supported by tabulation (see Table 5), which shows that the solution is stable when the executed minimum eigenvalues are positive and non-stable otherwise.

Conclusions
In this study, a three-dimensional Al 2 O 3 − Cu/water MHD hybrid nanofluid with heat transfer rotating flow on a linear stretched/shrunk surface underthe effect of radiation and Joule heating has been investigated through a bvp4c solver on the MATLAB computing platform. This research concentrates on the behavior of f (0), g (0), −θ (0), f (η), g(η), and θ(η) underthe effect of suction/injection, MHD, radiation, and Joule heating on