The Wall Properties Effect on Peristaltic Transport of Micropolar Non-Newtonian Fluid with Heat and Mass Transfer

The problem of the unsteady peristaltic mechanism with heat and mass transfer of an incompressible micropolar non-Newtonian fluid in a two-dimensional channel. The flow includes the viscoelastic wall properties and micropolar fluid parameters using the equations of the fluid as well as of the deformable boundaries. A perturbation solution is obtained, which satisfies the momentum, angular momentum, energy, and concentration equations for case of free pumping original stationary fluid . Numerical results for the stream function, temperature, and concentration distributions are obtained. Several graphs of physical interest are displayed and discussed.


Introduction
The peristaltic transport is traveling contraction wave along a tube-like structure, and it results physiologically from neuromuscular properties of any tubular smooth muscle. Peristalsis is now well known to physiologists to be one of the major mechanisms for fluid transport in many biological systems. In particular, a peristaltic mechanism may be involved in swallowing food through the esophagus, in urine transport from the kidney to the bladder through the ureter, movement of chyme in the gastrointestinal tract, in the transport of spermatozoa in the ducts efferentes of the male reproductive tracts and in the cervical canal, movement of ovum in the female fallopian tubes, the transport of lymph in the lymphatic vessels, and the vasomotion in small blood vessels such as arterioles, venues, and capillaries. In addition, peristaltic pumping occurs in many practical applications involving biomedical systems. It has now been accepted that most of the physiological fluids behave in general like suspensions of deformable or rigid particles in a Newtonian fluid. Blood, for example, is a suspension of red cells, white cells, and platelets in plasma. Another example is cervical mucus, which is a suspension of macromolecules in a water-like liquid 1 . However only a few studies have considered this aspect since the initial investigation by Raju and Devanathan 2, 3 . The biviscosity fluid can represent the behavior of blood in vessels of small diameter as the mean shear rate is about 20-150 s −1 4 . Several studies have been made to analyze both theoretical and experimental aspects of the peristaltic motion of non-Newtonian fluids in different situations. Böhme and Friedrich 5 studied the mechanism of peristaltic transport of an incompressible viscoelastic fluid by means of an infinite train of sinusoidal waves traveling along the wall of the duct in the case of plane flow. The effects of an Oldroyd-B fluid on the peristaltic mechanism are examined by Hayat et al. 6 , under the long wavelength assumption. They noticed that in the narrow part of the channel, the behavior of an Oldroyd-B fluid is much more different from that of a Newtonian fluid than in the wide part of the channel. El-Shehawey et al. 7 studied the peristaltic motion of an incompressible non-Newtonian fluid through a porous medium. They showed that the pressure rise increases as the permeability decreases and noted that both pressure rise and friction force do not depend on permeability parameter at a certain value of flow rate. Haroun 8 investigated the peristaltic flow of a third-order fluid in an asymmetric channel under the assumption of long wavelength approximation. He expanded the velocity components and pressure in a regular perturbation series in a small parameter Deborah number that contained the non-Newtonian coefficients appropriate to shear-thinning. The effects of slip boundary conditions on the dynamics of fluids in porous media is investigated by El-Shehawey et al. 9 . They studied the flow of non-Newtonian Maxwell fluids in an axisymmetric cylindrical tube pore , in which the flow is induced by traveling transversal waves on the tube wall peristaltic transport . Hakeem et al. 10 discussed the effect of an endoscope on the peristaltic mechanism of a generalized Newtonian fluid. Mekheimer 11 studied the effect of a uniform magnetic field on peristaltic transport of a blood in nonuniform two-dimensional channels, when blood is represented by a couplestress fluid. Tsiklauri and Beresnev 12 analyzed the effect of viscoelasticity on the dynamics of fluids in porous media by studying the peristaltic flow of a Maxwell fluid in circular tube, in which the flow is induced by a wave traveling on the tube wall.
The micropolar fluid represents fluids, which consist of rigid, randomly oriented or spherical particles suspended in a viscous medium where the deformation of the particles is ignored. The theory of thermo-micropolar fluid was developed by Eringen 13 . Agrawal and Dhanapal 14 considered the micropolar fluid to be the model of blood flow in small arteries, and the calculation of theoretical velocity profiles is observed, in good agreement with experimental data 15 . They have many applications in physiological and chemical engineering 16 . Srinivasachrya et al. 17 studied the problem of peristaltic transport of a micropolar fluid in a circular tube. Girija Devi and Devanathan 18 studied the peristaltic flow of a micropolar fluid in a cylindrical tube with a sinusoidal deformation of small amplitude travelling down its flexible wall for the case of low Reynolds number flow devoid of wall properties like tension and damping. However, the wall properties are essential to be taken into consideration in various real situations. The peristaltic motion of a simple microfluid which accounts for microrotation and microstretching of the particles contained in a small volume element is studied by Philip and Chandra 19 , using long wavelength approximation.
Combined heat and mass transfer problems are of importance in many processes and have, therefore, received a considerable amount of attention in recent years. In processes such as drying, evaporation at the surface of a water body, energy transfer in a wet cooling tower and the flow in a desert cooler, heat and mass transfer occurs simultaneously. Possible applications of this type of flow can be found in many industries. For example, in the power industry, among the methods of generating electric power is one in which electrical energy is extracted directly from a moving conducting fluid. Many practical diffusive operations involve the molecular diffusion of a species 20 . Radhakrishnamacharya and Radhakrishna 21 investigated the peristaltic flow with heat transfer in a nonuniform channel. The effect of elasticity of the flexible walls on peristaltic transport of an incompressible viscous fluid, with heat transfer, in a two-dimensional uniform channel under long wavelength approximation is explained by Radhakrishnamacharya and Srinivasulu 22 . El Dabe et al. 23 studied heat and mass transfer of a steady slow motion of a Rivlin-Ericksen fluid in tube of varying cross-section with suction. Agrawal et al. 24 obtained numerical solutions of flow and heat transfer of a micropolar fluid at a stagnation point on a porous stationary wall. They observed that the heat sources increase the velocity and temperature in the pipe while the heat sinks decrease them.
Muthu et al. 25 carried out a study of the peristaltic motion of an incompressible micropolar fluid in two-dimensional channel. They investigated the effects of viscoelastic wall properties and micropolar fluid parameters on the flow using the equations of the fluid as well as of the deformable boundaries. The non-Newtonian property and equations of heat and mass transfer were not taken into their consideration. Because of the wide range of practical importance of the heat and mass transfer, the present study considered the heat and mass transfer of an unsteady peristaltic motion of a micropolar non-Newtonian biviscosity fluid 4 . The following analysis includes the dynamic boundary condition. Analytical approximate solutions for the stream function, microrotation velocity, temperature, and concentration equations are obtained as a power series in terms of the small amplitude ratio. We have shown the relation between the different parameters of motion in order to investigate how to control the motion of the fluid by changing these parameters.

Formulation of the Problem
Consider a two-dimensional symmetric unsteady flow of an incompressible micropolar biviscosity fluid in an infinite channel of uniform thickness 2d, with heat and mass transfer. The walls of the channel are flexible membranes on which they are imposed traveling sinusoidal waves of moderate amplitude see Figure 1 .
We choose a rectangular coordinate system such that the axes x and y are in the directions of wave propagation and normal to the mean position of the membranes, respectively. The origin is located at the center line of the channel. It was shown by Wilson and Taylor 26 that the limiting process by which the biviscosity model is approached with the limiting process implicit in lubrication theory. Models of this type are much easier to handle mathematically than many models such as Oldroyd's and represent the experimental facts just as well, at least in many cases.
The biviscosity model which has been used in similar contexts by Nakayama and Sawada 4 can be written as 2π c e ij , π < π c .

4 Mathematical Problems in Engineering
The following quantity is introduced as a nondimensional parameter including π c where μ B is the plastic viscosity, p y is the yielding stress, π e ij e ij , where e ij is the i, j component of the deformation rate and the value of γ denotes the upper limit of apparent viscosity coefficient. The biviscosity model is approached in the limit γ → ∞ see Nakayama and Sawada 4 for a few further details .
The governing continuity, momentum, angular momentum, temperature, and concentration equations for this problem can be written as 13 where u x, y, t and υ x, y, t are the velocity components in the x and y directions, respectively, N θ x,y,t is the microrotation velocity component in the direction normal to both the x and y axes. Here ρ is the density of the fluid, P is the pressure, J is the microinertia constant, k is the vortex viscosity coefficient also known as the coefficient of gyroviscosity , and γ is the spin-gradient viscosity. Further, the material constants μ B and k satisfy the following inequalities 13 : Mathematical Problems in Engineering The temperature and concentration of fluid are T x, y, t and C x, y, t , C p is the specific heat at constant pressure, D is the coefficient of mass diffusivity, K T is the thermal diffusion ratio, T m is the mean fluid temperature, k c is the thermal conductivity and α c is the heat conduction parameter for a micropolar fluid. The last term in 2.8 signifies the thermaldiffusion effect.
We consider a symmetric motion of the flexible walls in the boundary conditions. Let the vertical displacements of the upper and lower walls be η and −η, respectively, where where a is the amplitude, λ is the wavelength, and c is the wave speed. We assume that the walls are inextensible so that only lateral motion takes place and the horizontal displacement of the wall is zero 25 .
Thus the no-slip boundary conditions for the velocity and microrotation are u 0, N θ 0 at y ± d η x, t .

2.11
The dynamic boundary conditions are imposed on the fluid by the symmetric motion of the flexible walls, which can be written as 18 Here L η is the pressure at the walls, is the tension in the membrane, m is the mass per unit area, and n is the coefficient of viscous damping force.

Mathematical Problems in Engineering
The plate temperature starts oscillating about a nonzero mean temperature. Under these physical conditions, the temperature and concentration at the upper wall give 2.14 where C m is the mean fluid concentration, T w and C w are the uniform temperature and concentration at the lower wall, that is, 3 allows the use of the stream function ψ x, y, t such that By introducing the following nondimensional quantities: Ec

2.17
Mathematical Problems in Engineering 7 substituting from 2.17 into 2.3 -2.8 , eliminating P between 2.4 and 2.5 , and dropping the accent mark for simplicity, the governing equations may be written as Also, the boundary conditions at y ± 1 η x, t , when η x, t ε cos α x − t , are

2.24
The parameter ε represents the amplitude ratio, α is wave number, Re is the Reynolds number, μ 1 is the ratio between the viscosity coefficient for micropolar fluids and the classical viscosity coefficient, and M is the micropolar fluid parameter characterizing spin-gradient viscosity. The parameter M can be thought of as a fluid property depending upon the size of the microstructure 18 . We note that R is the modified Reynolds number and involves the quantity J. In this paper, we considered that the effect of microinertia is neglected and R is taken to be zero 24 . The nondimensional quantities Pr are the Prandtl number, Ec is Eckert number, af is the dimensionless heat conduction coefficient of the micropolar fluid, Sc is the Schmidt number, Sr is the Soret number and represents the thermal-diffusion effect, K 2 and K 3 represent the dissipative and rigiditive feature of walls, and m 1 indicates the stiffness property of walls.

Method of Solution
Assume the amplitude ratio ε of the wave to be small. As given by Muthu et al. 25 , we seek for an approximate solution as a power series in terms of ε in the form

3.4
Here the asterisk denotes complex conjugate. The solutions for ψ 0 y and N 0 y are obtained as follows: where K Re / 2 1 γ −1 μ 1 dP/dx 0 is the poiseuille flow parameter for the micropolar fluid. For pure peristalsis, which means that the flow is generated by wall motion only. The Mathematical Problems in Engineering 9 pressure gradient dP/dx 0 0; this implies that K 0, which gives ψ 0 y 0 and N 0 y 0. The expressions of φ 1 and ξ 1 for the case of free pumping can be written as

3.8
To obtain the mass and heat transfer, we shall use the perturbation scheme 3.3 , 3.4 in 2.20 and 2.21 , and equating the zero-order terms on both sides, we obtain the following set of equations: In case of plane Poiseuille flow of micropolar fluid ε 0 , the governing 3.9 and 3.10 may be solved in the light of the boundary conditions 3.11 with the constant pressure gradient ∂P/∂x 0 . At this stage, the temperature and concentration distributions T 0 and C 0 , respectively, are given by Let us consider the case in which the pressure gradient ∂P/∂x 0 vanishes. In this case there will be no flow if the wall motion stops. Hence ψ 0 N 0 0 and T 0 C 0 1, and equating the first-order terms on both sides of 2.20 and 2.21 , we get the following set of equations: The boundary conditions are Mathematical Problems in Engineering

3.17
In this case, the governing equations for T 20 and C 20 reduce to with the boundary conditions The solutions of 3.18 are 3.20 where g 1 − g 3 , f 1 − f 10 K 3 and A 5 − A 8 are given in the appendix.

3.25
Hence g 4 − g 17 , f 11 − f 28 and A 9 − A 12 are given in the appendix. Using the expressions for Ψ 1 , T 1 , C 1 , T 20 , C 20 , T 22 , and C 22 in 3.1 , 3.3 , and 3.4 , we can determine the stream function Ψ x, y, t , temperature T x, y, t , and concentration C x, y, t of micropolar fluid for the case of free pumping as 3.27 3.28 where φ * 1 , T * 1 y , C * 1 y , T * 22 y , and C * 22 y are conjugate to 3.6 , 3.16 , 3.17 , 3.24 , and 3.25 , respectively.

Results and Discussion
To discuss the effect of the parameters of micropolar, non-Newtonian fluid, and viscoelastic membrane properties on the solutions of the considered problem. A numerical results are calculated from formula 3.6 , 3.27 , and 3.28 for the temperature T x, y, t , concentration C x, y, t and stream function ψ x, y, t and shown by Figures 2-24. Take into account the parameters t and x have the numerical value 2.0, while parameter α and amplitude ratio ε take the value 0.5 and 0.01, respectively, but α takes the value 2 for the stream function. Figures 2-4 illustrate the effects of μ 1 k/μ B , M 2d μ B /γ 0.5 , and af, in order, on the temperature distribution with taking into account that the elastic wall K 2 0 . The ratios k/μ B and γ/μ B measure the relative strengths of the vortex viscosity coefficient to the viscosity coefficient and the spin gradient viscosity coefficient to the viscosity coefficient, respectively. These coefficients μ B , k, and γ may be greater or equal to zero. For example, the blood 50% hematocrit has the values μ B 0.0029, k 0.000232, and γ 0.000001 27 . When the viscous effects are much larger than the spin gradient viscosity effects, γ/μ B may tend to be zero, and the parameter M will tend to infinity. Ahmadi 28 has stated that the parameter μ 1 depends on the shape and concentration of the microelements while the parameter M signifies the size of microstructure. In this context, small large value of μ 1 means low high concentration of the microelements. Similarly, a given small large value of M is related to small large size of the particles. It is noted from Figure 2 that the temperature distribution decreases with the increasing of μ 1 near the lower wall, while it increases as μ 1 increases when y ≥ −0.7. In Figure 3, it can be seen that the temperature distribution increases as M increases. Also, Figure 4 shows that the temperature increases as af increases but at y > 0.8, an opposite effect is noticed near the upper wall. The results in our study are consistent with those which are obtained by Agarwal et al. 24 in case of isothermal wall. The choice K 2 0 implies that the walls move up and down with no damping force on them and hence indicates that the membrane is treated as an elastic wall. The parameter K 2 depends upon the wall tension and represents the rigid nature of the walls. Figures 5, 6, and 7 illustrate the effect of the parameters of the viscoelastic wall properties K 2 , K 3 , and m 1 on the temperature distribution, respectively. It is found that the temperature distribution increases as K 2 increases as shown in Figure 5, but near the lower wall, it decreases with the increase of K 2 . The temperature distribution decreases with the increase of both K 3 and m 1 , but when y ≥ −0.7, it starts increasing as both K 3 and m 1 increase for the case of elastic walls K 2 0 , as shown in Figures 6 and 7. This phenomenon reflects the fact that the increasing of the values of tension or dissipative for the membrane leads to loss of heat from the walls and gives it to the fluid. Therefore, the temperature distribution of fluid will increase. In Figure 8, we see the effect of Reynolds number Re on the temperature distribution. It is observed that the temperature distribution behaves a dual role with the variation of Re. The variations of temperature distribution for various values of Prandtl number Pr and  Eckert number Ec are displayed in Figures 9 and 10, respectively. The graphical results of Figures 9 and 10 indicate that the temperature of fluid increases as Pr and Ec increase. While for |y| > 0.7, near the lower and upper wall, the temperature decreases with the increase of Pr and Ec. In the case when the temperature of the plate is periodic, Das et al. 29 showed that the temperature decreases with the increase of the Prandtle number of Newtonian fluid. In Figure 11, the temperature is plotted against y for various values of the upper limit apparent viscosity coefficient γ; it is observed that an increase of γ increases the temperature. Figures 12-18 reveal the influence of physical parameters entering in the problem on the concentration distribution of fluid. The effect of the coupling parameter μ 1 on concentration is indicated in Figure 12; we find that the concentration distribution of fluid decreases with the increase of μ 1 . It has been observed that the concentration distribution decreases with the increase of af, but an opposite effect is seen at y > 0.8. This is noticed in Figure 13. Also, in Figures 12 and 13, for small values of μ 1 and af, its effects are rather insignificant. In Figures 14, 15, and 16, it is found that the increase of the dissipative, rigid nature of the walls and stiffness in the walls causes a decrease of the concentration distribution. The behavior of the concentration distribution with the nondimensional distance y for various values of Soret number Sr is depicted in Figure 17. It is clear, from Figure 17, that the effect of Sr on the concentration can be neglected near the lower wall. When −0.8 < y < 0.5, the concentration decreases with the increase of Sr, but at y > 0.5, an increase in Sr leads to  Figure 17, at the upper wall, is in agreement with those obtained by El Dabe et al. 23 . Figure 9 illustrates the effects of the upper limit apparent viscosity coefficient γ on the concentration distribution. It is found that the concentration decreases near lower and upper walls, but it increases with the increase of γ at the middle of the channel. The effects of various physical parameters on the stream function ψ are indicated in  In these figures the stream function is plotted versus the dimensionless distance y, and it is clear that there is an inversely non-linear relation between y and ψ and that the obtained curves of all figures will intersect at the origin where y 0, ψ 0. Figure 19 shows the variation of ψ with y for different values of the coupling parameter μ 1 ; it is observed that, in the region −1 ≤ y < 0, the stream function decreases by increasing μ 1 , while in the region 0 < y ≤ 1, ψ increases by increasing μ 1 . Note that, when μ 1 ≥ 1, the peristaltic pumping region becomes slightly wider as μ 1 increases, and if μ 1 < 1, it is slightly narrow as μ 1 decreases. Figure 20 illustrates the stream function ψ for different values of M. From this figure, it is observed that the effect of M on ψ is similar to the effect of μ 1 on ψ illustrated in Figure 19. In Figures 21 and 22, the effects of Reynolds number Re and the parameter of dissipates of walls K 2 on the stream function ψ, respectively, are reverse to the effect of μ 1 and M on ψ, that is, in the region −1 ≤ y < 0, ψ increases by increasing both of Re and K 2 , while in the region 0 < y ≤ 1, ψ decreases by increasing both of Re and K 2 . Note that the peristaltic pumping region becomes slightly wider as Re and K 2 increase. Figures 23  and 24 show that the behaviors of the stream function ψ with the dimensionless distance y for different values of the parameter of rigidities feature of walls K 3 and upper limit apparent viscosity coefficient γ, respectively, are similar to the behaviors of ψ for different values of μ 1 and M illustrated in Figures 19 and 20, with the only difference that the peristaltic pumping will disappear in these cases.