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     

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.     

Stability of fluid flow past a membrane

The stability of the flow of a fluid past a solid membrane of infinitesimal thickness is investigated using a linear stability analysis. The system consists of two fluids of thicknesses R and H R and bounded by rigid walls moving with velocities and , and separated by a membrane of infinitesimal thickness which is flat in the unperturbed state. The fluids are described by the Navier-Stokes equations, while the constitutive equation for the membrane incorporates the surface tension, and the effect of curvature elasticity is also examined for a membrane with no surface tension. The stability of the system depends on the dimensionless strain rates and in the two fluids, which are defined as and for a membrane with surface tension , and and for a membrane with zero surface tension and curvature elasticity K. In the absence of fluid inertia, the perturbations are always stable. In the limit , the decay rate of the perturbations is O(k ³ ) smaller than the frequency of the fluctuations. The effect of fluid inertia in this limit is incorporated using a small wave number asymptotic analysis, and it is found that there is a correction of smaller than the leading order frequency due to inertial effects. This correction causes long wave fluctuations to be unstable for certain values of the ratio of strain rates and ratio of thicknesses H. The stability of the system at finite Reynolds number was calculated using numerical techniques for the case where the strain rate in one of the fluids is zero. The stability depends on the Reynolds number for the fluid with the non-zero strain rate, and the parameter , where is the surface tension of the membrane. It is found that the Reynolds number for the transition from stable to unstable modes, , first increases with , undergoes a turning point and a further increase in the results in a decrease in . This indicates that there are unstable perturbations only in a finite domain in the plane, and perturbations are always stable outside this domain. Received: 29 May 1997 / Revised: 9 October 1997 / Accepted: 26 November 1997     

Quasi-one-dimensional foam drainage

Foam drainage is considered in a froth flotation cell. Air flow through the foam is described by a simple two-dimensional deceleration flow, modelling the foam spilling over a weir. Foam microstructure is given in terms of the number of channels (Plateau borders) per unit area, which scales as the inverse square of bubble size. The Plateau border number density decreases with height in the foam, and also decreases horizontally as the weir is approached. Foam drainage equations, applicable in the dry foam limit, are described. These can be used to determine the average cross-sectional area of a Plateau border, denoted A, as a function of position in the foam. Quasi-one-dimensional solutions are available in which A only varies vertically, in spite of the two-dimensional nature of the air flow and Plateau border number density fields. For such situations the liquid drainage relative to the air flow is purely vertical. The parametric behaviour of the system is investigated with respect to a number of dimensionless parameters: K (the strength of capillary suction relative to gravity), α (the deceleration of the air flow), and n and h (respectively, the horizontal and vertical variations of the Plateau border number density). The parameter K is small, implying the existence of boundary layer solutions: capillary suction is negligible except in thin layers near the bottom boundary. The boundary layer thickness (when converted back to dimensional variables) is independent of the height of the foam. The deceleration parameter α affects the Plateau border area on the top boundary: weaker decelerations give larger Plateau border areas at the surface. For weak decelerations, there is rapid convergence of the boundary layer solutions at the bottom onto ones with negligible capillary suction higher up. For strong decelerations, two branches of solutions for A are possible in the K = 0 limit: one is smooth, and the other has a distinct kink. The full system, with small but non-zero capillary suction, lies relatively close to the kinked solution branch, but convergence from the lower boundary layer onto this branch is distinctly slow. Variations in the Plateau border number density (non-zero n and h) increase individual Plateau border areas relative to the case of uniformly sized bubbles. For strong decelerations and negligible capillarity, solutions closely follow the kinked solution branch if bubble sizes are only slightly non-uniform. As the extent of non-uniformity increases, the Plateau border area reaches a maximum corresponding to no net upward velocity of foam liquid. In the case of vertical variation of number density, liquid content profiles and Plateau border area profiles cease to be simply proportional to one another. Plateau border areas match at the top of the foam independent of h, implying a considerable difference in liquid content for foams which exhibit different number density profiles. Received 3 July 2001     

Asymptotic analysis of wall modes in a flexible tube

The stability of wall modes in a flexible tube of radius R surrounded by a viscoelastic material in the region R < r < H R in the high Reynolds number limit is studied using asymptotic techniques. The fluid is a Newtonian fluid, while the wall material is modeled as an incompressible visco-elastic solid. In the limit of high Reynolds number, the vorticity of the wall modes is confined to a region of thickness in the fluid near the wall of the tube, where the small parameter , and the Reynolds number is , and are the fluid density and viscosity, and V is the maximum fluid velocity. The regime is considered in the asymptotic analysis, where G is the shear modulus of the wall material. In this limit, the ratio of the normal stress and normal displacement in the wall, , is only a function of H and scaled wave number . There are multiple solutions for the growth rate which depend on the parameter .In the limit , which is equivalent to using a zero normal stress boundary condition for the fluid, all the roots have negative real parts, indicating that the wall modes are stable. In the limit , which corresponds to the flow in a rigid tube, the stable roots of previous studies on the flow in a rigid tube are recovered. In addition, there is one root in the limit which does not reduce to any of the rigid tube solutions determined previously. The decay rate of this solution decreases proportional to in the limit , and the frequency increases proportional to . Received: 5 November 1997 / Revised: 10 March 1998 / Accepted: 29 April 1998     

Effect of tangential interface motion on the viscous instability in fluid flow past flexible surfaces

The stability of linear shear flow of a Newtonian fluid past a flexible membrane is analysed in the limit of low Reynolds number as well as in the intermediate Reynolds number regime for two different membrane models. The objective of this paper is to demonstrate the importance of tangential motion in the membrane on the stability characteristics of the shear flow. The first model assumes the wall to be a “spring-backed” plate membrane, and the displacement of the wall is phenomenologically related in a linear manner to the change in the fluid stresses at the wall. In the second model, the membrane is assumed to be a two-dimensional compressible viscoelastic sheet of infinitesimal thickness, in which the constitutive relation for the shear stress contains an elastic part that depends on the local displacement field and a viscous component that depends on the local velocity in the membrane. The stability characteristics of the laminar flow in the limit of low are crucially dependent on the tangential motion in the membrane wall. In both cases, the flow is stable in the low Reynolds number limit in the absence of tangential motion in the membrane. However, the presence of tangential motion in the membrane destabilises the shear flow even in the absence of fluid inertia. In this case, the non-dimensional velocity (Λ_t) required for unstable fluctuations is proportional to the wavenumber k ( Λ _t∼k) in the plate membrane type of wall while it scales as k² in the viscoelastic membrane type of wall ( Λ _t∼k ²) in the limit k→ 0. The results of the low Reynolds number analysis are extended numerically to the intermediate Reynolds number regime for the case of a viscoelastic membrane. The numerical results show that for a given set of wall parameters, the flow is unstable only in a finite range of Reynolds number, and it is stable in the limit of large Reynolds number. Received 8 November 2000 and Received in final form 20 March 2001     

Asymptotic analysis of wall modes in a flexible tube revisited

The stability of wall modes in fluid flow through a flexible tube of radius R surrounded by a viscoelastic material in the region R < r < HR is analysed using a combination of asymptotic and numerical methods. The fluid is Newtonian, while the flexible wall is modelled as an incompressible viscoelastic solid. In the limit of high Reynolds number (Re), the vorticity of the wall modes is confined to a region of thickness O(Re ^-1/3) in the fluid near the wall of the tube. Previous numerical studies on the stability of Hagen-Poiseuille flow in a flexible tube to axisymmetric disturbances have shown that the flow could be unstable in the limit of high Re, while previous high Reynolds number asymptotic analyses have revealed only stable modes. To resolve this discrepancy, the present work re-examines the asymptotic analysis of wall modes in a flexible tube using a new set of scaling assumptions. It is shown that wall modes in Hagen-Poiseuille flow in a flexible tube are indeed unstable in the limit of high Re in the scaling regime Re∼Σ^3/4. Here Σ is a nondimensional parameter characterising the elasticity of the wall, and Σ≡ρGR ²/η², where ρ and η are the density and viscosity of the fluid, and G is the shear modulus of the wall medium. The results from the present asymptotic analysis are in excellent agreement with the previous numerical results. Importantly, the present work shows that the different types of unstable modes at high Reynolds number reported in previous numerical studies are qualitatively the same: they all belong to the class of unstable wall modes predicted in this paper. Received 12 June 2000 and Received in final form 8 November 2000     

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.     

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     

Monitoring the buckling threshold of drying colloidal droplets using water-ethanol mixtures

We visualize the drying of droplets of colloids suspended in a mixture of two miscible solvents, namely water and ethanol. After a period of isotropic shrinkage, droplets suddenly buckle like elastic shells. For a fixed colloid solid fraction, the buckling threshold evolves as a function of ethanol content, due to changes of the solvent mixture physical properties, such as viscosity and evaporation rate. A simplified model predicting the qualitative behavior of the buckling threshold as a function of the initial ethanol mass fraction has been developed that fits well experimental results.     

Gravitational compression of colloidal gels

We study the compression of depletion gels under the influence of a gravitational stress by monitoring the time evolution of the gel interface and the local volume fraction, , inside the gel. We find is not constant throughout the gel. Instead, there is a volume fraction gradient that develops and grows along the gel height as the compression process proceeds. Our results are correctly described by a non-linear poroelastic model that explicitly incorporates the -dependence of the gravitational, elastic and viscous stresses acting on the gel.     

Dissipative flows of 2D foams

We analyze the flow of a liquid foam between two plates separated by a gap of the order of the bubble size (2D foam). We concentrate on the salient features of the flow that are induced by the presence, in an otherwise monodisperse foam, of a single large bubble whose size is one order of magnitude larger than the average size. We describe a model suited for numerical simulations of flows of 2D foams made up of a large number of bubbles. The numerical results are successfully compared to analytical predictions based on scaling arguments and on continuum medium approximations. When the foam is pushed inside the cell at a controlled rate, two basically different regimes occur: a plug flow is observed at low flux whereas, above a threshold, the large bubble migrates faster than the mean flow. The detailed characterization of the relative velocity of the large bubble is the essential aim of the present paper. The relative velocity values, predicted both from numerical and from analytical calculations that are discussed here in great detail, are found to be in fair agreement with experimental results from the preprint Experimental evidence of flow destabilization in a 2D bidisperse foam by the present authors (2005).     

Uniqueness, stability and Hessian eigenvalues for two-dimensional bubble clusters

A recent conjecture on two-dimensional foams suggested that for fixed topology with given bubble areas there is a unique state of stable equilibrium. We present counter-examples, consisting of a ring of bubbles around a central one, which refute this conjecture. The discussion centres on a novel form of instability which causes symmetric clusters to become distorted. The stability of these bubble clusters is examined in terms of the Hessian of the energy. Received 8 November 2001     

Effects of the rate of evaporation and film thickness on nonuniform drying of film-forming concentrated colloidal suspensions

In this paper, we report on nonuniform distribution of film-forming waterborne colloidal suspensions above the critical concentration _c of the colloidal glass transition during drying. We found that colloidal suspension films dry nonuniformly when the initial rate of evaporation E and/or the initial thickness l₀ are high. We found that a Peclet number Pe, defined as Pe = El₀/D, where D is the diffusion coefficient of the colloids in the diluted suspensions, does not predict uniformity of drying of the concentrated suspensions, contrary to the reported work on drying of diluted suspensions. Since the colloidal particles are crowded and their diffusive motion is restricted in concentrated suspensions, we assumed that above _c water is transported to the drying surface by hydrodynamic flow along the osmotic pressure gradient. The permeability of water through channels between deforming particles is estimated by adapting the theory of foam drainage. We defined a new Peclet number Pe by substituting the transport coefficient of flow (defined as the permeability divided by the viscosity, multiplied by the osmotic pressure gradient) for the diffusion coefficient. This extended Peclet number predicted the nonuniform drying with a criterion of Pe > 1. These results indicate that the mechanism of water transport to the drying surface in concentrated suspensions is water permeation by osmotic pressure, which is faster than mutual diffusion between water and particles --that has been observed in diluted suspensions and discussed by Routh and Russel. The theory fits well the experimental drying curves for various thicknesses and rates of evaporation. The particle distribution in the drying films is also estimated and it is indicated that the latex distribution is nonuniform when Pe > 1.     

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.     

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.     

Tracer dispersion in power law fluids flow through porous media: Evidence of a cross-over from a logarithmic to a power law behaviour

An analytical model is presented to describe the dispersion of tracers in a power-law fluid flowing through a statistically homogeneous and isotropic porous medium. The model is an extension of Saffman's approach to the case of non-Newtonian fluids. It is shown that an effective viscosity depending on the pressure gradient and on the characteristics of the fluid, must be introduced to satisfy Darcy's law. An analytical expression of the longitudinal dispersivity is given as a function of the Peclet number Pe and of the power-law index n that characterizes the dependence of the viscosity on the shear rate . As the flow velocity increases the dispersivity obeys an asymptotic power law: . This asymptotic regime is achieved at moderate Peclet numbers with strongly non-Newtonian fluids and on the contrary at very large values when n goes to 1 ( for n=0.9). This reflects the cross-over from a scaling behaviour for towards a logarithmic behaviour predicted for Newtonian fluids (n=1). Received: 22 July 1997 / Revised and Accepted: 2 July 1998