3.1. General Bubble Behavior in CFMF
We study bubble formation as function of the transmembrane pressure (
). Bubble formation only occurs if the transmembrane pressure exceeds the activation pressure. This activation pressure equals the capillary pressure of the meniscus in the narrow pore, which is defined as
, where
is the surface tension between the air and the continuous phase,
is the diameter of the pore (opening), and
is the contact angle between the continuous phase and the membrane surface [
12]. At
as low as 50 mbar, a few pores form bubbles at an extremely low formation frequency (too low for quantitative analysis of
and corresponding
). The activation pressure for bubble formation is thus assumed to be approximately 50 mbar. For contact angles >90°, the air phase would flow out of the pores automatically, and the growing bubble would possibly spread out at the (hydrophobic) membrane surface, which would lead to extensive coalescence at the membrane surface [
11,
16], which is not observed in the present study. Thus, the macroscopic contact angle (120° as measured for an unused membrane) is not the actual contact angle in the system. Possible explanations are differences between micro- and macroscopic contact angle, which is known to occur on porous surfaces such as membranes, but more probably protein adsorption renders the membrane surface more hydrophilic (<90°).
We can estimate the effective diameter
of (circular) pores that are activated based on
= 50 mbar and
varying between equilibrium value (~50 mN/m) and that of pure air/water interface (~72 mN/m). For
we assumed a microscopic contact angle of about 80°. The estimated diameter of the activated pores is of the order of 10 μm, which corresponds with the largest observed pores. This means that at low
, only the relatively large surface pores are activated to form bubbles, while the smaller pores need a higher activation pressure. When increasing
above 50 mbar, bubbles are increasingly present in the membrane module. This behavior is ascribed to the increasing number of active pores (
Figure 3A1,2, indicated by the arrows) and the increasing frequency of the pores, with smaller pores having a lower frequency.
Along the system, going from membrane module (where
is estimated) to downstream in the flow cell (where
is derived), there is more substantial increase in bubble size at higher
(compare
Figure 3A from 2 to 4 versus
Figure 3A from 1 to 3). The two bubble diameters (estimated in the module and flow cell) and their corresponding formation frequencies are plotted to quantify the effects of
on bubble formation and stabilization (
Figure 3B,C). The bubble diameter
only slightly increases when
is increased to up to 300 mbar, showing a similar trend as that was reported in [
8]; and it increases relatively more at
mbar. In contrast to the overall slight variation in
, the coalesced bubble diameter
significantly increases with
. The difference indicates the co-existence of bubble formation and bubble coalescence in the 1 wt % whey protein system.
increases much more strongly than the rather constant
for increasing
(
Figure 3B), which is indicative of an increasing extent of bubble coalescence. Last but not least, with increasing
the frequency of coalesced bubbles was increasingly lower than the corresponding
(
Figure 3C). For a given continuous phase flow rate, bubble formation frequency is determined by both the
and the bubble size, which have opposite effects [
17]. As shown in
Figure 3C,
first increases and then decreases with increasing
. The decrease may be explained by the increased size of the bubbles, possibly in combination with more extensive coalescence (at high
), and this may lead to the much lower overall
.
To confirm whether coalescence mainly occurs in the membrane module or continues downstream, we measured
at two distinct positions along the flow path, namely at a distance of 10 cm and further of 100 cm downstream the membrane module (
Figure S2). The obtained
values are similar, which proves that bubbles coalesce mostly in the membrane module during or shortly after their formation, reaching a stable situation within a short term. This also allows us to limit ourselves to a flow cell position of 10 cm downstream for all the remaining experiments.
3.2. Bubble Formation at the Membrane Surface
We first introduce the forces which dictate bubble formation. As was reported for cross-flow membrane emulsification systems, the bubble is subjected to four forces while growing at the membrane surface, which are the shear force imposed by the continuous phase flow rate and viscosity, the buoyancy force, the inertia force and the interfacial tension force [
10,
17]. For bubble formation in the present system, the interfacial tension force is the holding force, and scales as either
or
, depending on the applicability of a force or torque balance (
is the radius of the pore and
is the in-line radius of bubble which grows towards
) [
18]. The shear force, which scales as
, is considered to be the only driving force among the rest of the forces since the other two forces are at least two orders of magnitude smaller. The bubble stays attached to the pore opening and snaps off when, for a certain bubble size, the shear force exceeds the surface tension force [
19,
20]. The force or torque balance leads to
, with
the capillary number of the flow and power-law exponent
whose value is between 0.5 and 1.
In our experiments, the continuous phase flow rate (
), protein concentration (
) and continuous phase viscosity (
) are varied, and bubble formation at the membrane surface is studied. In accordance with the definition of forces, our experimental parameters can influence this force or torque balance: by increasing either the flow rate or the viscosity of the continuous phase, the shear force increases; by increasing the protein concentration, faster protein adsorption can lower the (dynamic) surface tension faster and thus rapidly reduce the holding force (that is if protein adsorption can appreciably take place within the very short timescales for bubble formation). In both cases, bubbles can be detached earlier from the pore openings [
21,
22].
The bubble size is firstly investigated as function of continuous phase flow rate. Please note that bubbles growing at the membrane surface tend to be deformed by the flow of the continuous phase, particularly when a high continuous phase flow rate (
) is applied.
significantly decreases as
increases (see
Figure 4A), which is in line with what was reported for other membrane foaming systems [
8,
19] or membrane emulsification [
23,
24]. At lower
, the bubble stays attached to the pore for a longer time, thus obtaining a larger size, while at higher
,
reduces to about 50 µm for the highest
’s measured. The results agree with the proposed scaling with
with a fitted power-law exponent
0.8 ± 0.1.
almost shows no dependency on transmembrane pressure for
100–300 mbar for any given
(
Figure 4A). At
400 mbar, the small increase in
represents the volume of air added during the final detachment process of the bubble due to a higher air flow rate [
24]. The strong dependency of
on
indicates that the high shear flow dominates overall behavior; the effects of transmembrane pressure and fluid properties on bubble formation at the membrane surface can be neglected.
Additionally, only a minor decrease can be observed in
as function of protein concentration and continuous phase viscosity (see
Figure 4B,C), which are both evaluated at
200 mbar and
(
Figure 4A). By raising the continuous phase viscosity, the average bubble size
is reduced to about 50 µm due to effects on increasing the shear force and potentially also the (dynamic) surface tension force. The decreasing trend with viscosity collapses onto the same
behavior (assuming a constant surface tension,
Figure S3). For these experiments taken at larger
, the rather constant
can be ascribed to the fast bubble detachment process induced by shear force exerted by the continuous phase [
24]. When the protein concentration increases to 2.5 wt %, the
decreases also just slightly to approximately 51 µm. The decrease is ascribed to faster protein adsorption which further lowers the (dynamic) surface tension and advances the force balance for bubble snap-off. Moreover, because the surface tension only can be decreased to a limited extent, bounded by the equilibrium surface tension of the air/water interface stabilized by whey proteins, the influence of protein adsorption on bubble size is moderate. Hence, bubble formation at the membrane surface falls into a regime where bubble snap-off is dominated by shear force, and under the conditions studied here, mainly controlled by the continuous phase flow rate [
10].
The bubble size distribution varies with bubble size. The histogram plot indicates that bubbles with higher monodispersity are formed as
increases, with the main peak in the plot shifting from right to left (
Figure 4D). In addition to the effects of
, the bubble size distribution can further be narrowed by increasing the protein concentration and/or the continuous phase viscosity (
Figure 4E,F). Bubbles with diameters in the range of 40–60 µm account for more than 90% of the size distributions (
Figure 4E,F). The corresponding coefficient of variation (CV) is below 15% when the highest protein concentration and/or viscosity is used (
Figure 4E,F). However, the average
can be decreased only to a limited extent (to approximately 50 µm), with the smallest bubbles having a diameter of approximately 40 µm (corresponding to the left side of the histogram plot), which is independent of the experimental conditions; and the bubble size distribution cannot be infinitely narrowed. This can be ascribed to the characteristics of protein adsorption and the properties of the used membrane—the pore size and the pore size distribution. Specifically, when bubble formation is much faster than protein adsorption, the surface tension keeps constant as that of a pure interface, limited by the efficiency of protein adsorption; and moreover, the smallest pores that have a lower holding force, produce smaller bubbles thus widening the bubble size distribution [
17], and determine the lower boundary of the bubble size. In addition, interactions between bubbles can influence the bubble size distribution by detaching the bubbles early on [
25].
3.3. Bubble Coalescence
To thoroughly understand the CFMF system, the bubble coalescence is also evaluated as a function of the same experimental parameters. Firstly, within the resolution of our experimental results, coalescence of bubbles growing at adjacent active pores is rarely visually observed. Instead, coalescence often happens upon collision amongst flowing bubbles or a collision between a flowing bubble and a forming bubble. To demonstrate the extent of bubble coalescence, we report here the bubble size measured in the flow cell (
). In addition, we compare the volumes (
and
) of single and coalesced bubbles to estimate the number
of coalescence events a bubble has undergone:
. Because
may show some variation across different positions at the membrane surface (
Figures S4-1 and S4-2), due to the pressure drop in the continuous phase and the presence of bubbles (which further influence the local flow rate and viscosity of the continuous phase), the absolute values of the estimates for
are shown in
Figure S5, yet general conclusions can be drawn and will be discussed below.
For a given transmembrane pressure,
decreases as
increases (
Figure 5A) and is always larger than its corresponding
. Yet,
increases as a function of
, which possibly resulted from faster bubble formation (i.e., less protein adsorption) (
Figure S5). At any fixed
,
as well as
increase with
due to increasing bubble crowdedness (
Figure S5), and higher propensity to coalescence. Furthermore, as function of the protein concentration,
decreases strongly, much more strongly than
(see
Figure 5B). With the highest protein concentration tested (2.5 wt %),
is very similar to
, and the coalescence is effectively suppressed with
(
Figures S5 and S6). When a sufficiently high protein concentration is used, more proteins are able to adsorb to the bubble interface before bubble–bubble interactions take place, thus preventing the bubble coalescence [
26,
27]. Lastly,
also decreases and converges to
when a higher continuous phase viscosity is used (see
Figure 5C), which can be explained by the slower movement of bubbles [
5] and the slower drainage process of the liquid thin film (between bubbles) [
28], delaying (or diminishing) bubble coalescence. To summarize, the continuous phase properties (namely, protein concentration and phase viscosity) show vastly different influence on either bubble formation or coalescence: (1) they can control the size at bubble formation (
) only to a very limited extent, with
the dominating factor; (2) they significantly suppress coalescence and thus the formation of larger bubbles (
).
The coefficient of variation (CV) is used to characterize the bubble size distribution (see
Figure 5D–F). The CV is always higher for the coalesced bubbles than for the initially formed bubbles. The CV can only be reduced to a limited extent when we increase the continuous phase flow rate, and we obtain a minimum CV of 20% and 34% for
and
, respectively. The CV does decrease strongly down to 14% and 17%, respectively, when measures are taken to enhance stabilization of the freshly-created interfaces, namely by raising the protein concentration, or increasing the continuous phase viscosity. In both cases a smaller average bubble size is accompanied by a smaller CV as the bubble size distribution is narrowed on the upper end. These results indicate that to tightly control the properties of bubbles in the end product, irrespective of membrane properties, it is crucial to manipulate the operation conditions and the fluid properties to control bubble formation and, especially, to prevent bubble coalescence by sufficiently fast stabilization of the freshly-created interface and by slowing down the approaching bubbles, and thus the film drainage process.
3.4. Bubble Formation Dynamics—Timescales
Within our experimental resolution, we observe that bubble formation can be separated into two stages: a bubble growing stage and an interval stage. The latter stage corresponds to a period of time between the moment that one bubble detaches and that the next bubble appears. The two corresponding timescales are the bubble growing time (red symbols in
Figure 6A) and the interval time (grey symbols in
Figure 6A). For a given transmembrane pressure, the bubble growing time significantly decreases from 1.64 to 0.13 ms as
increases (
Figure 6A). This is because with increasing
the shear force increases, and the bubble is snapped off faster. The corresponding interval time only slightly decreases and converges towards the bubble growing time for
. We highlight that both timescales can decrease to hundreds of microseconds. Muijlwijk and co-authors (2017) [
29] reported that in a 1 wt %
β-lactoglobulin system, which is the main component of whey protein, bubbles are stable against coalescence only if a 100-millisecond duration (within the experimental resolution) is allowed for protein adsorption before bubble–bubble interactions occur. Additionally, in a 5 wt % whey protein system, we recently demonstrated that micrometer-sized bubbles can be sufficiently stabilized at a timescale larger than 1 millisecond; and at a timescale of ~0.01–1 millisecond, bubble formation co-exists with finite bubble coalescence [
30]. Therefore, the similar observations of bubble formation and controllable bubble coalescence in the current CFMF system can be explained by the intersecting timescales of bubble formation and protein adsorption [
30].
The final bubble size reflects a balance between bubble formation and bubble coalescence. The probability of bubble coalescence is steeply decreasing when a monolayer surface coverage is achieved. When the protein concentration or the continuous phase viscosity is manipulated, and both
and
are fixed, bubbles are likely to grow at a fixed surface expansion rate within a fixed period of time, ranging mostly from 0.13 up to 0.23 ms (see
Figure 6B). To explain bubble stabilization within the 2.5 wt % whey protein system, we introduce a dimensionless
Péclet number (
, defined as
, where
is the diffusion coefficient and
is a characteristic length), which describes the relative importance of convection and diffusion during transport of proteins. We first calculated the diffusion coefficient of
β-lactoglobulin (as the representative component of whey protein) using the Stoke–Einstein equation and then obtained
(
S7). Therefore, the above-mentioned surprisingly high bubble stability can be ascribed to high bulk concentration and bulk convection (enhanced mass transport of proteins). Additionally, given the bubble size encountered in this study, the highly curved surface can also accelerate the protein adsorption process and thus contribute to the high bubble stability [
31].