A poroelastic model for the perfusion of the lamina cribrosa in the optic nerve head

In this work we present a mathematical model for the coupling between biomechanics and hemodynamics in the lamina cribrosa, a thin porous tissue at the base of the optic nerve head which is thought to be the site of injury in ocular neurodegenerative diseases such as glaucoma. In this exploratory two-dimensional investigation, the lamina cribrosa is modeled as a poroelastic material where blood vessels are viewed as pores in a solid elastic matrix. The model is used to investigate the influence on the distributions of stress, blood volume fraction (or vascular porosity) and blood velocity within the lamina cribrosa due to the application of different levels of the intraocular pressure (IOP) and the enforcement of different mechanical constraints at the lamina's boundary. The model simulations suggest that the degree of fixity of the boundary constraint strongly influences the lamina's response to IOP elevation. Specifically, when the boundary is mechanically clamped, IOP elevation leads to an increase in stress close to the lamina's boundary, making it more susceptible to tissue damage. On the other hand, when rotations are allowed at the boundary, the most vulnerable region appears to be located at the lamina's central axis, in proximity of the eye globe, where increased stress and reduced vascular porosity and blood velocity are predicted for increased levels of IOP.