MHD Nanofluid Flow Past a Vertical Plate Embedded in a Rotating Porous Medium with Chemical Reaction ()
1. Introduction
Nanofluids are a relatively new technology consisting of nanoparticles mixed into a base fluid. Particles used in such a type of fluid are so small that each one is between 1 and 100 nanometres in size, which is made up of metal oxides, carbon nanotubes or other materials, while the base fluid is mainly a liquid such as water, oil or ethanol. This mixture forms a fluid with properties that are much more efficient than those of traditional fluids [1]. It was Choi [2] who first used the expression “nanofluid”. Nanofluids are found to have a higher thermal conductivity compared to the pure fluid [3]. In particular, the higher thermal conductivity allows nanofluids to be used as efficient coolants in electronics and other heat-generating devices. Nanofluids can also be engineered to have specific electrical conductivity and energy storage properties. In the medical field, nanofluids are also being used as drug delivery systems that can target specific areas of the body, and this type of delivery system can significantly improve the effectiveness of treatments and reduce the side effects of drugs. Earliest studies in the field of nanofluid by [4]-[10].
Many researchers investigate the field of nanofluid, among them Salem [11] who studied the effects of variable viscosity, viscous dissipation and chemical reaction on heat and mass transfer. Hammed and Pop [12] studied the unsteady free convection flow of a nanofluid past a permeable semi-infinite plate in the presence of an applied magnetic field, theoretically. They have succeeded in approximating an analytical solution to the boundary layer equations by means of the perturbation method. Loganathan and Sangeetha [13] used the implicit finite difference method to analyse the augmentation of thermal conductivity and mass transfer rate across a semi-infinite vertical plate due to nanoparticles suspended in Williamson fluid. Ullah et al. [14] numerically studied the two-dimensional Oldroyd-B flow of a nanofluid past a stretching sheet using the finite element method, under the effects of magnetic, electrical and thermal radiation. In an important study, Krishna and Chamkha [15] used the perturbation approximation approach to study the unsteady free convective flow of nanofluids at room temperature over a moving vertical plate. The plate is immersed in a homogeneous porous medium under thermal buoyancy with a constant heat source and convective boundary conditions. Furthermore, Krishna and Chamkha [16] studied unsteady MHD convective rotating flow of a nanofluid past a permeable moving vertical plate with a constant heat source using the perturbation method. Riasat et al. [17] used nanofluid with multi-walled carbon nanotubes between two rotating discs in a Darcy-Forchheimer permeable medium subject to an induced magnetic field to improve thermal conductivity in squeezing 3D magnetohydrodynamic thin film flow. Maneengam et al. [18] carried out numerical study of the flow of a hybrid nanofluid within a corrugated porous cavity in the presence of a magnetic field using the generalised finite element technique. Abumandour et al. [19] investigated the effects of slip, temperature conditions, and porosity on nanoparticle fluid suspension as well as the thermal nanofluid behavior as it flowed within a stenotic artery. Utilizing a catheter and vertical orientation, they explored the role of magnetic fields both with and without their presence. Using the homotopy analysis method, Ur Rasheed et al. [20], successfully evaluated the mixed convective MHD Jeffrey nanofluid flow. The study was conducted in the presence of a vertical stretchable cylinder with heat and radiation absorption effects. Hossain et al. [21] used the finite volume method to study the bio-nanofluid flow regimes in a real bifurcated artery. The study showed that the influence of a magnetic source induced gradual changes in the pressure and vorticity profiles. Several investigations on nanofluid properties and MHD flow are thoroughly detailed in the literature [22]-[29].
The coupling of heat and mass transfer through porous media has a wide range of applications in thermal engineering, physics, mechanical engineering and chemistry, such as the design and optimization of heating, cooling, ventilation and air conditioning systems in various buildings and industries. The rate of heat and mass transfer through porous materials is influenced by the physical properties of these materials, including thermal conductivity, density, volume, temperature, humidity, and others. Porous media refers to materials that possess interconnected pores, channels, or voids. Porous materials are characterized by their porosity, which is a measure of the interconnected voids in the material. The intricate dynamics of nanofluid flow in a rotating porous media is a challenging subject to comprehend. It often requires the application of advanced mathematical models and simulations to accurately predict the transport phenomena. Moreover, the interplay between the nanoparticles and the pore surfaces plays a crucial role in determining the flow behaviour of nanofluids in porous media. The use of porous media is very encouraging, as it provides new insights into optimising the flow of nanofluids, thereby improving the efficiency of heat transfer and facilitating the sustainable use of energy. Reddy et al. [30] studied an unsteady natural convective flow of a nanofluid through a porous medium in a rotating system with convective and diffusive boundary conditions. They have solved the governing equations analytically using the perturbation approach. Kumar & Sood [31] investigated the three-dimensional MHD nanofluid flow over a vertical surface embedded in a porous rotating system. The principles of homogeneous-heterogeneous chemical reaction are introduced simultaneously. The governing equations are converted to dimensionless form and solved explicitly using a finite difference approach. Reddy et al. [32] demonstrated MHD flow of Cu-water and Ag-water nanofluids with heat and mass transfer in the presence of chemical reaction, thermal radiation and partial slip through a porous medium. They used similarity variables to convert the governing equations into ordinary differential equations, which were then solved numerically using the finite element method. Rashid et al. [33] investigated the MHD hybrid nanofluid flow between two parallel spinning discs in present of joule heating and chemical reactions. They utilised the Matlab programme “bvp4c” to numerically solve the governing equations. Recently, Abd-Alla et al. [34] introduced heat and mass transfer analysis on MHD Jeffery fluid flow through a porous medium utilising a rotating frame with a chemical reaction. Recently, nanofluids flow through rotating porous media with or without chemical reaction investigated widely by [35]-[41].
Magnetohydrodynamics convective flows with simultaneous heat and mass transfer influenced by chemical reactions occur in many natural and artificial transport systems for instance in scientific and technological areas. Such applications of this phenomenon can be found in chemical vapour deposition on surfaces, nuclear reactor cooling, the petroleum industries, power and cooling sectors. In addition, convective flows occur naturally in temperature variations, concentration differences or a mixture of the two; for example, in air currents, concentration differences in the water where flow is affected by these concentration differences. Elsewhere this type of flow, which takes place through a porous plate, is of crucial importance in various scientific and mechanical applications. This prompted several scholars [42]-[48] to do extensive study on heat and mass transport of MHD free convective flow of fluid past a rotating plate with chemical reaction under different conditions. Prabhakar Reddy & Makinde [49] investigated a heat and mass transfer on MHD free convection flow under the influence of a chemical reaction. Also, Krishna et al. [50] studied numerically the problem of unsteady MHD natural convective rotating flow of an incompressible electrically conducting viscous fluid past an infinite vertical moving porous surface using cubic Bsplines collocation method. Very recently, Ayub et al. [51] investigated the problem of free convection on rotational flow of viscous fluid over an extended plate with heat and mass transfer in the presence of a constant-magnitude magnetic field through a porous medium.
Due to its importance in heat and mass transfer in some types of fluid flow, variable viscosity has been widely studied, recently. In a study by Ali et al. [52], temperature-dependent viscosity, convective boundary conditions, thermo-diffusion, and radiation effects were all taken into consideration to analysed the problem of unsteady MHD nanofluid flow across a radially nonlinear stretching sheet. They used the variational finite element method to solve the governing equations. A numerical study was carried out by Khan et al. [53] to analysed the impact of the variable viscosity on the flow of a Casson nanofluid over non-linear stretching sheet with slip and convective boundary conditions. With the motivation from all of the above efforts on a fact, specialised literature investigation tolerates that the effects of temperature-dependent viscosity, an induced magnetic field, nanofluid fluid, natural convection flow with heat, and mass transfer in the presence of a first-order chemical reaction through a semi-infinite vertical rotating plate embedded in a porous medium have not yet been discussed. Therefore, the current study claims to extend the Prabhaker Reddy and Makinde [49] for it is limitation of considering unsteady MHD free convective flow past through vertical plate in presence of chemical reaction such that the current study offered advancements of adding more specific consideration of variable viscosity, induced magnetic field with viscous dissipation and Joule heating, rotating plate, porous media saturated in fluidity with Copper(Cu) nanoparticles. Then to obtain the problem solution the novel BI-SRM numerical approach [54] [55] was chosen because it convincingly demonstrated the best in terms of accuracy and convergence speed, and computational cost is lacking as a direct comparison with other established numerical methods (e.g., finite element, finite difference).
2. Mathematical Formulation
Let us consider an unsteady MHD boundary layer flow of an incompressible nanofluid driven by a semi-infinite vertical rotating plate embedded in a porous medium subject to chemical reaction. A coordinate frame is selected in which x-axis runs vertically along the plate, the y-axis is normal to the plate and directed into the fluid zone, in which the z-axis is taken to be normal to the xy-plane. Assumed that the nanofluid is formed with non-magnetic nanoparticles and base liquid with the velocity
. Let
be the induced magnetic vector at the given point in the fluid. A homogeneous transverse magnetic field of intensity H is applied parallel to the y-axis and perpendicular to the plate. Together, the fluid and plate rotate about y-aixs with a constant angular velocity Ω in the anticlockwise direction. Initially (
), both the fluid and plate are kept at a constant temperature
where
is for a heated plate and
corresponding to a cooled plate. At a certain time
, the temperature of the plate is growing up to
here
is the characteristic time. Moreover, the distribution of species concentration over the plate is boosted to
. It is assumed that a homogeneous chemical reaction of first order with uniform rate
between the diffusing species and the ambient fluid exists. Figure 1 shows the geometry of the flow problem. The boundary layer equations controlling the flow include the continuity, motion, energy, concentration, and induction under the using Boussinesq’s standard approximation respectively:
![]()
Figure 1. Physical model.
(1)
(2)
(3)
(4)
(5)
(6)
(7)
where u and w represents the velocity components in the x and z directions, respectively, t is the time,
and
are the wall-dimensional temperatures and concentrations, respectively,
and
, respectively are the free stream’s starting temperature and concentration,
is the nanofluid density,
is the nanofluid dynamic viscosity,
is the nanofluid electrical conductivity,
is the nanofluid thermal diffusivity, T is the dimensional temperature, C is the dimensional mass concentration,
is the specific heat of the nanofluid at constant pressure, g is the gravitational acceleration,
and
are the nanofluid thermal and concentration expansion coefficients respectively,
and k are the magnetic and porous plate permeabilities respectively,
stands for the chemical reaction parameter.
The problem’s initial boundary conditions are similar to those of Prabhaker Reddy and Makinde [49], and are presented as follows
(8)
where Initial velocity and time are
and
, respectively. As a function of the nanoparticle, the following definitions for effective density, effective dynamic viscosity, thermal expansion coefficient, thermal diffusivity, heat capacitance, thermal conductivity, and electrical conductivity (see Minea [56], and Hamad & Pop [12]) are given, respectively:
(9)
where
denotes the volume of the nanoparticles that are made of solid material,
and
denote the densities of the fluid and solid particles, respectively, and
and
denote the coefficients of thermal expansion of the solid and fluid, respectively. The thermal conductivity of the base fluid and the solid are
and
, respectively, whereas the electrical conductivity of the fluid and the solid are
and
. Table 1 lists the thermo-physical properties of the base fluid (water) and copper that are utilized for code validation.
Table 1. Base fluid and nanoparticle thermophysical properties.
Thermo-physical properties |
Water |
Copper (Cu) |
Heat capacity
|
4179 |
385 |
Density
|
997.1 |
8933 |
Thermal conductivity
|
0.613 |
400 |
Thermal expansion
|
21 |
1.67 |
Electrical conductivity
|
0.05 |
5.96 × 107 |
The fluid dynamic viscosity
is considered to vary as an inverse function of temperature as given by Lai & Kulacki [57]:
(10)
The transformation variables in dimensionless form are introduced as follows:
(11)
Such that, Equations (2)-(7) can be written in the following dimensionless form (dropping asterisks):
(12)
(13)
(14)
(15)
(16)
(17)
Here
is the Rotation parameter,
is the Permeability parameter,
is the Reynolds number,
is the local temperature Grashof number,
is the local mass Grashof number,
is the Magnetic field parameter,
is the Prandtl number,
is the Eckert number,
is the Joule heating parameter,
is the Schmidt number,
is the Chemical reaction parameter,
is the Magnetic Reynolds number,
is the temperature difference and where,
The boundary conditions (8) in nondimensional form are
(18)
Local skin friction coefficient:
The skin friction coefficient is a dimensionless quantity defined by the shear stress on the wall:
(19)
where
and
are the values of the nanofluid density and of the longitudinal velocity near the boundary layer’s edge. Accordingly, at the plate with
, the primary and secondary shear stresses
and
are given by:
(20)
3. Numerical Method
This part is constructed to show the implementation of the Bivariate Spectral Relaxation Method (BI-SRM). The BI-SRM is a very effective tool that can be used in solving the highly nonlinear differential equations. The BI-SRM, a relatively new technique, improves upon the Spectral Relaxation Method (SRM)(see Motsa et al. [58]), which is widely used for the numerical solution of the partial differential equations involved in the modeling of fluid flow. BI-SRM is introduced in this section in order to solve the system of the nonlinear partial differential Equations (12)-(17) along with the boundary conditions in (18). In order to make the spectral collocation method applicable, BI-SRM linearises the non-linear PDEs (12)-(17) to have
(21)
(22)
(23)
(24)
(25)
(26)
subject to
(27)
The variable coefficients in above equations are determined by
where
and r denote the current and previous iterations, respectively. The Chebyshev spectral collocation is adopted to discretize the linear partial differential Equations (21)-(26) in both time (t) and space (y) directions. However, the spectral method is applicable only in the domain
, due to this, we have to convert the physical domain of the problem from
into
by suitable linear transformations. The linear partial differential Equations (21)-(26)) are decoupled, so each of them can be solved separately and their solutions can be approximated by a bivariate Lagrange polynomial standard form as follows:
(28)
(29)
(30)
(31)
(32)
(33)
To estimate variables in Equations (28)-(33) at predetermined points in the indicated directions
and
we introduce the Chebyshev-Gauss-Lobatto points
(34)
and
, respectively, denote the number of grid points in the
and
directions. Additionally, the discrete derivatives at the grid points are obtained by transforming the continuous time and space derivatives using Chebyshev-Gauss-Lobatto (CGL) (34) grid points. Moreover, the common Lagrange cardinal polynomials cardinal polynomials
and
are defined in the form
(35)
where
(36)
The definition of the function
is identical to that of the function
. Solution procedure starts by the determination of the derivative of
and
with respect to
and
, therefore, calculating the Chebyshev partial derivative of
with respect to
at any point
is given by
(37)
(38)
(39)
When
, the “Chebyshev Differentiation Matrix” [59] is utilized to produce the matrix’s elements
. Partial derivatives of
with respect to
are derived similarly. Similarly, at
CGL points, the partial derivative with respect to the space variable
is app
(40)
(41)
(42)
where i is the index nodes of
and p is the index basis of
, such that
(43)
where
denotes the matrix transpose. Consequently, the derivative of second order with respect to
is given by
(44)
Differentiation matrices of other dependent variables
are approximated by collocation, similar to Equations (39) and (42). Using the equations (39), (42) and (44) in Equation (24) we arrive at
(45)
According to the definitions of the dependent variables
and their derivatives, the Equations (21)-(26) can be written in the following operator form
(46)
(47)
(48)
(49)
(50)
(51)
where
Therefore, for
, the boundary conditions in Equation (49) are established as in the following
matrix system:
(52)
where
(53)
The boundary condition (27) is imposed similarly for the remaining variables, and the bivariate Chebyshev collocation technique is performed similarly to produce
(54)
(55)
(56)
(57)
(58)
The vectors
,
,
,
,
and
are the initial solutions derived from the supplied boundary conditions (27). After the boundary conditions are applied, Equations (54)-(58) are written as the following
matrices:
(59)
(60)
(61)
(62)
(63)
where
is the
unit matrix. The system of above matrices (59)-(63) is iteratively solved to get the numerical solution for
and
.
4. Results and Discussion
The results obtained by applying the Bivariate Spectral Relaxation Method (BI-SRM) to the problem of time-dependent MHD nanofluid flow driven by a semi-infinite vertical plate in a porous medium are carefully presented and thoroughly discussed in this section. As a check on the accuracy and convergent nature of the BI-SRM, we calculated the error of the method over orders see Figure 2 and Figure 3. From these figures, it is observed that the BI-SRM achieves convergence in only 12 iterations with a maximum error of order 10−10 for the velocity and temperature profiles and 10−11 for both
of the induction profiles. To understand the physics of the problem, we have shown in Figures 4-14 the effects of varying our governing parameters on velocity, temperature, concentration, induction fields as well as skin friction. Through all these figures, it is apparent that the mass volume fraction appearing leads to slowdown in the fluid motion. There are two potential explanations for this phenomenon. Firstly, it is worth
![]()
Figure 2. Residual error graphs.
Figure 3. Solution error graphs.
Figure 4. The velocity profiles for u and w against K and
.
Figure 5. The velocity profiles for u and w against Re and
.
Figure 6. The velocity profiles for u and w against R and
.
noting that larger nanoparticles undergo less Brownian motion, which increases resistance in the current and effectively reduces velocity. Secondly, raising nanoparticle concentration may increase fluid viscosity, thus causing a thicker momentum boundary layer, which could ultimately decelerate fluid motion.
The influence of permeability K of the porous media on the velocity profiles is displayed in Figure 4 for different values of nanoparticle volume fractions. It is observed that raising the permeability coefficient of the porous medium reduces the resistance between the fluid and medium, leading to a thinner momentum boundary layer, hence increasing the primary and secondary flow. Figure 5 shows the influence of the Reynolds number Re on the velocity components u and w. The Reynolds number relates the fluid momentum to the viscous shear stress in the fluid system. Increasing the Reynolds number decreases the ratio 1/Re, causing the momentum force to become negligible compared to the viscous forces. This enhances the thickness of the momentum boundary layer so that the velocity components u and w increase as 1/Re decreases. See Fig. 6 for the effect of the rotation parameter R on the primary velocity u and the secondary velocity w. It is evident that as the rotation parameter R increases, the primary velocity u decreases. This phenomenon occurs due to the dominance of the Coriolis force exerted by the fluid medium in the direction of rotation. For the same reason, it is also worth noting that the secondary velocity increases with the rotation parameter. Figure 7 illustrates how the magnetic field M affects the primary and secondary velocity, respectively.
It can be seen that the primary flow velocity decreases significantly as the magnetic parameter M increases. In a particular case, secondary flow velocity increases near the plate as M increases and then decreases away from the plate as M increases. In practice, this occurs due to the Lorentz force created by the magnetic field which acts against the flow. Figure 8 and Figure 9 show the variations of the primary and secondary velocities with respect to the thermal Grashof number Gr and the concentration Grashof number Gm, which represent the ratio of the thermal and concentration buoyancy forces to the viscous force, respectively. As it is known that the thermal and concentration Grashof numbers increase, Gr and Gm increase the effectiveness of the thermal and concentration buoyancy forces, resulting in an increase in the thermal Grashof number, the mass Grashof number, or both, which increases the thickness of the momentum boundary layers, thus increasing the velocities.
Figure 10 shows the velocity profiles of the primary and secondary components, u and w, respectively, for different values of the Prandtl number Pr. In fact, the ratio 1/Pr is a decreasing function of Pr, so in our fluid case, increasing Pr causes highly viscous flow. Therefore, both velocity components increase with increasing Pr. Figure 11 illustrates the effect of the chemical reaction parameter on the velocity components u and w. The results indicate that increasing the chemical reaction parameter leads to a reduction in both primary and secondary velocities, this is due to a decrease in fluid concentration.
Figure 7. The velocity profiles for u and w against M and
.
Figure 8. The velocity profiles for u and w against Gr and
.
Figure 9. The velocity profiles for u and w against Gm and
.
Figure 10. The velocity profiles for u and w against Pr and
.
Figure 11. The velocity profiles for u and w against Kc and
.
Figure 12. The velocity profiles for u and w against Sc and
.
Figure 12 shows the effect of Schmidt number Sc on the primary and secondary velocity profiles. As Sc increases, both u and w tend to decrease. This is because an increase in Sc leads to a decrease in molecular diffusivity, resulting in a thicker momentum boundary layer. The fluid velocity then decreases in both directions.
The primary and secondary velocity patterns for different values of the Reynolds magnetic number (Rem) which estimates the induction of magnetic diffusion by a magnetic field are plotted in Figure 13. Here is roughly estimated by the magnetic Reynolds number. As Rem increases, the primary velocity rises, while the secondary speed decreases. Physically, an increase in Rem leads to enhanced magnetic diffusion generates Lorentz drag, which interacts with fluid motion, resulting in a decrease in the primary velocity profile. Conversely, the secondary velocity exhibits the opposite trend for a while and then decreases due to the same indicated reasons. On the other hand, the side effects of the magnetic Reynolds number Rem on both axially induced magnetic fields, Hx and Hz, are shown in Figure 14. It can be seen that both Hx
Figure 13. The velocity profiles for u and w against Rem and
.
Figure 14. The magnetic profiles for Hx and Hz against Rem and
.
and Hz are decreasing functions in Rem. This is because as Rem decreases, the mass diffusivity increases and this accelerates the flow regime.
Figure 15 illustrates the temperature distribution for different values of Reynolds number Re and mass volume fraction
. The figure shows that as the Reynolds number increases, the heat flux cells become smaller and smaller in size. An increasing Reynolds number indicates that the inertial forces become dominant compared to the viscous force. This causes the fluid to accelerate, as a result, more heat is transported by convection. It also can be seen that increasing the amount of mass volume fraction has a positive effect on the thermal conductivity of the system. Figure 16 and Figure 17 are plotted to illustrate the impacts of Schmidt number Sc and Prandtl number Pr on temperature distribution, respectively. From Figure 16, it is very noticeable that iso-thermal cells become an increasing function of Sc. This means that by increasing Sc, the temperature is well distributed over the flow area. Conversely, it is clear that from Figure 17, the effects of the Prandtl number on the temperature distribution is to produce a completely opposite behavior. For different values of the Schemdit number Sc at Re = 2 and Re = 5, the distribution of the concentration was plotted in Figure 18. We found that the concentration distribution is strongly time dependent near the edge of the plate. When both Sc and Re are increasing, the concentration occupies a narrow region of the solution domain. Figure 19 illustrates the effect of the Eckert number Ec in the temperature distribution field. It is seen that increasing of Ec enhances the temperature distribution over the region domain. This demonstrates how the Ec number plays a significant role in the temperature distribution.
![]()
Figure 15. Temperature distribution for different values of Re and
.
Figure 16. Temperature distribution for different values of Re and Sc.
Figure 17. Temperature distribution for different values of Re and Pr.
Figure 20 and Figure 21 are plotted to illustrate the impact of the Reynolds magnetic number on both x-axis magnetic field Hx and y-axis magnetic field Hz. From these figures we can see that as the Reynolds magnetic number, the orientation and direction of the magnetic lines change from positive to negative or vice versa, see the first surface plot. The two separated groups of cells in the second contour plot are evidence of this. When the magnetic field is oriented in a
Figure 18. Concentration distribution for different values of Re and Sc.
Figure 19. Temperature distribution for different values of Re and Ec.
Figure 20. X-axis magnetic filed induction distribution for different values of Re and Rem
Figure 21. Z-axis magnetic filed induction distribution for different values of Re and Rem
direction opposite to the normal vector of the domain by changing the magnetic polarity, the flux is negative, or when an arrangement of magnets takes place.
Table 2 and Table 3 display the primary and secondary skin friction under the influence of M, K, Gr, Gm, R, Sc, Ec, Pr, Kc, Re, and Rem for various values of the
Table 2. Skin friction coefficient for primary velocity u.
M |
K |
Gr |
Gm |
R |
Ec |
Pr |
Sc |
Kc |
Re |
Rem |
|
|
|
10 |
0.5 |
10 |
10 |
1 |
0.5 |
1.7 |
0.5 |
1 |
0.5 |
0.5 |
4.3251 |
4.4641 |
4.7435 |
|
|
|
|
|
|
|
|
|
|
|
4.8302 |
4.9784 |
5.2754 |
|
|
|
|
|
|
|
|
|
|
|
5.3038 |
5.4617 |
5.7763 |
|
1 |
|
|
|
|
|
|
|
|
|
4.2670 |
4.3968 |
4.6637 |
|
1.5 |
|
|
|
|
|
|
|
|
|
4.2082 |
4.3287 |
4.5829 |
|
|
15 |
|
|
|
|
|
|
|
|
4.0050 |
4.1769 |
4.4831 |
|
|
20 |
|
|
|
|
|
|
|
|
3.6748 |
3.8826 |
4.2177 |
|
|
|
15 |
|
|
|
|
|
|
|
2.9961 |
3.0720 |
3.2836 |
|
|
|
20 |
|
|
|
|
|
|
|
1.6214 |
1.6329 |
1.7755 |
|
|
|
|
2 |
|
|
|
|
|
|
4.3351 |
4.4749 |
4.7541 |
|
|
|
|
3 |
|
|
|
|
|
|
4.3701 |
4.5122 |
4.7922 |
|
|
|
|
|
1 |
|
|
|
|
|
4.0937 |
4.2121 |
4.4572 |
|
|
|
|
|
1.5 |
|
|
|
|
|
3.8952 |
3.9971 |
4.2153 |
|
|
|
|
|
|
3 |
|
|
|
|
4.0543 |
4.1754 |
4.4241 |
|
|
|
|
|
|
4 |
|
|
|
|
3.8767 |
3.9834 |
4.2104 |
|
|
|
|
|
|
|
1 |
|
|
|
4.3400 |
4.4801 |
4.7607 |
|
|
|
|
|
|
|
1.5 |
|
|
|
4.3545 |
4.4958 |
4.7776 |
|
|
|
|
|
|
|
|
2 |
|
|
4.3397 |
4.4799 |
4.7604 |
|
|
|
|
|
|
|
|
3 |
|
|
4.3540 |
4.4952 |
4.7770 |
|
|
|
|
|
|
|
|
|
1 |
|
2.8157 |
2.9008 |
3.0787 |
|
|
|
|
|
|
|
|
|
1.5 |
|
2.1457 |
2.2164 |
2.3607 |
|
|
|
|
|
|
|
|
|
|
1 |
4.1546 |
4.2765 |
4.5344 |
|
|
|
|
|
|
|
|
|
|
1.5 |
4.0186 |
4.1321 |
4.3792 |
Table 3. Skin friction coefficient for secondary velocity w.
M |
K |
Gr |
Gm |
R |
Ec |
Pr |
Sc |
Kc |
Re |
Rem |
|
|
|
10 |
0.5 |
10 |
10 |
1 |
0.5 |
1.7 |
0.5 |
1 |
0.5 |
0.5 |
2.7835 |
2.8359 |
2.9276 |
|
|
|
|
|
|
|
|
|
|
|
2.8552 |
2.9123 |
3.0114 |
|
|
|
|
|
|
|
|
|
|
|
2.9255 |
2.9869 |
3.0928 |
|
1 |
|
|
|
|
|
|
|
|
|
2.8016 |
2.8563 |
2.9507 |
|
1.5 |
|
|
|
|
|
|
|
|
|
2.8201 |
2.8770 |
2.9743 |
|
|
15 |
|
|
|
|
|
|
|
|
2.8559 |
2.8993 |
2.9833 |
|
|
20 |
|
|
|
|
|
|
|
|
2.9314 |
2.9649 |
3.0406 |
|
|
|
15 |
|
|
|
|
|
|
|
3.0405 |
3.0965 |
3.1901 |
|
|
|
20 |
|
|
|
|
|
|
|
3.3038 |
3.3632 |
3.4589 |
|
|
|
|
2 |
|
|
|
|
|
|
3.1750 |
3.2497 |
3.3655 |
|
|
|
|
3 |
|
|
|
|
|
|
3.5540 |
3.6504 |
3.7901 |
|
|
|
|
|
1 |
|
|
|
|
|
2.8029 |
2.8539 |
2.9444 |
|
|
|
|
|
1.5 |
|
|
|
|
|
2.8231 |
2.8730 |
2.9624 |
|
|
|
|
|
|
3 |
|
|
|
|
2.8022 |
2.8520 |
2.9409 |
|
|
|
|
|
|
4 |
|
|
|
|
2.8138 |
2.8624 |
2.9497 |
|
|
|
|
|
|
|
1 |
|
|
|
2.7788 |
2.8310 |
2.9225 |
|
|
|
|
|
|
|
1.5 |
|
|
|
2.7742 |
2.8263 |
2.9175 |
|
|
|
|
|
|
|
|
2 |
|
|
2.7789 |
2.8312 |
2.9227 |
|
|
|
|
|
|
|
|
3 |
|
|
2.7745 |
2.8266 |
2.918 |
|
|
|
|
|
|
|
|
|
1 |
|
2.1384 |
2.1875 |
2.2688 |
|
|
|
|
|
|
|
|
|
1.5 |
|
1.8064 |
1.8479 |
1.9173 |
|
|
|
|
|
|
|
|
|
|
1 |
2.7337 |
2.7747 |
2.8518 |
|
|
|
|
|
|
|
|
|
|
1.5 |
2.6694 |
2.7028 |
2.7702 |
nanoparticle mass fraction. From Table 2, it is clear that increases in M, R, Sc, and Kc have been found to enhance the primary skin friction. On the other hand, the Primary skin friction reduces with increases in K, Gr, Gm, Ec, Pr, Re, and Rem. Secondary skin friction is enhanced by increases in M, K, Gr, Gm, R, Ec and Pr while it is diminished by increases in Sc, Kc, Re, and Rem. It is also observed that both the primary and secondary skin friction coefficient increases with an increase in the nanoparticle mass fraction
. Table 2 and Table 3 represent the skin friction for various leading parameter variations. This involved setting up the values R = 1, K = 0.5, Gr = 10, Gm = 10, M = 10, Pr = 0.71, Sc = 0.5, Je = 2, Kc = 1, Rem = 0.5, and Re = 0.5 which were then scattered around the range for computing reasons.
5. Conclusion
In this investigation, we carried out attributes of the convective flow of nanofluid flow with a vertical rotating plate imbedded in a porous medium by incorporating-induced magnetic field influence. The current model of the flow is utilized to stimulate the physical phenomena of the flow that occur inside the washing machine of the new generation, which is one of the frequent engineering utilizations in designing engineering products. The findings acquired are crucial in improving the performance of the new generation washing machine by adjusting the physical parameter values to get the best performance out of it. This paper uses the BI-SRM technique to numerically examine the effects of temperature-dependent viscosity, an induced magnetic field, nanofluid natural convection flow with heat and mass transfer, and a first-order chemical reaction through a semi-infinite, vertical rotating plate embedded in a porous media. Analysis of the momentum equations and induction equations, which incorporate variable viscosity and variable magnetic field, is the major objective. The problem investigates the primary physical parameters, which are then utilized in the governing equations. The most significant engineering outcomes are highlighted and well explained. The key findings of the current research lead to the following significant conclusions:
In case of varying the nanoparticle volume fraction, increases in Gr, Gm, K, and Re cause both the velocity components to accelerate, whereas increases in M, Kc, and Sc cause both velocity components to slow down.
An increase in Re causes a reduction in fluid temperature, while Re has a reverse tendency in the concentration profile.
Primary velocity slowed when R was raised and accelerated when Rem was increased.
Primary velocity increased when Pr increased in cases of low nanoparticle volume fraction concentration, but had the opposite tendency in cases of high nanoparticle volume fraction concentration.
Secondary velocity accelerates with an increase in R, whilst an increase of Rem decelerates the secondary velocity.
The axial components of the induced magnetic field decline as Rem increases.
When M, R, Sc, and Kc increase, primary skin friction
increases, but when K, Gr, Gm, Ec, Pr, Re, and Rem increase,
decrease.
Increases in M, K, Gr, Gm, R, Ec, and Pr were linked to increased secondary skin friction
, whereas increases in Sc, Kc, Re, and Rem were linked to decreased secondary skin friction
.
In both primary and secondary velocities, an increase in nanoparticle volume fraction
tends to increase skin friction.
In terms of future study, the authors recommend incorporating cross-diffusion effects, which are common in many real applications. Entropy generation analysis may also be used to improve the thermodynamic optimization of the model’s system.
Data Availability
All required data is available in the text of the paper.
Acknowledgments
The African Union and the Pan African University’s Institute for Basic Sciences, Technology, and Innovation (PAUSTI) were responsible for supporting this work, which the first author would like to acknowledge.