^{1}

^{2}

^{1}

^{1}

^{*}

This paper provides the exact solutions for the fully developed two layer pressure driven flows of incompressible Phan-Thien-Tanner fluids in a horizontal cylindrical pipe. Exact equations are formulated and solved for important kinematic properties, such as, velocity profiles, normal and shear stresses, total volume fluxes through a circular cross-section and average velocities. Graphical results are provided and discussed for the different flow parameters. A comparison of Upper Convected Maxwell (UCM), Linear Phan-Thien-Tanner (LPTT) and Exponential Phan-Thien-Tanner (EPTT) shows that UCM is a low viscosity fluid as compared to LPTT, and EPTT and LPTT is lighter than EPTT.

The Navier-Stokes equations governing the motion of the viscous fluid are non-linear and their exact solutions are rare or non-existent. The nonlinearity in these equations is caused by the presence of convected term and in the case of non-Newtonian fluids, the use of constitutive equations. Because of the great variety in physical structure of non-Newtonian fluids, it is not possible to describe their mechanical behavior by a single constitutive equation. For this reason, a great variety of constitutive equations have been proposed [

Recently, the interest in stratified laminar flow of two or more immiscible fluids has grown considerably because of their wide use in technological processes. There has been some theoretical and experimental work on the stratified laminar flow of two immiscible fluids in a horizontal pipe. For instance, Packham and Shall [

In this paper mathematics reveals a need for a system of four first-order linear ordinary differential equations associated with four inhomogeneous boundary conditions. We will solve this system of differential equations simultaneously, including the boundary conditions, in order to calculate the velocity profiles, total volume fluxes and the average velocities resulting from the flow UCM, LPTT and EPTT models.

The basic equations governing the isothermal two-layer flow of immiscible incompressible fluids are the continuity and momentum equations. Mathematically,

d i v V ( k ) = 0 ; k = 1 , 2 , (1)

ρ ( k ) D V ( k ) D t = − ∇ p + ∇ ⋅ T ( k ) + ρ ( k ) b ; k = 1 , 2 , (2)

where the superscript k denotes the two fluids, the number k = 1 represents the fluid in the core, while the fluid properties for the outer fluid along the wall are denoted by the superscript k = 2 . V ( k ) , ρ ( k ) , T ( k ) are the velocity, density and extra stress tensor of the k fluid respectively. D D t is the material derivative defined as

D ( ⋆ ) D t = ∂ ( ⋆ ) ∂ t + ( V ⋅ ∇ ) ( ⋆ ) .

The general form of the constitutive equation defining the PTT fluid is

f ( t r ( T ( k ) ) ) T ( k ) + λ ( k ) T ∇ ( k ) = η ( k ) A 1 ( k ) ; k = 1 , 2 , (3)

where η ( k ) , λ ( k ) and t r ( T ( k ) ) are the viscosity, relaxation time and trace of stress tensor of the k fluid and Rivlin-Ericksen tensor A 1 ( k ) given by

A 1 ( k ) = ∇ V ( k ) + ( ∇ V ( k ) ) T ; k = 1 , 2 , (4)

and upper convected derivative T ∇ ( k ) is defined as

T ∇ ( k ) = D T ( k ) D t − ( T ( k ) ( ∇ V ( k ) ) + ( ∇ V ( k ) ) T T ( k ) ) ; k = 1 , 2. (5)

The three PTT models, commonly used are

f ( t r T ( k ) ) = 1, ( UCMModel ) , (6)

f ( t r T ( k ) ) = 1 + ϵ ( k ) λ ( k ) η ( k ) t r ( T ( k ) ) , ( LPTTModel ) , (7)

f ( t r T ( k ) ) = e x p ( ϵ ( k ) λ ( k ) η ( k ) t r ( T ( k ) ) ) , ( EPTTModel ) , (8)

where ϵ ( k ) is a parameter related to the elongation behavior of the model. In the absence of the parameters ϵ ( k ) , the models (7) and (8) reduce to the well known Maxwell model (6).

We consider two-layer immiscible flow of incompressible PTT fluids in a horizontal cylindrical pipe of radius R_{2} as shown in _{1} from the axis of the pipe. To visualize the problem a cylindrical coordinate system is chosen, where the axis of the pipe will be denoted by z, while r will denote the radially outward measured distance from the z-axis. The flow is driven by the constant pressure gradient, assuming that the flow is steady, unidirectional and fully developed.

The velocities and extra stresses can be taken of the form

V ( k ) = [ 0 , 0 , w ( k ) ( r ) ] , T ( k ) = T ( k ) ( r ) ; k = 1 , 2 , (9)

where w is the axial velocity component. Using Equation (9) in Equations (3)-(5), non-zero stress components become

f ( t r T ( k ) ) T r z ( k ) = η ( k ) d w ( k ) d r ; k = 1,2, (10)

f ( t r T ( k ) ) T r r ( k ) = 2 λ ( k ) T r z ( k ) d w ( k ) d r ; k = 1 , 2. (11)

Solving Equation (10) and Equation (11) simultaneously, we obtain

T r r ( k ) = 2 λ ( k ) T r z ( k ) 2 η ( k ) ; k = 1 , 2. (12)

Trace of stress tensor of kth PTT fluid is given by

t r ( T ( k ) ) = T r r ( k ) + T θ θ ( k ) + T z z ( k ) = 2 λ ( k ) T r z ( k ) 2 η ( k ) ; k = 1 , 2. (13)

Using Equation (9) in Equation (1), the continuity equation is identically satisfied and the non-zero component of the momentum Equation (2) becomes

1 r ∂ ∂ r ( r T r z ( k ) ) = − C ; k = 1 , 2 , (14)

where ( − ∂ p ∂ z ) = C . The boundary conditions for the problem under consideration are given by

T r z ( 1 ) = 0 , at r = 0 , (15)

T r z ( 1 ) = T r z ( 2 ) , at r = R 1 , (16)

w ( 1 ) = w ( 2 ) , at r = R 1 , (17)

w ( 2 ) = 0 , at r = R 2 . (18)

We note here that these conditions are four in number and are sufficient to find the solution to our problem which consists of four first order differential equations, two first order differential equations represented by Equation (14) and two differential equations obtained for velocity by simplifying either Equation (10) or Equation (11). Hence, our problem is consistent and we attempt to find exact solutions to the problem.

Equation (14), upon integration with respect to r, yields

r T r z ( k ) = − C r 2 2 + A ( k ) , k = 1 , 2 (19)

where A ( k ) , k = 1 , 2 , are constants of integration. Using the condition (15) in the equation (19) when k = 1 , we find that A ( 1 ) = 0 , so that we obtain

T r z ( 1 ) = − C r 2 ; 0 ≤ r ≤ R 1 . (20)

Similarly, Equation (19) with the help of boundary condition (16) gives A ( 2 ) = 0 which yields

T r z ( 2 ) = − C r 2 ; R 1 ≤ r ≤ R 2 . (21)

The expressions (20) and (21) clearly indicate that the shear stresses of both the fluids are the same and vary linearly with the distance from the axis of the pipe, which implies

T r z ( k ) = − C r 2 ; k = 1 , 2 , (22)

for which we see that for fluid (1) and fluid (2) the regions are 0 ≤ r ≤ R 1 and R 1 ≤ r ≤ R 2 , respectively.

This last result when used in Equation (12) for normal stresses leads to

T r r ( k ) = λ ( k ) C 2 2 η ( k ) r 2 ; k = 1 , 2. (23)

Thus the normal stresses for both the fluids are different, depending on their corresponding material constants. Moreover, these normal stresses increase with the square of the distance from the axis of the tube. For Newtonian fluids, λ ( k ) = 0 and we find that the shear stresses remain unchanged while the normal stresses are zero.

Trace of extra stress tensor is given by

t r ( T ( k ) ) = λ ( k ) C 2 2 η ( k ) r 2 ; k = 1 , 2. (24)

Velocity profile can be calculated by rewriting in simplified form equation (10) resulting in the equation

d w ( k ) d r = f ( t r T ( k ) ) T r z ( k ) η ( k ) ; k = 1 , 2 (25)

which is a system of first order non-homogeneous linear ordinary differential equations for two immiscible PTT fluids.

By making use of expressions (6) and (22) in expression (25), we get

d w u ( k ) d r = − C 2 η ( k ) r , k = 1 , 2 , (26)

where subscript u stands for UCM model. By integrating Equation (26), we get velocity profiles for fluid (1) and fluid (2), with two constants of integration and those constants are evaluated by using boundary conditions (18) and (17). Velocity profiles for both immiscible fluid are given by

w u ( 1 ) = C ( R 1 2 − r 2 4 η ( 1 ) + R 2 2 − R 1 2 4 η ( 2 ) ) ; 0 ≤ r ≤ R 1 , w u ( 2 ) = C ( R 2 2 − r 2 4 η ( 2 ) ) ; R 1 ≤ r ≤ R 2 . (27)

Flow rates through a circular cross section of a pipe bounded in the region 0 ≤ r ≤ R 1 and R 1 ≤ r ≤ R 2 are given by

Q u ( 1 ) = 2 π ∫ 0 R 1 w u ( 1 ) ( r ) r d r , Q u ( 2 ) = 2 π ∫ R 1 R 2 w u ( 2 ) ( r ) r d r . (28)

We use the expression for velocity profiles (27) in formulas (28) and obtain

Q u ( 1 ) = π R 1 2 C ( R 2 2 − R 1 2 4 η ( 2 ) + R 1 2 8 η ( 1 ) ) , (29)

Q u ( 2 ) = π C ( R 2 4 + R 1 4 − 2 R 1 2 R 2 2 8 η ( 2 ) ) . (30)

Net volume flux through the circular pipe of radius R 2 is given by

Q u = π C 8 ( R 2 4 − R 1 4 η ( 2 ) + R 1 4 η ( 1 ) ) . (31)

Average velocity is given by

v ¯ u = Q u π R 2 2 . (32)

Substituting expression (31), average velocity becomes

v ¯ u = C 8 R 2 2 ( R 2 4 − R 1 4 η ( 2 ) + R 1 4 η ( 1 ) ) , (33)

Above results are same as for the case of two immiscible Newtonian fluids with viscosities η ( 1 ) and η ( 2 ) .

By using trace of extra stress tensor (24) in linear model (7) and then using that model and shear stresses (22), in differential Equation (25), we get

d w l ( k ) d r = − C 2 η ( k ) r − C 3 ϵ ( k ) λ ( k ) 2 4 η ( k ) 3 r 3 , k = 1 , 2 , (34)

where subscript l stands for LPTT fluid model. Integrating with respect to r we get

w l ( k ) = − C 4 η ( k ) r 2 − C 3 ϵ ( k ) λ ( k ) 2 16 η ( k ) 3 r 4 + B k , k = 1 , 2 , (35)

where B k , k = 1 , 2 are the constants of integration. For k = 2 , the Equation (34) with the boundary condition (18), yields

w l ( 2 ) = C 4 η ( 2 ) ( R 2 2 − r 2 ) + C 3 ϵ ( 2 ) λ ( 2 ) 2 16 η ( 2 ) 3 ( R 2 4 − r 4 ) ; R 1 ≤ r ≤ R 2 . (36)

with the help of the expression (36) for w l ( 2 ) and the boundary condition (17) the Equation (35) for k = 1 , gives

w l ( 1 ) = C 4 [ ( R 1 2 − r 2 ) η ( 1 ) + ( R 2 2 − R 1 2 ) η ( 2 ) ] + C 3 16 [ ϵ ( 1 ) λ ( 1 ) 2 η ( 1 ) 3 ( R 1 4 − r 4 ) + ϵ ( 2 ) λ ( 2 ) 2 η ( 2 ) 3 ( R 2 4 − R 1 4 ) ] ; 0 ≤ r ≤ R 1 . (37)

Flow rates for both the fluids are obtained by using Equation (36) and Equation (37) in Equation (28) as

Q l ( 1 ) = π C R 1 2 { R 2 2 − R 1 2 4 η ( 2 ) + R 1 2 8 η ( 1 ) } + π C 3 R 1 2 8 { ϵ ( 2 ) λ ( 2 ) 2 2 ( η ( 2 ) ) 3 ( R 2 4 − R 1 4 ) + ϵ ( 2 ) λ ( 1 ) 2 3 ( η ( 1 ) ) 3 R 1 4 } , (38)

Q l ( 2 ) = π C { R 1 4 + R 2 4 − 2 R 2 2 R 1 2 8 η 2 } + π C 3 ϵ ( 2 ) λ ( 2 ) 2 48 η ( 2 ) 3 { 2 R 2 6 + R 1 6 − 3 R 2 4 R 1 2 } . (39)

The total volume flux Q l through the circular cross section of the pipe of radius R 2 for two immiscible LPTT fluids is obtained as

Q l = Q l ( 1 ) + Q l ( 2 ) = π C 8 [ R 2 4 − R 1 4 η ( 2 ) + R 1 4 η ( 1 ) + C 2 ( ϵ ( 2 ) ( λ ( 2 ) ) 2 3 ( η ( 2 ) ) 3 ( R 2 6 − R 1 6 ) + ϵ ( 1 ) ( λ ( 1 ) ) 2 3 ( η ( 1 ) ) 3 R 1 6 ) ] , (40)

and

v ¯ l = Q l π R 2 2 = C 8 R 2 2 [ R 2 4 − R 1 4 η ( 2 ) + R 1 4 η ( 1 ) + C 2 ( ϵ ( 2 ) ( λ ( 2 ) ) 2 3 ( η ( 2 ) ) 3 ( R 2 6 − R 1 6 ) + ϵ ( 1 ) ( λ ( 1 ) ) 2 3 ( η ( 1 ) ) 3 R 1 6 ) ] , (41)

is the average velocity of two immiscible LPTT fluids.

Constitutive equation for EPTT fluid model (8) is

f ( t r T ( k ) ) = exp ( ϵ ( k ) λ ( k ) η ( k ) t r T ( k ) ) , k = 1 , 2 , (42)

with the help of expression (24), we get

f ( t r T ( k ) ) = exp ( C 2 ϵ ( k ) λ ( k ) 2 2 η ( k ) 2 r 2 ) , k = 1 , 2. (43)

Let

M ( k ) = ϵ ( k ) λ ( k ) 2 2 η ( k ) 2 k = 1 , 2 , (44)

then

f ( t r T ( k ) ) = exp ( C 2 M ( k ) r 2 ) , k = 1 , 2. (45)

Substituting function (42) and shear stresses (22) in differential Equation (25), we get

d w e ( k ) d r = − exp ( C 2 M ( k ) r 2 ) C r 2 η ( k ) , k = 1 , 2 , (46)

where subscript e is used for EPTT fluid model. Integration with respect to r gives

w e ( k ) = − 1 4 C M ( k ) η ( k ) exp ( C 2 M ( k ) r 2 ) + E ( k ) , k = 1 , 2 , (47)

where E ( k ) , k = 1 , 2 are the constants of integration, and can be evaluated by using boundary conditions (18) and (17) for k = 1 , 2 . Both the velocity profiles become

w e ( 1 ) = 1 4 C M ( 1 ) η ( 1 ) [ exp ( C 2 M ( 1 ) R 1 2 ) − exp ( C 2 M ( 1 ) r 2 ) ] + 1 4 C M ( 2 ) η ( 2 ) [ exp ( C 2 M ( 2 ) R 2 2 ) − exp ( C 2 M ( 2 ) R 1 2 ) ] ; 0 ≤ r ≤ R 1 (48)

w e ( 2 ) = 1 4 C M ( 2 ) η ( 2 ) [ exp ( C 2 M ( 2 ) R 2 2 ) − exp ( C 2 M ( 2 ) r 2 ) ] ; R 1 ≤ r ≤ R 2 . (49)

Flow rates of both the immiscible fluids are obtained as

Q e ( 1 ) = π 4 C M ( 1 ) η ( 1 ) [ R 1 2 exp ( C 2 M ( 1 ) R 1 2 ) − 1 C 2 M ( 1 ) { exp ( C 2 M ( 1 ) R 1 2 ) − 1 } ] + π R 1 2 4 C M ( 2 ) η ( 2 ) [ exp ( C 2 M ( 2 ) R 2 2 ) − exp ( C 2 M ( 2 ) R 1 2 ) ] , (50)

Q e ( 2 ) = π 4 C M ( 2 ) η ( 2 ) [ − 1 C 2 M ( 2 ) { exp ( C 2 M ( 2 ) R 2 2 ) − exp ( C 2 M ( 2 ) R 1 2 ) } + ( R 2 2 − R 1 2 ) exp ( C 2 M ( 2 ) R 2 2 ) ] . (51)

Total volume flux is given by

Q e = π 4 C [ 1 M ( 1 ) η ( 1 ) { 1 C 2 M ( 1 ) + ( R 1 2 − 1 C 2 M ( 1 ) ) exp ( C 2 M ( 1 ) R 1 2 ) } + 1 M ( 2 ) η ( 2 ) { ( R 2 2 − 1 C 2 M ( 2 ) ) exp ( C 2 M ( 2 ) R 2 2 ) − ( R 1 2 − 1 C 2 M ( 2 ) ) exp ( C 2 M ( 2 ) R 1 2 ) } ] . (52)

Average velocity v ¯ is given by

v ¯ = 1 4 C R 2 2 [ 1 M ( 1 ) η ( 1 ) { 1 C 2 M ( 1 ) + ( R 1 2 − 1 C 2 M ( 1 ) ) exp ( C 2 M ( 1 ) R 1 2 ) } + 1 M ( 2 ) η ( 2 ) { ( R 2 2 − 1 C 2 M ( 2 ) ) exp ( C 2 M ( 2 ) R 2 2 ) − ( R 1 2 − 1 C 2 M ( 2 ) ) exp ( C 2 M ( 2 ) R 1 2 ) } ] , (53)

for two immiscible EPTT fluids.

Case #1 If R 1 = 0 , R 2 = R , η ( k ) = η , η ( k ) = η , λ ( k ) = λ and M ( k ) = M for k = 1 , 2 , then the pressure driven two fluids reduce to a single fluid and flow becomes Poiseuile flow for PTT fluids, that is,

w u = C ( R 2 − r 2 ) 4 η ; 0 ≤ r ≤ R . (54)

w l = C ( R 2 − r 2 ) 4 η [ 1 + M 2 C 2 ( R 2 + r 2 ) ] ; 0 ≤ r ≤ R . (55)

w e = 1 4 C M η [ exp ( M C 2 R 2 ) − exp ( M C 2 r 2 ) ] ; 0 ≤ r ≤ R . (56)

where R represents the radius of the cylinder, η is the viscosity of the fluid, ϵ the parameter of elongation, λ the relaxation time, and M = ϵ λ 2 2 η 2 .

Here, w u , w l and w e given by Equations (54)-(56), are the velocity profiles for pressure driven flows of UC (same as Newtonian), LPTT and EPTT fluid models, respectively.

Case #2 If R 2 → R 1 then we obtain solutions for one layer UC, LPTT, and EPTT fluid models, same as given by Case # 1.

Case #3 If λ ( 1 ) = 0 then two immiscible fluids will reduce to the cases of Newtonian-Newtonian, Newtonian-LPTT, Newtonian-EPTT immiscible fluid models.

Case #4 If λ ( 2 ) = 0 then we will obtain the cases for Newtonian-Newtownian, LPTT-Newtonian, EPTT-Newtonian models, respectively.

Case #5 If ϵ ( 1 ) = 0 then all of the three cases will reduce to Maxwell-UC, Maxwell-LPTT, Maxwell-EPTT fluid models, respectively.

Case #6 If ϵ ( 2 ) = 0 then the three cases will reduce to UC-Maxwell, LPTT-Maxwell, EPTT-Maxwell fluid models, respectively.

Introducing non-dimensional parameters

r * = r R 2 , z * = z R 2 , δ = R 1 R 2 , w ( k ) * = w ( k ) V ,

λ ( k ) * = λ ( k ) R k / V , ℜ k = ρ V R k η ( k ) , C * = C ρ V 2 / R 2 , k = 1 , 2 , (57)

with the assumption ρ ( 1 ) = ρ ( 2 ) = ρ , Equation (27) take the form

w u ( 1 ) * = C * 4 [ ℜ 1 δ ( δ 2 − r * 2 ) + ℜ 2 ( 1 − δ 2 ) ] ; 0 ≤ r ≤ δ , w u ( 2 ) * = C * 4 ℜ 2 ( 1 − r * 2 ) ; δ ≤ r ≤ 1. (58)

Equation (37) and Equation (36) are transformed as

w l ( 1 ) * = C * 4 [ ℜ 1 δ ( δ 2 − r * 2 ) + ℜ 2 ( 1 − δ 2 ) ] + C * 3 16 [ ℜ 1 3 δ M ( 1 ) * ( δ 4 − r * 4 ) + ℜ 2 3 M ( 2 ) * ( 1 − δ 4 ) ] ; 0 ≤ r ≤ δ , (59)

w l ( 2 ) * = C * 4 ℜ 2 ( 1 − r * 2 ) [ 1 + C * 2 4 ℜ 2 2 M ( 2 ) * ( 1 + r * 2 ) ] ; δ ≤ r ≤ 1. (60)

Velocity components for EPTT fluid model becomes

w e ( 1 ) * = 1 2 C * M ( 1 ) * 1 δ ℜ 1 [ exp ( δ 2 ℜ 1 2 2 C * 2 M ( 1 ) * ) − exp ( ℜ 1 2 2 C * 2 M ( 1 ) * r * 2 ) ] + 1 2 C * M ( 2 ) * 1 ℜ 2 [ exp ( ℜ 2 2 2 C * 2 M ( 2 ) * ) − exp ( δ 2 ℜ 2 2 2 C * 2 M ( 2 ) * ) ] ; 0 ≤ r * ≤ δ , (61)

w e ( 2 ) * = 1 2 C * M ( 2 ) * 1 ℜ 2 [ exp ( ℜ 2 2 2 C * 2 M ( 2 ) * ) − exp ( ℜ 2 2 2 C * 2 M ( 2 ) * r * 2 ) ] ; δ ≤ r * ≤ 1, (62)

where M ( k ) * = M ( k ) R k 2 / η ( k ) 2 V 2 .

In this paper, the exact solution of the two layer non-Newtonian fluid through a horizontally placed cylindrical pipe is obtained and derived mathematical results are discussed and graphed.

In this section, graphs demonstrating the behavior of the various parameters in conjunction with the velocity profile are presented.

In

Observing

fluid model are showing thinning effect (

To summarize the results of our analysis of the pressure driven flow of two immiscible non-Newtonian fluids in a pipe by focusing our attention on three fluid models known as UCM, LPTT and EPTT, we present visual evidence reflecting the behavior of the various parameters in conjunction with the velocity profile.

・ For all the three fluid models, the viscosity of fluid k = (1), in the region 0 ≤ r < δ does not effect the flow of the fluid k = (2), in the region δ < r ≤ 1 . Thus, the flow parameters of fluid (2) are independent of the viscosity of fluid (1).

・ Viscosity of fluid (1) is affected by the flow parameters (velocities and flow rates) of fluid (1) and fluid (2).

・ Increase in pressure gradient, elongation and relaxation time parameters enhance the speed of flow, showing the thinning effect.

・ For fixed value of parameters we compare the velocities as follows

w u > w l > w e

where w u , w l and w e are the velocities of the UCM, LPTT and EPTT fluid models, respectively.

・ We recover the Newtonian case from LPTT and EPTT fluid models as λ → 0 .

The authors greatly acknowledge the financial support by the Penn State York Campus to publish this work.

The authors declare no conflicts of interest regarding the publication of this paper.

Siddiqui, A.M., Walait, A., Allison, T. and Haroon, T. (2018) Exact Solution for Pressure Driven Flow of Two Immiscible Phan-Thien-Tanner Fluids in a Pipe. Open Journal of Fluid Dynamics, 8, 378-391. https://doi.org/10.4236/ojfd.2018.84024