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.
Glaucoma is the leading cause of incurable blindness worldwide [1–5]. The most common clinical type in East Asia is primary angle-closure glaucoma (PACG) [6–8], which is characterized by a narrow anterior chamber angle (ACA) . ACA is considered as the outlet of aqueous humor (AH) in the eye. Once the narrowing of ACA occurred, the outflow of AH would be blocked, resulting in an obvious increase in intraocular pressure (IOP), which could cause progressive damage to the structure and function of the optic nerve . Therefore, the maintenance of IOP is critical to the stability of eyeball structure and eye health [11, 12].
IOP is regulated by the flow of aqueous humor (AH) . AH is secreted by the ciliary body continuously, flows from the posterior chamber (PC) into the anterior chamber (AC) through the pupil, and drains through the trabecular meshwork (Figure 1). Under the normal physiological condition, the rate of AH production is equal to that of drainage, maintaining a stable IOP in the eye. However, many factors can block AH flowing and the most common one is the pupillary block (PB) in PACG patients . Under PB physiological condition, iris-lens channel (ILC) became narrow, which increased the resistance in the flow pathway from PC to AC and hence resulted in a higher IOP in PC and consequently increased the pressure difference between PC and AC . Finally, the elevation of the differential pressure (Pd) would push iris anteriorly and further narrow ACA, which eventually caused a high IOP and further an irreversible damage to the optic nerve.
To quantify the effect of PB on the resistance of AH flow, pupil blocking force (PBF) has been proposed to represent that AH encountered finite resistance in passing through the narrow ILC. By assuming that iris was a type of linear elastic material, a formula based on Hooke’s law to calculate PBF has been deduced in early research [16, 17]. However, this method only included the limited physiological influencing factors of PBF and neglected the main resistance in ILC, i.e., AH flow resistance. A method based on computational fluid dynamics (CFD) has been successfully used to study AH flow under PB condition and visualize the flow field . CFD is a powerful tool to analyze the biofluid flow as well as mechanical distribution, and it has been widely used in biomechanics. Therefore, PBF could be assessed precisely by using CFD technology.
As a biomechanical phenomenon, PBF could be affected by a variety of factors. Wang et al. compared the PB degree under different lens’ relative position and revealed that lens’ position could affect PBF significantly . Huang and Barocas found that lens curvature also had an obvious effect on PBF by setting an active sphincter to simulate PB conditions . Besides, other factors such as AH flow through the iris-lens passage [20, 21] and lens size  were proved to be the effect factors of PBF. However, only one factor was investigated in most of the previous research, and a systematic study of various factors has not been reported yet. Moreover, despite high importance, a credible calculation method of PBF was still unknown. Therefore, it was necessary to investigate the comprehensive effect of several factors including pupil diameter, eye size, AC depth, lens size, and its location on PBF as well as AH flow and further construct a computational method to evaluate PBF.
The aim of this study was to investigate the effect of different optical conditions on PBF and construct a model to evaluate PBF. Firstly, twenty-one eye models with different pupil diameters and ILC were constructed and were conducted with numerical simulations. Then, PBF was calculated based on a modified Bernoulli’s equation. After that, response surface methodology was utilized to fit the relationship between PBF and pupil diameter as well as ILC distance. The effect of pupil diameter, ILC distance, and their interaction on PBF was also discussed. Finally, the AC morphology including variation in ACA and iris deformation under different PBF was investigated and a threshold value of PBF was given to guide clinical decisions.
2.1. Pupil Blocking Force
AH mainly flows from PC to AC through ILC. Assuming that AH was ideal fluid, its mechanical energy should remain constant at any position based on Bernoulli’s equation. However, due to complex components, the viscidity of AH could not be neglected . Considering the fluid viscosity, Bernoulli’s equation could be modified aswhere is the pressure at one position, h is its height, ρ is the fluid density, is the fluid velocity, is the acceleration of gravity, hf represents the flow resistance during AH flow pathway, and C is constant. When AH flows through ILC (Figure 2), the following equation could be induced to describe the conservation of energy:where u1 and u2 are the average velocities in the two surfaces (S1 and S2), p1 and p2 are the pressures in these surfaces, A1 and A2 are their areas, respectively, Z is the vertical distance from S1 to S2, and Hf (unit of J/kg) is the resistance (i.e., PBF) when AH flows from S1 to S2 [14, 24].
Known that AH is secreted by the ciliary body, its volumetric flux (represented by Q) in the surfaces S1 and S2 is identical and can be calculated as the multiplication by area and velocity; therefore, the following equation can be obtained:
2.2. Geometrical Model
Due to the hard-to-measure pressure and area in real eye, numerical simulation is used to simulate the AH flow and assess those parameters in equation (4). The geometric models for finite element analysis, including cornea, iris, lens, and AH (Figures 3(a)–3(c)), were developed from the finite model in our previous study . All the models were established with similar geometry parameters as typical ocular dimensions [25–27], with the assumption that the anterior segment was axisymmetric. Notably, different ILC distances as well as pupil diameters were used to construct these models, while the other geometric parameters remained the same. Detailed dimensions of the models are shown in Table 1. Considering the fact that obvious resistance in AH could be observed when the ILC distance was under 20 μm and PB was defined when the distance was less than 5 μm [28, 29], seven levels (20 μm, 15 μm, 10 μm, 5 μm, 4 μm, 3 μm, and 2 μm) of ILC distances were chosen to simulate healthy and different PB conditions. Three levels (2.4 mm, 3 mm, and 3.6 mm) of pupil diameters  were set to simulate different pupil sizes in PACG patients. Based on full factorial design, twenty-one models were constructed and are shown in Figures 3(a)–3(c). All of these models were conducted to generate the tetrahedron mesh using ICEM CFD and Mechanical APDL (Figures 3(d) and 3(e)). The mesh was refined, and the independency was checked. The computational grids were exported to ANSYS Workbench (ANSYS Inc., Canonsburg, PA, USA) for fluid-solid coupling simulation analysis.
2.3. Fluid-Solid Coupling Analysis
In order to calculate the Pd between surfaces S1-2 and investigate the effect of PB on eye shape with different pupil diameters, all the finite element models were conducted to simulate AH flow by fluid-solid coupling numerical simulation. To model the motion of AH steady flow from PC to AC, the Navier–Stokes equation representing the conservation of the momentum is written as
From left to right, the terms in equation (5) represent the change in momentum caused by convection, pressure gradients, viscous diffusion, and gravity. Considering AH is an incompressible Newtonian viscous fluid, the continuity equation is given as
For the boundary conditions, the inlet (the yellow surface in Figure 3(b)) was considered as constant flux and the flux rate was set at 3 μL/min, which was equal to the generation rate of AH in real eyes . The outlet (the green surface in Figure 3(b)) was set as constant pressure (scleral venous pressure) for AH fluid. The exterior surface of lens, iris root (the blue surfaces in Figure 3(c)), and pupil (Figure 3(c) line 1) were set as fixed positions. The interface between the cornea and AH and the interface between iris and AH (the red surfaces in Figures 3(b) and 3(c)) were chosen as fluid-solid interaction surfaces to simulate the interaction between AH fluid and solid bodies . Moreover, the bottom (the blue surface in Figure 3(b)) was applied to the wall boundary condition. As for the material properties, the liquid properties in the simulation were assumed to be those of pure water. The solid properties of cornea and iris were set according to the measurement of real eyes. A nonlinear elasticity model (2nd Ogden model) was used to simulate the mechanical properties of iris. Cornea was considered to be a linear elastic material and iris was a nonlinear elastic material. Detailed information about material properties is shown in Table 2.
2.4. Soft Measurement Method of PBF
Based on the result of fluid-solid coupling analysis, the Hf in each model was calculated by equation (4). In order to investigate the effect of pupil diameters, ILC distances, and their interaction on PBF, the response surface methodology (RSM) was selected to fit the experimental data. RSM was considered as a simple and fast empirical tool to study both the effect of each individual factor and their interactions on the response variable [38, 39]. By using SPSS Statistics software, a RSM model was constructed and was optimized with a sequential quadratic programming approach. After that, analysis of variance was used to verify the significance of the optimized model and its parameters. Finally, an empirical equation was derived to predict PBF based on the measurement of ILC distance and pupil diameter.
2.5. Analysis and Calculation
The areas of surfaces S1-2 (A1-2) were calculated using SolidWorks2014 software, and the pressures (P1-2) in those surfaces were derived from the results of fluid-solid coupling analysis by CFD-Post in ANSYS Workbench. Then, the magnitude of PBF was calculated by equation (4) according to the above parameters. The deformation of iris (), equivalent stress, and equivalent elastic strain of iris were also obtained by CFD-Post. Because of the irregular curved surface of the iris, the anterior chamber angle can be measured by averaging the ACA values in the results of iris deformation images of repeated manual measurements.
3. Results and Discussion
3.1. Effects of Iris-Lens Channel Distance and Pupil Diameter on AH Flow
With various ILC distances and pupil diameters, the velocity and pressure distributions of AH within eyes are shown in Figure 4, where gravity was opposite to the direction of y-axis, representing the supine position in clinical. It was obvious that the main flux of AH from PC to AC was at a low velocity. However, the fluid was accelerated rapidly and the velocity achieved peak when AH flowed through ILC. Besides, an obvious difference in pressure between PC and AC could be observed because of the resistance of ILC, according to a previous study . The maximum velocity of AH (Vmax) and the Pd between PC and AC were different with various ILC distances and pupil diameters. Despite the same pupil diameter (3.6 mm), the Vmax and Pd were 5.39e − 4 m/s and 1.60 Pa when the ILC distance was 20 μm (Figure 4(a)), while they became 1.72e − 3 m/s and 257.00 Pa once the distance decreased to 2 μm (Figure 4(b)). When the pupil diameter decreased from 3.6 mm to 2.4 mm, the Vmax and Pd increased to 2.13e − 3 m/s and 392.00 Pa further, indicating both the ILC distance and diameter had effects on the Vmax as well as Pd. The previous study has demonstrated that Vmax and Pd could influence AC morphology , which was also confirmed by the AC deformation in our study (Figure 4). Iris bomb and ACA blocking could be observed obviously when the ILC became narrow (Figures 4(a) and 4(b)), according to clinical trials . Notably, the degree of iris bomb and ACA blocking varied with pupil diameter, indicating that pupil diameter was also an important factor to assess the risk of PB .
The results of Vmax and Pd under different ILC distances and pupil diameters are shown in Figure 5. With the decrease of ILC distance, an exponential increase of Vmax and Pd could be observed under each level of pupil diameter (Figure 5(a)). Besides, both Vmax and Pd varied with pupil diameter, especially under narrow ILC distances (Figures 5(a) and 5(b)). Similarly, an exponential increase of ACA could be observed with the decrease of ILC distance (Figure 5(c)), indicating a great effect of PB on AC morphology. The effect of pupil diameter on iris deformation showed greater differences under different levels of ILC distances (Figure 5(c)). When the pupil diameter decreased from 3.6 mm to 2.4 mm, ACA increased 85.48% under the ILC distance of 2 μm while only 23.88% under the ILC distance of 20 μm. Furthermore, the curves of the relationships between Vmax and ILC distance were not parallel under each level of pupil diameter. Similar phenomena could be observed in Figures 5(b) and 5(c), suggesting a potential interaction between ILC and pupil diameter on Vmax, Pd, and iris deformation .
3.2. Soft Measurement Method of PBF
Considering the obvious effects of ILC as well as pupil diameter on the Pd and AH flow velocity, it could be induced that these two factors also had an impact on PBF. The results of PBF in the twenty-one models are shown in Figure 6(a). A distinct difference in PBF values could be observed under various ILC distances. It was only 1049.2 J/kg − 2342.4 J/kg under a wide ILC (20 μm) while it achieved 256449.2 J/kg − 472235.6 J/kg under an extreme narrow ILC (2 μm). This phenomenon verified the rationality of ILC as an indicator to assess PB risks. Moreover, similar trends could be observed in the relationship between the ILC distance and PBF under each level of pupil diameter, which appeared to be exponential growth curves. The PBF exploded when the ILC distance decreased to less than 5 μm, verifying the reliability of 5 μm as the threshold of PB in clinical [28, 29]. Furthermore, the PBF values increased in different degrees under various ILC distances with the same change of pupil diameter from 3.6 mm to 2.4 mm, showing a potential interaction relationship between ILC distance and pupil diameter on PBF.
Despite high accuracy to assess PBF using finite element analysis, it was not suitable for clinical application because of complex operation and high cost. Therefore, it was necessary to establish a soft measurement method to evaluate PBF. Considering the exponential distributions of PBF under each pupil diameter level, the logarithm of PBF value was used to develop a regression model. Besides, ILC ratio and diameter ratio, defined by the ratio of the measured ILC distance and pupil diameter to normal values (ILC distance: 20 μm; pupil diameter: 3 mm), were calculated to eliminate the discrepancy in magnitude. The relationships between logarithmic PBF and the ILC ratio as well as the diameter ratio are shown in Figure 6(b). By applying multiple regression analysis on these data, a second-order polynomial equation for the logarithmic PBF was obtained as follows:where y represents the logarithmic PBF and x1 and x2 are the ILC ratio and pupil diameter ratio, respectively. The result of variance analysis showed that the determination coefficient of this model achieved 0.973, implying approximately 98% of the variability in logarithmic PBF could be expressed by the model. Therefore, the empirical relationship between PBF value and the measured parameters (including ILC distance and pupil diameter) was quantified using the following equation:where Y represents the value of PBF (J/kg) and X1 and X2 are, respectively, the measured ILC distance (μm) and pupil diameter (mm).
3.3. Effect of PBF on AC Morphology
In order to testify the feasibility of PBF as an indicator to assess the risks of PB in clinical, the relationship between PBF and variation of ACA was investigated and is shown in Figure 7. ACA increased with PBF, and a linear relationship was observed between ACA and PBF. The determination coefficient achieved 0.92, indicating that more than 92% variation in ACA could be explained by PBF. Similar phenomenon could be observed in the relationship between deformation of iris and PBF. With the increase of PBF, the deformation of iris increased linearly. The determination coefficient of the fitted line was 0.95, suggesting the linear model was able to express more than 95% of the variability in iris deformation. Iris deformation and variation in ACA were considered as the pathogenic mechanism of PACG . Increased iris deformation could decrease AC depth and narrow ACA, resulting in a severe obstruction of AH outflow, which could cause high IOP in eye and induce acute attack of angle-closure glaucoma. Therefore, PBF could be utilized as an indicator to assess the effect of PB on AC morphology and evaluate its risks. Furthermore, obvious tendency of increase in equivalent elastic strain and equivalent stress of iris could be observed with the increase of PBF. The equivalent elastic strain of iris represented its deformability, which could be utilized to evaluate the risk of iris bomb. Increased equivalent stress of iris under high PBF could hasten the detachment of epithelial cells in iris, which might cause an obstruction in AC and rise the resistance of AH outflow , verifying the reliability of PBF to assess the risks of BP.
According to clinical experience, obvious deformation of iris and narrowing of ACA could be observed when ILC distance decreased to 5 μm with a 3 mm pupil diameter. The value of PBF under that pathology condition was 45835.6 J/kg. Considering the effect of pupil diameter on PBF, the value should fluctuate in the range from 27849.2 J/kg to 53843.4 J/kg. As mentioned in Section 3.2, a quadratic polynomial equation has been solved to explain the variation in PBF, and a second-order empirical model has been developed to provide a good estimate of PBF. Based on the linear phenomena between PBF and ACA as well as PBF and De, it was reasonable to use the PBF value of 27000 J/kg as PB risk threshold in clinical. Further study will focus on the optimization of the empirical model of PBF to enlarge the range of pupil diameter.
This study presented a simple and practical method to assess PBF by the measurement of ILC distance and pupil diameter, which were accessible parameters in clinical. The results demonstrated that the ILC distance had an obvious effect on PBF. In a specific range of pupil diameter from 2.4 mm to 3.6 mm, the ACA became narrow with the decrease of the diameter. The quantitative analysis of PBF in this study could be used to research PACG pathological correlation and its pathogenesis, so as to provide a threshold for clinical diagnosis of PACG.
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
This work was supported by the National Natural Science Foundation of China (nos. 10802053 and 31570952) and Beijing Natural Science Foundation (no. 7152022).
J. van der Schoot, The Complexity of Diagnosing and Monitoring Glaucoma, Erasmus University Rotterdam, Rotterdam, Netherlands, 2017.
B. S. Shastry, “Genetic susceptibility to primary angle closure glaucoma (PACG),” Discovery Medicine, vol. 15, no. 80, pp. 17–22, 2013.View at: Google Scholar
J. Caprioli, “The ciliary epithelia and aqueous humor,” in Adlers Physiology of the Eye, pp. 228–247, Mosby, St. Louis, MO, USA, 1992.View at: Google Scholar
T. Kondo and M. Miura, “A method of measuring pupil-blocking force in the human eye,” Graefe’s Archive for Clinical and Experimental Ophthalmology, vol. 225, no. 5, pp. 361–364, 1987.View at: Google Scholar
R. Millar, “Flow resistance,” in Encyclopedia of Earth Science, pp. 476–478, Springer, Berlin, Heidelberg, Germany, 1978.View at: Google Scholar
R. Amini and V. H. Barocas, “Anterior chamber angle opening during corneoscleral indentation: the mechanism of whole eye globe deformation and the importance of the Limbus (vol. 50, pg. 5288, 2009),” Investigative Ophthalmology & Visual Science, vol. 54, p. 7323, 2013.View at: Google Scholar
V. H. Barocas and E. C. Huang, “Aqueous humor and iris mechanics in pupillary block and primary angle–closure glaucoma,” Investigative Ophthalmology & Visual Science, vol. 45, no. 13, p. 5033, 2004.View at: Google Scholar
J. D. Stitzel, S. M. Duma, J. M. Cormier, and I. P. Herring, “A nonlinear finite element model of the eye with experimental validation for the prediction of globe rupture,” Stapp Car Crash Journal, vol. 46, pp. 81–102, 2002.View at: Google Scholar
G. E. P. Box and N. R. Draper, Empirical Model-Building and Response Surfaces, Wiley, New York, NY, USA, 1987.
M. Dvoriashyna, R. Repetto, M. R. Romano, and J. H. Tweedy, “Aqueous humour flow in the posterior chamber of the eye and its modifications due to pupillary block and iridotomy,” Mathematical Medicine & Biology A Journal of the IMA, vol. 35, no. 4, pp. 447–467, 2017.View at: Publisher Site | Google Scholar