No abstract

The mathematical model described appears to provide a framework for further development, capturing the essential features of mechanical anisotropy of the cornea. The tunnel incision simulation indicated the importance of the anisotropy in this case.

SUMMARYIn this paper a Galerkin least-squares (GLS) finite element method, in which residuals in least-squares form are added to the standard Galerkin variational equation, is developed to solve the Helmholtz equation in two dimensions. An important feature of GLS methods is the introduction of a local mesh parameter that may be designed to provide accurate solutions with relatively coarse meshes. Previous work has accomplished this for the one-dimensional Helmholtz equation using dispersion analysis. In this paper, the selection of the GLS mesh parameter for two dimensions is considered, and leads to elements that exhibit improved phase accuracy. For any given direction of wave propagation, an optimal GLS mesh parameter is determined using two-dimensional Fourier analysis. In general problems, the direction of wave propagation will not be known a priori. In this case, an optimal GLS parameter is found which reduces phase error for all possible wave vector orientations over elements. The optimal GLS parameters are derived for both consistent and lumped mass approximations. Several numerical examples are given and the results compared with those obtained from the Galerkin method. The extension of GLS to higher-order quadratic interpolations is also presented.

These results match our previous observation of a depth-dependent gradient in stromal collagen interconnectivity in the central cornea, and show that this gradient extends from the central cornea to the limbus. The lack of a preferred distribution of angled fibers with regard to corneal quadrant or radial position likely serves to evenly distribute loads and to avoid the formation of areas of stress concentration.

No abstract