Modelling of Applied Magnetic Field and Thermal Radiations Due to the Stretching of Cylinder

: In this study, a numerical approach was adopted in order to explore the analysis of magneto ﬂuid in the presence of thermal radiation combined with mixed convective and slip conditions. Using the similarity transformation, the axisymmetric three-dimensional boundary layer equations were reduced to a self-similar form. The shooting technique, combined with the Range–Kutta–Fehlberg method, was used to solve the resulting coupled nonlinear momentum and heat transfer equations numerically. When physically interpreting the data, some important observations were made. The novelty of the present study lies in ﬁnding help to control the rate of heat transfer and ﬂuid velocity in any industrial manufacturing processes (such as the cooling of metallic plates). The numerical results revealed that the Nusselt number decrease for larger Prandtl number, curvature, and convective parameters. At the same time, the skin friction coefﬁcient was enhanced with an increase in both slip velocity and convective parameter. The effect of emerging physical parameters on velocity and temperature proﬁles for a nonlinear stretching cylinder has been thoroughly studied and analyzed using plotted graphs and tables.


Introduction
Computational experiments in natural convection heat and mass transfer through thermal radiative and hydromagnetic flow have recently overwhelmed numerous scientists due to their importance in industrial technology [1]. This research has various industrial applications including the aerodynamic extrusion of plastic sheets, the cooling of metallic plates, energy storage devices, oil reservoirs, and micro magnetohydrodynamics (MHD) pumps [2].
To begin, [3] investigated the boundary layer flow (BLF) on continuous solid surfaces. They investigated whether the new boundary layer class had a different solution sustainability than the BLF on finite surfaces. In the presence of physical conditions, such as the stretching of surfaces, significant properties of flow and heat transfer have been investigated in the literature. A hybrid analytical and experimental investigation related to velocity-temperature fields in the boundary layer due to the stretching of surfaces has been conducted by [4]. Ref. [5] later applied this definition to the BLF in 1970. Crane studied fluid flow over the stretchable surface and found a similar solution for two-dimensional steady flow. Additionally, [6] examined the two-dimensional model and proposed similar solutions for the stretching of both the channel and cylinder. Ref. [7], on the other hand, found an exactly similar solution to the Navier-Stokes s equations concerning stagnation point flow against a stretching surface. Ref. [8] investigated hydromagnetic free convective heat transfer over a stretching surface embedded in a porous medium when heat generation or absorption was present. Meanwhile, [9] addressed the movement of viscoelastic flow near the stagnation point.
One of the most influential research fields in modern applied mathematics is fluid flow studies in mechanics due to nonlinear stretching [10]. Apart from air or water, the majority of fluids used in different processes have nonlinear viscosities [11]. Ref. [12] discussed stretching cylinders with both active and passive-controlled nanoparticle concentration conditions. They observed that the behavior of fluid momentum changes in the case of stretching cylinders and sheets was due to the nanoparticles' concentration. Ref. [13] studied heat transfer in the presence of hybrid nanofluid and found that these types of fluid transfer heat more effectively than conventional fluids. Ref. [13] unveiled the influence of the shapes of Al 2 O 3 nanoparticles in an aqueous medium flowing in the porous space between eccentric cylinders. Ref. [14] studied finite surfaces with localized necking as bifurcation from deformation with wave numbers. They built an approach and validated that the necking problem is treated as a bifurcation with zero wave numbers. Ref. [15] analyzed the MHD flow of a Casson fluid over a stretching cylinder numerically. They revealed the reverse behavior of velocity and thermal fields by increasing the Hartman number. Nonlinear liquid connections are useful for understanding the flow behavior of commonly found liquids in nature and industry. As a result, various liquid models have been proposed to account for the popular characteristics of nonlinear liquids. Considerable attention towards convective flow and heat transfer in has been seen in recent years due to their nature and potential engineering applications. Ref. [16] examined two-dimensional flow with non-uniform heat generation with slip conditions. It was observed that an increase of the Eckert number also increased the thermal boundary layer and temperature distribution. The available literature on convective flow over stretching surfaces shows that the combined effects of nonlinear stretching, MHD, thermal radiation, and slip effects for distributions has not been investigated for Newtonian fluids.
Many researchers have published their findings on the significance of slip conditions in regards to the phenomenon of the lack of fluid adherence to surface borders. Slip has an important impact on the thickness and momentum of the thermal boundary layer. Furthermore, partial slip conditions are necessary in some practical circumstances such as fluid flow in micro-electro-mechanical devices. Some researchers have published their findings on the significance of slip states. Ref. [17] investigated heat transfer under slip conditions and found that as the Brownian motion and thermophoresis parameters increased, so did the local Nusselt number. Ref. [18] had obtained an exact solution by considering the effect of a second slip on nanofluids peristaltic flow in an asymmetric channel. Ref. [19] had examined Hall current, Newtonian heating, and second-order slip effects on convective magneto-micropolar fluid flow over a sheet. Furthermore, increases in the slip decreased the horizontal velocity profile. Ref. [20] investigated second-order slip and heat transfer models over permeable stretching/shrinking surfaces. They discovered that there are multiple solutions for certain combinations of parameters, with an increase in slip value in the case of shrinking sheets. Ref. [21] investigated dual and triple solutions for an electrically conducting MHD slip flow and a non-Newtonian fluid past a shrinking sheet.
In natural flows, the MHD effects can be described by Navier-Stokes s equations for fluids. These equations describe the motion of a conducting fluid flowing under an external magnetic field s influence and, therefore, are of much importance. Ref. [22] investigated the steady laminar MHD mixed convection stagnation-point flow of an incompressible Processes 2021, 9, 1077 3 of 15 viscous fluid. They found that dual solutions exist in the presence of external magnetic fields both for assisting and opposing flows. Ref. [23] performed analyses for that heat and mass transfer for MHD viscoelastic flow with Ohmic dissipation. He had found that the magnetic parameter reduces heat transfer efficiency. Ref. [24] had studied second-grade MHD fluid in a porous medium with vorticity and Hall current effects. They concluded that the Newtonian flow field could be recovered by setting second-grade parameters equal to zero. Ref. [25] investigated the MHD flow of a micropolar liquid towards a nonlinear stretched surface using viscous dissipation, Joule heating, and a convective boundary state. They analyzed local Nusselt number enhances for the micropolar parameter, Biot, and Prandtl number. Ref. [26] used Buongiorno s model with nanoparticles for heat and mass transfer in the presence of MHD flow. They proposed using a passive control model to control the thickness of the boundary layer and the rate of mass transfer at the wall. The Buongiorno model was used to investigate MHD non-Newtonian Maxwell fluid with nanomaterials in an exponentially stretching surface [27]. They discovered that the thermophoresis parameter has a negligible effect on the temperature profile. Because of the curved stretching sheet, researchers investigated mass and heat transfer in MHD flow of an electrically conducting non-Newtonian micropolar fluid. Ref. [28] discovered that magnetic field interaction increased fluid temperature while decreasing microrotation and velocities at all domain flow points. Ref. [29] studied the thermal analysis of radiative bioconvection MHD flow comprising gyrotactic microorganism with activation energy. They concluded the density of motile microorganisms decreases with Peclet number and bioconvection Lewis number.
The current study investigates nonlinear extended flow for MHD viscous fluid flow over a cylinder in conjunction with nonlinear thermal radiations and slips effects for both velocity and temperature distributions. Therefore, an attempt is made to investigate the resulting nonlinear boundary layer problem numerically. The governing partial differential equations of the physical model have been reduced to ordinary differential equations (ODEs). The Runge-Kutta-Fehlberg method, in conjunction with shooting procedures, was used to numerically solve the resulting system of equations. Furthermore, graphical findings are shown and analyzed.

Mathematical Formulation
The behavior of two-dimensional incompressible boundary layer Newtonian fluid flow is apprehended in this article. In two-dimensional cylindrical flow, the radius of the cylinder is greater than the thickness of the boundary layer [30]. The following conditions were applied and considered from the applications point of view (i.e., MHD): velocitytemperature slips, nonlinear thermal radiations, and stretching. Nonlinear stretching cylinder of radius R in the presence of a non-uniform magnetic field B = B 0 (z) n−1 2 . The flow is confined in the region r > R, where z is the horizontal axis of the cylindrical coordinate system (and is taken along the cylinder s axis) while r and θ are radial and tangential components, respectively. The flow is generated by stretching the cylinder under the influence of two forces equal in magnitude and opposite in direction along the z-axis, keeping the origin fixed. We suppose that at the time, t < 0, the fluid is at rest. The cylinder is stretched with the velocity U(z) = U 0 z n , where U 0 and n are constant.
The assumed nonlinear temperature T w (z) = T ∞ + bz 2n−1 function of z. Free stream velocity is u e = u 0 z n where U 0 >> u 0 ambient fluid temperature T ∞ , thermal diffusivity α, thermal expansion coefficient β, and the magnetic field are applied traversal to the cylinder s axis. Flow description and coordinate axis are shown in Figure 1.
traversal to the cylinder′s axis. Flow description and coordinate axis are shown in Figure  1. The boundary layer′s assumptions are taken as mentioned [15,31,32]. Governing flow and heat transfer equations with the Boussinesq approximation are: Velocity slip L , thermal slip length S , and stretching index n are constants. The supposed temperature difference within the BLF is such that 4 T may be expanded in Taylor′s series about ∞ T and neglect the higher terms. Thus, we get, we acquired [33], The boundary layer s assumptions are taken as mentioned [15,31,32]. Governing flow and heat transfer equations with the Boussinesq approximation are: With corresponding boundary conditions (t = 0), Velocity slip L, thermal slip length S, and stretching index n are constants. The supposed temperature difference within the BLF is such that T 4 may be expanded in Taylor s series about T ∞ and neglect the higher terms. Thus, we get, and By substituting Equation (5) in Equation (6) we acquired [33], Here, the mean absorption coefficient k * , Stefan-Boltzmann constant σ * and Quinn Brewster heat flux q r given by [34] using Rosseland approximation.
Here, surface shear stress and heat flux of the stretching cylinder is denoted by τ w and q w , respectively. Applying the similarity transformation on Equations (12) and (13) achieves,

Numerical Method
The present model contains the nonlinear dimensionless coupled ODEs in Equations (9)-(13) with sufficient boundary conditions. Both Runge-Kutta-Fehlberg s fourthand fifth-order (RKF45) and the shooting scheme were used to solve the ODEs [13,35]. To solve these coupled nonlinear differential equations, we converted higher-order nonlinear ODEs into first-order ODEs. Following that, the applicable initial solutions (zeroth order) were chosen, which must validate the boundary conditions under consideration. The obtained solutions demonstrated that a variety of relevant parameters affect velocity and temperature distributions. The shooting approach was repeated in order to achieve the convergence criteria so that no numerical oscillations would occur. The ∆η = 0 .01 step size was taken, and the selection η ∞ varied from 3 to 5.
boundary conditions,

Model Validity
It is significant to validate and confirm the accuracy of the present numerical procedure. This model's validity is confirmed by comparing the special case with previous studies in the literature. In reference to other studies, our results agreed well with them (as mentioned in Tables 1 and 2). The present results of the local Nusselt number in the Newtonian fluid, compared with the existing data for different temperature exponent values, are present in Table 1. The skin friction coefficient numeric errors are also considered in Table 2. An agreement was found between both physical phenomena and the publicized data.

Results and Discussions
The non-linearity parameter n for the velocity profile has been taken into account in Figure 2. Here, the velocity profile decreased when the fluid exhibited heavy oscillating behavior. Furthermore, for greater values of n, the velocity profile became zero in the short-range.

Results and Discussions
The non-linearity parameter n for the velocity profile has been taken into account in Figure 2. Here, the velocity profile decreased when the fluid exhibited heavy oscillating behavior. Furthermore, for greater values of n , the velocity profile became zero in the short-range In Figure 3, an increase in the curvature made the cylinder cross-section smaller. It enhanced the fluid velocity due to the small surface area. A smaller surface means a smaller amount of friction production between the solid and liquid surface. In Figure 3, an increase in the curvature made the cylinder cross-section smaller. It enhanced the fluid velocity due to the small surface area. A smaller surface means a smaller amount of friction production between the solid and liquid surface. In Figure 3, an increase in the curvature made the cylinder cross-section smaller. It enhanced the fluid velocity due to the small surface area. A smaller surface means a smaller amount of friction production between the solid and liquid surface.     Figure 5 represents the velocity-slip effect on the velocity profile. The abrupt decrease in the momentum boundary layer was due to this effect. Here, the measured viscosity was significantly lower than the actual viscosity of the fluid.     Figure 5 represents the velocity-slip effect on the velocity profile. The abrupt decrease in the momentum boundary layer was due to this effect. Here, the measured viscosity was significantly lower than the actual viscosity of the fluid. In Figure 6, the velocity profile is portrayed in order to measure the influence of mixed convection. According to this graph, it was observed that the velocity profile was enhanced with an increase in the value of the mixed convection effect. This was because an increase in buoyancy caused an increase in the momentum of the fluid. In Figure 6, the velocity profile is portrayed in order to measure the influence of mixed convection. According to this graph, it was observed that the velocity profile was enhanced with an increase in the value of the mixed convection effect. This was because an increase in buoyancy caused an increase in the momentum of the fluid.  It reveals that an increment in n reduced the temperature profile. Nonlinear stretching caused a disturbance in the flow, which caused the excess particles of fluid to collide with each other.   It reveals that an increment in n reduced the temperature profile. Nonlinear stretching caused a disturbance in the flow, which caused the excess particles of fluid to collide with each other.   Figure 2). This caused the fluid to move in faster and have less interaction or collision with the fluid molecules. Since there was less collision between the molecules, the small amount of friction generated a lesser amount of heat. This reduced the overall temperature profile.   Figure 2). This caused the fluid to move in faster and have less interaction or collision with the fluid molecules. Since there was less collision between the molecules, the small amount of friction generated a lesser amount of heat. This reduced the overall temperature profile. The effect of the Prandtl number on the temperature profile is shown in Figure 9. The Prandtl number is the ratio of kinematic viscosity to thermal diffusivity that varies with temperature. Owing to this fact, different fluids can have different Prandtl numbers because of variations in their temperature. In this case, the fundamental law is that fulfilling a smaller value of Pr decreases the temperature profile. The effect of the Prandtl number on the temperature profile is shown in Figure 9. The Prandtl number is the ratio of kinematic viscosity to thermal diffusivity that varies with temperature. Owing to this fact, different fluids can have different Prandtl numbers because of variations in their temperature. In this case, the fundamental law is that fulfilling a smaller value of Pr decreases the temperature profile. The temperature profile, as shown in Figure 10, decreased as the value of the thermal slip increased. As the temperature rose, the heat transfer from the solid surface to the fluid decreased. As a result, the temperature boundary layer decreased. Furthermore, although the momentum equation is dependent on ( ) θ η , no important effect γ on velocity profiles was observed.  The temperature profile, as shown in Figure 10, decreased as the value of the thermal slip increased. As the temperature rose, the heat transfer from the solid surface to the fluid decreased. As a result, the temperature boundary layer decreased. Furthermore, although the momentum equation is dependent on θ(η), no important effect γ on velocity profiles was observed. The temperature profile, as shown in Figure 10, decreased as the value of the thermal slip increased. As the temperature rose, the heat transfer from the solid surface to the fluid decreased. As a result, the temperature boundary layer decreased. Furthermore, although the momentum equation is dependent on ( ) θ η , no important effect γ on velocity profiles was observed.   the heated fluid moves away from the heat source or heated surface. It makes the fluid have a lower temperature near the stretch surface, reducing the overall temperature profile. Figure 11 exhibits that temperature profiles decreased rapidly with the increment in mixed convection. It had a prominent effect due to coupled governing equations. In general, the mixed convective heat transfer phenomena, due to fluid motion, occurs when the heated fluid moves away from the heat source or heated surface. It makes the fluid have a lower temperature near the stretch surface, reducing the overall temperature profile. for different values of λ . Figure 12 shows the thermal radiations effect on the temperature profile. It was concluded that an increase in K also increased the temperature profile. It is a common fact that, in order to enhance thermal radiation, the temperature of the surrounding environment must also rise. Therefore, in this case, an excess amount of radiation enhanced the temperature and boundary thickness.  Figure 12 shows the thermal radiations effect on the temperature profile. It was concluded that an increase in K also increased the temperature profile. It is a common fact that, in order to enhance thermal radiation, the temperature of the surrounding environment must also rise. Therefore, in this case, an excess amount of radiation enhanced the temperature and boundary thickness. In Table 3, the skin friction In Table 3, the skin friction − f (0) coefficient showed decreasing behavior for the non-linearity exponent, curvature parameter, and magnetic parameter. The Nusselt number −θ (0) showed decreasing behavior for the non-linearity exponent, curvature parameter, mixed convection parameter, and Prandtl number.

Conclusions
In this study, we investigated the effects of velocity-temperature slip on steady flow and heat transfer over an axisymmetric cylinder with mixed convection. We also investigated both the traversal magnetic field and thermal radiation with non-linearity exponent behavior. By using the appropriate similarity transformations, the governing equations were reduced to ODEs. The shooting technique was used to reduce the ODEs and was solved numerically through RKF45. Graphs and tables were used to present and discuss the effect of main parameters on velocity and temperature profiles, as well as on the skin friction coefficient and Nusselt number. A stability analysis was also performed for the accuracy and validation of the results, and found that the present results are in good agreement with previous studies. In this work, we showed that:

•
The mixed convection parameter enhances the velocity profile.

•
The Prandtl number reduces the temperature distribution across the flow.

•
The velocity profile abruptly decreases in the presence of a velocity slip.

•
The temperature profile diminishes due to the thermal slip effect.