Temperature elevation in the eye of anatomically based human head models for plane-wave exposures

The influence of eyelashes on electric field distribution and absorbed power density in the cornea under millimeter-wave exposure

As millimeter wave (MMW) technology, particularly in fifth-generation (5G) devices, gains prominence, there is a crucial need for comprehensive electromagnetic (EM) models of ocular tissues to understand and characterize EM exposure conditions accurately. This study employs numerical modeling to investigate the interaction between MMW and the cornea, aiming to characterize EM field distributions and absorption within an anatomically accurate eye model while considering the influence of eyelashes. Using the finite-difference time-domain (FDTD) method, we conduct simulations of EM radiation interactions from 20.0 to 100.0 GHz with a human eye model. Moreover, we analyze the temperature distribution increase within the eye model using a thermal sensor in XFdtd, employing a scheme based on the finite difference (FD) method. Our findings reveal a nonuniform distribution of the EM field, particularly intensified in corneal regions adjacent to eyelashes and eyelids. Despite similar EM field patterns, the presence or absence of eyelashes has minimal impact on temperature differences. However, the study highlights increased radiation absorption by the eyelid's epidermis at 100.0 GHz, reducing the rise in the cornea's temperature.

The Effects of mmW and THz Radiation on Dry Eyes: A Finite-Difference Time-Domain (FDTD) Computational Simulation Using XFdtd

The importance of investigating the health effects of RF radiation on the cornea cannot be overstated. This study aimed to address this need by utilizing a mathematical simulation to examine the absorption of millimeter wave (mmW) and terahertz (THz) waves by the cornea, considering both normal and pathological conditions. The simulation incorporated variations in tear film thickness and hydration levels, as these factors play a crucial role in corneal health. To assess the impact of RF radiation on the cornea, the study calculated temperature rises, which indicate heating effects for both dry and normal eyes. XFdtd, a widely used commercial software based on the Finite-Difference Time Domain (FDTD) method, was employed to evaluate the radiation absorption and resulting temperature changes. The outcomes of this study demonstrated a crucial finding, i.e., that changes in the water ratio and thickness of the tear film, which are associated with an increased risk of dry eye syndrome, directly impact the absorption of mmW and THz waves by the cornea. This insight provides valuable evidence supporting the interconnection between tear film properties and the vulnerability of the cornea to RF radiation.

Assessment of incident power density in different shapes of averaging area for radio-frequency exposure above 6 GHz

Objective. The International Commission on Non-Ionizing Radiation Protection guidelines and IEEE C95.1-2019 standard for human protection from local electromagnetic field exposure above 6 GHz state that absorbed (or epithelial) power density (APD) and incident power density (IPD), averaged over a square area, are internal and external physical quantities, respectively, that set the exposure limit. Per exposure standards, the measurement procedure and evaluation of the IPD have been established in technical standards, where a circular averaging area is recommended only for non-planar surfaces in IEC/IEEE 63195-1 and -2. In this study, the effects of two averaging shapes on the APD and IPD are evaluated computationally to provide new insights from the viewpoint of exposure standards. Approach. The relation between the APD, IPD, and the steady-state temperature rise (heating factor) in rectangular and human models for exposure to a single dipole, dipole arrays, and the Gaussian beams is investigated computationally with finite-difference method. Main results. The maximum differences in the heating factor of the APD and IPD for square and circular averaging areas were 4.1% and 4.4% for the antenna–model distance >5 mm, respectively. These differences appear when the beam pattern on the model surface has an elliptical shape. For an antenna–model distance ≤5 mm and at frequencies ≤15 GHz, the heating factors for square averaging areas were not always conservative to those for circular ones (−7.8% for IPD), where only the antenna feed point are visible before beam formation. Significance. The heating factors of the APD and IPD for a circular averaging area are conservative for near-field exposure of canonical sources for frequencies up to 300 GHz, except for a beam with a significant major-to-minor axis ratio and an angle of 30°–60° to a square averaging area. This tendency would help bridge the gap between exposure and product standards.

Characterizing the Aging Process of the Human Eye: Tear Evaporation, Fluid Dynamics, Blood Flow, and Metabolism-Based Comparative Study

Eye temperature and intraocular pressure are two measurable parameters that can be monitored as a health index with aging. Deviations from the normal range of intraocular pressure and temperature lead to the formation of many diseases. This study has been carried out to evaluate the relations between the physiological and anatomical changes of the eye with aging using mathematical modeling. 2D computer-aided design of the human eye has been developed for two major groups: 21 to 30 years and 41 to 50 years. The computer simulation has been carried out to determine the effects of physiological changes of tear evaporation, fluid dynamics, blood flow, and metabolism of eye tissues with aging. The simulation has been carried out in the standing and the supine position of a human body. The rate of temperature change is – 0.0075 K per year in the standing position and – 0.007 K per year in the supine position because of the modeled anatomical and physiological effects. All the three simulation parameters of this study, the temperature of the human eye, the intraocular pressure, and the aqueous humor flow velocity, have been compared with the recent practical and simulation-based experiments to validate our results.

Assessment of absorbed power density and temperature rise for nonplanar body model under electromagnetic exposure above 6 GHz

The averaged absorbed power density (APD) and temperature rise in body models with nonplanar surfaces were computed for electromagnetic exposure above 6 GHz. Different calculation schemes for the averaged APD were investigated. Additionally, a novel compensation method for correcting the heat convection rate on the air/skin interface in voxel human models was proposed and validated. The compensation method can be easily incorporated into bioheat calculations and does not require information regarding the normal direction of the boundary voxels, in contrast to a previously proposed method. The APD and temperature rise were evaluated using models of a two-dimensional cylinder and a three-dimensional partial forearm. The heating factor, which was defined as the ratio of the temperature rise to the APD, was calculated using different APD averaging schemes. Our computational results revealed different frequency and curvature dependences. For body models with curvature radii of >30 mm and at frequencies of >20 GHz, the differences in the heating factors among the APD schemes were small.

Model-based approach for analyzing prevalence of nuclear cataracts in elderly residents

Recent epidemiological studies have hypothesized that the prevalence of cortical cataracts is closely related to ultraviolet radiation. However, the prevalence of nuclear cataracts is higher in elderly people in tropical areas than in temperate areas. The dominant factors inducing nuclear cataracts have been widely debated. In this study, the temperature increase in the lens due to exposure to ambient conditions was computationally quantified in subjects of 50-60 years of age in tropical and temperate areas, accounting for differences in thermoregulation. A thermoregulatory response model was extended to consider elderly people in tropical areas. The time course of lens temperature for different weather conditions in five cities in Asia was computed. The temperature was higher around the mid and posterior part of the lens, which coincides with the position of the nuclear cataract. The duration of higher temperatures in the lens varied, although the daily maximum temperatures were comparable. A strong correlation (adjusted R² > 0.85) was observed between the prevalence of nuclear cataract and the computed cumulative thermal dose in the lens. We propose the use of a cumulative thermal dose to assess the prevalence of nuclear cataracts. Cumulative wet-bulb globe temperature, a new metric computed from weather data, would be useful for practical assessment in different cities.

Comparison of Thermal Response for RF Exposure in Human and Rat Models

In the international guidelines/standards for human protection against electromagnetic fields, the specific absorption rate (SAR) is used as a metric for radio-frequency field exposure. For radio-frequency near-field exposure, the peak value of the SAR averaged over 10 g of tissue is treated as a surrogate of the local temperature elevation for frequencies up to 3⁻10 GHz. The limit of 10-g SAR is derived by extrapolating the thermal damage in animal experiments. However, no reports discussed the difference between the time constant of temperature elevation in small animals and humans for local exposure. This study computationally estimated the thermal time constants of temperature elevation in human head and rat models exposed to dipole antennas at 3⁻10 GHz. The peak temperature elevation in the human brain was lower than that in the rat model, mainly because of difference in depth from the scalp. Consequently, the thermal time constant of the rat brain was smaller than that of the human brain. Additionally, the thermal time constant in human skin decreased with increasing frequency, which was mainly characterized by the effective SAR volume, whereas it was almost frequency-independent in the human brain. These findings should be helpful for extrapolating animal studies to humans.

Temperature elevation in the human brain and skin with thermoregulation during exposure to RF energy

Background

Two international guidelines/standards for human protection from electromagnetic fields define the specific absorption rate (SAR) averaged over 10 g of tissue as a metric for protection against localized radio frequency field exposure due to portable devices operating below 3–10 GHz. Temperature elevation is suggested to be a dominant effect for exposure at frequencies higher than 100 kHz. No previous studies have evaluated temperature elevation in the human head for local exposure considering thermoregulation. This study aims to discuss the temperature elevation in a human head model considering vasodilation, to discuss the conservativeness of the current limit.

Methods

This study computes the temperature elevations in an anatomical human head model exposed to radiation from a dipole antenna and truncated plane waves at 300 MHz–10GHz. The SARs in the human model are first computed using a finite-difference time-domain method. The temperature elevation is calculated by solving the bioheat transfer equation by considering the thermoregulation that simulates the vasodilation.

Results

The maximum temperature elevation in the brain appeared around its periphery. At exposures with higher intensity, the temperature elevation became larger and reached around 40 °C at the peak SAR of 100 W/kg, and became lower at higher frequencies. The temperature elevation in the brain at the current limit of 10 W/kg is at most 0.93 °C. The effect of vasodilation became notable for tissue temperature elevations higher than 1–2 °C and for an SAR of 10 W/kg. The temperature at the periphery was below the basal brain temperature (37 °C).

Conclusions

The temperature elevation under the current guideline for occupational exposure is within the ranges of brain temperature variability for environmental changes in daily life. The effect of vasodilation is significant, especially at higher frequencies where skin temperature elevation is dominant.

Monte Carlo simulations of skin exposure to electromagnetic field from 10 GHz to 1 THz

In this study, we present an assessment of human-body exposure to an electromagnetic field at frequencies ranging from 10 GHz to 1 THz. The energy absorption and temperature elevation were assessed by solving boundary value problems of the one-dimensional Maxwell equations and a bioheat equation for a multilayer plane model. Dielectric properties were measured [Formula: see text] at frequencies of up to 1 THz at body temperature. A Monte Carlo simulation was conducted to assess variations of the transmittance into a skin surface and temperature elevation inside a body by considering the variation of the tissue thickness due to individual differences among human bodies. Furthermore, the impact of the dielectric properties of adipose tissue on temperature elevation, for which large discrepancies between our present measurement results and those in past works were observed, was also examined. We found that the dielectric properties of adipose tissue do not impact on temperature elevation at frequencies over 30 GHz. The potential risk of skin burn was discussed on the basis of the temperature elevation in millimeter-wave and terahertz-wave exposure. Furthermore, the consistency of the basic restrictions in the international guidelines set by ICNIRP was discussed.

Human exposure to pulsed fields in the frequency range from 6 to 100 GHz

Restrictions on human exposure to electromagnetic waves at frequencies higher than 3-10 GHz are defined in terms of the incident power density to prevent excessive temperature rise in superficial tissue. However, international standards and guidelines differ in their definitions of how the power density is interpreted for brief exposures. This study investigated how the temperature rise was affected by exposure duration at frequencies higher than 6 GHz. Far-field exposure of the human face to pulses shorter than 10 s at frequencies from 6 to 100 GHz was modelled using the finite-difference time-domain method. The bioheat transfer equation was used for thermal modelling. We investigated the effects of frequency, polarization, exposure duration, and depth below the skin surface on the temperature rise. The results indicated limitations in the current human exposure guidelines and showed that radiant exposure, i.e. energy absorption per unit area, can be used to limit temperature rise for pulsed exposure. The data are useful for the development of human exposure guidelines at frequencies higher than 6 GHz.

THE EFFECT OF GAZE ANGLE ON THE EVALUATIONS OF SAR AND TEMPERATURE RISE IN HUMAN EYE UNDER PLANE-WAVE EXPOSURES FROM 0.9 TO 10 GHZ

This article investigates the effect of gaze angle on the specific absorption rate (SAR) and temperature rise in human eye under electromagnetic exposures from 0.9 to 10 GHz. Eye models in different gaze angles are developed based on biometric data. The spatial-average SARs in eyes are investigated using the finite-difference time-domain method, and the corresponding maximum temperature rises in lens are calculated by the finite-difference method. It is found that the changes in the gaze angle produce a maximum variation of 35, 12 and 20 % in the eye-averaged SAR, peak 10 g average SAR and temperature rise, respectively. Results also reveal that the eye-averaged SAR is more sensitive to the changes in the gaze angle than peak 10 g average SAR, especially at higher frequencies.

Mathematical model of thermal effects of blinking in human eye

Blinking is regarded as the continuous interrupted eyelid closure or opening and its thermal effect will compromise between these two. During a blink, the heat loss via convection, radiation and tear evaporation from cornea is prevented, warm tear is layered across corneal surface and the vessels of the palpebral conjunctiva provide heat to anterior eye. In most of the thermal models in human eye that are found in literatures, effect of blinking is not included, simulation is carried out only in open eye. Thus, in this paper, thermal effects of blinking are investigated using one-dimensional finite element method in transient state case. The bio-heat transfer process is simulated during different blinking rates, lid closure and opening. The simulation is carried out using normal and extreme values of ambient temperatures, blood temperatures, evaporation rates, blood perfusion rates, and lens thermal conductivities. Blinking is found to increase corneal and lens temperature by 1.29°C and 0.78°C respectively when compared to open eye. The results obtained from this model are useful in predicting temperature distribution in different laser eye surgeries, hyperthermia and cryosurgery treatment of eyelid carcinoma, choroidal melanoma and can be used for diagnosing temperature-related diseases.

Dielectric property measurement of ocular tissues up to 110 GHz using 1 mm coaxial sensor

Measurement of the dielectric properties of ocular tissues up to 110 GHz was performed by the coaxial probe method. A coaxial sensor was fabricated to allow the measurement of small amounts of biological tissues. Four-standard calibration was applied in the dielectric property measurement to obtain more accurate data than that obtained with conventional three-standard calibration, especially at high frequencies. Novel data of the dielectric properties of several ocular tissues are presented and compared with data from the de facto database.

Generation of infant anatomical models for evaluating electromagnetic field exposures

Realistic anatomical modeling is essential in analyzing human exposure to electromagnetic fields. Infants have significant physical and anatomical differences compared with other age groups. However, few realistic infant models are available. In this work, we developed one 12-month-old male whole body model and one 17-month-old male head model from magnetic resonance images. The whole body and head models contained 28 and 30 tissues, respectively, at spatial resolution of 1 mm × 1 mm × 1 mm. Fewer identified tissues in the whole body model were a result of the low original image quality induced by the fast imaging sequence. The anatomical and physical parameters of the models were validated against findings in published literature (e.g., a maximum deviation as 18% in tissue mass was observed compared with the data from International Commission on Radiological Protection). Several typical exposure scenarios were realized for numerical simulation. Dosimetric comparison with various adult and child anatomical models was conducted. Significant differences in the physical and anatomical features between adult and child models demonstrated the importance of creating realistic infant models. Current safety guidelines for infant exposure to radiofrequency electromagnetic fields may not be conservative.

Numerical analysis of specific absorption rate in the human head due to a 13.56 MHz RFID-based intra-ocular pressure measurement system

A modern wireless intra-ocular pressure monitoring system, based on 13.56 MHz inductively coupled data transmission, was dosimetrically analyzed with respect to the specific absorption rate (SAR) induced inside the head and the eye due to the electromagnetic field exposure caused by the reader antenna of the transmission system. The analysis was based on numerical finite difference time domain computations using a high resolution anatomical eye model integrated in a modern commercially available anatomical model of a male head. Three different reader antenna configurations, a 7-turn elliptic (30 mm × 50 mm) antenna at 12 mm distance from the eye, a flexible circular antenna (60 mm diameter, 8 turns on 2 mm substrate) directly attached to the skin, and a circular 7-turn antenna (30 mm diameter at 12 mm distance to the eye) were analyzed, respectively. Possible influences of the eye-lid status (closed or opened) and the transponder antenna contained in a contact lens directly attached to the eye were taken into account. The results clearly demonstrated that for typical reader antenna currents required for proper data transmission, the SAR values remain far below the limits for localized exposure of the head, as defined by the International Commission for Non-Ionizing Radiation Protection. Particularly the induced SAR inside the eye was found to be substantially (orders of magnitudes for typical reader antenna currents in the order of 1 A turn) below values which have been reported to be critical with respect to thermally induced adverse health effects in eye tissues.

FDTD analysis of temperature elevation in the lens of human and rabbit models due to near-field and far-field exposures at 2.45 GHz

The eye is said to be one of the most sensitive organs to microwave heating. According to previous studies, the possibility of microwave-induced cataract formation has been experimentally investigated in rabbit and monkey eyes, but not for the human eye due to ethical reasons. In the present study, the temperature elevation in the lens, the skin around the eye and the core temperature of numerical human and rabbit models for far-field and near-field exposures at 2.45 GHz are investigated. The temperature elevations in the human and rabbit models were compared with the threshold temperatures for inducing cataracts, thermal pain in the skin and reversible health effects such as heat exhaustion or heat stroke. For plane-wave exposure, the core temperature elevation is shown to be essential both in the human and in the rabbit models as suggested in the international guidelines and standards. For localised exposure of the human eye, the temperature elevation of the skin was essential, and the lens temperature did not reach its threshold for thermal pain. On the other hand, the lens temperature elevation was found to be dominant for the rabbit eye.

Dominant factors affecting temperature rise in simulations of human thermoregulation during RF exposure

Numerical models of the human thermoregulatory system can be used together with realistic voxel models of the human anatomy to simulate the body temperature increases caused by the power absorption from radio-frequency electromagnetic fields. In this paper, the Pennes bioheat equation with a thermoregulatory model is used for calculating local peak temperatures as well as the body-core-temperature elevation in a realistic human body model for grounded plane-wave exposures at frequencies 39, 800 and 2400 MHz. The electromagnetic power loss is solved by the finite-difference time-domain (FDTD) method, and the discretized bioheat equation is solved by the geometric multigrid method. Human thermoregulatory models contain numerous thermophysiological and computational parameters--some of which may be subject to considerable uncertainty--that affect the simulated core and local temperature elevations. The goal of this paper is to find how greatly the computed temperature is influenced by changes in various modelling parameters, such as the skin blood flow rate, models for vasodilation and sweating, and clothing and air movement. The results show that the peak temperature rises are most strongly affected by the modelling of tissue blood flow and its temperature dependence, and mostly unaffected by the central control mechanism for vasodilation and sweating. Almost the opposite is true for the body-core-temperature rise, which is however typically greatly lower than the peak temperature rise. It also seems that ignoring the thermoregulation and the blood temperature increase is a good approximation when the local 10 g averaged specific absorption rate is smaller than 10 W kg(-1).

Children and adults exposed to electromagnetic fields at the ICNIRP reference levels: theoretical assessment of the induced peak temperature increase

To avoid potentially adverse health effects of electromagnetic fields (EMF), the International Commission on Non-Ionizing Radiation Protection (ICNIRP) has defined EMF reference levels. Restrictions on induced whole-body-averaged specific absorption rate (SAR(wb)) are provided to keep the whole-body temperature increase (T(body, incr)) under 1 °C during 30 min. Additional restrictions on the peak 10 g spatial-averaged SAR (SAR(10g)) are provided to prevent excessive localized tissue heating. The objective of this study is to assess the localized peak temperature increase (T(incr, max)) in children upon exposure at the reference levels. Finite-difference time-domain modeling was used to calculate T(incr, max) in six children and two adults exposed to orthogonal plane-wave configurations. We performed a sensitivity study and Monte Carlo analysis to assess the uncertainty of the results. Considering the uncertainties in the model parameters, we found that a peak temperature increase as high as 1 °C can occur for worst-case scenarios at the ICNIRP reference levels. Since the guidelines are deduced from temperature increase, we used T(incr, max) as being a better metric to prevent excessive localized tissue heating instead of localized peak SAR. However, we note that the exposure time should also be considered in future guidelines. Hence, we advise defining limits on T(incr, max) for specified durations of exposure.

Computation of temperature elevation in rabbit eye irradiated by 2.45-GHz microwaves with different field configurations

This study calculated the temperature elevation in the rabbit eye caused by 2.45-GHz near-field exposure systems. First, we calculated specific absorption rate distributions in the eye for different antennas and then compared them with those observed in previous studies. Next, we re-examined the temperature elevation in the rabbit eye due to a horizontally-polarized dipole antenna with a C-shaped director, which was used in a previous study. For our computational results, we found that decisive factors of the SAR distribution in the rabbit eye were the polarization of the electromagnetic wave and antenna aperture. Next, we quantified the eye average specific absorption rate as 67 W kg(-1) for the dipole antenna with an input power density at the eye surface of 150 mW cm(-2), which was specified in the previous work as the minimum cataractogenic power density. The effect of administrating anesthesia on the temperature elevation was 30% or so in the above case. Additionally, the position where maximum temperature in the lens appears is discussed due to different 2.45-GHz microwave systems. That position was found to appear around the posterior of the lens regardless of the exposure condition, which indicates that the original temperature distribution in the eye was the dominant factor.