This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

3D Placement and Orientation of mmWave-based UAVs for Guaranteed LoS Coverage

Javad Sabzehali, Vijay K. Shah, Harpreet S. Dhillon, and Jeffrey H. Reed The authors are with Wireless@VT, The Bradley Department of ECE at Virginia Tech, Blacksburg, VA (emails: {jsabzehali, vijays, hdhillon, reedjh}@vt.edu). This work is partially supported by NSF Grant CNS-1923807 and Commonwealth Cyber Initiative (CCI), an investment in the advancement of cyber R&D, innovation, and workforce development. For more information about CCI, visit www.cyberinitiative.org
Abstract

Unmanned aerial vehicles (UAVs), as aerial base stations, are a promising solution for providing wireless communications, thanks to their high flexibility and autonomy. Moreover, emerging services, such as extended reality, require high-capacity communications. To achieve this, millimeter wave (mmWave), and recently, terahertz bands have been considered for UAV communications. However, communication at these high frequencies requires a line-of-sight (LoS) to the terminals, which may be located in 3D space and may have extremely limited direct-line-of-view (LoV) due to blocking objects, like buildings and trees. In this paper, we investigate the problem of determining 3D placement and orientation of UAVs such that users have guaranteed LoS coverage by at least one UAV and the signal-to-noise ratio (SNR) between the UAV-user pairs are maximized. We formulate the problem as an integer linear programming(ILP) problem and prove its NP-hardness. Next, we propose a low-complexity geometry-based greedy algorithm to solve the problem efficiently. Our simulation results show that the proposed algorithm (almost) always guarantees LoS coverage to all users in all considered simulation settings.

Index Terms:
Unmanned Aerial Vehicles, mmWave Communications, Line-of-Sight coverage, Directional Antenna

I Introduction

Unmanned Aerial Vehicles (UAVs), also known as drones, have numerous applications in telecommunications, rescue operations, aerial sensing, to name a few. Significant research efforts have gone into understanding the role of UAVs as aerial base stations (BSs) to complement the coverage and capacity of terrestrial wireless networks. mmWave and future teraHertz communications with orders of magnitude larger bandwidths than the conventional sub-6 GHz systems are expected to play a pivotal role in supporting emerging applications, such as extended reality, in UAV-assisted communications networks [1]. However, characterized by very high path loss in non-line-of-sight (NLoS) conditions, it is critical that mmWave (and teraHertz) based UAV communications always have a line-of-sight (LoS) coverage to wireless devices  [2][3].

To provision high quality wireless LoS coverage to wireless devices (or users), UAVs should jointly adjust their coordinates and the orientations of their directional antennas. There exist some work  [4, 5, 6, 7, 8, 9] that study the problem of communication coverage for UAVs, however, there are several limitations: (i) Existing work consider wireless devices on the ground (i.e., horizontal plane) whereas our work considers wireless devices in 3D space. (ii) Existing work model air-to-ground path loss model as a probabilistic LoS and NLoS model in sub-6 GHz communications. However, as discussed in [3, 2], to fully take advantage of mmWave-based communication, especially for UAVs, the receiver must be located such that it can secure a LoS communication with the transmitter. Our 3D LoS coverage model attempts to always guarantee LoS, a key requirement of mmWave communications. (iii) And finally, most of existing works build UAV wireless coverage model as a disk. In this paper, we propose to utilize 3D directional coverage model [10] to realistically model mmWave-based UAV LoS and user direct line-of-view (LoV) coverage (See Fig. 1). Following this, we investigate the joint problem of 3D placement and orientation of UAVs such that (i) users with limited direct LoV region have guaranteed LoS coverage111Note that this work only considers LoS paths, however, they may be some strong NLoS paths and is left for future studies. by at least one UAV, (ii) the SNR between the UAV and its associated users are maximized, and (iii) the number of UAVs required is minimized. To our understanding, this is the first work which studies the guaranteed LoS coverage problem for UAV mmWave communications, where users are considered in 3D space and has limited directed LoV region.

Refer to caption
Figure 1: System Model. The red and blue cones respectively represent Line-of-view (LoV) region of users and Line-of-sight (LoS) region of UAVs.

The key contributions of this work are as follows:

  • We first present LoS and LoV coverage model for UAVs and users respectively, which realistically model mmWave based UAV systems.

  • We formulate the problem of determining the 3D placement and antenna orientations of UAVs as an integer linear programming optimization problem. Then, we prove that it is NP-Hard. Next, we propose a low-complexity geometric approach for solving the problem efficiently.

  • Our extensive experimental results depict that the proposed algorithm (almost) always guarantees LoS coverage to all users irrespective of its limited LoV region, compared to that of the baseline approach222The baseline algorithm is similar to the proposed algorithm in pseudocode  1 except that the antenna orientation is fixed.. Furthermore, we notice that compared to the baseline approach, the proposed algorithm incurs fewer number of UAVs and achieves much higher average SNR under scenarios where the user’s LoV are higher and both the approaches can provide LoS coverage to the majority of users.

The rest of the paper is organised as follows: In Section II, we model the considered UAV-based mmWave wireless system and formulate the problem. In Section III, we discuss the proposed solution approach. In Section IV, we present the experimental results, followed by conclusion in Section V.

II System Model and Problem Overview

Let us consider a mmWave-based UAV communication system with |C||C| UAVs and |U||U| users (i.e., wireless devices, such as, aerial sensing drones, mobile terminals etc.) both deployed in 3D space (See Fig. 1). We assume that each user u^jU\hat{u}_{j}\in U is located at a known 3D location uju_{j} in the region. Since, there may be buildings, trees or other objects on the side, top or bottom of any given user u^j\hat{u}_{j} (at location uju_{j}), it is quite likely u^j\hat{u}_{j} may not always have a LoV from all directions. Let vector direction \vvduj=(\vvxj,\vvyj,\vvzj)\vv{d}_{u_{j}}=(\vv{x_{j}},\vv{y_{j}},\vv{z_{j}}) with angle-of-view (AoV), θj\theta_{j} (where 0<θj1800^{\circ}<\theta_{j}\leq 180^{\circ}) denotes the user LoV for the considered user u^j\hat{u}_{j}, and is known a priori.

Let cic_{i} denote the 3D location of UAV c^iC\hat{c}_{i}\in C (in case of deployment in the 3D space). We assume that UAV c^i\hat{c}_{i} is equipped with a directional antenna. Furthermore, we consider that the antenna can face any direction \vvdci\vv{d}_{c_{i}} that can be defined by two angles in 3D space (2D azimuth and elevation angles [11]) or a 3D Cartesian vector (\vvx^i,\vvy^i,\vvz^i)(\vv{\hat{x}_{i}},\vv{\hat{y}_{i}},\vv{\hat{z}_{i}}). Moreover, UAV’s directional antenna half power beamwidth (HPBW) ϕ\phi is a prespecified value between (0,180)(0^{\circ},180^{\circ}). However, this limitation has no influence to the coverage orientation, because UAV can hover in the air and rotate itself to face any orientation \vvdci\vv{d}_{c_{i}}. So in other words, considering \vvdci\vv{d}_{c_{i}} as a 3D Cartesian vector (\vvx^i,\vvy^i,\vvz^i)(\vv{\hat{x}_{i}},\vv{\hat{y}_{i}},\vv{\hat{z}_{i}}) provides 360360^{\circ} rotation capability to a certain UAV in our system model.

Just like terrestrial cellular networks, we assume that neighboring UAVs transmit on orthogonal frequencies and hence do not interfere with each other. Also, our problem setup is general and is applicable to any multiple access scheme, such as time/frequency division multiple access or orthogonal frequency division multiple access.

Path Loss Model

As discussed earlier, for the communication to be successful, we need to ensure LoS between the UAVs and the users. Given an access link between user u^j\hat{u}_{j} and UAV c^i\hat{c}_{i} located respectively at 3D coordinates uj=(xj,yj,hj)u_{j}=(x_{j},y_{j},h_{j}) and ci=(xi,yi,hi)c_{i}=(x_{i},y_{i},h_{i}), the path loss of the channel (in dB) between user u^j\hat{u}_{j} and UAV c^i\hat{c}_{i} is modeled as follows:

φijL=20log(4πfcdijc).{\varphi_{ij}^{L}=20\log\left(\frac{{4\pi{f_{c}}{d_{ij}}}}{c}\right)}. (1)

where fcf_{c} is the carrier frequency, dij=[(xjxi)2+(yjyi)2+(hihj)2]12d_{ij}=[(x_{j}-x_{i})^{2}+(y_{j}-y_{i})^{2}+(h_{i}-h_{j})^{2}]^{\frac{1}{2}} is the 3D distance between user u^j\hat{u}_{j} and UAV c^i\hat{c}_{i}, and cc is the speed of light.

Let PiP_{i} be the transmission power and GiG_{i} be the antenna gain of UAV c^i\hat{c}_{i}, and δ2\delta^{2} as the noise power. The SNR between UAV c^i\hat{c}_{i} and user u^j\hat{u}_{j} hence can be expressed as

γij=10logPi+GiφijL10logδ2.\gamma_{ij}=10\log{P_{i}}+G_{i}-\varphi_{ij}^{L}-10\log{\delta^{2}}. (2)

Note that the SNR between a user u^j\hat{u}_{j} and a UAV c^i\hat{c}_{i} determines if the user u^j\hat{u}_{j} is covered by the corresponding UAV c^i\hat{c}_{i}. In other words, a user u^j\hat{u}_{j} is within the coverage of a UAV c^i\hat{c}_{i} if γij\gamma_{ij} meets the SNR threshold γ0\gamma_{0} (i.e., γijγ0\gamma_{ij}\geq\gamma_{0}).

UAV LoS (and user LoV) coverage model

Both UAV LoS coverage model and user LoV coverage model will be based on the 3D directional coverage model (inspired by a recent work [10]), as discussed in the following.

Refer to caption
(a)
Refer to caption
(b)
Figure 2: (a) Depiction of UAV LoS coverage model and user LoV coverage model and (b) Geometric approach for determining the optimal UAV orientation.
Definition 1 (3D Directional Coverage).

A user u^j\hat{u}_{j} at 3D location uju_{j}, with its vector direction \vvduj\vv{d}_{u_{j}}, is said to be in the directional coverage of UAV c^i\hat{c}_{i} at 3D location cic_{i}, with antenna orientation \vvdci\vv{d}_{c_{i}}, and HPBW ϕ\phi, if u^j\hat{u}_{j} is covered by c^i\hat{c}_{i} (i.e., γijγ0\gamma_{ij}\geq\gamma_{0}) and (\vvduj,\vvujci)<θj\angle(\vv{d}_{u_{j}},\vv{u_{j}c_{i}})<\theta_{j}.

UAV LoS Coverage Model. According to Definition 1, UAV LoS coverage model can be established as a spherical base cone as shown in Fig. 2. Let c^i\hat{c}_{i} be the UAV with HPBW ϕ\phi, and DmaxD_{max}333For simplicity, we consider a sectorized antenna model (also called flat top antenna model) [1], where the UAV antenna gain GiG_{i} is a positive constant amount over every direction inside of the formed spherical based cone and is zero outside of it. This gives a constant maximum 3D distance between UAV c^i\hat{c}_{i} and user u^j\hat{u}_{j} in every direction to form the mentioned spherical based cone. be the maximum 3D distance between UAV c^i\hat{c}_{i} and user u^j\hat{u}_{j} at which γijγ0\gamma_{ij}\geq\gamma_{0}, then, the UAV LoS coverage model reduces to:

{|\vvujci|Dmax(\vvciuj,\vvdci)ϕ/2.\small\begin{cases}|\vv{u_{j}c_{i}}|\leq D_{max}\\ \angle(\vv{c_{i}u_{j}},\vv{d}_{c_{i}})\leq\phi/2.\end{cases} (3)

Based on Equation 3, UAV c^i\hat{c}_{i} can efficiently cover a certain user u^j\hat{u}_{j} only when it positions itself in the spherical cone of user u^j\hat{u}_{j} where |\vvujci|Dmax|\vv{u_{j}c_{i}}|\leq D_{max} and (\vvciuj,\vvdci)ϕ/2\angle(\vv{c_{i}u_{j}},\vv{d}_{c_{i}})\leq\phi/2.

User LoV coverage model. Similar to UAV LoS coverage model, user LoV coverage model can also be established as a spherical base cone as shown in Fig. 1. For simplicity, in this work, we only consider downlink connection. For downlink connection from UAV c^i\hat{c}_{i} to user u^j\hat{u}_{j}, then we only need to ensure that user u^j\hat{u}_{j} is 3D directional covered by UAV c^i\hat{c}_{i}. Hence, in this case, we only need to ensure |\vvujci|=|\vvciuj|Dmax|\vv{u_{j}c_{i}}|=|\vv{c_{i}u_{j}}|\leq D_{max}. Let user u^j\hat{u}_{j} be the user, with AoV θj\theta_{j} and DmaxD_{max} as the maximum distance at which γjiγ0\gamma_{ji}\geq\gamma_{0}, then user LoV coverage model is as follows.

{|\vvciuj|Dmax(\vvujci,\vvduj)θj.\small\begin{cases}|\vv{c_{i}u_{j}}|\leq D_{max}\\ \angle(\vv{u_{j}c_{i}},\vv{d}_{u_{j}})\leq\theta_{j}.\end{cases} (4)

Combining Equations 3 and 4, the guaranteed LoS coverage function can be defined as follows:

Flos(ci,uj,\vvdci,\vvduj)={1,|\vvciuj|Dmax,(\vvujci,\vvduj)θj,and(\vvciuj,\vvdci)ϕ/20,otherwise.\small F_{los}(c_{i},u_{j},\vv{d}_{c_{i}},\vv{d}_{u_{j}})=\begin{cases}1,|\vv{c_{i}u_{j}}|\leq D_{max},\\ \angle(\vv{u_{j}c_{i}},\vv{d}_{u_{j}})\leq\theta_{j},\text{and}\\ \angle(\vv{c_{i}u_{j}},\vv{d}_{c_{i}})\leq\phi/2\\ 0,\text{otherwise}.\end{cases} (5)

Problem Formulation

Let fc^if_{\hat{c}_{i}} denote whether UAV c^i\hat{c}_{i} is located at 3D location ci=(xi,yi,hi)c_{i}=(x_{i},y_{i},h_{i}) or not, and \vvdci\vv{d}_{c_{i}} denote the directional antenna orientation of UAV c^i\hat{c}_{i} (at 3D location cic_{i}). Let βij\beta_{ij} denote that user u^j\hat{u}_{j} (at known 3D location uju_{j}) is associated with UAV c^i\hat{c}_{i}. Then, the key objective function is to determine a set of strategies, i.e., ciCci,\vvdci\cup_{c_{i}\in C}\langle c_{i},\vv{d}_{c_{i}}\rangle, that maximizes the SNR between each UAV and its associated users while minimizing the number of deployed UAVs.

P0:max<ci,\vvdci>c^iCu^jUγijβijc^iCfc^i\displaystyle\textbf{P0}:\mathop{\max_{<c_{i},\vv{d}_{c_{i}}>}}\hskip 2.84526pt\frac{\sum_{\hat{c}_{i}\in C}\sum_{\hat{u}_{j}\in U}\gamma_{ij}\beta_{ij}}{\sum_{\hat{c}_{i}\in C}f_{\hat{c}_{i}}} (6)
s.t.c^iCβij.Flos(ci,uj,\vvdci,\vvduj)=1\displaystyle s.t.\sum_{\hat{c}_{i}\in C}\beta_{ij}.F_{los}(c_{i},u_{j},\vv{d}_{c_{i}},\vv{d}_{u_{j}})=1 (7)

Constraint 7 ensures that a certain wireless user is associated with one unique UAV among all the UAVs which provides LoS coverage to the user.

Theorem 1.

The problem P0 is NP-Hard.

Proof.

Assume the antenna is an omnidirectional antenna, i.e., ϕ=360\phi=360^{\circ} and the user has direct line-of-view (LoV) from all directions, i.e., θ=180\theta=180^{\circ}. Moreover, also consider that users are placed on the ground (i.e., horizontal plane). Under these considerations, the problem P0 is transformed into determining the minimum number of UAVs and their optimal 3D coordinates that ensure 100%100\% ground user coverage. This is a case of geometric set cover problem, which is a classic NP-Hard problem [12]. Therefore, since the specific instance of P0 is NP-Hard, P0 itself is NP hard. ∎

Problem relaxation under practical considerations

The set of strategies, Π=ciCci,\vvdci\Pi=\cup_{c_{i}\in C}\langle c_{i},\vv{d}_{c_{i}}\rangle is an exponentially large solution space. This is because the 3D coordinate space is infinite and also, the orientation of directional antenna \vvdci\vv{d}_{c_{i}} is an infinite space. However, under practical scenario, we consider the following relaxations to the problem:

  • The 3D coordinate space is discretized in |𝕂||\mathbb{K}| equal-sized cubic grids (or discrete locations) as usually considered in the literature [13], and the 3D coordinate location for a certain UAV must belong to a certain grid. Let c^k\hat{c}_{k} denote a certain UAV (in set CC) located at grid k𝕂k\in\mathbb{K}.

  • By considering the information known about the users’ location and its LoV, we can use an efficient geometrical approach (discussed in the next subsection) to obtain a finite UAV antenna orientation space.

Geometrical Approach

For each user u^jUck\hat{u}_{j}\in U_{c_{k}} (where UckU_{c_{k}} is the set of all users that is LoS covered by a UAV placed at a grid location k𝕂k\in\mathbb{K}), we do the following: First, we generate an arbitrary orientation for deploying a UAV c^k\hat{c}_{k} at candidate location kk so that the selected user u^j\hat{u}_{j} lie on the surface of the UAV c^k\hat{c}_{k} LoS coverage region (See Fig. 2(b)). Then, by keeping the user u^j\hat{u}_{j} on the surface of the UAV c^k\hat{c}_{k} LoS coverage region, we will change the orientation of UAV c^k\hat{c}_{k}. It implies that the line ckujc_{k}u_{j} will remain on the surface of all UAV c^k\hat{c}_{k} LoS coverage regions. Since the antenna HPBW ϕ\phi is assumed to be constant in all cases, the orientations of all UAV placements will consist of the surface of a cone AA with radius rr and height |\vvckuj||\vv{c_{k}u_{j}}| (See Fig. 2(b)). To figure out these orientations, we need to figure out the equation of the circumference of the base of the cone AA (i.e. circle cc with center uju_{j} and radius rr). Given the location of grid kk, if we know all points lying on the circle cc, then the line passing kk and a point on circle cc can be uniquely determined.

By knowing the normal vector of a plane and a point on that plane, we can also uniquely figure out the equation of that plane. Hence, the equation of the plane that circle cc lies on can be determined by knowing the location of uj=(x^,y^,z^)u_{j}=(\hat{x},\hat{y},\hat{z}) and \vvujck\vv{u_{j}c_{k}}, which is the normal vector of the plane (i.e. \vvn=(u,v,w)\vv{n}=(u,v,w)). Therefore, the equation of the plane would be u.x+v.y+w.z=tu.x+v.y+w.z=t. By plugging uju_{j} in this equation, i.e. u.x^+v.y^+w.z^=tu.\hat{x}+v.\hat{y}+w.\hat{z}=t, tt will be found uniquely. Now we know the equation of the plane that circle cc lies on, we can define a random vector \vva\vv{a} on the plane such that passes through uju_{j} and its length be equal to one. Let \vvb\vv{b} be the cross product of \vvn\vv{n} and \vva\vv{a}, i.e. \vvb=\vva×\vvn\vv{b}=\vv{a}\times\vv{n}, that will be another vector on the plane perpendicular to both \vvn\vv{n} and \vva\vv{a}. Then c{r.(\vva.sin(ω)+\vvb.cos(ω))+uj}c\in\{r.(\vv{a}.sin(\omega)+\vv{b}.cos(\omega))+u_{j}\}, where 0ω<3600^{\circ}\leq\omega<360^{\circ}.

Observation 1.

Given ckc_{k} lies inside the LoV region of user u^j\hat{u}_{j} at 3D location uju_{j}, |\vvckuj|Dmax|\vv{c_{k}u_{j}}|\leq D_{max}. It means that placing a UAV c^k\hat{c}_{k} at grid location kk in some certain orientation like \vvdck=\vvckuj\vv{d}_{c_{k}}=\vv{c_{k}u_{j}}, user u^j\hat{u}_{j} will be LoS covered by c^k\hat{c}_{k} (See Eq. 5). Hence, if a UAV c^k\hat{c}_{k} placed at a certain location k𝕂k\in\mathbb{K} lies inside the LoV region (spherical cone base shape) of a user u^j\hat{u}_{j}, then there always exist a feasible orientation for a UAV c^k\hat{c}_{k} that provides LoS coverage to the user u^j\hat{u}_{j}.

III Proposed Algorithm

In this section, we propose a low-complexity greedy algorithm based on the aforementioned geometric approach to solve the problem efficiently.

As depicted in pseudocode 1, the algorithm computes the set of all users, UckU_{c_{k}}, that can be LoS covered by a UAV placed at a discrete candidate location k𝕂k\in\mathbb{K} with at least one feasible antenna orientation (\vvx^k,\vvy^k,\vvz^k)(\vv{\hat{x}_{k}},\vv{\hat{y}_{k}},\vv{\hat{z}_{k}}). Let 𝕂LoV𝕂\mathbb{K}_{LoV}\subseteq\mathbb{K} denote the set of such candidate grid locations (See Observation 1 and lines 3-7). Notice that, UAV c^k\hat{c}_{k} may not necessarily provide LoS coverage to all the users in UckU_{c_{k}} with any one unique orientation (among all possible orientations). Thus, the algorithm further computes the subset of users U<ck,dckui>UckU_{<c_{k},d^{u_{i}}_{c_{k}}>}\in U_{c_{k}} LoS covered by orientation dckuid^{u_{i}}_{c_{k}} that cover the largest possible subset of UckU_{c_{k}} while necessarily a user u^i\hat{u}_{i} is being covered (See lines 11 - 12). To compute U<ck,dckui>U_{<c_{k},d^{u_{i}}_{c_{k}}>} and dckuid^{u_{i}}_{c_{k}}, the algorithm uses the geometric approach described previously in subsection II-e.

By repeating this procedure for every single user u^iUck\hat{u}_{i}\in U_{c_{k}}, the algorithm selects the case in which the most number of users can be covered by grid location kk, i.e., subset UckmaxU^{max}_{c_{k}}, and corresponding orientation \vvdckmax\vv{d}^{max}_{c_{k}}. We also compute the average of SNR between the UAV c^k\hat{c}_{k} and each user u^jUckmax\hat{u}_{j}\in U^{max}_{c_{k}}, and term it γckavg\gamma^{avg}_{c_{k}} for a UAV at grid location kk and selected orientation \vvdckmax\vv{d}^{max}_{c_{k}} (See lines 9-14 of pseudocode 1). After repeating this procedure for all grid locations, among the cases that cover maximum number of users, we will pick the case which has the largest γckavg\gamma^{avg}_{c_{k}} (See lines 15-18 of pseudocode 1). We delete associated users to the selected grid location cselc_{sel}, which is UcselmaxU^{max}_{c_{sel}}, from all UkU_{k}s (See lines 19-21 of pseudocode 1). It means that each user’s coverage ensures only once in the proposed algorithm. These steps will be repeated until no users remains such that could be covered by grid candidate locations in the space or the total number of UAVs CC has been deployed. It is worth mentioning that |kUck||\cup_{k}{U_{c_{k}}}| will be a decreasing quantity after each iteration, and thus will reach zero in countable number of iterations. Therefore, the algorithm 1 is terminated in a polynomial-time running complexity as also analyzed in Theorem 2. In our work, we consider CC to be sufficiently large (i.e., |C|>𝕂|C|>\mathbb{K}), however, a certain network provider may consider lesser value of CC for enormous cost reasons. Nevertheless, the proposed algorithm will always provide a solution with minimum number of UAVs (|C|\leq|C|) with LoS coverage to almost all users and also, high SNR between UAV and users.

Input: UU users, CC UAVs, and 𝕂\mathbb{K} grid locations.
Output: Π\Pi set of strategies

1:Initialize strategy set Π=ϕ\Pi=\phi, Usofar=ϕU_{so-far}=\phi
2:Set of LoV candidate grids, 𝕂LoV=ϕ\mathbb{K}_{LoV}=\phi
3:for a user, u^iU\hat{u}_{i}\in U do
4:    for a UAV at grid location, k𝕂k\in\mathbb{K} do
5:       if kk lies inside LoV region of user u^i\hat{u}_{i} then
6:          Uck=Uck{u^i}U_{c_{k}}=U_{c_{k}}\cup\{\hat{u}_{i}\}
7:          𝕂LoV=𝕂LoVck\mathbb{K}_{LoV}=\mathbb{K}_{LoV}\cup c_{k}            
8:while |kUck|0|\cup_{k}{U_{c_{k}}}|\neq 0 and |kc^k||C||\cup_{k}\hat{c}_{k}|\leq|C| do
9:    for a UAV at grid location, k𝕂LoVk\in\mathbb{K}_{LoV} do
10:       for a user, u^iUck\hat{u}_{i}\in U_{c_{k}} do
11:          dckui=maxx𝒟ckuixd^{u_{i}}_{c_{k}}=\underset{x\in\mathcal{D}^{u_{i}}_{c_{k}}}{max}x, where 𝒟ckui\mathcal{D}^{u_{i}}_{c_{k}} is the set of all possible antenna orientations with user u^i\hat{u}_{i} covered by a UAV c^k\hat{c}_{k}
12:          Compute list of users U<ck,dckui>U_{<c_{k},d^{u_{i}}_{c_{k}}>} covered by UAV c^k\hat{c}_{k} with antenna orientation dckuid^{u_{i}}_{c_{k}}        
13:       <Uckmax,dckmax>={<Uckui,dckui>maxu^iUckU<ck,dckui>}<U^{max}_{c_{k}},d^{max}_{c_{k}}>=\{<U^{u_{i}}_{c_{k}},d^{u_{i}}_{c_{k}}>\mid\underset{\hat{u}_{i}\in U_{c_{k}}}{max}U_{<c_{k},d^{u_{i}}_{c_{k}}>}\}
14:       Average SNR, γckavg=u^jUckmaxγckuj|Uckmax|\gamma^{avg}_{c_{k}}=\frac{\underset{\hat{u}_{j}\in U^{max}_{c_{k}}}{\sum}\gamma_{c_{k}u_{j}}}{|U^{max}_{c_{k}}|}
15:       Usofar=Usofar<ck,dckmax,Uckmax,γckavg>U_{so-far}=U_{so-far}\cup<c_{k},d^{max}_{c_{k}},U^{max}_{c_{k}},\gamma^{avg}_{c_{k}}>     
16:    Sort UsofarU_{so-far} based on non-increasing order of UckmaxU^{max}_{c_{k}} followed by non-increasing order of γckavg\gamma^{avg}_{c_{k}}
17:    Select the first item in sorted UsofarU_{so-far} as the chosen strategy, i.e., <csel,dcsel>=<Usofar[0][0],Usofar[0][1]><c_{sel},d_{c_{sel}}>=<U_{so-far}[0][0],U_{so-far}[0][1]>
18:    Π=Π<csel,dcsel>\Pi=\Pi\cup<c_{sel},d_{c_{sel}}>
19:    for a UAV at grid location, k𝕂LoVk\in\mathbb{K}_{LoV} do
20:       for user, u^iUsofar[0][2]\hat{u}_{i}\in U_{so-far}[0][2] do
21:          Uck=Uck{u^i}U_{c_{k}}=U_{c_{k}}\setminus\{\hat{u}_{i}\}            
22:return Π\Pi
Algorithm 1 Proposed Algorithm
Theorem 2.

The time complexity of the proposed algorithm is O(|U|2.|𝕂|+|U|.|𝕂|.log(|𝕂|)O(|U|^{2}.|\mathbb{K}|+|U|.|\mathbb{K}|.\log(|\mathbb{K}|).

Proof.

The time complexity for lines 373-7 of the proposed algorithm is O(|U|.|𝕂|)O(|U|.|\mathbb{K}|). Lines 10-12 have time complexity of O(|U|)O(|U|) followed by the same time complexity for lines 131413-14, while line 15 takes O(1)O(1). Therefore, the running time of lines 9159-15 will be O(|U|.|𝕂|)O(|U|.|\mathbb{K}|). The cost of line 1616 is O(|𝕂|.log(|𝕂|)O(|\mathbb{K}|.\log(|\mathbb{K}|), followed by O(1)O(1) for lines 171817-18. The running time for lines 192119-21 is O(|U|.|𝕂|)O(|U|.|\mathbb{K}|). Hence, the runtime complexity of lines 8218-21 is O(|U|.|𝕂|.(|U|+log(|𝕂|))=O(|U|2.|𝕂|+|U|.|𝕂|.log(|𝕂|)O(|U|.|\mathbb{K}|.(|U|+log(|\mathbb{K}|))=O(|U|^{2}.|\mathbb{K}|+|U|.|\mathbb{K}|.\log(|\mathbb{K}|), which will also be the worst-case time complexity of the proposed algorithm. ∎

IV Simulation Results

In this section, we discuss the comparative analysis of proposed algorithm against a simple baseline approach (discussed below) in terms of key performance metrics, namely, percentage of users LoS covered, number of UAVs needed and average SNR achieved. For the baseline approach, we consider that the UAV antenna is always facing downward (i.e., azimuth angle to 00^{\circ} and elevation angle to 900-90^{0} or 3D cartesian vector is (0,0,1)(0,0,-1)), which is usually considered in the literature [14]. Notice that this baseline approach does not take into account the limited LoV region of each user, and thus, has limitations in providing guaranteed LoS coverage to all users as discussed later in this section.

Refer to caption
(a)
Refer to caption
(b)
Refer to caption
(c)
Figure 3: Varying user AoV: (a) Avg. perc. of users covered, (b) Avg. no. of UAVs per covered users, and (c) Avg. SNR achieved.

Simulation Setting. Unless otherwise stated, 100100 users are uniformly at random in clusters, where each cluster contains 10-15 devices, in a 3D simulation region of (length ×\times breadth ×\times altitude =1000×1000×100=1000\times 1000\times 100 m3), whose base corresponds to the ground plane. The simulation region is divided in 𝕂\mathbb{K} 3D cubic grids, where each cubic grid is 20×20×2020\times 20\times 20 m3. These cubic grids are the only feasible locations of a UAV in the considered 3D simulation region444We have defined extra grid locations outside of the considered 3D space region such that they are placed within a specific 5050 m outside the boundaries of the 3D region to cover the users near the boundaries which have a LoV region outside the considered 3D region.. We consider that each user’s angle-of-view (AoV) is uniformly at random in a unique interval from a set of intervals – {[15,45]\{[15^{\circ},45^{\circ}], [30,90][30^{\circ},90^{\circ}], [60,120][60^{\circ},120^{\circ}], [90,150][90^{\circ},150^{\circ}], [120,180][120^{\circ},180^{\circ}], [180[180^{\circ},180]}180^{\circ}]\} for each simulation setting. This represents the fact that the user’s direct LoV may be blocked from due to buildings, trees or other objects. Major simulation parameters are listed in Table I.

TABLE I: Simulation parameters
Parameters Value
Simulation 3D region 10001000 m ×\times 10001000 m ×\times 100 m
Number of Users 100100
Carrier frequency (fcf_{c}) 6060 GHz
UAV Transmit Power (PtP_{t}) 0.50.5 W
Minimum UAV Antenna Gain (GtG_{t}) 1010 dB
Bandwidth 1 GHz
Noise Power Spectral Density -174 dBm/Hz
SNR Threshold for user coverage (γ0\gamma_{0}) 1515 dB
UAV antenna HPBW (ϕ\phi) 4545^{\circ}

Note using the parameters listed in Table I, the SNR γij\gamma_{ij} between a UAV c^i\hat{c}_{i} and user u^j\hat{u}_{j} can be computed using Eq. 2. Now given γij\gamma_{ij}, known user u^j\hat{u}_{j}’s AoV (θj\theta_{j}), UAV antenna HPBW ϕ\phi, the maximum 3D distance (DmaxD_{max}) between UAV c^i\hat{c}_{i} and user u^j\hat{u}_{j} can be calculated, where DmaxD_{max} is 80\approx 80 m.

Experimental Results. Results are reported for percentage of covered users, number of UAVs needed, and average SNRs against varying user AoV. We also conducted experiments with varying number of users, region sizes, SNR thresholds, and antenna HPBWs, and the results obtained show similar behavior for proposed and baseline approaches as that of varying user AoVs. (Not shown due to space constraint.)

Average percentage of users LoS covered. As shown in Fig. 3(a), the percentage of LoS covered users increases with increasing user AoV. This is intuitive as the probability of a user to be in LoS coverage by a certain UAV increases with increasing user AoV, given the other parameters (i.e., SNR threshold and UAV HPBW) are constant (See Equation 5). It is important to note that the proposed approach significantly outperforms the baseline approach, particularly for limited user AoV (e.g., [15,45][15^{\circ},45^{\circ}]). This is because unlike the baseline approach, the proposed approach specifically accounts for users’ LoV, and chooses (sub)optimal 3D placement and antenna orientation of UAVs to guarantee LoS coverage to (almost) all users.

Average number of UAVs needed per number of users LoS covered. Fig 3(b) shows that, with increasing user AoV, the proposed algorithm requires fewer number of UAVs to provide guaranteed LoS coverage to all users. This is again because a certain UAV can cover larger number of users when user AoV is larger (See Eq. 5), and thus, on average fewer UAVs are required to cover all users in the considered 3D region. Besides providing LoS coverage to less number of users in the baseline approach compared to the proposed approach (as shown in Fig. 3(a)), the average number of UAVs needed per number of users LoS covered in the proposed algorithm is significantly less than the baseline approach, thus showcasing the superiority of the proposed algorithm.

Average SNRs achieved. As shown in Fig. 3(c), by considering an SNR of zero for the users out of the coverage (or in other words, users in NLoS region), our approach experiences significantly higher average SNR at users when compared to that of the baseline approach. This is expected as the proposed algorithm determines the suitable 3D placement and antenna orientation for UAVs such that the average SNR experienced by users are maximized (See Obj. function in Equation 6).

V Conclusion

In this letter, we investigated the design of joint 3D placement and orientation of high frequency band (e.g., mmWave or teraHertz) UAVs which require guaranteed line-of-sight (LoS) coverage to users – with limited direct line-of-view (LoV) due to various blocking objects around it. We first formulated the aforestated joint problem as an integer linear programming (ILP) optimization problem and proved its NP-Hardness. By utilizing the geometric properties of the setup, a low-computational complexity algorithm was presented to solve the problem. Our comprehensive evaluations showed the superiority of the proposed algorithm in providing LoS coverage for a near 100%100\% of the users while experiencing significantly higher SNRs, when compared to that of baseline approach, for all considered simulation setting.

References

  • [1] S. Tripathi, N. V. Sabu, A. K. Gupta, and H. S. Dhillon, “Millimeter-wave and Terahertz Spectrum for 6G Wireless,” 6G Mobile Wireless Networks. Springer, 2021.
  • [2] M. T. Dabiri, H. Safi, S. Parsaeefard, and W. Saad, “Analytical channel models for millimeter wave uav networks under hovering fluctuations,” IEEE Trans. on Wireless Communications, vol. 19, no. 4, pp. 2868–2883, 2020.
  • [3] Y. Wang, K. Venugopal, A. F. Molisch, and R. W. Heath, “Mmwave vehicle-to-infrastructure communication: Analysis of urban microcellular networks,” IEEE Transactions on Vehicular Technology, vol. 67, no. 8, pp. 7086–7100, 2018.
  • [4] E. Kalantari, H. Yanikomeroglu, and A. Yongacoglu, “On the number and 3d placement of drone base stations in wireless cellular networks,” in IEEE Vehicular Technology Conference (VTC-Fall), 2016, pp. 1–6.
  • [5] B. Perabathini, K. Tummuri, A. Agrawal, and V. S. Varma, “Efficient 3d placement of uavs with qos assurance in ad hoc wireless networks,” in 2019 28th International Conference on Computer Communication and Networks (ICCCN), 2019, pp. 1–6.
  • [6] R. I. Bor-Yaliniz, A. El-Keyi, and H. Yanikomeroglu, “Efficient 3-d placement of an aerial base station in next generation cellular networks,” in 2016 IEEE International Conference on Communications (ICC), 2016, pp. 1–5.
  • [7] M. Mozaffari, W. Saad, M. Bennis, and M. Debbah, “Efficient deployment of multiple unmanned aerial vehicles for optimal wireless coverage,” IEEE Comm. Letters, vol. 20, no. 8, pp. 1647–1650, 2016.
  • [8] W. Yi, Y. Liu, M. Elkashlan, and A. Nallanathan, “Modeling and coverage analysis of downlink uav networks with mmwave communications,” in 2019 IEEE International Conference on Communications Workshops (ICC Workshops), 2019, pp. 1–6.
  • [9] J. Lyu, Y. Zeng, R. Zhang, and T. J. Lim, “Placement optimization of uav-mounted mobile base stations,” IEEE Communications Letters, vol. 21, no. 3, pp. 604–607, 2017.
  • [10] W. Wang, H. Dai, C. Dong, X. Cheng, X. Wang, P. Yang, G. Chen, and W. Dou, “Placement of unmanned aerial vehicles for directional coverage in 3d space,” IEEE/ACM Trans. on Networking, vol. 28, no. 2, pp. 888–901, 2020.
  • [11] S. Garcia Sanchez, S. Mohanti, D. Jaisinghani, and K. R. Chowdhury, “Millimeter-wave base stations in the sky: An experimental study of uav-to-ground communications,” IEEE Trans. on Mobile Computing, 2020.
  • [12] V. E. Brimkov, A. Leach, J. Wu, and M. Mastroianni, “Approximation algorithms for a geometric set cover problem,” Discrete Applied Mathematics, vol. 160, no. 7, pp. 1039 – 1052, 2012.
  • [13] L. Bertizzolo, T. X. Tran, B. Amento, B. Balasubramanian, R. Jana, H. Purdy, Y. Zhou, and T. Melodia, “Live and let live: Flying uavs without affecting terrestrial ues,” in International Workshop on Mobile Computing Systems and Applications, 2020, pp. 21–26.
  • [14] P. Yu, W. Li, F. Zhou, L. Feng, M. Yin, S. Guo, Z. Gao, and X. Qiu, “Capacity enhancement for 5g networks using mmwave aerial base stations: Self-organizing architecture and approach,” IEEE Wireless Communications, vol. 25, no. 4, pp. 58–64, AUGUST 2018.