# Laminar, transitional and turbulent flow of Herschel--Bulkley fluids in concentric annulus.

INTRODUCTIONNon-Newtonian flow in annuli is encountered in many situations of importance in a variety of industries and particularly in oil-well drilling. For many non-Newtonian fluids, the two parameters models of Bingham plastic (Bingham, 1922) and of power law (Skelland, 1967; Govier and Aziz, 1972; Bourgoyne et al., 1991) are used most often because of their simplicity and the fair description of the fluid rheograms. The nonlinear three parameters model proposed by Herschel and Bulkley (1926) has not been used widely until very recently although it has been shown to fit much better rheological data of aqueous clay slurries and of drilling fluids (Fordham et al., 1991; Hemphil et al., 1993; Maglione and Ferrario, 1996; Kelessidis et al., 2005, 2007). Reasons for the non-frequent use were not only the complexity in derivation of the model's three parameters (Nguyen and Boger, 1987; Hemphil et al., 1993) but also the fact that analytical solutions for laminar or turbulent flow in annuli are not possible, requiring either graphical or trial-and-error solutions (Govier and Aziz, 1972; Hanks, 1979; Fordham et al., 1991). The on-line use, however, of personal computers enabled Maglione et al. (1999, 2000), Bailey and Peden (2000) and Becker et al. (2003) to utilize Herschel-Bulkley rheological model in fluid mechanics computations of drilling fluids.

Laminar flows of Bingham plastic and power law fluids in concentric annuli have been analyzed by Fredrickson and Bird (1958) while laminar flow of Herschel-Bulkley fluids in concentric annuli has been investigated by Hanks (1979) who provided a non-analytical solution, although with some errors (Buchtelova, 1988), together with several tables covering various annulus aspect ratios. Bird et al. (1983) provided an overview of solutions for the flow of yield-pseudoplastic fluids in various conduits. Analytical solutions for power-law and yield-pseudoplastic fluids, but not for Herschel-Bulkley fluids, in concentric annuli have been presented by Gucuyener and Mehmetoglu (1992), David and Filip (1995) and Filip and David (2003). Analytical and numerical codes covering laminar flow of non-Newtonian fluids in both concentric and eccentric annuli have been presented by Haciislamoglu and Langlinais (1990), Walton and Bittleston (1991), Szabo and Hassager (1992), Yang and Chukwu (1995), Hussain and Sharif (1997), Siginer and Bakhtiyarov (1998), Meuric et al. (1998), Fang et al. (1999) and Hussain and Sharif (2000), while Round and Yu (1993) have studied numerically laminar non-Newtonian flows in the entrance of a concentric annuli. Fordham et al. (1991) provided an analysis leading to a numerical code for solving the annulus situation aswell as the case of slot approximation for laminar flow, together with some experimental data. Chin (2001) has provided a general simulator covering laminar non-Newtonian flows in complex geometries. Measured and computed velocity profiles for laminar flow of yield pseudoplastic fluids in concentric and eccentric annuli have been given also by Escudier et al. (2002a,b).

No analytical solutions exist for turbulent flow of Herschel-Bulkley fluids in annuli and resort is made to experimental data with the use of the friction factor. Heywood and Cheng (1984) have reported variations of predictions of different proposed correlations of as much as [+ or -] 50% for turbulent flow of Herschel-Bulkley fluids in pipes. Harnett and Kostic (1990) reported that the best approach for turbulent flow of Herschel-Bulkley fluids in pipes was through the use of the Metzner and Reed (1955) graph. No such studies have been performed for the case of concentric annuli except the study by Reed and Pilehvari (1993) where a fairly complex model covering laminar, transitional and turbulent flow of Herschel-Bulkley fluids flowing in concentric annuli has been presented. Reed and Pilehvari (1993) proposed modifications to the diameters of the conduits to account for the non-Newtonian effect, using the generalized power law model of Metzner and Reed (1955) and the friction factor suggested for pipe flow by Dodge and Metzner (1959). Their model described very well their experimental data. Some confirmation of the Reed and Pilehvari model was also provided by Subramanian and Azar (2000) but there were cases for which predictions were at odds with the measurements. The approach of Reed and Pilehvari has been followed also by Merlo et al. (1995), Maglione et al. (2000) and Bailey and Peden (2000) to cover all flow regimes for flow of Herschel-Bulkley fluids and of generalized non-Newtonian fluids in concentric annuli.

A consistent and accurate methodology for predicting pressure drop for laminar, transitional and turbulent flow as well as the transition points from laminar to turbulent flow for the flow of Herschel-Bulkley fluids in concentric annuli, without resorting to heavy numerical computations, does not exist and this is the scope of the present paper. And as the slot model has been used often to approximate flows in concentric annuli of large diameter ratios, normally greater than 0.3 (Guillot and Dennis, 1988; Guillot, 1990; Bourgoyne et al., 1991; Fordham et al., 1991), the slot model will be used in the present approach.

[FIGURE 1 OMITTED]

THEORY

The Herschel-Bulkley model is given by:

[tau] = [[tau].sub.y] + K [(du / dy).sup.n], for [tau] > [[tau].sub.y] and du / dy > 0

[tau] = [[tau].sub.y] + K [(-du / dy).sup.n], for [tau] > [[tau].sub.y] and du / dy < 0 (1)

The geometry of the system is depicted in Figure 1. There is a central core of the fluid which moves as a rigid plug if the shear stress levels are smaller than the yield stress of the fluid. The flow rate for a slot of width, w, is then given by Grinchik and Kim (1974), Fordham et al. (1991) and Kelessidis et al. (2006):

q = [(d[p.sub.f]/dL / K).sup.m] 2w[(h/2).sup.m+2] [(1 - [xi]).sup.m+1] / (m + 1)(m + 2) [[xi] + (m + 1)] (2)

where m = 1/n and the dimensionless shear rate, [xi], is given by:

[xi] = [[tau].sub.y] / [[tau].sub.w]. (3)

The flow rate can then be given in terms of pressure drop and the annulus geometric parameters as:

q = [pi][([r.sup.2.sub.2] - [r.sup.2.sub.1]) ([r.sub.2] - [r.sub.1]).sup.1+m] [((1/K)(d[p.sub.f] / dL))].sup.m] [(1 - ([[tau].sub.y]/[([r.sub.2] - [r.sub.1])/2](d[p.sub.f]/dL))).sup.1+m] / [2.sup.m](m + 1)(2m + 4) x ([[tau].sub.y] / [([r.sub.2] - [r.sub.1])/2](d[p.sub.f]/dL) + m + 1) (4)

while the mean velocity, designated as, V, is given by:

V = [(d[p.sub.f]/dL / K).sup.m] [(h/2).sup.m+1] (1 - [xi]).sup.m+1] ([xi] + m + 1) / (m + 1)(m + 2). (5)

Using Equation (5), the velocity profiles can be recast in terms of the dimensionless shear rate as follows:

u/V = (m + 2) {[(1 - [xi]).sup.m+1] - [(1 - [xi] - (y/(h/2))).sup.m+1] / [(1 - [xi]).sub.m+1] ([xi] + m + 1)}; 0 [less than or equal to] y [less than or equal to] h / 2 (1 - [xi]) (6)

u/V = m + 2 / [xi] + m + 1 h/2 (1 - [xi]) [less than or equal to] y [less than or equal to] h/2 (1 + [xi]) (7)

u/V = (m + 2) {[(1 - [xi]).sup.m+1] - [(-1 + [xi] + (y/(h/2))).sup.m+1] / [(1 - [xi]).sup.m+1] ([xi] + m + 1)}; h/2 (1 + [xi]) [less than or equal to] y [less than or equal to] h. (8)

Laminar flow is amenable to analytical solution but for transition and turbulent flows, empirical equations must be used. Following Metzner and Reed (1955), use is made of the local power law parameters. The relevant Rabinowitsch (1929) equation for flow in an annulus, modelled as a slot, can be shown to be:

[[??].sub.Nw] = f ([[tau].sub.w]) = 12V / ([d.sub.2] - [d.sub.1]) = [[tau].sub.w]/3 d/d[[tau].sub.w] (12V/[d.sub.2] - [d.sub.1]) + 2/3 (12V / [d.sub.2] - [d.sub.1]). (9)

This equation shows that the shear stress at the wall, [[tau].sub.w], depends only on the Newtonian shear rate at the wall [[??].sub.Nw]. The local power law parameters, K', n', are then defined by:

n' = d ln ([[tau].sub.w])/d ln ([[??].sub.Nw]) (10)

and

[[tau].sub.w] = K'[([[??].sub.Nw]).sup.n'] (11)

The Herschel-Bulkley wall shear rate, [[??].sub.w] is then related to the Newtonian shear rate, [[??].sub.Nw], by:

[[??].sub.w] = 2n' + 1/3n' [[??].sub.Nw]. (12)

The generalized Reynolds number for the flow of Herschel-Bulkley fluid in an annulus then becomes, with [[mu].sub.e] = [tau]/[[??].sub.w]:

[Re.sub.MRa] = [rho]V ([d.sub.2] - [d.sub.1])/[[mu].sub.e] = [[rho][V.sup.2-n'] ([d.sub.2] - [d.sub.1]).sup.n'] / K'[(12).sup.n'-1 . (13)

The relationship between (n') and the Herschel-Bulkley parameters can be derived by combining the laminar flow solution (Equation (2)) with the definition of (n'), Equation (10), as it is done in the case of pipe flow, giving finally:

n' = n(1 - [xi])(n[xi] + n + 1) / 1 + n + 2n[xi] + 2[n.sup.2][[xi].sup.2] (14)

while the K' parameter can be shown to be:

K' = [[tau].sub.y] + K[((2n' + 1/3n')[[??].sub.Nw]).sup.n] / [([[??].sub.Nw]).sup.n'] (15)

Equations (14) and (15) fully define the generalized power law parameters of Herschel-Bulkley fluids flowing in an annulus, modelled as a slot, and are functions of the Herschel-Bulkley rheological parameters, [[tau].sub.y], K, n, and of the particular flow situations, through [[tau].sub.w].

The friction factor, f, is estimated, after adjusting the pipe flow equation proposed by Dodge and Metzner (1959) to flow in annuli, as:

1/[square root of (f)] = 4/[(n').sup.0.75] log [[Re.sub.MRa] [f.sup.1-n'/2]] - 0.395/[(n').sup.1.2]. (16)

Guillot and Dennis (1988) have reported good agreement of this equation with their data for flow of drilling fluids and of cement slurries in concentric annuli.

The pressure drop is then computed with:

d[p.sub.f]/dL = 2f[rho][V.sup.2]/[d.sub.2] - [d.sub.1]. (17)

A transition region, instead of a single transition point, is defined for a range of Reynolds numbers, with [Re.sub.1] denoting the beginning of transition and [Re.sub.2] the end of transition. Guillot (1990) has suggested to take these limits from the Dodge and Metzner (1959) graph, thus giving:

[Re.sub.1] = 3250 - 1150(n') (18)

[Re.sub.2] = 4150 - 1150(n'). (19)

The friction factor for the transition region, [f.sub.tr], is estimated at the particular Reynolds number, [Re.sub.MRa], with linear interpolation between the laminar friction factor at [Re.sub.1], [f.sub.l1], and the turbulent friction factor, [f.sub.t2], at [Re.sub.2]. Other interpolation approaches have not proven to provide better estimates. Hence, [f.sub.l1] can easily be shown to be:

[f.sub.l1] = 24/[Re.sub.1] (20)

while [f.sub.t2] is taken from the solution of Equation (16) with [Re.sub.MRa] = [Re.sub.2]. Finally, the transitional friction factor is estimated as:

[f.sub.tr] = [f.sub.l1] + ([Re.sub.MRa] - [Re.sub.1]) x ([f.sub.t2] - [f.sub.l1]) / ([Re.sub.2] - [Re.sub.1]). (21)

Derivation of pressure drop from the flow rate requires iterative solution. Similarly, derivation of the flow rate for a given pressure drop, normally used for design purposes, although gives a direct solution for laminar flow, it still requires iteration because the type of flow, laminar, turbulent or transition, is not known a priori and is known only after computing the pressure drop which gives the wall shear stress, [[tau].sub.w], on which the local power law parameter, n', depends.

The solution procedure for the former case will be given here below, as the latter case can be easily resolved. The current approach uses a bisection method and the flow diagram is given in Figure 2. Two values for [[tau].sub.w] are assumed, one close to the yield stress and the second at very high value. The third value is then the arithmetic mean of the two. The Newtonian shear rate at the wall is computed together with the dimensionless shear stress [xi] from Equations (3) and (9). The local power law parameters are then computed, n' from Equation (14) and K' from Equation (15). Then the generalized Reynolds number is estimated from Equation (13). The limits for laminar and turbulent flow can be computed and accordingly, the friction factor is computed. The flow rate is then estimated from:

q = A x [square root of (2 x [[tau].sub.w] / f x [rho])] (22)

The procedure converges when [absolute value of q - [q.sub.given]] [less than or equal to] [epsilon]. Following the above procedure, any flow problem can be solved without prior knowledge of the type of flow (laminar, transitional, or turbulent) and the procedure, based on theoretical and semi-theoretical grounds, has proven to work as it will be demonstrated later.

[FIGURE 2 OMITTED]

Many times is of great interest to know the start and end points of laminar to transition and transition to turbulent flow, that is, [Re.sub.1] and [Re.sub.2]. For accomplishing this, the procedure shown in Figure 3 is followed. Input data for the procedure are the geometry ([d.sub.2]-[d.sub.1]), fluid rheology and density ([rho]). The average velocity, V, is considered as variable parameter.

The first step consists in assuming that the local power law index equals the Herschel-Bulkley index, n' = n. Then, the Newtonian shear rate at the annulus wall is calculated from Equation (9). Following this, the Herschel-Bulkley shear rate at the wall, the shear stress at the wall, the dimensionless shear stress, [[xi].sub.1], and the local power-law flow-behaviour index [n'.sub.1], are calculated. The local consistency index K' and the generalized Reynolds number [Re.sub.MRa] are then calculated. Subsequently, the local wall shear rate and the shear stress, the dimensionless shear stress, [[xi].sub.2], and the local power-law flow-behaviour index [n'.sub.2], are again calculated. The procedure is repeated until the following conditions are satisfied, for the beginning of transition, [V.sub.1] and [Re.sub.MRa1] and for the end of transition for [V.sub.2] and [Re.sub.MRa2]:

[Re.sub.MRa1] - 3250 + 1150 x [n'.sub.j] [less than or equal to] [epsilon], (23)

[Re.sub.MRa2] - 4150 + 1150 x [n'.sub.j] [less than or equal to] [epsilon] (24)

The error in the calculation [epsilon], is chosen to be smaller than [10.sup.-4]. It has been seen that differences between computed ([[??].sub.wj-1] and [[??].sub.wj]) and ([[tau].sub.wj-1] and [[tau].sub.wj]) values become trivial after j = 4.

[FIGURE 3 OMITTED]

RESULTS

The accuracy of the predictions of the proposed model has been determined by comparing with experimental data and simulator results from works reported previously. More specifically, experimental data for laminar flow in concentric annulus, presented by Fordham et al. (1991), Okafor and Evers (1992) and Hansen et al. (1999) will be used for laminar flow. Similarly, controlled field data provided by Langlinais et al. (1983) will be used to cover all three flow regime types, laminar, transitional and turbulent flows. Comparison of the predictions of this work with predictions from full simulators (Hanks, 1979; Fordham et al., 1991) will also be performed.

Comparison of predictions from the approach presented in this work with experimental results of Okafor and Evers (1992) is shown in Figure 4 for an annulus aspect ratio of 0.62. There is excellent match of the predictions of this work with the experimental results. For a different fluid but the same annuli, the comparison with experimental results of Okafor and Evers (1992) is shown in Figure 5. There is excellent match with the data up to the velocity of 0.65 m/s, which corresponds to Reynolds number of 541. There is deviation for the last four points (maximum deviation is of the order of ~17%), having velocities between 1.0 and 1.25 m/s which correspond to Reynolds numbers 992 and 1275, respectively, thus the flow is still in the laminar regime. The slope is similar but there is constant under-prediction the cause of which is unknown.

In Figure 6, the comparison of predictions with experimental laminar flow results of Hansen et al. (1999) is shown. Very good match is seen for the laminar flow data with slight underprediction for the points close to transitional flow. In the graph, predictions in the transition region are also shown and the comparison with the data shows that probably the transition occurs earlier than is predicted by current methodology.

[FIGURE 4 OMITTED]

[FIGURE 5 OMITTED]

[FIGURE 6 OMITTED]

Langlinais et al. (1983) presented an extensive range of pressure drop measurements for single phase flow of different drilling fluids in an 1829 m vertical well. Predictions of this work with their results for the full range of flow regimes, laminar, transitional and turbulent flow is given in Figure 7. Very good match with the data at all flow types is seen. Similarly, the correct prediction of the points for laminar-transition and transition-turbulent flow regimes is observed.

Hanks (1979) provided the full solution to laminar flow of Herschel-Bulkley fluids in concentric annuli and presented his results in tabular and graphical form. In Figure 8, the predictions from this work are compared to the solution of Hanks, in the usual [DELTA]p-Q manner, where the curves of Hanks are taken from graphs of the original paper which were provided in dimensionless forms and for a range of annuli diameter ratio. One observes very close match of the predictions of this work with the predictions of the full numerical solution of Hanks (1979) for all aspect ratios, even down to an aspect ratio of 0.1.

[FIGURE 7 OMITTED]

[FIGURE 8 OMITTED]

Fordham et al. (1991) presented experimental results and compared with the predictions of their numerical solution for the flow of a 0.5% xanthan gum in water, in a 3 m long, 4 cm by 5 cm annuli for a range of flow rates. Comparison of the predictions of this work with their experimental data as well as their predictions (which matched their data) is shown in Figure 9. The data cover laminar flow only. There is excellent match of model predictions both with the experimental data and the predictions of the simulator of Fordham et al. (1991) which uses the full laminar flow solution in an annulus.

[FIGURE 9 OMITTED]

CONCLUSIONS

A consistent model which predicts the pressure drop for laminar, transitional and turbulent flow of Herschel-Bulkley fluids in concentric annuli, modelled as a slot, has been presented. For laminar flow, the model provides analytical solution for the pressure drop giving also full velocity profiles. For transitional and turbulent flows, semi-analytical solutions have been provided using local power law approach giving equations which relate the local power parameters, K', n', to the rheological parameters of Herschel-Bulkley model and to the annulus flow geometry. A transition region has also been defined which depends on (n').

The predictions for laminar flow have been compared to literature data containing experimental data and predictions of numerical simulators which model fully flow in annuli. The comparison ranges from good to very good. The predictions for transition and turbulent flow compare favourably well with the scarce data from literature for drilling fluids flowing in a model 1829 m well. More experimental data is needed which will span the transition and turbulent flows to further test the transition equations.

ACKNOWLEDGEMENT

The paper is dedicated to the memory of the first author, Mr. Konstantinos Founargiotakis, who tragically has left us in December 2006.

NOMENCLATURE A flow area ([m.sup.2]) [d.sub.1] diameter of inner tube of annulus (m) [d.sub.2] diameter of outer tube of annulus (m) d[p.sub.f]/dL pressure drop (Pa/m) f friction factor [f.sub.l1] friction factor at [Re.sub.1] [f.sub.t2] friction factor at [Re.sub.2] h slot height (m) K flow consistency index (Pa [s.sup.n]) K' flow consistency index for local power law parameters (Pa [s.sup.n]) L length (m) m inverse of flow behaviour index (=1/n) n flow behaviour index n' flow behaviour index for local power law parameters q flow rate ([m.sup.3]/s) [q.sub.given] given value of flow rate ([m.sup.3]/s) [r.sub.2] outer cylinder radius of annulus (m) [r.sub.1] inner cylinder radius of annulus (m) [Re.sub.1] Reynolds number for beginning transition laminar to transitional flow [Re.sub.2] Reynolds number for beginning transition transitional to turbulent flow [Re.sub.MRa] Metzner-Reed Reynolds number for annulus u velocity (m/s) V mean velocity (m/s) w slot width (m) y distance from bottom wall of slot (m) Greek Symbols [??] shear rate ([s.sup.-1]) [[??].sub.Nw] Newtonian shear rate on the wall ([s.sup.-1]) [[??].sub.w] wall shear rate ([s.sup.-1]) [[mu].sub.e] effective viscosity, (Pa s) [xi] dimensionless shear stress for annulus [rho] fluid density (kg/[m.sup.3]) [tau] shear stress (Pa) [[tau].sub.y] yield stress for Herschel-Bulkley fluid (Pa) [[tau].sub.w] wall shear stress (Pa)

Manuscript received November 17, 2006; revised manuscript received October 18, 2007; accepted for publication December 16, 2007.

REFERENCES

Bailey, W. J. and J. M. Peden, "A Generalized and Consistent Pressure Drop and Flow Regime Transition Model for Drilling Hydraulics," SPE Drill. Completion 15, 44-56 (2000).

Becker, T. E., R. G. Morgan, W. C. Chin and J. E. Griffith, "Improved Rheology Model and Hydraulic Analysis for Tomorrow's Wellbore Fluid Applications," Paper SPE 84215 presented at the SPE Productions and Operations Symposium, Oklahoma City, OK (2003).

Bingham, E. C., "Fluidity and Plasticity," McGraw-Hill, New York (1922).

Bird, R. B., G. C. Dai and B. J. Yarusso, "The Rheology and Flow of Viscoplastic Materials," Rev. Chem. Eng. 1, 1-70 (1983).

Bourgoyne, A. T., M. E. Chenevert, K. K. Millheim and F. S. Young, Jr., "Applied Drilling Engineering," SPE Textbook Series, Richardson TX (1991).

Buchtelova, M., Comments on "The Axial Laminar Flow of Yield-Pseudoplastic Fluids in a Concentric Annulus," Ind. Eng. Chem. Res. 27, 1557-1558 (1988).

Chin, W. C., "Computational Rheology for Pipeline and Annular Flow," Gulf. Pub. Co., Boston (2001).

David, J. and P. Filip, "Relationship of Annular and Parallel-Plate Poiseuille Flows for Power Law Fluids," Polym. Plast. Technol. Eng. J. 34, 947-960 (1995).

Dodge, D. W. and A. B. Metzner, "Turbulent Flow in Non-Newtonian Systems," AIChE J. 5, 189-204 (1959).

Escudier, M. P., P. J. Oliveira and F. T. Pinho, "Fully Developed Laminar Flow of Purely Viscous Non-Newtonian Liquids Through Annuli, Including the Effects of Eccentricity and Inner-Cylinder Rotation," Int. J. Heat Fluid Flow 23, 52-73 (2002a).

Escudier, M. P., P. J. Oliveira and F. T. Pinho, "Fully Developed Laminar Flow of Non-Newtonian Liquids Through Annuli: Comparison of Numerical Calculations with Experiments," Exp. Fluids 33, 101-111 (2002b).

Fang, P., R. M. Manglik and M. A. Jog, "Characteristics of Laminar Viscous Shear-Thinning Fluid Flows in Eccentric Annular Channels," J. Non-Newtonian Fluid Mech. 84, 1-17 (1999).

Filip, P. and J. David, "Axial Couette-Poiseuille Flow of Power-Law Viscoplastic Fluids in Concentric Annuli," J. Petr. Sci. Eng. 40, 111-119 (2003).

Fordham, E. J., S. H. Bittleston and M. A. Tehrani, "Viscoplastic Flow in Centered Annuli, Pipes and Slots," IEC Res. 29, 517-524 (1991).

Fredrickson, A. G. and R. B. Bird, "Flow of Non-Newtonian Fluids in Annuli," Ind. Eng. Chem. 50, 347-352 (1958).

Govier, G. W. and K. Aziz, "The Flow of Complex Mixtures in Pipes," Krieger, Malabar, FL (1972).

Grinchik, I. P. and A. K. Kim, "Axial Flow of a Non-Linear Viscoplastic Fluid Through Cylindrical Pipes," J. Eng. Phys. 23, 1039-1041 (1974).

Gucuyener, H. I. and T. Mehmetoglu, "Flow of Yield-Pseudoplastic Fluids Through a Concentric Annulus," AIChE J. 38, 1139-1143 (1992).

Guillot, D., "Rheology of Well Cementing Slurries," in E.B. Nelson, Ed., "Well Cementing," Schlumberger Educational Services, Houston (1990) pp. 4.01-37.

Guillot, D. and J. D. Dennis, "Prediction of Laminar and Turbulent Friction Pressures of Cement Slurries in Pipes and Centered Annuli," Paper SPE 18377 presented at the European Petroleum Conference, London, 18-19 October (1988).

Haciislamoglu, N. and J. Langlinais, "Non-Newtonian Flow in Eccentric Annuli," J. Energy Res. Tech. 112, 163-169 (1990).

Hanks, R. W., "The Axial Flow of Yield--Pseudoplastic Fluids in a Concentric Annulus," Ind. Eng. Chem. Process Des. Dev. 18, 488-493 (1979).

Hansen, S. A., R. Rommetveit, N. Sterri, B. Aas and A. Merlo, "A New Hydraulic Model for Slim Hole Drilling Applications," Paper SPE 57579 presented at the SPE/IADC Middle East Drilling Technology Conference, Abu Dhabi, 8-10 November (1999).

Harnett, J. P. and M. Kostic, "Turbulent Friction Factor Correlations for Power Law Fluids in Circular and Non-Circular Channel," Int. Commun. Heat Mass Trans. 17, 59-65 (1990).

Hemphil, T., W. Campos and M. A. Tehrani, "Yield Power Law Model Mode Accurately Predicts Mud Rheology," Oil Gas J. 91, 45-50 (1993).

Herschel, W. H. and R. Bulkley, "Konsistenzmessungen von Gummi-Benzollosungen," Kolloid Z. 39, 291-300 (1926).

Heywood, N. I. and D. Ch. Cheng, "Comparison of Methods for Predicting Head Loss in Turbulent Pipe Flow of Non-Newtonian Fluids," Trans. Inst. M.C. 6, 33-45 (1984).

Hussain, Q. E. and M. A. R. Sharif, "Viscoplastic Fluid Flow in Irregular Eccentric Annuli due to Axial Motion of the Inner Pipe," Can. J. Chem. Eng. 75, 1038-1045 (1997).

Hussain, Q. E. and M. A. R. Sharif, "Numerical Modeling of Helical Flow of Viscoplastic Fluids in Eccentric Annuli," AIChE J. 46, 1937-1946 (2000).

Kelessidis, V. C., A. Mihalakis and C. Tsamantaki, "Rheology and Rheological Parameter Determination of Bentonite-Water and Bentonite-Lignite-Water Mixtures at Low and High Temperatures," Proceedings of the 7th World Congress of Chem. Engr., Glasgow (2005).

Kelessidis, V. C., R. Maglione, C. Tsamantaki and Y. Aspirtakis, "Optimal Determination of Rheological Parameters for Herschel-Bulkley Drilling Fluids and Impact on Pressure Drop, Velocity Profiles and Penetration Rates during Drilling," J. Petrol. Sci. Eng. 53, 203-224 (2006).

Kelessidis, V. C., G. Christidis, P. Makri, V. Hadjistamou, C. Tsamantaki, A. Mihalakis, C. Papanikolaou and A. Foscolos, "Gelation of Water-Bentonite Suspensions at High Temperatures and Rheological Control with Lignite Addition," Appl. Clay Sci. 36, 221-231 (2007).

Langlinais, J. P., A. T. Bourgoyne and W. R. Holden, "Frictional Pressure Losses for the Flow of Drilling Mud and Mud/Gas Mixtures," Paper SPE 11993 presented at the 58th Annual Technical Conference and Exhibition, San Francisco, CA, 5-8 October (1983).

Maglione, R. and G. Ferrario, "Equations Determine Flow States for Yield-Pseudoplastic Drilling Fluids," Oil Gas J. 94, 63-66 (1996).

Maglione, R., A. Guarneri and G. Ferrari, "Rheologic and Hydraulic Parameter Integration Improves Drilling Operations," Oil Gas J. 97, 44-48 (1999).

Maglione, R., G. Robotti and R. Romagnoli, "In-Situ Rheological Characterization of Drilling Mud," SPE J. 5, 377-386 (2000).

Merlo, A., R. Maglione and C. Piatti, "An Innovative Model for Drilling Fluid Hydraulics," Paper SPE 29259 presented at the Asia Pacific Oil & Gas Conf., Kuala-Lumpur, Malaysia (1995).

Metzner, A. B. and J. C. Reed, "Flow of Non-Newtonian Fluids--Correlations of the Laminar, Transition and Turbulent-Flow Regions," AIChE J. 1, 434-440 (1955).

Meuric, O. F. J., R. J. Wakeman, T. W. Chiu and K. A. Fisher, "Numerical Flow Simulation of Viscoplastic Fluids in Annuli," Can. J. Chem. Eng. 76, 27-40 (1998).

Nguyen, Q. D. and D. V. Boger, "Measuring the Flow Properties of Yield Stress Liquids," Ann. Rev. Fluid Mech. 24, 47-88 (1987).

Okafor, M. N. and J. F. Evers, "Experimental Comparison of Rheology Models for Drilling Fluids," Paper SPE 24086 presented at the Western Regional Meeting, Bakersfield, OK, 30 March-11 April (1992).

Rabinowitch, B., "On the Viscosity and Elasticity of Sols," Z. Phisik. Chem. 1, 1-26 (1929).

Reed, T. D. and A. A. Pilehvari, "A New Model for Laminar, Transitional and Turbulent Flow of Drilling Muds," Paper SPE 25456 presented at the Productions Operations Symposium, Oklahoma City, OK, 21-23 March (1993).

Round, G. F. and S. Yu, "Entrance Laminar Flows of Viscoplastic Fluids in Concentric Annuli," Can. J. Chem. Eng. 71, 642-645 (1993).

Siginer, D. A. and S. L. Bakhtiyarov, "Flow of Drilling Fluids in Eccentric Annuli," J. Non-Newtonian Fluid Mech. 78, 119-132 (1998).

Skelland, A. H. P., "Non-Newtonian Flow and Heat Transfer," John Wiley, New York (1967).

Subramanian, R. and J. J. Azar, "Experimental Study on Friction Pressure Drop for Non-Newtonian Drilling Fluids in Pipes," Paper SPE 64647 presented at the SPE International Conference and Exhibition, Beijing, China, 7-10 November (2000).

Szabo, P. and O. Hassager, "Flow of Viscoplastic Fluids in Eccentric Annular Geometries," J. Non-Newton. Fluid Mech. 45, 149-169 (1992).

Walton, I. C. and S. H. Bittleston, "The Axial-Flow of a Bingham Plastic in a Narrow Eccentric Annulus," J. Fluid Mech. 222, 39-60 (1991).

Yang, L. and G. A. Chukwu, "A Simplified Couette-Flow Solution of Non-Newtonian Power-Law Fluids in Eccentric Annuli," Can. J. Chem. Eng. 73, 241-247 (1995).

K. Founargiotakis, (1) V. C. Kelessidis (1) * and R. Maglione (2)

(1.) Department of Mineral Resources Engineering, Technical University of Crete, Polytechnic City, 73100 Chania, Greece

(2.) Consultant, Vercelli, Italy

* Author to whom correspondence may be addressed. E-mail address: kelesidi@mred.tuc.gr

Printer friendly Cite/link Email Feedback | |

Author: | Founargiotakis, K.; Kelessidis, V.C.; Maglione, R. |
---|---|

Publication: | Canadian Journal of Chemical Engineering |

Date: | Aug 1, 2008 |

Words: | 4858 |

Previous Article: | Horizontal air--water flow in a square cross-section channel: pseudo-slug flow regime. |

Next Article: | Experimental and numerical study of heat transfer in horizontal concentric annulus containing phase change material. |