Rate of gas absorption on a slippery bubble mattress

We investigate the absorption of a pure gas into a liquid in laminar ﬂ ow past a superhydrophobic surface consisting of alternating solid walls and micro-bubbles. We experimentally measure and numerically estimate the dynamic mass transfer of gas absorption at stable gas – liquid interfaces for short contacting times. We study the net rate of gas absorption experimentally by in situ measurements of dissolved oxygen concentration pro ﬁ les in aqueous solutions ﬂ owing over oxygen bubbles by ﬂ uorescence lifetime imaging microscopy. We numerically analyze the dynamics of interfacial mass transfer of dissolved oxygen by considering (i) kinetic equilibrium conditions at bubble surfaces that are conventionally described by Henry's Law and (ii) non-equilibrium conditions at bubble surfaces using Statistical Rate Theory (SRT). Our experimental results show that kinetic equilibrium is not established for short contact times. Mass transfer of gas into liquid ﬂ ow past micro-bubbles can be well described by our simulations performed with the non-equilibrium theory for a short exposure time ( (cid:1) 180 m s) of liquid with a microbubble, deviating from the commonly accepted Henry's Law.


Introduction
2][3] Traditional gas-liquid contacting equipment such as distillation columns, packed towers, tray/plate columns and bubble columns provide direct contact of gas-liquid and phase equilibrium-based absorption/desorption processes at these interfaces.The mass transfer in such processes has been studied extensively and the mass transfer resistances have been reported to be determined by the transport of gaseous and dissolved species in the gas and liquid phases, respectively. 1,4,5nder equilibrium conditions, the concentrations in both gas and liquid phases at the gas-liquid boundary are proportional to their partial pressures, which are commonly described by Henry's Law. 61][12][13][14][15] The overall mass transfer process is reported to consist of four consecutive steps: (i) transport from the bulk gas phase to the outer surface of the membrane, (ii) diffusion through the membrane, (iii) dissolution of gas into liquid, and (iv) liquid phase transport.In these studies, the overall mass transfer coefficients are investigated.Analytical and experimental mass transfer correlations are obtained for developing/fully developed boundary layers using analogous heat transfer solutions of the Graetz-Lévèque type. 14,15The mass transfer resistance in the liquid phase, in relation to operating parameters, and the mass transfer resistance in the membrane, in relation to the wetting phenomena, have been investigated in detail. 12,13,15The mass transfer resistances in the gas phase and at the gas-liquid surfaces are commonly neglected due to the fast gas phase diffusion and the assumption of equilibrium at the gas-liquid boundaries.
7][18][19][20][21][22][23][24][25] The surface resistances have been investigated in many theoretical studies using both phase-equilibrium and non-equilibrium conditions.For example, the stagnant lm theory, 18 Higbie's penetration theory 19 and the lm-penetration theory 20 are among the widely used interface models under equilibrium conditions.The dynamics of non-equilibrium interface transport has also been investigated in previous analytical and theoretical studies based on kinetic and thermodynamic considerations. 22,24Moreover, the mass transfer resistance of gas dissolution at the gas-liquid interfaces has been reported as the rate limiting step in gas absorption/desorption processes for short contacting times. 22,25,26][29][30][31] However, to the best of our knowledge, the mass transfer resistances of gas dissolution into a liquid under pressuredriven laminar ow past slippery hydrophobic substrates with a slip velocity have not been investigated.To generate hydrodynamic slippage, hydrophobic substrates containing gas bubbles have been used in recent studies, [32][33][34] as the liquid ows over those shear-free bubble surfaces.3][34][35][36][37][38][39][40] Such uidic platforms with slippery interfaces have been shown to enhance liquid convection and have been suggested to increase the interfacial mass transport rates. 33,39,41In a recent theoretical study, 42 enhancements of up to two orders of magnitude in interfacially driven transport phenomena are predicted on slippery surfaces.These studies highlight the room for investigation of interfacial gas transport at bubble surfaces.In our previous study, 33 we presented a hydrophobic microuidic device that allows for the formation of stable and activegas pressure controlledmicrobubbles at the boundary of microchannels.We experimentally and numerically examined the effect of micro-bubble protrusion angle q and surface porosity on the slippage on our bubble mattresses.Our results revealed the dependency of hydrodynamic slip on the gas-liquid interface curvature and surface porosity.
In this study, we examine experimentally and numerically the oxygen absorption into a water stream in a pressure-driven ow over a superhydrophobic surface with transversely embedded oxygen gas micro-bubbles (Fig. 1).The fabricated devices allow for tunable gas phase pressure, bubble geometry, and liquid phase ow. 33We measure the detailed spatial dissolved oxygen gas concentration proles in the aqueous phase with frequency-domain uorescence lifetime imaging microscopy (FLIM).We numerically study the dynamics of interfacial mass transfer of dissolved oxygen in the liquid side microchannel embedded with curved oxygen microbubbles, using two models of interfacial gas concentration at the (i) kinetic equilibrium state and (ii) non-equilibrium state using Statistical Rate Theory (SRT).

Microuidic setup
Fig. 1a shows a representative scanning microscopy image of silicon-glass based microuidic devices.The microuidic devices consist of two main parallel microchannels for separate aqueous solution and oxygen gas streams, connected by an array of side channels. 33Silicon microchannels were fabricated by standard photolithography followed by a deep reactive ion etching process.The microchannels were sealed by anodic bonding to glass.To prevent the wetting of the connecting side channels, the original hydrophilic silicon microchannels were hydrophobized.We adopted the protocol described in ref. 43 for the hydrophobization via ultraviolet (UV) light illumination of silicone oil.The width of the main microchannels H was 50 mm or 100 mm.The width of the gas-lled side channels L g was kept constant at 20 mm and the length of the liquid-solid interface L s was 20 mm or 30 mm.The periodic bubble unit length is dened as L ¼ L g + L s and the surface porosity (shear-free fraction, 4 ¼ L g /L) is 0.38 or 0.54 for different L s (Fig. 1).The depth of the microchannels was 100 mm.
The ow rate of the liquid side was controlled using a syringe pump (HARVARD Apparatus PHD 2000) and the pressure of the oxygen gas stream was controlled using a pressure controller (Bronkhorst, EL-PRESS, P-602C).The side-channels were lled with oxygen gas due to the hydrophobicity of the microchannel walls and sufficiently large gas pressure P g .The protrusion angles q of the oxygen micro-bubbles into the aqueous stream were controlled by the active control of the applied oxygen gas pressure.The protrusion angles of the oxygen micro-bubbles were varied in a range from 2 to 45 by changing the applied gas pressure ($0.6 to $1.2 bar) for varying liquid ow rates from 9 ml min À1 to 45 ml min À1 .During the measurements the microbubbles were stable and symmetrically pinned at the sharp corners of the side channels.

FLIM experiments
The local oxygen concentrations dissolved in the liquid side channels were quantied by using a Ru-based oxygen sensitive luminescent dye, ruthenium tris(2,2 0 -dipyridyl) dichloride hexahydrate (RTDP) (Sigma Aldrich).The frequency domain uorescence lifetime imaging microscopy (FD-FLIM) technique was employed to measure the lifetime domains of RTDP aqueous solutions in the microchannels. 44,45With the lifetime measurements, the dissolved oxygen concentrations can be calculated. 44For the FLIM experiments, a LIFA-X system (Lambert Instruments) was used on a Zeiss Axio Observer inverted microscope attached to a modulated LED light source and an intensied CCD camera (LI2CAM) in gated mode.During the measurements, the intensier MCP (micro-channel plate) voltage was kept at 650 V.The optical sensing was carried out at a xed modulation frequency of 100 kHz.We used 10 mM uorescein (Sigma Aldrich) solution as a reference uorophore to test any instrumentation phase shi at the operating frequency.Fluorescein is suitable as a reference phase shi uorophore due to its short lifetime s < 6 ns compared to RTDP ($600 ns).The data acquired were analyzed with the LI-FLIM soware package to resolve the lifetime elds.
400 mM aqueous solutions of RTDP were prepared.Lifetimes of the RTDP were measured in oxygen-free (N 2 saturated), aerated and oxygen-saturated aqueous solutions.The solutions were saturated with a certain gas by bubbling that gas through the solution for sufficiently long times.The dissolved oxygen contents of the solutions were conrmed using a ber optic oxygen sensor (FIBOX 3, PreSens).For these measurements using the microuidic platform, the microbubbles were established by nitrogen, air and oxygen gases for corresponding oxygen-free, aerated and oxygen-saturated conditions.
During the mass transfer experiments, oxygen gas microbubbles were established at the boundary of the microchannels (Fig. 1c and 3).If otherwise not stated, deoxygenated RTDP aqueous solution was the working liquid owing past the transversely aligned oxygen bubbles.The experiments were performed at room temperature without further temperature regulation.

Governing equations
Numerically we study the mass transfer characteristics of oxygen gas dissolution at the surfaces of curved oxygen microbubbles aligned perpendicular to a pressure-driven microuidic laminar ow.We solve for coupled mass and momentum transport on our bubble mattress geometry (Comsol Multiphysics v4.3).The computational domain represents the experimental geometric parameters (Fig. 1b) and we solve for velocity and concentration proles in the liquid side microchannels with the same experimental operating settings.In our simulations, the gas side transport is not considered since the effects of mass transfer limitations in the gas phase and the effects of oxygen depletion during gas dissolution vanish due to the active supply of pure oxygen gas.
The Navier-Stokes equations and the conservation of species equation (convection-diffusion) are solved for a steady pressuredriven ow of water in a microchannel consisting of 60 successive bubble units at the bottom surface (y ¼ 0).In order to eliminate any developing ow effects and entrance/outlet effects, an entrance and outlet section of 2-bubble unit cell length was included in the front and end of the bubble mattress.The bubbles are non-deformed and pinned at the corners of side channels, as observed in the experiments.We approximate the bubble proles as circular arcs.We parametrize the bubble interface curvature by the protrusion angle q (Fig. 1b).
The pressure-driven laminar ow was produced by applying a mean velocity as the inlet condition.We varied the inlet mean velocities in the range of 0.04 to 0.2 m s À1 , consistent with our experiments.We applied shear-free boundary conditions (i.e.perfect slip) for the bubble surfaces and no-slip boundary conditions for the solid microchannel walls (at the bottom wall along each L s at y ¼ 0 and the upper wall at y ¼ H).The convective velocity (u) was obtained from the Navier-Stokes equations.
The steady-state conservation of species eqn (1) was solved for the dissolved oxygen transport.The oxygen dissolves at the microbubble surfaces.
Here D is the diffusion coefficient of oxygen in water (¼ 2.67 Â 10 À9 m 2 s À1 at 25 C (ref.25)).C O 2 is the dissolved oxygen concentration in water.The inlet boundary condition to the microchannel was set at a constant concentration value.If otherwise not stated, the inlet concentration is C O 2 | (x¼0) ¼ 0, representing deoxygenated solutions in the experiments.We applied no-ux boundary conditions (vC O 2 /vy ¼ 0) at the solid microchannel walls (at each L s at y ¼ 0 and the upper wall at y ¼ H), as the silicon microchannel walls are not permeable to oxygen.
The mass transfer boundary conditions at the bubble surfaces represent the gas dissolution mechanism at the gasliquid interface that determines the rate of gas absorption by the liquid.We examine two different boundary conditions along the gas protrusions as described below (Fig. 2).

Mass transfer boundary conditions
7][18][19][20][21][22][23][24][25] The surface resistances are commonly described under phaseequilibrium conditions.Higbie's penetration theory 19 and the two lm-penetration theory 20 are among the well-known mathematical models.It is important to note that these theories consider equilibrium conditions at the interface.The net exchange rate of absorbed and desorbed species at the surface determines the amount of accumulated species at the interface and thereby the interface concentration.Under equilibrium conditions, the amount of accumulated species is given by Henry's Law. 63][24][25] In kinetic interface models, the net exchange rate at the interface is described by the Hertz-Knudsen equation where mass accommodation (condensation) coefficients are used as tting parameters. 22,23esides the kinetic approaches, a statistical rate theory (SRT) is suggested to predict the instantaneous rate of interfacial molecular transport. 24,25In SRT the gas-liquid boundary is described to consist of two interfaces; the gas-surface interface and the surface-liquid interface, with the consideration of the Gibbs dividing surface 26 (Fig. 2b).The molecular distributions in each phase (bulk gas, surface, and bulk liquid) are related to thermodynamic properties using the Boltzmann denition of entropy.The transition probability of molecules from one phase to another is derived via a rst order perturbation analysis of the Schrödinger equation.The calculated transition probability of gaseous species from the bulk gas phase to the bulk liquid phase has been used in approximating the instantaneous nonequilibrium exchange rate of absorbed and desorbed molecules at the surfaces.For long contacting times, the SRT predicts equilibrium exchange rates and equilibrium interface concentrations at the saturation value. 26he theoretical models [16][17][18][19][20][21][22][23][24][25] for interfaces would suggest that on hybrid substrates with hydrodynamic slippage, the mass transfer resistance of gas dissolution at the phase boundary can become the rate limiting step for short contacting times of the liquid with bubbles.Therefore in our numerical simulations, in addition to the phase equilibrium state (Fig. 2a), we also study the dynamics of interfacial mass transfer at curved oxygen microbubbles at the non-equilibrium state (Fig. 2b).
Firstly, the saturation concentration of oxygen in water is applied as boundary conditions at the bubble surfaces in our phase-equilibrium based simulations (Fig. 2a).Under equilibrium conditions, the saturation value of the dissolved oxygen concentration at the interface C i O2 is proportional to its partial pressure P O 2 and described by Henry's Law 6 where K H is the Henry's law constant given as K H ¼ 7.7 Â 10 7 Pa l mol À1 at 25 C (ref. 1) and C sat O2 is the saturation concentration of oxygen in water.For pure oxygen bubbles, the dissolved oxygen concentration at the interface is calculated to be C i O2 ¼ C sat O2 ¼ 1:3 mol m À3 at 1 atm and 25 C. Secondly, the rate of molecular transport between the surface phase and the bulk liquid phase is predicted by SRT and applied as ux boundary conditions along the gas protrusions (coinciding each L g at y ¼ 0) in our non-equilibrium based simulations.According to SRT, this instantaneous molecular rate of transport between the phases is described as: where C sat O2 is the xed equilibrium concentration of the dissolved oxygen, hence the solubility at the given temperature and pressure.C i O2 is the interface concentration of oxygen in the liquid phase.At the bubble surface, concentration is a function of the axial x location and to some extent the bubble protrusion depth into liquid in the y direction (Fig. 2).K sl (moles m À2 s À1 ) is the equilibrium oxygen exchange rate between the bubble surface and the liquid phase and has a constant value for a system with given local thermodynamic properties.
where P g w is the partial pressure of water in the gaseous phase, which is the vapor pressure of the owing water.P g O2 is the partial pressure of oxygen in the gaseous phase and is approximated by neglecting the Laplace pressure as P g O2 ¼ P À P g w where P is the total pressure at the bubble surface.C l w is the concentration of water in the liquid phase, which is 55.2 M under the experimental conditions.Here, r and m are the radius and molecular mass of the oxygen gas molecule, k is the Boltzmann constant, and T is the temperature.K H 0 is the "pressure-based" Henry's law constant, given as 4.85 Â 10 9 Pa at 25 C. 1 Without any tting parameters, with a known value of K sl , we implemented the expression given in eqn (3) as boundary conditions dening the oxygen ux at the bubble surfaces.It is important to note here that for longer contacting times, the interface concentration approaches the saturation value ðC i O2 /C sat O2 Þ, and the rate of transport between phases ð J s O2 Þ vanishes, which is indicative of equilibrium conditions, i.e. a constant concentration, at the bubble surface.

FLIM measurements
Oxygen has been reported elsewhere 46 to be an effective collisional dynamic quencher of RTDP and its effect on RTDP uorescence decay is described by the Stern-Volmer equation ( 5) where I o and s o are the uorescent intensity and the lifetime of unquenched RTDP under oxygen-free conditions, respectively.I and s are the uorescent intensity and the lifetime of RTDP, respectively, for an arbitrary oxygen concentration of the solution.K q is the Stern-Volmer constant.In previous studies, the uorescence intensity and lifetime of RTDP were calibrated as a function of oxygen concentration and linear relationships were reported in agreement with the Stern-Volmer relationship. 2,44,45,47-50K q values were reported in the range of 2.1 to 2.7 Â 10 À3 mM À1 . 44,45,48In most of these studies, the oxygen concentrations were in the physiologically attractive range from oxygen-free to air saturated solutions.The linear response of RTDP, consistent with the Stern-Volmer equation, was also shown in the entire range of oxygen up to the oxygen saturation concentration. 47,51Based on our standard measurements for oxygen-free, aerated and oxygen-saturated aqueous solutions, we calculate the experimental K q to be 1.9 Â 10 À3 mM À1 , consistent with the literature ndings. 44,45,48

Rate of O 2 absorption
During the mass transfer experiments, oxygen gas microbubbles were established at the boundary of the microchannels (Fig. 3).We measured the lifetime of RTDP across the liquid side microchannel height (0 # y # H) at different axial locations x (Fig. 3).The bubble interface proles and locations were experimentally determined by locating the minimum lifetime data measured near the hybrid wall.The bubble proles were also calculated by circular arc estimation and depicted by the dashed lines in Fig. 3. Fig. 3 shows the successive lifetime elds measured with FLIM at different axial locations along the same microchannel embedded with microbubbles q ¼ 43 AE 2 .In Fig. 3c, the thickness of the diffusion boundary layer is $23% of the microchannel height H, and does not extend further into the microchannel due to a relatively large Reynolds number Re ¼ 7.5.
From the lifetime elds as shown in Fig. 3, the dissolved oxygen concentration gradients were obtained by correlating the lifetime data by eqn (5).When processing the data, no ltering or smoothing was applied to the measured lifetime elds.Due to the scattering of the experimental lifetime data, the obtained oxygen concentration values were averaged over the bubble surfaces (over $30 pixels in x, equivalent to L g ).The presented results are the average of independent experiments performed with the same operating uidic settings.The statistical error bars are calculated from these experiments.
Fig. 4a-c represent the local dissolved oxygen concentration proles obtained from the lifetime elds shown in Fig. 3a-c, which are successive measurements at different axial x positions in the microchannel during the same experiment.The local proles across the microchannel height, with respect to the y direction normal to the bubble surfaces, are obtained at x positions coinciding with the bubble surfaces that are indicated with dashed arrows in Fig. 3.Here the deoxygenated RTDP aqueous solution ows at a Reynolds number Re ¼ 7.5.The experimental results are compared with numerical results obtained for phase-equilibrium based and SRT-based gas absorption.
The bubble interface position is y $ 3.8 mm in Fig. 4 due to a large protrusion angle q ¼ 43 AE 2 of the bubbles into the water stream.The interface O 2 concentration predicted by phaseequilibrium based simulations converges to the saturation value 1.3 mol m À3 as imposed by Henry's law.Whereas the interfacial dissolved gas concentration predicted by SRT based simulations is $0.8 mol m À3 at x ¼ 1.05 mm and slightly increases for increasing x.When y T 3.8 mm, the oxygen concentration gradients are steep due to a relatively large Reynolds number.As can be clearly seen in Fig. 4a-c, the effect of the microbubbles on the concentration proles is damped when y T 25 mm at x ¼ 1.05 mm and when y T 35 mm at x ¼ 2.30 mm.Beyond the boundary layers, the oxygen concentration drops to the inlet bulk concentration value.Our ndings are consistent with previous studies 45,51 reporting the local oxygen gradients resolved by in situ FLIM measurements.In these studies, oxygen diffuses from an oxygen enriched liquid phase to an oxygen-free liquid phase.
The FLIM measurements in Fig. 4a-c reveal very low concentrations near bubble interfaces when y ( 10 mm.There is a sharp discrepancy between the measurements and both numerical results in the proximity of the bubble interfaces (3.8 mm # y ( 10 mm).The standard errors calculated close to the bubble surface are also larger compared to those of the remainder of the microchannel.The deviation can be related to the systematic error in our FLIM measurements close to the gas-liquid interfaces (additional measurements are provided in the Appendix).The deviations near the bubble surfaces do not interfere with the results presented here since we do not consider the experimental data obtained in the proximity of bubble interfaces.
Considering the oxygen concentration data when z10 mm # y # H, the measured oxygen concentration proles are in good agreement with the simulations performed with an incorporated SRT boundary condition.Whereas, the numerical results computed for a constant interface concentration at the solubility limit of oxygen in water predict steeper proles and faster gas absorption than those in the measurements.Due to a relatively large Re ¼ 7.5, the water residence time within the hold-up volume of the microchannel embedded with the bubbles is as small as 96 ms.The inset in Fig. 4d shows the numerically calculated ow velocity by solving the Navier-Stokes equations (in Section 3.1) using the same experimental parameters in these sets of experiments, revealing a slip velocity of $25 mm s À1 , consistent with our previous study. 33This slip velocity can induce a very short contact time of the water stream passing a single bubble interface.The exposure time of a slippery bubble surface to water is calculated to be $830 ms.These short exposure times are insufficient for a gas-liquid phase equilibrium.Our ndings are indicative of an additional mass transfer resistance induced by the gas dissolution at the bubble surfaces at short exposure times on a laminar ow over a slippery bubble mattress.
To further highlight our results indicative of interfacial mass transfer resistances, we present the local instantaneous convective ux distribution in Fig. 4d.The concentration prole shown in Fig. 4c and the numerical ow velocity prole (Fig. 4d, inset) are multiplied to calculate the local ux.The local ux near the bubble surfaces has low values up to y z 8 mm due to the low local velocity in this region even though the dissolved oxygen concentrations are at their maximum values.The ux then increases to a peak value when 8 mm ( y ( 11 mm.Above y z 11 mm, the local ux decreases with a slope mainly determined by the decrease in C O 2 concentration towards the end of the diffusion layer at y z 35 mm.The differences in the experimental and numerical slopes are mainly attributed to concentration proles, as the coupling velocity prole is unique at the given settings.Here the differences in the uxes obtained by the measurements and phase-equilibrium based simulations are more pronounced. The O 2 uxes are determined by integration of the local ux proles.In Fig. 5, we present the O 2 uxes obtained from the measurements, the simulations performed under phase equilibrium conditions on the bubble surfaces and the simulations incorporating an interfacial mass transfer resistance using SRT.The experimental and numerical total convective uxes are calculated for the entire microchannel height, i.e. for 3.8 mm # y # H. Whereas the area under the local ux curves when 11 mm ( y ( 35 mm contributes most to the total convective ux values.It is important to note here that in this region the experimental values are not hampered by the systematic error in the measurements (y ( 0.1H).The experimental oxygen uxes revealed in Fig. 5 are in a good quantitative agreement with the ones computed by considering an additional gas dissolution resistance using SRT.The ux values calculated with the assumption of a spontaneous phase-equilibrium state are higher, as in this case the interface dissolved oxygen concentration is at the maximum value.The phase equilibrium ux obtained at x ¼ 1.05 mm for Re ¼ 7.5 is $1.8 times higher than the ones obtained from the measurements and the SRT-based numerics.In Fig. 5, the oxygen uxes obtained at Re ¼ 10 are also presented.The deviation from the phase equilibrium mass transfer is more pronounced for a higher Reynolds number Re ¼ 10.The phase equilibrium ux obtained at x ¼ 1.05 mm is $2.2 times higher than those obtained from the measurements and the SRT-based numerics.An increase in Re from 7.5 to 10 results in a decrease in water residence time within the hold-up volume of bubble mattress from 96 ms to 48 ms, with an increase in mean velocity.Furthermore for a higher Reynolds number Re ¼ 10, the slip velocity at the hybrid wall (y ¼ 3 mm for q ¼ 35 AE 3 ) has a much higher value $112 mm s À1 , thereby a much shorter exposure time $180 ms of a slippery bubble surface to the water stream compared to the exposure time $830 ms for Re ¼ 7.5 (Fig. 5).These ndings show that these short exposure times lead to larger deviations between the equilibrium based boundary conditions (Henry's law) and the non-equilibrium based boundary conditions (SRT).
We compare the experimental and gas dissolution resistance incorporated-numerical O 2 uxes obtained at x ¼ 1.05 mm for Re ¼ 7.5 and Re ¼ 10.The experimental ux value for Re ¼ 10 is $2.3 fold higher than that for Re ¼ 7.5.In good agreement with the measurements, our simulations predict a $3 fold higher ux value for Re ¼ 10 compared to Re ¼ 7.5 (Fig. 5).This slight increase in the Reynolds number results in a larger decrease in exposure times due to a larger slip velocity near the bubble surface, and thereby a decrease in the dissolved oxygen concentration at the interface.Nevertheless, this enhanced hydrodynamics results in enhanced convection, and therefore an enhanced convective mass ux, as tripled in this case.Consistent with the literature ndings, 10,39,42,52 our ndings clearly demonstrate that the hydrodynamic slippage amplies the mass transport on hybrid substrates consisting of liquidsolid and liquid-gas interfaces.

Conclusions
We studied the interfacial mass transport accompanied by hydrodynamic slippage at gas-liquid interfaces in microuidic devices.The uidic conguration allows for the investigation of high-resolution gas absorption dynamics at curved gas-liquid interfaces.The local oxygen concentration elds resolved by FLIM measurements reveal slower gas absorption rates compared to those predicted by phase-equilibrium based simulations.The experimental results are in good agreement with the numerical results obtained from the simulations considering an additional mass transfer resistance at the gas-liquid boundary during gas dissolution.Our results indicate that the equilibrium state may not be established at short contacting times.Our ndings also reveal that even at such short exposure times, the hydrodynamic slippage enhances the total mass ux driven by the amplied ow throughput on the bubble mattresses.

Appendix Additional measurements and near-bubble deviations
We verify our FLIM results by the measurements of different inlet oxygen concentrations which were prescribed to certain values prior to microuidic mass transfer experiments by the aid of an external oxygen sensor during the bulk solution preparation.Fig. 6a and b show the dissolved oxygen concentration proles across the microchannel height, with respect to the y direction normal to the bubble surfaces for varying inlet concentrations of 0.5 and 0.2 mol m À3 .The experimental results are compared with numerical results obtained for phaseequilibrium based and SRT-based gas absorption.
The bubble interface position is y $ 0.3 mm as shown in Fig. 6a, as the bubbles are nearly at (q ¼ 2 AE 1 ).At the bubble interface (y $ 0.3 mm), the oxygen concentration is measured to be $0.9 mol m À3 , consistent with the SRT boundary condition numerical results.Whereas the interface O 2 concentration predicted by phase-equilibrium based simulations converges to a saturation value of 1.3 mol m À3 .The thickness of the diffusion boundary layer obtained by FLIM measurements extends to $18 mm, consistent with the numerical results.Beyond the thickness of the diffusion boundary layer, the dissolved oxygen concentration drops to an inlet bulk concentration of 0.5 mol m À3 .Similarly, the experimental oxygen concentration prole shown in Fig. 6b reaches an inlet bulk concentration of 0.2 mol m À3 beyond the diffusion boundary layer thickness at y $ 14 mm.In Fig. 6b, the measured O 2 prole is in good agreement with SRT based simulations when y T 4 mm.However there is a discrepancy between the measurements and both numerical results in the proximity of the bubble interfaces, as shown in Fig. 4. The standard errors close to the bubble surface are larger compared to those of the remainder of the microchannel.The discrepancy in the measured and calculated values when y ( 4 mm is still larger than the standard error margins.The near-bubble deviations can be related to the systematic error in our FLIM measurements close to the gas-liquid interfaces.In the FLIM measurements, the normalized uorescence lifetime signal decreases from 1 (oxygen-free state) to $0.3 (oxygen-saturated state), providing a signal level of $0.7 in the entire oxygen concentration range.We measure the noise level to be $0.1, hence a signal to noise ratio of $7. However we observed reduced signal-to-noise ratios near microchannel walls and bubble interfaces.Such deviations near gas-liquid interfaces have also been reported in previous reports [53][54][55][56] and are attributed to optical distortions and reection effects near interfaces.Besides, we observed differences in the range of z7-15% between the lifetime data calculated by the phase shi and that of the amplitude demodulation of the emitted signal near interfaces.Differences in the phase shi s p and demodulation s m lifetimes have been reported to be indicative of complex quenching mechanisms and are observed for long lifetime probes as Ru-based complexes in previous studies.

Fig. 1
Fig. 1 Microfluidic bubble mattress and oxygen dissolution at bubble surfaces.(a) Scanning electron microscopy image of the microfluidic device showing two main microchannels for oxygen gas (P O 2 ) and water (u x (y)) streams connected by oxygen-filled side channels.(b) Numerical results of the dissolved oxygen concentration in the water side microchannel solved with identical operating settings as in (c).(c) Lifetime field resolved by FLIM superimposed on the brightfield microscopy image showing bubbles protruding at 43 AE 2 into the water side microchannel with a height of H ¼ 100 mm.L g is the width of the oxygenfilled side channel (L g ¼ 20 mm), and L s is the width of the solid boundary (L s ¼ 32.6 mm).The color bars refer to the lifetime of RTDP aqueous solution which is given in ns (c) and the oxygen concentration which is given in mol m À3 in (b).In b and c, surface porosity 4 ¼ 0.38 and Q w ¼ 45 ml min À1 .

Fig. 2
Fig. 2 Schematic illustration of the applied boundary conditions and the effects of interface resistances on the concentration profiles in the presence of liquid flow.(a) The phase equilibrium state where the interface liquid concentration is fixed at the saturation value given by Henry's Law ( ).The time (t) evolution results in increasing boundary layer thicknesses.(b) The non-equilibrium state where the interface liquid concentrations ( ) are determined by the instantaneous exchange rates of absorbed and desorbed molecules at the surface liquid boundary (sl).The time (t) evolution results in increasing boundary layer thicknesses and increasing interface concentrations approaching the saturation value.In (a and b) the mass transfer resistances in the bulk gas phase are neglected.

Fig. 3
Fig. 3 Successive lifetime fields in axial position x.Quantitative visualization of the increasing boundary layer thickness along downstream flow with Q w ¼ 45 ml min À1 at (a) 0.7 ( x ( 1.1, (b) 1.4 ( x ( 1.8, and (c) 1.9 ( x ( 2.3.Here q ¼ 43 AE 2 and 4 ¼ 0.38.In (a-c) the flow direction is from left to right.The color bar refers to the lifetime, which is given in nanoseconds.The dashed arrows indicate the axial positions at which the local oxygen concentration profiles across the microchannel height are obtained.

Fig. 4
Fig. 4 Successive dissolved oxygen profiles at different axial x positions.(a) Local O 2 profile at x ¼ 1.05 mm.The inset is a zoom-in of the O 2 profile in (a).(b) Local O 2 profile at x ¼ 1.80 mm.The inset is a zoom-in of the O 2 profile in (b).(c) Local O 2 profile at x ¼ 2.30 mm.The inset is a zoom-in of the O 2 profile in (c).(d) Local convective flux J O 2 profile at x ¼ 2.30 mm.The inset shows the flow velocity profile obtained by simulations performed at the same experimental parameters.In (a-d) symbols represent the experimental results obtained by FLIM, the color solid lines represent the numerical results solved with SRT, and the color dashed lines represent the numerical results solved using Henry's law.Here the microchannel height H ¼ 100 mm, q ¼ 43 AE 2 , 4 ¼ 0.38, and Re ¼ 7.5.The black dashed lines in (a-c) depict the equilibrium saturation value of oxygen in water.

Fig. 5
Fig. 5 O 2 fluxes J O 2 as functions of axial x positions and Re obtained by FLIM measurements and numerical calculations.Red filled circles ( ), the red solid line ( ) and the red dashed line ( ) indicate the experimental, the SRT-based numerical and equilibrium-based numerical results, respectively, for Re ¼ 10 and q ¼ 35 AE 3 .Blue filled squares ( ), the blue solid line ( ) and the blue dotted line ( ) indicate the experimental, the SRT-based numerical and equilibriumbased numerical results, respectively, for Re ¼ 7.5 and q ¼ 43 AE 2 .In (a) and (b), 4 ¼ 0.38.

Fig. 6
Fig. 6 Experimental and numerical local dissolved oxygen profiles.(a) When the inlet oxygen concentration is 0.5 mol m À3 .Here Re ¼ 2.2.(b) When the inlet oxygen concentration is 0.2 mol m À3 .Here Re ¼ 4. In (a and b), symbols represent the experimental results obtained by FLIM, the colored solid lines represent the numerical results solved with SRT, and the colored dashed lines represent the numerical results solved using Henry's law.Here H ¼ 50 mm, q ¼ 2 AE 1 and 4 ¼ 0.54.The black dashed lines depict the equilibrium saturation value of oxygen in water.