3D simulation of the aqueous flow in the human eye

Drug delivery from a ring implant attached to intraocular lens: An in-silico investigation

Multiple iterations required to design ocular implants, which will last for the desired operational period of months or even years, necessitate the use of in-silico models for ocular drug delivery. In this study, we developed an in-silico model to simulate the flow of Aqueous Humor (AH) and drug delivery from an implant to the Trabecular Meshwork (TM). The implant, attached to the side of the intraocular lens (IOL), and the TM are treated as porous media, with their effects on AH flow accounted for using the Darcy equation. This model accurately predicts the physiological values of Intraocular Pressure (IOP) for both healthy individuals and glaucoma patients, as reported in the literature. Results reveal that the effective diffusivity of the drug within the implant is the critical parameter that can alter the bioavailability time period (BTP) from a few days to months. Intuitively, BTP should increase as effective diffusivity decreases. However, we discovered that with lower levels of initial drug loading, BTP declines when effective diffusivity falls below a specific threshold. Our findings further reveal that, while AH flow has a minimal effect on the drug release profile at the implant site, it significantly impacts drug availability at the TM.

Low Intraocular Pressure Induces Fibrotic Changes in the Trabecular Meshwork and Schlemm's Canal of Sprague Dawley Rats

Purpose

Continuous artificial aqueous humor drainage in the eyes of patients with glaucoma undergoing trabeculectomy likely exerts abnormal shear stress. However, it remains unknown how changes in intraocular pressure (IOP) can affect aqueous humor outflow (AHO).

Methods

Here, we induced and maintained low intraocular pressure (L-IOP) in healthy Sprague Dawley (SD) rats by puncturing their eyes using a tube (200-µm diameter) for 2 weeks. After the rats were euthanized, their eyes were removed, fixed, embedded, stained, and scanned to analyze the physiological and pathological changes in the trabecular meshwork (TM) and Schlemm's canal (SC). We measured SC parameters using ImageJ software and assessed the expression of various markers related to flow shear stress (KLF4), fibrosis (TGF-β1, TGF-β2, α-SMA, pSmad1/5, pSmad2/3, and fibronectin), cytoskeleton (integrin β1 and F-actin), diastolic function (nitric oxide synthase and endothelial nitric oxide synthase [eNOS]), apoptosis (cleaved caspase-3), and proliferation (Ki-67) using immunofluorescence or immunohistochemistry.

Results

L-IOP eyes showed a larger SC area, higher eNOS expression, and lower KLF4 and F-actin expression in the TM and SC (both P < 0.05) than control eyes. The aqueous humor of L-IOP eyes had a higher abundance of fibrotic proteins and apoptotic cells than that of control eyes, with significantly higher TGF-β1, α-SMA, fibronectin, and cleaved caspase-3 expression (all P < 0.05).

Conclusions

In conclusion, a persistence of L-IOP for 2 weeks may contribute to fibrosis in the TM and SC and might be detrimental to conventional AHO in SD rat eyes.

Translational Relevance

Clinicians should consider that aberrant shear force induced by aqueous humor fluctuation may damage AHO outflow channel when treating patients.

A computational fluid dynamics investigation of endothelial cell damage from glaucoma drainage devices

Glaucoma drainage devices (GDDs) are prosthetic-treatment devices for treating primary open-angle glaucoma. Despite their effectiveness in reducing intraocular pressures (IOP), endothelial cell damage (ECD) is a commonly known side-effect. There have been different hypotheses regarding the reasons for ECD with one being an induced increase in shear on the corneal wall. A computational fluid dynamics (CFD) model was used to investigate this hypothesis in silico. The Ahmed Glaucoma Valve (AGV) was selected as the subject of this study using an idealised 3D model of the anterior chamber with insertion angles and positions that are commonly used in clinical practice. It was found that a tube-cornea distance of 1.27 mm or greater does not result in a wall shear stress (WSS) above the limit where ECD could occur. Similarly, a tube-cornea angle of 45° or more was shown to be preferable. It was also found that the ECD region has an irregular shape, and the aqueous humour flow fluctuates at certain insertion angles and positions. This study shows that pathological amounts of WSS may occur as a result of certain GDD placements. Hence, it is imperative to consider the associated fluid force interactions when performing the GDD insertion procedure.

A Mechanistic Model of Aqueous Humor Flow to Study Effects of Angle Closure on Intraocular Pressure

Purpose

To study the relationship between the circumferential extent of angle closure and elevation in intraocular pressure (IOP) using a novel mechanistic model of aqueous humor (AH) flow.

Methods

AH flow through conventional and unconventional outflow pathways was modeled using the unified Stokes and Darcy equations, which were solved using the finite element method. The severity and circumferential extent of angle closure were modeled by lowering the permeability of the outflow pathways. The IOP predicted by the model was compared with biometric and IOP data from the Chinese American Eye Study, wherein the circumferential extent of angle closure was determined using anterior segment OCT measurements of angle opening distance.

Results

The mechanistic model predicted an initial linear rise in IOP with increasing extent of angle closure which became nonlinear when the extent of closure exceeded around one-half of the circumference. The nonlinear rise in IOP was associated with a nonlinear increase in AH outflow velocity in the open regions of the angle. These predictions were consistent with the nonlinear relationship between angle closure and IOP observed in the clinical data.

Conclusions

IOP increases rapidly when the circumferential extent of angle closure exceeds 180°. Residual AH outflow may explain why not all angle closure eyes develop elevated IOP when angle closure is extensive.

Translational Relevance

This study provides insight into the extent of angle closure that is clinically relevant and confers increased risk of elevated IOP. The proposed model can be utilized to study other mechanisms of impaired aqueous outflow.

Introduction and Clinical Validation of an Updated Biomechanically Corrected Intraocular Pressure bIOP (v2)

PURPOSE

To improve the stability of the Corvis ST biomechanically-corrected intraocular pressure measurements (bIOP) after refractive surgery and its independence of corneal biomechanics.

METHODS

A parametric study was carried out using numerical models simulating the behavior of the eye globe under the effects of IOP and Corvis ST external air pressure and used to develop a new algorithm for bIOP; bIOP(v2). It was tested on 528 healthy participants to evaluate correlations with CCT and age. Its ability to compensate for the geometrical changes was tested in 60 LASIK and 80 SMILE patients with six months follow up. The uncorrected Corvis ST IOP (CVS-IOP) and the two versions of biomechanically corrected IOP; bIOP(v1) and bIOP(v2), were compared.

RESULTS

In the healthy dataset, bIOP(v2) had weak and non-significant correlation with both CCT (R = -0.048, p = .266) and age (R = 0.085, p = .052). For bIOP(v1), the correlation was non-significant with CCT (R = -0.064, p = .139) but significant with age (R = -0.124, p < .05). In both LASIK and SMILE groups, the median change in bIOP(v2) following surgery was below 1 mmHg at follow-up stages and the interquartile range was smaller than both bIOP(v1) and CVS-IOP.

CONCLUSION

The bIOP(v2) algorithm performs better than bIOP(v1) and CVS-IOP in terms of correlation with CCT and age. The bIOP(v2) also demonstrated the smallest variation after LASIK and SMILE refractive surgeries indicating improved ability to compensate for geometrical changes.

An aqueous humour fluid dynamic study for normal and glaucomatous eye conditions

Glaucoma is the leading cause of irreversible blindness worldwide. Currently, the only treatable risk factor for glaucoma is elevated intraocular pressure (IOP). Glaucoma is commonly caused due to a decreased permeability of the trabecular meshwork, a porous structure at the eye outlet. This prevents the effective outflow of aqueous humour, increasing IOP. This study aims to simulate both normal and glaucomatous conditions of aqueous humour flow in the eye via computational fluid dynamics (CFD). Using clinical data, an idealised geometrical model of the eye was created. Darcy's law was employed to calculate the permeability values for various IOPs, which was then applied to the CFD model. Subsequently, verifiable and validated models for a normal and glaucomatous eye were achieved. Clinical Relevance- Glaucoma surgical treatments are often met with post-operative complications due to an insufficient or even excessive outflow of aqueous humour. The resulting glaucomatous eye model from this study can be used to test how different glaucoma filtration surgeries affect the efficacy of aqueous humour outflow. In turn, the most effective glaucoma surgical procedure may be identified for specific eye geometries according to race, age, gender, etc.

Rational Design of Microfluidic Glaucoma Stent

Glaucoma is a common, irreparable eye disease associated with high intraocular pressure. One treatment option is implantation of a stent to lower the intraocular pressure. A systematic approach to develop a microchannel stent meshwork that drains aqueous humor from the anterior chamber of the eye into the subconjunctival space is presented. The stent has a large number of outlets within its mesh structure that open into the subconjunctiva. The development approach includes a flow resistance model of the stent. Local adaption of the stent’s tubular dimensions allows for adjustment of the flow resistance. In this way, an evenly distributed outflow into the subconjunctiva is achieved. We anticipate that microblebs will form at the stent outlets. Their size is crucial for drainage and control of intraocular pressure. An analytical model for bleb drainage is developed based on the porous properties of the subconjunctival tissue. Both models—the stent flow resistance model and the bleb drainage model—are verified by numerical simulation. The models and numerical simulation are used to predict intraocular pressure after surgery. They allow for a systematic and personalized design of microchannel stents. Stents designed in this way can stabilize the intraocular pressure between an upper and lower limit.

Numerical study of aqueous humor flow and iris deformation with pupillary block and the efficacy of laser peripheral iridotomy

BACKGROUND

Disclosing the mechanism of primary angle closure glaucoma with pupillary block is important to the diagnosis as well as treatments, such as the laser peripheral iridotomy. Comparing with abundant clinical researches, there have been fewer quantitative studies of the aqueous humor flows with synechia iris configurations, and the efficacy of laser peripheral iridotomy in treating glaucoma.

METHODS

Based on the mathematical models of aqueous humor flow and iris deformation, the flow fields were simulated by computational fluid dynamics with normal and synechia iris configurations (iris-lens gap of 30, 5 and 2 μm, respectively), and through one-way fluid-structure coupling technique the deformations of the iris under the flow field pressure were calculated by finite element analysis. The efficacy of glaucoma treatment with different orifice sizes was also investigated.

FINDINGS

Results show that the pressure difference between anterior and posterior chambers and iris deformation increase dramatically with the iris-lens gap distance less than 5 μm, and when further decreasing this gap may lead the iris touch the cornea causing angle closure glaucoma with noticeable iris bombé. Laser peripheral iridotomy simulation results show that iridotomy size of 0.2 mm can effectively decrease the pressure difference across the iris and relieve iris bombé.

INTERPRETATION

This is a biomechanical numerical study, and the results are reasonable compare to those of published works. It may shed additional light on the diagnosis and treatment of angle closure glaucoma with pupillary block.

Ocular Fluid Mechanics and Drug Delivery: A Review of Mathematical and Computational Models

The human eye is a complex biomechanical structure with a range of biomechanical processes involved in various physiological as well as pathological conditions. Fluid flow inside different domains of the eye is one of the most significant biomechanical processes that tend to perform a wide variety of functions and when combined with other biophysical processes play a crucial role in ocular drug delivery. However, it is quite difficult to comprehend the effect of these processes on drug transport and associated treatment experimentally because of ethical constraints and economic feasibility. Computational modeling on the other hand is an excellent means to understand the associated complexity between these aforementioned processes and drug delivery. A wide range of computational models specific to different types of fluids present in different domains of the eye as well as varying drug delivery modes has been established to understand the fluid flow behavior and drug transport phenomenon in an insilico manner. These computational models have been used as a non-invasive tool to aid ophthalmologists in identifying the challenges associated with a particular drug delivery mode while treating particular eye diseases and to advance the understanding of the biomechanical behavior of the eye. In this regard, the author attempts to summarize the existing computational and mathematical approaches proposed in the last two decades for understanding the fluid mechanics and drug transport associated with different domains of the eye, together with their application to modify the existing treatment processes.

Dynamic volume-rendered optical coherence tomography pupillometry

PURPOSE

To assess intrapupillary space (IPS) changes in healthy subjects with regard to decreased iris motility in patients with pseudoexfoliation glaucoma (PEXG) or non-arteritic anterior ischaemic optic neuropathy (NAION) in a feasibility study in a clinical environment.

METHODS

Scotopic and photopic IPS measurements using three-dimensionally rendered swept-source optical coherence tomography (SS-OCT) data were obtained and compared for all subjects. Intrapupillary space (IPS) parameters were evaluated such as absolute volumetric differences, relative light response for volumetric ratios and pupillary ejection fraction (PEF) for functional contraction measurements.

RESULTS

From a total of 122 IPS from 66 subjects, 106 IPS were eligible for comparison providing values for 72 normal, 30 PEXG and 4 NAION eyes. In healthy, PEXG and NAION subjects, scotopic overall mean IPS was 8.90, 3.45 and 4.16 mm³ , and photopic overall mean IPS was 0.87, 0.74 and 1.13 mm³ , respectively. Three-dimensional contractility showed a mean absolute difference of 8.03 mm³ for normals (defined as 100% contractility), 2.72 mm³ for PEXG (33.88% of normal) and 3.03 mm³ for NAION (38.50% of normal) with a relative light response ratio between scotopic and photopic volumes of 10.26 (100%), 4.69 (45.70%) and 3.67 (35.78%), respectively. Pupillary ejection fraction (PEF) showed a contractile pupillary emptying of 88.11% for normals, 76.92% for PEXG and 70.91% for NAION patients.

CONCLUSION

This 3D pupillometry OCT assessment allows for quantitative measurements of pupil function, contractility and response to light. More specifically, PEF is presented as a potential (neuro)-pupillary outcome measure that could be useful in the monitoring of ophthalmic disorders that affect pupillary function.

Phacoemulsification and IOL-Implantation without Using Viscoelastics: Combined Modeling of Thermo Fluid Dynamics, Clinical Outcomes, and Endothelial Cell Density

Phacoemulsification is a widely used surgical method in cataract surgery with a high energy ultrasound source. The viscoelastic is considered to be tissue protective. The aim of this study is to investigate during surgery the impact of using viscoelastic versus no viscoelastic on clinical outcomes, potential complications and effect on endothelial cell density. The study group included 64 patients, who were subjected to phacoemulsification using balanced salt solution (BSS). Control group consisted of 62 patients, who underwent phacoemulsification using Hyaloronic acid 1% Healon 1%. Student’s t-test was applied for statistical analysis. The simulations of temperature changes during phacoemulsification were performed by COMSOL Multiphysics software. In the BSS group, a mean endothelial cell loss (ECL) of 4.5% was measured one month postoperatively, while in the Healon group ECL was 5.3%. Data analysis showed no significant difference in ECL between the groups (Student’s t-test, p = 0.8). No significant difference was observed in endothelial cell morphology and IOP between the two groups pre- and postoperatively (all p > 0.05). The modeling of thermo fluid dynamics showed that the heating of the cornea is slightly less when Healon was used as irrigation fluid. The phacoemulsification technique can be performed by an experienced surgeon with viscoelastics or continuous anterior chamber (AC) irrigation on the same level of safety regarding endothelial cell damage, providing equally satisfying clinical outcomes.

Numerical simulation of the flow of aqueous humor in the Schlemm's canal

Schlemm's canal (SC) is a ring structure responsible for maintaining fluid homeostasis in the anterior chamber of the eye by draining aqueous humor (AH) from the trabecular meshwork (TM) into the collecting channel (CC). Obstruction of the AH flow in this conventional pathway leads to an elevation of intraocular pressure in glaucoma, which is internationally recognized as the most common cause of blindness. Although there are obvious morphological differences between patients and healthy individuals, the AH drainage capacity of SC with specific parameters is difficult to quantitatively evaluate due to the limitations of analyses or experimental technologies. In this paper, the AH flow dynamic features under various SC morphological parameters are studied by numerical simulations of 3D models. The axis length of the SC cross-section and the level of openness are chosen as the main morphological factors according to the clinical data. The flow features, including the velocity, streamline and wall shear stress, are analyzed. This study aimed not only to explore the feasibility of 3D numerical modeling in evaluating the SC drainage capacity, but also to reveal the relationship between SC morphological parameters and the dynamic features of AH drainage.

Computational and experimental analysis of a Glaucoma flat drainage device

This paper presents a computational and experimental analysis of a glaucoma flat drainage device (FDD). The FDD consists of a metallic microplate placed into the eye sclerocorneal limbus, which creates a virtual path between the anterior chamber and its exterior, allowing the intraocular pressure (IOP) to be kept in a normal range. It also uses the surrounding tissue as a flow regulator in order to provide close values of IOP for a wide range of aqueous humor (AH) flow rates. The Neo Hookean hyperelastic model is used for the solid part, while the Reynolds thin film fluid model is used for the fluid part. On the other hand, a gravitational-driven flow test is implemented in order to validate the simulation process. An in vitro experiment evaluated the flow characteristics of the device implanted in fourteen extirpated pig eyes, giving as a result the best-fit for the Young modulus of the tissue surrounding the device. Finally, according to the resulting computational model, for a range of 1.4-3.1 μL/min, the device presents a pressure variation range of 6-7.5 mmHg.

Effect of aging on heat transfer, fluid flow and drug transport in anterior human eye: A computational study

There are a lot of geometrical and morphological changes that happen in the human eye with age. Primary open-angle glaucoma, which is caused by the increase in intraocular pressure inside the anterior chamber of the eye is also associated with the physiological aging of the eye. Therefore, it is crucial to understand the effects of aging on drug delivery in the human eye when applied topically. Consequently, a numerical model of topical drug delivery for an aging human eye has been developed using commercial software COMSOL Multiphysics in the current study. Three different age groups (young, middle and old) have been considered and the changes in geometrical and tissue properties of different domains of the eye with age have been included in the numerical model. The effect of aging on heat transfer, aqueous humor flow, intraocular pressure and drug concentration in different domains and orientations of the eye have been investigated. Additionally, an attempt has been made to predict the best class of anti-glaucomatic treatment in silico that should be preferred to treat primary open-angle glaucoma effectively. Results illustrate that there is a decrease in the average corneal temperature and an increase in the temperature deviation across the cornea with age. Further, there is a decrease in the aqueous humor flow magnitude in the anterior chamber of the eye and an increase in intraocular pressure in the anterior chamber of older age groups, which leads to primary open-angle glaucoma. The reduced aqueous humor flow leads to increased drug concentration in the anterior chamber as well as iris and reduced drug concentration in the trabecular mesh of the older age groups, thereby affecting the treatment efficacy. Additionally, our simulated results demonstrate that anti-glaucomatic treatments should be more focused on treating the trabecular mesh rather than the ciliary body of the eye.

Numerical model to predict and compare the hypotensive efficacy and safety of minimally invasive glaucoma surgery devices

Purpose

To predict and compare the hypotensive efficacy of three minimally-invasive glaucoma surgery (MIGS) implants through a numerical model.

Methods

Post-implant hypotensive efficacy was evaluated by using a numerical model and a computational fluid dynamics simulation. Three different devices were compared: the XEN 45 stent (tube diameter, 45 μm), the XEN 63 stent (63 μm) and the PreserFlo microshunt (70 μm). The influence of the filtration bleb pressure (Bp) and tube diameter, length, and position within the anterior chamber (AC) on intraocular pressure (IOP) were evaluated.

Results

Using baseline IOPs of 25, 30 and 50 mmHg, respectively, the corresponding computed post-implant IOPs for each device were as follows: XEN 45: 17 mmHg (29% decrease), 19 mmHg (45%) and 20 mmHg (59%) respectively; XEN 63: 13 mmHg (48%), 13 mmHg (62%), and 13 mmHg (73%); PreserFlo: 12 mmHg (59%), 13 mmHg (73%) and 13 mmHg (73%). At a baseline IOP of 35 mmHg with an increase in the outflow resistance within the Bp from 5 to 17 mmHg, the hypotensive efficacy for each device was reduced as follows: XEN45: 54% to 37%; XEN 63: 74% to 46%; and PreserFlo: 75% to 47%. The length and the position of the tube in the AC had only a minimal (non-significant) effect on IOP (<0.1 mmHg).

Conclusions

This hydrodynamic/numerical model showed that implant diameter and bleb pressure are the two most pertinent determinants of hypotensive efficacy. In distinction, tube length and position in the AC do not significantly influence IOP.

Numerical modeling of therapeutic lens drug delivery in the anterior human eye for the treatment of primary open-angle glaucoma

A numerical model of drug delivery from a therapeutic lens in the anterior portion of the human eye has been developed for a more effective treatment plan of primary open-angle glaucoma. The numerical model takes into account the drug diffusion through the therapeutic lens along with heat transfer and aqueous humor flow in different orientations of the human eye (supine (two-dimensional) as well as standing (three-dimensional)). Results illustrate that the drug diffuses through the therapeutic lens to the cornea and is convected into the anterior chamber of the eye due to the temperature gradient across the eye. In addition, eye orientation significantly affects drug delivery with supine orientation providing better and uniform drug exposure in different target regions of the eye as compared to standing in the case of the therapeutic lens. Furthermore, a comparison of the therapeutic efficacy of the therapeutic lens has been done with topical administration and the drug uptake results from both the drug delivery modes have been validated with the experimental data reported in the literature. The developed model may help ophthalmologists to comprehend the transport and retention of different drugs in different domains and orientations of the human eye when administered through a therapeutic lens.

Simulation of Air Puff Tonometry Test Using Arbitrary Lagrangian-Eulerian (ALE) Deforming Mesh for Corneal Material Characterisation

: Purpose: To improve numerical simulation of the non-contact tonometry test by using arbitrary Lagrangian-Eulerian deforming mesh in the coupling between computational fluid dynamics model of an air jet and finite element model of the human eye.

METHODS

Computational fluid dynamics model simulated impingement of the air puff and employed Spallart-Allmaras model to capture turbulence of the air jet. The time span of the jet was 30 ms and maximum Reynolds number was Re=2.3×104, with jet orifice diameter 2.4 mm and impinging distance 11 mm. The model of the human eye was analysed using finite element method with regional hyperelastic material variation and corneal patient-specific topography starting from stress-free configuration. The cornea was free to deform as a response to the air puff using an adaptive deforming mesh at every time step of the solution. Aqueous and vitreous humours were simulated as a fluid cavity filled with incompressible fluid with a density of 1000 kg/m3.

RESULTS

Using the adaptive deforming mesh in numerical simulation of the air puff test improved the traditional understanding of how pressure distribution on cornea changes with time of the test. There was a mean decrease in maximum pressure (at corneal apex) of 6.29 ± 2.2% and a development of negative pressure on a peripheral corneal region 2-4 mm away from cornea centre.

CONCLUSIONS

The study presented an improvement of numerical simulation of the air puff test, which will lead to more accurate intraocular pressure (IOP) and corneal material behaviour estimation. The parametric study showed that pressure of the air puff is different from one model to another, value-wise and distribution-wise, based on cornea biomechanical parameters.

Computational Study on the Biomechanics of Pupil Block Phenomenon

Pupil blocking force (PBF) can indicate the potential risk of pupil block (PB), which is considered as a main pathogenic factor of primary angle-closure glaucoma (PACG). However, the effect of PB on the PBF under different pupil diameters and iris-lens channel (ILC) distance was unknown. Besides, a simple and practical method to assess PBF has not been reported yet. In this study, 21 finite element models of eyes with various pupil diameters (2.4 mm–2.6 mm) and ILC (2 μm–20 μm) were constructed and were conducted to simulate aqueous humor flow by fluid-solid coupling numerical simulation. PBF in each model was calculated based on the numerical simulation results and was fitted using response surface methodology. The results demonstrated that ILC distance had a more significant effect than pupil diameter on PBF. With the decrease of ILC distance, the PBF increased exponentially. When the reduced distance was lower than 5 μm, the PBF exploded quickly, resulting in a high risk of iris bomb. The PBF also varied with pupil diameter, especially under the condition of narrow ILC. Both ILC distance and pupil diameter could explain more than 97% variation in PBF, and a second-order empirical model has been developed to be a good predictor of PBF. Based on the linear relationship between anterior chamber deformation and PBF, a threshold value of PBF was given to guide clinical decisions. This study could be used to investigate PACG pathological correlation and its pathogenesis, so as to provide a reference value for clinical diagnosis of PACG.

Determination of Corneal Biomechanical Behavior in-vivo for Healthy Eyes Using CorVis ST Tonometry: Stress-Strain Index

Purpose: This study aims to introduce and clinically validate a new algorithm that can determine the biomechanical properties of the human cornea in vivo.

Methods: A parametric study was conducted involving representative finite element models of human ocular globes with wide ranges of geometries and material biomechanical behavior. The models were subjected to different levels of intraocular pressure (IOP) and the action of external air puff produced by a non-contact tonometer. Predictions of dynamic corneal response under air pressure were analyzed to develop an algorithm that can predict the cornea's material behavior. The algorithm was assessed using clinical data obtained from 480 healthy participants where its predictions of material behavior were tested against variations in central corneal thickness (CCT), IOP and age, and compared against those obtained in earlier studies on ex-vivo human ocular tissue.

Results: The algorithm produced a material stiffness parameter (Stress-Strain Index or SSI) that showed no significant correlation with both CCT (p > 0.05) and IOP (p > 0.05), but was significantly correlated with age (p < 0.01). The stiffness estimates and their variation with age were also significantly correlated (p < 0.01) with stiffness estimates obtained earlier in studies on ex-vivo human tissue.

Conclusions: The study introduced and validated a new method for estimating the in vivo biomechanical behavior of healthy corneal tissue. The method can aid optimization of procedures that interfere mechanically with the cornea such as refractive surgeries and introduction of corneal implants.

A novel patient-oriented numerical procedure for glaucoma drainage devices

The present work analyses the performance of four glaucoma drainage devices, by means of a novel patient-oriented numerical procedure. The procedure is based on the three-dimensional geometry reconstruction from the stacks of tomographic images of a human eye, at different angles, on meshing and on thermo-fluid dynamics modelling activities, carried out on the reconstructed computational domain. The current three-dimensional eye model considers anterior chamber (AC), trabecular meshwork, Schlemm's canal, and collector channels, making use of generalised porous medium approach for modelling ocular porous tissue and cavities. The intraocular pressure (IOP) management inside AC of human eye is analysed, by comparing the results obtained for four drainage devices implanted in a human eye for glaucoma treatment, ie, ExPRESS shunt, iStent inject, SOLX gold micro shunt, and the novel silicon shunt device. The numerical results allow predicting the effects of the installation of these implants on human eyes, in terms of IOP decrease, aqueous humour velocity, pressure, friction coefficient, and local Nusselt number, pointing out the clear distinction between pre-operative and post-operative eye conditions for different glaucoma surgical techniques.

Computational simulation of aqueous humour dynamics in the presence of a posterior-chamber versus iris-fixed phakic intraocular lens

Purpose

To compare aqueous humour (AH) dynamics in the presence of a precrystalline (Implantable Collamer Lens^®; ICL) or iris-fixed (Artiflex^®) phakic intraocular lens (PIOL).

Methods

By computational fluid dynamics simulation, AH flow was modelled through a peripheral iridotomy (PI) or central lens hole (both 360 μm) in the presence of an Artiflex or ICL lens, respectively. The impacts of AH flow were then determined in terms of wall shear stress (WSS) produced on the endothelium or crystalline lens. Effects were also modelled for different scenarios of pupil diameter (PD 3.5 or 5.5 mm), ICL vault (100, 350, 800 μm) and number of Artiflex iridotomies (1 or 2) and location (12 or 6 o’clock).

Results

For a PD of 3.5 mm, AH volumes flowing from the posterior to the anterior chamber were 37.6% of total flow through the lens hole (ICL) and 84.2% through PI (Artiflex). For an enlarged PD (5.5 mm), corresponding values were 10.3% and 81.9% respectively, so PI constitutes a very efficient way of evacuating AH. Central endothelial WSS in Pa was lower for the large vault ICL and the Artiflex (1⁻⁰³ and 1.1⁻⁰³ respectively) compared to the PIOL-free eye (1.6⁻⁰³). Crystalline lens WSS was highest for the lowest vault ICL (1⁻⁰⁴).

Conclusions

AH flow varied according to the presence of a precrystalline or iris-fixed intraocular lens. Endothelial WSS was lower for an implanted ICL with large vault and Artiflex than in the PIOL-free eye, while highest crystalline WSS was recorded for the lowest vault ICL.

Automated segmentation of ophthalmological images by an optical based approach for early detection of eye tumor growing

PURPOSE

Iris neoplasm is a non-symptom cancer that causes a gradual loss of sight. The first purpose of this study was to present a novel and automatic method for segmenting the iris tumors and detecting the corresponding areas changing along time. The second aim of this work was to investigate several recently published methods after being applied for the iris tumors segmentation.

METHODS

Our approach consists firstly in segmenting the iris region by using the Vander Lugt correlator based active contour method. Secondly, by treating only the iris region, a K-means clustering model was used to assign the tumorous tissue to one pixel-cluster. This model is quite sensitive to the center initialization and to the choice of the distance measure. To solve these problems, a proportional probability based approach was introduced for the cluster center initialization, and the impact of several distance measure was investigated. The proposed method and the different comparative methods were evaluated on two databases: the Eye Cancer and the Miles Research.

RESULTS

Results reported using several performance metrics reveal that the first step assures the detection of all iris tumors with an accuracy of 100%. Additionally, the proposed method yields better performance compared to the recently published methods.

A generalised porous medium approach to study thermo-fluid dynamics in human eyes

The present work describes the application of the generalised porous medium model to study heat and fluid flow in healthy and glaucomatous eyes of different subject specimens, considering the presence of ocular cavities and porous tissues. The 2D computational model, implemented into the open-source software OpenFOAM, has been verified against benchmark data for mixed convection in domains partially filled with a porous medium. The verified model has been employed to simulate the thermo-fluid dynamic phenomena occurring in the anterior section of four patient-specific human eyes, considering the presence of anterior chamber (AC), trabecular meshwork (TM), Schlemm's canal (SC), and collector channels (CC). The computational domains of the eye are extracted from tomographic images. The dependence of TM porosity and permeability on intraocular pressure (IOP) has been analysed in detail, and the differences between healthy and glaucomatous eye conditions have been highlighted, proving that the different physiological conditions of patients have a significant influence on the thermo-fluid dynamic phenomena. The influence of different eye positions (supine and standing) on thermo-fluid dynamic variables has been also investigated: results are presented in terms of velocity, pressure, temperature, friction coefficient and local Nusselt number. The results clearly indicate that porosity and permeability of TM are two important parameters that affect eye pressure distribution. Graphical abstract Velocity contours and vectors for healthy eyes (top) and glaucomatous eyes (bottom) for standing position.

Development of a device to prevent the late consequence of non-treated/lately diagnosed glaucoma

People residing in rural and remote areas (worldwide) have substantially worse outcome in early detection and diagnosis of glaucoma than those living in metropolitan areas. This gap can be reduced by improved glaucoma diagnosis activities in primary care, but there is little empirical evidence regarding use of tonometry in rural settings, or the expertise associated with quality of eye care. This article describes a feasibility study of a novel through-eyelid tonometer based on the use of an instrumented form of indentation and applanation tonometry.

Microstructure visualization of conventional outflow pathway and finite element modeling analysis of trabecular meshwork

BACKGROUND

The intraocular pressure (IOP) is maintained through a dynamic equilibrium between the production and drainage of aqueous humor. Elevation of intraocular pressure is mainly caused by the blocking of aqueous humor outflow pathway. Therefore, it is particularly important to study the structure of drainage pathway and the effect of ocular hypertension at the process of aqueous humor outflow.

METHODS

Conventional drainage pathway of aqueous humor, including trabecular meshwork (TM), Schlemm's canal (SC) and aqueous vein, were imaged by using trans-scleral imaging method with lateral resolution of 2 μm. For quantitative assessment, the morphological parameters of the TM were measured with different IOP levels via a combination of measurements and simulations.

RESULTS

Images of the TM and the adjacent tissues were obtained. The porosity of TM with normal intraocular pressure varies from 0.63 to 0.74 as the depth increases, while in high IOP it is changed from 0.44 to 0.59. The diameter of aqueous vein varies from 32 to 43 μm, and is smaller than that of SC, which varies from 48 to 64.67 μm.

CONCLUSIONS

Our research provides a non-contact method to visualize the microstructure of tissue for clinical examination associated with the blocking of the outflow pathway of aqueous humor in humans. The three-dimensional (3D) microstructures of limbus and the results of finite element modeling analysis of the TM model will serve for the future evaluation of new glaucoma surgical techniques.

Development and validation of a correction equation for Corvis tonometry

PRIMARY OBJECTIVE

This study uses numerical analysis and validation against clinical data to develop a method to correct intraocular pressure (IOP) measurements obtained using the Corvis Tonometer for the effects of central corneal thickness (CCT), and age.

MATERIALS AND METHODS

Finite element analysis was conducted to simulate the effect of tonometric air pressure on the intact eye globe. The analyses considered eyes with wide variations in IOP (10-30 mm Hg), CCT (445-645 microns), R (7.2-8.4 mm), shape factor, P (0.6-1) and age (30-90 years). In each case, corneal deformation was predicted and used to estimate the IOP measurement by Corvis (CVS-IOP). Analysis of the results led to an algorithm relating estimates of true IOP as a function of CVS-IOP, CCT and age. All other parameters had negligible effect on CVS-IOP and have therefore been omitted from the algorithm. Predictions of corrected CVS-IOP, as obtained by applying the algorithm to a clinical data-set involving 634 eyes, were assessed for their association with the cornea stiffness parameters; CCT and age.

RESULTS

Analysis of CVS-IOP measurements within the 634-large clinical data-set showed strong correlation with CCT (3.06 mm Hg/100 microns, r(2) = 0.204) and weaker correlation with age (0.24 mm Hg/decade, r(2) = 0.009). Applying the algorithm to IOP measurements resulted in IOP estimations that became less correlated with both CCT (0.04 mm Hg/100 microns, r(2) = 0.005) and age (0.09 mm Hg/decade, r(2) = 0.002).

CONCLUSIONS

The IOP correction process developed in this study was successful in reducing reliance of IOP measurements on both corneal thickness and age in a healthy European population.

Lymphatic regulator PROX1 determines Schlemm's canal integrity and identity

Schlemm's canal (SC) is a specialized vascular structure in the eye that functions to drain aqueous humor from the intraocular chamber into systemic circulation. Dysfunction of SC has been proposed to underlie increased aqueous humor outflow (AHO) resistance, which leads to elevated ocular pressure, a factor for glaucoma development in humans. Here, using lymphatic and blood vasculature reporter mice, we determined that SC, which originates from blood vessels during the postnatal period, acquires lymphatic identity through upregulation of prospero homeobox protein 1 (PROX1), the master regulator of lymphatic development. SC expressed lymphatic valve markers FOXC2 and integrin α9 and exhibited continuous vascular endothelial-cadherin (VE-cadherin) junctions and basement membrane, similar to collecting lymphatics. SC notably lacked luminal valves and expression of the lymphatic endothelial cell markers podoplanin and lymphatic vessel endothelial hyaluronan receptor 1 (LYVE-1). Using an ocular puncture model, we determined that reduced AHO altered the fate of SC both during development and under pathologic conditions; however, alteration of VEGF-C/VEGFR3 signaling did not modulate SC integrity and identity. Intriguingly, PROX1 expression levels linearly correlated with SC functionality. For example, PROX1 expression was reduced or undetectable under pathogenic conditions and in deteriorated SCs. Collectively, our data indicate that PROX1 is an accurate and reliable biosensor of SC integrity and identity.

Experimental research on intraocular aqueous flow by PIV method

BACKGROUND

Aqueous humor flows regularly from posterior chamber to anterior chamber, and this flow much involves intraocular pressure, the eye tissue nutrition and metabolism.

PURPOSE

To visualize and measure the intraocular flow regular pattern of aqueous humor.

METHOD

Intraocular flow in the vitro eyeball is driven to simulate the physiological aqueous humor flow, and the flow field is measured by Particle Image Velocimetry(PIV). Fluorescent particle solution of a certain concentration was infused into the root of Posterior Chamber(PC) of vitro rabbit eye to simulate the generation of aqueous and was drained out at a certain hydrostatic pressure from the angle of Anterior Chamber(AC) to represent the drainage of aqueous. PIV method was used to record and calculate the flow on the midsagittal plane of the eyeball.

RESULTS

Velocity vector distribution in AC has been obtained, and the distribution shows symmetry feature to some extent. Fluorescent particle solution first fills the PC as it is continuously infused, then surges into AC through the pupil, flows upwards toward the central cornea, reflecting and scattering, and eventually converges along the inner cornea surface towards the outflow points at the periphery of the eyeball. Velocity values around the pupillary margin are within the range of 0.008-0.012 m/s, which are close to theoretical values of 0.0133 m/s, under the driving rate of 100 μl/min.

CONCLUSIONS

Flow field of aqueous humor can be measured by PIV method, which makes it possible to study the aqueous humor dynamics by experimental method. Our study provides a basis for experimental research on aqueous humor flow; further, it possibly helps to diagnose and treat eye diseases as shear force damage of ocular tissues and destructions on corneal endothelial cells from the point of intraocular flow field.

Stress free configuration of the human eye

Numerical simulations of eye globes often rely on topographies that have been measured in vivo using devices such as the Pentacam or OCT. The topographies, which represent the form of the already stressed eye under the existing intraocular pressure, introduce approximations in the analysis. The accuracy of the simulations could be improved if either the stress state of the eye under the effect of intraocular pressure is determined, or the stress-free form of the eye estimated prior to conducting the analysis. This study reviews earlier attempts to address this problem and assesses the performance of an iterative technique proposed by Pandolfi and Holzapfel [1], which is both simple to implement and promises high accuracy in estimating the eye's stress-free form. A parametric study has been conducted and demonstrated reliance of the error level on the level of flexibility of the eye model, especially in the cornea region. However, in all cases considered 3-4 analysis iterations were sufficient to produce a stress-free form with average errors in node location <10(-6)mm and a maximal error <10(-4)mm. This error level, which is similar to what has been achieved with other methods and orders of magnitude lower than the accuracy of current clinical topography systems, justifies the use of the technique as a pre-processing step in ocular numerical simulations.