Gravity-driven instability of a thin liquid film underneath a soft solid.

The gravity-driven instability of a thin liquid film located underneath a soft solid material is considered. The equations and boundary conditions governing the solid deformation are systematically converted from a Lagrangian representation to an Eulerian representation, which is the natural framework for describing the liquid motion. This systematic conversion reveals that the continuity-of-velocity boundary condition at the liquid-solid interface is more complicated than has previously been assumed, even in the small-strain limit. We then make clear the conditions under which the commonly used simplified version of this boundary condition is valid. The small-strain approximation, lubrication theory, and linear stability analysis are applied to derive an expression for the growth rate of small-amplitude perturbations. Asymptotic analysis reveals that the coupling between the liquid and solid manifests itself as a lower effective liquid-air interfacial tension that leads to larger instability growth rates. Although this suggests that it is more difficult to maintain a stable liquid coating underneath a soft solid, the effect is expected to be weak for cases of practical interest.

Tear film dynamics with evaporation, wetting, and time-dependent flux boundary condition on an eye-shaped domain.

We study tear film dynamics with evaporation on a wettable eye-shaped ocular surface using a lubrication model. The mathematical model has a time-dependent flux boundary condition that models the cycles of tear fluid supply and drainage; it mimics blinks on a stationary eye-shaped domain. We generate computational grids and solve the nonlinear governing equations using the OVERTURE computational framework. In vivo experimental results using fluorescent imaging are used to visualize the influx and redistribution of tears for an open eye. Results from the numerical simulations are compared with the experiment. The model captures the flow around the meniscus and other dynamic features of human tear film observed in vivo.

An overset grid method for the study of reflex tearing.

We present an overset grid method to simulate the evolution of human tear film thickness subject to reflex tearing. The free-surface evolution is governed by a single fourth-order non-linear equation derived from lubrication theory with specified film thickness and volume flux at each end. The model arises from considering the limiting case where the surfactant is strongly affecting the surface tension. In numerical simulations, the overset grid is composed of fine boundary grids near the upper and lower eyelids to capture localized capillary thinning referred to as 'black lines' and a Cartesian grid covers the remaining domain. Numerical studies are performed on a non-linear test problem to confirm the accuracy and convergence of the scheme. The computations on the tear film model show qualitative agreement with in vivo tear film thickness measurements. Furthermore, the role of the black lines in the presence of tear supply from the lid margins, reflex tearing, was found to be more subtle than a barrier to tear fluid flow between the anterior of the eye and the meniscus at the lid margin. During reflex tearing, tears may flow through the region normally containing the black line and drift down over the cornea under the influence of gravity.

Single-equation models for the tear film in a blink cycle: realistic lid motion.

We consider model problems for the tear film over multiple blink cycles that utilize a single equation for the tear film; the single non-linear partial differential equation that governs the film thickness arises from lubrication theory. The two models that we consider arise from considering the absence of naturally occurring surfactant and the case when the surfactant is strongly affecting the surface tension. The film is considered on a time-varying domain length with specified film thickness and volume flux at each end; only one end of the domain is moving, which is analogous to the upper eyelid moving with each blink. Realistic lid motion from observed blinks is included in the model with end fluxes specified to more closely match the blink cycle than those previously reported. Numerical computations show quantitative agreement with in vivo tear film thickness measurements under partial blink conditions. A transition between periodic and non-periodic solutions has been estimated as a function of closure fraction and this may be a criterion for what is effectively a full blink according to fluid dynamics.