Meniscal tear film fluid dynamics near Marx's line.

Extensive studies have explored the dynamics of the ocular surface fluid, though theoretical investigations are typically limited to the use of the lubrication approximation, which is not guaranteed to be uniformly valid a-priori throughout the tear meniscus. However, resolving tear film behaviour within the meniscus and especially its apices is required to characterise the flow dynamics where the tear film is especially thin, and thus most susceptible to evaporatively induced hyperosmolarity and subsequent epithelial damage. Hence, we have explored the accuracy of the standard lubrication approximation for the tear film by explicit comparisons with the 2D Navier-Stokes model, considering both stationary and moving eyelids. Our results demonstrate that the lubrication model is qualitatively accurate except in the vicinity of the eyelids. In particular, and in contrast to lubrication theory, the solution of the full Navier-Stokes equations predict a distinct absence of fluid flow, and thus convective mixing in the region adjacent to the tear film contact line. These observations not only support emergent hypotheses concerning the formation of Marx's line, a region of epithelial cell staining adjacent to the contact line on the eyelid, but also enhance our understanding of the pathophysiological consequences of the flow profile near the tear film contact line.

Kinetics of surfactant desorption at an air-solution interface.

The kinetics of re-equilibration of the anionic surfactant sodium dodecylbenzene sulfonate at the air-solution interface have been studied using neutron reflectivity. The experimental arrangement incorporates a novel flow cell in which the subphase can be exchanged (diluted) using a laminar flow while the surface region remains unaltered. The rate of the re-equilibration is relatively slow and occurs over many tens of minutes, which is comparable with the dilution time scale of approximately 10-30 min. A detailed mathematical model, in which the rate of the desorption is determined by transport through a near-surface diffusion layer into a diluted bulk solution below, is developed and provides a good description of the time-dependent adsorption data. A key parameter of the model is the ratio of the depth of the diffusion layer, H(c), to the depth of the fluid, H(f), and we find that this is related to the reduced Péclet number, Pe*, for the system, via H(c)/H(f) = C/Pe*(1/2). Although from a highly idealized experimental arrangement, the results provide an important insight into the "rinse mechanism", which is applicable to a wide variety of domestic and industrial circumstances.

Coupling fluid and solute dynamics within the ocular surface tear film: a modelling study of black line osmolarity.

We present a mathematical model describing the spatial distribution of tear film osmolarity across the ocular surface of a human eye during one blink cycle, incorporating detailed fluid and solute dynamics. Based on the lubrication approximation, our model comprises three coupled equations tracking the depth of the aqueous layer of the tear film, the concentration of the polar lipid, and the concentration of physiological salts contained in the aqueous layer. Diffusive boundary layers in the salt concentration occur at the thinnest regions of the tear film, the black lines. Thus, despite large Peclet numbers, diffusion ameliorates osmolarity around the black lines, but nonetheless is insufficient to eliminate the build-up of solute in these regions. More generally, a heterogeneous distribution of solute concentration is predicted across the ocular surface, indicating that measurements of lower meniscus osmolarity are not globally representative, especially in the presence of dry eye. Vertical saccadic eyelid motion can reduce osmolarity at the lower black line, raising the prospect that select eyeball motions more generally can assist in alleviating tear film hyperosmolarity. Finally, our results indicate that measured evaporative rates will induce excessive hyperosmolarity at the black lines, even for the healthy eye. This suggests that further evaporative retardation at the black lines, for instance due to the cellular glycocalyx at the ocular surface or increasing concentrations of mucus, will be important for controlling hyperosmolarity as the black line thins.

The effect of polar lipids on tear film dynamics.

In this paper, we present a mathematical model describing the effect of polar lipids, excreted by glands in the eyelid and present on the surface of the tear film, on the evolution of a pre-corneal tear film. We aim to explain the interesting experimentally observed phenomenon that the tear film continues to move upward even after the upper eyelid has become stationary. The polar lipid is an insoluble surface species that locally alters the surface tension of the tear film. In the lubrication limit, the model reduces to two coupled non-linear partial differential equations for the film thickness and the concentration of lipid. We solve the system numerically and observe that increasing the concentration of the lipid increases the flow of liquid up the eye. We further exploit the size of the parameters in the problem to explain the initial evolution of the system.

On the predictions and limitations of the Becker-Döring model for reaction kinetics in micellar surfactant solutions.

We investigate the breakdown of a system of micellar aggregates in a surfactant solution following an order-one dilution. We derive a mathematical model based on the Becker-Döring system of equations, using realistic expressions for the reaction constants fit to results from Molecular Dynamics simulations. We exploit the largeness of typical aggregation numbers to derive a continuum model, substituting a large system of ordinary differential equations for a partial differential equation in two independent variables: time and aggregate size. Numerical solutions demonstrate that re-equilibration occurs in two distinct stages over well-separated timescales, in agreement with experiment and with previous theories. We conclude by exposing a limitation in the Becker-Döring theory for re-equilibration of surfactant solutions.

The role of cell-cell interactions in a two-phase model for avascular tumour growth.

A two-phase model is presented to describe avascular tumour growth. Conservation of mass equations, including oxygen-dependent cell growth and death terms, are coupled with equations of momentum conservation. The cellular phase behaves as a viscous liquid, while the viscosity of the extracellular water manifests itself as an interphase drag. It is assumed that the cells become mechanically stressed if they are too densely packed and that the tumour will try to increase its volume in order to relieve such stress. By contrast, the overlapping filopodia of sparsely populated cells create short-range attractive effects. Finally, oxygen is consumed by the cells as it diffuses through the tumour. The resulting system of equations are reduced to three, which describe the evolution of the tumour cell volume fraction, the cell speed and the oxygen tension. Numerical simulations indicate that the tumour either evolves to a travelling wave profile, in which it expands at a constant rate, or it settles to a steady state, in which the net rates of cell proliferation and death balance. The impact of varying key model parameters such as cellular viscosity, interphase drag, and cellular tension are discussed. For example, tumours consisting of well-differentiated (i.e. viscous) cells are shown to grow more slowly than those consisting of poorly-differentiated (i.e. less viscous) cells. Analytical results for the case of oxygen-independent growth are also presented, and the effects of varying the key parameters determined (the results are in line with the numerical simulations of the full problem). The key results and their biological implications are then summarised and future model refinements discussed.