International Journal of Antennas and Propagation
Volume 2018 (2018), Article ID 2969237, 9 pages
https://doi.org/10.1155/2018/2969237

Research Article

Modeling Human Body Using Four-Pole Debye Model in Piecewise Linear Recursive Convolution FDTD Method for the SAR Calculation in the Case of Vehicular Antenna

Ammar Guellab and Qun Wu

Harbin Institute of Technology, Xidazhi Street, Harbin, China

Correspondence should be addressed to Ammar Guellab; ammarguellab400@hit.edu.cn and Qun Wu; qwu@hit.edu.cn

Received 16 September 2017; Revised 7 December 2017; Accepted 28 December 2017; Published 19 April 2018

Academic Editor: Rodolfo Araneo

Copyright © 2018 Ammar Guellab and Qun Wu. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We propose an efficient finite difference time domain (FDTD) method based on the piecewise linear recursive convolution (PLRC) technique to evaluate the human body exposure to electromagnetic (EM) radiation. The source of radiation considered in this study is a high-power antenna, mounted on a military vehicle, covering a broad band of frequency (100 MHz–3 GHz). The simulation is carried out using a nonhomogeneous human body model which takes into consideration most of the internal body tissues. The human tissues are modeled by a four-pole Debye model which is derived from experimental data by using particle swarm optimization (PSO). The human exposure to EM radiation is evaluated by computing the local and whole-body average specific absorption rate (SAR) for each occupant. The higher in-tissue electric field intensity points are localized, and the SAR values are compared with the crew safety standard recommendations. The accuracy of the proposed PLRC-FDTD approach and the matching of the Debye model with the experimental data are verified in this study.

1. Introduction

Nowadays, the battlefield is making more use of the electromagnetic spectrum to satisfy diverse operational needs that range from high-rate tactical links to broadband jammers. This calls for the use of vehicular antennas that transmit high powers, putting the life of crew personnel at stake through high electromagnetic (EM) exposure [1]. Such high-power radiation can occur in a broad frequency band (HF, VHF, and UHF). IEEE Technical Committee 95 (IEEE-TC95) proposed a standard for military workplaces whose purpose is to provide exposure limits to assure the personnel safety in a military workplace and provide protection against unfavorable effects of the electromagnetic radiation on the human body [2]. This standard expresses its recommendations as dosimetry reference limits (DRLs) which can be expressed by the within-tissue electric field strength or the specific absorption rate (SAR).

With the aim of protecting the crew against high-power EM radiation, much research is devoted to compute the level of the induced SAR from various EM sources. Accurate and efficient computational methods are sought after for this purpose [3–6]. This can be so challenging and complicated due to the human body and the vehicle structure complexity especially when a broad frequency band is considered.

A numerical simulation for the computation of the SAR in vehicle passengers due to an onboard 900 MHz transmission system is carried out in [3]. The FDTD method is used for the analysis of the SAR from a cell phone inside a vehicle in [4]. However, those studies are concerned with a single narrow band of frequency. Diao [5] uses the FDTD method for the estimation of the induced SAR from multiple communication devices, which occupy different frequencies inside the vehicle. In [6] the FDTD simulation of a complex structure vehicle is carried out to investigate the effect of the passenger’s number on the SAR level for many frequencies.

The previously cited works obtain a rigorous estimation for the SAR level induced by multiple sources. Nevertheless, they focus on isolated and narrow bands of frequency and use human body models that do not take into account the frequency dependency of electric characteristics of the human tissues.

A frequency-dependent model for the human tissue has been provided by Gabriel et al. [7], which is based on a four-pole Cole-Cole equation. This model leads to a highly precise evaluation of the complex electric characteristics (permittivity and conductivity) of the dispersive media constituting the human tissue. Nevertheless, the complexity of this model makes it disadvantageous for the implementation of time-domain numerical simulation.

Gabriel’s Cole-Cole model is reduced to a two-pole Debye model by using the least squares fitting technique in [8, 9]. The two-pole Debye tissue model is more suitable for the time domain methods but is less accurate than the Cole-Cole model. In [10], a fourth-order Debye model is used in the convolutional-based FDTD for the modeling of the electromagnetic waves propagation in the human head tissues. Based on this study, the fourth-order Debye model is more accurate than the two-pole Debye model and is simpler than the Cole-Cole model. Therefore, it combines the high accuracy of computations and the simplicity of implementation.

Despite the fact that the high-order Debye model matches perfectly with the experimental dielectric characteristics, it leads to computational complexity and an accumulative numerical dispersion. In the literature, many works are devoted to the development of highly accurate methods for frequency-dependent media. In [10], a general recursive convolution FDTD is applied for the simulation of the fourth-order Debye model. However, this method is less accurate. In [11] the alternating direction implicit finite difference time domain (ADI-FDTD) method is extended to the simulation of the Debye dispersive media. The advantage of the ADI-FDTD is its unconditional stability. However, it suffers from the loss of precision. In [12], a simple trapezoidal recursive (TRC) technique is used in the FDTD analysis of frequency dependent media, where it is proven that the TRC method is more accurate than the RC technique. For more accuracy, a piecewise linear recursive convolution (PLRC) technique is used in [13]. The PLRC technique is an efficacious method for dealing with the dispersive models (Debye, Lorentz, and Drude) [13].

In this work, a highly accurate four-pole Debye model is developed, which perfectly fits with the experimental data [14] of the whole-body tissues in a broad band of frequency ranging from 100 MHz to 3 GHz. For this purpose, we applied the particle swarm optimization (PSO) algorithm for the optimization role. Then, we proposed an implementation of a three-dimension FDTD based on the piecewise linear recursive convolution (PLRC) technique for the simulation of the four-pole Debye model. The agreement of the Debye model with the experimental dielectric properties of each tissue of the human body and the accuracy of the proposed FDTD method are verified.

After that, a typical military vehicle, commonly known as (Humvee) is considered. It is equipped with a large antenna used for radio-transmission applications. A reduced complexity model is used for the vehicle. Human bodies are present inside where they are exposed to high-power electromagnetic radiation issued from the antenna. We applied the developed PLRC-FDTD method for the study of the electromagnetic waves’ behavior in this case and to compute the human bodies’ exposure to those radiations through the evaluation of the specific absorption rate (SAR).

2. Formulations

2.1. Four-Pole Debye Model

With the aim of modeling the response of the human body issues to the electromagnetic radiation at a broad band of microwave frequencies, the measured permittivity and conductivity data for each tissue [14] are approximated by a fourth-order Debye dispersive model, as expressed in the following equation:

where is the angular frequency, is the free space permittivity, and are the real part and the imaginary part of the complex relative permittivity () of the Debye dispersive media where denotes the electric conductivity, is the permittivity at the infinite frequency, K is the number of Debye modes, and are, respectively, the magnitude and the relaxation time of the kth Debye dispersion mode, and is the static conductivity.

2.2. The Particle Swarm Optimization

In this paper, the Matlab toolbox based on the particle swarm optimization (PSO) algorithm is used to accurately fit the Debye parameterized model with a set of measured data [14], by minimizing the cost function ζ which is the mean of the relative errors between the data measurements and the Debye computed model at each measured frequency, as expressed in the following equation [10].

where , and are, respectively, the computed permittivity and conductivity using the Debye model. and are the measured permittivity and conductivity, respectively [14].

The application of the particle swarm algorithm implies the minimization of the cost function . The experimental data [14] of 19 tissues are fitted over the frequency ranging from 100 MHz to 3 GHz with the fourth-order Debye model.

The PSO algorithm’s flow involves first the creation of a random swarm of 2000 particles. Each particle has a position which is represented by a set of Debye parameters and a velocity which is the used to compute the next position. As a constraint, all the parameters are nonnegative numbers. At each iteration of the PSO algorithm, the cost function is evaluated at all particles to record the location of the best particle in the neighborhood of each particle and the best one over all particles. For each particle, the neighborhood is the first 500 nearest particles. The velocity of each particle is updated as a weighted sum of the previous velocity and the position of the neighborhood best particle and the global best particle. Then, the position of each particle is updated as a function of the previous position and the velocity. Over consecutive iterations, the swarm of particles converge toward the optimal solution with the lowest cost function. The maximum number of iterations is 10,000.

For the aim of comparison, the genetic algorithm is performed in such a manner to have the same resource consumption as the PSO algorithm. The population size is 2000. The maximum number of generations is 10,000.

Table 1 lists the four-pole Debye parameters of each tissue (due to the limited space, only a few of them are listed) and the evaluated errors of each model expressed by (2). The fitting results of some tissues are plotted in Figure 1.

Table 1: Parameters of the four-pole Debye model of human tissues.

Tissues
Blood
Bone
Cerebellum
Fat
Heart
Kidney
Liver
Lung
Muscle
White matter

1.96E + 01
7.11E + 00
1.58E + 01
2.07E + 00
1.80E + 01
1.80E + 01
1.46E + 01
6.68E + 00
2.11E + 01
1.81E + 01

1.76E + 04
5.38E + 04
2.74E + 01
1.74E + 04
2.07E + 06
1.03E + 05
1.35E + 04
9.21E + 06
6.67E + 05
1.03E + 04

4.59E + 03
2.52E + 02
3.40E + 01
1.09E + 01
1.57E + 03
1.32E + 03
4.42E + 02
4.62E + 02
1.95E + 03
1.93E + 02

1.17E + 02
2.99E + 01
2.05E + 02
2.56E + 00
1.74E + 02
3.72E + 01
3.10E + 01
6.85E + 01
4.70E + 01
1.96E + 01

3.70E + 01
1.16E + 01
1.22E + 03
1.51E + 00
4.05E + 01
2.55E + 02
5.07E + 01
1.42E + 01
3.77E + 01
4.11E + 01

τ₁ [rad/s]
8.65E − 04
1.49E − 05
9.98E − 04
6.79E − 06
5.48E − 05
3.17E − 06
5.22E − 07
9.17E − 04
1.13E − 05
7.68E − 07

τ₂ [rad/s]
9.73E − 08
4.97E − 08
1.37E − 11
9.92E − 09
6.04E − 08
4.81E − 08
2.13E − 08
4.37E − 08
1.14E − 07
1.35E − 08

τ₃ [rad/s]
5.86E − 09
1.76E − 09
3.50E − 09
7.51E − 12
4.45E − 09
1.56E − 11
1.64E − 11
4.17E − 09
1.37E − 09
2.23E − 11

τ₄ [rad/s]
1.41E − 11
1.68E − 11
4.56E − 08
8.82E − 10
1.52E − 11
4.20E − 09
1.72E − 09
1.41E − 11
1.68E − 11
1.27E − 09

σs [S/m]
1.064
0
1.13E − 01
0
0
0
0
0
0
0

Errors

Debye2 (PSO)
4.91E − 03
4.02E − 02
1.27E − 02
8.12E − 02
1.17E − 02
1.15E − 02
1.96E − 02
1.01E − 02
2.78E − 02
2.43E − 02

Debye4 (PSO)
1.30E − 03
6.05E − 03
2.04E − 03
4.10E − 03
2.98E − 03
2.57E − 03
2.55E − 03
2.01E − 03
3.41E − 03
4.56E − 03

Debye4 (GA)
1.61E − 02
8.58E − 03
4.52E − 02
2.52E − 02
2.79E − 02
2.66E − 01
7.21E − 02
4.01E − 03
2.61E − 02
3.54E − 02

Cole-Cole (PSO)
1.80E − 04
3.87E − 04
4.98E − 04
5.09E − 03
1.23E − 04
4.03E − 04
4.16E − 04
1.40E − 04
2.13E − 04
2.31E − 03

(a)

(b)

Figure 1: Properties of human body tissues using Debye model compared with the experimental data. (a) Relative permittivity. (b) Conductivity.

From Figure 1 and Table 1, it is shown that the Cole-Cole model offers the best agreement with the experimental data [7], then the four-pole Debye model computed by using the PSO algorithm offers almost the same performance as the Cole-Cole model, then the two-pole Debye offers lower fitting with the experimental data. In addition, if we consider the same model, the PSO algorithm offers more performance than the Genetic algorithm.

2.3. Formulation of the PLRC-FDTD

To investigate the electromagnetic wave behavior through a Debye dispersive media, the relationship relating the electric field intensity to the electric induction by the frequency-dependent permittivity is converted to the time domain by using the inverse Fourier transform.

where is the convolutional product, is the Dirac delta function and is the unit step function. The substitution of in (3) yields the following:

where represents the polarization density of the kth Debye mode. is conductive loss term.

In the discretized mode, the polarization density and the conductive loss become

and in (8) can be computed by a recursive formula:

The PLRC method is based on assuming that the electric field follows a piecewise linear behavior within each discretization interval, then it can be expressed as follows.

By the substitution of (10) into the (9), the expression of the polarization density and the conductive loss become as follows.

where

By solving (5) after substitution of (11), the updated equation of the electric field intensity is resumed as follows:

where

The advantage of using PLRC technique compared to the convolution-based FDTD method [10] resides in increasing the accuracy, by assuming that the electric field follows a piecewise linear function of the time, whereas the convolution-based method considers it as a constant in every discretization interval.

2.4. Validation

For the validation issues, a simple EM problem for which the analytical solution is well known is solved using the present method. Then, the results are compared with the analytical solution. This EM problem consists of a TEM Gaussian-pulse plane wave which is propagating along the x-direction and reflected by an air-dispersive media plane interface, where the chosen Debye parameters are those characterizing the fat in Table 1. The computational environment is surrounded by a PEC boundary condition from the upside and the downside, and a PMC boundary from the left and right sides, and ten cells of perfectly matched layer (PML) from the front and back sides, as shown in Figure 2. The Gaussian pulse is parameterized in such a way to cover the frequency band of interest (100 MHz–3 GHz). A uniform Yee grid is used with a space increment of 4 mm, which provides a minimum resolution of 25 cells per wavelength (at 3 GHz). Figure 2 shows the incident and the reflected wave from the interface in the time domain. The reflection coefficient computed by the PLRC, TRC, and RC methods is compared with the analytical solution in the frequency domain as shown in Figure 3(a). The relative error calculated using (15), which evaluates the accuracy of the used numerical methods, is shown in Figure 3(b).

where is the analytically computed reflection coefficient and is that computed numerically by using one of the techniques: PRLC, TRC, or RC.

Figure 2: Plane wave normal incident on a vacuum/Debye model interface.

(a)

(b)

Figure 3: Comparison of the computed reflection coefficient by RC, TRC, and PLRC techniques with the theoretical values. (a) Magnitude and phase. (b) Relative errors.

Figure 2 shows that the three methods present an excellent agreement with the analytical solution. However, both PLRC and TRC methods demonstrate a lower error compared to the RC method along the frequency band of interest as shown in Figure 3. This is explained by the accuracy order of each method, where PLRC and TRC methods are of second-order accuracy, and the RC is of first order.

3. Simulation Results

A military vehicle (Humvee) with a mounted antenna on its rear part is modeled by using a reduced complexity model. The complexity is reduced by eliminating the geometrical details which are far smaller than the smallest studied wavelength. Then, a crew composed of one, two, or four people is added. The geometrical characteristics of the human bodies are acquired from a large dataset of human tissue model [15] which is evaluated and accepted by the IEEE International Committee on Electromagnetic Safety for the calculation of the specific absorption rates. All the geometric entities are converted from vertex data to voxel model which is appropriate to the FDTD grid as shown in Figure 4. Each crew member is represented by a heterogeneous phantom which is formed by taking into consideration the geometrical characteristics of all the organs in [15]. Then, the phantoms are postured in a seating position by using geometrical transformation for the limbs.

Figure 4: 3D model of the Humvee with four human bodies inside.

The radiating antenna is fed with a wideband Gaussian pulse covering a broad frequency band from 0.1 GHz to 3 GHz. The computational domain is  m in length,  m in width, and  m in height. In order to allow the simulation for an open space, a perfectly matched layer (PML) of 20 cells surrounds the computational space. The simulation is carried out for the sufficient time allowing the total attenuation of the electromagnetic waves. The space discretization increment is  mm, which provides a minimum resolution of 25 points per wavelength, which is sufficient for obtaining accurate results as proven in Section 2.4.

The simulation is performed by using the Matlab parallel computing toolbox. The electric and magnetic flux densities are represented by 1360 × 715 × 875 sized three-dimensional variables covering all the computational space, whereas the electric and the magnetic field intensity are computed just in the nonvacuum space. The simulation requires 41 Go of memory with an average of 55 octets by cell. The FDTD algorithm is run for 10,000 iterations, the equivalent time of 50 ns which is the necessary time for the electromagnetic waves to vanish throughout the PML absorbing boundary. In order to showcase the electric field distribution at various frequencies, we applied the discrete Fourier transformation for each cell in the computational environment.

First, we consider the case of a crew composed of four people. The human body tissues exposure to electromagnetic radiation is evaluated by the within-tissue electric field strength as reported in Figure 5 which shows the computed distribution of the electric field intensity in an arbitrary chosen vertical slice passing through two passengers at the frequency of 200 MHz. It is clear that the electric field intensity in the vacuum is higher than that inside the human body. This fact is due to the conductor aspect of the human tissues.

Figure 5: Electric field distribution at the frequency 200 MHz.

Then, to perform an electromagnetic resonance analysis of the human body, particular locations in the human body have been chosen as detection points (head, heart, and ankle) as marked in Figure 5 by white marks. Figure 6 shows the spectrum of the electric field strength at the marked points. We notice the existence of resonance modes at the interval 100 MHz–1 GHz, depending on the geometry of each part of the body and the dielectric characteristics of its tissue. For the brain, the main resonance frequency is observed at 0.25 GHz; for the heart and the ankle, the main resonant frequency is 0.6 GHz.

Figure 6: Electric field strength at the head, the heart, and the ankle of the 3rd occupant.

Second, the simulation is repeated for the cases of a crew composed of one and two occupants. The computed electromagnetic dosimetry is compared with the recommended dosimetry reference limits (DRLs) established by the IEEE standard for the military workspaces [2], where different local SAR limits are identified for some regions of the body. In [2], the environments are classified into two classes: unrestricted zones, where the access is allowed for all personnel and the DRLs are strict; and restricted zones, where only the informed personnel are authorized to access and the DRLs are more relaxed than the unrestricted areas. For unrestricted environments, the standard [2] defines the limits for the whole-body average SAR as 0.08 W/kg, the 10 g local SAR in the head and the trunk as 2 W/kg, and the local SAR for the arms and legs as 4 W/kg.

Table 2 summarizes the maximum values of the computed SAR for a normalized radiated power of 1 W CW, and the frequencies for which those maximums are detected, for one, two, and four members of the crew and for each occupant in the crew (1, 2, 3, and 4 denote the driver, front passenger, left rear passenger, and rear right passenger, resp.).

Table 2: Computed SAR For AT 1 W CW.

SAR measure
Occupant 1
Occupant 2
Occupant 3
Occupant 4

One occupant

Whole-body average
Maximum (W/kg)
4.174E − 05

Frequency (MHz)
600

Local SAR for head and trunk
Maximum (W/kg)
1.236E − 03

Frequency (MHz)
600

Local SAR for limbs
Maximum (W/kg)
1.362E − 03

Frequency (MHz)
600

Two occupants

Whole-body average
Maximum (W/kg)
2.904E − 05
2.062E − 05

Frequency (MHz)
600
600

Local SAR for head and trunk
Maximum (W/kg)
1.015E − 03
8.671E − 04

Frequency (MHz)
600
600

Local SAR for limbs
Maximum (W/kg)
1.260E − 03
5.757E − 04

Frequency (MHz)
600
600

Four occupants

Whole-body average
Maximum (W/kg)
1.59E − 05
8.82E − 06
2.92E − 05
8.30E − 06

Frequency (MHz)
600
600
600
600

Local SAR for head and trunk
Maximum (W/kg)
4.674E − 4
7.417E − 4
8.606E − 4
3.277E − 4

Frequency (MHz)
650
600
600
675

Local SAR for limbs
Maximum (W/kg)
7.433E − 4
5.775E − 4
1.207E − 4
4.388E − 4

Frequency (MHz)
675
600
600
600

The obtained whole-body average and local SAR values for a radiation power of 1 W are significantly lower than the recommended limits. Since the computed SAR values are proportional to the radiated power, those values can be scaled to a higher level without exceeding the recommended limits.

From the obtained results, it is clear that the SAR level depends on the total number of occupants. For example, if we consider the first occupant, the SAR level is most important in the case of a single occupant. Also, we notice that in the case of four occupant crew, the left rear occupant is more affected by the antenna radiation than the other occupants. This can be explained by the fact that this occupant is nearest to the antenna.

According to those results, the maximum radiation power for which the SAR limits might not be attained is 1.62 KW, which is commonly greater than the usually used power for military operations on this type of vehicles.

4. Conclusion

In this paper, an efficient 3D PLRC-FDTD method is developed to provide an accurate simulation of the electromagnetic wave propagation in the human body to better assess the specific absorption rates in environments characterized by high levels of radiation.

The PLRC-FDTD was used for the simulation of a 3D voxel model of the military vehicle (Humvee) with a mounted transmission antenna in the presence of human bodies. The particle swarm optimization was used to derive four-pole Debye models for the whole-body tissues in the frequency band from 0.1 GHz to 3 G Hz. In addition to its accurate matching with the experimental electric characteristics of the human tissues, the Debye model also offers an ease of implementation.

This allows the computation of the in-tissue electromagnetic field’s intensity and consequently the calculation of the specific absorbing rate than the localization of regions with a high exposure to electromagnetic radiation in a broad band of frequency. The PLRC method, which is second-order accuracy, fulfills considerably improved precision over the RC method.

The computed SAR levels for 1 W CW radiated power are significantly lower than the recommended limits. These results also demonstrate that those limits cannot be attended unless the radiated power exceeds 1.62 kW, which is uncommon regarding the typically used equipment in the military vehicle.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Acknowledgments

This project is partly supported by the National Natural Science Foundation of China (Grant no. 61371044).

References

D. Garrido, M. Ignatenko, and D. S. Filipovic, “Computational study of electromagnetic exposure of military personnel in a Humvee,” in 2014 United States National Committee of URSI National Radio Science Meeting (USNC-URSI NRSM), p. 1, Boulder, CO, USA, 2014.

“IEEE standard for military workplaces--force health protection regarding personnel exposure to electric, magnetic, and electromagnetic fields, 0 Hz to 300 GHz,” in IEEE Std C95.1-2345-2014, pp. 1–69, New York, NY, USA, 2014, IEEE.

A. R. Ruddle, “Computed SAR levels in vehicle occupants due to on-board transmissions at 900 MHz,” in 2009 Loughborough Antennas & Propagation Conference, pp. 137–140, Loughborough, UK, 2009, IEEE.

G. Anzaldi, F. Silva, M. Fernandez, M. Quilez, and P. J. Riu, “Initial analysis of SAR from a cell phone inside a vehicle by numerical computation,” IEEE Transactions on Biomedical Engineering, vol. 54, no. 5, pp. 921–930, 2007.

Y. Diao, “SAR evaluation for multiple wireless communication devices inside a vehicle,” in 2013 International Symposium on Electromagnetic Theory, pp. 626–629, Hiroshima, Japan, 2013.

A. R. Ruddle, L. Low, J. M. Rigelsford, and R. J. Langley, “Variation of computed in-vehicle SAR with number and location of occupants at commonly used communications frequencies,” in 10th International Symposium on Electromagnetic Compatibility, pp. 756–761, York, UK, 2011.

S. Gabriel, R. W. Lau, and C. Gabriel, “The dielectric properties of biological tissues: III. Parametric models for the dielectric spectrum of tissues,” Physics in Medicine & Biology, vol. 41, no. 11, pp. 2271–2293, 1996.

T. Wuren, T. Takai, M. Fujii, and I. Sakagami, “Effective 2-Debye-pole FDTD model of electromagnetic interaction between whole human body and UWB radiation,” IEEE Microwave and Wireless Components Letters, vol. 17, no. 7, pp. 483–485, 2007.

M. Fujii, R. Fujii, R. Yotsuki, T. Wuren, T. Takai, and I. Sakagami, “Exploration of whole human body and UWB radiation interaction by efficient and accurate two-Debye-pole tissue models,” IEEE Transactions on Antennas and Propagation, vol. 58, no. 2, pp. 515–524, 2010.

S. Mustafa, A. M. Abbosh, and P. T. Nguyen, “Modeling human head tissues using fourth-order Debye model in convolution-based three-dimensional finite-difference time-domain,” IEEE Transactions on Antennas and Propagation, vol. 62, no. 3, pp. 1354–1361, 2014.

S. G. Garcia, R. G. Rubio, A. R. Bretones, and R. G. Martin, “Extension of the ADI-FDTD method to debye media,” IEEE Transactions on Antennas and Propagation, vol. 51, no. 11, pp. 3183–3186, 2003.

J. Shibayama, R. Ando, A. Nomura, J. Yamauchi, and H. Nakano, “Simple trapezoidal recursive convolution technique for the frequency-dependent FDTD analysis of a Drude–Lorentz model,” IEEE Photonics Technology Letters, vol. 21, no. 2, pp. 100–102, 2009.

I. Giannakis and A. Giannopoulos, “A novel piecewise linear recursive convolution approach for dispersive media using the finite-difference time-domain method,” IEEE Transactions on Antennas and Propagation, vol. 62, no. 5, pp. 2669–2678, 2014.

C. Gabriel, Compilation of the Dielectric Properties of Body Tissues at RF and Microwave Frequencies, Radiofrequency Radiation Division, Brooks Air Force Base, San Antonio, TX, USA, 1996, http://www.dtic.mil/cgi-bin/GetTRDoc?AD=ADA309764.

S. Makarov, G. Noetscher, and A. Nazarian, Low-frequency electromagnetic modeling for electrical and biological systems using MATLAB, Wiley, Hoboken, NJ, USA, 2016.