Imposed pressure driven flow in peristalsis

Abstract

Peristaltic transport phenomena are of great significance in biological sciences. The physiological transport of fluid takes place under the action of peristalsis generated as a pressure gradient. The peristaltic waves generate a pressure gradient which is responsible for the fluid flow in the forward direction. The further properties of this phenomena can be seen if an imposed pressure gradient is applied in addition to the one appearing due to peristaltic waves. This situation has not been discussed in the literature that needs further attention. The effects of the wavy boundaries and imposed pressure on the velocity of the flow field are analyzed. Here we impose a question: what happens if an imposed pressure gradient is also applied? This question of physical importance has not been addressed; and thus, remains the topic of this study. In previous papers of peristaltic motion, the flow generated by peristaltic waves only has been examined while in this study we will discuss the contribution of imposed pressure gradient on velocity field. The analytical results for the velocity field are obtained using the boundary perturbation method. The study shows that the impact of the wavy boundaries on the flow increases with the increase in corrugation parameter and imposed pressure.

Introduction

This article is concerned with the fluid mechanics of peristaltic pumping under conditions for which the length of peristaltic wave is large compared with the channel width and the appropriate Reynolds number is sufficiently small for the flow to be considered inertia- free. The basic principle of peristaltic motion is based on the sinusoidal wave boundary which is produced due to the traveling of sinusoidal wave along flexible walls of the channel. There are other muscular structures which are responsible for pumping the fluids by means of peristaltic motion. Technologically, peristaltic devices are commonly used for slurries, pumping of blood through human body, transporting corrosive fluids etc. In bio-mechanics; the peristaltic phenomena have great relevance in arterioles, heart-lung machine, intestines, ureter and many other. The procedures are adopted by pumping of blood, food in small intestines, passing of fluid from kidney to ureter. Peristalsis in human physiological systems have been the reasons of many recent investigations

Mathematically, the peristaltic flows induced by the waves traveling along a flexible channel have been treated under various geometrical situations and various fluid models. It must be kept in mind that the fluid is transported by the pressure gradient developed in the peristalsis mechanism. Before going further, we would like to refer few papers of major interest.

A literature review gives a general idea of the literature and theories related to the subject of the present problem. With the help of literature review, historical aspects of theories and literature is reached which lays importance forthrightly to the subject of the research. Even though peristalsis existed very well in physiology, its relevance came about mainly through the work of Engelmann. ¹ He was the first to investigate the phenomena of peristaltic waves in human ureter. He performed an analysis on physiological or methodological issues which opened new ways for further investigation of physiological mechanisms. After that, Latham ² investigated the mechanism of augmented fluid mixing within a peristaltic pump. After Latham several investigations have been made on the functioning of ureter following the peristaltic motion. Impacts of long wavelengths at low Reynolds number are assumed by Shapiro et al. ³ He modeled the problem and found all elegant analytical results for channel and tube. A number of physically impact and phenomena of bolus, augmented flow and backward flow are discussed. An analysis of peristaltic transport is addressed by considering nonlinear convective acceleration and no slip condition on wavy wall by Fung. ⁴ Hayat et al. ⁵ adopted perturbation method to investigate the peristaltic transport in a planar channel. Shapiro et al. ⁶ modeled a problem to present the simple physical picture of working of a peristaltic pump by means of peristaltic waves in the laboratory frame of reference. Importance of investigations in medical and industrial sciences based on peristaltic motion cannot be ignored. Duane ⁷ examined the impact of peristaltic motion in urinary tract and provide an overview of the basic anatomy and physiology of the urinary tract with an emphasis on their specific roles in host defense. Srivastava et al. ⁸ modeled a problem to investigate the attributes of peristaltic fluid suspension in transport of fluid and revealed that the reversal flow increases with increasing particle concentration and Reynolds number. Yin et al. ⁹ employed a relation between theory and experiments in peristaltic transport. When theoretical and experimental results are compared, they are in good accord and provide confirmation to the theoretical analysis. Zien and Ostrach ¹⁰ investigated the peristaltic transport for a long wave approximation and indicated that mean pressure gradient is associated with the motion. Leal gave a brief study of fluid flow in curved tubes. ¹¹ Hanin ¹² investigated the peristalsis in nano-fluid under endoscopic effects.

Pressure driven flows has demonstrated the ability to control fluids pressure in modern investigation like peristaltic pumping, cell sorting and microfluid injections. In microfluidics, it is often desirable to minimize the circulating volume of fluid. Imposed pressure can regulate the desired volume of fluid flow. Pressure driven flows has been analyzed by.^16,17,18

We observe that in these cases, the flow is generated by the peristaltic waves only through the pressure gradient developed in the channel or tube. However, the question remains; what happens if an external pressure gradient is imposed in addition to the one generated by the peristalsis. We believe that such situation is of much interest both from mathematical and physical points of view that has not been discussed so far.

The imposed pressure gradient is superimposed, and the peristalsis problem is remodeled. The solution is obtained using perturbation technique. The effects and the consequences of the imposed pressure gradient are investigated analytically. Hopefully, this will be a step further in the discussion of peristalsis flows.

Mathematical formulation

Consider a steady two-dimensional boundary layer flow of a viscous and in compressible fluid with imposed pressure gradient along a horizontal channel bounded by two periodic sinusoidal boundaries having separation $2h$ . The fluid is of density $\mathrm{\rho }$ , dynamic viscosity $\mathrm{\mu }$ and with imposed pressure gradient $G$ . The walls of the channel are flexible, on which are traveling, sinusoidal wave of small amplitude. In figure 1, ${\mathrm{y}}_{\mathrm{t}}=1+\epsilon \sin \left(2\pi \overline{x}\right)$ represents the top wavy boundary and $y_B=-(1+\epsilon \sin \left(2\pi \overline{x}\right))$ represents bottom wavy boundary.

Figure 1.

Geometry of fluid flow.

The above figure shows a long channel in which a peristaltic wave of contraction is produced by moving walls to the right at speed $c$ . In peristalsis the sinusoidal wave form is assumed that propagates in one direction with a constant velocity. It is important to mention that the variable boundary give rise to an additional dimension in the modeling of the problem. For the present case, this makes the problem two dimensional instead of one dimensional.

We assume that the pressure gradient $G$ is a nonzero constant and that the upper boundary moves in the same direction as the pressure gradient with a constant velocity $c$ . As the consequence of the variations in the channel width with $x$ , the fluid must accelerate and decelerate, and the pressure gradient will therefore be a function of position. However, we can still impose a mean pressure gradient in the $x$ direction, which we denote as $-G$ . We can then write

$$\frac{\partial p}{\partial x}=-G+\frac{\partial p\prime }{\partial x},$$

(1)

$$\frac{\partial p}{\partial y}=\frac{\partial p\prime }{\partial y}.$$

(2)

Law of Mass conservation:

$$\frac{\partial \overline{u}}{\partial \overline{x}}+\frac{\partial \overline{v}}{\partial \overline{y}}=0,$$

(3)

Law of Momentum:

Momentum equation with imposed pressure gradient G is presented as

In x-direction

$$\mathrm{\rho }(u\frac{\partial u}{\partial x}+v\frac{\partial u}{\partial y})=G-\frac{\partial p\prime }{\partial x}+\mathrm{\mu }(\frac{{\partial }^{2}u}{\partial x^2}+\frac{{\partial }^{2}u}{\partial y^2}),$$

(4)

in y-direction

$$\mathrm{\rho }(u\frac{\partial v}{\partial x}+v\frac{\partial v}{\partial y})=-\frac{\partial p\prime }{\partial y}+\mathrm{\mu }(\frac{{\partial }^{2}v}{\partial x^2}+\frac{{\partial }^{2}v}{\partial y^2}).$$

(5)

In the above equations $u$ and $v$ are velocities in $x$ and $y$ directions respectively and $2h$ is the channel’s width. Here the channel (of plane geometry) has been taken to be infinite in the $y$ plane. The boundaries are located at $y=1/2$ and $y=-1/2$ .

For the present model the relevant boundary conditions are given as

$$\frac{\partial u}{\partial y}=0\mathrm{at}y=0,$$

(6)

$$u=-c,v=0\mathrm{at}y=\pm h.$$

(7)

We introduce the following dimensionless variable transformations

$$\overline{u}=\frac{u}{c_x},\overline{v}=\frac{v}{c_y},$$

(8)

$$\overline{x}=\frac{x}{L},\overline{y}=\frac{y}{d},$$

(9)

$$c_y=c_x\left(\frac{d}{L}\right),\overline{p}=\frac{p\prime }{\mathrm{GL}}.$$

(10)

Using the above transformation, the equation of continuity is comparatively satisfied and the remaining equations (4) and (5) along with boundary conditions (6) and (7) are reduced into the following non-dimensional partial differential equations.

$$\frac{\partial \overline{u}}{\partial \overline{x}}+\frac{\partial \overline{v}}{\partial \overline{y}}=0,$$

(11)

$$\begin{array}{c}\multicolumn{1}{c}{\mathrm{\rho }\prime \frac{c_{x}d}{\mathrm{\mu }}\frac{d}{L}(\overline{u}\frac{\partial \overline{u}}{\partial \overline{x}}+\overline{v}\frac{\partial \overline{v}}{\partial \overline{y}})=\frac{G{d}^2}{\mathrm{\mu }{c}_x}(1-\frac{\partial \overline{p}}{\partial \overline{x}})}\\ +(\frac{d^2}{L^2}\frac{{\partial }^2\overline{u}}{\partial {\overline{x}}^2}+\frac{{\partial }^2\overline{u}}{\partial \overline{y^2}}),\hfill \end{array}$$

(12)

$$\frac{\mathrm{\rho }\prime c_{x}d}{\mathrm{\mu }}\frac{d}{L}(\overline{u}\frac{\partial \overline{v}}{\partial \overline{x}}+\overline{v}\frac{\partial \overline{v}}{\partial \overline{y}})=-\frac{G{L}^2}{\mathrm{\mu }{c}_x}\frac{\partial \overline{p}}{\partial \overline{y}}+\frac{d^2}{L^2}\frac{{\partial }^2\overline{v}}{\partial {\overline{x}}^2}+\frac{{\partial }^2\overline{v}}{\partial \overline{y^2}}.$$

(13)

The notation $\frac{\mathrm{\rho }\prime c_{x}d}{\mathrm{\mu }}$ the right-hand side of the (12) and (13) emphasize that this is a dimensionless parameter known as Reynolds number, multiplies by the geometric parameter $d/L.$

However, we do assume that the Reynolds number is vanishingly small compared even with the geometric amplitude parameter $\epsilon$ , which we shall also assume is small $\epsilon \ll 1$ , showing that the inertial effects can be neglected and that the flow is dominated by viscous and pressure gradient effects. The resulting problem is still of some practical importance.

For $\mathrm{Re}<<1$

$$0=\frac{G{d}^2}{\mathrm{\mu }{c}_x}(1-\frac{p_c}{\mathrm{GL}}\frac{\partial \overline{p}}{\partial \overline{x}})+(\frac{d^2}{L^2}\frac{\partial {\overline{u}}^2}{\partial {\overline{x}}^2}+\frac{{\partial }^2\overline{u}}{\partial \overline{y^2}}),$$

(14)

$$0=-\frac{G{L}^2}{\mathrm{\mu }{c}_x}\frac{\partial \overline{p}}{\partial \overline{y}}+\frac{d^2}{L^2}\frac{{\partial }^2\overline{v}}{\partial {\overline{x}}^2}+\frac{{\partial }^2\overline{v}}{\partial \overline{y^2}},$$

(15)

The transformed boundary conditions are

$$\frac{\partial \overline{u}}{\partial \overline{y}}=0\mathrm{at}\overline{y}=0,$$

(16)

$$\overline{u}=-1,\overline{v}=0\mathrm{at}\overline{y}=\pm h.$$

(17)

Where $h=1+\epsilon \sin \left(2\pi \overline{x}\right)$ .

We see that the dimensionless problem (and its solution) depends on a single dimensionless parameter, $\frac{G{L}^2}{\mathrm{\mu }{c}_x}$ .

This parameter is just the ratio of the two possible velocity scales $\frac{G{L}^2}{\mathrm{\mu }}\mathrm{and}{c}_x$ , one characterized by the magnitude of the pressure gradient and the other by the magnitude of the boundary velocity. The form of the velocity distribution (usually called the velocity profile) depends on the magnitude of the dimensionless ratio of velocities. When $\frac{G{d}^2}{\mathrm{\mu }U}<<1$ , in this case, the fluid motion is dominated by the motion of the boundary and the velocity profile reduces to a linear (simple) shear flow. When $\frac{G{d}^2}{\mathrm{\mu }U}>>1$ , on the other hand, that the quadratic contribution dominates the velocity profile.

Here we consider the domain perturbation method. In this section we seek the solution for limiting case $ϵ<<1$ by means of asymptotic expansions of the form

$$\overline{u}=u_0+\epsilon u_1+O\left({\epsilon }^2\right),$$

(18a)

$$\overline{v}=v_0+\epsilon v_1+O\left({\epsilon }^2\right),$$

(18b)

$$\nabla \overline{p}=\epsilon [G_1(\left\{\overline{x}\right\},\left\{\overline{y}\right\})i+G_2(\left\{\overline{x}\right\},\left\{\overline{y}\right\})j]+O\left({\epsilon }^2\right).$$

(18c)

For this limiting case, the boundary conditions, can also be approximated in terms of asymptotically equivalent boundary conditions applied at $y=\pm 1/2$ .

Applying Taylor’s expansion on the boundary we have:

$$\begin{array}{c}\multicolumn{1}{c}{{u|}_{y=1+\epsilon \sin \left(2\pi \overline{x}\right)}={(u_0+\epsilon u_1+{\epsilon }^2{u}_2)}_{y=\frac{1}{2}}}\\ +\epsilon \sin (2\pi \overline{x})\frac{\partial }{\partial y}{(u_0+\epsilon u_1)}_{y=\frac{1}{2}}\hfill \\ +\frac{1}{2}{(\epsilon \sin (2\pi \overline{x}))}^2{\left(\frac{{\mathrm{d}}^2{u}_0}{\mathrm{d}{y}^2}\right)}_{y=\frac{1}{2}}+\cdots \hfill \end{array}$$

(19)

Using the above relation in equations and comparing the coefficients of $\epsilon$ up to first order, we obtain following set of equations:

Zeroth order equations

$$\frac{\partial u_0}{\partial \overline{x}}+\frac{\partial v_0}{\partial \overline{y}}=0,$$

(20)

$$-\frac{G{d}^2}{\mathrm{\mu }{c}_x}=\frac{{\partial }^2{u}_0}{\partial \overline{y^2}},$$

(21)

$$0=\frac{{\partial }^2{v}_0}{\partial \overline{y^2}}+\frac{d^2}{L^2}\frac{{\partial }^2{v}_0}{\partial \overline{x^2}},$$

(22)

$$u_0=-1,v_0=0\mathrm{at}\overline{y}=\pm \frac{1}{2}.$$

(23)

First order equations

$$\frac{\partial u_1}{\partial \overline{x}}+\frac{\partial v_1}{\partial \overline{y}}=0,$$

(24)

$$\frac{G{d}^2}{\mathrm{\mu }{c}_x}{G}_1(\overline{x},\overline{y})=\frac{{\partial }^2{u}_1}{\partial \overline{y^2}}+\frac{d^2}{L^2}\frac{{\partial }^2{u}_1}{\partial \overline{x^2}},$$

(25)

$$\frac{G{L}^2}{\mathrm{\mu }{c}_x}{G}_2(\overline{x},\overline{y})=\frac{{\partial }^2{v}_1}{\partial \overline{y^2}}+\frac{d^2}{L^2}\frac{{\partial }^2{v}_1}{\partial \overline{x^2}},$$

(26)

$$\begin{array}{c}{u}_1=-\left(\frac{\partial u_0}{\partial \overline{y}}\right)\frac{1}{2}\sin \left(2\pi \overline{x}\right),\hfill \\ |v_1=-\left(\frac{\partial v_0}{\partial \overline{y}}\right)\frac{1}{2}\sin \left(2\pi \overline{x}\right)\left|\mathrm{at}\right|\overline{y}=\pm \frac{1}{2}.\hfill \end{array}$$

(27)

The Reynolds number is $\mathrm{Re}=\mathrm{\rho }\prime \frac{c_{x}d}{\mathrm{\mu }}.$

For the O(1) problem, the boundaries are flat; hence $v_0=0.$

Perturbative solution

In this article we provide the solution of the formulated mathematical model by perturbation method. The focus has been on velocity along the channel. This method applies to problems in which the geometry of the flow or transport domain is irregular in the sense that the boundaries do not correspond to coordinate surfaces of any known analytic coordinate system but are nevertheless near to such coordinate surfaces. The basic idea is to replace the exact boundary condition with the approximate boundary condition that is asymptotically equivalent for $\epsilon <<1$ but now applied at the coordinate surface $y=\pm \frac{1}{2}$ .

Zeroth order solution

$$u_0=-1+\frac{1}{2}\frac{G{d}^2}{\mathrm{\mu }{c}_x}(-\overline{y^2}+\frac{1}{4}).$$

(28)

First order solution

This solution at $O\left(ϵ\right)$ is the first correction for the geometry of sinusoidal wave and $G_1$ must be

$$G_1=\mathrm{Csin}\left(2\pi \overline{x}\right),$$

(29)

To solve this problem, we note that the boundary condition imposes a specific $\overline{x}$ dependence on the solution, which will be preserved throughout the domain. This suggests that we should a solution in the form

$$u_1=F\left(\overline{y}\right)\sin \left(2\pi \overline{x}\right).$$

(30)

It is clear from (28) that the $x$ component of the pressure gradient must have the same dependence on $\overline{x}$ , where we have assumed that the coefficient C is a constant to be determined. We will, of course, need to verify the assumption $C=\mathrm{const}$ is sufficiently general by showing that all of the equations and the boundary conditions for $u_1\mathrm{and}{v}_1$ can be satisfied. The function $F\left(\overline{y}\right)$ must then satisfy the ODE

If we substitute this proposed solution form into equation (24) we have

$$\frac{{\partial }^{2}F\left(\overline{y}\right)\sin \left(2\pi \overline{x}\right)}{\partial \overline{y^2}}+\frac{d^2}{L^2}\frac{{\partial }^{2}F\left(\overline{y}\right)\sin \left(2\pi \overline{x}\right)}{\partial \overline{x^2}}=\frac{G{d}^2}{\mathrm{\mu }{c}_x}\mathrm{Csin}\left(2\pi \overline{x}\right),$$

(31)

$$\begin{array}{c}\sin \left(2\pi \overline{x}\right)\frac{{\partial }^{2}F\left(\overline{y}\right)}{\partial \overline{y^2}}+2\pi \frac{d^2}{L^2}\frac{\partial F\left(\overline{y}\right)\cos \left(2\pi \overline{x}\right)}{\partial \overline{x}}\hfill \\ =\frac{G{d}^2}{\mathrm{\mu }{c}_x}\mathrm{Csin}\left(2\pi \overline{x}\right).\hfill \end{array}$$

(32)

Simplifying equation (31) we obtain

$$\frac{{\partial }^{2}F\left(\overline{y}\right)}{\partial \overline{y^2}}-{\left(2\pi \right)}^2\frac{d^2}{L^2}F\left(\overline{y}\right)=\frac{G{d}^2}{\mathrm{\mu }{c}_x}C,$$

(33)

with boundary conditions

$$\frac{\mathrm{dF}}{d\overline{y}}=0\mathrm{at}\overline{y}=0\mathrm{and}F\left(\overline{y}\right)=\frac{1}{4}\frac{G{d}^2}{\mathrm{\mu }{c}_x}\mathrm{at}\overline{y}=\pm \frac{1}{2}.$$

(34)

The general solution that satisfy the symmetry is

$$F\left(\overline{y}\right)=A\exp \left(\frac{2\pi d\overline{y}}{L}\right)+B\exp \left(\frac{-2\pi d\overline{y}}{L}\right)+C^{*}.\frac{G{d}^2}{\mathrm{\mu }{c}_x},$$

(35)

where $C^{*}=\frac{-C{L}^2}{{\left(2\pi d\right)}^2}$ .

The boundary condition then requires that

$$F=\frac{1}{4}\mathrm{at}\overline{y}=\frac{1}{2}.$$

Hence,

The boundary condition then requires that $F=\frac{1}{4}\mathrm{at}\overline{y}=\frac{1}{2}$ .

Hence,

$$F\left(\overline{y}\right)=\frac{G{d}^2}{\mathrm{\mu }{c}_x}[(\frac{1}{4}-C^{*})\frac{\cosh \left(\frac{2\pi d\overline{y}}{L}\right)}{\cosh (\frac{2\pi d}{L}.\frac{1}{2})}+C^{*}].$$

(36)

The constant $C^{*}$ can be determined from the “mass conservation” constraint that the volume flux must be independent of $\overline{x}$ . This means that

$${\int }_0^{y}12[u_0+\epsilon u_1+O\left({\epsilon }^2\right)]d\overline{y}=\mathrm{const},$$

(37)

that is independent of $\overline{x}$ , because $u_1=F\left(\overline{y}\right)\sin \left[\left(\frac{2\pi d}{L}\right)\overline{x}\right]$ , this implies

$${\int }_0^{\frac{1}{2}}F\left(\overline{y}\right)d\left(\overline{y}\right)=0.$$

(38)

To satisfy this condition, we must determine a specific value for $C^{*}$ such that

$$\overline{u}=u_0+\epsilon u_1+O\left({\epsilon }^2\right),$$

(39)

$$\begin{array}{c}\multicolumn{1}{c}{\overline{u}=-1+\frac{1}{2}\frac{G{d}^2}{\mathrm{\mu }{c}_x}(-\overline{y^2}+\frac{1}{4})+ϵ\left(\sin \left(2\pi \overline{x}\right)\right)}\\ \frac{G{d}^2}{\mathrm{\mu }{c}_x}[(\frac{1}{4}-C^{*})\frac{\cosh \left(\frac{2\pi d\overline{y}}{L}\right)}{\cosh (\frac{2\pi d}{L}.\frac{1}{2})}+C^{*}].\hfill \end{array}$$

(40)

The only remaining issue is to determine the $O\left(ϵ\right)$ contribution to $\overline{v}$ to satisfy the continuity equation.

To solve for $v_1$ we first note from (23) that

$$\frac{\partial \overline{v_1}}{\partial \overline{y}}=-\frac{\partial \overline{u_1}}{\partial \overline{x}},$$

(41)

$$\frac{\partial \overline{v_1}}{\partial \overline{y}}=-\frac{G{d}^2}{\mathrm{\mu }{c}_x}[(\frac{1}{4}-C^{*})\frac{\cosh \left(\frac{2\pi d\overline{y}}{L}\right)}{\cosh (\frac{2\pi d}{L}.\frac{1}{2})}+C^{*}]2\pi \cos \left(2\pi \overline{x}\right).$$

(42)

Hence

$$\begin{array}{c}\overline{v_1}=-\frac{G{d}^2}{\mathrm{\mu }{c}_x}2\pi \cos \left(2\pi \overline{x}\right)\hfill \\ {\int }_0^{\overline{y}}[(\frac{1}{4}-C^{*})\frac{\cosh \left(\frac{2\pi d\overline{y}}{L}\right)}{\cosh (\frac{2\pi d}{L}.\frac{1}{2})}+C^{*}]\mathrm{dy}.\hfill \end{array}$$

(43)

Results

This section deals with the physical interpretation of the analytical results obtained in Section 3. The velocity of the fluid flow is analyzed. The influence of the imposed pressure on the viscous fluid flow in the corrugated channel is particularly of interest. The 3D-graphical illustration of the modification of the axial velocity in the channel is provided in Figures 2 to 8. These graphs describe the behavior of imposed pressure gradient and the velocity profile of peristaltic flow of a Newtonian fluid with different values of $\epsilon$ . Figures 2 and 3 depicts the velocity $u_0$ and $u_1$ . Figures 4 to 8 represents velocity component $u$ in $x$ -direction with different values of perturbation parameter $\epsilon$ . Figure 9 presents velocity $v_1$ of fluid in $y$ direction. Whereas $v_0$ is zero.

Figure 2.

Velocity profile of $u_0$ .

Figure 3.

Velocity profile of $u_1$ .

Figure 4.

Velocity profile $u$ for $ϵ=0$ .

Figure 5.

Velocity profile $u$ for $ϵ=0.025$ .

Figure 6.

Velocity profile u for $ϵ=0.05$ .

Figure 7.

Velocity profile u for $ϵ=0.075$ .

Figure 8.

Velocity profile u for $ϵ=0.1$ .

Figure 9.

Velocity profile $v_1$ in y direction.

The solution that we have generated is quite complicated. The fluid velocity in the x-direction, $u_0+\epsilon u_1$ , decreases where the channel expands and increases where it contracts to maintain constant-volume flux along the channel. For $ϵ=0$ , $v_0=0.$ The flow resistance increases with the corrugation wavenumber and consequently decreases the flow. There is an increase in velocity by increasing external pressure. In general, these results could have a potential application in enhancing fluid transport in physical, biological and engineering processes.