# The Effect of Boundary Slip on the Transient Pulsatile Flow of a Modified Second-Grade Fluid

Authors: N. Khajohnsaksumeth; B. Wiwatanapataphee; Y. H. Wu
Journal: Abstract and Applied Analysis
 (2013)

---

## Abstract

We investigate the effect of boundary slip on the transient pulsatile fluid flow through a vessel with body acceleration. The Fahraeus-Lindqvist effect, expressing the fluid behavior near the wall by the Newtonian fluid while in the core by a non-Newtonian fluid, is also taken into account. To describe the non-Newtonian behavior, we use the modified second-grade fluid model in which the viscosity and the normal stresses are represented in terms of the shear rate. The complete set of equations are then established and formulated in a dimensionless form. For a special case of the material parameter, we derive an analytical solution for the problem, while for the general case, we solve the problem numerically. Our subsequent analytical and numerical results show that the slip parameter has a very significant influence on the velocity profile and also on the convergence rate of the numerical solutions.


---

## Body

## 1. Introduction

In this paper, we study a fluid-structure interaction problem, namely, the effect of boundary slip on the flow of a non-Newtonian fluid through microchannels. This problem has many applications, and in this paper we particularly focus on blood flow in the cardiovascular system.For the study of blood flow in arteries, two major types of constitutive models have been used. The first type of models is based on the microcontinuum or the structured continuum theories [1–6] in which the balance laws are used to determine the characteristics of blood motion. In the other type of models, blood is considered as a suspension, and its flow is modeled by the non-Newtonian fluid mechanics. Due to the red blood cells (RBC) migration as shown experimentally, blood has been modeled as a two-stage fluid by many researchers [7–9]. The first stage is a peripheral layer which is modeled as a Newtonian viscous fluid, while the other one is a centre core which is modeled as a non-Newtonian fluid. The effect of body acceleration and pulsatile conditions were taken into account under the same problem by Majhi et al. [7, 10]. Later, Massoudi and Phuoc [11] used the (generalized) second-grade fluid constitutive model to describe the shear thinning and normal stress effect, and the behavior of blood flow near the wall is modeled by the Newtonian fluid model, while the behavior of the blood flow at the core is described by the second-grade fluid model.In all of the above mentioned models, the so-called no-slip boundary condition is used; namely, the velocity of flow relative to the solid is zero on the fluid-solid interface [12]. Although the no-slip condition is supported by many experimental results, the existence of slip of a fluid on the solid surface was also observed by many other researches [13–20]. The Navier slip condition has been used by various researchers to describe boundary slip and is a more general boundary condition, in which the fluid velocity component tangential to the solid surface, relative to the solid surface, is proportional to the shear stress on the fluid-solid interface and the slip length. The surface characteristics constant, slip length, describes the “slipperiness” of the surface. Recently, we and many other researchers have investigated various flow problems of Newtonian fluids with the traditional no-slip and the Navier slip boundary conditions [12, 20–30], and it is found that the boundary slip and the slip parameter have significant influence on the flow of Newtonian fluids through microchannels and tubes.Motivated by the above mentioned work, we extend previous work on slip flows of Newtonian fluids [21, 22] to the case involving both Newtonian and non-Newtonian fluid flow in the flow region. The new feature and contribution of this work include establishment of the underlying boundary value problem for the problem, the derivation of an exact solution for a special case, and demonstration of the influence of the slip parameter on the flow profile and flow behavior. The rest of the paper is organized as follows. In Section 2, we present the underlying boundary value problem for the problem in dimensionless form. Then in Section 3, we derive an exact solution for a special case. In Section 4, we investigate numerically the effect of the slip parameter for the general case. Finally, a conclusion is given in Section 5.


## 2. Mathematical Formulation

The flow of a fluid with no thermochemical and electromagnetic effects can be described by the conservation equations of mass and linear momentum; namely,(1)∂ρ∂t+div(ρv)=0,ρ(∂v∂t+v·∇v)=divT+ρb,
where ρ is the density of the fluid, ∂/∂t is the partial derivative with respect to time, v is the velocity vector,b is the body force vector, and T is the stress tensor.The stress tensor is related to the velocity gradient by the constitutive equations. For a modified (generalized) second-grade fluid [11, 31, 32], the constitutive equations can be expressed by
(2)T=-pI+Πm/2(μA1+α1A2+α2A12,),
where m is a material parameter, Π=(1/2)trA12 is the second invariant of A1,p is the fluid pressure, μ is the coefficient of viscosity, αi are material moduli (the normal stress coefficients), and Ai are the kinematical tensors given by
(3)A1=L+LT,A2=∂A1∂t+[grad(A1)]v+A1L+(L)TA1,
in which L is gradv and the superscript T refers to matrix transposition.For the axially symmetrical blood flow through a circular tube of radiusb, we can assume that  v=v(r,t)ez, where z is the axial direction and r is the radial direction. Under the periodic body acceleration and a unsteady pulsatile pressure gradient [7, 10], the momentum equation in the z-direction in the cylindrical polar coordinate (r,θ,z) is
(4)ρ∂v∂t=-∂p∂z+ρG+1r∂∂r(rTrz).The shear stressTrz for a generalized second-grade fluid can be expressed by
(5)Trz={μ1|∂v1∂r|m∂v1∂r0≤r≤a,μ2∂v2∂ra≤r≤b.The approximate periodic form of the pressure gradient generated by the heart can be described by(6)-∂p∂z=A0+A1cosωpt,
where A0,A1,ωp=2πfp, and fp are the constant component of the pressure gradient, the amplitude of the pressure fluctuation (establishing the systolic and diastolic pressures), the circular frequency, and the frequency of pulse rate, respectively.The body accelerationG can be approximated by
(7)G=Agcos(ωbt+ϕ),
where Ag is the amplitude, fb=ωb/2π is the frequency, and ϕ is the lead angle of G with respect to the action of the heart.Substituting (5)–(7) into (4), the blood flow equation for a modified second-grade fluid in the z-direction, in the inner and outer core, becomes
(8)ρ1∂v1∂t=A0+A1cosωpt+ρAgcos(ωbt+ϕ)+1r∂∂r(rμ1|∂v1∂r|m∂v1∂r),for0≤r≤a,ρ2∂v2∂t=A0+A1cosωpt+ρAgcos(ωbt+ϕ)+1r∂∂r(rμ2∂v2∂r),fora≤r≤b.In order to completely define the problem, boundary and initial conditions are required. In this work, the Navier slip condition is applied. That is, on the solid-fluid interfacer=b, the axial fluid velocity, relative to the solid surface, is proportional to the shear stress on the interface. As the fluid layer near the wall is modeled as a Newtonian fluid in our model, the shear stress on the boundary is related to the shear strain rate by σrz=μ2(∂v/∂z). Thus, the Navier slip condition can be written as
(9)v2(b,t)+l∂v2∂t(b,t)=0,
where l is the slip parameter. Moreover, we assume that the slip parameter does not change along the axial direction.Onr=0, the symmetry condition is introduced:
(10)∂v1∂r(0,t)=0.On the interface between two different fluids, for continuous and smooth behavior of the velocity and shear stresses, we require(11)v1(a,t)=v2(a,t),[μ1|∂v1∂r|m∂v1∂r](a,t)=[μ2∂v2∂r](a,t).The initial conditions are set to(12)v1(r,0)=0=v2(r,0),
which is essential for the numerical scheme adopted to estimate the time at which the pulsatile steady state is achieved.To simplify the equations, we introduce the following nondimensional variables and parameters:(2)r¯=rb,v¯=vv0,t¯=ωp2πt,u0=A0b2μ2,e=A1A0,ωr=ωbωp,r0=ab,mu¯=μ(u0b)m,ρ*=ρ1ρ2,μ*=μ2μ¯,C1=A0b2μ¯u0,C2=ρ1Agb2μ¯u0=ρ1AgA0B1,α=ρ1ωpb22πμ¯,γ=ρ2ωpb22πμ¯μ*=ρ2ωpb2ρ12πμ¯μ*ρ1=αρ*μ*,C^1=A0b2μ¯u0μ*=C1μ*,C^2=ρ2Agb2ρ1μ¯u0μ*ρ1=C2ρ*μ*.In terms of the nondimensional variables and parameters, (8)–(12) can be written in the form of
(14)α∂v-1∂t¯=C1(1+ecos2πt¯)+C2cos(2πωrt¯+ϕ)+1r¯∂∂r¯[r¯|∂v-1∂r¯|m∂v-1∂r¯],for0≤r≤r0,γ∂v-2∂t¯=C-1(1+ecos2πt¯)+C-2cos(2πωrt¯+ϕ)+1r¯∂∂r¯[r¯∂v-2∂r¯],forr0≤r≤1.The boundary conditions and initial conditions, in dimensionless form, can be expressed by(15)∂v-1∂r¯(0,t¯)=0,(16)bv-2(1,t)+l∂v-2∂r¯(1,t¯)=0,(17)v-1(r0,t¯)=v-2(r0,t¯),(18)[|∂v-1∂r¯|m∂v-1∂r¯](r0,t¯)=[μ*∂v-2∂r¯](r0,t¯),(19)v-1(r,0)=0=v-2(r,0).


## 3. Analytical Solution

Form=0, the model reduces to the linear model with different viscosity in the peripheral layer and the centre core. In this case, (14) have the same form:
(20)L(v)=β∂v∂t-1r∂v∂r-∂2v∂r2=B1(1+ecos(2πt))+B2cos(2πωrt+ϕ).By the superposition principle, ifv0,v1, and v2 are the solution of L(v)=f(t), respectively, for f(t)=B1e0ti,B1ae2πti, and B2e(2πωrt+ϕ)i, then the complete solution of (20) is v=∑n=02Re(vn).To determinevn, we solve
(21)β∂vn∂t=Dnegn(t)i+1r∂vn∂r+∂2vn∂r2,
where g0(t)=0,g1(t)=2πt,g2(t)=2πωrt+ϕ,D0=B1,D1=aB1, and D2=B2. As (21) admits solutions of the formvn=fn(r)egn(t)i, we have from (21) that
(22)βgn′(t)fn(r)egn(t)ii=Dnegn(t)i+1rfn′(r)egn(t)i+fn′′(r)egn(t)i.
Dividing by egn(t)i on both sides of (22), we obtain
(23)βgn′(t)fn(r)i=Dn+1rfn′(r)+fn′′(r).Forn=0, we get
(24)f0′′(r)+1rf0′(r)=-B1,
which has the general solution: f0(r)=(c1+c2lnr)-(B1/4)r2.Forn=1, we have
(25)f1′′(r)+1rf1′(r)-2πβif1(r)=-eB1.Letβ-12=-2πβi; then,
(26)1β-12f1′′(r)+1β-12rf1′(r)+f1(r)=-eB1β-12.Letr^=β-1r; we have
(27)r^2f1′′(r^)+r^f1′(r^)+r^2f1(r^)=-eB1β-12r^2.The general solution of (27) is
(28)f1(r)=d1J0(β-1r)+e1Y0(β-1r)-eB12πβi,
where d1 and e1 are integration constants and J0 and Y0 denote the zero-order Bessel functions of the first kind and the second kind, respectively.Similarly, forn=2, we have
(29)f2′′(r)+1rf2′(r)-2βπωrf2(r)i=-B2,
and the general solution is
(30)f2=d2J0(β-2r)+e2Y0(β-2r)-B22βωrπi,
where β-22=-2πβωri.Because the boundness ofv1,v2, c2,e1, and e2 are set to zero, hence, from (14) and the solutions for (20), we have
(31)v-1=Re{c1-C14r¯2+[d1J0(β1r¯)-eC12παi]e2πt¯i+[d2J0(β2r¯)-C22πωrαi]e(2πωrt¯+ϕ)i[d1J0(β1r¯)-eC12παi]},v-2=Re{c^1+c^2lnr¯-C-14r¯2+[d^1J0(β^1r¯)+e^1Y0(β^1r¯)-eC^12πγi]e2πt¯i+[d^2J0(β^2r¯)+e^2Y0(β^2r¯)-c-2i2πωrγ]×e(2πωrt¯+ϕ)iC-14},
where β^12=-2πγi, β^22=-2πωrγi,β12=-2παi, and β22=-2πωrαi.AsdJ0(x)/dx=-J1(x) and dY0(x)/dx=-Y0(x), we have
(32)∂v-1∂r¯=Re(-C12r¯-d1β1J1(β1r¯)e2πt¯i-d2β2J1(β2r¯)e(2πωrt¯+ϕ)iC12r¯),∂v-2∂r¯=Re(c^21r¯-C-12r¯+[-d^1β^1J1(β^1r¯)-e^1β^1Y1(β^1r¯)]e2πt¯i+[-d^2β^2J1(β^2r¯)-e^2β^2Y1(β^2r¯)]e(2πωrt¯+ϕ)iC-121r¯).
Obviously, v1 satisfies the boundary condition (15) automatically. We now consider the boundary condition (16); namely,
(33)Re[(bc^1+lc^2-(l+b2)C-12)+(ebC-1i2πγ(bJ0(β^1)-lβ^1J1(β^1))d^1+(bY0(β^1)-lβ^1Y1(β^1))e^1-ebC-1i2πγ)e2πt¯i+(bC2¯2πωrγid^2(bJ0(β^2)-lβ^2J1(β^2))+e^2(bY0(β^2)-lβ^2Y1(β^2))-bC-22πωrγi)e(2πωrt¯+ϕ)i]=0.Further, from boundary conditions (17) and (18), we have
(34)Re[(c1-c^1-c^2lnr0-(C1-C^1)r024)+(ei2παγd1J0(β1r0)-d^1J0(β^1r0)-e^1Y0(β^1r0)-(γC1-αC-1)ei2παγ)e2πt¯i+(i2πωrγαd2J0(β2r0)-d^2J0(β^2r0)-e^2Y0(β^2r0)-(γC2-αC-2)i2πωrγα)×e(2πωrt¯+ϕ)i(c1-c^1-c^2lnr0-(C1-C^1)r024)]=0,Re[((μ*C-1-C1)r02-μ*c^2r0)+((β^1r0)-d1β1J1(β1r0)+d^1μ*β^1J1(β^1r0)+e^1μ*β^1Y1(β^1r0))e2πt¯i+(-d2β2J1(β2r0)+d^2μ*β^2J1(β^2r0)+e^2μ*β^2Y1(β^2r0))×e(2πωrt¯+ϕ)i((μ*C-1-C1)r02-μ*c^2r0)]=0.As (33)-(34) must be satisfied for any instant of time t, we require that the constant terms and the coefficients of the exponential terms all vanish; namely,
(35)bc^1+lc^2-(l+b2)C-12=0,c1-c^1-c^2lnr0-(C1-C-1)r024=0,(μ*C-1-C1)r02-μ*c^2r0=0,d^1(bJ0(β^1)-lβ^1J1(β^1))+e^1(bY0(β^1)-lβ^1Y1(β^1))-ebC-12πγi=0,d1J0(β1r0)-d^1J0(β^1r0)-e^1Y0(β^1r0)-eC12παi+eC-12πγi=0,-d1β1J1(β1r0)+d^1μ*β^1J1(β^1r0)+e^1μ*β^1Y1(β^1r0)=0,d^2(bJ0(β^2)-lβ^2J1(β^2))+e^2(bY0(β^2)-lβ^2Y1(β^2))-bC-22πωrγi=0,d2J0(β2r0)-d^2J0(β^2r0)-e^2Y0(β^2r0)-C12πωrαi+C-22πωrγi=0,-d2β2J1(β2r0)+d^2μ*β^2J1(β^2r0)+e^2μ*β^2Y1(β^2r0)=0.Solving the above system of equations yields(36)c1=(lnr0-lb)((μ*C-1-C1)r022μ*)+(lb+1-r022)C-12+C1r024,c^1=-lb((μ*C-1-C1)r022μ*)+(lb+12)C-12,c^2=(μ*C-1-C1)r022μ*,d1=μ*[β^1ei2πγα(J1(β^1r0)Y0(β^1r0)-J0(β^1r0)Y1(β^1r0))×ebβ^1C-1i2πγ+(γC1-αC-1)×[J1(β^1r0)(bY0(β^1)-lβ^1Y1(β^1))-Y1(β^1r0)(bJ0(β^1)-lβ^1J1(β^1))]×β^1ei2πγα]/(bJ0(β^1)-lβ^1J1(β^1))×(β1J1(β1r0)Y0(β^1r0)-μ*β^1Y1(β^1r0)J0(β1r0))+(bY0(β^1)-lβ^1Y1(β^1))×(μ*β^1J0(β1r0)J1(β^1r0)-β1J1(β1r0)J0(β^1r0)),d^1=[(β1J1(β1r0)Y0(β^1r0)-μ*β^1Y1(β^1r0)J0(β1r0))×ebC-12πγi+(γC1-αC-1)2πγα×(bY0(β^1)-lβ^1Y1(β^1))×eβ1J1(β1r0)i(β1J1(β1r0)Y0(β^1r0)-μ*β^1Y1(β^1r0)J0(β1r0))]/(bJ0(β^1)-lβ^1J1(β^1))×(β1J1(β1r0)Y0(β^1r0)-μ*β^1Y1(β^1r0)J0(β1r0))+(bY0(β^1)-lβ^1Y1(β^1))×(μ*β^1J0(β1r0)J1(β^1r0)-β1J1(β1r0)J0(β^1r0)),e^1=[β1ei2παγ(μ*β^1J0(β1r0)J1(β^1r0)-β1J0(β^1r0)J1(β1r0))×ebC-1i2πγ-(bJ0(β^1)-lβ^1J1(β^1))×(γC1-αC-1)×J1(β1r0)β1ei2παγ]/(bJ0(β^1)-lβ^1J1(β^1))×(β1J1(β1r0)Y0(β^1r0)-μ*β^1Y1(β^1r0)J0(β1r0))+(bY0(β^1)-lβ^1Y1(β^1))×(μ*β^1J0(β1r0)J1(β^1r0)-β1J1(β1r0)J0(β^1r0)),d2=μ*[β^2i2πωrγα(J1(β^2r0)Y0(β^2r0)-J0(β^2r0)Y1(β^2r0))×bβ^2C-1i2πωrγ+(γC2-αC-2)×[J1(β^2r0)(bY0(β^2)-lβ^2Y1(β^2))-Y1(β^2r0)(bJ0(β^2)-lβ^2J1(β^2))]×β^2i2πωrγα]/(bJ0(β^2)-lβ^2J1(β^2))×((β^2r0)β2J1(β2r0)Y0(β^2r0)-μ*β^2Y1(β^2r0)J0(β2r0))+(bY0(β^2)-lβ^2Y1(β^2))×(μ*β^2J0(β2r0)J1(β^2r0)-β2J1(β2r0)J0(β^2r0)),d^2=[(β2J1(β2r0)Y0(β^2r0)-μ*β^2Y1(β^2r0)J0(β2r0))×bC-12πωrγi+(γC1-αC-1)2πγαωr×(bY0(β^2)-lβ^2Y1(β^2))×β2J1(β2r0)i(β2J1(β2r0)Y0(β^2r0)-μ*β^2Y1(β^2r0)J0(β2r0))]/(bJ0(β^2)-lβ^2J1(β^2))×(β2J1(β2r0)Y0(β^2r0)-μ*β^2Y1(β^2r0)J0(β2r0))+(bY0(β^2)-lβ^2Y1(β^2))×(μ*β^2J0(β2r0)J1(β^2r0)-β2J1(β2r0)J0(β^2r0)),e^2=[(μ*β^2J0(β2r0)J1(β^2r0)-β2J0(β^2r0)J1(β2r0))×bC-1i2πωrγ-(bJ0(β^2)-lβ^2J1(β^2))×(γC1-αC-1)J1(β2r0)×β2i2πωrαγ]/(bJ0(β^2)-lβ^2J1(β^2))×((β^2r0)β2J1(β2r0)Y0(β^2r0)-μ*β^2Y1(β^2r0)J0(β2r0))+(bY0(β^2)-lβ^2Y1(β^2))×(μ*β^2J0(β2r0)J1(β^2r0)-β2J1(β2r0)J0(β^2r0)).To show the flow behavior and the effect of the slip parameter, we investigate the velocity profiles in the arteries with different values of the slip parameter under various different conditions. In the first example of investigation, the radius of the artery is taken asr=b=0.15 cm, and the other parameters are set to A0=698.65 dyne/cm3, Ag=0.5g, fb=fp=1.2,ϕ=0,C1=6.6,C2=4.64,A1=1.2A0, and ρ1/ρ2=1. Figure 1 shows the 3-dimensional velocity profile as a function of time and location and the 2-dimensional velocity profile as a function of time at three different radial locations for two different slip parameters l=0 (no-slip) and l=2. The results show that boundary slip has a very dramatical effect on the fluid flow in the artery. It affects not only the magnitude of the flow velocity significantly, but also the flow pattern and velocity profile on the cross-section of the artery. For the no-slip flow (lb=0), the pulsatile flow nature gradually disappears toward the arterial wall, while with boundary slip, the flow near the arterial wall also displays a pulsatile nature.The velocity profile in the small artery with radius 0.15 cm under two different slip parameter values: (a)lb=0; (b) lb=2. In the figure, the 3D graphs show the variation of velocity as a function of time and location, while the 2D graphs show the variation of velocity with time at three radial locations including the artery centre (r=0), the interface of inner-outer layer (r=0.6), and the arterial wall (r=1).
(a)
(b)We then investigate whether the above observed flow phenomena associated with boundary slip are affected or not by the radius of the artery, and for this purpose, we consider the fluid flow through an artery with a larger radiusr=0.5 cm. The constant pressure gradient is set to A0=32 dyne/cm3 in order to achieve a mean velocity magnitude approximately equal to that in the smaller artery, while all other parameters are set to the same values as those used for the smaller radius. Figure 2 shows the velocity profile in the artery for two different slip parameter values including lb=0 (no-slip) and lb=2. The 3-dimensional graph shows the variation of the flow velocity with time and radial position, while the 2-dimensional graphs demonstrate the variation of the flow velocity with time at three different radial locations including r=0 (centre), r=0.6 (inner-outer layers interface), and r=1 (arterial wall). From Figures 1 and 2, it is clear that the boundary slip related flow phenomena and behavior observed for the smaller artery also appear in the artery with a larger radius, and further, a more significant pulsatile nature of fluid flow is observed for the larger artery.The velocity profile in the large artery with radius 0.50 cm under two different slip parameter values: (a)lb=0; (b) lb=2. In the figure, the 3D graphs show the variation of velocity as a function of time and location, while the 2D graphs show the variation of velocity with time at three radial locations including the artery centre (r=0), the interface of inner-outer layer (r=0.6), and the arterial wall (r=1).
(a)
(b)To further investigate the effect of the slip parameter on the velocity profile near the artery wall, we show in Figure3 the velocity of fluid on the artery wall for four different values of the slip parameter including lb=0, 2, 4, 6, and 8. The results clearly demonstrate that the slip parameter has a very significant effect on the near-wall velocity and that the magnitude of the average wall velocity is proportional to the slip parameter.Diagrams showing the velocity profile on the arterial wall with five different slip parameterslb for two different artery radii (a) r=0.15 cm; (b) r=0.5 cm.
(a)
(b)


## 4. Numerical Investigation

A numerical scheme, based on the finite different method, is established to solve the underlying boundary value problem for the general casem≠0, consisting of (14) and boundary condition (15)–(19). To validate the numerical technique, we apply the numerical scheme to generate a series of numerical solutions for the case m=0 and then compare the numerical results with the exact solution derived in Section 3.Figure4 presents the velocity profile in the small and large arteries for two different slip parameters lb=0 (no-slip) and lb=2 obtained by the numerical technique. The numerical errors between the exact solution and the numerical solution, Er=V-U, are presented in Figure 5 in which V is the exact solution and U is the numerical solution. The results clearly indicate that the numerical solution converges to the exact solution. This shows that a larger slip length has a lower convergence rate.Velocity profiles in arteries with different radiir: (a) r=0.15 cm; (b) r=0.5 cm. In the figure, the graphs on the left column correspond to lb=0, while the graphs on the right column correspond to lb=2.
(a)Diagrams showing the convergence of numerical solutions for different slip parameters and artery radii: (a)r=0.15 cm; (b) r=0.50 cm.
(a)
(b)We then investigate the flow phenomena for the general casem≠0, and here we consider m=-1/4 in the investigation. Figure 6 gives the 3D graph showing the convergence of the transient velocity field to a steady state pulsatile velocity field and also demonstrating the substantial influence of boundary slip on the steady state velocity profile in both magnitude and flow pattern. Figure 7 shows the variations of velocities with time at three arterial locations for different slip parameters and artery radii and also clearly demonstrates the significant effect of boundary slip on the flow through the artery. Figure 8 shows the variation of fluid velocity along the artery wall under different slip parameters and artery radii. The results show that as the slip parameter increases, the time required for achieving convergence results increases, and the magnitude of the average steady state velocity also increases.Velocity profiles in arteries with different slip parameterslb and radii r: (a) r=0.15 cm; (b) r=0.50 cm. In the Figure, the graphs on the left column correspond to lb=0, while the graphs on the right column correspond to lb=2.
(a)
(b)Velocity profiles at three arterial locations(r1,r2,r3): for m=-1/4 and under different slip parameters lb and artery radii (a) r=0.15 cm; (b) r=0.50 cm. In the Figure, the graphs on the left column correspond to lb=0, while the graphs on the right column correspond to lb=2.
(a)
(b)Diagrams showing the convergence of numerical results of the fluid velocity on the wall to the steady state pulsatile velocity field under various slip parameterslb for two different artery radii: (a) r=0.15 cm; (b) r=0.50 cm.
(a)
(b)


## 5. Conclusion

In this paper, a mathematical model for the transient pulsatile flow of fluids through vessels, taking into account boundary slip and the Fahraeus-Lindqvist effect, is established. For a special case of the underlying boundary value problem, an exact solution for the velocity field has been derived in explicit form, which provides one with an exact analytical method for investigating the flow phenomena under the special case and also a mean for validating the subsequently developed numerical scheme for generating numerical results for the general case. Our analytical and numerical studies show that for the flow of fluids with the Fahraeus-Lindqvist effect, boundary slip has a very significant influence on the magnitude of the mean flow velocity and on the flow pattern and velocity profile on the cross-section. With boundary slip, the boundary layer near the wall also displays significant pulsatile flow nature. The results also show that as the boundary slip length increases, the convergence rate of numerical results to the exact solutions decreases and the time required to achieve the steady state pulsatile flow increases.



---

*Source: 858597-2013-09-25.xml*