Atomistic Description of Interdroplet Ice-Bridge Formation during Condensation Frosting.

The propagation of frost in an assembly of supercooled dew droplets takes place by the formation of ice protrusions that bridge ice particles and still-liquid droplets. In this work, we develop a Kinetic Monte Carlo (KMC) model to study the formation kinetics of the ice protrusions. The KMC simulations reproduce well the experimental results reported in the literature. The elongation speed of the ice protrusions does not depend on the droplet size but increases when the interdroplet distance decreases, the temperature increases, or the substrate wettability increases. While 2D diffusion of the water molecules on the substrate surface is sufficient to explain the process kinetics, high 3D (vapor) water-molecule concentration can lead to the development of 3D lateral branches on the ice protrusions. A 1D analytical model based on the water-molecule concentration gradient between a droplet and a nearby ice particle reproduces well the simulation results and highlights the relation between the protrusion elongation kinetics and parameters like the interdroplet distance, the water diffusivity, and the concentration gradient. The bridge-formation time has a quadratic dependence on the droplet-ice distance. Comparing the simulations, the analytical model, and the experimental results of the literature, we conclude that the propagation of frost on a flat substrate in an assembly of supercooled dew droplets with interdroplet spacing larger than about 1 µm is limited by water-molecule diffusivity.

Three-dimensional liquid foam flow through a hopper resolved by fast X-ray microtomography.

We probe the complex rheological behaviour of liquid foams flowing through a conical constriction. With fast X-ray tomographic microscopy we measure in situ the displacement and deformation of up to fifty thousand bubbles at any single time instance while varying systematically the foam liquid fraction, the bubble size and the flow direction - convergent vs. divergent. The large statistics and high spatio-temporal resolution allows to observe and quantify the deviations from a purely viscous flow. We indeed reveal an asymmetry between the convergent and divergent flows associated to the emergence of elastic stresses in the latter case, and enhanced as the liquid fraction is reduced. Such effect is related to the reorientation of the deformed bubbles flowing out of the constriction, from a prolate to an oblate shape in average, while they pass through the hopper waist.

Dynamics of Frost Propagation on Breath Figures.

We investigate the process of condensation frosting on flat surfaces using thermal imaging microscopy. This method is particularly well-suited to characterize the frosting of polydisperse assemblies of dew droplets, also called breath figures, that transform into ice droplets by the propagation of frost fronts. The front propagation speed is found to be a nonmonotonous function of the characteristic droplet size of the breath figure. In our experimental conditions, the propagation speed is maximum around 70 µm s-1 for a characteristic droplet radius of around 300 µm. We mainly show that the frost propagation speed is governed by the competition between two characteristic time scales. The first one is the freezing time of individual droplets, and the other one is the formation time of interdroplet ice bridges that grow from frozen to liquid droplets. In addition, the experiments reveal that the mean ice bridge speed is constant regardless of the characteristic radius of the liquid droplets in the breath figure. A theoretical mean-field analysis without any adjustable parameters recovers all of the features of the front propagation observed in experiments.

Collapse of a hemicatenoid bounded by a solid wall: instability and dynamics driven by surface Plateau border friction.

The collapse of a catenoidal soap film when the rings supporting it are moved beyond a critical separation is a classic problem in interface motion in which there is a balance between surface tension and the inertia of the surrounding air, with film viscosity playing only a minor role. Recently [Goldstein et al., Phys. Rev. E, 2021, 104, 035105], we introduced a variant of this problem in which the catenoid is bisected by a glass plate located in a plane of symmetry perpendicular to the rings, producing two identical hemicatenoids, each with a surface Plateau border (SPB) on the glass plate. Beyond the critical ring separation, the hemicatenoids collapse in a manner qualitatively similar to the bulk problem, but their motion is governed by the frictional forces arising from viscous dissipation in the SPBs. We present numerical studies of a model that includes classical laws in which the frictional force fv for SPB motion on wet surfaces is of the form fv ∼ Can, where Ca is the capillary number. Our experimental data on the temporal evolution of this process confirms the expected value n = 2/3 for mobile surfactants and stress-free interfaces. This study can help explain the fragmentation of bubbles inside very confined geometries such as porous materials or microfluidic devices.

Effect of gravity on the orientation and detachment of cubic particles adsorbed at soap film or liquid interfaces.

We investigate the interaction that occurs between a light solid cube falling under gravity and a horizontal soap film that is pinned to a circular ring. We observe in both experiments and quasi-static simulations that the final orientation of a cube that becomes entrapped by a soap film is strongly dependent on the Bond number. A cube is rotated by a soap film into one of three main orientations in a process that is driven by energy minimisation. The likelihood of observing each of these final orientations is shown to depend on the Bond number, and the most energetically favourable orientation depends on the terminal height reached by the cube. We also find a critical value for the Bond number, above which a cube is no longer supported by a soap film and detachment occurs, to be less than one.

Stress-Induced Dinoflagellate Bioluminescence at the Single Cell Level.

One of the characteristic features of many marine dinoflagellates is their bioluminescence, which lights up nighttime breaking waves or seawater sliced by a ship's prow. While the internal biochemistry of light production by these microorganisms is well established, the manner by which fluid shear or mechanical forces trigger bioluminescence is still poorly understood. We report controlled measurements of the relation between mechanical stress and light production at the single cell level, using high-speed imaging of micropipette-held cells of the marine dinoflagellate Pyrocystis lunula subjected to localized fluid flows or direct indentation. We find a viscoelastic response in which light intensity depends on both the amplitude and rate of deformation, consistent with the action of stretch-activated ion channels. A phenomenological model captures the experimental observations.

Ferrofluid Leidenfrost droplets.

We experimentally investigate the behavior of ferrofluid Leidenfrost droplets subject to a static magnetic field gradient. The droplets are deposited on a hot substrate and trapped over the vertical axis of a permanent magnet placed at a distance d above the substrate. Several effects are observed. Firstly, the droplet evaporation rate is strongly influenced by the distance d. Secondly, the droplet takes off from the substrate when its radius decreases to a critical value. The introduction of an effective gravity, which accounts for the magnetic force, allows a successful description of these effects. Finally, we observe an instability for which the droplet starts bouncing with irregular amplitudes. This behavior is qualitatively interpreted by introducing the synchronization of the free fall time between successive bounces with the period of the fundamental vibration mode of the droplet.

Experimental Characterization of Droplet Adhesion: The Ejection Test Method (ETM) Applied to Surfaces with Various Hydrophobicity.

We study the wetting and the adhesive behavior of substrates made by electropolymerization of copolymers of pyrene substituted with fluoroalkyl and adamantyl groups. The hydrophobicity and water adhesion properties can be tuned by the molar percentage (mol %) of each pyrene monomer so that the substrate properties can vary from superhydrophobic to parahydrophobic, with respectively low and high water adhesion. The ejection test method (ETM) is proposed as an original tool to discriminate and characterize such substrates. Using a catapult-like apparatus, a droplet initially at rest on the surface is subject to a large acceleration and is subsequently ejected. Depending on the surface properties and initial catapult acceleration, the ejection is more or less efficient and occurs with or without fragmentation of the droplet. The ETM is shown to be a complementary test to the lateral adhesion and hysteresis classical measurements. This work is of importance for the understanding of adhesion phenomena on various surfaces and for a better quantitative characterization of their adhesive properties.

Melde's Experiment on a Vibrating Liquid Foam Microchannel.

We subject a single Plateau border channel to a transverse harmonic excitation, in an experiment reminiscent of the historical one by Melde on vibrating strings, to study foam stability and wave properties. At low driving amplitudes, the liquid string exhibits regular oscillations. At large ones, a nonlinear regime appears and the acoustic radiation splits the channel into two zones of different cross section area, vibration amplitude, and phase difference with the neighboring soap films. The channel experiences an inertial dilatancy that is accounted for by a new Bernoulli-like relation.

Superpropulsion of Droplets and Soft Elastic Solids.

We investigate the behavior of droplets and soft elastic objects propelled with a catapult. Experiments show that the ejection velocity depends on both the projectile deformation and the catapult acceleration dynamics. With a subtle matching given by a peculiar value of the projectile/catapult frequency ratio, a 250% kinetic energy gain is obtained as compared to the propulsion of a rigid projectile with the same engine. This superpropulsion has strong potentialities: actuation of droplets, sorting of objects according to their elastic properties, and energy saving for propulsion engines.

Quantifying order in breath figure patterns through Voronoi entropy.

In this study, we simulate breath figures that are evolving two-dimensional assemblies of droplets on a substrate. We focus on the Voronoi/Shannon entropy of these figures, which quantifies the order related to the coordination number of droplets. We show that the Voronoi entropy of the complete breath figure pattern converges to a value that is the one of a randomly distributed point system. Conversely, the subset containing exclusively large droplets of the breath figure exhibits significantly lower entropy than that obtained for all droplets. Using molecular dynamics simulations, we show that coalescence events in breath figures induce the same Voronoi entropy as that caused by repulsive interactions in a bidimensional atomic system.

Scaling the tail beat frequency and swimming speed in underwater undulatory swimming.

Undulatory swimming is the predominant form of locomotion in aquatic vertebrates. A myriad of animals of different species and sizes oscillate their bodies to propel themselves in aquatic environments with swimming speed scaling as the product of the animal length by the oscillation frequency. Although frequency tuning is the primary means by which a swimmer selects its speed, there is no consensus on the mechanisms involved. In this article, we propose scaling laws for undulatory swimmers that relate oscillation frequency to length by taking into account both the biological characteristics of the muscles and the interaction of the moving swimmer with its environment. Results are supported by an extensive literature review including approximately 1200 individuals of different species, sizes and swimming environments. We highlight a crossover in size around 0.5-1 m. Below this value, the frequency can be tuned between 2-20 Hz due to biological constraints and the interplay between slow and fast muscles. Above this value, the fluid-swimmer interaction must be taken into account and the frequency is inversely proportional to the length of the animal. This approach predicts a maximum swimming speed around 5-10 m.s-1 for large swimmers, consistent with the threshold to prevent bubble cavitation.

Impact of composite soft-rigid elastic projectiles: A numerical study.

We study by numerical simulation the impact of a one-dimensional composite projectile, composed of two superposed homogeneous parts, on an infinitely rigid and massive wall. The coefficient of restitution and the contact time are systematically measured as functions of the contrasts of mass and stiffness between the two parts. For purely elastic parts, these quantities show complex trends associated with different dynamics of the deformation waves propagating inside the projectile. A significant portion of the initial kinetic energy can be trapped in the deformation modes: the coefficient of restitution is lowest, about 0.2, when there is a strong stiffness contrast between the two parts and the stiff and soft parts are at the leading and trailing edges of the projectile respectively. In this case, we highlight the presence of multiple bounces, whose number increases as the proportion of the soft part increases. Finally, viscoelastic parts can be implemented in the same numerical framework to successfully recover the results obtained in real composite projectile impact experiments [D'Angelo et al., Phys. Rev. E 103, 053005 (2021)2470-004510.1103/PhysRevE.103.053005].

Reinforcement learning approach to control an inverted pendulum: A general framework for educational purposes.

Machine learning is often cited as a new paradigm in control theory, but is also often viewed as empirical and less intuitive for students than classical model-based methods. This is particularly the case for reinforcement learning, an approach that does not require any mathematical model to drive a system inside an unknown environment. This lack of intuition can be an obstacle to design experiments and implement this approach. Reversely there is a need to gain experience and intuition from experiments. In this article, we propose a general framework to reproduce successful experiments and simulations based on the inverted pendulum, a classic problem often used as a benchmark to evaluate control strategies. Two algorithms (basic Q-Learning and Deep Q-Networks (DQN)) are introduced, both in experiments and in simulation with a virtual environment, to give a comprehensive understanding of the approach and discuss its implementation on real systems. In experiments, we show that learning over a few hours is enough to control the pendulum with high accuracy. Simulations provide insights about the effect of each physical parameter and tests the feasibility and robustness of the approach.

Use of compliant actuators for throwing rigid projectiles.

Muscles and tendons, actuators in robotics, and various sports implements are examples that exploit elasticity to accelerate objects. Tuning the mechanical properties of elastic elements connecting objects can greatly enhance the transfer of mechanical energy between the objects. Here, we study experimentally the throw of rigid projectiles by an actuator, which has a soft elastic element added to the distal end. We vary the thickness of the elastic layer and suggest a simple mass-spring chain model to find the properties of the elastic layer, which will maximize the energy transfer from the actuator to the projectile. The insertion of a soft layer, impedance matched to the ejection frequency of the projectile mass, can increase the throwing efficiency by over 400%. Finally, we identify that very thick and very soft compliant layers could potentially lead to high efficiency and flexibility simultaneously.

Impact dynamics of composite elastorigid projectiles onto solid surfaces.

We investigate the impact of composite objects. They consist of a soft layer on top of a rigid part with a hemispherical impacting end. The coefficient of restitution (e) of such objects is studied systematically as a function of the mass ratio and of the nature of the materials. For rather elastic materials, the coefficient of restitution is a nonmonotonic function of the mass ratio and exhibits important variations. The dynamics of the impact can be characterized by several bounces depending on the ratios between the four timescales at play. These include the duration of contact of the rigid part with the substrate and the time for the elastic waves to travel back and forth in the soft layer. In that sense, describing these projectiles requires one to take into account both the Hertzian theory of contact and the elastic waves described by Saint-Venant's approach.

Geometry of catenoidal soap film collapse induced by boundary deformation.

Experimental and theoretical work reported here on the collapse of catenoidal soap films of various viscosities reveal the existence of a robust geometric feature that appears not to have been analyzed previously; prior to the ultimate pinchoff event on the central axis, which is associated with the formation of a well-studied local double-cone structure folded back on itself, the film transiently consists of two acute-angle cones connected to the supporting rings, joined by a central quasicylindrical region. As the cylindrical region becomes unstable and pinches, the opening angle of those cones is found to be universal, independent of film viscosity. Moreover, that same opening angle at pinching is found when the transition occurs in a hemicatenoid bounded by a surface. The approach to the conical structure is found to obey classical Keller-Miksis scaling of the minimum radius as a function of time, down to very small but finite radii. While there is a large body of work on the detailed structure of the singularities associated with ultimate pinchoff events, these large-scale features have not been addressed. Here we study these geometrical aspects of film collapse by several distinct approaches, including a systematic analysis of the linear and weakly nonlinear dynamics in the neighborhood of the saddle node bifurcation leading to collapse, both within mean curvature flow and the physically realistic Euler flow associated with the incompressible dynamics of the surrounding air. These analyses are used to show how much of the geometry of collapsing catenoids is accurately captured by a few active modes triggered by boundary deformation. A separate analysis based on a mathematical sequence of shapes progressing from the critical catenoid towards the Goldschmidt solution is shown to predict accurately the cone angle at pinching. We suggest that the approach to the conical structures can be viewed as passage close to an unstable fixed point of conical similarity solutions. The overall analysis provides the basis for the systematic study of more complex problems of surface instabilities triggered by deformations of the supporting boundaries.

Bubble dynamics inside an outgassing hydrogel confined in a Hele-Shaw cell.

We report an experimental study of bubble dynamics in a non-Newtonian fluid subjected to a pressure decrease. The fluid is a hydrogel, composed of water and a synthetic clay, prepared and sandwiched between two glass plates in a Hele-Shaw geometry. The rheological properties of the material can be tuned by the clay concentration. As the imposed pressure decreases, the gas initially dissolved in the hydrogel triggers bubble formation. Different stages of the process are observed: bubble nucleation, growth, interaction, and creation of domains by bubble contact or coalescence. Initially bubble behave independently. They are trapped and advected by the mean deformation of the hydrogel, and the bubble growth is mainly driven by the diffusion of the dissolved gas through the hydrogel and its outgassing at the reactive-advected hydrogel-bubble interface. In this regime, the rheology of the fluid does not play a significant role on the bubble growth. A model is proposed and gives a simple scaling that relates the bubble growth rate and the imposed pressure. Carbon dioxide is shown to be the gas at play, and the hydrogel is degassing at the millimeter scale as a water solution does at a smaller scale. Later, bubbles are not independent anymore. The growth rate decreases, and the morphology becomes more anisotropic as bubbles interact because they are separated by a distance smaller than the individual stress field extension. Our measurements show that the interaction distance scales with the bubbles' size.

Two-dimensional flow of foam around an obstacle: force measurements.

A Stokes experiment for foams is proposed. It consists of a two-dimensional flow of a foam, confined between a water subphase and a top plate, around a fixed circular obstacle. We present systematic measurements of the drag exerted by the flowing foam on the obstacle versus various separately controlled parameters: flow rate, bubble volume, bulk viscosity, obstacle size, shape, and boundary conditions. We separate the drag into two contributions: an elastic one (yield drag) at vanishing flow rate and a fluid one (viscous coefficient) increasing with flow rate. We quantify the influence of each control parameter on the drag. The results exhibit in particular a power-law dependence of the drag as a function of the bulk viscosity and the flow rate with two different exponents. Moreover, we show that the drag decreases with bubble size and increases proportionally to the obstacle size. We quantify the effect of shape through a dimensional drag coefficient, and we show that the effect of boundary conditions is small.

Drop coalescence and liquid flow in a single Plateau border.

We report a comprehensive study of the flow of liquid triggered by injecting a droplet into a liquid foam microchannel, also called a Plateau border. This drop-injected experiment reveals an intricate dynamics for the liquid redistribution, with two contrasting regimes observed, ruled either by inertia or viscosity. We devoted a previous study [A. Cohen et al., Phys. Rev. Lett. 112, 218303 (2014)] to the inertial imbibition regime, unexpected at such small length scales. Here we report other features of interest of the drop-injected experiment, related to the coalescence of the droplet with the liquid microchannel, to both the inertial and viscous regimes, and to the occurrence of liquid flow through the soap films as well as effects of the interfacial rheology. The transition between the two regimes is investigated and qualitatively accounted for. The relevance of our results to liquid foam drainage is tackled by considering the flow of liquid at the nodes of the network of interconnected microchannels. Extensions of our study to liquid foams are discussed.