Computational Dynamics of Stagnation Point Flow of Micropolar Fluid Past Vertical Porous Plates

This work examines the flow of a micropolar fluid over a vertical porous plate at the MHD stagnation point under viscous dissipation, convective boundary conditions, and thermal radiation. The governing partial differential equations and a set of similarity parameters were used to transform them into ordinary differential equations. The Runge-Kutta fourth-order algorithm is used in conjunction with the Newton Raphson shooting technique to numerically solve the generated self-similar equations. Results were tabulated both numerically and graphically, and examples for different controlling factors are quantitatively analyzed. According to the study, the vortex viscosity parameter (k) causes the velocity profiles to rise while the magnetic parameter, suction parameter, and radiation parameter cause them to fall. In contrast, as the flow’s suction and prandtl values rise, so do the magnetic parameter, radiation, and vortex viscosity, while the thickness of the thermal boundary layer decreases.


Introduction
Fluids having microstructure are known as micropolar fluids.They fall under the category of fluids known as polar fluids, which have nonsymmetric stress tensors.Erigen [1] first proposed the fundamental continuum theory for this group of fluids and has been a well-liked subject of study.The movement of fluids with suspensions, colloidal fluids, polymer, bodily fluids, blood, and liquid crystals are all explained by this hypothesis.Incompressible micropolar boundary layer flow over a semi-infinite plate was researched by Ahmadi [2].In their investigation of the heat transfer on continuously rotating plates in micropolar fluids, Soundalgekar and Takhar [3] discovered the impact of surface temperature on the fluid dynamics.Similar to this, Hayat and Ali [4] explored the peristaltic flow of micropolar fluid in an asymmetric channel as well as the nature of the endoscope.Rees and Pop [5] studied the flow of a micropolar fluid on a continuously moving plate, whereas Sajid et al. [6] looked at the homotopy analysis for boundary layer flow of micropolar fluid via porous channels.The mixed convection flow of micropolar fluid across a non-linearly expanding surface was further investigated by Hayat et al. [7].
In many technical applications, such as the design of radial diffusers, thrust bearing, thermal oil recovery and transpiration cooling, stagnation point flow is crucial.Hiemenz [8] employed a similarity transformation method to convert the Navier-Stokes equations to nonlinear ordinary differential equations and made the discovery of stagnation point flow.Homann [9] expanded on this issue by incorporating the axisymmetric stagnation point scenario in both two and three dimensional instances.In addition to making a substantial contribution to the stagnation-point flow of micropolar fluid towards a stretched surface, Nazar et al. [10] also added to our understanding of the flow dynamics of stagnation point flow.
The design of several inventive energy conversion devices that run at high temperature has a significant problem from thermal radiation.The emissions from heated walls and working fluid are principally caused by the effects of thermal radiation.Several research investigations, like those by Zhu et al. [11] and Pop et al. [12], have discussed the effects of radiation.While Seini and Makinde [13] looked at the impact of radiation on chemically reacting MHD boundary layer flow via a vertical porous plate, Christian et al. [14] examined MHD stagnation flow with chemical reaction and radiation toward a heated shrinking porous surface.T. Ayando [15] also looked into the Blasius flow of MHD micropolar fluid caused by heat radiation through a permeable plate.Viscous dissipation's function is to alter the temperature distribution by acting as energy source, which impacts heat transfer rates.Whether the sheet is being heated or cooled determines the benefits of the viscous dissipation effect.
The impact of viscous dissipation in a heated vertical plate with natural convection was examined by Pantokratoras [16].The effects of suction and viscous dissipation on MHD boundary layer flow in a porous material over a moving vertical plate were also explored by Lakshmi et al. [17].The effects of viscous and ohmic dissipation on heat transfer and viscoelastic MHD flow across a stretched sheet were next investigated by Subhas et al. [18].Imoro et al. [19] investigated the presence of viscous dissipation and nth order chemical reaction on heat and mass transfer over a vertical surface with convective boundary conditions, whereas Arthur et al. [20] examined the presence of radiation with viscous dissipa-tion and convective boundary condition chemically reacting hydromagnetic flow over a flat surface.Makinde [21] investigated internal heat generation and convective boundary conditions on a moving vertical plate with natural convection using the similarity solution method.
The stagnation point flow of an MHD micropolar fluid in the presence of melting processes and heat absorption was studied by Mamta et al. [22], who found that the heat transfer rate reduces with melting and heat absorption.He said that at the fluid-solid boundary, heat production parameters greatly rise, while Lok et al. [23] looked the mixed convection flow of micropolar fluid at the stagnation point on a vertical surface and Ramachandra et al. [24] studied the mixed convection stagnation point flows close to a vertical porous surface.
Ghasemi et al. [25] investigated the effects of solar radiation on MHD stagnation point flow and heat transfer across a stretched sheet.
Lund et al. [26] investigated the dual similarity solution of MHD stagnation point flow of casson fluid with thermal radiation and viscous dissipation effects.
Hsio [27] studied the stagnation electrical MHD nanofluid mixed convection with slip boundary on stretching sheet, while Bilal [28] examined the micropolar flow of electrical MHD nanofluid with nonlinear thermal radiation and slip effects.Both researchers discovered that increasing the magnetic parameter (M) or electrical parameter (E) results in an increase in temperature distribution at a specific point of the flow region.Furthermore, Shahzada et al. [29] also expanded on Hsiao [27] research by looking at the stagnation point flow of an EMHD micropolar nanofluid with mixed convection and slip boundary.Recent studies by Ishak et al. [30] and Olanrewaju et al. [31] examined the effects of thermal radiation on magnetohydrodynamic (MHD) flow of micropolar fluid towards a stagnation point on a vertical surface and found that thermal radiation and absorption has greater impact on the velocity, angular velocity and temperature field.
To the best of the authors' knowledge, there is no documentation of the interaction of viscous dissipation, thermal radiation, and convective boundary conditions of stagnation point flow of micropolar fluid via vertical porous plate in the literature that is currently accessible.
This study investigates the impact of convective boundary condition, viscous dissipation, and thermal radiation on stagnation point flow of micro polar fluid through a vertical porous plate.The practical applications of this research include the extraction of polymers in melt-spinning processes, the cooling of nuclear reactors during emergency shutdown, and the cooling of electronic devices, serve as the driving force behind it.

Mathematical Model
On a heated vertical surface, a constant laminar two-dimensional flow of a viscous incompressible electrically conducting micropolar fluid through a porous media has been taken into consideration.Our research took into account the tangential and normal velocity components, as well as the x-axis running parallel to the wall in the direction of the flow motion and y-axis perpendicular to it.We also failed to consider the generated magnetic field that the movement of the electrically conducting fluid caused.The following steps are taken to get the boundary layer equations representing the flow: With boundary conditions: where u and v are the x and y-axis velocity components, respectively, g is the acceleration caused by gravity, T is the fluid temperature in the boundary layer, k is the vortex viscosity, μ is the dynamic viscosity, κ is the thermal conductivity, ρ is the fluid density, β is the thermal expansion coefficient, j is the microinertia density, H is the microrotation vector.
We further assume that K k µ = is the vortex viscosity parameter where the field of equations predicts proper behavior of the vortex when the microstructure effects are insignificant.
Using Rosseland approximation for radiative, simplified the radiated heat flux to: where K' and σ * are mean absorption coefficient and Stefan-Boltzmann con- stant.The term T 4 may be expressed as a linear function of temperature difference within the flow.Hence, in a Taylor series expansion about T ∞ and neglecting higher order terms, we get; Introducing the following dimensionless quantities: where ψ, θ, h and υ are respectively stream function, temperature, dimensionless microrotation and kinematic viscosity.Equations ( 2)-( 4) and the boundary conditions in ( 5) are transformed into non-linear higher order differential equations in the form: ( ) ( ) The transformed boundary conditions are; where j a υ = is the characteristic length,

Numerical Procedure
The higher order ordinary differential Equations ( 9) through (11) and the associated transformed boundary conditions (12) are reduced to a coupled first order system of ODEs using the standard Newton-Raphson shooting method and the fourth-order Runge-Kutta integration algorithm to arrive at the numerical solution.We employ numerical shooting technique where the two ending boundary conditions are utilized to produce two unknown initial conditions at at η = 0.In this calculation, the step size ∆η = 0.001 was used while obtaining the numerical solution with η max = 10 and six-decimal (10 −6 ) accuracy as the criterion for convergence.
The following reduction steps are allowed in order: Therefore, equations ( 9) -( 11) can be reduced as first order system as follows: For different prandtl number (pr) values, the studies of Ramachandran et al., Lok   3.According to the results, the shear stress and the couple stress rise with increasing values of fw, Bi, M, Ra, Br and λ (λ > 0), whereas they decrease with rising values of K, Pr and λ (λ < 0).At the plate surface, the shear stress and couple stresses are reduced by momentum and angular momentum diffusion, while they are increased by the magnetic field intensity, suction, viscous dissipation, convective heat transfer and thermal radiation.
The correlation between rate of heat transfer and shear stress is seen in Table 4.
The impact of buoyant forces results in an increase in shear stress and a decrease in the rate of heat transfer at the surface, as shown numerically in the table.Table 1.An example of f"(0) results for various values of Pr with λ = 1, Br = K = Ra = M = 0, and Bi x = 10

Microrotation Profiles
As seen in Figure 8 & Figure 9, the microrotation profiles decrease as the values of magnetic field (m) and suction (fw) parameter rise.The Lorenz force, which tends to work against the direction of the fluid flow when there is a magnetic field present, is a drag-like force.As the fluid temperature, magnetic field, and suction rise as result, the fluid velocity and microrotation decrease.

Conclusions
1) The vortex viscosity parameter (K) causes a rise in the velocity profiles, while other parameters such as the magnetic parameter (M), suction (fw), Brinkmann number (Br), Biot (Bi), and radiation parameter (Ra) cause a reduction.In  contrast, the thickness of the thermal boundary layer rises with the magnetic parameter (M), Brinkmann number (Br), Biot number (Bi), radiation (Ra) and vortex viscosity (K), whereas it decreases with increasing suction and prandtl.
2) The microrotation profiles decreases with increasing magnetic parameter

A.
Timothy et al.DOI: 10.4236/jamp.2023.11112213488 Journal of Applied Mathematics and Physics Noting the usual relationship of the velocity components u y The local Nusselt number, couple stress, and plate surface shear stress, which are each represented by the numerical values f"(0), h'(0) and −θ'(0) are computed numerically and results are shown in tabular form.

Figures 1 -Figure 1 .
Figures 1-7 shows the velocity curve with various regulating parameters.The

Figure 2
Figure 2 shows how the velocity field and suction (fw) values vary over boundary layer.Velocity decreases with increasing suction.Suction is a medium that creates flow resistance, slowing the flow velocity.As shown in Figures 3-5, the velocity boundary layer thickness decreases as the values of thermal radiation, Brinkmann number (Br) and Biot number (Bi) on the velocity increase, respectively, due to viscous dissipation and convective heat transfer.Additionally, as shown in Figure 6 & Figure 7, the thermal boundary layer thickness rises for varying values of the vortex viscosity parameter (k), improving the flow kinematics of fluid flows due to an increase in angular momentum diffusion.

Figure 10 &
Journal of Applied Mathematics and Physics

Figure 11 show
Figure 11 show the microrotation profiles for the Brinkmann (Br) and Biot (Bi) numbers.Due to internal heat production from convective heat transfer and viscous dissipation, the microrotation profiles at plate surface decrease as Br and Bi number rise. Figure 12 & Figure 13 depict the microrotation profiles for the vortex viscosity parameter (k) for aiding and opposing flows.As the angular

Figures 14 -
Figures 14-21 show the changes in temperature profiles.The fluids temperaturereaches its peak at the plate surface before decreasing to meet the boundary requirements at free stream temperature.With an increase in the values of magnetic parameter, a drop in temperature profiles is seen.According to Figure14

Figure 15 .
Figure 15.Profiles of temperature for the suction parameter.

Figure 16 .
Figure 16.Profiles of temperature for the Prandtl number.

Figure 17 .
Figure 17.Profiles of temperature for the Brinkman number.

Figure 18 .
Figure 18.Profiles of temperature for the Biot number.

Figure 19 .
Figure 19.Profiles of temperature for the radiation parameter.

Table 1 and
et al., Ishak et al., Olanrewaju et al. and the present study are presented in Table 2.The current study is in line with earlier published efforts in the field.The Nusselt number, couple stress, and shear stress values are shown numerically in Table

Table 4 .
Shear stress and Nusselt number numerical findings.