In-situ observation of collective bubble collapse dynamics in a quasi-two-dimensional foam.

The stability of foams is an important subject not only for fundamental science, but for applications in daily life. The most destructive phenomenon underpinning foam collapse is a collective bubble collapse, yet the mechanism behind this is unclear. In this study, we clarify the dynamics of the collective bubble collapse in a quasi-two-dimensional foam by in-situ observation with a high speed camera. We find two modes for collective bubble collapse: one is the propagation of liquid film breakage via impact with the stream of another broken liquid film. The other is breakage of a distant liquid film due to penetration by a liquid droplet, emitted by impact with the flow of a broken liquid film. As the liquid fraction increases, the velocity of liquid droplets decreases. Instead of penetration, the liquid droplet bounces like a billiard ball or it is absorbed into other films.

Introduction

A liquid-gas foam is a non-equilibrium jammed state system composed of gas bubbles and liquid films. Foam, which has a cell structure of gas bubbles, is widely used in daily life e.g. in food and cosmetic applications (ice-cream, food mousse and whipped cream, shaving and styling foams), pharmaceutical products (topical and rectal foams), bio-technologies (foam in bio-fermenters), chemical technology (froth flotation) amongst many other applications[1-13]. Therefore, it is a matter of importance for a variety of fields to understand and control foam stability.

The foam is a metastable non-equilibrium state and eventually collapses as time elapses[5,12]. The dynamics of foams is governed by three different processes: these are coarsening, drainage, and collapse in order of evolution rate. Coarsening, where the number of bubbles in the foam decreases and the average bubble size grows, occurs due to different pressures between the bubbles. Drainage leads to thinning of the liquid films due to gravity. Drainage and coarsening are well studied and understood from both experimental and theoretical perspectives. Bubble collapse has also been investigated[9,14-16]. A previous study showed that bubble collapse occurred when a lower coalescence limit was reached in liquid volume fraction[9]. A relationship between bubble collapse and T1 rearrangement has also been studied by injecting air into a system[14]. In addition, when bubbles burst in the foam, a crackling sound is emitted. By measuring and analyzing the acoustic emissions from foams, it has been suggested that a collective bubble collapse (CBC) phenomenon occurs during the collapse of foams[5,12]. This is sometimes called “cascade collapse”. A similar effect has also been observed in emulsions which are even more difficult to analyze due to smaller droplet size. Although CBC plays a role in foam stability, the mechanism behind CBC is still unclear; a framework and perspective for understanding CBC is yet to be established.

In this article, we investigate the mechanism behind CBC. Having prepared a quasi-two dimensional foam, a single bubble is popped at the edge; details are given in the Materials and Methods section. We then directly observe the CBC phenomenon by using a high speed camera. We identify two mechanisms behind the propagation of a rupture event; “propagating” and “penetrating” modes. As increases, it becomes harder for liquid droplets to penetrate liquid films, leading to droplets bouncing off films due to their elasticity, or being absorbed into the film.

Results

Firstly, we perform in-situ observation of the CBC process for different . Figure 1(a) shows enlarged successive images of the CBC process from time t = 0 ms to 3.12 ms for  = 0.0099 [see Supplementary Movie 1]. We puncture a bubble out side the foam using a capillary glass needle with a very small amount of silicon grease on the tip. Collapse of the bubbles is firstly observed at an outer edges since rapture of the foam film proceeds on the outer edge. It can be called a surface effect. Then we find two modes for CBC collapse in bulk: one is a ‘propagating’ mode, the other is a ‘penetrating’ mode. When a liquid film breaks, the broken liquid film is rapidly absorbed by the Plateau border. Due to the strong absorbing impact, another liquid film sharing the same Plateau border is also broken, indicated by the red circle in Fig. 1(a). This is the propagating mode of liquid film breakage [see Fig. 1(b)]. At the same time, when the broken liquid film is absorbed by the Plateau border, a liquid droplet is emitted, indicated by the blue and green circles in Fig. 1(a). The droplets run into the inside of the foam as shown by the arrows in Fig. 1(a), and penetrate a distant liquid film. The velocity of the liquid droplets V is found to be approximately 3 m/s. This penetration impact induces breaking the distant liquid film. This is the penetrating mode [see Fig. 1(b)]. The broken liquid film induces another broken film by both propagating and penetrating modes. The CBC collapse occurs as a result of the repeated breakage of liquid films via the two modes.

Figure 1

(a) Enlarged successive images of the CBC process for  = 0.0099 until t = 3.12 ms. The white bar corresponds to 5 mm. Red dashed circles show the breaking of the liquid film by the propagating mode. Blue and green dashed circles show liquid droplets emitted by the impact of the stream of the broken liquid. The arrows show trajectories of those liquid droplets. (b) Schematic image for CBC. The impact of the stream of the broken liquid induces the breaking of the liquid film (propagating mode) and the emission of a liquid droplet (penetrating mode).

As increases ( ≥ 0.015), the probability of a liquid droplet being emitted decreases. The velocity of liquid droplets also decreases; it becomes harder for liquid droplets to penetrate distant liquid films. Instead of penetration, the liquid droplet bounces like a billiard ball [see Supplementary Movie 2]. Figure 2 shows how the liquid droplet bounces over 30 ms; the dashed line shown in the top indicates the trajectory it follows. Measuring V after each bounce, we can plot V as a function of the number of impacts n, as shown in Fig. 3. We find that V decreases with each impact. We define the coefficient of restitution of the liquid film as where is a droplet velocity after the i-th bounce and we obtain that e = 0.50~0.74. After the droplet bounces several times, it is absorbed into the liquid film, as indicated by the filled symbol in Fig. 3. We note that a bouncing and absorption of the drop is similar to results in the soap film[17].

Figure 2

Enlarged successive images of the bouncing process for  = 0.015 until t = 25.9 ms. The white bar corresponds to 5 mm. Red circles show the liquid droplet emitted by the propagating mode. Yellow dashed line shows the trajectory of the liquid droplet.

Figure 3

The velocity of the liquid droplet V as a function of a number of collisions with the liquid film i. V decreases after each impact. The coefficient of restitution of the liquid film is 0.50~0.74. The droplet is finally absorbed into the liquid film, indicated by the filled symbol.

As is increased further ( > 0.022), the liquid film is successfully absorbed into the Plateau border, and liquid droplets are not emitted. The bubbles at the outer edges of the foam are broken by the surface effect even with greater ; however, fewer bubbles are broken, and the collapse stops quickly. As a result, a CBC does not occur. This suggests that the penetrating mode is crucial for CBC occurrence.

Next, we investigate the dependence of the number of collapsing bubbles quantitatively. A number of total collapsing bubbles N includes the bubble at the outer edge of the foam, which are collapsed by the surface effect, not only the propagation mode and the penetration mode in bulk. Since we focus on the CBC in bulk, we count a number of the collapsing bubbles inside the foam N, which excludes the surface effect. Thus, N corresponding to the number of bubbles broken by both propagating mode and penetrating mode in bulk. We count the number of broken bubbles, from the first to last bubble broken; this occurs over approximately 0.04 seconds. This number includes double collapsing events[18]; the time scale of the CBC in our experiment is consistent with that reported in ref.[18]. Figure 4 shows N (red) and N (blue) as a function of . Triangles, circles, and squares correspond to N or N at N ~ 200 for glycerol concentrations of 9.4 wt%, 17.8 wt% and 29 wt%, respectively. N is the number of bubbles in the whole foam. As shown, both N and N decreases as increases. From the power law fit, we find that with γ = 2.3 ± 0.36. Note that N is quite small for  > 0.015; CBC collapse seems to occur in dry states, where the boundary between dry and wet foams is in a quasi-two-dimensional foam[13,19]. We also find that both N and N is independent of the glycerol concentration below 29 wt%, that is, with a change in viscosity. This suggests that the behavior of the liquid film is primarily inertial rather than viscous. It becomes difficult to puncture a bubble when the glycerol concentration is more than 40 wt% and CBC is not observed. We also note that N (or N) collapses into the power law curve although the size distribution of the bubbles varies in each experiment. We also investigated the number of collapsing bubbles with a larger foam (N ~ 500), as shown by the diamond symbols in Fig. 4. This suggests that N and N are independent of N. We also considered the fact that we used silicone grease in order to induce the first bubble breakage; the grease may also affect the collapse of the other bubbles in the CBC process[20-22]. Thus, we measured N and N without using silicone grease, that is, N and N when the CBC phenomenon spontaneously occurs. This is shown by the filled symbols in Fig. 4. We find that the number of collapsing bubbles without using silicone grease collapses onto the existing curve, suggesting that the CBC phenomenon is independent of the use of the silicone grease.

Figure 4

dependence of N (red symbols) and N (blue symbols). N and N decreases with ; fit lines reveal that γ is 2.3 ± 0.36. Triangles, circles, and squares correspond to N and N for glycerol concentrations of 9.4 wt%, 17.8 wt% and 29 wt%, respectively. Thin diamond symbols correspond to N and N for a circular foam, with a glycerol concentration of 17.8 wt%. Diamond symbols correspond to N and N for a large foam. The filled symbols correspond to N and N when the CBC phenomenon occurs spontaneously before a bubble is popped with a needle. This means that the CBC phenomenon is independent of the use of silicone grease.

When is small, it is likely that the shape of each bubble is anisotropic and that the distorted bubbles reconnected in a chain-like form[13]. Bubbles with larger size than a mean size or anisotropic shape have large excess surface energy and one might think that such larger or anisotropic bubbles tend to collapse easily. Thus, we investigate the relationship between CBC and the size distribution and between CBC and the shapes of the bubbles. We use a parameter λ to characterize the anisotropy of bubble i[13]. λ is defined aswhere j is a pixel at the edge of the bubble, n is a total number of the pixel j, r is a distance between the center of bubble i and the pixel j and is a mean distance of r. λ = 0 when the bubble i is circular, whereas λ > 0 if the bubble is anisotropic. Figure 5(a) is a foam with  = 0.0086 before the CBC. The bubble color is set to black when λ is larger, and closer to white when λ is smaller. Red dots indicate the bubbles collapsed during the CBC event. We find that the bubbles located at the left of the foam are uniformly collapsed. Figure 5(b) shows a probability distribution as a function of a mean bubble diameter of a bubble i d before and after the CBC for the bubbles. d is computed by averaging a distance between a center and an interface. The shape of the probability distribution after the CBC is the same as that before the CBC. Figure 5(c) shows that d after CBC d as a function of d before CBC d. It is found that d = d as shown a line in Fig. 5(c) and it means that each bubble size is unchanged during CBC. We also shows the probability of λ before and after the CBCs in Fig. 5(d). As shown in Fig. 5(d), the shape of the probability distribution after the CBC is also the same as that before the CBC. We compute λ before CBC λ and after CBC λ and we find that for almost all bubbles shown as in Fig. 5(e). It means that the shape of the bubble is not relaxed after the collapsing occurs and it is consistent with a previous report[13]. This suggests that the bubbles collapsed during the CBC are also independent of the shape of the bubbles. In addition, we investigate the dependence of N when the shape of the whole foam is circular. These are shown by the thin diamond symbols in Fig. 4. According to a previous report, the shape of each bubble corresponds to the shape of the whole foam[13]. Thus, this result supports the hypothesis that the CBC phenomenon is independent of the size distribution and the shape of each bubble. It is also reported in ref.[18] that the collapsing bubbles along the air/foam interfaces have no particular topological selection i.e. the bubbles popped are independent of bubble shape. This is consistent with the propagating and penetrating modes being induced by the mechanical impact of the broken liquid film rather than excess surface energy due to the larger and anisotropic shape. Here we should note that the mean bubble size should affect the CBC process. We will discuss it later.

Figure 5

(a) Relationship between bubbles broken during CBC and deformed bubbles in the foam, at  = 0.0086. The bubble color is darker for larger deformations, lighter for smaller deformations. Red dots indicate bubbles broken during the CBC. (b) Probability of a diameter d before (blue line) and after CBC (red line). (c) A relation between bubble diameters before CBC d and after CBC d. It means that CBC is independent of the size distribution. (d) Probability of λ before CBC λ and after CBC λ. (e) A relation between λ before and after CBC. It is found that the shape is not related with the CBC process.

In order to evaluate the mechanical impact of the broken liquid film, we measure the velocity of the broken liquid film V in the foam with 17.8 wt% glycerol when it is absorbed. We define V as l/t, where l is the length of the liquid film and t is the time it takes for the liquid film to be absorbed, from start to finish. Figure 6(a) shows the dependence of V as a log-log plot. We find that V is approximately 1~10 m/s. The order of V is consistent with that in inertial liquid films; note that the velocity in viscous liquid films is ~mm/s[23,24]. From the fit, we obtain V = V0 × −, where V0 = 0.083 ± 0.018 m/s and α = 0.76 ± 0.07. We also obtain V in a foam with 29 wt% glycerol, as shown by the square symbols in Fig. 6. We found that V with 29 wt% glycerol is almost the same as V with 17.8 wt% glycerol. This result suggests that the behavior of the liquid film is inertial rather than viscous. The rapid decrease of V as a function of is consistent with our observations as shown in Fig. 1. With increasing , the emitted droplet is bounced or absorbed due to decreasing of V. Here we compare the Taylor-Culick velocity V. V can be described as where σ and are the surface tension and the density of the liquid, respectively, and h is the thickness of the film[25,26]. α should be 0.25 if V is the Taylor-Culick velocity[25,26] since a local liquid film thickness h is proportional to [27]. It is not consistent with our results (α = 0.76). We also investigate the relation between the velocity and osmotic pressure . The relationship between and has been studied previously[28-30]. According to ref.[28], for 2 dimensional foam becomeswhere σ is a surface tension and R is a mean radius of the bubble and is a jamming point and  = 0.16 in 2 dimension. In our experiment, we used σ = 37 mN/m and R = 1.7 mm. In the Hele-Shaw cell, a thin liquid layer exists on the wall. Thus we should subtract a fraction of the wetting layer from the liquid fraction in order to use the equation. Then we estimate with assuming that a thickness of the thin wetting layer is 1 μm. Figure 6(b) shows V as a function of ; it is found that V is proportional to . We confirm the linear relation between V and even though we estimate with changing the thickness of the wetting layer from 0 μm to 2 μm. From the linear relation, we consider that the driving force of the absorption is negative pressure in the foam film. Furthermore, we also compute dependence of as shown in Fig. S1 in the Supplementary. We find that for 0.005 <  < 0.05. Then we obtain and it is consistent with our observation.

Figure 6

(a) Stream velocity V as a function of , shown as a log-log plot. We find that V is proportional to . Circles and squares correspond to V for glycerol concentrations of 17.8 wt% and 29 wt%, respectively. (b) V as a function of osmotic pressure . is calculated by Eq. (2). It is found that V is proportional to .

Furthermore, we also look at the V dependence of N in Fig. 7. We find that N drastically increases with increasing V. This means that CBC is strongly correlated with the stream velocity of the liquid film. Again, it is consistent with the notion that CBC is induced by the mechanical impact of the liquid film. We also obtain N = N0 × (V/V1) where we set V1 = 1 m/s. The variable parameters are found from the fit, N0 = 1.54 ± 0.56 and β = 2.9 ± 0.40. N steeply increases with increasing V. It supports that the penetration mode is crucial for the CBC. From the relation , we compute the power law relationship where γ = 2.2 ± 0.50. The calculated γ is consistent with γ (2.3); this gives us more assurance as to the reliability of these experimental results.

Figure 7

N dependence of V, shown as a log-log plot. We find N = 1.54 (V/V1)2.9 where we set V1 = 1 m/s.

Discussion

Here, we discuss the linear relationship between V and . We start from the Navier-Stokes equation for the liquid film. We set x axis as a direction of absorption and we neglect velocity in a perpendicular direction to absorption.where P is a pressure and , ν is a kinematic viscosity, and k is a coefficient. −kV is a term for resistance force between a glass substrate and the liquid film. From our observations, V seems almost constant with time. We cannot observe a rim during an absorption, which differs from an absorption of a soap film (Taylor-Culick case). It suggests that local volume is unchanged and then V is almost homogeneous in x direction in the liquid film. Then we also obtain that and it is consistent with our observation, which V is independent of ν. Then we obtain that

Then , where l is a length of the liquid film. Thus, we find that

Our observation of can be explained by a balance between the negative pressure and the resistance between the liquid film and the glass substrate. It is similar that a velocity of a falling liquid droplet is determined by a balance between a gravity force and air resistance. Thus we consider that V is determined by the osmotic pressure, rather than . We note that is related with the bubble size R according to Eq. (2). Thus V depends on R and then the number of the collapsing bubbles should be depend on R. In near future, we investigate the CBC with changing the mean bubble size, not only the bubble size distribution.

Here we compare other experiments, which the drop penetrates the film without breaking. We introduce the weber number We = , where is density. According to previous reports, the drop penetrates the film without breaking when We > 20[31,32]. We calculate We in the penetration mode of the foam using  = 1000 kg/m3, R = 100 μm, V = 3 m/s, and σ = 37 mN/m. It is found that We ~ 30 and it is above the critical We. it seems inconsistent with the experiment with the soap film. The liquid film in the foam is quite thin (~100 nm) due to the negative osmotic pressure due to packing in the dry foam[28,29], while the soap film has several μm[31,32]. When the emitted droplet impacts the liquid film in the foam, the liquid film should be distorted. We consider two possibilities. One is that the liquid film is broken due to thinning by the distortion. The other is that the bending modulus is large due to resistance for thinning. Then, a large hole is created by an emitted droplet, while the soap film is easily bended. Thus we consider that the rapid droplet can break the liquid film in the foam. We should investigate how the droplet penetrates and the liquid film is broken in future.

Next, we discuss the other exponents, β and γ. To establish a theoretical approach, we consider the probability of a breakage being induced by one penetrating droplet. Since this depends on the droplet velocity and the local liquid fraction, this probability can be expressed as . Meanwhile, N is not the number of bubbles broken in the first collapse process but a summation over the entire process. Also note that the absorption of the broken liquid film should change the local liquid fraction, since the diffusion of the liquid in the foam is much slower than the CBC dynamics. Thus, we find that N is not a suitable parameter for a theoretical approach and we cannot relate to N. To reveal the origin of the exponents β and γ, we may need to investigate the first breakage process distinctly, separated from secondary and higher order processes to obtain an expression for .

Summary

To summarize, we perform in-situ observation of collective bubble collapse in a quasi-two-dimensional liquid-gas foam by using a high speed camera. It was found that there were two crucial modes for the CBC. One is the propagating mode, the liquid-film breaking due to the impact of absorption. The other mode is the penetrating mode, when liquid droplets are emitted due to the impact of absorption and penetrate through other, distant liquid films. We also investigate the CBC phenomenon for different liquid fractions. As the liquid fraction increases, the velocity of the liquid droplet decreases, making it harder to penetrate liquid films. Instead of penetration, the liquid droplets bounce due to the elasticity of the liquid film. When the velocity becomes slower, the liquid droplet is absorbed into the other liquid film. We find that the absorption velocity is proportional to the osmotic pressure; an explanation for the scaling of other parameters is left for future work, along with an examination of the local dynamics of the first breakage to further clarify the collective collapse process.

Materials and Methods

We used a 14 wt% solution of TTAB (tetradecyl trimethyl ammonium bromide) in glycerol and deionized water. In this article, we mainly show the results of foams with 17 wt% glycerol. The foams are created by using a capillary glass tube equipped with an air pump. We put the foam at the center of a glass plate and cover it with another glass plate. The sample thickness is controlled to 2.2 mm by using a spacer. This sample cell is horizontal; the foam is regarded as quasi two dimensional. The diameter of the bubbles is about 3~5 mm, and the number of the bubbles N is about 150~200 or approximately 500. We punctured a bubble outside the foam using a capillary glass needle with a very small amount of silicon grease on the tip to induce CBC collapse[20-22]. We took images during the collapsing process of the bubbles using a high speed camera (KEYENCE, VW-9000) with a frame rate between 4000 and 10000 fps. The temperature was set to around 20 °C using an air conditioner. Here, we determine the liquid fraction in quasi-two-dimensional foams by the method described in ref.[19]. It takes only a few minutes to perform the experiment, thus evaporation of the liquid is negligible. We also neglect drainage and coarsening.

Supplementary information

Movie for CBC process

Movie for bouncing process