Analysis of tear film rupture: effect of non-Newtonian rheology

Abstract

We investigate the rupture mechanism of a precorneal thin mucus coating sandwiched between the aqueous tear film and the corneal epithelial surface with a monolayer of surfactant overlying the aqueous layer. The Ostwald constitutive relation is employed to model mucus and a linear equation of state describing the relationship between surface tension and surfactant concentration is adopted. Three nonlinear coupled evolution equations governing the transport of surfactant, mucus, and total liquid layer thicknesses, based on lubrication theory and a perturbation expansion technique, have been derived. The resulting equations are solved numerically in order to explore the influence of the rheological properties of mucus, aqueous–mucus thickness ratio, aqueous–mucus interfacial tension, Marangoni number, and surfactant concentration on both the onset of instability and tear film evolution in the presence of van der Waals interactions, which could rupture the tear film. Our results reveal that the influence of rheological properties, aqueous–mucus thickness ratio, and interfacial tension on the time required for film rupture can be significant and varies considerably, depending on the magnitude of the Hamaker constants governing the strength of the van der Waals forces.