# Circular Antenna Array Synthesis Using Multiverse Optimizer

Authors: Ahmet Emre Taser; Kerim Guney; Erhan Kurt
Journal: International Journal of Antennas and Propagation
 (2020)

---

## Abstract

Antenna array synthesis is one of the most popular topics in the electromagnetic field. Since achieving a desired antenna radiation pattern is a mathematical problem, in the literature, there are various optimization algorithms applied to the synthesis process of different kinds of antenna arrays. In this study, Multiverse Optimizer (MVO) and modified MVO (MMVO) are used to perform circular antenna array (CAA) synthesis. During the exploration, exploitation, and local search phases of calculation, MVO uses three concepts in cosmology; white hole, black hole, and wormhole. Convergence capability of this nature-inspired algorithm is employed for finding optimum amplitude and position values of CAA elements in order to achieve an array pattern with low maximum sidelobe level (MSL) and minimum circumference. The performance of MVO and MMVO was tested on five design examples of pattern synthesis, and the obtained results were compared with ten different algorithms. The simulation results show that MVO and MMVO provide low MSLs with small circumferences.


---

## Body

## 1. Introduction

In recent days, with the developing technology, devices have started to host more wireless components. The application of wireless technology in all areas of data transfer and communication has accelerated the spread of its use. However, increment in use of this technology brought some drawbacks such as interference and electromagnetic pollution. To overcome this issue and use frequency spectrum much more efficiently, RF components are improved year by year. Antennas, which constitute one of the most important elements of wireless technology, also took part in this improvement. As a result, different antenna types appeared such as wire, horn, dish, and array.The antenna arrays consist of a combination of multiple antenna elements and they are widely used in communication. Electromagnetic radiation of each element can be constructively or destructively added. Thus, whole array can radiate towards desired direction while minimizing the interference on undesired angles [1]. Thanks to their useful and flexible structure, antenna arrays can provide high gain, high directivity, and high beam steering capability with less interfering effect. These features of antenna arrays make them popular both in military and civil applications. They are widely used in sonar, radar, navigation systems, cellular systems, and satellite communications.Antenna arrays take their names according to the geometric shape of overall array structure. Array elements can be positioned on a single line, called linear antenna arrays, or distributed over a plane, called planar antenna arrays. Array structures can be in form of circle and ellipse too. Moreover, by creating nested structures from circular or elliptical arrays, concentric circular and concentric elliptical antenna arrays are achieved.In this paper, circular antenna array (CAA) is considered for the array synthesis. Unlike the linear antenna arrays, the main lobe of CAA can be directed between 0 and 360 degrees. According to the changes in position, amplitude, and phase of array elements, or changes in circle radius, overall radiation pattern is shaped. However, achieving a desired antenna pattern is a challenging task since the mentioned array parameters have to be calculated. It is seen from the literature that the optimization algorithms can be used to overcome this kind of pattern synthesis problems. The following can be given as examples of these studies: linear antenna array synthesis using Tabu Search Algorithm (TSA) [2], Harmony Search Algorithm (HSA) [3], and Touring Ant Colony Optimization (TACO) [4]; planar antenna array synthesis using Taguchi Method (TM) [5] and Genetic Algorithm (GA) [6]; elliptical antenna array synthesis using Backtracking Search Optimization Algorithm (BSA) [7]; circular antenna array synthesis using Differential Evolution (DE) Algorithm [8]; and concentric CAA synthesis using BSA [9].Multiverse Optimizer (MVO), which is a nature-inspired algorithm, was proposed in 2016 [10] and used to solve practical optimization problems in engineering such as welded beam design, gear train design, three-bar truss design, pressure vessel design, and cantilever beam design. In [10], the performance of MVO was compared with GA, Particle Swarm Optimization (PSO), Grey Wolf Optimizer (GWO), and Gravitational Search Algorithm (GSA). The majority of test results show that MVO outperforms the popular algorithms in the literature. MVO was also used for multiobjective optimal antenna synthesis of microwave wireless power transmission [11].The main objective of this paper is to compute optimally the amplitudes and positions of CAA elements to achieve low MSL and small circumference using MVO and modified MVO (MMVO). Synthesis process was applied for CAAs with 8, 10, 12, and 20 elements. The obtained results were compared with the results of ten different optimization algorithms: hybrid algorithm (HA) [12] based on Antlion Optimization (ALO) and Grasshopper Optimization Algorithm (GOA), ALO [12], Sequential Quadratic Programming (SQP) [13], Self-Adaptive Differential Evolution (SADE) [13], Biogeography Based Optimization (BBO) [13], GA [14], PSO [15], Firefly Algorithm (FA) [16], Seeker Optimization Algorithm (SOA) [17], and Simulated Annealing (SA) [18].In the following section, MVO and MMVO are explained in detail. Section3 includes the CAA geometry and equation of array factor. Section 4 presents results and discussions. Section 5 covers the conclusion part.


## 2. Multiverse Optimizer

Multiverse theory was proposed by physicists, and MVO is inspired by this theory. It is believed that more than one big bang occurred, and each one of them is resulted with the birth of a universe [19, 20]. The idea behind this theory is the opposite of universe, which states the presence of other universes in addition to ours. Besides, according to this model, these universes may interact and/or collide with each other.MVO is based on three concepts in the cosmology: white holes, black holes, and wormholes [10]. Physicists state that the white hole, which is also considered as big bang, is the main component of the birth of a universe [21]. It appeared after the collision of two parallel universes [22]. Black holes are the white holes in reverse and, because of their extremely high gravitational force, they attract everything including light beams [23, 24]. Wormholes can be described as time/space highways where objects are able to travel instantly between two different points of a universe or two different universes too [25].Basically, MVO converges the target by sending and receiving the objects (variables) between the universes (solutions). This sending/receiving process depends on the inflation rates (fitness values) of the universes calculated over iterations. If a universe has high inflation rate, it is highly probable to have white holes and tends to send object towards other universes. If a universe has low inflation rate, most probably it has a black hole and tends to receive objects. This mentioned process is the exploration step of the algorithm and is represented mathematically as follows [10]:(1)U=x11x12⋯x1px21x22⋯x2p⋮⋮⋮⋮xn1xn2⋯xnp,xij=xkj,r1<NIUi,xij,otherwise,where Ui is the ith universe, p is the number of parameters (variables), n is the number of universes (solutions), xij is the jth parameter of ith universe, r1 is a random number in [0, 1], NI (Ui) is normalized inflation rate of the ith universe, and xkj is the jth parameter of kth universe selected by a roulette wheel selection mechanism.In the exploitation phase of the algorithm, it is assumed that each universe is using wormholes to transfer the objects through space randomly regardless of inflation rate. Besides, wormholes are used to perform local updates in each universe by establishing a connection to the best universe formed so far. This mechanism is formulated as follows [10]:(2)xij=Xj+TDR×ubj−lbj×r4+lbjr3<0.5Xj−TDR×ubj−lbj×r4+lbjr3≥0.5,r2<WEP,xij,r2≥WEP,where Xj is the jth parameter of best universe formed so far; r2, r3, and r4 are the random numbers in [0, 1]; and ubj and lbj are used to define upper and lower bound of jth variable, respectively. WEP and TDR are coefficients which are explained below.WEP stands for the wormhole’s existence probability in the universe. In the algorithm, this coefficient is increased linearly over repetitions to emphasize exploitation. WEP is calculated as follows [10]:(3)WEP=min+l×max−minL,where, by default, min is defined as 0.2 and max is defined as 1. L indicates the maximum iteration value and l represents the current step of the iteration.Traveling distance rate (TDR) is also a coefficient. It is used while sending objects through a wormhole to the best universe obtained so far. The formula for TDR is as follows [10]:(4)TDR=1−l1/pL1/p,where p indicates the exploitation accuracy and is equal to 6 as a default value.The pseudocode of the MVO is shown in Figure1.Figure 1
Pseudocode of MVO.


### 2.1. Modified MVO (MMVO)

As is stated in No Free Lunch (NFL) Theorem, there is no single algorithm to solve every problem set [26]. Besides, it is realized that the performance of an algorithm can be improved by modifying it and thus the algorithm can be used to solve problems in different fields. Even though the performance of original MVO is good enough to optimize the problems, structure of the algorithm itself is suitable for modifications to improve the convergence. Therefore, in this paper, the exploitation phase of the MVO is changed as follows:(5)xij=Xj+TDR×Xj−xij×r4+xijr3<0.5Xj−TDR×Xj−xij×r4+xijr3≥0.5,r2<WEP,xij,r2≥WEP,and this modification aims to perform the solution updates locally using the best universe and current universe instead of upper and lower bounds of problem set. Thus, by shortening the exploitation steps, the algorithm is improved and forced to find better amplitude and position values for array elements.


## 2.1. Modified MVO (MMVO)

As is stated in No Free Lunch (NFL) Theorem, there is no single algorithm to solve every problem set [26]. Besides, it is realized that the performance of an algorithm can be improved by modifying it and thus the algorithm can be used to solve problems in different fields. Even though the performance of original MVO is good enough to optimize the problems, structure of the algorithm itself is suitable for modifications to improve the convergence. Therefore, in this paper, the exploitation phase of the MVO is changed as follows:(5)xij=Xj+TDR×Xj−xij×r4+xijr3<0.5Xj−TDR×Xj−xij×r4+xijr3≥0.5,r2<WEP,xij,r2≥WEP,and this modification aims to perform the solution updates locally using the best universe and current universe instead of upper and lower bounds of problem set. Thus, by shortening the exploitation steps, the algorithm is improved and forced to find better amplitude and position values for array elements.


## 3. Circular Antenna Array

In Figure2, geometry of a nonuniform CAA is depicted. N isotropic radiators (antenna elements) are placed nonuniformly on a circle of radius a in the x-y plane. The array factor of such structure can be mathematically described as follows [1]:(6)AFθ,ϕ=∑n=1NInexpjkasinθcosϕ−ϕn+αn,where(7)ka=2πλa=∑i=1Ndi,ϕn=2π∑i=1ndika.Figure 2
Geometry of a nonuniform CAA.In order to focus the main beam towards the (θo, ϕo) direction, the excitation phase is written as follows [1]:(8)αn=−kasinθocosϕo−ϕn.In the above,In and αn indicate the excitation amplitude and excitation phase of the nth antenna element, respectively. d is used to define the arc distance between the adjacent antenna elements. ϕn is the angular position of the nth radiator on x-y plane. θ represents the measured elevation angle from the positive z-axis, and ϕ is used to represent the measured azimuth angle from the positive x-axis. Since the antenna array is placed in the x-y plane, θ is assumed to be 90o. (θo, ϕo) shows the direction of the main beam of antenna array. In order to simplify the calculations, θo is taken as 90o and ϕo is taken as 0o. In other words, the main beam is focused along the positive x-axis.


## 4. Results and Discussion

In this section, the overall performance of MVO and MMVO on CAA synthesis is evaluated in detail. The purpose of performing array synthesis is to achieve a desired antenna pattern. As mentioned in the Introduction, different synthesis processes were applied to different kinds of antenna arrays with different algorithms. The main aim of those works is to get low sidelobe levels, place nulls in desired directions, and increase the gain or directivity of the antenna.In this paper, MVO and MMVO algorithms were used to find optimum amplitude and position values of CAA elements in order to achieve low MSL and minimum circumference. The obtained results from both MVO and MMVO were comprehensively compared with the results of ten different algorithms which are GA [14], PSO [15], FA [16], SQP [13], SADE [13], HA [12], ALO [12], BBO [13], SOA [17], and SA [18]. During the optimization, the following fitness function was used:(9)Fitness=max20logAFϕMSLAFmax+∑i=1Ndi−C,where ϕMSL is the angular interval of the sidelobe region that is minimized, AFmax is the peak of main beam, and C is the desired circumference. In each iteration step, the difference between the desired circumference and the arc distance summation is taken into account and the algorithm reduces it while keeping the sidelobe level to a minimum. The fitness function is used for CAA with 8, 10, 12, and 20 elements. In order to observe the real performance of MVO and MMVO and to make a fair comparison, 20 independent runs were performed under the default parameters. All simulations were done by setting the algorithm to 300 iterations with population size of 50 on MATLAB R2016b with a computer that has 2.4 GHz i7 processor and 8 GB RAM.In the first example, the optimum amplitude and position values of CAA with 8 elements are calculated by MVO and MMVO to get the minimum sidelobe level in the regionϕ = [−180°, −34°] and ϕ = [34°, 180°]. The resultant radiation patterns of MVO and MMVO are shown in Figure 3. Table 1 illustrates MSL and circumference values of the radiation patterns obtained using MVO and MMVO. Table 1 also shows the comparison results of GA [14], PSO [15], FA [16], SQP [13], SADE [13], HA [12], ALO [12], BBO [13], SOA [17], and SA [18]. It can be clearly observed from Table 1 that MSL and circumference values of MVO and MMVO are better than those of GA [14], PSO [15], FA [16], SQP [13], SADE [13], HA [12], ALO [12], BBO [13], SOA [17], and SA [18]. The convergence curves of this example can be seen in Figure 4.Figure 3
Radiation pattern of the optimized 8-element CAA.Table 1
MSL and circumference comparison of algorithms for 8-, 10-, and 12-element CAA.
Present studyGA [14]PSO [15]FA [16]SQP [13]SADE [13]HA [12]ALO [12]BBO [13]SOA [17]SA [18]MMVOMVON = 8MSL (dB)−15.08−15.04−9.81−10.80−13.00−13.16−12.70−15.00−13.71−12.18−11.11−12.00∑ (d) (λ)4.38104.39364.40944.49314.65564.65214.65054.53844.54164.58144.40705.8750N = 10MSL (dB)−14.71−14.49−11.03−12.31−13.30−13.02−13.43−14.20−13.52−12.72−12.50−13.00∑ (d) (λ)5.86005.86536.08865.90296.07006.03196.05325.88215.96866.05656.10428.0214N = 12MSL (dB)−14.70−14.60−11.80−13.68−14.21−13.41−13.19−14.12−14.56−14.01−14.11−13.91∑ (d) (λ)7.14747.14797.77247.15017.21457.08417.17937.21537.26617.21967.14997.9523Figure 4
Convergence curve of 8-element CAA.In the second example, the number of array elements is increased from 8 to 10 and the sidelobe region is taken asϕ = [−180°, −27°] and ϕ = [27°, 180°]. The algorithm has been run 20 times. The patterns of this design example can be seen in Figure 5. The MSL and circumference results of MVO and MMVO are listed in Table 1 in comparison with GA [14], PSO [15], FA [16], SQP [13], SADE [13], HA [12], ALO [12], BBO [13], SOA [17], and SA [18]. As can be seen from this table, the MSL and circumference values achieved by MVO and MMVO are better than those by the other compared algorithms.Figure 5
Radiation pattern of the optimized 10-element CAA.The 12-element CAA is the third design example of this part. The best MSL and circumference results are obtained after 20 runs in theϕ = [−180°, −23°] and ϕ = [23°, 180°] angular region. The radiation patterns are given in Figure 6. In Table 1, the MSL and circumference values of MVO and MMVO are compared with those of GA [14], PSO [15], FA [16], SQP [13], SADE [13], HA [12], ALO [12], BBO [13], SOA [17], and SA [18]. It is clearly seen from Table 1 that MSL and circumference values of the methods proposed in this study are better than those of the other algorithms except that the circumference value of SQP [13] is slightly better than that of our methods.Figure 6
Radiation pattern of the optimized 12-element CAA.In the fourth example, the number of array elements is set to 20 and the sidelobe suppression region is expanded toϕ = [−180°, −14°] and ϕ = [14°, 180°]. Result comparison of this synthesis example is made in Table 2, and array patterns are shown in Figure 7. According to Table 2, it is evidently seen that the MSL and circumference values obtained by MVO and MMVO are better than those of BBO [13], SADE [13], SQP [13], HA [12], and ALO [12].Table 2
MSL and circumference comparison of algorithms for 20-element CAA.
Varying circumferenceConstant circumferencePresent studyBBO [13]SADE [13]SQP [13]HA [12]ALO [12]Present studyBBO [13]SADE [13]SQP [13]HA [12]ALO [12]UniformMMVOMVOMMVOMVOMSL (dB)−15.40−15.11−13.84−13.95−14.87−14.98−14.39−12.62−12.51−10.70−11.30−11.81−12.35−12.15−6.08∑ (d) (λ)10.602410.657112.397912.185112.325810.742013.89481010101010101010Figure 7
Radiation pattern of the optimized 20-element CAA.In the last example, while suppressing the MSL, the circumference value is fixed to the circumference value of 20-element uniform CAA which is 10λ (di = 0.5 λ). After 20 independent runs, −12.51 dB and −12.62 dB MSL are achieved for MVO and MMVO, respectively. The obtained patterns are shown in Figure 8. In Table 2, the MSL results of patterns produced by MVO and MMVO are compared with the results of BBO [13], SADE [13], SQP [13], HA [12], and ALO [12]. It is apparent from Table 2 that MSL values of MVO and MMVO are better than the MSL values of the other compared algorithms.Figure 8
Radiation pattern of the optimized 20-element CAA with 10λ circumference.The amplitude and position values for all design examples are shown in Table3. Besides, box-and-whisker plots are depicted in Figure 9 to exhibit the overall statistical view of 20 independent runs for each design case. From this figure, improvement in results can be seen as an outcome of modification in the algorithm.Table 3
Amplitude and position values of array elements for design examples.
[d1, d2, d3, …, dN] in λ′s[I1, I2, I3, …, IN]Figure3MVO0.3285, 0.6531, 0.2366, 0.6674, 0.6277, 0.7766, 0.7634, 0.34030.6649, 0.2609, 0.2650, 0.8362, 1, 0.2546, 0.8236, 0.2229MMVO0.3233, 0.7422, 0.8153, 0.6103, 0.6470, 0.2177, 0.7008, 0.32441, 0.3581, 0.8332, 0.7690, 0.1326, 0.3330, 0.9277, 0.1107Figure5MVO0.3012, 0.2732, 0.6967, 0.4202, 0.9751, 0.5266, 0.9490, 0.8571, 0.5480, 0.31820.8811, 0.1196, 0.2380, 0.4808, 0.6397, 0.7434, 0.4745, 0.2438, 1, 0.1999MMVO0.1903, 0.1571, 0.9053, 0.4085, 0.9862, 0.5522, 0.9654, 0.4765, 0.8877, 0.33080.3621, 0.7931, 0.4929, 0.2722, 0.9545, 0.8882, 0.3011, 0.5697, 1, 0.2856Figure6MVO0.3204, 0.7554, 0.7360, 0.6310, 0.9008, 0.4192, 0.9589, 0.1723, 0.7213, 0.3601, 0.9386, 0.23390.9903, 0.4331, 0.6986, 0.4033, 1, 0.9938, 0.7948, 0.1214, 0.2272, 0.5110, 0.9336, 0.5467MMVO0.3159, 0.7274, 0.7454, 0.6109, 0.1259, 0.7763, 0.4936, 0.9420, 0.8175, 0.4464, 0.8833, 0.26280.6145, 0.3664, 0.6626, 0.1591, 0.1449, 0.9212, 1, 0.7296, 0.2799, 0.4902, 0.6086, 0.6184Figure7MVO0.2479, 0.4334, 0.9496, 0.7325, 0.6853, 0.5666, 0.9143, 0.5344, 0.1350, 0.3199, 0.1926, 0.4758, 0.7640, 0.2852, 0.4629, 0.7318, 0.5492, 0.9961, 0.5618, 0.11880.5082, 0.7995, 0.7677, 0.6978, 0.5223, 0.7467, 0.6694, 0.7362, 0.7616, 0.8147, 0.6043, 0.3239, 0.7624, 0.2968, 1, 0.8843, 0.9138, 0.6642, 0.9445, 0.8043MMVO0.3013, 0.4719, 0.8251, 0.1155, 0.6402, 0.8880, 0.3553, 0.9541, 0.5113, 0.1000, 0.2801, 0.4103, 0.8738, 0.2270, 0.6347, 0.7450, 0.6127, 0.9866, 0.4472, 0.22230.7599, 0.2832, 0.1555, 0.6541, 0.3102, 0.5869, 0.9399, 0.3185, 1, 0.1036, 0.5827, 0.7808, 0.1381, 0.4406, 0.6412, 0.7885, 0.5329, 0.5145, 0.9146, 0.4414Figure8MVO0.1008, 0.1130, 0.4527, 0.4436, 0.9314, 0.2575, 0.2372, 0.3995, 0.8275, 0.9908, 0.3814, 1.0813, 0.9903, 0.3910, 0.3395, 0.8914, 0.4817, 0.4714, 0.1111, 0.10690.5859, 0.5800, 0.9237, 0.6956, 0.7435, 0.8154, 0.3572, 0.7263, 0.9135, 0.7137, 0.9515, 0.7478, 0.7961, 0.7528, 1, 0.6600, 0.9821, 0.3981, 0.7731, 0.5785MMVO0.6364, 0.3793, 0.9334, 0.5274, 0.3513, 0.9207, 0.4730, 0.1859, 0.5125, 0.3305, 0.5517, 0.3355, 0.8594, 0.5929, 0.3201, 0.6699, 0.2335, 0.5339, 0.3883, 0.26440.7025, 0.7391, 0.2558, 0.3866, 0.8487, 0.5303, 0.2365, 0.7463, 1, 0.8262, 0.1349, 0.9108, 0.1730, 0.7486, 0.6675, 0.2566, 0.7451, 0.2617, 0.6425, 0.6143Figure 9
Box-and-whisker plots of design examples 1, 2, 3, 4, and 5 for 20 runs.


## 5. Conclusions

In this paper, MVO and MMVO algorithms were presented for the synthesis of CAA. These algorithms were used to calculate the optimum amplitude and position values of the isotropic radiators in the CAA in order to achieve low MSL and minimum circumference. Five different design examples were used to see the overall performance of the MVO and MMVO against ten different optimization algorithms. It was observed that MVO and MMVO have offered minimum sidelobe levels with competitive circumference values. MVO and MMVO can be a good alternative while solving the optimization problems in antenna and electromagnetics.



---

*Source: 3149826-2020-08-24.xml*