Mobile QR Code

1. (Faculty of Radio-Electronic Engineering, Le Quy Don Technical University, Hanoi, Vietnam )

Phased array antenna, Cylindrical arrays, Sparse arrays, Sonar arrays, Amplitude distribution

1. Introduction

Cylindrical arrays have been used in many applications, such as radar, sonar, navigation, and communication, due to the capability of 360$^{\circ}$ coverage using an omnidirectional beam or multiple beams, or a narrow beam that can be steered over 360$^{\circ}$ [1,2]. In fish-finding sonar, the cylindrical sonar array is usually mounted below the vessel hull to generate an umbrella-shaped beam pattern in transmit mode and single beams simultaneously steering over 360$^{\circ}$ in receiver mode [3].

As a cylindrical array, the sparse cylindrical sonar array (SCSA) (also called cylindrical arrays with triangular grid [4]) is advantageous in mitigating the number of elements while maintaining the aperture and image quality [3]. With this advantage, SCSA has been used widely in sonar applications to reduce the complexity of both hardware and operation [1,3].

With a complicated geometry, it is difficult for designers to reduce the side-lobe level (SLL) and half-power beamwidth (HPBW) in the beam pattern of SCSA when the main beam is steered to any desired direction in the azimuth plane. Some studies used an Annealing Algorithm [3,5] to control the SLL and HPBW in the array, which was implemented by adding and removing the non-active elements in the array. With this solution, the computation complexity in determining the beam pattern was increased due to the increased number of calculations. In addition, describing the amplitude distribution mathematically is challenging when controlling SLL by this solution. Another solution for investigating the beam pattern in the azimuth plane used only two consecutive circles [6]. These solutions in [3, 5, 6] generated high SLLs > -16 dB and difficult to control.

This paper proposes a determination solution of the amplitude distribution generating a low SLL by exploiting the mathematical expression of the beam pattern and analyzing simulation results in two cases with isotropic elements and directional elements. Using the proposed solution, the optimal amplitude distribution was determined to control the SLL to be lower than a required value and obtain the narrowest HPBW corresponding to the SLL. The effectiveness of the proposed solution was evaluated by comparing the beam patterns derived from different amplitude distributions and releasing the optimal distribution satisfying the particular requirements.

2. Geometry Model of SCSA

Fig. 1(a) presents a geometry model of a full SCSA. In this figure, $\textit{R}$, $\textit{N}$, and $\textit{P}$ denote the radius of each circle, number of elements in a circle, and number of circles in the full SCSA, respectively. The angle between two consecutive elements and the distance between two adjacent circles are denoted by $\Delta \theta =2\pi /N$and $\textit{h}$, respectively. Therefore, the total number of elements in the full SCSA is $\textit{N}$${\times}$$\textit{P}$. When the beam is steered to the desired direction in the azimuth plane (usually in receiver mode), the angle of the active sector in each circle might be chosen as 60$^{\circ}$, 90$^{\circ}$, or 120$^{\circ}$ [1,2], and the number of active elements in each circle is $\textit{Q}$ ($\textit{Q}$ ${\leq}$ $\textit{N}$) (Fig. 1(b)).

(4)
$\psi _{n}\left(\theta _{0},\varphi _{0}\right)=-kR\left(\cos \left(\left(q-1+\frac{1}{2}\left\lceil \frac{p}{2}-\left\lfloor \frac{p}{2}\right\rfloor \right\rceil \right)\Delta \theta \right)-1\right)\cos \theta _{0}\cos \varphi _{0}\\ -kR\sin \left(\left(q-1+\frac{1}{2}\left\lceil \frac{p}{2}-\left\lfloor \frac{p}{2}\right\rfloor \right\rceil \right)\Delta \theta \right)\sin \theta _{0}\cos \varphi _{0}-k\left(p-1\right)h\sin \varphi _{0}$

where $k=\frac{2\pi }{\lambda }$ is the wavenumber, and $\lambda$ is the wavelength.

The far-field continuous-wave (CW) beam pattern of the array with $\textit{P}$${\times}$$\textit{Q}$ isotropic elements in direction (${\Theta}$, ${φ}$), which is called the array factor (AF), is determined to be [1,4]

(5)
A F ( θ , φ ) = n = 1 P × Q a n exp j k R cos q 1 + 1 2 p 2 p 2 Δ θ 1 cos θ cos φ cos θ 0 cos φ 0 × exp j k R sin q 1 + 1 2 p 2 p 2 Δ θ sin θ cos φ sin θ 0 cos φ 0 × exp j k ( p 1 ) h sin φ sin φ 0

In the azimuth plane, to control the SLL and HPBW in the beam pattern of SCSA, the amplitude distributions on row $\textit{a}$$_{q} were chosen to be a uniform window, Dolph-Chebyshev window, Gaussian window, Hamming window, Taylor window, and Kaiser window. A raised cosine-squared weighting (A cosine squared on a pedestal distribution) was also chosen to control the SLL and HPBW in the beam pattern of SCSA [9,10]. This raised cosine-squared distribution (aperture distribution) was determined to be [10] (9) a_{q}\left(q\right)=\Delta +\left(1-\Delta \right)\cos ^{2}\left(\frac{q-\frac{Q+1}{2}}{Q-1}\right) where {\Delta} denotes the pedestal level, which determines the SLL and HPBW in the beam pattern [10]. When {\Delta} = 1, the aperture distribution becomes a uniform distribution. With {\Delta} = 0.08, the aperture distribution becomes a Hamming distribution [9]. With the selected distributions, the beam pattern of SCSA was analyzed in the azimuth plane by the simulation tools to determine the optimal amplitude distribution satisfying the particular requirements of a lower SLL than a required value and the narrowest HPBW. Based on this solution, the optimal amplitude will be determined as section 4 in two cases with isotropic elements and directional elements. 4. Simulation Results To determine the optimal amplitude distribution satisfying the requirements of the SLL and the narrowest HPBW, this study considered an example of SCSA with 54 elements (\textit{N} = 48) on a circle and 16 circles (\textit{P} = 16). Choosing the active sector as 120^{\circ}, the number of active elements on a circle was 19 (\textit{Q} = 19), and the total number of active elements in the array was 16{\times}19 = 304 elements. The simulation was carried out at frequency \textit{f} = 30 kHz, with the velocity of sound in seawater being \textit{c} = 1500 m/s. The distances between two consecutive circles and between two consecutive elements on a circle were chosen as \textit{h} = 0.5{\times}$${\lambda}$ =2.5 (cm) and $\textit{d}$ = 0.8${\times}$${λ} = 4 (cm), respectively. As a result, the radius of a circle in the SCSA was \textit{R} = 34.46 (cm). The desired steering angles in the azimuth plane and elevation plane were assumed to be \theta _{0}=60^{\circ}and \varphi _{0}=0^{\circ}, respectively. As a ULAA in the vertical direction, a Dolph-Chebyshev with the side-lobe attenuation (SLA) -25 dB was chosen for \textit{a}$$_{p}$ to control the SLL and HPBW. To determine the optimal amplitude distribution providing a lower SLL than a required value and the narrowest HPBW in the horizontal direction, the amplitude distributions were selected; the beam patterns were simulated according to the selected amplitude distributions, and the simulation results were compared. The simulations were implemented for two cases with isotropic elements and directional elements of the form 1+cos(${\theta}$$_{en}). To generate a beam pattern with a lower SLL than -23 dB when isotropic elements are used, the amplitude distributions \textit{a}$$_{q}$ in the azimuth plane were chosen, which included: Dolph-Chebyshev window $\textit{SLA}$ = -45 dB, Gaussian window with the standard deviation ${σ}$ = 2.09, Hamming window, Taylor window with nearly constant-level side-lobes adjacent to the mainlobe $\overline{n}$= 5 and a maximum side-lobe level ASLL = -45.5 dB [8], Kaiser window with ${β}$ = 5, and raised cosine-squared weighting with ${\Delta}$ = 0.16. In addition, a uniform window was also chosen to compare the parameters of the beam patterns. Based on MATLAB software and reference (9), these distributions can be generated, and the beam patterns with isotropic elements and directional elements were also generated to compare the results.

Figs. 2(a) and (b) shows the simulation results of the beam pattern of SCSA in the azimuth plane in the two cases. In these figures, the beam patterns obtained from uniform distribution are denoted by curves with black squares The beam patterns derived from the Dolph-Chebyshev distribution, Gaussian distribution, Hamming distribution, Taylor distribution, Kaiser distribution, and raised cosine-squared distribution (${\Delta}$ = 0.16) are denoted by solid blue curves, dashed red curves, dashdot brown curves, dotted cyan curves, curves with green triangles, and curves with magenta circles, respectively. The parameters of beam pattern, including SLL and HPBW, are depicted more clearly in Table 1.

From Fig. 2(a) and Table 1, a uniform distribution provides the narrowest HPBW but generates an uncontrollably high SLL (-10.20 dB), which is higher than in ULAA (approximately -13 dB [10]). In the investigated amplitude distributions generating an SLL < -23 dB in the case with isotropic elements, the Hamming distribution provided the lowest SLL (-25.06 dB) but generated a wide HPBW (5.44$^{\circ}$) and an uncontrollable SLL. Compared to the amplitude distributions providing an SLL < -23 dB, the raised cosine-squared distribution (${\Delta}$ = 0.16) can provide both an SLL < -23 dB (-23.99 dB) and the narrowest HPBW (5.04$^{\circ}$). The others, including the Dolph-Chebyshev distribution (SLA =-45 dB), Taylor distribution (ASLL = -45.5 dB), Kaiser distribution (${β}$=5), and Gaussian distribution (${σ}$ = 2.09) generated both a higher SLL and a wider HPBW than the above aperture distribution.

Considering the element patterns in the azimuth plane of the form 1+cos(${\theta}$$_{en}), the SLLs decreased significantly compared to when neglecting the beam pattern of each element, as shown in Fig. 2(b) and Table 1. In the amplitude distributions mitigating the SLL, except for the Kaiser distribution ({β} = 5) generating an SLL of -24.89 dB, the others provided SLLs < -25 dB. In contrast to the SLLs, HPBWs increased slightly compared to the case with isotropic elements. By comparing the HPBWs, when element patterns in azimuth plane were chosen in the form 1+cos({\theta}$$_{en}$), the raised cosine-squared distribution (${\Delta}$ = 0.16) also provided the narrowest HPBW (5.24$^{\circ}$) according to SLL < -25 dB (-26.57 dB). Therefore, the aperture distribution (${\Delta}$ = 0.16) is the optimal amplitude distribution for SCSA according to the SLL and HPBW in the two cases with isotropic elements and directional elements.

To investigate the change in SLL and HPBW according to parameter ${\Delta}$ of the aperture distribution, this study considered the beam pattern of SCSA with directional elements of the form 1+cos(${\theta}$$_{en}) when {\Delta} changes from 0 to 0.22 with step 0.02. Table 2 lists the parameters of beam pattern, including SLL and HPBW. A larger parameter {\Delta} indicated a lower HPBW (Table 2). On the other hand, SLL did not increase monotonously according to {\Delta}. When {\Delta} changed from 0 to 0.1, the SLL of the beam pattern decreased gradually and achieved the minimum value (-27.20 dB) at {\Delta} = 0.1. With pedestal levels larger than 0.1, the SLLs rose along with the increases in pedestal levels ({\Delta}). Based on expressions (8) and (9), Fig. 3 presents the beam patterns of SCSA in four cases with {\Delta} = 0.2, {\Delta} = 0.4, {\Delta} = 0.6, and {\Delta} = 0.8 when steering the main beams to (60^{\circ}, 0^{\circ}). In this figure, solid red curve, dashed blue curve, dashdot green curve, and dotted black curve express the beam patterns in cases {\Delta} = 0.2, {\Delta} = 0.4, {\Delta} = 0.6, and {\Delta} = 0.8, respectively. According to Table 2 and Fig. 3, parameter {\Delta} was determined to achieve beam patterns with SLLs decreasing to -27.20 dB and the narrowest HPBW. To achieve SLLs < -27.20 dB, a Gaussian distribution with a large {σ} can be chosen as the optimal amplitude distribution because, except for a raised cosine-squared distribution, the Gaussian distribution provides the narrowest HPBW according to the SLL less than requirement value, as listed in Table 1. For example, with {σ} = 2.5, a beam pattern with SLL -33.24 dB and HPBW 5.93^{\circ} for the above configured SCSA based on (8) were obtained. With the Gaussian distribution, a larger {σ} indicates a lower SLL. In contrast to the SLL, HPBW decreased with increasing {σ}. Therefore, parameter {σ} should be larger than 2.5 to reduce SLL to less than -33.24 dB. With the explicit expression of the beam pattern and the analysis of the simulation results, the parameters of the optimal distributions (the aperture distribution in cases of SLL decreasing to -27.20 dB and Gaussian distribution in the other cases) were also determined when changing the required SLL value and requiring the narrowest HPBW. Therefore, it is more flexible and simpler for designers to choose the optimal amplitude distribution for SCSA using the proposed solution compared to using the annealing algorithm. To evaluate the merits of the proposed solution more clearly, the beam pattern in the azimuth plane of SCSA raised by the amplitude distribution in [11] was compared with that derived by the optimal amplitude distribution from the proposed solution. The comparison was implemented based on expression (8) with the same conditions of SCSA configured as above, and the element patterns of each element in the form 1+cos({\theta}$$_{en}$). For simplicity, the amplitude distribution on the row in this reference was determined by considering the cylindrical array as a simpler circular array, whose individual elements represent the columns of the cylindrical array [11]. To reduce SLL, the amplitude distribution in [11] was chosen as a cosine on a pedestal distribution with a pedestal value equal to 0.4 [9]. With this distribution, the beam pattern in the azimuth plane had an SLL of -15.94 dB and HPBW of 4.43$^{\circ}$. Based on the proposed solution, the optimal amplitude distribution on the row for the above SCSA was the raised cosine-squared distribution (${\Delta}$ = 0.57), which provided a beam pattern in the azimuth plane with SLL = -16.07 dB and HPBW = 4.30$^{\circ}$. Fig 4 presents these beam patterns in the azimuth plane. The amplitude distribution from the solution generated both a lower SLL and narrower HPBW than the amplitude distribution in [11].

The proposed solution can be used for SCSA in sonar systems owing to its flexibility and simplicity. On the other hand, experiments have not been carried out using the proposed solution because of the lack of measurement equipment. This work will be implemented in the future.

Table 1. Parameters of the SCSA Beam Patterns with the Amplitude Distributions.
 Distribution for aq Isotropic elements Directional elements SLL (dB) HPBW(°) SLL (dB) HPBW(°) Uniform -10.20 3.76 -11.57 3.92 Dolph-Chebyshev (SLA = -45 dB) -23.16 5.3 -25.13 5.5 Gaussian (σ = 2.09) -23.16 5.1 -25.54 5.30 Hamming -25.06 5.44 -26.98 5.62 Taylor (ASLL= -45.5 dB) -23.08 5.27 -25.07 5.46 Kaiser (β = 5) -23.09 5.49 -24.89 5.66 Raised cosine-squared distribution (Δ = 0.16) -23.99 5.04 -26.57 5.24
Table 2. Parameters of the SCSA Beam Patterns with the Values of Pedestal Levels (Δ).
 Pedestal level (Δ) Directional elements Pedestal level (Δ) Directional elements SLL (dB) HPBW(°) SLL (dB) HPBW(°) 0 -24.26 6.16 0.12 -27.19 5.42 0.02 -25.11 6.01 0.14 -26.96 5.32 0.04 -25.88 5.87 0.16 -26.57 5.24 0.06 -26.52 5.81 0.18 -26.06 5.16 0.08 -26.98 5.62 0.20 -25.46 5.09 0.10 -27.20 5.52 0.22 -24.82 5.02

5. Conclusion

This paper proposed a new solution for determining the optimal amplitude distribution in SCSA that reduces the SLL to < -23 dB in cases with isotropic elements and reduces the SLL to < -25 dB in cases with directional elements, and generated the narrowest HPBW in each case. Based on the proposed solution, SLL and HPBW in the beam pattern of SCSA have been controlled to ensure the particular requirements for SLL and HPBW in two cases with isotropic elements and directional elements. In addition, by analyzing the simulation results, the distribution parameters were defined when the requirements of SLL and HPBW were changed.

REFERENCES

1
Tinh N. D., Khanh T. D., Feb. 2020, An Investigation on Amplitude Distribution for Controlling Side-lobe Level of Sparse Cylindrical Sonar Arrays, in Proc. of ICGHIT2020, pp. 45-46
2
Josefsson L., Persson P., 2006, Conformal Array Antenna Theory and Design, John Wiley & Sons, New Jersey, USA
3
Kirkebø J. E., Austeng A., Apr. 2008, Sparse Cylindrical Sonar Arrays, IEEE Journal of Oceanic Engineering, Vol. 33, pp. 224-231
4
Tinh N. D., et al. , Oct. 2018, A General Coordinate Formula for Designing Phased Array Antennas in Cylindrical Shape with Triangular Grid, Journal of Science and Technology, Le Quy Don Technical University, No. 193, pp. 64-74
5
Xie P. H., Chen K. S., He Z. S., 2009, Synthesis of Sparse Cylindrical Arrays Using Simulated Annealing Algorithm, Progress In Electromagnetics Research Letters, Vol. 9, pp. 147-156
6
Kirkebø J. E., Austeng A., Holm S., Nov. 2004, Layout-optimized cylindrical sonar arrays, Proc. IEEE OCEANS, Vol. 2, pp. 598-602
7
Rahim T., Davies D. E. N., Feb. 1982, Effect of directional elements on the directional response of circular antenna arrays, in IEE Proc. H - Microwaves, Optics and Antennas, Vol. 129, No. 1, pp. 18-22
8
Stutzman W. L., Thiele G. A., 2013, Antenna Theory and Design, John Wiley & Sons, West Virginia, USA, pp. 446-459
9
Harry L. Van Trees , 2002, Optimum Array Processing: Part IV of Detection, Estimation, and Modulation Theory, John Wiley & Sons, pp. 99-101
10
Thomas A. Milligan , 2005, Modern Antenna Design, John Wiley & Sons, Hoboken, New Jersey, pp. 140-144
11
Jacek Marszal , 2014, Digital signal processing applied to the modernization of Polish Navy sonars, Polish Maritime Research, pp. 65-75

Author

Nguyen Dinh Tinh

Nguyen Dinh Tinh received his B.E. degree in Electronics -Telecommuni-cations, M.E. degree in Radar Navi-gation Engineering from Le Quy Don Technical University, Vietnam, in 2008 and 2012, respectively. He has been working as a lecturer at Le Quy Don Technical University since 2009. His research interests include antenna, signal processing, synthetic aperture sonar, and sonar engineering.

Trinh Dang Khanh

Trinh Dang Khanh received his B.E. degree in Electronics -Telecommuni-cations and Ph.D. degree in Electronic Engineering from Le Quy Don Technical University, Vietnam, in 1983 and 2000, respectively. He has been working as an Associate Professor since 2011 at Le Quy Don Technical University. His research interests include microwave circuit design, signal processing, inertial navigation system, and sonar engineering.