Mechanical instabilities of bubble clusters between parallel walls

We have carried out a systematic study of buckling-like mechanical instabilities in simple two- (2D) and three-dimensional (3D) symmetric foam clusters sandwiched between parallel planar walls. These instabilities occur when the wall separation w is reduced below a critical value, w^*, for which the foam surface energy E reaches its minimum, E^*. The clusters under investigation consist of either a single bubble, or of twin bubbles of fixed equal sizes (areas A in 2D or volumes V in 3D), which are either free to slide or pinned at the confining walls. We have numerically obtained w^* for both free and pinned 2D and 3D clusters. Furthermore, we have calculated the buckled configurations of 2D twin bubbles, either free or pinned, and of 3D free twin bubbles, whose energy is independent of w and equal to the minimum energy E^* of the unbuckled state. Finally, we have also predicted the critical w_t^* at which the terminal configurations under extension of 2D and 3D single and twin bubbles are realised. Experimental illustrations of these transitions under compression and extension are presented. Our results, together with others from the literature, suggest that a bubble cluster bounded by two parallel walls is stable only if the normal force it exerts on the walls is attractive, i.e., if dE/dw > 0; clusters that cause repulsion between the walls are unstable. We correlate this with the distribution of film orientations: films in a stable cluster cannot be too parallel to the confining walls; rather, their average tilt must be larger than for a random distribution of film orientations.     

Symmetry-breaking transitions and dissociation of two-dimensional Plateau borders

We have experimentally studied the dissociation/coalescence of internal Plateau borders (PBs) in simple monolayer bubble clusters, as a result of changing the liquid fraction. At large liquid content, the clusters consist of n bubbles of the same size, symmetrically placed around an internal n-sided PB (n-PB). On decreasing the liquid fraction we observed symmetry-breaking transitions in the 4- and 5-bubble clusters (but not in the 3-bubble cluster), followed by dissociation of the PBs. We used the Surface Evolver to determine the various equilibrium configurations of the corresponding two-dimensional wet clusters and their surface energies. The sequence of 4-bubble cluster configurations observed on varying the liquid fraction correlates qualitatively with that predicted on the basis of Surface Evolver calculations. The same is not true of the 5-bubble cluster.     

Foam imbibition in microgravity

We report an experimental study of aqueous foam imbibition in microgravity with strict mass conservation. The foam is in a Hele-Shaw cell. The bubble edge width ℓ is measured by image analysis. The penetration of the liquid in the foam, the foam imbibition, the foam inflation, and the rigidity loss are shown all to obey strict diffusion processes. The motion of bubbles needed for the foam inflation is a slow two-dimensional process with respect to the one-dimensional capillary rise of liquid. The foam is found to imbibes faster than it inflates. Received 20 May 2002 / Received in final form 21 January 2003 Published online 23 May 2003 RID="a" ID="a"e-mail: herve.caps@ulg.ac.be     

Decorated vertices with 3-edged cells in 2D foams: exact solutions and properties

The energy, area and excess energy of a decorated vertex in a 2D foam are calculated. The general shape of the vertex and its decoration are described analytically by a reference pattern mapped by a parametric Moebius transformation. A single parameter of control allows to describe, in a common framework, different types of decorations, by liquid triangles or 3-sided bubbles, and other non-conventional cells. A solution is proposed to explain the stability threshold in the flower problem.     

Instabilities in foam clusters related to the pressures in the bubbles

We reinterpret an instability in a two-dimensional free foam cluster previously discussed by Weaire et al. (Eur. Phys. J. E 7, 123 (2002)) in terms of the excess pressure in the bubbles. We conclude that in a free foam cluster no bubble can have a pressure below that of the surrounding gas. Received: 20 September 2002 / Accepted: 5 May 2003 / Published online: 27 May 2003 RID="a" ID="a"e-mail: fatima.vaz@ist.utl.pt     

Discrete rearranging disordered patterns, part II: 2D plasticity, elasticity and flow of a foam

The plastic flow of a foam results from bubble rearrangements. We study their occurrence in experiments where a foam is forced to flow in 2D: around an obstacle; through a narrow hole; or sheared between rotating disks. We describe their orientation and frequency using a topological matrix defined in the companion paper (F. Graner, B. Dollet, C. Raufaste, and P. Marmottant, this issue, 25 (2008) DOI 10.1140/epje/i2007-10298-8), which links them with continuous plasticity at large scale. We then suggest a phenomenological equation to predict the plastic strain rate: its orientation is determined from the foam's local elastic strain; and its rate is determined from the foam's local elongation rate. We obtain a good agreement with statistical measurements. This enables us to describe the foam as a continuous medium with fluid, elastic and plastic properties. We derive its constitutive equation, then test several of its terms and predictions.     

Mixing and sorting of bidisperse two-dimensional bubbles

We have examined a number of candidates for the minimum-surface-energy arrangement of two-dimensional clusters composed of N bubbles of area 1 and N bubbles of area λ ( λ≤1). These include hexagonal bubbles sorted into two monodisperse honeycomb tilings, and various mixed periodic tilings with at most four bubbles per unit cell. We identify, as a function of λ, the minimal configuration for N → ∞. For finite N, the energy of the external (i.e., cluster-gas) boundary and that of the interface between honeycombs in “phase-separated” clusters have to be taken into account. We estimate these contributions and find the lowest total energy configuration for each pair (N,λ). As λ is varied, this alternates between a circular cluster of one of the mixed tilings, and “partial wetting” of the monodisperse honeycomb of bubble area 1 by the monodisperse honeycomb of bubble area λ. Received 1 August 2002 RID="a" ID="a"e-mail: paulo@ist.utl.pt     

An elastic,plastic, viscous model for slow shear of a liquid foam

We suggest a scalar model for deformation and flow of an amorphous material such as a foam or an emulsion. To describe elastic, plastic and viscous behaviours, we use three scalar variables: elastic deformation, plastic deformation rate and total deformation rate; and three material-specific parameters: shear modulus, yield deformation and viscosity. We obtain equations valid for different types of deformations and flows slower than the relaxation rate towards mechanical equilibrium. In particular, they are valid both in transient or steady flow regimes, even at large elastic deformation. We discuss why viscosity can be relevant even in this slow shear (often called “quasi-static”) limit. Predictions of the storage and loss moduli agree with the experimental literature, and explain with simple arguments the non-linear large amplitude trends.     

The response of 2D foams to continuous applied shear in a Couette rheometer

The continuum model that has reproduced the experimental observation of exponential shear localisation for straight-edge boundary conditions is adapted to the case of circular geometry. Essentially the same effect is found. However, the scenario of possible velocity profiles is much richer. Our calculations elucidate many recent experiments qualitatively and suggest further extensions of them. Various limits are analysed. In particular, the localisation length vanishes as the inner-boundary velocity tends to zero.     

Yield drag in a two-dimensional foam flow around a circular obstacle: Effect of liquid fraction

We study the two-dimensional flow of foams around a circular obstacle within a long channel. In experiments, we confine the foam between liquid and glass surfaces. In simulations, we use a deterministic software, the Surface Evolver, for bubble details and a stochastic one, the extended Potts model, for statistics. We adopt a coherent definition of liquid fraction for all studied systems. We vary it in both experiments and simulations, and determine the yield drag of the foam, that is, the force exerted on the obstacle by the foam flowing at very low velocity. We find that the yield drag is linear over a large range of the ratio of obstacle to bubble size, and is independent of the channel width over a large range. Decreasing the liquid fraction, however, strongly increases the yield drag; we discuss and interpret this dependence.     

Dry microfoams: formation and flow in a confined channel

We present an experimental investigation of the agglomeration of microbubbles into a 2D microfoam and its flow in a rectangular microchannel. Using a flow-focusing method, we produce the foam in situ on a microfluidic chip for a large range of liquid fractions, down to a few percent in liquid. We can monitor the transition from separated bubbles to the desired microfoam, in which bubbles are closely packed and separated by thin films. We find that bubble formation frequency is limited by the liquid flow rate, whatever the gas pressure. The formation frequency creates a modulation of the foam flow, rapidly damped along the channel. The average foam flow rate depends non-linearly on the applied gas pressure, displaying a threshold pressure due to capillarity. Strong discontinuities in the flow rate appear when the number of bubbles in the channel width changes, reflecting the discrete nature of the foam topology. We also produce an ultra flat foam, reducing the channel height from 250 μm to 8 μm, resulting in a height to diameter ratio of 0.02; we notice a marked change in bubble shape during the flow.     

Foams in contact with solid boundaries: Equilibrium conditions and conformal invariance

A liquid foam in contact with a solid surface forms a two-dimensional foam on the surface. We derive the equilibrium equations for this 2D foam when the solid surface is curved and smooth, generalising the standard case of flat Hele-Shaw cells. The equilibrium conditions at the vertices in 2D, at the edges in 3D, are invariant by conformal transformations. Regarding the films, conformal invariance only holds with restrictions, which we explicit for 3D and flat 2D foams. Considering foams confined in thin interstices between two non-parallel plates, normal incidence and Laplace’s law lead to an approximate equation relating the plate profile to the conformal map. Solutions are given for the logarithm and power laws in the case of constant pressure. The paper concludes on a comparison with available experimental data.     

An elasto-visco-plastic model for immortal foams or emulsions

A variety of complex fluids consists in soft, round objects (foams, emulsions, assemblies of copolymer micelles or of multilamellar vesicles--also known as onions). Their dense packing induces a slight deviation from their prefered circular or spherical shape. As a frustrated assembly of interacting bodies, such a material evolves from one conformation to another through a succession of discrete, topological events driven by finite external forces. As a result, the material exhibits a finite yield threshold. The individual objects usually evolve spontaneously (colloidal diffusion, object coalescence, molecular diffusion), and the material properties under low or vanishing stress may alter with time, a phenomenon known as aging. We neglect such effects to address the simpler behaviour of (uncommon) immortal fluids: we construct a minimal, fully tensorial, rheological model, equivalent to the (scalar) Bingham model. Importantly, the model consistently describes the ability of such soft materials to deform substantially in the elastic regime (be it compressible or not) before they undergo (incompressible) plastic creep--or viscous flow under even higher stresses.     

Soft Dynamics simulation. 1. Normal approach of two deformable particles in a viscous fluid and optimal-approach strategy

Discrete simulation methods are efficient tools to investigate the behaviors of complex fluids such as dry granular materials or dilute suspensions of hard particles. By contrast, materials made of soft and/or concentrated units (emulsions, foams, vesicles, dense suspensions) can exhibit both significant elastic particle deflections (Hertz-like response) and strong viscous forces (squeezed liquid). We point out that the gap between two particles is then not determined solely by the positions of their centers, but rather exhibits its own dynamics. We provide the first ingredients of a new discrete numerical method, named Soft Dynamics, to simulate the combined dynamics of particles and contacts. As an illustration, we present the results for the approach of two particles. We recover the scaling behaviors expected in three limits: the Stokes limit for very large gaps, the Poiseuille-lubricated limit for small gaps and even smaller surface deflections, and the Hertz limit for significant surface deflections. We find that for each gap value, an optimal force achieves the fastest approach velocity. The principle of larger-scale simulations with this new method is provided. They will consitute a promising tool for investigating the collective behaviors of many complex materials.     

Shear modulus of two-dimensional foams: The effect of area dispersity and disorder

We use the Surface Evolver to determine the shear modulus G of a dry 2D foam of 2500 bubbles, using both extensional and simple shear. We examine G for a range of monodisperse, bidisperse and polydisperse foams, and relate it to various measures of the structural disorder of each foam. In all cases, the shear modulus of a foam decreases with increasing disorder.     

Viscosity effects in foam drainage: Newtonian and non-newtonian foaming fluids

We have studied the drainage of foams made from Newtonian and non-Newtonian solutions of different viscosities. Forced-drainage experiments first show that the behavior of Newtonian solutions and of shear-thinning ones (foaming solutions containing either Carbopol or Xanthan) are identical, provided one considers the actual viscosity corresponding to the shear rate found inside the foam. Second, for these fluids, a drainage regime transition occurs as the bulk viscosity is increased, illustrating a coupling between surface and bulk flow in the channels between bubbles. The properties of this transition appear different from the ones observed in previous works in which the interfacial viscoelasticity was varied. Finally, we show that foams made of solutions containing long flexible PolyEthylene Oxide (PEO) molecules counter-intuitively drain faster than foams made with Newtonian solutions of the same viscosity. Complementary experiments made with fluids having all the same viscosity but different responses to elongational stresses (PEO-based Boger fluids) suggest an important role of the elastic properties of the PEO solutions on the faster drainage.     

Uniform foam production by turbulent mixing: new results on free drainage vs. liquid content

An apparatus is described for rapidly producing large quantities of foam via turbulent mixing of gas with a narrow jet of a surfactant solution inside a delivery tube. By controlling relative flow rates, the gas volume fraction in the resulting foam may be easily varied across . Using such foams, we present a comprehensive set of data for free drainage as a systematic function of gas fraction and sample geometry. The qualitative behavior can be understood in terms of simple theoretical considerations, emphasizing the importance of controlling the initial foam conditions. Quantitative features are compared with two approximate versions of the drainage equation, highlighting the crucial role of capillarity for very dry foams and small samples. Received 15 February 1999     

Physicochemical approach to the theory of foam drainage

We have investigated theoretically the effect of surface viscoelasticity on the drainage of an aqueous foam. Former theories consider that the flow in Plateau borders is either Poiseuille flow or plug-flow. In the last case, the dissipation is attributed to flow in the nodes connecting Plateau borders. Although we do not include this dissipation in our model, we obtain a drainage equation that includes terms equivalent to those of the earlier models. We show that when the water solubility of the surfactant stabilizing the foam is low, the control parameter M for the transition between the two regimes is the ratio , where μ is the bulk viscosity, D _s the surface diffusion coefficient, R the radius of curvature of the Plateau border and ɛ the surface elasticity. When the surfactant is more soluble M is rather related to the bulk diffusion coefficient. Within the frame of this approach and in view of the estimated M values, we show that the flow in Plateau borders is Poiseuille-like. Received 26 June 2001