Investigation of Taylor bubble dynamics in annular conduits with counter-current flow

被引:12
|
作者
Liu, Yaxin [1 ]
Mitchell, Travis [2 ]
Upchurch, Eric R. [3 ]
Ozbayoglu, Evren M. [1 ]
Baldino, Silvio [1 ]
机构
[1] Univ Tulsa, Tulsa Univ Drilling Res Project, Tulsa, OK 74104 USA
[2] Univ Queensland, Sch Mech & Min Engn, St Lucia, Qld 4072, Australia
[3] Chevron Technol Co, Houston, TX 77002 USA
关键词
Taylor bubble; Counter-current flow; Translational velocity; Distribution parameter; Inclined eccentric annuli; SLUG LENGTH DISTRIBUTION; LIQUID SLUG; 2-PHASE FLOW; REGIME TRANSITION; CHURN TRANSITION; MULTIPHASE FLOW; GAS; MODEL; VELOCITY; MOTION;
D O I
10.1016/j.ijmultiphaseflow.2023.104626
中图分类号
O3 [力学];
学科分类号
08 ; 0801 ;
摘要
This study numerically investigates counter-current slug flow by considering the motion of a Taylor bubble in annular conduits with downward-flowing liquids using the Volume-of-Fluid method implemented in the commercial computational fluid dynamics software ANSYS Fluent (Release 19.2). The translational velocity of a counter-current ascending or co-current descending Taylor bubble in vertical concentric annuli and the corresponding distribution parameter (C-0) are analyzed. The latter is correlated in terms of Eotvos number (Eo) and inverse viscosity number (Nf) within the range of Eo between 40 and 400 and Nf between 40 and 320. The proposed correlation provides an accurate fit to the numerical data with an average error of 2.64%, and is successfully compared with published numerical findings. In general, the smooth and stable shape of the bubble is disrupted as the counter-current flow velocity (Fr-l) increases above a critical value, leading to the formation of surface waves and the displacement of the bubble tip away from the annular gap center and towards the outer pipe. C-0 increases with Eo and Nf, plateauing at high values of Eo. The effects of annulus inclination (theta) and eccentricity (epsilon) on bubble rise velocity are examined within the common range of theta and epsilon encountered during the drilling of oil, gas, or geothermal wells, i.e. 0 degrees <= theta <= 60 degrees and 0 <= epsilon <= 0.7, and their impact on the C-0. The increasing Fr-l and theta lead to a streamlined bubble with pointed nose and thus a reduction in the wrap angle (theta(wrap)), ultimately leading to reduced drag compared to the vertical annulus case and a decrease in C-0. As the epsilon increases, which is accompanied by an increase in the degree of bubble eccentricity, the corresponding C-0 decreases. For a constant Eo = 100 and Nf = 160 with inclination angle of theta = 40 degrees and eccentricity of epsilon = 0.5, the C-0 < 1 is observed.
引用
收藏
页数:16
相关论文
共 50 条