Lattice Boltzmann study of chemically-driven self-propelled droplets

We numerically study the behavior of self-propelled liquid droplets whose motion is triggered by a Marangoni-like flow. This latter is generated by variations of surfactant concentration which affect the droplet surface tension promoting its motion. In the present paper a model for droplets with a third amphiphilic component is adopted. The dynamics is described by Navier-Stokes and convection-diffusion equations, solved by lattice Boltzmann method coupled with finite-difference schemes. We focus on two cases. First the study of self-propulsion of an isolated droplet is carried on and, then, the interaction of two self-propelled droplets is investigated. In both cases, when the surfactant migrates towards the interface, a quadrupolar vortex of the velocity field forms inside the droplet and causes the motion. A weaker dipolar field emerges instead when the surfactant is mainly diluted in the bulk. The dynamics of two interacting droplets is more complex and strongly depends on their reciprocal distance. If, in a head-on collision, droplets are close enough, the velocity field initially attracts them until a motionless steady state is achieved. If the droplets are vertically shifted, the hydrodynamic field leads to an initial reciprocal attraction followed by a scattering along opposite directions. This hydrodynamic interaction acts on a separation of some droplet radii otherwise it becomes negligible and droplets motion is only driven by Marangoni effect. Finally, if one of the droplets is passive, this latter is generally advected by the fluid flow generated by the active one.Comment: 14 pages, 9 figures. In press on EPJ

Hybrid lattice Boltzmann model for binary fluid mixtures

A hybrid lattice Boltzmann method (LBM) for binary mixtures based on the free-energy approach is proposed. Non-ideal terms of the pressure tensor are included as a body force in the LBM kinetic equations, used to simulate the continuity and Navier-Stokes equations. The convection-diffusion equation is studied by finite difference methods. Differential operators are discretized in order to reduce the magnitude of spurious velocities. The algorithm has been shown to be stable and reproducing the correct equilibrium behavior in simple test configurations and to be Galilean invariant. Spurious velocities can be reduced of about an order of magnitude with respect to standard discretization procedure.Comment: Final version, to appear in Phys. Rev.

Switching dynamics in cholesteric blue phases

Blue phases are networks of disclination lines, which occur in cholesteric liquid crystals near the transition to the isotropic phase. They have recently been used for the new generation of fast switching liquid crystal displays. Here we study numerically the steady states and switching hydrodynamics of blue phase I (BPI) and blue phase II (BPII) cells subjected to an electric field. When the field is on, there are three regimes: for very weak fields (and strong anchoring at the boundaries) the blue phases are almost unaffected, for intermediate fields the disclinations twist (for BPI) and unzip (for BPII), whereas for very large voltages the network dissolves in the bulk of the cell. Interestingly, we find that a BPII cell can recover its original structure when the field is switched off, whereas a BPI cell is found to be trapped more easily into metastable configurations. The kinetic pathways followed during switching on and off entails dramatic reorganisation of the disclination networks. We also discuss the effect of changing the director field anchoring at the boundary planes and of varying the direction of the applied field.Comment: 17 pages, 11 figure

Lattice Boltzmann modeling of cholesteric liquid crystal droplets under an oscillatory electric field

We numerically study the dynamics of quasi-two dimensional cholesteric liquid crystal droplets in the presence of a time-dependent electric field, rotating at constant angular velocity. A surfactant sitting at droplet interface is also introduced to prevent droplet coalescence. The dynamics is modeled following a hybrid numerical approach, where a standard lattice Boltzmann technique solves the Navier-Stokes equation and a finite difference scheme integrates the evolution equations of liquid crystal and surfactant. Our results show that, once the field is turned on, the liquid crystal rotates coherently triggering a concurrent orbital motion of both droplets around each other, an effect due to the momentum transfer to the surrounding fluid. In addition the topological defects, resulting from the conflict orientation of the liquid crystal within the drops, exhibit a chaotic-like motion in cholesterics with a high pitch, in contrast with a regular one occurring along circular trajectories observed in nematics drops. Such behavior is found to depend on magnitude and frequency of the applied field as well as on the anchoring of the liquid crystal at the droplet interface. These findings are quantitatively evaluated by measuring the angular velocity of fluid and drops for various frequencies of the applied field.Comment: 16 pages, 13 figures. Accepted on Communications in Computational Physic

Bistable defect structures in blue phase devices

Blue phases (BPs) are liquid crystals made up by networks of defects, or disclination lines. While existing phase diagrams show a striking variety of competing metastable topologies for these networks, very little is known as to how to kinetically reach a target structure, or how to switch from one to the other, which is of paramount importance for devices. We theoretically identify two confined blue phase I systems in which by applying an appropriate series of electric field it is possible to select one of two bistable defect patterns. Our results may be used to realise new generation and fast switching energy-saving bistable devices in ultrathin surface treated BPI wafers.Comment: 4 pages, 3 figures. Accepted for publication in Phys. Rev. Let

Playing with Casimir in the vacuum sandbox

The Casimir effect continues to be a subject of discussion regarding its relationship, or the lack of it, with the vacuum energy of fluctuating quantum fields. In this note, we propose a Gedankenexperiment considering an imaginary process similar to a vacuum fluctuation in a typical static Casimir set up. The thought experiment leads to intriguing conclusions regarding the minimum distance between the plates when approaching the Planck scale. More specifically, it is found that distance between the plates cannot reach a value below $(L/L_P)^{2/3}$ Planck lengths, being $L_P$ the Planck length and $L$ the typical lateral extension of the plates. Additional findings allow the conclusion that the approach between the two plates towards this minimum separation distance is asymptotic

Phase separation of binary fluids with dynamic temperature

Phase separation of binary fluids quenched by contact with cold external walls is considered. Navier-Stokes, convection-diffusion, and energy equations are solved by lattice Boltzmann method coupled with finite-difference schemes. At high viscosity, different morphologies are observed by varying the thermal diffusivity. In the range of thermal diffusivities with domains growing parallel to the walls, temperature and phase separation fronts propagate towards the inner of the system with power-law behavior. At low viscosity hydrodynamics favors rounded shapes, and complex patterns with different lengthscales appear. Off-symmetrical systems behave similarly but with more ordered configurations.Comment: Accepted for publication in Phys. Rev. E, 11 figures, best quality figures available on reques

Soft channel formation and symmetry breaking in exotic active emulsions

We use computer simulations to study the morphology and rheological properties of a bidimensional emulsion resulting from a mixture of a passive isotropic fluid and an active contractile polar gel, in the presence of a surfactant that favours the emulsification of the two phases. By varying the intensity of the contractile activity and of an externally imposed shear flow, we find three possible morphologies. For low shear rates, a simple lamellar state is obtained. For intermediate activity and shear rate, an asymmetric state emerges, which is characterized by shear and concentration banding at the polar/isotropic interface. A further increment in the active forcing leads to the self-assembly of a soft channel where an isotropic fluid flows between two layers of active material. We characterize the stability of this state by performing a dynamical test varying the intensity of the active forcing and shear rate. Finally, we address the rheological properties of the system by measuring the effective shear viscosity, finding that this increases as active forcing is increased—so that the fluid thickens with activity