Multi-User AF Relay Networks with Power Allocation and Transfer: A Joint Approach

: The Internet-of-Things (IoT) framework has been considered as an enabler of the smart world where all devices will be deployed with extra-sensory power in order to sense the world as well as communicate with other sensor nodes. As a result, smart devices require more energy. Therefore, energy harvesting (EH) and wireless power transfer (WPT) emerge as a remedy for relieving the battery limitations of wireless devices. In this work, we consider a multi-user amplify-and-forward (AF)-assisted network, wherein multiple source nodes communicate with destination nodes with the help of a relay node. All the source nodes and the relay node have the capability of EH. In addition, to cope with a single point of failure i.e., failure of the relay node due to the lack of transmit power, we consider the WPT from the source nodes to the relay node. For WPT, a dedicated energy control channel is utilized by the source nodes. To maximize the sum rate using a deadline, we adopt a joint approach of power allocation and WPT and formulate an optimization problem under the constraints of the battery as well as energy causality. The formulated problem is non-convex and intractable. In order to make the problem solvable, we utilize a successive convex approximation method. Furthermore, an iterative algorithm based on the dual decomposition technique is investigated to get the optimal power allocation and transfer. Numerical examples are used to illustrate the performance of the proposed iterative algorithm.


Introduction
In today's smart world, wireless communication services are growing in demand exponentially. The Internet-of-Things (IoT) framework has been viewed as an enabler of the smart world. The key requirements to enable IoT services are prolonging the lifetime of operating nodes and extending the coverage range [1,2]. Since relaying schemes are paramount for improving coverage, capacity, and reliability of wireless networks, the lifetime of wireless networks can be enhanced by deploying intermediate energy harvesting (EH) relay nodes, thus, we don't need to replace the battery sporadically [3][4][5].
As a result of the aforementioned advantages, EH cooperative communication can easily become the backbone of modern, green, smart wireless networks. However, the design challenges of EH wireless networks lie in the fact that solar, wind, and thermoelectric energy, which are common environmental energy sources, are intermittent and random in nature. These common energy sources are not always available. Hence, for transmission reliability, these uncertainties become more challenging for designing EH communication systems. In such a scenario, a promising solution has emerged as regards the wireless power transfer (WPT) in an EH relay-assisted network, wherein In this work, we consider a multi-user AF assisted network and focus on maximizing the sum rate of the designed framework by jointly optimizing optimal power allocation and transfer. We assume that each transmitter (sources/relay) has the capability of EH and that SWIPT takes place from the source nodes to the relay node only.
In order to improve the reliability of the relay node, the source nodes use a dedicated energy control channel to transfer the harvested energy to the relay. Notice that the dedicated energy control channel is assumed to be orthogonal to the data transfer channel. Since the relay node operates in AF mode, the relay-assisted network suffers from a noise magnification problem. This presents an interesting question: how can we balance energy disbursement and energy transfer between EH nodes. In respect to this, we formulate an optimization problem for maximizing the sum rate by considering both energy and battery causality constraints. The formulated problem is non-convex and intractable in nature. In order to make the problem solvable, we utilize a successive convex approximation method. Furthermore, an iterative algorithm based on the dual decomposition technique is proposed to achieve the optimal strategy.
The rest of this work is divided into the following sections. The system model, EH and WPT model, and the formulation of an optimization problem are described in Section 2. The proposed strategy is illustrated in Section 3. The numerical results are given in Section 4. We conclude the paper in Section 5.

System Model
Consider an EH dual-hop AF-assisted relay network consisting of N EH source nodes (S k , for k = 1, . . . , N), one relay node (R) and N destination nodes (D k , for k = 1, . . . , N), as depicted in Figure 1. A single antenna is deployed at each node in the network. As a result of long distance and deep fading, there is no direct link between S k and D k , ∀k. The nodes S k , ∀k, and R can harvest energy from ambient energy sources and store them in batteries. We adopt a time-slotted model with a unit time duration for each time slot. The total time slot length is N T . Furthermore, the battery capacity i.e., energy queues, at each S k , ∀k, and R are B s k,max and B r,max units of energy, respectively. At the end of each slot, the status of batteries is updated. We further assume frequency flat fading channels between all links. In each time slot, the relay R works in half-duplex (HD) fashion. Suppose i , for i = 1, 2, . . . , N T , indicates the time slot index. In the first phase, the signal received by the relay node in the i th time slot is written as where P k,i represents the transmission power of the node S k , s k,i is the transmitted message with E |s k,i | 2 = 1, h k,i is the channel coefficient from S k node to R, n r,i ∈ C 0, σ 2 r is AWGN at node R. The signal transmitted by node R after amplification is given by , and W r,i is the transmission power of node R. In the second phase, the received signal at the j th destination node in the i th time slot is expressed as y j,i = α i P j,i g j,i h j,i s j,i + N ∑ k=1,k =j P k,i g j,i h k,i s k,i + g j,i n r,i + n d j ,i , where g j,i , n d j ,i and σ 2 d j are defined similar to h k,i , n r,i and σ 2 r for the link from R to the D j .
S N Relay (R)   Using (2), we can compute the sum rate over N T consecutive time slots as where P s ∈ C N×N T , W r ∈ C 1×N T , and signal-to-interference-plus-noise ratio (SINR) at node D j is obtained as with

Energy Harvesting and Power Transfer Model
At the i th time slot, E k,i ≥ 0, k = 1 . . . N and E r,i ≥ 0 are the amounts of harvested energy at S k and R nodes, respectively, and we assume the energy profiles are non-casually known before data transmission [14]. Furthermore, during the i th slot, each source node S k , transfers δ k,i ≥ 0 amount of energy to R with an efficiency of β k , 0 ≤ β k ≤ 1 by utilizing a dedicated control channel with channel coefficients f k,i . Thus, when node S k has enough energy, it can share a part of its harvested energy with node R to extend the lifetime of the relay-assisted network. The power allocation and transfer depends on the following two energy causality constraints: whereδ k,i = f k,i 2 δ k,i , and f k,i indicates the energy control channel from node S k to node R at the i th time slot. Furthermore, the battery storage constraints are required to control the battery overflow at S k , ∀k and R nodes, respectively, and are given by where P c k,i and P c r,i are the fixed processing power consumed at the k th source for energy transfer and at the relay node in energy conversion, respectively.

Subproblem Solution
The values of power allocation and transfer at the i th time slot for the S k node and node R can be computed iteratively using the Karush-Kuhn-Tucker (K.K.T.) conditions. TheP are obtained by taking the partial derivative with respect toP k,i ,W r,i , and δ k,i and by equating the result to zero. The detailed derivation which is straight forward is omitted here.
With the updated multipliers, we computeP . Next, ρ i and β i using (26) and (27) are recomputed in order to further enhance the lower bound performance in (24). We repeat this procedure until convergence of the algorithm. The proposed algorithm is outlined in Algorithm 1 for the fixed values of Ω and Ω r . Algorithm 1: Proposed algorithm for given values of Ω s and Ω r .
Theorem 1. The optimal P s ,W r , δ of the problem (P4) must fulfill the following two conditions: Proof. This theorem can be proved by contradiction. Let P s ,W r , δ , Ω s , Ω r be the optimal solution of (P4). First suppose ∑ , Ω s , Ω r . Then, we can ameliorate eP k,T , ∀k to achieve an increased sum rate off R (P s ,W r ), directly implied by utilizing the SINR obtained in (4), without conflicting with any other constraints in P4. Thus, (T.1) contradicts with the optimality of P s ,W r , δ , Ω s , Ω r .
k=1 β kδ k,i is fulfilled with strict inequality, for P s ,W r , δ , Ω s , Ω r . By reducingδ k,T , ∀k and increasing eP k,T and Ω k,T−1 , we can get an ameliorated sum rate, while fulfilling other conditions in (P4). Thus, (T.2) directly contradicts with the optimality of P s ,W r , δ , Ω s , Ω r . Hence, it is proved. Now, we update the values of Ω s and Ω r with the help of the following lemma.

Lemma 1.
To obtain the maximum sum rate, the auxiliary variables Ω s and Ω r should be updated as Proof. To maximize the sum rate, the auxiliary variables Ω s and Ω r in the outer loop can be illustrated as [Ω s , Ω r ] = arg max wheref LB (Ω s , Ω r ) denotes the sum rate achieved by the Algorithm 1. It is very hard to exhaustively search over all values of Ω s and Ω r due to the very high computational complexity. Therefore, we propose a two-step technique to find the values of Ω s and Ω r . As the first step, relax the battery storage constraints in (P.4) by assuming Ω r,m = ∞, for m = 1, . . . , N T − 1, and k = 1, . . . , N, and find the optimal P (0) s ,W (0) r , δ (0) using the Algorithm 1. In the next step, we recompute the optimal solution by updating Ω

Numerical Results
The path loss models given by 25.17 + 20 × log 10 (d) dB (d: distance in km) [27,28] are adopted for the data transfer and the dedicated power transfer channels from S k , ∀k to R, while the 3GPP path loss model, given by 131.1 + 42.8 × log 10 (d) dB, is used for the channels from R to D k , ∀k. The distance from S k , ∀k to R is considered to be 1 m [28,29]. The thermal noise density and the channel bandwidth are considered as −174 dBm/Hz and 1 MHz. The value of interference suppression factor is set to 0.01. E L = [1,15] mJ and E H = [1, 100] mJ are uniformly distributed EH profiles considered for the simulation. For simplicity, we assume that the EH profiles of all the source nodes are E s , whilst the EH profiles for node R is E r . The battery capacity at each S k , ∀k, and R nodes is assumed to be E s,max = E r,max = 10 mJ. We consider two scenarios: (1) E s = E r ; (2) E s ≥ E r , for T = 15. Furthermore, we set α = 0.7. The distances from S k , ∀k to R and from R to D k , ∀k are indicated by d sr and d rd , respectively. In addition, we set L max = 10 and i = 0.01. For better performance comparison, we also simulated the without WPT transmission. Figure 2 validates the convergence behavior of the algorithm for E s = E H and E r = E L . Note that Figure 2 is plotted for a single channel realization. From this figure, we find that the sum rate increases monotonically when the number of iterations increases. In addition, it is noticed that the algorithm converges within four iterations which reflects the effectiveness of the algorithm.   Figures 3 and 4 illustrate the throughput performance of the network in terms of sum rate for each time index in scenario 1 and scenario 2. In this example, we also include the average sum rate of the network without WPT. In both scenarios, it can be observed that the sum rate performance of the algorithm with energy transfer is much better than that without energy transfer. It can also be observed in Figure 4 that for better EH conditions at the S k , ∀k nodes compared to the condition at the R node, the sum rate performance is remarkably enhanced at the end of the time.    Figure 5 shows the achievable sum rate performance of the system when E s = E L < E r = E H for d sr = 1 m and d rd = 50 m. Since the EH rate at the source nodes is the worst, the effect of power transfer from the source nodes to the relay node is almost zero and therefore the performance of the system with WPT is identical to that without WPT. Figure 6 illustrates the impact of the relay's position on the average sum rate performance. The parameter settings are as follows: We set N = 1 and d sr + d rd = 2 m. The distance from S k to R and from R to D k are represented b d SR and d RD , respectively, and a distance ratio is defined as r d = d SR /(d SR + d RD ). From Figure 6, we can observe that when d sr + d rd = 2 m, the average sum rate performance is dominated by the energy transfer rate, δ i ; ∀i, and the highest sum rate performance can be achieved when the relay node is placed closer to the source node.

Conclusions
We studied a joint design of power allocation and transfer for multi-user AF-assisted networks with EH. We focused on maximizing the achievable sum rate using a deadline, subject to the battery and energy causality constraints at the S k , ∀k, and R nodes, and proposed an iterative algorithm for obtaining the near-optimal solution using successive convex approximation and dual decomposition techniques. Through numerical simulation, we found that the multi-user AF-assisted network with energy transfer achieves a significant sum rate enhancement over the one without WPT.