The minimal perimeter for N confined deformable bubbles of equal area

Candidates to the least perimeter partition of various polygonal shapes into N planar connected equal-area regions are calculated for N � 42, compared to partitions of the disc, and discussed in the context of the energetic groundstate of a two-dimensional monodisperse foam. The total perimeter and the number of peripheral regions are presented, and the patterns classified according to the number and position of the topological defects, that is non-hexagonal regions (bubbles). The optimal partitions of an equilateral triangle are found to follow a pattern based on the position of no more than one defect pair, and this pattern is repeated for many of the candidate partitions of a hexagon. Partitions of a square and a pentagon show greater disorder. Candidates to the least perimeter partition of the surface of the sphere into N connected equal-area regions are also calculated. For small N these can be related to simple polyhedra and for N � 14 they consist of 12 pentagons and N −12 hexagons.

Partitions of Minimal Length on Manifolds

We study partitions on three dimensional manifolds which minimize the total geodesic perimeter. We propose a relaxed framework based on a $\Gamma$-convergence result and we show some numerical results. We compare our results to those already present in the literature in the case of the sphere. For general surfaces we provide an optimization algorithm on meshes which can give a good approximation of the optimal cost, starting from the results obtained using the relaxed formulation

Curvature driven motion of a bubble in a toroidal Hele-Shaw cell

We investigate the equilibrium properties of a single area-minimizing bubble trapped between two narrowly separated parallel curved plates. We begin with the case of a bubble trapped between concentric spherical plates. We develop a model which shows that the surface energy of the bubble is lower when confined between spherical plates than between flat plates. We confirm our findings by comparing against Surface Evolver simulations. We then derive a simple model for a bubble between arbitrarily curved parallel plates. The energy is found to be higher when the local Gaussian curvature of the plates is negative and lower when the curvature is positive. To check the validity of the model, we consider a bubble trapped between concentric tori. In the toroidal case, we find that the sensitivity of the bubble’s energy to the local curvature acts as a geometric potential capable of driving bubbles from regions with negative to positive curvature

Boundary elements method for microfluidic two-phase flows in shallow channels

In the following work we apply the boundary element method to two-phase flows in shallow microchannels, where one phase is dispersed and does not wet the channel walls. These kinds of flows are often encountered in microfluidic Lab-on-a-Chip devices and characterized by low Reynolds and low capillary numbers. Assuming that these channels are homogeneous in height and have a large aspect ratio, we use depth-averaged equations to describe these two-phase flows using the Brinkman equation, which constitutes a refinement of Darcy's law. These partial differential equations are discretized and solved numerically using the boundary element method, where a stabilization scheme is applied to the surface tension terms, allowing for a less restrictive time step at low capillary numbers. The convergence of the numerical algorithm is checked against a static analytical solution and on a dynamic test case. Finally the algorithm is applied to the non-linear development of the Saffman-Taylor instability and compared to experimental studies of droplet deformation in expanding flows.Comment: accepted for publication, Computers and Fluids 201

Statistical mechanics of two-dimensional shuffled foams: Geometry-topology correlation in small or large disorder limits

Bubble monolayers are model systems for experiments and simulations of two-dimensional packing problems of deformable objects. We explore the relation between the distributions of the number of bubble sides (topology) and the bubble areas (geometry) in the low liquid fraction limit. We use a statistical model [M. Durand, Europhys. Lett. 90, 60002 (2010)] which takes into account Plateau laws. We predict the correlation between geometrical disorder (bubble size dispersity) and topological disorder (width of bubble side number distribution) over an extended range of bubble size dispersities. Extensive data sets arising from shuffled foam experiments, surface evolver simulations, and cellular Potts model simulations all collapse surprisingly well and coincide with the model predictions, even at extremely high size dispersity. At moderate size dispersity, we recover our earlier approximate predictions [M. Durand, J. Kafer, C. Quilliet, S. Cox, S. A. Talebi, and F. Graner, Phys. Rev. Lett. 107, 168304 (2011)]. At extremely low dispersity, when approaching the perfectly regular honeycomb pattern, we study how both geometrical and topological disorders vanish. We identify a crystallization mechanism and explore it quantitatively in the case of bidisperse foams. Due to the deformability of the bubbles, foams can crystallize over a larger range of size dispersities than hard disks. The model predicts that the crystallization transition occurs when the ratio of largest to smallest bubble radii is 1.4