Research Article | Open Access
Relationship between Corneal Temperature and Intraocular Pressure in Healthy Individuals: A Clinical Thermographic Analysis
Purpose. To study the geographical distribution of corneal temperature (CT) and its influence on the intraocular pressure (IOP) of healthy human volunteers. Materials and Methods. Fifteen subjects (7 M, 8 F), years old, were enrolled in this pilot, cross-sectional study. Measurements of CT were taken after one hour with closed eyelids (CET) or closed eyelids with a cooling mask (cm-CET) and compared to baseline. Results. If compared to baseline, after CET, average CT significantly increased by 0.56°C in the RE and by 0.48°C in the LE () and IOP concomitantly significantly increased by 1.13 mmHg and 1.46 mmHg, respectively, in each eye (). After cm-CET, average CT significantly decreased by 0.11°C and 0.20°C, respectively, in the RE and LE (RE ; LE ), followed by a significant IOP decrease of 2.19 mmHg and 1.54 mmHg, respectively, in each eye (RE ; LE ). Conclusion. Significant variations of CT occurred after CET and cm-CET and were directly correlated with significant differences of IOP. It can be speculated that both oxidative stress and sympathetic nerve fiber stimulation by temperature oscillations may affect the regulation of AH vortex flow and turnover, thus influencing IOP values.
Temperature is one of the fundamental regulators of tissue metabolism [1, 2]. Interest in the temperature of the eye spans almost 130 years and the ability to measure the temperature of the eye, driven by prevailing technologies, has potential importance in both research and clinical situations, including the study of ocular physiology and pathology [3–9].
New infrared ocular thermographs allow a noncontact and nonintrusive characterization of the thermal profile across the ocular surface [10–17]. Applications have included dry eye, wearing contact lens, corneal sensitivity, and ophthalmic surgery [18–24]. In addition, some studies showed a correlation between ocular surface temperature and ocular blood flow. An increase of intraocular pressure (IOP) was found to be related to a contemporary decrease of ocular perfusion pressure and ocular temperature in monkeys . A recent study on humans has also reported that eyes with ischemic central venous retinal occlusion (CRVO) have lower ocular surface temperatures than nonischemic ones . Moreover, in carotid artery stenosis, the eye on the affected side has been found to have an impairment in retrobulbar hemodynamics along with a reduction in corneal temperature (CT) . Thermography has also been applied to explore the role of vascular factors in the physiopathology of glaucoma and Galassi et al. have recently defined ocular surface temperature as a marker of impaired retrobulbar hemodynamics in patients with glaucoma .
However, to date, little work has been undertaken to determine the relationship between IOP and CT [3, 22, 28]. The variations in IOP following closed eyelid test (CET) both under normal conditions and after the administration of antioxidants have recently been investigated, leading to the conclusion that CET-induced ocular hypertension could be a response to mixed stress—oxidative and thermic—with degenerative effects on the trabecular meshwork (TM) [22, 29]. Moreover, given the known influence of ambient pressure and temperature on IOP, an underwater mask has been proposed as a provocative test in the diagnosis of primary open angle glaucoma (POAG) .
POAG is a multifactorial and not yet well-understood pathology [31, 32]. Ocular hypertension is one known critical risk factor for glaucomatous optic neuropathy. IOP is generated and maintained via the aqueous humor circulation system in the anterior chamber (AC) of the eye. The major factor controlling IOP is the dynamic balance between aqueous humor production in the ciliary body and its draining through so-called conventional—TM and Schlemm’s canal (SC)—and uveoscleral outflow pathways [33–36].
The goal of the study here reported was to map the CT in different areas of the cornea of healthy human volunteers and follow its variations after CET at room temperature or with a cooling mask (cm-CET) and correlate such variations to IOP values. The application of a physical model was then used to validate the experimental results.
2. Patients and Methods
This pilot, prospective, cross-sectional study was conducted following the tenets of the Declaration of Helsinki. Fifteen healthy Caucasian human volunteers (8 females and 7 males) were enrolled. The mean age was years (range 18–72). After signing an informed consent, each study participant underwent a complete ophthalmological examination, including a medical history review, best-corrected visual acuity measurement (BCVA), slit-lamp biomicroscopy, dry eye tests, and dilated fundus examination. All subjects were free from ocular and systemic diseases with no history of previous ocular surgery. Patients enrolled in the study were accepted if they had no signs or symptoms of ocular dryness: Ocular Protection Index (OPI: calculated as the ratio between BUT and the blinking frequency per minute) score ≥ 1 and OSDI score ≤ 12. In addition, since corneal pachymetry may influence temperature fluctuations and their determination, patients’ central corneal thickness measured by ultrasonic pachymetry (Pacline compact multifeature pachymeter, Optikon 2000, Rome, Italy) had to be within normality limits, that is, microns.
Further inclusion criteria were as follows:(i)IOP ≤ 21 mmHg, without any treatment.(ii)Refraction values between −4 and +4 spheric diopters.(iii)BCVA for far distance equal to 10/10.(iv)Normal angle structure at gonioscopy.(v)Normal C/D ratio at slit-lamp examination.(vi)Normal SAP (30-2 SITA standard program).
2.1. Corneal Temperature Measurement
Precise spatiotemporal measurements of CT were captured through the latest generation infrared thermometer (Sola Electro-Optics, Shanghai, China). Since it is known that there is an uneven distribution of temperature in the cornea, the average CT was calculated based on the recordings of CTs in five different areas of the cornea (nasal, temporal, superior, inferior, and central). All patients were acclimatized to the clinical environment for at least 15 min. The room temperature was specifically set and controlled at all times at 25°C (≈77°F), humidity was maintained at 42.0%, and the average indoor illumination was maintained at 300 lux, considered a standard indoor level of illumination. Any air drafts were avoided. The measurements were performed between 9:00 and 11:00 a.m. in a seated position and under the conditions described by Mori et al. : the subject blinked normally, then closed both eyes for 5 s, and then kept the eyes open for more than 10 s. The thermography device was set up 20 cm in front of the eye and the head was held steady with a frame. During that time, the subject was asked not to blink. If the subject blinked, a new measurement was performed.
Three measurements were taken consecutively during a single session for each eye and the average value was recorded. The measurements were repeated after one hour of CET and after the same test following the application of a cooling mask on the volunteers’ eyelids. To avoid any operator related bias, all thermographic measurements were made by one single examiner.
2.2. IOP Measurements
IOP is reported as the mean value of three consecutive readings of each eye by Goldmann applanation tonometer (AT-900, Haag Streit Diagnostics, Switzerland) registered between 9:00 and 11:00 a.m. to minimize the effect of daily variations. To avoid interexaminer and intertonometer variances, all IOP measurements were taken by the same trained resident.
2.3. Physical Model
To simulate the temperature distribution of the human eye, we adopted the physical model developed by Karampatzakis and Samaras  that solves the Pennes bioheat transfer equation  coupled with the incompressible Navier-Stokes equation of fluid dynamics.
According to such a physical model, the eye can be modeled by seven regions with different thermal properties—the cornea, the anterior chamber, the trabecular meshwork, the iris, the lens, the vitreous humor, and the sclera. The basal metabolic heat generation, as well as the blood perfusion, mainly occurs in the iris and in the sclera. This model also includes secretory inflow, drainage, and circulation of the AH.
2.4. Statistical Analysis
The results obtained were statistically analyzed by Student’s -test for paired samples. values below 0.05 were considered significant and denoted by one () or two () asterisks (SPSS V.19, IBM SPSS Statistics, USA).
We analyzed 30 eyes of 15 healthy volunteers. Figure 1 shows the distribution of basal values across the cornea. In each eye, the lowest temperature was observed at the temporal side, the highest temperature was observed at the nasal side, and intermediate values were observed along the corneal longitudinal axis (superior, central, and inferior). The same distribution was also maintained after CET and cm-CET.
CT mean basal values were °C in the right eye (RE) and °C in the left eye (LE). After CET, all the temperatures tended to increase with respect to the basal values (mean values: °C in the RE and °C in the LE) whereas after cm-CET all the temperatures showed a tendency to decrease (mean values: °C in the RE and °C in the LE). On average, after CET a highly significant increase of the corneal temperature of 0.56°C for the RE and 0.48°C for the LE ( for both eyes) was observed and after cm-CET there was a significant decrease of 0.11°C for the RE and 0.20°C in the LE (RE and LE for both eyes) (Figure 2(a)).
Correspondingly, we observed significant variations in IOP, which increased after CET by 1.13 mmHg in the RE and 1.46 mmHg in the LE (basal mean values: mmHg in the RE and mmHg in the LE; after CET mean values: mmHg in the RE and mmHg in the LE; for both eyes), whereas it significantly decreased after cm-CET by 2.19 mmHg in the RE and by 1.54 mmHg in the LE (basal mean values: mmHg in the RE and mmHg in the LE; after cm-CET mean values: mmHg in the RE and mmHg in the LE; RE and LE ) (Figure 2(b)).
Figure 3 shows that there is indeed a direct correlation between average CT and IOP, with a correlation coefficient scoring 0.74 in the RE and 0.94 in the LE.
Finally, the application of the physical simulation model by Karampatzakis and Samaras  and Pennes bioheat transfer equation  (Figure 4) shows the following:(a)The mean temperature of the eye is influenced by the ambient temperature: it increases by about 1.0°C when the ambient temperature rises from 15°C to 30°C. The data that we obtained at an environmental temperature of 25°C were 36.3°C with open eyelids and 36.8°C after CET, in agreement with the theoretical prediction.(b)The evaporation rate of the thin tear film on the cornea decreases its temperature, similarly to sweat evaporation on the skin; therefore, blocking the evaporation increases the overall surface temperature, as less heat is being exchanged between the cornea and the environment. Mathematical simulations estimate such a temperature rise at about 0.4°C, in agreement with the experimental data after CET, when the evaporation was stopped for some time, and consequently the CT increased by 0.56°C (RE) and 0.48°C (LE).
This study shows that, in normal eyes, following CET and cm-CET, IOP significantly changes in response to variations of CT. In agreement with previous studies [34–37], we also observed an overall increase in CT after CET. In addition, our study is the first to report an overall decrease of CT after application of a cooling mask on the lid surface (cm-CET). Moreover, a direct correlation was found between CT and IOP, with IOP values following the increase or the decrease of average surface CT values. To date, the physiological link between CT and IOP still remains somewhat inconclusive. Previous investigations demonstrated that eyelid closure, or the use of an underwater mask with open eyelids, can raise IOP because of an increase of local temperature . Further studies confirmed the relationship between local temperature, variation in IOP, and anterior chamber oxidative stress following CET [29, 39–41].
As previously reported, our data confirm that CT varies in response to tear film evaporation rate, which in turn is influenced by environmental temperature and blinking rate [19, 42–44]. However, our results also show a peculiar pattern of CT distribution across the corneal surface, not in line with previous findings [12, 14, 15, 45]. In fact, in all experimental settings (basal, after CET, and cm-CET), we detected the coolest area at the corneal temporal region while the warmest area was at the nasal side. Along the corneal vertical axis, the temperatures showed intermediate values. To explain and validate our data, the physical model by Karampatzakis and Samaras  and Pennes bioheat transfer equation  was applied and allowed us to conclude that AH flow depends on CT distribution. Numerical simulations show that the circulation of AH follows the temperature difference between the temporal and nasal side of the cornea and that an increase of AH flow from stagnation to fast vortices correlates with the increasing asymmetry between temporal and nasal temperatures, thus influencing the balance between production and outflow, and finally IOP values. Accordingly, CT differences between nasal and temporal sides were higher after cm-CET (0.80 and 0.96 for the RE and LE, resp.) than after CET (0.52 and 0.44 for the RE and LE, resp.). This may indeed suggest a faster flow with lower temperatures, favoring a higher discharge of AH, and therefore a lower IOP.
IOP depends on the dynamics of AH turnover, which in turn is a balance between secretion and excretion. Three mechanisms are involved in AH formation by the ciliary body: diffusion, ultrafiltration, and active secretion, the last being the major contributor to its formation. AH can be considered analogous to a blood surrogate as it provides nutrition, removes excretory products from metabolism, transports neurotransmitters, stabilizes the ocular structure, and contributes to the regulation of the homeostasis of the surrounding ocular tissues [33–36]. AH, as part of a vascular circulatory loop, returns to the blood flow by passing through the TM, a uniquely modified vessel wall interposed between the anterior chamber and SC or through the uveoscleral pathway [34, 36].
Our results suggest that temperature may affect AH balance through the regulation of its secretion, excretion, and flow dynamics. Secretion is enhanced by increasing temperatures, because of increased blood flow in the anterior segment of the eye due to vasodilation and upregulation of metabolic processes in the ciliary body. The opposite effect is therefore expected after decreasing the temperature, since it is known that low temperatures stimulate the sympathetic system inducing vasoconstriction. Excretion is also expected to be regulated by temperature. At high temperatures, there is an increase of oxidative stress and inflammatory processes due to increased metabolic activities. These events may negatively affect the draining through the TM, since it has been shown that high temperature increases the cellularity of the anterior chamber angle structures [33, 34]. The opposite is expected with low temperatures.
The application of the physical model by Karampatzakis and Samaras  confirmed the direct correlation between IOP and CT and the influence of tear film evaporation and environmental temperature on CT. Moreover, and more interestingly, such numerical simulation suggested that AH responds to increasing surface temperature generating vortices that might contribute to the variations in IOP observed in our experiments.
Finally, the presence of several transient receptor potential (TRP) channel isotypes that are responsive to temperature variations in the three corneal layers and in corneal nerve fibers  may suggest an involvement of this receptor type also in controlling the aqueous humor dynamics, although no reports exist yet indicating such kind of interaction.
IOP is the only treatable risk factor in glaucoma, one of the world’s leading causes of blindness. The existing correlation between CT and IOP suggests that CT may affect short-term IOP control and its fluctuations. Such IOP fluctuations are related to faster glaucoma progression. Therefore, controlling IOP and its fluctuations is a main therapeutic goal in glaucoma treatment. Temperature and oxidative stress remain additional targets in the control of glaucoma progression [29, 31, 34, 36, 40]. It will be interesting to extend these findings studying the influence of temperature on IOP not only in the seated position but also in the supine position and on a larger population. Moreover, it could be interesting to see whether the same correlation persists in the presence of ocular hypertension. In such case, the use of a cooling mask during the night, when peaks of IOP are mostly expected with negative effects on POAG progression, might be advisable.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
English proofreading by Dr. Antony Bridgewood is greatly appreciated.
- E. F. Ring, “Progress in the measurement of human body temperature,” IEEE Engineering in Medicine and Biology Magazine, vol. 17, no. 4, pp. 19–24, 1998.
- B. F. Jones, “A reappraisal of the use of infrared thermal image analysis in medicine,” IEEE Transactions on Medical Imaging, vol. 17, no. 6, pp. 1019–1027, 1998.
- C. Purslow and J. Wolffsohn, “The relation between physical properties of the anterior eye and ocular surface temperature,” Optometry and Vision Science, vol. 84, no. 3, pp. 197–201, 2007.
- M. Sniegowski, M. Erlanger, R. Velez-Montoya, and J. L. Olson, “Difference in ocular surface temperature by infrared thermography in phakic and pseudophakic patients,” Clinical Ophthalmology, vol. 9, pp. 461–466, 2015.
- J. Alio and M. Padron, “Normal variations in the thermographic pattern of the orbito-ocular region,” Diagnostic Imaging, vol. 51, no. 2, pp. 93–98, 1982.
- J. Alio and M. Padron, “Influence of age on the temperature of the anterior segment of the eye. Measurements by infrared thermometry,” Ophthalmic Research, vol. 14, no. 3, pp. 153–159, 1982.
- S.-R. Shih, H.-Y. Li, Y.-L. Hsiao, and T.-C. Chang, “The application of temperature measurement of the eyes by digital infrared thermal imaging as a prognostic factor of methylprednisolone pulse therapy for Graves' ophthalmopathy,” Acta Ophthalmologica, vol. 88, no. 5, pp. e154–e159, 2010.
- T.-C. Chang, Y.-L. Hsiao, and S.-L. Liao, “Application of digital infrared thermal imaging in determining inflammatory state and follow-up effect of methylprednisolone pulse therapy in patients with Graves' ophthalmopathy,” Graefe's Archive for Clinical and Experimental Ophthalmology, vol. 246, no. 1, pp. 45–49, 2008.
- U. R. Acharya, E. Y. K. Ng, G. C. Yee, T. J. Hua, and M. Kagathi, “Analysis of normal human eye with different age groups using infrared images,” Journal of Medical Systems, vol. 33, no. 3, pp. 207–213, 2009.
- C. Purslow and J. S. Wolffsohn, “Ocular surface temperature: a review,” Eye and Contact Lens, vol. 31, no. 3, pp. 117–123, 2005.
- P. B. Morgan, M. P. Soh, and N. Efron, “Corneal surface temperature decreases with age,” Contact Lens and Anterior Eye, vol. 22, no. 1, pp. 11–13, 1999.
- N. Efron, G. Young, and N. A. Brennan, “Ocular surface temperature,” Current Eye Research, vol. 8, no. 9, pp. 901–906, 1989.
- M. K. J. Klamann, A.-K. B. Maier, J. Gonnermann, J. P. Klein, and U. Pleyer, “Measurement of dynamic ocular surface temperature in healthy subjects using a new thermography device,” Current Eye Research, vol. 37, no. 8, pp. 678–683, 2012.
- M. Pattmöller, J. Wang, J. Pattmöller et al., “Interobserver and intraobserver reliability of corneal surface temperature measurements with the TG-1000 thermograph in normal eyes,” Ophthalmologe, vol. 112, no. 9, pp. 746–751, 2015 (German).
- S. Moussa, T. Eppig, J. Pattmöller et al., “Diurnal and zonal analysis of corneal surface temperature in young healthy adults,” European Journal of Ophthalmology, vol. 23, no. 5, pp. 641–646, 2013.
- R. Mapstone, “Measurement of corneal temperature,” Experimental Eye Research, vol. 7, no. 2, 237-IN29 pages, 1968.
- L. Tan, Z.-Q. Cai, and N.-S. Lai, “Accuracy and sensitivity of the dynamic ocular thermography and inter-subjects ocular surface temperature (OST) in Chinese young adults,” Contact Lens and Anterior Eye, vol. 32, no. 2, pp. 78–83, 2009.
- A. Mori, Y. Oguchi, Y. Okusawa, M. Ono, H. Fujishima, and K. Tsubota, “Use of high-speed, high-resolution thermography to evaluate the tear film layer,” American Journal of Ophthalmology, vol. 124, no. 6, pp. 729–735, 1997.
- T. Y. Su, W. T. Ho, C. Y. Lu, S. W. Chang, and H. K. Chiang, “Correlations among ocular surface temperature difference value, the tear meniscus height, Schirmer's test and fluorescein tear film break up time,” British Journal of Ophthalmology, vol. 99, no. 4, pp. 482–487, 2015.
- S. Buschschlüter, C. Koch, J. von Eicken, and H. Höh, “Computation of the temperature rise at the corneal endothelium during cataract surgery by modeling of heat generation inside the anterior chamber,” Ultrasound in Medicine and Biology, vol. 40, no. 10, pp. 2431–2444, 2014.
- R. Mencucci, C. Mazzotta, F. Rossi et al., “Riboflavin and ultraviolet A collagen crosslinking: in vivo thermographic analysis of the corneal surface,” Journal of Cataract and Refractive Surgery, vol. 33, no. 6, pp. 1005–1008, 2007.
- F. Galassi, B. Giambene, A. Corvi, G. Falaschi, and U. Menchini, “Retrobulbar hemodynamics and corneal surface temperature in glaucoma surgery,” International Ophthalmology, vol. 28, no. 6, pp. 399–405, 2008.
- G. Singh and H. S. Bhinder, “Comparison of noncontact infrared and remote sensor thermometry in normal and dry eye patients,” European Journal of Ophthalmology, vol. 15, no. 6, pp. 668–673, 2005.
- D. K. Martin and I. Fatt, “The presence of a contact lens induces a very small increase in the anterior corneal surface temperature,” Acta Ophthalmologica, vol. 64, no. 5, pp. 512–518, 1986.
- C. R. Auker, L. M. Parver, T. Doyle, and D. O. Carpenter, “Choroidal blood flow. I. Ocular tissue temperature as a measure of flow,” Archives of Ophthalmology, vol. 100, no. 8, pp. 1323–1326, 1982.
- A. Sodi, B. Giambene, G. Falaschi et al., “Ocular surface temperature in central retinal vein occlusion: preliminary data,” European Journal of Ophthalmology, vol. 17, no. 5, pp. 755–759, 2007.
- F. Galassi, B. Giambene, A. Corvi, and G. Falaschi, “Evaluation of ocular surface temperature and retrobulbar haemodynamics by infrared thermography and colour Doppler imaging in patients with glaucoma,” British Journal of Ophthalmology, vol. 91, no. 7, pp. 878–881, 2007.
- A. Bill, “Autonomic nervous control of uveal blood flow,” Acta Physiologica Scandinavica, vol. 56, pp. 70–81, 1962.
- N. Pescosolido, R. Malagola, G. Scarsella, F. Lenarduzzi, L. Dapoto, and M. Nebbioso, “Oxidative stress in the closed-eyelid test: management of glaucoma,” European Review for Medical and Pharmacological Sciences, vol. 16, no. 10, pp. 1453–1457, 2012.
- N. Pescosolido, “La maschera subacquea e' veramente innocua? Effetti sulla pressione endoculare ed utilita' come test di provocazione per la diagnosi nel glaucoma cronico semplice,” Annali di Ottalmologia e Clinica Oculistica, vol. 107, no. 9, pp. 659–666, 1981 (Italian).
- S. C. Saccà, A. Izzotti, P. Rossi, and C. Traverso, “Glaucomatous outflow pathway and oxidative stress,” Experimental Eye Research, vol. 84, no. 3, pp. 389–399, 2007.
- H. Quigley and A. T. Broman, “The number of people with glaucoma worldwide in 2010 and 2020,” British Journal of Ophthalmology, vol. 90, no. 3, pp. 262–267, 2006.
- S. C. Saccà, M. Centofanti, and A. Izzotti, “New proteins as vascular biomarkers in primary open angle glaucomatous aqueous humor,” Investigative Ophthalmology & Visual Science, vol. 53, no. 7, pp. 4242–4253, 2012.
- S. C. Saccà, A. Pulliero, and A. Izzotti, “The dysfunction of the trabecular meshwork during glaucoma course,” Journal of Cellular Physiology, vol. 230, no. 3, pp. 510–525, 2015.
- A. Izzotti, A. Bagnis, and S. C. Saccà, “The role of oxidative stress in glaucoma,” Mutation Research—Reviews in Mutation Research, vol. 612, no. 2, pp. 105–114, 2006.
- A. Bagnis, A. Izzotti, M. Centofanti, and S. C. Saccà, “Aqueous humor oxidative stress proteomic levels in primary open angle glaucoma,” Experimental Eye Research, vol. 103, pp. 55–62, 2012.
- A. Karampatzakis and T. Samaras, “Numerical model of heat transfer in the human eye with consideration of fluid dynamics of the aqueous humour,” Physics in Medicine and Biology, vol. 55, no. 19, pp. 5653–5665, 2010.
- H. H. Pennes, “Analysis of skin, muscle and brachial arterial blood temperatures in the resting normal human forearm,” The American Journal of the Medical Sciences, vol. 215, no. 3, article 354, 1948.
- N. Pescosolido, G. Belcaro, D. Rusciano, R. D. Steigerwalt Jr., and M. Nebbioso, “Retrospective study of glaucoma and closed-eyelid test: long-term outcomes in an Italian native population,” Panminerva Medica, vol. 55, pp. 1–7, 2013.
- N. Pescosolido, “‘The oxydative-stress’ role of free radicals in oxygen in hypertension which follow the ‘closed eyelid test’,” Annali di Ottalmologia e Clinica Oculistica, vol. 118, no. 8, pp. 867–876, 1992 (Italian).
- M. G. Bucci and N. Pescosolido, “The closed eyelid test in the management of glaucoma,” Glaucoma, vol. 5, no. 2, pp. 84–89, 1983.
- R. D. Freeman and I. Fatt, “Environmental influences on ocular temperature,” Investigative Ophthalmology, vol. 12, no. 8, pp. 596–602, 1973.
- J.-H. Tan, E. Y. K. Ng, U. R. Acharya, and C. Chee, “Infrared thermography on ocular surface temperature: a review,” Infrared Physics and Technology, vol. 52, no. 4, pp. 97–108, 2009.
- R. Mapstone, “Determinants of corneal temperature,” British Journal of Ophthalmology, vol. 52, no. 10, pp. 729–741, 1968.
- E. Y. K. Ng and E. H. Ooi, “Ocular surface temperature: a 3D FEM prediction using bioheat equation,” Computers in Biology and Medicine, vol. 37, no. 6, pp. 829–835, 2007.
- P. S. Reinach, W. Chen, and S. Mergler, “Polymodal roles of transient receptor potential channels in the control of ocular function,” Eye and Vision, vol. 2, no. 1, article 5, 10 pages, 2015.
Copyright © 2016 Claudia Fabiani et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.