Numerical simulation of terrain-induced severe slugging coupled by hydrodynamic slugs in a pipeline–riser system

A numerical study based on a one-dimensional two-fluid model is carried out to describe the transient hydrodynamic slugging and terrain-induced severe slugging in a pipeline–riser system. The system of equations is rendered well-posed by interfacial pressure model for the riser. The selected flow conditions are restricted in the well-posed region for the horizontal and the downward inclined pipes to ensure the hydrodynamic slug characteristics are predicted correctly. The validity of the model is examined by water faucet problem and horizontal slug flow experiments. Simulations with and without slug capturing are conducted to address the effect of hydrodynamic slugs on severe slugging. It has been found that more accurate predictions are obtained by taking hydrodynamic slugs into account. At low superficial gas velocity, the simulation without slug capturing tends to overestimate the severe slugging period. When hydrodynamic slugs are captured, the upstream gas expansion is suppressed by the hydrodynamic slugs. At relatively high superficial gas velocity, the simulation without slug capturing tends to underestimate the severe slugging period. When hydrodynamic slugs are captured, the upstream compressible volume is greatly enlarged by the blowout of the hydrodynamic slugs. In both situations, the influences of the hydrodynamic slugs can reduce the errors of the predicted severe slugging characteristics.     

Numerical modeling of slug flow initiation in a horizontal channels using a two-fluid model

This paper presents a methodology for modeling slug initiation and growth in horizontal ducts. Transient two-fluid equations are solved numerically using a class of high-resolution shock capturing methods. The advantage of this method is that slug formation and growth in a stratified regime can be calculated directly from the solutions to the flow field differential equations. In addition, by using high-resolution shock capturing methods that do not contain numerical diffusion, the discontinuity generated by slugging in the flow field can be modeled with good accuracy. The two-fluid model is shown to be well-posed mathematically only under certain conditions. Under these circumstances, the two-fluid model is capable of correctly predicting and modeling the flow physics. When ill-posed, an unbounded instability occurs in the flow field solution, and the instability amplitude increases exponentially with decreasing mesh sizes. This work shows that there are three zones associated with slug formation. In addition, long wavelength slugs are shown to initiate from short wavelength waves. These short waves are generated at the interface of the two phases by the Kelvin-Helmholtz hydrodynamic instability. The results obtained through numerical modeling show good agreement with experimental results.     

Numerical evaluation of turbulence models for dense to dilute gas–solid flows in vertical conveyor

A two-fluid model (TFM) of multiphase flows based on the kinetic theory and small frictional limit boundary condition of granular flow was used to study the behavior of dense to dilute gas–solid flows in vertical pneumatic conveyor. An axisymmetric 2-dimensional, vertical pipe with 5.6 m length and 0.01 m internal diameter was chosen as the computation domain, same to that used for experimentation in the literature. The chosen particles are spherical, of diameter 1.91 mm and density 2500 kg/m³. Turbulence interaction between the gas and particle phases was investigated by Simonin's and Ahmadi's models and their numerical results were validated for dilute to dense conveying of particles. Flow regimes transition and pressure drop were predicted. Voidage and velocity profiles of each phase were calculated in radial direction at different lengths of the conveying pipe. It was found that the voidage has a minimum, and gas and solid velocities have maximum values along the center line of the conveying pipe and pressure drop has a minimum value in transition from dense slugging to dilute stable flow regime. Slug length and pressure fluctuation reduction were predicted with increasing gas velocity, too. It is shown that solid phase turbulence plays a significant role in numerical prediction of hydrodynamics of conveyor and the capability of particles turbulence models depends on tuning parameters of slip-wall boundary condition.     

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.     

Positive frictional pressure gradient in vertical gas-high viscosity oil slug flow

The present study describes the wall shear stress and the falling liquid film behavior in upward vertical slug flow of air and high viscosity oil. The frictional pressure gradient is directly related to the wall shear stress, and it is usually negative (opposite to the overall flow direction). However, in vertical slug flow, the average total wall shear stress of a slug unit may be negative (in the same direction of the overall flow), resulting in a positive frictional pressure gradient. However, this does not mean, by any way, generation of additional energy or violation of the second law of thermodynamics.The positive frictional pressure gradient phenomenon, reasons and required conditions were explained in this paper. A simplified model was developed and validated against recent experimental data of air-high viscosity oil slug flow in a 50.8 mm ID vertical pipe. The oil viscosity was in the range of 127 mPa s to 580 mPa s. Positive frictional pressure gradient appears when the liquid film wall shear stress supersede the wall shear stress in the slug body. The rate of increase of both wall shear stresses (with respect to the mixture Reynolds number) depend, not only, on the mixture Reynolds number but also, highly, on the liquid viscosity.     

Computational investigation of hydrodynamics and cracking reaction in a heavy oil riser reactor

This paper presents a computational investigation of hydrodynamics, heat transfer and cracking reaction in a heavy oil riser operated in a novel operating mode of low temperature contact and high catalyst-to-oil ratio. Through incorporating feedstock vaporization and a 12-lump cracking kinetics model, a validated gas–solid flow model has been extended to the analysis of the hydrodynamic and reaction behavior in an industrial riser. The results indicate that the hydrodynamics, temperature and species concentration exhibit significantly nonuniform behavior inside the riser, especially in the atomization nozzle region. The lump concentration profiles along the riser height provide useful information for riser optimization. Compared to conventional fluid catalytic cracking (FCC) process, feedstock conversion and gasoline yield are respectively increased by 1.9 units and 1.0 unit in the new FCC process, the yield of liquefied petroleum gas is increased by about 1.0 unit while dry gas yield is reduced by about 0.3 unit.     

The elimination of severe slugging-experiments and modeling

Severe slugging can occur in a pipeline-riser system operating at low liquid and gas rates. The flow of gas into the riser can be blocked by liquid accumulation at the base of the riser. This can cause formation of liquid slugs of a length equal to or longer than the height of the riser. A cyclic process results in which a period of no liquid production into the separator occurs, followed by a period of very high liquid production. This study is an experimental and theoretical investigation of two methods for eliminating this undesirable phenomenon, using choking and gas lift. Choking was found to effectively eliminate or reduce the severity of the slugging. However, the system pressure might increase to some extent. Gas lift can also eliminate severe slugging. While choking reduces the velocities in the riser, gas lift increases the velocities, approaching annular flow. It was found that a relatively large amount of gas was needed before gas injection would completely stabilize the flow through the riser. However, gas injection reduces the slug length and cycle time, causing a more continuous production and a lower system pressure. Theoretical models for the elimination of severe slugging by gas lift and choking have been developed. The models enable the prediction of the flow behavior in the riser. One model is capable of predicting the unstable flow conditions for severe slugging based on a static force balance. The second method is a simplified transient model based on the assumption of a quasi-equilibrium force balance. This model can be used to estimate the characteristics of the flow, such as slug length and cycle time. The models were tested against new severe slugging data acquired in this study. An excellent agreement between the experimental data and the theoretical models was found.     

Simulation of large amplitude waves in a slug tracking scheme compared to roll wave experiments at high pressure

Due to the similarities between large amplitude roll waves and slug flow in two-phase gas–liquid pipe flow, a slug tracking scheme is presented with the addition of a simplified model for roll waves. The waves are treated in a similar way to slugs, modelled as objects moving at the wave velocity and with a pressure variation across them. The two-fluid model is solved on a stationary staggered grid in stratified sections between moving waves and slugs. The model is dynamic meaning that the growth and decay of waves and slugs can be simulated. The wave model implementation within the tracking scheme is discussed and demonstrated in comparison to existing experimental data on wave velocities and averaged pressure drops. The results from the tracking scheme compared well to the experiments when waves were initiated with the experimental frequency. Wave initiation remains as a modelling challenge.     

Challenges in slug modeling and control for offshore oil and gas productions: A review study

The upstream offshore multi-phase well-pipeline-riser installations are facing huge challenges related to slugging flow: An unstable flow regime where the flow rates, pressures and temperatures oscillate in the multi-phase pipelines. One typical severe slug is induced by vertical wells or risers causing the pressure to build up and hence originates the oscillating pressure and flow. There exist many negative consequences related to the severe slugging flow and thus lots of investments and effort have been put into reducing or completely eliminating the severe slug. This paper reviews in details the state-of-the-art related to analysis, detection, dynamical modeling and elimination of the slug within the offshore oil & gas Exploration and Production (E&P) processes. Modeling of slugging flow has been used to investigate the slug characteristics and for design of anti-slug control as well, however most models require specific facility and operating data which, unfortunately, often is not available from most offshore installations. Anti-slug control have been investigated for several decades in oil & gas industry, but many of these existing methods suffer the consequent risk of simultaneously reducing the oil & gas production. This paper concludes that slug is a well defined phenomenon, but even though it has been investigated for several decades the current anti-slug control methods still have problems related to robustness. It is predicted that slug-induced challenges will be even more severe as a consequence of the longer vertical risers caused by deep-water E&P in the future.     

Developing structure of two-phase flow in a large diameter pipe at low liquid flow rate

In order to develop the interfacial area transport equation for the interfacial transfer terms in the two-fluid model, accurate data sets on axial development of local parameters such as void fraction, interfacial area concentration, interfacial gas velocity and Sauter mean diameter are indispensable to verify the modeled source and sink terms in the interfacial area transport equation. From this point of view, local measurements of both group 1 spherical/distorted bubbles and group 2 cap/slug bubbles in vertical upward air–water two-phase flow in a large diameter pipe with 200 mm in inner diameter and 26 m in height were performed at three axial locations of z/D = 41.5, 82.8 and 113 as well as 11 radial locations from r/R = 0–0.95 by using four-sensor probe method. Here, z, r, D and R are the axial distance from the inlet, radial distance from the pipe center, pipe diameter and pipe radius, respectively. The liquid flow rate and the void fraction ranged from 0.0505 m/s to 0.312 m/s and from 1.98% to 32.6%, respectively in the present experiment. The flow condition covered extensive region of bubbly flow, cap turbulent flow as well as their transition. The extensive analysis on the radial profiles of local flow parameters and their axial developments demonstrate the development of interfacial structures along the flow direction due to the bubble coalescence and breakup and the gas expansion. The significant decrease in void faction and interfacial area concentration and the increase in Sauter mean diameter and interfacial velocity were observed when the gradual flow regime transition occurred. Finally, the net change in the interfacial area concentration due to the bubble coalescence and breakup was quantitatively investigated in the present paper to reflect the true transfer mechanisms in observed two-phase flows.     

CFD study: Effect of pulsating flow on gas–solid hydrodynamics in FCC riser

Gas–solid flow in a fluid catalytic cracking (FCC) riser exhibits poor mixing in the form of a core–annulus flow pattern and a dense bottom/dilute top distribution of solids. To enhance gas–solid mixing, studies on dense fluidized beds have suggested using a pulsating flow of gas. The present study investigates the effect of pulsating flow on gas–solid hydrodynamics inside the FCC riser employing computational fluid dynamics. Two flow conditions are investigated: a cold flow of air-FCC catalyst in a pilot-scale riser and a reactive flow in an industrial-scale FCC riser. In the cold-flow riser, pulsating flows cause the slug flow of solids and thus increase the average solid accumulation in the flow domain and solid segregation towards the wall. In the industrial FCC riser, pulsating flows produce radial profiles that are more homogeneous. Pulsating flows further improve the conversion and yield in the initial few metres of height. At 7 m, the conversion from pulsating flow is 59%, compared with 44% in without pulsating flow. The results and analysis presented here will help optimize flow conditions in the circulating fluidized bed riser, in not only FCC but also applications such as fast pyrolysis and combustion.     

Fluid dynamics of horizontal air–water slug flows through a dividing T-junction

Experimental data from horizontal air–water slug flows were obtained in a test facility which was a 34 mm internal diameter, 10 m long Plexiglas pipe connected to the 90° branch arms from a T-junction. The test points were located on the flow pattern map in the proximity of the transition lines which separates different flow patterns. Capacitive probes with helical and concave plate sensors were used to quantify the dynamic liquid holdup in each branch. They were combined with Venturi nozzles + differential pressure transmitters in each outlet branch for measuring the two-phase mass flow rates. The dynamic characteristics of the slug flow splitting in a T-junction were studied from the acquired signals. Diaphragm straight-through type valves were used in the run and in the lateral branch arms to imitate equipments consuming the two-phase flow after the T-junction. This assembly can also be used as a gas–liquid separation system. The results showed different mechanisms acting on the slug flow division phenomenon. Liquid accumulation into the run branch, between the TJ and the control valve, caused more gas to come to the lateral branch.     

Evaluation of subgrid-scale models in large-eddy simulation of flow past a two-dimensional block

Large-eddy simulations of flow past a two-dimensional (2D) block were performed to evaluate four subgrid-scale (SGS) models: (i) the traditional Smagorinsky model, (ii) the Lagrangian dynamic model, (iii) the Lagrangian scale-dependent dynamic model, and (iv) the modulated gradient model. An immersed boundary method was employed to simulate the 2D block boundaries on a uniform Cartesian grid. The sensitivity of the simulation results to grid refinement was investigated by using four different grid resolutions. The velocity streamlines and the vertical profiles of the mean velocities and variances were compared with experimental results. The modulated gradient model shows the best overall agreement with the experimental results among the four SGS models. In particular, the flow recirculation, the reattachment position and the vertical profiles are accurately reproduced with a relative coarse grid resolution of (N_x × N_y × N_z=) 160 × 40 × 160 (n_x × n_z = 13 × 16 covering the block). Besides the modulated gradient model, the Lagrangian scale-dependent dynamic model is also able to give reasonable prediction of the flow statistics with some discrepancies compared with the experimental results. Relatively poor performance by the Lagrangian dynamic model and the Smagorinsky model is observed, with simulated recirculating patterns that differ from the measured ones. Analysis of the turbulence kinetic energy (TKE) budget in this flow shows evidence of a strong production of TKE in the shear layer that forms as the flow is deflected around the block.     

Void fraction measurement of bubble and slug flow in a small channel using the multivision technique

Using the multivision technique, a new void fraction measurement method was developed for bubble and slug flow in a small channel. The multivision system was developed to obtain images of the two-phase flow in two perpendicular directions. The obtained images were processed—using image segmentation, image subtraction, Canny edge detection, binarization, and hole filling—to extract the phase boundaries and information about the bubble or slug parameters. With the extracted information, a new void fraction measurement model was developed and used to determine the void fraction of the two-phase flow. The proposed method was validated experimentally in horizontal and vertical channels with different inner diameters of 2.1, 2.9, and 4.0 mm. The proposed method of measuring the void fraction has better performance than the methods that use images acquired in only one direction, with a maximum absolute difference between the measured and reference values of less than 6%.     

Quantifying the intensity of vehicle impacts using vehicle tracking devices during live training exercises

The operation of off-road vehicles during military training exercises can affect the environmental conditions of training lands by removing or disturbing vegetation. To quantify the impact of vehicle based military training, global positioning system (GPS)-based vehicle tracking systems were used to characterize the movement of vehicles during live training exercises. Methods were developed to spatially estimate the tracking intensity (number of vehicle passes per area) resulting from the training exercises. This method was then combined with previous developed methods that identified off-road trail formation and vehicle dynamic properties to quantify the overall training mission impacts of specific training events on installation resources. This approach to characterizing training impacts results in mission impact profiles that more accurately quantify live training mission impacts.Search radius and output grid size are important parameters of the proposed traffic intensity approximation method. Traffic intensities estimated using a variety of search radii and grid sizes were compared. Results indicated that a 10 m search radius and a 10-by-10 m output grid size worked the best for the study dataset. Approximately, 89% accuracy was found for traffic intensity (number of passes) estimation when using a 10 m search radius and a 10-by-10 m output grid size.     

Rotational asymmetry of reactant concentration and its evolution in a circulating fluidized bed riser

Rotational asymmetric distribution of reactant (ozone) concentration and its evolution along with the gas-solid reactive flow were studied in a 76 mm i.d., 10.2 m high circulating fluidized bed (CFB) riser reactor. The superficial gas velocity ranged from 3 to 5 m/s and the solids circulation rates were 50 and 100 kg/(m² s). Experimental results show that the asymmetry of reactant distribution can extend to a height close to the length of flow developing zone of the CFB riser reactor and then disappears. Based on the hydrodynamics of the gas and solid phases in the solids entrance region, this asymmetry can be attributed to the effect of the solids entrance structure.     

Numerical study on mitigating severe slugging in pipeline/riser system with wavy pipe

The gas/liquid two-phase flow in pipeline/wavy-pipe/riser systems was investigated numerically with CFD. A CFD model of the pipeline/wavy-pipe/riser system was obtained by adding a wavy pipe to the model of the pipeline/riser system verified by the experimental data previously. The effects of the geometrical parameters and location of the wavy pipe on its performance of slug mitigation and flow characteristics in pipeline/wavy-pipe/riser systems were examined through the CFD models. With the increase of the amplitude or length of the wavy pipe, the slug in the pipeline/riser system becomes shorter. The optimum location of the wavy pipe in the pipeline exists for a pipeline/riser system and a wavy pipe at given operating conditions. The CFD modelling provides a feasible and flexible way to investigate the effectiveness of the wavy pipes on mitigating severe slugging in pipeline/riser systems.     

Flow of mono-dispersed particles through horizontal bend

Pipeline slurry flow of mono-dispersed particles through horizontal bend is numerically simulated by implementing Eulerian two-phase model in FLUENT software. A hexagonal shape and Cooper type non-uniform three-dimensional grid is chosen to discretize the entire computational domain, and a control volume finite difference method is used to solve the governing equations. The modeling results are compared with the experimental data collected in 53.0 mm diameter horizontal bend with radius of 148.4 mm for concentration profiles and pressure drops. Experiments are performed on narrow-sized silica sand with mean diameter of 450 μm and for flow velocity up to 3.56 m/s (namely, 1.78, 2.67 and 3.56 m/s) and four efflux concentrations up to 16.28% (namely, 0%, 3.94%, 8.82% and 16.28%) by volume for each velocity. Eulerian model gives fairly accurate predictions for both the pressure drop and concentration profiles at all efflux concentrations and flow velocities.     

CFD based mini- vs. micro-system delineation in elongated bubble flow regime

Delineation of mini- and micro-scale channels with respect to two-phase flow has been the subject of many research papers. There is no consensus on when the small channel can be characterized as a mini-channel or micro-channel. The idea proposed by this paper is to use the normalized bubble nose radius, liquid film thickness top over bottom ratio, and bubble shape contour, which are found under normal gravity conditions in slug flow through a horizontal adiabatic channel, as the delineation criteria. The input parameters are bubble nose radius and bubble nose velocity as the characteristic length scale and characteristic velocity scale respectively. 3D numerical simulation with ANSYS FLUENT was used to obtain the necessary data. Following CFD practice, a mesh independence study and a numerical model validation against published experimental data were both conducted. Analysis of the numerical simulation results showed that channels with D ⩽ 100 μm can be characterized as a micro-system, while channels with D ⩾ 400 μm belong to mini-systems. The region 200 μm ⩽ D ⩽ 300 μm represents a transition from the micro-scale to mini-scale.     

A 5‐equation,transient, hyperbolic, 1‐dimensional model for slug capturing in pipes

A novel numerical scheme for slug capturing in pipes using a 1‐dimensional transient hyperbolic 5‐equation 2‐fluid model is presented. Previous work has shown that 1‐dimensional 2‐fluid models are able to capture slug flow automatically. In this work, a similar approach is further developed using a new numerical scheme, applied to a hyperbolic 5‐equation 2‐fluid model. Starting from a finite volume discretisation of a 5‐equation 2‐fluid hyperbolic model and adding appropriate closure relations, a second‐order code is implemented and applied to air‐water flows in horizontal pipes, simulating the 2‐phase to 1‐phase flow process. The code is evaluated in some common standard test cases. A slug capturing application is also discussed. We show, in an air/water horizontal pipe, slug initiation, growth, and development. Moreover, a grid refinement analysis is performed showing that the method is grid independent and we show the code capability to take into account eventual surface tension effects, through the instantaneous pressure relaxation process. Finally, a prediction of flow regime transitions is shown and compared with a well‐known theoretical flow pattern map in addition to a preliminary comparison of computed slug characteristics against well‐known empirical correlations.