We study localized modes in uniform one-dimensional chains of tightly packed and uniaxially compressed elastic beads in the presence of one or two light-mass impurities. For chains composed of beads of the same type, the intrinsic nonlinearity, which is caused by the Hertzian interaction of the beads, appears not to support localized, breathing modes. Consequently, the inclusion of light-mass impurities is crucial for their appearance. By analyzing the problem's linear limit, we identify the system's eigenfrequencies and the linear defect modes. Using continuation techniques, we find the solutions that bifurcate from their linear counterparts and study their linear stability in detail. We observe that the nonlinearity leads to a frequency dependence in the amplitude of the oscillations, a static mutual displacement of the parts of the chain separated by a defect, and for chains with two defects that are not in contact, it induces symmetry-breaking bifurcations.

Periodic pillarlike microstructures can be created from initially flat polymer films via the electrohydrodynamic instabilities. Those patterns, however, are metastable. Our experimental observations show that the average pillar size increases slowly after linear growth. Major coarsening events then take place over times several orders of magnitude longer than the linear growth time. For all fill ratios, a logarithmic time dependence of the average pillar size can be identified, i.e., proportional to ln t. Thicker films, however, have faster coarsening rates than thinner films. Linear stability analysis of the pseudosteady states reveals two major coarsening mechanisms, collision and Ostwald ripening, which can also be identified from experimental images. We then reduce the original partial differential equation (PDE) into a pair of ODEs, which govern the interaction between pillars due to the above two coarsening mechanisms. From this, a logarithm scaling law is obtained for both low and high fill ratios and the coarsening rate is slower for lower fill ratios, consistent with experimental observations. We also find that arrays with more uniform sizes tend to start coarsening later, but they coarsen faster than more "disperse" arrays, which could be possibly utilized in experiments for controlling the onset and speed of coarsening. The logarithm scaling in the electrohydrodynamic coarsening phenomenon, which differs from coarsening in spinodal decomposition and dewetting of thin liquid films, is due to the significant nonlinear effect of Maxwell stresses and geometric confinement on the disjoining pressure at both top and bottom electrodes.

Proper roughness design is important in realizing surfaces with fully tunable wetting properties. Engineering surface roughness boils down to an energy barrier optimization problem, in which the geometric features of roughness serve as the optimization parameters. Computations of energy barriers, separating admissible equilibrium wetting states on patterned surfaces, have been demonstrated utilizing fine-scale simulators (e.g., lattice-Boltzmann for mesoscale and molecular dynamics for microscale simulations), however with substantial computational requirements. Here, by solving an augmented Young–Laplace equation with a disjoining pressure term, we demonstrate accurate and efficient computations of equilibrium shapes of entire millimeter sized droplets on patterned surfaces. In particular, by adopting a natural parameterization of the Young–Laplace equation along the liquid/air and liquid/solid interfaces, the tedious implementation of the Young's contact angle boundary condition at multiple three phase contact lines is bypassed. We, thus, enable the computation of wetting transition energy barriers, separating the well-known Cassie–Baxter and Wenzel states, as well as intermediate states, but with negligible computational cost. We demonstrate the method's efficiency by computing the equilibrium of droplets on stripe-patterned surfaces, and compare the results with mesoscopic lattice Boltzmann simulations. Our computationally efficient continuum-level analysis can be readily applied to patterned surfaces with increased and unstructured geometric complexity, and straightforwardly coupled with shape optimizers towards the design of surfaces with desirable wetting behavior.

The complicated dynamics of the contact line of a moving droplet on a solid substrate often hamper the efficient modeling of microfluidic systems. In particular, the selection of the effective boundary conditions, specifying the contact line motion, is a controversial issue since the microscopic physics that gives rise to this displacement is still unknown. Here, a sharp interface, continuum-level, novel modeling approach, accounting for liquid/solid micro-scale interactions assembled in a disjoining pressure term, is presented. By following a unified conception (the model applies both to the liquid/solid and the liquid/ambient interfaces), the friction forces at the contact line, arXiv:1601.06540v1 [cond-mat.soft]

Aiming to illuminate mechanisms of wetting transitions on geometrically patterned surfaces induced by the electrowetting phenomenon, we present a novel modeling approach that goes beyond the limitations of the Lippmann equation and is even relieved from the implementation of the Young contact angle boundary condition. We employ the equations of the capillary electrohydrostatics augmented by a disjoining pressure term derived from an effective interface potential accounting for solid/liquid interactions. Proper parametrization of the liquid surface profile enables efficient simulation of multiple and reconfigurable three-phase contact lines (TPL) appearing when entire droplets undergo wetting transitions on patterned surfaces. The liquid/ambient and the liquid/solid interfaces are treated in a unified context tackling the assumption that the liquid profile is wedge-shaped at any three-phase contact line. In this way, electric field singularities are bypassed, allowing for accurate electric field and liquid surface profile computation, especially in the vicinity of TPLs. We found that the invariance of the microscopic contact angle in electrowetting systems is valid only for thick dielectrics, supporting published experiments. By applying our methodology to patterned dielectrics, we computed all admissible droplet equilibrium profiles, including Cassie-Baxter, Wenzel, and mixed wetting states. Mixed wetting states are computed for the first time in electrowetting systems, and their relative stability is presented in a clear and instructive way.

We present a model based on the lattice Boltzmann equation that is suitable for the simulation of dynamic wetting. The model is capable of exhibiting fundamental interfacial phenomena such as weak adsorption of fluid on the solid substrate and the presence of a thin surface film within which a disjoining pressure acts. Dynamics in this surface film, tightly coupled with hydrodynamics in the fluid bulk, determine macroscopic properties of primary interest: the hydrodynamic slip; the equilibrium contact angle; and the static and dynamic hysteresis of the contact angles. The pseudopotentials employed for fluid-solid interactions are composed of a repulsive core and an attractive tail that can be independently adjusted. This enables effective modification of the functional form of the disjoining pressure so that one can vary the static and dynamic hysteresis on surfaces that exhibit the same equilibrium contact angle. The modeled solid-fluid interface is diffuse, represented by a wall probability function which ultimately controls the momentum exchange between solid and fluid phases. This approach allows us to effectively vary the slip length for a given wettability (i.e. the static contact angle) of the solid substrate.

Micro-or nano-structurally roughened solid surfaces exhibit a rich variety of wetting behavior types, ranging from superhydro-or superoleophobicity to superhydro-or superoleophilicity. Depending on their material chemistry, the scale and morphology of their roughness or even the application of external electric fields, their apparent wettability can be significantly modified giving rise to challenging technological applications by exploiting the associated capillary phenomena at the micrometer scale. Certain applications, however, are limited by hysteretic wetting transitions, which inhibit spontaneous switching between wetting states, requiring external stimuli or actuation like thermal heating. The presence of surface roughness, necessary for the manifestation of the superhydrophobicity, induces multiplicity of wetting states and the inevitable hysteresis appears due to considerable energy barriers separating the equilibrium states. Here, by using continuum as well as mesoscopic computational analysis we perform a systems level study of the mechanisms of wetting transitions on model structured solid surfaces. By tracing entire equilibrium solution families and determining their relative stability we are able to illuminate mechanisms of wetting transitions and compute the corresponding energy barriers. The implementation of our analysis to 'real world' structured or unstructured surfaces is straightforward, rendering our computational tools valuable not only for the realization of surfaces with addressable wettability through roughness design, but also for the design of suitable actuation for optimal switching between wetting states.

Collapse (Cassie to Wenzel) wetting transitions impede the electrostatically induced reversible modification of wettability on superhydrophobic surfaces, unless a strong external actuation (e.g., substrate heating) is applied. Here we show that collapse transitions can be prevented (the droplet remains suspended on the solid roughness protrusions) when the electrostatic force, responsible for the wetting modification, is smoothly distributed along the droplet surface. The above argument is initially established theoretically and then verified experimentally.