energies
Review A Review of Experiments and Modeling of Gas-Liquid Flow in Electrical Submersible Pumps
Jianjun Zhu * and Hong-Quan Zhang
McDougall School of Petroleum Engineering, The University of Tulsa, Tulsa, OK 74104, USA; [email protected] * Correspondence: [email protected]; Tel.: +1-918-760-6666
Received: 18 October 2017; Accepted: 3 January 2018; Published: 11 January 2018
Abstract: As the second most widely used artificial lift method in petroleum production (and first in produced amount), electrical submersible pump (ESP) maintains or increases flow rate by converting kinetic energy to hydraulic pressure of hydrocarbon fluids. To facilitate its optimal working conditions, an ESP has to be operated within a narrow application window. Issues like gas involvement, changing production rate and high oil viscosity, greatly impede ESP boosting pressure. Previous experimental studies showed that the presence of gas would cause ESP hydraulic head degradation. The flow behaviors inside ESPs under gassy conditions, such as pressure surging and gas pockets, further deteriorate ESP pressure boosting ability. Therefore, it is important to know what parameters govern the gas-liquid flow structure inside a rotating ESP and how it can be modeled. This paper presents a comprehensive review on the key factors that affect ESP performance under gassy flow conditions. Furthermore, the empirical and mechanistic models for predicting ESP pressure increment are discussed. The computational fluid dynamics (CFD)-based modeling approach for studying the multiphase flow in a rotating ESP is explained as well. The closure relationships that are critical to both mechanistic and numerical models are reviewed, which are helpful for further development of more accurate models for predicting ESP gas-liquid flow behaviors.
Keywords: electrical submersible pump (ESP); multiphase flow; mechanistic modeling; CFD; flow pattern; viscosity effect; drag coefficient; bubble size; in-situ gas void fraction
1. Introduction The electrical submersible pump (ESP), as a high-efficiency downhole equipment for converting kinetic energy to hydraulic pressure head, has improved significantly since it was invented in the 1910’s by the Russian engineer Arutunoff. Globally, there are more than 100,000 ESP installations, making ESP the second most widely used artificial lift method in oil production, but the first based on the produced amount [1,2]. ESPs excel in producing crude oils at very high flow rate, but they have to be operated within a narrow application window. Gas involvement, changing production rate and high oil viscosity can greatly affect the ESP performance [3,4]. Previous studies showed that the presence of gas would cause ESP hydraulic head degradation. The flow behaviors inside ESPs under gassy conditions, such as pressure surging and gas pockets, further deteriorate ESP boosting pressure. Surging may result in vibrations and short service life, and gas pockets can severely limit liquid production rates [5]. Although handling gas-liquid mixture has gradually become common for ESPs, the physical mechanism of two-phase flow affecting ESP hydraulic performance is not well understood. The gas bubble formation, coalescence and breakup mechanisms inside ESPs, which affect the two-phase flow characteristics, are still unclear. Due to the compact and complex geometries of multistage ESPs, the visualization of internal flow structures and bubble movement is very difficult.
Energies 2018, 11, 180; doi:10.3390/en11010180 www.mdpi.com/journal/energies Energies 2018, 11, 180 2 of 41 compactEnergies 2018and, complex11, 180 geometries of multistage ESPs, the visualization of internal flow structures and2 of 41 bubble movement is very difficult. The main objective of this paper is to provide a comprehensive review of literature concerning experimentalThe main studies objective and of modeling this paper approaches is to provide fo ar comprehensivepredicting ESP review pressure of literature increment concerning under gas-liquidexperimental flow studies conditions. and modeling The review approaches helps better for predicting understand ESP the pressure multiphase increment flow undercharacteristics gas-liquid insideflow a conditions. rotating ESP, The in order review to helpsdevelop better mechanis understandtic model the to multiphase predict ESP flow boosting characteristics pressure under inside gassya rotating flow conditions. ESP, in order Section to develop 2 presents mechanistic the fundamental model to predictconcepts ESP and boosting definitions pressure of a undermultistage gassy centrifugalflow conditions. pump. Section Section2 presents3 discusses the fundamentalthe previous conceptsknowledge and of definitions ESP hydraulic of a multistage performance centrifugal and flowpump. structures Section 3under discusses gassy the conditions previous knowledgeobserved from of ESP experiments. hydraulic performance Section 4 is anddivided flow into structures four subsections.under gassy It conditionsbegins with observed the empirical from experiments.correlations, followed Section4 byis divided one-dimensional into four subsections.two-fluid modeling It begins approach.with the empiricalThen, the correlations, CFD-based followed(Computational by one-dimensional Fluid Dynami two-fluidcs) modeling modeling and approach. mechanistic Then, approachesthe CFD-based are discussed. (Computational Finally, Fluid the Dynamics)closure rela modelingtionships andto make mechanistic models approaches solvable are are evaluated discussed. inFinally, Section the 5. closure relationships to make models solvable are evaluated in Section5.
2.2. Boosting Boosting Pressure Pressure of of Centrifugal Centrifugal Pump Pump ESPsESPs are are widely widely used used in in petroleum petroleum industry industry to toincrease increase the the hydrocarbon hydrocarbon fluid fluid production production rates, rates, especiallyespecially for for off-shore off-shore deep-water deep-water oil oil fields. fields. Normally, Normally, multistage multistage ESP ESP is isassembled assembled in inseries series (Figure (Figure 1a)1a) withwith each each stage stage comprising comprising the the rotating rotating impeller andand stationarystationary diffuserdiffuser (Figure (Figure1b,c). 1b,c). In In oil oil fields, fields, the theESP ESP system system usually usually contains contains hundreds hundreds of stagesof stages tomeet to meet the the boosting boosting pressure pressure requirement requirement due todue the togreat the great depth depth of reservoir. of reservoir. The The impeller impeller is the is dominantthe dominant part part of ESP of ESP since sinc ite forces it forces the the fluid fluid flow flow and andadds adds kinetic kinetic energy energy to to the the fluids fluids by by spinning spinning the th blades.e blades. At At the the outlet outlet of of impeller, impeller, the the main main part part of ofkinetic kinetic energyenergy inin thethe fluidfluid is converted to pressu pressurere potential potential by by diffuser diffuser vanes, vanes, which which are are fixed fixed and and workwork as as the the guiding guiding channels channels for for the the impeller impeller of of the the next next stage. stage.
(b)
(c)
(a)
Figure 1. Main components of an electrical submersible pump (ESP). (a) Multistage assembly (figure Figure 1. Main components of an electrical submersible pump (ESP). (a) Multistage assembly (figure courtesy of http://www.franklinwater.com); (b) impeller components; (c) diffuser components [3]. courtesy of http://www.franklinwater.com); (b) impeller components; (c) diffuser components [3].
Energies 2018, 11, 180 3 of 41 Energies 2018, 11, 180 3 of 41 As a type of centrifugal pumps, ESPs can be classified into three different categories: radial, axial and mixed types based on the dimensionless specific speed (NS). The non-dimensional NS is given by: As a type of centrifugal pumps, ESPs can be classified into three different categories: radial, axial 1 2 and mixed types based on the dimensionless specificQ speed (N ). The non-dimensional N is given by: N S S (1) S gH 3 4 ΩQ1/2 NS = 3 (1) where Ω is the rotational speed (rad/s), Q is liquid (flowgH) 3/4rate (m /s), g is local gravitational acceleration (m/s2), H is pump head (m). The pump industry uses a more practical expression for the specific 3 speedwhere NΩS isas the below: rotational speed (rad/s), Q is liquid flow rate (m /s), g is local gravitational acceleration (m/s2), H is pump head (m). The pump industry uses a more practical expression for the specific N q speed NS as below: N √ (2) q Nh0.75q N = (2) q h0.75 where N, q and h are rotational speed (rpm), flow rate (gpm) and pump head (ft), respectively. Based N q h onwhere Equation, and (2), arecentrifugal rotational pumps speed (rpm),are categorized flow rate (gpm)into radial, and pump mixed head and (ft), axial respectively. types. The Based radial on Equation (2), centrifugal pumps are categorized into radial, mixed and axial types. The radial pumps pumps usually fall in the range 500 < NS < 1800, while the mixed pumps can reach a maximum usually fall in the range 500 < NS < 1800, while the mixed pumps can reach a maximum NS = 4500. NS = 4500. ForFor a centrifugal pump, there are three important variables to characterize its hydraulic and H mechanicalmechanical performance, performance, namely: namely: pumping pumping head head ( (H),), efficiency efficiency ( (ηη)) and and brake brake horsepower horsepower (BHP). (BHP). FigureFigure 22 belowbelow showsshows thethe typical pump performanceperformance curves provided inin the product brochure, whichwhich areare obtained experimentally using tap-water.
Figure 2. Typical ESP pump performance curves (figure courtesy of Wood Group ESP, Inc., Oklahoma Figure 2. Typical ESP pump performance curves (figure courtesy of Wood Group ESP,Inc., Oklahoma City, City, OK, USA). OK, USA). With the liquid flow rate increase, the ESP performance curves exhibit different trends. The horsepowerWith the and liquid pumping flow head rate render increase, a monotonic the ESP in performancecreasing and decreasing curves exhibit trend versus different flow trends. rate, respectively,The horsepower whereas and the pumping efficiency head is in render a semi-q a monotonicuadratic relationship increasing andwith decreasingthe pump capacity. trend versus The bestflow efficiency rate, respectively, point (BEP) whereas presented the efficiencywith the dash is ined a line, semi-quadratic corresponds relationship to the highest with efficiency the pump of 68.9%capacity. at Q The = 2700 best bpd, efficiency N = 3500 point rpm (BEP) in Figure presented 2. It withis an theimportant dashed parameter line, corresponds in characterizing to the highest the ESPefficiency overall of performance. 68.9% at Q =In 2700this study, bpd, N the= main 3500 ob rpmjective in Figure is to review2. It is different an important models parameter in predicting in ESPcharacterizing hydraulic boosting the ESP overallpressure performance. under gas-liquid In this flow study, conditions. the main objective is to review different models in predicting ESP hydraulic boosting pressure under gas-liquid flow conditions. 2.1. Euler Head 2.1. Euler Head For an ESP impeller, the velocity triangles at the intake and discharge are shown in Figure 3. All For an ESP impeller, the velocity triangles at the intake and discharge are shown in Figure3. the variables are in International System of units (SI). The absolute velocity C can be decoupled into All the variables are in International System of units (SI). The absolute velocity C can be decoupled into two components: relative velocity W and peripheral velocity U, which is calculated by =Ω × .
Energies 2018, 11, 180 4 of 41
→ → → twoEnergies components: 2018, 11, 180 relative velocity W and peripheral velocity U, which is calculated by U = Ω4× of r41. The subscripts 1 and 2 denote the inlet and outlet of impeller, respectively. Here, W is relative to the The subscripts 1 and 2 denote the inlet and outlet of impeller, respectively. Here, W is relative to the impeller and C is relative to the global coordinate and is equal to the vector summation of U and W, impeller→ and→ C→ is relative to the global coordinate and is equal to the vector summation of U and W, e.g., C = W + U. In Figure3c, C’2 is an ideal absolute velocity assuming infinite number of impeller e.g., = + . In Figure 3c, C’2 is an ideal absolute velocity assuming infinite number of impeller blades, while C2 is the real absolute velocity at the impeller outlet. blades, while C2 is the real absolute velocity at the impeller outlet.
(a) (b)(c)
Figure 3. Velocity triangles. (a) Schematic of impeller flow channel; (b) inlet velocity triangle; (c) outlet Figure 3. Velocity triangles. (a) Schematic of impeller flow channel; (b) inlet velocity triangle; (c) outlet velocity triangle. velocity triangle.
According to the ideal conservation law of angular momentum in rotating centrifugal pump, According to the ideal conservation law of angular momentum in rotating centrifugal pump, the Euler head (HE) can be expressed as [6]: the Euler head (HE) can be expressed as [6]: C U C U U C U C H → 2 → 2 → 1 → 1 2 2U 1 1U (3) E C · U − C · U U C − U C H = 2 2 g 1 1 = 2 2U g 1 1U (3) E g g where C1U and C2U are the projection of absolute velocities at the impeller inlet and outlet to the wheredirectionC1 Uof andperipheralC2U are velocities. the projection Applying of absolutethe velocity velocities triangle atrelationship the impeller in Figure inlet and 3, one outlet can obtain: to the direction of peripheral velocities. Applying the velocity triangle relationship in Figure3, one can obtain: U 2 U 2 W 2 W 2 C 2 C 2 H 2 1 1 2 2 1 (4) E U2 −2gU2 W2 2−gW2 C2 −2gC2 H = 2 1 + 1 2 + 2 1 (4) E 2g 2g 2g and: and: 2 2 2 2 22 Ω rr2−rr1 QQΩ 11 11 HH = 2 1 − − (5)(5) EE g 2πgh tan β tan β g 2gh tan 2 2 tan 1 1 where r is the radius of impeller, h is the channel height, β is the blade angle. If the fluids enter the where r is the radius of impeller, h is the channel height, β is the blade angle. If the fluids enter the impeller without pre-rotation, Equation (5) can be written as: impeller without pre-rotation, Equation (5) can be written as: 2 2 Ω 2r22 QΩ HE = r2 − Q (6) HE g 2πgh tan β2 (6) g 2ghtan2 2.2. Head Loss Mechanisms
2.2. HeadAs shown Loss Mechanisms in Equation (6), the ideal Euler head HE is in linear relationship with liquid flow rate. In reality, pressure losses in the impeller, the diffuser, and losses from the interaction between As shown in Equation (6), the ideal Euler head HE is in linear relationship with liquid flow rate. them must be deducted from the ideal Euler head, including friction, shock and recirculation losses. In reality, pressure losses in the impeller, the diffuser, and losses from the interaction between them Therefore, the actual pump head can be calculated by: must be deducted from the ideal Euler head, including friction, shock and recirculation losses. Therefore, the actual pump head can be calculated by: H = HE − Hf rition − Hshock − Hleakage − Hrecirculation − Hdi f f user − Hdisk (7)
H H E H frition H shock H leakage H recirculation H diffuser H disk (7)
Figure 4a shows the change of the ideal Euler head with the outlet blade angle β. Figure 4b shows the losses to the ideal head and the final H-Q curve of a centrifugal pump according to Equation (7). Friction loss becomes prominent at high flow rates. In contrast, leakage loss is more at relatively low
Energies 2018, 11, 180 5 of 41
Figure4a shows the change of the ideal Euler head with the outlet blade angle β. Figure4b shows the losses to the ideal head and the final H-Q curve of a centrifugal pump according to Equation (7). Energies 2018, 11, 180 5 of 41 Friction loss becomes prominent at high flow rates. In contrast, leakage loss is more at relatively lowflow flow rates. rates. Shock Shock loss loss takes takes place place when when the the liquid liquid flow flow rate rate differs differs from from the the designed designed flow flow rate. rate. TablesTables1–5 1–5 below below discuss discuss each each pressure pressure loss mechanismloss mechanism by presenting by presenting existing existing models models that are that available are inavailable literature. in literature.
(a) (b)
Figure 4. Schematic of head curves. (a) Euler heads with different outlet blade angles; (b) actual head Figure 4. Schematic of head curves. (a) Euler heads with different outlet blade angles; (b) actual head after deducting losses [7]. after deducting losses [7]. Table 1 summarizes the studies on the friction losses inside a centrifugal pump impeller. Takacs [3] pointedTable 1out summarizes that friction the losses studies progressively on the friction increase losses with inside liquid a centrifugal rate and are pump due to impeller. fluid friction Takacs in [ 3] pointedthe impeller. out that friction losses progressively increase with liquid rate and are due to fluid friction in the impeller. Table 1. Friction loss models in literature. Table 1. Friction loss models in literature. Reference Model Coefficient Ito [8], Jones [9], Churchill [10], − Reference Model Coefficient= ℎ = Shah [11], Sun [12] 8 Ito [8], Jones [9], Churchill [10], f Q2 Wiesner [13], Sun and Prado [14], ( γβω − )( r2 −+ r1 ) = h f riction = 2 2 3 fγβω Fγ Fβ Fω f Shah [11], Sun [12] ℎ = 8gDH π bmsin βm r1r2 b2—constant 8 Thin et al. [15] 2 Wiesner [13], Sun and Prado [14], (D2−D1)(W1+W 2) h f riction = b2 + b2—constant Thin et al. [15] 8gsinβ2rH Ito and Nanbu [16], Bing et al. [17] ℎ = W2+W4 2 Ito and Nanbu [16], Bing et al. [17] h = N f s 1 2 f riction a γβω DH 4g 2 Zhu [1] ℎ = VI LI fFI—empirical constant Zhu [1] h f riction = fFI fFI—empirical constant 2gD2I
TheThe mechanism mechanism of of shock shock losses losses inin centrifugalcentrifugal pump is is still still not not well well studied. studied. Thus, Thus, only only empirical empirical correlations are available in the literature. Table 2 lists the existing prediction models of shock loss as correlations are available in the literature. Table2 lists the existing prediction models of shock loss as implemented in the prediction models. Shock losses are negligible at the best efficiency point (BEP) implemented in the prediction models. Shock losses are negligible at the best efficiency point (BEP) of of the pump, but increase at lower and higher liquid flow rates. They occur at the entrance and the the pump, but increase at lower and higher liquid flow rates. They occur at the entrance and the exit of exit of the impeller and are caused by sudden changes in the direction of flow [3]. the impeller and are caused by sudden changes in the direction of flow [3]. Table 2. Shock loss models in literature. Table 2. Shock loss models in literature. Reference Model Coefficient
StepanoffReference [6], Amaral [18] Model Coefficient ℎ = ( − ) kshock—empirical constant Thin et al. [15] Stepanoff [6], Amaral [18] 2 hshock = kshock(Q − QBEP) kshock—empirical constant ThinWiesner et al. [ 15[13]] − WiesnerSun [13] and Sun Prado and Prado [14] [ 14] ℎ = 2 h =2 kshock Q− Q BEP U ThinThin et et al. al. [15 [15]] shock 2g QBEP 1 V 2 ZhuZhu [1 [1]] ℎh == f I fTI—empiricalf —empirical constant constant TI TI2 2g TI
Energies 2018, 11, 180 6 of 41
As suggested by Tackas [3], the leakage losses represent the losses through the clearances between the rotating and stationary parts of the pump stage (at the impeller eye, through balancing holes, etc.), which decrease with increase of liquid rates. Table3 below shows the calculation models available in literature.
Table 3. Leakage loss models in literature.
Reference Model Coefficient
Qc1Uc1U2 Bing et al. [17] hleakage = 2Qg 2 2 3V V LG Zhu [1] h = L + f L fLK—empirical constant leakage 4g LK 2gSL
In a rotating centrifugal pump, recirculation loss is interpreted as the dissipated fluid energy which is continuously obtained from the shaft due to adverse pressure gradient from the impeller inlet to outlet. As liquid flow rate decreases, the adverse pressure gradient increases, resulting more severe head loss due to fluid recirculation. Currently, the agreement has not been reached as for the definition of recirculation, which in turn a fully understanding in terms of recirculation loss mechanism is unavailable by far. However, the empirical correlations have been proposed in literature, as summarized in Table4 below.
Table 4. Recirculation loss models in literature.
Reference Model Coefficient
Gulich [19,20] hrecirculation = 0 3 2 2 hrecirculation = krecsinh 3.5α2 D f U2 /g α = Q 2 arcsin 2+ 2− −5 Sun and Prado [14], Bing et al. [17] 2πr2b2(U2 W2 2U2W2cosβ2) k = 8 × 10 ! rec D = + W2 0.75gHt 1 − f 1 W 2 N D D 1 1t U2 ( a ) 1− 1t + 1t π D2 D2
hrecirculation = 0 Q > QBEP 2.5 Tuzson [21], Thin et al. [15] 3 2 Q h = krecω D Q ≤ Q recirculation 1 QBEP BEP 2 2 C2 −C2E Zhu [1] hrecirculation = 2g C2E—effective outlet velocity
Table 5. Diffuser loss models in literature.
Reference Model Coefficient
2 Fγ Fβ f Q r3di f −r2di f Ito [8], Jones [9], Churchill [10], Shah [11], Sun [12] hdi f f user = − 2 2 3 8gDH π bmsin βm r3di f r2di f (V2 +V2) Sun and Prado [14], Bing et al. [17] h = N F F f S 2d 3 di f f user a γ β DH 4g (V2 −V2) 2 2d 3 V3 − hdi f f user = − Cp (p3 p2) 2g 2g Cp = 2 Amaral [18] Q 0.5ρV3 V = 2d 2 2 V3 = V1 π r2di f f −r1di f f
Sun and Prado [12,14,22] claimed that the diffuser loss is caused mainly by the friction on the diffuser walls. Table5 presents the calculation formulae of diffuser losses that were implemented in the prediction models. Most equations are similar to the friction loss formulae in Table1. Disk friction losses, as Stepanoff [6] defined, are due to the contact between a rotating disk and fluid. Thus, additional power is consumed to keep the disk rotating since the viscous shear forces act on the disk surfaces [7]. Table6 shows some of existing models of calculating disk friction losses in literature. By deducting the pressure losses from Euler head in Equation (7), the real pressure increment of a centrifugal pump is obtained for single-phase flow. As can be seen in the Tables1–6, different pressure loss models are available in literature. A proper combination of these models may output an optimized Energies 2018, 11, 180 7 of 41 estimation of pump performance. Yet, due to the empirical nature of these models, the universal validity of such combination to deduce the most appropriate prediction model is questioned. Thus, more and more researchers resort to numerical approaches to calculate the hydraulic performance and understand the complicated flow structures inside rotating ESPs.
Table 6. Disk loss models in literature.
Reference Model Coefficient
Sun and Prado [14] (V02−V2 ) = 2 2d BEP Amaral [18] hdisk 2g f ρω3r5 Thin et al. [15] = disk 2 hdisk 109 Q 3 5 Cmω r2 Van Esch [23] hdisk = 2gQ cm—empirical constant k f ω3r5 5 h = RR th 2 1 − D1 disk gQ D2 Gulich [19,20] πr2 0.02 1+h/r2 = kRR = + 0.2 h 0.05r2 Ladouani [24] 2Re h Re 1+h/(2r2) 1.34 = − × −5 v fth exp 2 10 10−6
3. Experimental Studies The ESP boosting pressure depends on many hydraulic factors including pump geometry, fluid properties and multi-phase flow conditions. In recent years, with more and more installations of ESP in oil production systems, the effects of fluid viscosity and gas entrainment on ESP’s pressure boosting ability have been the focal areas of research interests.
3.1. Single-Phase Tests For single-phase flow, ESP performance is sensitive to the fluid viscosity. With the increase of oil viscosity, ESP boosting pressure becomes lower corresponding to the same flow rate. By the same token, the flow rate through the ESP decreases corresponding to the same boosting pressure in ESPs. Ippen [25] conducted over 200 performance tests for oil viscosities up to 10,000 Saybolt Second Universal (SSU) on four variants of centrifugal pumps. The experimental results were summarized by plotting the ratio of oil head to water head, brake horsepower and efficiency against a Reynolds-type dimensionless number, based on which the general correction factors for specific speeds from 800 to 2200 were proposed to correlate pump’s boosting pressure under viscous fluid flow. Hydraulic Institute provided a typical empirical approach with correction factors to estimate the conventional centrifugal pump boosting pressure for viscous liquid flow if the water performance were known. However, the accuracy of this approach was questioned by Gülich [19,20] and Li [26] since the experiments carried out by Hydraulic Institute were within a narrow range of the pump specific speeds. Unreasonable errors were found if extrapolation was beyond that range. Stepanoff [6] proposed a similar Reynolds-type number by using only one correction factor to get the new H-Q curves if the water performance were known. A more general model based on the evaluation of viscous dissipation for disk and hydraulic frictions to predict the boosting pressure of centrifugal pumps was proposed by Gülich. The friction losses on disk and in flow passage, as the author claimed, were dominating factors impairing centrifugal pump’s ability to handle high viscosity fluids. Compared with available data, Gülich also pointed out that friction losses were affected significantly by pump geometry, fluid properties and thermal conditions. More recent experimental studies conducted by Amaral et al. [27] and Solano [28] further revealed that the Hydraulic Institute charts and empirical correlations available in literature were unable to give appropriate correction factors to predict the ESP boosting pressure for viscous oil flow. The discrepancies between experimental measurements and predictions by existing correlations or charts are discussed in Barrios et al. [29] and Banjar et al. [30]. Energies 2018, 11, 180 8 of 41
3.2. Two-Phase Tests The pioneering experimental studies conducted by Murakami and Minemura [31,32] investigated centrifugal pump performance with gas entrainment. By employing a semi-opened impeller pump with a transparent casing, they experimentally observed the behavior of entrained air bubbles. The decreasing total head of the pump caused by air admission and the work consumed for air delivery were reported. Since then, investigators have conducted more experimental measurements and mathematical modeling on centrifugal pump performance under gassy flow conditions. Experimental researches of ESP performance under gas-liquid flow conditions have been carried out extensively at the Tulsa University Artificial Lift Project (TUALP) in the past decades. Table7 lists a summary of these studies.
Table 7. Experimental studies on ESP two-phase performance conducted at TUALP.
Study Content Pump Fluid GN4000 Cirilo [33] Compare two-phase flow performance of three different ESPs Air/water GN7000 ESP gas-liquid performance with an advanced gas handler Romero [34] GN4000 Air/water installed upstream Pessoa [35] Measure stage-by-stage pump pressure increment of a multistage ESP GC6100 Air/water Beltur [36] Investigate pressure surging in ESP and affecting factors GC6100 Air/water Duran [37] Correlate experimental data of ESP two-phase performance GC6100 Air/water Investigate pump rotational speed effect on ESP Zapata [38] GC6100 Air/water two-phase performance Visualize the internal flow of a 2nd stage ESP under gas/liquid Barrios [39] GC6100 Air/water flow conditions Visualize ESP two-phase flow pattern using a similar pump Gamboa [40] GC6100 Air/water prototype as Barrios [37] Air/oil Trevisan [41] Visualize ESP internal flow under air/viscous-liquid flow GC6100 visualization Banjar [30] Investigate ESP performance with air/oil flow DN1750 Air/oil Salehi [42] Investigate ESP gas/liquid performance with various flow conditions TE2700 Air/water Oil/water Croce [43] Investigate ESP performance with water/oil emulsion flow DN1750 Emulsion Investigate ESP gas/liquid flow performance with/without Air/water Zhu [1] TE2700 surfactant injections Surfactant
As can been seen in Table7, the previous experimental studies, which were conducted by Cirilo [33], Pessoa [35], Beltur [36], Duran [37], Zapata [38], Barrios [39], Gamboa [40], and Salehi [42] among others, covers a wide range of flow conditions. Cirilo built the experimental flow loops for testing ESP performance at the TUALP. Using water and air as the working fluids, he measured both the water and air-water performances of three different types of ESPs as a function of GVFs (the volumetric fraction of gas phase at the ESP inlet), intake pressure and rotational speeds. The obtained data indicated that the mixed type pumps were able to handle much higher GVFs (>30%) than radial type pumps (<10%). With necessary modifications of the same testing flow loops, Pessoa conducted experimental investigations of two-phase performance of ESP using a 22-stage GC6100 pump. By monitoring the stage-by-stage pressure increment, his results revealed that the ESP average behavior was significantly different from that observed per stage. Also, phenomena like surging and gas locking were observed and their boundaries were mapped. Additionally, a second region after pressure surging was observed in mapping test curves, where the slope of pressure increment versus flow rate changed again. Using the same experimental flow loop of Pessoa, Beltur, Duran, Zapata and Salehi conducted extensive experimental measurements of ESP performance under both liquid- and gas-liquid flow conditions. Beltur focused on ESP performance deterioration with the presence of gas for varying intake GVFs and pressure. Data analysis revealed that the intake GVF is the most important factor affecting ESP boosting pressure under gassy flow. A higher deterioration of pumping head occurs at GVFs EnergiesEnergies2018 2018, 11, ,11 180, 180 9 of9 41 of 41
intake GVFs and pressure. Data analysis revealed that the intake GVF is the most important factor aboveaffecting 6%. Duran ESP boosting and Zapata pressure developed under gassy empirical flow. A correlations higher deterioration for predicting of pumping the pressure head occurs increment at acrossGVFs the above stage 6%. and Duran the flow and regime Zapata boundaries. developed Zapata empirical also correlations carried out for further predicting measurements the pressure with a wideincrement range of across rotational the speedstage and to study the itsflow effects regime on theboundaries. average efficiencyZapata also of ESP.carried Salehi out did further similar measurementsmeasurements using with a a 14-stage wide range TE2700 of rotational ESP, on whichspeed to the study effects itsof effects the stage on the number, average intake efficiency pressure of andESP. fluid Salehi properties did similar were measurements investigated. It using was founda 14-stag thate TE2700 ESP two-phase ESP, on which boosting the pressure effects of varied the stage stage bynumber, stage only intake when pressure the GVF and exceeded fluid properties a certain were value, investigated. below which It was the found deterioration that ESP was two-phase mild and independentboosting pressure of stage varied number. stage However, by stage theonly degradation when the GVF of pumpingexceeded heada certain became value, more below prominent which andthe affected deterioration by the was stage mild number and independent if the GVF of reached stage number. a critical However, value, atthe which degradation the pump of pumping boosting pressurehead became of the downstreammore prominent stage and is betteraffected than by the upstream stage number ones. if the GVF reached a critical value, at which the pump boosting pressure of the downstream stage is better than upstream ones. A typical testing result of ESP two-phase performance with pure water-air flow is shown in A typical testing result of ESP two-phase performance with pure water-air flow is shown in Figure5[ 40], including surging (Figure5a) and mapping tests (Figure5b). For surging test, the liquid Figure 5 [40], including surging (Figure 5a) and mapping tests (Figure 5b). For surging test, the liquid flow rate is fixed, but gas flow rate or inlet gas volumetric fraction (GVF) varies. On the contrary, flow rate is fixed, but gas flow rate or inlet gas volumetric fraction (GVF) varies. On the contrary, the the liquid flow rate changes while the gas flow rate is fixed in mapping test. In Figure5a, different liquid flow rate changes while the gas flow rate is fixed in mapping test. In Figure 5a, different curves curvescorrespond correspond to different to different pump pumpintake intakepressures. pressures. Obviously, Obviously, the pump the pressure pump increment pressure incrementdeclines declinesnon-linearly non-linearly with the with increase the increase of inlet of inletGVF. GVF. For Formapping mapping test testresults results in Figure in Figure 5b,5 b,each each curve curve correspondscorresponds to to a differenta different gas gas flow flow rate, rate, which which increasesincreases from the the top top to to bottom bottom curves. curves. With With a fixed a fixed gasgas flow flow rate, rate, a suddena sudden drop drop of of pump pump stage stage pressurepressure increment is is obse observedrved when when flow flow rate rate is isreduced reduced to ato certain a certain value. value.
(a) (b)
Figure 5. ESP two-phase pressure increment [40]. (a) Surging test; (b) mapping test. Figure 5. ESP two-phase pressure increment [40]. (a) Surging test; (b) mapping test.
A more recent experimental study regarding surfactant effect on ESP gas-handling ability and boostingA more pressure recent experimentalunder gassy flow study was regarding conducted surfactant by Zhu et effectal. [44]. on Surfactants ESP gas-handling are molecules ability with and boostinga hydrophobic pressure and under a hydrophilic gassy flow part, was and conducted therefore bythey Zhu preferentially et al. [44]. adsorb Surfactants at the are interface molecules of withcontinuous/dispersed a hydrophobic and phases. a hydrophilic In the part,process, and they therefore reduce they the preferentiallyinterfacial tension adsorb [45]. at theStudies interface of of continuous/dispersedsurfactant effect on bubble/droplet phases. In formation the process, in air/water they reduce or water/oil the interfacial binary tensionimmiscible [45 ].two-phase Studies of surfactantflows were effect carried on bubble/droplet out by Hu et al. formation [46], Omer in and air/water Pal [47]. or Surfactant water/oil effect binary on liquid immiscible loading two-phase in gas flowswell were has carriedbeen studied out by by Hu van et al.Nimwegen [46], Omer et andal. [48,49], Pal [47 Ajani]. Surfactant et al. [50,51]. effect Both on liquid reported loading a lower in gas wellcritical has been gas studiedvelocity by at vanwhich Nimwegen liquid loading et al. [48 occurs,,49], Ajani meaning et al. [that50,51 the]. Both surfactant reported presence a lower helps critical gasprevent velocity liquid at which loading liquid in wells. loading occurs, meaning that the surfactant presence helps prevent liquid loadingHowever, in wells. the studies on the surfactant effect on centrifugal pump performance under two-phase flowHowever, conditions the are studies few. Zhu on the et al. surfactant [44] measured effect the on ESP centrifugal two-phase pump pressure performance increment under with/without two-phase flowsurfactant conditions injection are few. in the Zhu flow et al. loop. [44] The measured surfactant the ESPused two-phase in their study pressure is isopropanol increment alcohol with/without (IPA). surfactantThe decline injection trend inof the normalize flow loop. pressure The surfactant (Np) in Figure used in6a theirversus study the intake is isopropanol GVF improves alcohol with (IPA). surfactant injections, and the pressure surging phenomenon depicted by a sudden drop of ESP The decline trend of normalize pressure (Np) in Figure6a versus the intake GVF improves with boosting pressure disappears. Figure 6b compares ESP H-Q performance curves under air-water flow surfactant injections, and the pressure surging phenomenon depicted by a sudden drop of ESP with different IPA volumetric concentrations. A clear difference can be seen at low QL in terms of the boosting pressure disappears. Figure6b compares ESP H-Q performance curves under air-water flow stage pressure increment. Without IPA injection, the pressure increment drops to zero if QL becomes with different IPA volumetric concentrations. A clear difference can be seen at low Q in terms of the very low. But, the ESP performance improves significantly with surfactant presence. L stage pressure increment. Without IPA injection, the pressure increment drops to zero if QL becomes very low. But, the ESP performance improves significantly with surfactant presence. Energies 2018, 11, 180 10 of 41 Energies 2018, 11, 180 10 of 41
Energies 2018, 11, 180 10 of 41
(a) (b)
Figure 6. Comparison of ESP two-phase pressure increment with/without surfactant injection [44]. Figure 6. Comparison of ESP two-phase pressure increment with/without surfactant injection [44]. (a) Constant liquid flow(a rate) testing; (b) constant gas flow rate testing. (b) (a) Constant liquid flow rate testing; (b) constant gas flow rate testing. Figure 6. Comparison of ESP two-phase pressure increment with/without surfactant injection [44]. 3.3. Flow(a) PatternConstant Recognition liquid flow rate testing; (b) constant gas flow rate testing. 3.3. Flow Pattern Recognition In multiphase pipe flow, flow pattern varies due to phase interaction and interfacial momentum 3.3. Flow Pattern Recognition transferIn multiphase changes, pipewhich flow, depend flow on pattern fluid properties varies due, flow to phase rates interactionand pipe geometry. and interfacial Many studies momentum on transferthe recognition changes,In multiphase of which flow pipe dependpatterns flow, flow onand fluidpattern their properties, transition varies due boundaries flow to phase rates interaction andare available pipe and geometry. interfacialin literature Many momentum [52–54]. studies onFlow thetransfer recognition pattern changes, is important of which flow patterns dependto accurately on and fluid theirpredict properties transition the liquid, flow boundaries holdup rates andand are pipepressure available geometry. gradient in literatureMany in pipe studies [flow.52– on54 ]. FlowTakethe pattern recognitionthe upward is important ofgas/liquid flow topatterns accurately flow inand a vertical predicttheir transition pipe the liquid for boundariesexample, holdup the and are main pressureavailable flow patterns gradientin literature are in pipebubble, [52–54]. flow. Takeslug,Flow the churn upwardpattern and is gas/liquid annularimportant flows to flow accurately as in illustrated a vertical predict in pipe theFigure forliquid example,7a. holdup Correspondingly, theand main pressure flow the gradient patternsflow pattern in arepipe bubble,map flow. slug,showingTake churn the the andupward transition annular gas/liquid boundaries flows flow as illustrated inbetween a vertical different in pipe Figure for flow7 a.example, patterns Correspondingly, the is mainplotted flow in the Figurepatterns flow pattern7b are with bubble, the map showingsuperficialslug, thechurn transitionvelocities and annular of boundaries gas flowsand liquid as between illustrated as its coordinates. different in Figure flow 7a. patterns Correspondingly, is plotted the in Figureflow pattern7b with map the superficialshowingThe velocitiesdispersed the transition ofbubble gas boundaries and flow liquid prevails between as itsif the coordinates. gasdifferent flow rateflow is patterns low and is the plotted liquid inflow Figure rate 7bis high, with asthe shownsuperficialThe dispersedin Figure velocities 7a. bubble For of this gas flow flow and prevails pattern, liquid ifas thethe its turbul coordinates. gas flowence rate forces is loware high and enough the liquid to breakup flow rate gas-phase is high, as showninto small inThe Figure bubblesdispersed7a. For and bubble this disperse flow flow pattern, them prevails in thethe if turbulencethecontinuous gas flow forcesliquid-phase rate is are low high and [52]. enoughthe liquid to flow breakup rate gas-phaseis high, as intoshown smallWhen bubblesin Figureboth andgas 7a. flowFor disperse this rate flow themand pattern, liquid in the theflow continuous turbul rate enceare liquid-phaselow, forces gas are bubbles high [52 enough]. become to breakuplarger and gas-phase move upward in a zigzag path. The bubble flow regime may or may not exist depending on the pipe intoWhen small both bubbles gas flow and disperse rate and them liquid in the flow continuous rate arelow, liquid-phase gas bubbles [52]. become larger and move diameter.When Confined both gas by flowmaximum rate and lattice liquid packing, flow ratethe in-situare low, gas gas void bubbles fraction become of bubble larger flow and should move upward in a zigzag path. The bubble flow regime may or may not exist depending on the pipe beupward below 0.52in a [55]. zigzag However, path. The for abubble practical flow application, regime may the or maximum may not gasexist void depending fraction belowon the 0.25 pipe diameter. Confined by maximum lattice packing, the in-situ gas void fraction of bubble flow should be isdiameter. commonly Confined used for by vertical maximum co-current lattice pipe packing, flow. the in-situ gas void fraction of bubble flow should below 0.52 [55]. However, for a practical application, the maximum gas void fraction below 0.25 is be below 0.52 [55]. However, for a practical application, the maximum gas void fraction below 0.25 commonly used for vertical co-current pipe flow. is commonly used for vertical co-current pipe flow.
(a) (b)
Figure 7. Typical flow patterns and transition boundaries in vertical pipe, (a) main observed flow patterns [54]; (b) transition(a) boundaries and flow pattern map [52]. (b) Figure 7. Typical flow patterns and transition boundaries in vertical pipe, (a) main observed flow Figure 7. Typical flow patterns and transition boundaries in vertical pipe, (a) main observed flow patterns [54]; (b) transition boundaries and flow pattern map [52]. patterns [54]; (b) transition boundaries and flow pattern map [52].
Energies 2018, 11, 180 11 of 41
Slug flow regime corresponds to the steady alternating flow of gas pockets (Taylor bubble) and liquid slugs [53] at even higher gas flow rates (see Figure7b), where Taylor bubbles occupy the entire cross section of the pipe except for a thin liquid film on the wall. Slug flow can be considered as steady-state flow pattern since its slug length, frequency, translational velocity etc. are time-independent. However, a further increase of gas flow rates triggers the collapse of slugs and lead to an unstable flow that is described as the churn flow regime due to the large degree of turbulence in the flow [54]. Annular flow is characterized by extremely high gas flow rates and continuous gas core flowing upward with thin liquid film on the pipe wall. Due to high interfacial shear and drag force effects, the liquid droplets from the wavy interface between gas core and liquid film can be swept into the gas core, and pushed upward by the high-speed gas phase. Several visualization studies have been conducted to observe multiphase flow structures/patterns inside a rotating centrifugal pump. Table8 below lists the visualization experiments on flow patterns inside a centrifugal pump in literature. Most of these studies used transparent casing to visually observe the pump internal flow patterns. Starting with the original volute-type centrifugal pumps [31,32,56–70], they cover several research topics including flow pattern recognition, bubble movement visualization, mapping transition boundaries etc.
Table 8. Experimental studies on two-phase flow patterns inside a centrifugal pump.
Study Approach Observation Murakami and Minemura 1st study that associated pump perfor-mance with Transparent casing [31,32,56] gas-liquid flow pattern inside an impeller Small bubble flow and stationary large bubble coinciding Patel & Runstadler [57] Transparent casing with significant reduction of pump head Transparent casing The evolvement of dispersed bubble to slug flow and a large Sekoguchi et al. [58] electric resistivity probe gas pocket Kim et al. [59] Transparent casing Phase slippage leads to bubble agglomeration Sato et al. [60] Transparent casing Bubble flow to separated flow as gas flow rate increases Takemura et al. [61] Transparent casing Bubble accumulation coincided with large pressure gradient Gas filled space starts to develop at the suction pipe of Chisely [62] Transparent casing the pump Suryawijaya et al. [63] Transparent casing Bubble accumulates at pressure side of impeller blade Estevam [64] Transparent casing 1st study of flow patterns in ESP Poullikkas [65] Transparent casing Gas bubbles, gas pocket and blockage of pump flow Larger accumulation of bubbles is observed at pressure side Thum et al. [66] Transparent casing than suction side, four flow patterns: bubble, plug, slug and stratified flow Barrios [39,67] Transparent casing Bubbly flow and segregated flow are observed Gamboa [40] Transparent casing 1st study of flow pattern map in ESP impeller Measurement of gas void distribution in multiphase Schäfer et al. [68] HireCT, non-intrusive centrifugal pump Measurement of gas void distribution in multiphase Neumann et al. [69] HireCT, non-intrusive centrifugal pump Flow patterns in ESP: bubble, agglomerated bubble, gas Verde et al. [70] Transparent casing pocket, segregated flow
Murakami and Minemura [31,32] experimental study first investigated the gas entrainment effect on the hydraulic pump head of a centrifugal pump visually, and associated the degradation of pump performance with gas-liquid flow pattern inside the rotor. Using a transparent casing, Estevam [62] conducted the first visualization experiment of flow pattern recognition inside a rotating ESP impeller. Later visualization experimental work, conducted by Gamboa [40] presented similar observations Energies 2018, 11, 180 12 of 41 Energies 2018, 11, 180 12 of 41
Later visualization experimental work, conducted by Gamboa [40] presented similar observations as as Estevam’s study, based on which the first flow pattern map with transition boundaries between Estevam’s study, based on which the first flow pattern map with transition boundaries between different flow patterns inside an ESP was proposed. His work offers significant insight regarding how different flow patterns inside an ESP was proposed. His work offers significant insight regarding to characterize the flow patterns in an ESP quantitatively. how to characterize the flow patterns in an ESP quantitatively. Among previous visualization experiments, two recent studies from Schäfer et al. [68] and Among previous visualization experiments, two recent studies from Schäfer et al. [68] and Neumann et al. [69] are different from the others. They both used a new non-intrusive technique, the Neumann et al. [69] are different from the others. They both used a new non-intrusive technique, so called high resolution computed tomography (HireCT) to accurately measure the distribution of the so called high resolution computed tomography (HireCT) to accurately measure the distribution in-situ gas void fractions. Compared to the intrusive technique, the new measurement approach, as of in-situ gas void fractions. Compared to the intrusive technique, the new measurement approach, the authors claimed, can provide better image resolution of the gas void distribution and thus a higher as the authors claimed, can provide better image resolution of the gas void distribution and thus a measurement accuracy. higher measurement accuracy. By visualizing the ESP internal flow, Barrios [39], Barrios and Prado [67] observed that bubbles By visualizing the ESP internal flow, Barrios [39], Barrios and Prado [67] observed that bubbles enlarged when inlet GVFs increased and pump rotational speeds decreased. Such enlargement enlarged when inlet GVFs increased and pump rotational speeds decreased. Such enlargement corresponding to the poorer pump performance indicated that bubble behaviors played a significant corresponding to the poorer pump performance indicated that bubble behaviors played a significant role in ESP’s ability of handling gas-liquid mixtures. In addition, visualization experiments also role in ESP’s ability of handling gas-liquid mixtures. In addition, visualization experiments also showed different flow patterns prevailing inside ESP channels at higher GVFs. showed different flow patterns prevailing inside ESP channels at higher GVFs. In Figure8, λG denotes the intake GVF. It is evident that flow behaviors inside ESP impeller In Figure 8, λG denotes the intake GVF. It is evident that flow behaviors inside ESP impeller altered significantly with flow conditions. From Figure8a,b, the GVF increase caused formation of altered significantly with flow conditions. From Figure 8a,b, the GVF increase caused formation of larger bubbles and gas pockets, which in turn choked the flow passage for liquids and decreased the larger bubbles and gas pockets, which in turn choked the flow passage for liquids and decreased the pump hydraulic head. Moreover, the flow patterns prevailing in ESP impeller at a specific value of λG, pump hydraulic head. Moreover, the flow patterns prevailing in ESP impeller at a specific value of are comparable to that of gas-liquid two-phase flow in pipelines. λG, are comparable to that of gas-liquid two-phase flow in pipelines.
(a) (b)
(c) (d)
Figure 8. Impeller channel flow behavior at N = 600 rpm and QL = 174 bpd with various gas flow rates. Figure 8. Impeller channel flow behavior at N = 600 rpm and QL = 174 bpd with various gas flow rates. (a) λG = 0.15%; (b) λG = 0.23%; (c) λG = 0.39%; (d) λG = 1.05% [39]. (a) λG = 0.15%; (b) λG = 0.23%; (c) λG = 0.39%; (d) λG = 1.05% [39].
Energies 2018, 11, 180 13 of 41 Energies 2018, 11, 180 13 of 41
Similar as thethe gas/liquidgas/liquid pipe pipe flow, flow, Estevam [64] [64] and Estevam et al. [[71]71] categorizedcategorized the flowflow patternspatterns inin a a rotating rotating ESP ESP impeller impeller as bubblyas bubbly flow, fl transitionow, transition flow, andflow, elongated-bubble and elongated-bubble flow. A moreflow. thoroughA more thorough experimental experimental study on stud flowy on pattern flow pattern recognition recognition of two-phase of two-phase flow inflow ESP in impellerESP impeller was carriedwas carried out by out Gamboa by Gamboa [40], Gamboa[40], Gamboa and Pradoand Prado [72]. Under[72]. Under a given a given flow condition,flow condition, the typical the typical flow patternflow pattern mapfrom map their from visualization their visualization experiments experi isments shown is in shown Figure 9in, where Figure the 9, curveswhere denotedthe curves by differentdenoted by symbols different represent symbols the represent transition the boundaries transition between boundaries flow between patterns, flow including patterns, homogenous including flow,homogenous bubbly flow, flow, gas-pocket bubbly flow, formation, gas-pocket and formation, segregated and flow. segregated Figure9 is flow. important Figure to 9 understand is important the to two-phaseunderstand flow the two-phase mechanisms flow in amechanisms rotating ESP. in From a rotating the perspective ESP. From ofthe mechanistic perspective modeling, of mechanistic each flowmodeling, regime each corresponds flow regime to thecorresponds specific governing to the spec equationsific governing for flow equations characteristics, for flow suchcharacteristics, as bubble sizesuch ( dasb), bubble in-situ size gas ( voiddb), in-situ fraction gas (α voidG), slippage fraction velocity (αG), slippage (Vslip) velocity between ( gasVslip) and between liquid gas phases. and liquid With thephases. flow With pattern the prevailing flow pattern inside prevailing a rotating inside ESP a determined, rotating ESP the determined, governing equationsthe governing based equations on mass andbased momentum on mass and conservations momentum can cons beervations simplified. can be simplified.
1.2
1 Homogenous Boundary
0.8 Homogenous Bubbly Flow Regime Regime Unstable Gas Pocket 0.6
0.4 Stable Gas Pocket Regime 0.2 Normalized Liquid Rate[fraction] Segregated Flow Regime 0 0 0.05 0.1 0.15 0.2 Normalized Gas Rate [fraction]
Figure 9. Flow pattern map for GC-6100 ESP at 2400 rpm andand 150150 psigpsig [[39].39].
Verde et al. [70] conducted visualization experiments on flow pattern recognition inside a Verde et al. [70] conducted visualization experiments on flow pattern recognition inside a rotating rotating ESP impeller using high speed and resolution imaging technique. As shown in Figure 10, ESP impeller using high speed and resolution imaging technique. As shown in Figure 10, four four different flow patterns were classified, including bubble flow (Figure 10a), agglomerated bubble different flow patterns were classified, including bubble flow (Figure 10a), agglomerated bubble flow flow (Figure 10b), gas pocket flow (Figure 10c) and segregated flow (Figure 10d). They observed that (Figure 10b), gas pocket flow (Figure 10c) and segregated flow (Figure 10d). They observed that the intensity of pump performance degradation is directly influenced by the flow pattern within the the intensity of pump performance degradation is directly influenced by the flow pattern within impeller. The occurrence of the gas pocket flow pattern is linked to the intensification of the the impeller. The occurrence of the gas pocket flow pattern is linked to the intensification of deterioration of pump performance and the appearance of operation instabilities. Moreover, the the deterioration of pump performance and the appearance of operation instabilities. Moreover, segregated flow patterns correspond to the severe performance degradation which makes the pump the segregated flow patterns correspond to the severe performance degradation which makes the incapable of generating pressure. pump incapable of generating pressure. Summarized by Verde et al. [70], the schematic representations of each flow pattern are shown Summarized by Verde et al. [70], the schematic representations of each flow pattern are shown in in Figure 11 with intake GVF increasing from left to right. Due to the relatively small intake GVF, the Figure 11 with intake GVF increasing from left to right. Due to the relatively small intake GVF, the homogenous flow regime is featured by tiny and evenly-dispersed bubbles inside impeller channels, homogenous flow regime is featured by tiny and evenly-dispersed bubbles inside impeller channels, as shown in Figure 11a. Under this regime, the bubbles are deemed to move together with liquid as shown in Figure 11a. Under this regime, the bubbles are deemed to move together with liquid phase. Slippage between gas and liquid is small, meaning the in-situ αG is almost the same as λG. phase. Slippage between gas and liquid is small, meaning the in-situ α is almost the same as λ . As the intake GVF increases, the tiny bubbles are prone to collideG and aggregate to form biggerG As the intake GVF increases, the tiny bubbles are prone to collide and aggregate to form bigger ones, resulting in bubbly flow regime. In contrast to homogenous flow regime, the phase slippage ones, resulting in bubbly flow regime. In contrast to homogenous flow regime, the phase slippage between gas and liquid, shown by Figure 11b can no longer be neglected. Thus, depending on the intake GVF, the in-situ αG under bubbly flow becomes higher than λG. A further increase of GVF induces more severe collision and aggregation of bubbles so that the large gas pocket forms. This
Energies 2018, 11, 180 14 of 41 between gas and liquid, shown by Figure 11b can no longer be neglected. Thus, depending on the intakeEnergies GVF, 2018, the11, 180 in-situ αG under bubbly flow becomes higher than λG. A further increase14 of of GVF41 inducesEnergies more 2018 severe, 11, 180 collision and aggregation of bubbles so that the large gas pocket forms.14 This of 41 flow patternflow pattern is similar is similar as slug/churn as slug/churn flow flow patterns patterns in pipelines,in pipelines, which which are are featuredfeatured by a a gas gas core/pocket core/pocket flow pattern is similar as slug/churn flow patterns in pipelines, which are featured by a gas core/pocket followedfollowed by by a liquida liquid slug. slug. As As shown shown inin FigureFigure 1111c,c, the Taylor-bubble-like Taylor-bubble-like gas gas pocket pocket forms forms near near the the followed by a liquid slug. As shown in Figure 11c, the Taylor-bubble-like gas pocket forms near the suctionsuction side side of of ESP ESP impeller, impeller, which which occupies occupies aa significantsignificant portion of of the the impeller impeller flow flow passage. passage. suction side of ESP impeller, which occupies a significant portion of the impeller flow passage.
(a()a ) (b(b) )
(c) (d) (c) (d) Figure 10. Visualized flow patterns in a rotating ESP impeller [68]. (a) Bubble flow pattern; Figure 10. Visualized flow patterns in a rotating ESP impeller [68]. (a) Bubble flow pattern; Figure(b) agglomerated 10. Visualized bubble flow flow patterns pattern; in (c ) agas-pocket rotating flowESP pattern;impeller (d )[68]. segregated (a) Bubble flow pattern. flow pattern; (b) agglomerated bubble flow pattern; (c) gas-pocket flow pattern; (d) segregated flow pattern. (b) agglomerated bubble flow pattern; (c) gas-pocket flow pattern; (d) segregated flow pattern.
(a) (b)(c)(d)
Figure(a) 11. Schematic representations(b)( of flow patterns inside a rotatingc)( ESP impeller [70], (a)d bubble) flow; (b) agglomerated bubble flow; (c) gas-pocket flow; (d) segregated flow. Figure 11. Schematic representations of flow patterns inside a rotating ESP impeller [70], (a) bubble Figure 11. Schematic representations of flow patterns inside a rotating ESP impeller [70], (a) bubble flow;A (relativelyb) agglomerated high gas bubble flow flow; rate (andc) gas-pocket low liquid flow; flow (d) ratesegregated lead to flow. the segregated flow pattern, flow; (b) agglomerated bubble flow; (c) gas-pocket flow; (d) segregated flow. similar to the annular/stratified flows in pipe. As it can be seen in Figure 11d, the elongated bubble expandsA relatively and occupies high gas the flow entire rate impeller and low length. liquid For flow practical rate lead applications, to the segregated the segregated flow pattern, flow similarApattern relatively to thecorresponds annular/stratified high gas to gas-locking, flow rateflows andan in operational pipe. low liquidAs it failure can flow be which seen rate gives leadin Figure tonull the pump11d, segregated the head elongated and flownull bubbleflow pattern, similarexpandsrate to [70]. the and annular/stratified occupies the entire flows impeller in pipe. length. As it For can practical be seen inapplications, Figure 11d, the the segregated elongated bubbleflow expandspattern and Ascorresponds a occupies direct observation theto gas-locking, entire of impeller flow an patterns operational length. inside For failure practicalthe ESP which impeller, applications, gives visualization null pump the segregated head experiments and flow null can patternflow correspondsratehelp [70]. reveal to gas-locking, gas-liquid flow an operationalbehaviors intuitively. failure which However, gives the null experimental pump head facility and null needs flow special rate [70]. designsAs a direct associated observation with necessary of flow patterns modifications inside onthe pumpESP impeller, geometries, visualization such as theexperiments removal ofcan help reveal gas-liquid flow behaviors intuitively. However, the experimental facility needs special designs associated with necessary modifications on pump geometries, such as the removal of
Energies 2018, 11, 180 15 of 41
As a direct observation of flow patterns inside the ESP impeller, visualization experiments can help reveal gas-liquid flow behaviors intuitively. However, the experimental facility needs special Energies 2018, 11, 180 15 of 41 designs associated with necessary modifications on pump geometries, such as the removal of impeller impellerhub and hub attaching and attaching the Pyrex the glass Pyrex on glass its top on for its visualized top for visualized observation. observation. Barrios [Barrios39] pointed [39] pointed out that outthe that visualization the visualization of two-phase of two-phase flow structures flow structures in a multistage in a multistage ESP assembled ESP assembled in series in series was much was muchmore more difficult. difficult. Thus, Thus, how how to characterize to characterize the the multiphase multiphase flow flow in in ESPs ESPs has has becomebecome a challenging challenging topic.topic. Although Although several several technologies technologies [68,69,73–75] [68,69,73–75] to to visualize visualize the the internal internal flows flows inside inside a a centrifugal centrifugal pumppump were were discussed, discussed, they they required required some some modifi modificationscations on on pump pump geometry geometry and and the the implementation implementation waswas difficult difficult to to carry carry out. out. It It required required mounting mounting Hi HireCTreCT into into the the apparatus apparatus so so that that the the internal internal flow flow structuresstructures could could be be visualized visualized regardless regardless of of the the opaque opaque pump pump casing casing or or volute. volute. In In addition, addition, the the data data processingprocessing involved involved time-averaged time-averaged rotation-synch rotation-synchronizedronized CT CT scanning scanning techniques, techniques, adding adding further further complexitycomplexity in in analyzing analyzing the the obtained obtained experimental experimental results. results. FigureFigure 1212 showsshows the the measurements measurements using using the theHireCT HireCT technique technique by Schäfer by Schäfer et al. et [68]. al. The [68]. The horizontal axis denotes the inlet GVF and the vertical axis is the volumetric averaged in-situ α in horizontal axis denotes the inlet GVF and the vertical axis is the volumetric averaged in-situ αG in Gthe the rotating pump impeller. Clearly, there is a sharp α increase at a GVF corresponding to the severe rotating pump impeller. Clearly, there is a sharp αG increaseG at a GVF corresponding to the severe gas-pocketgas-pocket formation formation and and pumping pumping head head degradation. degradation. The The step step change change is is between between inlet inlet GVF GVF = = 2.5% 2.5% andand 3.0%. 3.0%. As As confirmed confirmed by by Schäfer Schäfer et et al. al. observation, observation, such such change change is is due due to to the the rapid rapid flow flow pattern pattern changechange from from bubbly bubbly flow flow to to intermittent intermittent flow. flow.
Figure 12. In-situ gas void fractions at various intake GVFs measured by HireCT [68]. Figure 12. In-situ gas void fractions at various intake GVFs measured by HireCT [68].
The aforementioned experimental studies on ESP gas-liquid performance used tap water as the workingThe fluid, aforementioned while the experimentalgas phase was studies compressed on ESP air gas-liquid or nitrogen. performance Several usedrecent tap experimental water as the studiesworking focusing fluid, whileon ESP the gas-handling gas phase ability was compressed under viscous air fluid or nitrogen. flow were Several conducted recent by Trevisan experimental [41], Trevisanstudies focusing and Prado on ESP[76], gas-handling Banjar et al. ability[30], and under Paternost viscous et fluid al. [77]. flow Using were conducteda visualization by Trevisan prototype [41], builtTrevisan from and original Prado ESP [76], components Banjar et al. [and30], with and Paternostminimal etgeometrical al. [77]. Using modifications, a visualization Trevisan prototype [41], Trevisanbuilt from and original Prado ESP[76] componentsconducted experiments and with minimal to investigate geometrical the viscosity modifications, effect on Trevisan liquid/gas [41], two-phaseTrevisan and flow Prado through [76] ESP.conducted The authors experiments identified to four investigate liquid/air the flow viscosity patterns effect inside on the liquid/gas impeller channels:two-phase agglomerated flow through ESP.bubble, The gas authors pocket, identified segreg fourated liquid/air gas and flowintermittent patterns gas inside flows. the impellerIt was concludedchannels: agglomeratedthat the agglomerated bubble, gasbubble pocket, flow segregatedis responsible gas andfor pressure intermittent surging gas flows.phenomenon, It was whichconcluded is the that initiation the agglomerated of pump head bubble deterioration flow is responsible due to gas for entrainment. pressure surging The authors phenomenon, also observed which thatis the the initiation increase of pumpin viscosity head deteriorationcaused surging due to gasoccur entrainment. at relatively The lower authors inlet also GVFs. observed Similar that experimentalthe increase in observations viscosity caused were surgingmade by to Banjar occur atet relativelyal. [30]. Paternost lower inlet et al. GVFs. [77] Similarfurther experimentalinvestigated theobservations performance were of made a centrifugal by Banjar etpump al. [30 handling]. Paternost single-phase et al. [77] further viscous investigated liquids and the analyzed performance the impactof a centrifugal of free gas pump entrainment. handling They single-phase concluded viscous that the liquids degradation and analyzed of pump the head impact was ofdue free to the gas stagnation of large gas-pocket formation, which became worse with the increase of liquid viscosity. An empirical correlation similar to the dimension analysis procedure in Solano [28] was proposed to correlate ESP pressure increment under two-phase flow conditions accounting the inlet GVF and fluid viscosity.
Energies 2018, 11, 180 16 of 41 entrainment. They concluded that the degradation of pump head was due to the stagnation of large gas-pocket formation, which became worse with the increase of liquid viscosity. An empirical correlation similar to the dimension analysis procedure in Solano [28] was proposed to correlate ESP pressure increment under two-phase flow conditions accounting the inlet GVF and fluid viscosity.
4. Modeling Approaches In this section, several empirical correlations for calculating centrifugal pump performance under gas/liquid two-phase flow are reviewed. Although empirical correlations are much easier to implement compared to other numerical or mechanistic models, the application range is rather limited. With this consideration, the more sophisticated one-dimensional two fluid models coupled with interfacial momentum transfer mechanisms are discussed. The 3D CFD simulations, based on the conservation laws of mass and momentum, can provide detailed numerical results in terms of flow structure, phase distribution and slippage etc. However, CFD simulations for multiphase flow, which usually demand high computational cost, are not as reliable or robust as for single-phase flow. Furthermore, both one-dimensional two-fluid model and CFD-based numerical model are extremely sensitive to the closure relationships, such as bubble size, drag force model and other interphase momentum transfer terms, etc. Thus, the mechanistic models rooted in physics, including flow pattern, mass and momentum conservation equations, have become a new direction for developing the best prediction model of boosting pressure in ESPs with gas involvement.
4.1. Emperical Correlations The first empirical correlation based on experimental data of ESP pressure degradation due to free gas entrainment was proposed by Turpin et al. [78]: ! !! qg 346430 qg 410 = − − Hm H exp 2 (8) ql (C1 pin) ql C1 pin where, Hm is the head of ESP with gas/liquid flow (ft); H is the single-phase pump head (ft); qg and ql are gas and liquid flow rates (bpd), respectively; pin is pump intake pressure (psi); C1 = 0.145 is a unit conversion constant. As the authors claimed, a critical constraint of stable flow condition inside ESPs should be applied to this model so as to achieve a reasonable prediction accuracy. Thus, another empirical correlation to distinguish the stable/unstable flow boundaries are given as below: qg/ql φ = 2000 (9) 3pin
The stable flow condition corresponds to φ < 1, which is the recommended application range of Equation (8). Romero [34] established a multiphase head model for the type of mixed flow pumps based on Cirilo’s [33] experimental data: ! H q q 2 q m = 1 − ld a ld + ld + 1 (10) Hmax qdmax qdmax qdmax where Hmax is the shut-in pump head in ESP; qld is the dimensionless liquid flow rate, a ratio of the intake liquid flow rate to the open flow rate (OFR). The application range for Equation (10), as pointed out by the authors, is dispersed bubble flow or low degradation of ESP performance:
a = 2.902λG + 0.2751 (11)
qdmax = 1 − 2.2035λG (12) Energies 2018, 11, 180 17 of 41
To extend the application range of above two empirical correlations, Duran and Prado [79] presented models for the mild and severe head degradations in ESPs, which correspond to bubbly and elongated bubble flow regimes:
( ql qg (1 − α)ρl H 1−α + αρg H α ∆pm = q (1−α) (13) −0.47075 − 0.21626 ln g qmax where ∆pm is the pressure increment of single stage ESP with gas/liquid flow (psi), α is the in-situ gas void fraction. Equation (13) is applied to bubbly flow and elongated bubble flow, respectively. A closure relationship is presented below to solve for α in above equations:
1.622 ρm ql (1−α) qg(1 − α) −0.843 ρ + 0.85 q = l q max (14) qmax 0.435 ql (1−α) 1.6213qmax
Similarly, Equation (14) corresponds to Equation (13), which can be used to calculate the in-situ gas void fraction and then the pressure increment in ESP with gas presence can be obtained accordingly. Although the application range of Equations (13) has been extended compared to Equation (10), its formulation is only based on the experimental data from a single ESP geometry without any validation against different pump models. Thus, its general validity is questioned due to its empirical nature.
4.2. One-Dimension Two-Fluid Modeling Zakem [80] first developed a mathematical model using a one-dimensional control volume method to analyze gas bubbles and liquid interaction for straight blade impellers. Furuya [81] developed a similar analytical model by incorporating the pump geometry, void fraction, flow slippage, and the flow regime but neglecting the compressibility and condensation effects. The basic formula in Furuya study is given by:
. 2 . 2 m ml 1 1 dα − 1 d(dn) + g 1−α 1 dα + 1 d(dn) ρ (dn)2 1−α 1−α ds dn ds ρ (dn)2 α2 α ds dn ds l " g # . 2 . 2 ρ mg d(dn) m d(dn) + l 1 dα + 1 + l 1 dα − 1 (15) 2 ρgαdn α ds dn ds ρl (1−α)dn 1−α ds dn ds
. . 2 2 3 Cd ml mg = ρl − ρg (1 − α)rω sin β cos γ − ρl − 8 rb ρl (1−α)dn ρgαdn
. . ml and mg are the mass flow rates of liquid and gas. s is the streamline coordinate. n is the coordinate normal to the streamline coordinate s. β is the angle between relative flow velocity and circumferential direction. γ is the angle between the radial direction and the stream surface. Cd is the drag coefficient. rb is the bubble radius. Comparing with experimental data in literature, the model predictions were within the relative average error band of ±30% for GVF <20%, and ±50% for GVF >30%. Sachdeva et al. [82,83] conducted a comprehensive investigation of two-phase flow in ESPs with air/water and diesel/CO2 mixtures. A dynamic five-equation, one-dimensional, two-fluid model accounting for pump geometry, intake pressure and GVF, fluid properties, was developed and verified to predict ESP boosting pressure. The basic equations of Sachdeva model are presented below. The mass balance equations are:
d(α ρ W r sin β) G G G = 0 (16) ds
d((1 − α )ρ W r sin β) G L L = 0 (17) ds where r is radius in ESP, WG and WL are the mass flow rates of gas and liquid, respectively. Energies 2018, 11, 180 18 of 41
The momentum balance equations are:
∂V ∂P ∂r α ρ V R,G = −α + α ρ rω2 − F + F + F (18) G G R,G ∂s G ∂s G G ∂s W,G i v ∂V ∂P ∂r (1 − α )ρ V R,L = −(1 − α ) + (1 − α )ρ rω2 − F − F − F (19) G L R,L ∂s G ∂s G L ∂s W,L i v where VR,G and VR,L are the radial components of gas and liquid absolute velocities, respectively. FW,G and FW,L are the friction forces of gas and liquid against the channel walls per unit volume of fluids. Fi is the interfacial friction force between gas and liquid per unit volume. Fv is the virtual mass force per unit volume. Minemura et al. [84] also studied the performance of centrifugal pumps in the nuclear industry under air-water two phase flow conditions with a low inlet GVF (<10%). Based on energy change in the flow from the rotating impeller to the stationary volute, a 1D, two-fluid model considering fluid viscosity and air-phase compressibility in a rotating impeller was proposed. This model can be solved numerically with a prediction error of ±20% of the related flow capacity. However, both Sachdeva et al. and Minemura et al. models are valid only under a narrow application range or specific experimental conditions. Compared to Sachdeva model, the major difference in Minemura et al. model is that the momentum balance equations are based on the relative velocities of gas/liquid rather than the radial components of the absolute velocities, which are given as:
∂W ∂P ∂r α ρ W G = −α + α ρ rω2 − F + F + F (20) G G G ∂s G ∂s G G ∂s W,G i v ∂W ∂P ∂r (1 − α )ρ W L = −(1 − α ) + (1 − α )ρ rω2 − F − F − F (21) G L L ∂s G ∂s G L ∂s W,L i v Based on Sachdeva et al. and Minemura et al. one-dimensional two-fluid models, Sun [12,22] developed a new two-phase model including a set of one-dimensional mass and momentum balance equations for predicting ESP performance. He also improved analytical models for wall frictional losses and shock loss, as well as new correlations for the drag coefficient. The general momentum balance equation along the streamlines is given by:
" # dp dWp 2 dp ds Mp,s ds = −ρpWp + ρpΩ r + − (22) dr dr ds dr αp dr streamline f ,p streamline where the subscript p is g or l for gas or liquid phase, respectively. And Mp,s is the interfacial momentum transfer term. Eliminating the pressure increment term at the left hand side in Equation (24), the combined momentum balance equation can be expressed as
! 2 dp dp 1 ρlWldWl − ρgWgdWg = ρl − ρg Ω rdr + − ds − Ml,s ds (23) ds f ,l ds f ,g αlαg
Applying a finite differencing algorithm, a preferred numerical scheme [85], to solving this model, a good agreement of pump performance curve, αG distribution, surging and gas lock conditions against correspondent experimental measurements was obtained. Although many studies for modeling ESP performance under gassy conditions were conducted, mechanistic modeling is still needed due to the over-simplification and assumptions or narrow application ranges of the existing models.
4.3. Three-Dimension Computational Fluid Dynamics (CFD) With the advances of computer technology, computational fluid dynamics (CFD) becomes a more and more powerful tool to study centrifugal-pump performance under single-phase and multiphase flow conditions. Due to complicated ESP geometries, it is difficult to investigate the internal velocity Energies 2018, 11, 180 19 of 41 and pressure fields experimentally. However, CFD offers an alternative way to simulate the complex internal flow structures. A centrifugal pump consists of an impeller rotating at a set angular velocity and a volute which is stationary. For an ESP, the rotating and stationary parts are the impeller and diffuser, which are accommodated in the rotating and stationary computational domains in CFD, respectively. Using a 3D CFD code with the frozen-rotor interface model, the internal flow inside centrifugal pumps can be simulated, including velocity and pressure fields [86–88] as well as flow recirculation and separation [89,90]. The frozen-rotor model is considered as steady-state simulation because it holds the rotating and stationary parts in two separate reference frames. Some transient simulations were conducted using the sliding-mesh technique to investigate the dynamic flow structures in centrifugal pumps [91–94]. Gonzalez et al. [91] performed numerical simulations of unsteady flow in a single-phase centrifugal pump, considering impeller-volute interaction. By solving viscous, incompressible Navier-Stokes (N-S) equations with the sliding mesh technique, the unsteady flow behavior inside a centrifugal pump due to impeller-volute interaction was captured. A relationship between the global variables, such as torque, as a function of impeller relative position, secondary flow in volute, etc., was obtained numerically. This approach was successful in predicting the dynamic interaction between impeller and volute. Huang et al. [93] studied unsteady flow and pressure fluctuations due to interaction between impeller and diffuser vanes by the sliding mesh technique. His study confirmed that the global variables are primarily affected by impeller blade passing frequency. In addition to designing turbomachinery, CFD can help engineers study the viscosity effects on centrifugal pump performance. Shojaeefard et al. [95,96] conducted both experimental study and numerical simulation of a centrifugal pump handling viscous fluids. The authors stated that a good agreement between simulation and experimental data was obtained by solving the steady state Reynolds-averaged Navier–Stokes (RANS) equations with Shear Stress Transport (SST) k-ω turbulence model. Using the same pump geometry, Sirino et al. [97] and Stel et al. [98] performed numerical investigations of viscosity effects on single-stage and three-stage ESPs, respectively. Similar numerical methodologies were used in their work including SST turbulence model with transient rotor-stator techniques. Both simulation results matched experiments well under a wide range of fluid viscosity. Stel et al. [98,99] pointed out that CFD simulated boosting pressure with multistage ESP geometry agreed with experimental results better than that with a single-stage geometry. The phenomenon of rising head with moderate increase of fluid viscosity was studied by Li [100]. By implementing the standard k-ε turbulence model and non-equilibrium wall function into RANS equations, the author confirmed that the rising pump head was due to the transition from hydraulically rough regime to hydraulically smooth regime. Zhu et al. [84] solved a set of 3D, steady-state RANS equations with standard SST turbulence model using the commercial CFD solver ANSYS CFX by employing the frozen-rotor technique. Their results matched correspondent experimental data well. CFD has also been used to simulate pump performance with gas involvement, such as cavitation phenomenon [101,102], and free-gas entrainment flow [103–105]. Unlike single-phase simulations, the two-phase simulation requires the solution of conservation equations of mass, momentum and energy for the continuous and dispersed phases. Meanwhile, another set of constitutive equations, for describing the interphase interactions, such as interfacial momentum/heat transfer, need to be solved simultaneously. Minemura and Uchiyama [106] solved 3D momentum equations with a finite-element method to predict gas/liquid flow behavior in a rotating centrifugal pump impeller. Tremante et al. [107] numerically studied gas/liquid flow through a cascade axial pump by CFD simulation. Coupled with a modified k-ε turbulence model by considering the viscosity of the liquid phase and the compressibility of the gas phase, the gas/liquid distributions versus different attack angles were obtained. Caridad et al. [108,109] studied ESP impeller performance handling water/air mixtures using CFD simulation. Applying two fluid models in 3D CFD simulations, the pressure and velocity fields, as Energies 2018, 11, 180 20 of 41 energy for the continuous and dispersed phases. Meanwhile, another set of constitutive equations, for describing the interphase interactions, such as interfacial momentum/heat transfer, need to be solved simultaneously. Minemura and Uchiyama [106] solved 3D momentum equations with a finite-element method to predict gas/liquid flow behavior in a rotating centrifugal pump impeller. Tremante et al. [107] numerically studied gas/liquid flow through a cascade axial pump by CFD simulation. Coupled with a modified k-ε turbulence model by considering the viscosity of the liquid phase and the compressibilityEnergies 2018, 11, 180 of the gas phase, the gas/liquid distributions versus different attack angles 20were of 41 obtained. Caridad et al. [108,109] studied ESP impeller performance handling water/air mixtures using CFD simulation. Applying two fluid models in 3D CFD simulations, the pressure and velocity fields,well as as the well gas as phasethe gas distributions phase distributions (Figure (Figure 13) were 13) obtained. were obtained. The gas The pocket gas pocket near thenear impeller the impeller blade bladewas also was identified also identified and compared and compared with experimental with experimental observations. observations.
FigureFigure 13. 13. TypicalTypical multiphase CFD simulation resu resultslts of gas/liquid flow flow in in ESPs ESPs [108]. [108].
The sensitivity analysis on GVF and bubble diameter indicated that ESP performance The sensitivity analysis on GVF and bubble diameter indicated that ESP performance deteriorated deteriorated with GVF and bubble diameter increase. Barrios [39], Barrios et al. [103] conducted with GVF and bubble diameter increase. Barrios [39], Barrios et al. [103] conducted multiphase multiphase CFD simulations on a single-stage ESP impeller with the new models of bubble size and CFD simulations on a single-stage ESP impeller with the new models of bubble size and drag drag coefficient predictions. Their simulations agreed with laboratory visualization images of coefficient predictions. Their simulations agreed with laboratory visualization images of streamlines streamlines and gas-accumulation zones. Qi et al. [110] designed ESPs for geothermal application and gas-accumulation zones. Qi et al. [110] designed ESPs for geothermal application with with high-temperature gas-liquid two-phase flow. Using CFD simulations, the designed mix-type high-temperature gas-liquid two-phase flow. Using CFD simulations, the designed mix-type centrifugal impeller and diffuser were optimized for better gas-handling and higher efficiency within centrifugal impeller and diffuser were optimized for better gas-handling and higher efficiency within a wide range of production rate. Zhu and Zhang [104] conducted multiphase CFD simulations on a a wide range of production rate. Zhu and Zhang [104] conducted multiphase CFD simulations three-stage ESP with each stage comprising of an impeller and a diffuser. Comparing to experimental on a three-stage ESP with each stage comprising of an impeller and a diffuser. Comparing to measurements, their work revealed that bubble size was a critical factor affecting ESP performance experimental measurements, their work revealed that bubble size was a critical factor affecting ESP under gassy conditions. A new bubble size prediction model was later proposed in Zhu and Zhang [105] performance under gassy conditions. A new bubble size prediction model was later proposed in Zhu based on CFD simulations. They also predicted the in-situ gas void fraction (αG) with a mechanistic and Zhang [105] based on CFD simulations. They also predicted the in-situ gas void fraction (α ) with model and validated the results with numerically simulated values (Zhu and Zhang [90]). A typicalG a mechanistic model and validated the results with numerically simulated values (Zhu and Zhang [90]). CFD simulation result of gas/liquid flow inside a rotating ESP from their studies is shown in A typical CFD simulation result of gas/liquid flow inside a rotating ESP from their studies is shown in Figure 14 below. Figure 14 below. Multiphase flow phenomena in ESP are transient in nature, such as fluctuations of local pressure, Multiphase flow phenomena in ESP are transient in nature, such as fluctuations of local pressure, gas void fractions, and breakup or coalescence of bubbles. To better simulate the hydrodynamics of gas void fractions, and breakup or coalescence of bubbles. To better simulate the hydrodynamics of gas-liquid two-phase flow in a rotating centrifugal pump, the unsteady CFD simulation code coupled gas-liquid two-phase flow in a rotating centrifugal pump, the unsteady CFD simulation code coupled with multiphase flow model and transient rotor-stator algorithm to account for the interactions with multiphase flow model and transient rotor-stator algorithm to account for the interactions between impeller and diffuser should be adopted. However, the computational cost for transient between impeller and diffuser should be adopted. However, the computational cost for transient CFD CFD simulations is far more than that for steady-state simulations. Marsis et al. [111] performed simulations is far more than that for steady-state simulations. Marsis et al. [111] performed transient transient two-phase CFD simulations on eight multi-vane ESP designs. The predicted pump two-phase CFD simulations on eight multi-vane ESP designs. The predicted pump performance was performance was confirmed by experimental results. The final design was achieved with the stage confirmed by experimental results. The final design was achieved with the stage pressure increased by 4% for single-phase flow and 23% for two-phase flow at inlet GVF = 20% by optimizing the meridional profile and number of blades. Yu et al. [112,113] conducted unsteady numerical simulation on gas-liquid flow in a multiphase centrifugal pump. Considering multiple interfacial momentum transfer components, such as drag force, lift force, virtual mass force (VMF) and turbulent dispersion force (TDF), they concluded that the two-fluid multiphase model was able to capture the transport process more accurately than the homogeneous model. Compared to turbulent dispersion force, the drag force plays a more dominant role, as shown in Figure 15a. In Figure 15, z is in the streamwise direction, while L is the maximum Energies 2018, 11, 180 21 of 41 Energies 2018, 11, 180 21 of 41 pressure increased by 4% for single-phase flow and 23% for two-phase flow at inlet GVF = 20% by optimizingpressure increased the meridional by 4% pr forofile single-phase and number flow of blades.and 23% for two-phase flow at inlet GVF = 20% by optimizingYu et al. the[112,113] meridional conducted profile unsteady and number numerical of blades. simulation on gas-liquid flow in a multiphase centrifugalYu et pump. al. [112,113] Considering conducted multiple unsteady interfacial numerical momentum simulation transfer on gas-liquid components, flow in such a multiphase as drag force,centrifugal lift force, pump. virtual Considering mass force multiple(VMF) and interfacial turbulent momentum dispersion transfer force (TDF), components, they concluded such as that drag force, lift force, virtual mass force (VMF) and turbulent dispersion force (TDF), they concluded that theEnergies two-fluid2018, 11 ,multiphase 180 model was able to capture the transport process more accurately than21 ofthe 41 homogeneousthe two-fluid model. multiphase Compared model to was turbulent able to dispersi captureon the force, transport the drag process force moreplays accuratelya more dominant than the role,homogeneous as shown in model. Figure Compared 15a. In Figure to turbulent 15, z is indispersi the streamwiseon force, the direction, drag force while plays L is a themore maximum dominant streamlinerole, as shown asas thethe in fluidfluid Figure flows flows 15a. from from In theFigure the intake intake 15, toz to dischargeis indischarge the streamwise of theof the computational computational direction, while domain. domain. L is As the bubbleAs maximum bubble size streamline as the fluid flows from the intake to discharge of the computational domain. As bubble size(Db) ( mayDb) may vary vary in a rotatingin a rotating centrifugal centrifugal pump, pump, its effect its effect on the on local the draglocal force drag isforce limited, is limited, which which can be canverifiedsize be ( Dverified inb) may Figure invary Figure15b. in a 15b.rotating centrifugal pump, its effect on the local drag force is limited, which can be verified in Figure 15b.
(a) (b) (a) (b) Figure 14. CFDCFD predicted predicted gas gas void void fraction fraction distribution distribution and and the the corresponding corresponding pressure pressure contour contour [90], [90], Figure 14. CFD predicted gas void fraction distribution and the corresponding pressure contour [90], (a) gas void fraction color map; (b) pressure contourcontour onon halfhalf impellerimpeller span.span. (a) gas void fraction color map; (b) pressure contour on half impeller span.
(a) (b) (a) (b) Figure 15. Magnitude analysis of interface forces by CFD simulations [112], (a) magnitude ratios of non-dragFigureFigure 15. 15. forcesMagnitude Magnitude to drag analysis force;analysis (b of) ofbubble interface interface size forces foeffectrces by onby CFD dragCFD simulations forcesimulations magnitude. [112 [112],], (a ()a magnitude) magnitude ratios ratios of of non-dragnon-drag forces forces to to drag drag force; force; (b ()b bubble) bubble size size effect effect on on drag drag force force magnitude. magnitude.
Huang et al. [94] investigated the transient inhaling process of gas initially filling a section of pipe into centrifugal pump impeller by sliding mesh and transient-frozen-rotor methods. The phase distribution, pressure and velocity versus time were computed and analyzed. Pineda et al. [114] presented an alternative approach to obtain the distribution and concentration of the dispersed phase in a rotating centrifugal pump by solving the realizable k-ε turbulence model coupled with the volume of fluid (VOF) multiphase model. They observed that the numerically simulated results of in-situ αG could be correlated by the Lockhart-Martinelli parameter. Zhang et al. [115,116] compared the Energies 2018, 11, 180 22 of 41 unsteady CFD simulation results of multiphase flow patterns in a three-stage centrifugal pump with the corresponding visualization experiments. Good agreement on the positions and shapes of the gas pocket was observed. Ye et al. [117] combined the advanced transient CFD multiphase simulation with finite-element-analysis (FEA) to optimize the 3D metal printing process of hybrid stage prototype and ESPs for high-gas wells. A design-validation tool was developed and a prototype ESP was manufactured which can pump up to 75% GVF gas/liquid mixture without gas locking.
4.4. Mechanistic Modeling Unlike multiphase pipe flow, whose flow structures and characteristics are easier to visualize and quantify, the mechanistic modeling of ESP performance under gas/liquid flow is more difficult due to its highly twisted flow channel and compact assembly. At present, only limited studies are available in the literature related to this topic, such as Barrios [39], Gamboa [40] and Zhu [1]. Similar to the mechanistic model of multiphase pipe flow, the flow patterns in ESP two-phase flow should be identified first, based on which the corresponding flow models can be developed. Based on Sachdeva [82] and Sun [12] momentum equations in Equation (24) along the streamline inside a rotating centrifugal pump, Zhu [1] derived the combined momentum equation of gas and liquid phases by assuming the Taylor bubble flow pattern prevailing in rotating ESP two-phase flow:
− − − − − ρL(vT vF)(vS vF) ρC(vT vC)(vS vC) − τFSF − lF HLF A 1 1 2 (24) τI SI + + (ρL − ρC)Ω RI = 0. HLF A (1−HLF)A streamline where the subscripts L, F, S, C, LF, and I denote liquid, film, slug, gas core, liquid film and interface, respectively. vT is the translational velocity. Compared to the combined momentum equation for slug flow in pipe [115,116,118–120], the difference in Equation (24) is the body force term. All the velocities are the relative velocities to the ESP channel. By applying the flow pattern and necessary assumptions, Equation (24) can be reduced to the simplified flow model, under which the flow characteristics and pump performance can be solved.
4.4.1. Flow Pattern Map and Transition Boundary As aforementioned, the visualization of flow patterns inside ESP with gas/liquid flow either involves difficult experimentation to install transparent glasses so as to enable the visualization of internal flow fields [39,70] or requires the complicated and expensive non-intrusive instrumentation [66]. However, the identification of flow patterns and transition boundaries inside a rotating ESP impeller is still not well understood. Gamboa and Prado [72] proposed an alternative approach to recognize the flow patterns of ESP multiphase flow by relating the inflection points on pump performance curves. Figure 16 below illustrates the identification process schematically. Gamboa and Prado [72] ascribed the inflections on each performance curve to the flow pattern transitions inside an ESP with gas/liquid flow. As the authors pointed out that such analogy can merely achieve a rough accordance, it did offer significant convenience to obtain the transition boundary among flow patterns. In Figure 16, the horizontal axis is normalized liquid flow rate, and the vertical axis is the measured pump pressure increment. Compared to single-phase performance curve with no inflections, it is found that the two-phase performance curve suffers from significant pressure degradation with inflections at varying loci as qld drops. Different markers representing the inflection points on each mapping test performance curve correspond to different flow pattern transitions. The white circles denote the transition boundary between dispersed bubble flow and bubbly flow; green diamond shows the transition from bubbly flow to intermittent flow; green square corresponds to the transition from intermittent flow to segregated flow. Energies 2018, 11, 180 23 of 41 Energies 2018, 11, 180 23 of 41
Figure 16. FlowFlow pattern pattern recognition by inte interpretingrpreting performance curves [[72].72].
TransitionGamboa from and Dispersed Prado [72] Bubble ascribed Flow the to inflections Bubbly Flow on (1steach Transition performance Boundary) curve to the flow pattern transitions inside an ESP with gas/liquid flow. As the authors pointed out that such analogy can merelyAs achieve discussed a rough above, accordance, the dispersed it bubbledid offer flow significant is considered convenience as no-slip to flow.obtain On the the transition contrary, boundarybubbly flow, among featured flow bypatterns. phase In slippage, Figure 16, corresponds the horizontal to the axis uneven is normalized dispersion liquid of gasflow bubbles rate, and in theliquids. vertical For axis a rotating is the measured ESP, the initiationpump pressure of pressure increment. surging Compared is closely to related single-phase to the 1stperformance transition curveboundary with [44no, 72inflections,]. Due to theit is lackfound of athat valid the mechanistictwo-phase performance model to predict curve the suffers initiation from ofsignificant pressure surging inside ESPs, the empirical correlations have been developed from multi-variable regression of pressure degradation with inflections at varying loci as qld drops. Different markers representing the inflectionexperimental points data. on Table each9 belowmapping summarizes test perfor severalmance existing curve pertinentcorrespond studies. to different flow pattern transitions. The white circles denote the transition boundary between dispersed bubble flow and Table 9. Correlations of surging initiation (1st transition boundary) in literature. bubbly flow; green diamond shows the transition from bubbly flow to intermittent flow; green square correspondsStudy to the transition from intermittent Correlation flow to segregated flow. Application Range Transition from Dispersed Bubble Flow to Bubbly2000 FlowQG (1st Transition Boundary)φ—empirical constant Turpin et al. [78] φ = 3P Q i L Pi < 2.8 mPa As discussed above, the dispersed bubble flow is considered as no-slip flow. On the contrary, 0.4342 Pi < 3.4 mPa bubbly flow,Cirilo featured [33] by phase slippage,λC corresponds= 0.0187Pi to the uneven dispersion of gas3 bubbles in QL < 0.02 m /s liquids. For a rotating ESP, the initiation of pressure surging is closely related to the 1st transition Estevam [64] λC = 31.92 − 32.15(1 − λG) boundary [44,72]. Due to the lack of a valid mechanistic model to predict the initiation of pressure P < 2.4 mPa 1.421 i surging insideDuran ESPs, [37] the empiricalQ correlationsG ρG have beenQ Ldeveloped from multi-variableQ < 0.02 m3 /sregression Q = 5.58 ρ + 0.098 Q G max L max 3 of experimental data. Table 9 below summarizes several existing pertinent studies.QL < 0.013 m /s
Figure 17 shows the comparison of surging initiationQ models in ESPs. The sameP < 1.4 flow mPa conditions 0.027 L i QG Qmax 3 Zapata [38] = Q sep QG < 0.02 m /s with rotational speed N = 3500 rpm andQmax separator0.9001− pressureL at P = 1034 kpa are applied to all Qmax 3 calculations. As can be seen, the predictions of ESP surging initiation from QexistingL < 0.016 mmodels/s vary 0.2 2 0.4 4.4682 QG ρG ΩD QL significantly.Gamboa and Turpin Prado [et72 ]al. [78] and= Cirilo [33] suggested0.102 exp that pressure surging shouldPi < 1.7 be mPa independent Qmax ρL ν Qmax on liquid flow rates (QL). Other studies [37,38,72] claimed1 that the critical gas volumetric fraction (λC) 2 2 0.4σ at which ESP pressure surging initiates is in a close(ρ −ρ )relationshipΩ2 R with QL. Duran [37] and Zapata [38] Zhu et al. [44] = L G I λC 3 − 2 1 PQ ρ correlations predict λC in monotonic-increase10.056 trendσ 5 ∆ withL 5respectL 5 to QL. Gamboa and Prado correlated ρL ρLV ρG λC and QL with a concave quadratic function, and claimed that this model was validated for