Abstract

Radiation is an important branch of thermal engineering which includes geophysical thermal insulation, ground water pollution, food processing, cooling of electronic components, oil recovery processes etc. An analysis of unsteady magneto-convective heat-mass transport by micropolar binary mixture of fluid passing a continuous permeable surface with thermal radiation effect has been introduced in this paper. The governing equations are transformed into coupled ordinary differential equations along with Boussinesq approximation by imposing the similarity analysis. Applying the shooting technique, the obtained non-linear coupled similarity equations are solved numerically with the help of “ODE45 MATLAB” software. The results of the numerical solutions to the problem involving velocity, temperature, concentration and micro-rotation are presented graphically for different dimensionless parameters and numbers encountered. With an increase of suction parameter, the velocity distributions very closed to the inclined permeable wall decrease slightly where . But for the uplifting values of sunction, both micro-rotation profile and species concentration enhance through the boundary layer. The skin-friction coefficient increases about 61%, 13%, 27% for rising values of Prandtl number (0.71-7), radiation effect (0 - 1) and thermal Grashof number (5-10), respectively, but an adverse effect is observed for magnetic field (1 - 4), inclined angle and Schmidt number (0.22 - 0.75). Heat transfer and mass transfer reduce about 82%, 53%, respectively, in increasing of Pr (0.71-1) and 36%, 11%, respectively, in increasing of thermal radiation (0 - 1). The surface couple stress increases about 26%, 49%, 64% and 30% with the increasing values of magnetic field (1-4), inclination angle , suction (0-1) and Schmidt number (0.22-0.75), respectively. Finally, the present study has been compared with the earlier published results. It is observed that the comparison bears a good agreement.

1. Introduction

Some non-Newtonian fluids are common to flow along a vertical surface in industrial engineering processes. These non-Newtonian fluids include molten polymers, fossil fuels, animal blood, pulps and fluids which contain certain additives, etc. which are occasionally available in industrial production. Earlier the researchers of fluid dynamics have played a great role in explaining the flow of these non-Newtonian fluids through the stretching surfaces in vertical as well as horizontal positions. In the case of continuous moving surfaces, many practical applications of boundary layer flow are remarkable in industrial and engineering processes. The boundary layer along a liquid film in condensation processes, crystal growth, the moving of heat-treated material between feed and wind-up rolls, aerodynamic extrusion of plastic sheets, many metallurgical processes, such as drawing of continuous filaments through steady fluids, cooling of an infinite metallic plate in a cooling path, which may be an electrolyte, etc., are some notable examples of such continuous moving surfaces. The properties of steady natural convective fluid flow, heat and mass transfer passing an isothermal vertical surface embedded in a permeable medium were predicted by the non-Newtonian fluid flow model [14]. Considering heat absorption in a permeable medium, the effects of chemical reaction were discussed by Rao et al. [5]. Additionally, they studied the effects of a semi-infinite vertical surface on an unstable magnetohydrodynamic free convection fluid flow. Ali et al. [6] solved the governing equations for an unsteady MHD free convective-viscous fluid flow along an inclined surface taking changeable heat and mass transport. They considered optically thick fluid instead of optically thin and the surface makes a fixed angle with the vertical axis. Implications of variable viscosity and chemical reaction on MHD natural convective radiating flow through an inclined porous plate was investigated by Ali et al. [7]. The MHD flow of a micropolar fluid was studied over a vertical plate with heat generation and suction by Ferdous et al. [8] which leads to look for similarity solutions of the problem by considering the variable wall temperature. In recent time, Arifuzzaman et al. [9] investigated chemically reactive and naturally convective unsteady micropolar fluid flow including nanoparticles passing a vertical permeable surface with mass diffusion, radiation absorption, thermal radiation, heat source and MHD. The familiar boundary layer approximations helped them to develop a flow model.

Magneto-hydrodynamics (MHD) deals with the behavior of fluids which are electrically conducting. It also deals with the properties of magnetism. Some mentionable examples of such magnetofluids are liquid metals, salt water, electrolytes, Plasmas etc. Maxwell’s equations of electromagnetism and Navier–Stokes equations of fluid dynamics jointly represent the MHD. Many authors are interested to study MHD for its applications to engineering and geophysics. In presence of suction, Al Mahbub et al. [10] introduced a mathematical model to analyze the theory of micro-rotation of MHD fluids over a stretching plate. Ali et al. [11] found similarity solution of time dependent MHD border layer heat transport and fluid movement over a moving wedge applying the Buongiorno’s nanofluid model. The impacts of chemical reaction on the MHD nanofluid over a vertical surface embedded in a permeable medium were presented by Swain et al. [12]. Through the study of Alam and Huda [13], a novel concept was created due to local similarity solutions of an unsteady MHD free convective heat transfer flow. They took into account fluid movement via a porous flat plate without taking internal heat generation into account. In their extension of earlier research, Alam et al. [14] considered the impact of internal heat generation or absorption for free convective heat transfer flow crossing a vertical permeable flat surface.

MHD border layer nanofluid movement and thermal transport through a stretching surface was analyzed by Ali et al. [15]. The effects of chemical reaction on an unsteady MHD free convection fluid flow were demonstrated by Reddy et al. [16]. The flow through a semi-infinite vertical surface embedded in a permeable material with heat absorption was taken into consideration. Thakur and Hazarika [17] showed the influences of thermal conductivity and variable viscosity on unsteady free convective heat - mass transport MHD flow. They considered the micropolar fluid passing a permeable surface in presence constant heat flux. They also considered viscosity and thermal conductivity as the inverse linear functions of temperature. The term “micropolar fluids” refers to complicated fluids having microstructure. The applications of micropolar fluids are extensively found in the field of geophysics, biomechanics, energy systems and industrial engineering etc. In presence of non-uniform heat source and chemical reaction, a numerical study of MHD micropolar fluid was presented by [1820] along a stretching sheet embedded in permeable surface. Applying micropolar fluids model and hydrodynamic characteristics with default function in FLUENT was analyzed by Wang et al. [21]. Considering the heat transfer by radiative influence, Baruah and Hazarika [22] discussed the impacts of variable viscosity as well as thermal conductivity on unsteady micropolar fluid around a porous cylinder. Haque [23] deeply studied the impacts of induced magnetic field on the transient heat-mass transport and MHD micropolar fluid flow. He thought of fluid flow crossing a vertical surface that is semi-infinite and surrounded by a permeable medium and a steady heat sink.

The migration of tiny particles from a hot surface to a cool one is known as thermophoresis. It plays a great role to micro-particle transport caused by surrounding temperature gradient. The application of thermophoresis is mostly observed in aerosol technology. Darcy mixed convective heat and mass transfer that travel through a permeable wedge along with chemical reaction have been explored previously for the effects of thermophoretic influences. By passing a semi-infinite porous inclined surface with chemical reaction, Mondal et al. [24] observed the Soret-Dufour influences on MHD convective flow. They also observed the effects of thermophoresis in the same environment. Thermal radiation is a process which makes the variation of temperatures in order to exchange energy. Radiation is transmitted from ultraviolet into far-field infrared according as the material temperature. It is an important branch of thermal engineering which includes geophysical thermal insulation, ground water pollution, food processing, cooling of electronic components, oil recovery processes etc. Khader and Megahed [25] presented a numerical simulation applying FDM (Finite Difference Method). They studied the heat transfer and fluid flow along an unsteady stretching sheet embedded in a porous surface. Shamshuddin et al. [26] propagated the research work of Babu et al. [27] in presence of thermal radiation as well as heat source/sink. The consequences of first order chemical reaction and FDM (Finite Difference Method) numerical solutions were also presented for radiative-convection micropolar fluid for an inclined permeable surface. Nazir et al. [28] considered the effects of viscous dissipation, Joule heating, and heat generation on the transport of momentum in a rotating cone under generalized ohm’s law and heat. Also, the heat and mass transport past a stretched surface having pores in a pseudo-plastic model were investigated by Wang et al. [29]. They used the MAPLE software in their simulations. Bhatti et al. [30] investigated the electro-magnetohydrodynamic (EMHD) dissipative free convection in a micro-channel containing a permeable medium saturated with a viscoelastic fluid. Rasool et al. [31] introduced the concept of entropy optimization in fluid flow through a Darcy channel with a non-linearly stretching surface. Colak et al. [32] demonstrated the efficiency of the ANN procedure for Buongiorno’s mathematical model using Powell Eyring fluid. This included slip, Darcy-Forchheimer, heat source phenomena, and convective boundary conditions toward a radiative stretching surface.

The innovation of this research work is to explore the thermal radiation effect of unsteady magneto-convective heat and mass transfer by micropolar binary mixture of fluid passing a continuous permeable surface. The main novelty of this research is to compare our numerical results with the published results. We have also enhanced this paper further by considering the induced magnetic force with the inclined angle which has not been explained yet. The current research is to extend the works of Animasaun [33] by considering the effects of the magnetic force parameter and the inclined angle. Three surfaces (horizontal, vertical and inclined surfaces) have been considered to study three positions of the permeable surface at a time to compare the radiative impacts. In this study, we have carefully described the velocity, temperature, concentration and microrotation profiles. The impacts of the thermo-physical parameters on the flow and heat-mass transport characteristics are discussed in detail. Additionally, tabular forms are supplied for the thermophysical parameters of practical importance, including the Nusselt number, Sherwood number, surface couple stress, and coefficient of skin friction.

2. Mathematical Analysis

Figure 1 displays the considered physical geometry. A micropolar binary fluid mixture traveling over a continuous permeable surface is modeled for an unsteady magneto-convective heat-mass transport. The flow is taken into consideration across the x-axis. This is also expected along the semi-infinite surface and y-axis is normal to it. An optically thin fluid is considered with absorption coefficient (α) satisfying the condition α < <1. According to Cheng [34], the radiative heat flux is approximately represented by the fourth power of temperature in the energy balance equation. The temperature gradient normal to the surface is significantly higher in boundary layer theory than the temperature gradient along the surface, i.e.

For this reason, the thermophoretic velocity component along the surface is ignorable with comparison to velocity component normal to the surface. It is expected that the wall temperature is less than the free stream temperature . With the help of Mathematical model presented by Sastry et al. [35], all presumptions expressed above and Boussinesq’s approximation; the governing equations can be written as:

Here a uniform transverse magnetic field is set along y-axis i.e., along horizontal direction. Also, the induced magnetic field is ignored which makes magnetic Reynolds as small. Further, there exists no electric field and the Hall effect of MHDs is disregarded.

Here v can be considered either as a function of time or as constant using continuity equation (2). Makinde [36] has defined the velocity component along y-axis as follows- where c <0 indicates the injection parameter and c >0 represents the suction parameter. Abdel-Rahman [37] and Qasim et al. [38] have introduced micro-inertia per unit mass and spin gradient viscosity as follows:

Talbot et al. [39] also introduced the thermophoretic parameter in equation (5) and later by Tsai [40] as- where KTh represents the thermophoretic coefficient. Batchelor and Shen [41] considered the values of KTh from 0.2 to 1.2. The fluid properties are considered to be same aside from dynamic viscosity and the effect of density difference with concentration and temperature in the body force term (Boussinesq’s approximation).

In 1903, Joseph Boussinesq considered buoyancy-driven flow; he investigated acceleration due to gravity and density differences corresponding to buoyancy driven flow. The density model for slight temperature variation between the wall and free stream layer is given by- where is known as buoyancy term. Applying the above approximation, the buoyancy term is given by-

The buoyancy force as well as pressure term is now expressed as . Furthermore, Boussinesq approximation for combined heat and mass transfer convection can be expressed as follows:

This principle is known as Boussinesq approximation. According to Boussinesq [42], this approximation indicates the density variation is too small to be ignored.

Then buoyancy model (11) can be modified as-

Batchelor [43] has studied the temperature dependent viscosity model and later on it has been accepted by Mukhopadhyay [44]. Animasaun and Anselm [45] have developed which is valid only for . This was also changed to the form-

This satisfies the condition

Following Raptis [46], the radiative heat flux is formulated as-

Here where indicates mean absorption coefficient as well as represents Stefan-Boltzmann constant. The temperature variations within the flow have been considered as small through the study of Raptis [47]. For this reason, the fluid flow can be written as a linear function of temperature. Neglecting higher order terms in application of Taylor series in about , we get

Williams Landel-Ferry model, Masuko-Magill model, Arrhenius model, exponential model and Batchelor model are the examples of some models for shear viscosity. The analysts in fluid dynamics are yet practicing on empirical data which reveals the impact of temperature on both vortex and dynamic viscosity of micropolar fluid. Using the relations (13) and (14) in equation (3), we obtain-

(3) implies that,

Using the equation (13), the equation (5) becomes-

Using the equation (8), the equation (6) becomes-

Equations (17)–(20) are subjected to the following boundary conditions:

The above boundary conditions are valid when and . In Equation (22) when =0, we have N (0, t) =0 which indicates no-spin condition. This means that the microelements in a concentrated substance of the fluid adjacent to the wall as if flows cannot spin.

We introduce the following dimensionless vinary differential equations as foll-(19):

Finally, we have the dimensionless non-linear ordinary differential equations as follows:

Also, the following dimensionless boundary conditions are obtained for the above problem: where is the dimensionless velocity function, is the dimensionless temperature function, is the dimensionless concentration function, is the dimensionless micro-rotation function, is the thermal Grashof number, is the variable viscosity parameter, is the micro-rotation parameter, is known as the solutal Grashof number, is the radiation parameter, is the Prandtl number, is defined as the time dependent micro-rotation parameter, is the thermophoretic parameter, is the magnetic field parameter, is the Schmidt number, is the inclined angle of permeable surface with -axis and η is a similarity variable.

3. Numerical Results and Discussion

An analysis of unsteady magneto-convective heat-mass transport by micropolar binary mixture of fluid passing a continuous permeable surface with thermal radiation effect has been introduced in this research work. Firstly, the higher order nonlinear partial differential equations (2)–(6) has been transformed into second order simultaneous linear ordinary differential equations using modified Boussinesq’s approximation. Applying the similarity technique, then it has also been converted into initial value problem. Finally, we have solved numerically the coupled non-linear ordinary differential equations (25)–(28) with boundary conditions (29)–(30). Here we have used the shooting technique through “ODE45 MATLAB” software.

Applying the aforesaid numerical procedure, the non-dimensional velocity, temperature, concentration and microrotation profiles are determined against the dimensionless coordinate η with the changes of different thermophysical parameters controlling the phenomena of fluid flow. The numerical results are displayed through graphs 2 to 10. The thermo-physical quantities of practical interest such as Nusselt number, coefficient of skin friction and Sherwood number are also calculated and are displayed in Tables 17 for analyzing the internal properties of the fluid flow. The comparison between present numerical results and a published paper has been revealed in Tables 8 and 9. We have included here the parameters or numbers as magnetic parameter (M), thermal Grashof number (Gr), Prandtl number (Pr), solutal Grashof number (Gc), Schmidt number (Sc), thermophoretic parameter (), thermal radiation parameter (𝑅), micro-rotation parameter (), temperature dependent variable viscosity parameter (ξ), time dependent micro-rotation parameter (), suction parameter (c) and an inclined angle (). Here it is noted that we have chosen the changes of some selected numbers or parameters while overlooking the impact of one number or parameter on the field variables. The rest of the numbers or parameters are taken as constants. We have taken the fixed values of parameters or numbers as: , Pr = 0.71, R = 0.5, M = 2, c = 0.5,  = 0.33, = 1, Sc = 0.22, ξ = 3,  = 0.5, Gr = 10, Gc = 10. Furthermore, the boundary conditions have been considered at infinity at a finite point of η = 5.

3.1. Effect of Magnetic Field

The impacts of magnetic field (M) on fluid velocity and micro-rotation of fluid particles for various values of M have been presented in Figures 2(a) and 2(b). We observe from the figures that the increasing M leads to decrease the velocity fields and also micro-rotation profiles. Because a transversely applied magnetic field generates a Lorentz force. This force acts like a drag force which reduces the fluid motion. Therefore, as the values of M increase, resistive force is also increased resulting a reduction of the velocity and micro-rotation profiles.

3.2. Effect of Prandtl Number

The effects of Prandtl number (𝑃𝑟) on the velocity, temperature and micro-rotation profiles are presented in Figures 3(a)3(c). The graph 3(a) depicts that with the rising value of Pr, the velocity reduces. Physically, the kinematic viscosity of the fluid promotes owing to increasing values of Prandtl number. This makes the fluid much thicker which reduces the fluid velocity. Figure 3(b) shows that for rising values of Prandtl number, the temperature diminishes. Physically, a higher Prandtl number contains relatively low thermal conductivity. This leads to reduce heat conduction and as a result temperature decreases. Hence for the increasing values of 𝑃𝑟, the rate of heat transfer enhances which reduces the temperature profiles. But an opposite trend is noticed for the micro-rotation profiles with the variation of Prandtl number. The micro-rotation profiles are found to increase noticeably adjacent to the surface 0 ≤ η ≤ 2.2 and further than its effect is negligible as η increases.

3.3. Effects of Thermal Grashof Number and Solutal Grashof Number

The characteristics of thermal Grashof number and solutal Grashof number on velocity and micro-rotation profiles are presented in Figures 4(a)4(d). We observed from the Figures 4(a) and 4(b) that the infuences of thermal Grashof number and solutal Grashof number on velocity are identical. For increasing values of Gr and Gc, velocity enhances. But a small variety is noticed for the maximum value of f(η) which are found as 4.25 for Gr and 2.75 for Gc at η = 0.8. Physically, it is well established that the increase of buoyancy parameters (Gr and Gc) enhances the fluid flow. It is also seen from both the figures that thermal Grashof number and solutal Grashof number have no influence as the fluid move away from the bounding surface where 3.2 ≤ η ≤ 5.

The impacts of local thermal and solutal Grashof number on the microrotation distributions are found same in Figures 4(c) and 4(d). As they increase from zero to positive values, the microrotation of the fluid particles decreases from negative side near the wall in the region 0 ≤ η ≤ 2.4. Thereafter, there is no influence of Gr and Gc on the microrotation of fluid particles and finally, asymptotically approach 0 as η → 5. Hence the micro-rotation profiles are highly influenced by fluids within the domain. The behaviors of local thermal and solutal Grashof number on the microrotation distributions are found same in Figures 4(c) and 4(d). As they increase from zero to positive values, the microrotation of the fluid particles decreases from negative side near the wall in the region 0 ≤ η ≤ 2.4. Thereafter, there is no influence of Gr and Gc on the microrotation of fluid particles and finally, asymptotically approach 0 as η → 5.

3.4. Effect of Schmidt Number

Figures 5(a)5(c) describe the velocity, concentration and micro-rotation profiles for several values of the Schmidt number. It is observed from the Figures 5(a)5(b) that velocity and concentration reduce with an increase of Sc value. Since Sc defines the ratio of kinematic viscosity (momentum diffusivity) with the diffusivity of the particle (mass diffusivity), so an increase of Sc corresponds to a reduction of mass diffusivity where the viscosity of the medium enhances and consequently, the velocity of the fluid is reduced. From the figure it is seen that for increasing values of Sc, concentration boundary layer becomes thinner. This leads to reduce the concentration profiles. In a physical manner, molecular diffusivity decreases due to an increase in magnitude of Sc. For that reason, the species concentration is lower for large values of Schmidt number and higher for small values of the Schmidt number. Figure 5(c) is sketched to observe the effects of Schmidt number (Sc) on the micro-rotation profiles. The higher value of Sc corresponds to a big increase of the micro-rotation profiles in a closed proximity to the wall 0 ≤ η ≤ 2.2 and thereafter no change of the behavior is seen and finally approach zero asymptotically all together.

3.5. Effect of Thermophoresis

Figures 6(a)6(c) illustrate the velocity, concentration, and micro-rotation profiles for various values of thermophoretic parameter . Since the increasing value of thermophoretic parameter reduces the momentum boundary layer thickness and hence the velocity diminishes. An opposite behavior is observed for concentration and micro-rotation of the fluid particles as increases. The fluid velocity increases rapidly due to a small change in the thermophoresis parameter. This creates overplus heat energy which leads to rapid increase in the concentration distribution. We also overlooked that an enhancement in the thermophoretic parameter holds relatively large impacts on the micro-rotation of fluid particles adjacent to the wall (0 ≤ η ≤ 2.2). After that, like Sc the micro-rotation profiles have no effect and tend to zero asymptotically when η → 5.

3.6. Effect of Thermal Radiation

Figures 7(a)7(c) describe the impact of thermal radiation parameter (R) on the velocity, temperature, and micro-rotation profiles. Since the radiation parameter is the ratio of conduction heat transfer and thermal radiation transfer, rising values of R enhance temperature. As a result, fluid viscosity decreases and velocity upturns within the boundary layer.

Contrary, the thermal boundary layer always found to thicken in presence of radiation, since radiation gives an additional means to disseminate energy. Consequently, the temperature lessens with increasing values of radiation parameter (R). We also noticed that temperature increases for uplifting values of η with radiation parameter. In the range of 0 ≤ η ≤ 2.6, the micro-rotation profiles decrease with increasing R and after that no effect is observed with increasing η and asymptotically approaches zero as η → 5.

3.7. Effect of Variable Viscosity

Figures 8(a) and 8(b) display the velocity and micro-rotation profiles for several values of temperature dependent variable viscosity parameter (ξ). As ξ enhances, the velocity profiles very close to the wall are negligible where 0 ≤ η ≤0.3 and thereafter tend to decrease (0.3 ≤ η ≤ 1.2) and finally approach 1 asymptotically as η → 5. From Figure 8(b), it is obvious that the micro-rotation of the fluid particles rises averagely with the uplifting value of ξ, very near to the wall 0 ≤ η ≤ 0.1; but an alternate behavior is observed when 0.1 ≤ η ≤ 1.4 thereafter. Here the micro-rotation profiles diminish slightly with the upturning value of ξ and lead to zero asymptotically as η → 5.

3.8. Effect of Suction

From Figures 9(a)9(c), we observed the impacts of suction parameter (c) on the velocity of the fluids, the chemical species concentration and micro-rotation field. As we expect, if the temperature at the free stream is less than the wall temperature, then rising values of suction parameter (c) correspond to a reduction of fluid velocity from the wall to the free stream (for the case of cooling of the surface, i.e. for 𝐺𝑟 > 0). This is because of suction removes the low momentum fluid around the heated wall. For this reason, both separation and transition are delayed. Since the fluid adjacent to the heated wall (η = 0) is pushed towards the region where the buoyancy forces reduce the fluid motion owing to the large effect of viscosity. We also observed that with an enhancement of suction (c >0), the velocity distributions very closed to the inclined permeable wall decreased negligibility where 0 ≤ η ≤0.3. We further noticed that fluid velocity is found to increase with η and are attained to its maximum peak within 0.5 ≤ η ≤ 1.4 and after that tend to 1 asymptotically as η → 5. From the Figure 9(b), it is noticed that an increase in the magnitude of c, the species concentration through the boundary layer also increased. As the species transfer from lower concentration region to the region of higher concentration, increase in suction will remove thin boundary layer passing the inclined permeable surface and as a result, the high dense fluid will concentrate near the surface quickly. We also investigated that the micro-rotation of fluid particles enhances for rising value of c, where 0 ≤ η ≤ 2.6. Thereafter (2.6 ≤ η ≤ 5) the micro-rotation profiles are not influenced with the uplifting value of c.

3.9. Effect of Inclination

Figures 10(a) and 10(b) depict the impact of inclined angle () on the velocity and micro-rotation profiles.

We investigated that the velocity leads to reduce for the uplifting values of inclined angle (). For the increasing values of , the impact of the buoyancy force reduces due to the multiplication factor cos and consequently the velocity decreases rapidly for the increasing value of . For horizontal surface ( = ), we also see that velocity profile is a strictly monotonic increasing function of . But an opposed behavior is observed from the Figure 10(b) for the increasing values of . Because for increasing value of , the angular momentum increases within the boundary layer which leads to enhance the micro-rotation profile.

4. Skin Friction, Nusselt Number, Sherwood Number and Surface Couple Stress

The thermo-physical quantities of practical importance such as skin friction coefficient which corresponds to f/(0), Nusselt number (Nu) which corresponds to 𝜃/(0), Sherwood number (Sh) which corresponds to 𝜑/(0) and the surface couple stress h/(0) are also tabulated to explain the internal behavior of the fluid flow:

Table 1 displays how different Prandtl number values affect the types of skin-friction coefficient, Nusselt number, Sherwood number, and surface couple stress. Here, it is noted that only the skin-friction coefficient increases with growing Pr values, whereas all other thermophysical parameters (including the Nusselt number, Sherwood number, and surface couple stress) decrease as Pr values rise. The skin-friction coefficient increases about 61% for uplifting values of the Prandtl number (0.71-7). But Nusselt number, Sherwood number and surface couple stress reduce about 82%, 53% and 3817%, respectively, in increasing of Pr (0.71-1).

Table 2 displays the features of thermal Grashof number (Gr) on skin-friction coefficient, Nusselt number, Sherwood number and. We observe that the coefficient of skin-friction increases in case of mode of cooling (Gr>0) and reduces in case of mode of heating (Gr<0) but a reverse effect is noticed for the surface couple stress. It is also evident that no effect of Gr is seen on Nusselt number and Sherwood number in both cases (Gr>0 and Gr<0). It is also revealed that the skin-friction coefficient upturns about 27% whereas the surface couple stress diminishes about 26% with the increasing values of Gr (5-10).

The effects of the radiation parameter (R) on the surface couple stress, Nusselt number, Sherwood number, and skin-friction coefficient are shown in Table 3. Here, it is discovered that only the skin-friction coefficient increases with growing values of R, but all other thermo-physical variables decrease with rising R values. The skin-friction coefficient promotes about 13% for uplifting values of the radiation parameter (0 - 1). But Nusselt number, Sherwood number and the surface couple stress reduces about 36%, 11% and 126%, respectively, in increasing of R (0-1).

Table 4 illustrates how the magnetic parameter (M) affects the skin-friction coefficient, Nusselt number, Sherwood number, and surface couple stress. Here we see that the coefficient of skin-friction reduces for rising values of M but an adverse effect is observed for the surface couple stress. It is also evident that there is no influence of M on Nusselt number and Sherwood number for uplifting values of M. We also reveal that the skin-friction coefficient decreases by about 22% whereas the surface couple stress increases by about 26% with the increasing values of M (1.0-4.0).

The features of the inclined angle () on the surface couple stress, Nusselt number, and Sherwood number are shown in Table 5. We discovered that when values of increase, the coefficient of skin-friction decreases, while the surface couple stress has the opposite effect. It is also evident that the upturning values of do not affect the Nusselt number or the Sherwood number. Additionally, we observe that as the values of increase, the skin-friction coefficient decreases by around 49% while the surface couple stress increases by about 49%.

Table 6 details how the suction parameter (c) affects the skin-friction coefficient, Nusselt number, Sherwood number, and surface couple stress. We see that when the values of c rise, all physical quantities get better. Additionally, we observe that with increasing values of c (0-1), the skin-friction coefficient, Nusselt number, Sherwood number, and surface couple stress all rise by around 2%, 75%, 14%, and 64%, respectively.

The effects of the Schmidt number (Sc) on the skin-friction coefficient, the Nusselt number, the Sherwood number, and the surface couple stress are shown in Table 7. From Table 7, we are clear that the coefficient of skin-friction reduces with an increase in the magnitude of Sc whereas reverse impacts are observed for both Sherwood number and surface couple stress. It is also evident that there is no influence of Sc only on Nusselt number for increasing values of Sc. We also reveal that the skin-friction coefficient decreases by about 24% but the Sherwood number and surface couple stress increase by about 149% and 30%, respectively, with increasing values of Sc (0.22-0.75).

5. Comparison

We have compared the numerical results of the present study with Animasaun [33]. Tables 8 and 9 show the comparison of the Nusselt number for thermal radiation parameter (R) and the Sherwood number for Schmidt number (Sc), respectively, considering present numerical results with mentioned published paper. The current numerical results shown in these Tables 8 and 9 are found to be good in agreement with that of Animasaun [33].

6. Conclusion

The insertion of a continuous permeable surface significantly influences the motion and micro-rotation of the flow field. The fluid velocity enhances within the boundary layer whereas both temperature and micro-rotation of fluid particles reduce for rising values of thermal radiation. The velocity of fluid within the boundary layer increases as the radiation rises, whereas the temperature and the microrotation of fluid particles decrease as the radiation rises. The main findings are as follows: (i)With the movement of a continuous permeable surface, the velocity decreases, and the angular momentum increases.(ii)For the increasing value of thermophoretic parameter (), the velocity reduces but the species concentration of the fluid and the micro-rotation of particles enhance.(iii)The micro-rotation of fluid particles increases on average with the increase of ξ, very near to the wall 0 ≤ η ≤ 0.1. But thereafter, an alternate behavior is observed when 0.1 ≤ η ≤ 1.4.(iv)The velocity profile very close to the wall is negligible where 0 ≤ η ≤0.3 and thereafter tends to decrease (0.3 ≤ η ≤ 1.2) for increasing value of ξ.(v)The increase of buoyancy parameters (Gr and Gc) enhances the fluid flow. But the microrotation of the fluid particles decreases significantly from the negative side near the wall in the region 0 2.4.(vi)The velocity and temperature reduction for the rising value of Prandtl number, but the micro-rotation profiles are found to increase significantly adjacent to the surface (0 2.2).(vii)The velocity and species concentration reduce, and the micro-rotation profiles enhance with the uplifting value of Schmidt number (Sc).(viii)An increase in radiation parameter (R) leads to an increase in fluid velocity within the boundary layer whereas both temperature and micro-rotation of fluid particles reduce for upturning values of radiation parameter (R).(ix)With an enhancement of suction parameter (c >0), the velocity distributions very close to the inclined permeable wall decrease slightly where 0 ≤ η ≤0.3. But for the rising values of c, both micro-rotation profile and species concentration enhance through the boundary layer.(x)Suction parameter, radiation parameter, thermophoretic parameter, Schmidt number, and Prandtl number highly affect the fluid properties.(xi)The skin friction coefficient enhances but all other thermo-physical quantities reduce with increasing values of both Pr and R.(xii)In cases of cooling (Gr>0) and heating (Gr<0), the coefficient of skin-friction and the surface couple stress play opposing roles.(xiii)Only the coefficient of skin-friction reduces for rising values of M but an adverse effect is observed for the surface couple stress.(xiv)With changes in inclined angle (), neither heat transfer nor mass transfer is affected.(xv)Suction parameter improvements have a positive impact on all physical quantities.(xvi)The increasing values of the Schmidt number (Sc) result in a drop in the skin-friction coefficient, but an increase in mass transfer and surface couple stress.

Nomenclature

MHD:magneto-hydrodynamics
Sc:Schmidt number
:density of current
:free stream temperature
:volumetric expansion coefficient with temperature
:fluid concentration, kg m-3
:free stream concentration
:velocity component along the x-axis, ms-1
:uniform surface velocity
Gr:thermal Grashof number
λ:thermophoretic parameter
τ:vortex viscosity
:specific heat at constant pressure
:thermal diffusion ratio
:similarity parameter
:acceleration due to gravity, ms-2
:magnetic parameter
ξ:Variable viscosity parameter
:local skin friction coefficient
:Sherwood number
:dimensionless temperature
:radiative heat flux
:thermal radiation parameter
:suction parameter
:uniform magnetic field, Am-1
:dimensionless concentration
:temperature of fluid, k-1
:wall temperature, k-1
:volumetric expansion coefficient with concentration
:fluid density, kg m-3
:wall concentration, kg m-3
:velocity component along the y-axis, ms-1
:suction velocity
:solutal Grashof number
:concentration susceptibility
η:similarity variable
:fluid mean temperature
:mass diffusivity coefficient
:kinematic viscosity, m-2s-1
:suction and blowing
:Prandtl number
t:Time, s
:Nusselt number
:dimensionless velocity
:Chemical reaction rate of concentration
к:thermal conductivity
:chemical reaction parameter
µ:Dynamic viscosity.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The author (s) declare (s) that they have no conflicts of interest.

Acknowledgments

The research work is done in the Department of Mathematics, Bangladesh University of Engineering and Technology (BUET), Dhaka, Bangladesh. This article is a part of the first author’s PH.D. thesis.