MHD Mixed Convection Flow of Hybrid Ferroﬂuid through Stagnation-Point over the Nonlinearly Moving Surface with Convective Boundary Condition, Viscous Dissipation, and Joule Heating Effects

: This paper discusses a numerical study performed in analysing the performance regarding the magnetic effect on the mixed convection stagnation-point ﬂow of hybrid ferroﬂuid, examining the inﬂuence of viscous dissipation, convective boundary condition as well as Joule heating across a nonlinearly moving surface. Additionally, the hybrid ferroﬂuid exhibits an asymmetric ﬂow pattern due to the buoyancy force affecting the ﬂow. Water ( H 2 O ) is employed as the base ﬂuid collectively with the mixtures of nanoparticles containing magnetite ( Fe 3 O 4 ) and cobalt ferrite ( CoFe 2 O 4 ) , forming a hybrid ferroﬂuid. The partial differential equation’s complexity is reduced by similarity transformation into a system of ordinary differential equations, which are then numerically solved by applying the MATLAB function bvp4c for a speciﬁc range of values regarding the governing parameters. Dual solutions were identiﬁed under both opposing and assisting ﬂow conditions, and the stability analysis identiﬁed that the ﬁrst solution was stable. Furthermore, it was also revealed that the addition of 1% CoFe 2 O 4 in hybrid ferroﬂuid led to a higher skin friction coefﬁcient between 3.35% and 7.18% for both assisting and opposing ﬂow regions. Additionally, the growth of magnetic ﬁelds results in a reduced heat transfer rate between 8.75% to 10.65%, whilst the presence of the suction parameter expands the range of solutions, which then delays the boundary layer separation. With the Eckert number included, the heat transfer rate continuously declined between 7.27% to 10.24%. However, it increased by about 280.64% until 280.98% as the Biot number increased.


Introduction
Convective transport has played a vital role in various real-world applications, for instance, energy-distribution networks, energy-generating plants, as well as environmental concerns.Convective transport notions and theories have been incorporated in several fluid models that contain both Newtonian and non-Newtonian fluid types owing to the breadth and depth of their theory as well as concepts.Here, the findings of the optimal thermal fluid result in further advancements in fluid flow modelling.Due to its economic significance, nanotechnology is now more prevalent across a broad range of industries.The application of nanotechnology, which involves the combination of base fluid and nanoparticles, has brought about a modification to heat transfer fluids.Beyond that, Choi and Eastman [1] used the term "nanofluid" to describe a brand-new kind of fluid.Magnetic nanoparticles are disseminated in a non-magnetic base fluid in a magnetic nanofluid, commonly known as ferrofluid [2].In the nanotechnology sector, ferrofluid is employed in various applications, for example, coolants in thermal management systems, heat exchangers, and biomedical fields, for the treatment and diagnosis concerning diseases, such as cancer, rheumatoid arthritis, as well as angiocardiopathy [3].Ferrofluid research in areas with severe heat transfer has thus received a lot of interest [4,5].It may also encourage cross-disciplinary research in the fields of physics, chemistry, biology, as well as material science.Due to these concerns, various researchers and academicians have performed numerical and experimental studies on nano/ferrofluid behaviour concerning heat transfer issues [6][7][8][9][10][11].
The qualities of nano/ferrofluid are further enhanced by introducing a hybrid nanofluid, which combines two different nanomaterials to enhance thermal and rheological characteristics.Here, the occurrence of two or more nanoparticle types in a hybrid nanofluid improves its heat-transfer efficiency, thermophysical properties, as well as stability in comparison to those single-nanoparticle nanofluids [12].The study of hybrid nanofluid has become more important in the heat transmission field considering renewable energy, electromechanical engineering, as well as automotive cooling [13].This results from the wide variety of applications that can be found for the technology.In contrast, hybrid nanofluid transmits heat more quickly than other materials.Multiple review articles [14][15][16][17] have provided in-depth descriptions of the hybrid nanofluid.Hybrid nanoparticles were deployed in the earlier investigations by Turcu et al. [18] and Jana et al. [19] in their analysis.Aminuddin et al. [20] conducted a study on the fluid flow and heat transfer problem on a hybrid nanofluid containing graphene oxide (GO) and iron dioxide (Fe 2 O 4 ) dispersed in ethylene glycol (EG) as the base fluid.Their study on hybrid nanofluid includes several factors, such as magnetic field, thermal radiation, and viscous dissipation along with the presence of velocity slip.Interestingly, they found out that the heat transfer performance of hybrid nanofluid was enhanced with the influence of thermal radiation and velocity slip.Beyond that, Kho et al. [21] investigated the fluid flow and heat transfer problem of Ag − TiO 2 hybrid nanofluid on different geometry, which is on a permeable wedge displaying several effects, such as a magnetic field, thermal radiation, and viscous dissipation.In their findings, the heat transfer rate for Ag − TiO 2 hybrid nanofluid was enhanced with a higher concentration of TiO 2 .Recently, Atashafrooz et al. [22] conducted a study on Al 2 O 3 − CuO/water hybrid nanofluid for fluid flow and heat transfer problems inside an open trapezoidal enclosure with the impacts of magnetic field and thermal radiation.The fact that many researchers have evaluated numerical as well as experimental studies employing the hybrid nanofluid concept after that is worth noting.Scholars frequently employ Takabi and Salehi's [23] as well as Devi and Devi's [24,25] thermophysical correlations for hybrid nanofluids when undertaking numerical studies driven by a volume fraction-based model for hybrid nanofluids.The recent studies on thermophysical correlations by Devi and Devi [26][27][28][29] and [30,31] on thermophysical correlations by Takabi and Salehi.
Hybrid ferrofluid, which features a range of nanoparticles scattered throughout the ferrofluid, has been discovered by researchers to build a more effective fluid.Hybrid ferrofluids have caught the interest of numerous academicians in finding solutions to practical heat transmission issues because of several significant applications, which include heat dissipation, damping, dynamic sealing, etc.The utilisation of hybrid ferrofluid as a seed material to cure acid mine drainage (AMD) is an application that can be highlighted as especially significant.A long-lasting option for wastewater treatment, especially for AMD, will be provided using a hybrid ferrofluid.This is because AMD will affect animal life, aquatic life, human life, plant life, and other forms of life, as well as inflict serious environmental damage [32].Mathematical models have proliferated rapidly alongside the expansion of theoretical work.However, a review of the existing research reveals that only a small number of studies have concentrated on heat and flow transfer regarding the hybrid ferrofluid.Kumar et al. [33] proposed further study of MHD hybrid ferrofluid under the influence of non-uniform heat sources and sank on the transport of radiative thin films.They found that a hybrid ferrofluid is significantly more efficient in heat conduction than ferrofluid itself.The subsequent study of an MHD hybrid ferrofluid subjected to asymmetrical heat fall/rise was conducted by Tlili et al. [34].In other research, Manh Symmetry 2023, 15, 878 3 of 29 et al. [35] investigated how a hybrid ferrofluid submerged in porous media reacted to magnetic and radiation fields in terms of heat transmission.Recent research on this topic found that Zainodin et al. [36] examined the heat source/sink effect as well as velocity slip on exponentially deformable sheets around stagnation-point in a mixed convection flow on a hybrid ferrofluid (Fe 3 O 4 -CoFe 2 O 4 /water).In addition to the works already mentioned, references [37][38][39] contain other noteworthy publications on hybrid ferrofluid.
Due to the high level of practical interest, research has been performed on the flow properties and heat transmission concerning mixed convection.There are a wide variety of technological and industrial applications for mixed convection flow, including electronic gadgets, pipeline transportation, and nuclear reactors.The term "mixed convection" refers to a process in which both forced and free convection occurs.Convection, in general, is a method of heat transmission in which fluid moves from a hotter substance to a colder substance.It happened when both forced and natural convection mechanisms were operating at once to support flow and heat transmission.Sparrow et al. [40] conducted an earlier investigation involving mixed convection flow on the boundary layer by adding the assisting and opposing flow in their analysis.Concurrently, Jamaludin et al. [41] consider the suction effect along with heat source/sink for mixed convection flow in a nanofluid.Following then, interest in mixed convection among scholars increased, leading to some noteworthy mention [42][43][44].Since then, scholars have taken MHD into account in mixed convective flow to significantly improve boundary layer and heat transmission features.Since it depends on the heat observing variable, the MHD term is especially crucial for verifying heat transfer while producing significant items (MHD power generation, nuclear reactor, MHD generator) [45].Khan and Rasheed [46] investigated the heat transfer performance of nanofluid in MHD mixed convection flow with thermal radiation by accounting for the magnetic field impact in the mixed convection flow.They identified that the Lorentz force reduced fluid flow velocity, as well as reduced skin friction following the resistance it caused in flows through MHD.The effect of suction and thermal radiation on MHD mixed convection flow for ferrofluid over the nonlinearly moving surface was then addressed by Jamaludin et al. [8].It is appealing that the existence of MHD and suction expanded the solution range and delayed boundary layer separation.In addition, Wahid et al.'s [47] examination of the heat transfer effectiveness of a hybrid nanofluid established that the efficiency rose as the magnetic values rose.After decades of study, researchers have been fascinated by the analysis of MHD mixed convection flow, which has led to several important mentions in studies of boundary layers and heat transfer rate [27,[48][49][50].
In several fields of engineering and business, the Joule heating mechanism is being used more frequently, including the wiring of electrical and electronic devices.In theory, heat is generated because of resistive losses during the Joule heating process, which transforms electrical energy into thermal energy.The Joule heating characteristics are often described by the Eckert number and magnetic parameter in combination, while the viscous dissipation effects are usually described by the Eckert number alone.According to Reddy and Reddy [51], this control parameter is employed to raise the nanofluid temperature from the boundary layer flow perspective.Here, Khashi'ie et al. [52] examined the influence of Joule heating in a Cu-Al 2 O 3 /water hybrid nanofluid by taking into consideration the effect of MHD and suction over a moving plate.Their finding was that heat transfer performance is reduced by the Eckert number caused by Joule heating, while it has zero impact on the boundary layer separation.Concurrently, Yashkun et al. [53] assessed the MHD mixed convection hybrid nanofluid flow regarding Joule heating passing through the exponentially moving surface.Their research revealed the presence of a dual solution for the shrinking region within a particular range, and the impact of Ec widened the range of solutions.Recently, Zainal et al. [54] sought to conduct research on hybrid nanofluids toward stagnation-point flow for the unsteady flow problem while considering the influence of Joule heating together with viscous dissipation.The examination of boundary layers and heat transport in Joule heating has been a fascinating topic addressed by researchers over the years, and it has produced several noteworthy mentions [55][56][57][58].
Additionally, the attention in boundary layer flow exploration also includes the convective boundary condition.Convective boundary conditions are presumed to exist when heat convection and conduction at a material's surface both flow in the same direction.Physically speaking, it is possible to conceptualise it as a state in which the heat conduction rate at the surface equals the heat convection rate.This kind of scenario offers broad and accurate forecasting, especially in a variety of engineering and industrial tasks like the conjugate heat transfer around fins, material drying, heat exchangers, the transpiration cooling process, and so forth.Similarly, high-temperature processes are required for convective heat transfers.Such circumstances occur in nuclear power plants, material drying, laser therapy, haemodialysis, oxygenation, sanitary fluid transmission, and gas turbines [59,60].The effect of convective boundary conditions on MHD flow regarding hybrid nanofluid over the permeable moving surface is discussed by Aly and Pop [61] in their paper.They found that the temperature profile changed because of the rise in the Biot number, which caused an increase in the thickness of the thermal boundary layer in both instances of stretching and shrinking.Jusoh et al. [62] analyse the heat transfer performance for Ag-Cu hybrid nanofluid flow over the stretching/shrinking surface with the influence of convective boundary conditions along with viscous dissipation.They determined that as the Biot number rose, the temperature profile steadily decreased because the thermal reversal was caused by viscous dissipation.Recently, Khashi'ie et al. [63] evaluated the heat transfer performance regarding the MHD stagnation-point of hybrid nanofluid flow on a shrinking disk, including the influence of convective boundary conditions, Joule heating as well as viscous dissipation.Notably, a few intriguing articles [64][65][66][67] investigated convective boundary conditions.
As far as the authors are concerned, based on the literature reviewed and experience, no prior study has examined the boundary layer flow as well as heat transfer analysis, including the identification of the dual solutions and the stability of the system on the effects of Joule heating, viscous dissipation, along with the convective boundary condition by integrating mixed convection hybrid ferrofluid stagnation-point flow across a nonlinearly permeable moving surface.Moreover, a review of the existing research reveals that a limited number of studies have concentrated on heat and flow transfer regarding the similarity variable proposed by Shen et al. [68].Hence, the current work is motivated by the aforementioned literature and aims to fill the void by broadening the work of Jamaludin et al. [8] as well as Shen et al. [68] in developing new mathematical models following the Tiwari and Das's model [69].This is performed via the nonlinearly permeable moving surface by including new control parameters, such as viscous dissipation, Joule heating, including convective boundary conditions.Beyond that, this current work also contributed with a distinct type of fluid that incorporates an advanced type of nanofluid known as Fe 3 O 4 -CoFe 2 O 4 /water hybrid ferrofluid into the MHD mixed convection stagnation-point flow compared to the previous studies.Noteworthy to mention, the novel aspect of this work is the unforeseen complexity of the flow and heat transfer behaviour of Fe 3 O 4 -CoFe 2 O 4 /water hybrid ferrofluid comes from the interaction of the relevant regulating parameters, such as convective boundary conditions, Joule heating, viscous dissipation, as well as effects, such as a magnetic field and suction across a nonlinearly permeable moving surface.Via the similarity transformations, one can acquire the fluctuation in a mathematical model that may subsequently be solved numerically via the use of the MATLAB program bvp4c.An additional interesting aspect that will be explored is the possibility of dual solutions (first as well as second solutions).As proposed by Merkin [70] as well as Weidman and Turner [71], the dual solutions will go through the stability analysis procedure.The stability between the two different solutions is an important consideration that must be considered.It provides a technique for determining which solution is stable, as other scholars [27,42,54] have already done in their outstanding study on stability analysis.

Mathematical Formulations
This paper covered the steady, two-dimensional mixed convection flow of a viscous and incompressible electrically conducting fluid and heat transfer through a penetrable nonlinearly stretching/shrinking sheet.The coordinate frame was described as illustrated in Figure 1, where the stretched x-axis by the surface of the sheet and the measured orthogonal of the y-axis to the surface.Additionally, the sheet was positioned vertically, and the flow was along the x-axis.Moreover, the surface velocity u w (x) = bx n caused a stretching (b > 0) or a shrinking (b < 0) of the wall surface in the x-direction.The flow was affected by Joule heating and a transverse magnetic field of intensity B y (x) = B 0 x (n−1)/2 that was employed along the positive y-axis orthogonal to the surface being shrunk or stretched.It was hypothesised that the produced magnetic field was inconsequential regarding the results of the current research.Additionally, the flow near the stagnation points free-stream velocity is expressed by the symbol u e (x), whereas the velocity of the wall's mass flux is defined by v w (x) = − n+1 2 √ av f x (n−1)/2 S, in which v w < 0 in the case of injection as well as v w > 0 in the case of suction.
This paper covered the steady, two-dimensional mixed convection flow of a viscous and incompressible electrically conducting fluid and heat transfer through a penetrable nonlinearly stretching/shrinking sheet.The coordinate frame was described as illustrated in Figure 1, where the stretched -axis by the surface of the sheet and the measured orthogonal of the -axis to the surface.Additionally, the sheet was positioned vertically, and the flow was along the -axis.Moreover, the surface velocity   () =   caused a stretching ( > 0) or a shrinking ( < 0) of the wall surface in the -direction.The flow was affected by Joule heating and a transverse magnetic field of intensity   () =  0  (−1)/2 that was employed along the positive -axis orthogonal to the surface being shrunk or stretched.It was hypothesised that the produced magnetic field was inconsequential regarding the results of the current research.Additionally, the flow near the stagnation points free-stream velocity is expressed by the symbol   (), whereas the velocity of the wall's mass flux is defined by   () = − +1 2 √    (−1) 2 ⁄ , in which   < 0 in the case of injection as well as   > 0 in the case of suction.Next, convection from the hot fluid is used to maintain the sheet's temperature   () =  ∞ +  0  2−1 , where ℎ  = ℎ 0  (−1)/2 (see Fatunmbi et al. [66]) is the coefficient of heat transfer.The ambient temperature at this point was  ∞ .The critical assumption of this study is that the stagnation point happens when the plate temperature equals the ambient air temperature, denoted by   () =  ∞ .In the case of a mixed convection fluid flow, buoyancy force can arise due to density differences in the fluid.The temperature gradient difference in the fluid causes density variations that lead to buoyancy forces.As a result, the fluid can flow asymmetrically due to the buoyancy force.Next, convection from the hot fluid is used to maintain the sheet's temperature T w (x) = T ∞ + T 0 x 2n−1 , where h f = h 0 x (n−1)/2 (see Fatunmbi et al. [66]) is the coefficient of heat transfer.The ambient temperature at this point was T ∞ .The critical assumption of this study is that the stagnation point happens when the plate temperature equals the ambient air temperature, denoted by T w (x) = T ∞ .In the case of a mixed convection fluid flow, buoyancy force can arise due to density differences in the fluid.The temperature gradient difference in the fluid causes density variations that lead to buoyancy forces.As a result, the fluid can flow asymmetrically due to the buoyancy force.
The ferroparticles in the working fluid are a mix of magnetite Although the agglomeration effects are ignored, the hybrid ferrofluid takes on a stable spherical shape that is nanoparticle-sized throughout.The nanofluid model suggested by Tiwari and Das [69] is implemented in this instance.It is essential to acknowledge the single-phase approach employed in this nanofluid model as well as the assumption that all the particles are uniform in shape and size.In addition, the interactions between the particles and the fluid with which they come into contact are not considered [72,73].Under the assumption of the boundary layer and Boussinesq approximations mentioned before, the governing equations may be expressed (Jamaludin et al. [8]; Khashi'ie et al. [52]; Tiwari and Das [69]; Lund et al. [74]) ∂u ∂x with respect to the boundary conditions (Jamaludin et al. [8]; Zainal et al. [65]) At this point, the x-axis velocity component is denoted by u, while the y-axis velocity component is expressed by v. Additionally, the acceleration resulting from gravity is defined as g, the temperature of the hybrid ferrofluid is referred to by T, T 0 denotes the hybrid nanofluid temperature characteristic with T 0 > 0 signifies assisting flow, whereas T 0 < 0 signifies opposing flow, and the mass flux velocity S is categorised into S > 0 for suction while S < 0 for injection.Moreover, a, b, and n represent a positive constant.Furthermore, , and σ hn f illustrate the heat capacitance, thermal expansion coefficient, density, dynamic viscosity, as well as the electrical conductivity of the hybrid ferrofluid (Fe 3 O 4 -CoFe 2 O 4 /water), respectively.Furthermore, the thermal diffusivity of the hybrid ferrofluid is defined as , where k hn f expresses the thermal conductivity of the hybrid ferrofluid.In addition, Table 1 lists the theoretical models that were employed to describe the thermophysical properties of nanofluid as well as hybrid nanofluid, and Table 2 compiles the thermophysical properties of the base fluid as well as suspended nanoparticles.It is worth mentioning that subscripts n f , f , as well as hn f symbolise the nanofluid, fluid, and hybrid ferrofluid, respectively.

Properties Nanofluid Hybrid Nanofluid
Viscosity Electrical Conductivity Thermal Conductivity Following Jamaludin et al. [8] and Shen et al. [68], this study now presents the relevant similarity variables, which are defined as follows: In this context, "prime" denotes differentiation conditioned on η.When the similarity variables from Equation ( 5) are substituted into the governing equations from Equations ( 1)−( 3) with boundary condition Equation ( 4), the ODEs are obtained as follows: subject to in which M describes the magnetic parameter, Pr resembles the Prandtl number, Ec expresses the Eckert number, c denotes the stretching/shrinking parameter for c > 0 (shrinking) and c < 0 (stretching), λ represents the mixed convection parameter having conditions that λ < 0 (opposing flow) or λ > 0 (assisting flow), S indicates as suction/injection parameter, where S > 0 refers to for the suction and S < 0 represents the injection, while Bi refers to the Biot number.The governing parameter stated above may be expressed as given below: in which the local Grashof number, Gr x , as well as the local Reynold number, Re x , are described as: The physical quantities regarding practical interest are the skin friction coefficient, C f with the local Nusselt number, Nu x expressed as given below: Symmetry 2023, 15, 878 8 of 29 Replacing Equation ( 5) with Equation ( 10) yields

Stability Analysis
According to the numerical results in Equations ( 6)−( 8) and their respective governing parameters, it appears that there is a dual solution.As a direct result of this, a stability analysis was conducted to examine unstable and stable solutions.Equation ( 1) was maintained before the beginning of the stability analysis, but Equations ( 2)−(4) were deemed unsteady as follows: where t represents the time.Following the similarity solutions (Equation ( 5)), a new similarity transformation will be introduced for the unsteady issue as follows: As a result, the modified Equations ( 13) and ( 14) can be expressed as follows: 1 Pr subject to Weidman and Turner's [71] work established the foundation for our investigation into the dual solutions' stability, which may be summarised as follows: in which γ refers to the unknown eigenvalue determining the stability of the solution.Concurrently, G(η) and F(η) are relatively opposed to θ 0 (η) and f 0 (η), respectively.The exponential disturbance demonstrates the disruption's rapid growth or onset.Thus, by employing the time-dependent solutions of Equation (19) in Equations ( 16)−( 18), the subsequent expressions can be obtained: Symmetry 2023, 15, 878 9 of 29 1 Pr subject to The infinite set of eigenvalues γ 1 < γ 2 < γ 3 < . . .were given by Equations ( 20)−( 22), in which γ 1 represents the smallest eigenvalue which is crucial in determining the stability with respect to the dual solution.As per Weidman and Turner [71], either one of the boundary conditions F 0 (∞) = 0 as well as G 0 (∞) = 0 will be relaxed to demonstrate the infinite set of eigenvalues γ 1 < γ 2 < γ 3 < . . . .Furthermore, the boundary condition F 0 (∞) = 0 was taken into consideration to be relaxed and then substituted with a new boundary condition F 0 (0) = 1 in solving the linearised equation stated before.

Results and Discussion
The boundary value problem that was provided, in addition to the boundary conditions in Equation ( 8), was numerically assessed using Equations ( 6) and ( 7) and the MATLAB bvp4c solver.Several governing parameters are considered in this work, including the volume percentage of nanoparticles φ 1 , φ 2 , magnetic parameter M, mixed convection parameter λ, Biot number Bi, Eckert number Ec, stretching/shrinking parameter c, as well as suction parameter S. Values φ 1 = 0.01, φ 2 = 0.01, Ec = 0.1, Bi = 0.2, M = 1, and S = 1 were implemented in MATLAB during numerical simulations and may be altered relying on a variety of governing parameters.Here, the Prandtl number was set to 6.2 (water) while the surface was shrunk for c = −0.5 and n was set to 3 (nonlinear) for all calculations in this work.The findings for nanofluid, Al 2 O 3 /water, between skin friction coefficient, Re 1/2 x C f , including Nusselt number, Re −1/2 x Nu x , were contrasted with those from Jamaludin et al. [8] as well as Nazar et al. [76] to confirm the existing model.To fulfil this validation, a necessary change was made by replacing the boundary condition in Equation (8) with θ(0) = 1 at η = 0.These were calculated in accordance with Table 3, which displayed a favourable comparison result.For the case of shrinking sheet (c < 0) with the absence of mixed convection (λ = 0), magnetic field (M = 0), suction (S = 0), Eckert number (Ec = 0) and Biot number (Bi = 0), it was found that the value of Re 1/2 x C f and Re −1/2 x Nu x are reduced as the strength of the shrinking case drop from −1.1 to −1.2 as presented in Table 3.Moreover, this current model for the hybrid ferrofluid problem is reliable in producing dual solutions as this current model can produce and validate the dual solution with Jamaludin et al. [8] as well as Nazar et al. [76].Additionally, the present results in Table 3 were extended until 6 decimal places for the accuracy of the numerical result.
The volume fraction of the nanoparticles is a critical physical metric in examining how well nanoparticles impact the rate of heat transfer and fluid flow.Tables 4 and 5 present the influence of nanoparticle volume fraction on Re 1/2 x C f and Re −1/2 x Nu x , respectively, for water (φ 1 = φ 2 = 0), Fe 3 O 4 -CoFe 2 O 4 /water hybrid ferrofluid (φ 1 = 0.01, φ 2 = 0.01) and Fe 3 O 4 -CoFe 2 O 4 /water hybrid ferrofluid (φ 1 = 0.02, φ 2 = 0.02).For both cases of assisting (λ > 0) and opposing (λ < 0) flow regions, it can be seen in Table 4 that Fe 3 O 4 -CoFe 2 O 4 /water hybrid ferrofluid (φ 1 = 0.02, φ 2 = 0.02) has the highest value of Re 1/2 x C f when compared with water (φ 1 = φ 2 = 0) and Fe 3 O 4 -CoFe 2 O 4 /water hybrid ferrofluid (φ 1 = 0.01, φ 2 = 0.01).This signature represents a colloid consisting of two solid nanoparticles suspended in a base fluid, which results in a higher value for Re 1/2  x C f .Physically speaking, the fluid viscosity was enhanced due to the increment value of the nanoparticle volume fraction.Hence, by boosting the value of nanoparticle volume fraction, it seems the shear stress also enhanced and contributes to the surge up values of Re 1/2  x C f .It is also observable that the value of Re 1/2 x C f reduced in opposing flow region (λ < 0) but enhanced in assisting flow region (λ > 0).This phenomenon happens may be due to the direction in which fluid is flowing along with the direction of the sheet being stretchable (c = −0.5)give an impact on the behaviour of Re 1/2 x C f .Beyond that, the heat transfer measured by the values of Re −1/2 x Nu x , as presented in Table 5 is observed to reduce as the nanoparticle volume fraction increase from 0 to 0.02 for both φ 1 and φ 2 .Physically, the viscosity enhanced with the increment value of nanoparticle volume fraction due to the strong interaction between the nanoparticle and the fluid molecules.This phenomenon can create a more thermal boundary layer around a heated surface and thicken the thermal boundary layer thickness.Consequently, the temperature gradient will drop, reduce the heat flux, and decrease the heat transfer rate.x C f for both ranges in assisting (λ > 0) as well as opposing (λ < 0) flow region compared to Fe 3 O 4 /water ferrofluid and water.According to the findings, the presence of one or more solid nanoparticles is likely the primary contributor to the observed disparity in Re 1/2  x C f values between the Fe  4. Additionally, in the opposing flow region (λ < 0), Figure 4 demonstrates a decrease in values of Re 1/2 x C f with the increment of M until it reached λ = λ c for M = 1, 2, 3, 4. The occurrence of a magnetic field in a flow of magnetohydrodynamic (MHD) fluid provides a countervailing force called the Lorentz force, which endures the motion of fluid particles and slows the flow.This occurs because the force acts to push the fluid away from the surface, effectively creating a resistance that reduces the interaction between the fluid and the surface.Physically, buoyancy acts in contrast to the fluid flow in an opposing flow region.This finding points to the magnetic field's role in reducing fluid flow velocity by bolstering buoyant force with resistance built by Lorentz force.Because of this, the surface has less shear stress and reduced momentum in the boundary layer.Concurrently, contrary results were obtained in the assisting flow region, in which the surge-up values of M improved the Re 1/2 x C f values.In the physical aspect, buoyancy acts in the same direction as the flow of a fluid.Hence, it will boost the fluid's flow velocity in opposition to the Lorentz force's resistance in MHD fluid flow and assist in increasing the surface shear stress.Further observations prove the rise of M slightly lowered the Re −1/2 x Nu x values, as illustrated in Figure 5. Here, it may be recognised that the occurrence of M decreased the heat transfer rate due to the increase thermal boundary layer thickness and contributed to less heat flux produced.Physically, energy is released from the frictional force between the fluid flow and the resistance form the Lorentz force.Hence, it may contribute to enhancing the thickness of the thermal boundary layer with a load of energy released from the frictional force.Additionally, the presence of M also induced electrical currents in the fluid, which can dissipate energy and reduce the temperature gradient between the fluid and the surface, which eventually reduced the value of Re −1/2 x Nu x and drops the heat transfer rate.Moreover, the critical value λ c decreases with the growth of M as visualised in Figures 4 and 5. Based on this result, it seems that incorporating a magnetic field into the boundary layer flow may be able to accelerate boundary layer separation as well as widen the range of solutions.energy released from the frictional force.Additionally, the presence of  also induced electrical currents in the fluid, which can dissipate energy and reduce the temperature gradient between the fluid and the surface, which eventually reduced the value of Re  −1/2   and drops the heat transfer rate.Moreover, the critical value   decreases with the growth of  as visualised in Figures 4 and 5. Based on this result, it seems that incorporating a magnetic field into the boundary layer flow may be able to accelerate boundary layer separation as well as widen the range of solutions.x C f is shown to decrease slightly as the value of Bi rises for the first solution in the opposing flow region (λ < 0), but the trend illustrated for the assisting flow region (λ > 0) is shown to grow slightly.This suggests the influence of Bi in opposing flow regions may slightly increase the momentum at the boundary layer and adjourn the fluid flow.In contrast, the fluid flow may slightly enhance due to slight decreases in momentum boundary layer thickness.In this instance, it is likely that the presence of Bi only possesses a tiny impact on the momentum at the boundary layer, with only modest variations occurring due to the rise in the value of Bi in the flow.Physically, heat is conducted more quickly within a body than it is dissipated to the surrounding fluid when the Biot number is low (Bi < 1).In this instance, the temperature inside the plate is almost uniform, and there is also a minimal temperature difference across the fluid border layer.As a result, the behaviour of Re 1/2 x C f is not considerably impacted by heat transfer, and the heat transfer coefficient at the fluid boundary layer is essentially constant.Nonetheless, the outcome demonstrated in Figure 7 proves the main impact of Bi on the heat transfer rate, as the Re −1/2 x Nu x values grow with the spike-up value of Bi.Physically, since it is directly admissible to the heat transfer coefficient, h f , the Biot number has an opposite relation with thermal resistance.Due to this, the surface heat transfer will rise since the heat resistance will decrease as Bi increases, as demonstrated in Figure 7.Moreover, the heat flux was enhanced since the presence of Bi reduced the thickness concerning the thermal boundary layer.Additionally, regarding the duration for boundary layer separation, both Figures 6 and 7 show that the critical value λ c rises corresponding to the increasing Bi value, which indicates that the boundary layer quickly separates as Bi increases.
Apart from that, Figures 8 and 9 delineate the Eckert number, Ec's impact on Re 1/2 x C f and Re −1/2 x Nu x , respectively, with λ.The rising value of Ec illustrates a reduced value of Re 1/2  x C f in the opposing flow region (λ < 0) but improvements in assisting the flow region (λ > 0), as illustrated in Figure 8.In essence, the Eckert number may be referred to as the possible ratio of the advective transport and the heat dissipation.Physically, the combined actions of the Joule heating, viscous dissipation process, as well as the magnetic field accelerate the Eckert number.Moreover, energy may be stored in the fluid area because of dissipation owing to viscosity and elastic deformation, leading to frictional heating as the Eckert number increases.Generally, buoyancy acts counter to fluid flow in an opposing flow region.This proposes that the heat energy produced from frictional heating is capable of supporting buoyant force, which contributes to slowing down the fluid flow velocity.Thus, this consequence reduced the momentum at the boundary layer and reduced the shear stress, as depicted in Figure 8. Concurrently, a contrary result was obtained in assisting the flow region because, in terms of the physical aspect, the buoyancy effect follows the fluid flow's direction.This recommends that the heat energy produced from frictional heating may support buoyant force, which contributes to speeding up the fluid flow velocity.Because of this effect, as seen in Figure 8, the shear stress and momentum in the boundary layer were both increased.Beyond that, the findings demonstrated in Figure 9 depict that the value of Re −1/2 x Nu x decreases as Ec increases.Physically, when an electric current source is passed through a conducting moving plate, Joule heating (which includes magnetic fields) produces heat on the moving plate.Moreover, a higher Ec value produces greater heat because of fluid-particle friction forces.This suggests that more heat was built up during the process and increased the thermal boundary layer thickness, which results in reduced convective heat transfer and leads to a lower value of Re −1/2 x Nu x .Therefore, in this specific situation, the Ec value tends to reduce the effectiveness of heat transfer as Ec grows in value.Both Figures 8 and 9 reveal that the critical value λ c expands with the increased value of Ec, which indicates that the boundary layer separates more quickly as Ec increases and reduces the range of solution.Ec value produces greater heat because of fluid-particle friction forces.This suggests that more heat was built up during the process and increased the thermal boundary layer thickness, which results in reduced convective heat transfer and leads to a lower value of Re  −1/2   .Therefore, in this specific situation, the Ec value tends to reduce the effectiveness of heat transfer as Ec grows in value.Both Figures 8 and 9 reveal that the critical value λ c expands with the increased value of Ec, which indicates that the boundary layer separates more quickly as Ec increases and reduces the range of solution.x Nu x towards dimensionless λ with various values of S. The appearance of the dual solution is recognised within the range of λ > λ c in opposing and assisting flow regions, as displayed in Figures 10 and 11, accordingly.From the observation in Figure 10, it can be seen the value of Re 1/2  x C f rose with the boost in the value of S in both opposing and assisting flow regions.Physically, the nanofluid's motion slows down and the surface's velocity gradient increases because of the suction's effect at the boundary.Thus, the heated fluid is drawn toward the wall under the action of suction, in which the buoyant forces are reduced by the strong effect in viscosity, resulting in an abrupt gradient of velocity.Hence, this phenomenon enhanced the shear stress at the wall and increased the momentum at the boundary layer.Moreover, Figure 10 illustrates the critical value λ c increased as S grows.This indicates the growth strength of suction may delay the boundary layer or flow separation.By encasing the slowly moving molecules within the shrinking sheet, the permeable sheet facilitates the maintenance of the laminar flow.Apart from that, the value of Re −1/2 x Nu x is also enhanced with the increment strength of S as portrayed in Figure 11.Physically, a higher rate of heat transmission is achieved due to suction's role in inducing fluid motion, which simultaneously aids heated particles in moving closer to the wall.More ferrofluid particles can diffuse through the sheet as the S value rises, enhancing the permeability of the sheet.Hence, this phenomenon contributes more heat flux and reduces the thickness of the thermal boundary layer.Hence, growing S accelerates the rate of heat transmission.The profiles of velocity f (η) as well as temperature θ(η) in Figures 12 and 13, accordingly, highlight the contrast between Fe 3 O 4 /water ferrofluid (φ 1 = 0.01, φ 2 = 0), Fe 3 O 4 -CoFe 2 O 4 /water hybrid ferrofluid (φ 1 = 0.01, φ 2 = 0.01), as well as water (φ 1 = φ 2 = 0) when λ = −70 (opposing flow zone).Figure 12 impressively portrayed that Fe 3 O 4 -CoFe 2 O 4 /water hybrid ferrofluid has a slightly higher f (η) in comparison to Fe 3 O 4 /water ferrofluid and water.This shows that the Fe 3 O 4 -CoFe 2 O 4 /water hybrid ferrofluid minimised the momentum boundary layer thickness and contributed to a higher fluid flow, which enhanced the shear stress at the surface compared with Fe 3 O 4 /water ferrofluid and water.This finding agrees with those in Figure 2        In the region of opposing flow (λ = −70), the behaviours of f (η) and θ(η) as a consequence of adjusting the magnetic parameter, M are demonstrated in Figures 14 and 15, accordingly.Interestingly, it is observed that the values of f (η) increase as M grows in Figure 14.This illustrates the fluid flow velocity was enhanced since the momentum boundary layer thickness was diminished with the increasing M values.Physically, the Lorentz force causes the fluid to flow in a particular direction and at a certain velocity.As the value of M grows, the Lorentz force acting on the fluid will become stronger, which will cause the fluid to flow with a greater velocity.This phenomenon raised the shear stress at the wall surface due to the speed-up of fluid flow velocity.Moreover, this result is aligned with the finding in Figure 4 as M grows at λ = −70 (opposing flow region), the value of skin friction also increases.Apart from that, a build-up of temperature in the thermal boundary layer contributes to an increase in temperature profile with an increment value of M, as illustrated in Figure 15.Physically speaking, the effect of M along with the several factors in this fluid flow problem, such as Joule heating, convective boundary condition, and viscous dissipation, boost up the amount of heat energy, which thickens more the thermal boundary layer thickness.Contributing to that, the heat flux was dropped, and the heat transfer rate was lowered, and the proof may be seen in the result from Figure 5, which depicts the value of Re −1/2 x Nu x reduced as M increased at λ = −70 (opposing flow region).Next, the profiles of velocity f (η) as well as temperature θ(η) in Figures 16 and 17, accordingly, highlight the influence of Biot number, Bi in opposing flow regions (λ = −70).As illustrated in Figure 16, the surge of Bi slightly reduced the velocity profile, which then decreased the thickness of the momentum boundary layer.This phenomenon slightly impedes the fluid flow velocity and contributes to less shear stress at the surface wall.Beyond that, this suggests the presence of Bi does not significantly affect the velocity profile, where this outcome agrees with the finding in Figure 6, indicating the increment of Bi increased the skin friction value slightly in the opposing flow region (λ = −70).In contrast, the significance of Bi may be emphasised in Figure 17, in which the surge-up Bi value greatly influences the temperature profile from η = 0 until η → ∞ .Generally, the Biot number is a measurement of the relationship between convection at the surface and conduction within a solid.It can be seen as Bi grows, the temperature profile also increases.Physically, when the Biot number rises, the surface's thermal resistance gradually decreases.The increase in convection is what leads to the greater surface temperature, including the built thermal boundary layer thickness.Moreover, the heat flux was diminished, contributing to the thickening of the thermal boundary layer, which reduced the heat transfer performance as Bi grows.Next, the profiles of velocity ′() as well as temperature () in Figures 16 and 17, accordingly, highlight the influence of Biot number, Bi in opposing flow regions ( = −70) .As illustrated in Figure 16, the surge of Bi slightly reduced the velocity profile, which then decreased the thickness of the momentum boundary layer.This phenomenon slightly the fluid flow velocity and contributes to less shear stress at the surface wall.Beyond that, this suggests the presence of Bi does not significantly affect velocity profile, where this outcome agrees with the finding in Figure 6, indicating the increment of Bi increased the skin friction value slightly in the opposing flow region ( = −70).In contrast, the significance of Bi may be emphasised in Figure 17, in which the surge-up Bi value greatly influences the temperature profile from  = 0 until  → ∞.Generally, the Biot number is a measurement of the relationship between convection at the surface and conduction within a solid.It can be seen as Bi grows, the temperature profile also increases.Physically, when the Biot number rises, the surface's thermal resistance gradually decreases.The increase in convection is what leads to the greater surface temperature, including the built thermal boundary layer thickness.Moreover, the heat flux was diminished, contributing to the thickening of the thermal boundary layer, which reduced the heat transfer performance as Bi grows.Physically, the existence of the magnetic field in Joule heating combined with viscous dissipation built up resistance from the Lorentz force in a flow.This suggests the presence of resistance from the growth of Ec reduced the fluid flow velocity, which contributes to diminishing the momentum boundary layer thickness, as portrayed in Figure 18.Therefore, this phenomenon reduced the shear stress at the wall surface.As illustrated in Figure 19, the surge-up value of Ec enhances the () in the opposing flow region ( = −70).Physically speaking, Eckert number properties are associated with the dissipation as well as Joule heating terms in the energy equation.The magnetic resistance built from the Lorentz force as well as the internal friction of the fluid, transforms the mechanical with electrical energy into heat, accordingly.Therefore, as Ec is increased, the () is also Figure 18 illustrates the behaviour of f (η) with Ec adjusted in the opposing flow region (λ = −70).It may be observed in the trend of f (η) broaden with the surge of Ec.Physically, the existence of the magnetic field in Joule heating combined with viscous dissipation built up resistance from the Lorentz force in a flow.This suggests the presence of resistance from the growth of Ec reduced the fluid flow velocity, which contributes to diminishing the momentum boundary layer thickness, as portrayed in Figure 18.Therefore, this phenomenon reduced the shear stress at the wall surface.As illustrated in Figure 19, the surge-up value of Ec enhances the θ(η) in the opposing flow region (λ = −70).Physically speaking, Eckert number properties are associated with the dissipation as well as Joule heating terms in the energy equation.The magnetic resistance built from the Lorentz force as well as the internal friction of the fluid, transforms the mechanical with electrical energy into heat, accordingly.Therefore, as Ec is increased, the θ(η) is also increased, resulting in a thicker thermal boundary layer and lowered heat flux.This phenomenon reduced the performance of heat transfer as Ec grows, and this result aligns with the outcome in Figure 9, in which the Nusselt number was reduced with the acceleration of Ec in the opposing flow region (λ = −70).Figures 20 and 21 exemplify the behaviour of f (η) as well as θ(η), accordingly, displaying a variety of values regarding the suction parameter S in the opposing flow region (λ = −70).This may be observed from Figure 20, where the value of f (η) significantly enhanced with the surge-up value of S. Physically, an increase in S leads to a higher rate of fluid extraction from the boundary, which creates a stronger pressure gradient within the boundary layer.This, in turn, causes an increase in the velocity near the boundary, and, thus, an increase in the velocity gradient or slope of the velocity profile.This phenomenon led to a decrease in momentum boundary layer thickness and improved the shear stress at the wall.Moreover, this result aligns with the finding in Figure 10, where skin friction increases as S grows in the opposing flow region (λ = −70).Additionally, the thermal boundary layer thickness as well as the temperature profile, also drop with the increased strength of S, as depicted in Figure 21.Physically speaking, a decrement in the thickness of the thermal boundary layer will assist in producing more heat flux.Because of this phenomenon, heat transfer performance was also enhanced.This finding is parallel with the outcome and observation in Figure 11, in which the Nusselt number is raised with the enhanced strength of S. enhanced with the surge-up value of .Physically, an increase in  leads to a higher rate of fluid extraction from the boundary, which creates a stronger pressure gradient within the boundary layer.This, in turn, causes an increase in the velocity near the boundary, and, thus, an increase in the velocity gradient or slope of the velocity profile.This phenomenon led to a decrease in momentum boundary layer thickness and improved the shear stress at the wall.Moreover, this result aligns with the finding in Figure 10, where skin friction increases as  grows in the opposing flow region ( = −70).Additionally, the thermal boundary layer thickness as well as the temperature profile, also drop with the increased strength of , as depicted in Figure 21.Physically speaking, a decrement in the thickness of the thermal boundary layer will assist in producing more heat flux.Because of this phenomenon, heat transfer performance was also enhanced.This finding is parallel with the outcome and observation in Figure 11, in which the Nusselt number is raised with the enhanced strength of .It is possible to see the presence of a non-unique solution, also known as a dual solution, in the region represented by Figures 2− 21 that contains both the opposing flow ( < 0) as well as the assisting flow ( > 0).In multiple cases, the dual answer will represent itself somewhere within the range of  >   .The critical value   is a unique solution that occurs at the intersection of the first as well as second solutions.Bifurcation points of the first as well as the second solutions can be said to have occurred at the point denoted by  =   In the boundary layer, separation happens when  =   , and when utilising the boundary layer approximation, one cannot achieve a solution if the range is between  <   .Therefore, an analysis of the flow's stability is essential.The same numerical technique that was addressed in the prior section was implemented to solve the linearised Equations ( 20)-( 22) numerically to execute the stability analysis.In addition, to validate the results of stability analysis, it is necessary to collect eigenvalues such as  that are relatively modest.Here, Table 6 was created using hybrid ferrofluid ( 1 = 0.01,  2 = 0.01), to illustrate how small eigenvalues  behave for different Bi = 0.2, 0.4 when  = 1,  = 1,  = 3, Ec = 0.1, and  = −0.5.The first solution provides a positive number concerning the eigenvalue  1 , which, in accordance with Merkin et al. [66] and Weidman and Turner [67], resembles the initial drop of disturbance for steady and genuine flow through time.At the same time, the second solution's negative eigenvalue  2 indicates an early increase in disturbance for an unrealistic and unstable flow as time progresses.Table 6 makes it clear that when  moves closer to the critical values represented by   , the values of the small eigenvalues expressed by  1 as well as  2 grow closer to zero.As a result, it is possible to conclude that the second solution is unstable while the first solution is stable.It is possible to see the presence of a non-unique solution, also known as a dual solution, in the region represented by Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11, Figure 12, Figure 13, Figure 14, Figure 15, Figure 16, Figure 17, Figure 18, Figure 19, Figure 20, Figure 21 that contains both the opposing flow (λ < 0) as well as the assisting flow (λ > 0).In multiple cases, the dual answer will represent itself somewhere within the range of λ > λ c .The critical value λ c is a unique solution that occurs at the intersection of the first as well as second solutions.Bifurcation points of the first as well as the second solutions can be said to have occurred at the point denoted by λ = λ c .In the boundary layer, separation happens when λ = λ c , and when utilising the boundary layer approximation, one cannot achieve a solution if the range is between λ < λ c .Therefore, an analysis of the flow's stability is essential.The same numerical technique that was addressed in the prior section was implemented to solve the linearised Equations ( 20)-( 22) numerically to execute the stability analysis.In addition, to validate the results of stability analysis, it is necessary to collect eigenvalues such as γ that are relatively modest.Here, Table 6 was created using hybrid ferrofluid (φ 1 = 0.01, φ 2 = 0.01), to illustrate how small eigenvalues γ behave for different Bi = 0.2, 0.4 when M = 1, S = 1, n = 3, Ec = 0.1, and c = −0.5.The first solution provides a positive number concerning the eigenvalue γ 1 , which, in accordance with Merkin et al. [66] and Weidman and Turner [67], resembles the initial drop of disturbance for steady and genuine flow through time.At the same time, the second solution's negative eigenvalue γ 2 indicates an early increase in disturbance for an unrealistic and unstable flow as time progresses.Table 6 makes it clear that when λ moves closer to the critical values represented by λ c , the values of the small eigenvalues expressed by γ 1 as well as γ 2 grow closer to zero.As a result, it is possible to conclude that the second solution is unstable while the first solution is stable.

Conclusions
Investigations are being performed into the steady stagnation point of an MHD mixed convection hybrid ferrofluid flow performed by the implantation of a nonlinearly permeable moving surface.The self-similar equations produced by the similarity transformation are then numerically solved by MATLAB's built-in solver (bvp4c).We examine in detail the relevance of the Eckert number, the magnetic parameter, the Biot number, the suction parameter, along with the nanoparticle volume fraction.The following are the study's primary findings:

•
Stability analysis results imply the first solution is stable.On the other hand, the second solution is unstable.

•
Varying nanoparticle volume fraction from φ 1 , φ 2 = 0 until φ 1 , φ 2 = 0.02 enhanced the skin friction by about 14.44% × 100% = 3.35% and 3.37%, respectively.Hybrid ferrofluid temperature characteristic (constant) T w (x) The temperature at the surface (K) T ∞ Ambient temperature (K) u, v Velocity component in the xand ydirections ms −1 u e Free stream velocity ms −1 u w Surface velocity for stretching/shrinking cases ms −1 v w The velocity of the wall's mass flux ms
Fe 3 O 4 and cobalt ferrite CoFe 2 O 4 , making it a hybrid ferrofluid Fe 3 O 4 -CoFe 2 O 4 with water as its base fluid.Here, the first ferroparticle volume fraction parameter, φ 1 in this case, is represented by Fe 3 O 4 , while the second ferroparticle volume fraction parameter φ 2 is represented by CoFe 2 O 4 .

Figures 2 and 3 2 x
Figures 2 and 3 depict the behaviour of Re 1/2 x C f and Re −1/2 x Nu x , respectively, against λ with the comparison between Fe 3 O 4 /water ferrofluid (φ 1 = 0.01, φ 2 = 0), Fe 3 O 4 -CoFe 2 O 4 /water hybrid ferrofluid (φ 1 = 0.01, φ 2 = 0.01), as well as water (φ 1 = φ 2 = 0).It can be observed in Figure 2 that Fe 3 O 4 -CoFe 2 O 4 /water hybrid ferrofluid has a slightly better value of Re 1/2x C f for both ranges in assisting (λ > 0) as well as opposing (λ < 0) flow region compared to Fe 3 O 4 /water ferrofluid and water.According to the findings, the presence of one or more solid nanoparticles is likely the primary contributor to the observed disparity in Re1/2  x C f values between the Fe 3 O 4 -CoFe 2 O 4 /water hybrid ferrofluid,

Figure 3 .
Figure 3. Variation of Re  −1/2   for water ( 1 = 0,  2 = 0) , Fe 3 O 4 / water ( 1 = 0.01,  2 = 0) and Fe 3 O 4 -CoFe 2 O 4 /water ( 1 = 0.01,  2 = 0.01).The magnetic parameter,  effect on Re  1/2   of hybrid ferrofluid ( Fe 3 O 4 -CoFe 2 O 4 /water) can be seen in Figure 4. Additionally, in the opposing flow region ( < 0), Figure 4 demonstrates a decrease in values of Re  1/2   with the increment of  until it reached  =   for  = 1, 2, 3, 4. The occurrence of a magnetic field in a flow of magnetohydrodynamic (MHD) fluid provides a countervailing force called the Lorentz force, which endures the motion of fluid particles and slows the flow.This occurs because the force acts to push the fluid away from the surface, effectively creating a resistance that reduces the

Figure 3 .
Figure 3. Variation of Re  −1/2   for water ( 1 = 0,  2 = 0) , Fe 3 O 4 / water ( 1 = 0.01,  2 = 0) and Fe 3 O 4 -CoFe 2 O 4 /water ( 1 = 0.01,  2 = 0.01).The magnetic parameter,  effect on Re  1/2   of hybrid ferrofluid ( Fe 3 O 4 -CoFe 2 O 4 /water) can be seen in Figure 4. Additionally, in the opposing flow region ( < 0), Figure 4 demonstrates a decrease in values of Re  1/2   with the increment of  until it reached  =   for  = 1, 2, 3, 4. The occurrence of a magnetic field in a flow of magnetohydrodynamic (MHD) fluid provides a countervailing force called the Lorentz force, which endures the motion of fluid particles and slows the flow.This occurs because the force acts to push the fluid away from the surface, effectively creating a resistance that reduces the

Figure 4 .
Figure 4. Variation of Re  1/2   against  for some  values.Figure 4. Variation of Re 1/2x C f against λ for some M values.

Figure 4 .Figure 5 .
Figure 4. Variation of Re  1/2   against  for some  values.Figure 4. Variation of Re 1/2x C f against λ for some M values.Symmetry 2023, 15, x FOR PEER REVIEW

Figures 6 and 7
Figures 6 and 7 illustrate the plots of Re  1/2   and Re  −1/2   , accordingly the mixed convection parameter,  for a broad range of values concerning the B ber, Bi.The trend of Re  1/2   is shown to decrease slightly as the value of Bi rise first solution in the opposing flow region ( < 0), but the trend illustrated for the flow region ( > 0) is shown to grow slightly.This suggests the influence of Bi i ing flow regions may slightly increase the momentum at the boundary layer and the fluid flow.In contrast, the fluid flow may slightly enhance due to slight dec momentum boundary layer thickness.In this instance, it is likely that the prese only possesses a tiny impact on the momentum at the boundary layer, with only variations occurring due to the rise in the value of Bi in the flow.Physically, he ducted more quickly within a body than it is dissipated to the surrounding fluid w

Figure 5 . 2 x
Figure 5. Variation of Re −1/2 x Nu x against λ for some M values.

Figures 6 and 7 2 x
Figures 6 and 7 illustrate the plots of Re 1/2x C f and Re −1/2x Nu x , accordingly, against the mixed convection parameter, λ for a broad range of values concerning the Biot number, Bi.The trend of Re 1/2x C f is shown to decrease slightly as the value of Bi rises for the first solution in the opposing flow region (λ < 0), but the trend illustrated for the assisting flow region (λ > 0) is shown to grow slightly.This suggests the influence of Bi in opposing flow regions may slightly increase the momentum at the boundary layer and adjourn the fluid flow.In contrast, the fluid flow may slightly enhance due to slight decreases in momentum boundary layer thickness.In this instance, it is likely that the presence

Figure 6 .
Figure 6.Variation of Re  1/2   against  for some Bi values.

Figure 7 .
Figure 7. Variation of Re  −1/2   against  for some values of Bi.Apart from that, Figures 8 and 9 delineate the Eckert number, Ec's impact on Re  1/2   and Re  −1/2   , respectively, with .The rising value of Ec illustrates a reduced value of Re  1/2   in the opposing flow region ( < 0) but improvements in assisting the flow region ( > 0), as illustrated in Figure 8.In essence, the Eckert number may be referred to as the possible ratio of the advective transport and the heat dissipation.Physically, the

Figure 7 .
Figure 7. Variation of Re  −1/2   against  for some values of Bi.Apart from that, Figures 8 and 9 delineate the Eckert number, Ec's impact on Re  1/2   and Re  −1/2   , respectively, with .The rising value of Ec illustrates a reduced value of Re  1/2   in the opposing flow region ( < 0) but improvements in assisting the flow region ( > 0), as illustrated in Figure 8.In essence, the Eckert number may be referred to as the possible ratio of the advective transport and the heat dissipation.Physically, the

Figure 7 . 2 x
Figure 7. Variation of Re −1/2 x Nu x against λ for some values of Bi.

Figure 8 .
Figure 8. Variation of Re  1/2   against  for some Ec values.Figure 8. Variation of Re 1/2x C f against λ for some Ec values.

Figure 8 .
Figure 8. Variation of Re  1/2   against  for some Ec values.Figure 8. Variation of Re 1/2x C f against λ for some Ec values.Symmetry 2023, 15, x FOR PEER REVIEW 17 of 31

Figures 10 and 11
Figures 10 and 11 highlight the distribution of Re  1/2   and Re  −1/2   towards dimensionless  with various values of .The appearance of the dual solution is recognised within the range of  >   in opposing and assisting flow regions, as displayed in Figures 10 and 11, accordingly.From the observation in Figure 10, it can be seen the value of Re  1/2   rose with the boost in the value of  in both opposing and assisting flow regions.Physically, the nanofluid's motion slows down and the surface's velocity gradient increases because of the suction's effect at the boundary.Thus, the heated fluid is drawn toward the wall under the action of suction, in which the buoyant forces are reduced by the strong effect in viscosity, resulting in an abrupt gradient of velocity.Hence, this phenomenon enhanced the shear stress at the wall and increased the momentum at the boundary layer.Moreover, Figure 10 illustrates the critical value   increased as  grows.This indicates the growth strength of suction may delay the boundary layer or flow separation.By encasing the slowly moving molecules within the shrinking sheet, the per-

Figure 9 . 2 x
Figure 9. Variation of Re −1/2 x Nu x against λ for some Ec values.

Figures 10 and 11
Figures 10 and 11 highlight the distribution of Re 1/2x C f and Re −1/2

Figure 11 .
Figure 11.Variation of Re −1/2 x Nu x against λ for some S values.

31 Fe 3 O 4 /
when Fe 3 O 4 -CoFe 2 O 4 /water hybrid ferrofluid has a higher Re 1/2 x C f than Fe 3 O 4 /water ferrofluid and water in λ = −70 (opposing flow region).Concerning the θ(η) behaviour illustrated in Figure 13, water can be shown to have a slightly reduced θ(η) than both Fe 3 O 4 /water ferrofluid and Fe 3 O 4 -CoFe 2 O 4 /water hybrid ferrofluid.This finding concludes that the thickness of the thermal boundary layer concerning water is marginally small in comparison to Fe 3 O 4 /water ferrofluid and Fe 3 O 4 -CoFe 2 O 4 /water hybrid ferrofluid, which contributes more heat flux and improves the heat transfer performance.Symmetry 2023, 15, x FOR PEER REVIEW 19 of water ferrofluid and Fe 3 O 4 -CoFe 2 O 4 /water hybrid ferrofluid.This finding concludes that the thickness of the thermal boundary layer concerning water is marginally small in comparison to Fe 3 O 4 / water ferrofluid and Fe 3 O 4 -CoFe 2 O 4 /water hybrid ferrofluid, which contributes more heat flux and improves the heat transfer performance.

Symmetry 2023 ,
15,  x FOR PEER REVIEW 20 of 3115, accordingly.Interestingly, it is observed that the values of  ′ () increase as  grows in Figure14.This illustrates the fluid flow velocity was enhanced since the momentum boundary layer thickness was diminished with the increasing  values.Physically, the Lorentz force causes the fluid to flow in a particular direction and at a certain velocity.As the value of  grows, the Lorentz force acting on the fluid will become stronger, which will cause the fluid to flow with a greater velocity.This phenomenon raised the shear stress at the wall surface due to the speed-up of fluid flow velocity.Moreover, this result is aligned with the finding in Figure4as  grows at  = −70 (opposing flow region), the value of skin friction also increases.Apart from that, a build-up of temperature in the thermal boundary layer contributes to an increase in temperature profile with an increment value of , as illustrated in Figure15.Physically speaking, the effect of  along with the several factors in this fluid flow problem, such as Joule heating, convective boundary condition, and viscous dissipation, boost up the amount of heat energy, which thickens more the thermal boundary layer thickness.Contributing to that, the heat flux was dropped, and the heat transfer rate was lowered, and the proof may be seen in the result from Figure5, which depicts the value of Re  −1/2   reduced as  increased at  = −70 (opposing flow region).

Figure 14 .
Figure 14.Variation of  ′ () against  for some  values.Figure 14.Variation of f (η) against λ for some M values.

Figure 14 .
Figure 14.Variation of  ′ () against  for some  values.Figure 14.Variation of f (η) against λ for some M values.

Figure 17 .
Figure 17.Variation of () against  for some Bi values.

Figure 18
Figure 18 illustrates the behaviour of  ′ () with Ec adjusted in the opposing flow region ( = −70).It may be observed in the trend of  ′ () broaden with the surge of Ec.Physically, the existence of the magnetic field in Joule heating combined with viscous dissipation built up resistance from the Lorentz force in a flow.This suggests the presence of resistance from the growth of Ec reduced the fluid flow velocity, which contributes to diminishing the momentum boundary layer thickness, as portrayed in Figure18.Therefore, this phenomenon reduced the shear stress at the wall surface.As illustrated in Figure19, the surge-up value of Ec enhances the () in the opposing flow region ( = −70).Physically speaking, Eckert number properties are associated with the dissipation as well as Joule heating terms in the energy equation.The magnetic resistance built from the Lorentz force as well as the internal friction of the fluid, transforms the mechanical with electrical energy into heat, accordingly.Therefore, as Ec is increased, the () is also

Symmetry 2023 ,
15, x FOR PEER REVIEW 23 of 31increased, resulting in a thicker thermal boundary layer and lowered heat flux.This phenomenon reduced the performance of heat transfer as Ec grows, and this result aligns with the outcome in Figure9, in which the Nusselt number was reduced with the acceleration of Ec in the opposing flow region ( = −70).

Figure 19 .
Figure 19.Variation of () against  for some Ec values.

Figures 20 and 21
Figures 20 and 21 exemplify the behaviour of ′() as well as (), accordingly, displaying a variety of values regarding the suction parameter  in the opposing flow region ( = −70).This may be observed from Figure 20, where the value of ′() significantly

Figure 20 .
Figure 20.Variation of f (η) against λ for some S values.

Figure 21 .
Figure 21.Variation of () against  for some  values.

Table 3 .
Comparison values of Re 1/2 x C f and Re −1/2 x Nu x for Al 2 O 3 /water nanofluid with various values of c for φ

x Nazar et al. [76] Jamaludin et al. [8] Present Result Nazar et al. [76]
The results that are denoted by parentheses ( ) pertain to the second solution.
x Nu
3 O 4 -CoFe 2 O 4 /water hybrid ferrofluid, Fe 3 O 4 /water ferrofluid and water systems.This pattern unambiguously depicts a colloid possessing two solid nanoparticles dispersed in the base fluid (Fe 3 O 4 -CoFe 2 O 4 /water hybrid ferrofluid), which exhibits greater Re 1/2x C f , values than the water and mono nanoparticle (Fe 3 O 4 /water ferrofluid).Higher values for the nanoparticle volume fraction proved to physically enhance fluid viscosity performance.As a result, the momentum boundary layer's thickness is decreased, which boosts fluid flow and increases surface shear stress.Subsequently, the graphical finding in Figure3illustrated that water enhanced heat transfer rate in comparison to the Fe 3 O 4 -CoFe 2 O 4 /water hybrid ferrofluid as well as Fe 3 O 4 /water ferrofluid for a higher value in Re −1/2 x Nu x .This indicates water can enhance the heat transfer rate given to reduced thickness in the thermal boundary layer that sparks greater heat flux compared to the Fe 3 O 4 -CoFe 2 O 4 /water hybrid ferrofluid and Fe 3 O 4 /water ferrofluid.Furthermore, this finding proposes that raising the temperature gradient can occur from the existence of other effects or factors included in the flow.Figures 2 and 3 also illustrate the critical value λ c for Fe 3 O 4 -CoFe 2 O 4 /water hybrid ferrofluid is lower than Fe 3 O 4 /water ferrofluid and water.When compared to Fe 3 O 4 /water ferrofluid and water, this suggests that hybrid ferrofluid Fe 3 O 4 -CoFe 2 O 4 /water hybrid ferrofluid can postpone the boundary layer separation as well as improve the range of solution.The magnetic parameter, M effect on Re 1/2 x C f of hybrid ferrofluid (Fe 3 O 4 -CoFe 2 O 4 /water) can be seen in Figure

•
In terms of heat transfer rate, Fe 3 O 4 -CoFe 2 O 4 /water hybrid ferrofluid has a lower Re −1/2 x Nu x compared with water and Fe 3 O 4 /water ferrofluid due to an external heat Symmetry 2023, 15, 878Malaysia would like to acknowledge the research grant from Ministry of Higher Education (grant number: FRGS/1/2020/STG06/UPNM/03/1) and the author from the Universiti Kebangsaan Malaysia (UKM) would like to express her sincere thanks to the Universiti Kebangsaan Malaysia (UKM) for funding this research work from the grant ST-2022-015.The authors declare no conflict of interest.Specific heat at constant pressure Jkg −1 K −1 Thermal conductivity of the fluid Wm −1 K −1 Ferrofluid/Nanofluid (Fe 3 O 4 /water) hn f Hybrid Ferrofluid/Hybrid Nanofluid (Fe 3 O 4 -CoFe 2 O 4 /water) s 1 Solid component for Fe 3 O 4 (Magnetite) s 2 Solid component for CoFe 2 O 4 (Cobalt Ferrite)