Maintaining Connectivity of MANETs through Multiple Unmanned Aerial Vehicles

Recently, Unmanned Aerial Vehicles (UAVs) have emerged as relay platforms to maintain the connectivity of ground mobile ad hoc networks (MANETs). However, when deploying UAVs, existing methods have not consider one situation that there are already someUAVs deployed in the field. In this paper, we study a problem jointing themotion control of existingUAVs and the deployment of new UAVs so that the number of new deployed UAVs to maintain the connectivity of ground MANETs can be minimized. We firstly formulate the problem as a Minimum Steiner Tree problem with Existing Mobile Steiner points under Edge Length Bound constraints (MST-EMSELB) and prove the NP completeness of this problem. Then we propose three Existing UAVs Aware (EUA) approximate algorithms for the MST-EMSELB problem: Deploy-Before-Movement (DBM), Move-Before-Deployment (MBD), and Deploy-Across-Movement (DAM) algorithms. Both DBM and MBD algorithm decouple the joint problem and solve the deployment and movement problem one after another, while DAM algorithm optimizes the deployment and motion control problem crosswise and solves these two problems simultaneously. Simulation results demonstrate that all EUA algorithms have better performance than non-EUA algorithm. The DAM algorithm has better performance in all scenarios than MBD and DBM ones. Compared with DBM algorithm, the DAM algorithm can reduce at most 70% of the new UAVs number.


Introduction
Unmanned Aerial Vehicles (UAVs) have emerged as promising relay platforms to improve networking performance (such as connectivity and throughput) for ground mobile ad hoc networks (MANETs) [1].UAVs have some unique characteristics suitable for the relaying task.Firstly, the motion flexibility of UAVs can expand the scope of ground-based networks especially in scenarios with obstacles.Secondly, UAVs can communicate with ground nodes in the line of sight, which can improve the link capacity between ground nodes.Last but not least, UAVs are often integrated with Communication, Computation, and Control (3C) system and various sensors, so that UAVs can be aware of the environment and control their motion adaptively.The adaptability of UAVs makes them suitable in providing relay service for MANETs that has dynamic network topology.
A variety of efforts have been made to explore the benefits of using UAVs as communication relays for ground MANETs.Some work wants to optimize the deployment of UAVs to improve the connectivity of ground nodes.Chandrashekar et al. presented a method on deploying minimum number of UAVs to connect a disconnected MANET [2].Some works study the motion control of UAVs to improve the link capacity of ground nodes.Jiang and Swindlehurst proposed a UAV's heading control algorithm that can maximize the link capacity of ground-to-air uplink channel using a multiantenna UAV [3].Dixon and Frew proposed a motion control method for using chains of UAVs to improve the link capacity between two isolated ground nodes [4].Some work considers both the deployment and motion control of UAVs.For example, Han et al. considered both the deployment and motion control of a single UAV to improve the connectivity of ground MANETs [1].
However, existing works on the deployment of UAVs have not considered a situation that some UAVs have already been deployed in the field.With the movement of ground MANETs, existing UAVs may fail to connect all ground nodes.New UAVs need to be supplied to maintain the connectivity of ground MANETs.In order to minimize the number of new added UAVs, both the movement of existing UAVs and the deployment of new UAVs need to be considered.This is a joint optimization problem that optimizes both the deployment and motion control of multiple UAVs.
We consider the usage of existing UAVs by moving them to proper positions so that the number of new UAVs that is needed can be reduced.The existing UAVs have a limited motion range that depends on the speed of UAVs.In order to support bidirectional communication between UAVs and ground nodes, we assume that UAVs have the same communication range as ground nodes.Figure 1 shows an example of how motion control of existing UAVs may reduce the number of needed new UAVs.Suppose that there are two ground nodes and two existing UAVs deployed in the field.Since the distance between two ground nodes is larger than their communication range , the ground MANET is separated into two parts which is shown in Figure 1(a).
In order to maintain the connectivity of ground MANETs, methods that do not consider existing UAVs such as [2] will add new UAVs to connect partitioned parts as shown in Figure 1(b).Here a new UAV is added and deployed in the middle of two ground nodes.So these two ground nodes can now communicate with each other with the help of the new UAV.If we do not consider existing UAVs, at least one additional UAV needs to be deployed to maintain the connectivity of ground MANETs.
To reduce capital expenditure, users would try their best to reduce the number of newly deployed UAVs.In other words, they will exploit the existing UAVs instead of ignoring them.By moving existing UAVs to proper positions and using existing UAVs as relays, the connectivity of ground MANETs can be improved.Figure 1(c) shows one proper movement.Move these two existing UAVs directly towards the line that consisted of two ground nodes until the distance between existing UAVs and one ground node is less than the communication range .Then a communication link is set up between two ground nodes.Thus the connectivity of ground MANETs is maintained and no new UAVs are needed to be deployed.
In this paper, we study the joint optimization problem of deployment and motion control of multiple UAVs so that the number of new added UAVs can be minimized.We firstly formulate this problem as a Minimum Steiner Tree problem with Existing Mobile Steiner nodes under Edge Length Bound constraints (MST-EMSELB) and prove the NP completeness of the problem.Then we present non-Existing UAVs Aware (non-EUA) algorithm and propose three Existing UAVs Aware (EUA) polynomial time approximation algorithms: Deploy-Before-Movement (DBM), Move-Before-Deployment (MBD), and Deploy-Across-Movement (DAM).The first two algorithms decouple the joint problem into deployment problem of new UAVs and motion control problem of existing UAVs.DBM algorithm optimizes the deployment of new UAVs before movement of existing UAVs and the MBD algorithm solves the problem contrarily.DAM algorithm is a mixed algorithm that solves the movement and deployment problem crossly.Simulation experiments show that all EUA algorithms have better performance in terms of new UAVs' number than non-EUA algorithm.DAM algorithm is always better than DBM and MBD algorithms and can improve the performance to 70% at most comparing with DBM algorithm.
The main contributions of this paper are as follows: ( (3) We compare the performance of proposed algorithms with non-EUA algorithms in simulation environment and demonstrate the effectiveness of proposed algorithms.
The rest of this paper is organized as follows.Section 2 presents some related work.We present the system model in Section 3 and formulate the problem in Section 4. In Section 5, we present three polynomial time approximation algorithms for the problem.We demonstrate the performance of proposed algorithm through simulation in Section 6. Section 7 concludes this work.

Related Work
Related works lie in two research fields: static relay deployment problem and mobile relay motion control problem.

Static Relay Deployment Problem.
Static relay deployment problems have been widely studied in wireless sensor networks (WSN).In WSN, it is hard to recharge the sensors after deployment.Due to the unbalanced load of message routing, some sensors will run out of energy before others.Then network partition happens and the whole network may be not available even though some sensors are still alive.
To extend the lifetime of WSN, some researches propose to deploy some static relay nodes in the field.Lloyd and Xue studied the problem of deploying minimum number of relay nodes so that, for each pair of sensor nodes, there is a connecting path consisting of relay and/or sensor nodes [5].They proved that the problem is NP complete and propose a polynomial time 7-approximation algorithm for the problem.Zhang et al. studied four related fault-tolerant relay node placement problems, discussed their computational complexity, and present a polynomial time O(1)-approximation algorithm with a small approximation ratio [6].Based on previous research, Misra et al. studied the constrained versions of the relay node placement problem, where relay nodes can only be placed at a set of candidate locations [7].They also proposed a polynomial time O(1)-approximation algorithm for the problem.Lee and Younis extend the usage of relay node to federating disjoint segments of WSN and proposed a distributed cell-based optimized relay node placement (CORP) algorithm [8].Marinho et al. studied the use of UAVs and cooperative multiple input multiple output (MIMO) techniques to keep the WSN connected [9].

Mobile Relay Motion Control
Problem.Similar to static relays, mobile relays were also proposed to extend the lifetime of WSN.But mobile relays often have rich resources and are able to move in the field.Thus mobile relays are more flexible than static relays and can be reused in different positions.
For WSN, the main purpose of using mobile relays is to save energy consumption of sensors and to extend the lifetime of the network.Wang et al. first studied the performance of a large dense network with one mobile relay and proposed a joint mobility and routing algorithm which can yield a network lifetime close to the upper bound [10].Venkateswaran et al. proposed a novel relay deployment framework that utilizes mobility prediction of MANET nodes to optimally define the movement of the relay nodes [11].Their simulation results indicate significant energy savings.El-Moukaddem et al. studied the problem of using mobile relays in dataintensive WSN to save the energy consumption of the whole network [12].They consider the energy consumption of both mobility and wireless transmissions.
As airborne platform of mobile relay, UAVs have been introduced to WSN, MANET, and other kinds of ground networks to improve the connectivity or link capacity between ground nodes.Since UAVs have relatively high mobility, sensibility, and self-controllability, they are especially suitable for MANET that has dynamic topology.
Chandrashekar et al. considered the problem of providing full connectivity to disconnected ground MANET nodes by dynamically placing Unmanned Aerial Vehicles (UAVs) to act as relay nodes [2].But they did not consider the mobility control problem of deployed UAVs.Jiang and Swindlehurst considered using a multiantenna UAV to connect a collection of single-antenna ground nodes [3].By dynamically adjusting the UAV heading they can maximize the approximate sum rate of the ground-to-air uplink channel.Han et al. considered using one UAV to improve the connection of ground-based wireless ad hoc networks [1].The location and movement of the UAV are optimized to improve four types of network connectivity including global message connectivity, worst case connectivity, network bisection connectivity, and -connectivity.Both of these two works just consider using a single UAV.
With the development of UAV's manufacturing technology, the size of UAV becomes smaller and the price also becomes cheaper.Thus it is possible to use a team of UAVs to provide network connection for ground nodes or improve their link capacity.Zhan et al. investigated a communication system in which UAVs are used as relays between groundbased terminals and a network base station [13].They developed an algorithm for optimizing the performance of the ground-to-relay links through control of the UAV heading angle.Cetin and Zagli studied UAVs' motion control to achieve continuous long-range communication relay infrastructure [14].They have proposed a novel dynamic approach to maintain the communication between vehicles.Besides dynamically keeping vehicles in range and appropriate position to maintain communication relay, artificial potential field based path planning also provides collision avoidance system.Ponda et al. presented a cooperative distributed planning algorithm that ensures network connectivity for a team of heterogeneous agents operating in dynamic and communication limited environments [15].The algorithm predicts the network topology and proposes relay tasks to repair connectivity violations.Dixon and Frew considered using chains of UAVs to improve the connectivity for two isolated ground nodes [4].The mobility of UAVs is controlled to maximize the communication link capacity for the end-toend connection.But they just assume that there are two nodes in ground MANETs.

System Model
We assume a scenario that multiple UAVs are used to maintain the connectivity of ground MANETs and some UAVs have already been deployed in the field.However, due to the movement of ground nodes and limited communication range, existing UAVs are not able to connect all ground nodes.Thus we need to add new UAVs to maintain the full connectivity of MANETs.The system model is shown in Figure 2. In the figure, vehicles represent ground nodes.As we can see, there are already two UAVs deployed in the field, but, due to the long distance between vehicles, the ground network is partitioned into two parts.In order to keep the full connectivity, a new UAV is added as relay that connects these two partitioned subnetworks.

Mobility Model of UAVs.
We assume the UAVs used in this paper are small four-rotor UAVs.The four-rotor UAVs can stay in a constant position and fly directly up and down.It can also spin 360 degrees around itself with a zero radius.To simplify the system, we assume all UAVs fly in different altitude so that collision avoidance of UAVs needs not to be considered in this paper.
Since the four-rotor UAV is small and always uses battery as energy, the velocity of the UAVs is limited.Because ground nodes are continuously moving, the task of motion control of existing UAVs and deployment of new UAVs must be finished within a given deadline.This requirement is especially important for some military scenarios.So we assume in this paper the mobility of existing UAVs is constrained.They can move towards any direction.However, the distance between new positions and current positions must not be more than a constant length.And we call this constant length motion range in this paper.

Communication Model.
Connectivity represents the communication capability between nodes in a network.Here, the link capacity is used to represent the connectivity.The meaning of maintaining the connectivity of ground MANETs is to keep the link capacity higher than a given threshold.When the link capacity between two nodes is higher than the threshold, these two nodes are connected.Otherwise, they are disconnected.
Link capacity is the upper bound of data rate when transmitting data.According to Shannon equation the link capacity can be computed using (1).Here,  is the link capacity in b/s;  is the bandwidth and SNR is the signal noise ratio Signal noise ratio (SNR) is the ratio of receiving power of signal with power of noise signal.The SNR of node  from node  is defined as (2) [16].Here,   is the sending power,  , is the distance between node  and node , and  is the path loss. 0 is the power density of noise and here we consider the Gaussian white noise, so  0 = 4 * 10 −21 Hz.  is the bandwidth According to (1) and ( 2), we can find that there are two factors that affect the connectivity.One factor is the distance between two nodes.The larger the distance is, the weaker the link capacity will be.When distance is large enough, the link capacity will be smaller than the given threshold and these two nodes will be disconnected.The other factor is the path loss, which reflects the signal interference of the environment.Different environment has different path loss.Signal loses more power when it is transmitted near the ground than that in the air.So when two nodes have the same distance, the link capacity between these two nodes is higher when these two nodes are in the air than that when they are on the ground.
In this paper, we define the connectivity between any two nodes as a binary variable.If the distance between two nodes is not more than a constant length, we assume these two nodes are connected.Otherwise, these two nodes are disconnected.Here we name the constant length communication range.The communication range of ground nodes is smaller than the communication range between two UAVs or between one UAV and one ground node.

Problem Formulation
In our definitions, we assume that all current positions of ground nodes and existing UAVs are known.We also assume there are no obstacles that affect the mobility of UAVs or transmissions.Our problem can be described as follows: given a set of ground nodes and a set of existing UAVs, we want to find new positions for existing UAVs and positions for new added UAVs to form tree spanning all ground nodes so that the number of new added UAVs is minimized.There are two constraints in this problem.One is the distance between new position and current position of each existing UAV that is not more than a given motion range.The other is that length of each edge in the tree is no more than a given communication range.Given the number of new added Steiner points and a topology which specifies the edges in the final tree, a bottleneck tree under this given topology and edge bound constraint can be computed in polynomial time [19] and the mobility constraint of mobile Steiner points can be checked in polynomial time.Therefore MST-EMSELB problem belongs to the class NP.It follows from this remark and Theorem 1 that we have proved Theorem 2. Detailed proofing of Theorem 1 can be found in our previously published paper in [20].

Heuristic Solution
As previously mentioned, the MST-EMSELB problem belongs to NP complete problem; thus we try to find polynomial time approximation algorithms for this problem.In this section, we firstly present existing methods and then propose three heuristic algorithms for MST-EMSELB problem.
5.1.Non-EUA Algorithm.Currently, there are no particular algorithms designed for MST-EMSELB problem.The most related problem of MST-EMSELB is MST-EMS problem.When we crystallize the MST-EMS problem in the scenarios of using UAVs to maintain the connectivity of MANETs, a specific problem of MST-EMS is using new UAVs to maintain the connectivity of MANETs without considering existing UAVs.
Lin and Xue presented a minimum spanning tree (MST) based heuristic algorithm for MST-EMS problem whose worst case approximation ratio is 4 [18].The MST heuristic algorithm firstly generates a minimum spanning tree over .It then divides each edge  in the tree into small pieces of length at most  by inserting ⌈()/⌉ − 1 degree-2 Steiner points so that all pieces in edge  have equal length.Here () is the Euclidean length of edge .
Since no mobile Steiner points are considered, the MST heuristic algorithm cannot be directly used for MST-EMSELB problem.Here we just take the Lin and Xue methods as a comparative method.Since this method has not considered the existing UAVs, we call this method non-Existing UAVs Aware (non-EUA) algorithm.Non-EUA algorithm just computes minimum number of new UAVs that is needed to connect all ground nodes.None of the existing UAVs will be reused for connecting ground MANETs.So the number of needed new UAVs computed by non-EUA should be an upper bound of other Existing UAVs Aware algorithms.The non-EUA algorithm is shown in Algorithm 1.

DBM Algorithm.
Since this is a joint optimization problem and there are two variables that need to be optimized, one variable is the new position of existing UAVs and the other variable is the position of new added UAVs.So we decouple the MST-EMSELB problem into two subproblems: the movement control problem of existing UAVs and the deployment problem of new added UAVs.In order to optimize the joint problem, we solve these two subproblems one by one.The first algorithm we proposed is Deploy-Before-Movement (DBM) algorithm that firstly optimizes the deployment of new UAVs and then optimizes the movement control of existing UAVs.
The main idea of DBM algorithm is shown as follows.Firstly we use the non-EUA algorithm to generate candidate positions of new added UAVs without considering existing UAVs.Then we match existing UAVs with candidate positions of new added UAVs.A match between an existing UAV and a candidate position of a new added UAV means that the new added UAV will be replaced by the existing UAV by moving this existing UAV to the candidate position.Since the motion range of existing UAVs is limited, the number of matches is also constrained.Here, we use Hungary algorithm to find maximum matches so that the number of needed new UAVs can be minimized.The DBM algorithm is shown in Algorithm 2.
Although DBM algorithm can utilize the mobility of existing UAVs, the candidate position for moved existing UAVs is limited which is the position computed by non-EUA algorithm.Due to the motion range limitation, existing UAVs may not be able to move to these candidate positions.Figure 3 shows one worst case for DBM algorithm.The worst case scenario is similar to the example showed in Figure 1.But the motion range of existing UAVs is constrained in this scenario.In Figure 3(a), the radius of motion range of existing UAVs is .When using DBM algorithm, the candidate position of new UAVs is the middle point of two ground nodes as shown in Figure 3(b).The distance from one existing UAV to the candidate position is  and in this case  is larger than .Thus, none of these two existing UAVs can be matched to the candidate position.So the minimum number of needed new UAVs computed by DBM algorithm is 1.
However, the best solution for this case is shown in Figure 3(c).The best movement of existing UAVs is that they all move directly towards the line that consisted of two ground nodes until meeting the motion range.A communication link will be set up for two ground nodes using two existing UAVs as relays and no new UAVs are needed in this solution.The worst case indicates that DBM algorithm cannot make best use of existing UAVs due to the less optimal motion control.

MBD Algorithm.
As previously mentioned, we decouple the joint optimization problem into two subproblems: the deployment of new UAVs problem and motion control problem of existing UAVs.DBM algorithm firstly solves the deployment problem and then solves the motion control problem.However, due to the less optimal motion control, DBM algorithm encounters some worst cases that none of existing UAVs can be reused.So we reverse the solution and propose Move-Before-Deployment (MBD) algorithm that firstly solves the movement problem and then solves the deployment problem.
The main idea of MBD algorithm is as follows.Firstly, we use a heuristic function to generate new positions of existing UAVs   .Then, we merge the set of ground nodes  and the set of existing UAVs with new positions   into a big node set  ∪   .After that, we generate a minimum spanning tree Input: , , , ,  //  is the position of ground nodes,  is the position of existing UAVs,  is the motion-range of existing UAVs,  is the communication-range between ground nodes and  is the communication-range between ground nodes and UAVs.Output:   ,  //  is the new position of existing UAVs and  is the position of new added UAVs (1) Initially:   ← ;  ← Φ; (2) generate candidate positions  of new UAVs using non-EUA(, , ) algorithm; (3) cost ← Φ; (4) for each pair   ∈  and   ∈  do (5) if distance(  ,   ) ≤  then (6) cost [𝑖][] = distance(  ,   ); (7) else (8) cost [𝑖][] = +∞; (9) end if (10) end for (11)    over set  ∪   and then existing UAVs cut process will be used to cut all 1-degree existing UAVs in the tree  until all existing UAVs in the tree have at least 2 degrees.For the rest subtree   of , new UAVs will be added to edges of   that has larger length than .
The heuristic function mentioned in previous paragraph to generate new positions of existing UAVs might affect the performance of the whole MBD algorithm.We find that the DBM algorithm will be a perfect heuristic function to generate new positions for existing UAVs if the motion constraints of UAVs are released.This is because when the motion constraints of existing UAVs are released, they can move to any positions.So all of existing UAVs can match any candidate positions of new UAVs, since these candidate positions are computed by non-EUA algorithm that can minimize the number of new UAVs to maintain the connectivity of ground MANETs without considering existing UAVs.Thus, DBM algorithm can find the best positions for existing UAVs.The MBD algorithm is shown in Algorithm 3.
In MBD algorithm, the existing UAVs cut process is an important and necessary process to minimize the number of new added UAVs, because when constructing the minimum Input: , , , ,  //  is the position of ground nodes,  is the position of existing UAVs,  is the motion-range of existing UAVs,  is the communication-range between ground nodes and  is the communication-range between ground nodes and UAVs.Output:   ,  //  is the new position of existing UAVs and  is the position of new added UAVs (1) Initially:   ← ;  ← Φ; (2) compute new positions of existing UAVs   using DBM(, , , , ) algorithm; (3) generate a complete graph (, ) over  ∪   ; (4) compute a minimum spanning tree (,   ) based on (, ); (5) while true do (6)  ← Φ,  ← Φ; // is the set of 1-degree UAVs and  is the set of edge connecting vertexes in .(7) for each V  ∈  do (8) if degree(V  ) == 1, V  ∈   then (9)  ← V  ,  ←  , ; (10) end if (11) end for (12) if sizeof() > 0 then (13)  ← ( − ),   ← (  − ); ( 14) else (15) break; (16) end if (17) end while (18)   spanning tree , all existing UAVs have been considered as terminals.However, existing UAV that has 1 degree in  is the leaf of the tree and should not be considered in adding new UAVs in next step.So we recursively cut 1-degree existing UAVs until all existing UAVs have at least 2 degrees and then the number of new added UAVs can be minimized.
The existing UAVs cut process is shown in Figure 4.There are three ground nodes and four existing UAVs in this scenario.The generated minimum spanning tree over all ground nodes and existing UAVs is shown in Figure 4(a).We can find that the degree of existing UAV  is 1.So the existing UAVs cut process will firstly delete UAV  and generate a subtree in Figure 4(b).Again, we find the degree of existing UAV  is 1.So UAV  is also cut and generates a subtree in Figure 4(c).Now, we find that there are not any 1-degree existing UAVs in the tree and all existing UAVs have at least 2 degrees.So we will add new UAVs to edges of the subtree that has larger length than  which is shown in Figure 4(d).

DAM Algorithm.
Although MBD algorithm can minimize the number of new added UAVs using the existing UAVs cut process, the new positions of existing UAVs may not be optimal.This is because the new position of existing UAVs in MBD algorithm is generated using DBM algorithm which is based on the release of motion range constraints.So, when the motion range of existing UAVs is constrained, the performance of MBD algorithm will be depressed.
Since the deployment of new UAVs and the motion control of existing UAVs affect each other, we think that if we can solve the deployment and motion control problem crosswise, then the solution of the joint problem may be optimal.So here we propose a Deploy-Across-Movement (DAM) algorithm that solves these two problems simultaneously.
The main idea of DAM algorithm is as follows.Firstly, we generate a complete graph (, ) on ground nodes and sort all edges of  in length increasing order.Then, we consider all edges  , in the set  that the length of the edge is no more than  and vertexes of the edge belong to different components.After this step, we will get several components that consist of connected ground nodes.Now we will recursively move existing UAVs and add new UAVs to connect partitioned components until all partitioned components are connected into one component.In each loop, we will try to connect all vertexes pairs V  and V  that belong to different components using two different methods.One method uses existing UAVs to set up a communication chain between V  and V  by moving UAVs to certain positions.New UAVs will be added to edges of the chain which has larger length than .The other method does not consider existing UAVs and just try to set up a communication chain between V  and V  by adding new UAVs.The number of new added UAVs using these two methods will be compared and the less one will be recorded as the minimum number of new UAVs (MNN) for connecting V  and V  .The vertexes pair that has minimum MNN will be selected to connect two partitioned components in this loop.New positions of existing UAVs and positions of new added UAVs generated to connect this vertexes pair would also be recorded as part of the final result.The DAM algorithm is shown in Algorithm 4.
The New-UAV-Chain method in DAM algorithm is just the non-EUA algorithm in two-ground node case.The Existing-UAV-Chain (EUC) method is an important part of DAM algorithm since it controls the movement of existing UAVs.The details of Existing-UAV-Chain method are shown in Algorithm 5. Given two ground nodes  1 ,  2 and all existing UAVs , it will firstly generate a minimum spanning tree  over set  1 ∪  ∪  2 and get the existing UAVs chain  from  1 to  2 .For each UAV   in the chain , it will compute a new position that depends on position of its left node   and right node   in the chain.For the first UAV, the left node is  1 and for the last UAV the right node is  2 .There are three candidate new positions for   with different priority.The first candidate position with highest priority is the middle point of nodes  1 and  2 .The second candidate position with less priority is the shadow point from   to the line that consists of nodes  1 and  2 .The last candidate position with least priority is the point of all reachable points from   , nearest to the line that consists of nodes  1 and  2 when distance between   and the line is larger than the motion range .The new position of UAV   will be set to one of these three candidate positions with highest priority under motion range constraint.Afterwards, new UAVs will be added to edges in chain  that has larger length than .
Figure 5 shows the three candidate positions in EUC method.There are two ground nodes  1 ,  2 and five existing UAVs in this scenario.The minimum spanning tree generated is shown in Figure 5(a).The UAV chain from  1 to  2 is  →  → .Then EUC method will compute candidate Due to the motion range of , the new position it will choose is  3 in this situation.
5.5.Discussion.In the proposed algorithms, we do not consider the power consumption of UAVs.In fact, the battery power of each UAV is limited and the limited power will certainly affect the deployment and movement of UAVs.
Here, the power consumption is modeled as a constant factor.Existing UAVs that have very low battery power should land and recharge; thus it will be excluded from the set of existing UAVs.So all existing UAVs have adequate battery power.
The battery power of UAVs can also be modeled as dynamic factor; thus existing UAVs may have different battery power due to duration time of their task.Then the deployment and movement UAVs should consider the power cost and optimize the survival time of the whole network.However, that is beyond the focus of this paper and we may consider it in our future work.

Simulation Experiment
A simulation environment was set up to test the performance of the proposed algorithm.Simulation parameters are shown in Table 1.
To demonstrate the effectiveness of proposed Existing UAVs Aware (EUA) algorithms, we compare the performance in terms of new added UAV numbers with non-EUA algorithm and the CBBA (Consensus-Based Bundle Algorithm) proposed by Ponda et al. [15].Different from non-EUA algorithm, CBBA is an algorithm that considers existing UAVs in deploying new UAVs.It uses distributed planning strategies to control the movement of existing UAVs so that network connectivity can be ensured.We also compare performance among three proposed EUA algorithms: the DBM, MBD, and DAM algorithms.We change the testing scenario by varying 5 simulation parameters including the field size, number of ground nodes, number of existing UAVs, communication range, and motion range.For each scenario, we carry out the simulation on 100 randomly generated topologies and take the average performance as the final performance for this scenario.Our simulation scenarios can be classified into five cases.We will analyze our simulation results one by one.
6.1.Varying Field Size.We use a square as testing field.We set the number of ground nodes as 50, existing UAVs as 5, communication range as 500 m, and motion range as 50 m.In this scenario, we vary the edge of testing field from 1 km to 10 km by increments of 1 km and randomly generate 100 different topologies for both ground nodes and existing UAVs.
Figure 6 shows the average number of new added UAVs to connect all ground nodes.It is obvious that the number of new added UAVs computed using the EUA algorithms is smaller than the one using non-EUA algorithm.And it is possible to observe an increase in the number of new added UAVs with the increasing size of field.This is because the larger field size is the larger distance between ground nodes since the number of ground nodes is constant.So more relays are needed to maintain the connectivity of ground MANETs.The DAM method has best performance, while DBM method, MBD method, and CBBA method have similar performance.

Varying Number of Ground Nodes.
In this scenario, we set the edge of test field as 5 km and keep other three parameters the same while varying the number of ground nodes from 10 to 100.
Figure 7 shows the average number of new added UAVs to connect all ground nodes.It is obvious that the number of new added UAVs computed using the EUA algorithms is also smaller than the one using non-EUA algorithm.And the performance of DAM algorithm is always the best one of all algorithms.At the beginning, the DBM algorithm is little better than MBD and CBBA algorithm, but, with the increase of ground nodes number, the performance of MBD algorithm is better than DBM algorithm.There is a trend in the simulation result that the number of new added UAVs firstly increases but after it reaches a peak, it decreases very quickly, because the increasing number of ground nodes requires more UAVs to maintain its connectivity at the beginning.But afterwards, sufficient number of ground nodes will increase the connectivity of MANETs so less UAVs are needed.

Varying Number of Existing UAVs.
In this case, we set number of ground nodes as 50 and kept the field size, number of ground nodes, communication range, and motion range static.We vary the number of existing UAVs from 2 to 20, by increments of 2.
Figure 8 shows the average number of new added UAVs when given different number of existing UAVs.As we expected, the number of new added UAVs computed using the EUA method is smaller than the one using non-EUA method in all scenarios.Another observation can also emulate the importance of considering existing UAVs.The observation is that there is an obvious decrease in the number of new added UAVs with the increase of the number of existing UAVs when using EUA algorithms, while the number of new added UAVs using non-EUA algorithm almost keeps the same.We can also find that CBBA and MBD algorithm are better than DBM algorithm.DAM algorithm always has the best performance.From Figure 8, we can see that the MBD algorithm can averagely reduce about 30% of new UAVs number compared with CBBA and DBM algorithm.DAM algorithm can averagely reduce about 70% of new UAVs number compared with DBM algorithm.

Varying Motion Range.
In this case, we set existing UAVs as 5 and keep the field size, number of ground nodes, number of existing UAVs, and communication range static.We vary the motion range from 10 m to 100 m, by increments of 10 m, and generate 100 different topologies for each motion range.Figure 9 shows the average number of new added UAVs under different motion range of existing UAVs.EUA algorithms have better performance than non-EUA algorithms and DAM algorithm has best performance.CBBA method firstly has better performance than DBM and MBD method, but when motion range is larger than 50 m the performance of CBBA is worse than BDM and MBD method.The observation of BDM and MBD method is similar to results under varying number of existing UAVs.However, the difference between DBM algorithm and MBD algorithm is decreasing with the increase of motion range constraint.So when the motion range increases, the difference between DBM and DBM algorithm will vanish.

Varying Communication Range.
In this case, we set motion range as 50 m while keeping the number of ground nodes, number of existing UAVs, and motion range static and vary the communication range of ground nodes from 50 m to 100 m, by increments of 50 m.The communication range between ground nodes and UAVs is set as twice communication range of ground nodes.
Figure 10 shows the average number of new added UAVs under different communication ranges.We can observe that, for all algorithms, the number of new UAVs is decreasing with the increase of communication range.This is due to the fact that larger communication range will lead to better network connectivity.Thus fewer relays are needed to maintain the connectivity of the network.So when the communication range increases to certain value, the difference between all algorithms will vanish.The DAM method has the best performance.CBBA, DBM, and MBD methods have similar performance.

Conclusion
This paper studies the problem of using UAVs to maintain the connectivity of ground MANETs.Different from existing works, this paper considered a condition that some UAVs have already been deployed in the field.Due to the movement of ground MANETs and limited communication range, existing UAVs are not able to connect all ground nodes, so new UAVs need to be deployed to maintain the connectivity.

2 MathematicalFigure 1 :
Figure 1: An example that illustrates the importance of existing UAVs in maintaining the connectivity of ground MANETs.

Figure 4 :Figure 5 :
Figure 4: Existing UAVs cut process and new UAVs adding process of MBD algorithm.

Figure 6 :
Figure 6: Comparison with varying field size.

Figure 7 :
Figure 7: Comparison with varying numbers of ground nodes.

Figure 8 :
Figure 8: Comparison with varying numbers of existing UAVs.

Figure 9 :
Figure 9: Comparison with varying motion range of existing UAVs.

MathematicalFigure 10 :
Figure9shows the average number of new added UAVs under different motion range of existing UAVs.EUA algorithms have better performance than non-EUA algorithms and DAM algorithm has best performance.CBBA method firstly has better performance than DBM and MBD method, but when motion range is larger than 50 m the performance of CBBA is worse than BDM and MBD method.The observation of BDM and MBD method is similar to results under varying number of existing UAVs.However, the difference between DBM algorithm and MBD algorithm is decreasing with the increase of motion range and is almost equal to each other at the right end of the figure.This is because MBD algorithm uses DBM algorithm to generate new positions of existing UAVs.And it is based on the release of motion Since this problem is similar to Steiner Tree Problem with minimum number of Steiner points, we formulate this problem as a Minimum Steiner Tree problem with Existing Mobile Steiner points under Edge Length Bound constraints (MST-EMSELB).The Steiner points here stand for UAVs and the Edge Length Bound is the communication range.The formal definition of MST-EMSELB problem is shown as follows.Given.There are a set of ground nodes , a set of existing UAVs , a motion range , ground node communication range , and ground-air communication range .Here,  < ,  = { 1 ,  2 , . . .,   } ,  = { 1 ,  2 , . . .,   } .Output.There are new positions of exiting UAVs , positions of new added UAVs , and a tree  spanning : .2.2.Decision Version of MST-EMSELB Problem.Given a set  of  terminal points and a set  of  mobile Steiner points in the two-dimensional Euclidean plane R 2 , a nonnegative constant , a positive constant , and a nonnegative integer , the MST-EMSELB problem firstly asks whether there exists a motion of mobile Steiner points in set  such that distance between new position and old position of each mobile Steiner point is no more than .Then, the MST-EMSELB problem will check if there exists a tree spanning a point set  ⊇  such that the length of each edge in the tree is no more than  and the number of new added Steiner points (points in  ∩ ¬( ∪ )) is no more than .

Table 1 :
Simulation parameter.UAVs in the chain one by one.Take the first UAV  as an example.The three candidate positions for  are  1 ,  2 , and  3 shown in Figure5(b). 1 is the middle point of  1 and  2 . 2 is the shadow point of  to lines  1 ,  2 .And  3 is the point in the motion range of that is nearest to lines  1 ,  2 .The priority of these three candidates' position is  1 >  2 >  3 .