Experimentally‐based constitutive relations for co‐current gas‐liquid flow in randomly packed beds

Experimental Observations on average pulse velocity and frequency in concurrent gas‐liquid (down) flow through randomly packed beds are used to extract constitutive relations for the gas‐liquid interaction and mean curvature terms that appear in a recently proposed volume‐averaged two‐fluid model for bubbly flow. The proposed closures lead to a reasonably quantitative prediction of the average pressure drop and liquid saturation under bubbly flow conditions and in the near pulse regime. In addition, the proposed closures provide realistic estimates for the location of the bubble‐to‐pulse transition in microgravity and in 1g down‐flow and predict the disappearance of the bubbly flow pattern at low liquid fluxes in 1g down‐flow. © 2016 American Institute of Chemical Engineers AIChE J, 63: 812–822, 2017     

A new mechanistic model to predict gas–liquid interface shape of gas–liquid flow through pipes with low liquid loading

Devising a new mechanistic method to predict gas–liquid interface shape in horizontal pipes is concerned in this article. An experiment was conducted to find the pressure gradients of air–water flow through a 1‐in. pipe diameter. Comparing results of model with some experimental data available in the literature demonstrates that the model provides quite better predictions than existed models do. This model also predicts flow regime transition from stratified to annular flow better than Apparent Rough Surface and Modified Apparent Rough Surface models for both 1‐ and 2‐in. pipe diameters. The model also leads to reliable predictions of wetted wall fraction experimental data. Although one parameter of new model was evaluated based on air–water flow pressure loss experimental data for 1 in. pipe, it was considerably successful to predict pressure drop, liquid holdup, stratified‐annular transition and wetted wall fraction for other gas–liquid systems and pipe diameters. © 2014 American Institute of Chemical Engineers AIChE J, 61: 1043–1053, 2015     

Equilibrium and non‐equilibrium gas–liquid two phase flow in long and short pipelines following a rupture

The two‐phase flow following the blowdown of pipeline carrying flashing liquid is numerically investigated by using thermodynamic equilibrium and non‐equilibrium models. Model equations are solved numerically by the finite volume method. The values of fluxes at cell boundaries are obtained by AUSM+‐up. To obtain proper values for the coefficients of dissipation, both single phase liquid and two phase shock tube problems are investigated. The transient release from the pressurized pipeline is studied for two cases of long and short pipes. Comparison of the predictions against experimental data reveals non‐equilibrium model performs a little better than equilibrium model in the prediction of temporal variations of pressure and void fraction of the long pipe. However, equilibrium model totally overestimates pressure and void fraction of the short pipe. The relative error of equilibrium model in the prediction of pressure variation with time exceeds 50% and it is 20% for non‐equilibrium model. © 2017 American Institute of Chemical Engineers AIChE J, 63: 3214–3223, 2017     

Liquid phase residence time distribution for gas–liquid flow in Kenics static mixer

The residence time distribution (RTD) of the liquid phase for co-current gas–liquid upflow in a Kenics static mixer (KSM) with air/water and air/non-Newtonian fluid systems was investigated. The effect of liquid and gas superficial velocities on liquid holdup and Peclet number was studied. Experiments were conducted in three KSMs of diameter 2.54 cm with 16 elements and 5.08 cm diameter with 8 and 16 elements, respectively, of constant L_e/D_e = 1.5 for different liquid and gas velocities. A correlation was developed for Peclet number, in terms of generalized liquid Reynolds number, gas Froude number and liquid Galileo number, where as for liquid holdup, a correlation was developed as a function of gas Reynolds number. The axial dispersion model was found to be in good agreement with the experimental data.     

Sensor‐based flow pattern detection—gas–liquid–liquid upflow through a vertical pipe

Flow distribution during gas–liquid–liquid upflow through a vertical pipe is investigated. The optical probe technique has been adopted for an objective identification of flow patterns. The probability density function (PDF) analysis of the probe signals has been used to identify the range of existence of the different patterns. Dispersed and slug flow have been identified from the nature of the PDF, which is bimodal for slug flow and unimodal for dispersed flow. The water continuous, oil continuous, and emulsion type flow distributions are distinguished on the basis of the PDF moments. The method is particularly useful at high flow rates where visualization techniques fail. Based on this, a flow pattern detection algorithm has been presented. Two different representations of flow pattern maps have been suggested for gas–liquid–liquid three phase flow. © 2014 American Institute of Chemical Engineers AIChE J, 60: 3362–3375, 2014     

Liquid‐bridge flow in the channel of helical string and its application to gas–liquid contacting process

To solve the problems of the traditional packings, such as high pressure drop, mal‐distribution and short liquid residence time, a helical flow structured packings was proposed. Two different flow patterns, liquid‐bridge flow and liquid‐drop flow were identified when the width of the channel of the helical string was adjusted. Moreover, the characteristics of the helical liquid‐bridge flow including maximum liquid loading, mean thickness of liquid film, mean residence time and effective specific surface area, were examined. And the separation efficiency was studied by the lab‐scale distillation column. In comparison, the effective specific surface area of the helical flow type packings is almost as large as the traditional B1‐350Y structured packings, but with thinner liquid film, longer liquid residence time and finally higher separation efficiency. © 2018 American Institute of Chemical Engineers AIChE J, 64: 3360–3368, 2018     

Flow pattern transition in gas‐liquid downflow through narrow vertical tubes

The present report studies on the flow pattern transitions during vertical air water downflow through millichannels (0.83 ≤ Eötvös no. ≤ 20.63). Four basic flow patterns namely falling film flow, slug flow, bubbly flow, and annular flow are observed in the range of experimental conditions studied and their range of existence has been noted to vary with tube diameter and phase velocities. Based on experimental observations, phenomenological models are proposed to predict the transition boundaries between adjacent patterns. These have been validated with experimental flow pattern maps from the present experiments. Thus the study formalizes procedure for developing a generalized flow pattern map for gas‐liquid downflow in narrow tubes. © 2016 American Institute of Chemical Engineers AIChE J, 63: 792–800, 2017     

Viscous co‐current downward Taylor flow in a square mini‐channel

This article presents a computational study of the co‐current downward Taylor flow of gas bubbles in a viscous liquid within a square channel of 1 mm hydraulic diameter. The three‐dimensional numerical simulations are performed with an in‐house computer code, which is based on the volume‐of‐fluid method with interface reconstruction. The computed (always axi‐symmetric) bubble shapes are validated by experimental flow visualizations for varying capillary number. The evaluation of the numerical results for a series of simulations reveals the dependence of the bubble diameter and the interfacial area per unit volume on the capillary number. Correlations between bubble velocity and total superficial velocity are also provided. The present results are useful to estimate the values of the bubble diameter, the liquid film thickness and the interfacial area per unit volume from given values of the gas and liquid superficial velocities. © 2009 American Institute of Chemical Engineers AIChE J, 2010     

Chemical conversion and liquid–liquid extraction of 5‐hydroxymethylfurfural from fructose by slug flow microreactor

A dehydration of fructose in the water/methyl isobuthyl ketone (MIBK) biphasic system can yield 5‐hydroxymethylfurfural (HMF) to be successfully extracted into the organic MIBK phase. The HMF production and yield in MIBK phase was discussed by using a simplified model taking into consideration of the slug flow. The extraction resistance of HMF across the interface between water and MIBK depended on the line velocity and the flow rate ratio. It was likely that the velocity field generated in the slug flow contributed to an increase in the mass transfer of HMF. © 2016 American Institute of Chemical Engineers AIChE J, 62: 2135–2143, 2016     

Phase separation of gas–liquid two‐phase stratified and plug flows in multitube T‐junction separators

Using air and water as the working fluids, phase separation phenomena for stratified and plug flows at inlet were investigated experimentally, at a simple T‐junction and specifically designed multitube T‐junction separators with two or three layers. The results show that for these two flow patterns the separation efficiency of the two phases for any multitube T‐junction separator is much higher than that of the simple T‐junction. Increasing the number of connecting tubes in the multitube T‐junction separator can increase the separation efficiency. Generally, for stratified flow, complete separation of the two phases can be achieved by the two‐layer multitube T‐junction separator with five or more connecting tubes and by the three‐layer separator; increasing the gas flow rate, the liquid flow rate, or the mixture velocity under plug flow is detrimental to phase separation with a drop in peak separation efficiency. © 2016 American Institute of Chemical Engineers AIChE J, 63: 2285–2292, 2017     

Numerical simulation of dilute turbulent gas‐particle flow with turbulence modulation

A numerical study of a dilute turbulent gas‐particle flow with inelastic collisions and turbulence modulation in an Eulerian framework is described. A new interpretation is provided for the interaction/coupling terms, based on a fluctuating energy transfer mechanism. This interpretation provides for a new robust closure model for the interaction terms with the ability to predict the turbulence dampening as well as the turbulence enhancement phenomenon. Further, the model developed herein is investigated along with a variety of other published closure models used for the interaction/coupling terms, particle drag, and solid stress. The models are evaluated against several sets of benchmark experiments for fully‐developed, turbulent gas‐solid flow in a vertical pipe. © 2011 American Institute of Chemical Engineers AIChE J, 2012     

Bubble splitting under gas–liquid–liquid three‐phase flow in a double T‐junction microchannel

Gas–aqueous liquid–oil three‐phase flow was generated in a microchannel with a double T‐junction. Under the squeezing of the dispersed aqueous phase at the second T‐junction (T2), the splitting of bubbles generated from the first T‐junction (T1) was investigated. During the bubble splitting process, the upstream gas–oil two‐phase flow and the aqueous phase flow at T2 fluctuate in opposite phases, resulting in either independent or synchronous relationship between the instantaneous downstream and upstream bubble velocities depending on the operating conditions. Compared with two‐phase flow, the modified capillary number and the ratio of the upstream velocity to the aqueous phase velocity were introduced to predict the bubble breakup time. The critical bubble breakup length and size laws of daughter bubbles/slugs were thereby proposed. These results provide an important guideline for designing microchannel structures for a precise manipulation of gas–liquid–liquid three‐phase flow which finds potential applications among others in chemical synthesis. © 2017 American Institute of Chemical Engineers AIChE J, 63: 376–388, 2018     

Counter‐current gas‐liquid wavy film flow between the vertical plates analyzed using the Navier‐Stokes equations

The article is devoted to a theoretical analysis of counter‐current gas‐liquid wavy film flow between vertical plates. We consider two‐dimensional nonlinear waves on the interface over a wide variation of parameters. The main interest is to analyse the wave structure at the parameter values corresponding to the onset of flooding observed in experiments. We use the Navier‐Stokes equations in their full statement to describe the liquid phase hydrodynamics. For the gas phase equations, we use two models: (1) the Navier‐Stokes system and (2) the simplified Benjamin‐Miles approach where the liquid phase is a small disturbance for the laminar or turbulent gas flow. With the superficial gas velocity increasing and starting from some value of the velocity, the waves demonstrate a rapid decreasing of both the minimal film thickness and the phase wave velocity. We obtain a region of the gas velocity where we have two solutions at one set of the problem parameters and where the flooding takes place. Both the phase wave velocity and the minimal film thickness are positive numbers at such values of the velocity. We calculate the flooding point dependences on the liquid Reynolds number for two different liquids. The wave regime corresponding to the flooding point demonstrates negative u‐velocities in the neighbourhood of the interface near the film thickness maximum. At smaller values of the superficial gas velocity, the negative u‐velocities take place in the neighbourhood of the film thickness minimum. © 2009 American Institute of Chemical Engineers AIChE J, 2010     

An analysis of liquid CO2 drop formation with and without hydrate formation in static mixers

The formation process of CO₂ drops in various types of Kenics Static Mixers was analyzed from the perspective of energy dissipation in the mixer, focusing on the formation of drop surfaces. Experimental studies on CO₂ drop formation were conducted under varying temperatures, pressure, and flow rates, with and without hydrate formation. Analysis of the CO₂ drop size and distribution at several locations within the static mixer was conducted, as of pressure drop in the mixer, to determine dissipation energies. In all the experimental conditions, by considering the surface energy for hydrate formation, the energy required for the formation of CO₂ drops correlated well with total energy dissipation by mixer flow, which is represented by a pressure drop along the mixer. This process has important applications to the formation of liquid CO₂ for ocean disposal as a countermeasure to global warming. © 2010 American Institute of Chemical Engineers AIChE J, 2010     

Hydrodynamics and mass transfer of gas–liquid flow in a falling film microreactor

In this article, flow pattern of liquid film and flooding phenomena of a falling film microreactor (FFMR) were investigated using high‐speed CCD camera. Three flow regimes were identified as “corner rivulet flow,” “falling film flow with dry patches,” and “complete falling film flow” when liquid flow rate increased gradually. Besides liquid film flow in microchannels, a flooding presented as the flow of liquid along the side wall of gas chamber in FFMR was found at high liquid flow rate. Moreover, the flooding could be initiated at lower flow rate with the reduction of the depth of the gas chamber. CO₂ absorption was then investigated under the complete falling flow regime in FFMR, where the effects of liquid viscosity and surface tension on mass transfer were demonstrated. The experimental results indicate that k_L is in the range of 5.83 to 13.4 × 10^?5 m s^?1 and an empirical correlation was proposed to predict k_L in FFMR. © 2009 American Institute of Chemical Engineers AIChE J, 2009     

Spectral characterization of static mixers. The S‐shaped micromixer as a case study

We investigate the steady‐state performance of a planar micromixer composed of several S‐shaped units. Mixing efficiency is quantified by the decay of the scalar variance downstream the device for generic feeding conditions. We discuss how this decay is controlled by the spectral properties of the advection‐diffusion Floquet operator, $ \cal F $ , that maps a generic scalar profile at the inlet of a single unit into the corresponding profile at the unit outlet section. Two advantages characterize the Floquet operator approach —(i) it allows to analyze an arbitrarily long device and (ii) it provides a quantitative assessment of mixing efficiency that is independent of the feeding conditions and that depends solely on the interaction between advection and diffusion. © 2009 American Institute of Chemical Engineers AIChE J, 2010     

Multistage rotor‐stator spinning disc reactor

The scale up of a rotor‐stator spinning disc reactor by stacking single stage rotor‐stator units in series is demonstrated. The gas‐liquid mass transfer per stage is equal to the mass transfer in a single stage spinning disc reactor. The pressure drop per stage increases with increasing rotational disc speed and liquid flow rate. The pressure drop is more than a factor 2 higher for gas‐liquid flow than for liquid flow only, and is up to 0.64 bar at 459 rad s^?1. The high mass and heat transfer coefficients in the (multistage) rotor‐stator spinning disc reactor make it especially suitable for reactions with dangerous reactants, highly exothermic reactions and reactions where selectivity issues can be solved by high mass transfer rates. Additionally, the multistage rotor‐stator spinning disc reactor mimics plug flow behavior, which is beneficial for most processes. © 2011 American Institute of Chemical Engineers AIChE J, 2012     

Numerical study of turbulent liquid‐liquid dispersions

A numerical approach is developed to gain fundamental insight in liquid‐liquid dispersion formation under well‐controlled turbulent conditions. The approach is based on a free energy lattice Boltzmann equation method, and relies on detailed resolution of the interaction of the dispersed and continuous phase at the microscopic level, including drop breakup and coalescence. The capability of the numerical technique to perform direct numerical simulations of turbulently agitated liquid‐liquid dispersions is assessed. Three‐dimensional simulations are carried out in fully periodic cubic domains with grids of size . The liquids are of equal density. Viscosity ratios (dispersed phase over continuous phase) are in the range 0.3–1.0. The dispersed phase volume fraction varies from 0.001 to 0.2. The process of dispersion formation is followed and visualized. The size of each drop in the dispersion is measured in‐line with no disturbance of the flow. However, the numerical method is plagued by numerical dissolution of drops that are smaller than 10 times the lattice spacing. It is shown that to mitigate this effect it is necessary to increase the resolution of the Kolmogorov scales, such as to have a minimum drop size in the range 20–30 lattice units [lu]. Four levels of Kolmogorov length scale resolution have been considered , 2.5, 5, and 10 [lu]. In addition, the numerical dissolution reduces if the concentration of the dispersed phase is increased. © 2015 American Institute of Chemical Engineers AIChE J, 61: 2618–2633, 2015     

Analysis of roll wave characteristics under low liquid loading two‐phase flow conditions

An experimental study is conducted using a 0.152‐m ID facility to investigate the wave characteristics of two‐phase stratified wavy flow in horizontal pipelines. The experiments are conducted under low liquid loading condition, which is very commonly observed in wet gas pipelines. The experiments are conducted with water as the liquid phase, and repeated with 51 wt % of monoethylene glycol (MEG) in the aqueous phase to analyze the effects of MEG presence on wave characteristics. The experimental range of this study covers superficial gas velocity, v_Sg, values of 9–23 m/s and superficial liquid velocity, v_SL, values of 0.01–0.02 m/s. Similar test matrices are completed for the cases with and without MEG in the aqueous phase. A conductivity probe system is used to measure the wave characteristics at the liquid–gas interface. These characteristics include the wave celerity, frequency, amplitude, length, and liquid film thickness. The experimental oil–air wave characteristics data of Gawas et al. (Int J Multiphase Flow. 2014;63:93–104) is also used for comparison purposes. The trends in the resulting wave characteristics with respect to input parameters are investigated, for oil, water, or MEG–water mixture as the liquid phase. Common predictive methods for interfacial wave celerity, including shallow water theory, Watson (Proceedings of the 4th International Conference in Multi‐Phase Flows, Nice, France. 1989:495–512), Paras et al. (Int J Multiphase Flow. 1994;20(5):939–956), Al‐Sarkhi et al. (AIChE J. 2012;58(4):1018–1029), and Gawas et al. (Int J Multiphase Flow. 2014;63:93–104) are evaluated in comparison with the experimental data. The results of the wave frequency correlation of Al‐Sarkhi et al. (AIChE J. 2012;58(4):1018–1029) are also compared with the experimental wave frequency data. Lastly, a correlation is developed to predict the relative wave amplitude, as a function of superficial gas Weber number and liquid velocity number. Most of the commonly used two‐phase stratified flow models are developed with the assumption of steady‐state conditions, and neglect the transient wave effects. This study provides valuable experimental results on wave characteristics of stratified wavy flow for different types of liquid phase. Moreover, a comprehensive analysis of the parameters affecting the wave characteristics of stratified wavy flow is presented. © 2017 American Institute of Chemical Engineers AIChE J, 63: 3177–3186, 2017