Research Article  Open Access
JiHyun Jung, SeYun Kim, "Enhancement of CrossBorehole Pulse Radar Signature on a Partially WaterFilled Tunnel", International Journal of Antennas and Propagation, vol. 2014, Article ID 396835, 9 pages, 2014. https://doi.org/10.1155/2014/396835
Enhancement of CrossBorehole Pulse Radar Signature on a Partially WaterFilled Tunnel
Abstract
Crossborehole pulse radar has been employed to detect a deeply located empty tunnel. In this paper, effects of underground water collected in the bottom of an empty tunnel on crossborehole pulse radar signatures are analyzed numerically. Bscan images, stacks of received pulses, are calculated by applying the finitedifference timedomain (FDTD) method for 6 different heights of water from the bottom to the half height inside an empty tunnel. The most important features of an empty tunnel, the fastest time of peak (TOP) and time of arrival (TOA) extracted from the Bscan images, are slowed considerably depending on the increased height of water inside the tunnel. To compensate the weak TOP like that of an empty tunnel, a relation curve is formulated only utilizing measurable parameters of the fastest TOP and the fastest TOA. Then, a unified curve including the effects of two granites with the low and high dielectric properties is derived to cover widely varied dielectric properties of underground rocks. Based on the fastest TOP of an empty tunnel, the average difference between the fastest TOP of an empty tunnel and that of a partially waterfilled tunnel decreases from 22.92% to 2.59% after enhancement.
1. Introduction
Two types of crossborehole radars have been operated for detecting subversive intrusion tunnels deeply located in Korea since 1980s. One type is steppedfrequency continuouswave (CW) radar characterized by the doubledip pattern at some specific frequencies in received signals [1–6]. Crossborehole CW radar played a key role in detecting the Fourth Tunnel in Korea [2, 4]. The other type is pulse radar [7–15], characterized by fast arrival in received pulse near the depth of an empty tunnel. Recently, crossborehole pulse radar was developed in our group and its detectability of an empty tunnel was confirmed successfully in a wellsuited tunnel test site in Korea [13–15]. However, underground water oozed out of underground rock can be collected in the bottom of an empty tunnel owing to highly weathered, jointed, and fractured underground rock of Korea. To detect such a partially waterfilled tunnel, it should be analyzed properly whether conventional signatures of an empty tunnel can be employed immediately or not.
The pulse waveforms scattered by a fully waterfilled cylindrical tunnel were already calculated analytically using 3dimensional Green’s function [16]. As expected, the first positive peak arrived latest in the central position of the tunnel. TE and TM scattering by a tunnel halffull of water were calculated in the presence of an airearth interface using the 2dimensional method of moment [17]. Imaging of scatterer buried in a lossy half space was developed using a regularizing procedure based on twostep inexactNewton method [18]. A method to detect buried multilayer dielectric cylinders with elliptic crosssections was analyzed using a global optimization procedure called memetic algorithm [19]. Until now, however, crossborehole pulse radar signatures considering various height of collected water inside an empty tunnel are not available. Hence, experiments at a real test site in Korea can be considered utilizing our crossborehole pulse radar. However, the safety of many visitors will render it nearly impossible to fill the insite tunnel with water. Technically, the precise control of the water level in the tunnel is also very difficult since water is continuously oozed out or permeated into the underground rock. It leads us to consider numerical simulations on a partially waterfilled tunnel instead of field experiments. A numerical simulator developed by our group [20] using the finitedifference timedomain (FDTD) [21, 22] is modified.
In this paper, the crossborehole pulse radar operated in a wellsuited tunnel test site is represented by 3dimensional model and the finitedifference timedomain (FDTD) method is applied to the model for numerical analysis. Two separate waterfilled boreholes and an empty tunnel with crosssection of 2.1 m by 2.1 m are surrounded by underground rock, lossy and dispersive granite. Two dipole antennas immersed in the bottom of two waterfilled boreholes are moving upwardly with a uniform distance step. At every depth, the radiated pulse from a transmitting antenna is propagated through the granite and picked by a receiving antenna. By stacking picked pulses at every depth continually, Bscan images are obtained. Numerical simulations are repeated by changing the height of collected water from the bottom to the half height of the tunnel with a uniform interval.
Conventional signatures, the time of peak (TOP) and the time of arrival (TOA), are extracted from calculated Bscan images. Fast arrival of the TOP and TOA, the feature of an empty tunnel, is delayed gradually as the height of water inside the tunnel increases. However, the weak TOP and TOA may not be considered as tunnel signatures because those signatures do not meet the accumulating experience of field experiments related with an empty tunnel. Thus, one may be confused whether the delayed TOP and TOA are tunnel signatures or not. To enhance slowed TOP like the TOP of an empty tunnel, a relation curve is derived using measurable parameters of the fastest TOP and the fastest TOA. By adding the calculated enhancement amount obtained from the relation curve, delayed TOP can be improved similar to the TOP of an empty tunnel. Finally, the relation curve is unified by considering two additional granites which cover the low and high dielectric properties of granite.
In the following section, the setup of 3dimensional numerical model representing the operation of crossborehole pulse radar in a wellsuited test site in Korea is described. In Section 3, numerical simulations are performed and Bscan images of partially waterfilled tunnels are obtained. Then, conventional signatures extracted from Bscan images are interpreted. In Section 4, a relation curve to enhance weak conventional signature is derived. The effects of two additional granites with the low and high dielectric properties on the relation curve are treated properly. Finally, conclusions are summarized in the last section.
2. Numerical Modeling
Exploration of a deeply located subversive intrusion tunnel using crossborehole pulse radar in a wellsuited test site is depicted in Figure 1. An empty tunnel with crosssection about 2.1 m by 2.1 m is located at 73 m deep. Two waterfilled boreholes are separated by 15.9 m. A transmitting antenna, inserted near the bottom of the left borehole, radiates an electromagnetic pulse. The pulse waveform distorted by the tunnel and underground rock is picked out by a receiving antenna at the same depth in the right borehole. Described measurement is performed until more than 5 m below the air and water surface in the boreholes by pulling up both antennas simultaneously.
Threedimensional numerical model of the crossborehole pulse radar in a wellsuited test site is shown in Figure 2(a). In our numerical model, the spatial resolution () is uniformly taken by 0.03 m by considering the electrical property of water and the pulse width of the excited pulse, 10 ns. The relative complex permittivity of water is expressed by Debye’s formula [23] as where is the high frequency permittivity, is the zero frequency permittivity, is the relaxation time, is the conductivity, and is the permittivity in free space. The corresponding parameters of water are ps, and S/m [24].
(a)
(b)
The center of an empty tunnel perpendicularly crossing the vertical plane connecting two boreholes is located at a depth of 73 m. The archshaped upper boundary of the tunnel is approximated into a trapezoid. The width and height of the rectangle are 2.1 m () and 1.5 m (), respectively. And the height of the trapezoid is 0.6 m (). Hence, the total width and height of the tunnel crosssection are 2.1 m () and 2.1 m (), respectively. Proportional to the total height of the tunnel, the height of collected water inside the tunnel is expressed using the height ratio of water inside a tunnel as as illustrated in Figure 2(b). The tunnel is infinitely stretched along the axis. Two boreholes and an empty tunnel are surrounded by a lossy and dispersive rock, granite. The relative complex permittivity of granite is given by Debye’s formula with ps, and S/m [7].
Two boreholes with diameter of 0.18 m () are separated with the distance of 15.9 m () and filled with water. Two dipole antennas separately enclosing transmitter and receiver are immersed in two different waterfilled boreholes. The radius and height of the antennas are 0.03 m () and 3 m (), respectively. And the center of the antennas, feeder region, is enclosed by Teflon (). Each antenna is suspended by a cable which connects the antenna with a ground controller consisting of a computer, a winch, and a power supply. To reduce computational burden, two antennas are winding up only from the depth of 77.75 m to 68.25 m with uniform distance interval of 0.15 m. In each depth, the transmitting antenna excites a Gaussianshaped electromagnetic pulse with the width of 10 ns towards underground rock by employing the onecell gap model [25]. And the receiving antenna in the other borehole collects the propagated pulse in the same depth. The time step () is chosen by 51.9615 ps to satisfy the stability condition [21, 22]. The width, length, and height of the 3dimensional model are 18.9 m (), 3.0 m (), and 9.3 m (), respectively. Finally, as an absorbing boundary condition, the convolutional perfectly matched layer (CPML) [26] is set in the outer region of the 3dimensional model.
3. Conventional Tunnel Signatures
3.1. Time of Peak
FDTD simulations are performed for no tunnel, an empty tunnel, and a partially waterfilled tunnel, respectively. In general, an empty tunnel is not fully filled with water since natural drainage is considered during the construction of manmade intrusive tunnel. Hence, half height of the tunnel () is selected as the height of collected water inside an empty tunnel. Calculated 3 different Bscan images are displayed in Figure 3. The received waveforms of no tunnel in Figure 3(a) are uniform. Thus, the first positive peaks commonly arriving at 280.4086 ns are used as a reference. The Bscan image of the empty tunnel in Figure 3(b) illustrates faster arrival of the first positive peak near the depth of 73 m than 280.4086 ns because excited pulse propagates relatively faster through the air region inside a tunnel than the underground rock. The Bscan image of a partially waterfilled tunnel in Figure 3(c) also displays faster arrival of the first positive peak near the depth of 73 m than 280.4086 ns. However, the fast arrival of the first positive peak slows down considerably compared with that of the empty tunnel. It implies that the signature of an empty tunnel, the fast arrival of the first positive peaks, is affected by collected water inside an empty tunnel.
(a)
(b)
(c)
To find the effects of the stagnant water inside the tunnel on tunnel signatures precisely, additional FDTD simulations are performed for 4 different height ratios of water inside the tunnel (), 0.1, 0.2, 0.3, and 0.4, respectively. Then, the most important tunnel signature, the time of peak (TOP), is interpreted. At every depth, the arrival time of the first positive peak is extracted from 6 different Bscan images and illustrated in Figure 4. In Figure 4, the fastest arrival of TOP profiles marked by a black dot is gradually delayed depending on the increased height ratio of water in a tunnel . The fastest arrival of each TOP profile corresponding from to is 270.9008 ns, 272.0438 ns, 272.9790 ns, 273.6025 ns, 274.1220 ns, and 274.8494 ns, respectively. The cause of the delay can be explained by the growing effect of the water in the tunnel which renders the propagation of pulse to be slowed down. In practice, however, the weak TOP which does not meet the accumulating experience of field experiments related with an empty tunnel may not be considered as a tunnel signature. The delayed TOP means the existence of a small tunnel but the crosssection of manmade tunnel is about 2.1 m by 2.1 m in general owing to the movement of people through the tunnel. Moreover, information related with the inner state of the tunnel is not provided in advance. Thus, the chance to explore a tunnel may be decreased if the tunnel is partially filled with water. As a detour, another conventional signature which may not be affected severely by collected water in a tunnel needs to be interpreted.
3.2. Time of Arrival
On behalf of the TOP, the time of arrival (TOA) is adopted since the TOA is only related with pulse dominantly propagated through the fastest route inside the tunnel, air [13, 14]. To find the optimum amplitude level of the TOA at first, the Bscan image of a half waterfilled tunnel is considered again. At every depth, time corresponding to the 10%, 1%, and 0.1% of the amplitude of the first positive peak is selected as the TOA and illustrated in Figure 5, respectively. When the amplitude level is 0.1% of the amplitude of the first positive peak, the fast arrival of the TOA is the most accelerated compared with two amplitude levels of 10% and 1%. Despite the described advantage, the low amplitude level may be difficult to detect in practice since the amplitude may be smaller than the detectable range of receiving system. The TOA profile with 0.5% amplitude level of the measured first positive peak utilizing the developed crossborehole pulse radar in our group was seriously contaminated by noise. However, the much clear TOA profiles were obtained when the amplitude level was increased to 1%, 5%, and 10%, respectively [14]. Thus, 10% of the amplitude level of the first positive peak is considered next. The amplitude level seems to be enough to detect using the receiving system but the pattern of the TOA profile approaches that of the TOP profile. By considering both the fast arrival and the detectable range of receiving system, 1% of the amplitude of the first positive peak is selected as the optimum amplitude level. The TOA profiles extracted from 6 different Bscan images with the optimum amplitude level of 1% are illustrated in Figure 6. As expected, the TOA profiles also arrive fast near the depth of 73 m. The fastest arrival of each TOA profile corresponding from to is 238.5329 ns, 238.5849 ns, 238.8966 ns, 239.2603 ns, 239.7799 ns, and 240.5072 ns, respectively.
To quantify fast arrivals, the time advance is introduced by subtracting 6 different TOA profiles from the TOA profile of no tunnel, respectively. Then, the maximum time advance of the TOA profile () is extracted and illustrated in Figure 7 depending on the height ratio of water in a tunnel . By applying the same method to the TOP profile, the maximum time advance of the TOP profile () is chosen and also plotted in Figure 7 for comparison. In Figure 7, the maximum time advance is larger than the even though the height ratio of water inside the tunnel is the same. In spite of the described advantage, the maximum time advance also slows down depending on the increased height ratio of water in the tunnel . The cause of the slowed TOA can be explained as the decreased portion of the air region inside the tunnel. Even though the degree of decrement is small, the same problem with the TOP has also arisen from the TOA.
4. Enhanced Tunnel Signature
4.1. Enhanced Time of Peak
Instead of finding another conventional signature which is robust to the effect of collected water inside an empty tunnel, a method to enhance the weak TOP is investigated. In Figure 7, the delayed maximum time advance owing to the collected water inside the tunnel can be improved like corresponding value of an empty tunnel () if decreased time advance expressed as can be obtained immediately. Then, the enhanced can be calculated easily by adding the decreased time advance to the measured maximum time advance . However, the decreased time advance cannot be determined immediately using measurable parameter of the maximum time advance since both the maximum time advance of the empty tunnel and the height ratio of water in the tunnel are unknown in practice. Hence, it is very important to find the delayed time advance by employing measurable parameters to enhance the weak TOP.
Thus, a conversion formula is derived to calculate the delayed time advance . During the process of formulating, the delayed time advance divided by the maxim time advance , the enhancement rate , is utilized to exclude the effect of the dielectric property of underground rock. The enhancement rate expressed as increases gradually depending on the height ratio of water in the tunnel since the delayed time advance increases depending on the decreased maximum time advance . If the height ratio of water in the tunnel can be changed into a measurable parameter, the enhancement rate can be calculated immediately. Thus, measurable rate is defined as The measurable rate also increases monotonically depending on the height ratio of water in the tunnel since decrease in the maximum time advance is larger than the maximum time advance . Thus, onetoone correspondence between the measurable rate and the enhancement rate is illustrated in Figure 8 instead of the height ratio of water in the tunnel . The relation between the and displayed in Figure 8 is approximated by a curve of the thirdorder polynomials as Finally, the enhanced maximum time advance can be calculated utilizing the developed relation curve in (5) and measurable parameters and as
The possibility of the proposed cubic curve is confirmed by comparing the maximum time advances , , and . To compare the degree of variation relatively, the maximum time advances and are divided by corresponding value of an empty tunnel and and plotted in Figure 9. The enhanced maximum time advance is divided by the value of and also illustrated in Figure 9. Based on the relative time advance of 1.0, the computed average differences corresponding to the maximum time advances , , and are 22.86%, 5.38%, and 0.16%. Different from gradually decreased and , keeps similar relative time advance to the value of an empty tunnel even though the increases.
4.2. Effects of a Wide Variety of Granites on the Enhanced Time of Peak
One may suppose that developed curve model to enhance the maximum time advance may be affected by measurement errors caused during field experiments since the relation curve is introduced from simulation results. However, measurement error caused by receiving system can be ignored since the amplitude of the TOP is large enough compared to the receiver noise [14]. In contrast, measurement error that has arisen from the inhomogeneous underground rock of granite may not be ignored since the dielectric property of granite in Korea suffers from severe variations [15]. It means that developed formula in (5) may be affected dominantly by the various dielectric properties of granite during field experiments in real sites. To cover the effect of inhomogeneous underground rock, simulations are repeated by considering the lower and higher dielectric properties of granite [15] than the medium dielectric properties of granite employed in Section 2. Parameters related with frequency dependent low and high dielectric properties of granite are summarized in Table 1.

The maximum time advances and corresponding to the high and low dielectric properties of granites are extracted from calculated Bscan images and listed in Table 2. Then, the measurable rate and the enhancement rate are calculated using the values in Table 2. By considering the low, medium, and high dielectric properties of granite together, relations between the measurable rate and the enhancement rate are plotted in Figure 10. As a result, a unified relation curve can also be approximated by a cubic relation as

The possibility of the relation curve in (7) is confirmed in Figure 11 by comparing the relative maximum time advances , , and again. The maximum time advances and are divided by corresponding value of an empty tunnel and . And the enhanced maximum time advance is divided by the value of . In Figure 11, small variations are observed at enhanced maximum time advance corresponding to the medium dielectric granite. In the cases of the high and low dielectric properties of granite, relatively larger variation seems to be observed at enhanced maximum time compared to the maximum time when the height ratio of water in tunnel is smaller than 0.3. Thus, average differences between relative time advance of 1.0 and the maximum time advances , , and are computed. The computed average differences corresponding to the maximum time advances , , and are 22.92%, 5.82%, and 2.59% based on the relative time advance of 1.0. It leads us to conclude that developed curve model in (7) can be utilized during field experiments to improve weak due to stagnant water in tunnel.
5. Conclusion
FDTD simulations were performed for 6 different water heights from the bottom to the half height of an empty tunnel to interpret the effects of collected water on conventional signatures used to find an empty tunnel. The fastest arrivals of the TOP and TOA were slowed gradually depending on the increased height of water inside a tunnel. To enhance the weak TOP similar to the TOP of an empty tunnel, a cubic relation was formulated using the fastest TOA divided by the fastest TOP. Then, a unified cubic relation including the effect of two granites with the low and high dielectric properties was derived and the possibility was confirmed. Hence, the developed enhancement process will be helpful for detecting a tunnel in real sites even though the tunnel is partially filled with water.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
This work was supported by the Korea Institute of Science and Technology Institutional Program under Grant 2E24790.
References
 R. J. Lytle, E. F. Laine, D. L. Lager, and D. T. Davis, “Crossborehole electromagnetic probing to locate highcontrast anomalies,” Geophysics, vol. 44, no. 10, pp. 1667–1676, 1979. View at: Publisher Site  Google Scholar
 T.K. Lee, S.O. Park, J.W. Ra, and S.Y. Kim, “Nearfield diffraction pattern by an underground void of circular cylinder,” Microwave and Optical Technology Letters, vol. 2, no. 5, pp. 179–183, 1989. View at: Publisher Site  Google Scholar
 J. Schneider, J. Brew, and I. C. Peden, “Electromagnetic detection of buried dielectric targets,” IEEE Transactions on Geoscience and Remote Sensing, vol. 29, no. 4, pp. 555–562, 1991. View at: Publisher Site  Google Scholar
 S.Y. Kim and J.W. Ra, “The role of cross borehole radar in the discovery of fourth tunnel at Korea DMZ,” in Proceedings of 4th Tunnel Detection Symposium on Subsurface Exploration Technology, pp. 253–257, Golden, Colo, USA, 1993. View at: Google Scholar
 H.K. Choi and J.W. Ra, “Detection and identification of a tunnel by iterative inversion from crossborehole CW measurements,” Microwave and Optical Technology Letters, vol. 21, no. 6, pp. 458–465, 1999. View at: Publisher Site  Google Scholar
 I. C. Peden and J. Brew, “A laboratory scale model for the study of subsurface scattering in lowloss media with applications to ground penetrating radar,” Journal of Applied Geophysics, vol. 33, no. 13, pp. 109–118, 1995. View at: Publisher Site  Google Scholar
 G. R. Olhoeft, “Interpretation of holetohole radar measurements,” in Proceedings of 3rd Technical Symposium on Tunnel Detection, pp. 616–628, Golden, Colo, USA, 1988. View at: Google Scholar
 M. L. Moran and R. J. Greenfield, “Radar signature of a 2.5D tunnel,” Geophysics, vol. 58, no. 11, pp. 1573–1587, 1993. View at: Publisher Site  Google Scholar
 H. Zhou and M. Sato, “Subsurface cavity imaging by crosshole borehole radar measurements,” IEEE Transactions on Geoscience and Remote Sensing, vol. 42, no. 2, pp. 335–341, 2004. View at: Publisher Site  Google Scholar
 S. Ebihara, “Directional borehole radar with dipole antenna array using optical modulators,” IEEE Transactions on Geoscience and Remote Sensing, vol. 42, no. 1, pp. 45–58, 2004. View at: Publisher Site  Google Scholar
 K. Takahashi and M. Sato, “Parametric inversion technique for location of cylindrical Structures by crosshole measurements,” IEEE Transactions on Geoscience and Remote Sensing, vol. 44, no. 11, pp. 3348–3355, 2006. View at: Publisher Site  Google Scholar
 H. Zhou, M. Sato, T. Takenaka, and G. Li, “Reconstruction from antennatransformed radar data using a timedomain reconstruction method,” IEEE Transactions on Geoscience and Remote Sensing, vol. 45, no. 3, pp. 689–696, 2007. View at: Publisher Site  Google Scholar
 S.W. Kim and S.Y. Kim, “Analysis of crossborehole pulse radar signatures measured at various tunnel angles,” Exploration Geophysics, vol. 41, no. 1, pp. 96–101, 2010. View at: Publisher Site  Google Scholar
 S.W. Kim, S.Y. Kim, and S. Nam, “Estimation of the penetration angle of a manmade tunnel using time of arrival measured by shortpulse crossborehole radar,” Geophysics, vol. 75, no. 3, pp. J11–J18, 2010. View at: Publisher Site  Google Scholar
 S.W. Kim and S.Y. Kim, “Estimation of the strike angle of an empty tunnel in different rocks using a crossborehole pulse radar,” Radio Science, vol. 47, no. 1, pp. 1–11, 2012. View at: Publisher Site  Google Scholar
 M. L. Moran and R. J. Greenfield, “Radar waveforms from a three dimensional tunnel,” in Proceedings of the 4th Tunnel Detection Symposium on Subsurface Exploration Technology, pp. 71–84, Golden, Colo, USA, 1993. View at: Google Scholar
 G. A. Ellis and I. C. Peden, “Analysis technique for buried inhomogeneous dielectric objects in the presence of an airearth interface,” IEEE Transactions on Geoscience and Remote Sensing, vol. 33, no. 3, pp. 535–540, 1995. View at: Publisher Site  Google Scholar
 G. Bozza, C. Estatico, M. Pastorino, and A. Randazzo, “Application of an inexactNewton method within the secondorder Born approximation to buried objects,” IEEE Geoscience and Remote Sensing Letters, vol. 4, no. 1, pp. 51–55, 2007. View at: Publisher Site  Google Scholar
 S. Caorsi, A. Massa, M. Pastorino, M. Raffetto, and A. Randazzo, “Detection of buried inhomogeneous elliptic cylinders by a memetic algorithm,” IEEE Transactions on Antennas and Propagation, vol. 51, no. 10, pp. 2878–2884, 2003. View at: Publisher Site  Google Scholar
 J.H. Jung, S.W. Kim, Y.S. Kim, and S.Y. Kim, “Electromagnetic propagation from the intestineingested source in the human body model,” IEEE Transactions on Antennas and Propagation, vol. 58, no. 5, pp. 1683–1688, 2010. View at: Publisher Site  Google Scholar
 K. S. Yee, “Numerical solution of initial boundary value problems involving Maxwell’s equations in isotropic media,” IEEE Transactions on Antennas and Propagation, vol. 14, no. 3, pp. 302–307, 1966. View at: Publisher Site  Google Scholar
 A. Taflove and S. C. Hagness, Computational Electrodynamics: FiniteDifference TimeDomain Method, Artech House, Norwood, NJ, USA, 2000. View at: MathSciNet
 R. Luebbers, F. P. Hunsberger, K. S. Kunz, R. B. Standler, and M. P. Schneider, “A frequencydependent finitedifference timedomain formulation for dispersive materials,” IEEE Transactions on Electromagnetic Compatibility, vol. 32, no. 3, pp. 222–227, 1990. View at: Publisher Site  Google Scholar
 D. M. Hagl, D. Popovic, S. C. Hagness, J. H. Booske, and M. Okoniewski, “Sensing volume of openended coaxial probes for dielectric characterization of breast tissue at microwave frequencies,” IEEE Transactions on Microwave Theory and Techniques, vol. 51, no. 4, pp. 1194–1206, 2003. View at: Publisher Site  Google Scholar
 S.Y. Hyun, Electromagnetic characterization of impulse groundpenetrating radar [Ph. D. dissertation], Department of Radio Sciences and Engineering, Korea University, Seoul, South Korea, 2008.
 J. A. Roden and S. D. Gedney, “Convolutional PML (CPML): an efficient FDTD implementation of the CFSPML for arbitrary media,” Microwave and Optical Technology Letters, vol. 27, pp. 334–339, 2000. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 JiHyun Jung and SeYun Kim. 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.