FEM Approach for Transient Heat Transfer in Human Eye

In this paper, a bio-heat transfer model of temperature distribution in human eye is discussed using appropriate boundary conditions for cornea and sclera. Variational finite element method with Crank-Nicolson scheme is used to calculate the transient temperature distribution in normal human eye. The temperature with and without the effect of blood perfusion and metabolism on retina is simulated and compared for various ambient temperatures, evaporation rates and lens thermal conductivities. The obtained results are compared with experimental results and past results found in literatures. The results show that the steady state corneal temperature is achieved in around 31 and 45 minute of exposure at ambient temperatures 10 ̊C and 50 ̊C respectively. Steady state eye temperature is achieved earlier at higher evaporation rate. Similar result is achieved for higher lens thermal conductivity and also for lower ambient temperature.


Introduction
In a human body the internal body temperature almost remains constant despite the fluctuation of environmental temperatures up to certain limits.The main organ that keeps core temperature constant is dermal part [1].There is no skin layer to keep the core temperature constant in case of human eye.The skin layer (eyelid) covers eye surface (cornea) for 3 seconds in a minute (in an average).But for 57 seconds in a minute,ocular surface (cornea) has to manage thermal stress of an environment.The human eye is relatively a small and complex organ, consists of several sub domains with different material properties and having complex geometry.
The calculation of the temperature distribution in human eye when it is heated or cooled is an important aspect of the development of infrared and radiofrequency safety guidelines and for hyperthermia and thermo-therapy treatments of various ocular diseases [2].The severity of the physiological effect produced by small temperature increases can cause eyesight to worsen.Actually, a small temperature increase in the eye of 3˚C -5˚C leads to induce cataracts formation [3].Some researchers believe that thermal effect can induce cataracts; other believes that it is the result of other biological and genetic issues.One of the early theory suggested that heat exchanges within the anterior eye caused the cataract.
Verhoeff and Bell argued that cataract formed on the posterior surface of the lens because the anterior surface was cooled by circulation of the aqueous humor and so the cornea was air cooled [4].Investigation in Germany showed that the cataracts were due to the raised temperature induced indirectly through heat absorbed by the iris, where a rich blood supply would be consistent with a high degree of heating.At the same time Salil noticed the rise in cataracts one year after a very hot, dry summer in Iowa, highlighting likely environmental causes of cataracts [4].
Due to convective heat transport of the blood vessels; the blood picks up energy from hot areas and deposits this at cooler area or vice versa.The difficulty of modeling in the eye is due to the impact of blood flow on the heat transfer; however, incorporating the impact of blood flow in heat transport calculation is very important [2].The temperature inside the human body depends on the degree of temperature, duration of exposure and the environmental conditions which cause heat gain/loss from tissues [1].Hence, blood flow and time are the main factors that affect temperature distribution in human eye.
Lagendijk [5] 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 [6] utilized finite element method to obtain the temperature profile based on heat conduction using various heat transfer coefficients given by Lagendijk.Amara [7] presented a numerical thermal model of laser-ocular media interaction.Ng and Ooi [8] studied the effect of aqueous humor hydrodynamics on heat transfer within human eye by neglecting blood perfusion and metabolism.They neglected the effects because they assumed that perfusion in iris/cilliary is sufficient to maintain their temperature at 37˚C.Li et al. [9] studied the bio-heat transfer in the human eye, neglecting the effect of perfusion and metabolism, using 3D alpha finite element method.They assumed the contribution due to perfusion and metabolism is very small because they occur only on small part in the human eye.Cvetkovic et al. [10] developed a thermal model and studied effects of pulsed laser in human eye.Narasimhan et al. [11] developed transient model to study heat transfer in human eye undergoing laser surgery.Flyckt et al. [2] studied the impact of choroidal blood flow and scleral convection heat transfer coefficient in human eye.
All the previously developed models have neglected the effects of blood perfusion and metabolism on retina/iris/cilliary body.The significance of blood perfusion and metabolism on the temperature distribution in the eye is debatable, since they take place only on retina, choroid, iris and cilliary body, which constitute very small part of human eye.The blood flow in the iris/sclera part plays significant role to adjust eye temperature with the rest of the body [3].In [8] it is clearly mentioned that "Due to lack of literature data, the perfusion term is neglected.The effect of this assumption on the accuracy of the model however remaining unknown".The retina has perhaps the highest oxygen consumption rate (metabolism) of any tissue in the body [12].Hence it is necessary to investigate the effects of blood perfusion and metabolism in order to obtain a more accurate result.The objective of this paper is to determine the transient temperature distribution in human eye.Hence, this model provides the results for transient temperature distribution including the effect of blood perfusion and metabolism.

Model Formulation
The eye is assumed a perfectly bonded solid structure with each component homogeneous.The eye is considered having six major components: cornea, aqueous humor, lens, vitreous humor, retina, and sclera.As sketched on Figure 1, the thickness of cornea, aqueous humor, lens, vitreous humor, retina and sclera have been considered as 1   A one dimensional finite element human eye model has been developed to simulate its thermal unsteady state conditions.The governing differential equation used for heat flow in the human eye due to Pennes [13] is given by: where, ρ = tissue density (kg/m 3 ), c = tissue specific heat (J/kg˚C), ρ b = blood density (kg/m 3 ), c b = blood specific heat (J/kg˚C), k = tissue thermal conductivity (W/m˚C), ω = volumetric blood perfusion rate per unit volume (s −1 ), T b = blood temperature (˚C), T = tissue temperature (˚C), t = time (s), Q m = heat generation due to metabolism (W/m 3 ).
The three terms on the right-hand side of bio-heat Equation ( 1) are blood conduction, blood perfusion, and metabolism.In this model, the effect of blood perfusion and metabolism is analyzed only on retina.
Boundary conditions for the system can be defined as follows [9,14]: 1) In the back of the eye, heat is transferred from blood via ophthalmic vessels to the sclera: where  is the normal direction to the surface boundary, k s is the thermal conductivity of sclera, h b is the heat transfer coefficient between blood and eye (W/m 2 ˚C), and T b is blood temperature (˚C).
2) At the cornea, heat loss from the eye occurs through convection, radiation, and tear evaporation: where h  represents the convection heat transfer coeffi-cient between the cornea and ambient environment (W/m 2 ˚C), T  is the ambient room temperature (˚C), σ is the Stefan Boltzmann constant (5.67 × 10 −8 W/m 2 ˚C4 ), ε is the emissivity of the cornea, and 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:   where, T represents an initial guess of temperature.The iteration is completed when the convergent condition is satisfied: where  is iteration tolerance.
The inner body core temperature T b is assumed to be 37˚C.Therefore, the initial boundary condition is The partial differential Equation (1) together with boundary conditions (2) and ( 5) in one dimensional variational form is: We write I separately for the six layers in the following equation: Further, to optimize I, we differentiate I partially with respect to T i and equating to zero, we get, 0, 0,1, 2, , 5 The system of Equations ( 11) can be written in matrix form where and [K] are 6 6   matrices called capacity and conductivity matrices respectively.Now we apply Crank-Nicolson method to solve the system (12) with respect to time, we get the following relation where t  is time interval.The temperature increases from outer surface of cornea towards eye core when ambient temperature is less than 37˚C and vice versa.Hence, we consider the temperature increases/decreases in linear order towards body core with regard to thickness.For initial nodal temperatures   0 T at time 0 t  , we assume the following initial condition where   0, 0 T  20˚C and r  constant to be deter- mined.The Equation ( 13) is repeatedly solved to get the required nodal temperatures.

Results and Discussion
To solve equation ( 13), the following values of parameters are considered [14]: .The numerical calculation for unsteady state temperature distribution is carried out for different parts of human eye with and without taking blood perfusion and metabolism on the retina.The effects of different ambient temperatures, tear evaporation rates and lens thermal conductivities are studied and compared.The results obtained are compared with experimental and other results found in literatures.The overall thermal behavior of human eye is observed for 3600 seconds using one second time step size.

Effect of Blood Perfusion and Metabolism
The transient temperature distribution of several parts of human eye for It can be seen from Figures 2 and 3 that the anterior and posterior part of cornea, aqueous and lens temperatures begin to stabilize around 1870 seconds (approxi-mately 31.17minutes) and 2740 seconds(approximately 45.67 minutes) respectively.In human eye, heat gain occurs through conduction, perfusion, metabolism, blinking, tear flow, evaporation, and convection but heat loss occurs only through conduction, evaporation, convection and radiation.More factors are involved in heating eye components than cooling.Hence, eye is more vulnerable when it is exposed to high temperatures (high ambient temperatures, hyperthermia treatment, laser surgery etc.) than low (low ambient temperatures, cryosurgery treatment etc).The temperature difference obtained between anterior and posterior parts of cornea are 2.2˚C and 0.65˚C respectively at 10 C T    and T ∞ = 50˚C.As well the temperature differences between anterior and posterior parts of lens are 0.2˚C and 0.05˚C.
The corneal temperature distribution with and without blood perfusion and metabolism in retina at T ∞ = 10˚C and T ∞ = 50˚C are shown in Figures 4 and 5 below.
It can be observed from Figures 4 and 5 that when corneal surface temperature is attained steady state in both cases, the temperature differences with and without blood perfusion and metabolism obtained are 0.23˚C and   0.06˚C respectively.Steady state corneal surface temperature is reached earlier in case of having blood perfusion and metabolism in retina.Accordingly, steady state corneal surface temperature is reached earlier at T ∞ = 10˚C compared to that at T ∞ = 50˚C.This is due to better heating mechanism than cooling in eye.
For further investigation, we now put the effect with blood perfusion and metabolism as Case I and without them as Case II.

Effect of Tear Evaporation
The corneal surface contains a three-layered structure-a mucoid layer, a thick aqueous layer and thin oily layer.The function of oily layer is to prevent evaporation of tear from the corneal surface.When the oily layer is destroyed, the evaporation rate increases dramatically [15].For different evaporation rates of 20 W/m 2 , 40 W/m 2 , 100 W/m 2 , 180 W/m 2 , the corneal surface temperature differences between Case I and Case II are obtained as 0.10˚C, 0.12˚C, 0.18˚C, and 0.27˚C respectively.The differences show that when evaporation rate increases the nodal temperature decreases in both cases but decreasing rate is slow in Case I than Case II.It means, steady state corneal surface temperature is reached earlier in Case I than Case II.This is due to the effect of perfusion and metabolism on retina.The corneal surface temperature is dropped in by 4.78˚C and 4.95˚C in Cases I and II respectively, when evaporation rate is reached from 20 W/m 2 to 180 W/m 2 .Also, steady state corneal temperature is achieved earlier at E = 180 W/m 2 than at 20 W/m 2 .The sudden decrease in corneal temperature at E = 180 W/m 2 increases the temperature difference between cornea and aqueous.It is because the greater the temperature difference between two surfaces, the faster is the rate of transfer of thermal energy.Hence steady state corneal temperature is reached earlier at higher evaporation rate.

Effect of Ambient Temperatures
Heat losses occur due to convection, radiation and tear evaporation at cornea.This loss is strongly related to ambient temperature.In addition, ambient temperature is one of the factors affecting the amount of tear in the eyes [9].Four sets of ambient temperatures 10˚C, 25˚C, 40˚C, and 50˚C are considered for analysis.The numerical results are presented in Figures 8 and 9.
It can be seen from Figures 8 and 9 that, increase in ambient temperature from 10˚C to 50˚C increases corneal temperature from 27.69˚C to 39.69˚C in Case I and from 27.41˚C to 39.78˚C in Case II.It shows that increase in ambient temperature increases eye temperatures lower in case I compared to Case II and vice versa.This is due to the cooling effect of perfusion at retina.Also, the steady state corneal temperature is reached earlier at T ∞ = 10˚C than at T ∞ = 50˚C.This is because; more factors are involved in heating eye components than cooling.

Effect of Lens Thermal Conductivities
It is well known that the water content of lens decreases as age increases.Decrease in water level of lens increases its hardness.This process changes the thermal conductivity of lens due to age [14].Four sets of values  It is found that, increase in lens thermal conductivity from 0.21 W/m ˚C to 0.54 W/m ˚C increases corneal temperature by 0.87˚C in Case I and by 0.85˚C in Case II but decreases posterior lens temperature by 0.23˚C in Case I and by 0.25˚C in Case II.More heat transfer occurs via conduction from the posterior region to anterior region, when thermal conductivity of lens increases.As a result corneal surface temperature is increased.The steady state corneal temperature is achieved earlier at K 3 = 0.54 W/m ˚C than at 0.21 W/m ˚C but the process is reversed in case of posterior part of lens.Similarly, steady state eye temperature is reached earlier in Case I compared to Case II due to the effect of perfusion and metabolism on retina.

Conclusions
In this model, a comparative study of temperature distribution with and without considering the effect of blood perfusion and metabolism on retinal part of human eye is  presented.The steady state corneal temperature is found to be 32.17˚C in Case I and 32.05˚C in Case II respectively.Many authors [1,3,10,13,15] reported similar temperature distribution along the eye pupillary axis.The reported values range from 30.92˚C to 33.7˚C on corneal surface.Hence, our steady state results are in accordance with past results.In earlier studies [1,3,10,13,15], the effect of blood perfusion and metabolism is assumed negligible on retinal part.We are able to show that blood perfusion and metabolism play an important role in maintaining the eye temperature.
The surface temperature of the cornea was measured using bolometer by Mepastone [16], the mean temperature variation on the cornea surface was obtained to be 0.8˚C under 33.2˚C -36˚C.Kessel et al. [17] found that steady state corneal temperature was achieved between 33˚C to 35˚C, when ambient temperature was increased.Also they found that 20˚C increase in ambient temperature, from 2˚C to 22˚C, was required to increase corneal temperature by 3˚C.In our model, the steady state corneal temperature is achieved between 32.17˚C to 35.58˚C, when ambient temperature is increased from 20˚C to 40˚C.Also, a 20˚C increase in ambient temperature is required to increase corneal temperature by 3.38˚C.Hence, the temperature obtained from our model agrees with the experimental results obtained by Mepastone [16] and Kessel et al. [17].The slight differences whatsoever obtained may be due to consideration of various parameter values at the layers of human eye.
This model may help to understand the thermal behavior of eye, which can be very crucial in ocular diseases like corneal pain, presbyopia, cataracts etc.This may also prove to be valuable in several eye therapies and surgeries like hyperthermia using microwave, heating due to laser etc.Thus, the model could be useful for the researchers to study the effects of heat flux inside the human eye and medical scientists to improve the diagnosis and treatment.
. Similarly, 0 1 2 3 4 5 , , , , , T T T T T T and 6 b T T  (body core temperature) are the nodal temperatures at a dis-

Figure 1 .
Figure 1.Finite element sketch of human eye.
Four sets of data values for E equals 20 W/m 2 , 40 W/m 2 , 100 W/m 2 , 180 W/m 2 are used in this investigation.The temperature variations for different evaporation rates at 25 C T    , 3 0.4 W/m C K   are shown in Figures 6 and 7.

Figure 6 .
Figure 6.Corneal surface temperature for different evaporation rates (Case I).

Figure 7 .
Figure 7. Corneal surface temperature for different evaporation rates (Case II).

Figure 8 .
Figure 8. Corneal surface temperature for different ambient temperatures (Case I).

Figure 9 .
Figure 9. Corneal surface temperature for different ambient temperatures (Case II).

Figure 10 .
Figure 10.Anterior corneal and posterior lens temperature for different Lens thermal conductivities (Case I).

Figure 11 .
Figure 11.Anterior corneal and posterior lens temperature for different Lens thermal conductivities (Case II).