Investigation on the void fraction of gas–liquid two-phase flows in vertically-downward pipes

A special experimental loop is designed and constructed to study the characteristics of the void fraction of gas–liquid two-phase flow in vertically-downward pipes. The test section is made of transparent pipe with a length of 6 m and an internal diameter of 25 mm. The void fraction ranging from 0.1 to 0.98 widely is measured using quick-closing valve method. It is found that the range of the void fraction could be divided into three regions with different flow patterns and different relationships between the void fraction and the gas–liquid volumetric flow rate ratio. Moreover, 39 correlations for calculating the void fraction collected from present literature, are classified, and evaluated using the experimental data obtained in this study. The prediction of correlations in the literature needs to be improved when the void fraction is small.     

Heat transfer model for gas–liquid slug flows under constant flux

This paper investigates the mechanisms leading to enhanced heat and/or mass transfer rates in two-phase non-boiling slug flows. The problem is analyzed in a minichannel geometry subjected to a constant heat flux boundary. Local Nusselt numbers, obtained using Infrared thermography are analyzed in both entrance and fully developed flow regions. These novel measurements highlight the physics governing slug-flow heat transfer and results indicate that optimized slug geometries can yield up to an order of magnitude heat transfer enhancement. Finally, based on the physics identified, a heat transfer model is developed which is also applicable to similar mass transfer problems.     

Physical modelling and advanced simulations of gas–liquid two-phase jet flows in atomization and sprays

This review attempts to summarize the physical models and advanced methods used in computational studies of gas–liquid two-phase jet flows encountered in atomization and spray processes. In traditional computational fluid dynamics (CFD) based on Reynolds-averaged Navier–Stokes (RANS) approach, physical modelling of atomization and sprays is an essential part of the two-phase flow computation. In more advanced CFD such as direct numerical simulation (DNS) and large-eddy simulation (LES), physical modelling of atomization and sprays is still inevitable. For multiphase flows, there is no model-free DNS since the interactions between different phases need to be modelled. DNS of multiphase flows based on the one-fluid formalism coupled with interface tracking algorithms seems to be a promising way forward, due to the advantageous lower costs compared with a multi-fluid approach. In LES of gas–liquid two-phase jet flows, subgrid-scale (SGS) models for complex multiphase flows are very immature. There is a lack of well-established SGS models to account for the interactions between the different phases. In this paper, physical modelling of atomization and sprays in the context of CFD is reviewed with modelling assumptions and limitations discussed. In addition, numerical methods used in advanced CFD of atomization and sprays are discussed, including high-order numerical schemes. Other relevant issues of modelling and simulation of atomization and sprays such as nozzle internal flow, dense spray, and multiscale modelling are also briefly reviewed.     

Thermocapillary instabilities in liquid columns under co – and counter-current gas flows

The thermocapillary flows in an infinite liquid column surrounded by an annular channel of gas with the axial temperature gradient are investigated. The gas is pumped through the channel parallel to the interface with a given flow rate. The solution for stationary motion is derived and possible flow regimes are analyzed depending on the Marangoni number and dimensionless gas flow rate. Linear stability analysis of these regimes is performed. It is shown that gas pumping in the same direction with respect to the thermocapillary motion on the interface has a destabilizing effect on the system. Gas pumping in the opposite direction can be stabilizing or destabilizing depending on the gas flow rate. It is shown that the motion of liquid can be completely vanished by the gas flow with a specified flow rate. The obtained results demonstrate the possibility of controlling thermocapillary flows and their stability in liquid columns (liquid bridges) by gas flows.     

Estimation of wettability in gas–liquid–rock systems

The wettability in gas–liquid–rock systems, including geothermal (steam–water–rock) systems, is central to defining the relative permeability and capillary pressure, which govern the forecasts of future production and reinjection performance. Since the frequently used Amott and USBM approaches are appropriate only for specific liquid–liquid–rock systems, we developed a model to evaluate the wettability in gas–liquid–rock systems. The method can be applied to determine the wettability at specific wetting-phase saturation if the capillary pressure and relative permeability at this fluid saturation are known. The proposed model can also be used in liquid–liquid–rock systems. The validity of the technique was tested qualitatively and quantitatively using experimental data from different rock–fluid systems. The calculated wettability indices in gas–liquid–rock systems were higher than in liquid–liquid–rock systems, as expected. Results also show that the wettability index calculated using the new method may or may not be independent of fluid saturation in the cases studied. In addition, the wettability index in drainage was seen to be greater than that in imbibition, a wettability difference that is particularly important in steam–water–rock (geothermal) systems because of the impact on reinjection.     

Three-dimensional measurement of gas dissolution process in gas–liquid microchannel flow

The three-dimensional CO₂ dissolution process through a gas–liquid interface in microfluidic devices was investigated experimentally, for the precise control of CO₂ dissolution. The gas dissolution was evaluated by using confocal micron-resolution particle image velocimetry (micro-PIV) combined with laser induced fluorescence (LIF), which has the ability to measure the velocity and dissolved CO₂ concentration distribution in a liquid flow field. The measurement system is based on the confocal microscope, which has excellent depth resolution and enables visualization of the three-dimensional distributions of velocity and dissolved CO₂ concentration by rendering two-dimensional data. The device is comprised of a polydimethylsiloxane chip, whose microchannels were fabricated by using a cryogenic micromachining system. The width and depth of the liquid flow channel are larger than those of the gas flow channel. This is due to the need for decreasing the width of the gas–liquid interface and increasing the hydraulic diameter of the liquid channel, whose conditions generate a static gas–liquid interface. The experiments were performed for three different liquid flow conditions corresponding to Reynolds numbers of 1.0 × 10^?2, 1.2 × 10^?2 and 1.7 × 10^?2, and the gas flow rate was set to be constant at 150 μL/min. The LIF measurements indicate that an increase in the Reynolds number yields a decrease in dissolved gas in the spanwise directions. Furthermore, molar fluxes by convection and diffusion were evaluated from the experimental data. The molar fluxes in the streamwise direction were at least 20 times as large as those in the spanwise and depthwise directions. This reveals that an increase in momentum transport in the spanwise and depthwise directions is an important factor for enhancing mass transfer in the gas–liquid microchannel flow.     

Analysis of mixed convection—Radiation interaction in a vertical channel: Entropy generation

This work examines analytically the effects of radiation heat transfer on mixed convection through a vertical channel in the presence of transverse magnetic field. Both First and Second Laws of thermodynamics are applied to analyze the problem. Special focus is given on the entropy generation characteristics and its dependency on the various dimensionless parameters, i.e. Hartmann number (Ha), Plank number (Pl), Richardson number (Ri), group parameter (Br/Π) etc. A steady-laminar flow of an incompressible-viscous fluid is assumed through the channel with negligible inertia effect. Fluid is further considered as an optically thin gas and electrically conducting. Governing equations in Cartesian coordinate are solved analytically after reasonable simplifications. Expressions for velocity, temperature, local, and average entropy generation rate are derived and presented graphically.     

Thermal stochastic collision model in turbulent gas–solid pipe flows

New thermal stochastic particle collision model in gas–solid flow in a riser is developed. The simulation is based on four-way coupling of phases considering inter-particle collision and heat transfer. It is shown that the limitation of excessive computational time in Eulerian–Lagrangian simulation of gas–solid flows for the high loading ratios is eliminated by using the stochastic particle collision model. The simulation results demonstrate that the predictions of the developed thermal stochastic particle collision modem are in good agreement with those obtained by the direct particle collision model and the available experimental data. The new stochastic modeling is used and nearly dense gas–solid flow is simulated for high loading ratios up to eight and the results are presented and discussed.     

Prandtl and capillary effects on heat transfer performance within laminar liquid–gas slug flows

This paper investigates a two-phase non-boiling slug flow regime for the purposes of enhancing heat transfer in microchannel heat sinks or compact heat exchangers. The primary focus is upon understanding the mechanisms leading to enhanced heat transfer and the effect of using different Prandtl number fluids, leading to variations in Capillary number. Experimental work was conducted using Infrared thermography and results are presented in the form of Graetz solution, spanning both the thermal entrance and fully developed flow regions. Nusselt numbers enhancements were observed throughout when data was reduced to account for void fraction. However, the gaseous void was also noted to demonstrate an artificial increase with greater thicknesses of the liquid film, due to higher Capillary numbers. Up to 600% enhancement in heat transfer rates were observed over conventional Poiseuille flow. This was verified through Nusselt number measurements over inverse Graetz number ranges from 10^?4 to 1 and slug length to channel diameter ratios from 0.88 to 32. Varying Prandtl and Capillary numbers caused notable effects in the transition region between entrance and fully developed flows. Significant Nu oscillations were observed for low Pr fluids due to internal circulation within the slug. However, these oscillations are observed to be damped out when higher Prandtl number fluids are employed. The thickness of the liquid film surrounding the gas bubbles is shown to have a significant influence on heat transfer performance. Overall, this study provides a greater understanding of the mechanisms leading to significant enhancements in heat exchange devices employing two-phase gas–liquid flows without boiling.     

Heat transfer of gas–liquid mixture in micro-channel heat sink

The main objective of the present investigation is to study heat transfer in parallel micro-channels of 0.1 mm in size. Comparison of the results of this study to the ones obtained for two-phase flow in “conventional” size channels provides information on the complex phenomena associated with heat transfer in micro-channel heat sinks. Two-phase flow in parallel micro-channels, feeding from a common manifold shows that different flow patterns occur simultaneously in the different micro-channels: liquid alone (or single-phase flow), bubbly flow, slug flow, and annular flow (gas core with a thin liquid film, and a gas core with a thick liquid film). Although the gas core may occupy almost the entire cross-section of the triangular channel, making the side walls partially dry, the liquid phase always remained continuous due to the liquid, which is drawn into the triangular corners by surface tension. With increasing superficial gas velocity, a gas core with a thin liquid film is observed. The visual observation showed that as the air velocity increased, the liquid droplets entrained in the gas core disappeared such that the flow became annular. The probability of appearance of different flow patterns should be taken into account for developing flow pattern maps. The dependence of the Nusselt number, on liquid and gas Reynolds numbers, based on liquid and gas superficial velocity, respectively, was determined in the range of Re_LS = 4–56 and Re_GS = 4.7–270. It was shown that an increase in the superficial liquid velocity involves an increase in heat transfer (Nu_L). This effect is reduced with increasing superficial gas velocity, in contrast to the results reported on two-phase heat transfer in “conventional size” channels.     

Flow and heat transfer in a micro-cylindrical gas–liquid Couette flow

Analytic solutions for the gas and liquid velocity and temperature distribution are determined for steady state one-dimensional microchannel cylindrical Couette flow between a shaft and a concentric cylinder. The solution is based on the continuum model and takes into consideration the velocity slip and temperature jump in the gaseous phase defined by the Knudsen number range of 0.001 < Kn < 0.1. The two fluids are assumed immiscible. The gas layer is adjacent to the shaft which rotates with angular velocity ω_s and is thermally insulated. The outer cylinder rotates with angular velocity ω_o and is maintained at uniform temperature. The governing parameters are identified and the effects of the Knudsen number and accommodation coefficients on the velocity and temperature profiles, reduction in the overall temperature rise due to the gas layer, the Nusselt number and shear reduction are examined. It was found that the required torque to rotate the liquid in the annular space is significantly reduced by introducing a thin gas layer adjacent to the shaft. Also, reduction in shaft temperature is enhanced through a combination of high energy accommodation coefficient and low momentum accommodation coefficients. Results also indicate that the gas layer becomes more effective in reducing the shaft temperature when the housing angular velocity is much larger than the shaft angular velocity.     

VOF modelling of gas–liquid flow in PEM water electrolysis cell micro-channels

In this study, the gas–liquid flow through an interdigitated anode flow field of a PEM water electrolysis cell (PEMEC) is analysed using a three-dimensional, transient, computational fluid dynamics (CFD) model. To account for two-phase flow, the volume of fluid (VOF) method in ANSYS Fluent 17.2 is used. The modelled geometry consists of the anode channels and the anode transport layer (ATL). To reduce the complexity of the phenomena governing PEMEC operation, the dependence upon electro-chemistry is disregarded. Instead, a fixed source of the gas is applied at the interface between the ATL and the catalyst layer. An important phenomenon that the model is able to capture is the gas–liquid contact angle on both the channel wall and ATL-channel interface. Particularly, the latter interface is crucial in capturing bubble entrainment into the channel. To validate the numerical simulation, photos taken of the gas–liquid flow in a transparent micro-channel, are qualitative compared against the simulation results. The experimental observations confirm the models prediction of long Taylor bubbles with small bubbles in between. From the simulation results, further intriguing details of the flow are revealed. From the bottom to the top of the outgoing channel, the film thickness gradually increases from zero to 200 μm. This increase in the film thickness is due to the particular superficial velocity field that develops in an interdigitated flow. Here both the superficial velocities change along the length of the channel. The model is capable of revealing effect of different bubble shapes/lengths in the outgoing channel. Shape and the sequence of the bubbles affect the water flow distribution in the ATL. The model presented in this work is the first step in the development of a comprehensive CFD model that comprises multiphase flow in porous media and micro-channel, electro-chemistry in catalyst layers, ion transport in membrane, hydrogen evolution, etc. The model can aid in the study of gas–liquid flow and its impact on the performance of a PEMEC.     

Direct numerical simulation of interphase heat and mass transfer in multicomponent vapour–liquid flows

A Volume-of-Fluid methodology for direct numerical simulation of interface dynamics and simultaneous interphase heat and mass transfer in systems with multiple chemical species is presented. This approach is broadly applicable to many industrially important applications, where coupled interphase heat and mass transfer occurs, including distillation. Volume-of-Fluid interface tracking allows investigation of systems with arbitrarily complex interface dynamics. Further, the present method incorporates the full interface species and energy jump conditions for vapour–liquid interphase heat and mass transfer, thus, making it applicable to systems with multiple phase changing species. The model was validated using the ethanol–water system for the cases of wetted-wall vapour–liquid contacting and vapour flow over a smooth, stationary liquid. Good agreement was observed between empirical correlations, experimental data and numerical predictions for vapour and liquid phase mass transfer coefficients. Direct numerical simulation of interphase heat and mass transfer offers the clear advantage of providing detailed information about local heat and mass transfer rates. This local information can be used to develop accurate heat and mass transfer models that may be integrated into large scale process simulation tools and used for equipment design and optimization.     

Experimental investigation of reaction kinetics of gas–liquid–solid Bunsen reaction in iodine–sulfur process

Iodine–sulfur (IS) cycle is the most promising thermochemical water-splitting process for nuclear hydrogen production. The Bunsen reaction, which produces sulfuric and hydriodic acid for the two decomposition reactions, plays a crucial role for the continuous stable operation of the IS cycle. Insufficient kinetics studies on Bunsen reaction, particularly under the gas–liquid–solid heterogeneous conditions, have caused difficulties for the design of Bunsen reactor, as well as the optimization and improvement of the efficiency of the process. In this work, the reaction kinetics of gas–liquid–solid Bunsen reaction denoted in the pressure drop of SO₂ was experimentally investigated, and the effluences of the main factors, including the initial SO₂ pressure, molar ratio of I₂ to H₂O, temperature, and stirring rate, were studied. In addition, a kinetics model for simulating the heterogeneous reaction was proposed and verified by the experimental data obtained under the three-phase Bunsen reaction conditions.     

Determinants of global natural gas consumption and import–export flows

Based on the data of the BP Statistical Review of World Energy, this paper constructs the consumption and import–export of natural gas identities. It discusses the drivers of changes in global natural gas consumption and trade flows from 2008 to 2015 using the extended logarithmic mean Divisia index. The results show that differences in the natural gas supply and demand across countries or regions, as well as the distribution of energy between the domestic and international markets, can be better explained when natural gas trade movements are considered. By comparing the supply and consumption increment of natural gas, this study finds that only the energy intensity, economic growth, and demographic effects are consistent with each other. The changes in the impact of other effects mainly depend on storage variations and statistical errors. In addition, the primary drivers of the incremental changes in natural gas consumption vary in different countries. They include production scale, import scale, export scale, consumption structure proportion, energy intensity, economic growth, and population and balance effects. Finally, the consumption competitiveness of the liquefied natural gas significantly improved over the examined period.     

Experimental investigation of opposing mixed convection hear transfer in the vertical flat channel in a laminar–turbulent transition region

In this paper we present the results on experimental investigation of the local opposing mixed convection heat transfer in the vertical flat channel with symmetrical heating in a laminar–turbulent transition region. The experiments were performed in airflow (p = 0.1–1.0 MPa) in the range of Re from 1.5 × 10³ to 6.6 × 10⁴ and Gr_q up to 1 × 10¹¹ at the limiting condition q_w1 = q_w2 = const. The analysis of the results revealed significant increase in the heat transfer with increasing of air pressure (Gr number). Also sharp increase in heat transfer was noticed in the region with vortex flow in comparison with the turbulent flow region.     

Modeling of Navier–Stokes equations for high Knudsen number gas flows

The possibility of modeling the Navier–Stokes equations and together with the conventional second order slip boundary condition at high Knudsen numbers is explored in this paper by incorporating the Knudsen diffusion phenomenon in rarefied gases. An effective mean free path (MFP) model is augmented to the governing equation and the slip boundary condition, as gas transport properties can be related to the MFP. This simple modification is shown to implicitly take care of the complexities associated in the transitional flow regime, without necessitating dependency of the slip coefficients on the Knudsen number. Unique analytical model with fixed values of slip coefficients is proposed and rigorous comparisons with the experimental and simulation data for pressure driven and thermally driven rarefied gas flows support this conjecture. First and second order slip coefficients have been proposed as 1.1466 and 0.9756 for rectangular channels and 1.1466 and 0.14 for the capillaries, from the continuum to the transition flow regime. The current work is significant from the numerical simulation point of view because simulation tools are better developed for Navier–Stokes equations.     

Experimental investigation of water droplet–air flow interaction in a non-reacting PEM fuel cell channel

It has been well documented that water production in PEM fuel cells occurs in discrete locations, resulting in the formation and growth of discrete droplets on the gas diffusion layer (GDL) surface within the gas flow channels (GFCs). This research uses a simulated fuel cell GFC with three transparent walls in conjunction with a high speed fluorescence photometry system to capture videos of dynamically deforming droplets. Such videos clearly show that the droplets undergo oscillatory deformation patterns. Although many authors have previously investigated the air flow induced droplet detachment, none of them have studied these oscillatory modes. The novelty of this work is to process and analyze the recorded videos to gather information on the droplets induced oscillation. Plots are formulated to indicate the dominant horizontal and vertical deformation frequency components over the range of sizes of droplets from formation to detachment. The system is also used to characterize droplet detachment size at a variety of channel air velocities. A simplified model to explain the droplet oscillation mechanism is provided as well.     

The processes of dissolution and hydrate forming behind the shock wave in the gas–liquid medium with gas mixture bubbles

The processes of dissolution and hydrate forming behind the front of the shock wave, propagating in liquid with bubbles of carbon dioxide and nitrogen mixture, have been experimentally studied under various initial pressures, medium temperatures and surfactant concentrations in water. The theoretical model of this process is presented for the gas–liquid medium with bubbles of multicomponent gas mixture with consideration of accompanying heat effects. Numerical solutions to this problem were found. Close fit of calculation results and experimental data was achieved. Experimental data on profiles of gas content, times of dissolution and hydrate forming were generalized on the basis of the suggested model.