MHD boundary layer micropolar fluid flow over a stretching wedge surface Thermophoresis and brownian
Journal of Thermal Engineering 2024, Vol. 10, Issue 2, pp. 330-349; doi.org/10.18186/thermal.1448609
Abstract
Keywords: Chemical Reaction; Micropolar Fluid; Micro-rotation; MHD; Radiation; Shooting Method; Stretching Wedge Surface
Introduction
The micropolar fluid attracted a great deal of interest from the researcher due to its importance in various fluid fields. The micropolar fluid is the sub-class of micro fluids that exhibits the effects and inertia of micro-rotation. Couple stress and body couples are supported by the micropolar fluids. The bar-like elements containing fluids can effectively symbolize the micropolar fluids. Micropolar fluid is a type of anisotropic fluid like liquid crystals composed of dumbbell molecules. Animal blood is also this type of fluid. Eringen [1] has introduced a fluid model which demonstrates several microscopic effects that arise from the micro-motion as well as the local formation of fluid constituents. The fluids which are consisting of inflexible, indiscriminately oriented, or spherical species suspended in the viscous medium are symbolized by the micropolar fluid model and have ignored the particles’ deformation. “Micropolar fluids have been shown to accurately simulate the flow characteristics of polymeric additives, geo-morphological sediments, colloidal suspensions, hematological suspensions, liquid crystals, lubricants, etc.”[2]. Various authors have investigated heat and mass transfer boundary layer flow considering micropolar fluid in different physical situations. Srinivasacharya et al. [3] studied the boundary layer free convective flow considering variable wall temperature as well as a concentration in stratified micropolar fluid on a vertical plate. It was noticed that when stratification is present, the sign of microrotation changes within the boundary layer region from negative to positive values. Khalid et al. [4] have investigated natural convection heat and mass transfer boundary layer flow in a micropolar fluid. For this, they have considered the unsteady flow over an infinite oscillating vertical plate and found that the velocity of Newtonian fluid is larger than the velocity of the micropolar fluid. Many researchers have a great interest in the study of nanofluids because of their several applications in different fields like automotive, electronics, computing technology, communications, optical devices, material processing, lasers, medicine, and material amalgamation, etc., while heat dissipation is required. In 1995, Choi and Eastman [5] first invented the nanofluids notion and they have observed high thermal conductivities demonstrated by nanofluids owing to the existence of solid particles which are nanometer (nm) sized. Some authors [6 -12] have studied mixed convective heat and mass transfer considering different types of nanofluids with diverse conditions as well as their various applications such as Jilte et al.[13] have tested heat sink cooling concert and focused on the geometry of four channels liquid-cooled heat sink configured in an alternate flow passage concentric shape as well as found that the sink of nanofluid cooled heat compared with water rejected a higher heat rate, on the other hand, Jilte et al. [14] have presented two systems of organizing cooling component which is liquid circulated battery cooling systems (LcBS) and liquid-filled battery cooling systems
(LfBS) containing water-based nanofluid (Alumina 40nm diameter nanoparticles) and the results showed that the production and implementation of Li-ion batteries are easier in electric vehicles for the modular nature of both battery cooling system and LfBS uncomplicated in operation. While the turbulent effect is ignorable, the thermophoresis, as well as Brownian diffusion, play an important role in the fluid flow and Buongiorno [15] has recommended a model considering the consequence of thermophoresis and Brownian diffusion and found that an essential decline of viscosity for a heated fluid in the boundary layer which enhanced heat transfer. A few researchers [16-18] investigated the convective boundary layer flow in diverse fields employing the Buongiorno model. Due to its use in a variety of sectors, including MHD power generators, plasma studies, nuclear reactor cooling, the petroleum industry, crystal development, and boundary layer control in aerodynamics, magnetohydrodynamics (MHD) plays a significant role in the study of free convection flow. Pavlov [19] has investigated the MHD boundary layer flow over a stretching surface. In his study, he considered a uniform magnetic field and obtained exact similarity solutions. The viscoelastic fluid flow past a stretching surface has been investigated considering a transverse magnetic field [20]. In this regard, he has obtained a set of the exact solution and found that a transverse magnetic field has a similar effect as the viscoelasticity on the fluid flow. On the other hand, Kasiviswanathan and Gandhi [21] have determined a set of the exact solution for the flow of micropolar fluid between infinite non-coaxial revolving disks in the presence of magnetic field effect and observed that the velocity profile is deformed as well as the magnitude of micro-rotation components vector is diminished for the enhance in Hartmann number. Chamkha et al. [22] have investigated the MHD dominated Cu-water nanofluid flow in a porous enclosure and found that for all values of heat generation parameter inside the porous cavity the entropy generation is decreasing for increasing the nanoparticles volume fraction and miniature volume fraction of nanofluid has a significant consequence on the thermal rendering. The Falkner-Skan equation (F-SE) plays a substantial role in the augmentation of fluid dynamics. “The Falkner– Skan equation, a third-order ordinary differential equation over a semi-infinite domain was originally developed in 1931 by Falkner and Skan [23] from a similarity analysis on the steady, two-dimensional, boundary-layer flow of an incompressible, viscous fluid over a wedge [24]”. Several researchers [25-28] have reviewed flow over a wedge. The flow with a triangular shape is called a wedge flow. In the procedure under consideration, it can be utilized to pull up one object while keeping another in a plane. The wedge transforms the lateral force into a crosswise splitting force. Later, several authors [29-31] investigated the analytical and numerical solutions of the flow past a wedge considering different flow properties. In high-temperature operating systems radiation has a great impact on the heat and mass
transfer of the fluid flows. Many processes are carried out at a high temperature in the fluid of engineering. In such circumstances, the investigation of radiative heat transfer plays a vital role in the choice of suitable equipment. There aren’t many industries that produce gas turbines, radioactive and nuclear power plants, manmade satellites, aircraft, weapons, wind turbines, and other things. Dogonchi et al. [32] have investigated the consequence of radiation and heat source on the heat transfer of nanofluid flow and found significant effects of nanofluid volume fraction and radiation parameter on the Nusselt number. Zheng et al. [33] have found a significant consequence of radiation parameters on the heat transfer of nanofluid flow in their investigation. Ahmed et al.[34] have investigated MHD mixed convection water-based micropolar nanofluids flow in trapezoidal areas numerically and found that average Nusselt number enhances for decreasing values of Richardson number as well as dimensionless viscosity, heat source location, Hartmann number have an essential consequence on the average Nusselt number. Bagh et al.[35] have studied the g-jitter effects on mixed convection for the two colloidal liquids flow which are the alone micropolar nanofluid such as copper/water and micropolar hybrid nanofluid such as alumina-copper/water and observed the distinguished increased performance of hybrid nanofluids in the upper space of inclined sheet g-jitter influenced flow. They also found that the variation of skin friction coefficient and heat flux gradient enhanced with larger values of inflection amplitude. Kasmani et al.[36] investigated the double-diffusive convective boundary layer nanofluid flow past a moving wedge with the effects of Soret and Dufour in the presence of suction and found that the temperature enhances with increasing effect of the Soret parameter as well as declines for the increasing effect of Dufour parameter and the opposite trend arises in case of the mass transfer rate. Ashraf et al.[37] have studied the steady MHD convective tangent hyperbolic nanofluid flows across an elastic
surface with the impact of thermal radiation and thermal conductivity. And they observed that for the increasing effects of the Weissengberg number and power-law index, the velocity profile decreased but the temperature profile increased. Ahmed et al. [38] have identified the consequences of thermal conductivity as well as dynamic viscosity on the boundary layer flow of nanofluid and characteristics of heat transfer over a stretching permeable tube with a heat source/sink. The skin friction coefficient decreases, on the other hand, the local Nusselt number increases for the effects of the Reynolds number and suction/injection parameter. The new approach of the present study contributes significantly to the understanding of numerical solutions to nonlinear differential equations in fluid mechanics and micropolar fluid flow. Micropolar fluids are becoming even more of a focus due to the desire for engineering applications in various fields of medical, mechanical engineering, and chemical processing. This study looks into how thermophoresis and Brownian diffusion affect convective heat and mass transfer boundary layer micropolar fluid flow over a stretching wedge-shaped surface in the presence of thermal radiation and chemical reaction.
Mathematical Model With FLOW Configuration
The two-dimensional steady convective boundary layer flow of viscous, incompressible, and electrically conducting fluid and heat transfer of micropolar fluid over a stretching wedge-shaped surface considering chemical reaction and radiation effects are demonstrated in this study. Coordinate system and configuration of flow as shown in Figure 1. The Cartesian coordinate x- and y-axes are determined along and normal to the surface of the wedge, respectively. Fluid velocities such as u and v act along x- and y- axes, correspondingly. Uw(x) = axm, where a > 0, is the stretching
wedge velocity and U(x) = bxm is the free stream velocity outside the boundary layer of the micropolar fluid. In Figure 1, Ω = πβ is the wedge angle where β is the Hartree pressure gradient. Under the above assumptions, the basic equations are considered as [39]:
where is the Stefan-Boltzmann constant as well as k* the mean absorption coefficient, and T 4 stand for the linear sum of the temperature. Applying Taylor series expansion T 4 can be expanded T∞ as:
The temperature differences are considered very small, so ignoring the higher order terms involving (T − Too , Eq. (7) one can be written as:
To transfer the set of PDEs (Eqs. (2), (3), (5), and (10)) into a set of dimensionless ODEs introduce the following similarity transformations:
1. At the wedge surface (i.e. at y = 0), N = 0, while s = 0.
It represents the flow of concentrated particles in which the microelements closed to the wall surface are unable to rotate. This case is also known as a strong concentration of microelements. The antisymmetric part of the stress tensor will vanish when s = 0.5 consequently concentration of microelements becomes weak. Whereas, the case s = 1 is used for the modeling of turbulent boundary layer flows [41]”. Considering Rossel and approximation the radiative heat flux can be written as (6)
where Cw is the concentration at the stretching surface and Tw is convective fluid temperature. Using Eq. (11) in Eqs. (2), (3), (5), and (10), one can be determined the following ODEs: (12)
for the diverse values of concerned physical parameters employing the Runge-Kutta fourth-order method by exploiting symbolic software MATLAB. In this regard Eqs. (12) - (15) along with boundary conditions are converted into a set of first-order ODEs and then are situated into the MATLAB numerical software code to find the solutions. Let us consider into account
In the case of A > 0 and A < 0 the wedge moves in the same and opposite directions, respectively. On the other hand, the wedge will be static while A = 0 [42]. The skin friction coefficient, Nusselt number, Sherwood number, and couple stress coefficient are the physical quantities of engineering interest. One can be defined as:
(21) The numerical results of velocity (f '), temperature (θ), and nanoparticle volume fraction (ϕ), angular velocity (h) profiles are calculated for diverse values of obtained physical parameters and presented graphically, as well as discussed physical viewpoint.
Results And Discussion
, and indicate the shear stress, the transfer rate of heat, and mass per unit area of the surface, respectively. The three quantities are defined as (16) With the aid of the above equation (Eq. (16)), one can determine the one-dimensional form of Cf, Nu and Sh, Cr which are proportional to and respectively. Solution Procedure It is impossible to find the exact solution of higher-order nonlinear coupled ordinary differential equations (CODEs). Therefore the numerical solution of the CODEs Eqs. (12) - (15) with boundary conditions are obtained
The Investigations are conducted on the impacts of the relevant physical parameters on the velocity, temperature, and concentration profiles. Equations (11) through (14) are dimensionless coupled nonlinear ordinary differential equations having boundary conditions (Eq.21) that are numerically solved using the fourth-order Runge-Kutta (RK) technique and the shooting approach. This results in the conversion of Eqs. (11) through (14) into nine first-order ordinary differential equations including boundary conditions, described in Section 3. The typical values of the relevant physical parameters (especially the Prandtl number and Grashof number are taken for the water and the cool surfaces, respectively) are taken into consideration for numerical solutions as coupling constant parameter (0.01 ≤ B1 ≤ 0.05), magnetic parameter (1.0 ≤ M ≤ 15.0), Grashof number (0.3 ≤ Gr ≤ 0.9), modified Grashof number (0.3 ≤ Gm ≤ 0.8), micropolar parameter (2.0 ≤ G2 ≤ 7.5) vortex viscosity constraint (0.02 ≤ G1 ≤ 0.2), Prandtl number (7.0 ≤ Pr ≤ 15.0), thermal radiation parameter (0.25 ≤ R ≤ 0.50), Brownian
motion parameter (0.2 ≤ Nb ≤ 0.62), thermophoresis parameter (0.04 ≤ Nt ≤ 0.10), heat generation parameter (0.1 ≤ Q ≤ 0.5), Biot number (0.65 ≤ Bi ≤ 1.0), stretching parameter (0.2 ≤ A ≤ 0.5), Lewis number (3.0 ≤ Le ≤ 7.0), and chemical reaction parameter (0.2 ≤ K ≤ 0.7). In this section, the impacts of the relevant parameters on velocity, temperature, concentration, and angular velocity are discussed in detail. Effect of the Magnetic Parameter (M) Figure 2 depicts the velocity (f'), temperature (θ), concentration (ϕ), and angular velocity (h) profiles for the effect of the magnetic parameter (M). For higher values of the magnetic parameter, a substantial resistive force called the Lorentz
force opposes the fluid motion, and as a result, the rate of heat transfer within the boundary layer increases. It indicates that one of the sources of fluid ruggedness and heat transport is symbolized by the magnetic field. Once higher values of M the velocity (concentration) profile enhance in the range 0 < η <
2.5. (0 < η < 1) depict in Figures 2a and c, respectively, while the
temperature (concentration) profile decreases in the range 0 ≤ η < 2 (1 ≤ η< 2) for enhancement of M elucidating Figures 2b and c, respectively. Figures 2a and b also show that, due to the influence of the magnetic parameter, the momentum boundary layer is growing and the thermal boundary layer is shrinking. On the other hand, weak Lorentz force (excited for the effect of magnetic field) guides to deprecation in temperature
Figure 2. Influence of M on (a) f '(η), (b) θ(η), (c) ϕ(η), and (d) h(η) profiles, taking remaining parameters as m = 0.1, Gr = 0.3, Gm = 0.3 , β = 0.1, Le = 5, R = 0.4, Pr = 10, Nt = 0.1, Nb = 0.2, Q = 0.5, K = 0.5, Bi = 1.0, A = 0.2, B1 = 0.01, G1 = 0.1, G2 = 2.0, s = 0.1.
and development in velocity profile with prosperous values of the magnetic parameter. Figure 2d elucidates the angular velocity profile (h) for the increasing effect of M. It is seen from the figure that the h is decreasing for increasing values of M in the limit 0 ≤ η < 4 and then static. Effect of Lewis Number (Le) Figure 3a exhibits the effect of Lewis number (Le) on the velocity profile (f ') and found that Le has an inconsequential increasing effect of f '. With increasing values of Le the temperature, the profile is increasing in the range 0 ≤ η < 1.8 and then static, which depicts in Figure 3b. Figure 3c
reveals the concentration profile (ϕ) for the effect of Le and it is observed that ϕ is decreasing near the wedge surface (0 ≤ η < 1), but a few distances from the wedge’s surface, it is rising in the range 0 ≤ η < 1.8 before becoming static due to the improvement of Le. The Lewis number (Le) is the ratio of thermal and mass diffusivities, so while the values of Le enhance, then the rate of thermal diffusion goes ahead to the rate of mass diffusion, and in the fullness of time, the concentration profile declines illustrates in Figure 3c. The effect of Le on the angular velocity profile predicts in Figure 3d and it is seen from the figure that the angular velocity
Figure 3. Influence of Le on (a) f '(η), (b) θ(η), (c) ϕ(η), and (d) h(η) profiles, taking remaining parameters as M = 1.0, m = 0.1, Gr = 0.3, Gm = 0.3, β = 0.1, R = 0.4, Pr = 10, Nt = 0.1, Nb = 0.2, Q = 0.5, K = 0.5, Bi = 1.0, A = 0.2, B1 = 0.01, G1 = 0.1, G2 = 2.0, s = 0.1.
decreases in the range 1 ≤ η < 2.4 and then there is no effect of Le on the angular velocity profile. Effect of Prandtl Number (Pr) The impact of the Prandtl number (Pr) on the velocity, temperature, concentration, and angular velocity profiles is depicted in Figures 4a, b, 4c, and d, respectively. Figure 4a demonstrates that the Prandtl number has a slight rising effect on the velocity profile (f') at a small distance from the wedge surface. On the other hand, as the Pr increases, the temperature profile rises close to the wedge surface in the range 0 ≤ η < 1.5, and there is no longer an effect, as seen in Figure 4b. This phenomenon physically ascribes as, due
to the larger value of Pr increases the viscous forces of the micropolar fluid with large values of η, consequently the temperature profile increases in the vicinity of the wedge surface. The variation of the concentration profile detects in Figure 4c for the different values Pr and the result indicates that the concentration profile declines close to the wedge surface (in the range 0 ≤ η < 0.9 ) while it is increasing far from the wedge surface (1 ≤ η < 2) and then static with enhancesof Pr. These phenomena explain as, for the raising values of the Prandtl number the nanoparticles accumulated into the wedge surface, as a result, the nanoparticle volume fraction profile decrease at the wedge surface and increase far from the wedge surface. The effect of Pr on the angular velocity
Figure 4. Influence of Pr on (a) f '(η), (b) θ(η), (c) ϕ(η), and (d) h(η) profiles, taking remaining parameters as M = 1.0, m = 0.1, Gr = 0.3, Gm = 0.3, β = 0.1, Le = 5, R = 0.4, Nt = 0.1, Nb = 0.2, Q = 0.5, K = 0.5, Bi = 1.0, A = 0.2, B1 = 0.01, G1 = 0.1, G2 = 2.0, s = 0.1.
profile predicts Figure 4d and the graphic illustrates how the angular velocity declines in the range 0 ≤ η < 2.2 and then has no further impact on the angular velocity profile. Effect of Brownian Motion Parameter (Nb) Figure 5a reveals the effects of the Brownian motion parameter (Nb) on the velocity profile (f ') and it is seen that the velocity profile is found to decrease for increasing levels of Nb and progressively increase over η. On the other hand, it is observed from Figure 5b that in the vicinity of the wedge surface, the temperature profile is decreasing for the effect of Nb. It is also observed that for small values of Nb the increasing rate of temperature is perceptible. The
temperature profile of the fluid and the related thermal boundary layer decline because smaller nanoparticles maintain a correlation with higher amounts of Nb. For the effect of Nb, the concentration profile increases in the vicinity of the wedge surface, while it is decreasing far from the wedge surface which is observed in Figure 5c. The motion of the nanoparticles is raised by using higher values of Nb, which causes the concentration profile to grow close to the wedge surface while decreasing farther away. Figure 5d divulges the angular velocity profiles (h) for the increasing effect of Nb. It is seen that the h is increasing for increasing values
Figure 5. Influence of Nb on (a) f '(η), (b) θ(η), (c) ϕ(η), and (d) h(η) profiles, taking remaining parameters as M = 1.0, m = 0.1, Gr = 0.3, Gm = 0.3, β = 0.1, Le = 5, R = 0.4, Pr = 10, Nt = 0.1, Q = 0.5, K = 0.5, Bi = 1.0, A = 0.2, B1 = 0.01, G1 = 0.1, G2 = 2.0, s = 0.1.
of Nb in the limit 0 ≤ η < 2 and then there is no effect and finally, it will be static for η > 4. Effect of Thermophoresis Parameter (Nt) Figure 6 demonstrates the velocity, temperature, concentration, and angular velocity profiles for the effect of the thermophoresis parameter, Nt. As can be observed from Figure 6a, the velocity profile is slightly increased by the thermophoresis parameter. The flow steams warm up by thermophoresis, aggravating nanoparticle deprivation in the process. This emphasizes the fluid velocity and speeds up the motion of the nanoparticles. When it comes to the analysis of the temperature and concentration profiles for
nanofluid flow, the thermophoresis parameter is essential. For the consequences of Nt the temperature and concentration profiles are presented in Figures 6b and c, respectively. For the increasing effect of Nt, the temperature profile increases in the vicinity of the surface of the wedge, and away from the wedge surface, the temperature profile becomes static as η → ∞ depicted in Figure 6b. This means that for the effect of the thermophoresis parameter the thermal boundary layer is thicker in the surrounding area of the wedge surface and away from the wedge surface the boundary layer thickness is consistent. Therefore, because of the collision between elevated thermal energy and lowgrade energy particles the heat transfer rate increases. From
Figure 6. Influence of Nt on (a) f '(η), (b) θ(η), (c) ϕ(η), and (d) h(η) profiles, taking the remaining parameter as M = 1.0, m = 0.1, Gr = 0.3, Gm = 0.3, β = 0.1, Le = 5, R = 0.4, Pr = 10, Nb = 0.2, Q = 0.5, K = 0.5, Bi = 1.0, A = 0.2, B1 = 0.01, G1 = 0.1, G2 = 2.0, s = 0.1.
Figure 6c it is found that in the case of the concentration profile, there is a point of intersection and the nanoparticle volume fraction profile is decreasing before the point, the concentration profile increasing after the point for the increasing effect of Nt. Figure 6d elucidates the effect of the thermophoresis parameter, Nt on the angular velocity (h) and found that in the range 0 ≤ η < 2.5 the angular velocity decreases with an increase of Nt and there is no effect in the range η > 2.5, finally, it becomes static in the range η > 4. Effect of Biot Number (Bi) Figure 7a depicts the velocity profile due to the increasing effect of Biot number (Bi) and observed that the velocity profile, as well as momentum boundary thickness,
are increasing. The effect of Bi in the temperature and nanoparticle volume fraction profiles elucidates Figure 7b and c, respectively. From Figure 7b it is seen that the temperature profile is increasing in the range 0 < η < 1.75 with the enhancement of Bi, besides from Figure 7c it is observed that the concentration profile is decreasing in the interval 0 < η < 1 and then increasing in the range 1 ≤ η <
1.9. for the increasing effect of Bi. In the range 0 < η < 1 , the
concentration boundary layer thickness decline, while that of enhanced for the effect of Bi in the range 1 ≤ η < 1.9 . It is clear that as the Biot number rises, the buoyant force also rises, which causes the fluid’s thermal energy to rise more quickly. The increase in fluid density and temperature, as
Figure 7. Influence of Bi on (a) f '(η), (b) θ(η), (c) ϕ(η), and (d) h(η) profiles, taking remaining parameters as M = 1.0, m = 0.1, Gr = 0.3, Gm = 0.3, β = 0.1, Le = 5, R = 0.4, Pr = 10, Nt = 0.1, Nb = 0.2, Q = 0.5, K = 0.5, A = 0.2, B1 = 0.01, G1 = 0.1, G2 = 2.0, s = 0.1.
well as the enhancement of the fluid’s flow by the buoyancy force, are necessary for the thermal energy’s extensibility. The Biot number is a measure of how the body’s internal conductive resistance compares to its surface-level convective resistance. Therefore, as the Biot number increases, the surface’s heat resistance decreases and its temperature rises as a result of enhanced convection. Figure 7d describes the angular velocity (h) for the increasing effect of Bi and found that the h is decreasing in the range 0 ≤ η < 2.2 , while in the range 2.2 ≤ η < 4.2 there is no effect of Bi on h, and then the angular velocity becomes static.
Effect of Chemical Reaction Parameter (K) Figures 8a and b display the velocity and temperature profiles for the increasing effect of chemical reaction parameter (K) and found that in the vicinity of the surface, the velocity profile decreases in the range 0 < η < 2.7 and then becomes static in the limit η ≥ 2.7, and for the effect of K the temperature profile also decreases in the interval 0 < η < 1.7 and then becomes static. It is also observed that on the surface the temperature profile increases and then gradually decreases with time. Consequently, the momentum, as well as thermal boundary layer thickness, is decreasing
Figure 8. Influence of K on (a) f '(η), (b) θ(η), (c) ϕ(η), and (d) h(η) profiles, taking remaining parameters as M = 1.0, m = 0.1, Gr = 0.3, Gm = 0.3, β = 0.1, Le = 5, R = 0.4, Pr = 10, Nt = 0.1, Nb = 0.2, Q = 0.5, Bi = 1.0, A = 0.2, B1 = 0.01, G1 = 0.1, G2 = 2.0, s = 0.1.
with enhances of K. The chemical reaction insertion in the flow scheme is a well-organized method to control the mass and heat transfer. Therefore, the chemical reaction guides the amendment of the flow characteristics. On the other hand, the application of chemical reactions is interconnected with the fields of medicine, manufacturing, and so on. Therefore, the influence of chemical reaction (K) on the concentration profile is displayed in Figure 8c and found that the concentration profile is increasing (decreasing) in the interval 0 ≤ η < 0.9 (0.9 ≤ η < 1.9) with the enhancement of K. Consequently, the concentration boundary layer thickness grows up and declines due to the enhancement
of K the neighborhood of the wedge surface and far from the wedge surface, respectively. For the effect of chemical reaction (K) the angular velocity increases in the range 0 ≤ η < 2.5 and there is no effect of K on angular velocity in the range 2.5 ≤ η < 4 , finally, the angular velocity follow the static condition for any η (η ≥ 4) depicts in Figure 8d. Effect of Radiation Parameter (R) The effect of radiation parameter (R) on velocity, temperature, concentration, and angular velocity profiles reveal in Figures 9a, b, c, and d, respectively. It is observed from Figure 9a that R has an increasing effect on the velocity profile in the range 0 ≤ η < 2.2 and then
Figure 9. Influence of R on (a) f '(η), (b) θ(η), (c) ϕ(η), and (d) h(η) profiles, taking remaining parameters taking as M = 1.0, m = 0.1, Gr = 0.3, Gm = 0.3, β = 0.1, Le = 5, Pr = 10, Nt = 0.1, Nb = 0.2, Q = 0.5, K = 0.5, Bi = 1.0, A = 0.2, B1 = 0.01, G1 = 0.1, G2 = 2.0, s = 0.1.
becomes static. Figures 9b and c expound the temperature and concentration profiles, respectively for the variation of radiation parameters, R. From Figure 9b it is observed that the temperature profile is increasing for the effect of R in the range 0 ≤ η < 1.8, consequentl, the thermal boundary layer increases as well. The temperature profile and the thermal boundary layer rise because more energy enters the system due to the increased values of the radiation parameter. It is seen from Figure 9c that there exists a point of intersection and the concentration profile is decreasing before the point and after the point, the nanoparticle volume fraction profile is increasing for enhancement of R. Consequently, the concentration
boundary layer is decreasing as well as increasing for the effect of R in the surrounding area and away from the wedge surface, respectively. In the vicinity of the wedge surface, the angular velocity profile decreases with the increasing effect of R depicted in Figure 9d. Effect of Vortex Viscosity Constraint (G1) Figure 10 demonstrates the velocity, temperature, concentration, and angular velocity profiles for G1 = 0.020, 0.065, 0.070, and 0.100, taking the remaining parameters are constant. Figure 11a displays the effect of G1 on the velocity profile. Velocity profile as well as momentum boundary layer thickness increase in the range 0 ≤ η < 2.5 and then follow the static condition for large values of η (η
Figure 10. Influence of G1 on (a) f '(η), (b) θ(η), (c) ϕ(η), and (d) h(η) profiles, taking remaining parameters as M = 1.0, m = 0.1, Gr = 0.3, Gm = 0.3, β = 0.1, Le = 5, R = 0.4, Pr = 10, Nt = 0.1, Nb = 0.2, Q = 0.5, K = 0.5, Bi = 1.0, A = 0.2, B1 = 0.01, G2 = 2.0, s = 0.1.
≥ 2.5). In the vicinity of the wedge surface, the temperature (concentration) profiles decrease which reveals Figure 10b (Figure 10c). From Figure 10c, it is observed that being far away from the wedge surface G1 leads to increasing the concentration profile. Figure 10d exposes the effect of G1 on the angular velocity profile and found that in the range 0 ≤ η < 2.2 the angular velocity profile decreases, and at a certain range 2.2 ≤ η < 4.2 there is no effect of G1. Effect of Micropolar Parameter (G2) Figure 11a demonstrates the velocity profile G2 = 2.0, 2.5, 4.5, and 7.5 taking the remaining parameters as
constant. From Figure 11a it is seen that the velocity profile decreases in the range 0 ≤ η < 2.2 , while it becomes static in the range η ≥ 2.2. Figure 11b displays the effect of G2 the temperature profile and found that temperature profile as well as thermal boundary layer thickness decrease in the vicinity of the wedge surface (0 ≤ η < 1.45) and then follows the static condition for large values of η (η ≥ 1.45). Near the wedge surface, the concentration profile increases in the range 0 ≤ η < 1 revealed in Figure 11c. From Figure 11c, it is observed that being far away from the wedge surface G2 leads to decreasing the concentration profile in the range 1 ≤ η < 2. Figure 11d exposes
Figure 11. Influence of G2 on (a) f '(η), (b) θ(η), (c) ϕ(η), and (d) h(η) profiles, taking remaining parameters as M = 1.0, m = 0.1, Gr = 0.3, Gm = 0.3, β = 0.1, Le = 5, R = 0.4, Pr = 10, Nt = 0.1, Nb = 0.2, Q = 0.5, K = 0.5, Bi = 1.0, A = 0.2, B1 = 0.01, G1 = 0.1, s = 0.1.
the effect of G2 on the angular velocity profile and found that in the range 0 ≤ η < 1.5 the angular velocity profile increases, and at a certain range 1.5 ≤ η < 5.0 the angular velocity decreases for the effect of G2. Effect of Coupling Constant Parameter (B1) Figure 12a demonstrates the velocity profile B1 = 0.1,0.2,0.4, and 0.5, taking the remaining parameters as constant. From Figure 12a it is observed that B1 has an insignificant decreasing effect on the velocity profile in the range 0 ≤ η < 2.5 , while it becomes static in the range η ≥ 2.5. Figure 12b displays the temperature profile for B1 = 0.1,0.2,0.4, and 0.5 and and found that temperature profile
as well as thermal boundary layer thickness decrease in the vicinity of wedge surface (0 ≤ η < 1.5) and then follow the static condition for large values of η (η ≥ 1.5). Near the wedge surface, the concentration profile increases in the range 0 ≤ η < 0.8 revealed in Figure 12c. From Figure 12c, it is observed that being far away from the wedge surface B1 leads to decreasing the concentration profile in the range
0.8. ≤ η < 1.7. Figure 12d exposes the effect of B1 on the
angular velocity profile and found that in the range 0 ≤ η <
1.5. ≤ η < 5.0 there is no effect of B1.
Figure 12. Influence of B1 on (a) f '(η), (b) θ(η), (c) ϕ(η), and (d) h(η) profiles, taking remaining parameters as M = 1.0, m = 0.1, Gr = 0.3, Gm = 0.3, β = 0.1, Le = 5, R = 0.4, Pr = 10, Nt = 0.1, Nb = 0.2, Q = 0.5, K = 0.5, Bi = 1.0, A = 0.2, G1 = 0.1, G2 = 2.0, s = 0.1.
Table 1. Skin friction coefficient (−f ''(0)), Nasselt number (−θ'(0)), Sherwood number (−ϕ'(0)), and Couple stress coefficient (−h'(0)) for distinct values of M, β, Le, Nt, Nb, K, and A, taking remaining parameters as: m = 0.1, Gr = 0.3, Gm = 0.3, R = 0.4, Pr = 10, Q = 0.5, Bi = 1.0, B1 = 0.01, G1 = 0.1, G2 = 2.0, and s = 0.1. M
Table 2. The correlation coefficient for the mentioned parameters of the skin friction f ''(0), the Nusselt number θ'(0), the Sherwood number ϕ'(0), and the couple stress coefficient h'(0) by taking 2nd order coefficient M
Correlation Analysis From Table 2 it is seen that β, Nb, K, and A are positively correlated with a velocity gradient but M, Le, and Nt are negatively correlated with a velocity gradient. Higher values of one variable tend to result in higher values of another
variable whenever there is a positive correlation, but the opposite is true when there is a positive correlation. If the numerical values are 0.25 or above, the correlations between variables are statistically significant. So, from Table 2 it is also observed that the temperature gradient is positively
correlated with M, Nb, and K but negatively correlated with β, Le, Nt, and A respectively. The concentration is positively correlated with β, Le, and Nt but negatively correlated with M, Nb, K and A. Again, the couple stress coefficient is positively correlated with β, Le, M, Nb, K, and A but negatively correlated with Nt.
Conclusion
Under consideration of the impact of the concerned physical parameters, a computational examination of the MHD boundary layer of micropolar fluid flow across a stretched wedge is carried out. The Buongiorno model is taken into account for this. The significant findings of this investigation are listed as : • The magnetic parameter act as a controlling parameter in this study because fluid velocity and the rate of temperature enhance but the rate of mass transfer decreases within the boundary layer region consequence of the magnetic parameter. • The Brownian motion parameter acts as a source of cooling but the thermophoresis parameter acts as a heat source parameter because, within the boundary layer region, the temperature of the fluid decreases for the effect of the Brownian motion parameter on the other hand the reverse results arises for thermophoresis parameter. • The stretching or velocity ratio parameters act as a formation of the boundary layer. While the surface velocity is less than the free stream velocity, the structure of the boundary layer will be inverted but the opposite shape is formed when the surface velocity is greater than the free stream velocity. In the study, one can find the reverse inverted shape boundary layer.
Nomenclature
Latin alphabet Specific heat at constant pressure [Jkg−1 K −1] Cp Strength of transverse magnetic field [kg s−2 A−1] B0 u Velocity component of micropolar along x- direction [ms-1] v Velocity component of micropolar along y- direction [ms-1] Ω Wedge angle β The gradient of Hartree pressure [Nm-3] σ* Stefan-Boltzmann constant [Wm-2K-4] k* Mean absorption coefficient The stretching velocity of the surface [ms-1] Uw U Free stream velocity of the fluid [ms-1] f (η) Dimensionless stream function f '(η) Dimensionless velocity f ''(η) Dimensionless velocity gradient a Stretching rate constant b The free stream rate constant Heat generation Q0
Chemical reaction Volumetric coefficient of thermal expansion [K−1] Volumetric coefficient expansion of nanoparticle volume fraction Acceleration due to gravity vector [ms-2] The dimensional temperature of the fluid [K] Convective fluid temperature [K] Ambient nanofluid temperature [K] The dimensional nanoparticle volume fraction [kg m-3] Nanoparticle volume fraction at the stretching surface [kg m-3] Free stream nanoparticle volume fraction [kg m-3] Brownian diffusion coefficient [cm2 s -1] Thermophoresis diffusion coefficient Radiative heat flux [Wm-2] Coefficient of heat transfer Coupling constant micro-rotation or angular velocity [rads-1] Thermal conductivity [Wm-1K] Microinertia density
Greek alphabet γ Spin gradient τ The ratio of nanoparticle effective heat capacity and base fluid heat capacity τw Shear stress [Nm-2] σ Electrical conductivity [Sm-1] ρ The density of the base fluid [kgm-3] (ρCp)f Micropolar heat capacity α Thermal diffusivity of base fluid [m2s-1] ν Kinematic viscosity of the base fluid [kg m-3] J Stream function [m-2s-1] θ(η) Dimensionless temperature θ'(η) Rate of temperature ϕ(η) The dimensionless nanoparticle volume fraction ϕ'(η) Rate of nanoparticle volume fraction Dimensionless numbers/parameters M Magnetic parameter Gr Grashof number Gm Modified Grashof number Pr Prandtl number R Thermal radiation parameter Nb Brownian motion parameter Nt Thermophoresis parameter Q Heat generation parameter Sc Schmidt number Bi Biot number A Stretching parameter Le Lewis number K Chemical reaction parameter B1 Coupling constant parameter G1 Vortex viscosity constraint G2 Micropolar parameter
Transfer rate of heat per unit area of the surface Transfer rate of mass per unit area of the surface skin friction coefficient Nusselt number Sherwood number
Abbreviations
BL Boundary layer MHD Magnetohydrodynamic PDE’s Partial Differential Equations ODEs Ordinary Differential Equations CODEs Coupled Ordinary Differential Equations BCs Boundary conditions
Data Availability Statement
The authors confirm that the data that supports the findings of this study are available within the article. Raw data that support the finding of this study are available from the corresponding author, upon reasonable request.
Conflict Of Interest
The author declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Ethics
There are no ethical issues with the publication of this manuscript.
References
- [1] Eringen AC. Theory of micropolar fluids. J Math Mech 1966;6:1–18. [CrossRef]
- [2] Reddy GVR, Krishna YH. Soret and Dufour effects on MHD micropolar fluid flow over a linearly stretching sheet through a non-darcy porous medium. Int J Appl Mech Eng 2018;23:485–502. [CrossRef]
- [3] Srinivasacharya D, Upendra M. Effect of double stratification on MHD free convection in a micropolar fluid. J Egypt Math Soc 2013;21:370–378. [CrossRef]
- [4] Khalid A, Khan I, Khan A, Shafie S. Conjugate transfer of heat and mass in the unsteady flow of a micropolar fluid with wall couple stress. AIP Adv 2015;5:127125. [CrossRef]
- [5] Choi SUS, Eastman JA. Enhancing thermal conductivity of fluids with nanoparticles. Proc Int Mech Eng Congress and Exhibition. 1995, San Francisco, CA(United States), Nov. 12–17. [CrossRef]
- [6] Mebarek-Oudina F, Chaban I. Review on nano-fluids applications and heat transfer enhancement techniques in different enclosures. J Nanofluids 2022;11:155–168. [CrossRef]
- [7] Gul T, Usman M, Khan I, Nasir S, Saeed A, Khan A, Ishaq M. Magneto hydrodynamic and dissipated nanofluid flow over an unsteady turning disk. Adv Mech Eng 2021;13:1–11. [CrossRef]
- [8] Mahdy A, Chamkha AJ. Unsteady MHD boundary layer flow of tangent hyperbolic two-phase nanofluid of moving stretched porous wedge. Int J Numerical Methods Heat Fluid Flow 2018;28:2567–2580. [CrossRef]
- [9] Abderrahmane A, Qasem NAA, Younis O, Marzouki R, Mourad A, Shah NA, Chung JD. MHD hybrid nanofluid mixed convection heat transfer and entropy generation in a 3-D triangular porous cavity with zigzag wall and rotating cylinder. Math 2022;10:769. [CrossRef]
- [10] Islam T, Yavuz M, Parveen N, Fayz-Al-Asad Md. Impact of non-uniform periodic magnetic field on unsteady natural convection flow of nanofluids in a square enclosure. Fractal Fract 2022;6:101.
- [11] Vaidya H, Rajashekhar C, Mebarek-Oudina F, Prasad KV, Vajravelu K, Ramesh BB. Examination of chemical reaction on three-dimensional mixed convective magnetohydrodynamic Jeffrey nanofluid over a stretching sheet. J Nanofluids 2022;11:113–124. [CrossRef]
- [12] Abbas Z, Hussain S, Naveed M, Nadeem A, Ali A. Analysis of joule heating and generalized slip flow in ferromagnetic nanoparticles in a curved channel using Cattaneo-Christof heat flux theory. Therm Sci 2022;26:437–448. [CrossRef]
- [13] Jilte R, Mohammad HA, Ravinder K, Vilas K, Amirhosein M. Cooling performance of a novel circulatory flow concentric multi-channel heat sink with nanofluids. Nanomaterials 2020;10:647. [CrossRef]
- [14] Jilte R, Ravindra K, Mohammad HA. Cooling performance of nanofluid submerged vs. nanofluid circulated battery thermal management systems. J Cleaner Prod 2019;240:118131. [CrossRef]
- [15] Buongiorno J. Convective Transport in Nanofluids. ASME J Heat Transf 2006;128:241–250. [CrossRef]
- [16] Abdullah A, Ibrahim F, Chamkha A. Nonsimilar solution of unsteady mixed convection flow near the stagnation point of a heated vertical plate in a porous medium saturated with a nanofluid. J Porous Media 2018;21:363–388. [CrossRef]
- [17] Mabooda F, Mastroberardino A. Melting heat transfer on MHD convective flow of a nanofluid over a stretching sheet with viscous dissipation and second-order slip. J Taiwan Int Chem Eng 2015;57:62–68. [CrossRef]
- [18] Mahdy A, Chamkha AJ. Heat transfer and fluid flow of a non-Newtonian nanofluid over an unsteady contracting cylinder employing Buongiorno’s model. Int J Numerical Methods Heat Fluid Flow 2015;25:703–723. [CrossRef]
- [19] Pavlov KB. The magnetohydrodynamic flow of an incompressible viscous fluid is caused by the deformation of a plane surface. Magn Gidrodin 1974;4:146–147.
- [20] Andersson HI. MHD flow of a viscoelastic fluid past a stretching surface. Acta Mech 1992;95:227–230. [CrossRef]
- [21] Kasiviswanathan SR, Gandhi MV. A class of exact solutions for the MHD flow of micropolar fluid. Int J Eng Sci 1992;30:409–417. [CrossRef]
- [22] Chamkha AJ, Rashad AM, Armaghani T, Mansour MA. Effects of partial slip on entropy generation and MHD combined convection in a lid-driven porous enclosure saturated with a Cu–water nanofluid. J Therm Anal 2018;132:1291–1306. [CrossRef]
- [23] Falkner VM, Skan SW. Some approximate Solutions of the boundary-layer equations. Lond Edinb Dublin Philos Mag J Sci 1931;12:865–896. [CrossRef]
- [24] Schlichting H, Gersten K. Boundary-Layer Theory. Berlin, Germany: Springer; 2016. [CrossRef]
- [25] Falkner VM, Skan SW. Some approximate solutions of the boundary layer for flow past a stretching boundary. SIAM J Appl Math 1931;49:1350–1358. [CrossRef]
- [26] Goyal KP, Kassoy DR. Heat and mass transfer in a saturated porous wedge with impermeable boundaries. Int J Heat Mass Transf 1979;22:1577–1585. [CrossRef]
- [27] Stewartson K. Further solutions of the Falkner-Skan equation. Math Procs Cambridge Philos Soc. 1954;50:454–465. [CrossRef]
- [28] Riley, Weidman PD. Multiple solutions of the Falkner-Skan equation for flow past a stretching boundary. SIAM J Appl Math 1989;49:1350–1358. [CrossRef]
- [29] Dimian MF. Effect of the magnetic field on forced convection flows along a wedge with variable viscosity. Mech Mech Eng 2004;7:107–117.
- [30] Ishak A, Nazar R, Pop I. Falkner-Skan equation for flow past a moving wedge with suction or injection. J Appl Math Comput 2007;25:67–83. [CrossRef]
- [31] Fang T, Zhang J. An exact analytical solution of the Falkner-Skan equation with mass transfer and wall stretching. Int J Non-Linear Mech 2008;43:1000–1006. [CrossRef]
- [32] Dogonchi AS, Chamkha AJ, Seyyedi SM, Ganji DD. Radiative nanofluid flow and heat transfer between parallel disks with penetrable and stretchable walls considering Cattaneo–Christov heat flux model. Heat Transf Asian Res 2018;47:735–753. [CrossRef]
- [33] Zheng L, Zhang C, Zhang X, Zhang J. Flow and radiation heat transfer of a nanofluid over a stretching sheet with velocity slip and temperature jump in a porous medium. J Franklin Inst 2013;350:990–1007. [CrossRef]
- [34] Ahmed SE, Hussein AK, Mansour MA, Raizah ZAS, Zhang X. MHD mixed convection in trapezoidal enclosures filled with micropolar nanofluids. Nanosci Tech 2018;9:343–372. [CrossRef]
- [35] Bagh A, Khan SA, Hussein AK, Thumma T, Hessain S. Hybrid nanofluids: Significance of gravity modulation, heat source/sink, and magnetohydrodynamic on dynamics of a micropolar fluid over an inclined surface via finite element simulation. Appl Math Comp 2022;419:126878. [CrossRef]
- [36] Kasmani R, Sivasankaran S, Bhuvaneswari M, Hussein AK. Analytical and numerical study on convection of nanofluid past a moving wedge with Soret and Dufour effects. Int J Numerical methods Heat Fluid Flow 2017;27:2333–2354. [CrossRef]
- [37] Ashraf M, Rehman S, Farid S, Hussein AK, Ali B, Shah N, Weera W. Insight into the significance of bioconvection on MHD tangent hyperbolic nanofluid flow of irregular thickness across a slender elastic surface. Math 2022;10:2592–2608. [CrossRef]
- [38] Ahmed S, Hussein AK, Mohammed H, Sivasankaran S. Boundary layer flow and heat transfer due to permeable stretching tube in the presence of heat source/sink utilizing nanofluids. Appl Math Comp 2014;238:149–162. [CrossRef]
- [39] Sudarsana RP, Chamkha AJ. Soret and Dufour effects on MHD heat and mass transfer flow of a micropolar fluid with thermophoresis particle deposition. J Naval Arch Marine Eng 2016;13:1813–8235. [CrossRef]
- [40] Narender G, Govardhan K, Sarma S. Viscous dissipation and thermal radiation effects on the flow of maxwell nanofluid over a stretching surface. J Nonlinear Anal Appl 2021;12:1267–1287.
- [41] Mishra SR, Khan I, Al-Dallal QM, Asifa T. Free convective micropolar fluid flow and heat transfer over a shrinking sheet with a heat source. Case studies Therm Eng 2018;11:113–119. [CrossRef]
Share and Cite
HANI, U.; ALI, M.; ALAM, M.S. MHD boundary layer micropolar fluid flow over a stretching wedge surface Thermophoresis and brownian. Journal of Thermal Engineering 2024, Vol. 10, pp. 330-349. https://doi.org/10.18186/thermal.1448609

