Buoyancy-driven bubbly flows: scaling of velocities in bubble columns operated in the heterogeneous regime

Abstract The hydrodynamics of bubble columns in the heterogeneous regime is revisited. Focusing on air–water systems at large aspect ratio, we show from dimensional analysis that buoyancy equilibrates inertia, and that velocities scale as $(gD\varepsilon )^{1/2}$, where $D$ is the bubble column diameter, $\varepsilon$ the void fraction and $g$ the gravitational acceleration. From new experiments in a $0.4$ m diameter column with ${{O}}(10^3)$ particle Reynolds number bubbles and from a detailed analysis of published data, we confirm the self-organization prevailing in the heterogeneous regime, and that the liquid flow rate is only set by the column diameter $D$. Besides, direct liquid and gas velocity measurements demonstrate that the relative velocity increases above the terminal velocity $U_T$ in the heterogeneous regime, and that it tends to ${\sim }2.4 U_T$ at very large gas superficial velocities $V_{sg}$. The proposed velocity scaling is shown to hold for liquid and gas mean velocities and for their standard deviations. Furthermore, it is found to be valid over a wide range of conditions, corresponding to Froude numbers $Fr=V_{sg}/(gD)^{1/2}$ from 0.02 to 0.5. Then, the relevance of this scaling for coalescing media is discussed. Moreover, following the successful prediction of the void fraction with a Zuber & Findlay approach at the beginning of the heterogeneous regime, we show how the void fraction is correlated with $Fr$. Further investigations are finally suggested to connect the increase in relative velocity with meso-scale structures known to exist in the heterogeneous regime.


Introduction
In bubble columns gas is injected at the bottom of an initially stagnant liquid contained in a vertical cylinder.Such systems are often used in industry as reactors (chemical and biochemical transformations), in separation techniques (flotation), to promote agitation and mixing (metallurgy). . .At low inlet gas flow rate, bubbles homogeneously rise over the column cross-section.In consequence, mixing and turbulence are mainly generated by interactions of bubble wakes, and the void fraction is linearly increasing with the gas superficial velocity V sg (V sg is defined as the ratio of the injected gas flow rate to the entire column cross-section).When the inlet gas flow rate is increased above some threshold, the flow loses its spatial uniformity and its unsteadiness growths.In this so-called heterogeneous regime, complex flow structures appear as secondary motions are superimposed on the mean recirculation arising at the reactor scale.These motions are said to be 'chaotic' by Noël De Nevers who noticed in 1968: 'In unbaffled systems these circulations are unstable and chaotically change in size, shape, and orientation.These chaotic circulations provide the principal mode of vertical bubble transport in bubble columns over a wide range of operating conditions'.In addition, the increase of the void fraction with V sg turns to become nonlinear (Joshi et al. 1998;Ruzicka et al. 2001;Sharaf et al. 2016).
In industrial applications, superficial velocities are typically in the range 5 to 30 cm s −1 and the gas volume fraction evolves between 10 % and 40 %: most of these systems are thus operated in the heterogeneous regime that corresponds to values of gas hold up above 20 %.A key issue in this regime is the extrapolation of results acquired in small-scale experiments to full-size reactors whose diameters D can be as large as many metres.Indeed, in spite of a sustained scientific production over more than 70 years, the hydrodynamics of bubble columns operated in the heterogeneous regime is not yet fully understood.In particular, there is still no consensus on the scaling of key variables such as void fraction, mean liquid velocity, velocity fluctuations. . .As an illustration of that situation, more than twenty different correlations are currently proposed to evaluate the void fraction (see the reviews by Deckwer & Field 1992;Joshi et al. 1998;Kantarci, Borak & Ulgen 2005;Rollbusch et al. 2015;Kikukawa 2017;Besagni, Inzoli & Ziegenhein 2018).
Similarly, and despite progresses in two-fluid modelling, simulations of bubble columns based on two-fluid approaches have not yet reached a fully predictive status (Shu et al. 2019).Indeed, ad hoc adjustments are still required to reach some agreement with experiments.Concerning the momentum transferred from bubbles to the liquid, in earlier attempts it was the bubble size that was adjusted to modulate the relative velocity (Ekambara, Dhotre & Joshi 2005).However, the current approach consists in multiplying the drag force by an ad hoc coefficient function of the local void fraction in order to represent an effective momentum exchange between phases.That correction combines an increase of the drag at low void fractions -the well-known hindering effect (e.g.Ishii & Zuber (1979) that was originally identified in solid suspensions by Richardson & Zaki 1954), with a neat decrease at larger void fractions -with the so-called swarm effect (Ishii & Zuber 1979;Simonnet et al. 2007).This swarm effect represents the impact of neighbour bubbles on the motion of a test inclusion: it notably arises from the entrainment of bubbles in the wake of larger inclusions.In air/water systems, drag reduction has been observed to start at void fractions ∼15 %, and to reach a factor 5 at a void fraction of approximately 30 %. Various empirical expressions have been proposed for the swarm coefficient: let us quote Roghair et al. (2011);McClure et al. (2017); Gemello et al. (2018) and Yang et al. (2018). . .Despite slight variations (on the value of the critical void fraction, on some extra dependency on the mean bubble size), all these proposals are quite similar and they all mainly depend on the local void fraction.When used in bubble column simulations, such swarm correction leads to a strong increase of the apparent relative velocity in the heterogeneous regime, a trend that is in qualitative agreement with experiments (Krishna, Wilkinson & Van Dierendonck 1991;Raimundo et al. 2019).Provided that the bubble size is known beforehand (implying that coalescence/breakup mechanisms -whose modelling are also important issues -are not active in the flow conditions considered), such corrections ensure a reasonable agreement with air/water experiments with deviations up to 15 % in void fractions (global and local) and up to 40 % in the liquid velocity on the axis (Ertekin et al. 2021).Interesting, these figures hold over a significant range of conditions, namely for gas superficial velocities up to 25 cm s −1 and for column diameters from 0.2 to 3 m, indicating that the swarm effect is a key feature of the heterogeneous regime.
Other uncertainties in the modelling of bubble column hydrodynamics arise in the description of the turbulence in the liquid phase for which different approaches have been attempted (Khan, Bhusare & Joshi 2017) and a new production model has been proposed (Panicker, Passalacqua & Fox 2020).This issue also concerns the description of the agitation induced by bubbles as a variety of mechanisms are able to generate velocity fluctuations in the continuous phase (Risso 2018) including cluster induced turbulence (Capecelatro, Desjardins & Fox 2015;Shu et al. 2019;Panicker et al. 2020;Baker et al. 2020) due to meso-scale structures.Furthermore, the presence of such structures in bubble columns operated in the heterogeneous regime has been demonstrated in recent experiments (Raimundo et al. 2019).
Hence, despite the limitations of current turbulence and agitation models, the introduction of an ad hoc swarm coefficient in simulations happens to provide a somewhat reasonable agreement with experiments in terms of void fraction.Its seems thus that some robust physics takes place in the complex flows prevailing in the heterogeneous regime.To identify that physics, we revisit in this paper the hydrodynamics of bubble columns.Our starting point is that recent experimental results support the idea that a strong analogy exists between a bubble column in the heterogeneous regime and turbulent buoyancy-driven flows in confined channels with zero mean flow.This prompted us to hypothesize that a dynamical equilibrium between inertia and buoyancy holds in the heterogeneous regime: such an equilibrium leads to a liquid velocity that scales as V liquid ∼ (gDε) 2 , where g is the gravitational acceleration, D the column diameter and ε stands for the void fraction (Cartellier 2019).In this paper, that scaling proposal is analysed and tested against previous experimental data reported in the literature.It is also tested against a new experimental dataset collected in a D = 0.4 m column that includes gas phase velocity statistics measured with a newly developed Doppler optical probe (Lefebvre et al. 2019(Lefebvre et al. , 2022)).We will notably show that, for bubble columns operated in the heterogeneous regime, both mean and fluctuating axial velocities of the liquid phase and of the gas phase closely follow the proposed scaling.Moreover, and to complement the velocities scalings, an empirical expression will be proposed for the void fraction that is backed up by a dimensional analysis, and by the Zuber & Findlay (1965) one-dimensional model relating void fraction and gas flow rate fraction.The comparison of that model with experiments shows that the axial evolution of the flow is significant, and that it deserves to be investigated to capture the flow structuration in the heterogeneous regime.
3 to 35 cm s −1 ) while keeping a fixed bubble size.More precisely, the Sauter mean horizontal diameter was the same within the range ±1 mm in all columns at a given value of V sg : it increased from 6 to 8 mm over the range of V sg studied with a mean eccentricity always close to 0.7.Accordingly, the equivalent diameter evolved from 4.7 to 7 mm, and thus, the terminal velocity was nearly constant, equal to 0.21 ± 0.01 cm s −1 whatever the flow conditions.Coalescence being avoided (or at least, it was too weak to influence the flow behaviour), we have been able to clarify some key features of bubble columns operated in the heterogeneous regime. Notably: (i) The homogeneous-heterogeneous transition was observed with a fixed bubble size, meaning that, contrary to the current belief, coalescence is not necessary for such a transition to occur.(ii) Local void fraction fluctuations happen to be very significant in the heterogeneous regime: they evolve between one tenth and ten times the average gas hold-up when quantified by a one-dimensional Voronoï analysis (Raimundo 2015;Mezui, Cartellier & Obligado 2018;Raimundo et al. 2019).
The co-existence of 'dense' regions corresponding to clusters of bubbles with regions almost 'free of bubbles' called 'voids' (Raimundo et al. 2019) induces strong differences in local velocities: the bubble transport is controlled by these clusters/voids meso-scale structures.This could be the origin of the observed increase in the apparent relative velocity of the gas in the heterogeneous regime, an effect usually accounted by way of a swarm coefficient.Moreover, the presence of clusters and voids in the mixture induces strong local shear rates as well as intense three-dimensional (3-D) vortical structures that are expected to significantly contribute to turbulence production.Clearly, the fluctuations in the mixture density induce strong spatial and temporal fluctuations in buoyancy (see figure 1 and movies in the supplementary movies): they are thus reminiscent of convective instabilities arising in turbulent buoyancy-driven flows with zero mean flow.
Some analogies also hold on the global flow structure.Let us first discuss the existence of a quasi-fully developed region in bubble columns operated in the heterogeneous regime.Such a region appears in systems that are not too strongly confined: according to Wilkinson, Spek & van Dierendonck (1992), a minimum internal diameter of 0.15 m is a necessary condition (possibly, that condition is helpful to avoid the development of a slug flow regime).In addition, the aspect ratio should be large enough so that end effects do not affect the flow organization in the central portion of the column: Wilkinson et al. (1992) argue that the dynamic height H D of the mixture should exceed 5 diameters while Forret (2003) established that H D /D = 3 is a sufficient condition for large (namely D = 0.4 m and 1 m) columns.When introducing the static liquid height H 0 , these conditions transform into H 0 /D 3.8 or 2.3, respectively, indicating that the bubble column should not be operated in the shallow water limit for a quasi-fully developed region to exist.Moreover, when the above conditions are fulfilled, the way gas injection is performed has no impact on the flow organization outside the entrance region.That conclusion has been ascertained in air-water systems with a gas injection evenly distributed over the column cross-section, and for large enough injection orifices (orifice diameter above 1 mm according to Wilkinson et al. (1992), above 0.5 mm according to Sharaf et al. 2016).
When one considers the so-called 'pure' heterogeneous regime, i.e. flow conditions such that the void fraction vs the superficial velocity is concave and that are far enough from the homogeneous/heterogeneous transition (Ruzicka 2013;Sharaf et al. 2016), that conclusion seems to hold also when changing the coalescence efficiency by way of surfactants or of Images of the flow in the vicinity of column walls between approximately 0.8 and 2 m above gas injection, in the homogeneous regime (a-d) and in the heterogeneous regime (e-h).As side lightning was used, the grey level is an indication of the presence of bubbles: liquid structures comprising a few bubbles appear as dark zones while clusters of bubbles correspond to bright zones.Results correspond to an air-water bubble column with D = 0.4 m and static liquid height H 0 = 2 m.The unsteadiness of these structures can be appreciated from the movies included in the supplementary movies available at https://doi.org/10.1017/jfm.2022.833.The time increment between images is 1/30 s.
water purity (Sharaf et al. 2016).The precise extent of the quasi-fully developed region is not exactly known: it is said to range from the end of the entrance region whose extent is approximately one (Forret et al. 2006) or two (see Guan et al. (2016) and references therein) column diameters, up to typically one column diameter below the free surface (Forret et al. 2006).Within that quasi-fully developed region, the self-similarity of the flow structure in the heterogeneous regime -in terms of transverse profiles of void fraction and liquid mean as well as fluctuating velocities -was shown to hold for diameters ranging from 0.15 up to 3 m, and for superficial velocities spanning almost a decade, that is from the transition that arises for V sg of approximately a few cm s −1 up to around 35-40 cm s −1 (Forret et al. 2006;Raimundo et al. 2019).
Let us now consider the structure of the flow organization.Globally, there is no imposed external pressure gradient.Instead, the buoyancy due to bubbles injected at the bottom forces the liquid upward.As there is no liquid flux entering or leaving the column, the liquid must flow downward in some regions of the system.In the homogeneous regime, that occurs in between bubbles everywhere in a given cross-section.In the heterogeneous regime, a stable global recirculation takes place with an upflow region in the centre and a downflow region near the walls.Indeed, the mean liquid velocity profiles consistently exhibit an inversion of the velocity direction at a distance from the column axis equal to 0.7R where R is the column radius.
Such large-scale organization is reminiscent of turbulent buoyancy-driven flows in confined channels with zero mean flow where density gradients are due to solute concentration (Cholemari & Arakeri 2009) or to temperature (Castaing et al. 2017).As for long channels or pipes, the translational invariance observed along the bubble column axis implies that the only characteristic length is the column diameter, provided that the column is high enough.Besides, and owing to that translational invariance, a uniform density gradient develops along the channel (Cholemari & Arakeri 2009;Castaing et al. 2017).That expectation is corroborated in bubble columns.Indeed, in the quasi-fully developed region, the local void fraction measured on the column axis ε axis exhibits a linear growth with the vertical distance above the gas injector, and the slope dε axis /d(H/R) is proportional to the mean void fraction in the column (see figure 26 in Lefebvre et al. 2022).Neglecting the density of the gas compared with that of the liquid (all the experiments discussed hereafter were performed at ambient pressure conditions), the local density ρ of the mixture becomes ρ ∼ (1 − ε)ρ L (with ρ L the density of the liquid phase).As the void fraction increases with the height above the injector, the density decreases with the height.Hence, the vertical stratification observed on the bubble column axis is stable.This contrasts with most investigations made on zero-mean-flow buoyancy-driven turbulence in pipes.For the later, the boundary conditions imposed at the top and bottom ends of the vertical pipe or channel are fixed temperature or solute concentration corresponding to unstable situations, that may lead to intermittent up and down flows (Gibert et al. 2009;Rusaouen et al. 2014).
While in such studies the output is to evaluate the vertical flux of heat or of solute concentration, in bubble columns the boundary conditions are different since it is the gas volumetric flux through the system that is imposed, and the void fraction is the unknown.Owing to the local flow structure evoked above, injected bubbles are entrained in the central portion of the column, and most of them disappear at the free surface.A small fraction of these bubbles (approximately a few per cent, see Lefebvre et al. 2022) recirculate along the walls.These transport mechanisms lead thus to a strong transverse gradient in void fraction that induces radially distributed axial buoyancy forcing.Therefore, in bubble columns, the flow destabilization arises from the radial distribution of the two phases (instead of an unstable vertical stratification).
In the liquid momentum balance, fluid inertia terms (namely ρ∂ t v Li , ρv Lj ∂ j v Li and ∂ i p) equilibrate buoyancy.Within a Boussinesq approximation, the later equals g i ρ, where g is the gravitation acceleration and ρ the difference in density at the origin of buoyancy forces.Hence, the velocity scale for the liquid obeys where ρ/ρ is evaluated at a large length scale.In turbulent buoyancy-driven flows in confined channels, the constant axial gradient of the mean density is used to evaluate ρ (Cholemari & Arakeri 2009;Castaing et al. 2017).In bubble columns, as the flow destabilization arises from lateral differences in density and thus in buoyancy, we sought a relevant scale from the radial void fraction profile.The void fraction typically evolves between ε axis on the column axis and nearly 0 in the wall zone.Thus, the magnitude of the radial difference in density ρ over a length scale of the order of the column diameter D is ρ ∼ D∂ρ/∂r ∼ ρ L D∂ε/∂r ∼ ρ L ε axis .Therefore, the void fraction on the axis ε axis is a measure of the radial density gradient ρ/ρ L .The self-similarity of the radial void fraction profiles ε/ε axis = f (r/R) also supports that result.Such self-similarity is found empirically and therefore the functional form of f (r/R) has not been deduced from first principles yet.Nevertheless, different fits of f (r/R) have been proposed in the literature (see Forret et al. (2006) and discussions therein).
Hence, as ε axis is proportional to the global void fraction (Raimundo et al. 2016), one can use any characteristic gas fraction ε in the system to estimate the magnitude of the radial density gradient ρ/ρ L .Consequently, the scaling from (2.1) becomes (2.2) The Boussinesq approximation is not mandatory for the derivation of (2.2).Indeed, for a bubbly flow, the dynamical equilibrium for the liquid phase balances at first-order inertia terms (i.e.ρ L DV L /Dt) with the momentum transfer between phases.The latter, homogeneous to a force per cubic metre, can be estimated as the void fraction times the force F exerted by a single bubble on the fluid divided by the bubble volume V. Hence, the momentum equation for the liquid writes at first order where the pressure gradient term p includes the hydrostatic contribution.The last term in (2.3) represents the momentum source for the liquid phase due to the presence of bubbles: it is written here with the force −F on a bubble because of the neat scale separation between the disturbance fields (evolving over at the scale close to the bubble size), and the undisturbed fields V L and P (that change over a scale of order D).Given the dynamic equilibrium of the dispersed phase, and since there is no mean vertical acceleration of the continuous phase in the quasi-fully developed region, the force F along a vertical corresponds to the buoyancy on a bubble, i.e. to ρ L g i V.The momentum transfer amounts then to ερ L g i , and the (2.2) scaling is recovered by balancing inertia and buoyancy without using the Boussinesq approximation, that is without constraints on the relative velocity between phases.Former experimental results support the velocity scaling proposed in (2.2).Indeed, in Raimundo et al. (2019), we evaluated the liquid flow rate Q Lup in the core region of the flow that is in the zone where the mean liquid velocity is upward directed.Owing to the self-organization of the flow occurring in the heterogeneous regime, that region extends from the column axis up to a radial distance of 0.7R (that 0.7R limit was also identified by Kawase & Moo-Young 1986).We have shown that, in the heterogeneous regime, Q Lup is independent of the gas superficial velocity.Instead, Q Lup only depends on the column diameter D and it scales as D 5/2 : consequently, (gD) 1/2 was identified as the proper velocity scale for the mean flow circulation.Equation (2.2) is consistent with that result since the void fraction is known to be weakly sensitive (if any) to the column diameter.In the following, we test the relevance of (2.2) for bubble columns in the heterogeneous regime by examining a number of experimental results relative to the liquid and to the gas velocities, including mean and fluctuating components.In § 3, we consider new experiments in which we succeeded to gather reliable statistics on bubble velocity.In § 4, we examine datasets extracted from the literature.patented by A2 Photonic Sensors.The probe design, the signal processing and the sensor qualification are detailed in Lefebvre et al. (2022), where mean bubble velocity profiles in a D = 0.4 m bubble column are also presented.In the following, we exploit further that probe to examine how bubble velocity statistics evolve with the gas superficial velocity.

Test
In parallel, classical Pavlov tubes are also used to access the liquid velocity.Let us first summarize the experimental conditions.

Experimental conditions
The experiment consisted in a 3 m high and D = 0.4 m internal diameter bubble column functioning with air and water.The gas injector was a 10 mm thick Plexiglas plate perforated by 352 orifices of 1 mm internal diameter.These orifices were uniformly distributed over the column cross-section.The column was filled with tap water at an initial height H 0 = 2.02 m.The surface tension of the tap water used was 67 mNm −1 at 25 • C, its pH evolved in the interval [7.7, 7.9] and its conductivity varied within the range 330-450 μS cm −1 , indicating the presence of a significant solid content.All the data presented here were gathered at H = 1.45 m above injection, that is at H/D = 3.625, a position well within the quasi-fully developed region.Besides, and owing to the large ratio H 0 /D = 5.05, the information collected in this zone is not sensitive to the static liquid height H 0 .Experiments were performed for superficial gas velocities V sg ranging from 0.6 cm s −1 to 26 cm s −1 ; the maximum global volume fraction was approximately 35 %.
Information relative to bubbles was acquired with the Doppler probe.For each bubble detected, the probe gives access to the gas residence time t G , that is the time spent by the probe tip inside the bubble.The void fraction is given by the sum of residence times divided by the measuring duration.For the present flow conditions in terms of fluids, bubble size, absolute velocity of bubbles and probe dimension, the uncertainty on void fraction measurements is less than 1 % according to the study of Vejražka et al. (2010) in air-water systems.Besides, a Doppler signal is recorded from the rear interface (that is at the gas to liquid transition) for some bubble signatures, and its analysis provides the bubble velocity V b projected along the fibre axis.When both the gas residence time t Gi and the bubble velocity V bi are available for the ith bubble, one can infer the gas chord C i = V bi t Gi cut by the probe through that bubble.The typical reproducibility on chord length measurements is of order 3 % on the mean value, and less than 20 % on the standard deviation (Lefebvre et al. 2022).
For the liquid phase, velocity statistics were measured with a Pavlov device made of two parallel tubes (external diameter 6 mm, internal diameter 5 mm), each drilled with a 0.5 mm in diameter hole.These two orifices faced opposite directions: they were aligned along a vertical, and the vertical distance between them was 12 mm.The pressure transducer was a Rosemount 2051 CD2 with a dynamics of ±15 000 Pa, a resolution of ±9.75 Pa and a response time of 130 ms.The differential pressure was transformed into the local liquid velocity using v 2 L (t) = ±2|p(t)|/ρ L (no correction dependent on void fraction was considered) with the appropriate sign.Hence, the range in velocity was ±5.48 m s −1 and the resolution ±0.14 m s −1 .
Liquid and gas velocity probability density functions (p.d.f.s) measured with these sensors in the centre of the column are illustrated figure 2. By construction, the Pavlov tube detects both positive (upward directed, i.e. against gravity) and negative (downward directed, i.e. along gravity) velocities.For bubbles, as the Doppler probe detects only inclusions approaching it head on, the p.d.f.s were built by cumulating the information gathered over the same measuring duration and at the same position with an upward directed probe and with a downward directed probe.More details and discussion concerning these measurements are presented in Lefebvre et al. (2022).The reproducibility on mean bubble velocity and on its standard deviation is better than ±5 %.
Concerning the bubble size, the analysis of the axial evolution of chords distributions along the column indicates that coalescence was absent (or at least extremely weak) in our experimental conditions (Lefebvre et al. 2022), most probably because of the partial contamination of the tap water used.Over the investigated range of superficial velocities, the Sauter mean vertical diameter of bubbles remained in the interval [6.2 mm; 6.7 mm] while their Sauter mean horizontal diameter measured with the correlation technique (Raimundo et al. 2016) increased with V sg from 6.6 to 7.8 mm.Overall, the mean equivalent bubble diameter remained in the interval [6.62 mm; 7.35 mm]: that corresponds to a terminal velocity from 21 cm s −1 to 23 cm s −1 (Maxworthy et al. 1996) and to particle Reynolds numbers in the range 1450-1550.
The local void fraction measured on the column axis ε axis is plotted vs the gas superficial velocity in figure 3(a).The homogeneous regime ends for V sg between ∼4 cm s −1 and ∼5 cm s −1 , while the 'pure' heterogeneous regime starts at V sg ∼ 6.5 cm s −1 .Following Krishna et al. (1991), these data are plotted as V sg /ε axis vs V sg in figure 3(b): they exhibit a constant rise velocity, close to the bubble terminal velocity U T , up to V sg ∼ 4 cm s −1 , that is within the homogeneous regime.Beyond that, the apparent rise velocity (called 'rise velocity of swarm' by Krishna et al. 1991), monotonically increases with the gas superficial velocity.It reaches a magnitude of approximately 3 times U T at the largest V sg investigated here (namely 24.7 cm s −1 ).That increase is the signature of the heterogeneous regime.Note also that the latter correspond to a void fraction on the axis that exceeds 20 %.In the following, the transition will be represented by a vertical dash line at V sg = 5 cm s −1 in figures as a guide to the eye.

Local gas and liquid velocity on the column axis vs V sg
Figure 4 provides the mean vertical velocities of bubbles V G and of the liquid V L on the column axis as well as the standard deviations V G for the gas phase and V L for the liquid phase.Two datasets are presented for the mean velocities.
For bubbles, a set named 'up flow' corresponds to measurements achieved with a Doppler probe pointing downwards that collects only upward directed (i.e.positive) vertical velocities.A second dataset named 'up and down flow' was obtained by gathering direct (i.e.without interpolation) velocity measurements from a probe pointing downward, with direct (i.e.without interpolation) velocity measurements from the same probe pointing upward.In this process, the measuring duration was the same for the two probe orientations.In the flow conditions considered here, the mean velocities from these two Buoyancy-driven bubbly flows: velocity scaling sets are close, with a difference of at most 4 % (Lefebvre et al. 2022).Similarly, the difference on bubble velocity standard deviations from these two sets is at most 8 %.
For the liquid velocities, two datasets are also presented: one corresponds to moments evaluated over the entire distribution (named 'up and down flow') while the other concerns positive velocities only (named 'up flow').In the heterogeneous regime, the difference between the two sets is at most 3.6 % for the mean value and 18 % for the standard deviation.Oddly, larger liquid velocity deviations between 'up flow' and 'up and down flow' statistics appear in the homogeneous regime.The difference is especially pronounced for V sg below 3 cm s −1 .These deviations are related with the unexpected apparition of a significant negative tail in the liquid velocity p.d.f.s when V sg becomes small, a defect that may possibly be due to the flow perturbation induced by the rather large probe holder used in our experimental set-up.
Except for these low V sg cases, the differences between the 'up flow' and the 'up and down flow' datasets remain weak.These small differences partly come from the fact that these data are all collected on the column axis, where the probability of occurrence of absolute negative velocities in the laboratory frame remains small.Indeed, in our experiments, the probability to observe a downward directed liquid velocity on the axis was less than 3 % for any V sg in the heterogeneous regime.Similarly, Xue et al. (2008) found a probability for observing negative bubble velocities on the column axis between 4 % and 4.5 % at V sg = 14 cm s −1 and approximately 6 % at V sg = 60 cm s −1 .For the gas phase, and because of the positive (upward directed) relative velocity, these probabilities should be lower than the above figures.Hence, considering either 'up flow' or 'up and down flow' datasets does not change the conclusions proposed hereafter.Yet, the distinction between the two series is worth to be kept in mind in the perspective of analysing other radial positions where the probability of occurrence of downward flow increases.
From the data presented figure 4, a series of comments and conclusions emerge: (i) The local relative velocity V G − V L remains nearly constant in the homogeneous regime.It is approximately 27 cm s −1 , a magnitude comparable to the 'mean' terminal velocity U T of the bubbles generated in the column.(ii) In the heterogeneous regime, the relative velocity becomes larger than U T .
Furthermore, it seems to monotonically increase with V sg (see the trend indicated by the green and red dashed lines in figure 4).At approximately V sg = 19.5 m s −1 , the measured relative velocity amounts to 54.6 cm s −1 that is 2.5 times the terminal velocity U T .Thus, these direct velocity measurements are consistent with the behaviour of the 'rise velocity of swarm' V sg /ε axis shown in figure 3(b).They also confirm the conclusion we previously obtained (see Raimundo et al. 2019) by analysing the flow in the core region of the bubble column: the apparent relative velocity was indeed found to range between 2 and 8 times U T using a conservative evaluation of the gas liquid flow rate in the core region of the bubble column.(iii) The relative fluctuations in velocity V /V are nearly constant in the heterogeneous regime (as shown figure 5).For the liquid phase, the average of V L /V L over the data collected in the heterogeneous regime equals 36.5 %, in agreement with previous findings (see the discussion in Raimundo et al. 2019).For the gas phase, the average of V G /V G is even higher; it equals 52 %-53 % when considering positive velocities only, and it rises up to 56 %-57 % when combining positive and negative velocities measurements (Lefebvre et al. 2022).These strong figures confirm that intense turbulent motions take place in heterogeneous conditions.(iv) From a closer examination of figures 4 and 5, two different behaviours could possibly be distinguished in the heterogeneous regime.From the homogeneous/heterogeneous transition up to V sg approximately 13-15 cm s −1 , the relative velocity V G − V L and also to some extent the fluctuation V G /V G exhibit a clear monotonic increase with the gas superficial velocity.Above V sg ∼ 13-15 cm s −1 , these two quantities seem to become constant.In particular, the increase in the bubble velocity illustrated by the black dashed line in figure 4 becomes nearly parallel to that of the mean liquid velocity (dashed red line in figure 4): accordingly, the relative velocity seems to stabilize at a value approximately 2.3-2.4UT at large V sg .With regard to flow dynamics and scaling laws, it would be worthwhile to clarify whether the relative velocity reaches an asymptote, or if it continues to grow with V sg : more experimental data covering an enlarged range of gas superficial velocities are required for that.
To test the scaling proposed in (2.2) on these experimental data, we used local void fraction and velocities measured on the axis and at the same height H/D above the injector.Since the two sets 'up flow' and 'up and down flow' are close, let us consider only 'up and down flow' velocity statistics for the analysis.The mean velocities as well as the standard deviations scaled by (gDε) 1/2 are represented for both phases in figure 6.Clearly, all these quantities remain nearly constant for all flow conditions pertaining to the heterogeneous regime.For the mean bubble velocity on the column axis, one gets while for the mean liquid velocity on the column axis According to these results, the relative velocity on the column axis scales as U R ∼ 0.41(gDε) 1/2 , i.e. it monotonically increases with the void fraction ε.Such an increase of the relative velocity with the void fraction has been identified in bubble columns using indirect arguments (see for e.g.Raimundo et al. 2019).It has also been observed in other  gas-liquid systems.Such a behaviour is sometimes represented by a swarm coefficient that quantifies the decrease of the drag force acting on a bubble with the void fraction (Ishii & Zuber 1979;Simonnet et al. 2007).Nowadays, ad hoc swarm coefficients are routinely introduced in numerical simulations based on two-fluid models (e.g.McClure et al. 2017;Gemello et al. 2018).We bring here direct experimental evidence of the increase of the relative velocity with the void fraction in a bubble column operated in the heterogeneous regime.
Concerning the standard deviation of velocities, the standard deviation being proportional to the mean (see figure 4), they also follow the same scaling with V G /(gDε) 1/2 ∼ 0.6 ± 0.02 and V L /(gDε) 1/2 ∼ 0.22 ± 0.02.Hence, all the above results obtained on the mean and on the fluctuating components of bubbles and liquid velocities confirm the soundness of the scaling proposed in (2.2) with respect to void fraction.
The scaling resulting from an inertia-buoyancy equilibrium proposed above also includes an increase of velocities with the square root of the bubble column diameter.As the experiments presented in this section concern only one bubble column diameter, the dependency with the column diameter cannot be tested.New experimental data gathered in bubble columns of variable diameters are needed to test the validity of the proposed scale.In that perspective, available experiments from the literature and relative to bubble columns of variable diameter will be analysed in the next section.
Nevertheless, we already have accumulated strong experimental evidence of the relevance of the square root of the bubble column diameter as a scaling factor for the mean liquid velocity (Raimundo et al. 2019).Indeed, the neat upward liquid flow rate Q Lup in the column, where Q Lup is obtained by integrating the liquid flux (1 − ε)V L over the core region (i.e. from the axis up to 0.7-0.71R),happens not to depend on V sg in the heterogeneous regime.In the present experiments, we also found that Q Lup is independent on V sg in the heterogeneous regime.That conclusion was confirmed from data collected at three heights above injection, namely H/D = 2.625, H/D = 3.625 and H/D = 4.875.The result is Q Lup = 0.0122 m 3 s −1 , with variations between +0.001 m 3 s −1 and 0.0007 m 3 s −1 depending on the set of data considered to evaluate the mean.Moreover, we have previously shown (Raimundo et al. 2019) that Q Lup is proportional to D 2 (gD) 1/2 over a significant range of column diameters (from D = 0.15 m to 3 m) and of gas superficial velocities (from V sg = 9 to 25 cm s −1 ).As shown figure 7(a), the present experiments confirm that finding in a D = 0.4 m column, for V sg between 6.5 cm s −1 and 22.7 cm s −1 and for 2.625 H/D 4.875.The dependency of Q Lup with D is further illustrated in figure 7(b) where we have reported the present data, the data collected by Raimundo et al. (2019), as well as one data produced by Guan et al. (2015) in the following flow condition: V sg = 47 cm s −1 in a D = 0.8 m column and for 2.75 H/D 4.625.All these data fall onto the same curve.Overall, the observed liquid flow rate -column diameter relationship writes This result can also be expressed as a Froude number based on the average liquid velocity Q Lup /S core in the core region.Here, the cross-section S core of the ascending flow zone is evaluated as πD 2 /8 (the mean liquid becomes zero at a radial position between 0.7R and 0.71R: that limit is well approximated as 2 1/2 /2R = 0.707R).Hence, Fr L = (Q Lup /S core )/(gD) 1/2 = 0.098.(There is a typo error in (13) of Raimundo et al. (2019): the coefficient 0.024 should be replaced by 0.098.)All the above mentioned experiments bring clear evidence that the velocity of the mean liquid circulation in a bubble column operated in the heterogeneous regime scales as (gD) 1/2 with the column diameter.

Confrontation of the proposed scaling with experimental data from the literature
In this section, we examine whether available experimental data on local liquid or gas/bubble velocities obey the scaling proposed above when considering a larger range of flow conditions, and in particular variable column diameters.To this end, we target meaningful experiments in bubble columns in the sense that we seek conditions such that the flow dynamics is controlled by the same mechanisms as those discussed in § 2.More precisely, we select experiments pertaining to the 'pure' heterogeneous regime, a regime that occurs well beyond the transition region, and for which the global void fraction is an increasing and concave function of the gas superficial velocity.In addition, we consider only data collected in the quasi-fully developed region of the bubble column: as discussed in § 2, that condition implies a minimum column diameter, a minimum static liquid height as well as a proper range of measuring heights.
Besides, and although it is known that in the heterogeneous regime the flow organization is weakly sensitive to the injector design, we also select experiments such that the gas injection is pretty well uniformly distributed over the column cross-section to avoid forcing of large-scale instabilities by an uneven gas repartition at injection and/or to avoid strong asymmetry of the mean flow such as the one observed by Chen et al. (2003) in their largest column.
Another constraint in the search of relevant data is that both local velocity and local void fraction should be simultaneously available.In the following, we focus on local data gathered on the column axis.The sets of data fulfilling the above-mentioned constraints are presented in tables 1 and 2 for the liquid phase and in tables 3 and 4 for the gas phase.Note that almost all experimental conditions in tables 1, 2, 3 and 4 correspond to air-water systems and to 'large' bubbles, i.e. with an equivalent diameter between 3 and 10 mm.Their terminal velocity typically ranges between 21 and 27 cm s −1 (Maxworthy et al. 1996), so that all these flow conditions involve bubble dynamics at high (from 800 to 2100) particle Reynolds number.
Yet, these experiments remain difficult to quantify with respect to coalescence.We qualitatively estimated the coalescence efficiency based on measurements of the axial evolution of the bubble size when such data were available.When such information was absent, we considered how the bubble size changes with the superficial velocity: a strong increase of the latter from homogeneous to heterogeneous conditions could be (but this is not certain) the mark of a neat coalescence.In addition, let us underline that, when a significant coalescence is present, the flow regime may continuously evolve with the height above injection so that a quasi-fully developed region may not exist or may require column heights much larger than the ones available in standard experiments.All the information collected is summarized in tables 1, 2, 3 and 4, where the situation with respect to coalescence has been classified into three main categories: none or weak coalescence, medium coalescence, strong coalescence, and the situation is said to be unclear when information was insufficient to draw a conclusion.

Mean liquid velocity on the column axis
We start this section by briefly detailing the datasets from the literature that will be discussed.Tables 1 and 2 lists the experiments that obey the above constraints and that provide the liquid velocity and the local void fraction on the bubble column axis.Some choices were made to exploit the data.For Hills (1974), we considered only the data acquired with the 'plate B' that corresponds to a uniform gas injection over the column cross-section.In Forret et al. (2006), the only quantitative data on gas hold-up are global void fractions deduced from static and dynamic liquid heights.We transformed the global void fraction into a local void fraction on the axis by multiplying it by 1.5 as done by the authors (see their (4)), but this factor could be inappropriate.The local void fractions for the data of Vial et al. (2001) were collected in Camarasa et al. (1999).For Yao et al. (1991), we present the data they collected at various heights with H/D from 2.6 up to 12, and we use an extrapolation to evaluate the void fraction at V sg = 10 cm s −1 .Note also that these authors imposed a forced liquid motion but the mean velocity of 1 cm s Buoyancy-driven bubbly flows: velocity scaling negligible compared with the measured liquid and gas velocities on the column axis: these data are therefore believed to be representative of a bubble column operated in a batch mode and they have been kept in the analysis.Finally, all the measurements mentioned in tables 1 and 2 considered positive and negative velocity realizations, although some (hard to evaluate) bias may be present notably with Pavlov tubes, as evoked by Hills (1974).To be consistent, we compared with our data series named 'up and down flow' (see § 3).
Figure 8 provides the quantity V L /(gDε) 1/2 vs the superficial gas velocity with V L and ε measured on the axis.The figure includes all the contributions listed in tables 1 and 2, and that have been collected in the homogeneous regime as well as in the heterogeneous regime.The data concern column diameters from 0.1 to 3 m and superficial velocities varying from 1.2 to 62 cm s −1 .In the limit of small superficial gas velocities, V L /(gDε) 1/2 evolves between 0.01 and 1.02.As V sg increases, the range of variation of V L /(gDε) 1/2 smoothly diminishes.Above V sg ∼ 15 cm s −1 , V L /(gDε) 1/2 evolves within the interval [0.48; 0.77].Note that this interval corresponds to bubble column diameters ranging from 0.1 m to 3 m.Above V sg ∼ 20 cm s −1 , that interval further narrows and the quantity V L /(gDε) 1/2 tends to be constant in the pure heterogeneous regime: that feature supports the scaling argument presented in § 2. Although experimental data are lacking at very large V sg to precisely define the asymptotic behaviour, the latter is estimated as as shown by the inset in figure 8 that provides V L vs (gDε) 1/2 for all available data at V sg above 8 cm s −1 .The proportionality factor equals 0.5774 (and the correlation coefficient is 0.867).Note that the same plot using all data available at V sg above 6 cm s −1 also provides a linear behaviour with a proportionality factor equals to 0.5786 (and a correlation coefficient of 0.87).Alternately, when the quantity V L /(gD) is plotted vs the local void fraction for all the data shown in figure 8, V L /(gD) is found to evolve as ε 0.50±0.01: this is fully consistent with the ε 1/2 dependency predicted.It is worth discussing the uncertainty on these figures as the data presented in figure 8 come from different operators and from various measuring techniques.The typical dispersion in the measurements can be appreciated by comparing the data collected in identical bubble columns.For example, in the D = 0.15 m column and at V sg ∼ 20 cm s −1 , there is a 0.2 difference in V L /(gDε) 1/2 between the data from Forret et al. (2006) and those from Raimundo et al. (2019).Similarly, in the D = 0.4 m column and at V sg ∼ 20 cm s −1 , the values of V L /(gDε) 1/2 deduced from the data of Forret et al. (2006), from those of Raimundo et al. (2019) and from the present work all fall within a 0.1 band.These are quite reasonable dispersions especially when considering that the water quality was not always the same, so that the coalescence efficiency varied.Finally, as far as we can judge from the available information given in articles, all the experimental conditions in figure 8 correspond to no or weak coalescence.The monotonic allure of the evolution of the local void fraction on the axis with V sg shown in figure 9 also supports that statement.

Mean gas velocity on the column axis
Experiments providing statistics on the bubble velocity are not common.This is due to the lack of reliable measuring techniques giving access to bubble velocity in the difficult conditions encountered in the heterogeneous regime, in particular with respect to high void fractions, flow unsteadiness and 'chaotic' 3-D trajectories of bubbles.Each bubble velocity technique has its own limitations, and their respective uncertainty and resolution in such flow conditions are not well known.For example, for phase detection techniques based on immersed probes either single, double or multiple, it is well known that erroneous velocity data are collected in heterogeneous conditions because of the unsteady 3-D motions of these two-phase flows.Yet, average quantities relative to velocity or size seem to be meaningful (Chaumat et al. 2007;Raimundo et al. 2016).Moreover, while some bias may alter the statistics and when the latter occurs, it is usually hard to quantify.For all these reasons, we will not discuss further the respective merits of each measuring technique: instead, we report all available raw data as given in the original papers, keeping in mind that some issues on resolution, accuracy or bias remain as an open question.V sg (m s -1 ) 0.4 0.5 0.6 0.7 0.8 Figure 10.Evolution of V G /(gDε) 1/2 where V L and ε are measured on the column axis vs the superficial gas velocity V sg for all gas or bubble velocity measurements from the articles quoted in tables 3 and 4.
Tables 3 and 4 summarize the set of experiments used to analyse the gas velocity that corresponds here to the velocity of bubbles.The corresponding data on the mean bubble velocity scaled by (gDε) 1/2 , where both the velocity and the void fraction were measured on the column axis, are reported vs V sg in figure 10 irrespective of the flow regime.
We will now provide some information on how the data were exploited for the datasets used in this section.First, when different injectors were tested, we always selected the data acquired with multiple-orifice injectors distributed over the entire column cross-section.Second, original data were sometimes interpolated to estimate missing local void fraction data: that process was used only when the interpolation process was safe.If some extrapolation was required, the corresponding data were discarded unless otherwise specifically stated in the text in the legend.It is worth to notice that almost all data correspond to tap water and air (under ambient pressure and temperature conditions), with two exceptions: Yao et al. (1991) used deionized water, and one set of data from Camarasa et al. (1999) was gathered in an aqueous solution of alcohol (water and pentanol at a concentration 4 × 10 −4 mol l −1 ).
We also mention some specific choices we made when extracting the data.For Yao et al. (1991), only the bubble velocity detected by the ultrasound technique is plotted because these authors found comparable results with their five-point conductivity probe.In Camarasa et al. (1999), the bubble velocity was measured by an ultrasonic Doppler technique for single orifice and porous plate gas injection but not for the multiple orifice nozzle that provides the injection conditions we are looking for here.For that sparger, the bubble velocity was measured by a DGD (dynamic gas disengagement) technique: Camarasa et al. (1999) report the velocity of 'large' bubbles and of 'small bubbles' but they do not explicitly specify how 'large' bubbles are distinguished from 'small' ones.It happens that, for a given flow condition, the mean velocity of 'small bubbles' measured by Camarasa et al. (1999) is 2 to 4 times (in water) and 2 to 2.5 times (in aqueous solution of alcohol) lower than the mean velocity measured for 'large bubbles'.Presumably, 'small' bubbles correspond to the 1-2 mm bubbles at the lower end of the bubble size distributions they provide, while 'large' bubble can be as large as 8-10 mm in water (see their figure 10 and 6-8 mm in water plus pentanol (see their figure 20).According to these comments, only the data for 'large' bubbles are reported in figure 10.
Let us finally mention that, for the data of Camarasa et al. (1999) gathered in the aqueous solution of pentanol, no data are given on the local void fraction and we used the global void fraction instead.Accordingly, the values of V G /(gDε) 1/2 are overestimated for that series.Chen et al. (2003) performed measurements in a D = 0.4 m bubble column at H/D = 2.6, and in a D = 0.8 m bubble column at H/D = 1.3: the latter case, for which they observed non-symmetrical flows, was discarded because the data were not collected in the quasi-fully developed region.
Concerning the experiments by Xue et al. (2008), all the data they collected in the heterogeneous regime at H/D = 5.1 correspond to strongly asymmetrical flows.Indeed, the velocity difference between the upward motion on one side of the column and the downward motion on the opposite side ranges between 20 cm s −1 and 40 cm s −1 : these figures are therefore quite significant compared with the mean bubble motion on the column axis that evolved between 40 and 90 cm s −1 .That asymmetry was further confirmed by void fraction and bubble detection frequency profiles.The only symmetrical bubble velocity profile reported by Xue et al. (2008) was collected at a larger distance from injection (namely H/D = 8.5) and at V sg = 30 cm s −1 .Unfortunately, they do not provide the void fraction for these conditions.Despite these shortcomings, all the data of Xue et al. (2008) collected at H/D = 5.1 have been integrated in our analysis.Let us also underline that these authors are among the very few who have explored large gas superficial velocities.
Among the contributions listed in tables 3 and 4, almost all gathered (at least in principle) positive and negative bubble velocities, except for Raimundo (2015) and Raimundo et al. (2016) who collected positive velocities only since they exploited the dewetting of a single fibre tip.Yet, as the sources of bias are usually not analysed for bubble velocity measurement techniques, it is difficult to ascertain that the information collected was indeed the faithful assembly of positive and negative realizations.An indication of these difficulties is that the standard deviation of bubble size distributions are never provided, nor discussed, except by Yao et al. (1991) who measured bubble velocity fluctuations with an ultrasound technique.As for the liquid phase, and to be consistent, we thus consider our data series from § 3 named 'up and down flow' for the comparison.
The data presented figure 10 concern column diameters from 0.1 to 1 m and gas superficial velocities varying over two decades from 0.6 to 60 cm s −1 .At low gas superficial velocities, say below a few cm s −1 , that is within the homogeneous regime, the quantity V G /(gDε) 1/2 evolves from 0.2 up to 1.8.That ratio significantly varies from one experiment to the other.For a given data series, the ratio V G /(gDε) 1/2 tends to become somewhat constant when moving towards large superficial velocities.For V sg above approximately 10 cm s −1 , that is well within the heterogeneous regime, the dispersion of the data significantly diminishes and V G /(gDε) 1/2 evolves inside a narrower band between 0.6 and 1.4.Note that these last figures encompass bubble column diameters ranging from 0.15 to 1 m.The trends are therefore the same as those detected when analysing the liquid velocity.Yet, the fluctuations observed from one series to another are larger than those observed for the liquid velocity.This is probably because of the stronger uncertainties of gas velocity measuring techniques.Also, and compared with the mean liquid velocity presented figure 8, it is more difficult to estimate the asymptotic value of V G /(gDε) 1/2 .According to Xue et al. (2008) and to some runs of Raimundo et al. (2019), that limit is near 0.8 while from the present data, as well as from those of Yao et al. (1991), the limit is possibly closer to 1. Again, experimental data at very large V sg are required to more accurately determine the asymptotic behaviour of V G /(gDε) 1/2 .Despite the limitations on available data, the trends observed on V G /(gDε) 1/2 for column diameters between 0.1 and 1 m are consistent with the scaling argument proposed in § 2.
Figure 11 provides the evolution of the void fraction on the axis with V sg for all the experiments quoted in tables 3 and 4. As in figure 9, all these evolutions are monotonic, so that they presumably all correspond to no or weak coalescence.Yet, in terms of local void fraction, one data series happens to be neatly above all others: this is the one from Xue et al. (2008).Contrary to all other experiments presented in figure 10 (see also tables 3 and 4), the coalescence was probably significant in their experiments.Indeed, they observe an increase of the mean bubble chord length on the axis from 2-3 mm in the homogeneous regime up to 6 mm in the heterogeneous regime while the standard deviation of the size distribution growths from approximately 1 mm up to 10 mm.Also, their bubble chord distributions indicate that bubbles up to 15 mm are detected, and bubbles on the axis are significantly larger than near walls where their mean chord is less than 3 mm.Let us also underline that the local void fractions measured by Xue et al. (2008) in the heterogeneous regime are among the largest of all data series presented figures 9 and 11.Xue et al. (2008) used injectors made of orifices 0.5 mm in diameter, and small orifices are known to produce smaller bubbles and thus to increase the maximum void fraction reached in the homogeneous regime (Joshi et al. 1998).Chaumat, Billet & Delmas (2006) obtained similar large values of the local hold-up when using the same small orifices: yet, the values of V G /(gDε) 1/2 deduced from their measurements remain comparable to other contributions for V sg below 0.13 m s −1 (figure 10).These comments indicate that a characterization of flow conditions with respect to coalescence is not easy.In addition, experimental data on bubble velocity are missing to evaluate how much the magnitude of V G /(gDε) 1/2 may vary with the coalescence efficiency.
So far, the scaling of the liquid velocity has been discussed based on data collected on the axis.For the liquid, it is known that, in the heterogeneous regime, both the transverse liquid velocity and the local void fraction profiles assume self-similar shapes when scaled by their respective value on the axis (Forret et al. 2006).Hence, all the above findings are expected to remain valid at any radial position in the column provided one remains in the quasi-fully developed region.For the gas phase, we have shown in Lefebvre et al. (2022) that the bubble velocity profiles collected in the quasi-fully developed region when in heterogeneous conditions also happen to be self-similar when scaled by the bubble velocity on the column axis.Hence, the proposed scaling is also expected to remain valid at any radial position for the gas phase.

Liquid and gas velocity fluctuations
Experimental data on velocity fluctuations collected in bubble columns are scarce, and this is particularly true in the heterogeneous regime.For the liquid phase, Menzel et al. (1990) provide two profiles of the axial liquid velocity fluctuations (quantified here by the standard deviation V of the velocity distribution).Nevertheless, these datasets were gathered in a 80 wt% glycerol/water mixture and, unfortunately, the authors do not indicate the corresponding mean velocities and void fractions for that fluid.Otherwise, data for the liquid phase are available in Yao et al. (1991), Vial et al. (2001), Forret (2003), Raimundo et al. (2019) and from the present contribution.All these data concern deionized or tap water.
For the gas phase, single tip or multiple tips probes that exploit a transit time technique for velocity measurements are common, but, with these techniques, the measured velocity distributions are strongly biased by the detection of erroneous, large velocities (see e.g.Chaumat et al. 2007;Raimundo et al. 2016).It happens that the mean velocity is significant, but the standard deviation is not reliable.Also, some authors (such as Chen et al. 2003;Xue et al. 2008;Guan & Yang 2017) provide velocity distributions but the standard deviations are not quantified.For these reasons, one is left only with the data from Yao et al. (1991) acquired from an ultrasound technique in a D = 0.29 m column with deionized water/air as fluids, and the data we collected with the Doppler probe in a D = 0.4 m column with tap water/air as fluids (see Lefebvre et al. 2022).
The experiments of Menzel et al. (1990) indicate that, in the heterogeneous regime, the radial profiles of liquid velocity fluctuations remain self-similar when normalized by their maximum.One can also use the velocity fluctuations evaluated on the axis of the bubble column for that normalization.Hence, in the following, we focus our discussion on velocity fluctuations measured on the axis.The relative fluctuation V L /V L in the liquid phase (respectively V G /V G for the gas phase) measured on the column axis is presented figure 12 (respectively figure 13) vs the superficial gas velocity.All these measurements have been done in the quasi-fully developed region.Despite the limited number of independent data, each of the quantities V L /V L and V G /V G tends towards a constant value when moving inside the heterogeneous regime.That feature is well established for the liquid phase since the data come from different operators and from different sensors.Remarkably, the asymptotic behaviour is the same whatever the column diameter ranging from 0.15 to 3 m.The only series exhibiting a different trend are those from Yao et al. (1991) for the liquid phase.These data were obtained from hot-film probes, a technique that could be delicate to exploit in bubbly flows.The difficulties are expected to be even stronger in the conditions encountered in bubble columns at high gas superficial velocities.Unfortunately, Yao et al. (1991) do not comment on the signals they collected, nor on the signal processing they develop: it is therefore difficult to evaluate the reliability of their measurements.
Figures 12 and 13 show that, in the heterogeneous regime, the velocity fluctuations in the liquid as well as in the gas remain proportional to the mean velocity.Hence, all the findings on the scaling of mean velocities in the core region of the bubble column also apply to velocity fluctuations.One may be puzzled by such result, but physical arguments V sg (m s -1 ) . Evolution of the relative fluctuation V G /V G of the velocity of bubbles measured on the column axis vs the superficial gas velocity.
similar to those evoked in § 2 can explain that feature.The idea is that, in the heterogeneous regime, the contributions to velocity fluctuations in the liquid phase by the relative motion at the bubble scale, i.e. the so-called bubble-induced turbulence, which is at the origin of large ratio V L /V L observed figure 12 in the homogeneous regime (see Risso 2018), is not the leading mechanism.Instead, in the heterogeneous regime, the agitation in the liquid is due to the presence of meso-scale structures.The later have been put in evidence and quantified with a 1-D Voronoï analysis performed on the phase indicator function delivered by optical probes.With such Voronoï tessellations, we have shown that in the heterogeneous regime these flows as organized in clusters (high void fraction regions) and voids (low void fraction regions) (Raimundo et al. 2019).The variations in void fraction from one meso-scale structure to the other induce buoyancy forces that spatially fluctuate, and hence a velocity field that changes from one structure to the other.These meso-scale structures are transported by the mean flow and they could also form and disappear.Y. Mezui, M. Obligado and A. Cartellier Hence, at a fixed point in space, the passage of successive structures induces the velocity fluctuations that are precisely those detected by an Eulerian measuring technique, i.e. they are the quantities V L and V G measured with local probes.

Velocity scaling: further considerations on the void fraction prediction
We have shown that, in the heterogeneous regime, all velocities, namely the mean velocity and its standard deviation in both phases as well as relative velocity between phases, scale with (gDε) 1/2 as expected from an inertia-buoyancy equilibrium.We have shown that this scaling holds in a D = 0.4 m column for all the flow conditions pertaining to the heterogeneous regime that we have investigated.The analysis of literature data has confirmed the validity of the proposed scaling for the mean and for the standard deviation of the liquid velocity in columns of diameter between 0.1 and 3 m, for the mean gas velocity in columns of diameter between 0.1 and 1 m, and for the standard deviation of the gas velocity in columns of diameter D = 0.29 and 0.4 m.For the relative velocity, the only data available concern the D = 0.4 m bubble column considered here.Let us recall that almost all flow conditions analysed correspond to air-water systems with large bubbles at high particle Reynolds numbers ( § 4).
To reach a fully predictive status for velocities, one also needs the gas hold-up ε as a function of flow parameters.However, there is no consensus on the void fraction prediction in bubble columns operated in the heterogeneous regime as tens of different correlations involving various sets of parameters are proposed in the literature (as notably shown by the reviews of Joshi et al. 1998;Kantarci et al. 2005;Kikukawa 2017;Besagni et al. 2018).There is even no clear consensus on the set of non-dimensional parameters governing the response of the system.We propose in Appendix A a dimensional analysis dedicated to the heterogeneous regime that is restricted to high aspect ratio bubble columns, to systems far from critical conditions, and without or with weak coalescence.We identify five independent non-dimensional parameters, and a possible choice could be: (i) The Archimedes number Ar = gD 3 /ν 2 L , Ar is the square of a Reynolds number based on the column diameter D, on the liquid viscosity and on the velocity scale (gD) 1/2 .(ii) The Froude number Fr = V sg /(gD) 1/2 , that quantifies the injected gas flow rate.(iii) The Eötvös number Eo = ρ L gd 2 /σ , that measures the mean bubble size d relative to the capillary length.(iv) The Morton number M = gμ 4 L /(ρ L σ 3 ), that involves the physical properties of the couple of fluids selected and the gravitational acceleration.(v) A non-dimensional parameter quantifying the degree of polydispersity in the system defined as the standard deviation of the bubble size distribution std(d) divided by the mean bubble diameter d.
Among these, the Eötvös and Morton numbers completely define the dynamics, that is the shape, the trajectory and the relative velocity, of an isolated bubble having the mean equivalent diameter d immersed in the stagnant liquid and for the given gravitation field (strictly speaking, this only holds for clean interfaces).
All the experimental conditions analysed here (see § 4) involve the same couple of fluids (i.e.air and water in ambient thermodynamic conditions) and Earth's gravity so the M parameter is the same (∼10 −11 ).The Eo parameter evolves in a rather narrow range, roughly from 1 to 10.For these M and Eo values, the bubbles are in the so-called wobbling regime, and they have a similar dynamics with O( 10 Buoyancy-driven bubbly flows: velocity scaling bubble size distributions, it does not change much (the minimum bubble size is typically of the order of 0.5-1 mm while the maximum bubble size never exceeds ∼10 mm).Hence, in the experiments quoted in tables 1, 2, 3 and 4, only the two parameters Ar and Fr have been significantly changed.The Ar number evolves between 9.8 × 10 9 and 2.6 × 10 14 .
Furthermore, to be sure to analyse data pertaining to the heterogeneous regime (and thus to escape from the transition zone), let us consider gas superficial velocities above 7 cm s −1 or 9 cm s −1 .The corresponding Froude numbers span the range [0.016; 0.475].According to figures 7 and 9, the void fraction seems to be mainly driven by the superficial gas velocity V sg .To check that, we first attempted a correlation with both Fr and Ar, and an exponent as low as 0.047 was found for the Ar number.We therefore examined how the void fraction evolves with the Froude number alone.As shown figure 14, the local void fraction correlates well with Fr as one gets ε axis = 0.853Fr 0.389 , V sg 7 cm s −1 , ε axis = 0.838Fr 0.377 , V sg 9 cm s −1 , (5.1) with correlation coefficients of approximately 0.8.The maximum deviation of these fits from measurements is ±30 % except for two data collected in a D = 3 m column at V sg = 16 cm s −1 and 25 cm s −1 for which the deviation reaches 35 %.The measurements in large columns are not easy (probably due to vibration of the probe holder, and/or to flow perturbation induced by the latter when it is too large).If these two data points are discarded, the correlation becomes ε axis = 0.897Fr 0.415 , ( with a correlation coefficient of 0.887.Equation (5.2) holds for V sg 7 cm s −1 as well as for V sg 9 cm s −1 .The deviation remains then within ±22 % for all data.The correlation coefficients as well as the maximum deviations found here appear as acceptable, especially if one accounts for the fact that the data considered were collected at different heights above injection (see tables 1, 2, 3 and 4).Indeed, we have shown in Lefebvre et al. (2022) that the local void fraction on the axis linearly increases with the height in the quasi-fully developed region of the flow.The impact of the measuring height on the local void fraction is illustrated in figure 14 by two datasets (closed symbols) gathered on the axis of the D = 0.4 m column: the upper dataset corresponds to H/D = 6.37 and the lower dataset corresponds to H/D = 2.85.Clearly, the distance between these series is comparable to the dispersion.The disappearance of the Archimedes number in the above empirical expressions for the void fraction is not unexpected.Indeed, as discussed in Appendix A, the heterogeneous regime in a bubble column corresponds to a turbulent regime in free thermal convection.Thus, the Archimedes number somewhat controls the transition to that regime (the critical Rayleigh number introduced by Ruzicka & Thomas (2003) to identify the homogeneous-heterogeneous transition is proportional to Ar).However, once the turbulent regime is installed, buoyancy forcing overwhelms viscous effects, and the precise value of Ar is no longer relevant for setting the dynamical equilibrium: this is why its outcome, i.e. the void fraction, is no longer dependent on Ar.
There is also a debate as whether the void fraction should depend or not on the bubble column diameter D. For example, according to Besagni et al. (2018), the correlation proposed by Akita & Yoshida (1973) should be considered as the state of the art for determining the global gas hold-up.That correlation does not include any dependency of the void fraction on D. Yet, among the experiments quoted in tables 1, 2, 3 and 4 differences (with factors much larger than one) appear between gas hold-up measurements and predictions using the Akita and Yoshida's correlation.Also, Ruzicka et al. (2001) unambiguously demonstrated that an increase in the column diameter advances the transition.Therefore, one expects some dependency of the void fraction on the bubble column diameter.According to the fits proposed here, ε axis evolves as ∼V 0.4 sg and as ∼D −0.2 .
(5.3) Equation ( 19) leads to V L ∼ 1.09V 0.19 sg D 0.3 , which is close to the empirical fit proposed by Raimundo et al. (2019), who wrote V L ∼ 1.35V 0.16 sg D 0.4 .Hence, we recover a formula similar to that empirical fit by exploiting the inertia-buoyancy argument leading to the (gDε) 1/2 velocity scaling, combined with an empirical relationship between the local void fraction and the Froude number.These findings are therefore consistent with each other.
Let us compare (5.3) with available results from the literature.The majority of the correlations proposed for the mean liquid velocity are not in a dimensionless form.Among those that are dimensionally consistent, only one proposal uses (gD) 1/2 as a velocity scale.The latter writes V L /(gD) 1/2 = 0.737Fr 1/3 .
(5.4) Equation ( 5 the void fraction and the mean liquid velocity are the fits proposed by Forret et al. (2006) and the mean gas velocity is the fit proposed by Lefebvre et al. (2022).For the gas, the profile only concerns upward directed velocities, the downward directed gas velocity between x/R = 0.85 and x/R = 1 is not represented.
same prediction.Equation (5.4) was derived using modelling considerations: it is based on a simplified axial momentum equilibrium and on a global energy balance where the dissipation in the column is estimated from mean liquid velocity profiles using a mixing length approach.According to Kawase & Moo-Young (1986), (5.4) is reliable for column diameters from 0.1 to 1 m, and for Newtonian fluids with a dynamic viscosity between 10 −3 and 2 × 10 −2 Pa s: the range of Fr is not indicated.Besides, Kawase & Moo-Young (1987) correlated the global gas hold-up R G3 for Newtonian fluids as Let us first underline that the exponents of Fr appearing in (5.4) and (5.5) are compatible with the scaling (see (2.2)) derived from inertia-buoyancy equilibrium.However, the 1/3 exponent of Fr in (5.4) is significantly larger than the empirical value 0.2 found here (see (5.3)).Similarly, concerning the dependency of the void fraction on the Froude number, the exponent 2/3 found by Kawase & Moo-Young (1987), is somewhat larger than the empirical value (∼0.38-0.4)found here (see figure 14).Let us finally underline that Kawase & Moo-Young (1987) found (5.5) valid for D between 0.1 and 1.07 m and for Fr from 0.005 to approximately 0.05.It happens that these flow conditions mostly correspond to, or/and are very close to, the homogeneous regime (see Kawase, Halard & Moo-Young 1987), so that the comparison of their proposal with the results we obtained in a pure heterogeneous regime may not be entirely relevant.
Further investigations are therefore required to accurately determine how the void fraction evolves with non-dimensional parameters, and in fine to predict absolute as well as relative velocities in bubble columns operated in the heterogeneous regime.The above proposal is believed to be a rather robust first step in that direction.
At this stage, it is relevant to briefly discuss what would be necessary ingredients to elaborate a predictive model.In Appendix B, we present an elementary model for estimating the void fraction in the heterogeneous regime.Following Zuber & Findlay (1965), we assume a 1-D mean vertical flow in the central part of the column.That central part is considered as an inner tube fed by the liquid flow rate Q Lup set by (3.3) and by a gas flow rate Q Lup = (1 + c)πR 2 V sg .Here, c denotes the fraction of the injected gas flow rate πR 2 V sg that is recirculated.The value of c has been estimated to vary between 0 and 0.2 (Lefebvre et al. 2022).The gas flow rate fraction β injected in that tube is then Q Gup /(Q Gup + Q Lup ).Using Zuber & Findlay's approach, we derive a relationship between the gas concentration ε in that tube and β.The difference between ε and β is controlled by the ratio of the mean gas velocity to the mean liquid velocity, where mean values correspond here to an average over the tube cross-section.Introducing known velocities and void fraction profiles in such a model, the predicted mean void fraction ε happens to be in reasonable agreement (within 20 %) with experiments in the D = 0.4 m column at the beginning of the heterogeneous regime that is up to V sg ∼ 9 cm s −1 .This is a quite encouraging result, and it success lies on the knowledge of the circulated liquid flow rate, and on direct measurements of the actual relative velocity between phases.In addition, as β is only a function of Fr (see (B5) in Appendix B), that elementary model brings some support to the correlation ε(Fr) discussed in this section.
However, the deviation between predictions and measurements happens to monotonically increase with V sg : it reaches 80 % at V sg ∼ 25 cm s −1 in the D = 0.4 m column (see figure 16).A similar trend arises from data extracted from literature for columns diameters between 0.13 and 3 m: a reasonable agreement holds at the beginning of the heterogeneous regime but deviations reach a factor of approximately 2 and above at V sg 40 cm s −1 (see figure 17).
Our understanding of the reason for these deviations is the following.First, we stress that the Zuber & Findlay approach holds for a truly fully developed flow.Second, although self-similarity holds, a few measurements (Lefebvre et al. 2022) indicate that, in the heterogeneous regime, the void fraction on the axis of the column significantly increases with the distance to injection.These measurements indicate also that the axial void fraction gradient increases with V sg .We therefore suspect that such axial evolutions are at the origin of the deviations of the model from reality at large V sg .In other words, the 1-D assumption of the model needs to be relaxed so that an extra dependency of the dynamics on V sg could be identified.It could well be that the downward directed gas flow rate is not fully recirculated at the bottom but that it continuously feeds the central tube.That would render the gas flow rate fraction but also the mean phasic velocities as well as the void fraction dependent on the height.A consequence of that discussion is that much more attention should be paid in future investigations to characterize and to analyse the axial evolutions of key variables.

Conclusions and perspectives
We revisited the hydrodynamics of bubble columns operated in the heterogeneous regime.
Conditions with large enough aspect ratio were selected to ensure the presence of a quasi-fully developed region where transverse profiles of void fraction, liquid and bubble mean velocities remain self-similar.We also focused the analysis on air-water systems in ambient thermodynamic conditions involving bubbles in the wobbling regime with large, O(10 3 ) particle Reynolds numbers.
We have shown that the dynamical equilibrium in these gravity-driven bubbly flows balances liquid inertia with buoyancy.Contrary to thermal convection in tubes that involve an unstable vertical stratification, the vertical density gradient on the axis is stable in bubble columns.Instead, the transport of bubbles induces a transverse gradient in void fraction, so that the driving force of the main motion in bubble columns arises from the radial density distribution.The resulting scaling for velocities is V ∼ (gDε) 1/2 , where D is the diameter of the bubble column, ε the void fraction and g the gravitational acceleration.
Using new experiments performed in a D = 0.4 bubble column, as well as data extracted from the literature, this scaling proposal has been shown to hold for the large-scale motion of both liquid and gas phases.This proposal happens to be valid over a wide range of flow conditions, namely for D from 0.1 to 3 m and for gas superficial velocities V sg up to 60 cm s −1 : the corresponding Froude number Fr = V sg /(gD) 1/2 spans a range from ∼0.02 to ∼0.5.The same scaling applies to velocity fluctuations, as the latter were found proportional to the mean velocity both in the gas and in the liquid.Moreover, we also confirmed the finding of Raimundo et al. (2019) that the recirculating liquid flow rate Q Lup in the heterogeneous regime is uniquely set by the column diameter, namely Q Lup = 0.098S core (gD) 1/2 , where S core is the cross-section of the upward mean-flow region centred on the column axis.That result also supports the fact that (gD) 1/2 is a natural velocity scale for the vertical transport.
Direct measurements of unconditional mean phasic velocities in the D = 0.4 m column show that the relative velocity levels off at the homogeneous-heterogeneous transition, increases with the gas superficial velocity, and seems to asymptote to ∼2.4 times the terminal velocity of bubbles at large V sg .Such an analysis deserves to be pursued in order to connect these findings with the known presence of meso-scale structures in the heterogeneous regime.In particular gas velocity measurements conditioned by the local concentration which are now accessible, will help identifying and quantifying collective effects leading to enhanced relative velocities.
As a prediction of the gas hold-up is still crucially needed, an empirical proposal has been made for air-water systems involving wobbling bubbles, in which the local void fraction ε is a function of the Froude number alone.We also attempted a void fraction prediction using a Zuber & Findlay approach: the latter proves successful just after the homogeneous-heterogeneous transition but it fails at larger gas superficial velocities.That feature indicates that the evolutions of key variables such as void fraction, phasic velocities along a vertical are significant even in the quasi-fully developed region of bubble columns.Further efforts should therefore be dedicated to characterizing the axial evolutions of variables, and to understanding the origin of the global self-organization prevailing in the heterogeneous regime.
Finally, a question of importance concerns the impact of coalescence on the above findings.We anticipate that the flow dynamics discussed here and the proposed scalings would remain valid as far as bubbles do not become too large.This statement is already supported by the experiments we analysed as the latter cover different situations in terms of coalescence efficiency.Another argument is the following: the size of a bubble needed for its terminal velocity to equal the asymptotic relative velocity measured in the D = 0.4 m column (approximately 0.7 m s −1 ) would be 5 cm.Hence, up to that size limit for bubbles, we do not expect the flow dynamics to drastically change because of coalescence.
We have seen that the response of the system does not depend on the static liquid height Ho when the aspect ratio Ho/D is large enough.We are thus left with 5 independent non-dimensional parameters, namely Ar, Fr, Eo, M and std(d)/d.
Within that list, the Eötvös and the Morton numbers control the dynamics of an isolated bubble in a quiescent fluid (Clift, Grace & Weber 2005): in particular, they control the shape of the bubble and its terminal velocity U T in the selected fluids and gravity field.Almost all experimental data mentioned in tables 1, 2, 3 and 4 concern large (say 3 to 10 mm) air bubbles in water for which the particulate Reynolds number is quite high (∼800-2100): they thus all correspond to bubbles in the same regime.The polydispersity is also significant in the experiments in tables tables 1, 2, 3 and 4: the parameter std(d)/d is not often quantified, but available measured size distributions indicate that this parameter keeps the same magnitude even though coalescence efficiency varies.In particular, let us underline that the flow conditions in tables 1, 2, 3 and 4 never concern bubbles whose size becomes of the order of the bubble column diameter (in other words, the flow conditions never correspond to slug flow).
Thus, over the conditions mentioned in tables 1, 2, 3 and 4 that almost exclusively concern high aspect ratio bubble columns operated with a few millimetres in size air bubbles in water under ambient T, P conditions, only two non-dimensional parameters have been significantly varied, namely: Note that the Archimedes number equals Re 2 , where Re is the Reynolds number based on the velocity scale (gD) 1/2 , on the size D of the column and on the viscosity ν L of the liquid.For the data shown figure 8, Ar ranges from 9.8 × 10 9 to 2.6 × 10 14 when in the heterogeneous regime.The Archimedes number Ar is the equivalent of the Grashof number used in thermal convection where changes in density arise from differences in temperature and from fluid dilation instead of differences in local void fraction.In free thermal convection, the transition from laminar to turbulent regime corresponds to a Grashof number of approximately 10 9 (Metais & Eckert 1964).The magnitude of the Archimedes number in the heterogeneous regime discussed above exceeds indeed that critical Grashof limit: the heterogeneous regime in a bubble column can be seen as equivalent to the turbulent regime in thermal free convection.
Owing to that observation, one could be tempted to associate the homogeneous/ heterogeneous transition with a critical Grashof or Archimedes number.However, the above analysis was achieved in the asymptotic limit of a large aspect ratio.As far as the transition is concerned, both the column height (more precisely the static liquid height) and its diameter do affect the transition, as demonstrated by Ruzicka et al. (2001).Consequently, the parameter Ho should be accounted for when discussing the transition, and the Grashof/Archimedes numbers definition should be adapted accordingly.Under such conditions, and as shown by Ruzicka & Thomas (2003), the homogeneous/heterogeneous transition can be seen as an equivalent of thermal layers instability those transition is driven by a Rayleigh number.Experiments indicate that this region has a radius between 0.7R (Kawase & Moo-Young 1986) and 0.71R (Forret et al. 2006).Raimundo et al. (2019) suggested that the actual limit may be the radius that equalizes the cross-section area of upflow and downflow regions i.e. (2 1/2 /2)R.As we consider distances to injection such that the mean flow is steady and quasi-fully developed (see the discussion in § 2), let us assume that there is no radial exchange with the downward directed flow near walls so that the mean flow is purely one-dimensional.
As the core zone forms a vertical inner tube of constant cross-section where the mean flow is one-dimensional, one can rely on kinematic models (Zuber & Findlay 1965) to relate the void fraction to the gas flow rate fraction β.The latter is defined as β = Q Gup /(Q Gup + Q Lup ) where Q Lup (respectively Q Gup ) is the liquid (respectively gas) flow rate flowing upward in the tube.By definition, Q Gup = εU G S core and Q Lup = (1 − ε)U L S core , where S core is the cross-section of the inner tube, and where the brackets .denote the spatial average over the cross-section (i.e.f equals the spatial integral of f over the cross-section divided by the area of that cross-section).Following Zuber & Findlay (1965) (B2) This equation states that the difference between the gas concentration ε and the gas flow rate fraction β is set by the ratio of the relative velocity of the gas averaged over the cross-section of the inner tube U R , to the liquid velocity U L averaged over the same cross-section.As expected, a positive relative velocity leads to a concentration lower than the gas flow rate fraction.
From the known empirical profiles for void fraction and for mean liquid and gas velocities in the heterogeneous regime (see 15), one gets a = 1.02 and b = 0.98.These figures change by less than 1 % when the limit of the inner tube is varied from 0.7R to 0.71R.As both coefficients a and b are very close to unity, and as they do not change with flow conditions, it is the ratio U G / U L = 1 + U R / U L that directly controls the proportionality between ε and β .The ratio U L /U L (0) and U G /U G (0) evaluated from the known velocity profiles are given in table 5.
To exploit (B2) or (B3), one needs to know the gas flow rate fraction β relative to the inner tube, namely β = Q Gup /(Q Gup + Q Lup ) where both liquid Q Lup and gas Q Gup flow rates are evaluated within the inner tube.
(i) For the liquid phase, the total upward liquid flow rate Q Lup is given by (3.3), a result valid over a large range of flow conditions.(ii) For the gas phase, we argue that the gas flow rate flowing upward in the bubble column is equal to the flow rate injected over the entire column cross-section, that is Q G = πR 2 V sg , plus some gas flow rate due to the global recirculation that re-entrain bubbles from the top to the bottom.We assume that the latter is a fraction c of the injected gas flow rate, so that the total upward directed gas flow rate is (1 + c)Q G .(iii) Besides, the mean bubble velocity and void fraction profiles (figure 15) indicate that the gas flows upwards between the axis and a radius equal to ∼0.85R.That region is larger than the inner tube into which the liquid flows upward, and only the gas flux within the inner tube of radius 0.7R to 0.71R must be counted when exploiting (B4).
According to known transverse profiles, the fraction of the gas flow rate flowing upward in the corona between 0.7R or 0.71R and 0.85R represents between 8.2 % and 9.2 % of the total upward gas flow rate.Hence, the actual gas flow rate Q Gup flowing through the inner tube is 91 % ± 0.5 % of the total gas flow rate flowing upward.Therefore where the coefficient 1.83 is known within ±0.02 and S core denotes the cross-section of a tube of radius (2 1/2 /2)R.Therefore, for a bubble column operated the heterogeneous regime (with the mentioned restrictions on flow conditions), the actual gas flow rate fraction in the inner tube writes Note that a direct quantification of c is not accessible because reliable bubble velocity measurements in near wall regions are lacking.Yet, according to the estimations made in Lefebvre et al. (2022), c should vary between 0 and 0.2.This result is expected to hold for all the flow conditions considered here (see § § 3 and 4).
Let us first reanalyse the data collected in the D = 0.4m column and in the heterogeneous regime.From the void fraction ε axis measured on the axis and for H/D = 3.625 (figure 3a), the average void fraction ε in the core region is estimated as 0.873ε axis according to the void fraction profile proposed by Forret et al. (2006).The ratio U G (0)/U L (0) deduced from (B3) and (B5) is plotted vs Fr in figure 16.Three values of c, namely 0; 0.1 and 0.2 have been considered: the results remain identical within 20 %.
From the scalings given by (3.1) and (3.2), one expects U G (0)/U L (0) ∼ 1.65.This value is indeed recovered in figure 16 at the lowest Froude numbers shown, that is for superficial velocities at the beginning of the heterogeneous regime.A second observation is that U G (0)/U L (0) is monotonically increasing with V sg : it growths from 1.5 up to approximately 3 at the largest Fr considered that is Fr ∼ 0.12.This result contrasts with the direct measurements of the V L /(gDε) 1/2 and of V G /(gDε) 1/2 presented in figure 6 that show that the ratio V G /V L does not evolve with V sg .A plausible reason for this apparent discrepancy is that U G (0)/U L (0) is deduced from a 1-D model that represents a strong idealization of the actual flow.In particular, the Zuber & Findlay (1965) approach exploited here requires the flow to be fully developed.However, we have shown (Lefebvre et al. 2022) that in the quasi-fully developed region of a bubble column, the void fraction linearly increases with the height above the injector, and that this increase is larger than the hydrostatic contribution (the latter changes the volume of bubbles and thus the void fraction).Even more, the slope dε axis /dH is zero in the homogeneous regime, while this slope linearly increases with V sg in the heterogeneous regime (see figure 26 in Lefebvre et al. 2022).These observations indicate that the flow is evolving along a vertical, that is, even if self-similarity holds, the values of phasic velocities on the column axis are expected to change with height.In addition, and as shown by the vertical profiles of void fraction, these axial evolutions are significantly sensitive to the amount of gas injected.
(i) Consequently, an extra dependency on V sg , which is not present in the Zuber & Findlay (1965) approach by construction, is thus expected to intervene in the relationship between β and ε .(ii) Another consequence is that not enough attention has been paid so far to the axial evolutions of key variables, and new data are needed to characterize the axial changes.Such information will also be helpful to more precisely identify the experimental trends (see the discussion on dispersion related with figure 14).
Similar conclusions arise from a re-analyse of data from the literature.Figure 17 provides the ratio U G (0)/U L (0) vs Fr as deduced from heterogeneous regime data extracted from the articles quoted in tables 1, 2, 3 and 4 (that is for D between 0.138 and 3 m, and for V sg from 4 cm s −1 to 60 cm s −1 ).The trend is the same as the one identified for the D = 0.4 m column in figure 16 .Let us also notice that the dispersion resulting from the selection of c in the interval [0; 0.2] is less than 20 %: this magnitude is comparable to the dispersion of void measurements observed when considering various heights above injection (see figure 14).The precise value of c is therefore not critical.
The increase of the ratio U G (0)/U L (0) vs Fr observed both in 16 and in figure 17 means that the pure 1-D assumption needs to be relaxed.In particular, it is probable that the downward bubble flux along the walls is not entrained down to the bottom of the column in its entirety.Instead, some bubbles moving downward possibly feed the upflow motion all along the column height.That would continuously feed the upflow region at all altitudes with an extra bubble flux.If so, the gas flow rate fraction becomes a function of the altitude, and so does the void fraction ε axis as well as the transport velocities V L and V G on the axis.New experimental data are needed to test and characterize these options.

Figure 2 .
Figure 2. Velocity p.d.f.s for the liquid (blue lines) and for the bubbles (red lines) measured on the column axis at H/D = 3.625 for v sg = 13 cm s −1 (a), v sg = 16.25 cm s −1 (b) and v sg = 22.75 cm s −1 (c).

952Figure 3 Figure 4 .
Figure 3. (a) Evolution of the local void fraction ε axis on the column axis with the gas superficial velocity V sg .(b) Plot of the apparent rise velocity estimated as V sg /ε axis vs V sg .Measurements with a downward directed Doppler probe at a height H/D = 3.625 above injection.

Figure 5 .
Figure 5. Vertical velocity fluctuations V /V of liquid and gas phases vs the gas superficial velocity V sg .Measurements performed in a D = 0.4 m column, on the column axis at H/D = 3.625.

Figure 6 .
Figure 6.Evolution of phasic velocities ('up and down flow' velocity statistics) scaled by (gDε) 1/2 vs the superficial velocity V sg .The mean (V) and fluctuating (V ) components of the bubble and the liquid vertical velocities as well as the void fraction ε are local quantities measured in a D = 0.4 m column, on the column axis at H/D = 3.625.

952
Figure 7. (a) Upward directed liquid flux Q Lup /[D 2 (gD) 1/2 ] vs V sg measured at different heights above injection in D = 0.4 m columns.(b) Evolution of Q Lup with the bubble column diameter from Raimundo et al. (2019), from Guan et al. (2015) and from present data.

Figure 8 .
Figure8.Evolution of V L /(gDε) 1/2 where V L and ε are measured on the column axis vs the superficial gas velocity from the contributions quoted in tables 1 and 2. The inset plots V L vs (gDε) 1/2 for all data collected on the column axis in the heterogeneous regime for V sg 8 cm s −1 and for 0.1 m D 3 m: the dashed line in the inset corresponds to the fit V L = 0.577(gDε) 1/2 .

952
Figure 9. Evolution of the local void fraction on the column axis vs the superficial gas velocity for all the contributions quoted in tables 1 and 2 and exploited in figure 8.For Forret et al. (2006), the local void fraction has been estimated as the global void fraction divided by 1.5.

952
Figure 11.Evolution of the local void fraction on the column axis vs the superficial gas velocity for contributions quoted in tables 3 and 4 and exploited in figure 10.

952
Figure 12.Evolution of the relative fluctuation in the liquid phase velocity V L /V L measured on the column axis vs the superficial gas velocity.

Figure 14 .
Figure 14.Correlation between the local void fraction on the axis of the column and the Froude number Fr = V sg /(gD) 1/2 for all experiments quoted in tables 1 and 3. Measurements performed for 1.3 H/D 12.To illustrate the impact of H/D, the data with closed symbols correspond to measurements in the D = 0.4 m column performed at moderate (H/D = 2.85, triangles) and at large (H/D = 6.37, squares) distances from injection.
Figure 15.(a) Sketch of the quasi-1-D two-phase flow in the core region of a bubble column operated in the heterogeneous regime.(b) Transverses profiles of flow parameters normalized by their value on the axis vs r/R: the void fraction and the mean liquid velocity are the fits proposed by Forret et al. (2006) and the mean gas velocity is the fit proposed byLefebvre et al. (2022).For the gas, the profile only concerns upward directed velocities, the downward directed gas velocity between x/R = 0.85 and x/R = 1 is not represented.

952FrFigure 16 .
Figure 16.Plot of the ratio U G (0)/U L (0) deduced from (B3) and (B5) vs Fr for the data collected in the D = 0.4 m column at H/D = 3.625 and in the heterogeneous regime.

Figure 17 .
Figure 17.Plot of the ratio U G (0)/U L (0) deduced from (B3) and (B5) vs Fr for datasets from the literature mentioned in tables 1, 2, 3 and 4 in the heterogeneous regime (data within the transition region have been discarded and c = 0 assumed).Data in red are those from the D = 0.4 m column.The fit (black solid line) for the datasets for the literature is 6.91Fr 0.47 .The red solid line corresponds to the fit to the present work's data: 8.20Fr 0.50 .

Table 1 .
List of references and flow conditions exploited to extract liquid velocity and local void fraction on the column axis.Further information is provided in table 2.

Table 2 .
For caption see the next page.

Table 2 (
cntd).List of references and flow conditions exploited to extract liquid velocity and local void fraction on the column axis (this table complements the information given in table1). https://doi.org/10.1017/jfm.2022.833

Table 3 .
For caption see next page.

Table 3 (
cntd).List of references and flow conditions exploited to extract the gas velocity and local void fraction measured on the column axis.Further information is provided in table 4. https://doi.org/10.1017/jfm.2022.833

Table 4 (
cntd).List of references and flow conditions exploited to extract the gas velocity and local void fraction measured on the column axis (this table complements the information given in table3).OD, external (outer) diameter. https://doi.org/10.1017/jfm.2022.833

Table 5 .
Coefficients U L /U L (0) and U G /U G (0) in the heterogeneous regime.