Detailed Analysis of Amplitude and Slope Diffraction Coefficients for knife-edge structure in S-UTD-CH Model

Online: 25 March, 2017 In urban, rural and indoor applications, diffraction mechanism is very important to predict the field strength and calculate the coverage accurately. The diffraction mechanism takes place on NLOS (non-line-of-sight) cases like rooftop, vertex, corner, edge and sharp surfaces. S-UTD-CH model computes three type of electromagnetic wave incidence such as direct, reflected and diffracted waves, respectively. As obstacles in diffraction geometry are in the same or closer height, contribution of the diffraction mechanism is dominant. To predict the diffracted fields accurately, amplitude and slope diffraction coefficients and the derivative of these coefficients have to be taken correctly. In this paper, all the derivations about diffraction coefficients are made for knife edge type structures and extensive simulations are performed in order to analyze the amplitude and diffraction coefficients. In plane angle diffraction, contributions of amplitude and slope diffraction coefficient are maxima.


Introduction
Diffraction means bending of electromagnetic waves from a sharp edge, wedge, corner and narrow slits. The diffraction takes place in the case of relatively higher wavelength of electromagnetic wave from obstacle or slit aperture. This paper is an extension of work originally presented in Electromagnetic wave theory symposium [1]. In case that wavelength of electromagnetic wave is less than the height of obstacle; high frequency asymptotic techniques are used [2]. Geometrical optics model runs for refraction, reflection, direct propagation and have been used for a long time in wave propagation in lit region. Geometrical optics model cannot explain diffraction concept [3], and thus cannot give accurate field prediction particularly behind an obstacle. In order to solve the diffraction problem in the field prediction, geometrical theory of diffraction (GTD) model is introduced [4]. However, GTD model gives inaccurate outcomes in optical boundaries [5]. GTD gives inaccurate results in the situation that source, diffraction and observation points are in the same line [6]. GTD model is an extension to geometric optic model by adding diffracted fields terms and used in high frequency spectrum [7]. UTD model is introduced in order to overcome singularity problem of GTD model [8]. UTD model gives appropriate and exact outcomes in case of single diffraction [5]. UTD model is also an extension to geometric optic model and reflected waves are also considered in field prediction [6]. UTD model gives inaccurate outcomes in field prediction in the situation of closer building height. In order to predict the field strength more exactly, the second order diffracted field terms are added to total field at the receiver [9]. Slope diffraction model is introduced in [10,11]. However, there are some singularities in field prediction at shadow boundary points [12]. Higher order diffraction coefficient and diffracted fields are developed to predict the relative path loss at the receiver for wedge and knife edge structures [13][14][15]. Amplitude and slope diffraction coefficients are developed via providing amplitude, slope and phase continuities in optical boundaries [16]. Slope UTD model applied multi-shaped cascaded obstacles [17] and adapted into ASTESJ ISSN: 2415-6698 time domain and all formulas are derived in time domain [18]. A new GTD model is introduced for equal heights case [19]. After 10 diffractions slope UTD model loses the accuracy and requires more computation time. S-UTD-CH model is proposed to predict the field strength accurately after 10 diffractions for knife edge and wedge structures [20,21]. Amplitude and slope diffraction coefficients are used in S-UTD-CH model developed in MATLAB environment. Firstly location and height data of transmitter, obstacle and receiver is entered to program. And then a ray tracing is made and all contribution to total path loss at the receiver is determined. Afterwards, by using amplitude and slope diffraction coefficient path loss is calculated. As an example a simple test scenario is used in order to show that effect of diffraction coefficient is maximum in case of plane angle diffraction. All the simulations are performed in MATLAB.

Amplitude Diffraction Coefficient
To calculate the total field behind an obstacle formula [8] is given by where, Ei is incoming field, D is amplitude diffraction coefficient. ( ) stands for spreading factor. s and k represent travelled distance and wave number, respectively. The amplitude diffraction coefficient for the knife edge type structure is given [12] by where, α is diffraction angle ( α ≡ φ -φ´) and illustrated in Fig. 1. k is wave number and F[x] refers to transition function.

Fig.1. Multiple Knife edge Diffraction Scenario
Amplitude diffraction coefficient used in calculation of the diffracted fields behind the obstacle. The amplitude diffraction coefficient changes with respect to diffraction angle. The diffraction angle is measured between incoming and diffracted fields as depicted in Fig. 1. Changing of the amplitude diffraction coefficient with regard to diffraction angle is shown in Fig. 2.

Fig. 2. Amplitude diffraction coefficient vs diffraction angle
As it is seen in Fig. 2, the diffraction angle variation is indicated on horizontal axis. Variation of the amplitude diffraction coefficient is indicated on vertical axis. Also it can be seen in Fig. 2, the magnitude of amplitude diffraction coefficient has the maximum value at π and -π radian (180̊ and -180). Plane angle diffraction takes place in the case that source, diffracting and observation points are in the same line.

Slope Diffraction Coefficient
In order to calculate the diffracted field behind multiple obstacles, slope diffraction coefficient have to be used to predict the field strength accurately. The slope diffraction coefficient can be calculated by the formula [12] expressed by where, k is wave number, D(α) is amplitude diffraction coefficient and j is complex number. The slope diffraction coefficient obtained by partial derivation of amplitude diffraction coefficient with respect to diffraction angle [12] is given by, where, k is wave number, Ls is distance parameter for slope diffraction coefficient and F[x] refers to transition function.
Slope diffraction coefficient used in calculation of the diffracted fields behind multiple obstacles. The slope diffraction coefficient changes according to the diffraction angle. Again the diffraction angle is measured between incoming and diffracted fields as depicted in Fig. 1. UTD model fails in predicting of the diffracted field in case of equal or closer obstacle height. In order to remove the discontinuities in the predicted field, derivatives of the amplitude diffracted fields have to be added on the receiver. Changing of the slope diffraction coefficient with respect to diffraction angle is shown in Fig. 3. As it is seen in Fig. 3, the diffraction angle variation is indicated on x-axis. Variation of the amplitude diffraction coefficient is indicated on y-axis. Also it can be seen in Fig. 3, the magnitude of the slope diffraction coefficient has the maximum value at π and -π radian (180̊ and -180). By adding the slope diffracted fields, the deficiency of UTD model in equal height case is removed [22][23].
If there are 3 or more obstacles in the propagation scenario, normal derivative of the slope diffraction coefficient have to be used in prediction [12] is given by, where, k is wave number, s is travelling distances, is Ls is distance parameter for slope diffraction coefficient and F[x] refers to transition function. Also F՛[x] is derivative of transition function.
Derivative of the slope diffraction coefficient used in calculation of the diffracted fields behind multiple obstacles, when the obstacle number is greater than 2. The derivative of the slope diffraction coefficient changes with the diffraction angle. Again the diffraction angle is measured between incoming and diffracted fields as depicted in Fig. 1. Changing of the derivative of the slope diffraction coefficient with respect to diffraction angle is shown in Fig. 4.   Fig. 4. Derivative of slope diffraction coefficient vs diffraction angle As it is seen in Fig. 4, the diffraction angle variation is indicated on abscissa. Variation of the amplitude diffraction coefficient is indicated on ordinate. Also it can be seen in Fig. 4, the magnitude of the slope diffraction coefficient has the maximum value at π and -π radian (180̊ and -180).

Results and Discussion
In order to verify the contribution of diffraction is maxima in the plane angle diffraction, test scenario in Fig. 5 is used.

Fig. 5. Plane angle test scenario
As it is depicted in Fig. 5, the transmitter antenna height is 24 m, obstacle is 35 m apart from the transmitter and whose height is 24 m. Operating frequency is 300 MHz. Path loss variation for direct and diffracted fields with respect to the receiver height just behind the obstacle (36 m) is illustrated in Fig. 6   Fig. 6. Path loss of direct and diffracted fields together As can be seen in Fig. 6, diffraction effect is maxima in the plane angle diffraction case. As the receiver height is 24 m, the angle between transmitter and receiver becomes 180. After this height of receiver, direct field is dominated.
Path loss variation for only diffracted fields with respect to the receiver height just behind the obstacle (36 m) is illustrated in Fig. 7. As can be seen in Fig. 7, the diffraction effect is maxima in the plane angle diffraction case. As the receiver height is 24 m, the angle between transmitter and receiver becomes 180.
Path loss variation for reflected and diffracted fields with respect to the receiver height just behind the obstacle (36 m) is illustrated in Fig. 8.

Fig. 8. Path loss of reflected and diffracted fields together
As can be seen in Fig. 8, diffraction effect is very much on the receiver position. In that position ground reflected fields make some contribution to total field. If the receiver height is increased the diffracted field contribution is decreased. Fig. 9. Coverage Map Furthermore, a coverage map is generated by using the simple test scenario in Fig. 5. By using only diffracted field, a coverage map is indicated in Fig. 9.
As it is seen in Fig. 9, the diffraction contribution is maxima in the case of plane angle diffraction. Plane angle diffraction takes place in the case that the receiver height is 24 m.

Conclusions
The diffraction mechanism is significant to predict the field strength and it is used for coverage prediction in urban, rural and indoor environment behind an obstacle like rooftop, edge, corner, and vertex. This field generated by direct, reflected and diffracted waves. The diffracted fields take place in the case of that wavelength of the electromagnetic wave is close or greater than the slit. In order to solve the diffraction problem in shadow region, a lot of electromagnetic wave propagation models are presented, such as GTD, UTD, S-UTD and S-UTD-CH. Amplitude and slope diffraction coefficient directly affect the total electromagnetic fields at the receiver. As the diffraction angle is close to plane angle, the contribution of slope diffracted field is so high. In the case of diffraction angle π radian; in other words observation, diffraction and field points are in the same line, contribution of amplitude and slope diffraction coefficient is maximum. In future, amplitude and slope diffraction coefficient will be developed for rectangular type structure.