Independent variables in foam clusters

We find, by counting the degrees of freedom of two-dimensional bubble clusters (finite or periodic) of given topology and bubble areas, that the Plateau laws determine a unique configuration of a finite free cluster, but allow an infinite number of configurations of a periodic cluster. Each of these configurations is associated with a particular strain (stress) state of the cluster; there is in general one unstrained configuration, which corresponds to the minimum of the (surface) energy. Configurations of given topology that satisfy Plateau's laws may only exist in certain ranges of bubble area ratios and/or strains. Received 31 May 2001 and Received in final form 12 September 2001     

Minimum perimeter partitions of the plane into equal numbers of regions of two different areas

We identify the minimum-perimeter periodic tilings of the plane by equal numbers of regions (cells) of areas 1 and λ (minimal tilings), with at most two cells of each area per period and for which all cells of the same area are topologically equivalent. For λ close to 1 the minimal tiling is hexagonal. For smaller values of λ the minimal tilings contain pairs of 5/7, 4/8 and 3/9 cells, the cells with fewer sides having smaller area. The correlation between area fraction and number of sides in the minimal tilings is approximately linear and consistent with Lewis' law. Received 27 June 2001 and Received in final form 29 August 2001     

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.     

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.     

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     

Randomness and a step-like distribution of pile heights in avalanche models

This paper considers a one-parameter family of sand-piles. The family exhibits the crossover between the models with deterministic and stochastic relaxation. The mean pile height is used to describe the crossover. The height densities corresponding to the models with relaxation of both types approach one another as the parameter increases. Relaxation is supposed to deal with the local losses of grains by a fixed amount. In that case the densities show a step-like behaviour in contrast to the peaked shape found in the models with the local loss of grains down to a fixed level [S. Lübeck, Phys. Rev. E 62, 6149 (2000)]. A spectral approach based on the long-run properties of the pile height considers the models with deterministic and random relaxation more accurately and distinguishes between the two cases for admissible parameter values.     

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.     

A foam film propagating in a confined geometry: Analysis via the viscous froth model

A single film (typical of a film in a foam) moving in a confined geometry (i.e. confined between closely spaced top and bottom plates) is analysed via the viscous froth model. In the first instance the film is considered to be straight (as viewed from above the top plate) but is not flat. Instead it is curved (with a circular arc cross-section) in the direction across the confining plates. This curvature leads to a maximal possible steady propagation velocity for the film, which is characterised by the curved film meeting the top and bottom plates tangentially. Next the film is considered to propagate in a channel (i.e. between top and bottom plates and sidewalls, with the sidewall separation exceeding that of the top and bottom plates). The film is now curved along as well as across the top and bottom plates. Curvature along the plates arises from viscous drag forces on the channel sidewall boundaries. The maximum steady propagation velocity is unchanged, but can now also be associated with films meeting channel sidewalls tangentially, a situation which should be readily observable if the film is viewed from above the top plate. Observed from above, however, the film need not appear as an arc of a circle. Instead the film may be relatively straight along much of its length, with curvature pushed into boundary layers at the sidewalls.     

Anomalous scaling in the Zhang model

We apply the moment analysis technique to analyze large scale simulations of the Zhang sandpile model. We find that this model shows different scaling behavior depending on the update mechanism used. With the standard parallel updating, the Zhang model violates the finite-size scaling hypothesis, and it also appears to be incompatible with the more general multifractal scaling form. This makes impossible its affiliation to any one of the known universality classes of sandpile models. With sequential updating, it shows scaling for the size and area distribution. The introduction of stochasticity into the toppling rules of the parallel Zhang model leads to a scaling behavior compatible with the Manna universality class. Received 8 August 2000 and Received in final form 4 October 2000     

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.     

Synchronization and structure in an adaptive oscillator network

We analyze the interplay of synchronization and structure evolution in an evolving network of phase oscillators. An initially random network is adaptively rewired according to the dynamical coherence of the oscillators, in order to enhance their mutual synchronization. We show that the evolving network reaches a small-world structure. Its clustering coefficient attains a maximum for an intermediate intensity of the coupling between oscillators, where a rich diversity of synchronized oscillator groups is observed. In the stationary state, these synchronized groups are directly associated with network clusters.     

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.     

Statics and dynamics of adhesion between two soap bubbles

An original set-up is used to study the adhesive properties of two hemispherical soap bubbles put into contact. The contact angle at the line connecting the three films is extracted by image analysis of the bubbles profiles. After the initial contact, the angle rapidly reaches a static value slightly larger than the standard 120^° angle expected from Plateau rule. This deviation is consistent with previous experimental and theoretical studies: it can be quantitatively predicted by taking into account the finite size of the Plateau border (the liquid volume trapped at the vertex) in the free energy minimization. The visco-elastic adhesion properties of the bubbles are further explored by measuring the deviation Δθ_d(t) of the contact angle from the static value as the distance between the two bubbles supports is sinusoidally modulated. It is found to linearly increase with Δr _c/r _c , where r_c is the radius of the central film and Δr _c the amplitude of modulation of this length induced by the displacement of the supports. The in-phase and out-of-phase components of Δθ_d(t) with the imposed modulation frequency are systematically probed, which reveals a transition from a viscous to an elastic response of the system with a crossover pulsation of the order 1rad · s^-1. Independent interfacial rheological measurements, obtained from an oscillating bubble experiment, allow us to develop a model of dynamic adhesion which is confronted to our experimental results. The relevance of such adhesive dynamic properties to the rheology of foams is briefly discussed using a perturbative approach to the Princen 2D model of foams.     

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.     

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.     

Fractal formations in the Lattice Limit Cycle model

We examine the fractal patterns arising in the Lattice Limit Cycle model, when it is restricted on square and fractal lattices. We show that, for processes taking place on regular 2d substrates, the fractal dimensions depend on the kinetic constants and we have observed a sharp phase-transition from uniform 2d spatial distributions (d_f=2) when the kinetic parameters are near the Hopf bifurcation point, to a inside the limit cycle region. For processes taking place on substrates which contain inactive sites, we observe nucleation of homologous species around inactive regions leading to poisoning, when the active sites are distributed in a fractal manner on the substrate. This is less frequent in cases where the active sites are distributed uniformly and randomly on the lattice leading, normally, to non-trivial steady states.     

The emergence of coordination in public good games

In physical models it is well understood that the aggregate behaviour of a system is not in one to one correspondence with the behaviour of the average individual element of that system. Yet, in many economic models the behaviour of aggregates is thought of as corresponding to that of an individual. A typical example is that of public goods experiments. A systematic feature of such experiments is that, with repetition, people contribute less to public goods. A typical explanation is that people “learn to play Nash” or something approaching it. To justify such an explanation, an individual learning model is tested on average or aggregate data. In this paper we will examine this idea by analysing average and individual behaviour in a series of public goods experiments. We analyse data from a series of games of contributions to public goods and as is usual, we test a learning model on the average data. We then look at individual data, examine the changes that this produces and see if some general model such as the EWA (Expected Weighted Attraction) with varying parameters can account for individual behaviour. We find that once we disaggregate data such models have poor explanatory power. Groups do not learn as supposed, their behaviour differs markedly from one group to another, and the behaviour of the individuals who make up the groups also varies within groups. The decline in aggregate contributions cannot be explained by resorting to a uniform model of individual behaviour. However, the Nash equilibrium of such a game is a total payment for all the individuals and there is some convergence of the group in this respect. Yet the individual contributions do not converge. How the individuals “self-organsise” to coordinate, even in this limited way remains to be explained.     

Collapse dynamics of smectic-A bubbles

The collapse dynamics of smectic-A bubbles are analyzed experimentally and theoretically. Each bubble is expanded from a flat film stretched at the end of a hollow cylinder and deflated through a pressure release by means of a capillary tube. Its total collapse time can be varied between 0.1s and 20s by suitably choosing the length and the internal diameter of the capillary. This experiment allowed us to show that the collapse takes place in two steps: an initial one, which lasts a fraction of a second, where the meniscus destabilizes and fills up with focal conics, followed by a much longer period during which the bubble collapses and exchanges material with the meniscus. By measuring simultaneously the Laplace pressure and the internal pressure inside the bubble, we were able to fully characterize the shear-thinning behavior of the smectic phase within the meniscus. We emphasize that this method is generic and could be applied as well to other systems such as soap bubbles, on condition that inertial effects are negligible.     

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.     

On some natural and model 2D bimodal random cellular structures

The topological and metric properties of a few natural 2D random cellular structures, namely an armadillo shell structure and young soap froths, which are formed from two classes of cells, large and small, have been characterized. The topological properties of a model generated from a Kagome tiling, which mimics such random binary structures, have also been exactly calculated. The distribution of the number of cell sides is bimodal for the structures investigated. In contrast to the classical Aboav-Weaire law for homogeneous 2D random cellular structures, nm(n), the mean total number of edges of neighbouring cells of cells with n sides does not vary linearly with n. Only the nm(i, n) (i=1,2) determined separately for every class of cells are linear in n for all investigated structures. Topological properties and correlations between metric and topological properties are finally compared with the predictions of various literature models. Received: 24 December 1997 / Revised: 7 April 1998 / Accepted: 20 April 1998