## Numerical Solution of Power-law Fluid Flow through Eccentric Annular Geometry

**Nuha Hussein Ebrahim**^{1,}, **Noaman El-Khatib**^{1}, **Mariyamni Awang**^{1}

^{1}Petroluem Engineering, Universiti Teknologi Petronas, Perak Darul Ridzuan, Malaysia

### Abstract

Cuttings transport modeling in inclined and horizontal wellbores is complicated due to the eccentricity of the annulus. Development of a model for cuttings transport requires a deep understanding of the drilling mud flow behavior in the eccentric annular geometry. In this paper, bipolar coordinates system is used to solve for the eccentric annular geometry due to irregular shape of the boundaries. Finite difference method is used to obtain the velocity profile of Power-law non-Newtonian fluids through eccentric annuli. The discretized dimensionless Equation of flow using the finite difference method is solved iteratively using Point Successive Over Relaxation (S.O.R.) method. The results for Newtonian eccentric annular flow with 0.0001 eccentricity are in agreement with the Newtonian concentric annular flow. The Power-law eccentric annular flow results with flow index of 1.0 are verified with Newtonian fluid eccentric annular flow results. The parametric effects of flow index, Pipe/hole radius ratio, and eccentricity are investigated. We expect the development of a new model for flow in eccentric annular geometries to be an important new tool for application in oil - well drilling and production.

### At a glance: Figures

**Keywords:** bipolar coordinate, eccentric annulus, Power-law fluid, Point Successive over Relaxation

*American Journal of Numerical Analysis*, 2013 1 (1),
pp 1-7.

DOI: 10.12691/ajna-1-1-1

Received December 08, 2012; Revised January 16, 2013; Accepted March 03, 2013

**Copyright**© 2014 Science and Education Publishing. All Rights Reserved.

### Cite this article:

- Ebrahim, Nuha Hussein, Noaman El-Khatib, and Mariyamni Awang. "Numerical Solution of Power-law Fluid Flow through Eccentric Annular Geometry."
*American Journal of Numerical Analysis*1.1 (2013): 1-7.

- Ebrahim, N. H. , El-Khatib, N. , & Awang, M. (2013). Numerical Solution of Power-law Fluid Flow through Eccentric Annular Geometry.
*American Journal of Numerical Analysis*,*1*(1), 1-7.

- Ebrahim, Nuha Hussein, Noaman El-Khatib, and Mariyamni Awang. "Numerical Solution of Power-law Fluid Flow through Eccentric Annular Geometry."
*American Journal of Numerical Analysis*1, no. 1 (2013): 1-7.

Import into BibTeX | Import into EndNote | Import into RefMan | Import into RefWorks |

### 1. Introduction

Annular flow analysis is important in drilling and production in deviated and horizontal wells. In drilling, the annular velocity is one of two major variables in the process of cleaning solids (drill cuttings)* *from the wellbore. By maintaining the annular velocity above a certain rate (speed) in conjunction with the rheological properties of the drilling fluid, the wellbore can be kept clean of the drill cutting to prevent them from settling and causing drilling problems. Many investigators assumed the annulus to be “almost concentric,” however, the drill pipe is usually positioned eccentrically in the wellbore where a drill pipe has strong tendency to offset toward the low side because of gravitational effect. Eccentric annular flows are more complicated, and have been modeled under limiting assumptions ^{[1]}. This paper describes a numerical solution of fully developed laminar flow of Power-law non-Newtonian fluid in an eccentric annulus using bipolar coordinates transformation. The resulting model will be used in modeling of cuttings transport during drilling operation.

### 2. Background

The main function of drilling fluid is to carry the drilled particles “cuttings” generated by the drill bit to the surface through the annulus. Most drilling fluids exhibit non-Newtonian properties ^{[2, 3]}. Non-Newtonian behavior is characterized by nonlinear relationship between the shear stress and shear strain. Most drilling fluids are pseudoplastic as their shear thinning nature is very desirable in drilling.

The practical complications of flow of non-Newtonian fluid through eccentric annulus are: the eccentric geometry of the annulus (Figure 1); and the non linearity of the equations of motion of non-Newtonian fluids due to the non linear relationship between shear stress and shear rate.

**Fig**

**ure**

**1**

**.**Annular velocity Profile for a non-Newtonian Fluid with Eccentric Drillpipe [4]

Therefore, it is difficult to find exact analytical solutions, numerical methods are usually used to solve this problem. A number of studies have focused on the problem of flow of non-Newtonian fluids flow through eccentric annulus ^{[9, 10, 13]}. Two alternative approaches have been used in the previous studies. The first one uses the bi-polar coordinates to describe the eccentric annular geometry; the second one treats an eccentric annulus as a slot of variable height. Using a bipolar coordinate system and Green’s function, Heyda ^{[5]} presented analytical solutions for Newtonian fluid flow in an eccentric annulus in the form of an infinite series. He showed that the velocity profile of Newtonian fluid in laminar flow regime would differ dramatically from that of an eccentric annulus. Redberger and Charles ^{[6]} reported their numerical evaluations of the velocity profile and the volumetric flowrate for the Newtonian eccentric annular flow and concluded that the fluid velocity is higher in the enlarged region than in the reduced region of an eccentric annulus and the volumetric flowrate increases with increase of the offset of the inner tube for a given pressure gradient. Their velocity profile agreed with Heyda analytical solution.

Snyder and Goldstein ^{[7]} presented an analysis of fully developed laminar flow of Newtonian fluid in an eccentric annulus. An exact solution for the velocity distribution was presented in terms of infinite series summation. From this solution expressions were obtained for local shear stresses on the inner and outer surfaces of the annulus, friction factors based on the inner and outer surfaces, and the overall friction factor. Curves of these parameters were presented covering a range of eccentricity values and radius ratio values. Mitsuishi and Aoyagi ^{[8]} extended the bipolar coordinate method to the eccentric annular flow of non-Newtonian fluids using the Sutterby model and similar conclusions as those obtained by Redberger and Charles were reached. They validated the calculated velocities and the volumetric flowrates with the experimental data. Guckes ^{[9]} presented a systematic procedure for calculating the volumetric flowrate by using the bipolar coordinate approach and developed a series of dimensionless curves for the laminar eccentric annular flow of power-law fluids. However, after a lengthy treatment, Tosun ^{[10]} developed expressions for the volumetric flow-rate of Newtonian fluids through an eccentric annulus. Using the expressions for volumetric flow-rate in a concentric annuli developed by Bird et al. ^{[11]}, he presented the ratio of the volumetric flow rates in eccentric and concentric annuli in terms of the diameter ratio and the eccentricity, for a given pressure gradient. Their results showed that this ratio increases as the diameter ratio and the eccentricity increases. Haciislamoglu and Langlinais ^{[12]} used an iterative finite difference method to model Bingham fluids, but they did not provide information on computing times and numerical stability or code portability. However, limitations on the mapping used means that the methodology cannot* *be extended to handle boreholes with cuttings beds ^{[1]}.

Examining the analytical solution that approximate an eccentric annulus with equivalent slot ^{[13, 14, 15, 16, 17]}, it is can be realized that bipolar coordinates system should be used to obtain an exact solution for the eccentric annular geometry rather than approximate slot solution.

### 3. Methodology

The fluid flow in an annulus is idealized as steady, isothermal, fully developed laminar flow through a straight annulus consisting of an outer cylindrical casing () and an inner cylinder which is offset (i.e. eccentric). The geometry of the system is depicted in Figure 2.

The flow domain under consideration is not easy to describe in any of the classical coordinate system due to the irregular shape of the boundaries. Therefore a bipolar coordinate system is suggested.

**Fig**

**ure**

**2.**Physical domain in cylindrical coordinate

In bipolar coordinates system, the two cylindrical boundaries of the annulus coincide with two surfaces having constant values of (and , which can be expressed in terms of the annulus radius ratio k and the dimensionless eccentricity ξ. The other coordinate () represents a set of eccentric cylinders whose centers lie on the y-axis and which intersect orthogonally the boundaries of the annulus (Figure 3).

**Fig**

**ure**

**3.**Eccentric Annulus in Bipolar Coordinates [18]

The relationships needed to transform from Cartesian coordinates to bipolar coordinates are given by Speigel ^{[19]}

(1) |

(2) |

Where a is the pole of the bipolar coordinate system on the x-axis (a>0), which is expressed as

To achieve the scope of this paper, the work done was categorized into two sections

I. Development of model for Newtonian fluid flow in eccentric annulus.

II. Extend the developed model for non-Newtonian eccentric annular flow.

**3.1 Newtonian Eccentric Annular Flow**

For axial, laminar, steady flow of incompressible fluid, the z-direction momentum equation in bipolar coordinate system can be expressed in terms of velocity as:

(3) |

Here denotes the Newtonian viscosity and v the velocity

According to Guckes ^{[9]}, the outer and inner walls of an eccentric annulus can be expressed respectively in a bipolar coordinate system as:

(4) |

(5) |

Where is the non-dimensional eccentricity and is the radius ratio given by:

(6) |

Equation (3) is made dimensionless by introducing dimensionless velocity

(7) |

After the dimensionless analysis, equation (3) becomes

(8) |

This is so-called Poisson’s equation in mathematics. In an eccentric annulus, the axis that goes through the center points of the inner and outer pipes is the line-of-symmetry which divides the velocity profile into two identical parts. Therefore we shall solve for the velocity profile only in one half of the eccentric annulus ^{[20]}.

To build a system of linear equations involving the unknowns, we apply the above differential equation at the (i,j) grid point and approximate the second partial derivatives using central difference, setting

(9) |

(10) |

These approximations allow us to convert the partial differential equation to the finite difference equation choosing mesh spacing on a rectangular domain

(11) |

Where

Equation (11) can be rearranged to yield

(12) |

For a uniform grid, the final formulation of the finite-difference equation is:

(13) |

Where A, C, E, F, D, and B are given as:

And

These C, A, F**, **E, and D expressions refer to the coefficients associated with any set of left, right, top, bottom and central mesh points as shown in Figure 4. Because of this point configuration, equation (12) is called a five-point formula

**Fig**

**ure**

**4.**Eccentric annulus in transformed coordinates

The boundary conditions of the region are

v=0 at (at the pipe wall)

v=0 at (at the hole wall)

at

The problem is then to find values of the dependent variable V at internal locations. The matrix form of Equation (13) is pentadiagonal (most entries are zero, and the five diagonals contain all the nonzero entries).** **The basic idea when storing sparse matrices is to store only the non-zero entries as opposed to storing all entries. Fig. (5) demonstrated the final matrix form after applying boundary conditions. For a 2D-grid of mxn interior points, **V** has m x n grid points, but the** **coefficient matrix (A) has (m x n) x (m x n) entries. Hence, this yields linear system of equations as

(14) |

**Fig**

**ure**

**5.**Structure of pentadiagonal matrix coefficients A

In this paper point SOR method (Fig.6) is used to solve Equation (9). It is an iterative scheme that uses a relaxation parameter ω and is a generalization of the Gauss-Seidel method in the special case ω = 1.

**Fig**

**ure**

**6.**Flow chart of S.O.R. method

A FORTRAN computer program was written to obtain the velocity profile for Newtonian eccentric annular flow. To verify our numerical model, the results of Newtonian eccentric annular flow with 0.0001 eccentricity are compared with the concentric annular flow results. Results are shown in Figure 9. It shows that this program can be used for modeling Newtonian eccentric annular flow.

**3.2 Power-law Eccentric Annular Flow**

For axial, laminar, steady flow of incompressible fluid, the z-direction momentum equation in Cartesian coordinate system is:

(15) |

For Power-law model, the relationship between the shear stress and shear rate is

(16) |

Where K denotes the fluid consistency index

Then the governing equation of motion for Power-law fluids in eccentric annulus can be expressed in terms of velocity as:

(17) |

The z-direction momentum equation in bipolar coordinate system

(18) |

The dimensionless form of Equation (18) is

(19) |

Where is the dimensionless velocity given as:

(20) |

Equation (19) can be written as:

(21) |

After discretizing, Equation (21) will be:

(22) |

After grouping the terms of Equation (22), the discreitized differential equation constants will be:

(23) |

Equation (23) is strongly nonlinear, the direct solutions are not numerically possible, but exact solutions can be obtained iteratively.

### 4. Results and Discussion

The eccentricity of an annular geometry has a strong influence on the laminar flow velocity distribution. This is seen in Figure 7, where the dimensionless axial velocity profiles of Newtonian fluid, are plotted versus x direction (21) for 0.625 Pipe/hole radius ratio and at 0.6 dimensionless eccentricity. Due to the eccentricity, the velocity in the narrowest part of the annulusis reduced because the resistance to flow is increased as the gap between two pipes decreases. While at the widest gapthe fluid reaches its maximum velocity.

**Fig**

**ure**

**7.**Velocity distribution of Newtonian fluid flow in the widest and narrowest gabs of eccentric annulus

To verify the proposed model, firstly, the velocities at widest gap of an eccentric annulus with 0.0001 eccentricity were plotted with the velocity values at narrowest (Figure 8), it was found that they were identical that means eccentric annulus with small eccentricity resembles the concentric annular geometry. Secondly, the results obtained for Newtonian fluid flow through eccentric annulus with 0.0001 eccentricity were compared with the results of Newtonian fluid flow in concentric annulus (e=0.0). As seen in Figure 9, the results of eccentric annulus with 0.0001 eccentricity are in agreement with the concentric annulus (e=0.0) solution.

**Fig**

**ure**

**8.**Relationship between velocity profiles in the widest and narrowest gap of annulus with 0.0001 dimensionless eccentricity

**Fig**

**ure**

**9.**Velocity profiles of Newtonian fluid flow in concentric annulus and eccentric annulus with small eccentricity

The investigated model was extended to be used for the non-Newtonian fluid flow through eccentric annular geometry. To validate non-Newtonian fluid model, the results of Power-law flow through eccentric annulus with 0.0001eccentricity are compared with the experimental data by Nouri and Whitelaw (1994) for Power-law fluid flow in concentric annulus.

The agreement between the experimental data and the proposed model are clearly satisfactory (Figure 10).

**Fig**

**ure**

**10.**Comparision between proposed model normalized velocity data and Nouri and Whitelaw (1994) normalized velocity data

The parametric effect of eccentricity and pipe/hole radius ratio on the laminar flow velocity distribution was investigated. The eccentricity of an annular geometry has a strong influence on the laminar flow velocity distribution. This is seen in Figure 11, where the dimensionless axial velocity profiles of Power-law non-Newtonian fluid, are presented for 0.625 Pipe/hole radius ratio and at four different dimensionless eccentricities (0.05, 0.2, 0.5, and 0.6). For zero eccentricity, the velocity profiles on both sides of the inner cylinder are identical. With the increase of the eccentricity, the difference between the width of the wide and narrow gaps increases. The velocity increases in the wide and decreases in the narrow gap with increase of eccentricity.

Figure 12 illustrates the effect of radius ratio on the velocity distribution, results obtained for three different radius ratios (k=0.5, 0.625, and 0.7) at 0.6 dimensionless eccentricity. Larger is the radius ratio, lesser is the velocity.

**Fig**

**ure**

**11.**Effect of eccentricity on the velocity distribution of Power-law non-Newtonian fluid in the wide and narrow gap of annulus

**Fig**

**ure**

**12.**Effect of radius ratio on the velocity distribution of Power-law non-Newtonian fluid in the wide and narrow gap of annulus

### 5. Conclusion

The equation of motion in axial direction has been transformed into bipolar coordinate system. First, Newtonian model has been incorporated into the equation of motion. A numerical solution of the discretized equation of motion has been developed using Successive Over Relaxation method to obtain the velocity profile in eccentric annular geometry. Comparison between Newtonian concentric and eccentric annular flow with small eccentricity showed the validity of the numerical model. Second, the numerical model for Newtonian model has been extended to obtain the Power-law velocity profile.

The results of Power-law flow through eccentric annulus with small eccentricity were in agreement with the experimental data published in literature for Power-law fluid flow in concentric annulus. The velocity profiles have shown higher velocities in the wide gap of the annulus and lower velocities in the narrow gap of the annulus where the cuttings bed will be formed. Therefore, pipe eccentricity contributes to hole cleaning problems. It was found that the velocity increases with the increase of the flow index and decreases with the increase of the radius ratio.

### Nomenclature

*A*, *B*, *C, D, F, and G* = Coefficients of Pentadiagonal system Equation

*a*= is the pole of the bipolar coordinate system on the x-axis

*dP/dz* = Pressure gradient, Pa/m [Psi/ft]

*e*= offset distance between the centers of the inner and outer pipes m [in.]

*e** =dimensionless eccentricity

*h* =scale factor of the bipolar coordinate system

*K*=Flow consistency index, Pa. s ^{n} [lbf.s^{n}/ft^{2}]

*k* = Annulus aspect ratio

*n *= Flow behavior index, dimensionless

*R*_{i}= Drillpipe radius, m [in.]

*R*_{o}= Outer radius of annulus, m [in.]

*r *= Radius coordinate, m [in.]

*r*^{*}= radius ratio of the annulus

*v*= Axial velocity, m/s [ft/min]

*v*^{*} = Dimensionless velocity

x = the first transverse direction in the Cartesian

coordinate system

*Z= *axial coordinate, m [ft]

*Greek Symbols*

*η*_{i}=constant variable used to describe bipolar coordinate system

*θ** *= angular coordinate, (deg.)

*τ** *= Shear stress, Pa [Psi]

*ρ** *= fluid density, kg/m^{3} [lbm/gal]

*γ*= Shear rate, (*s*^{-1})

*μ* = Viscosity, Pa.s [cp]

*Subscript*

*i* =outer surface of inner cylinder

*o*= inner surface of outer cylinder

### Acknowledgment

Financial support from PETRONAS NASIONAL BERHAD gratefully acknowledged.

### References

[1] | Chin, W., Computational rheology for pipelines and annular flow, Gulf Publishing Company, Houston, Texasm, 2001. | ||

In article | |||

[2] | Alderman, N. Ram Babu, D. Hughes, T. and Maitland, G. “The rheological properties of oilwell drilling fluids,” Proc Xth Int. Cong. Rheology, Sydney, 140-142, 1998. | ||

In article | |||

[3] | Hamed, S. and Belhadri, M. “Rheological properties of biopolymers drilling fluids,” Journal of Petroleum Science and Engineering, 67 (3-4), 84-90, 2009. | ||

In article | CrossRef | ||

[4] | Carden, R. and Grace, R., Horizontal and directional drilling, PetroSkills. OGCI Co., Oklahoma, USA, 2007. | ||

In article | PubMed | ||

[5] | J. Heyda, “A green’s function solution for the case of laminar incompressible flow between non-concentric circular cylinders,” Journal of Franklin Inst. 267, 25-34, 1959. | ||

In article | CrossRef | ||

[6] | Redberger, P. and Charles, M., “Axial laminar flow in a circular pipe containing a fixed eccentric core,” Canadian J. of Chem. Eng., 40(4), 148-151, 1962. | ||

In article | CrossRef | ||

[7] | Snyder, W. and Goldstein, G., “An Analysis of Fully Developed Laminar Flow in an Eccentric Annuli,” AIChE Journal, 11, 462-467, 1965. | ||

In article | CrossRef | ||

[8] | Mitsuishi, N. and Aoyagi, Y., “Non-Newtonian Fluid Flow in an Eccentric Annulus,” Journal of Chemical Eng., 6(5), 402-08, 1973. | ||

In article | |||

[9] | Guckes, T.,“Laminar flow of non-Newtonian Fluids in an eccentric annulus,” Trans. ASME, 97(2), 498-506, 1974. | ||

In article | |||

[10] | Tosun, I., “Axial laminar Flow in an Eccentric Annulus: An Approximate Solution,” AICHE Journal, 30(5), 877-78, 1984. | ||

In article | CrossRef | ||

[11] | Bird, R. Stewart, W. and Lightfoot, E., Transport Phenomena, John Wiley & Sons, New York City, 1960 | ||

In article | |||

[12] | Haciislamoglu, M. and Langlinais, J., “Non-Newtonian Flow in Eccentric Annuli,” J. Energy Res. Tech. 112(3), 163-169, 1990. | ||

In article | CrossRef | ||

[13] | Tao, L. and Donovan, W., “Through Flow in Concentric and Eccentric Annuli of Fine Clearance With and Without Relative Motion of the Boundaries,” Trans., ASME, 77, 1291-1301, 1995. | ||

In article | |||

[14] | Vaughn, R., “Axial Laminar flow of Non-Newtonian Fluids in Narrow Eccentric Annuli,” SPEJ, 277-80, Trans., AIME, 234, 1965. | ||

In article | |||

[15] | Iyoho, A. and Azar, J., “An Accurate Slot Model for Non-Newtonian Flow through Eccentric Annuli,” SPEJ (Oct. 1981) 565-72, 1981. | ||

In article | |||

[16] | Uner, D. Ozden, C., and Tosun, I., “An Approximate Solution for Non-Newtonian Flow in Eccentric Annuli,” Ind. And eng. Chem., 27(4), 698-701, 1988. | ||

In article | |||

[17] | Luo, Y. and Peden, J., “Flow of Non-Newtonian Fluids Through Eccentric Annuli,” SPE Production Engineering, 91-96, 1990. | ||

In article | |||

[18] | Pillutla, J., “Laminar non-Newtonian flows in eccentric annuli with inner cylinder rotation,” Msc diss., University of Indhra, 2001. | ||

In article | |||

[19] | Speigel, M. R., Mathematical Handbook of Formulas and Tables, McGraw-Hill Book Co., pp. 125-128, 1968. | ||

In article | |||

[20] | Haciislamoglu, M. and Langlinais, J., “Effect of pipe Eccentricity on Surge Pressures,” Journal of Energy Resources Technology, 113, 157-160, 1991. | ||

In article | CrossRef | ||

[21] | Adariani, Y. H.: “Numerical simulation of laminar flow of non-Newtonian fluids in eccentric annuli,” Msc diss., University of Tulsa, 2005. | ||

In article | |||

[22] | Nouri, J.M., Whitelaw, J.H., “Flow of Newtonian and non-Newtonian fluids in a concentric annulus with rotation of the inner cylinder,” Journal of Fluids Engineering, 116, 1994, 821-827. | ||

In article | CrossRef | ||