# Efficient Evaluation of the External Mutual Coupling in Dielectric-Covered Waveguide Slot Arrays

Authors: Zusheng Jin; Giorgio Montisci; Giovanni Andrea Casula; Hu Yang; Junqi Lu
Journal: International Journal of Antennas and Propagation
 (2012)

---

## Abstract

An accurate and efficient procedure is devised to evaluate the mutual coupling in dielectric-covered planar arrays of longitudinal slots. This approach takes full advantage of the discrete complex image method to cast the spatial-domain Green's functions into closed forms, and hence a direct calculation of mutual coupling in the spatial-domain is available. The computation time reduces significantly compared to the previous spectral-domain procedure, without any loss in the accuracy, rendering this approach very attractive for the design of large dielectric-covered planar arrays.


---

## Body

## 1. Introduction

Broadwall waveguide slots are widely employed as radiating elements in arrays for remote sensing, communication, and radar systems up to the millimeter wave frequency band. The advantages of these antennas are the high efficiency, the high gain, a good polarization purity [1], and a considerable mechanical strength.The analysis and synthesis procedures available in the open literature [2–12] allow to achieve a very good accuracy in the design of waveguide slot arrays. Moreover, the recent advances in manufacturing technology assure high precision in the realization of these antennas even in the millimeter wave range [13]. Different kinds of broadwall radiating slots have been proposed in the past literature: longitudinal slots [1], inclined slots for linear [8], or circular polarization [14], U-shaped slots for circular polarization [15]. Among them, longitudinal slots are the most popular solution, since they allow the higher polarization purity than other waveguide slots. Therefore, their use, started in the 40s, is still continuing nowadays.The design of planar arrays of longitudinal slots reached a state of maturity with the work by Elliott and his team, and the proposed iterative synthesis procedure has become a standard for the design of these antennas [1, 2].In many practical applications, waveguide slot arrays need to be covered by a dielectric material for protection, aerodynamic purposes, or pressurization of the radiating waveguides [5, 10, 16]. The extension of the Elliott design procedure for such case has been recently proposed by Casula and Montisci [5] for a small planar array.As it is well known, the external mutual coupling has a significant effect on the performance of a slot array, causing the slots to be detuned at the design frequency. Therefore, an accurate evaluation of the external mutual coupling is required in order to comply with the design specifications on the far-field pattern and on the impedance matching. Moreover, since the iterative synthesis procedure for large arrays is computationally expensive, efficient computation of the external mutual coupling is also required.In [5], the external mutual coupling is evaluated resorting to the spectral-domain by using the Parseval identity and through numerical integration. A specially tailored numerical technique, the so-called Weighted Average Algorithm (WAA) [17], has been devised, and the integration of the highly oscillatory integrand has been accelerated by using the Shank's transform [18]. However, this approach cannot be considered entirely satisfactory, in particular for the designing of large arrays. Actually, the convergence of the integral in the mutual coupling expression deteriorates rapidly as the distance between the slots increases. As a result, its evaluation becomes time-consuming, and it is hard to ensure the accuracy for large distances. Moreover, the WAA requires to extensively investigate the behavior of the integrand, rendering the spectral-domain approach inefficient and not robust. In fact, the oscillation period depends in a complicated manner on the spacing between the slots, and the WAA must be tailored to this behavior.These drawbacks can be avoided if the calculation is performed directly in the spatial-domain, which is the aim of this paper. An alternative spatial-domain approach for the computation of the external mutual coupling has been proposed with the help of the discrete complex image method (DCIM). Numerical examples show that the computational time required for an accurate evaluation of the mutual coupling is only a small fraction of that required by the previous spectral-domain procedure. Moreover, the computational time and the accuracy are independent of the distance between the slots, and no special treatments are needed for large distances, as required by the procedure presented in [5].The spatial-domain approach is described in detail in Section2. Section 3 provides several numerical examples to validate the efficiency and accuracy of this procedure, as compared to the spectral-domain approach of [5].


## 2. Evaluation of the External Mutual Coupling

The external mutual coupling between two arbitrary slots in a dielectric-covered planar array, as depicted in Figure1, is proportional to the reaction integral between the magnetic field of one slot (slot n) and the equivalent current distribution of another (slot m) [2] which can be represented as
(1)gnm=A∫slotmH_d(M_n)⋅M_mdS,
where A is a suitable normalization coefficient; M_m and M_n are the normalized equivalent magnetic currents on the slot m and n, respectively; H_d is the magnetic field in the dielectric cover due to the equivalent magnetic current M_n and can be expressed through the vector potential  F_d:(2)H_d(M_n)=-jωεF_d(M_n)+1jωμ∇∇⋅F_d(M_n),
wherein ω is the angular frequency; ε and μ denote the permittivity and permeability of the dielectric cover, respectively.Two arbitrary slots in a planar array (a) top view and (b) sectional view.
(a)
(b)Substituting (2) into (1), after an integration by parts of the second term, we obtain
(3)gnm=A(-jωεI1,nm+1jωμI2,nm)
with
(4)I1,nm=∫slotmF_d(M_n)⋅M_mdS,(5)I2,nm=-∫slotm∇⋅F_d(M_n)∇⋅M_mdS.Using the well-known assumption of a monomodal electric field distribution on a longitudinal slot [2], we get(6)I1,nm=VnsVms×∫-lm/2lm/2∫-ln/2ln/2G(ρnm)cos⁡(πynln)×cos⁡(πymlm)dyndym,I2,nm=πVnsVmslm×∫-lm/2lm/2∫-ln/2ln/2∇⋅[G(ρnm)cos⁡(πynln)]×sin⁡(πymlm)dyndym,
wherein Vns and Vms are the slot voltages; ln and lm are the slot lengths (Figure 1); G(ρ) is the Green’s function of F_d, where ρnm denotes the distance between the source point and the field point. Once the spatial-domain Green’s function G(ρ) has been obtained, the double integrals in (6) can easily be solved via a Gaussian quadrature.The spatial-domain Green’s functions can be obtained from their spectral-domain counterparts by the following Sommerfeld integral(7)G(ρ)=14π∫-∞∞dkρkρH0(2)(kρρ)G~(kρ),
where H0(2) is the zero-order-second-kind Hankel function, and G~(kρ) is the spectral-domain Green’s function, which can be obtained analytically both in the free space and in planar media [19]. For the free-space case, G~(kρ)=jexp(-jkz0z)/kz0 with kz02=k02-kρ2, and therefore G(ρ)=exp(-jk0r)/(2πr) with r2=ρ2+z2. For the dielectric-covered case, G~(kρ) is given by
(8)G~(kρ)=1+jεr(kz0/kz1)tan⁡(kz1h)DTM
with DTM=-jεrkz0+kz1tan(kz1h), kzi2=ki2-kρ2,i=0,1, where εr is the relative permittivity of the dielectric cover, and h is the dielectric thickness. The constants k0 and k1 are the wave numbers of the free space and of the dielectric medium, respectively. In this case, integrals (6) cannot be solved in closed form, and the direct numerical integration of the Sommerfeld integral is very time-consuming due to the highly oscillating and slow converging nature of its integrand.In order to avoid the time-consuming integration typical of the spectral-domain approach, in this work, the computation of the mutual coupling is performed in the spatial-domain by means of a three-level DCIM [20].The DCIM was first proposed by Fang et al. in 1988 [21] and is under development till now [22–24]. The basic idea of DCIM is to approximate the spectral-domain Green's functions in terms of complex exponentials and then transform them into the spatial-domain analytically with the help of the Sommerfeld identity.Before the implementation of the approximation, the surface wave poles (SWPs) must be extracted fromG~(kρ) and their contribution will be analytically added after the approximation. This extraction is necessary in order to ensure the accuracy when the field point is far from the source point. Since the spectral function is an even function of kρ, the poles occur in complex conjugate pairs and the surface wave term in the spectral-domain can be represented mathematically as
(9)G~sw(kρ)=∑i2kρiResikρ2-kρi2,
where kρi and Resi are the ith SWP and the corresponding residue. Substituting (9) into (7), the SWP contributions are obtained as
(10)Gsw(ρ)=-j2∑iH0(2)(kρiρ)kρiResi.The remaining partG~(kρ)-G~sw(kρ), denoted as G~1(kρ), can be approximated by complex exponentials via a three-level generalized pencil of function (GPOF) method [20, 25]. The so-called three-level GPOF, an extension of the previous two-level version, is capable of approximating the lateral waves by introducing an additional path around the branch point. There may be many different parameterizations of the paths, and the one applied here, first proposed and demonstrated by Alparslan et al. [20, Equation 13], provides very accurate closed-form Green's functions over all ranges of distances.By performing the three-level GPOF,G~1(kρ) can be approximated in terms of complex exponentials. Then, using the well-known Sommerfeld identity, the spatial-domain counterpart of G~1(kρ), denoted as G1(ρ), can be written as
(11)G1(ρ)=12π∑q=13∑p=1PqApqe-jk0rpqrpq,
where rpq=ρ2-Bpq2 is a complex function of ρ. The complex coefficients Apq and the complex exponents Bpq are obtained through the GPOF procedure; Pq denotes the number of exponentials used in the approximation in the qth level.Finally, combining all the approximating exponentials, together with the SWP contributions, the spatial-domain Green’s function can be expressed in closed formG(ρ)=Gsw(ρ)+G1(ρ).Once the spatial-domain Green’s function is available, the computation of the mutual coupling (6) is straightforward. A good accuracy is expected for all slots spacings, since the spatial-domain Green’s functions can be calculated accurately in all ranges via the three-level DCIM [20]. Moreover, the coefficients and exponents in the closed-form expression of the Green’s function need to be generated for only once, leading to a very efficient calculation.


## 3. Results

In all the simulations presented in this section, we consider two longitudinal slots with lengthsln=lm=15.5 mm on a standard WR90 waveguide (22.86 mm × 10.16 mm). Let λg be the guided wavelength of the radiating guide fundamental mode and θ the angle between the waveguide axis and the line connecting the centers of the slots (Figure 1).In order to assess the procedure of this paper, first of all, we test its accuracy against the expression of the mutual coupling obtained by Elliott in  of [2, Equation 30] for an array radiating in free space, that is, for εr=1 (see Figure 1). The operating frequency is 9 GHz and θ=45°. The results are reported in Figures 2 and 3, where we show the mutual coupling coefficient as a function of the normalized distance between the slots dnm/λg. An excellent agreement can be observed.Figure 2
Real part of the mutual coupling coefficientgnm, θ=45°.Figure 3
Imaginary part of the mutual coupling coefficientgnm, θ=45°.Then, the spatial-domain procedure has been compared against the spectral-domain approach, which has been validated for a small array in [5]. A number of simulations have been performed (on an Intel Core I7 CPU at 2.80 GHz with 8 GB RAM) at the operating frequency of 9.4 GHz. The dielectric cover thickness is h=0.165 mm, and its permittivity is εr=2.05.In Figures4 and 5, we show the comparison between the spatial-domain procedure proposed in this paper and the procedure used in [5] for θ=0. The curves have been divided in two figures with different scales in order to better show the behavior of the mutual coupling for large values of dnm. The mutual coupling coefficient computed either with the DCIM-based procedure or with the spectral-domain approach of [5] virtually overlap until dnm=12.5λg. For larger values, the results of the procedure used in [5] are not satisfactory,  according to what we have pointed out in the Introduction.Real part of the mutual coupling coefficientgnm, θ=0.
(a)
(b)Imaginary part of the mutual coupling coefficientgnm, θ=0.
(a)
(b)For small values ofθ (near zero), the results of the comparison are similar to those of Figures 4 and 5, though the convergence of [5, Equation 6]  is achieved for a slightly larger value of the maximum distance between the slots. In Figures 6 and 7, we show the comparison for θ=45°. In this case, a good agreement can be observed until dnm=16λg. For θ  around 90°, the procedure of [5] is more stable. In Figures 8 and 9, the comparison for θ=90° is shown: in this case, the results of the procedure in [5] are accurate enough until dnm around 18λg.Figure 6
Real part of the mutual coupling coefficientgnm, θ=45°. The results for  dnm<10λg are not shown since the curves virtually overlap.Figure 7
Imaginary part of the mutual coupling coefficientgnm, θ=45°. The results for  dnm<10λg are not shown since the curves virtually  overlap.Figure 8
Real part of the mutual coupling coefficientgnm,θ=90°. The results for  dnm<15λg are not shown since the curves virtually overlap.Figure 9
Imaginary part of the mutual coupling coefficientgnm, θ=90°. The results for  dnm<20λg are not shown since the curves virtually overlap.As apparent in Figures4–9, the deterioration of the mutual coupling coefficient computed with the spectral-domain approach of [5] is mainly due to the contribution of its real part.In [5], the mutual coupling coefficient has been expressed in the spectral-domain by means of the Parseval identity and the resulting integral has been evaluated by dividing the integration domain into three intervals (see [5, Section 2]). We have found that the integral in the unbounded region, wherein the integrand is highly oscillatory and convergence is very slow, is always real, and this is the reason for the strong instability of the real part of the mutual coupling coefficient.The proposed spatial-domain procedure has been tested and compared with the spectral-domain approach of [5] also for different dielectric cover permittivity and thickness. The same behavior as for the case investigated above has been observed.We have found that the DCIM-based spatial-domain approach is not only stable but also efficient. Actually, the convergence of the results is assured regardless both of the distance between the slotsdnm and of the angle θ. Here are the parameters of the DCIM for the related simulations presented in this section: T01=200, T02=2, T03=99, P1=6, P2=2, P3=2, and N1=200, N2=100, N3=100. Pq(q=1,2,3) denotes the number of exponentials used in the approximation in the qth level and is defined in (11). According to [20], T01,T02, and T03 are the truncation points for the parametrization of the paths (Cap1, Cap2, and  Cap3); N1, N2, and  N3 are the number of samples of the GPOF approximation in the paths (Cap1, Cap2, and  Cap3), respectively. With this choice, the spatial-domain procedure takes about 0.05 seconds to compute the mutual coupling coefficient, and, as pointed out before, this computational time does not depend on the relative position between the slots.On the other hand, if the mutual coupling is calculated in the spectral-domain by using  [5, Equation 6], the integration parameters strongly depend on dnm and θ. According to [5], we have selected here the parameters of the WAA in order to obtain the convergence of the integrations for each θ  and a maximum value of dnm equal to 12λg. In particular, 15 iterations have been used in the WAA. In this case, the computational time required for the evaluation of the mutual coupling is always about tens of times slower than the spatial-domain procedure. If a dnmlarger than 12λg is required, in order to achieve a good accuracy, the integration parameters in the spectral-domain approach of [5] must be suitably recomputed. However, since the convergence of the integral in the mutual coupling expression deteriorates rapidly as the distance between the slots increases, in this case, its evaluation becomes very time-consuming.To rate the improvements in the computational time of the approach proposed in this paper, we take a generalθ=0 case (Figures 4 and 5) as an example. The results are shown in Table 1. The larger the maximum distance between two arbitrary slots, the larger the computational time required for the convergence of the integration technique of [5]. On the contrary, the approach presented in this paper, based on the DCIM, takes about 0.05 seconds to compute the mutual coupling with the same accuracy of [5], and this time is independent of the distance between the slots. Therefore, the technique proposed in this paper is significantly faster than that of [5] and the improvement increases for larger arrays.Table 1
Comparison between the computational times needed for the evaluation of the mutual coupling,θ=0.
d
n
m
/
λ
g
Procedure in [5]
Procedure of this paper
4
2.1 sec
0.05 sec
6
3.2 sec
0.05 sec
8
4.2 sec
0.05 sec
10
5.3 sec
0.05 sec
12
6.3 sec
0.05 sec


## 4. Conclusion

A spatial-domain approach has been proposed to evaluate the external mutual coupling in dielectric-covered waveguide slot planar arrays. The presented procedure is based on the use of the DCIM for the approximation of the spatial-domain Green’s function. The computational time and the accuracy are independent of the distance between the slots, leading to an accurate and efficient evaluation for all ranges of distances. Numerical examples show that the computational time reduces significantly with respect to the spectral-domain approach without any loss in the accuracy.



---

*Source: 491242-2012-07-05.xml*