
Citation: Luís J. del Valle, Lourdes Franco, Ramaz Katsarava, Jordi Puiggalí. Electrospun biodegradable polymers loaded with bactericide agents[J]. AIMS Molecular Science, 2016, 3(1): 52-87. doi: 10.3934/molsci.2016.1.52
[1] | Lin Chen, Mingjun Wang, Zhanyu Wu, Jinbo Sun, Jianglong Li, Chun Chen, Chuan Ye . Advancements in health informatics: finite element insights into medial open-wedge high tibial osteotomy and lateral meniscal tears. Mathematical Biosciences and Engineering, 2024, 21(4): 5394-5410. doi: 10.3934/mbe.2024237 |
[2] | Ke Li, Qiaolin Zhu, Jianzhang Wu, Juntao Ding, Bo Liu, Xixi Zhu, Shishi Lin, Wentao Yan, Wulan Li . DCT-Net: An effective method to diagnose retinal tears from B-scan ultrasound images. Mathematical Biosciences and Engineering, 2024, 21(1): 1110-1124. doi: 10.3934/mbe.2024046 |
[3] | Christophe Prud'homme, Lorenzo Sala, Marcela Szopos . Uncertainty propagation and sensitivity analysis: results from the Ocular Mathematical Virtual Simulator. Mathematical Biosciences and Engineering, 2021, 18(3): 2010-2032. doi: 10.3934/mbe.2021105 |
[4] | Ziad Khan, Hari Mohan Srivastava, Pshtiwan Othman Mohammed, Muhammad Jawad, Rashid Jan, Kamsing Nonlaopon . Thermal boundary layer analysis of MHD nanofluids across a thin needle using non-linear thermal radiation. Mathematical Biosciences and Engineering, 2022, 19(12): 14116-14141. doi: 10.3934/mbe.2022658 |
[5] | Giovanna Guidoboni, Alon Harris, Lucia Carichino, Yoel Arieli, Brent A. Siesky . Effect of intraocular pressure on the hemodynamics of the central retinal artery: A mathematical model. Mathematical Biosciences and Engineering, 2014, 11(3): 523-546. doi: 10.3934/mbe.2014.11.523 |
[6] | H. T. Banks, Robert Baraldi, Karissa Cross, Kevin Flores, Christina McChesney, Laura Poag, Emma Thorpe . Uncertainty quantification in modeling HIV viral mechanics. Mathematical Biosciences and Engineering, 2015, 12(5): 937-964. doi: 10.3934/mbe.2015.12.937 |
[7] | Mohamed Abbas, Ali Algahtani, Amir Kessentini, Hassen Loukil, Muneer Parayangat, Thafasal Ijyas, Abdul Wase Mohammed . Qualitative theoretical modeling to study the possibility of detecting multi-virus in blood flow using Nano-quartz crystal microbalance. Mathematical Biosciences and Engineering, 2020, 17(5): 4563-4577. doi: 10.3934/mbe.2020252 |
[8] | Xiujuan Zheng, Zhanheng Li, Xinyi Chun, Xiaomei Yang, Kai Liu . A model-based method with geometric solutions for gaze correction in eye-tracking. Mathematical Biosciences and Engineering, 2020, 17(2): 1396-1412. doi: 10.3934/mbe.2020071 |
[9] | Li Wang, Changyuan Wang, Yu Zhang, Lina Gao . An integrated neural network model for eye-tracking during human-computer interaction. Mathematical Biosciences and Engineering, 2023, 20(8): 13974-13988. doi: 10.3934/mbe.2023622 |
[10] | Antonio Fasano, Marco Gabrielli, Alberto Gandolfi . Investigating the steady state of multicellular spheroids by revisiting the two-fluid model. Mathematical Biosciences and Engineering, 2011, 8(2): 239-252. doi: 10.3934/mbe.2011.8.239 |
The tear film is critical for nourishment and protection for the ocular surface [1,2] as well as providing a smooth optical surface essential for good vision. [3,4]. A sketch of the lacrimal supply and drainage system is shown in Figure 1. The tear film is painted onto the ocular surface during the upstroke of the upper eyelid during the blink; it is a thin liquid film with multiple layers that, in health, establishes itself rapidly after each blink.
At the anterior interface with air is an oily lipid layer that primarily retards evaporation [7,8], which helps to retain a smooth, well-functioning tear film [9]. Posterior to the lipid layer is the aqueous layer, which consists mostly of water [10]. At the corneal surface, there is a region with transmembrane mucins protruding from the epithelial cells of the cornea or conjunctiva. This forest of glycosolated mucins is called the glycocalyx. It is generally agreed that the presence of the hydrophilic glycocalyx on the ocular surface prevents the tear film from dewetting [11,12,13]. The overall thickness of the tear film is a few microns [14], while the average thickness of the lipid layer is on the order of tens to 100 nanometers [9,15,16,17,18] and the thickness of the glycocalyx is a few tenths of a micron [2]. This overall structure is reformed on the order of a second after each blink in a healthy tear film.
The aqueous part of tear fluid is primarily supplied from the lacrimal gland near the temporal canthus (corner) and the excess is drained through the puncta near the nasal canthus. Doane [5] proposed the mechanism of tear drainage in vivo whereby tear fluid is drained into the canaliculi through the puncta during the opening interblink phase. Water lost from the tear film due to evaporation into air is an important process also [7,19,20]. We believe that this is the primary mechanism by which the osmolarity (essentially, the total concentration of ions on a molar basis) is increased in the tear film [18]. Some water is supplied from the ocular epithelia via osmosis when the osmolarity increases above the isotonic level (the so-called hyperosmolar state) [18,21,22]. According to the Dry Eye WorkShop (DEWS) reports, the osmolarity is important because it plays an important etiological role in the development of dry eye disease [23], which affects millions [24,25]. An accepted view in these DEWS reports is that repeated exposure to sufficient hyperosmolarity in tears, as may occur in tear breakup, can result in inflammation and ocular surface damage in a vicious cycle that may can lead to dry eye disease. There is still considerable ongoing research to more fully understand the causes and etiology of dry eye disease; some mathematical modeling results are discussed in this paper help to quantify the relationship between tear breakup and hyperosmolarity.
The characteristic ranges of osmolarity on the ocular surface may be summarized as in [26] and [27]; the following closely follows those sources. The healthy tear film maintains homeostasis with the blood in the range 296-302 mOsM (or, equivalently, mOsm/L or Osm/m$ ^3 $) [28,29,30], while healthy blood is in the range 285-295 mOsM [31]. In dry eye disease, the lacrimal system may be unable to maintain this homeostasis and osmolarity values in the meniscus rise to 316-360 mOsM [29,32,33,34], and even higher values may be attained over the cornea. Using in vivo experiment and sensory feedback, Liu et al. [35] estimated peak values of 800-900 mOsM in tear break up (TBU). Similar or higher values were computed from mathematical models of TBU in [18,36], or for models of the whole ocular surface [27]. The estimates from these models take on added significance since the osmolarity in TBU or in most locations cannot be measured directly; the only available clinical method used in the clinic measures the osmolarity in the inferior meniscus [28].
The ongoing supply and drainage of tear fluid affects the distribution and flow of the tear film. A number of methods have been used to visualize and/or measure tear film thickness and flow, including interferometry [14,37,38], optical coherence tomography [39], fluorescence imaging [40,41,42,43] and many others. There are different ways that fluorescence may be used to visualize the tear film [44]. In the dilute limit, the fluorescein concentration is below the critical concentration, and the intensity of the fluorescence from the tear film is proportional to its thickness. In the concentrated (or self-quenching) limit, the intensity drops as the tear film thins in response to evaporation, and the thickness is roughly proportional to the square root of the intensity for a spatially uniform (flat) tear film [45,46]. We use the shorthand FL to indicate fluorescence or fluorescein below.
Recently, the DEWS reports [25] and others (summarized in [18]) have argued that TBU is driven, at least in many cases, by evaporation. Evaporation causes increased osmolarity in the tear film, leading to irritation, inflammation and damage to the ocular epithelium. Currently, the authors are unaware of any directly measurements of the osmolarity in regions of TBU over the cornea. Theoretical efforts have responded to this situation by creating models that incorporate evaporation and osmosis, as well as other effects, into the mathematical models to better understand the dynamics of the process at a small scale. These models emphasize understanding the overall dynamics of the thickness, flow, and transport of solutes inside the film. These predictions can produce insights beyond what direct measurements can currently provide, and include significant new results about the distribution of osmolarity in TBU. We believe that these results impact the understanding of osmolarity dynamics in the tear film as well as the understanding of FL imaging of TBU dynamics.
A variety of mathematical models have incorporated various important effects of tear film dynamics as has been recently reviewed [18,21,47]. The most common assumptions for these models are a Newtonian tear fluid and a flat cornea [48,49]. This paper emphasizes local partial differential equation (PDE) models for tear film dynamics in TBU. Beyond previous reviews, we include some recent results using non-Newtonian fluids and non-uniform properties of the ocular surface. We also review new results for fitting local thin film models to in vivo imaging data; these are reviewed in the Results section. Opportunities for improvements are suggested in the Discussion section. Conclusions follow thereafter.
This paper reviews mathematical models and theories for human tear film dynamics. We typically refer the reader to the original works for derivations and detailed discussions of solution methods. Generally, multiple scale perturbation methods are applied to the equations of fluid dynamics and transport of mass and energy to derive the so-called thin film or lubrication eqautions for the quantities of interest. Reviews of the sizable literature on thin films include [50,51,52]. Recent reviews of tear film models and dynamics include [18,21,47,53]. The resulting systems are typically dissipative, in one- or two-dimensional domains, and may be solved via a variety of numerical methods; a good summary is given in Maki et al. ([54], p. 364).
Quite a few parameters appear in the equations. Table 1 gives typical dimensional values that are needed in these models (after [55]). The thickness of the tear film ($ d $) is about 100 times less the length scale along the tear film ($ \ell $); this small ratio enables the use of the thin film approximation to derive governing equations for the dependent variables.
Dimensional parameters | |||
Parameter | Description | Value | Reference |
$ \mu $ | Viscosity | $ 1.3 \times 10^{-3} $ Pa $ \cdot $ s | [58] |
$ \sigma_0 $ | Surface tension | 0.045 N $ \cdot $ m$ ^{-1} $ | [59] |
$ \rho $ | Density | $ 10^3 $ kg $ \cdot $ m$ ^{-3} $ | Water |
$ A^* $ | Hamaker constant | $ 6 \pi \times 3.5 \times 10^{-19} $ s $ \cdot $ m$ ^{-1} $ | [60] |
$ d $ | Initial TF thickness | $ 3-6 \times 10^{-6} $ m | Calculated; see [61] |
$ \ell $ | $ (\sigma_0/\mu/v_{\max})^{1/4}d $ | 0.426-0.694 mm | Calculated |
$ v_{\max} $ | Peak thinning rate | 0.5-25 $ \mu $m/min | [56] |
$ V_w $ | Molar volume of water | 1.8 $ \times 10^{-5} $ m$ ^3 \cdot $ mol$ ^{-1} $ | Water |
$ h_{\text{eq}} $ | Equilibrium thickness | 0.25 $ \mu $m | Estimated, [13] |
$ D_f $ | Diffusivity of fluorescein | 0.39 $ \times 10^{-9} $ m$ ^2 $/s | [62] |
$ D_o $ | Diffusivity of salt | 1.6 $ \times 10^{-9} $ m$ ^2 $/s | [63] |
$ c_0 $ | Isotonic osmolarity | 300 Osm/m$ ^3 $ | [28] |
$ P_o $ | Permeability of cornea | $ 12.1 \times 10^{-6} $ m/s | [18] |
$ \epsilon_f $ | Napierian extinction coefficient | $ 1.75 \times 10^7 $ m$ {-1} $ M$ ^{-1} $ | [64] |
Overall flows over a two-dimensional (2D) eye-shaped domain with stationary boundaries have been able to capture thinning near the lid margins due to the presence of a meniscus there, and flow within the meniscus along each lid margin. Tear film models are often formulated on a one-dimensional (1D) domain oriented vertically through the center of the cornea with stationary ends corresponding to the eyelid margins. We refer to models on this kind of domain as 1D models. Surface tension, viscosity, gravity and evaporation are often incorporated into 1D models; wetting forces have been included as well. Due to space considerations we refer the reader to a recent review [21].
In this section we briefly discuss two-dimensional (2D) results over the open eye shape. For a stationary boundary (no blinking), thin film fluid (lubrication) models have been developed for the dynamics of the tear film. A series of papers developing these models have used a piecewise quartic polynomial fit to the lid margins in an image of an eye as the boundary. The first papers used time-independent pressure [65] and flux [54] boundary conditions. The thickness of the tear film and the flows inside it were computed. Thinning occurred near the lid margin, and flow around the menisci was observed. Subsequent work added time-dependent to the boundary conditions [66], evaporation and osmolarity transport [27], and fluorescein transport and fluorescence [67]. Evaporation and osmolarity transport revealed that hyperosmolarity occurred near the eyelid margins (boundaries) in the so-called "black line" where the tear film is thinnest. Adding fluorescein and fluorescence captured elements of a commonly used procedure for imaging the tear film in the clinic, and set the stage for quantitative comparison with these imaging techniques.
A sample of papers that review blinking includes [68,69]. For a recent review of tear film dynamics and blinking observations, together with some mathematical modeling, the reader is referred to [18]. Of particular interest with regard to visualizing tear film dynamics are high-resolution imaging of the lipid layer in both space and time during blinks, simultaneous imaging with fluorescence, and a survey of methods for imaging the tear film. The latter work also includes first efforts to closely compare computed and observed FL intensities in TBU.
Recently, an approach that incorporates blinking of a 2D eye-shaped domain has been developed. A model problem treats the eyelid margins as sectors of circles and uses conformal mapping methods to map the eye shape to more convenient domains for computing and moving the upper lid [70]. That work solved a model second-order nonlinear diffusion problem on the blinking domain to high accuracy. Combining the moving domain with thin film fluid flow equations is in progress (TA Driscoll and RJ Braun, in preparation). Preliminary results for one model, corresponding to a model for pure water without evaporation, appeared in [47].
The situation for tear film dynamics on moving 1D domains representing the opening of the eyelid and the complete blink cycle has been treated in a number of papers. Some models that combine the eyelid upstroke and interblink include [71,72,73], and a two layer model with a dynamic lipid layer [74]. Several models have included full blink cycles, including fluid dynamics only [75]; flow and heat transfer [76,77]; and flow, evaporation and osmolarity [78]. For a recent summary of these models, the reader is referred to [47].
We briefly discuss two papers that included the lid upstroke and interblink but simplified the lipid layer to its polar lipid component by treating it as an insoluble surfactant. In Aydemir et al. [71], the boundary conditions on the insoluble surfactant were of no-flux type; the aqueous flux at the ends was specified. This choice fixed the total amount of surfactant on the tear film surface. They were able to analyze some limits, e.g., the early time dynamics, via a perturbation approach. They obtained good agreement between the analysis and the numerical solutions of the full model. Maki et al. [79] used Dirchlet boundary conditions on the insoluble surfactant in an effort to study a hypothesized lipid reservoir at the lid margin [1]. They also studied different conditions on the aqueous layer. The model could be analyzed in some limits, yielding good agreement between theory and numerical solutions for parts of the dynamics. The agreement with expectations of tear film dynamics appears to be better for the no-flux conditions in most aspects of the results; however, taken together with the results for the surfactant distributions in two-layer models, there are still aspects of lipid layer dynamics that warrant further study.
We note that Zubkov et al. [80] studied flow near the lid margin. They compared flows in the meniscus using 2D Navier-Stokes equations and lubrication approximations near the domain boundary corresponding to the lids. The two approaches showed qualitative agreement there, but the lubrication theory was unable to capture some fine details that the Navier-Stokes model exhibited. Thus, models that rely on lubrication theory for fine detail of meniscus flows should be viewed as qualitative or model problems. We also note that no models to our knowledge have taken the fine details of the geometry of the cross section of the lid margin into account.
We now turn to solving local problems for the failure of the tear film via TBU. One solves for the thickness, osmolarity and fluorescein concentration that vary with both time and one space dimension. There are no menisci in these local problems; they are applied to a small region of interest within the open eye domain. The shape of TBU regions observed using FL imaging may be complex, and they are dynamic. PDE models of TBU have been solved only for simple shapes so far. Axisymmetric circular shapes, corresponding to "spots, " have been solved in an axisymmetric polar coordinate system. Linear shapes, corresponding to "streaks, " have been solved in Cartesian coordinates. We limit our discussion to these simplified geometries here. Such models can approximate and produce comparable data for simple TBU regions seen in vivo [18].
According to driving mechanism, there are at least two types of TBU: evaporative and rapid glob-driven. Figure 2 shows sketches of the different effects that are present and their relative importance in TBU when these respective mechanisms are dominant. The evaporative mechanism is thought to be from defects in the lipid layer; a local thin region of the lipid layer is sufficient to drive faster localized thinning and subsequent TBU [82]. The Marangoni effect is caused by differences in lipid layer composition, in particular for the surface active components at the lipid/aqueous layer interface. The Marangoni effect is known to be active just after a blink [48,83], where the formation of new area during eye opening causes a lower surfactant concentration near the upper lid. It has also been been hypothesized to be the driving force in some instances of rapid TBU [84,85], where uneven lipid spreading after the blink can cause divergent tangential flow, and associated local thinning, of the tear film.
A third mechanism is thought to be important as well: dewetting. The ocular surface may be defective, causing it to be unwettable [86,87]; the idea has been discussed in the context of TBU types [88] and for dry eye disease diagnosis [89]. In this case, a pressure gradient is developed inside the tear film due to van der Waals forces pushing fluids from between the air/tear and tear/cornea interfaces [90]. This mechanism typically causes accelerated decrease of thickness as the film becomes thinner (e.g., [52]); however, the level of detail seen in carefully controlled physical systems may not be readily achieved in the relatively messy tear film system. This could make observation of accelerated thinning difficult in humans, though animal models may be more amenable to modification for detailed observation of dewetting-driven TBU.
We first consider evaporative TBU, then proceed to rapid TBU.
The thin film approximation, or lubrication theory results in a simplified system of equations to solve for the tear film dynamics. The following equations apply to axisymmetric spot TBU; the derivation and the Cartesian versions of the equations may be found in [91].
One must solve
$ ∂th=1r∂r(rhˉu)−Je+Jo, $
|
(3.1) |
$ ˉu=−h212r∂rp, $
|
(3.2) |
$ p=−1r∂r(r∂rh)−Ah−3, $
|
(3.3) |
$ h(∂tm+ˉu∂rm)=Pe−1m1r∂r(rh∂rm)+(Je−Jo)m, m=c,f, $
|
(3.4) |
$ Jo=Pc(c−1). $
|
(3.5) |
These equations are solved for the tear film thickness $ h(r, t) $, the pressure $ p(r, t) $, the osmolarity $ c(r, t) $, and the fluorescein concentration $ f(r, t) $. The domain is $ 0 < r < R_O $, where $ R_O $ is the outer radius of the domain. The nondimensional parameters are summarized in Table 2, with some typical values for the scalings used in [91].
Nondimensional parameters with typical values | ||
Parameter | Expression | Value |
$ \epsilon $ | $ \frac{d}{\ell} $ | $ 1\times10^{-2} $ |
$ S $ | $ \frac{\sigma_0 \epsilon^4}{\mu v_{max}} $ | 1 |
$ A $ | $ \frac{A^*}{\mu v_{max}~~ \ell^2} $ | $ 0.9 \times 10^{-3} $ |
$ \alpha $ | $ \frac{\alpha_0 \mu}{\rho \ell \epsilon^3} $ | 0.0368 |
$ P_c $ | $ \frac{P_o V_m~~ c_0}{v_{max}} $ | 0.1960 |
$ Pe_{f} $ | $ \frac{v_{max}~~ \ell}{\epsilon D_f} $ | 30.5 |
$ Pe_{c} $ | $ \frac{v_{max}~~ \ell}{\epsilon D_o} $ | 7.43 |
$ \phi $ | $ \epsilon_f f_{cr} d $ | $ 0.326 $ |
These equations apply to axisymmetric spot TBU; those for streak TBU are similar. The pressure is treated as a dependent variable to be found as part of the numerical treatment. Treating the pressure this way requires only first order derivatives in the boundary conditions, which is favorable for the numerical methods that we employ [36,65,91,92,93].
The boundary conditions are symmetry at the origin, homogenous Neumann at $ R_O $, and the initial conditions are uniform in space:
$ h(r,0)=c(r,0)=1,f(r,0)=f0,p(r,0)=−A. $
|
(3.6) |
The initial FL concentration $ f_0 $ is relative to the critical FL concentration $ f_{cr} = 0.2 \% $.
The dynamics are driven by the evaporative flux $ J_e $, which is scaled with $ v_{max} $, the maximum volume flux measured from experiment. The evaporative flux is assumed to be composed of a time-independent distribution at the initial time, which is modified by the film pressure. The latter includes capillary and wetting terms and are motivated by one-sided models of evaporation [60]. Empirical evidence for a fixed distribution may be found in Figure 6 of [82], where a very thin region of lipid forms, then remains roughly stationary, leading to TBU by 20 seconds post blink. For in vivo experiments with controlled humidity and little airflow, measurements yield constant thinning rates over a significant time interval [56]. Thus, we expect that our assumed form of evaporation rate applies best for comparison with data from controlled laboratory conditions. Nondimensionally, we have
$ Je(r,t)=J(r)+αp. $
|
(3.7) |
We show results for $ J(r) $, the evaporation rate with a Gaussian distribution inside a disk, namely
$ J(r)=vb+(1−vb)e−(r/rw)2/2. $
|
(3.8) |
The standard deviation $ r_w $ will be used to indicate the radius of the evaporation distribution. Here $ v_b = v_{min}/v_{max} $ is the nondimensional version of the background evaporation rate $ v_{min} $, and the scaling is based the peak evaporation rate $ v_{max} $.
The FL intensity $ I $ is computed from the thickness $ h $ and the fluorescein concentration $ f $ via the nondimensional equation [44,45]
$ I=I01−e−ϕfh1+f2. $
|
(3.9) |
Here $ \phi $ is the nondimensional Napierian extinction coefficient, and $ I_0 $ is a normalization constant that takes the optical system into account. As mentioned previously, for fixed thickness $ h $ the FL intensity decreases linearly at small concentration in the so-called dilute regime. Sufficiently above a critical concentration $ f_{cr} = 0.2\% $, the intensity decreases quadratically with increasing $ f $, which we call the self-quenching regime [45]. For a tear film thinning by evaporation, the situation is slightly more complicated. For a flat tear film, the product $ hf $ is constant, and so in the dilute regime, $ I $ is constant during thinning, while in the self-quenching regime, the quadratic decrease for increasing $ f $ still holds [46].
The numerical method used to solve this model (and the other results shown herein) is an application of the method of lines that used a Chebyshev spectral discretization in space [94]. The resulting differential-algebraic system being solved using $\texttt{ode15s}$ in Matlab (MathWorks, Natick, MA, USA). This solver uses adaptive BDF time stepping with variable size time steps and variable order of 1 to 5; this solver is well suited to these stiff, high-order problems. With sufficient resolution of the spatial discretization, we found that we did not need to use positivity-preserving schemes that have met with success elsewhere, e.g., [95]. We note that finding TBU times was enabled by using event detection capability built into the solver.
In this case, we show results from [91] for the problem with PDEs given by Eqs (3.1)-(3.4) and (3.7) subject to the homogeneous Neumann boundary conditions and uniform initial conditions Eq (3.6). For illustration, the Gaussian distribution has peak thinning rate $ v_{max} = 20\; \mu $m/min, background rate $ v_{min} = 1\; \mu $m/min, and width $ r_w = 1 $; this width is the same as the length scale along the film ($ 0.35 $ mm in this case). The time scale is $ d/v_{max} = 10.5 $ s. This approximates a hole in the lipid layer allowing localized fast evaporation compared to a slow background rate. In this example, wetting forces are active ($ A\neq 0 $) and there is a pressure correction to the evaporation rate ($ \alpha\neq 0 $). The thickness is scaled with typical film thickness $ d $, the evaporative flux with the maximum evaporation rate $ \rho v_{max} $, the osmolarity $ c $ is scaled with the isotonic concentration $ c_0 $, and the fluorescein concentration $ f $ with its critical concentration.
Figure 3 shows tear film thickness $ h $, evaporative flux $ J_e $, osmolarity $ c $ and FL concentration $ f $ as a function of space for different times. We see that the tear film thickness decreases to a low level that we interpret as TBU and that this takes slightly longer than one unit of time; this confirms the suitability of our choice of time scale for a thinning rate distribution of this width. When the thickness of the tear film gets small enough, the evaporative flux is reduced according to the model. The osmolarity and FL concentrations develop elevated values in the TBU region, with the peak in $ f $ larger than the peak in $ c $ due to its smaller diffusivity. However, these peak values of osmolarity are not sufficient to stop the thinning, which progresses to TBU [36]. The evaporative flux $ J_e $ decreases at the thinnest $ h $ values due to the presence of the wetting term $ Ah^{-3} $. This prevents $ h $ from decreasing below about $ 0.25\ \mu $m to mimic the thickness of the glycocalyx [96]. If desired, we could compute solutions beyond the initial appearance of TBU, but we do not treat the spreading of TBU in this model because it is limited by the assumed evaporation distribution. While the TBU time is often used as a measure of tear film instability, the progression of TBU as measured by increasing area of TBU during the interblink is an important clinical aspect [97]. We note that when different evaporation models are used, the thinning may stop at much smaller values than we use here [36,98]. For a more complete discussion, see [91].
In Figure 4, we show TBU time as a function of spot size $ r_w $ for different peak thinning rates $ v_{max} $. The TBU time is defined as when the tear film thickness reaches twice the wetting equilibrium thickness $ h_{eq} = (A\alpha)^{1/3} $ where evaporation stops due to the balance with the wetting term. In Figure 4, we use $ h_{eq} = 0.25\; \mu $m, so the TBU time was taken to be when the tear film reached $ 0.5\; \mu $m. This size was chosen because the mathematical equilibrium is only slowly approached, which artificially extends the TBU time. In these results, osmosis is active with $ P_o = 12.1\; \mu $m/s permeability [18].
The results show that for relatively large spots with $ r_w > 0.8 $ and $ w_0\ge 10\mu $m/min, the TBU time is about one unit of time. This indicates that the dimensional time to get to TBU is about $ d/v_{max} $ which is the time for a flat film to evaporate away. If the peak thinning rate is small, such as $ v_{max} = 5\; \mu $m/min, then the thinning is slow and the TBU time increases because there is sufficient time for capillary flow and osmosis to supply fluid to the thinning region [91].
In the range of moderate spot sizes $ 0.4 \le r_w \le 0.8 $ and with $ v_{max} \ge 10\; \mu $m/min, the TBU time increases as $ r_w $ decreases. In this regime, capillary flow is increasing, which brings more fluid into the TBU region. However, osmolarity diffuses out rapidly enough that the peak value of $ c $ does not approach the high levels that occur for the flat film case [36], and thinning can proceed to TBU in similar times to the large spots.
Finally, in the range of small spots with $ r_w\le 0.4 $, the TBU time is very long because capillary flow into the TBU region (healing flow due to capillarity) is fast relative to the diffusion of solutes out of the TBU region. Increased advection causes osmolarity to reach values near the flat film case, and then osmosis and capillarity cooperate to slow thinning dramatically. However, the theoretical TBU times are long compared to in vivo results that we have fit with TBU models. This suggests that a different or additional mechanism may be present in small TBU spots that is not currently in the models.
Besides affecting TBU time via flow and osmolarity distributions, these dynamics have important consequences for imaging the tear film from the distributions of $ f $ in TBU.
Figure 5 shows results at different values of spot size and initial fluorescein concentration for evaporative TBU. It is clear that the best agreement between thickness and a function of intensity, for this mechanism, is for the lower left part of the array for evaporative TBU. This area is where $ f_0 $ is at the critical concentration or higher, and the spot sizes are near or larger than the length scale $ \ell $ (about 0.35 mm in this case). For further details about variation with these parameters, the interested reader is referred to [91].
In [55], the evaporation-driven model above was optimized to match intensity data extracted from in vivo experiments on subjects with a healthy tear film. The instances of TBU selected had simple spot (circular) and streak (linear) shapes for directly matching to the mathematical models. The TBU instances were also selected to be relatively slow to develop, as is expected for evaporation-driven TBU (e.g., [82]). The details of the instances used and the method for extracting the intensity data are given in their paper. The theoretical model is solved computationally, the intensity evaluated, and the process repeated until optimal parameters for matching the experimental intensity are found. Matlab's $\texttt{lsqnonlin}$ was used with either Levenburg-Marquardt or trust region reflective options (MathWorks, Natick, MA, USA). Detailed results of the fits can be found in [55].
Here we focus on the thinning rate found from the best fits and compare those results with experimentally determined thinning rates. Results are shown in Figure 6.
The fitting to the experimental intensity determines the peak evaporation rates $ v_{max} $ and the background evaporation rate $ v_{min} $. The resulting computed thinning rate (denoted $ dh/dt $) corresponds to the experimentally measured thinning rate from Nichols et al. [56]. We note that the experimental data is a point measurement from the central cornea and did not specifically target TBU instances. The results for thinning rate fall within the experimental range, and their maximum evaporation rate can fall above the maximum experimentially-measured thinning rate. Since flow within the film can reduce the thinning rate, particularly for evaporative thinning, this is an expected result.
Sharma and Ruckenstein developed mathematical models with van der Waals driven rupture using linear [99] and nonlinear [100,101] theories. They included a separate (fluid) mucus layer between the aqueous layer and the epithelium, and there was an insoluble surfactant on the interface between the two layers. These two-layer models were shear-dominated in both layers, and a surfactant transport governed the dynamics of the surfactant on the interface between them. The expressions for the fluid and surfactant fluxes are complicated; see those papers for details. The mucus layer was unstable to van der Waals forces in the model, and the models could all give reasonable TBU time ranges.
The two-layer film theory was generalized to include van der Waals forces in both mucus and aqueous layers, and surfactant transport on the aqueous-air interface [102,103]. The mucus layer was treated as a power-law fluid with a fit to experimental data for whole tears [104]. The aqueous layer was assumed to be Newtonian, and the lipid layer was simplified to the transport of an insoluble surfactant. They also found that the tear film could be unstable, with TBU driven by van der Waals forces. Thinner mucus layers led to reduced TBU times, while increased Marangoni effect led to increased TBU times. Related papers, which may apply to eyes via analogy with lung surfactants, include [105]. For a comprehensive review of related bilayer work, see [52].
Dey et al [106,107] posed a model with continuous viscosity variation through the depth of the tear film. It replaced separate aqeuous and mucin layers with variable mucin concentration and concentration-dependent viscosity. The concentration is governed by a transport equation. They computed linear and nonlinear stabilty results for the model, and found that for the fixed average velocity, a mucin distribution that caused higher viscosity near the ocular surface stabilized the tear film by increasing TBU time.
The TBU model by Peng et al. [36] adds a stationary non-polar lipid layer atop a dynamic aqueous layer. The evaporation rate is determined by two factors. It is inversely proportional to the lipid layer thickness and it incorporates the resistance to water vapor transport through the boundary layer in air just outside the tear film. That model was the first to find that the osmolarity rises dramatically in TBU, up to about 10 times the isotonic value. This result supports the notion that hyperosmolarity and evaporative TBU can contribute to ocular surface inflammation and dry eye. This paper also found that the osmolarity diffuses out of the TBU region, and this prevents osmosis from stopping TBU when there is spatial variation in the tear film quantities. The paper also investigated various conditions in the air outside the tear film.
Stapf et al. [98] extended [36] to include a mobile nonpolar lipid layer and an insoluble surfactant representing polar lipid at the lipid/aqueous interface. It was based on an earlier model of Bruna and Breward [74] that was used for one dimensional blinking. Both the nonpolar lipid and aqueous layer were assumed to be Newtonian; the lipid layer is dominated by extension, and the aqueous layer is dominated by shear. A nonlinear least squares fit was carried out for the dependence of the evaporation on the thickness the lipid layer from measured data [108]. They found that evaporative TBU can occur within this model with concentrated surfactant; the surfactant spread, causing a local thin region in the lipid layer, which then enhances evaporation there, leading to localized thinning of the tear film. The results presented in [98] were a TBU process dominated by evaporation; the early thinning of the lipid layer was very rapid, but did not drive TBU. Beyond that motion driven by the Marangoni effect, the lipid layer dynamics were rather slow compared to the other processes such as evaporation and aqueous flux due to capillarity. This result is consonant with simultaneous imaging experiments [82]: those observations contained cases where the lipid layer is essentially quiescent while thinning occurs due to evaporation through holes in it. The fixed lipid layer model of [36] was motivated by such observations, but the model with two mobile layers could capture the rapid development of the lipid imperfections that lead to increased evaporation processes that are slower.
The above two-layer models generate tangential flows that help drive thinning toward TBU. In this section, we turn to single-layer models which may have either an insoluble surfactant or a boundary condition playing the role of the lipid layer.
A model for TBU with an insoluble surfactant was developed by Siddique et al. (2015). In that model, the surfactant is assumed to be a monolayer. The evaporation rate linearly decreases with the surface concentration of a surfactant that represented the polar lipid layer of tear film. The nondimensional evaporation rate was given by
$ Je=1+δp−βΓˉK+h, $
|
(3.10) |
where $ \delta $, $ \beta $ and $ \bar{K} $ are constants. The results from that model showed initial rapid thinning occurs, but then the tear film could develop a competition between capillarity-driven healing flow and an outward flow from the Marangoni effect. The upshot of the model was that for the parameters tested, TBU occurred in too long a time for comparison with experiments of interest.
Results for a simplified model are now presented in this review. The evaporation rate is assumed independent of surfactant concentration, and entire tear film surface is mobile with an insoluble surfactant. Flow is driven by an initial "pile" of excess surfactant at the center of the film, and the surfactant is governed by a transport equation on the surface of the film. The low surface tension there causes a rapid divergent tangntial flow that thins the tear film in the center of the domain. Typically such thinning is much faster than the evaporative thinning [82].
Here, we give the equations in Cartesian coordinates for the streak model. Nondimensionally, the system is given by
$ ∂th+Je−Pc(c−1)+∂x(hˉu)=0, $
|
(3.11) |
$ p=−∂2xh−Ah−3, $
|
(3.12) |
$ h(∂tc+ˉu∂xc)=Pe−1c∂2xc+∂xh∂xc−Pc(c−1)c+Jec, $
|
(3.13) |
$ h(∂tf+ˉu∂xf)=Pe−1f∂2xf+∂xh∂xf−Pc(c−1)f+Jef. $
|
(3.14) |
$ ∂tΓ=Pe−1s∂2xΓ−∂x(usΓ), $
|
(3.15) |
here,
$ ˉu=−13h2∂xp−12Mh∂xΓ, $
|
(3.16) |
$ us=−12h2∂xp−Mh∂xΓ. $
|
(3.17) |
In addition to the dependent variables in the evaporative case, we must also find the insoluble surfactant concentration $ \Gamma(r, t) $. The surface velocity $ u_s $ also appears in the surfactant transport equation. The problem is solved on $ 0 < x < X_0 $, where $ X_0 $ is the outer end of the domain. Homogeneous Neumann boundary conditions are applied at $ x = 0 $ and $ x = X_0 $. The initial conditions are the same as for the spot case with $ r $ and $ R_0 $ replaced by $ x $ and $ X_0 $. We give results here for a uniform evaporation rate (constant in space and time). Nondimensional parameters arising from scaling the problem are given in Table 3.
Non-dimensional parameters with typical values | ||
Parameter | Description | Expression |
$ \epsilon $ | Aspect ratio | $ d / \ell $ |
$ S $ | Contribution of surface tension | $ (\sigma_0 \epsilon^3)/(\mu U) $ |
$ M $ | Contribution of Marangoni effect | $ \epsilon (\Delta \sigma)_0 \left[t_s^3/(\sigma_0 \mu^3 d^3) \right]^{1/4} $ |
$ A $ | Non-dimensional Hamaker constant | $ A^*/\left[\epsilon (\Delta \sigma)_0 d \ell \right] $ |
$ P_c $ | Permeability of cornea | $ (P_o V_w c_0)/(\epsilon U) $ |
Pe$ _f $ | Péclet number for FL diffusion | $ U \ell/D_f $ |
Pe$ _c $ | Péclet number for salt iron diffusion | $ U \ell / D_o $ |
Pe$ _s $ | Péclet number for surface diffusion | $ \epsilon (\Delta \sigma)_0 \ell/(\mu D_s) $ |
$ \phi $ | Non-dimensional Napierian extinction coefficient | $ \epsilon_f f_{\text{cr}} d$ |
The initial condition for $ \Gamma $ is a high concentration of unity near the origin and a low concentration outside $ X_I $, using the transition function $ B $. The initial condition uses the transition function $ B(x) $:
$ Γ(x,0)=1⋅[1−B(x)]+0.1⋅B(x), $
|
(3.18) |
where
$ B(x;XI,XW)=12+12tanh(x−XIXW). $
|
(3.19) |
Here, $ X_I $ defines the half-width of the region of higher concentration of lipid, and $ X_W = 0.1 $ is the (fixed) transition width.
Computational solution of this problem yields dynamics where the thickness decreases rapidly to a limited degree due to the outward divergent flow. That flow is driven by surfactant concentration gradients and the Marangoni effect, but the surfactant distribution rapidly relaxes toward a more uniform distribution, and this lessens outward tangential flow as time increases. This slows thinning and results in longer TBU times when compared to experiment or other TBU models to follow. We find that the thickness and intensity lag experimental results later in the trial; this model has difficulty keeping up with experimentally observed thinning. A version of these results for a representative set of dimensional parameter values are shown in Figure 7. In the upper left, we see that the thickness and intensity are quite similar in this case, with the intensity slightly lagging the thickness. The osmolarity and FL concentrations (upper right), are also quite similar. The pressure variation decreases with increasing time (middle left), and the surfactant distribution rapidly flattens (middle right). The depth-averaged and surface velocities (bottom row) are both primarily outward and their magnitude decreases steadily as the surface concentration of surfactant flattens.
Thus in both models described in this subsection, the insoluble monolayer with a linearized equation of state appears to relax too rapidly to drive the tear film to TBU. We now turn to a model which addresses this issue via a hypothesized "glob" of excess lipid that is tangentially immobile on the tear film surface.
Zhong et al. [84] developed a model that hypothesized that a localized region of excess lipid, a so-called lipid glob, could supply polar lipid (surfactant) for a period of time long enough to drive rapid TBU. Experimental evidence suggests that this occurs in some kinds of rapid TBU [82]. Their model has a region under the glob that is deformable, is tangentially immobile and has constant high concentration of surfactant. Outside the glob, the surface is mobile, and a standard surfactant transport equation conserves surfactant. No flux conditions apply at the center of the glob and at the edge of the film away from the glob.
These different conditions on the surface are combined by weighting them with a smooth blend function $ B(r) $ (for a spot) that has a narrow transition at the edge of the glob. For the streak geometry, the Cartesian model applies, and the glob edge is at the fixed location $ r = R_I $. The nondimensional blend function is given by
$ B(r;RI,RW)=0.5+0.5tanh(r−RIRW). $
|
(3.20) |
One may choose $ R_w $ to be small so that the transition width of $ B(R; R_I, R_W) $ is very narrow compared to the domain size. $ B(r; R_I, R_W) $ then approximates a step function which is zero to the left of $ R_I $ but one to the right. The equations on the two subdomains are then combined into one continuous equation over the whole domain.
Once the boundary conditions are blended using this function, then lubrication theory can be applied to the thin film throughout the domain to obtain simplified equations for the dependent variables. The reduced equations in this case govern the tear film thickness $ h(r, t) $, pressure $ p(r, t) $ the solute concentrations $ c(r, t) $ and $ f(r, t) $ and surfactant surface concentration $ \Gamma(r, t) $:
$ ∂th+1r∂r(rhˉu)=−Je+Jo, $
|
(3.21) |
$ p=−1r∂r(r∂rh)−Ah−3, $
|
(3.22) |
and
$ h(∂tm+ˉu∂rm)=Pe−1s1r∂r(rh∂rm)+[Je−Jo]m, m=c,f, $
|
(3.23) |
where the osmotic flux is given by
$ Jo=Pc(c−1), $
|
(3.24) |
and the depth-averaged velocity component along the film is given by
$ ˉu=−13∂rph2[B+(1−B)(14h)]−12hB∂rΓB+(1−B)h. $
|
(3.25) |
One also has a surfactant transport equation governing the surface concentration of $ \Gamma $,
$ ∂tΓ=[Pe _s^{-1} (1r∂r(r∂rΓ))−1r∂r(rΓur)]B, $
|
(3.26) |
where the surface velocity component along the film is given by
$ ur(r,h,t)=12h2∂rp+−h(∂rp)[(12h)(1−B)+B]−(∂rΓ)BB+(1−B)hh. $
|
(3.27) |
The boundary conditions are as follow. At the glob center, $ r = 0 $, we enforce symmetry with homogeneous Neumann conditions to obtain the following equations for use in numerical solution:
$ ∂th=−h∂rˉu−Je, 0=p+∂2rh+Ah−3, ∂rc=∂rf=0, and∂rΓ=0. $
|
(3.28) |
At $ r = R_O $, no-flux conditions require
$ ∂rh=∂rc=∂rf=0,∂rp=0,and∂rΓ=0. $
|
(3.29) |
Finally, we choose uniform initial conditions except for the surfactant concentration:
$ h(r,0)=c(r,0)=1,f(r,0)=f0,andΓ(r,0)=1⋅(1−B)+0.1⋅B. $
|
(3.30) |
We also have the consistent pressure initial condition $ p = -A $. For the initial $ \Gamma $, $ \Gamma = 1 $ in the glob and 0.1 outside the glob with the smooth, narrow transition between them. All dynamics are driven by the surfactant distribution and any nonzero evaporative flux $ J_e $. The variable step size and variable order solver ($\texttt{ode15 s}$) that we used for the spatially-discretized system can readily handle the computations on short time intervals that we find for this model.
Different choices for the evaporative flux $ J_e(x) $ were made in Zhong et al. [84,85]. We only give a sample of those results here using a constant evaporation rate; the reader is referred to the original papers for more complete results.
Figure 8 portrays a spot TBU result for our default parameters from Tables 1 and with no evaporation ($ J_e = 0 $). In this case, TBU occurs and is induced only by the different composition of lipid inside the glob. The locally-elevated surfactant concentration there spreads out immediately over the tear/air interface, and the surfactant concentration variation generates a significant shear stress, which in turn drives a strong divergent tangential flow. A significant depression of the tear/air interface forms at 0.15 s. In 0.59 s, the aqueous layer thins from 3.5 to 0.25 $ \mu $m, the thickness which we defined as TBU. This TBU time is near unity, which indicates that the time scales we chose work well. The Marangoni flow dominates the thinning process in this case, which leads to TBU in less than a second. In some in vivo observations with simultaneous fluorescence and lipid layer interferometry [82], the dark spot has already appeared by 0.14 s (practically instantaneously); the bright circular glob in the lipid layer spreads rapidly but its center remains relatively bright and thus relatively thick. This experimental observation corresponds well with the model.
Within an appropriate range of glob size, $ R_I $, TBU can happen in about a second or less, successfully capturing the short time scale in rapid TBU. Figure 8 shows that $ h $ thins to 0.25 $ \mu $m at $ 0.57 $ s, which is the threshold for TBU we defined for this model. At the same time, $ \Gamma $ spreads outside the glob (to the right), so the glob acts as a source that drives tangential flow away from the glob. Due to the short time and our assumed uniform evaporation profile, $ c $ and $ f $ increased only slightly under the glob. The solute concentration increases were quite small compared to the evaporative TBU cases discussed above.
Comparing the first and last columns of Figure 8 shows that $ I $ is a good approximation to $ h $ in this case. This is because the intensity change is all from the thickness change since $ f $ changes very little. Detailed parametric study showed that the best results for imaging the tear film with rapid TBU appeared to be around $ f_0 = 0.7 $ (or 0.14$ \% $) for this fast thinning case. This value of FL concentration is below the critical concentration [85,109]. For related cases, a value slightly lower could also work best. Thus, for rapid thinning, using a more dilute initial fluorescein concentration, together with the intensity itself as indicator of thickness, appeared to betst approximate thickness of the tear film. See [85] for a more detailed explanation.
Equations (3.25) and (3.27) imply that the Marangoni contribution to flow is directly proportional to $ \partial_r \Gamma $. As the surfactant spreads out from the glob, $ \partial_r \Gamma $ on the tear/air interface decreases, which results in a weaker Marangoni flow. Meanwhile, the contribution to aqueous flow from capillarity increases due to the deformation of the tear film. In some cases, the aqueous layer thins around the edge of the glob at short times, but then thickens at longer times. When that happens, the divergent flow from the Marangoni effect is insufficient to lead to TBU because the capillary contribution to flow dominates the dynamics at later times. This phenomenon matches our expectations for the model and some experimental observations where TBU does not occur.
Recently, Luke et al. [57] chose to explore parameters and scalings that intentionally combined both evaporation and tangential flow at comparable levels. The models are rescaled versions of the models developed by Zhong et al. [85] but then wrapped in optimization to determine parameter values that best match experimental observations in a least squares sense. The models were solved computationally, and the computed FL intensity was calculated from the resulting $ h $ and $ f $ distributions as functions of space and time. The 2-norm of the difference between computed ($ I_{th} $) and experimentally measured ($ I_{ex} $) FL intensities was minimized over the glob size ($ R_I $), Marangoni effect ($ \Delta \sigma_0 $), and evaporation parameters ($ J_e $). In terms of a continuous optimization, they solved
$ argminRI,Je,Δσ0||Ith−Iex||22, $
|
(3.31) |
where the argmin designates the minimization over the parameters. Eight instances of TBU were fit this way, and tangential flow was found to contribute in all cases they selected. In a small number of cases, evaporation only played a small role, while there was a fairly even mix in five cases. We consider this to be strong evidence that divergent tangential flow is important in TBU instances that are significantly under 10 s for TBU time.
The value of the osmolarity in TBU is of particular interest to clinicians due to its potential for causing pain and inflammation [110,111,112]. Results are shown for both evaporative and mixed-mechanism TBU in Figure 9.
This figure shows that rapid TBU can still reach significant osmolarity values even when evaporation is moderate. Those results argue against eliminating evaporation from TBU modeling except for the extremely rapid cases (under 1 s for TBU time). The figure also shows that the osmolarity in the fits does not continue to grow with time as the trial duration grows. The evaporative cases peaked in the 650-900 mOsM range, with a typical value around 700 mOsM. These results suggest that the high osmolarities obtained from ODE models for spatially uniform tear films [18,47] should be viewed as upper bounds that are not very tight when compared to PDE models. We also note that these fits are from a limited sample, and so the range in vivo may be larger than we have found here.
The tear film is a complex fluid structure that is highly dynamic [46]. There are recent studies that solve model problems for flows all over the exposed ocular surface. The problems use thin film approximations, even in the meniscus; they have also used simplified material properties. The models still capture a number of aspects of observed tear film dynamics in vivo including flow around the eyelid margins. However, the models do not yet closely link fine scale dynamics for phenomena such as tear break up to the overall flow or blinking dynamics. There is clearly room for progress in this area.
Blinking dynamics and tear film dynamics have been studied more extensively in one-dimensional domains. The situation is similar in that a number of aspects of the overall flows have been described and explained by the models. However, linking small scale dynamics of TBU, for example, have not been linked to the overall flow. And simplified material parameters have been used in these models as well, but more complex models have been studied that included two-layer films. There is room for including more realistic parameters here also.
TBU has been studied in local thin film models. The models have differed in complexity, but they have had some successes. Evaporative models have predicted high levels of osmolarity [36,91] that comport well with clinical estimates based on subject feedback [35]. More recently, parameters of TBU have been estimated by fitting theoretical models to in vivo experiments on normal subjest using FL imaging [55,57]. The thinning rates found by parameter estimation have agreed well with in vivo measurements [56] that did not specifically target TBU. Additional insights were gained about thinning rates versus evaporation rates. The relationship between the of TBU and osmolarity are complex, however. Very fast instances of TBU may have a small increase in osmolarity, and more slowly developing instances caused by evaporation may have quite high osmolarity. But there are intermediate cases where both outward tangential flow and evaporation contribute to TBU, and the resulting hyperosmolarity takes on an intermediate value. Developing fast, automated methods to fit the TBU data may lead to inferring the hyperosmolarity in TBU and to quantitative understanding of the hyperosmolarity in regions of TBU. These kinds of models offer the possibility of inferring the mechanism in specific instances of TBU as well; we believe that this close connection between theory and experiment is crucial to better understanding of the tear film in general, and of TBU and dry eye disease in particular.
To date, it has been easier to image the aqueous and lipid layers in the tear film [81]. This has contributed substantially to the modeling efforts described in this work. Imaging the glycocalyx and the ocular surface remains more challenging than the other parts of the tear film [2,113]. It may be that advances in technology for that purpose lead to advances in modeling and pinpointing TBU thought to be driven primarily by dewetting or other effects chiefly due to the epithelial surface [114]. These remain areas for advancement.
While mathematical modeling of tear film dynamics has made progress in explaining and predicting some aspects of the tear film, there remain many open problems and challenges. Better models in both 1D and 2D domains with blinking would be welcome, as well as the use of realistic material properties. We also believe that it is crucial to work closely with clinicians to better understand the tear film. Nothing is simple about the tear film, yet there is a need to simplify the mathematical descriptions of the tear film to make progress. Making appropriate modeling decisions based on insight from data and from clinicians is essential. Testing the results on data leads to better understanding for all involved.
This work was supported by National Science Foundation grant DMS 1909846. The content is solely the responsibility of the authors and does not necessarily represent the official views of the funding source. The authors gratefully acknowledge P.E. King-Smith for helpful conversations and sharing of data.
The authors declare that they have no conflicts of interest.
[1] | Costerton JW, Stewart PS, Greenberg EP (1999) Bacterial Biofilms: A Common Cause of Persistent Infections. Science 27: 1318-1322. |
[2] |
Davies D (2003) Understanding biofilm resistance to antibacterial agents. Natur Rev Drug Discov 2: 114-122. doi: 10.1038/nrd1008
![]() |
[3] | Sun G (2014) Prevention of hospital and community acquired infections by using antibacterial textiles and clothing. In: Polymeric materials with antimicrobial activity, Muñoz-Bonilla A, Cerrada ML, Fernández-García M., Eds., RSC Polymer Chemistry Series, Ch 6, 139-154. |
[4] |
Chen L, Bromberg L, Hatton TA, et al. (2008) Electrospun cellulose acetate fibers containing chlorhexidine as a bactericide. Polymer 49: 1266-1275. doi: 10.1016/j.polymer.2008.01.003
![]() |
[5] |
DeLeo FR, Chambers HF (2009) Reemergence of antibiotic-resistant Staphylococcus aureus in the genomics era. J Clin Invest 119: 2464-2474. doi: 10.1172/JCI38226
![]() |
[6] | Melaiye A, Sun Z, Hindi K, et al. (2005) Silver(I)−Imidazole Cyclophane gem-Diol Complexes Encapsulated by Electrospun Tecophilic Nanofibers: Formation of Nanosilver Particles and Antimicrobial Activity. J Am Chem Soc 8: 2285-2291. |
[7] |
Baley GJ, Peck GE, Banker GS (1977) Bactericidal properties of quaternary ammonium compounds in dispersed systems. J Pharm Sci 66: 696-699. doi: 10.1002/jps.2600660522
![]() |
[8] | Ortenzio LF, Stuart LS (1959) The behavior of chlorine-bearing organic compounds in the AOAC available chlorine gemicidal equivalent concentration test. J Assoc Off Ana Chem 42: 630-633. |
[9] |
Chopra I, Roberts M (2001) Tetracycline Antibiotics: Mode of Action, Applications, Molecular Biology, and Epidemiology of Bacterial Resistance. Microbiol Mol Biol Rev 65: 232-260. doi: 10.1128/MMBR.65.2.232-260.2001
![]() |
[10] |
Danese PN (2002) Antibiofilm Approaches: Prevention of Catheter Colonization. Chem Biol 9: 873-880. doi: 10.1016/S1074-5521(02)00192-8
![]() |
[11] |
Jarvis WR (1996) Selected aspects of the socioeconomic impact of nosocomial infections: morbidity, mortality, cost, and prevention. Infect Cont Hosp Ep 17: 552-557. doi: 10.2307/30141291
![]() |
[12] |
McGowan JE (2001) Economic impact of antimicrobial resistance. Emerg Infect Dis 7: 286-292. doi: 10.3201/eid0702.010228
![]() |
[13] | Reneker DH, Chun I (1996) Nanometre diameter fibres of polymer, produced by electrospinning, Nanotechnology 7: 216-223. |
[14] |
Li D, Xia Y (2004) Electrospinning of nanofibers: Reinventing the wheel? Adv Mater 16: 1151-1170. doi: 10.1002/adma.200400719
![]() |
[15] |
Deitzel, JM, Kleinmeyer J, Harris D et al. (2001) The effect of processing variables on the morphology of electrospun nanofibers and textiles. Polymer 42: 261-272. doi: 10.1016/S0032-3861(00)00250-0
![]() |
[16] |
Daels N, De Vrieze S, Sampers I, et al (2011) Potential of a functionalised nanofibre microfiltration membrane as an antibacterial water filter. Desalination 275: 285-290. doi: 10.1016/j.desal.2011.03.012
![]() |
[17] |
Zodrow K, Brunet L, Mahendra S (2009) Polysulfone ultrafiltration membranes impregnated with silver nanoparticles show imporved biofouling resistance and virus removal. Water Res 43: 715-723. doi: 10.1016/j.watres.2008.11.014
![]() |
[18] |
Payne RJH, Jansen VAA (2003) Pharmacokinetic principles of bacteriobacteriophage therapy. Clin Pharmacokinet 42: 315-325. doi: 10.2165/00003088-200342040-00002
![]() |
[19] | Gao Y, Truong YB, Zhu Y, et al. (2014) Electrospun antibacterial nanofibers: production, activity, and in vivo applications. J Appl Polym Sci 131: 40797. |
[20] | Wang S, Zheng F, Huang Y, et al. (2012) Encapsulation of amoxicillin within laponite-doped poly(lactic-co-glycolic acid) nanofibers: preparation, characterization, and antibacterial activity. Acs Appl Mater Interfaces 4: 6396-6401. |
[21] |
Özdemir I, Özcan EO, Günal S (2010) Synthesis and Antimicrobial Activity of Novel Ag-N-Heterocyclic Carbene Complexes. Molecules 15: 2499-2508. doi: 10.3390/molecules15042499
![]() |
[22] |
Clement JL, Jarret PS (1994) Antibacterial Silver. Met Based Drugs 1: 467-482. doi: 10.1155/MBD.1994.467
![]() |
[23] |
Tambe SM, Sampath L, Modak SM (2001) In-vitro evaluation of the risk of developing bacterial resistance to antiseptics and antibiotics used in medical devices. J Antimicrob Chemoth 47: 589-598. doi: 10.1093/jac/47.5.589
![]() |
[24] |
Liu JJ, Galettis P, Farr A, et al. (2008) In vitro antitumour and hepatotoxicity profiles of Au(I) and Ag(I) bidentate pyridyl phosphine complexes and relationships to cellular uptake. J Inorg Biochem 102: 303-310. doi: 10.1016/j.jinorgbio.2007.09.003
![]() |
[25] | Jeong L, Kim MH, Jung JY, et al. (2014) Effect of silk fibroin nanofibers containing silver sulfadiazine on wound healing. Int J Nanomed 9: 5277-5287. |
[26] |
Burd A, Kwok CH, Hung SC, et al. (2007) A comparative study of the cyto¬toxicity of silver-based dressings in monolayer cell, tissue explant, and animal model. Wound Repair Regen 15: 94-104. doi: 10.1111/j.1524-475X.2006.00190.x
![]() |
[27] |
Lok CN, Ho CM, Chen R, et al. (2006) Proteomic analysis of the mode of antibacterial action of silver nanoparticles. J Proteome Res 5: 916-924. doi: 10.1021/pr0504079
![]() |
[28] | Lansdown AB (2005) A guide to the properties and uses of silver dressings in wound care. Prof Nurse 20: 41-43. |
[29] |
Vepari C, Kaplan DL (2007) Silk as a biomaterial. Prog Polym Sci 32: 991-1007. doi: 10.1016/j.progpolymsci.2007.05.013
![]() |
[30] |
Kasoju N, Bora U (2012) Silk fibroin in tissue engineering. Adv Healthc Mater 1: 393-412. doi: 10.1002/adhm.201200097
![]() |
[31] |
Dong G, Xiao X, Liu X, et al. (2009) Functional Ag porous films prepared by electrospinning. J Appl Surf Sci 255: 7623-7626. doi: 10.1016/j.apsusc.2009.04.039
![]() |
[32] |
Xu X, Yang Q, Wang Y, el al. (2006) Biodegradable electrospun poly(l-lactide) fibers containing antibacterial silver nanoparticles. Eur Polym J 42: 2081-2087. doi: 10.1016/j.eurpolymj.2006.03.032
![]() |
[33] |
Son WK, Youk JH, Park WH (2006) Antimicrobial cellulose acetate nanofibers containing silver nanoparticles. Carbohyd Polym 65: 430-434. doi: 10.1016/j.carbpol.2006.01.037
![]() |
[34] |
Lala NL, Ramaseshan R, Bojun L, et al. (2007) Fabrication of nanofibers with antimicrobial functionality used as filters: protection against bacterial contaminants. Biotechnol Bioeng 97: 1357-1365. doi: 10.1002/bit.21351
![]() |
[35] | Yang QB, Li DM, Hong YL, el al. (2003) Preparation and characterization of a PAN nanofibre containing Ag nanoparticles via electrospinning Proceedings of the 2002 International Conference on Science and Technology of Synthetic Metals, Elsevier SA, Lausanne, Switzerland, 973-974. |
[36] |
Sheikh FA, Barakat NAM, Kanjwal MA, el al. (2009) Electrospun antimicrobial polyurethane nanofibers containing silver nanoparticles for biotechnological applications. Macromol Res 17: 688-696. doi: 10.1007/BF03218929
![]() |
[37] | Gliscinska E, Gutarowska B, Brycki B, el al. (2012) Electrospun Polyacrylonitrile Nanofibers Modified by Quaternary Ammonium Salts. J Appl Polym Sci 128: 767-775. |
[38] |
Yarin AL, Zussman E (2004) Upward needleless electrospinning of multiple nanofibers. Polymer 45: 2977-2980. doi: 10.1016/j.polymer.2004.02.066
![]() |
[39] |
Nomiya K, Tsuda K, Sudoh T, el al. (1997) Ag(I)-N bond-containing compound showing wide spectra in effective antimicrobial activities: Polymeric silver(I) imidazolate. J Inorg Biochem 68: 39-44. doi: 10.1016/S0162-0134(97)00006-8
![]() |
[40] | Nomiya K, Noguchi R, Oda M (2000) Synthesis and crystal structure of coinage metal(I) complexes with tetrazole (Htetz) and triphenylphosphine ligands, and their antimicrobial activities. A helical polymer of silver(I) complex [Ag(tetz)(PPh3)2]n and a monomeric gold(I) complex [Au(tetz)(PPh3)]. Inorg Chim Acta 298: 24-32. |
[41] |
Atiyeh BS, Ioannovich J, Al-Amm CA, el al. (2002) Management of acute and chronic open wounds: the importance of moist environment in opimal wound healing. Curr Pharm Biotechnol 3: 179-195. doi: 10.2174/1389201023378283
![]() |
[42] |
Akhavan O, Ghaderi E (2009) Enhancement of antibacterial properties of Ag nanorods by electric field. Sci Technol Adv Mat 10: 015003. doi: 10.1088/1468-6996/10/1/015003
![]() |
[43] |
Tijing LD, Amarjargal A, Jiang Z (2013) Antibacterial tourmaline nanoparticles/polyurethane hybrid mat decorated with silver nanoparticles prepared by electrospinning and UV photoreduction. Curr Appl Phys 13: 205-210. doi: 10.1016/j.cap.2012.07.011
![]() |
[44] |
Raffi M, Mehrwan S, Bhatti TM, et al. (2010) Investigations into the antibacterial behavior of copper nanoparticles against Escherichia coli. Ann Microbiol 60: 75-80. doi: 10.1007/s13213-010-0015-6
![]() |
[45] |
Zhong W, Zishena W, Zhenhuana Y, et al. (1994) Synthesis, characterization and antifungal activity of copper (II), zinc (II), cobalt (II) and nickel (II) complexes derived from 2-chlorobenzaldehyde and glycine. Synth React Inorg Met-Org Chem 24: 1453-1460. doi: 10.1080/00945719408002572
![]() |
[46] |
Sheikh FA, Kanjwal MA, Saran S, et al. (2011) Polyurethane nanofibers containing copper nanoparticles as future materials. Appl Surf Sci 257: 3020-3026. doi: 10.1016/j.apsusc.2010.10.110
![]() |
[47] |
Hwang SH, Song J, Jung Y, et al. (2011) Electrospun ZnO/TiO2 composite nanofibers as a bactericidal agent. J Chem Commun 47: 9164-9166. doi: 10.1039/c1cc12872h
![]() |
[48] |
Lee S (2009) Multifunctionality of layered fabric systems based on electrospun polyurethane/zinc oxide nanocomposite fibers. J Appl Polym Sci 114: 3652-3658. doi: 10.1002/app.30778
![]() |
[49] |
Lee K, Lee S (2012) Multifunctionality of poly(vinyl alcohol) nanofiber webs containing titanium dioxide. J Appl Polym Sci 124: 4038-4046. doi: 10.1002/app.34929
![]() |
[50] | Yu W, Lan CH, Wang SJ, et al. (2010), Influence of Zinc oxide nanoparticles on the crystallization behavior of electrospun poly(3-hydroxybutyrate-co-3-hydroxyvalerate) nanofibers. Polymer 51: 2403. |
[51] |
Drew C, Liu X, Ziegler D, et al. (2003) Metal oxide-coated polymer nanofibers. Nano Lett 3: 143-147. doi: 10.1021/nl025850m
![]() |
[52] |
Horzum N, Mari M, Wagner M, et al. (2015) Controlled surface mineralization of metal oxides on nanofibers. RSC Adv 5: 37340-37345. doi: 10.1039/C5RA02140E
![]() |
[53] |
Klabunde KJ, Stark J, Koper O, et al. (1996) Nanocrystals as stoichiometric reagents with unique surface chemistry. J Phys Chem 100: 12142-12153. doi: 10.1021/jp960224x
![]() |
[54] |
Stoimenov PK, Klinger RL, Marchin GL (2002) Metal Oxide Nanoparticles as Bactericidal Agents. Langmuir 18: 6679-6686. doi: 10.1021/la0202374
![]() |
[55] | Koper O, Klabunde J, Marchin G, et al. (2012) Nanoscale Powders and Formulations with Biocidal Activity Toward Spores and Vegetative Cells of Bacillus Species, Viruses, and Toxins. Curr Microbiol 44: 49-55. |
[56] | McDonnell G, Russell AD (1999) Antiseptics and disinfectants: activity, action, and resistance. Clin Microbiol Rev 12: 147-179. |
[57] | Kügler R, Bouloussa O, Rondelez F (2005) Evidence of a charge-density threshold for optimum Synthesis and antimicrobial activity of modified poly(glycidyl methacrylate-co-2-hydroxy -ethyl methacrylate) derivatives with quaternary ammonium and phosphonium salts. J Polym Sci Pol Chem 40: 2384-2393. |
[58] | Kenawy ER, Abdel-Hay FI, El-Shanshoury A, et al. (2002) Biologically active polymers. V. Synthesis and antimicrobial activity of modified poly(glycidyl methacrylate-co-2-hydroxyethyl methacrylate) derivatives with quaternary ammonium and phosphonium salts. J Polym Sci Part A: Polym Chem 40: 2384-2393. |
[59] |
Kenawy ER, Mahmoud YAG (2003) Biological active polymers. Macromol Biosci 3: 107-116. doi: 10.1002/mabi.200390016
![]() |
[60] |
Lundin JG, Coneski PN, Fulmer PA, et al. (2014) Relationship between surface concentration of amphiphilic quaternary ammonium biocides in electrospun polymer fibers and biocidal activity. React Funct Polym 77: 39-46. doi: 10.1016/j.reactfunctpolym.2014.02.004
![]() |
[61] |
Arumugam GK, Khan S, Heiden PA (2009) Comparison of the Effects of an Ionic Liquid and Other Salts on the Properties of Electrospun Fibers, 2-Poly(vinyl alcohol). Macromol Mater Eng 294: 45-53. doi: 10.1002/mame.200800199
![]() |
[62] |
You Y, Lee SJ, Min BM, et al. (2006) Effect of solution properties on nanofibrous structure of electrospun poly(lactic-co-glycolic acid). J Appl Polym Sci 99: 1214-1221. doi: 10.1002/app.22602
![]() |
[63] |
Tong HW, Wang M (2011) Electrospinning of poly(hydroxybutyrate-co-hydroxyvalerate) fibrous scaffolds for tissue engineering applications: effects of electrospinning parameters and solution properties. J Macromol Sci 50: 1535-1558. doi: 10.1080/00222348.2010.541008
![]() |
[64] |
Park JA, Kim SB (2015) Preparation and characterization of antimicrobial electrospun poly(vinyl alcohol) nanofibers containing benzyl triethylammonium chloride. React Funct Polym 93: 30-37. doi: 10.1016/j.reactfunctpolym.2015.05.008
![]() |
[65] |
Kim SJ, Nam YS, Rhee DM, et al. (2007) Preparation and characterization of antimicrobial polycarbonate nanofibrous membrane. Eur Polym J 43: 3146-3152. doi: 10.1016/j.eurpolymj.2007.04.046
![]() |
[66] |
Nicosia A, Gieparda W, Foksowicz-Flaczyk J, et al. (2015) Air filtration and antimicrobial capabilities of electrospun PLA/PHB containing ionic liquid. Sep Purif Technol 154: 154-160. doi: 10.1016/j.seppur.2015.09.037
![]() |
[67] |
Buruiana EC, Buruiana T (2002) Recent developments in polyurethane cationomers. Photoisomerization reactions in azoaromatic polycations. J Photoch Photobiol A: Chem 151: 237-252. doi: 10.1016/S1010-6030(02)00181-8
![]() |
[68] |
Uykun N, Ergal I, Kurt H, et al. (2014) Electrospun antibacterial nanofibrous polyvinylpyrrolidone/ cetyltrimethylammonium bromide membranes for biomedical applications. J Bioact Compat Pol 29: 382-397. doi: 10.1177/0883911514535153
![]() |
[69] |
Qun XL, Fanf Y, Shan Y, et al. (2010) Antibacterial Nanofibers of Self-quaternized Block Copolymers of 4-Vinyl Pyridine and Pentachlorophenyl Acrylate. High Perform Polym 22: 359-376. doi: 10.1177/0954008309104776
![]() |
[70] |
Cashion MP, Li X, Geng Y, at al. (2010) Gemini Surfactant Electrospun Membranes. Langmuir 26: 678-683. doi: 10.1021/la902287b
![]() |
[71] |
Singh G, Bittner AM, Loscher S, et al. (2008) Electrospinning of diphenylalanine nanotube. Adv Mater 20: 2332-2336. doi: 10.1002/adma.200702802
![]() |
[72] |
Worley SD, Williams DE (1988) Halamine Water Disinfectants. Crit Rev Environ Control 18: 133-175. doi: 10.1080/10643388809388345
![]() |
[73] |
Song J, Jang J (2014) Antimicrobial polymer nanostructures: Synthetic route, mechanism of action and perspective. Adv Colloid Interface Sci 203: 37-50. doi: 10.1016/j.cis.2013.11.007
![]() |
[74] |
Cerkez I, Worley SD, Broughton RM, et al. Huang (2013) Rechargeable antimicrobial coatings for poly(lactic acid) nonwoven fabrics. Polymer 54: 536-541. doi: 10.1016/j.polymer.2012.11.049
![]() |
[75] |
Chen Z, Sun Y (2005) N-Chloro-Hindered Amines as Multifunctional Polymer Additives. Macromolecules 38: 8116-8119. doi: 10.1021/ma050874b
![]() |
[76] | Chen Z, Sun Y (2005) Antimicrobial Functions of N-Chloro-Hindered Amines. Polym Preprint 46: 835-836. |
[77] |
Dong A, Zhang Q, Wang T, et al. (2010) Immobilization of cyclic N-halamine on polystyrene-functionalized silica nanoparticles: synthesis, characterization, and biocidal activity. J Phys Chem C 114: 17298-17303. doi: 10.1021/jp104083h
![]() |
[78] |
Kang J, Han J, Gao Y, et al. (2015) Unexpected Enhancement in Antibacterial Activity of N-Halamine Polymers from Spheres to Fibers. ACS Appl Mater Interfaces 7: 17516-17526. doi: 10.1021/acsami.5b05429
![]() |
[79] |
Sun X, Zhang L, Cao Z, et al. (2010) Electrospun Composite Nanofiber Fabrics Containing Uniformly Dispersed Antimicrobial Agents As an Innovative Type of Polymeric Materials with Superior Antimicrobial Efficacy. ACS Appl Mater Interfaces 2: 952-956. doi: 10.1021/am100018k
![]() |
[80] |
Tan K, Obendorf SK (2007) Fabrication and evaluation of electrospun antimicrobial nylon 6 membranes. J Membrane Sci 305: 287-298. doi: 10.1016/j.memsci.2007.08.015
![]() |
[81] | Dickerson MB, Sierra AA, Bedford NM, et al. (2013) Keratin-based antimicrobial textiles, films, and nanofibers. J Mater Chem B 1: 5505-5514. |
[82] |
Fan X, Jiang Q, Sun Z, et al. (2015) Preparation and Characterization of Electrospun Antimicrobial Fibrous Membranes Based on Polyhydroxybutyrate (PHB). Fiber Polym 16: 1751-1758. doi: 10.1007/s12221-015-5108-1
![]() |
[83] |
Ignatova M, Rashkov I, Manolova N (2013) Drug-loaded electrospun materials in wound-dressing applications and in local cancer treatment. Expert Opin Drug Del 10: 469-483. doi: 10.1517/17425247.2013.758103
![]() |
[84] |
Kim K, Luu YK, Chang C et al. (2004) Incorporation and controlled release of a hydrophilic antibiotic using poly(lactide-co-glycolide)-based electrospun nanofibrous scaffolds. J Control Release 98: 47-56. doi: 10.1016/j.jconrel.2004.04.009
![]() |
[85] |
Yoo HS, Kim TG, Park TG (2009) Surface-functionalized electrospun nanofibers for tissue engineering and drug delivery. Adv Drug Deliver Rev 61: 1033-1042. doi: 10.1016/j.addr.2009.07.007
![]() |
[86] |
Zong XH, Li S, Chen E, at al. (2004) Prevention of postsurgery-induced abdominal adhesions by electrospun bioabsorbable nanofibrous poly(lactide-co-glycolide)-based membranes. Ann Surg 240: 910-915. doi: 10.1097/01.sla.0000143302.48223.7e
![]() |
[87] |
Hong Y, Fujimoto K, Hashizume R, et al. (2008) Generating elastic, biodegradable polyurethane/poly(lactide-co-glycolide) fibrous sheets with controlled antibiotic release via two-stream electrospinning. Biomacromolecules 9: 1200-1207. doi: 10.1021/bm701201w
![]() |
[88] |
Thakur RA, Florek CA, Kohn J, et al. (2008) Electrospun nanofibrous polymeric scaffold with targeted drug release profiles for potential application as wound dressing. Int J Pharm 364: 87-93. doi: 10.1016/j.ijpharm.2008.07.033
![]() |
[89] |
Gilchrist SE, Lange D, Letchford DK, et al. (2013) Fusidic acid and rifampicin co-loaded PLGA nanofibers for the prevention of orthopedic implant associated infections. J Control Release 170: 64-73. doi: 10.1016/j.jconrel.2013.04.012
![]() |
[90] |
Qi R, Guo R, Shen M, et al. (2010) Electrospun poly(lactic-co-glycolic acid)/halloysite nanotube composite nanofibers for drug encapsulation and sustained release. J Mater Chem 20: 10622-10629. doi: 10.1039/c0jm01328e
![]() |
[91] |
Moghe AK, Gupta BS (2008) Co‐axial Electrospinning for Nanofiber Structures: Preparation and Applications. Polym Rev 48: 353-377. doi: 10.1080/15583720802022257
![]() |
[92] |
Viseras C, Cerezo P, Sanchez R, et al. (2010) Current challenges in clay minerals for drug delivery. Appl Clay Sci 48: 291-295. doi: 10.1016/j.clay.2010.01.007
![]() |
[93] |
McMurry LM, Oethinger M, Levy SB (1998) Triclosan targets lipid synthesis. Nature 394: 531-532. doi: 10.1038/28970
![]() |
[94] | Green JBD, Fulghum T, Nordhaus MA (2011, Immobilized antimicrobial agents: a critical perspective, In: Science Against Microbial Pathogens: Communicating Current Research and Technological Advances, Formatex Microbiology Books Series, Méndez-Vila E, Ed., Badajoz (Spain). |
[95] |
Kaehn K (2010) Polihexanide: A Safe and Highly Effective Biocide. Skin Pharmacol Physiol 23: 7-16. doi: 10.1159/000318237
![]() |
[96] |
Gilbert P, Pemberton D, Wilkinson DE (1990) Synergism within polyhexamethylene biguanide biocide formulations. J Appl Bacteriol 69: 593-598. doi: 10.1111/j.1365-2672.1990.tb01553.x
![]() |
[97] |
del Valle LJ, Camps R, Díaz A, et al. (2011) Electrospinning of polylactide and polycaprolactone mixtures for preparation of materials with tunable drug release properties. J Polym Res 18: 1903-1917. doi: 10.1007/s10965-011-9597-3
![]() |
[98] |
del Valle LJ, Díaz A, Royo M, et al. (2012) Biodegradable polyesters reinforced with triclosan loaded polylactide micro/nanofibers: Properties, release and biocompatibility. Express Polym Lett 6: 266-282. doi: 10.3144/expresspolymlett.2012.30
![]() |
[99] |
Llorens E, Bellmunt S, del Valle LJ, et al. (2014) Scaffolds constituted by mixed polylactide and poly(ethylene glycol) electrospun microfibers. J Polym Res 21: 603. doi: 10.1007/s10965-014-0603-4
![]() |
[100] |
Llorens E, Ibañez H, del Valle LJ, et al. (2015) Biocompatibility and drug release behavior of scaffolds prepared by coaxial electrospinning of poly(butylene succinate) and polyethylene glycol. Mater Sci Eng C 49: 472-484. doi: 10.1016/j.msec.2015.01.039
![]() |
[101] |
Llorens E, del Valle LJ, Puiggalí J (2015) Electrospun scaffolds of polylactide with a different enantiomeric content and loaded with anti-inflammatory and antibacterial drugs. Macromol Res 23: 636-648. doi: 10.1007/s13233-015-3082-5
![]() |
[102] |
Veiga M, Merino M, Cirri M. et al. (2005) Comparative study on triclosan interactions in solution and in the solid state with natural and chemically modified cyclodextrins. J Incl Phenom Macrocycl Chem 53: 77-83. doi: 10.1007/s10847-005-1047-6
![]() |
[103] |
Kayaci F, Umu OCO, Tekinay T, et al. (2013) Antibacterial electrospun Poly(lactic acid) (PLA) nanofibrous webs incorporating triclosan/cyclodextrin inclusion Complexes. J Agric Food Chem 61: 3901-3908. doi: 10.1021/jf400440b
![]() |
[104] |
Celebioglu A, Umu OCO, Tekinay T, et al. (2014) Antibacterial electrospun nanofibers from triclosan/cyclodextrin inclusion complexes. Colloid Surface B 116: 612-619. doi: 10.1016/j.colsurfb.2013.10.029
![]() |
[105] |
Chen L, Bromberg L, Hatton TA (2008) Electrospun cellulose acetate fibers containing chlorhexidine as a bactericide. Polymer 49: 1266-1275. doi: 10.1016/j.polymer.2008.01.003
![]() |
[106] |
Ignatova M, Stoilova O, N Manolova N, et al. (2010) Electrospun mats from styrene/maleic anhydride copolymers: Modification with amines and assessment of antimicrobial activity. Macromol Biosci 10: 944-954. doi: 10.1002/mabi.200900433
![]() |
[107] |
Fernandes JG, Correia DM, Botelho G, et al. (2014) PHB-PEO electrospun fiber membranes containing chlorhexidine for drug delivery applications. Polym Test 34: 64-71. doi: 10.1016/j.polymertesting.2013.12.007
![]() |
[108] | del Valle LJ, Roa M, Díaz A, et al. (2012) Electrospun nanofibers of a degradable poly(ester amide). Scaffolds loaded with antimicrobial agents. J Polym Res 19: 9792. |
[109] |
Díaz A, Katsarava R, Puiggalí J (2014) Synthesis, properties and applications of biodegradable polymers derived from diols and dicarboxylic Acids: From Polyesters to poly(ester amide)s. Int J Mol Sci 15: 7064-7123. doi: 10.3390/ijms15057064
![]() |
[110] | Rodríguez-Galán A, Franco L, Puiggalí J (2011) Degradable poly(ester amide)s for biomedical applications. Polymers 3: 65-99. |
[111] |
Murase SK, del Valle LJ, Kobauri S, et al. (2015) Electrospun fibrous mats from a L-phenylalanine based poly(ester amide): Drug delivery and accelerated degradation by loading enzymes. Polym Degrad Stab 119: 275-287. doi: 10.1016/j.polymdegradstab.2015.05.018
![]() |
[112] |
Díaz A, del Valle LJ, Tugushi D, et al. (2015) New poly(ester urea) derived from L-leucine: Electrospun scaffolds loaded with antibacterial drugs and enzymes. Mater Sci Eng C 46: 450-462. doi: 10.1016/j.msec.2014.10.055
![]() |
[113] | Saha K, Butola BS, Joshi M (2014) Drug-loaded polyurethane/clay nanocomposite nanofibers for topical drug-delivery application. J Appl Polym Sci 131: 40230. |
[114] |
Llorens E, Calderón S, del Valle LJ, et al. (2015) Polybiguanide (PHMB) loaded in PLA scaffolds displaying high hydrophobic, biocompatibility and antibacterial properties. Mater Sci Eng C 50: 74-84. doi: 10.1016/j.msec.2015.01.100
![]() |
[115] | Dash M, Chiellini F, Ottenbrite RM, et al. (2007) Chitosan-a versatile semi-synthetic polymer in biomedical applications. Prog Polym Sci 36: 981-1014. |
[116] | Fischer TH, Bode AP, Demcheva M, et al. (2007) Hemostatic properties of glucosamine-based materials. J Biomed Mater Res A 80: 167-174. |
[117] | Muzzarelli RAA, Belcher R, Freisers H, In: Chitosan in natural chelating polymer; alginic acid, chitin and chitosan, Pergamon Press, Oxford, Pergamon Press, Oxford, 1973, 144-176. |
[118] |
Geng X, Kwon OH, Jang J (2005) Electrospinning of chitosan dissolved in concentrated acetic acid solution. Biomaterials 26: 5427-5432. doi: 10.1016/j.biomaterials.2005.01.066
![]() |
[119] | Spasova M, Manolova N, Paneva D, et al (2004) Preparation of chitosan-containing nanofibres by electrospinning of chitosan/poly(ethylene oxide) blend solutions. e-Polymers 4: 624-635. |
[120] |
Bhattaraia B, Edmondson D, Veiseh O, et al. (2005) Electrospun chitosan-based nanofibers and their cellular compatibility. Biomaterials 26: 6176-6184. doi: 10.1016/j.biomaterials.2005.03.027
![]() |
[121] |
Ignatova M, Starbova K, Markova N, et al. (2006) Electrospun nano-fibre mats with antibacterial properties from quaternised chitosan and poly(vinyl alcohol). Carbohyd Res 341: 2098-2107. doi: 10.1016/j.carres.2006.05.006
![]() |
[122] |
Min BM, Lee SW, Lim JN, et al. (2004) Chitin and chitosan nanofibers: Electrospinning of chitin and deacetylation of chitin nanofibers. Polymer 45: 7137-7142. doi: 10.1016/j.polymer.2004.08.048
![]() |
[123] |
Ohkawa K, Cha D, Kim H, et al. (2004) Electrospinning of chitosan. Macromol Rapid Comm 25: 1600-1605. doi: 10.1002/marc.200400253
![]() |
[124] |
Neamnark A, Rujiravaniti R, Supaphol P (2006) Electrospinning of hexanoyl chitosan. Carbohydrate. Polymers 66: 298-305. doi: 10.1016/j.carbpol.2006.03.015
![]() |
[125] |
Homayoni H, Ravandi SH, Valizadeh M (2009) Electrospinning of chitosan nanofibers: Processing optimization. Carbohyd Polym 77: 656-661. doi: 10.1016/j.carbpol.2009.02.008
![]() |
[126] |
Zong X, Kim K, Fang D, et al. (2002) Structure and process relationship on electrospun bioabsorbable nanofiber membranes. Polymer 43: 4403-4412. doi: 10.1016/S0032-3861(02)00275-6
![]() |
[127] |
Au HT, Pham LN, Vu THT, et al. (2012) Fabrication of an antibacterial non-woven mat of a poly(lactic acid)/chitosan blend by electrospinning. Macromol Res 20: 51-58. doi: 10.1007/s13233-012-0010-9
![]() |
[128] |
Shalumon KT, Anulekha KH, Girish CM, et al. (2010) Single step electrospinning of chitosan/poly(caprolactone) nanofibers using formic acid/acetone solvent mixture. Carbohyd Polym 80: 413-417. doi: 10.1016/j.carbpol.2009.11.039
![]() |
[129] |
Cooper A, Oldinski R, Ma H, et al. (2013) Chitosan-based nanofibrous membranes for antibacterial filter applications. Carbohyd Polym 92: 254-259. doi: 10.1016/j.carbpol.2012.08.114
![]() |
[130] |
Zheng H, Du Y, Yu J, et al. (2001) Preparation and characterization of chitosan/poly(vinyl alcohol) blend fibers. J Appl Polym Sci 80: 2558-2565. doi: 10.1002/app.1365
![]() |
[131] |
Chuang WY, Young TH, Yao CH, et al. (1999) Properties of the poly(vinyl alcohol)/chitosan blend and its effect on the culture of fibroblast in vitro. Biomaterials 20: 1479-1487. doi: 10.1016/S0142-9612(99)00054-X
![]() |
[132] |
Jia YT, Gong J, Gu XH, et al. (2007) Fabrication and characterization of poly(vinyl alcohol)/chitosan blend nanofibers produced by electrospinning method. Carbohyd Polym 67: 403-409. doi: 10.1016/j.carbpol.2006.06.010
![]() |
[133] |
Zhang H, Li S, White CJB, et al. (2009) Studies on electrospun nylon-6/chitosan complex nanofiber interactions. Electrochim Acta 54: 5739-5745. doi: 10.1016/j.electacta.2009.05.021
![]() |
[134] |
Lin SJ, Hsiao WU, Jee SH, et al. (2006) Study on the effects of nylon–chitosan-blended membranes on the spheroid-forming activity of human melanocytes. Biomaterials 27: 5079-5088. doi: 10.1016/j.biomaterials.2006.05.035
![]() |
[135] |
Ma Z, Kotaki M, Yong T, et al. (2005) Surface engineering of electrospun polyethylene terephthalate (PET) nanofibers towards development of a new material for blood vessel engineering. Biomaterials 26: 2527-2536. doi: 10.1016/j.biomaterials.2004.07.026
![]() |
[136] |
Jung KH, Huh MW, Meng W, et al. (2007) Preparation and antibacterial activity of PET/chitosan nanofibrous mats using an electrospinning technique. J Appl Polym Sci 105: 2816-2823. doi: 10.1002/app.25594
![]() |
[137] |
An J, Zhang H, Zhang JT, et al. (2009) Preparation and antibacterial activity of electrospun chitosan/poly(ethylene oxide) membranes containing silver nanoparticles. Colloid Polym Sci 287: 1425-1434. doi: 10.1007/s00396-009-2108-y
![]() |
[138] |
Zhuang XP, Cheng BW, Kang WM, et al. (2010) Electrospun chitosan/gelatin nanofibers containing silver nanoparticles. Carbohyd Polym 82: 524-527. doi: 10.1016/j.carbpol.2010.04.085
![]() |
[139] |
Cai ZX, Mo XM, Zhang KH, et al. (2010) Fabrication of chitosan/silk fibroin composite nanofibers for wound dressing applications. Int J Mol Sci 11: 3529-3539. doi: 10.3390/ijms11093529
![]() |
[140] |
Alipour SM, Nouri M, Mokhtari J, et al. (2009) Electrospinning of poly(vinyl alcohol)–water-soluble quaternized chitosan derivative blend. Carbohyd Res 344: 2496-2501. doi: 10.1016/j.carres.2009.10.004
![]() |
[141] |
Fu GD, Yao F, Li Z, et al. (2008) Solvent-resistant antibacterial microfibers of self-quaternized block copolymers from atom transfer radical polymerization and electrospinning. J Mater Chem 18: 859-867. doi: 10.1039/b716127a
![]() |
[142] |
Parfitt T (2005) Georgia: an unlikely stronghold for bacteriophage therapy. The Lancet 365: 2166-2167. doi: 10.1016/S0140-6736(05)66759-1
![]() |
[143] | Sulakvelidze A, Kutter E (2005) Bacteriophage therapy in humans, In: Bacteriophages: Biology and Applications, Kutter E, Sulakvelidze A, Eds., CRC Press, Boca Raton, FL, 381-436. |
[144] |
Kutter E, De Vos D, Gvasalia G, et al. (2010) Phage therapy in clinical practice: treatment of human infections. Curr Pharm Biotechnol 11: 69-86. doi: 10.2174/138920110790725401
![]() |
[145] |
Abedon S, Kuhl S, Blasdel B, et al. (2011) Phage treatment of human infections. Bacteriophage 1: 66-85. doi: 10.4161/bact.1.2.15845
![]() |
[146] | FDA (2006) FDA approval of Listeria-specific bacteriophage preparation on ready-to-eat (RTE) meat and poultry products. FDA, Washington, DC. |
[147] |
Frykberg RG, Armstrong DG, Giurini J, et al. (2000) Diabetic foot disorders: a clinical practice guideline. J Foot Ankle Surg 39: S1-S60. doi: 10.1016/S1067-2516(00)80057-5
![]() |
[148] |
Anany H, Chen W, Pelton R (2011) Biocontrol of Listeria monocytogenes and Escherichia coli O157:H7 in Meat by Using Phages Immobilized on Modified Cellulose Membranes. Appl Environ Microb 77: 6379-6387. doi: 10.1128/AEM.05493-11
![]() |
[149] |
Markoishvili K, Tsitlanadze G, Katsarava R, et al. (2002) A novel sustained-release matrix based on biodegradable poly(ester amide)s and impregnated with bacteriophages and an antibiotic shows promise in management of infected venous stasis ulcers and other poorly healing wounds. Int J Dermatol 41: 453-458. doi: 10.1046/j.1365-4362.2002.01451.x
![]() |
[150] | Katsarava R, Alavidze Z (2004) Polymer Blends as Biodegradable Matrices for Preparing Biocomposites, US Patent 6,703,040 (Assigned to Intralytix, Inc.). |
[151] | Katsarava R, Gomurashvili Z (2011) Biodegradable Polymers Composed of Naturally Occurring α-Amino Acids, In: Handbook of Biodegradable Polymers—Isolation, Synthesis, Characterization and Applications. Lendlein A, Sisson A, Eds., Wiley-VCH, Verlag GmbH & Co KGaA, Ch. 5. |
[152] |
Jikia D, Chkhaidze N, Imedashvili E, et al. (2005) The use of a novel biodegradable preparation capable of the sustained release of bacteriophages and ciprofloxacin, in the complex treatment of multidrug - resistant Staphylococcus aureus-infected local radiation injuries caused by exposure to Sr90. Clin Exp Dermatol 30: 23-26. doi: 10.1111/j.1365-2230.2004.01600.x
![]() |
[153] |
Puapermpoonsiri U, Spencer J, van der Walle CF (2009) A freeze-dried formulation of bacteriophage encapsulated in biodegradable microspheres. Eur J Pharm Biopharm 72: 26-33. doi: 10.1016/j.ejpb.2008.12.001
![]() |
[154] |
Gervais L, Gel M, Allain B, et al. (2007) Immobilization of biotinylated bacteriophages on biosensor surfaces. Sensor Actuat B-Chem 125: 615-621. doi: 10.1016/j.snb.2007.03.007
![]() |
[155] |
Nanduri V, Balasubramanian S, Sista S, et al. (2007) Highly sensitive phage-based biosensor for the detection of β-galactosidase. Anal Chim Acta 589: 166-172. doi: 10.1016/j.aca.2007.02.071
![]() |
[156] | Tawil N, Sacher E, Rioux D, et al. (2015) Surface Chemistry of Bacteriophage and Laser Ablated Nanoparticle Complexes for Pathogen Detection. J Phys Chem C 119: 14375-14382. |
[157] | Sun W, Brovko L, Griffiths M (2001) Food-borne pathogens. Use of bioluminescent salmonella for assessing the efficiency of Constructed phage-based biosorbent. J Ind Microbiol Biot 27: 126-128. |
[158] |
Cademartiri R, Anany H, Gross I, et al. (2010). Immobilization of bacteriophages on modified silica particles. Biomaterials 31: 1904-1910. doi: 10.1016/j.biomaterials.2009.11.029
![]() |
[159] |
Bennett AR, Davids FGC, Vlahodimou S, et al. (1997) The use of bacteriophage-based systems for the separation and concentration of Salmonella. J Appl Microbiol 83: 259-265. doi: 10.1046/j.1365-2672.1997.00257.x
![]() |
[160] |
Pearson HA, Sahukhal GS, Elasri MO, et al. (2013) Phage-Bacterium War on Polymeric Surfaces: Can Surface-Anchored Bacteriophages Eliminate Microbial Infections? Biomacromolecules 14: 1257-1261. doi: 10.1021/bm400290u
![]() |
[161] |
Dai M, Senecal A, Nugen SR (2014) Electrospun water-soluble polymer nanofibers for the dehydration and storage of sensitive reagents. Nanotechnology 25: 225101. doi: 10.1088/0957-4484/25/22/225101
![]() |
[162] | Korehei R, Kadla J (2011) Incorporating pre-encapsulated bacteriophage in electrospun fibres, 16th International Symposium on Wood, Fiber and Pulping Chemistry, Tianjin, People’s Republic of China. Proceedings ISWFPC 2: 1302-1306. |
[163] |
Lee SW, Belcher AM (2004) Virus-based fabrication of micro- and nanofibers using electrospinning. J Nano Lett 4: 387-390. doi: 10.1021/nl034911t
![]() |
[164] |
Salalha W, Kuhn J, Dror Y, et al. (2006) Encapsulation of bacteria and viruses in electrospun nanofibers. Nanotechnology 17: 4675-4681. doi: 10.1088/0957-4484/17/18/025
![]() |
[165] |
Korehei R, Kadla J (2013) Incorporation of T4 bacteriophage in electrospun fibres. J Appl Microbiol 114: 1425-1434. doi: 10.1111/jam.12158
![]() |
1. | M.Y. Thanuja, Sudhir H. Ranganath, Joseph A. Bonanno, Sangly P. Srinivas, Nanoliposomes for Sensing Local Osmolarity of the Tear Film on the Corneal Surface, 2022, 38, 1080-7683, 549, 10.1089/jop.2022.0044 | |
2. | Nannan Wang, Yunsen Zhang, Wei Wang, Zhuyifan Ye, Hongyu Chen, Guanghui Hu, Defang Ouyang, How can machine learning and multiscale modeling benefit ocular drug development?, 2023, 0169409X, 114772, 10.1016/j.addr.2023.114772 | |
3. | Jiaxin Han, Zihong Xie, Xingwei Zhu, Wanting Ruan, Meng Lin, Zhiqiang Xu, Lijie Miao, Junyang Zhong, Fan Lu, Liang Hu, The effects of eyelash extensions on the ocular surface, 2024, 47, 13670484, 102109, 10.1016/j.clae.2023.102109 |
Dimensional parameters | |||
Parameter | Description | Value | Reference |
$ \mu $ | Viscosity | $ 1.3 \times 10^{-3} $ Pa $ \cdot $ s | [58] |
$ \sigma_0 $ | Surface tension | 0.045 N $ \cdot $ m$ ^{-1} $ | [59] |
$ \rho $ | Density | $ 10^3 $ kg $ \cdot $ m$ ^{-3} $ | Water |
$ A^* $ | Hamaker constant | $ 6 \pi \times 3.5 \times 10^{-19} $ s $ \cdot $ m$ ^{-1} $ | [60] |
$ d $ | Initial TF thickness | $ 3-6 \times 10^{-6} $ m | Calculated; see [61] |
$ \ell $ | $ (\sigma_0/\mu/v_{\max})^{1/4}d $ | 0.426-0.694 mm | Calculated |
$ v_{\max} $ | Peak thinning rate | 0.5-25 $ \mu $m/min | [56] |
$ V_w $ | Molar volume of water | 1.8 $ \times 10^{-5} $ m$ ^3 \cdot $ mol$ ^{-1} $ | Water |
$ h_{\text{eq}} $ | Equilibrium thickness | 0.25 $ \mu $m | Estimated, [13] |
$ D_f $ | Diffusivity of fluorescein | 0.39 $ \times 10^{-9} $ m$ ^2 $/s | [62] |
$ D_o $ | Diffusivity of salt | 1.6 $ \times 10^{-9} $ m$ ^2 $/s | [63] |
$ c_0 $ | Isotonic osmolarity | 300 Osm/m$ ^3 $ | [28] |
$ P_o $ | Permeability of cornea | $ 12.1 \times 10^{-6} $ m/s | [18] |
$ \epsilon_f $ | Napierian extinction coefficient | $ 1.75 \times 10^7 $ m$ {-1} $ M$ ^{-1} $ | [64] |
Nondimensional parameters with typical values | ||
Parameter | Expression | Value |
$ \epsilon $ | $ \frac{d}{\ell} $ | $ 1\times10^{-2} $ |
$ S $ | $ \frac{\sigma_0 \epsilon^4}{\mu v_{max}} $ | 1 |
$ A $ | $ \frac{A^*}{\mu v_{max}~~ \ell^2} $ | $ 0.9 \times 10^{-3} $ |
$ \alpha $ | $ \frac{\alpha_0 \mu}{\rho \ell \epsilon^3} $ | 0.0368 |
$ P_c $ | $ \frac{P_o V_m~~ c_0}{v_{max}} $ | 0.1960 |
$ Pe_{f} $ | $ \frac{v_{max}~~ \ell}{\epsilon D_f} $ | 30.5 |
$ Pe_{c} $ | $ \frac{v_{max}~~ \ell}{\epsilon D_o} $ | 7.43 |
$ \phi $ | $ \epsilon_f f_{cr} d $ | $ 0.326 $ |
Non-dimensional parameters with typical values | ||
Parameter | Description | Expression |
$ \epsilon $ | Aspect ratio | $ d / \ell $ |
$ S $ | Contribution of surface tension | $ (\sigma_0 \epsilon^3)/(\mu U) $ |
$ M $ | Contribution of Marangoni effect | $ \epsilon (\Delta \sigma)_0 \left[t_s^3/(\sigma_0 \mu^3 d^3) \right]^{1/4} $ |
$ A $ | Non-dimensional Hamaker constant | $ A^*/\left[\epsilon (\Delta \sigma)_0 d \ell \right] $ |
$ P_c $ | Permeability of cornea | $ (P_o V_w c_0)/(\epsilon U) $ |
Pe$ _f $ | Péclet number for FL diffusion | $ U \ell/D_f $ |
Pe$ _c $ | Péclet number for salt iron diffusion | $ U \ell / D_o $ |
Pe$ _s $ | Péclet number for surface diffusion | $ \epsilon (\Delta \sigma)_0 \ell/(\mu D_s) $ |
$ \phi $ | Non-dimensional Napierian extinction coefficient | $ \epsilon_f f_{\text{cr}} d$ |
Dimensional parameters | |||
Parameter | Description | Value | Reference |
$ \mu $ | Viscosity | $ 1.3 \times 10^{-3} $ Pa $ \cdot $ s | [58] |
$ \sigma_0 $ | Surface tension | 0.045 N $ \cdot $ m$ ^{-1} $ | [59] |
$ \rho $ | Density | $ 10^3 $ kg $ \cdot $ m$ ^{-3} $ | Water |
$ A^* $ | Hamaker constant | $ 6 \pi \times 3.5 \times 10^{-19} $ s $ \cdot $ m$ ^{-1} $ | [60] |
$ d $ | Initial TF thickness | $ 3-6 \times 10^{-6} $ m | Calculated; see [61] |
$ \ell $ | $ (\sigma_0/\mu/v_{\max})^{1/4}d $ | 0.426-0.694 mm | Calculated |
$ v_{\max} $ | Peak thinning rate | 0.5-25 $ \mu $m/min | [56] |
$ V_w $ | Molar volume of water | 1.8 $ \times 10^{-5} $ m$ ^3 \cdot $ mol$ ^{-1} $ | Water |
$ h_{\text{eq}} $ | Equilibrium thickness | 0.25 $ \mu $m | Estimated, [13] |
$ D_f $ | Diffusivity of fluorescein | 0.39 $ \times 10^{-9} $ m$ ^2 $/s | [62] |
$ D_o $ | Diffusivity of salt | 1.6 $ \times 10^{-9} $ m$ ^2 $/s | [63] |
$ c_0 $ | Isotonic osmolarity | 300 Osm/m$ ^3 $ | [28] |
$ P_o $ | Permeability of cornea | $ 12.1 \times 10^{-6} $ m/s | [18] |
$ \epsilon_f $ | Napierian extinction coefficient | $ 1.75 \times 10^7 $ m$ {-1} $ M$ ^{-1} $ | [64] |
Nondimensional parameters with typical values | ||
Parameter | Expression | Value |
$ \epsilon $ | $ \frac{d}{\ell} $ | $ 1\times10^{-2} $ |
$ S $ | $ \frac{\sigma_0 \epsilon^4}{\mu v_{max}} $ | 1 |
$ A $ | $ \frac{A^*}{\mu v_{max}~~ \ell^2} $ | $ 0.9 \times 10^{-3} $ |
$ \alpha $ | $ \frac{\alpha_0 \mu}{\rho \ell \epsilon^3} $ | 0.0368 |
$ P_c $ | $ \frac{P_o V_m~~ c_0}{v_{max}} $ | 0.1960 |
$ Pe_{f} $ | $ \frac{v_{max}~~ \ell}{\epsilon D_f} $ | 30.5 |
$ Pe_{c} $ | $ \frac{v_{max}~~ \ell}{\epsilon D_o} $ | 7.43 |
$ \phi $ | $ \epsilon_f f_{cr} d $ | $ 0.326 $ |
Non-dimensional parameters with typical values | ||
Parameter | Description | Expression |
$ \epsilon $ | Aspect ratio | $ d / \ell $ |
$ S $ | Contribution of surface tension | $ (\sigma_0 \epsilon^3)/(\mu U) $ |
$ M $ | Contribution of Marangoni effect | $ \epsilon (\Delta \sigma)_0 \left[t_s^3/(\sigma_0 \mu^3 d^3) \right]^{1/4} $ |
$ A $ | Non-dimensional Hamaker constant | $ A^*/\left[\epsilon (\Delta \sigma)_0 d \ell \right] $ |
$ P_c $ | Permeability of cornea | $ (P_o V_w c_0)/(\epsilon U) $ |
Pe$ _f $ | Péclet number for FL diffusion | $ U \ell/D_f $ |
Pe$ _c $ | Péclet number for salt iron diffusion | $ U \ell / D_o $ |
Pe$ _s $ | Péclet number for surface diffusion | $ \epsilon (\Delta \sigma)_0 \ell/(\mu D_s) $ |
$ \phi $ | Non-dimensional Napierian extinction coefficient | $ \epsilon_f f_{\text{cr}} d$ |