Research Article  Open Access
Titus Ntow Ofei, "Effect of Yield Power Law Fluid Rheological Properties on Cuttings Transport in Eccentric Horizontal Narrow Annulus", Journal of Fluids, vol. 2016, Article ID 4931426, 10 pages, 2016. https://doi.org/10.1155/2016/4931426
Effect of Yield Power Law Fluid Rheological Properties on Cuttings Transport in Eccentric Horizontal Narrow Annulus
Abstract
Narrow annular drilling such as casingwhiledrilling technique is gaining popularity due to its ability to mitigate nonproductive time during oil and gas drilling operations. However, very little is known about the flow dynamics in narrow annular drilling. In this study, the EulerianEulerian twofluid model was used to examine the influence of Yield Power Law fluid rheological properties on cuttings transport in eccentric horizontal narrow annulus. The flow was assumed as fully developed, laminar, and transient state. The present simulation model was validated against experimental data, where a mean percent error of −1.2% was recorded. Results revealed an increase in the radial distribution of cuttings transport velocity in the wide annular region as the consistency index, , and the flow behavior index, , increase. Nonetheless, increasing the yield stress, , had insignificant effect on the cuttings transport velocity. Threedimensional profiles showed how cuttings preferred to travel in less resistant flow area, whereas cuttings concentration builds up in the narrow annular region. Furthermore, annular frictional pressure losses also increased as , , and increased. This study serves as a guide to properly optimize drilling fluid rheological properties for efficient cuttings transport and equivalent circulating density (ECD) management in narrow annular drilling.
1. Introduction
Drilling fluid rheological properties are important parameters which contribute to effective hole cleaning. Adari et al. [1] indicated that drilling fluid rheological properties are highly influential on cuttings transport; hence, care must be taken in predicting the optimum parameters to enhance better hole cleaning. The HerschelBulkley (HB) viscosity model, also called Yield Power Law (YPL) model, is known to correlate better to the drilling fluid viscosity curve than most other rheological models [2, 3].
In literature, most experimental studies on drilling fluids observed that the drilling fluid rheological curves (rheogram) conformed best to that of YPL fluid model. Ahmed and Miska [4] conducted extensive experimental study with polymerbased fluids in concentric and eccentric horizontal annuli without the presence of cuttings. The fluid rheological properties were accurately characterized using YPL model. The effects of inner pipe rotation, eccentricity, and flow rates on frictional pressure losses were analyzed. Other experimental studies [5–8] also accurately fitted their fluid rheological parameters using the Yield Power Law model. The authors proposed pressure loss equations for YPL fluid flow in pipes and annuli which is an integration of analytical, semianalytical, and empirical equations for laminar, transitional, and turbulent flows. Taghipour et al. [9] were among the first researchers to use YPL fluids to perform cuttings transport experimental study in inclined and horizontal annulus of diameter ratio, . They measured annular pressure losses, solids bed height, and drill string torque.
Several experimental and numerical studies have also evaluated the effects of drilling fluid rheological properties on cuttings transport. Cho et al. [10] observed that a decrease in flow behavior index, , resulted in a decrease of stationary bed, whereas moving bed layer increases. Adari et al. [1] also observed that the increase in the ratio of flow behavior index to consistency index () reduced cuttings bed height. In addition, the drilling fluid yield stress (YS) and plastic viscosity (PV) also influenced cuttings removal. Valluri et al. [11] presented their experimental findings that an increase in YS at constant flow rate without drillpipe rotation resulted in negligible cuttings bed erosion. Furthermore, Mohammadsalehi and Malekzadeh [12] also observed that a reduction in PV and YS resulted in a better hole cleaning at reduced flow rates. Horton et al. [13] explored how rheological properties of YPL fluids, bottomhole temperature, and fluid configuration in the wellbore affect the temperature profile in deep offshore wells. It was concluded that the yield stress, , value is the most important parameter that significantly affects the shutin time for the well to cool to 68°F, while the flow behavior index, , and consistency index, , have minor effect on both oilbased and waterbased YPL fluids.
In addition, other authors have utilized CFD approach to model YPL fluid through pipes and annuli. Bui [14] conducted a CFD study of YPL fluid flowing through pipes and concentric and eccentric annuli with tool joints. The author examined the effects of drillpipe rotation and tool joint on annular pressure losses. Ofei et al. [15] also employed CFD technique to model YPL fluid flowing in both concentric and eccentric narrow annuli. The authors analyzed the effects of bulk velocity, diameter ratio, inner pipe rotation speed, and eccentricity on the radial distributions of axial and tangential flow profiles. Mokhtari et al. [16] carried out CFD study to simulate the effects of YPL fluid rheological parameters on pressure losses in eccentric annulus. The authors established that an increase in the fluid yield stress, , increases the pressure loss substantially and creates a wider plug zone velocity. Furthermore, they indicated that as the fluid becomes more shear thinning, pressure loss decreases considerably.
It is evident that very rare or few studies on cuttings transport using YPL fluids exist in literature. As the oil and gas industries are developing the interest in the use of YPL fluids in drilling operations, there is need to conduct more research to understand the behavior of this fluid of interest on cuttings transport.
Several studies [17–28] have adopted CFD techniques to model solidliquid flow in annular geometries with much accuracy. This is due to the many advantages CFD has over experimental setup and empirical correlations in handling complex multiphase flow problems with unlimited physical and operational conditions.
The present study adopts the inhomogeneous EulerianEulerian twofluid model to analyze the effect of fluid rheological properties on cuttingsYPL flow in eccentric narrow horizontal annulus. A finite volume method is used to solve the continuity and momentum equations. CFD methods have been proven to be very effective for multiphase flow problems due to their ability in handling unlimited number of physical and operational conditions as well as eliminating the need for expensive experimental and materials setups.
Within the context of this study, the radial distributions of cuttings transport velocity were observed in both wide and narrow annular regions. The wide region (sector AA) is the gap between the top of the drillpipe and hole, while the gap between the bottom of the drillpipe and hole represents the narrow region (sector BB) as shown in Figure 1.
2. Materials and Methods
The inhomogeneous (EulerianEulerian) twofluid model in ANSYS CFX 14.0, where both the liquid and solid phases are considered interpenetrating continua, is adopted in this study. Other models such as EulerianLagrangian model, however, simulate the solid phase as a discrete phase and allows for particle tracking. The EulerianEulerian model is preferred to the EulerianLagrangian model due to its ability to handle high solid volume fractions. Furthermore, it accounts for particleparticle interaction and includes turbulence automatically. A drawback of this model is, however, the need for complex closure relations. The commercial software package ANSYS CFX 14.0 consists of the following five specialized components: (a) DesignModeler for building the geometries, (b) CFX Mesh for mesh generation (c) CFXPre for flow model definition, (d) CFXSolver for solving the governing equations, and (e) CFXPost for analyzing the results. The following continuity and momentum equations representing the twophase flow model are described for the sake of brevity.
2.1. Governing Equations
The solidliquid flow is assumed as (a) isothermal and (b) laminar and transient state.
2.1.1. Continuity Equations
The governing continuity equations for both liquid and solid phases could be expressed, respectively, as [29]where the solid and liquid phase volume fractions sum up to unity as
2.1.2. Momentum Equations
The forces acting on each phase and interphase momentum transfer term that models the interaction between each phase are given below [29].
For liquid phase,Similarly, for solid phase,
2.2. Other Closure Models
2.2.1. Interphase Drag Force Model
Considering spherical particles, the drag force per unit volume is given asFor densely distributed solid particles, where the solid volume fraction , the drag coefficient model by WenYu [30] may be utilized. This model is modified and implemented in ANSYS CFX 14.0 to ensure the correct limiting behavior in the inertial regime aswhere and .
For large solid volume fraction, , the drag model by Gidaspow [31] may be used with the interphase drag force per unit volume defined as
2.2.2. Lift Force Model
For spherical solid particles, ANSYS CFX employs the lift force model by Saffman [32, 33] asThe authors correlated the lift force for low Reynolds number past a spherical solid particle, where and . For higher range of solid particle Reynolds number, Saffman’s correlation was generalized by Mei and Klausner [34] aswhere , andand ,
2.3. Solid Phase Viscosity Model
Solid particles suspended in a liquid phase are affected by shear rate redistribution. The solid phase viscosity, , is related to the apparent suspension viscosity, , and liquid apparent viscosity, , asThe suspension viscosity is expressed in terms of relative viscosity, , asThe relative viscosity of more concentrated suspension with particleparticle interactions is given by [35]The first two terms of (13) represent Einstein [36] equation.
Equations (11) to (13) are coded into CFX library as expression language to compute the solid viscosity.
2.4. Carrier Fluid Viscosity
There is singularity problem associated with the classical YPL viscosity model at vanishing shear rate. To alleviate this, the proposed YPL viscosity function by Mendes and Dutra [37] is implemented in the present CFD study asThis viscosity function is numerically stable and devoid of discontinuity. Equation (14) is currently not available in the library of variables of ANSYS CFX 14.0; hence, it is modeled as expression language and coupled with the momentum equations.
2.5. Boundary Conditions and Meshing
At the inlet, a mass flow rate is specified, while a zerogauge pressure is specified at the outlet boundary. At the pipe walls, different boundary conditions were used for both liquid and solids. The usual noslip condition was imposed at the walls for the liquid phase, while for the solid phase, the freeslip condition was assumed at the walls to prevent the solid phase from adhering to the walls. This is consistent with real flow behavior of solid particles flowing near a solid boundary. The solid volume fraction in the domain was specified at the beginning of each simulation to correspond to the desired solid loading.
The annular 3D geometry of diameter ratio was meshed into structured hexahedral grids of approximately elements. Grid independent study of numerical results of pressure losses were carried out until variations in results were insignificant. In this study, the optimum grid sizes used for the radial, circumferential, and axial directions are . Figure 1 shows a section of the 2D and 3D annular mesh, while Figure 2 depicts the results from the grid independence study.
2.6. Description of Simulation Study
The simulation of the twophase solidliquid flow was set up in three dimensions using ANSYS CFX 14.0 with the transport equations solved using CFXSolver. The geometry dimensions, fluid rheological properties, solid properties, and operating parameters are presented in Table 1. To achieve a fully developed flow, the annular length, , must be longer than the hydrodynamic entrance length of the flow. In singlephase Newtonian fluids flowing in pipes, the hydrodynamic entrance length is presented as [38]However, for a twophase flow in annular gap with a nonNewtonian fluid, such expression as in (15) does not exist in literature. As a rule of thumb, the author adopted (15) by replacing the pipe diameter with a hydraulic diameter . It should be noted that a much longer annular length would only result in a longer simulation run.

The CFXSolver is based on a finite volume method in which the flow equations are integrated over each control volume. The advection scheme is set to high resolution to satisfy both accuracy and boundedness, where the blend factor, , is computed locally to be close to 1 without resulting in nonphysical values. The algorithm by Rhie and Chow [39] is used to solve the pressurevelocity coupling, since it can overcome pressurevelocity oscillations.
It should be noted that the complexity of the transport equations could not permit numerical convergence under steady state. However, all simulations were run in transient state. It is usually recommended that simulation of such steady state nature should first be run under transient state when it is difficult to attain convergence [23]. A very small time step of s was used in all simulations to help the solution converge; hence, the solutions finally converged when the convergence criterion was met at a root mean square (RMS) value of . An average of 200 time steps was required to achieve convergence with 1–10 loop iterations at each time step.
2.7. Validation of Simulation Model
There is scarcity of data for cuttings transport study using YPL fluids. The most recent and only experimental cuttings transport study is conducted by Taghipour et al., where annular pressure losses were measured for cuttingsYPL fluid flow in inclined (30°) eccentric wellbore and are used to validate the twophase CFD model. To validate the present model, the author simulated the experimental condition as presented by Taghipour et al. [9]. Table 2 presents the experimental setup with operational and rheological parameters. It is worth noting that only a section of the experimental flow loop length of 12 m was simulated. A simulated length of was chosen based on the calculated hydrodynamic length in (15) to ensure a fully developed flow as well as saving computational time. The flow was assumed as laminar and isothermal, while the drag and lift forces were modelled using the models by WenYu and Saffman as presented by (5) and (8), respectively. Figure 3 depicts the comparison between experimental and simulated pressure loss data. The comparison shows good agreement with a mean percentage error (MPE) between the calculated and measured pressure loss data of −1.2%, thus confirming the validity of the current model setup.

3. Results and Discussion
The simulation results reported here include the effects of yield stress, consistency index, and flow behavior index on cuttings transport velocity in YPL fluid in eccentric horizontal narrow annulus. Measurements of the radial distributions of the cuttings velocity profiles were taken along sectors AA and BB or at the positions 0° and 180° (see Figure 1), representing the wide and narrow regions of the eccentric wellbore, respectively. To generate the cuttings velocity profiles in CFXPost, a plane was first created parallel to the symmetry of the wellbore. Secondly, a polyline was also created with a boundary intersection method, where the boundary list of wellbore outlet was intersected with the plane. Finally, a chart was chosen where the velocity profiles were plotted using the polyline.
Studies have shown that YPL fluids used in the field have a wide range of values of rheological properties. These rheological properties could include the following range of values [13]: yield stress, = 0.048–50.06 Pa, consistency index, = 0.043–10.27 Pa·s^{n}, and flow behavior index, = 0.314–0.978. The range of rheological parameters used in this study is within the aforementioned range of values usually encountered during drilling operation using YPL fluids.
3.1. Yield Stress, , Effect on Cuttings Transport Velocity
The yield stress value evaluates the ability of the drilling fluid to suspend drilled cuttings during circulation. A high yield stress fluid is believed to transport cuttings better than fluid with low yield stress having similar density values. Figures 4(a) and 4(b) depict the effect of yield stress on the radial distributions of cuttings transport velocity in both wide and narrow annular regions, respectively. Figure 4(a) shows that cuttings transport velocity increases with increasing yield stress in the wide annular region. This is an indication of better carrying capacity of high yield stress fluids. Nonetheless, this improvement in cuttings transport velocity is less significant especially in the core region as yield stress increases. In contrast, the narrow annular region provides greater flow resistance to high yield stress fluids, thus the least to transport cuttings, as shown in Figure 4(b). In the wide annular region (sector AA), increasing from 2 to 8 Pa led to ~0.8% increment in the maximum cuttings velocity, whereas ~5.8% reduction in cuttings velocity occurred as increases from 2 to 8 Pa in the narrow annular region (sector BB).
(a)
(b)
3.2. Consistency Index, , Effect on Cuttings Transport Velocity
The consistency index of a drilling fluid is a rheological property related to the cohesion of the individual particles of the fluid, its ability to deform, and its resistance to flow. Figures 5(a) and 5(b) illustrate the consistency index’s effect on cuttings transport velocity in both wide and narrow annular regions, respectively. It is observed in Figure 5(a), that is, in the wide annular region, that when all other parameters are kept constant, increasing results in an increase in cuttings velocity.
(a)
(b)
It should be noted that is directly proportional to the fluid’s effective viscosity. Therefore, cuttings will travel less in low fluids with low effective viscosity as the cuttings tend to settle faster at the bottom of the annulus due to gravity. Whereas in high fluids, having high effective viscosity, cuttings will suspend for a longer period and, hence, travel faster and farther.
On the contrary, due to the flow restriction induced by the narrow annular region, cuttings travelling in high fluids had the least velocity as opposed to those travelling in low fluids (see Figure 5(b)). Investigations on the range of values presented in this study show that increasing from 1.7 to 6.3 Pa·s^{n} leads to ~11.0% increment in the maximum cuttings velocity in the wide annular region (sector AA), whereas ~37.7% reduction in cuttings velocity was recorded as increases from 1.7 to 6.3 Pa·s^{n} in the narrow annular region (sector BB).
3.3. Flow Behavior Index, , Effect on Cuttings Transport Velocity
The flow behavior index, , is a measure of a fluid’s shearthinning performance in both pseudoplastic and yield stress fluids; that is, apparent viscosity decreases with increasing shear rate. Fluids with higher shearthinning properties, , exhibit a plug (wider flat) flow region in the central core. Figures 6(a) and 6(b) demonstrate the effect of on cuttings transport velocity in both wide and narrow annular regions, respectively. In the wide annular region, as shown in Figure 6(a), cuttings transport velocity in the fluid domain decreases as the fluid becomes more shearthinning, especially in the core region. In the vicinity of the outer and inner pipe boundaries, high shearthinning fluids exhibit a flat radial cuttings velocity profile. This increases the high velocity zone towards the boundaries with improved cuttings transport. Within the study range, increasing from 0.31 to 0.75 led to increment of ~36.3% in the maximum cuttings velocity within the wide annular region (sector AA) as shown in Figure 6(a). In the narrow annular region (sector BB), cuttings travel faster in the radial distribution in low fluids compared to high fluids as shown in Figure 6(b). It is also observed that increasing values from 0.31 to 0.75 led to a significant reduction in maximum cutting velocity by ~84.6%.
(a)
(b)
3.4. ThreeDimensional Cuttings Velocity and Concentration Profiles
Figure 7 presents the threedimensional distributions of cuttings travelling in the annular geometry. Skewed cuttings velocity and concentration profiles were observed in the eccentric annulus as shown in Figures 7(a) and 7(b), respectively. These irregular distributions were due to the pipehole eccentricity which modified the flow by creating a high velocity zone in the wide annular region and a low or no velocity zone in the narrow annular region. The narrow region restricted the movement of the trapped cuttings and hence showed a very low velocity (see Figure 7(a)). Furthermore, high accumulation of cuttings bed can be observed in the narrow region of the annulus due to the effects of gravity and eccentricity as shown in Figure 7(b).
(a)
(b)
3.5. Annular Pressure Losses
Accurate estimation of annular pressure losses is very vital when designing drilling hydraulic programs, particularly the equivalent circulating densities (ECD) required for efficient transport of drilled cuttings from the wellbore to the surface. Figure 8 examines the influence of , , and on annular pressure loss. Figure 8(a) indicates that as the fluid becomes more shearthinning, that is, , the pressure loss decreases. There is a gradual rise in pressure as loss from = 0.31 to 0.51; however, a dramatic percentage increase of 316.4% is computed as increased from 0.51 to 0.75. Furthermore, Figure 8(b) also reveals an increase in pressure loss as increases from 1.7 to 6.3 Pa·s^{n} with approximately 132.6% increase in pressure loss. Last but not the least, increase in from 2 to 8 Pa had the least influence on pressure with nearly 13.6% increase.
(a)
(b)
(c)
4. Conclusions
A study on the effect of rheological parameters on cuttings transport velocity in YPL fluid flowing in eccentric narrow horizontal annulus is analyzed using EulerianEulerian twofluid CFD model. The proposed viscosity model for YPL fluid by Mendes and Dutra [37] was adopted since it is capable of avoiding numerical singularity difficulties at vanishing shear rate. For twophase cuttingsYPL fluid flow, CFD calculated pressure loss data were validated using recent but rare experimental data from literature. The present model was used to analyze the effects of yield stress, consistency index, and flow behavior index of the carrier YPL fluid on cuttings transport velocity.
Within the context of this study, the radial distributions of cuttings transport velocity were observed in both wide and narrow annular regions. The increase in yield stress of the carrier fluid from 2 to 8 Pa did not have much influence on the cuttings transport velocity in the wide region especially in the core region. However, there was much improvement in the cuttings transport velocity at the vicinity of the walls, indicating less cuttings bed accumulation. In the narrow gap, cuttings travelled much faster in low yield stress fluids as a result of their low resistance to flow compared to high yield stress fluids.
The study also revealed that carrier fluids with high consistency index value enhanced more cuttings transport especially in the wide annular region due to the fluid’s cuttings lifting ability. At the vicinity of the walls, cuttings travelled faster in low consistency index fluids due to the fluid’s tendency to high shearing and, hence, resulted in less cuttings bed formation. In the narrow region, however, cuttings travelled faster in low consistency index fluids due to their less resistance to flow.
Increasing the flow behavior index of the carrier fluid also showed much improvement in the cuttings transport velocity, especially in the core region of the wide annular gap. Meanwhile, carrier fluid with low flow behavior index transported cuttings better at the vicinity of the walls, an indication of less formation of cuttings bed. A reverse trend was, however, observed in the narrow annular gap, where there was significant cuttings transport in the carrier fluid with low flow behavior index. Threedimensional flow distribution profiles have shown the actual dynamics of cuttings travelling in the eccentric annulus, where most of the cuttings were inclined to travel in the wide margin with less stress.
The YPL fluid model has been shown to fit much better rheological data of drilling fluids in the oil and gas industry compared to the Bingham plastic and power law fluid models. Furthermore, YPL fluids with high yield stress values reduce convective heat loss in most especially high temperature wellbores, thus possessing the very rheological properties which are important for an insulating fluid to perform well. It is noteworthy that YPL fluids have viscosities that increase significantly as shearstrain rate diminishes; hence, by increasing the viscosity of the fluid, the drilling engineer can gain partial control over convection heat loss. More importantly, YPL fluids tend to have relatively low viscosity at high shear rates (shearthinning), making them easier to place initially, to bleed off pressure that may build up in the annulus equipped with venting capability, and to displace the drilling fluid in the event of well intervention. This study provides a guide to the drilling engineer on the selection of YPL fluid rheological properties which would enhance efficient transport of drilled cuttings in narrow horizontal wellbores.
It was further observed that the rheology of YPL fluids has significant effect on the annular pressure losses. Proper caution must be taken in selecting the fluid rheology to ensure efficient cuttings transport while maintaining a bottomhole pressure which will not fracture the formation.
Nomenclature
:  Drag coefficient (—) 
:  Solid particle mean diameter (m) 
:  Outer diameter of inner pipe (m) 
:  Inner diameter of outer pipe (m) 
:  Hydraulic diameter, (m) 
:  Offset distance (m) 
:  Gravity vector (m/s^{2}) 
:  Liquid phase volume fraction (—) 
:  Solid phase volume fraction (—) 
:  Consistency index (Pa·s^{n}) 
:  Annular geometry length (m) 
:  Hydrodynamic length (m) 
:  Interphase momentum transfer 
:  Drag force per unit volume (N/m^{3}) 
:  Lift force per unit volume (N/m^{3}) 
:  Flow behavior index (—) 
:  Fluid Reynolds number (—) 
:  Solid particles Reynolds number (—) 
:  Vorticity Reynolds number (—) 
:  Solid particle pressure (Pa) 
:  Radial direction 
:  Normalized radial distance 
:  Rate of penetration (m/s) 
:  Cuttings velocity at any radial distance (m/s) 
:  Bulk fluid velocity (m/s) 
:  Fluid phase velocity vector (m/s) 
:  Solid phase velocity vector (m/s) 
:  Axial direction. 
:  Eccentricity 
:  Fluid phase density (kg/m^{3}) 
:  Solid phase density (kg/m^{3}) 
:  Viscous stress tensor (Pa) 
:  Diameter ratio 
:  Zero shear rate viscosity (Pa·s) 
:  Viscosity defined in (14) (Pa·s) 
:  Liquid viscosity (Pa·s) 
:  Apparent viscosity (Pa·s) 
:  Relative viscosity 
:  Solid viscosity (Pa·s) 
:  Suspension viscosity (Pa·s) 
:  Specific volume (m^{3}/kg) 
:  Angular velocity (1/min) 
:  Shear rate (1/s) 
:  Rotation vector (1/min) 
:  Circumferential direction. 
Competing Interests
The author declares that there are no competing interests regarding the publication of this paper.
References
 R. B. Adari, S. Miska, E. Kuru, P. Bern, and A. Saasen, “Selecting drilling fluid properties and flow rates for effective hole cleaning in highangle and horizontal wells,” in Proceedings of the SPE Annual Technical Conference and Exhibition, pp. 273–281, Dallas, Tex, USA, October 2000. View at: Publisher Site  Google Scholar
 A. Mihalakis, P. Makri, V. C. Kelessidis, G. Christidis, A. Foscolos, and K. Papanikolaou, “Improving rheological and filtration properties of drilling muds with addition of greek lignite,” in Proceedings of the 7th National Congress on Mechanics (HSTAM '04), Chania, Greece, 2004. View at: Google Scholar
 R. Maglione, G. Robotti, and R. Romagnoli, “Insitu rheological characterization of drilling mud,” SPE Journal, vol. 5, no. 4, pp. 377–386, 2000. View at: Publisher Site  Google Scholar
 R. Ahmed and S. Miska, Experimental Study and Modeling of Yield PowerLaw Fluid in Annuli with Drillpipe Rotation, Society of Petroleum Engineers, Richardson, Tex, USA, 2008. View at: Publisher Site
 V. C. Kelessidis, R. Maglione, C. Tsamantaki, and Y. Aspirtakis, “Optimal determination of rheological parameters for HerschelBulkley drilling fluids and impact on pressure drop, velocity profiles and penetration rates during drilling,” Journal of Petroleum Science and Engineering, vol. 53, no. 34, pp. 203–224, 2006. View at: Publisher Site  Google Scholar
 K. Founargiotakis, V. C. Kelessidis, and R. Maglione, “Laminar, transitional and turbulent flow of HerschelBulkley fluids in concentric annulus,” Canadian Journal of Chemical Engineering, vol. 86, no. 4, pp. 676–683, 2008. View at: Publisher Site  Google Scholar
 V. C. Kelessidis, P. Dalamarinis, and R. Maglione, “Experimental study and predictions of pressure losses of fluids modeled as HerschelBulkley in concentric and eccentric annuli in laminar, transitional and turbulent flows,” Journal of Petroleum Science and Engineering, vol. 77, no. 34, pp. 305–312, 2011. View at: Publisher Site  Google Scholar
 O. Erge, E. M. Ozbayoglu, S. Z. Miska et al., The Effects of Drillstring Eccentricity, Rotation, and Buckling Configurations on Annular Frictional Pressure Losses While Circulating Yield Power Law Fluids, Society of Petroleum Engineers, 2014. View at: Publisher Site
 A. Taghipour, B. Lund, J. D. Ytrehus et al., “Experimental study of friction and cutting transport in noncircular borehole geometry,” in Proceedings of the SPE/IADC Middle East Drilling Technology Conference & Exhibition, pp. 821–835, Society of Petroleum Engineers, Dubai, UAE, October 2013. View at: Publisher Site  Google Scholar
 H. Cho, S. N. Shah, and S. O. Osisanya, “A threesegment hydraulic model for cuttings transport in horizontal and deviated wells,” in Proceedings of the SPE/Petroleum Society of CIM International Conference on Horizontal Well Technology, SPE/PSCIM 65488, Calgary, Canada, 2000. View at: Google Scholar
 S. G. Valluri, S. Z. Miska, R. Ahmed, M. Yu, and N. Takach, “Experimental study of effective hole cleaning using ‘sweep’ in horizontal wellbores,” in Proceedings of the SPE Annual Technical Conference and Exhibition, SPE 101220, San Antonio, Tex, USA, 2006. View at: Google Scholar
 M. Mohammadsalehi and N. Malekzadeh, “Optimization of hole cleaning and cutting removal in vertical, deviated and horizontal wells,” in Proceedings of the Asia Pacific Oil and Gas Conference and Exhibition, SPE 143675, Society of Petroleum Engineers, Jakarta, Indonesia, September 2011. View at: Publisher Site  Google Scholar
 R. L. Horton, T. S. Froitland, W. E. Foxenberg, and D. A. Knox, “A new yield power law analysis tool improves insulating annular fluid design,” in Proceedings of the International Petroleum Technology Conference, Doha, Qatar, November 2005. View at: Publisher Site  Google Scholar
 B. Bui, “Modeling the effect of pipe rotation on pressure loss through tool joint,” in Proceedings of the SPETT Energy Conference and Exhibition, SPE157982MS, Society of Petroleum Engineers, PortofSpain, Trinidad, June 2012. View at: Publisher Site  Google Scholar
 T. N. Ofei, S. Irawan, W. Pao, and R. E. Osgouei, “Modified Yield PowerLaw fluid flow in narrow annuli with inner rotating pipe,” Canadian Journal of Chemical Engineering, vol. 93, no. 1, pp. 150–165, 2015. View at: Publisher Site  Google Scholar
 M. Mokhtari, M. Ermila, and A. N. Tutuncu, “Accurate bottomhole pressure for fracture gradient prediction and drilling fluid pressure programpart 1,” in Proceedings of the 46th U.S. Rock Mechanics/Geomechanics Symposium, American Rock Mechanics Association, Chicago, Ill, USA, June 2012. View at: Google Scholar
 T. N. Ofei, S. Irawan, and W. Pao, “CFD method for predicting annular pressure losses and cuttings concentration in eccentric horizontal wells,” Journal of Petroleum Engineering, vol. 2014, Article ID 486423, 16 pages, 2014. View at: Publisher Site  Google Scholar
 X.H. Zhu, C. Sun, and H. Tong, “Distribution features, transport mechanism and destruction of cuttings bed in horizontal well,” Journal of Hydrodynamics, vol. 25, no. 4, pp. 628–638, 2013. View at: Publisher Site  Google Scholar
 L.C. Qiu and C.Y. Wu, “A hybrid DEM/CFD approach for solidliquid flows,” Journal of Hydrodynamics, vol. 26, no. 1, pp. 19–25, 2014. View at: Publisher Site  Google Scholar
 S.M. Han, Y.K. Hwang, N.S. Woo, and Y.J. Kim, “Solidliquid hydrodynamics in a slim hole drilling annulus,” Journal of Petroleum Science and Engineering, vol. 70, no. 34, pp. 308–319, 2010. View at: Publisher Site  Google Scholar
 S. A. Hashemi, A. Sadighian, S. I. A. Shah, and R. S. Sanders, “Solid velocity and concentration fluctuations in highly concentrated liquidsolid (slurry) pipe flows,” International Journal of Multiphase Flow, vol. 66, pp. 46–61, 2014. View at: Publisher Site  Google Scholar
 J. Capecelatro and O. Desjardins, “EulerianLagrangian modeling of turbulent liquidsolid slurries in horizontal pipes,” International Journal of Multiphase Flow, vol. 55, pp. 64–79, 2013. View at: Publisher Site  Google Scholar
 M. Eesa and M. Barigou, “Horizontal laminar flow of coarse nearlyneutrally buoyant particles in nonNewtonian conveying fluids: CFD and PEPT experiments compared,” International Journal of Multiphase Flow, vol. 34, no. 11, pp. 997–1007, 2008. View at: Publisher Site  Google Scholar
 Z.M. Wang, X.L. Guo, M. Li, and Y.K. Hong, “Effect of drillpipe rotation on borehole cleaning for extended reach well,” Journal of Hydrodynamics, vol. 21, no. 3, pp. 366–372, 2009. View at: Publisher Site  Google Scholar
 K. Ekambara, R. S. Sanders, K. Nandakumar, and J. H. Masliyah, “Hydrodynamic simulation of horizontal slurry pipeline flow using ANSYSCFX,” Industrial and Engineering Chemistry Research, vol. 48, no. 17, pp. 8159–8171, 2009. View at: Publisher Site  Google Scholar
 X. Sun, K. Wang, T. Yan, S. Shao, and J. Jiao, “Effect of drillpipe rotation on cuttings transport using computational fluid dynamics (CFD) in complex structure wells,” Journal of Petroleum Exploration and Production Technology, vol. 4, no. 3, pp. 255–261, 2014. View at: Publisher Site  Google Scholar
 T. N. Ofei, S. Irawan, and W. Pao, “Modelling of pressure drop in eccentric narrow horizontal annuli with the presence of cuttings and rotating drillpipe,” International Journal of Oil, Gas and Coal Technology, vol. 9, no. 1, pp. 39–60, 2015. View at: Publisher Site  Google Scholar
 M. Sommerfeld and S. Lain, “Parameters influencing dilutephase pneumatic conveying through pipe systems: a computational study by the Euler/Lagrange approach,” Canadian Journal of Chemical Engineering, vol. 93, no. 1, pp. 1–17, 2015. View at: Publisher Site  Google Scholar
 B. G. M. van Wachem and A. E. Almstedt, “Methods for multiphase computational fluid dynamics,” Chemical Engineering Journal, vol. 96, no. 1–3, pp. 81–98, 2003. View at: Publisher Site  Google Scholar
 C. Y. Wen and Y. H. Yu, “Mechanics of fluidization,” Chemical Engineering Progress Symposium Series, vol. 62, pp. 100–111, 1966. View at: Google Scholar
 D. Gidaspow, Multiphase Flow and Fluidization, Academic Press, Boston, Mass, USA, 1994. View at: MathSciNet
 P. G. Saffman, “The lift on a small sphere in a slow shear flow,” Journal of Fluid Mechanics, vol. 22, no. 2, pp. 385–400, 1965. View at: Publisher Site  Google Scholar
 P. G. Saffman, “The lift on a small sphere in a slow shear flow—corrigendum,” Journal of Fluid Mechanics, vol. 31, no. 03, p. 624, 1968. View at: Publisher Site  Google Scholar
 R. Mei and J. F. Klausner, “Shear lift force on spherical bubbles,” International Journal of Heat and Fluid Flow, vol. 15, no. 1, pp. 62–65, 1994. View at: Publisher Site  Google Scholar
 D. G. Thomas, “Transport characteristics of suspension: VIII. A note on the viscosity of Newtonian suspensions of uniform spherical particles,” Journal of Colloid Science, vol. 20, no. 3, pp. 267–277, 1965. View at: Publisher Site  Google Scholar
 A. Einstein, “Eine neue Bestimmung der Moleküldimensionen [AdP 19, 289 (1906)],” Annalen der Physik, vol. 14, supplement 1, pp. 229–247, 2005. View at: Publisher Site  Google Scholar
 P. R. S. Mendes and S. S. E. Dutra, “A viscosity function for viscoplastic liquids,” Annual Transactions of the Nordic Rheology Society, vol. 12, pp. 183–188, 2004. View at: Google Scholar
 Y. Cengel and J. A. Cimbala, Fluid Mechanics: Fundamentals and Applications, McGrawHill Higher Education, 2004.
 C. M. Rhie and W. L. Chow, “Numerical study of the turbulent flow past an airfoil with trailing edge separation,” AIAA Journal, vol. 21, no. 11, pp. 1525–1532, 1983. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2016 Titus Ntow Ofei. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.