# Effect of the Velocity Second Slip Boundary Condition on the Peristaltic Flow of Nanofluids in an Asymmetric Channel: Exact Solution

Authors: Emad H. Aly; Abdelhalim Ebaid
Journal: Abstract and Applied Analysis
 (2014)

---

## Abstract

The problem of peristaltic nanofluid flow in an asymmetric channel in the presence of the second-order slip boundary condition was investigated in this paper. To the best of the authors’ knowledge, this parameter was here incorporated for the first time in such field of a peristaltic flow. The system governing the current flow was found as a set of nonlinear partial differential equations in the stream function, pressure gradient, nanoparticle concentration, and temperature distribution. Therefore, this system has been successfully solved exactly via a very effective procedure. These exact solutions were then proved to reduce to well-known results in the absence of second slip which were published very recently in the literature. Effect of the second slip parameter on the present physical parameters was discussed through graphs and it was found that this type of slip is a very important one to predict the investigated physical model. Moreover, the variation of many physical parameters such as amplitudes of the lower and upper waves, phase difference on the temperature distribution, nanoparticle concentration, pressure rise, velocity, and pressure gradient were also discussed. Finally, the present results may be viewed as an optimal choice for their dependence on the exact solutions which are obtained due to the highly complex nonlinear system.


---

## Body

## 1. Introduction

Peristalsis is an interest subject and has recently attracted much attention due to its importance in engineering and medical applications. In human body, peristalsis is found in the swallowing food through the esophagus, chyme motion in the gastrointestinal tract, vasomotion of small blood vessels such as venules, capillaries and arterioles, urine transport from kidney to bladder, and intrauterine fluid flow within the uterine cavity. Various experimental and theoretical studies have been conducted to understand peristaltic flow in asymmetric channel or an axisymmetric tube; see [1–7]. According to De Vries et al. [8], they observed that myometrial contractions are peristaltic-type motion and therefore these contractions of the uterine wall may occur in both symmetric and asymmetric directions. Due to the importance of this field of research, a great effort was devoted to study this type of flow for Newtonian and non-Newtonian fluids in an asymmetric channel [9–17].Recently, a few papers have been published in the field of nanofluid flows under the peristaltic action [18–23]. Nanofluids usually lead to the enhancement of the thermal conductivity of the base fluid [24]. The concept of nanofluids, analogous to that of nanoparticles, may be initiated by Choi et al. [25]. The fact that nanofluids have higher thermal conductivity than the other heat transfer fluids because of their nanostructure has attracted many engineers and theoretical scientists to investigate their behavior. This higher thermal conductivity may be very useful in the treatment of tumors by injecting the blood vessel nearest to the tumor with magnetic nanoparticles along with placing a magnet close to the tumor, where these particles act like heat sources, in the presence of the applied magnetic field of alternating nature. Hence, the undesirable tissues (cancer's tissues) are destroyed when the temperature reaches 42–45°C [6]. Besides, the drug may be placed on the magnetic nanoparticles and is injected near the tumor. The drug is then absorbed by the tumor through a high gradient magnetic field, which is concentrated near the tumor center [26]. In fact, the drug absorption due to high concentration of magnetic particles increases and magnetic force prevents uniform drug distribution in circulatory system. This approach reduces the side effect and allows using high dose of anticancer drug [27].Very recently, Ebaid and Aly [28] investigated the system of partial differential equations describing the peristaltic flow of a nanofluid with slip effect of the velocity, temperature, and concentration has been analytically solved. The obtained exact solutions have been applied to study effects of the slip parameter, thermophoresis, Brownian motion parameters, and many other parameters on the pressure rise, velocity profiles, temperature distribution, nanoparticle concentration, and pressure gradient. These exact solutions were proved to reduce to the results in the literature at special cases when β2=0. Moreover, on comparing the present results with those in [22], remarkable differences were noticed for behavior of the included previous physical phenomena.The main feature of the previous studies on the nanofluid flow under peristaltic action [18–23] is the ignorance of the second slip effect, while this effect has been recently discussed by many authors; see for example [29–33]. These papers showed clearly that the second-order slip flow model is necessary to predict the flow characteristics accurately. Therefore, as an extension of these studies and for the first time, in this paper we aim to investigate effect of the second slip on the peristaltic flow of nanofluids in an asymmetric channel. Hence, we focus here on the problem discussed by Akbar et al. [22] and Ebaid and Aly [28] with slight differences in the boundary conditions.


## 2. The Mathematical Model

Consider peristaltic transport of an incompressible Newtonian nanofluid in an asymmetric channel with flexible walls. The channel asymmetry is generated by propagation of waves on the channel walls traveling with different amplitudes and phases but with the same constant speedc. In the Cartesian coordinates system (X¯,Y¯) of the fixed frame, the upper and lower walls h1¯ and h2¯, respectively, are given by (see Figure 1)(1a)h1¯=d1+a1cos⁡(2πλ[X¯-ct¯]),(1b)h2¯=-d2-b1cos⁡(2πλ[X¯-ct¯]+ϕ),Figure 1
Schematic diagram of a two-dimensional asymmetric channel.wherea1 and b1 are amplitude of the waves, λ is the wave length, and d1+d2 is the width of the channel. The phase difference ϕ varies in the range 0≤ϕ≤π, where ϕ=0 and ϕ=π correspond to symmetric channel with waves out of the phase and in the phase, respectively. Further, a1, b1, d1, d2, and ϕ have to satisfy the following condition [34]:
(2)a12+b12+2a1b1cos⁡ϕ≤(d1+d2)2,
with the following nondimensional phenomena [22]:
(3)a=a1d1,b=a2d1,d=d2d1.
Furthermore, in the moving frame of references (x¯,y¯), we have
(4)x¯=X¯-ct¯,y¯=Y¯,
which have been used with the following nondimensional phenomena [22]:
(5)x=2πx¯λ,y=y¯d1,t=2πt¯λ,h1=h1¯d1,h2=h2¯d2,
to obtain the present physical model. On taking into account (i) heat transfer along with nanoparticle phenomena, (ii) long wavelength, and low Reynolds number approximation, recently, Akbar et al. [22] found that the flow is governed by the following system of partial differential equations in nondimensional form:
(6)∂4ψ∂y4+Gr∂θ∂y+Br∂σ∂y=0,(7)dpdx=∂∂y[∂2ψ∂y2+Grθ+Brσ],(8)∂2θ∂y2+Nb∂θ∂y∂σ∂y+Nt[∂θ∂y]2=0,(9)∂2σ∂y2+[NtNb]∂2θ∂y2=0,
where ψ, θ, σ, and p are the stream function, temperature distribution, nanoparticle concentration, and pressure gradient, respectively. In addition, Nb, Nt, Gr, and Br are the Brownian motion parameter, thermophoresis parameter, local temperature Grashof number, and nanoparticle Grashof number, respectively. The above system of PDEs. (6)–(9) has to be solved subject to the following boundary conditions:
(10)ψ=F2,∂ψ∂y=-β1∂2ψ∂y2-β2∂3ψ∂y3-1ath1=1+acos⁡(x),ψ=-F2,∂ψ∂y=β1∂2ψ∂y2+β2∂3ψ∂y3-1ath2=-d-bcos⁡(x+ϕ),θ+γ∂θ∂y=0,aty=h1,θ-γ∂θ∂y=1,aty=h2,σ+γ1∂σ∂y=0,aty=h1,σ-γ1∂σ∂y=1,aty=h2,
where β1, β2, γ, and γ1 represent the first-order slip parameter, second-order slip parameter, thermal slip parameter, and concentration slip parameter, respectively. It should be noted that the full formulation of the above model can be checked in [22], and there is no need to repeat it again here. However, as mentioned in Section 1 in the present work, we focused on obtaining the general closed form solution in the presence of the second-order slip parameter, as a very important one to predict the investigated physical model.


## 3. The General Closed Form Solution

As mentioned by Ebaid and Aly [28], we have the following exact solutions for the temperature distribution and nanoparticle concentration:
(11)θ(x,y)=f4e-Nbf1y+1Nbf3f1,σ(x,y)=-NtNbf4e-Nbf1y+f1y+f2-NtNb2f3f1,
where
(12)f2=NtNb2f3f1-(γ1f1-1Nb)Ntf4r1f1-(γ1+h1)f1,f3=-Nbf1(1-γNbf1)r1f1(1+γNbf1)r2f1-(1-γNbf1)r1f1,f4=1(1+γNbf1)r2f1-(1-γNbf1)r1f1,
where
(13)r1=e-Nbh1,r2=e-Nbh2.f1 can be obtained from the following implicit algebraic equation:
(14)NtNb[(γ1Nbf1-1)r1f1+(γ1Nbf1+1)r2f1(γNbf1-1)r1f1+(γNbf1+1)r2f1]+(2γ1+h1-h2)f1=-1.
It can be seen from (14) that it is of complex structure to be solved exactly in its general form. However, it has an exact solution for f1 when γ=γ1 and given by
(15)f1=-(1+Nt/Nb)(2γ1+h1-h2).
Regarding the stream function, it can be solved to give
(16)ψ=f8+f7y+12f6y2+16f5y3+g(y),
where
(17)g(y)=16Ω1y3-124Brf1y4-Ω2(Nbf1)3e-Nbf1y.
Applying the boundary conditions on the ψ-equation, we obtain the following system:(18a)f8+f7h1+12f6h12+16f5h13=R1,(18b)f8+f7h2+12f6h22+16f5h23=R2,(18c)f7+(β1+h1)f6+(12h12+β1h1)f5=S1,(18d)f7+(h2-β1)f6+(12h22-β1h2-β2)f5=S2,where(19a)R1=F2-g(h1),(19b)R2=-F2-g(h2),(19c)S1=-1-g′(h1)-β1g′′(h1)-β2g′′′(h1),(19d)S2=-1-g′(h2)+β1g′′(h2)+β2g′′′(h2).On solving the system (18a)–(18d), we get
(20)f5=6(-2R1+2R2+(h1-h2)(S1+S2))(h1-h2)2(6β1+h1-h2),f6=(2(-h13(S1+2S2)+6h12(R1-R2+S2(h2-2β1))+2h2(-h2(3R1-3R2+h2(2S1+S2))+6(R1-R2+h2S1)β1)+6h1(h22S1+2(R1-R2+h2(-S1+S2))β1)+12(2R1-2R2-(h1-h2)(S1-S2))β22(-h13(S1+2S2)+6h12(R1-R2+S2(h2-2β1))))×(((h1-h2)2(h1-h2+2β1)(h1-h2+6β1)))-1,f7=(+6h2((R1-R2)β1-(S1+S2)β2))h14S2+2h13(h2S1-(S1-2S2)β1)+2h1(+6(R1-R2)(β12-β2)h23S2+12h2(-R1+R2)β1+3h22(R1-R2-S2β1)+6(R1-R2)(β12-β2))+h12(-3h2(2R1-2R2+h2(S1+S2))+6(R1-R2+h2S1)β1+6(S1+S2)β2)+h2(+6h2((R1-R2)β1-(S1+S2)β2)h23S1+2h22(-2S1+S2)β1-12(R1-R2)×(β12+β2)+6h2((R1-R2)β1-(S1+S2)β2)))×((h1-h2)2(h1-h2+2β1)(h1-h2+6β1))-1,f8=(h14(R2-h2S2)-h13(4R2+h2(S1-S2))(h2-2β1)+h22R1(h2-6β1)(h2-2β1)+h12(h23(2S1-S2)+12R2β12+3h22(R1+R2+2(-S1+S2)β1)+6h2(-(R1+3R2)β1-(S1+S2)β2h23(2S1-S2)+12R2β12))+h1h2(-h23S1-12(R1+R2)β12-2h22(2R1+(-2S1+S2)β1)+12(R1-R2)β2+6h2((3R1+R2)β1+(S1+S2)β2)))×((h1-h2)2(h1-h2+2β1)(h1-h2+6β1))-1.
To get the pressure gradient dp/dx, we obtain from (7) that
(21)dpdx=Ω3-Brf1y+(1+Nbf1)Ω2e-Nbf1y,
where
(22)Ω3=Ω1+f5+Brf1,Ω2=(BrNtNb-Gr)f4,Ω1=(BrNtNb-Gr)1Nbf3f1-Brf2.
The pressure rise Δp in terms of the flow rate Q is given as follows:
(23)Δp=∫01‍(dpdx)dx.
Therefore, the exact expression for the pressure rise at the center of the channel is given as
(24)Δp=∫01[6[2g(h1)-2g(h2)+(h1-h2)(S1+S2)](h1-h2)2(6β1+h1-h2)+Ω1+Brf1(1-h1+h22)+(1+Nbf1)Ω2e-Nbf1((h1+h2)/2)6[2g(h1)-2g(h2)+(h1-h2)(S1+S2)](h1-h2)2(6β1+h1-h2)]dx-12(Q-1-d)∫01‍1(h1-h2)2(6β1+h1-h2)dx.


## 4. Results and Discussion

This section is devoted for investing the exact solutions which have been obtained in the previous section for the temperature distribution, nanoparticle concentration, velocity, pressure gradient, and pressure rise. These solutions are expected to provide us with the correct physical effect of the second slip for the five investigated physical phenomena. It is observed from (11) that there is no effect for the second slip parameter on the solutions of the temperature distribution and nanoparticle concentration. Although (11) have been used by Ebaid and Aly [28] to discuss the effect of many parameters on the temperature distribution and nanoparticle concentration for comparing with the results in [22], many other parameters were ignored. These important parameters are investigated in the next subsections.


### 4.1. The Temperature Distributionθ

Very recently in [28], Ebaid and Aly discussed the effect of Nt, γ, and γ1 on the temperature distribution. Here, we aim to discuss the effect of many other parameters on this phenomenon. In Figures 2(a) and 2(b), effect of the amplitude a on the temperature distribution is depicted at two different values for the Browning motion parameter and fixed values for the other investigated parameters. From these figures, it can be concluded that there are no remarkable influences for the variation of a on θ in the first part of the channel, while an increase in θ occurs in the upper part with increasing a. The variation of the amplitude b on the temperature distribution is presented in Figures 3(a) and 3(b). These figures show that very small decrease in θ occurs as b increases. This means that while the amplitude a has remarkable effect on θ, the influence of the amplitude b on θ may be ignored. In Figures 4(a) and 4(b), the variation of the phase difference ϕ on θ is introduced. Although the increase in ϕ between the waves of the channel increases the temperature distribution, this increase is very small. Instead of studying the variation of θ at separate values of Brownian motion parameter Nb, we present in Figure 5(a) the variation of θ against a continuous range for the values of Nb. This figure indicates a clear description for the variation of Brownian motion parameter Nb on θ. In addition, Figure 5(b) shows other descriptions for the variation of the thermophoresis parameter Nt on the temperature distribution. It can be concluded from the last two figures that θ always decreases with the increase in Nb and Nt.Variation of temperature profile at different values ofa when Nt=5, d=1, b=0.5, x=1, ϕ=0.2, γ=0.1, and γ1=0.1 for (a) Nb=0.8 and (b) Nb=3.
(a)
(b)Variation of temperature profile at different values ofb when Nt=5, d=1, a=0.5, x=1, ϕ=0.2, γ=0.3, and γ1=0.3 for (a) Nb=0.8 and (b) Nb=3.
(a)
(b)Variation of temperature profile at different values ofϕ when Nt=5, d=1, a=0.5, x=1, b=0.5, γ=0.7, and γ1=0.7 for (a) Nb=0.8 and (b) Nb=3.
(a)
(b)Variation of temperature profile at different values of (a)Nb for Nt=0.8 and (b) Nt for Nb=0.8; both cases when d=1, a=0.5, b=0.5, ϕ=0.2, x=1, γ=0.1, and γ1=0.1.
(a)
(b)


### 4.2. The Nanoparticle Concentrationσ

The nanoparticles concentrationσ is depicted in Figures 6(a) and 6(b). It is observed from these figures that σ reaches its highest value at the lower wall of the channel and then it decreases in a certain domain across the channel. A converse of this behavior occurs after that domain. It is also detected that at a higher value of Brownian motion parameter Nb the domain in which σ decreases becomes wider. In addition, the wide of such domain increases with increasing the amplitude of the upper wave. The variation of the lower wave amplitude on σ is presented in Figures 7(a) and 7(b). For all values of b, it is noticed that σ decreases until it reaches a certain point. After that point the curves become identical whatever the value of b becomes. This refers to that there is no effect for b on σ after this point.Variation of nanoparticle concentration at different values ofa when Nt=5, d=1, b=0.5, x=1, ϕ=0.2, γ=0.1, and γ1=0.1 for (a) Nb=0.8 and (b) Nb=3.
(a)
(b)Variation of nanoparticle concentration at different values ofb when Nt=5, d=1, a=0.5, x=1, ϕ=0.2, γ=0.3, and γ1=0.3 for (a) Nb=0.8 and (b) Nb=3.
(a)
(b)In Figures8(a) and 8(b), the variation of the phase difference ϕ on σ is displayed. Although the increase in the phase difference between the waves of the channel increases the nanoparticles concentration, this occurs in certain domain of the channel. However, after that domain the curves are identical whatever the value of the phase difference. Moreover, instead of studying the variation of σ at separate values of Brownian motion parameter Nb, the variation of σ against a continuous range for the values of Nb is displayed in Figure 9(a). In addition, Figure 9(b) shows other descriptions for the variation of the thermophoresis parameter Nt on the nanoparticle concentration. It can be concluded from Figure 9(a) that σ always decreases when Nb is in the range from 1 to 10 at fixed values for the other parameters. However, the situation is different for the variation of σ against y at certain values of Nt, where it is noticed that when 1<Nt<4 then σ decreases within the channel. In addition, at values of Nt higher than 4, we observe that σ decreases in certain domain of the channel and changes this behavior in the rest of channel.Variation of nanoparticle concentration at different values ofϕ when Nt=5, d=1, a=0.5, x=1, b=0.5, γ=0.7, and γ1=0.7 for (a) Nb=0.8 and (b) Nb=3.
(a)
(b)Variation of nanoparticle concentration at different values of (a)Nb for Nt=0.8 and q=0.5 and (b) Nt for Nb=0.8 and b=0.5; both cases when d=1, a=0.5, ϕ=0.2, x=1, γ=0.1, and γ1=0.1.
(a)
(b)


### 4.3. The Pressure RiseΔp

The exact expression for the pressure rise is given by (24) which is valid for any set of the physical parameters. On using this equation, we observe that the pressure rise is always a decreasing function in terms of the flow rate. This notice is clarified in Figures 10(a)–10(d). It is also observed from Figure 10(a) that the pressure rise increases with increasing the negativity of the second slip parameter β2 in the whole range of the flow rate Q. The effect of the thermophoresis parameter Nt on the pressure rise is depicted in Figure 10(b). It can be seen from this figure that the pressure rise decreases with the increase in the values of Nt in the whole region of Q.Variation of the pressure rise versus flow rate ata=0.3, d=1, and γ=γ1=0.5 for (a) b=0.5, Nb=0.8, Nt=0.5, Br=0.5, Gr=0.2, ϕ=0.2, and β1=0.1; (b) b=0.5, Nb=0.5, Br=0.5, Gr=0.5, ϕ=0.2, β1=0.8, and β2=-2; (c) Nb=0.8, Nt=0.5, Br=0.5, Gr=0.2, ϕ=0.5, β1=0.8, and β2=-2; and (d) b=0.5, Nb=0.8, Nt=0.5, Gr=0.2, ϕ=0.2, β1=0.8, and β2=-2.
(a)
(b)
(c)
(d)Figure10(c) shows that the pressure rise decreases with decreasing the amplitude b when Q lies in the region -2.5<Q<1, while a converse behavior occurs for Q>1. Figure 10(d) reveals that the pressure rise increases with increasing the local nanoparticle Grashof number Br. This behavior differs from the results obtained by Akbar et al. [22] in the absence of the second slip effect. Therefore, the present results reveal that the existence of this parameter has a significant effect on the pressure rise and leads to different behavior when compared with its ignorance.


### 4.4. The Velocityu(y)

Regarding the axial velocity, it was plotted in Figures11(a)–11(d) at several values of the present physical parameters. Variations of the second slip parameter β2, Brownian motion parameter Nb, flow rate Q, and local temperature Grashof number Gr on the velocity profile have been displayed in Figures 11(a)–11(d), respectively. It is found from Figure 11(a) that the velocity increases near the lower wall with increasing β2, while a converse behavior occurs near the upper wall; that is, the velocity decreases with increasing β2. Unlike the study made by Akbar et al. [22], the effect of Brownian motion parameter Nb leads to that there are three different regions in which alternating behavior for the velocity occurs. As shown in Figure 11(b), in a central region of the channel, that is, -0.6<y<0.8, the velocity increases with the increase in Nb, while in the other two regions of the channel, that is, h2<y<-0.6 and 0.8<y<h1, a different situation is detected. At the same three regions just mentioned, it can be seen from Figure 11(d) that effect of the local temperature Grashof number Gr on the velocity profiles happens in a converse manner regarding Brownian motion parameter Nb. However, it is clear from Figure 11(c) that an increase in the flow rate leads to an increase in the velocity and this is similar to the variation of the velocity profile against flow rate as in Akbar et al. [22].Variation of the velocity profile ata=0.1, b=0.5, d=1, Br=0.5, x=0.1, ϕ=0.2, γ=γ1=0.5, and β1=0.1 for (a) Nt=0.5, Nb=0.8, Gr=0.5, Q=2, and β2=-0.2,-0.4,-0.6,-0.8; (b) Nt=0.7, Gr=3, Q=2, β2=-0.6, and Nb=3,4,5,6; (c) Nt=0.5, Nb=0.8, Gr=0.5, β2=-0.2, and Q=0.2,0.4,0.6,0.8; and (d) Nt=0.5, Nb=0.8, Q=2, β2=-0.2, and Gr=1,2,3,4.
(a)
(b)
(c)
(d)


### 4.5. The Pressure Gradientdp/dx

Figures12(a)–12(d) indicate the pressure gradient for different values of β2, Nb, Gr, and a, respectively. From these figures, it is shown that magnitude of the pressure gradient increases in view of an increase in β2, Gr, and a, and it decreases with an increase in Nb. In addition, it is observed that our results on the effect of β2, Nb, and Gr on the pressure gradient are completely different comparing with those obtained in [22] for the effect of the first slip parameter only. This also may refer to the importance of including the second slip parameter β2 in such problems of peristaltic flow. However, the effect of the amplitude a on the pressure gradient remains the same even in the presence of β2, where Figure 12(d) shows that an increase in the amplitude a leads to an increase in the pressure gradient. Besides, one of the main features of these graphs is that the maximum pressure gradient occurs at the same point when x=0.5. In the absence of β2, this value was shifted slightly to the left; that is, it occurred at x=0.45 as in [22]. The final notice here is that the present numerical discussion is based on the obtained exact analytical solutions for the current physical model, unlike the approximate series solutions obtained in [22] whose disadvantages have been proved by Ebaid and Aly [28].Variation of the pressure gradient atb=0.5, d=1, Nt=0.5, Br=0.5, ϕ=0.2, Q=0.5, and γ=γ1=0.5 for (a) β1=0.1, Nb=0.8, Gr=0.5, a=0.3, and β2=-0.5,-1,-1.5,-2; (b) β1=0.8, β2=-2, Gr=0.5, a=0.3, and Nb=1,5,9,13; (c) β1=0.8, β2=-2, Nb=0.8, a=0.3, and Gr=0.5,0.6,0.7,0.8; and (d) β1=0.3, β2=-2, Nb=0.8, Gr=0.5, and a=0.2,0.3,0.4,0.5.
(a)
(b)
(c)
(d)


## 4.1. The Temperature Distributionθ

Very recently in [28], Ebaid and Aly discussed the effect of Nt, γ, and γ1 on the temperature distribution. Here, we aim to discuss the effect of many other parameters on this phenomenon. In Figures 2(a) and 2(b), effect of the amplitude a on the temperature distribution is depicted at two different values for the Browning motion parameter and fixed values for the other investigated parameters. From these figures, it can be concluded that there are no remarkable influences for the variation of a on θ in the first part of the channel, while an increase in θ occurs in the upper part with increasing a. The variation of the amplitude b on the temperature distribution is presented in Figures 3(a) and 3(b). These figures show that very small decrease in θ occurs as b increases. This means that while the amplitude a has remarkable effect on θ, the influence of the amplitude b on θ may be ignored. In Figures 4(a) and 4(b), the variation of the phase difference ϕ on θ is introduced. Although the increase in ϕ between the waves of the channel increases the temperature distribution, this increase is very small. Instead of studying the variation of θ at separate values of Brownian motion parameter Nb, we present in Figure 5(a) the variation of θ against a continuous range for the values of Nb. This figure indicates a clear description for the variation of Brownian motion parameter Nb on θ. In addition, Figure 5(b) shows other descriptions for the variation of the thermophoresis parameter Nt on the temperature distribution. It can be concluded from the last two figures that θ always decreases with the increase in Nb and Nt.Variation of temperature profile at different values ofa when Nt=5, d=1, b=0.5, x=1, ϕ=0.2, γ=0.1, and γ1=0.1 for (a) Nb=0.8 and (b) Nb=3.
(a)
(b)Variation of temperature profile at different values ofb when Nt=5, d=1, a=0.5, x=1, ϕ=0.2, γ=0.3, and γ1=0.3 for (a) Nb=0.8 and (b) Nb=3.
(a)
(b)Variation of temperature profile at different values ofϕ when Nt=5, d=1, a=0.5, x=1, b=0.5, γ=0.7, and γ1=0.7 for (a) Nb=0.8 and (b) Nb=3.
(a)
(b)Variation of temperature profile at different values of (a)Nb for Nt=0.8 and (b) Nt for Nb=0.8; both cases when d=1, a=0.5, b=0.5, ϕ=0.2, x=1, γ=0.1, and γ1=0.1.
(a)
(b)


## 4.2. The Nanoparticle Concentrationσ

The nanoparticles concentrationσ is depicted in Figures 6(a) and 6(b). It is observed from these figures that σ reaches its highest value at the lower wall of the channel and then it decreases in a certain domain across the channel. A converse of this behavior occurs after that domain. It is also detected that at a higher value of Brownian motion parameter Nb the domain in which σ decreases becomes wider. In addition, the wide of such domain increases with increasing the amplitude of the upper wave. The variation of the lower wave amplitude on σ is presented in Figures 7(a) and 7(b). For all values of b, it is noticed that σ decreases until it reaches a certain point. After that point the curves become identical whatever the value of b becomes. This refers to that there is no effect for b on σ after this point.Variation of nanoparticle concentration at different values ofa when Nt=5, d=1, b=0.5, x=1, ϕ=0.2, γ=0.1, and γ1=0.1 for (a) Nb=0.8 and (b) Nb=3.
(a)
(b)Variation of nanoparticle concentration at different values ofb when Nt=5, d=1, a=0.5, x=1, ϕ=0.2, γ=0.3, and γ1=0.3 for (a) Nb=0.8 and (b) Nb=3.
(a)
(b)In Figures8(a) and 8(b), the variation of the phase difference ϕ on σ is displayed. Although the increase in the phase difference between the waves of the channel increases the nanoparticles concentration, this occurs in certain domain of the channel. However, after that domain the curves are identical whatever the value of the phase difference. Moreover, instead of studying the variation of σ at separate values of Brownian motion parameter Nb, the variation of σ against a continuous range for the values of Nb is displayed in Figure 9(a). In addition, Figure 9(b) shows other descriptions for the variation of the thermophoresis parameter Nt on the nanoparticle concentration. It can be concluded from Figure 9(a) that σ always decreases when Nb is in the range from 1 to 10 at fixed values for the other parameters. However, the situation is different for the variation of σ against y at certain values of Nt, where it is noticed that when 1<Nt<4 then σ decreases within the channel. In addition, at values of Nt higher than 4, we observe that σ decreases in certain domain of the channel and changes this behavior in the rest of channel.Variation of nanoparticle concentration at different values ofϕ when Nt=5, d=1, a=0.5, x=1, b=0.5, γ=0.7, and γ1=0.7 for (a) Nb=0.8 and (b) Nb=3.
(a)
(b)Variation of nanoparticle concentration at different values of (a)Nb for Nt=0.8 and q=0.5 and (b) Nt for Nb=0.8 and b=0.5; both cases when d=1, a=0.5, ϕ=0.2, x=1, γ=0.1, and γ1=0.1.
(a)
(b)


## 4.3. The Pressure RiseΔp

The exact expression for the pressure rise is given by (24) which is valid for any set of the physical parameters. On using this equation, we observe that the pressure rise is always a decreasing function in terms of the flow rate. This notice is clarified in Figures 10(a)–10(d). It is also observed from Figure 10(a) that the pressure rise increases with increasing the negativity of the second slip parameter β2 in the whole range of the flow rate Q. The effect of the thermophoresis parameter Nt on the pressure rise is depicted in Figure 10(b). It can be seen from this figure that the pressure rise decreases with the increase in the values of Nt in the whole region of Q.Variation of the pressure rise versus flow rate ata=0.3, d=1, and γ=γ1=0.5 for (a) b=0.5, Nb=0.8, Nt=0.5, Br=0.5, Gr=0.2, ϕ=0.2, and β1=0.1; (b) b=0.5, Nb=0.5, Br=0.5, Gr=0.5, ϕ=0.2, β1=0.8, and β2=-2; (c) Nb=0.8, Nt=0.5, Br=0.5, Gr=0.2, ϕ=0.5, β1=0.8, and β2=-2; and (d) b=0.5, Nb=0.8, Nt=0.5, Gr=0.2, ϕ=0.2, β1=0.8, and β2=-2.
(a)
(b)
(c)
(d)Figure10(c) shows that the pressure rise decreases with decreasing the amplitude b when Q lies in the region -2.5<Q<1, while a converse behavior occurs for Q>1. Figure 10(d) reveals that the pressure rise increases with increasing the local nanoparticle Grashof number Br. This behavior differs from the results obtained by Akbar et al. [22] in the absence of the second slip effect. Therefore, the present results reveal that the existence of this parameter has a significant effect on the pressure rise and leads to different behavior when compared with its ignorance.


## 4.4. The Velocityu(y)

Regarding the axial velocity, it was plotted in Figures11(a)–11(d) at several values of the present physical parameters. Variations of the second slip parameter β2, Brownian motion parameter Nb, flow rate Q, and local temperature Grashof number Gr on the velocity profile have been displayed in Figures 11(a)–11(d), respectively. It is found from Figure 11(a) that the velocity increases near the lower wall with increasing β2, while a converse behavior occurs near the upper wall; that is, the velocity decreases with increasing β2. Unlike the study made by Akbar et al. [22], the effect of Brownian motion parameter Nb leads to that there are three different regions in which alternating behavior for the velocity occurs. As shown in Figure 11(b), in a central region of the channel, that is, -0.6<y<0.8, the velocity increases with the increase in Nb, while in the other two regions of the channel, that is, h2<y<-0.6 and 0.8<y<h1, a different situation is detected. At the same three regions just mentioned, it can be seen from Figure 11(d) that effect of the local temperature Grashof number Gr on the velocity profiles happens in a converse manner regarding Brownian motion parameter Nb. However, it is clear from Figure 11(c) that an increase in the flow rate leads to an increase in the velocity and this is similar to the variation of the velocity profile against flow rate as in Akbar et al. [22].Variation of the velocity profile ata=0.1, b=0.5, d=1, Br=0.5, x=0.1, ϕ=0.2, γ=γ1=0.5, and β1=0.1 for (a) Nt=0.5, Nb=0.8, Gr=0.5, Q=2, and β2=-0.2,-0.4,-0.6,-0.8; (b) Nt=0.7, Gr=3, Q=2, β2=-0.6, and Nb=3,4,5,6; (c) Nt=0.5, Nb=0.8, Gr=0.5, β2=-0.2, and Q=0.2,0.4,0.6,0.8; and (d) Nt=0.5, Nb=0.8, Q=2, β2=-0.2, and Gr=1,2,3,4.
(a)
(b)
(c)
(d)


## 4.5. The Pressure Gradientdp/dx

Figures12(a)–12(d) indicate the pressure gradient for different values of β2, Nb, Gr, and a, respectively. From these figures, it is shown that magnitude of the pressure gradient increases in view of an increase in β2, Gr, and a, and it decreases with an increase in Nb. In addition, it is observed that our results on the effect of β2, Nb, and Gr on the pressure gradient are completely different comparing with those obtained in [22] for the effect of the first slip parameter only. This also may refer to the importance of including the second slip parameter β2 in such problems of peristaltic flow. However, the effect of the amplitude a on the pressure gradient remains the same even in the presence of β2, where Figure 12(d) shows that an increase in the amplitude a leads to an increase in the pressure gradient. Besides, one of the main features of these graphs is that the maximum pressure gradient occurs at the same point when x=0.5. In the absence of β2, this value was shifted slightly to the left; that is, it occurred at x=0.45 as in [22]. The final notice here is that the present numerical discussion is based on the obtained exact analytical solutions for the current physical model, unlike the approximate series solutions obtained in [22] whose disadvantages have been proved by Ebaid and Aly [28].Variation of the pressure gradient atb=0.5, d=1, Nt=0.5, Br=0.5, ϕ=0.2, Q=0.5, and γ=γ1=0.5 for (a) β1=0.1, Nb=0.8, Gr=0.5, a=0.3, and β2=-0.5,-1,-1.5,-2; (b) β1=0.8, β2=-2, Gr=0.5, a=0.3, and Nb=1,5,9,13; (c) β1=0.8, β2=-2, Nb=0.8, a=0.3, and Gr=0.5,0.6,0.7,0.8; and (d) β1=0.3, β2=-2, Nb=0.8, Gr=0.5, and a=0.2,0.3,0.4,0.5.
(a)
(b)
(c)
(d)


## 5.  Conclusion

The effect of the second slip condition on the peristaltic transport of a nanofluid in an asymmetric channel was investigated for the first time in the present paper. The system describing the problem was exactly solved via a simple, but very effective, analytical procedure. Based on the obtained exact solutions, some important results were introduced through graphs to indicate effect of the second slip on the velocity profiles, pressure rise, and pressure gradient. Moreover, the variation of many physical parameters such as amplitudes of the lower and upper waves, phase difference on the temperature distribution, and nanoparticle concentration were also discussed for a new set of the presented parameters rather than those in the literature. It was found that the temperature distributionθ always decreases with the increase in Nb and Nt. In addition, the amplitude a has remarkable effect on θ, while the influence of the amplitude b on θ may be ignored. For all values of b, nanoparticles concentration σ decreases until it reaches a certain point; however, there is no effect for b after this point. In a certain domain of the channel, the phase difference between the waves of the channel increases σ; however, the curves are identical whatever value of the phase difference after that. Finally, it should be mentioned that the second slip condition had a significant effect on the various phenomena involved in the current model. Accordingly, the present suggested technique may be useful for obtaining exact solutions for many other similar problems. In addition, such type of the second slip should be included for modelling the peristaltic flow.



---

*Source: 191876-2014-03-23.xml*