An experimental investigation of gas-liquid two-phase flow in single microchannel contactors

This paper describes two-phase flow pattern and pressure drop characteristics during the absorption of CO₂ into water in three horizontal microchannel contactors which consist of Y-type rectangular microchannels having hydraulic diameters of 667, 400 and , respectively. With the help of a high-speed photography system, flow patterns such as bubbly flow, slug flow (including two sub-regimes, Taylor flow and unstable slug flow), slug-annular flow, churn flow and annular flow were observed in these microchannels. The applicability of the currently available correlations for describing flow pattern transitions in microchannels has been examined. Generally, the predicting performance of these correlations deteriorates as the channel diameter further reduces. Toward solving this discrepancy, an empirical correlation based on the superficial Weber numbers was developed to interpret the transition from Taylor flow to unstable slug flow in three microchannels. Taylor bubble formation process in microchannels was found to be in the squeezing regime at lower superficial liquid velocities (Ca ranging from 0.0019 to 0.029) while the transition to the dripping regime was observed at the highest superficial liquid velocity of 1.0 m/s. Lengths of Taylor bubbles formed in the squeezing regime can be well represented by the scaling relation proposed by Garstecki et al. [Formation of droplets and bubbles in a microfluidic T-junction—scaling and mechanism of break-up. Lab on a Chip, 6, 437-446]. For flow patterns including slug-annular flow, annular flow and churn flow, a simple analysis based on the separated flow model has been performed in order to reveal the observed effect of the superficial liquid velocity on two-phase frictional multiplier in the present microchannels. Then, reasonable correlations for the prediction of two-phase frictional pressure drop under these flow patterns were suggested.     

Hydrodynamics and mass transfer characteristics in gas-liquid flow through a rectangular microchannel

Researches on two-phase transfer and reaction processes in microchannnels are important to the design of multiphase microchemical systems. In the present work, hydrodynamics and mass transfer characteristics in cocurrent gas-liquid flow through a horizontal rectangular microchannel with a hydraulic diameter of have been investigated experimentally. Liquid side volumetric mass transfer coefficients were measured by absorbing pure CO₂ into water and a 0.3 M NaHCO₃ / 0.3 M Na₂CO₃ buffer solution. Interfacial areas were determined by absorbing pure CO₂ into a 1 M NaOH solution. Two-phase flow patterns and pressure drop data were also obtained and analyzed. This paper shows that two-phase frictional pressure drop in the microchannel can be well predicted by the Lockhart-Martinelli method if we use a new correlation of C value in the Chisholm's equation. Liquid side volumetric mass transfer coefficient and interfacial area as high as about and , respectively, can be achieved in the microchannel. Generally, liquid side volumetric mass transfer coefficient increases with the increasing superficial liquid or gas velocity, which can be described satisfactorily by the developed empirical correlations. A comparison of mass transfer performance among different gas-liquid contactors reveals that the gas-liquid microchannel contactor of this study can provide at least one or two orders of magnitude higher liquid side volumetric mass transfer coefficients and interfacial areas than the others.     

Effect of gas expansion on the velocity of individual Taylor bubbles rising in vertical columns with water: Experimental studies at atmospheric pressure and under vacuum

This study was designed to determine the effect of gas expansion on the velocity of Taylor bubbles rising individually in a vertical column of water. This experimental study was conducted at atmospheric pressure or under vacuum (33.3 and ) using three different acrylic columns with internal diameters of 0.022, 0.032, and 0.052 m, and more than 4.0 m high. A non-intrusive optical method was used to measure velocity and length of Taylor bubbles at five different locations along the columns. The operating conditions used correspond to inertial controlled regime.In experiments performed under vacuum, there is considerable gas expansion during the rise of Taylor bubbles, particularly when they approach the liquid free surface where the pressure drop (due to the hydrostatic pressure) is of the order of magnitude of the absolute pressure. The liquid ahead of the bubble is displaced upward by an amount proportional to the gas expansion resulting in increased bubble velocity. The calculated Reynolds number suggests a laminar regime in the liquid ahead of the bubble. However, the experimentally determined velocity coefficient C for each column was much smaller than 2, which would be expected for laminar flow. The value of C obtained ranges from 1.13±0.09, for the narrowest column, to 1.40±0.24, for the widest column. This suggests that a fully developed laminar flow in the liquid ahead of the bubble is never achieved due to continuous bubble expansion at a variable rate, regardless of column height.The velocity coefficient C can be used to calculate the contribution of liquid motion to bubble velocity. Subtracting this contribution from the measured bubble velocity defines a constant value which is nearly identical to the bubble rise velocity measured in the same column operated as a constant volume system (two ends closed) where gas expansion is absent.     

Induced pulsing in trickle beds—characteristics and attenuation of pulses

Gas/liquid down-flow in packed beds is studied, under periodic liquid feeding (at sufficiently high frequencies to be classified as “fast” mode of pulsing), in a range of mean liquid and gas flow rates within the steady “trickling flow regime”. The aim is to identify periodic feeding conditions resulting in improved fluid-mechanical characteristics (e.g. uniform fluids distribution) and possibly enhanced transport rates in this flow regime, which is common in industrial processes. From instantaneous, cross-sectionally averaged holdup measurements, at various locations along the packed bed, quantitative information is obtained on the axial propagation and attenuation of induced pulses. A phenomenological treatment of the pulse decay process facilitates data interpretation and leads to the determination of a characteristic attenuation factor for the various conditions tested. Key parameters of the process studied include, in addition to dynamic holdup, pressure drop, pulse celerity and intensity, as a function of fluid feed rates (G,L) and liquid cyclic frequency. Under the conditions of these tests, and for fixed mean rates G,L, the time averaged holdup and the pulse celerity are practically constant along the bed; furthermore, these quantities as well as the pressure drop do not seem to be affected by the imposed cyclic liquid feeding frequency. An expression to tentatively correlate pulse celerity data is recommended.The computed attenuation factors indicate that there is a rather narrow band of mean gas and liquid rates (along the so-called “pseudo-transition” boundary to pulsing flow) where pulse decay is at a minimum. Based on these results as well as on pulse intensity vs. bed length data, recommendations are made on preferred conditions for induced pulsing (from the fluid-mechanical standpoint) which would maximize expected benefits.     

An experimental study of air-water Taylor flow and mass transfer inside square microchannels

Flow and mass transfer properties under air-water Taylor flow have been investigated in two square microchannels with hydraulic diameters of 400 and 200 μm. Experimental data on Taylor bubble velocity, pressure drop and liquid side volumetric mass transfer coefficient (k_La) have been presented. It was shown that the measured Taylor bubble velocity in square microchannels could be well interpreted based upon an approximate measurement of the liquid film profile therein. Then, the obtained two-phase frictional pressure drop values in both microchannels were found to be significantly higher than the predictions of the correlation proposed by Kreutzer et al. [2005b. Inertial and interfacial effects on pressure drop of Taylor flow in capillaries. A.I.Ch.E. Journal 51, 2428-2440] when the liquid slug was very short, which can be explained by the inadequacy of their correlation to describe the excess pressure drop caused by the strong inner circulation in such short liquid slugs. An appropriate modification has been made to this correlation in order to improve its applicability in microchannels. Finally, the experimental (k_La) values in the microchannel with hydraulic diameter of 400 μm were found to be in poor agreement with those predicted by the existing correlations proposed for capillaries with diameters of several millimeters. The observed deviation was mainly due to the fact that mass transfer experiments in this microchannel actually corresponded to the case of short film contact time and rather poor mixing between the liquid film and the liquid slug, which was not in accordance with mass transfer assumptions associated with these correlations. A new empirical correlation has been proposed to describe mass transfer data in this microchannel.     

Prediction of pulsation frequency of pulsing flow in trickle-bed reactors using artificial neural network

Based on an extensive experimental database (946 measurements) set up from the literature published over past 30 years, a new correlation relying on artificial neural network (ANN) was proposed to predict the basic pulsation frequency of pulsing flow in the trickle-bed reactors. Seven dimensionless groups employed in the proposed correlation were liquid and gas Reynolds (Re_L,Re_G), liquid Weber (We_L), gas Froude (Fr_G), gas Stokes (St_G) and liquid Eötvös numbers and a bed correction factor (S_b). The performance comparisons of literature and present correlations showed that ANN correlation is significantly an improvement in predicting pulsation frequency with an AARE of 10% and a standard deviation less than 18%. The effects of the variables including the properties of fluid and bed, and flow rate of liquid and gas on pulsing frequency were investigated by ANN parametric simulations and the trends were compared with exiting experimental results that confirmed the coherence of the proposed method with the previous experiments.     

Heat transfer characterization of gas-liquid flows in a trickle-bed

Instantaneous local fluid-solid heat transfer coefficient (h_t) in a laboratory scale trickle-bed was measured using a constant-voltage anemometry technique. It was observed that convective heat transfer rate in the liquid-rich pulses was approximately 4 times that in gas-continuous bases for the air-water system. Time-averaged heat transfer rate was found to be positively influenced by both gas and liquid flow rates, with a stronger dependence on the latter. Heat removal efficiency, taking pressure drop penalty into account, suggested an optimum at intermediate liquid flow rate. Based on the measurements, a four-parameter heat transfer model featuring heat transfer coefficients in liquid-rich pulses (h_tp) and gas-continuous bases (h_tb), pulsing frequency and pulse fraction was developed to characterize transient h_t under various flow regimes. This model can be used in any trickle-bed reactor simulation that accounts for the dynamic interactions of catalytic reactions and heat transfer. It was found that while h_tp and h_tb correspond to liquid-solid and gas-solid heat transfer, respectively, and are determined mainly by the fluid properties, pulsing frequency and pulse fraction are the factors characterizing different flow regimes. Pulsing frequency, which can significantly impact reaction, may be tuned by selecting appropriate packing size, since smaller sizes generate higher frequency pulses. For example, a two-fold higher frequency was detected in packing as compared to that with packing. Flow regime evolution along the column axial location was identified visually, while the dispersed bubbling flow retreating to pulsing flow owing to gas bubble coalescence was evidenced by the heat transfer measurements.     

Bubble size effect on low liquid input drift-flux parameters

We investigated the effect of bubble size on the drift-flux parameters at low liquid flow conditions by measuring the radial profiles of void fraction and phase velocities in a vertical bubbly pipe flow of diameter and height . To study the effect of the bubble size we used two different types of bubble inlets. We measured the local bubble fraction and velocity U_g by using single and four-point-optical fibre probes, and we used Laser Doppler Anemometry to determine the liquid velocity U_l. The distribution parameter C₀ and the weighted mean drift velocity |U_drift| were directly computed from the local measurements at a height on our experimental set-up. Both parameters were influenced by the bubble size. Provided no liquid flow reversal occurred at the near wall region, the distribution parameter reached a below unity minimum plateau value of C₀=0.95 for wall peaking void fraction profiles. At low liquid input conditions both the liquid input and bubble size had an influence on the distribution parameter. Extreme values such as C₀>2 were measured. From these measurements we developed models for the drift-flux parameters to take into account the effect of bubble size and input-flow conditions for our intermediate pipe diameter value. These models were tested and validated with separately collected experimental data.     

Hydrodynamics of three-phase fluidized bed systems examined by statistical, fractal, chaos and wavelet analysis methods

Hydrodynamic studies were conducted in gas-liquid-solid systems (0.1 m ID, 2 m high) of 3.0 mm glass beads and of 2.1 mm polypropylene low-density particles, with particles densities of 2471 and 1290 kg/m³, respectively. Simultaneous measurement of differential pressure and bubble conductivity probe signals sampled at 500 Hz for 60 s enabled the investigation of the change in flow structure in relation to the flow regime transitions. Superficial gas velocities ranged between 0.010 and 0.052 m/s for polypropylene particles, and extended to 0.12 m/s for glass beads, while the superficial liquid velocities covered the ranges of 0.0007-0.045 m/s for polypropylene particles, and ranged up to 0.056 m/s for glass beads.Spectral analysis of the pressure fluctuations revealed a transition from dispersed to coalesced bubbling flow with decreasing liquid velocity for a given superficial gas velocity. The use of a conductivity probe facilitated characterization of the local flow structure in terms of bubble movement. The measurements were extensively analyzed using fractals and chaos, power spectra frequency analysis and wavelet decomposition in addition to the standard statistical analyses. The coefficient of variation of the bubble probe signals was found to be the most effective in deducing the transition velocity between coalesced and dispersed bubbling flow regimes, while wavelet energy confirmed the similarity in the distribution between two axial positions once operated in the dispersed flow regime. Comparison of the flow structure between glass beads and polypropylene particles showed that both the minimum liquid fluidization velocity and the transition velocity between the bubble flow regimes were much higher for the glass beads than for the lighter polypropylene particles. Furthermore, the standard deviations of the decomposed bubble probe signals through wavelet transformation successfully highlighted the difference between the two systems of particles.     

An experimental study of flow and heat transfer of supercritical carbon dioxide in multi-port mini channels under cooling conditions

This paper presents the fluid flow and heat transfer characteristics of supercritical CO₂ in a horizontal multi-port extruded aluminum test section consisting of 10 circular channels with an inner diameter of 1.31 mm. Both local and average pressure drop and heat transfer coefficients were measured as CO₂ was cooled in the multi-port circular channels with pressures ranging from 7.4 to 8.5 MPa, inlet fluid temperatures ranging from 22 to , and mass velocity ranging from 113.7 to 418.6 kg/m² s. The results indicate that the operating pressure, the mass velocity and the temperature of CO₂ had significant effects on fluid flow and heat transfer characteristics. The pressure drop and the average heat transfer coefficient increased greatly with increasing the average temperatures of CO₂ in the near-critical region; the average heat transfer coefficient attained a peak value near the corresponding pseudocritical temperature; and the maximum heat transfer coefficient decreased as the pressure increased. Both the pressure drop and the heat transfer coefficient increased with the mass velocity, but decreased with the operating pressure. The measured average heat transfer coefficients were compared with the experimental data reported in the literatures and a large discrepancy was observed. Based on the experimental data collected in the present work, a new correlation was developed for forced convection of supercritical CO₂ in horizontal multi-port mini channels under cooling conditions.     

EFFECT OF VISCOSITY ON HYDRODYNAMIC PROPERTIES OF PULSING FLOW IN TRICKLE BEDS

New data on pulsing flow onset, properties of pulses (frequency, celerity, length), liquid holdup and pressure drop are presented for aqueous glycerol solutions of viscosity 6.7 and 20.2 mPa s and compared with similar measurements from an air-water (1.0 mPas) system. With the exception of viscosity, all other physical properties of the liquid phase are kept constant and fairly close to those of water, thus allowing a direct assessment of the effect of viscosity. Pulse formation and propagation with viscous liquids is examined on the basis of time records from a conductance type technique. A striking effect due to increased liquid viscosity is the reduction of the pulsing flow regime; in particular, the pulsing-to-bubbling transition boundary is shifted towards higher gas flow rates. Pulse frequency and celerity appear to decrease only slightly with increasing liquid viscosity, whereas the two-phase pressure gradient increases significantly. Liquid holdup also tends to increase with viscosity. Moreover, holdup with viscous liquids tends to increase significantly with the liquid flow rate, whereas an insignificant effect is found for water. A new correlation for estimating liquid holdup is proposed, and a simple model for predicting pulsing flow characteristics is modified in order to take account of the aforementioned effects.     

Static mixers with a gas continuous phase

The aim of this work was to characterise hydrodynamics and mass transfer in a gas-liquid contactor containing static mixers (SMs). The originality of this study lies in the fact that these mixing organs are used with a gas continuous phase. Two types of SM were implemented in co-current flows, Statiflo and Lightnin. The pressure drop ΔP, the volumic interfacial area a and the volumic mass transfer coefficient k_La were measured in several configurations: horizontal flow, vertical up-flow and vertical down-flow. The influences of position and flow rates were studied in order to understand the behaviour of these contactors, and to optimise the operating conditions. As expected, the pressure drop was found to increase mainly with gas velocity but also with liquid velocity, and to reach 3300 Pa in the range of velocities studied (the gas flow rate varied between 4 and and the liquid flow rate between 0 and 100 L/h), far less than Sülzer SM. The volumic interfacial area and the volumic mass transfer coefficient showed the same changes, a varying between 100 and , and k_La reaching 0.07 L/s. This is interesting compared with other classical absorption processes: indeed, even if packing towers can provide the same range of values, the operating conditions are more drastic or the dimensions of the apparatuses are far larger than SM ones. The position was also found to have an influence on the hydrodynamic and mass transfer parameters (ΔP, a and k_La).     

Induced pulsing in trickle beds — Particle shape and size effects on pulse characteristics

Periodic liquid feeding of the ON-OFF type is investigated — at sufficiently high frequencies to be classified as “fast” mode of induced pulsing — in the range of mean gas and liquid flow rates corresponding to the steady “trickling flow” regime. Two of the most common types of catalyst-support particles, i.e. porous spherical and cylindrical extrudates, are employed to study the imposed pulse characteristics. Detailed information is obtained, on the axial propagation and attenuation of pulses, from instantaneous, cross-sectionally averaged holdup measurements. Key fluid-mechanical parameters studied include, aside from dynamic holdup and pressure drop, pulse celerity and intensity, as a function of fluid feed rates (G,L) and liquid cyclic frequency. Similar published data, for 6 mm glass spheres, are employed for comparison; it is shown that, for the particles examined, particle size has a pronounced effect, but not as significant as that of particle shape. For particles of comparable size, the cylindrical shape is associated with much greater global dynamic holdup and pressure drop, and with increased pulse attenuation. Moreover, packed extrudates exhibit a significant increase of holdup in the axial direction, recently also observed in steady trickling flow. For spherical particles, both time-average holdup and pulse celerity are practically constant along the bed for fixed L,G. Pressure drop, global holdup and pulse celerity are not affected by cyclic liquid feeding frequency, for both spherical and cylindrical extrudate particles. Based on the pulse attenuation characteristics, for the three particle types examined, recommendations are made on preferred conditions for induced pulsing (from the fluid dynamics point of view) which would maximize benefits. Overall, it appears that spherical packings hold significant advantages over cylindrical extrudates of comparable size. Finally, in view of the observed significant decay of imposed pulses along the bed, care should be exercised to properly interpret data obtained in short laboratory reactors (where pulse attenuation is limited) for scale-up of the much longer industrial beds.     

Local liquid saturation measurements inside a trickle bed reactor operating near the transition between pulsing and trickling flow

A wire-mesh tomography device was used to study the liquid saturation at 78 points covering the cross-sectional area of a 0.30 m diameter trickle bed reactor. Measurements in the pulsing flow near the transition were done for glass beads and alumina cylinders using air and water as fluids. Local liquid saturation measurements allow identification of flow regime and pulse frequency for each of the 78 points considered. To the best of the authors’ knowledge this is the first time that such local measurements are done for the complete sectional area of a trickle bed reactor. Flow rate conditions at which different flow regimes coexist at the same horizontal plane of the column were identified. A remarkable influence of the initial liquid distribution on the observed flow regime was also observed.     

Liquid mixing in gas-liquid two-phase flow by meandering microchannels

The influence of the channel radius on the mass transfer in rectangular meandering microchannels (width and height of ) has been investigated for gas-liquid flow. Laser induced velocimetry measurements have been compared with theoretical results. The symmetrical velocity profile, known from the straight channel, was found to change to an asymmetrical one for the meandering channel configuration. The changes in the secondary velocity profile lead to an enhanced radial mass transfer inside the liquid slug, resulting in a reduced mixing length. In the investigated experimental range (superficial gas velocity and superficial liquid velocity ) the mixing time was reduced eightfold solely due to changes in channel geometry. An experimental study on the liquid slug lengths, the pressure drop and their relation to the mass transfer have also been performed. Experimental results were validated by a simulation done in Comsol Multiphysics^®. To obtain information for higher velocity rates, simulations were performed up to . These velocity variations in the simulation indicate the occurrence of a different flow pattern for high velocities, leading to further mass transfer intensification.     

Comparison of normal phase operation and phase reversal studies in a pulsed sieve plate extraction column

The hydrodynamic characteristics of a pulsed sieve plate extraction column (PSPEC) was studied experimentally using two different liquid phase systems, namely water/kerosene and 30%TBP (tributyl phosphate) in NPH (normal paraffin hydrocarbon)/0.3 M HNO₃. The aqueous phase as the dispersed phase and the organic phase as the continuous phase (phase reversal) and vice versa (normal phase operation) studies in a pulsed sieve plate extraction column 0.076 m in diameter and 1 m height are presented in this paper. The hydrodynamic properties like drop size and holdup are characterized as a function of various operating parameters namely pulse velocity, dispersed phase and continuous phase velocity and duty cycle of pulsing. Flooding in the column was also investigated for the changes involving flow ratio of continuous phase to that of the dispersed phase for both insufficient and excessive pulsing. It has been observed that phase reversal mode of operation is not efficient as compared to normal phase operation for the PSPEC.     

Multiple hydrodynamic states in trickle flow: Quantifying the extent of pressure drop, liquid holdup and gas-liquid mass transfer variation

It is well established that pressure drop and liquid holdup under trickle flow conditions are functions of the flow history. However, the extent of possible variation of these and other critical hydrodynamic parameters has not been fully quantified. In this study, specifically defined prewetting procedures are used as limiting cases for hydrodynamic hysteresis. These are:

•

Non-prewetted.

•

Levec prewetted: the bed is flooded and drained and after residual holdup stabilisation the gas and liquid flows are introduced.

•

Kan_L prewetted: the bed is operated in the pulse flow regime (by increasing liquid velocity) after which liquid flow rate is reduced to the desired set point (all at the desired gas flow rate).

•

Kan_G prewetted: the bed is operated in the pulse flow regime (by increasing gas velocity) after which gas flow rate is reduced to the desired set point (all at the desired liquid flow rate).

•

Super prewetted: the bed is flooded and gas and liquid flows are introduced once draining commences.

It is shown that the upper limiting case for pressure drop is the Kan_L mode of operation. The lower limiting cases are the non-prewetted and Levec prewetted modes (these coincide). Pressure drop may vary by as much as 700% even for prewetted beds. Liquid holdup is different in all five prewetting modes. The upper limiting case is the Kan_G mode of operation, while the lower limiting case is the non-prewetted mode (Kan_G holdup is approximately 160% that of non-prewetted mode holdup at ). At low gas velocities the Kan_L holdup can be 400% of that of the non-prewetted beds. Importantly, the lower limiting case for prewetted beds is the Levec mode. Holdup in the Kan_G mode may be as much as 130% of the holdup in the Levec mode (at ).The effect of hydrodynamic multiplicity of the volumetric mass transfer coefficient is measured by the desorption of oxygen from water into nitrogen. In this case the different prewetting procedures result in three distinct regions, the upper region being the Kan and Super prewetted beds, the intermediate region being the Levec prewetted bed and the lower region being the dry bed. Mass transfer coefficients in the upper region can be as much as 600% of that of the lower region and 250% of that of the intermediate region. Evidently, prewetting (and even pulsing flow prewetting) does not guarantee that the bed is operating at the maximum values of pressure drop, holdup and mass transfer coefficient. Evidence of operation in between the limiting cases is presented. These non-limiting cases can be reached in multiple ways.     

Effects of liquid viscosity and surface tension on atomization in two-phase, gas/liquid fluid coker nozzles

This is an experimental study of the effects of liquid viscosity (μ_L) and surface tension (γ) on the size of droplets produced by an atomizer fed with a two-phase gas/liquid (TPGL) mixture. The nozzle is a one-quarter scale model of an industrial fluid coker nozzle used for bitumen upgrading. The primary objective of this study was to compare the Sauter mean diameter of the TPGL sprays for different μ_L and γ. Compressed air was used as the gas phase; the liquid phases were water, canola oil and a glycerine-water solution. The μ_L was varied from 1 mPa s to 67 mPa s, and γ was varied from 25 mN/m to 61 mN/m. The liquid flow rates were varied from 0.095 L/s to 0.105 L/s, and the gas-to-liquid ratio , by mass, was fixed at 2%, similar to the ratio used in commercial fluid coker nozzles. Fluid mixing pressures in the test were between 400 kPa and 700 kPa. The D₃₂ within the spray was measured using a Dantec 2D-Fibre mode Phase-Doppler-Particle-Anemometer (PDPA). With measurements performed at axial distances of 100 and 202 mm from the nozzle orifice and within a spray radius of 40 mm. The results showed that the 67-fold increase in μ_L results in 46 μm increase in D₃₂ which roughly extrapolates to a power law relationship exponent of 0.11 (based on two points). In contrast, the 2.4-fold increase in γ indicates roughly 42 μm increase in D₃₂, giving an exponent of −0.24. For both parameters the exponent is small, indicating weak dependence. The γ effects appear stronger, but the modest increase is barely exceeds the margin of error.     

Hydrodynamics and mass transfer in a novel multi-airlifting membrane bioreactor

A novel multiple-airlifting membrane bioreactor is built with four sintered stainless steel tubular filters as the risers and downcomers. This work investigates the hydrodynamics including gas holdup, liquid velocity, liquid circulation and mixing times by aerating different number of risers (one to three) at superficial gas velocities of 0.02-0.07 m/s The mass transfer phenomena, including oxygen mass transfer (k_La) and effective molecular diffusivity of lactic and acetic acids through the walls of tubular filters, are also investigated. It is found that gas holdup in individual risers increases linearly with the superficial gas velocity, and performs independently under multiple-airlifting conditions. The vessel-based gas holdup and liquid velocity in downcomer(s) increase with aeration rate of individual risers as well as the number of risers. The liquid velocity in downcomers reaches an upper limit (about 0.6 m/s), because of flow resistance or energy loss of liquid circulation. The oxygen mass transfer coefficient (k_La) is primarily affected by gas holdup and the number of risers, and to some extent influenced by liquid velocity. The novel airlifter configuration results in good liquid mixing in the bioreactor that quickly reaches new steady state in response to a sudden pH change from acid addition.