Interfacial forces and Marangoni flow on a nematic drop retracting in an isotropic fluid

Geometry-structure models for liquid crystal interfaces, drops and membranes: wrinkling, shape selection and dissipative shape evolution

We review our recent contributions to anisotropic soft matter models for liquid crystal interfaces, drops and membranes, emphasizing validations with experimental and biological data, and with related theory and simulation literature. The presentation aims to illustrate and characterize the rich output and future opportunities of using a methodology based on the liquid crystal-membrane shape equation applied to static and dynamic pattern formation phenomena. The geometry of static and kinetic shapes is usually described with dimensional curvatures that co-mingle shape and curvedness. In this review, we systematically show how the application of a novel decoupled shape-curvedness framework to practical and ubiquitous soft matter phenomena, such as the shape of drops and tactoids and bending of evolving membranes, leads to deeper quantitative insights than when using traditional dimensional mean and Gaussian curvatures. The review focuses only on (1) statics of wrinkling and shape selection in liquid crystal interfaces and membranes; (2) kinetics and dissipative dynamics of shape evolution in membranes; and (3) computational methods for shape selection and shape evolution; due to various limitations other important topics are excluded. Finally, the outlook follows a similar structure. The main results include: (1) single and multiple wavelength corrugations in liquid crystal interfaces appear naturally in the presence of surface splay and bend orientation distortions with scaling laws governed by ratios of anchoring-to-isotropic tension energy; adding membrane elasticity to liquid crystal anchoring generates multiple scales wrinkling as in tulips; drops of liquid crystals encapsulates in membranes can adopt, according to the ratios of anchoring/tension/bending, families of shapes as multilobal, tactoidal, and serrated as observed in biological cells. (2) Mapping the liquid crystal director to a membrane unit normal. The dissipative shape evolution model with irreversible thermodynamics for flows dominated by bending rates, yields new insights. The model explains the kinetic stability of cylinders, while spheres and saddles are attractors. The model also adds to the evolving understanding of outer hair cells in the inner ear. (3) Computational soft matter geometry includes solving shape equations, trajectories on energy and orientation landscapes, and shape-curvedness evolutions on entropy production landscape with efficient numerical methods and adaptive approaches.

Phase-field model for elastocapillary flows of liquid crystals

We propose a phase-field model to study interfacial flows of nematic liquid crystals that couple the capillary forces on the interface with the elastic stresses in the nematic phase. The theoretical model has two key ingredients: A tensor order parameter that provides a consistent description of the molecular and distortional elasticity, and a phase-field formalism that accurately represents the interfacial tension and the nematic anchoring stress by approximating a sharp-interface limit. Using this model, we carry out finite-element simulations of drop retraction in a surrounding fluid, with either component being nematic. The results are summarized by eight representative steady-state solutions in planar and axisymmetric geometries, each featuring a distinct configuration for the drop and the defects. The dynamics is dominated by the competition between the interfacial tension and the distortional elasticity in the nematic phase, mediated by the anchoring condition on the drop surface. As consequences of this competition, the steady-state drop deformation and the clearance between the defects and the drop surface both depend linearly on the elastocapillary number.

Mechanogeometry of nanowrinkling in cholesteric liquid crystal surfaces

Biological plywoods are multifunctional fibrous composites materials, ubiquitous in nature. The chiral fibrous organization is found in chitin (insects), cellulosics (plants), and collagen I (cornea and bone of mammals) and is a solid analog of that of cholesteric liquid crystals. The surface and interfaces of plywoods are distinguished by hierarchical topographies and nanowrinkling. In this paper, we present a theory to model the emergence of these surfaces and interfaces using liquid crystal-based shape equations that directly connect material properties with geometric wrinkling. The model applies to liquid crystal precursors of the plywood solid analoges. We focus on wrinkling geometry, wrinkling mechanics, and the mechanogeometry relationships that underlie multifunctionality ubiquitous in biological surfaces. Scaling wrinkling laws that connect mechanical pressures and stresses to folding and bending are formulated and quantified. A synthesis of the connections between mechanics and geometry is achieved using the topology of stress curves and curvature of the wrinkles. Taken together the results show that anchoring is a versatile surface morphing mechanism with a rich surface bending stress field, two ingredients behind many potential multifunctionalities.

On the migration of a droplet on an incline

A liquid droplet placed on a nonuniformly heated solid surface will migrate from a high temperature region to a low temperature region. The present study reports the results of an experimental investigation on the migration behavior of mineral oil droplets subjected to a thermal gradient on an inclined plane. A particular attention is paid to the relationship between the critical inclination angle and thermal gradients. It is shown that there exists a critical inclination angle at which the droplet migration is halted. This critical inclination angle can be readily predicted using analytical expressions derived in this paper. This study puts forward the understanding of the interface phenomenon of thermocapillary migration on an incline. The knowledge of the critical inclination is important in applications where the migration on an incline must be obstructed to retain adequate lubrication in the desired location.

Modeling of coupled motion and growth interaction of equiaxed dendritic crystals in a binary alloy during solidification

Motion of growing dendrites is a common phenomenon during solidification but often neglected in numerical simulations because of the complicate underlying multiphysics. Here a phase-field model incorporating dendrite-melt two-phase flow is proposed for simulating the dynamically interacted process. The proposed model circumvents complexity to resolve dendritic growth, natural convection and solid motion simultaneously. Simulations are performed for single and multiple dendritic growth of an Al-based alloy in a gravity environment. Computing results of an isolated dendrite settling down in the convective supersaturated melt shows that solid motion is able to overwhelm solutal convection and causes a rather different growth morphology from the stationary dendrite that considers natural convection alone. The simulated tip growth dynamics are correlated with a modified boundary layer model in the presence of melt flow, which well accounts for the variation of tip velocity with flow direction. Polycrystalline simulations reveal that the motion of dendrites accelerates the occurrence of growth impingement which causes the behaviors of multiple dendrites are distinct from that of single dendrite, including growth dynamics, morphology evolution and movement path. These polycrystalline simulations provide a primary understanding of the sedimentation of crystals and resulting chemical homogeneity in industrial ingots.

Order parameter description of electrochemical-hydrodynamic interactions in nanochannels

A novel phase-field model is developed for the quantitative modeling of the complex electrochemical-hydrodynamic interactions in narrow fluidic confinements. Through an order parameter variation, this model captures the underlying excluded volume effects, solvation interactions, and preferential polarizabilities in a self-consistent fashion, without resorting to computationally prohibitive molecular dynamics simulations. Agreement with molecular dynamics predictions is found to be quantitative.

Simulation of Neutrophil Deformation and Transport in Capillaries using Newtonian and Viscoelastic Drop Models

It is well known that neutrophils take much longer to traverse the pulmonary capillary bed than erythrocytes, and this is likely due to differences in the structure and rheology of the cells. In this study, we simulate the transit of a neutrophil in a capillary using a Newtonian drop model and a viscoelastic drop model. The cell membrane is represented by an interface with isotropic and constant tension, and the cell motion and deformation are described by a phase-field method. The governing equations are solved using finite elements in an axisymmetric geometry, and the thin interfaces are resolved by mesh adaptivity. With a fixed pressure drop, the entry of a cell into a capillary consists of several stages in which the flow rate varies in distinct manners. The entrance time is consistent with experimental measurements. It decreases with the pressure drop, increases with the cell viscosity and generally decreases with the relaxation time of a viscoelastic cytoplasm. The capillary geometry has a strong effect on the entry and transit of a neutrophil. The entrance time increases sharply when the capillary diameter decreases or when the capillary is constricted by a pinch.

The influence of electric field-induced orientation on the retraction of a liquid crystal droplet immersed in a flexible polymer matrix

We measured the apparent interfacial tension between a liquid crystal and a flexible polymer by deformed droplet retraction method. An external electric field is applied to change the director orientation in liquid crystal droplet. The deformation and recovery of a single liquid crystal droplet dispersed in a polydimethylsiloxane (PDMS) matrix were realized by a transient shear flow and observed by polarized optical microscope. In order to control the director orientation in LC droplet, the electric field is applied perpendicular and parallel to the flow field, respectively. The different orientation induced by electric field in liquid crystal droplet has different behavior during droplet retraction and affect the apparent interfacial tension between liquid crystal and flexible polymer.