Numerical Analysis of Micropolar Nanofluid Flow near a Stagnation Point over an Inclined Stretching Surface
Abstract:
The stagnation point flow behavior of a micropolar nanofluid over an inclined stretching surface was numerically investigated. The formulation accounts for the combined effects of Brownian motion, thermophoresis, thermal radiation, velocity slip, and the presence of internal heat generation or absorption. The governing system of non-linear partial differential equations was transformed into a set of coupled ordinary differential equations through the application of appropriate similarity transformations. These transformed equations were solved numerically to analyze the behavior of the fluid near the stagnation region, where both the stretching velocity of the surface and the external free stream velocity are assumed to vary linearly with distance from the stagnation point. Special attention was paid to the influence of dimensionless parameters on key physical quantities, including skin friction coefficient, energy transfer, and Sherwood number. It was observed that increasing the stagnation point parameter leads to a reduction in skin friction, while the inclination angle demonstrates an opposing effect on heat and mass transfer rates. Data extracted from graphical results was tabulated to provide quantitative insights into the impact of varying parameters. The findings offer significant implications for microscale heat and mass transfer systems, particularly in processes involving inclined geometries and nanoparticle-enhanced fluids under magnetohydrodynamic (MHD) effects.1. Introduction
When it comes to the study of fluid dynamics and heat transfer, the analysis of flow at a stagnation point contributes significantly to many different parts of the field. Places where the velocity of a fluid in respect to a solid surface is equal to zero are referred to as points of stagnation. This arrangement of movement is applied in a variety of domains, including chemical engineering, the research of boundary layers, and aerodynamics, among others. In order to improve the efficiency of heat transfer procedures and optimize the structure of various engineering systems, it is vital to have a solid understanding of the patterns and characteristics of flows that occur at stagnation sites. Hiemenz [1] first investigated the dynamics of fluids around the point where the fluid velocity relative to the surface becomes zero. Mahapatra and Gupta [2] investigated the features of flow in the context of viscoelastic fluid when they were inspired by the applications that were discussed earlier. Nazar et al. [3] conducted an investigation into the flow analysis in relation to a micropolar fluid. A comprehensive investigation into the phenomenon of flow at a stagnation point in connection to a plate being stretched was carried out by Chiam [4]. The movement that occurs as a result of this structure is of utmost significance in the subject of fluid dynamics. It is also utilized in a wide range of engineering fields, including heat transmission, the management of boundary layers, and the processing of materials. An analysis of the numerical modeling of the flow of micropolar fluids towards an inclined surface was carried out by Shamshuddin and Thumma [5].
The notion of micropolar fluids was first introduced in recent years as a result of the necessity to provide a description of a fluid that contains micro-elements whirling continuously. Micropolar fluids are able to link the rotating movement of atoms with the very small speed field. The components of these fluids are solid particles that are suspended in a dense material. The blood of animals, foamy fluids, magnetic liquids, and other types of liquids are all examples of micropolar liquids. There are a number of intriguing applications for these fluids, including the creation of organic compounds, slippery substances, and polymers mixed together. The theory of micropolar fluids has been thoroughly investigated by a remarkable number of researchers. At first, Eringen and Suhubi [6] introduced the idea of micropolar fluids by analyzing the applications of non-Newtonian fluids. An investigation into the movement of micropolar fluid as it approached a sloped surface was carried out by Rahman et al. [7]. In addition, Uddin [8] conducted research on the thermal interaction of a micropolar fluid with an inclined porous surface. Bhargava et al. [9] investigated the micropolar transfer processes brought about by the non-linear stretching of a sheet by utilizing numerical methods. Two distinct methods were applied simultaneously: the finite element technique and the finite difference technique. Takhar et al. [10] studied the flow of micropolar fluids over a flexible surface, which looked into the combination of MHD and convection.
As a means of improving the efficiency with which solar collectors absorb solar energy, nanomaterials are utilized in these devices. Nanoparticles have recently emerged as adaptable tools in the field of medical science. They are utilized for a wide range of applications, including the elimination of cancer cells, the diagnostic imaging of medical conditions, the combating of microorganisms, and the administration of stem cell therapies. The concept of nanofluids is not a novel one by any means. In the course of their investigation on novel coolants and cooling techniques, Choi and Eastman [11] published the first report on them. The numerous uses found in nuclear reactor systems, heat exchangers, electronic cooling, boilers, and energy storage devices were the contributing factors that led to its meteoric growth in popularity. Nanoparticles or nanofibers are extremely small particles with a diameter of less than 100 nanometers, which are present in nanofluids in minuscule quantities. Blending nanoparticles or nanofibers with a variety of heat transfer base fluids, such as water, ethylene glycol, toluene, and motor oil, results in the production of nanofluids. A nanofluid that contains microorganisms was discussed by Al-Amri and Muthtamilselvan [12] in relation to the flow that occurs at the stagnation point. Following the transformation of the system of partial differential equations into a set of ordinary differential equations, the system was reorganized. Following this, numerical methods were utilized with the support of the bvp4c MATLAB software in order to solve these equations. The flow of micropolar nanofluids near stagnation points on inclined surfaces was investigated by Rafique et al. [13]. Ram et al. [14] studied the stagnation point flow of micropolar fluid over a stretching sheet under variance parameter. Soomro et al. [15] studied the influence of non-linear radiation and heat production or absorption on the movement of nanofluid at its point of stop along a stretching sheet.
The influence of magnetism on a fluid is referred to as MHD. This influence is affected by an external magnetic force, heat exchange within the fluid, and the currents that occur within it. When the fluid comes into contact with a more powerful magnetic field as a result of more experience, it encounters resistance that impedes its flow, which causes it to slow down. The movement of liquid over a surface that is constantly expanding and is impacted by a magnetic field has a significant impact on a variety of engineering fields, including the study of plasma and the extraction of geothermal energy. Investigations on the influence that MHD effects have on the flow of fluid over a stretching surface can be found in the sources that are now available. An analytical resolution for the MHD flow that is induced by a stretched sheet was reported by Ganji et al. [16]. This resolution was presented in a non-linear form. Olkha and Dadheech [17], [18] conducted an entropy analysis in the MHD flow of a variety of non-Newtonian fluids that were caused by a stretched sheet with a slip effect and a heat source. The notion of heat and mass transfer in MHD flow was investigated by Chen [19]. This flow was driven by natural convection from a stretching surface that was permeable and correctly inclined, and it had a wall with variable temperatures and concentrations. Aydin and Kaya [20] investigated the flow of heat transfer using mixed convection in MHD around a plate that was angled at a right angle.
Since the beginning of time, the effects of thermal radiation on the movement of heat and the exchange of heat have been a source of tremendous fascination for scientists and researchers. These astonishing applications can be found in the fields of engineering, physics, and space technology, and they also play an important part in operations that involve high temperatures. These effects can also be an essential component of the factors that regulate heat transfer in the industrial sector. In this sector, the quality of the final product is dependent, to a certain extent, on the heat-regulating elements. Research conducted by Ram et al. [21] addresses a non-linear mathematical model for micropolar stagnation point flow via a permeable stretchable device. This model takes into account the influence of changing responding species as well as heat sink/source effects. Additionally, this one-of-a-kind combination sheds light on complex fluid dynamics, providing fresh insights that go beyond conventional studies in the field. In order to simplify the complex workings of micropolar fluid dynamics, Kausar et al. [22] investigated the complexities of viscous dissipation and thermal radiation. Farooq et al. [23] investigated the flow of a thick, flexible nanofluid by taking into account the effects of non-linear radiation. In a recent study, Maleki et al. [24] investigated the mobility of nanofluid by taking into account the impact of radiation and convective boundary conditions. An intricate flow of micropolar nanofluid with non-linear radiation effects was addressed by a mathematical framework that was proposed by Guedri et al. [25]. An investigation of the convective movement of a micropolar nanofluid over a vertical surface was carried out by Ahlawat and Sharma [26].
The combined effects of heat transmission, mass diffusion, and chemical reactions play an important part in a variety of processes, including the cooling of nuclear reactors, thermal insulation, and geothermal reservoirs, among others. In their study, Andersson et al. [27] investigated the phenomenon of chemically reactive species diffusing across a surface that was flat and elastic. An investigation into the movement and heat exchange of a non-Newtonian power law liquid with mass diffusion and chemical reaction around a revolving cylinder was carried out by Abo-Eldahab and Salem [28]. The influence of a magnetic field was taken into consideration in the investigation. Tripathy et al. [29] conducted an analysis of the quantitative evaluations of hydromagnetic micropolar fluids flowing over an extended surface that was placed in a porous channel. Additionally, the different heat sources and chemical reactions that were permitted were considered.
Several processes, including cooling nuclear reactors, thermal insulation, and geothermal reservoirs, are dependent on the combined effects of heat transmission, mass diffusion, and chemical reactions. These effects play an important role in the overall process. Yao [30] studied the analytical solution technique for diffusion of chemically reaction species over a stretching sheet. Hypothesized by Patil et al. [31] that the Eyring-Powell fluid has considerable effects when it flows over a stretched surface in the presence of a magnetic field and chemical reactions. An investigation into this kind of flow was carried out by Nadeem et al. [32] in the context of a viscous nanofluid flowing across a curved surface. The objective of the study was to analyze the behavior of the nanofluid when it comes into contact with a curved surface while being subject to the influence of a magnetic field and viscous effects.
On the other hand, many different devices are used for manipulating fluids, and one of the most important factors in these devices is the temperature difference between the surface and the fluid layer that is close to it. The fluctuation in temperature is the factor that is responsible for either the generation or absorption of heat, which has a significant impact on the parameters of heat exchange. This is especially striking in procedures in which the fluid being utilized is subject to chemical reactions that either release or absorb heat, as well as in the production of metal waste as a secondary result of the utilization of nuclear energy. The exploration of heat generation or absorption has become a primary topic for academics who are interested in examining difficulties related to fluid dynamics, which has sparked a great deal of interest among experts in the field. Vasanthakumari and Pondy [33] investigated the effect that the creation or absorption of heat has on the movement of nanofluid and its direction towards an inclined surface. An investigation into the change in heat transmission brought about by a micropolar nanofluid in the presence of a porous medium was carried out by Moatimid et al. [34]. An examination of the finite element associated with heat-mass transfer in hydromagnetic micropolar flow across a stretched surface was carried out by Kumar [35]. Emad et al. [36] investigated the effects of flowing or suction forces on the transport of heat through hydromagnetic systems. Specifically, the mixed convection on a surface that is continuously extending was studied by taking into consideration the creation or absorption of heat from within the surface.
One of the most distinctive factors that contributes to the reduction of resistance in a variety of technical applications, such as aerodynamics and marine vessels, is the participation of slip. Because slip can have an effect on the passage of fluids through porous rocks, slip effects are widely used in enhanced oil recovery and petroleum engineering to help improve oil recovery. The research conducted by Navier is connected to the concept of slip in the field of fluid dynamics. The analysis of fluid movement in viscous substances was the subject of a study by Fang et al. [37]. Ibrahim [38] investigated the comprehension of the slide phenomenon in non-Newtonian fluids and conducted a related analysis. The application of porous dissipation is what sets this work apart from others and provides a fresh perspective on the way fluid dynamics and thermal effects interact with one another. Providing experts and researchers [39], [40], [41] with useful views is one of the many ways in which these studies make significant contributions to their respective fields.
2. Mathematical Formulation
Within the context of a micropolar nanofluid, the primary purpose of this investigation is to investigate the influence that radiation and numerous slip effects have on the flow of MHD. This study uses the assumption that there is a magnetic field that is perpendicular to the inclined surface and has a strength of $B_0$. Furthermore, this study overlooks any magnetic field that is induced by the surface. The $x$‐axis is aligned along the inclined surface, while the $y$‐ axis is normal to the surface, as depicted in Figure 1. Stretching was accomplished by moving the wall with a velocity $U_w=a x$ (with $a>0$) along the -axis. Additionally, $\Omega$ is the angle made relative to the vertical axis of the sheet during the stretching process. The velocity was taken as $U_{\infty}=b x$ as $y \rightarrow \infty$. The temperature near and distant from the surface are denoted as $T_w$ and $T_{\infty}$, respectively, with corresponding concentrations of $C_w$ and $C_{\infty}$, respectively. Thermophoretic effects and Brownian motion were taken into account.

The governing equations are as follows [42]:
The Rosseland approximation is given by:
Using Taylor’s expansion, $T^4$ can be expressed as:
Combining Eq. (4) with Eqs. (6) and (7) leads to:
Boundary conditions are as follows [2]:
For the study at hand, the stream function $\psi=\psi(x, y)$ has the following form:
The definition of the similarity transformations is as follows:
Eqs. (2) to (5) in their transformed form can be obtained by applying Eq. (8) as follows:
where,
$\begin{aligned} & K=\frac{K_1^*}{\rho_f} M=\frac{\sigma_f B_0^2}{\rho_f a}, G r_x=\frac{g \beta_T\left(T_w-T_{\infty}\right) x^3}{\upsilon^2}, G r=\frac{G r_x}{\operatorname{Re}_x^2}, G c_x=\frac{g \beta_C\left(C_w-C_{\infty}\right) x^3}{\upsilon^2}, G c=\frac{G c_x}{\operatorname{Re}_x^2}, \\ & \operatorname{Re}_x=\frac{U_w x}{\upsilon}=\frac{a x^2}{\upsilon}, A=\frac{b}{a}, \operatorname{Pr}=\frac{\upsilon}{\alpha}, N b=\frac{\tau D_B\left(C_w-C_{\infty}\right)}{v}, N t=\frac{\tau D_T\left(T_w-T_{\infty}\right)}{\upsilon T_{\infty}}, R=\frac{4 \sigma^* T_{\infty}^3}{k^* k}, \\ & Q=\frac{Q_0}{a(\rho c)_f}, S c=\frac{\upsilon}{D_B}, \delta_1=\delta_1^* \sqrt{\frac{a}{\upsilon}}, \delta_2=\delta_2^* \sqrt{\frac{a}{\upsilon}}, \delta_3=\delta_3^* \sqrt{\frac{a}{\upsilon}}, S=-\frac{V_w}{\sqrt{a \upsilon}} \end{aligned}$
Conditions at those boundaries are as follows:
The relevant physical quantities—namely the skin friction coefficient, local Nusselt number, and Sherwood number—are defined as follows:
where, $\tau_w$ is the stress with the stretched surface, $q_w$ is the wall heat flux, $q_m$ is the mass heat flux, and $k$ is the thermal conductivity, which are stated as follows:
Using variables with no dimensions, the following can be obtained from Eq. (14):
3. Methodology for Solution
The boundary layer flow of micropolar nanofluids over an inclined stretched surface was investigated by taking into account the influence of thermal radiation and magnetic fields. The influences of the Soret effect and chemical reactions were also examined. For the purpose of determining the thermal efficiency of a fluid flow in the presence of Brownian motion and thermophoresis, the model developed by Buongiorno was used. Mathematica NDSolve was used to perform numerical resolution on the differential equations that are supported by the similarity transformation-converted differential equations. To validate the accuracy of the numerical scheme, the latest results of $-\theta^{\prime}(0)$ and $-\phi^{\prime}(0)$ were assessed against the outcomes of Khan and Pop [43] and Rafique et al. [42], as shown in Table 1.
$\boldsymbol{N b}$ | $\boldsymbol{N t}$ | Rafique et al. [42] | Khan and Pop [43] | Outcomes of This Study | |||
$-\theta^{\prime}(0)$ | $-\phi^{\prime}(0)$ | $-\theta^{\prime}(0)$ | $-\phi^{\prime}(0)$ | $-\theta^{\prime}(0)$ | $-\phi^{\prime}(0)$ | ||
0.1 | 0.1 | 0.9524 | 2.1294 | 0.9524 | 2.1294 | 0.952234 | 2.129302 |
0.2 | 0.2 | 0.3654 | 2.5152 | 0.3654 | 2.5152 | 0.365238 | 2.515209 |
0.3 | 0.3 | 0.1355 | 2.6088 | 0.1355 | 2.6088 | 0.135876 | 2.608783 |
0.4 | 0.4 | 0.0495 | 2.6038 | 0.0495 | 2.6038 | 0.049564 | 2.603893 |
0.5 | 0.5 | 0.0179 | 2.5731 | 0.0179 | 2.5731 | 0.0170987 | 2.572190 |
4. Results and Discussion
While most parameter values are varied across a range during the simulation, certain parameters are held constant throughout the process such as $ K=0.5$, $S=A=\delta_1=\delta_2=\delta_3=\delta_4=Q=G r=G c=M=N b=N t=0.1$, $R=0.5$, $L e=5.0$, $\Omega=60^{\circ}$, $\operatorname{Pr}=6.5$.
Figures 2, Figure 3, Figure 4, Figure 5 illustrate the effects of the material parameter ($K$) on the velocity $f^{\prime}(\zeta)$, microrotation $h(\zeta)$, thermal $\theta(\zeta)$ contour and concentration profiles. Whenever the values of $K$ increase, the momentum and thermal contours are improved, whereas the microrotation contour is reduced. Physically, in micropolar fluids, the material parameter that can affect the momentum contour is recognized as the micropolar fluidity parameter ($K$). When the micropolar fluidity parameter ($K$) increases, it implies that the microstructure or internal degrees of freedom have a stronger effect on the fluid movement. This can lead to an increase in the complexity of the flow patterns and the momentum contour. Figure 5 exemplifies the consequences of the material parameter ($K$) on the concentration profile. As the values of K increase, the concentration contours show a corresponding enhancement. Figure 6 displays the effect of the magnetic field parameter ($M$) on momentum contour. The values of $M$ increase, resulting in decay in the momentum contour $f^{\prime}(\zeta)$. Lorentz force came into existence when a magnetic field was imposed on the flow field. This force acts to resist the fluid motion by reducing its velocity. Hence, both the fluid flow velocity and the thickness of the momentum boundary layer decrease.





Figure 7 illustrates that as the value of $M$ increases, the microrotation, denoted as $h(\zeta)$, decays in proximity to the surface and rises as the distance from the surface upsurges. In other words, the microrotation tends to reduce close to the surface while it tends to rise farther from the surface when $M$ is enlarged. The effects of the magnetic field parameter on the flow temperature and concentration are shown in Figure 8, Figure 9. From the figures, it is observed that the magnetic field increases the temperature and concentration contours. This can be attributed to the Lorentz force—a fractional resistive force which opposes the fluid motion and produces heat. Therefore, the stronger the magnetic field, the thicker the temperature and concentration boundary layers.



The development of the non-dimensional velocity contour with varying $A$ values is presented in Figure 10. Within the edge film, the fluid velocity flows as it moves farther away from the wall and decays when it gets nearer to it. A smaller $A$ value means that a wall velocity is bigger than the free stream velocity. The fluid’s velocity initially declines as it approaches the wall because of the wall inertia, which interacts with the surrounding fluid in the edge film. The fluid does, however, experience a local acceleration at an exact distance from the wall, which may make the local velocity drop under the free stream velocity. The slower-moving fluid is then entrained by the free stream velocity, raising the velocity. As shown in Figure 10, the fluid does not feel the free stream’s dragging effect when $A$ equals 0. In fact, the free stream velocity is nearly zero at this point. For all values of $A > 0$, the boring effect of the free stream rushes the local velocity external to the boundary layer. Figure 11 depicts the difference of the thermal contour in a change in the values of the momentum ratio parameter $A$. It displays that as the velocity ratio parameter increases, the thermal boundary layer’s thickness decreases. Additionally, the surface temperature (in absolute value) increases as $A$ increases. As a result, the thermal contour declines. Figure 12 portrays the influence of velocity ratio parameter $A$ on the concentration graph. As the values of $A$ increase, the concentration boundary film’s thickness decreases. Moreover, it is known from the graph that the magnitude of the surface concentration upsurges as $A$ increases.



Figure 13 reveals the decreasing presentation of the momentum for increased values of the angle of inclination $\Omega$. The fluid flow becomes more resistant due to the Lorentz force, which lowers the fluid’s velocity. The effects of the velocity slip parameter on the momentum profiles are depicted in Figure 14. From the graph, it can be seen that as the value of the momentum slip upsurges, the velocity profiles decline. This is due to the fact that as the momentum slip increases, the velocity of the fluid decays because the pulling of the stretching sheet can transmit the fluid. The same trend can also be viewed in the case of the relationship between the thermal slip factor and the thermal profiles in Figure 15. The thermal curves for different values of the thermal radiation parameter are portrayed in Figure 16. From the graph, it is possible to observe that as the values of the thermal radiation parameter increase, both the thermal graph and the thermal boundary layer’s thickness increase. Figure 17 is included to examine the effect of the heat generation/absorption parameter on the thermal contours. From the figure, it can be seen that the thermal contours upsurge as the heat generation/absorption parameter increases because the thermal state of the fluid rises with an increase in the heat generation.





The effects of the Brownian motion parameter on the solutal contours are portrayed in Figure 18. It can be seen that as the values of the Brownian motion parameter increase, the solutal contours’ thickness decreases. Figure 19 portrays the impact of the thermophoresis parameter on the thermal contour. From the figures, it is perceived that when the thermophoresis parameter increases, there is an improvement in the thermal boundary layer’s thickness. The solutal contours are presented in Figure 20 for various values of the thermal slip parameter $\delta_3$. It is observed that as the values of the thermal slip parameter increase, the solutal contours decline. Figure 21 shows the variation of the skin-friction coefficient with respect to $M$ and $K$. It’s observed that as $M$ and $K$ increase, the skin-friction coefficient upsurges. In Figure 22, it is observed that the Nusselt number declines when both $K$ and $M$ increase. Figure 23 displays that the local Sherwood number decreases with increasing values of $M$ and $K$.






5. Conclusions
This study discusses the heat and mass exchange that occurs in the micropolar nanofluid flow adjacent to the stagnation point. This flow is characterized by magnetic effects, slip effects, radiation, and chemical reactions, as well as the generation or absorption of heat. The flow of the micropolar nanofluid, which is generated by the linear stretching inclined sheet, is a crucial component in the heat transfer procedures utilized in industrial and cooling systems. This study addresses the research gap related to this phenomenon.
The following are the key findings of this research:
$\bullet$ Larger Lorentz forces ($\digamma$) and porosity indices are beneficial to microrotation because they increase axial traction and fluid-medium interaction among the particles.
$\bullet$ Thermophoresis, magnetic effects, and Brownian motion are all factors that contribute to the heating of micropolar nanofluids, which ultimately results in an improvement in the thermal performance of these fluids.
$\bullet$ The utilization of enhanced Brownian motion and increased Lewis numbers has the potential to reduce concentration swings. This is accomplished by more effectively dispersing nanoparticles and reducing the amount of mass diffusion.
$\bullet$ When it comes to optimizing pharmaceutical components, microorganism systems, and manufacturing techniques, it is particularly important to understand the skin friction that occurs in micropolar nanofluids. With this newfound knowledge, it is possible to lessen the amount of drag, enhance the efficiency of heat transfer, and cut down on energy use.
$\bullet$ The findings of this study have the potential to be applied in the field of building envelope applications due to the heat transfer.
The future scope of this work includes the investigation of novel nanofluid compositions with the purpose of improving heat and mass transport, the construction of more accurate models to account for dynamic phenomena, and the verification of the results through experimental investigation. The incorporation of these ideas has the potential to improve efficiency, dependability, and energy conservation in a variety of domains, including industrial processes such as cooling mechanisms and pharmaceuticals.
The data used to support the findings of this study are available from the corresponding author upon request.
The authors declare that they have no conflicts of interest.
