# Network Coded Multicast over Multibeam Satellite Systems

Authors: R. Alegre-Godoy; M. A. Vazquez-Castro
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

We propose a multicast scheme for multibeam satellite
systems exploiting both the multiuser and spatial diversity inherent
in this type of systems while taking into account realistic physical
distributions of User Terminals (UTs) over the coverage. Our proposed
scheme makes use of the well-known Adaptive Coding and Modulation
(ACM) feature in Digital Video Broadcasting over Satellite, 2nd Generation
(DVB-S2) and Extension (DVB-S2X) standards but also incorporates a
set of innovative features. First, multilink reception, that is, receivers
that can syntonize different frequencies and/or polarizations, together
with Network Coding (NC) is used to enable decoding of signals from
adjacent beams (spatial diversity). Second, efficient and fair allocation
of resources is achieved through Proportionally Fair Multicast (PFM)
scheduling. Our results, obtained over realistic nonuniform UTs distributions,
show average system multicast throughput gains up to 88% with regard
to state-of-the-art multicast schemes. Furthermore, a complete cross-layer
architecture is proposed, fully compliant with the standard providing
Quality of Service (QoS) guarantees.


---

## Body

## 1. Introduction

Recent studies are showing an increasing demand for the efficient distribution of personalized contents in Internet based networks [1]. This has led to the deployment of satellite platforms delivering high throughputs (HTS systems) such as Ka-SAT [2] or constellations of communication satellites such as the O3b system [3]. Recent works even consider collocating two satellites covering the same region in order to cope with the user needs [4, 5]. Beyond the aforementioned satellite physical aspects, it is also possible to satisfy the user needs by improving the logical mechanisms delivering multimedia contents. Multicasting is one of the cornerstones for the effective dissemination and distribution of personalized multimedia contents in broadband networks and the focus of this paper. Applications such as audio/video streaming, online gaming, file distribution, and file downloading are based on multicast-like transmissions.In wireless networks, including multibeam satellite networks, the main challenge when multicasting is how to address the heterogeneous channels conditions of the User Terminals (UTs), that is, the presence of multiuser diversity. In clear-sky conditions, the difference in Signal to Interference plus Noise Ratio (SINR) between a UT located at the center of the beam and a UT located at the edge is typically 2-3 dB. When Adaptive Coding and Modulation (ACM) is adopted at the physical layer as in Digital Video Broadcasting over Satellite, 2nd Generation (DVB-S2) [6] and DVB-S2 Extension (DVB-S2X) [7] standards, this difference in SINR is translated into a spectral efficiency difference of 11%–25% [6, 8]. The differences in spectral efficiency can be much higher if we consider a beam partially affected by a rain event.Traditional colouring schemes in multibeam systems allow many opportunities for exploiting spatial diversity. A UT can potentially access a number of orthogonal transmissions from the adjacent beams [8]. In that case, UTs in the border of the beam would be more advantaged than those in the center of the beam. Since current UTs are syntonized at a single frequency, orthogonal transmissions are not exploited. A multilink receiver, that is, a receiver able to syntonize different frequencies and polarizations, could access and decode the signal meant for adjacent beams. As mentioned in [8], the design of this type of receivers is perfectly possible with the current technology. UTs distribution also plays a fundamental role. In real life, users are not uniformly distributed but concentrated in specific areas such as cities. The concentration of terminals in specific areas of the coverage also affects the performance of cellular based systems [9, 10].In this paper, we take into account these three aspects, multiuser and spatial diversity and UTs distribution to design a multicasting scheme for the efficient delivery of broadband contents.(a) Related Works on Satellite Multicast. Multicast in multibeam satellite systems has been little investigated from the scheduling and resources allocation point of view. The authors in [11] propose to choose a fixed Reed Solomon code and a fixed rate out of a set of possible rates in order to accomplish a certain degree of reliability. Following a similar approach to [11], works in [12, 13] propose to choose a modulation and codification (MODCOD) which ensures reception to a subset of the UTs in the multicast group. The rest could only decode with a certain probability. In [14], a traditional approach is adopted and information is multicasted according to the channel conditions of the worst UT in the multicast group. With respect to these works, our approach takes advantage of multiuser diversity in order to select in each time-slot the optimal MODCOD rather than assuming a fixed scheme. To the best of our knowledge, the authors in [15] provide the most similar approach to our work since the scheme they propose is based on ACM. In particular, authors propose a Network Utility Maximization (NUM) to trade delay and rate also accounting for Quality of Service (QoS) and multiuser diversity. With respect to this work, we introduce the novelty of multilink reception UTs together with Network Coding (NC) which enables decoding orthogonal transmissions.Furthermore, our paper breaks the traditional approach of assuming uniform UTs distributions and provides results and analysis for nonuniform distributions which are close to reality.(b) Related Works on NC for Multibeam Satellite Systems. In the past few years, a number of works have studied the implementation of NC in multibeam satellite systems. In [16], an overview of possible satellite scenarios where NC can be applied is provided. For the particular case of multibeam satellite systems, NC is proposed as a mechanism to reduce retransmissions. Works in [17, 18] take advantage of the orthogonal transmissions available using multilink reception. However, the focus is on unicast transmissions and NC is used to provide enhanced reliability and flexibility rather than increasing the throughput. Our previous works in [19, 20] assessed the feasibility of using NC for multicasting in multibeam satellite systems. As a result, it was identified that the multilink reception approach together with NC coding could bring important benefits subject to the location of the UTs. These papers describe a preliminary concept and lack of a(i)
method to decide when to use the multilink reception with NC feature,(ii)
scheduling policy selecting the optimal MODCOD for the multicast service.(c) Contributions of the Paper. This paper proposes a full multicasting scheme, that is, scheduling policy, packet scheduling architecture, and algorithm to decide if the multilink reception with NC feature must be used or not. Our work presents the following novel results with respect to use of NC technology:(i)
A technique for the joint use of multilink reception, NC, and Proportionally Fair Multicast (PFM).And it also presents the following novel results in the field of satellite multicasting: (i)
Introduce and adapt the PFM scheduling concept in [21]. More specifically, we provide MODCOD selection and use of the multilink reception with NC feature when suitable.(ii)
Its associated cross-layer packet scheduling architecture with respect to the Internet Engineering Task Force (IETF) differentiated services at IP level model (Diffserv).(iii)
A scheme providing multicast throughput gains employing the same resources as a traditional multicast scheme demonstrated via analysis and simulations over theoretical and realistic nonuniform UTs distributions.The rest of the paper is organized as follows: Section 2 introduces the multibeam satellite system model. Section 3 describes the proposed multicast scheme. The packet scheduling architecture is introduced in Section 4. Finally, Section 5 provides numerical evaluation of the system performance and Section 6 draws conclusions on the work done.


## 2. System Model

### 2.1. Multibeam Satellite System Model

We assume a multibeam and multigateway satellite system withNb beams, P polarizations, and frequency reuse factor fr. The number of colours of the system is Nc=Pfr. Forward link transmissions are based on DVB-S2/DVB-S2X with ACM. Each gateway (GW) is associated with a subset of the overall number of beams (or cluster). GWs receive channel state information (CSI) messages from the UTs through a feedback channel.Let a GW of the system serve the subset of beamsC={b1,b2,…,bC}. Each beam bj⊂C has assigned a number of UTs uj requesting the same multicast service. Let us derive the SINR for a UT i in bj, 1≤i≤uj. First, let the number of cochannel beams of the overall system be K=Nb/Nc. Now, we define H∈CK×K as the forward link channel matrix which can be decomposed as H=BG. Matrix B∈CK×K accounts for the atmospheric, propagation, space, and ground system effects and is defined as(1)B=diag⁡β1,β2,…,βK,where βj=OBOHPALsatLpropGUT, with OBOHPA being the output back-off of the satellite high power amplifier, Lsat the satellite repeater losses, Lprop the propagation losses, and GUT the UT antenna gain. Matrix G∈CK×K accounts for the square root of the satellite antennas gains towards the concrete position of the UT and is defined as(2)G=g11⋯g1K⋮⋱⋮gK1⋯gKK,where gvw stands for the square root of antenna gain for antenna v towards the location of the UT at beam w. Therefore, each element hvw∈H accounts for all the gains and losses from satellite antenna v towards UT location at beam w. The received signal yi,j at UT i∈bj can be expressed as(3)yi,jx,y=Psat⁡hjjx,ysj+Psat⁡∑l≠jhljx,ysl+ni,where Psat is the satellite transmitted power, ni is the Gaussian noise (zero mean complex circular noise of variance Ni), and sj and sl are the transmitted and interfering symbols, respectively. Assuming constant transmitting power, the SINR can be extracted directly from (3) and is given by(4)Γi,jx,y=Psat⁡hjjx,y2Psat⁡∑l≠jhlj2+Ni.Under the ACM specification of DVB-S2 and DVB-S2X, SINR values are mapped to spectral efficiencies (or equivalently MODCODs) as follows:(5)ηi,j=fDVB-S2Γi,j,ηi,j=fDVB-S2XΓi,j,where ηi,j is the spectral efficiency for UT i in bj and fDVB-S2, fDVB-S2X are mapping functions that relate SINRs and spectral efficiencies for DVB-S2 and DVB-S2X standards, respectively.


### 2.2. Multilink Reception System Model

The work in [8] introduces and models multilink multibeam systems. Such systems assume the use of multilink receivers, that is, receivers which can syntonize different frequencies or polarizations to simultaneously decode orthogonal transmissions from adjacent beams.The main concept is as follows. In a multibeam system withNc colours, UTs can potentially decode up to Nc transmissions, 1 transmission from the own beam and Nc-1 transmissions from adjacent beams in orthogonal frequencies and/or polarizations. This effect is produced because the antenna gain of each spot-beam is so high, that even UTs outside of the beam observe values of SINR that lie within the range of available MODCODs and can decode the signal. To do so,(i)
UTs must have multilink reception capabilities, for example, a terminal with a single antenna, one Low Noise Block downconverter (LNB), and multiple reception chains to detect and decode different polarizations and bands;(ii)
UTs must observe a value of SINR higher than or equal to the one required to decode the MODCOD transmitted in the orthogonal beam. Conversely, the GW can lower the MODCOD transmitted in a beam to let a number of UTs outside of the beam decode the signal.In our multicasting scheme, we assume that UTs can decode their own transmission and one out of the Nc-1 orthogonal transmissions, more specifically, the transmission with strongest SINR or equivalently the transmission from the closest adjacent beam. This enables an extra path to reach each UT and the opportunity to exploit spatial diversity. Let subscript a denote the adjacent beam whose signal intends to decode UT i in bj and let Γi,ja(x,y) be the SINR observed from this adjacent beam. Γi,ja(x,y) is obtained particularizing (4) with beam a, that is, substituting hjj with haj and computing the interference power from the cochannel beams of a. The spectral efficiency achievable from such adjacent beam is denoted as ηi,ja and obtained particularizing (5) with Γi,ja(x,y).Figure1 shows the MODCODs achievable for different locations (in U/V coordinates [22]) within a beam of a 70-beam system. More specifically, the top plot shows, per each location, the adjacent beam providing better SINR (the location and the adjacent beam with the strongest SINR are plotted in the same color and dashed lines separate the different areas). The mid and bottom plots show the achievable MODCODs from the own and determined adjacent beam transmissions, respectively. It can be observed that any point of the beam can decode an orthogonal transmission and that locations close to the edge of the beam and in the beam overlapping areas can decode it with a high order MODCOD. Locations in the center of the beam could only decode signals employing low spectral efficiency MODCODs.Figure 1
Example of achievable MODCODs for locations within a beam of interest. Results extracted from a 70-beam, 4-coloured system under DVB-S2.


### 2.3. Network Coding Model

In our multicast scheme, UTs can simultaneously decode their own transmission and one orthogonal transmission. The objective of NC is to avoid receiving duplicated versions of the same packet and ease the scheduling at the GW side. Therefore, whenever our proposed scheme sends packets over two paths, those packets will be coded employing NC. We perform NC at the IP layer for the following reasons:(i)
Our objective is to increase the average multicast throughput in a single multicast service. IP allows identifying the content of the packets and classifying them into different multicast services such that NC is performed in a intraservice manner. Working at lower layers would not allow this classification.(ii)
IP headers can be easily modified to include the encoding coefficients necessary for decoding [23].In particular, we adopt Random Linear Network Coding (RLNC) due to its simplicity and capacity-achieving performance, although any other network code would be, in principle, applicable. The scheme generates a random linear coded packet from Np native packets as follows. Let L denote the length in bits of a native IP packet. Then, the payload of each packet pt is split into blocks of ϵ bits. Let dtk denote the kth block of tth packet, with 1≤k≤L/ϵ. Next, we choose Np random coefficients ct from Fq, with q=2ϵ. The kth coded block of a coded packet is given by(6)Ck=∑t=1Npctdtk,1≤k≤Lϵ.The encoding coefficients c1,…,cNp are added in the header of the packet unless pseudorandom codes are used and then only the seed needs to be sent. In the receiver side, at least Np coded packets must be received to retrieve the original packets. The parameter ϵ is the size of the finite field and should be big enough to ensure that the probability of generating two linearly dependent (l.d.) coded packets is negligible. Sections 3 and 4 introduce the specific parameters and architectural details necessary to perform RLNC.


### 2.4. User Terminals Distribution

We employ the method in [24] to generate the nonuniform distributions of UTs. Specifically, we break the area of interest, that is, the cluster of beams C, in a number of bins Nz. Each bin z is assigned a different probability Pz of containing a UT according to a function ϕz. In this case, ϕz is a function of the coordinates xz, yz of the bin and given by a symmetrical 2D truncated Gaussian function [9]:(7)ϕzxz,yz=αexp⁡-αxz-x02+yz-y02π1-exp⁡-αR02,where x0,y0 is the UTs distribution center, R0 is the radius, and α is the location index. To express ϕz as a probability, it must be normalized by the total sum of ϕz over all the Nz bins considered:(8)Pzxz,yz=ϕzxz,yz∑z=1Nzϕzxz,yz.Then, we can model UTs distribution in a cluster C by the quadruple (α,R0,x0,y0) and analyze average multicast throughput results varying these parameters. Negative values of α concentrate UTs towards the edge of the cluster whilst positive values of α concentrate UTs around x0,y0. α=0 provides a uniform distribution.


## 2.1. Multibeam Satellite System Model

We assume a multibeam and multigateway satellite system withNb beams, P polarizations, and frequency reuse factor fr. The number of colours of the system is Nc=Pfr. Forward link transmissions are based on DVB-S2/DVB-S2X with ACM. Each gateway (GW) is associated with a subset of the overall number of beams (or cluster). GWs receive channel state information (CSI) messages from the UTs through a feedback channel.Let a GW of the system serve the subset of beamsC={b1,b2,…,bC}. Each beam bj⊂C has assigned a number of UTs uj requesting the same multicast service. Let us derive the SINR for a UT i in bj, 1≤i≤uj. First, let the number of cochannel beams of the overall system be K=Nb/Nc. Now, we define H∈CK×K as the forward link channel matrix which can be decomposed as H=BG. Matrix B∈CK×K accounts for the atmospheric, propagation, space, and ground system effects and is defined as(1)B=diag⁡β1,β2,…,βK,where βj=OBOHPALsatLpropGUT, with OBOHPA being the output back-off of the satellite high power amplifier, Lsat the satellite repeater losses, Lprop the propagation losses, and GUT the UT antenna gain. Matrix G∈CK×K accounts for the square root of the satellite antennas gains towards the concrete position of the UT and is defined as(2)G=g11⋯g1K⋮⋱⋮gK1⋯gKK,where gvw stands for the square root of antenna gain for antenna v towards the location of the UT at beam w. Therefore, each element hvw∈H accounts for all the gains and losses from satellite antenna v towards UT location at beam w. The received signal yi,j at UT i∈bj can be expressed as(3)yi,jx,y=Psat⁡hjjx,ysj+Psat⁡∑l≠jhljx,ysl+ni,where Psat is the satellite transmitted power, ni is the Gaussian noise (zero mean complex circular noise of variance Ni), and sj and sl are the transmitted and interfering symbols, respectively. Assuming constant transmitting power, the SINR can be extracted directly from (3) and is given by(4)Γi,jx,y=Psat⁡hjjx,y2Psat⁡∑l≠jhlj2+Ni.Under the ACM specification of DVB-S2 and DVB-S2X, SINR values are mapped to spectral efficiencies (or equivalently MODCODs) as follows:(5)ηi,j=fDVB-S2Γi,j,ηi,j=fDVB-S2XΓi,j,where ηi,j is the spectral efficiency for UT i in bj and fDVB-S2, fDVB-S2X are mapping functions that relate SINRs and spectral efficiencies for DVB-S2 and DVB-S2X standards, respectively.


## 2.2. Multilink Reception System Model

The work in [8] introduces and models multilink multibeam systems. Such systems assume the use of multilink receivers, that is, receivers which can syntonize different frequencies or polarizations to simultaneously decode orthogonal transmissions from adjacent beams.The main concept is as follows. In a multibeam system withNc colours, UTs can potentially decode up to Nc transmissions, 1 transmission from the own beam and Nc-1 transmissions from adjacent beams in orthogonal frequencies and/or polarizations. This effect is produced because the antenna gain of each spot-beam is so high, that even UTs outside of the beam observe values of SINR that lie within the range of available MODCODs and can decode the signal. To do so,(i)
UTs must have multilink reception capabilities, for example, a terminal with a single antenna, one Low Noise Block downconverter (LNB), and multiple reception chains to detect and decode different polarizations and bands;(ii)
UTs must observe a value of SINR higher than or equal to the one required to decode the MODCOD transmitted in the orthogonal beam. Conversely, the GW can lower the MODCOD transmitted in a beam to let a number of UTs outside of the beam decode the signal.In our multicasting scheme, we assume that UTs can decode their own transmission and one out of the Nc-1 orthogonal transmissions, more specifically, the transmission with strongest SINR or equivalently the transmission from the closest adjacent beam. This enables an extra path to reach each UT and the opportunity to exploit spatial diversity. Let subscript a denote the adjacent beam whose signal intends to decode UT i in bj and let Γi,ja(x,y) be the SINR observed from this adjacent beam. Γi,ja(x,y) is obtained particularizing (4) with beam a, that is, substituting hjj with haj and computing the interference power from the cochannel beams of a. The spectral efficiency achievable from such adjacent beam is denoted as ηi,ja and obtained particularizing (5) with Γi,ja(x,y).Figure1 shows the MODCODs achievable for different locations (in U/V coordinates [22]) within a beam of a 70-beam system. More specifically, the top plot shows, per each location, the adjacent beam providing better SINR (the location and the adjacent beam with the strongest SINR are plotted in the same color and dashed lines separate the different areas). The mid and bottom plots show the achievable MODCODs from the own and determined adjacent beam transmissions, respectively. It can be observed that any point of the beam can decode an orthogonal transmission and that locations close to the edge of the beam and in the beam overlapping areas can decode it with a high order MODCOD. Locations in the center of the beam could only decode signals employing low spectral efficiency MODCODs.Figure 1
Example of achievable MODCODs for locations within a beam of interest. Results extracted from a 70-beam, 4-coloured system under DVB-S2.


## 2.3. Network Coding Model

In our multicast scheme, UTs can simultaneously decode their own transmission and one orthogonal transmission. The objective of NC is to avoid receiving duplicated versions of the same packet and ease the scheduling at the GW side. Therefore, whenever our proposed scheme sends packets over two paths, those packets will be coded employing NC. We perform NC at the IP layer for the following reasons:(i)
Our objective is to increase the average multicast throughput in a single multicast service. IP allows identifying the content of the packets and classifying them into different multicast services such that NC is performed in a intraservice manner. Working at lower layers would not allow this classification.(ii)
IP headers can be easily modified to include the encoding coefficients necessary for decoding [23].In particular, we adopt Random Linear Network Coding (RLNC) due to its simplicity and capacity-achieving performance, although any other network code would be, in principle, applicable. The scheme generates a random linear coded packet from Np native packets as follows. Let L denote the length in bits of a native IP packet. Then, the payload of each packet pt is split into blocks of ϵ bits. Let dtk denote the kth block of tth packet, with 1≤k≤L/ϵ. Next, we choose Np random coefficients ct from Fq, with q=2ϵ. The kth coded block of a coded packet is given by(6)Ck=∑t=1Npctdtk,1≤k≤Lϵ.The encoding coefficients c1,…,cNp are added in the header of the packet unless pseudorandom codes are used and then only the seed needs to be sent. In the receiver side, at least Np coded packets must be received to retrieve the original packets. The parameter ϵ is the size of the finite field and should be big enough to ensure that the probability of generating two linearly dependent (l.d.) coded packets is negligible. Sections 3 and 4 introduce the specific parameters and architectural details necessary to perform RLNC.


## 2.4. User Terminals Distribution

We employ the method in [24] to generate the nonuniform distributions of UTs. Specifically, we break the area of interest, that is, the cluster of beams C, in a number of bins Nz. Each bin z is assigned a different probability Pz of containing a UT according to a function ϕz. In this case, ϕz is a function of the coordinates xz, yz of the bin and given by a symmetrical 2D truncated Gaussian function [9]:(7)ϕzxz,yz=αexp⁡-αxz-x02+yz-y02π1-exp⁡-αR02,where x0,y0 is the UTs distribution center, R0 is the radius, and α is the location index. To express ϕz as a probability, it must be normalized by the total sum of ϕz over all the Nz bins considered:(8)Pzxz,yz=ϕzxz,yz∑z=1Nzϕzxz,yz.Then, we can model UTs distribution in a cluster C by the quadruple (α,R0,x0,y0) and analyze average multicast throughput results varying these parameters. Negative values of α concentrate UTs towards the edge of the cluster whilst positive values of α concentrate UTs around x0,y0. α=0 provides a uniform distribution.


## 3. Multicast Proportional Fairness for Multibeam Satellite Systems

Figure2 shows a functional diagram of the proposed multicast scheme based on a cross-layer architecture. Novel blocks, at GW and UT side, are shadowed in green. We define multicast as the transmission of a data flow (or multicast service) from the GW to a number of UTs in the cluster C controlled by the GW. The UTs requesting the service are known as the multicast group. The scheme we propose applies to each multicast service individually (intraservice), and hence explanations of the multicast policies assume a single multicast service. Scheduling is performed on a slot-by-slot basis; that is, in the current slot we choose how packets are to be transmitted in the next slot. The scheme has been designed such that if follows the coherence time of the channel (i.e., time windows in the scheduling are not required because the satellite channel for fixed UTs shows only medium/long term variations even under rain conditions).Figure 2
Functional diagram of the proposed multicast scheme based on a cross-layer architecture.On the GW side, we introduce 3 novel blocks which carry out the Link Layer (LL) scheduling and substitute classical scheduling policies. Let us assume a system with ACM at the physical layer and letM be the number of MODCODs of the system. We denote by ηm the spectral efficiency of the mth MODCOD, with 1≤m≤M. The bitrate provided by ηm is given by(9)rm=Rsηm,where Rs denotes the system symbol rate. The objective of LL scheduling, formed by the 3 green blocks, is twofold: (1) Choose the most appropriate transmission mode in each time-slot; that is, decide whether the multilink reception capabilities of the UTs have to be exploited or not. (2) Choose the appropriate value of spectral efficiency (i.e., MODCOD) for the next time-slot according to the transmission mode selected and the CSI reports received from the UTs in the current slot. This is equivalent to say the multicast policy is scheduling a certain number of UTs (those able to decode the selected MODCOD). In particular, (i)
the block Proportionally Fair Multicast (PFM) computes per each beam of the clusterC (and per multicast service) which is the most suitable MODCOD according to the proportionally fair rule, and hence it exploits multiuser diversity;(ii)
the block PFM with NC (PFM-NC) computes for the entire clusterC a common MODCOD for all the beams of the cluster delivering that service; data is sent network coded using RLNC and UTs use its multilink reception capabilities to attempt decoding the signal from its own beam and its adjacent beam with higher SINR;(iii)
the block combined PFM (C-PFM) selects the most appropriate transmission mode, PFM (nonuse of multilink reception capabilities), or PFM-NC (use of multilink reception capabilities) in each time-slot by computing the mode obtaining the highest multicast rate;(iv)
the inputs required for the PFM, PFM-NC, and C-PFM blocks are the MODCODs achievable by the UTs in base to their CSI feedback. The LL Resource Manager (LLRM) provides this information to these 3 blocks. PFM, PFM-NC, and C-PFM blocks are described in detail in Sections 3.1, 3.2, and 3.3, respectively. The rest of the GW follows a traditional approach such that our proposed scheme can be implemented via software modifications with minimal intrusion. Incoming multicast traffic is classified per beam, per QoS following the IETF Diffserv model, and per multicast service before being forwarded to the LL schedulers. The specific packet scheduling architecture of the GW is described in detail in Section 4.The membership management of the multicast groups is supported by the Internet Group Management Protocol (IGMP) [25]. The LLRM collects information of the multicast groups available at the IGMP block and the CSI of the UTs. UTs not only send the CSI with respect to its own beam but also the CSI with respect to all its adjacent beams. The collaboration between the IGMP block and the LLRM provides the following information: (1) which beams are involved in the transmission of a certain multicast service; (2) the number of UTs per beam and its SINR with respect to the own and adjacent beams; (3) combining (1) and (2), identifying the orthogonal transmission with the strongest SINR.On the UT side, the RF chain must be able to syntonize two different bands (frequencies and/or polarizations) in order to obtain the packets from the own beam and from the strongest adjacent beam. Physical layer and link layer packets are deencapsulated in both bands up to the IP layer. At the IP layer, the UT checks if packets are network coded by looking at the IP headers:(i)
If not, packets from the adjacent beam are discarded and the packets from the own beam are deencapsulated as usual and forwarded to the upper layers.(ii)
If yes, decoding of the network coded packets is performed before deencapsulating and forwarding the packets to the upper layers.Although this may appear as a complex receiver, the multilink reception does only need a single antenna and a single LNB as explained in Section 2.2. Even if multiple reception chains to detect and decode different bands are necessary, such kind of terminals already exists and is defined for instance in Digital Video Broadcasting Satellite to Handheld (DVB-SH) for mode B devices [26]. DVB-S2X also provides tools to implement transmission of a single stream through several physical satellite transponders and how to carry out its reception through several DVB-S2X decoders in the receiver (the so-called channel bonding for multitune receivers [7]). With regard to the decoding of network coded packets, a number of works have successfully assessed the feasibility of the technology in broadband networks; see, for instance, [27–29].


### 3.1. Proportionally Fair Multicast (PFM)

In this subsection, we introduce the PFM block in Figure2. In a PF allocation, UTs with better channel state are scheduled more often than UTs with worse channel states. To this aim, the GW performs the scheduling in two steps. (1) Compute the optimal per-beam MODCOD for the next time-slot. (2) Update the scheduled rates of each UT according to the MODCOD selected in Step (1) for the next iteration.Step 1.
The MODCOD chosen for the next time-slot at beambj out of the M available MODCODs is given by(10)mjpft+1=arg⁡maxm⁡∏i∣ri,jt≥Rsηm1+RsηmRi,jt∀j∈1,C,where Ri,j(t) accumulates the scheduled rates of UT i in bj and ri,j(t) is the achievable bitrate by UT i in bj given by ri,j(t)=Rsηi,j(t).Step 2.
LetSjpf(t+1) denote the set of UTs in bj which will be scheduled at t+1 and given by(11)Sjpft+1=i∣ri,jt≥rmjpft+1,where rmjpf(t+1)=Rsηmjpf(t+1) is the bitrate provided by spectral efficiency ηmjpf(t+1) at beam bj. Subsequently, Ri,j is updated by(12)Ri,jt+1=Ri,jt+rmjpft+1,i∈Sjpft+1Ri,jt,elsewhere.The PFM scheme takes advantage of multiuser diversity because (10) selects the most appropriate MODCOD per time-slot according to the channel conditions across the UTs.


### 3.2.  Proportionally Fair Multicast with Network Coding (PFM-NC)

We now introduce the PF block employing NC (PFM-NC) in Figure2 to take advantage of the extra path provided by orthogonal transmissions. To this aim, the GW considers the cluster as one entity with a number of UTs U=∑j=iC‍uj and performs the scheduling in 3 steps. (1) Select for each UT the orthogonal transmission with the strongest SINR from which it will receive packets in addition to those sent through its own beam. (2) Compute the optimal MODCOD to use in the cluster for the next time-slot taking into account orthogonal transmissions. (3) Update the scheduled rates of each UT according to the MODCOD selected in Step (2) for the next iteration.Since UTs receive packets from two paths employing the same MODCOD, the min-cut/max-flow theorem states that up to 2 different packets can be received simultaneously. Due to the layout of the multibeam network, the only way to ensure that all scheduled UTs will receive both packets is to send packets network coded. For instance, Figure3(a) shows an example where one of the UTs would only receive one of the packets p1 and p2 independently of the packets sent by each queue. Instead, if every queue generated a single coded packet, all UTs would receive p1 and p2. Moreover, NC simplifies the scheduling because the system just transmits combined packets instead of looking at the topology of the network and decides which packets should be selected in each queue for transmission; for example, Figure 3(b) shows that all UTs can obtain p1 and p2 but previous knowledge of the topology is necessary to decide which packet forwards each queue. Hence, under PFM-NC, the per-beam queues of a specific multicast service will generate a single coded packet from the same Np=2 native packets as explained in Section 2.3. Then, all the UTs able to decode the chosen MODCOD will always get two different coded packets. Those UTs unable to decode orthogonal transmissions will get only one coded packet and will not be able to convey the original packets; that is, they are considered as nonscheduled UTs.Figure 3
Two examples of inefficient multicasting from a GW transmitting uncoded packets to three UTs in three beams. Continuous and dashed arrows represent transmissions from the own and best adjacent beam, respectively.
(a)
(b)Step 1.
Select for each UT an adjacent beam, that is, the orthogonal transmission with strongest SINR, in order to decode simultaneously from two different paths. Such adjacent beam is chosen as follows. LetAj={a1,…,aAj} be the set of adjacent beams to bj and let Γi,jae be the SINR of UT i in bj with respect to the adjacent beam ae. UT i in bj selects the adjacent beam a such that(13)e∗=arg⁡maxe⁡Γi,jae∣ae∈Aj,ae∈C,ue≠0.The related spectral efficiency of Γi,jae∗, ηi,jae∗, is obtained from (5). The conditions in the argmax in (13) restrict the selection of ae to those adjacent beams that belong to C and have at least one UT subscribed to the multicast service. Note that this differentiates our scheme from load balancing techniques since we do not use sparing resources from an adjacent beam but we take advantage of the resources already in use in the adjacent beam. If none of the adjacent beams accomplish the conditions in (13) or the SINR from the chosen adjacent beam lies without the MODCOD range, then ηi,jae∗=0.Step 2.
The optimal MODCOD in the PFM-NC scheme is given by(14)mnct+1=arg⁡maxm⁡∏i∣ri,jae∗t≥Rsηm1+2RsηmRi,jt,where ri,jae∗(t)=Rsηi,jae∗. Thus, given a value of Rsηm, the product in (14) takes into account only those UTs that can decode at such rate from the orthogonal transmission. Implicitly, this means that the UT can support that rate from the own transmission and therefore factor 2 in the numerator of (14).Step 3.
LetSnc(t+1) denote the set of UTs in C scheduled at t+1 and given by(15)Snct+1=i∣ri,jae∗t≥rmnct+1,where rmnc(t+1)=Rsηmnc(t+1) is the bitrate provided by spectral efficiency ηmnc in the cluster C. Subsequently, Ri,j is updated by(16)Ri,jt+1=Ri,jt+2rmnct+1,i∈Snct+1Ri,jt,elsewhere.When a UT is scheduled, we take into account that it receives information from two paths by adding a factor of 2 in the second term of the scheduled rates.


### 3.3. Combined PFM/PFM-NC (C-PFM)

In the two previous subsections, we have introduced the two blocks of Figure2 based on the PF rule, PFM, and PFM-NC. In short, (1) PFM exploits multiuser diversity. Under clear-sky conditions, it tends to favour UTs in the center of the beam. (2) PFM-NC exploits multiuser and spatial diversity. Under clear-sky conditions, it tends to favour UTs in the edge of the beam. Following, we introduce the control block to select the most appropriate scheme for transmitting in each time-slot, that is, the selection block named C-PFM. Let fjpf(ηm) denote the product in (10) and let fnc(ηm) denote the product in (14). fjpf(ηm) and fnc(ηm) evaluate the average multicast rates at a given ηm for the PFM and PFM-NC schemes, respectively. Hence, at a given time-slot, we compute mjpf(t+1) and mnc(t+1) and evaluate the scheme providing better average multicast rate by comparing the outputs of fjpf(ηmjpf(t+1)) and fnc(ηmnc(t+1)). Since the latter takes into account all the UTs in the cluster, we perform the comparison with the product of all fjpf(ηmjpf(t+1)). Then, scheduled rates must be updated according to the scheme selected using (12) or (16). Algorithm 1 shows the C-PFM multicasting scheme we propose.Algorithm 1:C-PFM.
(1)  At time-slott obtain mjpf(t+1)∀j and mnc(t+1) based on (10) and (14) respectively.
(2) Computefjpf(ηmjpf(t+1))∀j and fnc(ηmnc(t+1)).
(3)if  ∏j=1|C|fjpf(ηmjpf(t+1))≥fnc(ηmnc(t+1))  then
(4)Transmit in PFM mode with MODCOD mjpf(t+1)∀j
(5)Update rates according to (12)
(6)else
(7)Transmit in PFM-NC mode at MODCOD mnc(t+1)
(8)Update rates according to (16)
(9)end if
(10) Next time slot


## 3.1. Proportionally Fair Multicast (PFM)

In this subsection, we introduce the PFM block in Figure2. In a PF allocation, UTs with better channel state are scheduled more often than UTs with worse channel states. To this aim, the GW performs the scheduling in two steps. (1) Compute the optimal per-beam MODCOD for the next time-slot. (2) Update the scheduled rates of each UT according to the MODCOD selected in Step (1) for the next iteration.Step 1.
The MODCOD chosen for the next time-slot at beambj out of the M available MODCODs is given by(10)mjpft+1=arg⁡maxm⁡∏i∣ri,jt≥Rsηm1+RsηmRi,jt∀j∈1,C,where Ri,j(t) accumulates the scheduled rates of UT i in bj and ri,j(t) is the achievable bitrate by UT i in bj given by ri,j(t)=Rsηi,j(t).Step 2.
LetSjpf(t+1) denote the set of UTs in bj which will be scheduled at t+1 and given by(11)Sjpft+1=i∣ri,jt≥rmjpft+1,where rmjpf(t+1)=Rsηmjpf(t+1) is the bitrate provided by spectral efficiency ηmjpf(t+1) at beam bj. Subsequently, Ri,j is updated by(12)Ri,jt+1=Ri,jt+rmjpft+1,i∈Sjpft+1Ri,jt,elsewhere.The PFM scheme takes advantage of multiuser diversity because (10) selects the most appropriate MODCOD per time-slot according to the channel conditions across the UTs.


## 3.2.  Proportionally Fair Multicast with Network Coding (PFM-NC)

We now introduce the PF block employing NC (PFM-NC) in Figure2 to take advantage of the extra path provided by orthogonal transmissions. To this aim, the GW considers the cluster as one entity with a number of UTs U=∑j=iC‍uj and performs the scheduling in 3 steps. (1) Select for each UT the orthogonal transmission with the strongest SINR from which it will receive packets in addition to those sent through its own beam. (2) Compute the optimal MODCOD to use in the cluster for the next time-slot taking into account orthogonal transmissions. (3) Update the scheduled rates of each UT according to the MODCOD selected in Step (2) for the next iteration.Since UTs receive packets from two paths employing the same MODCOD, the min-cut/max-flow theorem states that up to 2 different packets can be received simultaneously. Due to the layout of the multibeam network, the only way to ensure that all scheduled UTs will receive both packets is to send packets network coded. For instance, Figure3(a) shows an example where one of the UTs would only receive one of the packets p1 and p2 independently of the packets sent by each queue. Instead, if every queue generated a single coded packet, all UTs would receive p1 and p2. Moreover, NC simplifies the scheduling because the system just transmits combined packets instead of looking at the topology of the network and decides which packets should be selected in each queue for transmission; for example, Figure 3(b) shows that all UTs can obtain p1 and p2 but previous knowledge of the topology is necessary to decide which packet forwards each queue. Hence, under PFM-NC, the per-beam queues of a specific multicast service will generate a single coded packet from the same Np=2 native packets as explained in Section 2.3. Then, all the UTs able to decode the chosen MODCOD will always get two different coded packets. Those UTs unable to decode orthogonal transmissions will get only one coded packet and will not be able to convey the original packets; that is, they are considered as nonscheduled UTs.Figure 3
Two examples of inefficient multicasting from a GW transmitting uncoded packets to three UTs in three beams. Continuous and dashed arrows represent transmissions from the own and best adjacent beam, respectively.
(a)
(b)Step 1.
Select for each UT an adjacent beam, that is, the orthogonal transmission with strongest SINR, in order to decode simultaneously from two different paths. Such adjacent beam is chosen as follows. LetAj={a1,…,aAj} be the set of adjacent beams to bj and let Γi,jae be the SINR of UT i in bj with respect to the adjacent beam ae. UT i in bj selects the adjacent beam a such that(13)e∗=arg⁡maxe⁡Γi,jae∣ae∈Aj,ae∈C,ue≠0.The related spectral efficiency of Γi,jae∗, ηi,jae∗, is obtained from (5). The conditions in the argmax in (13) restrict the selection of ae to those adjacent beams that belong to C and have at least one UT subscribed to the multicast service. Note that this differentiates our scheme from load balancing techniques since we do not use sparing resources from an adjacent beam but we take advantage of the resources already in use in the adjacent beam. If none of the adjacent beams accomplish the conditions in (13) or the SINR from the chosen adjacent beam lies without the MODCOD range, then ηi,jae∗=0.Step 2.
The optimal MODCOD in the PFM-NC scheme is given by(14)mnct+1=arg⁡maxm⁡∏i∣ri,jae∗t≥Rsηm1+2RsηmRi,jt,where ri,jae∗(t)=Rsηi,jae∗. Thus, given a value of Rsηm, the product in (14) takes into account only those UTs that can decode at such rate from the orthogonal transmission. Implicitly, this means that the UT can support that rate from the own transmission and therefore factor 2 in the numerator of (14).Step 3.
LetSnc(t+1) denote the set of UTs in C scheduled at t+1 and given by(15)Snct+1=i∣ri,jae∗t≥rmnct+1,where rmnc(t+1)=Rsηmnc(t+1) is the bitrate provided by spectral efficiency ηmnc in the cluster C. Subsequently, Ri,j is updated by(16)Ri,jt+1=Ri,jt+2rmnct+1,i∈Snct+1Ri,jt,elsewhere.When a UT is scheduled, we take into account that it receives information from two paths by adding a factor of 2 in the second term of the scheduled rates.


## 3.3. Combined PFM/PFM-NC (C-PFM)

In the two previous subsections, we have introduced the two blocks of Figure2 based on the PF rule, PFM, and PFM-NC. In short, (1) PFM exploits multiuser diversity. Under clear-sky conditions, it tends to favour UTs in the center of the beam. (2) PFM-NC exploits multiuser and spatial diversity. Under clear-sky conditions, it tends to favour UTs in the edge of the beam. Following, we introduce the control block to select the most appropriate scheme for transmitting in each time-slot, that is, the selection block named C-PFM. Let fjpf(ηm) denote the product in (10) and let fnc(ηm) denote the product in (14). fjpf(ηm) and fnc(ηm) evaluate the average multicast rates at a given ηm for the PFM and PFM-NC schemes, respectively. Hence, at a given time-slot, we compute mjpf(t+1) and mnc(t+1) and evaluate the scheme providing better average multicast rate by comparing the outputs of fjpf(ηmjpf(t+1)) and fnc(ηmnc(t+1)). Since the latter takes into account all the UTs in the cluster, we perform the comparison with the product of all fjpf(ηmjpf(t+1)). Then, scheduled rates must be updated according to the scheme selected using (12) or (16). Algorithm 1 shows the C-PFM multicasting scheme we propose.Algorithm 1:C-PFM.
(1)  At time-slott obtain mjpf(t+1)∀j and mnc(t+1) based on (10) and (14) respectively.
(2) Computefjpf(ηmjpf(t+1))∀j and fnc(ηmnc(t+1)).
(3)if  ∏j=1|C|fjpf(ηmjpf(t+1))≥fnc(ηmnc(t+1))  then
(4)Transmit in PFM mode with MODCOD mjpf(t+1)∀j
(5)Update rates according to (12)
(6)else
(7)Transmit in PFM-NC mode at MODCOD mnc(t+1)
(8)Update rates according to (16)
(9)end if
(10) Next time slot


## 4. Packet Scheduling Architecture

In this section, we describe the packet scheduling architecture for the proposed multicast scheme. Although our work focuses on how to improve the average multicast throughput of a single multicast service, we show how our scheme integrates with an architecture offering more services. We build upon the architecture provided in [30] for unicast transmissions and optimize it for multicast transmissions. Our scheme requires modifications only at the LL level of the scheduler and can be implemented by modifying the software/firmware of the GW. As in [30], QoS between the different multicast services is achieved adopting the IETF Diffserv model at IP level which is a well-known and prevailing satellite QoS provision model [31, 32]. Specifically, 3 Classes of Service (CoS) are defined: Expedited Forward (EF) for premium class services, Assured Forward (AF) for less constrained traffic, and Best Effort (BE) for services with very low or without QoS. Figure 4 shows the proposed packet scheduling architecture for the EF CoS. The structure would be analogous for the AF and BE CoS. The cross-layer control functionalities introduced in Section 3 (see Figure 2) are not shown in the figure to ease the understanding.Figure 4
Proposed packet scheduling architecture according to IETF Diffserv model at IP level.Let us assume that the GW of clusterC transmits σEF, σAF, and σBE multicast services per EF, AF, and BE CoS, respectively. First of all, incoming multicast traffic is classified per beam, per QoS, and per multicast service resulting in a number of IP queues σEF, σAF, and σBE per CoS. After that, per each CoS, a coordinated Round Robin (RR) stage pulls packets from a specific service within each CoS. Without loss of generality, let us assume that packets from EF CoS service 1 are pulled. At this point, C-PFM computes, from the UTs achievable MODCODs, which is the optimal transmission mode according to Algorithm 1:(i)
If PFM is selected, the optimal MODCOD per each beam of the clusterC is computed according to (10). Then, per each beam area, a packet from service 1 is pulled, encapsulated using Generic Stream Encapsulation (GSE) or Multiprotocol Encapsulation (MPE), and forwarded to the appropriate MODCOD queue.(ii)
If PFM-NC is selected, for the entire clusterC, a single optimal MODCOD is computed according to (14). Now, per each beam area, two packets from service 1 are pulled and coded into a single packet using RLNC. This packet is encapsulated using GSE or MPE and forwarded to the appropriate MODCOD queue.GSE/MPE packets are buffered in the MODCOD queues. The Tuneable Fairness Weighted Round Robin (TF-WRR) scheduler in [30] is substituted by a simple First come first served (FCFS) packet selector. Whilst the TF-WRR in [30] allows unicast sessions within a CoS to be fairly treated, this concept is already provided for each multicast session by our combination of PFM, PFM-NC, and C-PFM blocks. In the last stage before the DVB-S2/DVB-S2X modulator, the Adaptive Weighted Round Robin (AWRR) scheduler provides the QoS guarantees among the different CoS.


## 5. Simulation Results

In order to show results for the proposed scheme, we focus on clusters of 6/7 beams under the control of a single GW within a system of the characteristics of Table1. Note that this is a feasible approach since HTS systems require several GW feeder links in order to provide UTs with the required throughputs for broadband and multimedia applications [33, 34]. The UTs distributions in the cluster are generated as detailed in Section 2.4 and for realistic case in Spain. The main figure of merit we analyze is the average system multicast throughput defined as the time average bitrate sum offered to the UTs. Since our scheme applies to multicast services individually and with the purpose of not masking the actual results of the proposed scheme, we focus on the performance of a single multicast performance. The performance of our proposed scheme is compared to the Worst Case Multicast (WCM) scheme [14] and compared to a scheme with only the PFM implementation (Section 3.1) such that the improvements of the multilink reception with NC can be easily observed. In the following subsections, we justify the use of DVB-S2X for our particular scheme, provide simulation results for clear-sky and rainy conditions, and give practical insights on how to provide a Service Level Agreement (SLA) for UTs that are not so frequently scheduled.Table 1
System characteristics.
Parameter
Value
N
b
70
N
c
4
Frequency
19.95 GHz
System symbol rate
180 Mbauds
PHY protocol
DVB-S2X
LL protocol
GSE
Network layer protocol
IP
IP packets size
1344 bytes
NC overhead
10 bytes
Number of services
1


### 5.1. DVB-S2 versus DVB-S2X

In Figure1 we have shown, per each location within a beam, the achievable MODCODs when decoding signals from the own and from the adjacent beam providing better SINR for a DVB-S2 based system. Figure 5 shows the equivalent information for a DVB-S2X based system. It can be observed that at the same locations the UTs under DVB-S2X can be assigned MODCODs with higher spectral efficiencies (a combination of higher modulation order and higher code rate) from the own and adjacent beam providing better SINR and hence an improved throughput can be delivered.Figure 5
Example of achievable MODCODs for locations within a beam of interest. Results extracted from a 70-beam, 4-coloured system under DVB-S2X.Moreover, DVB-S2X incorporates natively the aforementioned channel-bonding tool for multitune receivers to split and uplink a single stream into several satellite transponders (i.e., different beams) and to receive it with a DVB-S2X receiver implementing several physical layers. Hence, DVB-S2X seems to be a more appropriate standard to support our scheme.


### 5.2. Clear-Sky Conditions

Figures6 and 7 show the average multicast throughput for 6 relevant UTs distributions, generated as in Section 2.4, for a cluster C of 7 beams.Figure 6
Average multicast throughput for 3-UT distribution within a cluster of beams;α=50, R0=RC, and x0, y0 = {beam center, close to the edge, overlapping area}.Figure 7
Average multicast throughput for 3-UT distribution within a cluster of beams;α=1000, R0=RC, and x0, y0 = {beam center, close to the edge, overlapping area}.When UTs tend to be distributed over a large area, as in Figure6, the proposed C-PFM scheme can achieve gains in average system multicast throughput up to the 27%. Such advantage is obtained when most of the UTs are located in, or around, the overlapping area of the beam because the scheme enables the transmission with multilink reception and NC. Note that the PFM scheme, which does not incorporate multilink reception, attains a similar throughput to the WCM scheme. As most of the UTs tend to be located further from the overlapping area, that is, close to the beam edge or in the beam center, the gain obtained decreases to 16% and 3%, respectively. This is mainly produced because the multilink reception capability with NC is not activated that often since it is not favourable to the UTs located in the center of the beam. Even if some MODCOD balancing is achieved due to the PF rule, it is not high enough to provide significant average system multicast throughput gains.When UTs are concentrated in smaller areas, as Figure7 shows, the gains tend to be larger. If the UTs are located in the beam overlapping area, the gain in average system multicast throughput is 88% and it gradually decreases to 64% and 5% when UTs are highly concentrated close to the edge of the beam or in the beam center, respectively. Ideally, if all the UTs were located in the beam overlapping area (e.g., in a regional satellite broadcast for a city), a 100% throughput gain could be achieved, except for the 10 bytes overhead introduced when including the encoding coefficients in the network coded IP packets.Hence, the proposed C-PFM scheme provides higher gains when UTs are located in the beam overlapping area or close to the edge of the beam concentrated in small (64%–88%) or large areas (16%–27%). When UTs are located in the center of the beam, the gain is marginal (3%–5%). Moreover, this gain is achieved without employing additional resources, by just enabling UTs accessing orthogonal transmissions already in place.


### 5.3. Service Level Agreement

The proposed C-PFM scheme optimizes the average multicast throughput of the system according to the PF rule but it does not guarantee a minimum throughput to each individual UT. This side effect of the PF rule is reflected in the right bottom plot of Figure7 where the worst UT is assigned a negligible throughput; that is, an SLA is not guaranteed to the UT. Many multicast applications are based on the delivery of video and audio such as video streaming, audio streaming, or online gaming. The SLA guarantees that a minimum amount of data is delivered to the UT such that these services are delivered properly in a similar way traditional Digital Subscriber Line (DSL) operators guarantee a minimum connection quality to users in disfavoured areas while users in better areas get better connection quality.The results in Figure8 consider a GW providing multicast services to Spain. The number of UTs subscribed to the multicast service is proportional to the number of inhabitants of the 5 largest Spanish metropolitan areas [35]. We consider a minimum SLA with the UT equivalent to half of the throughput achieved with the highest common achievable MODCOD from the own beam (see top plot of Figure 5); that is, in average each UT should be scheduled once every 2 slots with such MODCOD. As it can be observed in the right bottom plot of Figure 8, without an SLA mechanism such reference level is not reached (the blue bars mark the minimum throughput per UT and the worst UT under the C-PFM scheme does not overcome it). A simple mechanism to provide a minimum SLA consists in periodically introducing slots with a MODCOD decodable by all the UTs in the cluster. In particular, the throughput assigned to the worst UT is averaged for a past number of slots and in case it does not reach the proposed level, a MODCOD decodable by such UT is imposed. As it can be observed, now the worst UT is guaranteed the SLA but the average multicast throughput of the cluster has been affected since at certain slots the optimal MODCODs are not scheduled. In the case under study, the provision of SLA has reduced the average system multicast throughput gain compared to WCM from 32% to 21.5%.Figure 8
Average multicast throughput for realistic UTs distribution in Spain.The fact that each transmission is decodable by a different number of UTs due to the use of ACM should not affect the proper delivery of multicast services proved that proper techniques are employed at the application layer. For audio and video streaming and online gaming Scalable Video Coding (SVC) [36] and Scalable Audio Coding (SAC) [37] techniques can be used. For instance, Figure 9 shows the spectral efficiency of the scheduled MODCODs when the SLA is guaranteed. When the MODCOD decodable by all the UTs is transmitted, a video/audio base layer (BL) could be sent such that all the UTs retrieve it. When MODCODs decodable only by a subset of the UTs are transmitted, different levels of video/audio enhancement layers (ELs) can be sent.Figure 9
Spectral efficiency transmitted in each time-slot in the C-PFM scheme for UTs distribution in Spain and considering a minimum SLA. Green dots signal time-slots where the multilink with NC capability is not enabled. Red dots signal time-slots where the multilink with NC capability is enabled.For other multicast applications like file distribution and downloading which are less time sensitive, UTs able to decode all the MODCODs would retrieve the files faster and then abandon the multicast group. Then, the scheme would adapt the transmitted MODCODs according to this new set of UTs which could complete the download of the file a little later in time.


### 5.4. Rain Conditions

Finally, we show that our scheme is also advantageous under rainy conditions while it is able to maintain the SLA guarantee to the UT. To this aim, we assume that the southeast coast of Spain within the dotted ellipse in the left bottom plot of Figure8 is under rain event. Such event produces a 5 dB attenuation. As it can be observed in the right bottom plot of Figure 8, the C-PFM scheme still maintains a 15% gain in average multicast throughput compared to the WCM scheme while keeping the proposed SLA.


## 5.1. DVB-S2 versus DVB-S2X

In Figure1 we have shown, per each location within a beam, the achievable MODCODs when decoding signals from the own and from the adjacent beam providing better SINR for a DVB-S2 based system. Figure 5 shows the equivalent information for a DVB-S2X based system. It can be observed that at the same locations the UTs under DVB-S2X can be assigned MODCODs with higher spectral efficiencies (a combination of higher modulation order and higher code rate) from the own and adjacent beam providing better SINR and hence an improved throughput can be delivered.Figure 5
Example of achievable MODCODs for locations within a beam of interest. Results extracted from a 70-beam, 4-coloured system under DVB-S2X.Moreover, DVB-S2X incorporates natively the aforementioned channel-bonding tool for multitune receivers to split and uplink a single stream into several satellite transponders (i.e., different beams) and to receive it with a DVB-S2X receiver implementing several physical layers. Hence, DVB-S2X seems to be a more appropriate standard to support our scheme.


## 5.2. Clear-Sky Conditions

Figures6 and 7 show the average multicast throughput for 6 relevant UTs distributions, generated as in Section 2.4, for a cluster C of 7 beams.Figure 6
Average multicast throughput for 3-UT distribution within a cluster of beams;α=50, R0=RC, and x0, y0 = {beam center, close to the edge, overlapping area}.Figure 7
Average multicast throughput for 3-UT distribution within a cluster of beams;α=1000, R0=RC, and x0, y0 = {beam center, close to the edge, overlapping area}.When UTs tend to be distributed over a large area, as in Figure6, the proposed C-PFM scheme can achieve gains in average system multicast throughput up to the 27%. Such advantage is obtained when most of the UTs are located in, or around, the overlapping area of the beam because the scheme enables the transmission with multilink reception and NC. Note that the PFM scheme, which does not incorporate multilink reception, attains a similar throughput to the WCM scheme. As most of the UTs tend to be located further from the overlapping area, that is, close to the beam edge or in the beam center, the gain obtained decreases to 16% and 3%, respectively. This is mainly produced because the multilink reception capability with NC is not activated that often since it is not favourable to the UTs located in the center of the beam. Even if some MODCOD balancing is achieved due to the PF rule, it is not high enough to provide significant average system multicast throughput gains.When UTs are concentrated in smaller areas, as Figure7 shows, the gains tend to be larger. If the UTs are located in the beam overlapping area, the gain in average system multicast throughput is 88% and it gradually decreases to 64% and 5% when UTs are highly concentrated close to the edge of the beam or in the beam center, respectively. Ideally, if all the UTs were located in the beam overlapping area (e.g., in a regional satellite broadcast for a city), a 100% throughput gain could be achieved, except for the 10 bytes overhead introduced when including the encoding coefficients in the network coded IP packets.Hence, the proposed C-PFM scheme provides higher gains when UTs are located in the beam overlapping area or close to the edge of the beam concentrated in small (64%–88%) or large areas (16%–27%). When UTs are located in the center of the beam, the gain is marginal (3%–5%). Moreover, this gain is achieved without employing additional resources, by just enabling UTs accessing orthogonal transmissions already in place.


## 5.3. Service Level Agreement

The proposed C-PFM scheme optimizes the average multicast throughput of the system according to the PF rule but it does not guarantee a minimum throughput to each individual UT. This side effect of the PF rule is reflected in the right bottom plot of Figure7 where the worst UT is assigned a negligible throughput; that is, an SLA is not guaranteed to the UT. Many multicast applications are based on the delivery of video and audio such as video streaming, audio streaming, or online gaming. The SLA guarantees that a minimum amount of data is delivered to the UT such that these services are delivered properly in a similar way traditional Digital Subscriber Line (DSL) operators guarantee a minimum connection quality to users in disfavoured areas while users in better areas get better connection quality.The results in Figure8 consider a GW providing multicast services to Spain. The number of UTs subscribed to the multicast service is proportional to the number of inhabitants of the 5 largest Spanish metropolitan areas [35]. We consider a minimum SLA with the UT equivalent to half of the throughput achieved with the highest common achievable MODCOD from the own beam (see top plot of Figure 5); that is, in average each UT should be scheduled once every 2 slots with such MODCOD. As it can be observed in the right bottom plot of Figure 8, without an SLA mechanism such reference level is not reached (the blue bars mark the minimum throughput per UT and the worst UT under the C-PFM scheme does not overcome it). A simple mechanism to provide a minimum SLA consists in periodically introducing slots with a MODCOD decodable by all the UTs in the cluster. In particular, the throughput assigned to the worst UT is averaged for a past number of slots and in case it does not reach the proposed level, a MODCOD decodable by such UT is imposed. As it can be observed, now the worst UT is guaranteed the SLA but the average multicast throughput of the cluster has been affected since at certain slots the optimal MODCODs are not scheduled. In the case under study, the provision of SLA has reduced the average system multicast throughput gain compared to WCM from 32% to 21.5%.Figure 8
Average multicast throughput for realistic UTs distribution in Spain.The fact that each transmission is decodable by a different number of UTs due to the use of ACM should not affect the proper delivery of multicast services proved that proper techniques are employed at the application layer. For audio and video streaming and online gaming Scalable Video Coding (SVC) [36] and Scalable Audio Coding (SAC) [37] techniques can be used. For instance, Figure 9 shows the spectral efficiency of the scheduled MODCODs when the SLA is guaranteed. When the MODCOD decodable by all the UTs is transmitted, a video/audio base layer (BL) could be sent such that all the UTs retrieve it. When MODCODs decodable only by a subset of the UTs are transmitted, different levels of video/audio enhancement layers (ELs) can be sent.Figure 9
Spectral efficiency transmitted in each time-slot in the C-PFM scheme for UTs distribution in Spain and considering a minimum SLA. Green dots signal time-slots where the multilink with NC capability is not enabled. Red dots signal time-slots where the multilink with NC capability is enabled.For other multicast applications like file distribution and downloading which are less time sensitive, UTs able to decode all the MODCODs would retrieve the files faster and then abandon the multicast group. Then, the scheme would adapt the transmitted MODCODs according to this new set of UTs which could complete the download of the file a little later in time.


## 5.4. Rain Conditions

Finally, we show that our scheme is also advantageous under rainy conditions while it is able to maintain the SLA guarantee to the UT. To this aim, we assume that the southeast coast of Spain within the dotted ellipse in the left bottom plot of Figure8 is under rain event. Such event produces a 5 dB attenuation. As it can be observed in the right bottom plot of Figure 8, the C-PFM scheme still maintains a 15% gain in average multicast throughput compared to the WCM scheme while keeping the proposed SLA.


## 6. Conclusions

In this paper we have proposed a full multicasting scheme for multibeam satellite systems. Our scheme makes use of the ACM feature of DVB-S2 and DVB-S2X standards by employing a PF scheduling policy and opportunistically enables UTs, assumed to have multilink reception capabilities and a NC decoder, to retrieve information from the adjacent beam providing better SINR. In particular, we have introduced the functional architecture logic in the GW, based on a cross-layer design, and UT side. We have also detailed the algorithms employed in the novel blocks in the GW and we have designed the packet scheduling architecture respecting the IETF Diffserv model at the IP level.From the complexity point of view, the practical introduction of our scheme requires modifications at both the GW and the UT. While modifications in the GW can be implemented via software upgrades, the satellite UT requires further modifications. The RF chain must be able to syntonize two different bands in order to receive from the own and adjacent beam. Also, two separate decoding/deencapsulating chains are necessary up to the IP level. Finally, the UTs need a NC decoder. Although these elements may build up an apparent complex receiver, we have demonstrated that the standards and technology for its implementation are already in place.From the results point of view, the average system multicast throughput gain provided by our scheme highly depends on the UTs distribution. With the aim of providing fair and realistic results, our simulations have been carried out over nonuniform theoretical and realistic distributions within clusters of few beams where transmissions employ the recently appeared DVB-S2X standard. The results obtained show the following. When the UTs tend to be located in the centers of the beams, the gain provided by our scheme is marginal, that is, below 5%. However, when most of the UTs are located close to the edge of a beam or in the overlapping areas with another beam, even if the UTs spread over a large area of the cluster, the gain provided by our scheme can reach the 88%. This gain is mainly enabled by the multilink reception and NC.Moreover, a simple SLA mechanism has been introduced and simulated such that a minimum service is guaranteed to the UT while still providing significant gains (21%) even under rain conditions (15%). In addition, we have suggested, for some relevant multicast applications, proper techniques at the upper layers to deal with the different amounts of data delivered to each UT. Finally, the gains achieved by our scheme do not require the use of additional resources but only enable UTs accessing orthogonal transmissions already in use by UTs in other beams.As future line of research, we would like to investigate and analyze the performance of the proposed scheme in a mobile environment where the election of a suitable modulation and codification becomes more challenging since UTs are subject to fast channel variations. Moreover, our intraservice policy could be combined with interservice policies in order to meet the different requirements demanded by multicast services.



---

*Source: 364234-2015-07-09.xml*