×

3D low frequency electromagnetic modelling of the human eye with boundary elements: application to conductive keratoplasty. (English) Zbl 1244.78026

Summary: Conductive keratoplasty (CK) is a non-ablative surgical technique for the treatment of mild to moderate hyperopia (far-sightedness). In a CK session a thin electrode penetrates the cornea and delivers pulsed radio-frequency energy at 350 kHz to the surrounding tissue. The electromagnetic (EM) energy is dissipated into heat in the tissue surrounding the tip yielding thermally localised shrinkage and tightening of the collagen lamellae. When applied in a controlled way, this relatively new technique allows eye surgeons to correct the shape of the cornea and to treat common eye types of disease such as far-sightedness or astigmatism. The purpose of this work is to present a three-dimensional boundary element model of the human eye and its solution for the current density distribution in the different tissues when exposed to CK treatment. This is in order to examine in detail the induced currents that appear in the rest of the eye while delivering the EM signal, and to understand how the human eye behaves as a bulk imperfect conductor when a low frequency high voltage signal is locally applied. The Boundary element method is implemented with a multi-domain staggered approach in which each boundary element contains both continuous and discontinuous nodes. The former allocate electric potential only, while the latter allocate normal electric field only. In this way, continuity of the potential is ensured between adjacent elements and the complication of matching normal electric fields at corner points generated by elements of different orientations is avoided. A quantitative estimation of the EM energy absorbed in each tissue is presented, as well as the initial impedance of the electrodes.

MSC:

78M15 Boundary element methods applied to problems in optics and electromagnetic theory
78A70 Biological applications of optics and electromagnetic theory
92C99 Physiological, cellular and medical topics

Software:

LAPACK; GiD
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Poljak D. Human exposure to electromagnetic fields. Advances in electrical and electronic engineering. Southampton, UK: WIT Press; 2003.; Poljak D. Human exposure to electromagnetic fields. Advances in electrical and electronic engineering. Southampton, UK: WIT Press; 2003.
[2] Dodig, H.; Peratta, A.; Poljak, D., Analysis method for the heating of the human eye exposed to high frequency electromagnetic fields, J Commun Software Systems, 3, 1, 3-10 (2007)
[3] Jackson, J. D., Classical electrodynamics (1975), Wiley: Wiley NY, Chichester, Brisbane, Toronto, Singapore · Zbl 0997.78500
[4] Brebbia, C.; Telles, J.; Wrobel, L., Boundary elements techniques (1984), Springer: Springer Berlin, Heidelberg, New York, Tokio · Zbl 0556.73086
[5] Brebbia, C.; Dominguez, J., Boundary elements, an introductory course (1992), Computational Mechanics Publications, McGraw-Hill: Computational Mechanics Publications, McGraw-Hill New York, Colorado, San Francisco, Mexico, Toronto · Zbl 0780.73002
[6] Poljak, D.; Tham, C. Y.; Kovac, N., The assessment of human exposure to low frequency and high frequency electromagnetic fields using boundary element analysis, Eng Anal Boundary Elements, 27, 999-1007 (2003) · Zbl 1043.92003
[7] Gonzalez, M. C.; Peratta, A.; Poljak, D., Boundary element modeling of the realistic human body exposed to extremely-low-frequency (ELF) electric fields: computational and geometrical aspects, IEEE Trans Electromagn Compatibility, 49, 1, 153-162 (2007)
[8] Peratta, A.; Gonzalez, C.; Poljak, D., Current density induced in the human body due to power distributions lines using the boundary element method, J Commun Software Systems, 3, 1 (2007)
[9] Bottauscio O. Evaluation of induced currents in human body models exposed to ELF electromagnetic fields using BEM. In: Brebbia CA, editor. Boundary element research in Europe, European boundary element method symposium, vol. 2. Southampton, Boston: Computational Mechanics Publications; 1998. p. 125 [chapter 4].; Bottauscio O. Evaluation of induced currents in human body models exposed to ELF electromagnetic fields using BEM. In: Brebbia CA, editor. Boundary element research in Europe, European boundary element method symposium, vol. 2. Southampton, Boston: Computational Mechanics Publications; 1998. p. 125 [chapter 4]. · Zbl 0924.92001
[10] Gabriel, C.; Gabriel, S.; Corthout, E., The dielectric properties of biological tissues: I. Literature survey, Phys Med Biol, 41, 2231-2249 (1996)
[11] Gabriel, S.; Lau, R. W.; Gabriel, C., The dielectric properties of biological tissues: II. measurements in the frequency range 10Hz to 20gHz, Phys Med Biol, 41, 2251-2269 (1996)
[12] Peratta, A.; Popov, V., Modelling flow and solute transport in fractured porous media using the DRM multidomain technique, (Popov, V.; Power, H.; Škerget, L., Domain decomposition techniques for boundary elements (2007), WIT Press: WIT Press Southampton, Boston), 99-146, [chapter 4] · Zbl 1378.76112
[13] Taigbenu A. Simulation of coupled nonlinear electromagnetic heating with the green element method. In: Sunden B, Brebbia C, editors. Advanced computational methods in heat transfer IX, Engineering sciences, vol. 53. Wessex Institute of Technolgy, WIT Press; 2006.; Taigbenu A. Simulation of coupled nonlinear electromagnetic heating with the green element method. In: Sunden B, Brebbia C, editors. Advanced computational methods in heat transfer IX, Engineering sciences, vol. 53. Wessex Institute of Technolgy, WIT Press; 2006.
[14] Taigbenu, A., The green element method (1999), Kluwer: Kluwer Boston, USA · Zbl 0978.76061
[15] Mendez G, Mendez Noble A. Conductive keratoplasty for the correction of hyperopia, Baltimore: Williams and Wilkens; 1997. p. 163-71 [chapter 17].; Mendez G, Mendez Noble A. Conductive keratoplasty for the correction of hyperopia, Baltimore: Williams and Wilkens; 1997. p. 163-71 [chapter 17].
[16] Pallikaris, I.; Naoumidi, T.; Astyrakakis, N., Long-term results of conductive keratoplasty for low to moderate hyperopia, J Cataract Refractive Surg, 31 (2005)
[17] Lin, D.; Manche, E., Two-year results of conductive keratoplasty for the correction of low to moderate hyperopia, J Cataract Refract Surg, 29, 2339-2350 (2003)
[18] McDonald M, et al. Conductive keratoplasty for the correction of low to moderate hyperopia: US clinical trial 1-year results on 355 eyes, Am Acad Ophthalmol 2002;1978-89.; McDonald M, et al. Conductive keratoplasty for the correction of low to moderate hyperopia: US clinical trial 1-year results on 355 eyes, Am Acad Ophthalmol 2002;1978-89.
[19] Pearce J, Panescu D, Thomsen S. Simulation of diopter changes in radio frequency conductive keratoplasty. In: Ursino M, Brebbia C, Magosso E, editors. Modelling in medicine and biology VI, Advances in bioengineering, vol. 2. Southampton, UK: Wessex Institute of Technology, WIT Press; 2005. p. 469-77.; Pearce J, Panescu D, Thomsen S. Simulation of diopter changes in radio frequency conductive keratoplasty. In: Ursino M, Brebbia C, Magosso E, editors. Modelling in medicine and biology VI, Advances in bioengineering, vol. 2. Southampton, UK: Wessex Institute of Technology, WIT Press; 2005. p. 469-77.
[20] Berjano, E.; Alió, J.; Saiz, J., Modeling for radio-frequency conductive keratoplasty: implications for the maximum temperature reached in the cornea, Phys Meas, 26, 157-172 (2005)
[21] Smerdon, D., Anatomy of the eye and orbit, Curr Anaesthesia Crit Care, 11, 286-292 (2000)
[22] Hermans E, Dubbelman M, van der Heijde G, Heethaar R. Estimating the external force acting on the human eye lens during accommodation by finite element modelling. Vision research, 2006. Available online at \(\langle;\) www.sciencedirect.com \(\rangle;\); Hermans E, Dubbelman M, van der Heijde G, Heethaar R. Estimating the external force acting on the human eye lens during accommodation by finite element modelling. Vision research, 2006. Available online at \(\langle;\) www.sciencedirect.com \(\rangle;\)
[23] Lindemblatt, G.; Silny, J., A model of the electrical volume conductor in the region of the eye in the ELF range, Phys Med Biol, 46, 3051-3059 (2001)
[24] Barchanski, A.; Clemens, M.; De Gersem, H.; Weiland, T., Efficient calculation of current densities in the human body induced by arbitrarily shaped, low-frequency magnetic field sources, J Comput Phys, 214, 81-95 (2006) · Zbl 1086.92003
[25] Hirata, A.; Matsuyama, S.; Shiozawa, T., Temperature rises in the human eye exposed to EM waves in the frequency range 0.6-6GHz, IEEE Trans Electromagn Compatibility, 42, 4, 386-393 (2000)
[26] Pearce J, Panescu D. Radio frequency conductive keratoplasty in the cornea: prediction of diopter changes in numerical models. In: Proceedings of the 26th annual international conference of the IEEE EMBS, San Francisco, CA, USA, September 2004. p. 5426-8.; Pearce J, Panescu D. Radio frequency conductive keratoplasty in the cornea: prediction of diopter changes in numerical models. In: Proceedings of the 26th annual international conference of the IEEE EMBS, San Francisco, CA, USA, September 2004. p. 5426-8.
[27] GID resources. Website \(\langle;\) http://gid.cimne.upc.es \(\rangle;\); GID resources. Website \(\langle;\) http://gid.cimne.upc.es \(\rangle;\)
[28] CIMNE, International center for numerical methods in engineering, Barcelona, Spain. GID, The personal pre/postprocessor manual.; CIMNE, International center for numerical methods in engineering, Barcelona, Spain. GID, The personal pre/postprocessor manual.
[29] Ooi, E. H.; Ang, W. T.; Ng, E. Y.K., Bioheat transfer in the human eye: a boundary element approach, Eng Anal Boundary Elements, 31, 6, 494-500 (2007) · Zbl 1195.80031
[30] Anderson, E.; Bai, Z.; Bischof, C.; Blackford, S.; Demmel, J.; Dongarra, J.; Du Croz, J.; Greenbaum, A.; Hammarling, S.; McKenney, A.; Sorensen, D., LAPACK users’ guide (1999), SIAM: SIAM Philadelphia, PA · Zbl 0934.65030
[31] Choi, B.; Welch, A. J.; Pearce, J. A., Dynamic impedance measurements during radio-frequency heating of cornea, IEEE Trans Biomed Eng, 49, 6, 1610 (2002)
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.