# Simple and Exact Closed-Form Expression for Determination of a Penetrated Ray Path in a Ray Tracing

Authors: Hyun Wook Moon; Woojoong Kim; Sewoong Kwon; Jaeheung Kim; Young Joong Yoon
Journal: International Journal of Antennas and Propagation
 (2013)

---

## Abstract

A simple and exact closed-form equation to determine a penetrated ray path in a ray tracing is proposed for an accurate channel prediction in indoor environments. Whereas the penetrated ray path in a conventional ray tracing is treated as a straight line without refraction, the proposed method is able to consider refraction through the wall in the penetrated ray path. Hence, it improves the accuracy in ray tracing simulation. To verify the validation of the proposed method, the simulated results of conventional method, approximate method, and proposed method are compared with the measured results. The comparison shows that the proposed method is in better agreement with the measured results than the conventional method and approximate method, especially in high frequency bands.


---

## Body

## 1. Introduction

In recent years, a lot of indoor wireless applications have been proposed and researched to satisfy various demands such as high data rates, low power consumption [1]. In addition, the operating frequencies of these indoor applications are going up to millimeter-wave wireless personal area networks (WPAN) due to a number of merits and the lack of available frequency bands [1–3]. In order to deploy these applications successfully, exact prediction and detailed understanding of the wave propagation characteristics in indoor environments are necessary. Therefore, there have been many attempts to model indoor radio channel as accurately as possible by using ray tracing [4, 5].The ray tracing technique which is based on geometrical optics (GO) and uniform theory of diffraction (UTD) is capable of calculating an electromagnetic field at a receiving point by adding multipath waves coherently. These multipath waves with different amplitudes and phases depending on wave propagation mechanisms cause deep power fading and/or intersymbol interference (ISI) [6]. Hence, more accurate calculations of amplitude and phase for each multipath wave in ray tracing are very crucial to predict indoor radio channel exactly.In indoor environments, the penetration through the building wall is one of the important wave propagation mechanisms. In conventional ray tracing, the building wall has been treated as a semi-infinite half-space, and then the electric field for the penetrated ray path is calculated with the Fresnel transmission coefficients because finding the exact penetrated ray path in ray tracing is complex and time-consuming [7–9]. This conventional method makes the assumed penetrated ray path (L′) through the wall as an incorrect straight line although the actual penetrated ray path (L) is refracted through a building wall due to difference of refraction indices between the air and the wall and finite thickness of wall as shown Figure 1. Thus, the path error for the penetrated ray path is caused by incorrect two factors: an inaccurate incident angle (θi′), which leads to an error for transmission coefficient, and incorrect total path length from the transmitting point to the receiving point, which leads to a phase error that is proportional to the frequency. As a result, this path error can cause an inaccuracy in calculating the fading effect when the multipath waves are superposed after the building wall in indoor environments.Figure 1
Comparison between an assumed ray path in conventional ray tracing and an actual ray path for penetration through the wall.To overcome this path error within the ray tracing framework, a few studies were carried out. In [10, 11], it was demonstrated that the path error can be reduced when the number of source rays is increased during shooting and bouncing ray tracing (SBR). However, this approach is still not accurate enough and requires increasing simulation time. In [12], an approximate method for the image theory (IT) to determine the penetrated ray path with relatively small error was proposed. Although the error in this method is negligibly small in low frequency bands, it can cause significant phase error in high frequency bands because the phase error is growing proportionally to frequency. Therefore, this method could result in more error in the ray tracing calculation at a millimeter-wave communication environment. Moreover, the approximate method has some complex procedures since it is not given in a closed-form expression.In this paper, a closed-form equation to determine the penetrated ray path through the wall with finite-thickness in the ray tracing is proposed. This proposed equation is not only an exact solution for the penetrating ray path because it is derived from Fermat’s principle [13], but also simple because it is provided in a closed-form without any complex procedure. To validate the proposed method, fading characteristics at 700 MHz and 7 GHz are measured and then compared with simulation results of ray tracing.The paper is organized as follows. The proposed method is described in Section2 and an experimental setup to measure the fading characteristics in indoor environments is described in Section 3. The results in terms of fading effect are presented and discussed in Section 4 and conclusions are given in Section 5.


## 2. A Closed-Form Expression of the Penetrated Ray Path

Figure1 represents an assumed ray path in conventional ray tracing and an actual ray path for penetration through the wall. As mentioned before, since the wall is assumed to be a semi-infinite half-space in conventional ray tracing, the assumed penetrated ray path has two differences compared to the actual penetrated ray path: the incident angle (θi) and the total path length (L). These differences lead to an error in electromagnetic field computations. In order to correct the error and calculate rapidly, a simple and exact closed-form equation to determine the penetrating ray path is derived from Fermat’s principle.Figure2 represents a penetrating ray path through the wall. During ray tracing, the following information is given from geometric parameters: the locations of the transmitting/receiving points, the thickness of the wall, and the electromagnetic properties of air and the wall. Hence, refractive indices n1 and n2 and lengths h1, h2, h3, and d are known parameters. For a given situation, the goal is to determine the parameters d1, d2, and d3.Figure 2
Geometrical representation of penetrating ray path through the wall.According to Fermat’s principle, the derivative of the optical path length from the transmitting point (A) to the receiving point (D) must be zero to have the shortest path [13]. Using the geometry in Figure 2, the optical total path length (l) for the penetrating ray path, whose d2 and d3 are substituted by d1, is expressed as
(1)l=n1AB¯+n2BC¯+n1CD¯=n1h12+d12+n2h22+d22+n1h32+d32=n1h′h12+d12+n2h22+(d-h′d1)2,
where h′=1+h3/h1. To satisfy Fermat’s principle, the derivative of the optical path length (δl/δd1) should be zero to have the shortest path. Thus, from this condition, a fourth-order equation which agrees with Snell's law is obtained, as shown in the following equation:
(2)n1sinθi=n1d1h12+d12=n2(d-h′d1)h22+(d-h′d1)2.Then, four possible solutions ford1 from (2) can be expressed, as shown in the following equation:
(3)d11=0.5(dh′-2X1-X2+X3-X1+X2),d12=0.5(dh′+2X1-X2+X3-X1+X2),d13=0.5(dh′-2X1-X2+X3+X1+X2),d14=0.5(dh′+2X1-X2+X3+X1+X2),
where
(4)X1=1h′2(d2-2Y13Y2),X2=13Y2h′2(Y12(Y3+Y42)-1/3+(Y3+Y42)1/3),X3=(2d/h′)(Y1/Y2h′2-d2/h′2+2h12/Y2)X1+X2,Y1=(n122-1)d2-h12h′2+h22n122,Y2=n122-1,Y3=2Y13+108Y2d2h12h22h′2n122,Y4=Y32-4Y16,n12=n1n2.Among the above solutions,d12 is a unique solution because d1 must have a real value in the range [0,d]. When the value of d1 is fixed, d2 and d3 can be obtained using the relation between them by Snell’s law, and the penetrating ray path is finally calculated by the ray tracing. This result is not only an exact solution for the penetrating ray path because it is derived from Fermat’s principle, but also simple because it is provided in closed-form without any complex procedure.


## 3. Experimental Setup

To validate the proposed method, an experiment was performed at an empty office room on the sixth floor of the engineering building at Yonsei University. As shown in Figure3, the size of the office room was 8.00 m × 6.02 m. It was composed of concrete walls with 20 cm thickness, a window, and a metal door.Figure 3
Experimental site.To generate the fading including the penetrating ray path, the transmitter (Tx) was placed in the middle of a corridor outside the office room and the receiver (Rx) was placed in the office room and moved along the paths (1–8), which are 5 m in length and equally spaced at 50 cm intervals. The metal door was open resulting in reflection. Hence, the reflecting ray from metal door and various multi-path waves after penetration through the wall were superposed and caused the fading.Figure4 shows the measurement system to measure the received signal strength. At the transmitter, the continuous wave (CW) signal was generated from a signal generator (Anritsu MG3695A) and was radiated with 30 dBm transmitting power through a 30 dB gain amplifier and a reference dipole antenna at 1 m height. Then at the receiver, the signal was received through an identical reference dipole antenna at 1 m height and measured by a spectrum analyzer (Agilent 8566B) at the center frequency. This measured received signal strength was stored in a notebook PC through a General Purpose Interface Bus (GPIB) connection. During the measurement, it is very important that the received signal is measured at an exact point in the moving path of the receiver. Hence, four rectangular conductors were mounted on the wheel with equal spacing and a magnetic sensor was also mounted beside the wheel. This magnetic sensor generated trigger signal when it detected the conductors. Thus, the trigger signal which was generated with a regular time interval made the spectrum analyzer operate at the exact point. Since the diameter of the wheel was 7 cm, the spacing between the conductors was about 5.5 cm and total number of samples in each path was 91. In order to obtain the fading characteristics, the measurements at each path were repeated ten times and final received signal strength was calculated by their average value.Measurement system.
(a)
Transmitter
(b)
Receiver


## 4. Fading Effect Comparison

Figure5 shows measured received signal strength with unit dBm at 700 MHz and 7 GHz. The fading of the multi-path waves is shown and the fading effect is more severe and frequent in higher frequency bands because the phase for each multi-path wave changes rapidly in proportion to the frequency. In other words, the prediction of the exact phase for each multi-path wave is very important, especially in high frequency bands.Measured received signal strength (fading effect).
(a)
700 MHz
(b)
7 GHzFigures6 and 7 represent comparison of fading effect along path 2 for measured and simulated results at 700 MHz and 7 GHz, respectively. The reason that path 2 is selected for comparison is that the fading is the most severe on that path. The penetrated ray path in the simulated results using ray tracing is calculated by three methods: conventional method, approximate method, and proposed method. With the conventional method, the penetrated ray path is assumed as a straight line without refraction as shown in Figure 1. Hence, this method has the most inaccurate results. The approximate method in [12] which has some complex procedure and a more error in relation to the proposed method is compared together. At low frequency bands (700 MHz), the simulated results with three methods are matched well with measured results in terms of fading as shown in Figure 6. However, the degree of mismatch is becoming greater as the frequency increases because the inaccurate prediction of multipath waves including the penetrating ray path at high frequency bands due to rapid change of phase leads to an incorrect calculation when the multipath waves are added coherently. To investigate the degree of mismatch according to the methods for calculating the penetrating ray path, the root-mean-square (RMS) errors for three methods between measured results and simulated results at both frequencies are illustrated in Figure 8. With comparison between the proposed method and the conventional method, the proposed method reduces the mean RMS error by 1.38 dB at 700 MHz and 4.51 dB at 7 GHz, respectively. This means that the conventional assumption for the penetrating ray path which is treated as a straight line without refraction is only valid at low frequency bands. On the other hand, the approximate method has more RMS errors of 0.11 dB at 700 MHz and 0.89 dB at 7 GHz compared with the proposed method although it is also able to improve the accuracy compared to the conventional method. Therefore the proposed method is simple to calculate because it is provided in closed-form and also can improve the accuracy of ray tracing by modeling the penetrating ray path exactly, especially in high frequency bands.Figure 6
Comparison of fading effect along path 2 for measured and simulated results at 700 MHz.Figure 7
Comparison of fading effect along path 2 for measured and simulated results at 7 GHz.Comparison of RMS errors between measured results and simulated results.
(a)
700 MHz
(b)
7 GHz


## 5. Conclusion

In this paper, a new method to calculate the penetrating ray path with refraction in the wall exactly and simply is proposed. For this objective, a closed-form equation is derived from Fermat’s principle; hence, it has a simple form and can consider refraction when the ray penetrates through a wall. To validate the proposed method, the measured results and simulated results in terms of fading effect are compared. As a results, the significant improvement in accuracy for the proposed method is verified, especially in high frequency bands. Therefore, the proposed method is very helpful for channel analysis of indoor wireless communication using high frequency bands such as millimeter-wave WPAN application.



---

*Source: 545874-2013-12-11.xml*