(ProQuest: ... denotes non-US-ASCII text omitted.)
I. Ahmad 1 and N. Ali 2 and A. Abbasi 1 and W. Aziz 3 and M. Hussain 1 and M. Ahmad 1 and M. Taj 1 and Q. Zaman 1
Academic Editor:Takeshi Iwamoto
1, Department of Mathematics, University of Azad Jammu & Kashmir, Muzaffarabad 13100, Pakistan
2, Faculty of Basic and Applied Sciences, Department of Mathematics and Statistics, International Islamic University, Islamabad 44000, Pakistan
3, Department of Computer Sciences and Information Technology, University of Azad Jammu & Kashmir, Muzaffarabad 13100, Pakistan
Received 16 May 2014; Revised 11 August 2014; Accepted 11 August 2014; 28 October 2014
This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
Peristaltic pumping is a mean of fluid transport caused by waves of contraction along the boundary of a channel/tube. To analyze the properties of peristalsis theoretically, the flow caused by sinusoidal traveling waves of contraction along the boundary has been considered in an infinite channel/tube. In general, flow in a peristaltic channel is driven by two independent mechanisms: the movement of the walls and due to the mean pressure gradient along the channel. After the pioneering work of Shapiro et al. [1] the fluid flow due to prescribed wall motion has been described theoretically by Ebaid [2]. The approach of Fung and Yih [3] was used to analyze the peristaltic flow by Hayat et al. [4] and Haroun [5]. In the above mentioned studies, peristaltic motion was investigated under various assumptions. One of the important aspects in such studies is the nature of the fluid. Since, it is a fact that most of the biological and industrial fluids are non-Newtonian in nature, therefore, peristaltic flow of these fluids has also been studied by Hayat et al. [6] considering non-Newtonian Maxwell fluid in a channel. Srinivas et al. [7] studied the MHD peristaltic motion with slip and heat characteristics. The effects of wall flexibility on the peristaltic motion have also attracted the attention of various researchers in the field. Mittra and Prasad [8] initiated the study of peristaltic transport of viscous fluid with compliant walls. Their work has been generalized by Muthu et al. [9] for viscous incompressible fluid in circular tube and Abd Elnaby and Haroun [10] for two-dimensional channel. The work in [10] has been extended for Maxwell fluid by Ali et al. [11]. Peristaltic motion of Jeffrey fluid in a channel having compliant walls was also studied by Hayat et al. [12]. Recently the heat and mass transfer effects on the peristaltic flow of an Oldroyd-B fluid in a channel with complaint walls were investigated by Hayat et al. [13]. Motivated by the above mentioned investigations, the aim of the present paper is to extend the work of [13] by considering the fluid which obeys the constitutive equation of Burger's fluid model. The Burger's fluid is a rate type viscoelastic fluid model [14] which has been found quite successful in predicting the behavior of viscoelastic material, food products such as cheese and soil. The viscoelastic behavior of fine-grained polycrystalline olivine can also be predicted by Burger's fluid model [15]. Moreover, this model has also been used for calculation of transient creep proportion of the earth's mantle [16] and in the interpretation of postglacial uplift [17].
The paper is arranged in the following manner. Flow geometry is explained in Section 2 while governing equations are developed in Section 3. Series solution is provided in Section 4. Results and discussion presented in Section 5. Some conclusions are made at the end of article in Section 6.
2. Flow Geometry
Let us consider the flow of magneto hydrodynamic (MHD) Burger's fluid due to sinusoidal motion of the compliant walls in a channel. A constant magnetic field of strength B 0 is applied in the direction normal to the flow. The deformation of the walls of the channel is represented by the expression [figure omitted; refer to PDF] where a is the amplitude of the peristaltic wave, λ is the wavelength, t is the time, c is the speed of the wave, and d 1 is the half-width of the channel.
3. Governing Equations
Let u and v denote the velocity components in the axial and transverse direction, respectively, then we can write [figure omitted; refer to PDF] The law of conservation of mass and momentum for the flow of an incompressible fluid are [figure omitted; refer to PDF] where V is the velocity, ρ is the density, p is the pressure, and S is the extra stress tensor, which satisfies the following relation for a Burger's fluid: [figure omitted; refer to PDF] In the above equation μ is the dynamic viscosity, A 1 is the first Rivlin-Ericksen tensor, λ 1 and λ 2 are the relaxation and retardation times, respectively, β is Burger's fluid parameter, L is the gradient of velocity, and D / D t is the upper convective derivative and defined as [figure omitted; refer to PDF] Using the definition of velocity and D / D t in (3) we get the following equations: [figure omitted; refer to PDF] where the components of extra stress appearing in the above equations satisfy three equations that can be obtained directly from (4a).
The relevant boundary conditions are [8] [figure omitted; refer to PDF] where τ is the tension in the wall, m is the mass per unit area, d is the coefficient of viscous damping, and σ is the electrical conductivity of fluid. Defining the stream function ψ such that u = ∂ ψ / ∂ y , v = - ∂ ψ / ∂ x , and the following dimensionless variables: [figure omitted; refer to PDF] we find that the continuity equation (5) is identically satisfied and (8)-(9), after dropping the primes, reduce to [figure omitted; refer to PDF] In the above equations η = ± ( 1 + [straight epsilon] sin 2 π ( x - t ) ) , [straight epsilon] = a / d 1 is the dimensionless amplitude, δ = d 1 / λ denotes the characteristic ratio of transversal and axial scales of peristaltic motion, R = c d 1 / ν is the Reynolds number, M = σ / μ B 0 d 1 denotes the Hartman number, E 1 = - τ d 1 3 / λ 3 μ c , E 2 = m c d 1 3 / λ 3 μ , E 3 = d d 1 3 / λ 3 μ are the nondimensional elasticity parameters, and β * = β c 2 / d 1 2 is the dimensionless Burger's fluid parameter.
4. Solution of the Problem
Due to nonlinearity of (11) it seems impossible to get a closed form solution. Therefore, some plausible assumptions must be made to put these equations in tractable form. We would like to mention here that without such assumptions one cannot obtain solution of these equations even for Newtonian fluid. Following [1], we assume that the δ is small and proceed to write [figure omitted; refer to PDF] Substituting expressions (13) into (11) and (12) and collecting the terms of equal powers of δ , we get the following systems.
4.1. Zero-Order System
[figure omitted; refer to PDF]
4.2. First-Order System
[figure omitted; refer to PDF]
4.3. Second-Order System
[figure omitted; refer to PDF] The solution of these systems yields [figure omitted; refer to PDF] where [figure omitted; refer to PDF] and value of u 2 is given in the appendix.
The value of constants involved in (18) is given in the appendix.
5. Results and Discussion
The purpose of this section is to provide the graphical illustration of the solution obtained in previous section. Emphasis has been given to examine the influence of emerging parameters such as δ , R , M , and the rheological parameter of the Burger's fluid β on the axial velocity and the stream function.
In fact the novelty of the present analysis is based on the effects of Burger's fluid parameter β on various flow characteristics. This is because of the fact that the results in the absence of β are available already in the literature [13] for comparison. Particularly the velocity u is plotted for different values of the involved parameter including the Burger's fluid parameter β . Figure 1 interprets the effects of Hartman number M and β on u . It is noted from this figure that u is an increasing function of β , while it decreases with M . The decrease of u with M is a result of the fact that magnetic force acts in the transverse direction and hence resists the flow due to peristalsis.
Variation on u when E 1 = 1 , E 2 = 0.2 , E 3 = 0.5 , [straight epsilon] = 0.2 , δ = 0.01 , t = 0.2 , R = 0.01 , λ 1 = 0.1 , λ 2 = 0.05 , and x = 0.3 plotted for (a) for different values of M and (b) different Burger's fluid parameter β .
(a) [figure omitted; refer to PDF]
(b) [figure omitted; refer to PDF]
The effects of R and δ on u can be seen through Figure 2. This figure shows that u increases by increasing R and δ . The Reynolds number represents the relative magnitude of inertial to the viscous force and with the increase of Reynolds number, the viscous forces decrease resulting in increases in velocity. Similarly with an increase in δ the wavelength of wave decreases which causes the increase in the velocity. It is worth mentioning that these results are in accordance with those obtained in [11] for a Maxwell fluid. A maximum in the velocity at center of the channel and an increase of velocity with R and δ are also observed from Figure 2.
Variation on u when E 1 = 1 , E 2 = 0.2 , E 3 = 0.5 , [straight epsilon] = 0.2 , β = 0.05 , t = 0.2 , M = 2.5 , λ 1 = 0.1 , λ 2 = 0.0 , and x = 0.3 plotted for (a) for different values of R and (b) different values of δ .
(a) [figure omitted; refer to PDF]
(b) [figure omitted; refer to PDF]
Figure 3 is prepared to see the effects of complaint wall parameters E 1 , E 2 , and E 3 on longitudinal distribution of velocity u . This figure reveals that the velocity increases when both wall parameters E 1 and E 2 increase because of the increase in flexibility and elasticity of the wall. However, velocity decreases with E 3 due to the increase in rigidity of the wall. Note that we recover the results reported in [11] for λ 2 = β = 0 .
Variation on u when R = 0.01 , δ = 0.01 , [straight epsilon] = 0.2 , β = 0.01 , t = 0.2 , M = 2.6 , λ 1 = 0.1 , λ 2 = 0.0 , and x = 0.3 plotted for (a) for different values of E 1 , (b) for different values of E 2 , and (c) for different values of E 3 .
(a) [figure omitted; refer to PDF]
(b) [figure omitted; refer to PDF]
(c) [figure omitted; refer to PDF]
Figure 4 illustrates the effects of λ 1 and λ 2 on u . It is observed that when we increase the relaxation parameter the velocity u gets increased due to the increment in λ 1 . The effects of retardation parameter λ 2 on the velocity are quite opposite to that of λ 1 .
Variation on u when E 1 = 1 , E 2 = 0.2 , E 3 = 0.5 , [straight epsilon] = 0.2 , β = 0.05 , t = 0.2 , M = 1 , δ = 0.01 , R = 0.01 , and x = 0.3 plotted for (a) for different values of λ 1 and (b) for different values of λ 2 .
(a) [figure omitted; refer to PDF]
(b) [figure omitted; refer to PDF]
The purpose of Figures 5, 6, 7, and 8 is to capture the variation in size and circulation of the trapped bolus, in particular, the effects of E 1 , E 2 , β , and E 3 . It is observed from these figures that the size of bolus decreases by increasing E 1 and E 2 while it increases with E 3 and β .
Streamlines when E 1 = 0.2 , E 2 = 0.1 , E 3 = 0.1 , [straight epsilon] = 0.1 , δ = 0.01 , t = 0.0 , R = 0.01 , M = 1 , λ 1 = 0.1 , and λ 2 = 0.05 . (a) β = 0.01 and (b) β = 0.05 .
(a) [figure omitted; refer to PDF]
(b) [figure omitted; refer to PDF]
Streamlines when β = 0.05 , E 2 = 0.1 , E 3 = 0.1 , [straight epsilon] = 0.1 , δ = 0.01 , t = 0.0 , R = 0.01 , M = 1 , λ 1 = 0.1 , and λ 2 = 0.05 . (a) E 1 = 0.1 and (b) E 1 = 0.3 .
(a) [figure omitted; refer to PDF]
(b) [figure omitted; refer to PDF]
Streamlines when β = 0.05 , E 1 = 0.1 , E 3 = 0.1 , [straight epsilon] = 0.1 , δ = 0.01 , t = 0.0 , R = 0.01 , M = 1 , λ 1 = 0.1 , and λ 2 = 0.05 . (a) E 2 = 0.1 and (b) E 2 = 0.3 .
(a) [figure omitted; refer to PDF]
(b) [figure omitted; refer to PDF]
Streamlines when β = 0.05 , E 1 = 0.1 , E 2 = 0.1 , [straight epsilon] = 0.1 , δ = 0.01 , t = 0.0 , R = 0.01 , M = 1 , λ 1 = 0.1 , and λ 2 = 0.05 . (a) E 3 = 0.1 and (b) E 3 = 0.5 .
(a) [figure omitted; refer to PDF]
(b) [figure omitted; refer to PDF]
6. Conclusions
Peristaltic motion of a (MHD) non-Newtonian fluid with Burger's constitutive equation in channel with complaint walls is taken into account. The results for velocity and stream function are obtained and compared with the literature. It is observed that Reynolds number and δ cause an increase in the velocity while a decrease in velocity is noted by increasing Hartman number. In previous mathematical model [13], Burger's parameter is absent and therefore such models cannot accurately describe the peristaltic motion of food stuff like cheese and so forth. The physical impact of this parameter is quite interesting. Our mathematical model predicts an increase in velocity by increasing Burger's fluid parameter. The model further reveals that the size and circulation of the trapped bolus increase with Burger's fluid parameter.
Acknowledgments
The authors are grateful to the reviewers for their valuable suggestions which improved the quality of the paper.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
[1] A. H. Shapiro, M. Y. Jaffrin, S. L. Weinberg, "Peristaltic pumping with long wave length at low Reynolds number," Journal of Fluid Mechanics , vol. 37, no. 4, pp. 799-825, 1969.
[2] A. Ebaid, "Effects of magnetic field and wall slip conditions on the peristaltic transport of a Newtonian fluid in an asymmetric channel," Physics Letters A , vol. 372, no. 24, pp. 4493-4499, 2008.
[3] Y. C. Fung, C. S. Yih, "Peristaltic transport," Journal of Applied Mechanics , vol. 33, pp. 669-675, 1968.
[4] T. Hayat, N. Ali, S. Asghar, "An analysis of peristaltic transport for flow of a Jeffrey fluid," Acta Mechanica , vol. 193, no. 1-2, pp. 101-112, 2007.
[5] M. H. Haroun, "Effect of relaxation and retardation time on peristaltic transport of the Oldroydian viscoelastic fluid," Journal of Applied Mechanics and Technical Physics , vol. 46, no. 6, pp. 842-850, 2005.
[6] T. Hayat, N. Ali, "Peristaltic mechanism of a Maxwell fluid in an asymmetric channel," Nonlinear Analysis: Real World Applications , vol. 9, no. 4, pp. 1474-1490, 2008.
[7] S. Srinivas, R. Gayathri, M. Kothandapani, "The influence of slip conditions, wall properties and heat transfer on MHD peristaltic transport," Computer Physics Communications , vol. 180, no. 11, pp. 2115-2122, 2009.
[8] T. K. Mittra, S. N. Prasad, "On the influence of wall properties and Poiseuille flow in peristalsis," Journal of Biomechanics , vol. 6, no. 6, pp. 681-693, 1973.
[9] M. A. Muthu, K. Rathish, C. Peeyush, "Peristaltic motion in circular cylindrical tubes: effect of wall properties," Indian Journal of Pure and Applied Mathematics , vol. 32, no. 9, pp. 1317-1325, 2001.
[10] M. A. Abd Elnaby, M. H. Haroun, "A new model for study the effect of wall properties on peristaltic transport of a viscous fluid," Communications in Nonlinear Science and Numerical Simulation , vol. 13, no. 4, pp. 752-762, 2008.
[11] N. Ali, T. Hayat, S. Asghar, "Peristaltic flow of a Maxwell fluid in a channel with compliant walls," Chaos, Solitons and Fractals , vol. 39, no. 1, pp. 407-416, 2009.
[12] T. Hayat, M. Javed, N. Ali, "MHD peristaltic transport of a Jeffery fluid in a channel with compliant walls and porous space," Transport in Porous Media , vol. 74, no. 3, pp. 259-274, 2008.
[13] T. Hayat, S. Hina, A. A. Hendi, "Heat and mass transfer effects on peristaltic flow of an Oldroyd-B fluid in a channel with compliant walls," Heat Transfer--Asian Research , vol. 41, no. 1, pp. 63-83, 2012.
[14] J. M. Burgers, "Mechanical considerations-model systems-phenomenological theories of relaxation and of viscosity," First Report on Viscosity and Plasticity , Nordemann Publishing Company, New York, NY, USA, 1935.
[15] A. H. Gerritsen, C. A. P. M. van Gurp, J. P. J. van der Heide, A. A. A. Molenaar, A. C. Pronk, "Prediction and prevention of surface cracking in asphaltic pavements," in Proceeding of the 6th International Conference on Structural Design of Asphalt Pavements (ISAP '87), vol. 1, University of Michigan, Ann Arbor, Mich, USA, July 1987.
[16] P. N. Chopra, "High-temperature transient creep in olivine rocks," Tectonophysics , vol. 279, no. 1-4, pp. 93-111, 1997.
[17] B. H. Tan, I. Jackson, J. D. F. Gerald, "High-temperature viscoelasticity of fine-grained polycrystalline olivine," Physics and Chemistry of Minerals , vol. 28, no. 9, pp. 641-664, 2001.
Appendix
Here we provide the value of u 2 along with various constants appearing in it and in the expressions of u 0 and u 1 , as follows: [figure omitted; refer to PDF]
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Copyright © 2014 I. Ahmad et al. I. Ahmad et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Abstract
A theoretical analysis is presented for the peristaltic motion of a magneto-hydrodynamic (MHD) non-Newtonian fluid in channel with complaint walls. The fluid obeys viscoelastic non-Newtonian model with Burger's constitutive equation. The relevant equations are first developed and then solved using perturbation technique. Expressions of stream function and velocity components are constructed under the assumption that δ (characteristic ratio of transversal and axial scales of peristaltic motion) is small. The results indicate the strong effects of Burger's fluid parameter, Hartman number, Reynolds number, and complaint wall parameters on the velocity field and stream function. The obtained solutions are shown graphically for the different values of involved parameters.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer