May 2006
Volume 47, Issue 13
ARVO Annual Meeting Abstract  |   May 2006
Large–Scale Parallel Finite Element Simulations of the Monkey Lamina Cribrosa Microarchitecture
Author Affiliations & Notes
  • S. Kodiyalam
    Biological Sciences, Louisiana State University, Baton Rouge, LA
  • R.T. Hart
    Biomedical Engineering, Tulane University, New Orleans, LA
  • C.F. Burgoyne
    Devers Eye Institute, Portland, OR
  • M.D. Roberts
    Devers Eye Institute, Portland, OR
  • J.C. Downs
    Devers Eye Institute, Portland, OR
  • Footnotes
    Commercial Relationships  S. Kodiyalam, None; R.T. Hart, None; C.F. Burgoyne, None; M.D. Roberts, None; J.C. Downs, None.
  • Footnotes
    Support  NIH–BRIN/INBRE Grant P20 RR16456, NIH Grant NEI R01–EY11610
Investigative Ophthalmology & Visual Science May 2006, Vol.47, 1229. doi:
  • Views
  • Share
  • Tools
    • Alerts
      This feature is available to authenticated users only.
      Sign In or Create an Account ×
    • Get Citation

      S. Kodiyalam, R.T. Hart, C.F. Burgoyne, M.D. Roberts, J.C. Downs; Large–Scale Parallel Finite Element Simulations of the Monkey Lamina Cribrosa Microarchitecture . Invest. Ophthalmol. Vis. Sci. 2006;47(13):1229.

      Download citation file:

      © ARVO (1962-2015); The Authors (2016-present)

  • Supplements

Purpose: : To compare IOP–related stress and strain within large–scale voxel–based models of a monkey optic nerve head (ONH) computed on a parallel machine to those generated by tetrahedra–based models run with commercial software.

Methods: : We have developed novel parallel finite element code to handle large–scale, voxel–based models of 3–D reconstructions of the lamina cribrosa. Computations are carried out for voxel–based laminar microstructure from a small region of a monkey ONH containing ∼100,000 laminar voxels. These results are compared to results from a tetrahedral mesh for the same region obtained from commercial software. Another small region of the actual nerve head spanning 75x77x41 voxels in xyz size is used to construct cylindrical ONH–like structures containing ∼400,000 laminar voxels that are studied under different boundary conditions. In order to balance the load on the CPUs of the parallel machine an automated load–balancing scheme has been developed and tested with an entire monkey lamina containing ∼7 million laminar voxels.

Results: : Results from voxel–based computations are consistent with those from the tetrahedral–based mesh run using commercial software: the best–case average relative error in the displacement components is ∼3.0%. For this case, while the relative error in the stress–strain components is 11–26%, the spatial correlation is very strong: > 0.95 for all components. For the ONH–like configurations, a linear trend is observed in the anterior–to–posterior variation of average hydrostatic pressure within the lamina (translaminar pressure gradient) under IOP loading and fixed anterior–to–posterior displacement at the cylindrical lamina boundary. When applied to the entire lamina configuration (∼7 million laminar voxels), the automated load–balancing scheme results in a small load imbalance of 6% on 24 CPU's.

Conclusions: : Quantitative agreement with the results from tetrahedra–based calculation achieves a first validation of our voxel–based methodology. The observed translaminar pressure gradient in the ONH–like configurations is similar to the trend in experimental measurements of the pressure within the lamina [Morgan W.H. et al, IOVS 1995, Vol 36, No. 6, 1163–72]. With a linear, IOP–induced, translaminar pressure gradient as the loading input load balanced computations will allow us to simulate the biomechanical behavior of the entire laminar microstructure in a single model.

Keywords: optic disc • lamina cribrosa 

This PDF is available to Subscribers Only

Sign in or purchase a subscription to access this content. ×

You must be signed into an individual account to use this feature.