A continuum description of the buckling of a line of spheres in a transverse harmonic confining potential

A line of contacting hard spheres, placed in a transverse confining potential, buckles under compression or when tilted away from the horizontal, once a critical tilt angle is exceeded. This interesting nonlinear problem is enriched by the combined application of both compression and tilt. In a continuous formulation, the profile of transverse sphere displacement is well described by numerical solutions of a second-order differential equation (provided that buckling is not of large amplitude). Here we provide a detailed discussion of these solutions, which are approximated by analytic expressions in terms of Jacobi, Whittaker and Airy functions. The analysis in terms of Whittaker functions yields an exact result for the critical tilt for buckling without compression.

Analysis of the foam-forming of non-woven lightweight fibrous materials using X-ray tomography

Abstract

Foam-forming has in the past predominantly been used to create two-dimensional sheet-like fibrous materials. Allowing the foam to drain freely and decay under gravity, rather than applying a vacuum to remove it rapidly, we can produce lightweight three-dimensional fibrous structures from cellulose fibres, of potential use for thermal and acoustic insulation. \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu$$\end{document}μCT scanning of the fibrous materials enable us to determine both void size distributions and also distributions of fibre orientations. Through image analysis and uniaxial compression testing, we find that the orientation of the fibres, rather than the size of the voids, determine the compressive strength of the material. The fibrous samples display a layering of the fibres perpendicular to the direction of drainage of the precursor liquid foam. This leads to an anisotropy of the compressive behaviour of the samples. Varying the initial liquid fraction of the foam allows for tuning of the compressive strength. We show an increase in over seven times can be achieved for samples of the same density (13 kg.m^-3).

Graphic abstract

The energy of fcc and hcp foams

We present Surface Evolver evaluations of the difference in energy between face-centred cubic (fcc) and hexagonal close-packed (hcp) foams in the usual idealized model, for liquid fractions ranging from the dry to the wet limit. The difference vanishes in both limits, and favours hcp for all intermediate liquid fractions, as has been proven. The maximum relative energy difference is very small, of the order of 10-5. The asymptotic dependence on liquid fraction is non-analytic in both limits: we present explicit expressions in both cases, derived from first principles. They have been obtained from identifying node interactions (dry limit) and contact interactions (wet limit) as the respective sources for energy differences between fcc and hcp. The wet limit is well described by Morse-Witten theory which has proven to be very powerful for the analytic computation of the surface energy of slightly deformed bubbles.

Theory of rotational columnar structures of soft spheres

There is a growing interest in cylindrical structures of hard and soft particles. A promising new method to assemble such structures has recently been introduced by Lee et al. [Lee, Gizynski, and Grzybowski, Adv. Mater. 29, 1704274 (2017)ADVMEW0935-964810.1002/adma.201704274]. They used rapid rotations around a central axis to drive spheres of lower density than the surrounding fluid towards the axis. This resulted in different structures as the number of spheres is varied. Here, we present comprehensive analytic energy calculations for such self-assembled structures, based on a generic soft sphere model, from which we obtain a phase diagram. It displays interesting features, including peritectoid points. These analytic calculations are complemented by preliminary numerical simulations for finite sample sizes with soft spheres. A similar analytic approach could be used to study packings of spheres inside cylinders of fixed dimensions, but with a variation in the number of spheres.

Ideal wet two-dimensional foams and emulsions with finite contact angle

We present simulations that show that the equilibrium structure of an ideal two-dimensional foam with a finite contact angle develops an inhomogeneity for high liquid fraction φ. In liquid-liquid emulsions this inhomogeneity is known as flocculation. In the case of an ordered foam this requires a perturbation, but in a disordered foam inhomogeneity grows steadily and spontaneously with φ, as demonstrated in our simulations performed with the Surface Evolver.

Bubble-bubble interactions in a 2d foam, close to the wet limit

Following the general approach of Morse and Witten for the deformation of a bubble in contact with neighbouring bubbles, we develop a model for contacting bubbles in two dimensions which can be solved analytically. The force-displacement relations are derived by elementary methods; unlike the case of 3d, no logarithmic factors arise in two dimensions. We also discuss the case of a uniform compression of a symmetric foam structure; the (osmotic) compressibility depends on the number of contacts, as was shown in earlier work by Lacasse et al. Our model, which is based on first principles, without any free parameters, may be extended to simulate 2d foams.

Simulation and observation of line-slip structures in columnar structures of soft spheres

We present the computed phase diagram of columnar structures of soft spheres under pressure, of which the main feature is the appearance and disappearance of line slips, the shearing of adjacent spirals, as pressure is increased. A comparable experimental observation is made on a column of bubbles under forced drainage, clearly exhibiting the expected line slip.

Corrected Article: Simulation and observation of line-slip structures in columnar structures of soft spheres [Phys. Rev. E 96, 012610 (2017)]

This corrects the article DOI: 10.1103/PhysRevE.96.012610.

Analysis of liquid metal foams through X-ray radioscopy and microgravity experiments

Complex liquid structures such as metallic foams were produced in a furnace that allowed in situ X-ray monitoring of the evolution of the structure and distribution of the liquid in the foam. The experiments were carried out during parabolic flights which provided varying levels of gravity. The evolution of the characteristic liquid fraction profiles due to gravity induced drainage was measured and analysed in terms of the foam drainage equation, obtaining viscosity and surface tension by fitting solutions of the equation to the experimental data. The surface tension of the melt in the foam was decreased up to 40%. Effective viscosities of up to 139 times the viscosity of a pure bulk melt were observed. These effects could be attributed to the smaller influence of solid particles dispersed in the melt and the larger influence of the complex foam structure.

Theoretical analysis of the performance of a foam fractionation column

A model system for theory and experiment which is relevant to foam fractionation consists of a column of foam moving through an inverted U-tube between two pools of surfactant solution. The foam drainage equation is used for a detailed theoretical analysis of this process. In a previous paper, we focused on the case where the lengths of the two legs are large. In this work, we examine the approach to the limiting case (i.e. the effects of finite leg lengths) and how it affects the performance of the fractionation column. We also briefly discuss some alternative set-ups that are of interest in industry and experiment, with numerical and analytical results to support them. Our analysis is shown to be generally applicable to a range of fractionation columns.

A model system for foam fractionation

A model system for theory and experiment that is relevant to foam fractionation consists of a column of foam moving through an inverted U-tube between two pools of surfactant solution. The foam drainage equation and its variants are used for a theoretical analysis of this process. In the limit in which the lengths of the two legs is large , exact analytic formulae may be derived for the key properties of the system. Numerical computations and experiments support these results.

Dense packings of spheres in cylinders: simulations

We study the optimal packing of hard spheres in an infinitely long cylinder, using simulated annealing, and compare our results with the analogous problem of packing disks on the unrolled surface of a cylinder. The densest structures are described and tabulated in detail up to D/d=2.873 (ratio of cylinder and sphere diameters). This extends previous computations into the range of structures which include internal spheres that are not in contact with the cylinder.

Foam as a complex system

What is a 'complex system'? The two-dimensional foam, as originally popularized by Cyril Stanley Smith, provides an ideal context in which to explore this question.

Steady drainage in emulsions: corrections for surface Plateau borders and a model for high aqueous volume fraction

We compare extensive experimental results for the gravity-driven steady drainage of oil-in-water emulsions with two theoretical predictions, both based on the assumption of Poiseuille flow. The first is from standard foam drainage theory, applicable at low aqueous volume fractions, for which a correction is derived to account for the effects of the confinement of the emulsion. The second arises from considering the permeability of a model porous medium consisting of solid sphere packings, applicable at higher aqueous volume fractions. We find quantitative agreement between experiment and the foam drainage theory at low aqueous volume fractions. At higher aqueous volume fractions, the reduced flow rate calculated from the permeability theory approaches the master curve of the experimental data. Our experimental data demonstrates the analogy between the problem of electrical flow and liquid flow through foams and emulsions.

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.

Two-dimensional foam rheology with viscous drag

We formulate and apply a continuum model that incorporates elasticity, yield stress, plasticity, and viscous drag. It is motivated by the two-dimensional foam rheology experiments of Debregeas et al. [Phys. Rev. Lett. 87, 178305 (2001)10.1103/PhysRevLett.87.178305] and Wang et al. [Phys. Rev. E 73, 031401 (2006)10.1103/PhysRevE.73.031401], and is successful in exhibiting their principal features, which are an exponentially decaying velocity profile and strain localization. Transient effects are also identified.

The crystal structure of bubbles in the wet foam limit

We have observed a rich variety of three-dimensional crystal and defect structures spontaneously formed by small (diameter 200 µm) bubbles in a wet foam. The observations confirm and extend those made by Bragg and Nye in 1947. However, while their experiments with two-dimensional bubble rafts have stimulated many researchers, their work on assemblages does not appear to have been followed up. These ordered packings now pose intriguing questions for the physics of foams. The bubbles seem too large for conventional thermodynamics and kinetics to easily explain the high degree of ordering.

Drainage induced convection rolls in foams I. Convective bubble motion in a tilted tube

When liquid is added to a foam at sufficiently large flow rates, convective bubble motion will occur. Experiments are described in which the foam is confined in a tube which is tilted from the vertical. The theory of foam drainage is applied to this problem to show that the critical angle of tilt theta(c) at which convection occurs is related to the liquid flow-rate Q by theta(c) proportional to Q(-3/4).

Two-dimensional viscous froth model for foam dynamics

The two-dimensional viscous froth model is a simple tractable model for foam rheology and coarsening. It includes, but is not confined to, the quasistatic regime. Here we present a detailed analysis and implementation of the model, illustrated with various examples. With certain simplifying assumptions, it provides significant insight into strain-rate-dependent effects in foam rheology and elsewhere, particularly in relation to recent experiments.

The foam/emulsion analogy in structure and drainage

The often quoted analogy between foams and emulsions is experimentally tested by studying properties after settling and under forced drainage of oil-in-water emulsions of drop size similar as for bubbles generally used in foam experiments. Observations with regard to structure, water fraction and drainage wave properties confirm the expected similarity in the low flow rate range. However, while for foams a convective circulation on the scale of the container sets in for values of water fraction exceeding about 0.2, no such convection is found in emulsions. Here instabilities are only encountered at water fractions of about 0.4, close to the void fraction of random packings of spheres. These take on the form of descending pulses of increased water fraction and lead to the transition from a frozen to a locally agitated structure.