Transmission scheduling for wireless mesh networks with temporal reuse

Link-assigned transmission schedules with timeslot reuse by multiple links in both the space and time domains are investigated in this study for stationary multihop wireless mesh networks with both rate and power adaptivity. Specifically, cross-layer optimised schedules with proportionally fair end-to-end flow rates and network coding capability are constructed for networks operating under the physical interference model with single-path minimum hop routing. Extending transmission rights in a link-assigned schedule allows for network coding and temporal reuse, which increases timeslot usage efficiency when a scheduled link experiences packet depletion. The schedules that suffer from packet depletion are characterised, and a generic temporal reuse-aware achievable rate region is derived. Extensive computational experiments show improved schedule capacity, quality of service, power efficiency and benefit from network coding accrued with schedules optimised in the proposed temporal reuse-aware convex rate region.


Introduction
Wireless mesh networking maintains a strong position in future communications, especially since it can extend access from a base station or gateway. This prominent implementation scenario can also provide flexible last mile access for stationary wireless mesh clients, instead of more expensive wired options. The focus of this paper is thus on the optimisation of transmission scheduling for stationary wireless mesh networks that use link-assigned scheduling with extended transmission rights. Although properly researched, the transmission scheduling aspect of wireless network design should be re-evaluated in the light of recent advances in network coding (NC) and cross-layer network optimisation. Specifically, modifications are required to allow linkassigned schedules to perform opportunistic network coding and to improve timeslot usage efficiency when the schedule is predefined.
Client meshing in a wireless mesh network creates data transport flows between all node pairs, which causes a traffic pattern where opposing flows traverse the network links. This presents a high likelihood for a network coding opportunity, where a relay node can combine packets intended for exchange between two neighbouring nodes into one transmission. Wireless mesh networks are replete with network coding opportunities which have the potential of significantly improving network throughput, quality of service and power efficiency. Opportunistic network coding is simple to deploy in an existing network when there is no opportunistic listening involved. Link-assigned scheduling is designed in this paper with support for network coding in mind, without incurring a lesser spatial reuse penalty or increased difficulty of cross-layer optimisation.
The two predominant transmission schedule assignment methods are node-and link-assignment, and while nodeassigned schedules [1,2] allow for NC and have better usage efficiency, they lack the higher spatial reuse and finer granularity for optimisation achieved by link-assigned schedules [3][4][5]. Preference is given to link-assignment since it is more amenable to cross-layer optimisation, and it can achieve greater capacity owing to better spatial reuse. A novel hybrid node-link assignment method was presented in [6] (see also [7]), which keeps the high spatial reuse of link-assignment but integrates the NC capability and a constituent of the usage efficiency of node-assignment. This hybrid strategy was termed link-assignment with extended transmission rights (LET), where a node can attempt transmission on an alternative link if its primary scheduled link experiences packet depletion.
Packet depletion has been defined in [8] as the condition where a scheduled link does not have packets to transmit, resulting in a wasted timeslot with link-assigned scheduling. Conventional link-assignment only authorises one outgoing link per active node, and if that link is not utilised during its scheduled timeslot, then the opportunity and potential capacity are wasted. This lack of available packets for transmission on a link is a function of the network load, traffic pattern and transmission scheduling, and it is shown that this packet depletion occurs in most wireless mesh networks with client meshing where packets are added randomly to data flows. The hybrid assignment strategy provides alternative links for potential transmission, which can then make more efficient use of the available capacity in the case of primary link packet depletion. This beneficial aspect of some schedule assignment strategies in improving usage efficiency with multiple link options per node is analysed in this paper in terms of temporal reuse. The investigation of transmission scheduling with regard to temporal reuse is instructive in assisting with optimised schedule design, and it presents a novel viewpoint for the analysis of schedule assignment strategies.
The hybrid LET-schedules were cross-layer optimised in [8], according to the framework outlined in [9] (see also [10]). Specifically, cross-layer optimised schedules with proportionally fair end-to-end flow rates and network coding capability are constructed for networks operating under the physical interference model [11] with singlepath minimum hop routing. The adopted model does not capture the multi-carrier air interface, such as OFDM (orthogonal frequency-division multiplexing), but rather a single-channel physical interference model to place the focus on a capacity analysis for different transmission scheduling algorithms. The hybrid node-link assignment is redefined in terms of temporal reuse, which is proposed in this study as the time domain analogue of spatial reuse. The schedules that suffer from packet depletion and have the potential to benefit from temporal reuse are characterised, and an improved estimate of the achievable rate region is derived that incorporates the increased capacity afforded by temporal reuse.
The work in [8] is extended in this paper to include optimisation of networks with both variable transmission power and transmission at multiples of the base rate. When nodes can vary transmission power then the network spatial reuse potential and power efficiency can be increased, and when variable rates are allowed then network throughput and the benefit from network coding can be improved. Incorporating these changes into the temporal reuse-aware cross-layer optimisation is one of the unique contributions in this paper, where the capacity calculations, optimisation and complexity issues relating to variable power and rate control are dealt with.
Apart from the variable power and rate control, temporal reuse analysis and characterisation of networks that suffer from packet depletion, this paper further extends on the previous work of the authors in [8] by providing improved formulations for calculating the temporal reuse-aware rate region. The computationally feasible rate region is derived in a manner that enhances the understanding of the underlying factors that increase capacity with extended transmission rights. The capacity calculations and the optimisation program are also expanded to incorporate variable rate and power assignments, and the simulations include analyses of schedules with different rate and power objectives. In addition, a numerical analysis is done to determine the power efficiency of the various schedules considered in this paper. The main contribution is the calculation of the capacity region for wireless mesh networks that use predefined link-assigned transmission schedules with extended transmission rights and its use in a cross-layer optimisation framework similar to the one in [9].
Extensive computational experiments attest to the efficacy of the contributed generic scheduler in estimating the higher usage efficiency yielded by temporal reuse. Schedule capacity is measured as the constant number of packets per timeslot added to a random data flow such that a bounded mean end-to-end packet delay is conserved, and a greater capacity is achieved with the proposed temporal reuse-aware scheduler when compared to the non-aware scheduler for both schedules having temporal reuse. The capacity increases with NC are greater for the proposed scheduler, and in addition, flow rate fairness is improved for all combinations of rate/ power adaptivity and NC usage. Power efficiency, as related to the mean transmission power per packet, is also improved when the rate region used takes into account the benefits of temporal reuse.
The wireless networks and traffic models studied are defined in Section 2 and the achievable rate region and cross-layer optimisation framework for link-assigned schedules are reconsidered in Section 3. The role of temporal reuse in improving schedule performance is motivated in Section 4, where a new temporal reuse-aware rate region is also formulated. A generic temporal reuseaware schedule optimiser can be practically implemented with column generation according to the algorithm in Section 5, and the benefits of the new scheduler are determined in the computational experiments described in Section 6. fixed positions, have isotropic antennas with mutually exclusive reception and transmission capabilities, and have infinitely large buffers. A directed wireless link l (single-user Gaussian channel) is formed from a transmitting node tr(l) to a receiving node re(l) for a sufficiently high signal-to-interference and noise ratio (SINR) The transmit power of node n is P n (limited to 0 ≤ P n ≤ P max ), the thermal noise at receiver m is s m , and the effective power gain between transmitter n and receiver m is G nm . The gain is calculated according to a deterministic fading model G nm = K nm d −ρ nm , where K nm is a normalisation constant, d nm is the Euclidean distance between n and m, and r is a constant path loss exponent. The link SINR g l (P) is a function of the global power assignment (possible multiple transmitters owing to spatial reuse), which is denoted by P = (P n : n ∈ N ).
The target SINR for transmission at a discrete r-multiple of the base rate is γ (r) tgt , where r = 1 refers to the base transmission rate and r = 0 a zero transmission rate. In practice, the link rate is limited to c l = c tgt , where W is the system bandwidth in the Shannon capacity formulation. Links are activated during specific discrete timeslots, where all timeslot durations are equal and chosen such that one data packet can be forwarded in one timeslot. Perfect synchronisation between all nodes is assumed and every node operates according to a global transmission schedule.
The link gains are characterised in this paper only by the large-scale parameters, although realistically smallscale fading does exist even on the line-of-sight links between static nodes. CSI (channel state information) feedback is needed by the central scheduling agent in order to take the small-scale fading into account. However, since the focus is on centralised transmission scheduling, the main requirement is only that the network topology remains intact. By choosing a large SINR threshold (e.g. 10 dB) for the links in the topology, the small-scale fading analysis becomes less important. This allows for the study to focus on the capacity analysis of a fixed topology and not on the dynamics of the possible changes in topology due to small-scale fading.
Data traffic is modelled as F = N(N -1) transport flows f ∈ F between every possible source-sink node pair in the network, where data is generated at the source and delivered to the sink at flow rates s = (s f : f ∈ F ). Only single-path routing is considered in this study, and the link-flow incidence matrix R ℝ L×F couples the links with the flows so that 3 Link-assigned scheduling The purpose of a transmission schedule is to control link transmission powers and rates, regulate medium access and provide adequate capacity for all network traffic. Transmission schedules are predetermined for a limited number of timeslots, with the frame being repeated in time. The capacity provided during a specific timeslot is defined by the associated link rate vector, which gives the transmission rate for each active link in the spatial reuse group, or clique. The cross-layer optimisation framework in [9] is partly reviewed in this section to make the paper self-contained.

Link capacities
Each link can be assigned a total of r l, max + 1 different rates according to the link rate allocation policy tgt : r l ∈ [0, r l,max ] , where r l = 0 means an inactive link and r l = r l, max the maximum possible multiple of the base rate for the link l. A link rate vector is one possible network capacity assignment for a single timeslot, so that every link in the network is assigned a specific rate in the link rate vector. All K l∈L r l,max + 1 possible link rate vectors {v k , k [1, K]} are given by the set A schedule allocates a certain number of timeslots to each link rate vector, so that the long-term achievable link rate region can be defined as This is similar to the capacity region for wireless ad hoc networks defined in [12]. The convex hull conv() of the set V includes all achievable average link capacities owing to time-sharing of all the link rate vectors. This convex combination of link rate vectors allocates a time fraction of a k to every link vector v k to form the polyhedral rate region.

Schedule optimisation
Link capacities are represented by the transmission schedule, which can be optimised to maximise a global utility function. The level of satisfaction with a flow rate s f is given by a utility function u(s f ), which is concave and strictly increasing. Maximising the weighted network throughput f ∈F w f s f leads to grossly unfair flow rates [13], so a proportionally fair flow utility is decided on where u -∞ when s f 0 + . The flow rates s are proportionally fair if and only if u(s f ) = log(s f ) (see [14]), and this leads to a higher total network throughput than with maximisation of a uniform flow rate.
Network throughput depends on the available link capacities, which are consumed with relaying potentially multiple flows. The aggregate traffic load on the links is given by the vector Rs, which is due to the point-topoint data flows and the multihop nature of the mesh network. The achievable link rates are given by the capacity vector c ∈ C, so that the flow rates are limited by the primary constraint Rs ≼ c, where ≼ denotes a component-wise inequality. The purpose of the optimisation problem is to find the time-sharing coefficients (let K = [1, K])(let K = [1, K]) for a capacity vector c that maximises the global utility according to The concave global utility proposed in [9] is subject to convex constraints so it can be solved as a convex optimisation problem that produces a cross-layer optimised solution. The resultant schedule consists of link rate vectors {v k : k ∈ K} and the associated time-sharing coefficients {α k : k ∈ K}, which define transmission power and rate control in the physical layer and medium access control in the data link layer, respectively. The optimal solution involves at most L + 1 link rate vectors according to Carathéodory's theorem [15] (see [16]), so there are at most L + 1 non-zero a k values in the optimal schedule. The exact number of utilised link rate vectors depends on the specific optimal solution for a particular network. The greater the number of links, the more numerous the possible spatial reuse groups, which increases the number of used link rate vectors. For a proper TDMA (time-division multiple access) schedule, it can be seen that L link rate vectors would be needed for a feasible solution. When spatial reuse comes into play, less link rate vectors can describe the same capacity of TDMA, so the limit of the number of involved link rate vectors is not exceeded. The delivery route of each data flow is established in the network layer as represented by R, and the flow rates s are regulated in the transport layer.

Scheduling with temporal reuse
Temporal reuse is important in transmission scheduling, as it is a means to share capacity and improve schedule usage efficiency. This term is explained further in this section, and the networks that can have improved efficiency with temporal reuse are also characterised. The capacity region for link-assigned scheduling with extended transmission rights is formulated, and a computationally feasible suboptimum is proposed. This accurate temporal reuse-aware capacity region is then used in the subsequent sections to obtain the full schedule optimisation program.

Temporal reuse
Spatial reuse is the reuse of a timeslot by simultaneously activated links, with sufficient SINR, in the spatial domain. Likewise, temporal reuse is defined here as the reuse of a timeslot by a group of links in the time domain. A temporal reuse group is the set of outgoing links from a particular node, which may transmit in a certain timeslot. Since a transmission schedule is composed of a T-length frame that is repeated in time, every timeslot in the frame is repeated at frequency 1/T. This means that the links in a temporal reuse group will share the capacity of a timeslot over all the repetitions of the timeslot in which the temporal reuse group appears. The reuse of a timeslot in the time domain thus refers to a timeslot, at a particular location in the scheduling frame, that is reused by a group of links as the specific timeslot reoccurs with the repetition of the frame in time.
The purpose of temporal reuse is to increase schedule capacity usage efficiency in the case of link-assigned schedules, but it is also important in node-assignment to allow all links to transmit even though only nodes are scheduled. Node-assigned schedules [1] have to employ a form of temporal reuse to provide capacity for every outgoing link associated with a scheduled node. Having multiple link options per node results in higher capacity usage efficiency, at the expense of reduced spatial reuse. Alternatively, link-assigned schedules with extended transmission rights can achieve some measure of the same improved efficiency with better spatial reuse, which results in a higher schedule capacity.
In the case of node-assigned scheduling, each node forms its own temporal reuse group and secondary scheduling must be executed at each node to service each of the outgoing links of a node. However, temporal reuse groups are prioritised in link-assigned scheduling with extended transmission rights, which makes the schedule completely defined since the secondary scheduling is implied. In addition, by extending a link assignment to a prioritised group of links, the schedule allows for the implementation of opportunistic network coding, which further improves performance, as in the case of nodeassigned scheduling.

Schedule utilisation efficiency
A link must have a packet to transmit every timeslot that it is scheduled; otherwise, a usage inefficiency occurs due to packet depletion. The networks that can experience packet depletion and the associate usage inefficiencies are characterised in this subsection. Predefined schedule frames with a limited number of timeslots provide capacity π l for each link l, where π l packets per timeslot can be added to a data flow randomly chosen according to a uniform distribution, such that the schedule will be able to handle the resultant load on l.
Let p l be the probability that a packet added to a random flow will have to be routed on a link l, and let T l ≈ c −1 l be the average number of timeslots between the start of successive scheduled timeslots for link l with continual repetitions of the scheduling frame. The link scheduling period T l is approximately the inverse of the link rate c l , if the schedule is composed in a cyclical manner so that multiple schedulings of a link rate vector are as far removed in the scheduling frame as possible. Then, the average number of packets that have to be sent in a scheduled timeslot of link l is given by h l (π) in (6) for a global flow load π, which also represents the efficiency of the schedule in terms of link l, where at most one packet can be transmitted per timeslot.
This formulation sums the products of the probability that a certain number of i packets will be added to l in between successive schedulings of l with the number of packets i, to get the average number of added packets for l in each of its scheduled timeslots. This is a monotonically increasing function that relates a global load to a locally resultant load at a specific link l. The full efficiency load for link l is then the solution to h l (π l ) = 1, which equates to the applied global load that will fully utilise the specific link capacity provided by the schedule. The globally applied flow load is upper bounded by the link capacity that is saturated first; thus, the maximum packet load per timeslot is π = min l∈L π l . Thus, for every link with π m > π , the average utilisation efficiency is η m ( π) < 1, which means that the maximum load allowed is too small to make full use of the available capacity for link m.
Full schedule utilisation can thus only be achieved when η −1 l (1) = η −1 m (1), ∀l, m ∈ L; otherwise, for any inequality π l < π m , l ≠ m, it must mean that h l (π l ) <1. The global load may not exceed the minimum π l , because then link l will experience a load greater than its capacity, so only when the full efficiency loads for all the links are equal can every link be saturated without exceeding the capacity of any one link. For networks with spatial reuse and multihop flows, it becomes more difficult to calculate schedules that give equal full efficiency load π l for each link l ∈ L, and even if it is achieved, there is still no trivial guarantee that the schedule will achieve optimal schedule capacity and fairness. The optimisation problem is thus further constrained by , ∀l, m ∈ L, which mainly constrains the link scheduling periods T l , rather than flow-link probabilities p l which are determined by the network topology and routing.
By extending transmission rights to alternative links in the case of packet depletion, the unused timeslots of the unsaturated links can be used to increase the capacity for most of the links in the network to π m + ε m . Consequently, the maximum schedule capacity will be increased to π = min l∈L (π l + ε l ), and the schedule will have a higher possible usage efficiency due to the applied temporal reuse.

Link capacities with temporal reuse
Given a scheduled link rate vector v k ∈ V that experiences packet depletion for one of the associated links L(v k ) = l ∈ L : v k (l) > 0 , the scheduler can resort to transmission on an alternative link through temporal reuse. In such a timeslot where an alternate link is used, another link rate vector v m ≠ v k describes the capacity used for the temporal reuse instance. All possible link rate vectors that can describe temporal reuse for the primary vector v k are defined by the extended set An extended set of link rate vectors can be associated with any one specific link rate vector, such that spatial reuse is maintained but with increased temporal reuse. The node transmission powers for a link rate vector v k are given by P(v k ), and the power assignment should not be exceeded by any alternate link rate vector.

Temporal reuse-aware rate region
The average link capacities used for a scheduled extended set depend on its ordering policy, the network topology, radio propagation model and traffic load, as these factors influence the prevalence of packet depletion. For all repetitions of a timeslot where an extended set of v k is scheduled, a certain fraction km of repetitions will be represented by every link rate vector v m ∈ V(v k ) in the extended set. These fractions can only be determined by simulating the network, although the values will change when a schedule is used that is optimised using those values in the temporal reuse-aware rate region The aware rate region C e gives the achievable link capacities with the extended temporal reuse sets when taking the inefficiencies due to packet depletion into account. Optimising a schedule with more accurate representations of the link rate vectors leads to improved network performance.

Rate region estimation
A schedule optimisation algorithm should be independent of a specific network topology and radio propagation model, so several assumptions and simplifications are applied to create a generic optimiser.
(i) Each alternative link in the temporal reuse set must operate at the maximum possible transmission rate, within the power constraints.
(ii) The usage probabilities of links in a temporal reuse set are determined only on the basis of the order in which the links are queried for packets in the case of packet depletion, regardless of network topology or radio propagation model. The validity of this assumption is tested in the usage efficiency analysis in the simulation section, where it is shown that it is reasonable.
Since some links in a temporal reuse set could have sufficient SINR to transmit at different rates, only the highest rate (primary criterion) is chosen to maximise the capacity of each extended set. For schedules that use link rate vectors with variable power control, a minimum power criterion (secondary with a small ε >0) could also be used to differentiate between multiple link rate vectors with a maximum rate for a specific link. The temporal reuse set for a clique node n ∈ T (v k ) = n ∈ N : L(v k ) ∩ O(n) = ∅ , where O(n) are the outgoing links of node n, which is thus calculated as To calculate the rate region C e , a very large set of (V) = v k ∈V V(v k ) values particular to V is needed. A network-agnostic rate region can be rendered by focusing on the usage probabilities of alternative links in L n (v k ), instead of alternative link rate vectors as in the case of (V). With temporal reuse, when a primary scheduled link experiences packet depletion, the alternative links are queried in a certain order for available packets, according to the ordering policy, until a usable link is found.
Each link l in a temporal reuse set L n (v k ) is assigned an ordering value n (v k , l) ∈ [0, |L n (v k )| − 1], where a smaller value means a higher priority, to establish an ordering policy, so a link l ∈ L n (v k ) with a higher priority Θ n (v k , l) <Θ n (v k , m) than m ∈ L n (v k ) will have a higher usage probability U n (v k , l) > U n (v k , m). The number of usage probabilities is reduced with the assumption that the utilisation probabilities are approximately equal when the usage priorities are the same, regardless of differences in networks, the extended sets, the particular node temporal reuse sets or links in question. This simplification is made on the basis that it is the order in which links in the temporal reuse set are checked that determines the link usage probabilities. Let U be the simplified priority-based link usage probabilities, where the usage probability for a link with priority θ is U θ .
The primary outgoing link for a node, i.e. the active link that is originally implicated in a v k , has priority θ = 0 and lower priority reserve links have higher θ values. The simplified usage probabilities of outgoing links in the temporal reuse set for a specific node n active in v k are given by Although this form does not define a specific ordering policy, one based on relative traffic flow magnitudes is presented in the following subsection. Let I L be the identity matrix in ℝ L×L , where I L (l) is the jth column in I L , so that j is the index of link l in L. Temporal reuse link sets are mutually exclusive; thus, for any v k ∈ V, the relationship ∩ n∈T (v k ) L n (v k ) = ∅ holds, so the usage probability vector for a v k with temporal reuse is given by The full link rate template for v k with temporal reuse given by C(v k ) represents the non-achievable capacity before taking the actual usage limits specified by U(v k ) into account.
The Hadamard (Schur) product of the usage probability estimates U(v k ) and the rate template C(v k ) give the resultant achievable capacity for v k with full temporal reuse, so the simplification of C e gives the estimated achievable rate region for networks with spatial and temporal reuse. Apart from the time-sharing coefficients {α k : k ∈ K} and the link rate vector power assignments {P(v k ) : v k ∈ V}, each schedule has to include the prioritized temporal reuse link rate sets { (v k ), U(v k ), C(v k ) : v k ∈ V} associated with every active transmitter node n in the spatial reuse node set (clique) T (v k ) for every v k ∈ V. The schedule is composed by cycling between the maximum of L+1 involved link rate vectors until there are ⌈a k δ (L + 1)⌉ occurrences of each vector v k . A value of δ = 3 has been found to give good results for various network sizes.
U is independent of a particular network topology or radio propagation model and can be determined once with the averaged values of simulations of many different networks and used to optimise any network afterwards, as done in the numerical analysis of Section 6. The number of utilisation probability values needed to calculate C e is This estimate C e is thus greatly simplified while maintaining the effect of temporal reuse on the usage efficiency and achievable capacity of link-assigned schedules.

Ordering policy
The discriminator that is used to assign priorities to the alternative links in a temporal reuse set is taken as the number of data flows that traverse a link. The more flows that have to be routed by a link, the greater the need for capacity; thus, preference is given to links that relay more flows. This traffic flow-based ordering policy assigns higher priority to the higher traffic link as formulated by

Column generation approach
The global utility maximisation problem in (5) determines the time-sharing coefficients for each link rate vector in V. The number of link rate vectors is exponential in the number of network links, so calculating and storing them in order to solve (5) becomes a hurdle as the network grows. The solution also tends to be more inaccurate with computer solvers when more link rate vectors are involved.
For these reasons, the mathematical programming technique of column generation [17] is employed, where modifications specific to temporal reuse are made to the algorithm discussed in [9], which is reconsidered in this section. Optimisation can then be started with only a small number of link rate vectors, and more vectors that contribute most to increasing the global utility are added through column generation. The process of optimising a schedule with column generation, for networks that employ both spatial and temporal reuse, is explained in this section.

Main algorithm
The utility maximisation problem guides the column generation by providing information pertaining to capacity demands, so it is referred to as the master problem. The full master problem in (5) contains all K link rate vectors, but it is restricted to a subset J ⊆ K of vectors in the restricted master problem (16) used for column generation.
When solving the restricted master problem, the capacity demand l l for each link l can be determined, which indicates the relative utility improvement for a unit increase in link capacity. For u(·) = log(·), each link demand also equals the inverse flow rate of the flow that only uses the corresponding link [18]. These capacity demands are given by the Lagrangian duality for (16), which is as follows.
The optimal Lagrangian multipliers l* are used to update J by maximising the upper bound on the optimum (weak duality, see [9]) according to The lower bound (19) involves only the columns generated thus far, as indicated by J , to find the best supply vector.
Finding a new link rate vector from the K-element rate region C e that maximises capacity for the given demand solves the scheduling subproblem sup c∈ C e λ * T c . An extreme point in C e will maximise capacity [9] rather than a combination of link rate vectors. This generated column does not require that all K elements have to be known, but the new vector can be obtained by solving a mixed integer linear programming formulation discussed in the proceeding subsection.
Optimality is reached when the new link rate vector cannot supply a greater capacity than an existing vector. The process of solving the restricted master problem and expanding J is continued until the difference between the optimal upper bound u u and lower bound u l becomes smaller than a certain threshold ε, as shown by (20).
The main column generation algorithm that is used to optimise a schedule, by forming J to maximise the global utility, is as follows: 1) Initialise the index set J with the link rate vectors that correspond to a simple TDMA link-assigned schedule operating at base rates only with no spatial reuse.
2) Solve the restricted master problem to optimality with a primal-dual interior-point method [19] to obtain the optimal Lagrangian multipliers or capacity demands l*.
3) Add the new link rate vector, that solves the scheduling subproblem, to J . 4) Repeat the process until u uu l < ε, or when the generated column is already present in J .

Scheduling subproblem
The optimal Lagrangian multipliers λ * l give the sensitivity of the lower optimum bound u l at the optimal point c*. By satisfying the links with the greatest demand for capacity as indicated by λ * l , the lower bound can be improved most. By finding an extreme vertice in c ∈ C e to maximize λ * T l c, a new scheduling element is added to increase the network utility. This is termed the scheduling subproblem, and for networks with temporal reuse, it is formulated as finding A mixed integer linear program is used to determine which spatial reuse clique with extended transmission rights solves the scheduling subproblem. Let the binary variable x (r) l indicate whether link l is operating at r times the base rate.
The binary variable y n specifies whether node n is transmitting or not, so the following definition is included: Whenever a link is transmitting at a non-zero rate, the associated node also has to be activated, so this relationship is formulated as In the program, the instantaneous maximum transmission power of a node n is given by P n which is limited to a maximum as shown in (25): To ensure that all transmitting links have a sufficiently large SINR, the following linear constraint is imposed: By choosing a sufficiently large M (r) l , the above constraint will be met when the link is not transmitting, so the following value will suffice: For scheduling with fixed power, the node transmit powers are fixed at P n = P max , ∀n ∈ N , and for nonadaptive rate scheduling the rate is limited to r {0, 1}. Given the optimal Lagrange multipliers, the variable rate and power mixed integer linear program, which considers all possible temporal reuse priority policies, is thus formulated as By ensuring that the links with greater traffic flow loads are assigned higher rates in a temporal reuse link set, the complexity can be reduced. The complexity of the scheduling subproblem still remains NP-hard, however. The following constraint can be added to the linear program to achieve the complexity reduction:

Numerical analysis
Extensive simulation was conducted to evaluate the performance of C e -schedules optimised in the temporal reuseaware rate region as compared to C-schedules optimised in the default non-aware rate region. The simulation model is expounded on, and results are included that show the improvements that are gained with the addition of temporal reuse. The usage probability vector U is estimated through simulation, and the process is explained and values given for networks with differing rate/power adaptivity. Schedule performance is gauged according to schedule capacity, fairness and power efficiency, and in this section, the benefits of the temporal reuse-aware rate region are made clear. Finally, the relationship between spatial and temporal reuse for the various schedules is investigated. Multihop wireless mesh networks of 10 to 30 nodes are randomly generated in two dimensions where the maximum length of each dimension for 10, 15, 20, 25 and 30 node networks are set at 140, 160, 180, 200 and 220 m respectively. The various square dimensions have been experimentally chosen to produce fully connected networks with a good average spatial reuse potential that increases as the number of nodes increases. Singlepath minimum hop routing is used according to Dijkstra's algorithm, and only networks with a path between every node pair are considered. For each data point in the results, the averaged outcomes of 120 random networks are taken and each network size is analysed separately.

Simulation model
Schedules are determined centrally and shared with all nodes, so that with perfect synchronisation each node operates according to the schedule. Each data flow is communicated with packets where each packet can be relayed over one hop in one timeslot when transmitting at the base rate. The schedule capacity is determined as the maximum number of packets π that can be added per timeslot to a flow randomly chosen according to a uniform distribution, such that the mean end-to-end packet delays are bounded. The numerical analysis determines the mean end-to-end packet delays as well as the mean flow rates and their standard deviation and the mean transmission power per packet at the maximum schedule capacity.
The analyses with NC use opportunistic two-way network coding [20] where a relay node XORs two packets intended for exchange between two neighbouring nodes. Neighbouring nodes do not need to communicate which packets each received, neither is opportunistic listening required with two-way coding, and the throughput gain with n-way coding is maximised when n = 2. First-infirst-out (FIFO) transmit buffers are employed and with NC the buffer is checked in a FIFO manner to obtain the first two packets that can be coded and transmitted over the scheduled temporal reuse set.

Improvements with temporal reuse
The effect of temporal reuse on schedule capacity is determined for schedules optimised in the unaware rate region C with transmission rights extended afterwards. This C-optimisation does not consider the capacity changes that extended transmission rights cause and is thus temporal reuse-unaware. Every link rate vector v k that is activated in a C-schedule has transmission rights extended to L n (v k ) − 1 = q − 1 alternative links for q >0 to see what improvements in the schedule capacity results.
As the number of alternative links is increased, the schedule capacity π also increases, as shown in Figure 1a. So as transmission rights are extended, the added capacity min l∈L ε l (q + n) > min l∈L ε l (q) expands for n >0, where ε l is a function of q, with greater improvements seen for smaller network sizes owing to smaller spatial reuse and thus better temporal reuse. The data population sizes for the large-extreme end of the x-axis are much smaller, since not all networks can achieve such large temporal reuse. This data population size inequality and the trade-off between the schedule capacity and endto-end packet delays cause the reverse trend in the packet delay values at a large temporal reuse.
The schedule optimised in the convex rate region C is named CVXOPT-LET in Figure 1, where LET refers to link-assignment with extended transmission rights (full temporal reuse). The mean end-to-end packet delays for the schedule operating at capacity are given in Figure 1b and are reduced for most of the rights extension. Thus, both the schedule capacity and the average end-to-end packet delays can be improved with temporal reuse.

Estimation of usage efficiency
Temporal reuse-unaware C-optimised schedules with full rights extension (full temporal reuse) are analysed to estimate the usage probabilities U of links in temporal reuse sets based on their assigned priorities θ. The average usage probabilities of 10-30 node networks are given in Table 1 for networks operating with schedules that have no power-rate adaptivity, power adaptivity only, rate adaptivity only and both power and rate adaptivity.
The mean (for different power-rate adaptivities) of the standard deviations due to networks having different sizes is given in the column for SD A. The accuracy of the given usage probabilities was determined by comparing them to that of temporal reuse-aware C e -schedules and is given in the Acc. column. The validity of the assumption that different network realisations of the same size have similar usage probabilities is measured with the mean (for different power-rate adaptivities) of the mean (for different network sizes) of the standard deviations due to different network topologies as given in the column for SD B. The deviations are small compared to the means for the higher priority links, so the assumption is taken as sound.
Each column displays the effect of temporal reuse, which shares capacity among the temporal reuse set to increase timeslot usage efficiency. The primary links (corresponding to θ = 0) are not fully used, which indicates significant packet depletion occurring owing to unequal link capacities in proportionally fair cross-layer optimised schedules. These results thus show the important benefits that can be derived from temporal reuse.

Convergence of the master problem
The number of spatial reuse groups is limited for sizeconstrained fixed-topology networks, and the resulting link rate vectors are all determined. The accurate temporal reuse-aware link rate vectors are also determinable, and the number of vectors remains the same. The column generation must therefore end, since the number of link rate vectors is limited. Due to Carathéodory's theorem, there is a high probability that the master problem will converge before all the link rate vectors have been generated. The moment no new link rate vector can increase the capacity supply, the column generation converges.
Slater's condition is met for fully connected networks, since there is a strictly feasible solution for the master problem. One such solution is given by a simple linkassigned TDMA schedule that supplies a non-zero  positive capacityc l for each link l in the network. The network is fully connected so a non-zero positive flow ratesp can be achieved by every flow p in the network, and if a flow rate of min l c l / p r lp is assigned to each flow, then the convex constraint Rs c met. Since Slater's condition is met strong duality holds, and there is thus no duality gap. This means that the column generation must converge. The convergence of the column generation for the temporal reuse-aware C e -schedule with extended transmission rights is depicted in Figure 2. The reduction in the difference between the lower and upper bounds of the master problem is shown for different network sizes. The mean convergence rate was measured for 120 random network realisations, for each network size. Since larger networks have more spatial reuse groups and data flows, the convergence takes more iterations and the starting difference between the bounds is greater. The absence of the duality gap is confirmed, since all of the solutions converge with zero difference between the lower and upper bounds.

Schedule performance comparison
The practical superiority of schedule optimisation in the temporal reuse-aware rate region is proven through measurements of schedule capacity, flow rate fairness and transmit power efficiency.

Schedule capacity
Schedule capacity is taken as the primary schedule performance measure, and it is the maximum number of packets per timeslot that can be added to a random flow, so that the schedule maintains a limited mean end-to-end packet delay. The C-schedule (CVX) capacity is increased with extended transmission rights (CVX-LET) and when network coding is used (CVX-LET-NC), the capacity is increased even further, as shown in Figure 3a. The temporal reuse-aware C e -schedule with extended transmission rights (CVXAware) has greater capacity than with CVX-LET, and the capacity improvement with NC (CVXAware-NC) is also notably larger than with CVX-LET-NC. The temporal reuse-aware rate region gives greater capacity to link rate vectors with good temporal reuse, which provides more NC opportunities for the C e -schedule.
Reference graphs include those of TDMA, STDMA, traffic-adaptive STDMA with temporal reuse and NC (STDMA-Adaptive-LET-NC), and traffic-adaptive nodeassigned schedules with NC (Node-Assigned-NC). The STDMA schedule uses a packing heuristic [21], which packs as many links as have not yet transmitted into a timeslot, and the traffic-adaptive STDMA schedule is generated according to [22]. The node-assigned schedule lets a node choose a link with probability proportional to the number of flows that traverse the link, to achieve traffic adaptivity.
The mean end-to-end packet delays for the different schedules operating at their respective maximum capacities are given in Figure 3b. Node-assigned schedules have the lowest delay, although they do not have the capacity of the cross-layer optimised link-assigned schedules. The C-schedule (CVX) achieves lower delays when temporal reuse (CVX-LET) is enabled, even though the delays for CVX are taken at a smaller load. The timeslot duration is a function of the packet length, which affects the quality of service if the packet delays grow large. Although apart from the proportional fairness, the end-to-end packet delays are allowed to grow as large as possible while remaining bounded. This is warranted to obtain the  maximum schedule capacity, which is a more important metric in this study than a secondary quality of service concern.
The effect of rate and power adaptivity on the capacity of the cross-layer optimised schedules is shown in Figure 4. Schedules with power adaptivity (CVX-Power-LET, CVX-RatePower-LET, CVXAware-Power, CVXAware-RatePower) perform better because of more diverse spatial reuse considerations during optimisation. When NC is enabled for the temporal reuse-aware schedules (CVXAware), greater capacity increases are seen than with the C-schedules, as shown in Figure 4b.

Flow rate fairness
The second important schedule performance measure is flow rate fairness, which indicates the quality of service level that can be expected with a schedule. For the schedules operating at maximum capacity, the associated mean flow rates are compared in Figure 5a, where nodeassigned schedules have the highest rate and CVX the lowest. Node-assigned schedules achieve the highest mean flow rate at the expense of rate fairness, and the best fairness is achieved by CVX as measured by the flow rate standard deviations in Figure 5b. The sum of the logarithms of the achieved flow rates during network simulation can also be used to assess proportional fairness, but the standard deviation of the achieved flow rates is calculated instead to accompany the mean achieved flow rate. The temporal reuse-aware schedule (CVXAware) achieves a better fairness than the unaware (CVX-LET) schedule, and fairness is not compromised with the use of NC for the CVXAware schedule.
For the rate-power adaptive schedules, significantly lower flow rate deviation occurs with the temporal reuseaware schedules than with the unaware schedules, since the proportionally fair global utility is more accurately determined when taking the benefits of temporal reuse into account. The improved fairness with the C e -schedules is especially apparent with the power-adaptive schedules as shown in Figure 6b.
A network routes most of its data flows through the centre of the network, so by giving a centralised link more capacity, there can be more flows enjoying a higher transmission rate. This may increase the mean flow rate at the expense of reduced fairness, due to the reduced capacity experienced by flows routed across links on the edges of the network. A proportional fairness utility can balance the schedule capacity and fairness objectives explicitly, as opposed to the disregard of fairness by schedules such as node-assignment.

Transmit power efficiency
Power efficiency is used as the third schedule performance indicator and is related to the mean transmit    Figure 7. By performing network coding (CVX-LET-NC), one transmission can relay two packets, so only half the energy is used in that case, which improves the power efficiency. Significantly less power is used with the power-adaptive schedules, and when network coding is used, the temporal reuse-aware schedules achieve greater added efficiency than the nonaware schedules.

Spatial and temporal reuse
There is a relationship between spatial and temporal reuse, since greater spatial reuse implies greater interference and consequently smaller temporal reuse. These results are evident in Figure 8, where the power-adaptive schedules have the greatest spatial reuse (Figure 8a) but the smallest temporal reuse (Figure 8b). Nodeassigned schedules in turn achieve the highest temporal reuse, which reduces the spatial reuse possibilities, so that these also achieve the smallest spatial reuse. The power-adaptive C e -schedules give greater preference to temporal reuse when compared to the C-schedules, which have higher spatial reuse. Thus, with temporal reuse-aware optimisation a more beneficial balance is found between spatial and temporal reuse, which results in favourable network performance.

Conclusion
The research question that has been addressed in this study is as follows: What generic transmission scheduling algorithm can produce predefined link-assigned schedules that maximise schedule capacity with proportionally fair end-to-end flow rates for stationary multihop mesh networks with single-path minimum hop routing, rate and power adaptivity, both spatial and temporal reuse and network coding?
The following contributions were made in answering this question: 1) The hybrid node-link assignment strategy of [6] is re-evaluated in terms of a new concept that is defined here, namely temporal reuse, which is recognised as the time-domain analogue of spatial reuse. 2) A more specific characterisation of the schedules that could induce packet depletion is given than in the previous study of [8]. It is then shown that optimal schedules can potentially benefit from higher usage efficiency afforded by temporal reuse.
3) A new achievable rate region is formulated that provides better estimates of the increased link capacities due to temporal reuse. The rate region was derived to be generic, meaning that only information on the network topology and data flows are needed to calculate its capacity without having to simulate the network in question.

4) A cross-layer optimisation scheduling algorithm is
given that uses the newly proposed rate region and solves the scheduling subproblem with temporal reuse in mind.
Firstly, it was experimentally shown that by optimising schedules with temporal reuse in the temporal reuseaware rate region, the schedule capacity is increased when compared to optimisation in the non-aware rate region of [9]. In addition, flow rate fairness is also improved with the proposed scheduler and the benefit derived from network coding is also greater. For schedules with rate and power adaptivity, it was shown that the proposed scheduler achieves better power efficiency.
The relationship between spatial and temporal reuse was analysed, and the temporal reuse-aware scheduler managed a more effective balance for the different rate and power adaptivities. The trade-off between schedule capacity and fairness was affirmed, and it was shown that the proposed class of optimised link-assigned schedules outperforms node-assigned schedules in that respect.
The research in this paper thus includes a high-performance scheduler for modern wireless mesh networks that use recent enhancements, such as network coding and temporal reuse. Possible extensions of this work might include an analysis for networks with multi-path routing, cooperative coding and more complex n-way (n >2) network coding.