Abdelmadjid Recioui

Hamid Bentarzi

Fatma Zohra

Mohamed Amine Ouamri
email: ouamrimouhamedamine@gmail.com

Marius-Emil Oteşteanu
email: marius.otesteanu@upt.ro

Gordana Barb
email: gordana.barb@student.upt.ro

Cedric Gueguen
email: cedric.gueguen@irisa.fr

Coverage Analysis and Efficient Placement of Drone-BSs in 5G Networks †

Keywords: drone-Base Station, coverage, 5G networks, Grey Wolf Optimizer

The integration of drones as base stations has shown to be a potential approach for the future mobile communication systems. Hence, this emerging technology is currently being investigated within the 3GPP standardization community with the main objective of improving coverage and capacity in dense urban areas. Nevertheless, in order to provide adequate coverage for users, it is necessary to find the optimal location of the Drone-BS. This work proposes a novel approach for the Drone-BS in 5G communication systems, using the meta-heuristic algorithm. Firstly, we analyse the downlink coverage probability according to SINR by using stochastic geometry. Afterwards, we apply the Grey Wolf Optimizer algorithm in order to find the optimal Drone-BS placement under coverage probability constraint.

Introduction

The chronological progress and exponential growth of cellular networks has led the mobile communications community to suggest the use of drones as a Base Station (BS), with the incorporation of a transceiver, due to its massive potential in mobile communication systems [1]. This approach provides significant improved coverage, enhanced quality of service and increased capacity [2,[START_REF] Alzenad | 3-D Placement of an Unmanned Aerial Vehicle Base Station for Maximum Coverage of Users with Different QoS Requirements[END_REF]. Furthermore, the high position of the drone-BS encourages Line-of-Sight (LOS) communication to grounds users [START_REF] Mozaffari | Efficient Deployment of Multiple Unmanned Aerial Vehicles for Optimal Wireless Coverage[END_REF], and consequently the optimal drone-BS deployment is paramount for the proper functioning of the network. On the other hand, with the significant increase in traffic volume, 5G mobile networks will exploit the vast amount of available spectrum in the millimeter wave (mmWave) band [START_REF] Chen | The requirements, challenges, and technologies for 5G of terrestrial mobile telecommunication[END_REF]. Simultaneously, a drone equipped with a mmWave base station offers considerable throughput. Several research projects are paying more attention to the drone-BS deployment. For instance, authors in [START_REF] Li | Drone-assisted public safety wireless broadband network[END_REF] study the optimal Unmanned Aerial Vehicles (UAV) placement in order to analyse coverage with fixed UAV altitude. Similarly, in [START_REF] Kalantari | Backhaul-aware robust 3D drone placement in 5G+ wireless networks[END_REF], the authors investigate the optimal placement of drone-cells by taking into consideration the backhaul requirements. In addition, the authors propose the robustness solution of drone-BS placement under the user's movement. A heuristic approach based on Particle Swarm Optimization (PSO) is introduced in order to find the optimal 3D location of drone-cells in [START_REF] Kalantari | On the Number and 3D Placement of Drone Base Stations in Wireless Cellular Networks[END_REF]. The study aims to maximize coverage and limit interference while reducing the drone's altitudes. More recently, in [START_REF] Alzenad | 3-D Placement of an Unmanned Aerial Vehicle Base Station for Maximum Coverage of Users with Different QoS Requirements[END_REF], the authors propose a novel 3D UAV-BS placement using exhaustive search (ES) and maximal weighted area algorithms. They maximize coverage under the Eng. Proc. 2022, 14, 18 2 of 8 quality-of-service requirement constraint. The efficient UAVs deployment in order to provide coverage is suggested in [START_REF] Mozaffari | Efficient Deployment of Multiple Unmanned Aerial Vehicles for Optimal Wireless Coverage[END_REF]. The authors study the downlink coverage probability as function of altitude and antenna gain and apply circle packing theory in order to maximize the coverage area. Authors in [START_REF] Bor-Yaliniz | Efficient 3-D placement of an aerial base station in next generation cellular networks[END_REF] present UAVs-BS placement based on the location of users. In this paper we investigate the drone-BS placement by assuming random distribution and mmWave Path Loss model. We analyse the downlink coverage probability based on altitude and antenna gain. We apply a Grey Wolf Optimizer (GWO) algorithm [START_REF] Mirjalili | [END_REF] in order to find the optimal 3D position that satisfies coverage under interference constraint. The rest of this paper is organized as follows: Section 2 presents the system model adopted for our work, Section 3 describes the downlink coverage probability and Section 4 provides the proposed optimal drone-BSs placement algorithm using GWO. Simulation and numerical results are given in Section 5. Finally, conclusion and future work are discussed in Section 6.

System Model

In this section we present the system model adopted for our work. We consider an urban scenario where N users are distributed randomly in the area represented by (xj,j). The M drones-BSs are randomly distributed with fixed power transmission. Each drone is defined by its altitude h and location (xi,i). Geometrically, the coverage area of the drone is given by Cdrone = (R2 -h2), where R represents the radius of the transmitter drone [START_REF] Al-Turjman | Enhanced Deployment Strategy for the 5G Drone-BS Using Artificial Intelligence[END_REF]. Figure 1 depicts a possible scenario of drone-BSs deployment. Additionally, we assume that the propagation environments are both LOS and Non-Line-of-Sight (NLOS). On the other hand, the metric for satisfying quality of service is based on the determination of the Signal-to-Interference Plus Noise Ratio (SINR), which should be above a certain threshold SINR > SINRthr. According to [START_REF] Alzenad | 3-D Placement of an Unmanned Aerial Vehicle Base Station for Maximum Coverage of Users with Different QoS Requirements[END_REF], the probability of the mean path loss is given by: PL h, d ji = PL LOS P LOS + PL NLOS P NLOS (1) where P LOS represents the probability of the LOS link given by Equation (2) and P NLOS = 1 -P LOS . As discussed in [START_REF] Al-Hourani | Optimal LAP Altitude for Maximum Coverage[END_REF], the path loss between drone-BSs and users for the LOS and NLOS propagation in dB is given by Equation (3).

P LOS = 1 1 + aexp -b 180 π tan -1 h (xj-xi) 2 +(y j -y i ) 2 -a (2)
PL LOS = 20log 4π f c d ji c + η LOS PL NLOS = 20log 4π f c d ji c + η NLOS (3)
where f c represents the carrier frequency, d ji =

x jx i 2 + y jy i 2 + h 2 is the distance between the drone-BSs and user i. Furthermore, η LOS and η NLOS are the average additional losses for LOS and NLOS, respectively. The probability mean path loss can be further rewritten as:

PL h, d jt = (η LOS -η NLOS )P LOS + 10 log h 2 + x j -x i 2 + y j -y i 2 + B (4) with B = 20 log 4π f c c + η NLOS .
Eng. Proc. 2022, 14, 18

3 of 8 𝜔 = 1, 𝑖𝑓 𝑡ℎ𝑒 𝑑𝑟𝑜𝑛𝑒 -𝐵𝑆 𝑖𝑠 𝑑𝑒𝑝𝑙𝑜𝑦𝑒𝑑 𝑎𝑡 (𝑥 , 𝑦 , ℎ) 0, 𝑜𝑡ℎ𝑒𝑟𝑤𝑖𝑠𝑒 (5)
Then, for the second case (user's coverage):

𝜑 = 1, 𝑖𝑓 𝑢𝑠𝑒𝑟 𝑗 𝑖𝑠 𝑐𝑜𝑣𝑒𝑟𝑒𝑑 𝑏𝑦 𝑜𝑛𝑒 𝑑𝑟𝑜𝑛𝑒 -𝐵𝑆 𝑖 0, 𝑜𝑡ℎ𝑒𝑟𝑤𝑖𝑠𝑒 (6)

Downlink Coverage Probability Calculation

In this section we study the downlink coverage probability. When evaluating channel performance, a typical user in the area is associated with a drone-BS when the SINR is greater than a predefined threshold. According to [START_REF] Alzenad | 3-D Placement of an Unmanned Aerial Vehicle Base Station for Maximum Coverage of Users with Different QoS Requirements[END_REF], the received power at the user 𝑗 from its drone-BS is given by:

𝑝 = 𝑝 -𝑃𝐿 ℎ, 𝑑 -𝜗 (7)
where 𝒑 𝒕 is the power transmission and 𝜗 is the shadow fading with Gaussian random variable in dB for LOS and NLOS environments. The SINR ratio is given by the following formula:

𝑆𝐼𝑁𝑅 = 𝑝 ∑ 𝑝 + 𝑁 ∈ , (8)
ℙ 𝑆𝐼𝑁𝑅 > 𝑆𝐼𝑁𝑅 = ℙ 𝑝 ∑ 𝑝 + 𝑁 ∈ , ≥ 𝑆𝐼𝑁𝑅 (9)
Lemma 1. In the drone-BS network, the downlink coverage probability for a ground user located at a distance 𝒅 𝒋𝒊 from drone-BS can be expressed as:

𝒑 𝒄 = 𝑸 𝑷𝑳 𝒉, 𝒅 𝒋𝒊 + 𝑺𝑰𝑵𝑹 𝒕𝒉𝒓 ∑ 𝒑 𝒓 𝒋𝒊 + 𝑵 𝒊 ∈𝑫,𝒊 𝒊 -𝒑 𝒕 + 𝝁 𝝈 (10)
where 𝑄(. ) is the 𝑄 function. As shown in Figure 1, the blue circles represent the drone's coverage area, and the black points indicate the users in the area. At this point, there are two major binary decisions to be conducted. The first one is to determine the deployment of drone-BSs and the second one is to determine the users covered by drone-BSs. The procedure can be developed mathematically as follows: let D = {1 . . . M} represent the set of the candidate drone-BSs and U = {1 . . . N} the set of the users. For the first case, the binary decision can be formulated as:

ω i xyh = 1, i f the drone -BS is deployed at (x i , y i , h) 0, otherwise (5)
Then, for the second case (user's coverage):

ϕ i j = 1, i f user j is covered by one drone -BS i 0, otherwise (6)

Downlink Coverage Probability Calculation

In this section we study the downlink coverage probability. When evaluating channel performance, a typical user in the area is associated with a drone-BS when the SINR is greater than a predefined threshold. According to [START_REF] Alzenad | 3-D Placement of an Unmanned Aerial Vehicle Base Station for Maximum Coverage of Users with Different QoS Requirements[END_REF], the received power at the user j from its drone-BS is given by:

p ji r = p t -PL h, d ji -ϑ (7)
where p t is the power transmission and ϑ is the shadow fading with Gaussian random variable in dB for LOS and NLOS environments. The SINR ratio is given by the following formula:

SI NR = p ji r ∑ i ∈D, i =i p i j r + N (8)
P[SINR > SI NR thr ] = P   p ji r ∑ i ∈D, i =i p ji r + N ≥ SI NR thr   (9)
Lemma 1. In the drone-BS network, the downlink coverage probability for a ground user located at a distance d ji from drone-BS can be expressed as:

p c = Q    PL h, d ji +SINR thr ∑ i ∈D, i =i p ji r + N -p t + µ σ    (10)
where Q(.) is the Q function.

Eng. Proc. 2022, 14, 18 4 of 8

Proof.

p c = P p ji r ∑ i ∈D, i =i p ji r +N ≥ SINR thr = P p ji r ≥ SINR thr ∑ i ∈D, i =i p ji r + N = P p t -PL h, d ji -ϑ ≥ SINR thr ∑ i ∈D, i =i p ji r +N = P ϑ ≤ p t -PL h, d ji -SINR thr ∑ i ∈D, i =i p ji r + N Q   PL(h,d ji )+SINR thr ∑ i ∈D, i =i p ji r +N -p t +µ σ   (11)
where ϑ ∼ N µ LOS , σ 2 LOS and ϑ ∼ N µ NLOS , σ 2 NLOS represent the shadow fading with normal distribution, in dB, for LOS and NLOS cases, respectively. However, the Complementary Cumulative Distribution Function (CCDF) of a Gaussian random variable as show in Equation ( 12) below gives the downlink coverage probability for LOS environment.

p c = Q    PL LOS P LOS + SI NR thr ∑ i ∈D, i =i p ji r + N -p t + µ LOS σ LOS    (12)

Optimal Drone-BSs Placement

In this section we propose to solve the 3D deployment of drone-BSs using GWO [START_REF] Mirjalili | [END_REF], where the goal is to maximize coverage by increasing the number of users covered by drone-BSs. Generally, the drone-BSs placement constitutes a NP-hard problem. For this reason, evolutionary algorithms can be a suitable alternative. Mathematically, the coverage problem can be formulated using Mixed Integer Non-Linear Programming (MINLP), as follows:

max x i , y i , h, ϕ i j N ∑ j=1 M ∑ i=1 ϕ i j (13) Subject to: min M ∑ i=1 ω i xyh (14) h min ≤ h ≤ h max (15) N ∑ j=1 b j ϕ i j ≤ BW (16)
The first constraint assumes that the coverage must be met with a minimum number of drone-BSs deployed. Indeed, a considerable deployment of drones leads to a reduction in the distance between them and hence, the interference from the adjacent drone-BSs increases [START_REF] Mozaffari | Efficient Deployment of Multiple Unmanned Aerial Vehicles for Optimal Wireless Coverage[END_REF]. The variables h min , h max in the second constraint represent the minimum and maximum drone-BS's altitude, respectively. Nevertheless, as the altitude of drone-BS increases, the coverage radius also increases, and energy performance of drone-BSs decreases. The last proposed constraint indicates the total limiting available bandwidth, where bj is the bandwidth required by user j and BW denotes the total bandwidth of the drone-BS [START_REF] Kalantari | Backhaul-aware robust 3D drone placement in 5G+ wireless networks[END_REF]. Now we apply GWO in order to find the optimal 3D drone-BSs placement. Proposed by [START_REF] Mirjalili | [END_REF], GWO is a meta-heuristic algorithm inspired by the Grey Wolf hunting mechanism. The GWO algorithm takes into consideration the search, encircling and attacking prey by Grey Wolf. Notwithstanding, the algorithm classifies wolf into four groups to find leadership: α is the most optimal solution and β, are the second and third best solutions, respectively. X represents the alternative solutions. The process is mathematically modelled as follows:

X(t + 1) = 1 3 (X 1 (t) + X 2 (t) + X 3 (t)) (17)
where:

X 1 (t) = X α -A 1 × D α (18) X 2 (t) = X β -A 2 × D β (19) X 3 (t) = X θ -A 3 × D θ (20)
From the equations above, X 1 (t), X 2 (t), X 3 (t) are the position vectors of the prey, and X(t + 1) represents the grey wolf position vector. Grey wolves encircle prey during the hunt. The mathematically model-encircling behaviour is given by the following formula:

D = C × X p (t) -X(t) (21)
where: C = 2 × r 2 and A = 2a × r 2a are two coefficient vectors, and r 1 , r 2 are vectors randomly generated between [0, 1] . The Algorithm 1 of GWO is given by: Algorithm 1: GWO for drone-BS deployment.

Initializing the drone-BSs population X i = (i = 1, 2, . . . n) Initialize α, A and C for i = 1 : n do j = 1 : m do Calculate a distance, pathloss and evaluate coverage probability in each user end for Find the fitness of each search agent X α , X β , X θ while t < Max iter for each search agent X(t + 1) = 1 3 (X 1 (t) + X 2 (t) + X 3 (t)) end for Update α, A and C Calculate the fitness of all search agents Update X α , X β , X θ t = t + 1 end while return X α

Simulation Results

The simulation and numerical results for the downlink coverage probability and optimal drone-BSs placement are presented in this section. 200 users and 10 drone-BSs are randomly distributed within a 2 × 2 km 2 area. The simulation parameters are illustrated in Table 1. Our goal is to serve the maximum number of users based on a reduced number of drone-BSs using a mmWave carrier frequency. Figure 2 depicts the downlink coverage probability considering different altitudes for the drone-BSs. We observe that the coverage probability increases when the altitude of the drone-BSs also increases. However, as mentioned above, a very high altitude can influence the energy performance of the drone. Figure 3 shows the optimal drone-BSs placement and the minimum required number of drone-BSs to ensure that all users are covered. The results indicate that in order to satisfy 0.76 of downlink coverage probability, it is necessary to deploy five drone-BSs. On the other hand, a coverage probability can be improved by increasing the number of drones, although this would lead to a decrease in the distance between the drones and consequently the interference at user level would be higher.

𝑀 (number of drone-BSs) 10 𝑁 (number of users) 200 𝑝 30 dBm 𝑆𝐼𝑁𝑅 5 dBm 𝐵 20 MHz Figure 2 depicts the downlink coverage probability considering different altitudes for the drone-BSs. We observe that the coverage probability increases when the altitude of the drone-BSs also increases. However, as mentioned above, a very high altitude can influence the energy performance of the drone. Figure 3 shows the optimal drone-BSs placement and the minimum required number of drone-BSs to ensure that all users are covered. The results indicate that in order to satisfy 0.76 of downlink coverage probability, it is necessary to deploy five drone-BSs. On the other hand, a coverage probability can be improved by increasing the number of drones, although this would lead to a decrease in the distance between the drones and consequently the interference at user level would be higher. Figure 2 depicts the downlink coverage probability considering different altitudes for the drone-BSs. We observe that the coverage probability increases when the altitude of the drone-BSs also increases. However, as mentioned above, a very high altitude can influence the energy performance of the drone. Figure 3 shows the optimal drone-BSs placement and the minimum required number of drone-BSs to ensure that all users are covered. The results indicate that in order to satisfy 0.76 of downlink coverage probability, it is necessary to deploy five drone-BSs. On the other hand, a coverage probability can be improved by increasing the number of drones, although this would lead to a decrease in the distance between the drones and consequently the interference at user level would be higher. Furthermore, the number of users served by the deployed drone-BSs are illustrated in Figure 4. We observe a decrease in the number of users covered by the network when the number of drone-BSs is reduced. Moreover, the results obtained indicate that the drone-BSs are separated in distance, leaving uncovered areas. Besides, the optimal altitude selected for the drone-BSs can be improved to increase the coverage radius and get more users covered, but with reduced energy performance.

Furthermore, the number of users served by the deployed drone-BSs are illustrated in Figure 4. We observe a decrease in the number of users covered by the network when the number of drone-BSs is reduced. Moreover, the results obtained indicate that the drone-BSs are separated in distance, leaving uncovered areas. Besides, the optimal altitude selected for the drone-BSs can be improved to increase the coverage radius and get more users covered, but with reduced energy performance.

Conclusions

Drone-BSs deployment in 5G networks is currently an important subject with extreme potential in the mobile communication networks industry. Placement optimization is a complicated task that requires automation algorithms to find the optimal solution for the placement of drone-BSs. In this work, we studied a 3D deployment of drone-BSs in a 5G network using GWO algorithm. Firstly, we analysed the downlink coverage probability according to SINR and altitude. Next, we proposed a mathematical model in order to maximize coverage under a minimum number of drone-BSs, altitude and total limiting available bandwidth. The results obtained showed a minimization of the total number of drone-BSs deployed with a coverage reaching 82.5%. Our study does not take into consideration the blocking effect and overlap between drones that are important for the handoff, but these factors will be the focus of research in future work.

Conclusions

Drone-BSs deployment in 5G networks is currently an important subject with extreme potential in the mobile communication networks industry. Placement optimization is a complicated task that requires automation algorithms to find the optimal solution for the placement of drone-BSs. In this work, we studied a 3D deployment of drone-BSs in a 5G network using GWO algorithm. Firstly, we analysed the downlink coverage probability according to SINR and altitude. Next, we proposed a mathematical model in order to maximize coverage under a minimum number of drone-BSs, altitude and total limiting available bandwidth. The results obtained showed a minimization of the total number of drone-BSs deployed with a coverage reaching 82.5%. Our study does not take into consideration the blocking effect and overlap between drones that are important for the handoff, but these factors will be the focus of research in future work.

Figure 1 .

 1 Figure 1. System Model.

Figure 1 .

 1 Figure 1. System Model.

Figure 2 .

 2 Figure 2. Coverage probability versus SINR threshold, at different altitudes.

Figure 3 .

 3 Figure 3. Drone-BSs placement at the end of the GWO iteration.

Figure 2 .

 2 Figure 2. Coverage probability versus SINR threshold, at different altitudes.

Figure 2 .

 2 Figure 2. Coverage probability versus SINR threshold, at different altitudes.

Figure 3 .

 3 Figure 3. Drone-BSs placement at the end of the GWO iteration.Figure 3. Drone-BSs placement at the end of the GWO iteration.

Figure 3 .

 3 Figure 3. Drone-BSs placement at the end of the GWO iteration.Figure 3. Drone-BSs placement at the end of the GWO iteration.

Figure 4 .

 4 Figure 4. Number of covered users and drone-BSs deployment at each iteration.

Figure 4 .

 4 Figure 4. Number of covered users and drone-BSs deployment at each iteration.

Table 1 .

 1 Simulation parameters adopted for our work.

	Parameters	Value
	f c	28 GHz
	h min ≤ h ≤ h max	1000 ≤ h ≤ 3000 m
	M (number of drone-BSs)	10
	N (number of users)	200
	p t	30 dBm
	SI NR thr	5 dBm
	B	20 MHz

Acknowledgments:

The authors wish to acknowledge the technical support, infrastructure and documentation provided by Nokia Networks.

Acknowledgments:

The authors wish to acknowledge the technical support, infrastructure and documentation provided by Nokia Networks.

Funding: This Project is financed by Politehnica University Timisoara, Timisoara 300006, Romania Institutional Review Board Statement: Not applicable. Informed Consent Statement: Not applicable. Funding: This Project is financed by Politehnica University Timisoara, Timisoara 300006, Romania. Institutional Review Board Statement: Not applicable. Informed Consent Statement: Not applicable.

Data Availability Statement: Not applicable.

Data Availability Statement: Not applicable.

Conflicts of Interest:

The authors declare no conflict of interest.