Two-phase Flow Research Articles

Pore-scale fluid distribution and remaining oil during tertiary low-salinity waterflooding in a carbonate

The utilization of low-salinity waterflooding as a promising enhanced oil recovery method has exhibited exciting results in various experiments conducted at different scales. For carbonate rock, pore-scale understanding of the fluid distribution and remaining oil after low-salinity waterflooding is essential, especially the geometry and topology analysis of oil clusters. We performed the tertiary low-salinity waterflooding and employed X-ray micro-CT to probe the pore-scale displacement mechanism, fluid configuration, oil recovery, and remaining oil distribution. We found that the core becomes less oil-wet after low-salinity waterflooding. Furthermore, we analyzed the oil-rock and oil-brine interfacial areas to further support the wettability alteration. By comparing images after high-salinity waterflooding and low-salinity waterflooding, it is proven that wettability alteration has a significant impact on the behavior of the two-phase flow. Our research demonstrates that low-salinity waterflooding is an effective tertiary enhanced oil recovery technology in carbonate, which changes the wettability of rock and results in less film and singlet oil.

Numerical and experimental study on water-sediment flow in a lateral pumping station forebay

Multiple vortex flow patterns are commonly observed in lateral pumping station forebays, particularly in basins with high sediment concentrations. These patterns can lead to pump blockage, sediment deposition, and other issues that disrupt pump station operations. The water-sediment two-phase flow in lateral pumping station forebays is significantly influenced by the start-up combination, yet our understanding remains limited. To address this, the mixture multiphase flow theory is introduced to describe water-sediment dynamics, and the mathematical model is validated with experimental data. By analyzing the characteristics and formation mechanisms of large vortex and multiple small vortex regions in the original scheme, nine different start-up combination schemes were proposed. The research results indicate that, due to the narrow channel and slope effect in the lateral forebay, some of the gravitational potential energy of the water-sediment mixture is converted into kinetic energy upon entering the forebay, thereby increasing the velocity in the main flow area. Additionally, due to the friction and dissipation effects of the two sidewalls, a pressure difference is generated in the main flow area, resulting in the formation of multi-level vortices. Furthermore, the various types of proposed start-up combinations can optimize the flow patterns in the forebay to a certain extent. The preferred scheme improved the uniformity of flow velocity by 18.63% and increased the deviation angle by 17.62°, resulting in a 75.47% reduction in vortex area compared to the original scheme. These research results provide theoretical guidance for optimizing start-up combinations and reorganizing flow field structures to achieve hydrodynamic dredging effects.

An investigation of anisotropy in the bubbly turbulent flow via direct numerical simulations

We investigated the effects of bubble count, flow direction, and Eötvös number on deformable bubbles in turbulent channel flow. For a given shear Reynolds number Re = 180 and fixed bubble volume fractions (1.263% and 2.525%), we conducted a series of direct numerical simulations using a coupled level-set and volume-of-fluid solver to evaluate their impact on bubble volume fraction distribution, velocity fields, and turbulence characteristics. Each aspect was studied based on the microscopic equations of two-phase flow, and the accuracy of the modeling terms used in current Reynolds-averaged Navier–Stokes equation (RANS) models was assessed. The influence on the anisotropic state was analyzed using the Lumley triangle, and the anisotropy of Reynolds stresses was captured through the exact balance equations. The results indicate that in upward flow, bubbles tend to accumulate near the wall, with smaller Eötvös numbers leading to closer proximity to the wall and greater attenuation of the liquid-phase velocity. This distribution enhances energy dissipation and turbulence isotropy. In downward flow, bubbles cluster in the channel center, generating additional pseudo-turbulence and attenuating energy in the buffer layer. Moreover, the interfacial transfer of turbulent energy, as currently modeled in RANS, is found to be inadequate for upward flows.

The dynamic mechanism and optimal parameters of hydraulic jump induced disk floatation

Abstract The hydraulic jump induced disk floatation is one of the important and complicated kinetic issues. Based on the observation of the hydraulic jump phenomenon and disk floatation in experiment, the gas-liquid two-phase flow and free surface advection are simulated by COMSOL Multiphysics. The force and flow fields, as well as the movement mode of the disk are analyzed through numerical simulation. Meanwhile, a theoretical model is firstly proposed to investigate the motion equations from the perspectives of energy and force. The stable, metastable and instable motion modes are found out and analyzed in detail by variable-controlled and comparison methods. The critical parameters of stability and instability of the disk are obtained by continuous interpolation method.

A numerical model for solitary wave breaking based on the phase-field lattice Boltzmann method

This study presents a numerical investigation of a solitary wave breaking over a slope by using the phase-field lattice Boltzmann method. The incompressible two-phase flow equations are solved by using a velocity-based formulation of the two-phase lattice Boltzmann method with a central-moment collision model to accurately simulate wave breaking problems. For interface capture, a phase-field lattice Boltzmann method that ensures mass conservation is employed. The validity of the proposed method is confirmed through solitary wave propagation and transformation problems, and the obtained results are in good agreement with the experimental and calculated results. The proposed method is then employed to analyze wave breaking on a slope, demonstrating strong concordance with experimental data and existing computational findings. By analyzing the instantaneous flow characteristics and the temporal evolution of the variation in kinetic, potential, and total energy from deep to shallow water, the model can reveal the macroscopic characteristics of solitary wave breaking. Because the phase-field model effectively simulates wave breaking and air entrainment, it can depict wave energy dissipation more accurately than the single-phase lattice Boltzmann method with free surface tracking.

Direct numerical simulation of dispersion and mixing in gas–liquid Dean-Taylor flow with influence of a 90° bend

Gas-liquid capillary flow finds widespread applications in reaction engineering, owing to its ability to facilitate precise control and efficient mixing. Incorporating compact and regular design with Coiled Flow Inverter (CFI) enhances process efficiency due to improved mixing as well as heat and mass transfer leading to a narrow residence time distribution. The impact of Dean and Taylor flow phenomena on mixing and dispersion within these systems underscores their significance, but is still not yet fully understood. Direct numerical simulation based on finite element method enables full 3D resolution of the flow field and detailed examination of laminar flow profiles, providing valuable insights into flow dynamics. Notably, the deflection of flow velocity from the center axis contributes is followed by tracking of particle with defined starting positions, aiding in flow visualization and dispersion characterization. In this CFD study, the helical flow with the influence of the centrifugal force and pitch (Dean flow) as well as the capillary two-phase flow (Taylor bubble) is described and characterized by particle dispersion and related histograms. Future prospects in this field include advancements in imaging techniques to capture intricate flow patterns, as well as refined particle tracking methods to better understand complex flow behavior.

MRCNN: Multi-input residual convolution neural network for three-dimensional reconstruction of bubble flows from light field images

Accurate measurement of bubbles in air-water two-phase flows holds immense significance in the realm of thermal hydraulics assessments within nuclear reactors. Nevertheless, conventional bubble measurement techniques grapple with challenges encompassing system intricacy, limited real-time capabilities, and inaccuracies stemming from their inherent two-dimensional (2-D) nature. In response, we pioneered an innovative three-dimensional (3-D) analysis approach that leverages light field (LF) imaging diagnosis and deep learning algorithms. Unlike traditional 2-D reconstruction methods, our approach enables direct computation of bubble depth from LF images using digital refocusing technology. Following calibration, a seamless transformation is established between the camera coordinate system and the real-world coordinate system using a sharpness evaluation algorithm. This calibration process ensures precise alignment of refocused images with real-world positions. Subsequently, fully automated and highly accurate computations of bubble depth are realized from input images via the incorporation of a multi-input residual convolution neural network (MRCNN). The limitations of traditional two-dimensional imaging techniques are effectively addressed by this methodology, resulting in a reduction in measurement errors. The study confirms the feasibility of employing LF imaging diagnosis and deep learning algorithms for bubble measurements in an air-water two-phase flow, offering a significant improvement over traditional methods.

Computational fluid dynamics simulation of cryogenic vertical upflow boiling under Earth gravity

This study presents numerical simulations and their validation for flow boiling of liquid nitrogen (LN2) in a vertical upflow orientation, with a primary aim to understand the complex two-phase flow and heat transfer phenomena important to space applications. The computational fluid dynamics (CFD) model utilized the coupled level set volume of fluid (CLSVOF) method, incorporating additional source terms for bubble collision dispersion force and shear lift force in the momentum conservation equation to enhance simulation accuracy. The simulations were conducted for two mass velocities (G = 526 and 804 kg/m2s) and three different heat flux levels (approximately 10%, 30%, and 70% of critical heat flux (CHF) under Earth gravity. The model was validated against measured wall temperature data acquired from the authors’ previous experimental studies, demonstrating average deviations of less than 2.8 K across all operating conditions. The simulated two-phase flow contours illustrated various flow patterns, including bubbly, slug, churn, and annular. Both mass velocity and heat flux were observed to impact the onset of nucleate boiling (ONB), bubble nucleation, growth, and coalescence, and overall vapor structure. The simulations also offered insight into axial and radial void fraction and velocity profiles, revealing local flow acceleration trends synchronized with void fraction development. A comparison between predicted and measured bulk fluid temperature profiles showed excellent agreement, further validating the CFD model’s accuracy and practical usefulness for two-phase cryogenic flow boiling simulations in space applications.

Saturated/subcooled flow boiling heat transfer characteristics for R134a, HFE-7100 and Deionized water inside micro/mini-channels: Part I - Visualization of flow patterns and new diabatic two-phase flow pattern transition regimes

An experimental investigation on flow saturated/subcooled flow boiling heat transfer for R134a, HFE-7100 and Deionized water inside micro/mini-channels (L/dh = 40/53/70) was presented. Flow patterns and their transitions of three working fluids at various liquid subcoolings (0 ∼ 70 K) were visualized and analyzed. Experiments were carried out at a saturated pressure of 770 kPa for R134a and at atmospheric pressure for HFE-7100 and Deionized water. The mass flux varied from 400 to 800 kg/m2s, and heat fluxes up to 270.1 W/cm2 over the entire range of vapor quality (0 ∼ 1). The effects of mass flux, heat flux, liquid subcoolings and geometry parameters on flow pattern evaluation and transition were discussed. Bubble (isolated and confined bubbles), slug, churn, annular and dry-out flows were observed for the micro/mini-channels. The results indicated that the flow patterns have a significant difference on flow boiling under different liquid subcoolings. There is obvious flow pattern transition from slug flow to annular flow for nearly saturated boiling (ΔTsub ≤ 10 K). However, with the increase of inlet liquid subcooling, only obvious bubble flow and slug flow were observed during the entire boiling process, and the annular flow was not obvious before triggering the CHF. Besides, the flow pattern transition criterion is modified and the new flow pattern transition regimes including saturated/subcooled flow boiling for different working fluids are developed based on present experimental results.

Experimental investigation of different sub-regimes in horizontal stratified and intermittent gas-liquid two-phase flow: Flow map and analysis of pressure drop fluctuations

Horizontal stratified and intermittent gas-liquid two-phase flows exhibit several sub-regimes. Identifying these correctly would enable the development of more robust predictive models. This study reports on experimental investigation, based on observations and the collection of time series of pressure drop in a 40 mm ID pipe. Nine different sub-regimes (SS, 2D wave, 3D wave, RW, ED+RW, PS+RW, plug, LAS, HAS) were revealed. An original flow pattern map for this diameter is proposed including transitory regions between the sub-regime. Comparison of experimental observations with existing flow maps has highlighted the effect of pipe diameter on transition of sub-regimes. As reported in literature the sub-regimes, could be identified by direct visualization of the pressure drop time series obtained from differential pressure sensor a swell as using Probability Density Function. Furthermore, statistical analysis of standard deviation as function of gas superficial velocity enabled detection of the transition of various sub-regimes. Finally, a space feature based on standard deviation and mixture Froude number is proposed as an identification tool between the different sub-regimes investigated.

Ammonia adiabatic two-phase flow patterns in the horizontal tube

The manuscript reports the experimental study results dealing with the ammonia two-phase flow pattern observations in the transparent horizontal tube of ID 7.5 mm with a saturation temperature of 15 to 65°C, mass velocity of 50 to 160 kg·m−2s−1 and vapour quality of 0.1 to 0.8. Stratified-wavy, annular-wave and annular smooth patterns have been identified based on analyses of camera images. The annular-wavy pattern was not distinguished in the previous ammonia studies at low saturation temperatures (Lima et al., 2009; Gao et al., 2021). Transitions between stratified-wavy and annular-wavy and between annular-wavy and smooth annular patterns occur at constant interfacial vapour momentum flux. This observation aligns with Hewitt and Roberts's (1969) suggestions for vertical steam-water flows. It does not match the existing iterations of Taitel and Dukler's (1976) map, which accurately predicts the transition of HFC fluids and partly the low-temperature ammonia tests. Implementing newly developed transitions between the two-phase ammonia flow patterns in the horizontal tubes updated and validated the existing hydraulic and heat transfer correlations. The most significant among them is an update to the Müller-Steinhagen and Heck's (1986) model in predicting adequately the ammonia frictional pressure losses in the horizontal tube for the broad range of boundary conditions.

Numerical design of a flow restrictor for tanked liquid nitrogen undergoing reduced-gravity flights

NASA and the United States are designing multiple reduced-gravity cryogenic payloads to study various two-phase flow phenomena that will be ground-tested and subsequently flight-tested onboard parabolic flights. The setup will undergo 25+ parabolas per flight, during which liquid nitrogen will be cyclically or continuously demanded from a supply Dewar. The minimum gravity level during a parabolic maneuver is 0 ± 0.05 g, which can last around 20 seconds, and the maximum expected gravity is 2 g, which lasts for 40-50 seconds. This Computational Fluid Dynamics (CFD) study focuses on the design and optimization of a bi-directional perforated plate and a ring baffle to help ensure single-phase liquid nitrogen outflow during all flight phases, regardless of supply tank liquid level or gravity level. CFD analyses carried out using ANSYS FLUENT indicate that a double perforated plate with an open area percent equal or less than 0.63% can achieve high values of expulsion efficiency. Structural analyses performed using the ANSYS Static Structural solver determined the minimum plate thickness needed to withstand gravitational and hydrodynamic forces experienced during the flight.

Numerical study on flow instability in parallel helical tubes under ocean conditions

The superimposed coupling of flow oscillations generated by transient external forces under ocean conditions and flow oscillations in parallel helical tubes during two-phase flow instability may lead to more complex flow oscillation phenomena in parallel helical tubes. However, there is still a lack of understanding regarding the effect of ocean conditions on the two-phase flow instability within parallel helical tubes. The two-phase flow instability in parallel helical tubes under ocean conditions is numerically investigated in the present study. The effects of ocean conditions on the flow behavior and instability threshold in parallel helical tubes are studied with mass flux 200 kg/m2s and inlet subcooling within 5–70 °C. The results show that the mass flow rate decreases as the inclination angle increases. A larger rolling amplitude or shorter period increases the flow and temperature oscillation amplitude. Two characteristic frequencies are found in the natural circulation flow instability under the rolling motion. As the rolling amplitude increases, the compound flow oscillation amplitude is larger and more chaotic with a certain phase shift. The effect of stationary inclining on the flow instability threshold is rather significant than rolling effect. The results of the present study can serve as a significant reference for the operation and design of the helical-coiled once-through steam generator under ocean conditions.

Thermal transfer characteristics and thermoelasticity analysis of direct-steam-generation parabolic trough collector

Parabolic trough solar direct-steam-generation technology integrates clean energy with green carriers and is one of the favorable low-carbon technologies. However, the non-uniform heat flow distribution and flow stratification of collectors are susceptible to variations in climatic conditions, causing dry spots and raising the failure risk, and there are fewer studies on their influence laws based on three-dimensional two-fluid modeling. Therefore, the three-dimensional optical-thermal-stress-hydrodynamic model of the collector is built in this work based on Monte Carlo Ray Tracing method, finite volume method, finite element method and Eulerian two-fluid modeling method, and the model validation results agree with the experimental results. The impacts of direct solar radiation (I) and inlet pressures (Pin) on heat transfer performance and thermoelasticity characteristics of the three flow stages are analyzed. The results revealed that heat transfer performance in two-phase flow is better than single-phase flow at the same working conditions. The elevation of I effectively improves the Nusselt number and collector efficiency, but also dramatically raises the circumferential temperature difference (CTD) and thermal stresses. Especially, the CTD, maximum deformation, and maximum equivalent stresses in the superheated stage are as high as 54.37 K, 32.82 mm, and 38.23 MPa at Pin = 3 MPa and I = 1000 W m−2, similarly the maximum equivalent stress is also as high as 32.96 MPa in the stratified flow. The results serve as a reminder to pay particular attention to the superheated stage and the laminar flow, particularly at low Pin and high I. Furthermore, the lower the Pin, the higher the deformation and equivalent stresses in the stratified flow phase and the higher the risk of induced failure. Therefore, this implies that higher operating pressures are more appropriate for stratified flow, but this may be limited by high heat loss and low collector efficiency.

CFD-DEM analysis of internal flow characteristics in multistage deep-sea mining pumps

ABSTRACT The potential profits associated with deep-sea mineral deposits have encouraged the development of new mining techniques. In particular, the effectiveness of lift pumps is crucial to the success or failure of these mining operations. In this study, the interior pressure field, turbulence characteristics, and turbulent kinetic energy of deep-sea mining pumps were numerically simulated using the computational fluid dynamics/discrete element method (CFD-DEM). Turbulence in the solid–liquid two-phase flow in deep-sea mining lift pumps was also investigated using CFD-DEM, and the key factors affecting lift and efficiency were obtained. Data exchange among particles and between particles and fluids was fully considered to obtain reliable calculation results. The results show that lift pump operation improves when the speed increases. A real-world hydraulic test was conducted to verify the simulation result. That is, increasing the speed of the lift pump stabilises the flow of particles and fluids and gradually increases the lift and efficiency. The simulation method provides a reliable starting point for developing and optimising subsea lift pumps.

Research on the Erosion Law and Protective Measures of L360N Steel for Surface Pipelines Used in Shale Gas Extraction.

The erosion of surface pipelines induced by proppant flowback during shale gas production is significant. The surface pipelines in a shale gas field in the Sichuan Basin experienced perforation failures after only five months of service. To investigate the erosion features of L360N, coatings, and ceramics and optimize the selection of two protective materials, a gas-solid two-phase flow jet erosion experimental device was used to explore the erosion resistance of L360N, coatings, and ceramics under different impact velocities (15 m/s, 20 m/s, and 30 m/s). An energy dispersive spectroscope, a scanning electron microscope, and a laser confocal microscope were employed to analyze erosion morphologies. With the increase in flow velocity, the erosion depth and erosion rate of L360N, coating, and ceramic increased and peaked under an impact velocity of 30 m/s. The maximum erosion rate and maximum erosion depth of L360N were, respectively, 0.0350 mg/g and 37.5365 µm. Its primary material removal mechanism was the plowing of solid particles, and microcracks were distributed on the material surface under high flow velocities. The maximum erosion rate and maximum erosion depth of the coating were, respectively, 0.0217 mg/g and 18.9964 µm. The detachment of matrix caused by plowing is the main material removal mechanism. The maximum erosion rate and maximum erosion depth of ceramics were, respectively, 0.0108 mg/g and 12.4856 µm. The erosion mechanisms were micro-cutting and plowing. Under different particle impact velocities, different erosion morphologies were observed, but the primary erosion mechanism was the same. The erosion resistance of the ceramics was higher than that of the coatings. Therefore, ceramic lining materials could be used to protect the easily eroded parts, such as pipeline bends and tees, and reduce the failure rate by more than 93%. The study provided the data and theoretical basis for the theoretical study on oil and gas pipeline erosion and pipeline material selection.

A spring-mass-damper model based on separated phase flow mode for pulsating heat pipe with adjustive-structured channels

Pulsating heat pipe (PHP) is a kind of efficient passive phase-change cooling device. The pulsating behaviors of the two-phase flow inside PHP significantly affect the heat transfer performance for PHP, the investigation of which will greatly contribute to the optimal design of PHP for electronic heat dissipation in small space. In the present work, the heat transfer performance of PHP is optimized via structure analysis and modeling calculation. A new “spring-mass-damper” model in terms of separated phase flow mode is established, where the frictional pressure loss of the real two-phase flow pattern − slug flow in PHP is considered. Besides, a prototype of PHP with adjustive-structured channel (ASCPHP) is proposed. The heat transfer performance of ASCPHP is evaluated with the newly established model. With theoretical computation method, the frequency of ASCPHP the superiority of ASPHP is also confirmed by comparison with other types of PHPs.

Linearity of the Co-moving Velocity

AbstractThe co-moving velocity is a new variable in the description of immiscible two-phase flow in porous media. It is the saturation-weighted average over the derivatives of the seepage velocities of the two immiscible fluids with respect to saturation. Based on analysis of relative permeability data and computational modeling, it has been proposed that the co-moving velocity is linear when plotted against the derivative of the average seepage velocity with respect to the saturation, the flow derivative. I show here that it is enough to demand that the co-moving velocity is characterized by an additive parameter in addition to the flow derivative to be linear. This has profound consequences for relative permeability theory as it leads to a differential equation relating the two relative permeabilities describing the flow. I present this equation together with two solutions.

Investigation of Improved Energy Dissipation in Stepped Spillways Applying Bubble Image Velocimetry

This study investigates skimming flow regimes, two-phase air–water flow conditions, and simple measures to improve energy dissipation in stepped spillways. Experiments were conducted using two different scale physical models, 1:50 and 1:17, within separate rectangular flumes to define scale effects. Flow patterns were analyzed using the Bubble Image Velocimetry (BIV) technique, which tracks air bubbles. The introduction of splitters resulted in a 7% increase in relative energy dissipation. Additionally, the length of inception was reduced to Li/ks = 10, thereby decreasing the potential for subsequent cavitation. Beyond the BIV experiments, two experiments were conducted on the large-scale model using Acoustic Doppler Velocimetry (ADV), with and without splitters, to examine the impact of splitters on the velocity profile above the crest. In the experiment with splitters, the vertical velocity vector (v) contributed to turbulence by changing direction, thereby reducing average velocities both in front of and behind the ogee crest. This led to a reduction in energy on the downstream side of the spillway. Although the small-scale model appears unsuitable for studying two-phase flow, the change in relative energy dissipation from the baseline to the splitter configuration was practically identical for both scale models, thereby supporting the findings of the large-scale model.

Effect of fiber curvature on gas diffusion layer two-phase dynamics of proton exchange membrane fuel cells

Both straight and curved carbon fibers are widely used in various commercial gas diffusion layer (GDL) fabrications. The effect of the different carbon fiber curvatures on two-phase flow dynamics within the cathode GDLs of proton exchange membrane fuel cells remains unclear. In this study, we investigate liquid transport in three types of GDLs with varying fiber curvatures using the two-phase volume of fluid simulations in OpenFOAM. For the first time, a rod periodic surface model is combined with a layer-by-layer fiber stacking strategy, to stochastically reconstruct GDL structures while incorporating crucial parameters from physical (experimental) GDLs. A grid independence study and model validation are conducted. Following pore network analysis of pore size distribution and connectivity, we study the time-varying GDL total and local water saturation and capillary pressure. Despite maintaining similar layer and bulk porosity, increased fiber curvature enhances pore connectivity but raises water saturation and capillary pressure, increasing the risk of flooding. Additionally, droplets in gas channels with straight-fiber GDLs are larger and have slower movement than those in curved-fiber GDLs. Fiber curvature inversely affects drainage capacity in GDLs and connected channels. With comparable water saturation and capillary pressure, curved-fiber GDLs exhibit lower discrepancies, suggesting improved uniformity in water distribution.