Thermal Effects of Eyelid in Human Eye Temperature Model

J. Appl. Math. & Informatics Vol. 32(2014), No. 5 - 6, pp. 649 - 663 http://dx.doi.org/10.14317/jami.2014.649 THERMAL EFFECTS OF EYELID IN HUMAN EYE TEMPERATURE MODEL† K. C. GOKUL∗ , D. B. GURUNG AND P. R. ADHIKARY Abstract. Presence of eyelid on anterior ocular surface and its thermal effects play significant role in maintaining eye temperature. In most of the literatures of thermal modeling in human eye, the eyelid is not considered as an eye component. In this paper, finite element model is developed to investigate the thermal effects of eyelid closure and opening in human eye. Based on different properties and parameter values reported in literatures, the bio-heat transfer process is simulated and compared with experimental results in steady and transient state cases. The sensitivity analysis using various ambient temperatures, evaporation rates, blood temperatures and lens thermal conductivities is carried out. The temperature values so obtained in open eye show a good agreement with past results. The closure of eyelid is found to increase/decrease the eye temperature significantly than its opening, when the parameter values are considered to be at extreme. AMS Mathematics Subject Classification : 92B05, 80A20, 80M10. Key words and phrases : Bio-heat transfer, Eyelid, Finite Element Method. 1. Introduction Temperature is of vital importance for all biological functions. The rate of metabolic process increases with temperature. Temperature changes can affect tissues in several ways; it can kill cells, denature proteins, slow down or speed up metabolism, involve in pathological changes etc. Due to the lack of sufficient blood flow in the inner part and lack of skin as a protecting layer in anterior part, human eye is assumed to be the most vulnerable organ in human body, even for small thermal interactions. Exposure to elevated ambient temperature has been proposed to be a risk factor for presbyopia and cataract and low temperature for dry eye, eyelid spams, excess tearing etc[10]. Various environmental conditions like airflow, temperature, humidity and thermal radiation can influence on ocular temperature. Received November 22, 2013. Revised April 25, 2014. Accepted May 10, 2014. author. † This work was supported by the University Grants Commission Nepal. c 2014 Korean SIGCAM and KSCAM. ⃝ 649 ∗ Corresponding 650 K.C. Gokul, D.B. Gurung and P.R. Adhikary The human eye is relatively a small and complex organ, consists of several sub domains with different material properties and having complex geometry. Direct in vivo temperature monitoring in the interior of a live human eye is not yet completely possible. Invasive or direct contact techniques in measuring the eye temperature is now confined to animal experiments due to its damaging nature of test procedures[15]. Bolometer was used to measure corneal surface temperature in early 1960s. Mapstone[14] used bolometer to measure corneal temperature and obtained a mean temperature 34.8◦ C. Later on, the infrared radiation thermometry has been used in measuring corneal surface temperature ever since. Fujishima et al.[5] reported mean corneal surface temperature 34.2◦ C using infrared radiation technique. All these non invasive techniques measure corneal surface temperature only. The information of corneal surface temperature is not sufficient to predict intraocular temperature variations in case of medical treatment processes like hyperthermia, laser etc. Therefore, mathematical modeling is the alternative and appropriate approach to obtain the temperature distribution within the components of human eye. Lagendijk[11] used a finite difference method to calculate the temperature distribution in human and rabbit eyes during hyperthermia treatment. The heat transport from the sclera to the surrounding anatomy is described by a single heat transfer coefficient which includes the impact of blood flow in choroid and sclera. Scott[18] utilized finite element method to obtain the temperature profile based on heat conduction using various heat transfer coefficients given by Lagendijk. Flyckt et al.[4] studied the impact of choroidal blood flow and scleral convection on heat transfer coefficient in human eye. Ng and Ooi[15] presented a 2D finite element model and simulated ocular surface temperature using bioheat equation. Ooi et al.[16] studied the effect of aqueous humor hydrodynamics on heat transfer within human eye. Li et al.[13] studied the bio-heat transfer in the human eye using 3D alpha finite element method. Shafai and Vafai[19] studied the eye response to thermal disturbances and analyzed the role of primary thermal transport mechanisms on the eye subject to different conditions. All these past models have neglected the effects of eyelid in temperature distribution in human eye. Some authors[16, 18] acknowledged their deficiency of the model is excluding the eyelid as an eye component. Eyelid is closed while sleeping. During sleep, with the eyelid closed, the temperature of the anterior segment is increased about 2◦ C by the blood flow in the eyelid, although surface body temperature also declines with sleep[21]. Even in wake up, eyelid covers cornea surface for 7 seconds in a minute (on an average) while blinking, but it varies by individual and depends on different environmental and health conditions. When the eyes are focused on an object for long time, eyelid closure time decreases and when the eye is exposed on high speed wind flow, extreme hot/cold temperatures etc the eyelid closure time increases. Although, eyelid closure time is much shorter than opening, however, this change may contribute to the temperature of ocular surface[5]. Thermal Effects of Eyelid in Human Eye Temperature Model 651 There are four types of temperature effects on anterior corneal surface: heat loss/gain between blinks, heat transfer between cornea and environment, heating/cooling due to spread of tear across the surface of cornea, and the heating/cooling effect caused by the movement of the eyelid. Closure and opening of the eyelid show different variations in corneal temperature. During the opening of eyelid, the cornea surface temperature is controlled by convection, radiation and tear evaporation. When eyelid is closed these heat losses are prevented and the cornea is influenced by another thermal environment, heat conduction from vascular palpebral conjunctivamapstone1968. Therefore, modeling of heat transport is needed in order to completely explain the actual temperature variation in human eye in presence of eyelid. In the present paper, finite element method is used as a tool to predict the steady and transient temperature distribution in human eye in case of eyelid closure and opening. The eye temperature along pupillary axis using different values of evaporation rates, blood temperatures, ambient temperatures and lens thermal conductivities are studied. The temperature values obtained are also compared with past results. 2. Model Formulation 2.1. Eye geometry and properties. A schematic diagram of the human eye when eyelid is opened is presented in figure 1. For modeling purpose, the opened eye is divided into six regions: cornea, aqueous humor, lens, vitreous humor, retina, and sclera. In several literatures[20, 19, 18, 15, 16, 9] choroid/retina/sclera were modeled as a single layer as sclera. We separated choroid/retina from sclera and assumed as a different layer for three reasons: (1) Sclera is avascular but choroid/retina is vascular[3] (2) The blood perfusion rate in the choroid is the highest among any other tissue[9] (3) The metabolic rate is the highest in the retina among any other tissue[2] The diameter of the eye along pupillary axis is about 25.10mm[7]. The thickness of cornea, aqueous humor, lens, vitreous humor, retina and sclera have been considered as l1′ , l2′ − l1′ , l3′ − l2′ , l4′ − l3′ , l5′ − l4′ and l6′ − l5′ respectively. Similarly, T0′ , T1′ , T2′ , T3′ , T4′ and T5′ are the nodal temperatures at a distances x = 0, x = l1′ , x = l2′ , x = l3′ , x = l4′ and x = l5′ respectively and T6′ = Tc is the body core temperature at x = l6′ . Similarly, a schematic diagram of the eye when eyelid is closed is displayed in fig. 2. The eyelid consists of four layers: skin, orbicularis muscle, tarsus plate and palpebral conjunctiva. Palpebral conjunctiva is very thin compare to other layers, so we model palpebral conjunctiva and tarsal plate as a single layer. Hence, the eyelid is further divided into three layers: dermis (skin surface), orbicularis oculi muscle and tarsal plate. Orbicularis oculi is a kind of skeletal muscle that is responsible for closing of eyelid. Tarsal plate, considered to be skeleton of eyelid, consists not of cartilage, but of dense fibrous tissue that 652 K.C. Gokul, D.B. Gurung and P.R. Adhikary Figure 1. Schematic diagram of human eye when eyelid is opened Figure 2. Schematic diagram of human eye when eyelid is closed gives the eyelid its contour, and functions as stiffener. The eyelid thickness along pupillary axis is 4.25mm[12]. The eye in this case is divided into nine regions: dermis, orbicularis oculi, tarsal plate, cornea, aqueous humor, lens, vitreous humor, retina, and sclera. Therefore, the diameter of the closed eye along pupillary axis is about 29.35mm. The thickness of dermis, orbicularis oculi, tarsal plate, cornea, aqueous humor, lens, vitreous humor, retina and sclera have been considered as l1 , l2 − l1 , l3 − l2 , l4 − l3 , l5 − l4 , l6 − l5 , l7 − l6 , l8 − l7 and l9 − l8 respectively. Similarly, T0 , T1 , T2 , T3 , T4 , T5 , T6 , T7 , T8 and T9 = Tc (body core temperature) are the nodal temperatures at a distances x = 0, x = l1 , x = l2 , x = l3 , x = l4 , x = l5 , x = l6 , x = l7 , x = l8 and x = l9 . Due to the lack of appropriate data, the parameter values of cartilage are used for tarsal plate. The parameter values for different parts of eye are presented in table 1. 2.2. Governing equation and boundary conditions. The governing differential equation representing the bio-heat transfer in the human eye can be written by the well known Pennes equation addressing the effect of blood perfusion and metabolism[17] is given by: ρc ∂T = ∇.(K∇T ) + ωρb cb (Tb − T ) + Qm + Q ∂t (1) ◦ where, ρb = blood density (Kgm−3 ), cb = blood specific heat (JKg −1 C −1 ), k ◦ = tissue thermal conductivity (W m−1 C −1 ), ω = volumetric blood perfusion −1 rate per unit volume (s ), Tb = blood temperature (◦ C), T = tissue temperature (◦ C), Qm = heat generation due to metabolism (W m−3 ) and Q = heat generation due to external heat source(W m−3 ). Boundary conditions for the system can be defined as follows: Thermal Effects of Eyelid in Human Eye Temperature Model 653 Table 1. Thermal properties of human eye tissues Tissue Type Thermal Conductivity K o (W m−1 C −1 ) Dermis 0.34[4] Orbicularis 0.56[4] Tarsal 0.47[9] Cornea 0.58[15] Aqueous 0.58[15] Lens 0.40[15] Vitreous 0.603[15] Retina 0.565[3] Sclera 1.0042[15] Blood Metabolic Perfusion Rate ω Qm (s−1 ) (W m−3 ) 0.0087[4] 1620[9] 0.0034[4] 480[9] 0.0082[9] 1600[9] 0[4] 0[9] 0[4] 0[9] 0[4] 0[9] 0[4] 0[9] 0.0222[4] 22000[9] 0[3] 0[3] Density ρ (Kgm−3 ) 1070[4] 1050[4] 1250[9] 1050[15] 996[15] 1050[15] 1000[15] 1050[3] 1100[15] Specific heat C o (JKg −1 C −1 ) 3662[4] 3639[4] 3600[9] 4178[15] 3997[15] 3000[15] 4178[15] 3680[3] 3180[15] (1) On the outer surface of the sclera, the heat flows run into the eye with the complicated network of ophthalmic vessels which are located inside the choroidal layer acting as a heating source to the sclera. This heat exchange between the eye and the surrounding is modeled using the following convection boundary condition: ∂T = hb (T − Tb ) (2) Γ2 : −ks ∂η where η is the normal direction to the surface boundary, ks is the thermal conductivity of sclera, hb is the heat transfer coefficient between blood ◦ and eye (W m−2 C −1 ), and Tb is blood temperature (◦ C). (2) Since outer surface of the eye (cornea or skin) is exposed to the environment, the heat loss caused via convection, radiation, and evaporation. This loss is modeled using the following boundary condition : ∂T = ha (T − Ta ) + σϵ(T 4 − Ta4 ) + E (3) Γ1 : −kc ∂η { hc , When eyelid is opened where ha = , hc represents heat transhs , When eyelid is closed fer coefficient between environment and cornea and hs represents heat ◦ transfer coefficient between skin and environment (W m−2 C −1 ), Ta is the ambient temperature (◦ C), σ is the Stefan Boltzmann constant ◦ (5.67 × 10−8 W m−2 C −4 ), ϵ is the emissivity, E is evaporative heat loss (W m−2 ). The nonlinear radiation term in the boundary condition (3) is treated by using simple iterative procedure as follows: ∂T1 = [ha + σϵ(T1 + Ta )(T12 + Ta2 )](T1 − Ta ) + E (4) −kc ∂η 654 K.C. Gokul, D.B. Gurung and P.R. Adhikary −kc ∂T1m ∂η hcr hcr = hcr (T1m − Ta ) + E (5) = ha + σϵ(T1m−1 + Ta )((T1m−1 )2 + Ta2 ) (6) = hconvection + hradiation (7) where T1m are temperature sequences for m = 1, 2, 3, and T10 represents an initial guess of temperature. The iteration is completed when the convergent condition is satisfied: (8) ∥T1m − T1m−1 | < δ where δ is iteration tolerance. The inner body core temperature Tc is assumed to be 37◦ C. Therefore, the initial boundary condition is Tc = 37◦ C (9) The partial differential equation (1) together with boundary conditions (2) and (5) in one dimensional variational form is: ∫ ∂T 2 1 L dT ]dx I = [K( )2 + ωρb cb (Tb − T )2 − 2Qm T + ρc 2 0 dx ∂t 1 1 + hb (T − Tb )2 + hcr (T − Ta )2 + ET (10) 2 2 To optimize I, we differentiate I partially with respect to Ti and equating to zero as follows, ∂I =0 (11) ∂Ti Equation (11) is the system of linear equations which can be written in matrix form as [C]{Ṫ } + [K]{T } = {R} (12) i where {Ṫ } = { ∂T ∂t }, {T } = {Ti } and {R} = {Ri } are 6 × 1 and 9 × 1 vectors in case of eyelid opening and closure respectively. Similarly, [C] and [K] are 6 × 6 and 9 × 9 matrices in case of eyelid opening and closure respectively. Now, to get the steady state temperature distribution, we remove the time derivative from (12) and solve the following system of equation: [K]{T } = {R} (13) To get the transient temperature distribution, we apply Crank-Nicolson method to solve the system (12) with respect to time 1 1 1 1 [C] + [K]){T }n+1 = ( [C] − [K]){T }n + {R} (14) ∆t 2 ∆t 2 where ∆t is time interval. The temperature increases from outer surface of the eye towards core, when ambient temperature is less than 37◦ C and vice versa. Hence, we consider the temperature increases/decreases in linear order towards eye core with regard to ( Thermal Effects of Eyelid in Human Eye Temperature Model 655 thickness. For initial nodal temperatures {T }0 at time t = 0, we assume the following initial condition T (x = li , t = 0) = T (0, 0) + rli (15) Where i ranges from 1, 2, 3, · · · , 6 and 1, 2, 3, · · · , 9 in case of eyelid opening and closure respectively, T (0, 0) = 20◦ C and r = constant to be determined. The equation(13) and (14) are solved to get the required nodal temperatures in steady and transient state cases. 3. Results 3.1. Sensitivity analysis. In this section, the influences of various parameters on eye temperature are studied in case of eyelid closure and opening. These parameters are ambient temperatures, blood temperatures, evaporation rates and lens thermal conductivities. The aim is to determine the dominant affecting parameters on temperature distribution in human eye. For analysis, ◦ ◦ ◦ hs = 6.28W m−2 C −1 [18], hc = 10W m−2 C −1 [17], hb = 65W m−2 C −1 [17] and the parameter values in table 1 are used. 3.1.1. Effects of evaporation rates. When eyelids are opened the tear evaporation occurs from corneal surface. On the surface of corneal epithelium, there is usually a thin lipid layer produced by the Meibomian glands of the tarsal plate[18]. The function of this layer is to prevent evaporation of tear from the corneal surface[22]. When this layer is destroyed, the evaporation rate increases dramatically and can reach as high as 320W m−2 . In normal condition evaporation rate ranges from 20W m−2 to 100W m−2 [18]. Therefore, in this study E equals 0W m−2 and 120W m−2 are chosen for extreme cold and hot environment and 40W m−2 for room temperature respectively. When the eyelids are closed, the sweat evaporation occurs from eyelid skin surface. Although the skin layer in the eyelid is thinnest within whole body, same parameter values as in normal skin are used for analysis. A resting man doesn’t sweet when atmospheric temperature is below 20◦ C. A man at 21◦ C can lose 25% of his/her heat by sweat evaporation[8]. Thus, in this study E equals 0W m−2 and 192W m−2 are chosen for extreme cold and hot environment and 96W m−2 for room temperature respectively[1]. The temperature values so obtained are shown in figure 3. From figure 3 it is observed that the corneal surface temperature is dropped by 3.69◦ C and 4.75◦ C respectively when eyelids are opened and closed. The corneal surface temperature difference between eyelid closure and opening are 1.73, 0.67◦ C and 0.58◦ C respectively at low, average and high evaporation rates. The corneal temperature is dropped from 35.03◦ C to 30.28◦ C and from 33.30◦ C to 29.61◦ C, when the evaporation rate varies from low to high. The results show that when evaporation rate increases corneal temperature decreases but decreasing rate is slower in eyelid opening than closure. 656 K.C. Gokul, D.B. Gurung and P.R. Adhikary Figure 3. Temperature variation for different evaporation rates 3.1.2. Effects of ambient temperatures. Heat loss/gain from corneal surface depends on the temperature difference between the eye tissue (cornea or skin) and the surrounding air. This loss/gain occurs due to conduction, convection, radiation and evaporation. About 3%, 15%, 60% and 22% of total heat is lost by conduction, convection, radiation and evaporation respectively from human body[8]. In this study, the ambient temperatures 0◦ C and 50◦ C are used to study the thermal response to extreme cold and extreme hot environments and 25◦ C for normal. The numerical results are shown in figure 4. Figure 4. Temperature variation for different ambient temperatures Thermal Effects of Eyelid in Human Eye Temperature Model 657 It is observed from figure 4 that corneal temperature is increased from 24.37◦ C to 39.78◦ C and from 30.11◦ C to 38.76◦ C respectively in case of eyelid opening and closure. That is, an increase of 8.65◦ C and 15.41◦ C is observed at corneal surface when ambient temperature varies from 0◦ C to 50◦ C. The temperature difference at cornea between eyelid closure and opening are 5.74, 2.36◦ C and 1.02◦ C at ambient temperatures 0, 25◦ C and 50◦ C respectively. The eye is heated, when environmental temperature drops below normal level and vice versa by blood flow in eyelid. 3.1.3. Effects of blood temperatures. The blood flow in the choroid/retina/iris acts as the major source of heat for opened eye. In addition, blood flow in skin/orbicularis oculi/tarsal plate of eyelid acts as a heating source for anterior eye. Although, blood flow occurs only in a few parts, its temperature significantly affects the temperature throughout the human eye. The blood temperatures 35◦ C and 39◦ C are used to study the thermal response due to sickness and 37◦ Cfor normal case. The numerical result is presented in figure 5. Figure 5. Temperature variation for different blood temperatures It is observed from figure 5 that corneal temperature is increased from 31◦ C to 34.31◦ C and from 30.69◦ C to 33.45◦ C respectively in case of eyelid closure and opening. The corneal temperature is increased to 3.31◦ C and 2.76◦ C when blood temperature is increased from 35◦ C to 39◦ C. The corneal temperature differences are 0.31◦ C, 0.58◦ C and 0.86◦ C respectively between eyelid closure and opening at different blood temperatures 35◦ C, 37◦ C and 39◦ C respectively. Increase in blood temperature increases corneal temperature but the rate is higher in closed eye than opened eye. 658 K.C. Gokul, D.B. Gurung and P.R. Adhikary 3.1.4. Effects of lens thermal conductivities. The lens contains 65 percent of water. The water content of the lens decreases as age of individual increases[15]. Different water contents will produce different thermal conduc◦ tivities. In this study, the lens thermal conductivities 0.30W m−1 C −1 and −1 ◦ −1 −1 ◦ −1 0.54W m C are used as extreme values and 0.4W m C as a normal value. The numerical results are presented in figure 6. Figure 6. Temperature variation for different lens thermal conductivities It is observed that the corneal temperature increases from 32.46◦ C to 32.81◦ C and from 31.82◦ C to 32.27◦ C respectively in case of eyelid closure and open◦ ing, when the lens thermal conductivity is increased from 0.3W m−1 C −1 to −1 ◦ −1 0.54W m C . The corneal surface temperature is increased by 0.35◦ C and ◦ ◦ 0.45 C respectively, when lens thermal conductivity varies from 0.3W m−1 C −1 −1 ◦ −1 to 0.54W m C . The corneal temperature difference between eyelid closure and opening are 0.64◦ C, 0.58◦ C and 0.54◦ C at lens thermal conductivities 0.3, ◦ 0.4, 0.54W m−1 C −1 respectively. 3.2. Transient results. The transient temperature distribution cornea for t = 3600 seconds is presented in figure 7. The parameter values E = 0, 40 and 120W m−2 for cornea surface, E = 0, 96 and 192W m−2 for eyelid skin surface are used to find the temperature profile at Ta = 0, 25 and 50◦ C respectively. In eyelid opening, figure 7, the temperature rises rapidly over the first 1500 seconds, this increasing rate then slows down and steady state is achieved in around 1969, 2267 and 2767 seconds at Ta = 0, 25 and 50◦ C respectively. Similarly, the temperature rises rapidly over the first 8-9 min, this increasing rate then slows down and steady state is achieved in around 1326, 1374 and 1529 seconds at Ta = 0, 25 and 50◦ C respectively when eyelid is closed. The cornea surface temperature differences between eyelid opening and closure are 9.33, Thermal Effects of Eyelid in Human Eye Temperature Model 659 Figure 7. Transient temperature variations of cornea 3.62 and 0.9◦ C at ambient temperatures Ta = 0, 25 and 50◦ C respectively. The corneal surface temperature varies from 25.84◦ C to 37.21◦ C in case of eyelid opening and from 35.17◦ C to 36.31◦ C in case of eyelid closure respectively, when ambient temperature varies from 0◦ C to 50◦ C. 4. Discussion The steady state results show that when evaporation increases the eye temperature decreases. More heat is lost during closure of eyelid via evaporation. The skin is vascular where as cornea is avascular. Vascular skin has higher temperature than avascular cornea. Evaporative water molecules pick up high amount of heat energy from vascular skin than avascular cornea. Also the lipid layer of cornea prevents evaporation of tear film causing less heat transfer from cornea to environment. But eyelid skin doesn’t contain such layer to prevent evaporation. In addition, anterior part of eyelid contains glands of zeis and glands of moll, which are sweat glands that open into the hair of follicles and in a row near the base of the eyelashes. Hence more evaporation occurs from eyelid skin surface and causes rapid temperature drop in eyelid closure than opening. Increase in ambient temperature increases corneal temperature. Rapid increase in corneal temperature is observed in eyelid opening. The ambient temperature has direct effect on cornea in eyelid opening. At low ambient temperature, high amount of heat loss occurs from cornea to environment via convection and radiation, but the heating mechanism is only conduction from body core. Similarly at high ambient temperature, heat loss occurs from eye surface via radiation; tear evaporation and conduction process from cornea to body core. Hence at extreme hot/cold temperatures corneal temperature is increased /decreased rapidly during eyelid opening. However, in eyelid closure blood flow in 660 K.C. Gokul, D.B. Gurung and P.R. Adhikary vascular eyelid prevents cornea from high degree of heating/cooling in extreme environments. Hence, corneal surface temperature is less affected by ambient temperatures when eyelids are closed than opened. When blood temperature increases eye temperature also increases. The blood temperature on eyelid, when closed, has direct effect on anterior corneal surface temperature. The blood flow in anterior eyelid and posterior retina/choroid shielded the eyeball by blood temperature, when closed. So, anterior corneal surface temperature varies according to blood temperature in eyelid closure. But in eyelid opening, corneal surface temperature depends on posterior choroidal blood flow and ambient temperature. Conduction of heat from choroid/ retina is the only source of heating cornea. Hence increase in blood temperature increases corneal temperature but impact is less in case of eyelid opening. Increase in lens thermal conductivity increases corneal temperature and decreases posterior lens temperature in both cases. Higher thermal conductivities of the lens permits more heat transfer from the posterior region (high temperature region) to anterior region (low temperature region) of the eye, and thus causes the corneal surface temperature to increase. But there is no significant effect of lens thermal conductivities in case of eyelid closure and opening for corneal temperature. Hence the dominant parameters affecting eye surface temperature are blood temperatures, ambient temperatures and evaporation rates. In all these analysis the corneal center is noticed as the most sensitive point for thermal response. From figures 3, 4, 5 and 6, it is seen that the temperature changes are dominant in cornea and anterior part, very small in vitreous and almost not in sclera. The transient results show that the steady state temperature is achieved earlier in case of eyelid closure. In eyelid closure, the eye temperature increases rapidly for 600 seconds and then plateaus in around 1400 seconds. The variation in ambient temperatures has less effect on eye temperature when closed because the temperature starts to plateau from same time in all three cases and variation is only 1.14◦ C. In eyelid opening, temperature is increased mostly in the first 1500 seconds then plateaus in around 2300 seconds. The corneal temperature reaches in steady state earlier at low ambient temperature than high, when opened, and the variation 11.37◦ C is also high. High volume of blood flow in eyelid prevents rise/drop of corneal temperature in extreme ambient conditions and balance eye temperature faster in case of eyelid closure than opening. In this model the steady state corneal temperature is found to be 32.17◦ C and 34.43◦ C respectively in eyelid opening and closure. Ng and Ooi[15] summarized the corneal temperatures obtained from different experimental setup, which vary from 32◦ C to 36.6◦ C. Thus, our steady state results are in a good agreement with past experimental results. Mapstone[14] studied the effect of eyelid opening and closure experimentally using bolometer. He found a fall in temperature with a range of 0.6◦ C to 1.6◦ C when the eyelids are opened from closed position. Similarly, a rise in temperature with a range of 1.1◦ C to 2.2◦ C when the eyelids are closed from opened position. He obtained the results at 19 − 24◦ C variation Thermal Effects of Eyelid in Human Eye Temperature Model 661 of ambient temperatures. Our result at ambient temperature 25◦ C shows a rise in temperature of 2.36◦ C when the eyelids are closed from open position. Kessel et al.[10] found experimentally that steady state corneal temperature is achieved between 33◦ C to 35◦ C, when ambient temperature increased from 22◦ C to 28◦ C. Also they found that 20◦ C increase in ambient temperature, from 2◦ C to 22◦ C, increases corneal temperature by 3◦ C. In our model, the steady state corneal temperature is achieved 32.17◦ C at ambient temperature 25◦ C. Also, 20◦ C increase in ambient temperature increases corneal temperature by 3.38◦ C. The slight variation in temperature may be due to the consideration of different environmental temperatures and due to the consideration of various parameter values at the layers of human eye. 5. Conclusion We presented a one dimensional finite element model of human eye and simulate its steady and transient temperature distribution along pupillary axis in case of eyelid opening and closure. The blood flow in eyelid increases anterior temperature of human eye when eyelid is closed. The eyelid not only acts as heater/cooler of the human eye but also helps to maintain eye temperature constant. In sensitivity analysis, the effect of various blood temperatures, evaporation rates, ambient temperatures and lens thermal conductivities were analyzed. It is found that blood temperature, ambient temperature and evaporation rate are the dominant factors affecting corneal surface temperature strongly. The variation in the thermal conductivity of lens leads to minor changes in anterior corneal temperature. The transient analysis shows that, steady state temperature is achieved earlier and the variation in temperature is very less in eyelid closure than opening. The results obtained in this study are compared with the available experimental data and are in good agreement with them. The present model is one dimensional which underestimates the real temperature distribution in eye geometry. The iris is not included because it doesn’t lie on pupillary axis, but iris is assumed to be the heating/cooling source for anterior eye and it also prevents lens from extreme temperature conditions. The model did not include the effects of aqueous humor flow and tear flow in temperature distribution in human eye. Beside these disadvantages, we separated the retina/choroid region from sclera and the blood perfusion and metabolism effects were analyzed. We also modeled eyelid as a three layer structure and effects of those layers were analyzed. At last, we can conclude that eyelid is one of the significant eye component that affects eye temperature. References 1. S. Acharya, D.B. Gurung and V.P. Saxena, Effect of Metabolic Reactions on Thermoregulation in Human Males and Females Body, Appl. Math., 4(2013), 39-48. 662 K.C. Gokul, D.B. Gurung and P.R. Adhikary 2. B. Anderson, Ocular effects of changes in oxygen and carbon dioxide tension, Trans. American Ophthalmol. Society, 66 (1968), 423-474. 3. M. Cvetkovic, D. Poljak and A. Peratta, FETD Computation of the Temperature Distribution Induced into a Human Eye by a Pulsed Laser, Prog. Electromagnetics, 120(2011), 403-421. 4. V.M.M. Flycket, B.W. Roaymakers and J.J.W. Lagendijk, Modeling the impact of blood flow on the temperature distribution in the human eye and the orbit: fixed heat transfer coefficient versus the Pennes bioheat model versus discrete blood vessels, Phys. Med. Biol., 51(2006), 5007-5021. 5. H. Fujishima, I. Toda, M. Yamada, N. Sato and K. Tsubota, Corneal temperature in patients with dry eye evaluated by infrared radiation thermometry, British J. Ophthalmol., 80(1996), 29-32. 6. K.C. Gokul, D.B. Gurung and P.R. Adhikary, Effect of Blood Perfusion and Metabolism in Temperature Distribution in Human Eye, Adv. appl. math. Biosciences, 4(2013), 13-23. 7. K.C. Gokul, D.B. Gurung and P.R. Adhikary, FEM Approach for Transient Heat Transfer in Human Eye, Appl. math, 4(2013), 30-36. 8. C. Guyton and E. Hall, Text book of medical physiology, Elsevier, India, 2009. 9. A. Hirata, S. Watanabe, O. Fujiwara, M. Kojima, K. Sasaki and T. Shiozawa, Temperature elevation in the eye of anatomically based human head models for plane-wave exposures, Phys. Med. Biol., 52(2007), 6389-6399. 10. L. Kessel, L. Johnson, H. Aridsson and M. Larsen, The Relationship between Body and Ambient Temperature and Corneal Temperature, Invest. ophthalmol. & Vis. Science, 51(2010), 6593-6597. 11. J.J. Lagendijk, A mathematical model to calculate temperature distributions in human and rabbit eyes during hyperthermia treatment, Phys. Med. Biol., 27(1982), 1301-1311. 12. H.S. Lee, H. Lew and Y.S. Yun, Ultrasonographic Measurement of Upper Eyelid Thickness in Korean Children with Epicanthus, Korean J. Ophthalmol., 20 (2006), 79-81. 13. E. Li, G.R. Liu, V. Tan and Z.C. He, Modeling and simulation of bioheat transfer in the human eye using 3D alpha finite element method, International journal for numerical methods in biomedical engineering, 26 2010, 955-976. 14. R. Mapstone, Determinants of Corneal Temperature, British J. Ophthalmol., 52 (1968), 729-741. 15. E.Y.K. Ng and E.H. Ooi, FEM simulation of the eye structure with bio-heat analysis, Computational Methods Programs Biomed, 82 (2006), 268-276. 16. E.H. Ooi and E.Y.K. Ng, Simulation of aqueous humor hydrodynamics in human eye heat transfer, Comput. Biol. Med., 38 (2007), 252-262. 17. H.H. Pennes, Analysis of tissue and arterial blood temperatures in the resting human forearm, J. Appl. Physiol., 85 (1998), 5-34. 18. J.A. Scott, A finite element model of heat transport in the human eye, Phys. Med. Biol., 33 (1988), 227-241. 19. M. Shafahi and K. Vafai, Human Eye Response to Thermal Disturbances, J. Heat Transfer ASME, 133 (2011), Article ID 011009. 20. N. Sharon, P.Z.B. Yoseph, B. Bormusov and A. Dovrat, Simulation of heat exposure and damage to the eye lens in a neighborhood bakery, Exp. Eye Research, 87 (2008), 49-55. 21. D.H. Sliney, Physical factors in cataractogenesis: Ambient ultraviolet radiation and temperature, IOVS, 27 (1986), 781-790. 22. H. Wang and Q.H. Qin, FE approach with Green’s function as internal trial function for simulating bioheat transfer in the human eye, Arch. Mech., 62 (2010), 493-510. 23. T. Wessapan and P. Rattanadecho, Specific Absorbtion Rate and Temperature increase in Human eye subjected to Electromagnetic Fields at 900 MHZ, J. heat transfer, 134 (2012) 091101. Thermal Effects of Eyelid in Human Eye Temperature Model 663 K. C. Gokul received M.Phil in applied mathematics from Kathmandu University, Nepal. He has been teaching science and engineering at undergraduate and graduate levels at Kathmandu University since 2011. His research interests include computational mathematics and heat transfer modeling in human eye. Department of Natural Sciences(Mathematics), Kathmandu University, P.O.Box 6250, Kathmandu, Nepal. e-mail: [email protected] D. B. Gurung received Ph.D. in bio-mathematical modeling from Kathmandu University. He is currently an associate professor and graduate program coordinator in the department of natural sciences, school of science, Kathmandu University. His research interests are computational mathematics and mathematical modeling on various real life problems. Department of Natural Sciences(Mathematics), Kathmandu University, P.O.Box 6250, Kathmandu, Nepal. e-mail: db [email protected] P. R. Adhikary is a professor of mathematics at Kathamdnu University. His research interest is in applied mathematics. Department of Natural Sciences(Mathematics), Kathmandu University, P.O.Box 6250, Kathmandu, Nepal. e-mail: [email protected]