An official website of the United States government
Purpose: Little quantitative or mechanistic information about tear film breakup can be determined directly via current imaging techniques. In this paper, we present simplified mathematical models based on two proposed mechanisms of tear film breakup: evaporation of water from the tear film and tangential fluid flow within the tear film. We use our models to determine whether one or a combination of the two mechanisms causes tear film breakup in a variety of instances. In this study, we estimate related breakup parameters that cannot currently be measured in breakup during subject trials, such as tear film osmolarity and thinning rates. The present study validates our procedure against previous work.Methods: Five ordinary differential equation models for tear film thinning were designed that model evaporation, osmosis, and various types of tangential flow. Eight tear film breakup instances occurring within a time interval of 1–8 s postblink of five healthy subjects thatwere identified in fluorescence images in previous work were fit with these five models. The fitting procedure used a nonlinear least squares optimization that minimized the difference of the computed theoretical fluorescent intensity from the models and the experimental fluorescent intensity from the images. The optimization was conducted over the evaporation rate and up to three tangential flow rate parameters. The smallest norm of the difference was determined to correspond to the model that best explained the tear film dynamics.Results: All of the breakup instances were best fit by models with time-dependent tangential flow. Our optimal parameter values and thinning rate as well as tangential fluid flow profiles compare well with previous partial differential equation model results in most instances.Conclusion: Our fitting results suggest that a combination of tangential fluid flow and evaporation cause most of the breakup instances. Comparison with results from previous work suggests that the simplified models can capture the essential tear film dynamics in most cases, thereby validating this procedure for wider usage.
more »
« less
The influence of a lipid reservoir on the tear film formation
Abstract We present a mathematical model to study the influence of a lipid reservoir, seen experimentally, at the lid margin on the formation and relaxation of the tear film during a partial blink. Applying the lubrication limit, we derive two coupled non-linear partial differential equations characterizing the evolution of the aqueous tear fluid and the covering insoluble lipid concentration. Departing from prior works, we explore a new set of boundary conditions (BCs) enforcing hypothesized lipid concentration dynamics at the lid margins. Using both numerical and analytical approaches, we find that the lipid-focused BCs strongly impact tear film formation and thinning rates. Specifically, during the upstroke of the eyelid, we find specifying the lipid concentration at the lid margin accelerates thinning. Parameter regimes that cause tear film formation success or failure are identified. More importantly, this work expands our understanding of the consequences of lipid dynamics near the lid margins for tear film formation.
more »
« less
- Award ID(s):
- 1909846
- PAR ID:
- 10258072
- Date Published:
- Journal Name:
- Mathematical Medicine and Biology: A Journal of the IMA
- Volume:
- 37
- Issue:
- 3
- ISSN:
- 1477-8599
- Page Range / eLocation ID:
- 363 to 388
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
The human tear film is a multilayer structure in which the dynamics are often strongly affected by a floating lipid layer. That layer has liquid crystalline characteristics and plays important roles in the health of the tear film. Previous models have treated the lipid layer as a Newtonian fluid in extensional flow. Motivated to develop a more realistic treatment, we present a model for the extensional flow of thin sheets of nematic liquid crystal. The rod-like molecules of these substances impart an elastic contribution to the rheology. We rescale a weakly elastic model due to Cummings et al. [“Extensional flow of nematic liquid crystal with an applied electric field,” Eur. J. Appl. Math. 25, 397–423 (2014).] to describe a lipid layer of moderate elasticity. The resulting system of two nonlinear partial differential equations for sheet thickness and axial velocity is fourth order in space, but still represents a significant reduction of the full system. We analyze solutions arising from several different boundary conditions, motivated by the underlying application, with particular focus on dynamics and underlying mechanisms under stretching. We solve the system numerically, via collocation with either finite difference or Chebyshev spectral discretization in space, together with implicit time stepping. At early times, depending on the initial film shape, pressure either aids or opposes extensional flow, which changes the free surface dynamics of the sheet and can lead to patterns reminiscent of those observed in tear films. We contrast this finding with the cases of weak elasticity and Newtonian flow, where the sheet retains the same qualitative shape throughout time.more » « less
-
null (Ed.)Tear film plays a key role in protecting the cornea surface against contaminations and dry eye syndrome which can lead to symptoms of discomfort, visual trouble, and tear film instability with the potential to damage the ocular surface. In this paper, coupled nonlinear partial differential equations of the fourth order proposed by Aydemir et al. to describe the evolution of tear film dynamics are considered. These equations are of Benney type and known to suffer from unbounded behavior and lack of a global attractor. The objective here is to identify a reduced order modeling framework with the potential to be used as a basis for control in future work using smart tears with a surfactant that can modify the surface tension to prevent tear film breakup. Since the dynamics are infinite dimensional and nonlinear, a reduced order model based on the proper orthogonal decomposition (POD) is developed, analyzed, and compared to the full order model. Numerical simulations illustrate that only a small number of POD modes are required to accurately capture the tear film dynamics allowing for the full partial differential model to be represented as a low-dimensional set of coupled ordinary differential equations.more » « less
-
Abstract Despite significant study, when and how plate tectonics initiated on Earth remains contentious. Geologic evidence from some of Earth's earliest cratons has been interpreted as reflecting the formation of initial continental blocks by non‐subduction processes, which then trigger subduction initiation at their margins. Numerical models of mantle convection with a plastic yield stress rheology have shown this scenario is plausible. However, whether continents can trigger subduction initiation has not been tested with other rheologies. We, therefore, use numerical models of mantle convection with an imposed continental block to test whether continents facilitate subduction initiation with a grain‐damage mechanism, where weak shear zones form by grain size reduction. Our results show that continents modestly enhance stresses in the lithosphere, but not enough to significantly impact lithospheric damage or subduction initiation: continents have minimal influence on lithospheric damage or plate speed, nor does subduction preferentially initiate at the continental margin. A new regime diagram that includes continental blocks shows only a small shift in the boundary between the mobile‐lid and stagnant‐lid regimes when continents are added. However, as we do find that stresses are modestly enhanced at the continental margin in our models, we develop a scaling law for this stress enhancement to more fully test whether continents could trigger subduction initiation on early Earth. We find that lithospheric stresses supplied by continents are not sufficient to initiate subduction on the early Earth on their own with grain‐damage rheology; instead, additional factors would be required.more » « less
-
Abstract We describe elongate, wet, subglacial bedforms in the shear margins of the NE Greenland Ice Stream and place some constraints on their formation. Lateral shear margin moraines have been observed across the previously glaciated landscape, but little is known about the ice-flow conditions necessary to form these bedforms. Here we describe in situ sediment bedforms under the NE Greenland Ice Stream shear margins that are observed in active-source seismic and ground-penetrating radar surveys. We find bedforms in the shear margins that are ~500 m wide, ~50 m tall, and elongated nearly parallel to ice-flow, including what we believe to be the first subglacial observation of a shear margin moraine. Acoustic impedance analysis of the bedforms shows that they are composed of unconsolidated, deformable, water-saturated till. We use these geophysical observations to place constraints on the possible formation mechanism of these subglacial features.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1093/imammb/dqz018
