Power deposition in the head and neck of an anatomically based human body model for plane wave exposures

Electrical conductivity and permittivity maps of brain tissues derived from water content based on T1 -weighted acquisition

PURPOSE

To develop an electrical properties tomography (EPT) technique that can provide in vivo electrical conductivity and permittivity images of biological tissue without performing complex-valued radiofrequency field measurements.

THEORY AND METHODS

Electrical conductivity and permittivity images are modeled as a monotonic function of tissues' water content (W) under the principle of Maxwell's mixture theory. Water content maps are estimated from two spin-echo images having different repetition times (TRs). For the modeling functions, physically measured parameters (electrical properties, water content, and T₁ ) of brain cerebrospinal fluid (CSF), gray matter, and white matter are used as landmark literature references. The formulations are validated by a developed electrolyte-protein phantom and by human brain studies at 3 Tesla (T).

RESULTS

The electrical properties (EPs) of the phantom estimated by the proposed method match well with the values measured on the bench. The conductivity and permittivity maps from all experiments show uncompromised spatial resolution without boundary artifacts and higher contrast when compared with water content maps.

CONCLUSIONS

Human brain and phantom EP images suggest that water content is a dominating factor in determining the electrical properties of tissues. Despite possible literature inaccuracies, the proposed method offers EP maps that can provide complementary information to current approaches, to facilitate EPT scans in clinical applications. Magn Reson Med 77:1094-1103, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

An exponential growth of computational phantom research in radiation protection, imaging, and radiotherapy: a review of the fifty-year history

Radiation dose calculation using models of the human anatomy has been a subject of great interest to radiation protection, medical imaging, and radiotherapy. However, early pioneers of this field did not foresee the exponential growth of research activity as observed today. This review article walks the reader through the history of the research and development in this field of study which started some 50 years ago. This review identifies a clear progression of computational phantom complexity which can be denoted by three distinct generations. The first generation of stylized phantoms, representing a grouping of less than dozen models, was initially developed in the 1960s at Oak Ridge National Laboratory to calculate internal doses from nuclear medicine procedures. Despite their anatomical simplicity, these computational phantoms were the best tools available at the time for internal/external dosimetry, image evaluation, and treatment dose evaluations. A second generation of a large number of voxelized phantoms arose rapidly in the late 1980s as a result of the increased availability of tomographic medical imaging and computers. Surprisingly, the last decade saw the emergence of the third generation of phantoms which are based on advanced geometries called boundary representation (BREP) in the form of Non-Uniform Rational B-Splines (NURBS) or polygonal meshes. This new class of phantoms now consists of over 287 models including those used for non-ionizing radiation applications. This review article aims to provide the reader with a general understanding of how the field of computational phantoms came about and the technical challenges it faced at different times. This goal is achieved by defining basic geometry modeling techniques and by analyzing selected phantoms in terms of geometrical features and dosimetric problems to be solved. The rich historical information is summarized in four tables that are aided by highlights in the text on how some of the most well-known phantoms were developed and used in practice. Some of the information covered in this review has not been previously reported, for example, the CAM and CAF phantoms developed in 1970s for space radiation applications. The author also clarifies confusion about 'population-average' prospective dosimetry needed for radiological protection under the current ICRP radiation protection system and 'individualized' retrospective dosimetry often performed for medical physics studies. To illustrate the impact of computational phantoms, a section of this article is devoted to examples from the author's own research group. Finally the author explains an unexpected finding during the course of preparing for this article that the phantoms from the past 50 years followed a pattern of exponential growth. The review ends on a brief discussion of future research needs (a supplementary file '3DPhantoms.pdf' to figure 15 is available for download that will allow a reader to interactively visualize the phantoms in 3D).

Estimation of variability of specific absorption rate with physical description of children exposed to electromagnetic field in the VHF band

Recently, there has been an increasing concern regarding the effects of electromagnetic waves on the health of humans. The safety of radio frequency electromagnetic fields (RF-EMFs) is evaluated by the specific absorption rate (SAR). In recent years, SAR has been estimated by numerical simulation using fine-resolution and anatomically realistic reference whole-body voxel models of people of various ages. The variation in SAR with a change in the physical features of a real person is hardly studied, although every person has different physical features. In this study, in order to estimate the individual variability in SAR of persons, we obtained considerable 3D body shape data from actual three-year-old children and developed several homogeneous models of these children. The variability in SAR of the homogeneous models of three-year-old children for whole-body exposure to RF electromagnetic fields in the very high frequency (VHF) band calculated using the finite-difference time-domain method has been described.

Requirements for reliable worst-case assessment of human exposure to RF electromagnetic fields with known uncertainty

This paper gives an overview of the current state of science on numerical and experimental assessment of the exposure of the human body to radio frequency electromagnetic fields. Differences of near- and far-field exposure conditions are discussed with respect to the requirements on the correct representation of the human body and on the numerical or experimental technique. The general requirements for the application of these techniques on the assessment of human exposure are defined, and a combined numerical and experimental approach is proposed that allows the evaluation of the worst-case absorption considering the variability of exposure situations and anatomical properties with known uncertainty.

FDTD calculations of specific energy absorption rate in a seated voxel model of the human body from 10 MHz to 3 GHz

Finite-difference time-domain (FDTD) calculations have been performed to investigate the frequency dependence of the specific energy absorption rate (SAR) in a seated voxel model of the human body. The seated model was derived from NORMAN (NORmalized MAN), an anatomically realistic voxel phantom in the standing posture with arms to the side. Exposure conditions included both vertically and horizontally polarized plane wave electric fields between 10 MHz and 3 GHz. The resolution of the voxel model was 4 mm for frequencies up to 360 MHz and 2 mm for calculations in the higher frequency range. The reduction in voxel size permitted the calculation of SAR at these higher frequencies using the FDTD method. SAR values have been calculated for the seated adult phantom and scaled versions representing 10-, 5- and 1-year-old children under isolated and grounded conditions. These scaled models do not exactly reproduce the dimensions and anatomy of children, but represent good geometric information for a seated child. Results show that, when the field is vertically polarized, the sitting position causes a second, smaller resonance condition not seen in resonance curves for the phantom in the standing posture. This occurs at approximately 130 MHz for the adult model when grounded. Partial-body SAR calculations indicate that the upper and lower regions of the body have their own resonant frequency at approximately 120 MHz and approximately 160 MHz, respectively, when the grounded adult model is orientated in the sitting position. These combine to produce this second resonance peak in the whole-body averaged SAR values calculated. Two resonance peaks also occur for the sitting posture when the incident electric field is horizontally polarized. For the adult model, the peaks in the whole-body averaged SAR occur at approximately 180 and approximately 600 MHz. These peaks are due to resonance in the arms and feet, respectively. Layer absorption plots and colour images of SAR in individual voxels show the specific regions in which the seated human body absorbs the incident field. External electric field values required to produce the ICNIRP basic restrictions were derived from SAR calculations and compared with ICNIRP reference levels. This comparison shows that the reference levels provide a conservative estimate of the ICNIRP whole-body averaged SAR restriction, with the exception of the region above 1.4 GHz for the scaled 1-year-old model.

Resonance behaviour of whole-body averaged specific energy absorption rate (SAR) in the female voxel model, NAOMI

Finite-difference time-domain (FDTD) calculations have been performed of the whole-body averaged specific energy absorption rate (SAR) in a female voxel model, NAOMI, under isolated and grounded conditions from 10 MHz to 3 GHz. The 2 mm resolution voxel model, NAOMI, was scaled to a height of 1.63 m and a mass of 60 kg, the dimensions of the ICRP reference adult female. Comparison was made with SAR values from a reference male voxel model, NORMAN. A broad SAR resonance in the NAOMI values was found around 900 MHz and a resulting enhancement, up to 25%, over the values for the male voxel model, NORMAN. This latter result confirmed previously reported higher values in a female model. The effect of differences in anatomy was investigated by comparing values for 10-, 5- and 1-year-old phantoms rescaled to the ICRP reference values of height and mass which are the same for both sexes. The broad resonance in the NAOMI child values around 1 GHz is still a strong feature. A comparison has been made with ICNIRP guidelines. The ICNIRP occupational reference level provides a conservative estimate of the whole-body averaged SAR restriction. The linear scaling of the adult phantom using different factors in longitudinal and transverse directions, in order to match the ICRP stature and weight, does not exactly reproduce the anatomy of children. However, for public exposure the calculations with scaled child models indicate that the ICNIRP reference level may not provide a conservative estimate of the whole-body averaged SAR restriction, above 1.2 GHz for scaled 5- and 1-year-old female models, although any underestimate is by less than 20%.

Effects of posture on FDTD calculations of specific absorption rate in a voxel model of the human body

A change in the posture of the human body can significantly affect the way in which it absorbs radiofrequency electromagnetic radiation. To study this, an anatomically realistic model of the body has been modified to develop new voxel models in postures other than the standard standing position with arms to the side. These postures were sitting, arms stretched out horizontally to the side and vertically above the head. Finite-difference time-domain (FDTD) calculations of the whole-body averaged specific energy absorption rate (SAR) have been performed from 10 MHz to 300 MHz at a resolution of 4 mm. Calculations show that the effect of a raised arm above the head posture was to increase the value of the whole-body averaged SAR at resonance by up to 35% when compared to the standard, arms by the side position. SAR values, both whole-body averaged and localized in the ankle, were used to derive the external electric field values required to produce the SAR basic restrictions of the ICNIRP guidelines. It was found that, in certain postures, external electric field reference levels alone would not provide a conservative estimate of localized SAR exposure and it would be necessary to invoke secondary reference levels on limb currents to provide compliance with restrictions.

RF dosimetry: a comparison between power absorption of female and male numerical models from 0.1 to 4 GHz

Realistic numerical models of human subjects and their surrounding environment represent the basic points of radiofrequency (RF) electromagnetic dosimetry. This also involves differentiating the human models in men and women, possibly with different body shapes and postures. In this context, the aims of this paper are, firstly, to propose a female dielectric anatomical model (fDAM) and, secondly, to compare the power absorption distributions of a male and a female model from 0.1 to 4 GHz. For realizing the fDAM, a magnetic resonance imaging tomographer to acquire images and a recent technique which avoids the discrete segmentation of body tissues into different types have been used. Simulations have been performed with the FDTD method by using a novel filtering-based subgridding algorithm. The latter is applied here for the first time to dosimetry, allowing an abrupt mesh refinement by a factor of up to 7. The results show that the whole-body-averaged specific absorption rate (WBA-SAR) of the female model is higher than that of the male counterpart, mainly because of a thicker subcutaneous fat layer. In contrast, the maximum averaged SAR over 1 g (1gA-SAR) and 10 g (10gA-SAR) does not depend on gender, because it occurs in regions where no subcutaneous fat layer is present.

Facilitating arrhythmia simulation: the method of quantitative cellular automata modeling and parallel running

Background

Many arrhythmias are triggered by abnormal electrical activity at the ionic channel and cell level, and then evolve spatio-temporally within the heart. To understand arrhythmias better and to diagnose them more precisely by their ECG waveforms, a whole-heart model is required to explore the association between the massively parallel activities at the channel/cell level and the integrative electrophysiological phenomena at organ level.

Methods

We have developed a method to build large-scale electrophysiological models by using extended cellular automata, and to run such models on a cluster of shared memory machines. We describe here the method, including the extension of a language-based cellular automaton to implement quantitative computing, the building of a whole-heart model with Visible Human Project data, the parallelization of the model on a cluster of shared memory computers with OpenMP and MPI hybrid programming, and a simulation algorithm that links cellular activity with the ECG.

Results

We demonstrate that electrical activities at channel, cell, and organ levels can be traced and captured conveniently in our extended cellular automaton system. Examples of some ECG waveforms simulated with a 2-D slice are given to support the ECG simulation algorithm. A performance evaluation of the 3-D model on a four-node cluster is also given.

Conclusions

Quantitative multicellular modeling with extended cellular automata is a highly efficient and widely applicable method to weave experimental data at different levels into computational models. This process can be used to investigate complex and collective biological activities that can be described neither by their governing differentiation equations nor by discrete parallel computation. Transparent cluster computing is a convenient and effective method to make time-consuming simulation feasible. Arrhythmias, as a typical case, can be effectively simulated with the methods described.

A semi-automatic method for developing an anthropomorphic numerical model of dielectric anatomy by MRI

Complex permittivity values have a dominant role in the overall consideration of interaction between radiofrequency electromagnetic fields and living matter, and in related applications such as electromagnetic dosimetry. There are still some concerns about the accuracy of published data and about their variability due to the heterogeneous nature of biological tissues. The aim of this study is to provide an alternative semi-automatic method by which numerical dielectric human models for dosimetric studies can be obtained. Magnetic resonance imaging (MRI) tomography was used to acquire images. A new technique was employed to correct nonuniformities in the images and frequency-dependent transfer functions to correlate image intensity with complex permittivity were used. The proposed method provides frequency-dependent models in which permittivity and conductivity vary with continuity--even in the same tissue--reflecting the intrinsic realistic spatial dispersion of such parameters. The human model is tested with an FDTD (finite difference time domain) algorithm at different frequencies; the results of layer-averaged and whole-body-averaged SAR (specific absorption rate) are compared with published work, and reasonable agreement has been found. Due to the short time needed to obtain a whole body model, this semi-automatic method may be suitable for efficient study of various conditions that can determine large differences in the SAR distribution, such as body shape, posture, fat-to-muscle ratio, height and weight.

Specific absorption rate and temperature elevation in a subject exposed in the far-field of radio-frequency sources operating in the 10-900-MHz range

The exposure of a subject in the far field of radiofrequency sources operating in the 10-900-MHz range has been studied. The electromagnetic field inside an anatomical heterogeneous model of the human body has been computed by using the finite-difference time-domain method; the corresponding temperature increase has been evaluated through an explicit finite-difference formulation of the bio-heat equation. The thermal model used, which takes into account the thermoregulatory system of the human body, has been validated through a comparison with experimental data. The results show that the peak specific absorption rate (SAR) as averaged over 10 g has about a 25-fold increase in the trunk and a 50-fold increase in the limbs with respect to the whole body averaged SAR (SARWB). The peak SAR as averaged over 1 g, instead, has a 30- to 60-fold increase in the trunk, and up to 135-fold increase in the ankles, with respect to SARWB. With reference to temperature increases, at the body resonance frequency of 40 MHz, for the ICNIRP incident power density maximum permissible value, a temperature increase of about 0.7 degrees C is obtained in the ankles muscle. The presence of the thermoregulatory system strongly limits temperature elevations, particularly in the body core.

SARs for pocket-mounted mobile telephones at 835 and 1900 MHz

Increasingly, mobile telephones are becoming pocket-sized and are being left in the shirt pocket with a connection to the ear for hands-free operation. We have considered an anatomic model of the chest and a planar phantom recommended by US FCC to compare the peak 1 and 10 g SARs for four typical cellular telephones, two each at 835 and 1900 MHz. An agreement within +/- 10% is obtained between calculated and experimental 1 and 10 g SARs for various separations (2-8 mm) from the planar phantom used to represent different thicknesses of the clothing both for the antenna away from or turned back towards the body. Because of the closer placement of the antennas relative to the body, the peak 1 and 10 g SARs are considerably higher (by a factor of 2-7) for pocket-mounted telephones as compared to the SARs obtained using a 6 mm thick plastic ear head model--a procedure presently accepted both in the US and Europe. This implies that a telephone tested for SAR compliance against the model of the head may be severely out of compliance if it were placed in the shirt pocket.

Fine resolution calculations of SAR in the human body for frequencies up to 3 GHz

Finite-difference time-domain (FDTD) calculations of whole-body averaged specific energy absorption rate (SAR) have been performed from 100 MHz to 3 GHz at the basic 2 mm resolution of the voxel (volume pixel) model NORMAN without any rescaling to larger cell sizes. The reduction in the voxel size from previous work allows SAR to be calculated at higher frequencies. Additionally, the calculations have been extended down to 10 MHz, covering the whole-body resonance regions at a resolution of 4 mm. As well as for the adult phantom, SAR values are calculated for scaled versions representing 10-, 5- and 1-year-old children for both grounded and isolated conditions. External electric field levels are derived from limits of whole-body averaged SAR and localized SAR in the ankle, and compared with NRPB investigation levels and ICNIRP reference levels. The ICNIRP field reference levels alone would not provide a conservative estimate of the localized SAR exposure in the leg for grounded conditions. It would be necessary to invoke the secondary reference level on limb current to provide compliance with basic restrictions on localized SAR averaged over 10 g.

Head and neck resonance in a rhesus monkey--a comparison with results from a human model

The use of primates for examining the effects of electromagnetic radiation on behavioural patterns is well established. Rats have also been used for this purpose. However, the monkey is of greater interest as its physiological make-up is somewhat closer to that of the human. Since the behavioural effects are likely to occur at lower field strengths for resonant absorption conditions for the head and neck, the need for determination of resonance frequencies for this region is obvious. Numerical techniques are ideal for the prediction of coupling to each of the organs, and accurate anatomically based models can be used to pinpoint the conditions for maximum absorption in the head in order to focus the experiments. In this paper we use two models, one of a human male and the other of a rhesus monkey, and find the mass-averaged power absorption spectra for both. The frequencies at which highest absorption (i.e. resonance) occurs in both the whole body and the head and neck region are determined. The results from these two models are compared for both E-polarization and k-polarization, and are shown to obey basic electromagnetic scaling principles.