Abstract

The problem of energy shortage in sensor nodes caused by frequent data interactions is one of the major constraints on the development of twin workshops. A backscatter-assisted wireless powered communication network (BAWPCN) has been deemed a potential solution for addressing the problem of energy shortage in twin workshops. How to effectively ensure the link energy efficiency (EE) while satisfying the quality of servers for each user has been of high interest, while it has not been well studied in previous works. Inspired by this, in this paper, we propose a resource allocation scheme based on the max–min criterion, considering the user quality of service and energy-causality constraints. The optimization problem is formulated as a mixed-integer nonconvex fractional planning problem which is aimed at maximizing the minimum EE of each link. The generalized fractional theory is used to transform the nonconvex fractional planning problem into an equivalent mixed-integer nonconvex subtraction optimization problem, and then, the mixed-integer nonconvex subtractive optimization problem is transformed into an equivalent nonconvex optimization problem by introducing relaxation variables to eliminate the integer programming problem arising from the maximum-minimum function. Based on this, the block coordinate descent method is used to decompose the transformation problem into two convex subproblems, and an iterative algorithm is proposed to solve the transformation problem. Simulation results verify the fast convergence of the proposed iterative algorithm and show that the proposed resource allocation method can effectively guarantee the fairness of the energy efficiency of the system in twin workshops.

1. Introduction

Industry 4.0 refers to using the Cyberphysical System (CPS) to digitize and intellectualize the supply, manufacturing, and sales information in production and finally achieves fast, effective, and personalized product supply. The core of Industry 4.0 is the interconnection and intelligent operation of the physical and information worlds of manufacturing. The combination of intelligent automation (e.g., robotics) and a new generation of information technology (e.g., Internet of Things and artificial intelligence) produces the digital twin workshop, which is one of the current trends in manufacturing in the context of Industry 4.0 [1]. Digital twin (DT) is driven by the multidimensional virtual model and fused data and realizes monitoring, simulation, prediction, optimization, and other actual functional services and application requirements through virtual and real closed-loop interaction [2], as shown in Figure 1. Digital twin workshops (DTW), as an important enabling way to realize digital transformation, promote intelligent upgrading, and accelerate Industry 4.0, have moved from theoretical research to the practical application stage, among which the development of Internet of Things technology is the key driving force.

With the development of Internet of Things (IoT) technology and its wide application in the manufacturing industry, it greatly enhances the real-time data acquisition ability and the transmission ability of production factors in twin workshops. The realization of ubiquitous IoT requires the deployment of numerous low-power sensors, but the frequent data interaction will greatly consume sensor energy [4]. One way to solve the problem is replacing the battery, but the deployment of sensors in an industrial production environment cannot meet the demand to frequently replace the battery. At the same time, as the power consumption and volume of wireless sensor nodes become smaller and smaller, the advantages of small and flexible wireless sensor nodes will be limited if they only rely on their own power supply. Therefore, the energy shortage of sensor nodes caused by frequent data interaction is one of the important factors restricting the development of twin workshops. In recent years, thanks to the progress of science and technology, scholars have proposed wireless power communication network (WPCN) and backscatter communication technologies to solve the energy limitation problem of sensor nodes. The WPCN deploys dedicated energy stations to provide the energy resource to the sensor nodes through the wireless energy transfer (WET) technology, and the nodes leverage the harvested energy to transmit information. Therefore, the core of WPCN design lies in the joint allocation of energy and time resource. The authors in [5] proposed the harvest-then-transmit (HTT) protocol by incorporating WET into wireless communication networks, which is the main basis for the operation of current WPCNs. In [6], the half-duplex mode is extended to the full-duplex mode, and the base station contains two antennas so that it can simultaneously transmit energy signals and receive transmission data. However, the transmitter performs active transmission (AT) by the HTT protocol, which leads to high power consumptions. Thus, the transmitter has to allocate a long period to harvest sufficient energy but leaves a short period to transmit information. Backscatter communication (BackCom) technology adopts a relatively simple modulation method to load its information to external RF signals for transmission. Thus, BackCom avoids power-consuming components, providing a transmission strategy with low power consumption and low transmission rate for IoT scenarios [7]. The long period for HTT to harvest energy is available for BackCom to increase the throughput. Thus, BackCom can be combined with HTT to achieve different tradeoffs between energy harvesting (EH) and data transmission.

The authors of [8] take the wireless digital TV signal as an example to show that the sensor can maintain communication by using the existing wireless signal without batteries. In [9, 10], the authors combined BackCom with WPCN and proposed a backscatter-assisted wireless powered communication network (BAWPCN). For the BAWPCNs, the authors of [11] considered a network combining multimode backscatter communication and HTT and proposed a time slot resource allocation scheme to maximize subuser link capacity. In [12], the authors studied the optimal time allocation of energy harvesting, backscatter, and wireless transmission to maximize the throughput. The authors of [13] proposed a resource allocation scheme to maximize the throughput by jointly optimizing the transmission time, reflection coefficient, and transmit power in the full-duplex BackCom. The authors of [14] considered the scenario of nonorthogonal multiple access and maximized the throughput by jointly optimizing the backscatter time and reflection coefficient, subject to the constraints of energy collection threshold and signal-to-noise ratio. Energy efficiency (EE) is also one of the important indicators of wireless communication networks; the authors in [15] studied the resource allocation method to maximize the user EE in the BAWPCNs. Then, users’ EE was also studied in an unmanned aerial vehicle-based WPCN with backscatter communications [16]. Based on the analysis of existing literature, it can be seen that there are many researches on BAWPCNs from the evaluation of network outage capacity or spectral efficiency, while there are few researches from the perspective of energy efficiency. Meanwhile, the following three problems still exist in the research process: (1) there are only single nodes in BackCom; (2) only the time dimension is considered in the design of resource allocation, while the power resource allocation is ignored; and (3) the energy harvesting model adopts the linear mode. However, in the case of a multiuser network in the twin workshop, the communication network composed of multiple nodes should be considered. A large number of nonlinear components such as photoelectric sensors and Hall sensors are used in twin workshops, and the actual energy harvesting circuit presents nonlinear characteristics. At the same time, power resources also affect the level of energy efficiency. Inspired by the above three problems, aiming at the energy shortage caused by frequent data exchange of each sensor node in the twin workshop, this paper introduces a BAWPCN to study the time-power two-dimension resource allocation method for multiple users and optimize the energy efficiency of the communication link. Our main contributions are summarized as follows. (i)A joint energy efficiency/fairness time-power two-dimensional resource optimization model for the system based on the max–min criterion is proposed. The proposed optimization model not only ensures fair access to communication resources for users but also considers the optimization of the time slot, transmit power of the power beacon, and nodes simultaneously(ii)The joint optimization of time slot and power leads to coupling of multiple variables and a fractional form of the objective function. In addition, there is a max–min function of the objective function. The proposed optimization model is therefore a mixed-integer nonconvex fractional programming problem and cannot be solved directly using the existing tools such as CVX(iii)The proposed iterative algorithm is verified to converge quickly through simulation and is shown to achieve better fairness of the link EE when compared with similar algorithms

The rest of this paper is organized as follows. In Section 2, the system model for BAWPCN is introduced. In Section 3, we formulate the problem to maximize the minimum links’ EE and solve it for the EE fairness resource allocation scheme. In Section 4, the simulations are conducted to evaluate the performance of the proposed scheme. Section 5 concludes this paper.

2. System Model

Consider a BAWPCN as shown in Figure 2, the network consists of a gateway (GW) to receive information, a dedicated power beacon (PB) to provide energy, and power-limited users (EUs).

The EUs need to upload their data to the gateway within time slot . The PB was deployed to provide radio frequency (RF) signals to EUs who can use the received RF signals for backscatter communication and energy harvesting.

In order to avoid interference between EUs, time division multiple access is adopted to decompose the whole time slot into multiple small time slots, as shown in Figure 2. In , PB broadcasts the unmodulated RF signals, and all EUs can use the received RF signals for energy harvesting or backscatter information. Specifically, in , all EUs work in the energy harvesting mode. At the time slot (), EU transmits data to the GW through backscatter technology while the other EUs continue to harvest energy. In , the PB remains silent and EU transmits its data to the GW in slot by active communication. The signal received by the th EU in can be expressed as where represents the transmitted power of the PB, denotes RF signals transmitted by PB and , represents the receiving noise of th EU and follows a Gaussian distribution with the mean of 0 and variance of , and represents the channel coefficient from PB to th EU. Using the nonlinear energy harvesting model proposed in [17], the energy harvested by EU within time slot can be expressed as follows: where is the maximum harvested power of the energy collector, and are parameters of the nonlinear energy model, and their values can be obtained by fitting the actual measured data. Since the value of harvested energy is positive, the value of must meet the requirement . It should be noted that the time for EU to harvest energy is , not , because EU operates in the backscattering communication mode within the time slot. It should be pointed out that the energy collected in formula (2) will be used for th EU’s backscatter communication and HTT energy dissipation in and time slots. In , EU backscatters information to the gateway via the backscattering communication technique, where the instantaneous power of the reflected signal received by the gateway from th EU can be expressed as where is the scattering efficiency of the backscattering communication module, denotes the channel coefficient from th EU to the gateway, and and denote the reflection coefficient. The Fries transfer formula is used to model the channel gain, , , where denotes the wavelength; , , denote the antenna gain of the PB, the EUs, and the gateway, respectively; and and denote the distance from the PB to th EU and th EU to the gateway, respectively. According to Equation (3) and Shannon’s capacity theorem, the achievable throughput of th EU at can be calculated as follows: where denotes bandwidth. Since backscatter communication uses simple modulation, its channel capacity is smaller than that of conventional communication modes. In this paper, we use the same approach as in [18] to characterize this difference in channel capacity, i.e., multiplying the signal reception signal-to-noise ratio by a nonnegative real number with a factor less than 1 (). In , th EU transmits data to the gateway in time slot in the traditional communication mode, so the throughput that th EU can accomplish is expressed as where denotes the transmit power of th EU during time slot . Thus, the total throughput by th EU during the entire time slot can be expressed as

According to Equation (18) in [19], the EU EE is the ratio of the total throughput achieved by the user to the consumed energy. Equation (6) already gives the total throughput that can be completed by user throughout the time slot. Next, this section analyses the energy consumed by th EU throughout the time slot and the expression of EE.

During time slots and , th EU collects energy from a dedicated energy station. During time slots and , the th EU needs to consume energy. During time slot , there is only circuit loss as the th EU modulates its own information to the received RF signal (i.e., the RF signal emitted by the energy source) and does not need to generate its own carrier. During time slot , the th EU uses a conventional communication mode to transmit data, so its energy consumption consists of two components: the energy consumed by the transmit power and the circuit losses. In summary, the total system energy consumed by th EU is where denotes the circuit loss when th EU operates in the backscatter communication mode, denotes the circuit loss when th EU is operating in the conventional communication mode, and denotes the transmit power of th EU in time slot .

Combining Equations (6) and (7), the link EE of th EU can be expressed as

3. Resource Allocation Method for Energy Efficiency Fairness of the System

3.1. Problem Formulation

In this subsection, we propose a resource allocation scheme that can effectively guarantee the EE fairness of communication links. The max–min criterion is an effective way to guarantee user fairness [20], so that the resource allocation method for guaranteeing energy-efficient fairness of links among EUs can be achieved by solving the following optimization problem, namely,

C3 ensures that the harvested energy is greater than the energy consumed by user . C4 constrains the maximum transmit power of the th EU when working in the traditional information transmission mode. C5 indicates that the throughput achieved by th EU in the whole time slot cannot be less than the given minimum value, i.e., the communication service quality of user is guaranteed. C6 constrains the maximum transmit power of the PB.

As can be seen from the objective function, the EE of th link is a fraction function, which contains multiple coupled variables, such as the variables and , , , and . Secondly, since the max–min criterion is adopted, the integer variable needs to be optimized. Therefore, the optimization problem (9) is a mixed-integer nonconvex fractional programming optimization problem. Therefore, the existing optimization tools such as CVX cannot be directly used to solve the original problem.

3.2. Problem Transformation

For the mixed-integer nonconvex fractional programming problem in Equation (9), we divide three steps to obtain its optimal solution. Firstly, the nonconvex fractional programming problem is transformed into an equivalent mixed-integer nonconvex subtractive optimization problem by using the generalized fractional programming theory. Secondly, the integer programming problem caused by the max–min function is eliminated by introducing relaxation variables; that is, the mixed-integer nonconvex subtraction optimization problem is transformed into an equivalent nonconvex optimization problem. Finally, the block coordinate descent (BCD) technology was used to decompose the transformation problem into two convex subproblems [21], and then, an iterative algorithm was designed to solve the transformation problem.

The BCD algorithm encompasses a wider range of traditional alternating optimization algorithms as well as coordinate descent algorithms. The BCD algorithm is usually applied to optimization problems with nonconvex objective functions or feasible domains, but when the optimization variables are blocked, the objective functions and feasible domains of such optimization problems are convex with respect to the blocks of variables. In addition, the scope of the BCD algorithm can be extended, for example, by using approximate optimization for nonconvex objective functions and feasible domains after blocking the variables, where the BCD algorithm is also feasible [22].

The specific processing steps are described below. (1)Let the variable denote the value of the objective function of problem (9), i.e., the max–min EE. According to the generalized fractional programming theory, the sufficient condition for obtaining the optimal solution to the optimization problem (9) is that Equation (10) holds when the C1-C7 constraints are satisfied:

We can obtain the optimal solution to the original problem by solving the optimization problem (10). However, in practice, is often unknown, but we can obtain by continuously updating the value of , which is shown in Algorithm 1.

1: Solve the optimal solution of the optimization problem given any greater than zero
2: Substitute the optimal solution obtained in the first step into the objective function of the optimization problem (9) to update
3: If does not converge, let and repeat the first step; if not, let , and the optimal solution in the first step is the optimal solution of the original problem (9):

According to Algorithm 1, the core of solving the original problem (9) is to solve the optimization problem. Compared to the original problem (9), the objective function of optimization problem does not have a fractional form, but it is still a mixed-integer nonconvex optimization problem. (2)A relaxation variable is introduced to transform the mixed-integer nonconvex optimization problem into the following optimization problem, namely,

In the optimization problem (11), the optimization objective is a linear function and the constraints C1, C2, C4, C6, and C7 are linear constraints, but the remaining constraints C3, C5, and C8 are nonconvex, so the optimization problem (11) is a nonconvex problem. (3)The optimization problem (11) is transformed into two subproblems using the BCD technique as follows: (1)Given , solve for , , , , by the following problem:

In the optimization problem (12), constraints C3, C5, and C8 still exist in the coupled variable and are jointly nonconvex. Specifically, in C3 and C8, the variable and are coupled; in C5 and C8, the joint nonconvexity is present in . To solve the above problem, we introduce auxiliary variable and bring them into the optimization problem (13) to obtain the equivalent optimization problem as follows, i.e., (2)Given , , , , , we can obtain by solving the following optimization problem:

Lemma 1. The optimization problem (13) is a convex problem.

Proof. Please see the appendix.
According to Lemma 1, we can use the CVX tool to obtain the optimal solution to solve the optimization problem (13).

In the optimization problem (14), the objective function is linear and C1, C2, C4, and C7 are independent of the optimization variables. C6 is a linear constraint. Meanwhile, C3 and C5 are convex constraints. The left-hand side of C8 is shaped as , where , , , and are all constants and greater than zero. The second-order derivative of can be calculated as . , so C8 is a convex constraint. Therefore, the optimization problem (14) is convex, and we can use the CVX tool to obtain the optimal solution to the optimization problem (14).

Based on the above analysis, the optimization problem can be solved by an iterative BCD-based algorithm, as described in Algorithm 2.

1: Initialize the system parameters, given , the convergence accuracy , and the maximum iteration number
2: Repeat
3: For given , solve the problem (13) by using CVX to obtain , , , ,
4: For given , , , , , solve the problem (14) by using CVX to obtain
5: Until or
6: Return the optimal solutions: , , , , ,

As the optimization problems (13) and (14) are convex, the convergence of Algorithm 2 is guaranteed [23]. The overall algorithm of the link energy efficiency fairness resource allocation scheme can be described as Algorithm 3.

1: Given any greater than zero, use Algorithm 2 to obtain an optimal solution to the optimization problem (14)
2: Update by bringing the optimal solution obtained in the first step into the objective function of the optimization problem (9)
3: If does not converge, let and repeat the first step; if not, let , and the solution obtained in the first step is the solution to the original problem (9)

4. Simulation Results and Analysis

In this section, we provide simulation results to evaluate the performance of the proposed iterative algorithm. In addition, we illustrate the advantages of the proposed resource allocation scheme based on the max–min criterion to guarantee energy efficiency fairness by comparing with the total energy efficiency maximization scheme. Unless otherwise stated, the parameters listed in Table 1 are used in this section. In addition, the distance between the gateway and the PB is assumed to be 52 m; the distance between the three EUs and the PB is 1.8 m, 1.6 m, and 1.4 m, respectively. Based on the simulated channel gains and and the energy harvesting circuits, dedicated power stations produced by the powercast company, the antenna gains of the PB and gateways assume as 5 dBi, and the antenna gain of each EU is 2 dBi in this paper. Assume a carrier frequency of 915 MHz.

The convergence performance of the proposed resource allocation algorithm is shown in Figure 3. It can be seen that the proposed iterative algorithm converges to a certain constant after roughly three to four iterations, which validates the fast convergence of Algorithm 1. Secondly, we can see that has an important influence on the energy efficiency of the system link. As in Equation (4), the max–min EE should increase with the . And we can see that the larger the , the greater the max–min EE from the figure.

Figure 4 shows the achieved EE versus the performance gap . We compare the max–min link EE achieved by the proposed transmission strategy and the transmission strategies of the HTT mode and backscatter communication mode. As can be seen in Figure 4, the proposed transmission strategy achieves the max–min link EE no worse than these transmission strategies regardless of the value of . It is noted that the proposed hybrid transmission strategy can be degraded to these two transmission modes by adjusting parameters. It is consistent that the two transmission modes are special cases of the transmission strategy, and there is a tradeoff between the BackCom and HTT mode in the proposed resource allocation scheme as in theoretical analysis. Specifically, when is less than or equal to -18 dB, the max–min link EE achieved by the proposed transmission strategy is the same as the max–min link EE in the HTT mode. When is greater than -18 dB and less than -14 dB, the proposed transmission strategy makes EUs work in both backscatter and HTT modes in a full transmission time slot; when is greater than -14 dB, it can be seen that the max–min EE accomplished by the proposed transmission strategy is the same as the max–min link EE in the backscatter communication mode, which means that when is larger, the proposed transmission strategy degrades to the backscattered communication mode. In addition, compared to the throughput maximization resource allocation scheme, the EE by the proposed max–min fairness scheme presents better performance. The reason is that the throughput maximization scheme maximizing UE throughput in a BAWPCN does not take into account the factor of EE; the achieved EE is lower than that of the proposed max–min EE fairness scheme. The above analysis shows that the BAWPCNs can indeed combine the advantages of both backscatter and WPCNs, allowing them to adaptively adjust their parameters to meet the different communication goals in complex communication demands.

In Figure 5, we compare the total EE maximization resource allocation scheme with the proposed max–min link EE resource allocation scheme, under the case of large and small interuser channel differences. The major difference between the two resource allocation schemes is the objective function. According to Equation (9), the objective function is the minimum EE for the proposed max–min resource allocation scheme. For the total EE maximization resource allocation scheme, the objective function is the sum rates divided by the sum energy consumptions of EUs calculated as . It can be seen that the difference in EE between the best and worst EUs under the total EE maximization resource allocation scheme is significantly greater than that of the proposed max–min resource allocation scheme. When the channel difference between users is small, the average EE of users under the total EE maximization resource allocation scheme is slightly higher than that under the proposed max–min resource allocation scheme, but the difference between the EE of the best and worst users is significantly higher than the proposed resource allocation scheme. When the channel difference between users is large, it can be seen that the EE by the total EE maximization scheme results in the energy efficiency of the best user being 2.2 times that of the worst link, but the EE of the best link is about 1.4 times that of the worst link by the proposed resource allocation scheme, which effectively ensures fair access to resources between EUs. This is because the total EE maximization resource allocation scheme inclines the EUs with good channel status to maximize the EE, while the max–min scheme maximizes the EE of EUs with the worst link status to ensure fairness. As a result, the max–min scheme achieves much fairness at the sacrifice of less total EE.

5. Conclusion

This paper introduces a BAWPCN to address the energy shortage of sensor nodes due to frequent data interactions in digital twin workshops and proposes a resource allocation scheme to guarantee fairness in link EE. Considering the EUs’ minimum rate requirement and energy causality constraints, the max–min system link EE optimization problem was formulated as a mixed-integer nonconvex fractional programming problem with time-power two-dimension resource joint optimization. By introducing generalized fractional theory, relaxation variables, BCD, and auxiliary variables, an iterative algorithm is designed to solve the transformation problem to obtain the resource allocation scheme. Finally, the following three conclusions are verified by simulation: (1) the proposed iterative algorithm can quickly converge and (2) the proposed resource allocation scheme can effectively guarantee the fairness of link EE.

Appendix

In the optimization problem (14), the objective function, constraints C1-C4 and C7 are linear, so we only need to prove that constraints C5 and C8 are convex constraints. In C5 and C8, and are linear, so that a sufficient condition for the optimization problem (14) to be convex is that the Hessian matrix is seminegative definite. Construct the function , and its Hessian matrix can be expressed as follows:

The first-order determinant of the Hesse matrix shown in formula (A.1) is less than 0, and the second-order determinant is equal to 0, so the Hesse matrix is a seminegative definite matrix. Therefore, both constraints C5 and C8 are convex constraints. Based on the above analysis, Lemma 1 is proved to be right.

Data Availability

The simulation data used to support the findings of this study are included within the article. The code used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.