Planar dynamics of inclined curved flexible riser carrying slug liquid–gas flows

Flexible risers transporting hydrocarbon liquid–gas flows may be subject to internal dynamic fluctuations of multiphase densities, velocities and pressure changes. Previous studies have mostly focused on single-phase flows in oscillating pipes or multiphase flows in static pipes whereas understanding of multiphase flow effects on oscillating pipes with variable curvatures is still lacking. The present study aims to numerically investigate fundamental planar dynamics of a long flexible catenary riser carrying slug liquid–gas flows and to analyse the mechanical effects of slug flow characteristics including the slug unit length, translational velocity and fluctuation frequencies leading to resonances. A two-dimensional continuum model, describing the coupled horizontal and vertical motions of an inclined flexible/extensible curved riser subject to the space–time varying fluid weights, flow centrifugal momenta and Coriolis effects, is presented. Steady slug flows are considered and modelled by accounting for the mass–momentum balances of liquid–gas phases within an idealized slug unit cell comprising the slug liquid (containing small gas bubbles) and elongated gas bubble (interfacing with the liquid film) parts. A nonlinear hydrodynamic film profile is described, depending on the pipe diameter, inclination, liquid–gas phase properties, superficial velocities and empirical correlations. These enable the approximation of phase fractions, local velocities and pressure variations which are employed as the time-varying, distributed parameters leading to the slug flow-induced vibration (SIV) of catenary riser. Several key SIV features are numerically investigated, highlighting the slug flow-induced transient drifts due to the travelling masses, amplified mean displacements due to the combined slug weights and flow momenta, extensibility or tension changes due to a reconfiguration of pipe equilibrium, oscillation amplitudes and resonant frequencies. Single- and multi-modal patterns of riser dynamic profiles are determined, enabling the evaluation of associated bending/axial stresses. Parametric studies reveal the individual effect of the slug unit length and the translational velocity on SIV response regardless of the slug characteristic frequency being a function of these two parameters. This key observation is practically useful for the identification of critical maximum response.     

Gas/shear-thinning liquid flows through pipes: Modeling and experiments

In chemical and oil industry gas/shear-thinning liquid two-phase flows are frequently encountered. In this work, we investigate experimentally the flow characteristics of air/shear-thinning liquid systems in horizontal and slightly inclined smooth pipes. The experiments are performed in a 9-m-long glass pipe using air and three different carboxymethyl cellulose (CMC) solutions as test fluids. Flow pattern maps are built by visual observation using a high-speed camera. The observed flow patterns are stratified, plug, and slug flow. The effects of the pipe inclination and the rheology of the shear-thinning fluid in terms of flow pattern maps are presented. The predicted existence region of the stratified flow regime is compared with the experimental observation showing a good agreement. A mechanistic model valid for air/power-law slug flow is proposed and model predictions are compared to the experimental data showing a good agreement. Slug flow characteristics are investigated by the analysis of the signals of a capacitance probe: slug velocity, slug frequency, and slug lengths are measured. A new correlation for the slug frequency is proposed and the results are promising.     

An appraisal of liquid–liquid slug flow in different pipe orientations

Gas–liquid slug flow occurs over a wide range of phase flow rates and in a variety of practical applications during gas–liquid two-phase flows. The range of slug flow increases further in narrow pipes (<0.0254 m), undulated pipelines, riser tube, etc. On the other hand, the past literature shows that slug flow is rarely observed for liquid–liquid cases. In the present study, an interest was felt to investigate whether liquid–liquid slug flow occurs in situations known for excessive slugging in gas–liquid cases. For this, experiments have been performed in narrow (0.012 m ID) vertical and horizontal pipes and an undulated pipeline of 0.0254 m internal diameter where the V-shaped undulation comprises of an uphill and a downhill section between two horizontal pipes. The studies have been performed for both peak and valley orientation of the undulation. Kerosene and water have been selected as the test fluids and the optical probe technique has been used to supplement visual observations especially at higher flow rates. The studies have revealed the existence of the slug flow pattern over a wide range of phase flow rates in all the three geometries. Interestingly, it has been noted that the introduction of an undulation induces flow patterns which bear a closer resemblance to gas–liquid flows as compared to liquid–liquid flows through a horizontal pipe of 0.0254 m diameter.     

Multiphase drag reduction: Effect of eliminating slugs

Two-phase pressure drop measurements were taken for air/water mixtures in a 0.052-m diameter horizontal pipe with special focus on the superficial liquid velocity range of 0.03–1.2 m/s at superficial gas velocities of 3.8, 5.2, and 6.6 m/s. It was found that the addition of 400 ppm of sodium dodecyl sulfate (SDS) to the water reduced the pressure drop by 25–40% when compared to equal flow rates without SDS. The pressure drop reduction occurred where the SDS eliminated the occurrence of the intermittent flow present with water. It was also found that the same concentration of SDS had virtually no effect on single phase liquid pressure drop. The pressure drop reduction appears to be due solely to the suppression of intermittent flow patterns.     

Phase split of nitrogen/non-Newtonian fluid two-phase flow at a micro-T-junction

An experimental investigation has been undertaken to understand the phase split of nitrogen gas/non-Newtonian liquid two-phase flow passing through a 0.5 mm T-junction that oriented horizontally. Four different liquids, including water and aqueous solutions of carboxymethyl cellulose (CMC) with different mass concentrations of 0.1, 0.2 and 0.3 wt%, were employed. Rheology experiments showed that different from water, CMC solutions in this study are pseudoplastic non-Newtonian fluid whose viscosity decreases with increasing the shear rate. The inlet flow patterns were observed to be slug flow, slug–annular flow and annular flow. The fraction of liquid taken off at the side arm for nitrogen gas/non-Newtonian liquid systems is found to be higher than that for nitrogen gas/Newtonian liquid systems in all inlet flow patterns. In addition, with increasing the pseudoplasticity of the liquid phase, the side arm liquid taken off increases, but the increasing degree varies with each flow pattern. For annular flow, the increasing degree is much greater than those for slug and slug–annular flows.     

Single-beam gamma densitometry measurements of oil–water flow in horizontal and slightly inclined pipes

Gamma densitometry is a frequently used non-intrusive method for measuring component volume fractions in multiphase flow systems. The application of a single-beam gamma densitometer to investigate oil–water flow in horizontal and slightly inclined pipes is presented. The experiments are performed in a 15 m long, 56 mm diameter, inclinable stainless steel pipe using Exxsol D60 oil (viscosity 1.64 mPa s, density 790 kg/m³) and water (viscosity 1.0 mPa s, density 996 kg/m³) as test fluids. The test pipe inclination is changed in the range from 5° upward to 5° downward. Experimental measurements are reported at three different mixture velocities, 0.25, 0.50 and 1.00 m/s, and the inlet water cut is varied from 0 to 1. The gamma densitometer is composed of radioactive isotope of Am-241 with the emission energy of 59.5 keV, scintillation detector [NaI(Tl)] and signal processing system. The time averaged cross-sectional distributions of oil and water phases are measured by traversing the gamma densitometer along the vertical pipe diameter. Based on water volume fraction measurements, water hold-up and slip ratio are estimated. The total pressure drop over the test section is measured and frictional pressure drop is estimated based on water hold-up measurements. The measurement uncertainties associated with gamma densitometry are also discussed. The measured water hold-up and slip ratio profiles are strongly dependent on pipe inclination. In general, higher water hold-up values are observed in upwardly inclined pipes compared to the horizontal and downwardly inclined pipes. At low mixture velocities, the slip ratio decreases as the water cut increases. The decrease is more significant as the degree of inclination increases. The frictional pressure drop for upward flow is slightly higher than the horizontal flow. In general, there is a marginal difference in frictional pressure drop values for horizontal and downwardly inclined flows.     

Flow pattern transition,pressure gradient,hold-up predictions in gas/non-Newtonian power-law fluid stratified flow

This work focuses on gas/non-Newtonian power-law fluid stratified pipe flow. Two different theoretical approaches to obtain pressure gradient and hold-up predictions are presented: the steady fully developed two-fluid model and the pre-integrated model. The theoretical predictions are compared with experimental data available for horizontal and for slightly downward inclined air/shear thinning fluid stratified flow taken from literature. The predictions of the pre-integrated model are validated showing a good agreement when compared with experimental data. The criteria for the transition from the stratified flow pattern are applied to gas/non-Newtonian stratified flow. The neutral stability analysis (smooth/wavy stratified flow) and the well-posedness (existence region of stratified flow) of governing equations are carry out. The predicted transition boundaries are obtained using the steady fully developed two-fluid model and the pre-integrated model, where the shape factors and their derivatives are accounted for. A comparison between the predicted boundaries and experimental flow pattern maps is presented and shows a good agreement. A comment on the shear stress modeling by the pre-integrated model is provided.     

Pressure drop in two phase slug/bubble flows in mini scale capillaries

A segmented two phase slug/bubble flow occurs where a liquid and a gas are pumped into the same tube over a range of Reynolds numbers. This segmented two phase flow regime is accompanied by an increase in pressure drop relative to the single phase flow where only one fluid is flowing in a capillary. This work experimentally and theoretically examines the pressure drop encountered by the slug/bubble flow with varying slug lengths in mini channels. In the experimental work the dimensionless parameters of Reynolds number and Capillary number span over three orders of magnitude, and dimensionless slug length ranges over two orders of magnitude to represent flows typical of mini- and micro-scale systems. It is found, in agreement with previous work, that these dimensionless groups provide the correct scaling to represent the pressure drop in two phase slug/bubble flow, although the additional pressure drop caused by the interface regions was found to be ∼40% less than previously reported.     

Flow patterns of n-hexadecane–CO2 liquid–liquid two-phase flow in vertical pipes under high pressure

To promote a better understanding of liquid–liquid two-phase flow behavior, particularly under high pressure, flow patterns of n-hexadecane–CO₂ liquid–liquid two-phase upward flow in vertical stainless steel pipes were experimentally investigated. Observations were made in two 0.0015 m I.D. pipes of different lengths (0.068 m and 0.5 m) under high pressure varying from 10.3 to 29.6 MPa using a high pressure visualization system. The total flow rate was fixed at 2.0 × 10⁻⁶ m³/min, while the flow rate ratio (φ) varied from 0.05 to 19. Bubbly flow, plug flow, slug flow, annular flow, and near-one-phase flow regions were found in both pipes, while stratified flow was observed only in the 0.068 m pipe. Flow pattern maps were constructed in the flow rate ratio versus pressure graph, which demonstrates significant impacts of flow rate ratio, pipe length, and pressure on flow patterns. These impacts are discussed in detail. To the authors’ best knowledge, this work is the first attempt to observe complex liquid–liquid two-phase flow behavior with flow pattern transitions under high pressure, and contributes to a better understanding of liquid–liquid two-phase flow behavior.     

Effect of drag-reducing polymers on pseudo-slugs––interfacial drag and transition to slug flow

Drag-reducing polymers were added to air and water flowing in a stratified configuration in a horizontal 2.54 cm pipe. The interface was covered with large amplitude roll waves, that have been called pseudo-slugs, over a range of flow conditions. The damping of small wavelength waves causes a large decrease in the interfacial stress and, therefore, an increase in the liquid holdup. At superficial gas velocities greater than 4 m/s the transition to slug flow is delayed in that it occurs at larger liquid holdups. This observation is interpreted by assuming that turbulence in slugs is damped. This increases the shedding rate of a slug and, therefore, its stability. The pressure drop can increase or decrease when polymers are added. The increase in holdup is accompanied by an increase in gas velocity, which causes an increase in the pressure drop. The decrease in the interfacial stress has the opposite effect.     

Numerical modeling and experimental investigation of gas–liquid slug formation in a microchannel T-junction

Gas–liquid two-phase flow in a microfluidic T-junction with nearly square microchannels of 113 μm hydraulic diameter was investigated experimentally and numerically. Air and water superficial velocities were 0.018–0.791 m/s and 0.042–0.757 m/s, respectively. Three-dimensional modeling was performed with computational fluid dynamics (CFD) software FLUENT and the volume of fluid (VOF) model. Slug flow (snapping/breaking/jetting) and stratified flow were observed experimentally. Numerically predicted void fraction followed a linear relationship with the homogeneous void fraction, while experimental values depended on the superficial velocity ratio U_G/U_L. Higher experimental velocity slip caused by gas inlet pressure build-up and oscillation caused deviation from numerical predictions. Velocity slip was found to depend on the cross-sectional area coverage of the gas slug, the formation of a liquid film and the presence of liquid at the channel corners. Numerical modeling was found to require improvement to treat the contact angle and contact line slip, and could benefit from the use of a dynamic boundary condition to simulate the compressible gas phase inlet reservoir.     

Experimental investigation of air–oil slug flow using capacitance probes,hot-film anemometer,and image processing

Despite the importance of air–oil slug flows to many industrial applications, their available data reported in the literature are limited compared to air–water slug flows. The main objective of the present study is to explain how air–oil slug flow parameters can be experimentally investigated using hot-film anemometry, capacitance sensors and image processing. Experiments were performed using air–oil slug flow through a horizontal pipe for air superficial velocities ranged from 0.01 m/s to 0.65 m/s and oil superficial velocities ranged from 0.03 m/s to 2.3 m/s. The signal obtained from the hot-film anemometer was used to determine the time-averaged local void fraction and liquid velocity and turbulence intensity for air–oil slug flow. The capacitance signals along with the data obtained by image processing of the flow were used to determine the elongated bubble length and velocity. The measurements techniques used found to describe in detail the internal structure of the slug flow. Finally, the experimental results were compared to existing models and correlations.     

Frequency and development of slugs in a horizontal pipe at large liquid flows

The generation of slugs was studied for air–water flow in horizontal 0.0763 m and 0.095 m pipes. The emphasis was on high liquid rates (u_LS ? 0.5 m/s) for which slugs are formed close to the entry and the time intervals between slugs are stochastic. A “fully developed” slug flow is defined as consisting of slugs with different sizes interspersed in a stratified flow with a height slightly larger than the height, h₀, needed for a slug to be stable. Properties of this “fully developed” pattern are discussed. A correlation for the frequency of slugging is suggested, which describes our data as well as the data from other laboratories for a wide range of conditions. The possibility is explored that there is a further increase of slug length beyond the “fully developed” condition because slugs slowly overtake one another.     

Void fraction and pressure waves in a transient horizontal slug flow

The void fraction and the pressure waves in an air–water mixture flowing in the slug regime are experimentally investigated in a horizontal line. The test section is made of a transparent Plexiglas pipe with 26 mm ID and 26.24 m long, operating at ambient temperature and pressure. The flow induced transients are made by quickly changing the air or the water inlet velocity. The test grid has four operational points. This choice allows one to create expansion and compression waves due to the changes to the gas or to the liquid. Each experimental run is repeated 100 times to extract an ensemble average capable of filtering out the intrinsic flow intermittence and disclosing the void fraction and pressure waves’ features. The slug flow properties such as the bubble nose translational velocity, the lengths of liquid film underneath the bubble and the liquid slug are also measured. The objective of the work is two-fold: access the main characteristics of the void fraction and pressure waves and disclose the mechanics of the transient slug flow as described through the changes of the slug flow properties.     

The onset of gas entrainment from a flowing stratified gas–liquid regime in dual discharging branches: Part II: Critical conditions at low to moderate branch Froude numbers

This is the second of a two part investigation. Experiments were performed in a 50.8 mm diameter horizontal pipe with three 6.35 mm diameter branches located at the test section mid-span. The inlet length was 1.8 m, and three branch orientations were tested at 0° (side), 45° (inclined), and 90° (bottom) from horizontal. Water and air, operating at 206 kPa, were used and both fluids flowed co-currently within the inlet in the smooth-stratified regime. The inlet superficial velocity of the liquid phase ranged between 0.04 and 0.15 m/s while in the gas phase values of 0.3, 0.4, and 1 m/s were tested. Three different dual discharge combinations were tested and included side-inclined, side-bottom, and inclined-bottom. The tested branch flow Froude numbers were limited between low to moderate values which ranged between 1 and 23. Extensive experimental data are reported for the critical conditions at the onset of gas entrainment during dual discharge. A novel map was developed for the inclined-bottom branch configuration showing the relationship between the inlet superficial liquid velocity and branch Froude numbers. This map was used to quantify the three observed modes of gas entrainment during dual discharge. These modes were classified as onset of gas entrainment in the inclined branch only, in the bottom branch only, or both the inclined and bottom branches simultaneously. The critical height at the onset of gas entrainment results are compared to published models and data sets and poor agreement was found with studies conducted in stratified gas–liquid reservoirs.     

垂直上升管内油—气—水三相弹状流压力降的研究

给出一种垂直上升油－气－水三相弹状流压力降的计算模型。该模型考虑弹状流中Ｔａｙｌｏｒ气泡周围下降液膜的变化历程。通过油－气－水弹状流的实验研究发现，该模型的数值模拟结果与低压工况下的实验值符合得较好。本模型是计算垂直油－气－水三相弹状流中液相的连续相为水相时的压力降的有效方法。     

The formation of hydrodynamic slugs by the interaction of waves in gas–liquid two-phase pipe flow

This paper examines the effects of wave interaction on the formation of hydrodynamic slugs in two-phase pipe flow at relatively low gas and liquid superficial velocities. The experiments were conducted using a horizontal 31 m long, D = 10 cm internal diameter transparent pipe at atmospheric pressure. High resolution photography allowed the location of the gas–liquid interface to be measured with a high degree of accuracy at 5 Hz. Image analysis allowed individual waves to be tracked over a 14D section of the pipe. Regular waves having similar properties such as speed, amplitude and length were seen far from the region of slug formation. However, near the transition region, where hydrodynamic slugs were formed, significant differences between wave properties were observed which resulted in wave interaction leading to a type of sub-harmonic resonance and slug formation. The formation of hydrodynamic slugs due to wave interaction differs from predictions for slug formation using long wavelength stability theory. The properties of the waves were quantified which gave detailed information on the resonance mechanism found near the transition to slugging.     

Non-newtonian liquid-air stratified flow through horizontal tubes—ii

Non-Newtonian liquid gas stratified flow data were obtained using 0.052 and 0.025 m dia horizontal circular ducts. Unless the liquid velocity was very low, the flow pattern generally observed was non-uniform stratified flow having an interfacial level gradient between the two phases. The Heywood-Charles model is valid for predicting the pressure drop and liquid holdup in pseudoplastic (shear thinning) non-Newtonian liquid-gas uniform stratified flow. Two-phase drag reduction, which is predicted by the Heywood-Charles model did not occur because there was a transition to semi-slug flow before the model criteria were reached. Interfacial liquid and gas shear stresses were compared.