Unsteady Micropolar Fluid over a Permeable Curved Stretching Shrinking Surface

1Department of Mathematics, Faculty of Science, Universiti Putra Malaysia (UPM), 43400 Serdang, Selangor, Malaysia 2Institute For Mathematical Research, Universiti Putra Malaysia (UPM), 43400 Serdang, Selangor, Malaysia 3School of Mathematical Sciences, Faculty of Science and Technology, Universiti Kebangsaan Malaysia (UKM), 43600 Bangi, Selangor, Malaysia 4Department of Mathematics, Babeş-Bolyai University, 400084 Cluj-Napoca, Romania


Introduction
The fundamental nature of the theory of micropolar fluid flow lies in the extension of the constitutive equations for Newtonian fluids, so that more complex fluid such as particle suspensions, liquid crystals, animal blood, and turbulent shear flows can be expressed by this theory.The micropolar fluid theory which takes into account the inertial characteristics of the substructure particles which are permitted to endure rotation has been intended by Eringen [1] and was further discussed by Eringen [2].This theory has caused much attention and many aspects of flows being reexamined to verify the effect of the fluid microstructure.A detailed study on the applications of microcontinuum fluid mechanics has been accessed by Ariman et al. [3].Since then, many other studies related to the ideas of micropolar fluid have been done.Ahmadi [4] analyzed self-similar solution for incompressible micropolar flow over semi-infinite plate.A study on stagnation flows of micropolar fluids with strong and weak interactions was reported by Guram and Smith [5].Then, it seems Gorla [6] has obtained numerical results for the micropolar boundary layer flow at a stagnation point on moving wall using fourth-order Runge Kutta method.
During the past several years considerable interest has been verified in the study of steady and unsteady flows of a viscous incompressible fluid determined by a linearly stretching surface through a quiescent fluid.Such flow situations are encountered in a number of manufacturing practices, such as the cooling of metallic plate in a cooling bath, polymer sheet extrusion, and heat-treated materials that are on a conveyor belt.Crane [7] found a closedform exact solution for a two-dimensional laminar flow of an incompressible viscous fluid over a linearly stretching sheet, and Wang [8] discovered similarity solutions for the axisymmetric case.In another paper, Wang [9] has analyzed the viscous flow due to a stretching sheet with surface slip and suction.The combination of both stagnation point-flows past a stretching surface was solved numerically by Mahapatra and Gupta [10,11] and analytically by Mahapatra et al. [12].The effects of radiation with magnetic field on stretching sheet were obtained by Mat Yasin et al. [13].Also, in a very interesting recently published paper, Turkyilmazoglu [14] has 2 Mathematical Problems in Engineering mathematically resolved the domain of existence for magnetohydrodynamic mixed convection flow of a micropolar fluid past a heated or cooled stretching permeable surface by taking into account the heat generation and absorption effects.Furthermore, few studies have been done on the reaction of permeable stretching sheet in nanofluid which was conducted by Ibrahim and Shankar [15], Bachok et al. [16], and Zaimi et al. [17].Moreover, Turkyilmazoglu [18] analyzed the flow of a micropolar fluid due to a porous stretching sheet and heat transfer.
Furthermore, the work on boundary layer flow due to a shrinking sheet has also attracted much interest.It seems that Miklavčič and Wang [19] investigated the flow over a shrinking sheet and found an exact solution of the Navier-Stokes equations.It was discovered that mass suction is needed to sustain the flow over a shrinking sheet.From physical point of view, vorticity of the shrinking sheet is not restricted within a boundary layer, and the flow is improbable to subsist unless an adequate suction on the boundary is discussed [19].Further, the velocity of the shrinking sheet is confined in the boundary layer with stagnation flow was discovered by Wang [20].This new type of shrinking sheet flow is fundamentally a backward flow as discussed by Goldstein [21].On the other hand, Turkyilmazoglu [22] has analytically studied the steady flow of micropolar fluid and heat transfer over a shrinking sheet under porous wall conditions.He has verified mathematically the bounds of multiple existing solutions of purely exponential kind.It is worth mentioning that the importance and new results on the flow induced by a shrinking sheet in a various fluid and geometrical approach were recently presented by Rosali et al. [23], Naramgari and Sulochana [24], Merkin et al. [25], and Soid et al. [26].In addition, Ros ¸ca and Pop [27] investigated combined effects of micropolar fluid towards a permeable shrinking sheet in the presence of slip effect and Aurangzaib et al. [28] extend it with MHD mixed convection effects.Next, a study of the viscous flow over a shrinking sheet analytically solved using a second-order slip flow model was by Fang et al. [29] where they obtained that the result has two branches, or dual solutions, in a certain range of the parameters, while Bhattacharyya et al. [30] studied that the velocity and thermal boundary layer thicknesses for the second solutions are always larger than the first solutions.
Instead of plane boundary layer, there has been growing interest in the study of the effects of curvature by several investigators recently.Few papers have been discussed on this topic and found that its presence inside the boundary layer is no more negligible as in the case of a flat stretching sheet.Sajid et al. [33] considered stretching curved surface viscous fluid flow and discovered that the boundary layer thickness increases for a curved surface compared to flat surface.Moreover, the drag force to move the liquid on curved surface is less than on a flat surface.Then, they extend the problem to micropolar fluid (see Sajid et al. [34]).Abbas et al. [31] discussed the heat transfer analysis by considering the two heat processes, namely, prescribed surface temperature and described heat flux on MHD boundary layer flow on a curved stretching surface.Further, they continued with slip effect with generation and thermal radiation in Abbas et al. [35].Homogeneous-heterogeneous reactions in MHD flow due to an unsteady curved stretching surface were being observed by Imtiaz et al. [36].
It is worth mentioning that there are also many examples of unsteady flow and in fact, there is no real flow situation, natural or artificial, that does not occupy some unsteadiness.The flows in all manufacturing applications were randomly assumed to be steady.In many engineering instruments, unsteadiness is an integral part of the problem.The helicopter rotor, the cascades of turbomachinery blade, the ship propeller, and so forth usually work in an unsteady environment.Most of the elementary ideas of unsteady viscous flows are described by Telionis [37].Quite different flow behavior is observed for an unsteady shrinking sheet that is for an unsteady stretching sheet.Ros ¸ca and Pop [32] performed a study on the unsteady boundary layer flow over a permeable curved stretching/shrinking sheet.Boundary layer flow of a viscous fluid over an unsteady curved stretching sheet with magnetic field is discussed by Naveed et al. [38].
The aim of the present paper is to study the flow problem due to the unsteady plane, two-dimensional laminar flow caused by a permeable curved stretching/shrinking sheet in a micropolar fluid.The governing partial differential equations are first transformed into a system of ordinary differential equations, before being solved numerically.The expressions for the skin friction coefficients and couple stress coefficients, including velocity and microrotation velocity, are determined to understand the flow characteristics.It is hoped that the results obtained will not only provide useful information for applications, but also serve as a complement to the previous studies.To the best of our knowledge these results are new and original.

Basic Equations
Consider the unsteady two-dimensional boundary layer flow of a micropolar fluid over a permeable curved stretching/shrinking surface coiled in a circle of radius  about the curvilinear coordinates  and , where  is normal to tangent vector at any point on the surface and  is the arc length coordinate along the flow direction, as depicted in Figure 1(a), so that large values of  represent a slightly curved sheet.The geometry for the stretching and shrinking sheet is shown in Figures 1(b) and 1(c), respectively.A constant magnetic field,  0 , is applied in the -direction.By assuming that the magnetic Reynolds number is low, thus, we have neglected the effects of induced magnetic field.Also the applied electric field is assumed to be zero.Moreover, it is assumed that the surface is stretched/shrinked with the velocity   () along the -direction.It should be mentioned that V  (, ) is the mass flux velocity, where V  (, ) < 0 corresponds to suction and V  (, ) > 0 corresponds to injection, respectively, and  is the time.
Under these assumptions along with the boundary layer approximations, the governing equations for the flow are discussed by Abbas et al. [31] and Ros ¸ca and Pop [32], where V and  are the velocity components along and directions, respectively,  is the pressure,  is the density, and ] is the kinematic viscosity.Here,  is assumed to be where  =  * / is the micropolar or material parameter.The initial and boundary conditions applicable to the flow are  < 0: where  is a positive constant,  is the dimensionless constant stretching ( > 0) or shrinking ( < 0) parameter, respectively,  > 0 is for an accelerated surface and  < 0 is for a decelerated surface, respectively, and   () = .Further,  is a constant which varies in the range of 0 ≤  ≤ 1.It is also noted that the case  = 0 is called strong concentration and indicates  = 0 near the wall and represents concentrated particle flows in which the microelements close to the wall surface are unable to rotate.The case  = 0.5 indicates the vanishing of antisymmetric part of the stress tensor and denotes weak concentrations.
The case  = 1 is used for the modeling turbulent flows.
In the present paper, we are interested in the laminar flow configuration; therefore, we consider only strong and weak concentrations of microelements.We will introduce now the similarity variables defined as where prime denotes differentiation with respect to .Thus, we take where  is the constant wall mass transfer, with  > 0 for suction and  < 0 for injection, respectively.Substituting (4) into (1), we get the following system of ordinary differential equations: subject to the boundary conditions where  = √/](1 − ) is the curvature parameter, which depends on the variable t.However the parameter  should be constant, so that we assume  =  0 √ 1 − , where  0 is the characteristic radius of the curved surface.Thus,  =  0 √/] is the constant dimensionless curvature parameter and  = / is the unsteadiness parameter.For the current work, we assume a decelerating shrinking sheet with  ≤ 0.
Eliminating the pressure  from ( 6) and ( 7), we have subject to boundary condition The pressure  can now be determined from (7) as It is significant to point out that the result for flat surface can be obtained with  → ∞.The skin friction coefficient and couple stress coefficient are given as

Result and Discussion
The arrangement of ordinary differential equations (10) and (11) subject to the boundary condition (12) has been solved numerically using shooting method with the aid of shootlib function in Maple software for several values of , , , , , and .Table 1 shows the comparison of the skin friction coefficient for several values of  when  =  =  = 0 and  = 1 (stretching) with those reported by Abbas et al. [31] and Ros ¸ca and Pop [32].It can be observed from Table 1 that the results are in good agreement which verify the accuracy of the present method.In addition, Table 2  [32].In Table 3, we display a sample of our results for skin friction coefficient and couple stress coefficient corresponding to different values of  and  for strong and weak concentrations.Computations are carried out for the following values of the governing parameter both for the weak concentration ( = 0) and for the strong concentration ( = 0.5) cases:  = −0.1,0, 0.1;  = 0, 0.5, 1;  = 1.5, 2, 3; and  = −1, −2, −3.It will be seen later in the graphical result that dual solutions exist but only on certain range.
Variations of skin friction coefficient and couple stress coefficient with  for different values of K are presented in Figures 2 and 3 for  = 0 (strong concentration) while Figures 4 and 5 are for  = 0.5 (weak concentration).It can be seen from these figures that dual solutions exist for each value of considered , while for  equals a certain value, say  =   ( = 1, 2, 3), where   is the critical values of , there is unique solution, and when  <   , there is no solution.Dual solutions are obtained when  >   .The first solution is stable, while the second solution is not, according to the stability analysis performed by Weidman et al. [39]; that is, first solution is physically realizable (stable), while the second solution is not physically realizable (unstable).Also, this investigation confirms that the existence and uniqueness   of the solutions greatly depend on the stretching/shrinking velocity parameter  as also pointed out by Wang [20].In addition, we can see the effects of  from Figures 2-5.From these figures, similar influences have been detected for both weak and strong concentration where values of skin friction and couple stress coefficient decrease when  increase for shrinking case and increase with  for stretching case.
Next, the influence of  and  on the skin friction and couple stress coefficient when  = 0.5 can be examined in Figures 6 and 7, respectively.We noticed that dual solutions  exist when  >   ( = 1, 2, 3) and no solutions exist for  <   ( = 1, 2, 3), where   is the critical value of , a point where first and second solution meet.The values of   increase as  increase.This result was different for normal viscous fluid where the values of   decrease with  as reported by Ros ¸ca and Pop [32].Furthermore, Figures 6 and 7 explained that both skin friction and couple stress coefficient slightly increase when  increase.Figures 8 and 9 depict the variations of the skin friction and couple stress coefficient with the unsteadiness parameter, , for several values of  when  = 0 (strong concentration) and  = −0.1 (shrinking case).Both figures indicate the dual solutions exist for the range of  >   ( = 1, 2, 3), where  =   is the critical point where both solutions merge with one another.It also figured that no solution exists at the range of  <   ( = 1, 2, 3).Moreover, the figures show the values of skin friction and couple stress coefficient increase when  increase.Likewise, similar effect happens as the value of  increases.It is worth pointing out that as  → ∞, it represents the sheet which is not curved and hence, a bigger force is needed to move the fluid on the flat surface than the curved surface.The effects of each governing parameter to the velocity profile,   (), as well as microrotation velocity profile, (), are presented in Figures 10-19.It is worth mentioning that effects of  on velocity profiles are depicted in the same figure since they just have a slight difference compared to their effects on microrotation velocity profiles.Figure 10 shows the influence of  on velocity profile and it shows that the First solution Second solution   velocity decreases when  increase for both weak and strong concentration.An opposite observation has been described by Sajid et al. [34] on steady case, while microrotation velocity decreases with  for strong concentration but increases for weak concentration (see Figures 11(a) and 11(b)).Sajid et al. [34] reported different outcome in both concentrations for steady case.Figures 12 and 13 represent the effects of  on velocity and microrotation velocity profiles, respectively.First solution Second solution  The boundary layer thicknesses for   () and () increase with the increase of  as illustrated in both Figures 12 and  13. Figure 12 shows a slight decrease of velocity profiles for weak concentration compared to the strong concentration.Figure 13(a) demonstrates that by increasing the unsteadiness parameter  will increase () for strong concentration, but () decrease when  increase (Figure 13(b)).Then, we observed the effect of curvature parameter, , on velocity and microrotation velocity profiles, figured in Figures 14 and 15 compared to the second solution as  changed.Increasing  leads to an increasing of velocity profile in both weak and strong concentration.Similar result was observed by Sajid et al. [34] for steady case.Figure 15   is slightly higher than the strong concentration when other parameters are identical.Figure 19(a) illustrates clearly that () increase when  increases in strong concentration.Also, the boundary layer is getting thicker as  decreases.Besides, increasing  would decrease the microrotation velocity profile as shown in Figure 19(b).Generally, we noticed that the boundary layer thicknesses are larger for the second solutions than for the first solutions.All these profiles support the existence of dual nature of the solutions presented in Figures 2-9.They also satisfy the far field boundary conditions (12) asymptotically, which support the numerical results obtained.

Conclusions
An unsteady micropolar fluid flow over a permeable curved stretching and shrinking surface has been investigated.Numerical solutions of the developed mathematical problems are computed and discussed.It is noted that the boundary layer thickness increases for a curved sheet when compared to a flat sheet.The couple stress coefficients increase with  for strong concentration ( = 0) but decrease in weak concentration ( = 0.5).The skin friction decreases when  increase in both concentrations for shrinking case ( < 0) but increase for stretching case ( > 0).It is found that the solutions of the ordinary (similarity) differential equations

Nomenclature 𝑎 :
Constant   : Skin friction coefficient   : Couple stress coefficient : Microinertia density : Curvature parameter : Micropolar or material parameter : Ratio of the microrotation vector component and the fluid skin friction at the wall : Microrotation or angular velocity : Pressure Re  : Local Reynolds number , : Curvilinear coordinates : Radius/distance of the sheet from origin  0 : Characteristic radius : Constantmassflux : Time , V: Velocity component in the and -direction, respectively

Table 1 :
Comparison of results of the skin friction coefficient −  Re  1/2 for several values of  when  =  =  = 0 and  = 1 (stretching).
also provides the comparison values of   (0) for shrinking case by Ros ¸ca and Pop
describes that when  increase, () increases for strong concentration but decreases for weak concentration.Influence of stretching and shrinking parameter, , on   () and () is shown in Figures16 and 17, respectively.Both figures show similar behavior where increasing  enhanced both velocity and microrotation velocity profiles for each weak and strong concentration.Figures18 and 19described the effects of constant mass flux, , on the velocity and microrotation velocity profiles, respectively.Figure18explains that   () increase with  for both weak and strong concentration which happens because suction reduces drag force in order to avoid boundary layer separation.Furthermore, the velocity for weak concentration