Dual solutions of magnetized radiative flow of Casson Nanofluid over a stretching/shrinking cylinder: Stability analysis

The enhanced thermal efficiency exhibited by Casson nanofluids offers significant practical applications across various industrial and engineering sectors. This study focuses on the mathematical investigation of the steady magnetohydrodynamic (MHD) boundary layer flow of Casson nanofluid through a stretched/shrinking cylinder, taking into account the effects of suction and thermal radiation. The governing partial differential equations (PDEs) have been subjected to a similarity transformation, resulting in a set of nonlinear ordinary differential equations (ODEs). These ODEs were solved numerically utilizing the code of bvp4c in the software of Matlab which offers high accuracy (4th order). The employed nanofluid model incorporates the effects of Brownian motion and thermophoresis. The present study illustrates the graphical depictions of the impacts of different governing parameters, namely Hartmann (M) number, curvature (γ) parameter, Brownian motion (Nb) parameter, mass suction (S) parameter, thermal radiation (Rd) parameter, and thermophoresis (Nt) parameter, on heat transfer, flow, and mass transfer characteristics. Comprehensive determination and visual presentation of the coefficient of skin friction, local Nusselt number, and local Sherwood number were conducted for a range of estimates of applied parameters. Based on our examination, it has been determined that dual similarity solutions are present within a specific range of mass suction parameters. The relationship between the Casson parameter and various fluid dynamic properties, such as skin friction coefficient, heat transfer rate, and mass transfer rates, has been found to exhibit a decreasing trend. Furthermore, the stability analysis discovered that the first solution exhibits linear stability, whereas the second solution displays linear instability. Additionally, the motivation behind this study is to enhance industrial performance through the optimization of thermal power generation systems, thereby increasing their overall efficiency.


Introduction
Nanofluids exhibit enhanced characteristics, such as increased temporal stability, reduced flow passage obstruction, and enhanced conductivity of thermal, in contrast to simple base fluids due to nanoparticle inclusion within the fluid medium.Hence, these fluids exhibit a wide range of practical applications in various fields, including pharmacological delivery systems, solar collectors, diabetic treatment through peristaltic pumps, electronic cooling, pharmacological administration mechanisms, and atomic applications.The introduction of the nanofluid concept was pioneered by Choi and Eastman [1], who proposed the blending of a convection fluid with solid nano-particles of extremely small size.Meanwhile, Khan and Pop [2] were the early adopters of a mathematical way to examine the impacts of Brownian motion and thermophoresis of nanofluid on a stretched sheet.Furthermore, the examination of nanofluids in the context of flow over stretching/shrinking cylinders offers a captivating opportunity to examine the augmentation of mass and heat transfer phenomena.Buongiorno model [3] effectively represents the complex interaction among different forces occurring at the nanoscale.Nanoparticles disperse across the fluid medium due to Brownian motion caused by heat variations.In addition, nanoparticles move in response to thermal differences, a process known as thermophoresis, which happens when temperatures change.This movement helps to enhance heat transfer.Thermal control in electronics, solar collectors, and medical therapies are just a few of the many applications that engineers and researchers have investigated using the Buongiorno nanofluid model.Nanofluids, when added to these systems, may improve heat transfer rates and maximize efficiency.Applying the Buongiorno model to the phenomenon of expanding and contracting cylinders has been the subject of a large amount of scholarly work.For various solutions, Rana et al. [4] used Buongiorno's model and Nielsen's effective viscosity model together.They found that the Nusselt number and skin friction both rise noticeably with increasing nanoparticle concentration.Shatnawi et al. [5] extended the above concept for the hybrid nanofluid for Casson model.Zaimi et al. [6] expanded upon Buongiorno's model by applying it to a shrinking cylinder.Ros ¸ca et al. [7] employed the identical Buongiorno model to investigate the flow of nanofluid over a shrinking/stretching cylinder.Shatnawi et al. [8] and Abbas et al. [9] used the model for non-Newtonian flow.Usman et al. [10] utilized the Buongiorno model on an inclined plane and discovered solution.
Fluids can be primarily categorized into two groups: Newtonian and Non-Newtonian.Examples of Newtonian liquids contain water and engine oil, which exhibit a linear stress-strain relationship.On the other hand, non-Newtonian fluids such as honey, paste, and paint demonstrate a nonlinear stress-strain relationship.We frequently encounter these types of fluids in our everyday experiences.Therefore, the examination of these fluids holds significant importance.The current research endeavors will center on the examination of non-Newtonian fluid.Several mathematical models have been put forth to describe the phenomena of flow, heat transfer, and mass transfer in non-Newtonian fluids.Several fluid models have been proposed in the literature, including the Prandtl fluid model [11,12], the Casson fluid model [13], the Carreau fluid model [14], and the Maxwell fluid model [15].Furthermore, it is worth noting that the blood of humans is classified as a type of Casson fluid.The study conducted by Shah et al. [16] examined the flow behavior of Casson nano-fluid over the stretching sheet, taking into consideration the influence of activation energy and chemical reaction.The investigation conducted by Oyelakin et al. [17] examined the behavior of gyrotactic micro-organisms within the flow of Casson nanofluid.Authors Reza et al. [18] employed analysis of finite difference to investigate the unsteady MHD flow of a Casson fluid.The slip conditions of unsteady flow Casson nanofluid on a surface were examined in the studies undertaken by Ghadikolaei et al. [19] and Jamshed et al. [20].Abbas et al. [21] investigated the Casson hybrid nanofluid model numerically where it has been observed that "the skin friction of Casson nanofluid increased and temperature gradient declined due to the increment of solid nanoparticle concentration".Ali et al. [22] employed a non-Newtonian model of micropolar nanofluid and observed that an increase in material parameters leads to a decrease in the compactness, concentration, and temperature.In their research articles, Dero et al. [23,24] investigated non-Newtonian models, revealing the presence of dual solutions within specific ranges of the applied parameters.In their investigation of the Casson fluid flow model, Lund et al. [25] identified dual solutions.
In the year 1942, Haanas Alfren made a significant contribution to the field of MHD by introducing relevant terminology.Numerous scholars have conducted extensive research to comprehend the characteristics of MHD and examine its effects on various aspects of nanofluid.Presently, this phenomenon finds applications in diverse domains including astrophysics, medical science, geography, and numerous others.MHD is a field of study that explores the interplay between fluid dynamics and magnetic fields.Analogies of this nature are conducted for fluids that possess electrical conductivity and lack magnetic properties.These fluids include hot ionized gases, strong electrolytes, and liquid metals.The analysis of the flow of MHD has extensive utilizations in various fields including engineering, industry, and technology.These applications encompass generators of MHD, thin magnetic fluid film stabilization, electric motors, transformers, mass spectrometers for designing cooling systems, and more.The manipulation of the magnetic field allows for the modification of numerous physical properties exhibited by such materials.The study conducted by Hayat et al. [26] examines the influence of the activation energy of Arrhenius on a nonlinear stretched sheet in the existence of convective third-grade nanofluid in MHD flow.The treatment of nanomaterials, irrespective of the presence of MHD streamline, has been examined in the context of a melting sheet by Dinh et al. [27].The examination accomplished by Rashid et al. [28] studied the flow of MHD nanofluid on the formation of porous with shrunk walls of entropy.The research carried out by Pavar et al. [29] inspected the MHD current of nanofluid flowing semi-infinitely through a vertical permeable plate.In their study, Turkyilmazoglu [30] investigated the MHD flow that arises from a rotating disc undergoing a shrinking process.The flow behavior of nanofluid in the existence of the field of magnetic was investigated by Nagaraju et al. [31] via the ramped surface model.Abbas and Shatanawi [32] investigated the effects of Brownian motion and thermophoresis on non-Newtonian fluids.The same Buongiorno model was extended for temperature-dependent non--Newtonian fluid flow over the Riga sheet [33].In their study, Kumar et al. [34] investigated the influence of MHD effect on Casson nanofluid, incorporating considerations for thermal radiation and Hall current.In their study, Azmi et al. [35] expanded upon the MHD concept to encompass Casson fluid, investigating the implications of first-order slip conditions.In their study, Mahabaleshwar A.M. Soomro et al. et al. [36] explored the axisymmetric flow of a non-Newtonian fluid, examining the influence of MHD on the process.Their findings revealed a positive impact of MHD on heat transfer within the system.In their study, Lund et al. [37] incorporated the MHD effect while investigating the behavior of Casson fluid over a vertical surface.The findings revealed that the Hartman number exhibited opposite directions in both solutions.
The latest and intriguing article provides a compilation of significant references about MHD flow on diverse flow geometries [34][35][36][37].
The significance of convective heat transfer in various industrial and environmental applications, such as gas turbines, energy storage, photovoltaic panels, geothermal reservoirs, nuclear plants, and rocket propulsion, has generated considerable attention.The convective boundary condition (BC) has garnered attention in the literature, often being modeled using a Biot number within the thermal BC of the wall.Furthermore, a Biot number less than one signifies that the conduction resistance within the fluid is significantly smaller in comparison to the convection resistance at the surface.The Biot number is defined as unity, indicating that both conduction within the fluid and convection at the surface exhibit similar levels of significance.A Biot number greater than one directs that conduction resistance within an object is substantially higher than convection resistance at a sheet.In the latest examination undertaken by Shahzad et al. [38], authors examined the similarity solutions of the flow and heat transfer characteristics of power law fluid on the stretched cylinder.The analysis also incorporated the effect of a convective BC.In their study, Rashad et al. [39] communicated a similarity solution for convective surface BC of the Williamson hybrid nanofluid.Aziz [40] conducted a study on the characteristics of thermal and hydrodynamic slip flow boundary layers, employing a temperature BC.In their study, Gupta et al. [41] employed a variational finite element method to model the behavior of radiative mixed convective micropolar flow on shrunk surface.The simulation incorporated a convective BC.Analysis undertaken by Bég et al. [42] examined the influence of buoyancy and Biot number effects on MHD slip flows.
This research paper focuses on investigating the MHD flow of a Casson nanofluid around a permeable cylinder that undergoes stretching or shrinking.Enhanced thermal conductivity, reduced friction, and exceptional stability are just a few of the benefits of using Casson nanofluids for heat transfer and viscosity controlling.This research has real-world implications in many different areas, such as heat exchangers for industrial use, solar thermal systems, automobiles, and electronic cooling.Because these allow for controlled thermal treatments and medical imaging, these fluids are vital to biomedical devices, which are ideal for many uses for controlled and efficient heat transfer because of their impressionability in rheological characteristics.Current research mathematical model also depends totally on Buongiorno's [3] nanofluid.It is important to note that Hayat et al. [43], Bakar et al. [44], and Ros ¸ca et al. [45] have investigated the effect of radiation on nanofluid over the permeable cylinder during stretching and shrinking.It should be emphasized that, as far as we are aware, the investigation of stretching/shrinking cylinder in a magnetized Casson nanofluid for multiple solutions with stability analysis, utilizing the mathematical nanofluid model introduced by Buongiorno, was not previously explored.Hence, the prime aims of the current investigation endeavor are to examine the influence of MHD and Biot number on Casson nanofluid, anticipate the presence of multiple solutions, and conduct a stability analysis to ascertain the solution that exhibits stability.

Mathematical modeling
This study focuses on analyzing the characteristics of a specific fluid flow under following assumptions.
1. Flow is the incompressible, steady flow of Casson nanofluid.2. Surface is stretching/shrinking cylinder of radius R along with velocity U w = cx L , where c is continuous with c < 0 corresponding to shrunk constant and c > 0 corresponding to stretched constant, x is coordinate with cylinder, and L is the characteristic length.

Mass suction is assumed as
where S > 0 denotes suction and S < 0 denotes injection.
4. Temperature of surface is T f , and T ∞ is ambient temperature, where T f > T ∞ . 5. The magnetic field is oriented perpendicular to the direction of the fluid's movement (see Fig. 1 (a, b)).
∂C ∂r ∂T ∂r where u and v are the associated components of velocity, and r and x are the coordinates measured along the radial direction and cylinder's surface, respectively.Next, where denotes electrical conductivity, B 0 denotes the applied magnetic field, T denotes the boundary layer temperature, C denotes the concentration of nanoparticles therein, ϑ denotes kinematic viscosity, D B denotes Brownian diffusion, and D T denotes thermophoresis diffusion.In addition, α = k , where ρ f is the fluid density, ρ p is the particle density, and c is the coefficient of volumetric expansion.As for the BCs (refer to Hayat et al. [43] and Ros ¸ca et al. [45]), Under the stretching/shrinking cylinder, a flow with a coefficient of convective heat transfer h f and a temperature T f is present.By applying the similarity transformations, we search for solutions in similarity to Eqs. ( 1)-( 4) subject to BCs (5).It is important to acknowledge that within the field of mathematical modeling, transformations are employed to convert a system into an approximate system.Similarly, similarity transformations are utilized to decrease the number of independent and dependent variables in a given system of equations.In this context, we seek a suitable similarity transformation as proposed by Ros ¸ca et al. [45].
where η is the similarity variable.The BCs at r = a are transformed into BCs at η = 0, which simplifies the numerical calculations.
The following nonlinear ODEs are obtained by putting Eq. ( 6) into Eqs ( 2)-( 4). ( put under the influence of (5) BCs, The skin friction coefficient C f , the local Nusselt number Nu x , and the local Sherwood number Sh x are the physical quantities of interest, and they are specified as where the shear stress on the wall is denoted by τ 1 , the local heat flow is denoted by q w , and the local mass flux is denoted by q m .
Above quantities are calculated using the similarity variables (6).
where Re x = U∞x ϑ is the local Reynolds number.

Stability analysis
In cases when there are more than two possible solutions, a stability analysis must be conducted.Dual solutions have been obtained in this research.In order to ensure the reliability of the solution, it is necessary to perform a stability analysis of Equations ( 7)-( 9) with respect to the boundary conditions (10).We have modified Equations ( 2)-( 4) regulating the unsteady flow by adding a time-dependent variable, denoted by τ = c L .t.

Numerical method
The system represented by equations ( 7)-( 9) is solved through the utilization of the bvp4c procedure, employing the three-stage collocation formula.The collocation polynomial generated yields a consistently accurate fourth-order solution, maintaining C 1 continuity throughout the integration interval [44].The collocation method employs a mesh to discretize the interval into smaller subintervals.The solver ensures the provision of a solution for the governing system.To evaluate the accuracy of the solution, a solver is applied to each subinterval to predict the error as explained in Ref. [45].If the required tolerance is not met, the procedure is iterated with mesh adjustment.This method entails converting BVP into a corresponding IVP.As a result, we formed.
along with BCs { The system of first-order ODEs (28)(29)(30)(31) is numerically solved by specifying initial points on a mesh and appropriately adjusting the initial approximations.

Findings and analysis
The governing ODEs (7-9) and corresponding BCs (10) were solved using the three-stage Lobatto III-A formula implemented in bvp4c code within Matlab software which offers high accuracy (4th order).Results obtained were subsequently compared with the A.M. Soomro et al. findings reported in a previously published paper by Khashi'ie et al. [47] to assess the reliability of the methodology employed.The comparison revealed a significant similarity, as indicated in Table 1.This paper focuses on the characterization of velocity f ′ (η), temperature θ(η), and concentration ∅(η) profiles in the context of a shrinking cylinder, as it is important to note that multiple solutions can be observed in this scenario.In current consideration, we managed to obtain two different solutions within the given parameter ranges: 5 ≤ β ≤ ∞; S c ≤ S ≤ 3; 0 ≤ γ ≤ 0.2; 0.5 ≤ Bi ≤ 2; 0.1 ≤ N t ≤ 0.7; 0.1 ≤ N b ≤ 0.7; 0 ≤ Rd ≤ 1; and 0.07 ≤ M ≤ 0.15.It should be noted that while the first solution is easily found by assigning actual boundary conditions, the trial-and-error basis technique was used to obtain the second solution in boundary conditions by assuming different guesses, subject to all boundary conditions being met.
According to Wang [48] and other researchers, it has been suggested in their analyses that the attainment of a similarity solution for problems of fluid flow on the shrunk sheet is feasible through the application of adequate wall mass suction.The behavior of Newtonian fluid is distinct from that of non-Newtonian fluid.In the case of Casson nanofluid, it has been observed that a significant mass suction is necessary to achieve a solution when the value of the Casson parameter β decreases.This study has revealed the presence of two distinct regions for similarity solutions, specifically dual solutions and no solution, which are contingent upon the mass suction parameter.For the values of β equal to 5, 10, and ∞ (representing the Newtonian fluid), a set of dual solutions is observed when S is greater than or equal to S c1 , S c2 , and S c3 , respectively.Conversely, a single similarity solution is present when S is less than S c1 , S c2 , and S c3 , as depicted in Figs.2-4.When the level of suction is increased, skin friction coefficient (Re 1/2 x C f ) exhibits an increase in the first solution, while in the second solution, it initially experiences a decrease followed by an increase.Figs. 3 and 4 were generated in order to investigate the impact of suction S and Casson parameter β on Nusselt number (Re − 1/2 x Nu x ) and Sherwood number (Re − 1/2 x Sh x ), correspondingly.In all graphical representations, the augmentation of suction results in an enhancement of both the Nusselt number and Sherwood number across all solutions.
As, is already stated the skin friction coefficient (Re 1/2 x C f ), Nusselt number (Re − 1/2 x Nu x ), and Sherwood number (Re − 1/2 x Sh x ) are physical quantities that are of great interest due to their significant engineering applications.According to the observations made in Fig. 5, it is followed that coefficient of skin friction tends to surge as the parameter S increases.Conversely, the value of Re 1/2 x C f reduces as curvature γ parameter increases in the first solution.This can be attributed to the fact that greater values of curvature γ parameter negatively impact the shear stress.The physical phenomenon occurs as a result of the decrease in parameter S, indicating to a decline in fluid velocity and a corresponding increase in friction force at the wall.Conversely, an inverse pattern is observed initially for the second solution.Conversely, it is observed that the ranges of dual solutions and the no solution can also be identified.Based on the analysis of Figs. 6 and 7, it can be inferred that for γ values of 0, 0.1, and 0.2, there exists a range of dual solutions when the value of S is greater than or equal to S c1 , S c2 , and S c3 , respectively.Conversely, no similarity solution is observed when S is less than S c1 , S c2 , and S c3 , respectively.Fig. 7 illustrates the fluctuations in the gradient of concentration ∅(η) profile across several values of S and γ.It is observed that rate of concentration transfer for both solutions exhibits an increase as the value of S is incremented.
Fig. 8 illustrates the fluctuations observed in the Biot number (Bi) across the temperature profile θ(η).The results clearly demonstrate that a increase in Biot number outcomes in an augmentation of both θ(η) distribution and thermal boundary layer (TBL) for both solutions.In practical terms, the feasibility arises from the relationship involving the Biot number, which quantifies the relationship between internal and external thermal resistances.As Biot number improves, the external thermal resistance decreases, leading to a thicker TBL.Fig. 9 depicts the concentration profile ∅(η) corresponding to varying Biot numbers.It has been observed that there is an increase in concentration in both the first and second solutions.From a physical standpoint, it can be argued that an elevated Biot number has a positive impact on convective heat transfer, facilitating improved dispersion and mixing of nanoparticles in the nanofluid.Consequently, this results in the increased ∅(η) profile.
Fig. 10 illustrates a positive correlation between θ(η) profile and thermophoresis parameter N t in both first and second solutions.The phenomenon that has been noticed can be explained by to N t , and this is what propels the nanoparticles away from the surface and into the surrounding flow.Consequently, this leads to an increase in TBL.Fig. 11 illustrates the ∅(η) profile as N t increases.The concentration of the first and second solutions is observed to increase as the value of N t increases.The physical explanation for this phenomenon is that a large number of nanoparticles are moving toward the surrounding flow, which raises the concentration of these particles and causes the concentration boundary layer to grow.Fig. 12 illustrates the variation of the Brownian motion N b with respect to the temperature profile θ(η).The relationship between the values of N b and fluid temperature in the first and second solutions exhibits a clear positive correlation.Nevertheless, the

Table 1
Contrast of f ″ (0) for different magnitudes of S and γ when Bi = Rd = 0.5, λ = − 1, M = 0.1, Pr = 2, Sc = 2, N b = 0.5 = N t , and β = ∞.γ Khashi'ie et al. [47] Current Results         an increase as the enhancement of Rd occurs in both solutions.Nevertheless, as Rd raises, one sees a matching rise in the thickness of TBL in both solutions.The observed variations in the θ(η) distribution can be attributed to the influence of radiative heat transfer, specifically the net energy flux resulting from thermal radiation.According to Fig. 15, it can be observed that there is an increase in ∅ (η) for solutions with larger values of Rd parameter in both the first and second solutions cases.The velocity profile f ′ (η) for various Hartman number M is depicted in Fig. 16.f ′ (η) is noticed to reduce as M magnitude rises in the first solution.The improvement in the intensity of the magnetic field results in the generation of a resistant force recognized as the Lorentz force, which signifies the deceleration of f ′ (η).Velocity and momentum boundary layer's thickness both exhibit an increase as the Hartmann number (M) is augmented in the second solution.The temperature profile θ(η), as depicted in Fig. 17   number M is increased.There is a lack of significant observable alteration.In contrast, the ∅(η) profile exhibits an increase as the Hartmann number M is enhanced in the second solution.This behavior can be attributed to the higher Hartmann numbers inducing more pronounced magnetic effects, which in turn alter the movement of nanoparticles and fluid flow, consequently modifying the ∅ (η) profile of the nanofluid.

S
Table 2 presents the smallest eigenvalues (ε 1 ) for both solutions, considering various values of S and curvature parameter (γ).The elementary objective of our study is to assess the stability of the solutions in order to determine the stable solution among the two solutions.It is imperative to acknowledge that the physical interpretation of a solution is contingent upon its stability.From a mathematical perspective, it is essential to consider all potential multiple solutions, if they exist, as they are integral to the overall fluid flow problem and satisfy all necessary solution criteria.It is important to acknowledge that the stability of solution is contingent on the smallest eigenvalue sign.Findings presented in Table 2 demonstrate that ε 1 exhibits negative (positive) values.These values signify an initial growth (decay) of the disturbance, which in turn disrupts (resumes) the separation of the boundary layer and influences the stability (instability) of the flow.It is noteworthy to mention that the first solution, known as the stable solution, consistently yields a  meaningful physical interpretation that can be observed.

Conclusion
This article exhibits an examination of the steady flow, heat transfer, and mass transfer characteristics of a magnetized Casson nanofluid over the stretched/Shrunk cylinder.The Buongiorno model has been modified to examine the influence of Rd, N t and N b effects.Viable transformations are employed in order to derive a set of simplified ODEs.The bvp4c solver function in Matlab software is employed to investigate the solutions of these ODEs.The following key findings are enumerated: 1. Dual solutions are observed exclusively on a shrinking cylinder.2. There exists a range of dual solutions when the value of S is greater than or equal to S c1 , S c2 , and S c3 .3. The inclusion of γ has a significant effect on profiles of f ′ (η), θ(η), and ∅(η).

4.
Local Nusselt number exhibits an expanding trend with an upsurge in S, while it demonstrates a decreasing trend with an increase in γ. 5. First solution demonstrates an inverse correlation between the f ′ (η) profile and magnetic M parameter.6. Profiles of temperature and concentration increase in both solutions for the rising effects of biot number.7. The temperature and concentration profiles of both solutions exhibit an increase due to the escalating influence of thermal radiation.8. Sherwood number exhibits a rising trend as the suction parameter is expanded in both the second and first solutions.9. Stable solution refers to the first solution that satisfies the physical constraints.
The findings of current paper hold significant associations in the domain of heat and mass transfer, specifically in relation to the utilization of non-Newtonian Casson nanofluids.These results also present several promising avenues for further research.Casson nanofluids exhibit distinct characteristics that render them suitable for a wide range of applications, such as medical utilization in the       context of blood flow, enhancement of thermal energy technologies within the energy sector, and utilization in solar energy systems to power aircraft and ships.Potential areas of investigation for future studies may encompass the examination of entropy generation, multi-scale analysis, experimental validation, micro-morphic constitutive, and hybrid nanofluid model.

Fig. 1 .
Fig. 1. (a, b): Physical model consists of two cylinders: one that stretches (a) and one that shrinks (b).
A.M.Soomro et al.

Fig. 2 .
Fig. 2. Changes to Re 1/2x C f for various S and β values.

Fig. 5 .
Fig. 5. Changes to Re 1/2x C f for various S and γ values.

A
.M.Soomro et al.   heightened activity of nanoparticles results in a rise in thickness of TBL.Consequently, these nanoparticles have a substantial effect on behavior of the boundary layer and induce alterations in the θ(η) distribution in far proximity to the surface.Fig.13illustrates the impact of N b on the concentration profile ∅(η).A decline in behavior has been noted in the initial and subsequent solutions as the values of N b increase.This phenomenon can be attributed to the correlation between the increased value of N b and the heightened random motion occurring within the boundary layer.Consequently, this augmented motion does not result in a proportional increase in resistance, thereby leading to a reduction in the concentration of nanoparticles.The impact of thermal radiation (Rd) on the dimensionless temperature profile θ(η) is demonstrated in Fig.14.θ(η) profile exhibits
, is presented in relation to the Hartmann number M. Second solution exhibits a direct correlation between the Hartmann number and both TBL thickness and the temperature.The increase in Hartmann number leads to heightened resistance within the boundary layer, resulting in an elevation of the nanofluid's temperature.The observed phenomenon in the first solution is the declining trend exhibited by θ(η) profile and TBL.The analysis of Fig. 18 reveals that the thickness of concentration decreases in the first solution as the Hartmann