Foam Micromechanics

Foam evokes many different images: waves breaking at the seashore, the head on a pint of Guinness, an elegant dessert, shaving, the comfortable cushion on which you may be seated... From the mundane to the high tech, foams, emulsions, and cellular solids encompass a broad range of materials and applications. Soap suds, mayonnaise, and foamed polymers provide practical motivation and only hint at the variety of materials at issue. Typical of mukiphase materiaIs, the rheoIogy or mechanical behavior of foams is more complicated than that of the constituent phases alone, which may be gas, liquid, or solid. For example, a soap froth exhibits a static shear modulus-a hallmark of an elastic solid-even though it is composed primarily of two Newtonian fluids (water and air), which have no shear modulus. This apparent paradox is easily resolved. Soap froth contains a small amount of surfactant that stabilizes the delicate network of thin liq- uid films against rupture. The soap-film network deforms in response to a macroscopic strain; this increases interracial area and the corresponding sur- face energy, and provides the strain energy of classical elasticity theory [1]. This physical mechanism is easily imagined but very challenging to quantify for a realistic three-dimensional soap froth in view of its complex geome- try. Foam micromechanics addresses the connection between constituent properties, cell-level structure, and macroscopic mechanical behavior. This article is a survey of micromechanics applied to gas-liquid foams, liquid-liquid emulsions, and cellular solids. We will focus on static response where the foam deformation is very slow and rate-dependent phenomena such as viscous flow can be neglected. This includes nonlinear elasticity when deformations are large but reversible. We will also discuss elastic- plastic behavior, which involves yield phenomena. Foam structures based on polyhedra packed to fill space provide a unify- ing geometrical theme. Because a two-dimensional situation is always easier to visualize and usually easier to analyze, the roots of foam micromechanics lie in the plane packed with polygons. There are striking similarities as well as obvious differences between 2D and 3D

Worst case and probabilistic analysis of the 2-Opt algorithm for the TSP

2-Opt is probably the most basic local search heuristic for the TSP. This heuristic achieves amazingly good results on “real world” Euclidean instances both with respect to running time and approximation ratio. There are numerous experimental studies on the performance of 2-Opt. However, the theoretical knowledge about this heuristic is still very limited. Not even its worst case running time on 2-dimensional Euclidean instances was known so far. We clarify this issue by presenting, for every p∈N , a family of L p instances on which 2-Opt can take an exponential number of steps. Previous probabilistic analyses were restricted to instances in which n points are placed uniformly at random in the unit square [0,1]2, where it was shown that the expected number of steps is bounded by O~(n10) for Euclidean instances. We consider a more advanced model of probabilistic instances in which the points can be placed independently according to general distributions on [0,1] d , for an arbitrary d≥2. In particular, we allow different distributions for different points. We study the expected number of local improvements in terms of the number n of points and the maximal density ϕ of the probability distributions. We show an upper bound on the expected length of any 2-Opt improvement path of O~(n4+1/3⋅ϕ8/3) . When starting with an initial tour computed by an insertion heuristic, the upper bound on the expected number of steps improves even to O~(n4+1/3−1/d⋅ϕ8/3) . If the distances are measured according to the Manhattan metric, then the expected number of steps is bounded by O~(n4−1/d⋅ϕ) . In addition, we prove an upper bound of O(ϕ√d) on the expected approximation factor with respect to all L p metrics. Let us remark that our probabilistic analysis covers as special cases the uniform input model with ϕ=1 and a smoothed analysis with Gaussian perturbations of standard deviation σ with ϕ∼1/σ d

The Saffman-Taylor problem on a sphere

The Saffman-Taylor problem addresses the morphological instability of an interface separating two immiscible, viscous fluids when they move in a narrow gap between two flat parallel plates (Hele-Shaw cell). In this work, we extend the classic Saffman-Taylor situation, by considering the flow between two curved, closely spaced, concentric spheres (spherical Hele-Shaw cell). We derive the mode-coupling differential equation for the interface perturbation amplitudes and study both linear and nonlinear flow regimes. The effect of the spherical cell (positive) spatial curvature on the shape of the interfacial patterns is investigated. We show that stability properties of the fluid-fluid interface are sensitive to the curvature of the surface. In particular, it is found that positive spatial curvature inhibits finger tip-splitting. Hele-Shaw flow on weakly negative, curved surfaces is briefly discussed.Comment: 26 pages, 4 figures, RevTex, accepted for publication in Phys. Rev.

Flow of foam through a convergent channel

International audienceWe study experimentally the flow of a foam confined as a bubble monolayer between two plates through a convergent channel. We quantify the velocity, the distribution and orientation of plastic events, and the elastic stress, using image analysis. We use two different soap solutions: a sodium dodecyl sulfate (SDS) solution, with a negligible wall friction between the bubbles and the confining plates, and a mixture containing a fatty acid, giving a large wall friction. We show that for SDS solutions, the velocity profile obeys a self-similar form which results from the superposition of plastic events, and the elastic deformation is uniform. For the other solution, the velocity field differs and the elastic deformation increases towards the exit of the channel. We discuss and quantify the role of wall friction on the velocity profile, the elastic deformation, and the rate of plastic events

On the mechanisms governing gas penetration into a tokamak plasma during a massive gas injection

A new 1D radial fluid code, IMAGINE, is used to simulate the penetration of gas into a tokamak plasma during a massive gas injection (MGI). The main result is that the gas is in general strongly braked as it reaches the plasma, due to mechanisms related to charge exchange and (to a smaller extent) recombination. As a result, only a fraction of the gas penetrates into the plasma. Also, a shock wave is created in the gas which propagates away from the plasma, braking and compressing the incoming gas. Simulation results are quantitatively consistent, at least in terms of orders of magnitude, with experimental data for a D 2 MGI into a JET Ohmic plasma. Simulations of MGI into the background plasma surrounding a runaway electron beam show that if the background electron density is too high, the gas may not penetrate, suggesting a possible explanation for the recent results of Reux et al in JET (2015 Nucl. Fusion 55 093013)

Velocity-space sensitivity of the time-of-flight neutron spectrometer at JET

The velocity-space sensitivities of fast-ion diagnostics are often described by so-called weight functions. Recently, we formulated weight functions showing the velocity-space sensitivity of the often dominant beam-target part of neutron energy spectra. These weight functions for neutron emission spectrometry (NES) are independent of the particular NES diagnostic. Here we apply these NES weight functions to the time-of-flight spectrometer TOFOR at JET. By taking the instrumental response function of TOFOR into account, we calculate time-of-flight NES weight functions that enable us to directly determine the velocity-space sensitivity of a given part of a measured time-of-flight spectrum from TOFOR

Relationship of edge localized mode burst times with divertor flux loop signal phase in JET

A phase relationship is identified between sequential edge localized modes (ELMs) occurrence times in a set of H-mode tokamak plasmas to the voltage measured in full flux azimuthal loops in the divertor region. We focus on plasmas in the Joint European Torus where a steady H-mode is sustained over several seconds, during which ELMs are observed in the Be II emission at the divertor. The ELMs analysed arise from intrinsic ELMing, in that there is no deliberate intent to control the ELMing process by external means. We use ELM timings derived from the Be II signal to perform direct time domain analysis of the full flux loop VLD2 and VLD3 signals, which provide a high cadence global measurement proportional to the voltage induced by changes in poloidal magnetic flux. Specifically, we examine how the time interval between pairs of successive ELMs is linked to the time-evolving phase of the full flux loop signals. Each ELM produces a clear early pulse in the full flux loop signals, whose peak time is used to condition our analysis. The arrival time of the following ELM, relative to this pulse, is found to fall into one of two categories: (i) prompt ELMs, which are directly paced by the initial response seen in the flux loop signals; and (ii) all other ELMs, which occur after the initial response of the full flux loop signals has decayed in amplitude. The times at which ELMs in category (ii) occur, relative to the first ELM of the pair, are clustered at times when the instantaneous phase of the full flux loop signal is close to its value at the time of the first ELM

The microrheology of wet forms

The Kelvin cell is the only known topology for stable, perfectly ordered, dry foams. During topological transitions (T1s) associated with large elastic-plastic deformations, these cells switch neighbors and some faces gain or lose two sides, but the resulting bubbles with different shape are still Kelvin cells. The bubbles in a stable, perfectly ordered. wet foam are not limited to one topology (or even the two described here). The topological transitions considered here result in gain or loss of two dry films per bubble. The transition from Kelvin to RD topology is triggered by films shrinking in area, as in the dry case. However, the reverse transition from RD to Kelvin topology involves a different mechanism--opposite interfaces of an eight-way vertex touch and a new film grows from the point of contact as the foam is compressed. Microrheological analysis based on 2D models of foam structure has been useful preparation for 3D, despite obvious differences between 2D and 3D. Linear elastic behavior is anisotropic for perfectly ordered 3D foams--nonlinear elastic behavior is isotropic for 2D foams with polydisperse hexagonal structure. The shear moduli of a wet Kelvin foam decrease with increasing {phi}--the shear modulus of a wet 2D foam (with three-way Plateau borders) does not depend on {phi} at all. The effective isotropic shear moduli G of perfectly ordered wet foams tend to decrease with increasing {phi} but do not exhibit linear dependence, which may stem from the disorder of real systems

The structure of foam cells: isotropic plateau polyhedra

A mean-field theory for the geometry and diffusive growth rate of soap bubbles in dry 3D foams is presented. Idealized foam cells called isotropic Plateau polyhedra (IPPs), with F identical spherical-cap faces, are introduced. The geometric properties (e.g., surface area S, curvature R, edge length L, volume V) and growth rate of the cells are obtained as analytical functions of F, the sole variable. IPPs accurately represent average foam bubble geometry for arbitrary F ≥ 4, even though they are only constructible for F = 4,6,12. While R/V1/3, L/V1/3 and exhibit F1/2 behavior, the specific surface area S/V2/3 is virtually independent of F. The results are contrasted with those for convex isotropic polyhedra with flat faces