Search

Two-phase flow with low liquid loads is common in high-pressure natural gas offshore gathering and transmission pipelines. During gas production slowdowns or shutdowns, an accumulation of liquid in the lower sections of subsea pipelines may occur. This phenomenon is observed in jumpers that connect different units in deep-water subsea gas production facilities. The displacement of the accumulated liquid during production ramp-up induces temporal variations in pressure drop across the jumper and forces on its elbows, resulting in flow-induced vibrations (FIV) that pose potential risks to the structural integrity of the jumper. To bridge the gap between laboratory experiments and field conditions, transient 3D numerical simulations were conducted using the OpenFOAM software. These simulations facilitated the development of a mechanistic model to elucidate the factors contributing to increased pressure and forces during the liquid purging process. The study examined the influence of gas pressure level, pipe diameter, initially accumulated liquid amount, liquid properties, and gas mass flow rate on the transient pressure drop and the forces acting on the jumper's elbows. The critical gas production rate required for complete liquid removal of the accumulated liquid was determined, and scaling rules were proposed to predict the effects of gas pressure and pipe diameter on this critical value. The dominant frequencies of pressure and force fluctuations were identified, with low-pressure systems exhibiting frequencies associated with two-phase flow phenomena and high-pressure systems showing frequencies attributed to acoustic waves.

This work deals with stability of two-phase stratified air-water flows in horizontal circular pipes. For this purpose, we performed a linear stability analysis, which considers all possible three-dimensional infinitesimal disturbances and takes into account deformations of the air-water interface. The main results are presented in form of stability maps, which compare well with the available experimental data. The neutral stability curves are accompanied by the corresponding wavenumbers and wave speeds of the critical perturbations, as well as by spatial patterns of their velocity components. Accordingly, several modes of the critical perturbation are revealed. Long waves are found to be the critical perturbation over part of the stability boundary, and they are affected by the surface tension due to the confinement effect of the lateral direction. Exploring the effect of pipe diameter on the stability boundary and critical perturbations shows that for small water holdups (i.e., thin water film) the scaling of the critical gas velocity by the gas Froude number is valid for pipe diameters larger than about 0.1m, where the surface tension effects due to the lateral confinement become negligible. Comparing results obtained in pipe, square duct, and two-plate geometries, we show that there are cases where the simplified geometry of two parallel plates can be employed to model the realistic geometry reasonably well., Comment: 19 pages, 20 figures

We address the problem of steady laminar flow of a shear-thinning fluid in rectangular ducts, which is encountered in many systems, in particular, in microfluidic and biomedical devices. However, an exact solution for the flow of non-Newtonian fluids that considers a realistic shear-thinning rheological behavior is not available in the literature. In this study, an accurate solution for the case of Carreau fluid is obtained and investigated numerically. The numerical solution allows us to analyze the effects of the fluid rheology and the aspect ratio of the rectangular duct on the velocity field and pressure gradient that drives the flow. The relationship between the pressure gradient and the Carreau number is found to follow the rheological curve of the shear-thinning fluid. The analysis shows that the fluid rheology and the aspect ratio have independent contributions to the integral flow characteristics. Moreover, separate consideration of these contributions allows us to arrive at universal scaling and general formulae for the pressure gradient and friction factor for various rheological parameters of the fluid and aspect ratios of rectangular ducts., Comment: 13 pages, 9 figures

Removal of the accumulated liquid from jumpers of subsea gas production systems is essential to avoid possible hydrate creation and further damage to the pipeline. However, the displacement of high amounts of accumulated liquid during the production start-up leads to a gas pressure rise. Liquid plugs formed during the liquid displacement impact the structure's elbows. This in addition to cyclic pressure/forces fluctuations, may lead to harmful flow-induced vibrations (FIV). These flow phenomena that may endanger the jumper structure were explored in air-water experiments performed in a lab-scale jumper. The critical (minimal) gas velocity needed to purge the accumulated liquid was determined and the pressure and forces variations during the liquid removal were measured. In addition, the effects of the gas velocity, initial liquid amount, and gas flow ramp-up on the air-water flow phenomena were documented. Results of 3D and 2D numerical simulations (using OpenFOAM) were verified against the experimental data. The effects of employing different RANS turbulence models on the predictions were tested and demonstrated. A simple mechanistic model was established to predict the pressure and force variation during liquid displacement. The model enables inspecting the variation with the operational conditions of each pressure component (i.e., hydrostatic, friction, and acceleration) and examining their significance and contribution to the pressure rise.

A numerical framework for rigorous linear stability analysis of two-phase stratified flows of two immiscible fluids in horizontal circular pipes is presented. For the first time, three-dimensional disturbances, including those at the interface between two fluids, are considered. The proposed numerical framework is based on a finite-volume method and allows solving the problem numerically in bipolar cylindrical coordinates. In these coordinates, both the pipe wall and the unperturbed interface (of a constant curvature, e.g., plane interface, as considered in this work) coincide with the coordinate surfaces. Thereby, the no-slip as well as the interfacial boundary conditions can be imposed easily. It also enables investigation of the local behavior of the flow field and shear stresses in the vicinity of the triple points, where the interface contacts the pipe wall. The results obtained in the bipolar coordinates are verified by an independent numerical solution based on the problem formulation in Cartesian coordinates, where the pipe wall is treated by the immersed boundary method. Two representative examples of gas-liquid and liquid-liquid flows are included to demonstrate the applicability of the proposed numerical technique for analyzing the flow stability., Comment: 37 pages, 11 figures

The linear stability of stratified two-phase flows in rectangular ducts is studied numerically. The linear stability analysis takes into account all possible infinitesimal three-dimensional disturbances and is carried out by solution of the associated eigenproblem. The neutral stability boundary and the corresponding critical wave number are obtained for liquid - liquid and air - water systems. Depending on the problem parameters, the instability sets in owing to short, intermediate, of long wave most unstable perturbations. Patterns of the most unstable disturbances are reported and discussed. It is shown that the instability arises due to shear, or interfacial mechanisms. Effects of the surface tension and of width/height aspect ratio are also studied. The results support the premise that the stability analysis of stratified two-phase flow in the simpler geometry of two-infinite plates can provide a reasonable estimation of the conditions for which this flow pattern can be considered to be linearly stable.

The non-modal transient growth of perturbations in horizontal and inclined channel flows of two immiscible fluids is studied. 3D perturbations are examined in order to find the optimal perturbations that attain the maximum amplification of perturbation energy at relatively short times. Definition of the energy norm is extended to account for the gravitational potential energy along with the kinetic energy and interfacial capillary energy. Contrarily to the fastest exponential growth, which is reached by essentially 2D perturbations, the maximal non-modal energy growth is attained mostly by three-dimensional spanwise perturbations. Significant transient energy growth is found to occur in linearly stable flow configurations, which, similarly to single phase shear flows, may trigger non-linear destabilizing mechanisms within one of the phases. It is shown that the transient energy growth in linearly stable cases can be accompanied by noticeable interface deformations. Therefore, flow pattern transition due to non-modal transient growth and reduction of the range of operational conditions for which stratified-smooth flow remains stable cannot be ruled out.

Linear stability of horizontal and inclined stratified channel flows of Newtonian/non-Newtonian shear-thinning fluids is investigated with respect to all wavelength perturbations. The Carreau model has been chosen for the modeling of the rheology of a shear-thinning fluid, owing to its capability to describe properly the constant viscosity limits (Newtonian behavior) at low and high shear rates. The results are presented in the form of stability boundaries on flow pattern maps (with the phases' superficial velocities as coordinates) for several practically important gas-liquid and liquid-liquid systems. The stability maps are accompanied by spatial profiles of the critical perturbations, along with the distributions of the effective and tangent viscosities in the non-Newtonian layer, to show the influence of the complex rheological behavior of shear-thinning liquids on the mechanisms responsible for triggering instability. Due to the complexity of the considered problem, a working methodology is proposed to alleviate the search for the stability boundary. Implementation of the proposed methodology helps to reveal that in many cases the investigation of the simpler Newtonian problem is sufficient for the prediction of the exact (non-Newtonian) stability boundary of smooth stratified flow (i.e., in case of horizontal gas-liquid flow). Therefore, the knowledge gained from the stability analysis of Newtonian fluids is applicable to those (usually highly viscous) non-Newtonian systems. Since the stability of stratified flow involving highly viscous Newtonian liquids has not been researched in the literature, interesting findings on the viscosity effects are also obtained.

Exact solutions for laminar stratified flows of Newtonian/non-Newtonian shear-thinning fluids in horizontal and inclined channels are presented. An iterative algorithm is proposed to compute the laminar solution for the general case of a Carreau non-Newtonian fluid. The exact solution is used to study the effect of the rheology of the shear-thinning liquid on two-phase flow characteristics considering both gas/liquid and liquid/liquid systems. Concurrent and counter-current inclined systems are investigated, including the mapping of multiple solution boundaries. Aspects relevant to practical applications are discussed, such as the insitu hold-up, or lubrication effects achieved by adding a less viscous phase. A characteristic of this family of systems is that, even if the liquid has a complex rheology (Carreau fluid), the two-phase stratified flow can behave like the liquid is Newtonian for a wide range of operational conditions. The capability of the two-fluid model to yield satisfactory predictions in the presence of shear-thinning liquids is tested, and an algorithm is proposed to a priori predict if the Newtonian (zero shear rate viscosity) behaviour arises for a given operational conditions in order to avoid large errors in the predictions of flow characteristics when the power-law is considered for modelling the shear-thinning behaviour. Two-fluid model closures implied by the exact solution and the effect of a turbulent gas layer are also addressed., Comment: 36 pages, 27 Figures

Linear stability of stratified gas-liquid and liquid-liquid plane-parallel flows in inclined channels is studied with respect to all wavenumber perturbations. The main objective is to predict parameter regions in which stable stratified configuration in inclined channels exists. Up to three distinct base states with different holdups exist in inclined flows, so that the stability analysis has to be carried out for each branch separately. Special attention is paid to the multiple solution regions to reveal the feasibility of non-unique stable stratified configurations in inclined channels. The stability boundaries of each branch of steady state solutions are presented on the flow pattern map and are accompanied by critical wavenumbers and spatial profiles of the most unstable perturbations. Instabilities of different nature are visualized by streamlines of the neutrally stable perturbed flows, consisting of the critical perturbation superimposed on the base flow. The present analysis confirms the existence of two stable stratified flow configurations in a region of low flow rates in countercurrent liquid-liquid flows. These configurations become unstable with respect to shear mode of instability. It was revealed that in slightly upward inclined flows the lower and middle solutions for the holdup are stable in a part of the triple solution region, while the upper solution is always unstable. In the case of downward flows, in the triple solution region, none of the solutions are stable with respect to short-wave perturbations. These flows are stable only in the single solution region at low flow rates of the heavy phase, where long-wave perturbations are the most unstable ones.

Linear stability of stratified two-phase flows in horizontal channels to arbitrary wavenumber disturbances is studied. The problem is reduced to Orr-Sommerfeld equations for the stream function disturbances, defined in each sublayer and coupled via boundary conditions that account also for possible interface deformation and capillary forces. Applying the Chebyshev collocation method, the equations and interface boundary conditions are reduced to the generalized eigenvalue problems solved by standard means of numerical linear algebra for the entire spectrum of eigenvalues and the associated eigenvectors. Some additional conclusions concerning the instability nature are derived from the most unstable perturbation patterns. The results are summarized in the form of stability maps showing the operational conditions at which a stratified-smooth flow pattern is stable. It is found that for gas-liquid and liquid-liquid systems the stratified flow with a smooth interface is stable only in confined zone of relatively low flow rates, which is in agreement with experiments, but is not predicted by long-wave analysis. Depending on the flow conditions, the critical perturbations can originate mainly at the interface (so-called "interfacial modes of instability") or in the bulk of one of the phases (i.e., "shear modes"). The present analysis revealed that there is no definite correlation between the type of instability and the perturbation wavelength.