In all scenarios we first evaluate the case where only BSs are active, i.e., RIS are not active. In this case the coverage area of the BSs must include almost the whole scenario area. Large blue circles are illustrating this situation, for every scenario.
The total power transmitted by each BS or device depends on the scenario and the frequency carrier. At 100 GHz the power was set to 30mW (14.8dBm) for the Indoor Factory, whereas 3.16W (35dBm) or 250mW (24dBm) are transmitted from each BS at 28 GHz, in the case of Stadium. For the Smart City scenario, the power transmitted by IoT devices was set as 10mW (10dBm) both for 3.6 GHz and 28 GHz. Two types of links are considered, namely, a direct link between BS and devices or users and an indirect link through the RIS. For all scenarios, because of C-RAN there are always double-links, where devices or users are served simultaneously by a RIS and a BS. The noise power is N0= -85 dBm for the bandwidth B=0.4 GHz and N0= -94 dBm for the bandwidth B=50 MHz. The receiver noise figure F= 3dB. The spacing between each element of the RIS panel is dRIS= λ/2=42 mm (3.6GHz), dRIS= λ/2=5.4 mm (28GHz) and dRIS=1.5 mm (100GHz), resulting respectively in areas of A=1764 mm2, A=29.16 mm2 and A=2.25 mm2 per element. The gains of the individual antenna elements of the arrays are 0 dBi for both the transmitter and receiver.
4.2. System-Level Simulations
The signal-to-noise ratio (SNR) in dB considered in the system-level simulations is obtained from SNR = (Es/N0) + 10log(Rs/B) dB, where Rs is the total transmitted symbol rate per antenna and user, B is the total bandwidth (we considered 50 MHz at 3.6 GHz and 28 GHz for Smart City with numerology 0, and 400 MHz at 28 GHz and 100 GHz for the other scenarios with numerology 3), and Es/N0 is the ratio of symbol energy to noise spectral density in dB. Values of Es/N0 are obtained from the link-level BER results. In this work we only considered QPSK modulated symbols. 5G NR frame structure has both frequency division duplex (FDD), used in the paired spectrum, and time division duplex (TDD), used for the unpaired spectrum. We chose TDD in this work for all scenarios.
The subcarrier spacing, transmission time interval, cyclic prefix (CP), and the number of symbols per slot, are all defined by the proposed 5G NR scalable OFDM numerology [
15]. Our scenarios were simulated considering numerologies 0 and 3. Numerology 0 (4G) was chosen for Smart City scenario because the most important key performance indicator (KPI) is the coverage of IoT devices, not their throughput, and numerology 3 for the Stadium with Pitch and Indoor Factory scenarios, in mmWave and THz bands, respectively.
For Smart City scenario, we considered different UPA antenna arrays with
Nrx= 144, 256 up to 576 elements at the receiver side (BS). For the device (transmitter side)
Ntx=1,
Ns=1 and
Nc=1. For numerology 0, Δf=15KHz, every PRB has 12 carriers (B
prb=180KHz) with 14 symbols per subframe duration Ts = 1 ms, 168 (12 × 14) subcarriers. In real cellular networks with massive MIMO, there is a coherence interval where channel information does not change much, and pilot reference symbols must be retransmitted. The coherence interval is
where coherence time
Tc = 50 ms is well fitted for nomadic and static low-power IoT devices [
32] and
Bc=1.697MHz (9.42857x180KHz) is adequate to an urban city scenario such as Smart city. In
τc the maximum number of subcarriers is
Nsc=79200 (168 x 50 x 9.42857). Half of the coherence interval can be used for uplink data transmission, and the remaining half of the coherence interval can be used for downlink transmission. We need orthogonal reference pilots in uplink that must be subtracted from uplink data symbols. For uplink data and reference pilots, we must have
Nsc/4=79200/4 = 19800 subcarriers. Assuming that we divide 19800 into twelve groups (scenario has 4 BS x 3 sectors), then we obtain 19800/12=1650 subcarriers, one for each device transmitting in each sector at the same time without any intra cell interference.
Table 4 presents additional simulation parameters of the Smart City scenario. In real situation the number of active IoT devices per sector can be much higher because of the sleeping time of devices. Sensors do not need to transmit packets every 1 ms.
For the Smart City scenario, we consider only the uplink where devices are transmitting to BS and RIS at the same time. Uplink power control is appropriately applied so that the received powers are the same for almost all devices. However, data sent from devices that are too far from the receivers, i.e., BSs and/or RISs, might not be well received. There is a minimum received power that must be above the receiver sensitivity so that packets are correctly received. We take as reference the BER ≤ 5x10
−4 to decide, which packets are received without bit error. For numerology 0 of 5G NR, Δf=15KHz and 168 subcarriers are being transmitted in every PRB during 1 ms. Minimum throughput of each device is 2 Kbps, assuming transmission in one carrier with QPSK symbols, and minimum spectral efficiency per device is ℇ
u=2/15 bps/Hz. There is a theoretical equation for the maximum throughput, assuming perfect power control and single cell with no inter-cell interference [
33], which is:
where τ
p is fraction of time for pilot symbols transmission and B is the total bandwidth. SINR is the signal to interference plus noise ratio.
In
Figure 7, it is shown the throughput performance versus number of users, for different number of antennas at BS in each sector of scenario, and for different frequency carriers. The throughput curves of
Figure 7 were obtained with a maximum transmitted power of 10 mW. The settings are described in
Table 1. In this scenario we used
Ns=1,
Nc=1,
Ntx=1 and
Nrx=144, 256 or 576 antennas. As expected, the frequency carrier of 3.6 GHz provides the highest throughput compared to 28 GHz due to its lower path attenuation. Increasing the number of UPA antennas allows serving effectively more devices at the same time. For almost all curves, as
Nrx increases, the throughput increases up to a certain point, and then it decreases. Increasing
Nrx provides higher throughput. The decrease of throughput after a certain
Nu has to do with higher number of devices at the border of the cells with SINR below the reference to avoid errors in packet. Due to the C-RAN operation mode almost no reference pilot symbols contamination occurs. For 3.6 GHz the maximum simulated spectral efficiency occurs for
Nu=19800 and
Nrx=576, corresponding to ℇ=508/24.75=20.53 bps/Hz, where 508 Mbps is the aggregate throughput of 19800 devices and B=24.75 MHz is the occupied bandwidth. Theoretical maximum spectral efficiency is ℇ
max=554/24.75=22.40 bps/Hz.
In
Figure 8, it is shown the throughput performance versus
Nu, for different number of antennas at BS in each sector, for carrier frequency of 28GHz, without and with RIS panels
Nris = 576. We compare four standard communication curves without any RIS, only with BSs and devices, with three curves with RISs panels spread around the scenario. The former curves consist only of direct links between the devices and the receivers. The other cases consist of a combination of direct link connections and RIS-aided connections. They are represented as a percentage of devices that are transmitting signals to BS plus the percentage of devices also with RIS connections, namely as %BS + %RIS. For example, 100%BS+0%RIS, means that all devices are attached to the BSs, whereas 89%BS+11%RIS, represents 89% of spread out devices linked to nearest BS and the 11% remaining devices are linked to nearest RISs. C-RAN processes data received from the two best links involving each device. As expected, the introduction of RIS panels at cell borders provides higher throughput compared to BS only connections. As expected, increasing the number of receiving UPA antennas results in higher throughput due to their higher multiple access capability. The throughput gain introduced by RIS panels with
Nris = 576 increases with decreasing number of UPA antennas
Nrx = 576, 256 and 144. The gain is more noticeable for increasing
Nu because more users become close to where RIS panels are localized at the cell borders. Only for 28GHz we have considered an UPA with
Nrx = 1024 antennas because of its size. As expected, the maximum throughput is achieved for
Nrx = 1024 due to its higher multiple access capability.
Figure 9 presents the coverage versus transmitted power for
Nu=19800 with the same conditions of
Figure 8. The comparison between
Figure 9 and
Figure 8 shows that there is a direct correspondence between the throughput performance and the associated coverage. For the maximum transmitted power of 10 mW, the smallest coverage of 44% is achieved by the curve 100%BS+0%RIS,
Nrx = 144, whereas the highest coverage is 93% which is obtained by the curve 100%BS+0%RIS,
Nrx = 1024. Thus, the maximum coverage gain is 111%. The next highest coverage of 90% corresponds to the curve 89%BS+11%RIS,
Nrx = 576 and represents a coverage gain of 104% that includes a reduction of the number of antennas but the introduction of RIS. Still, the main contribution to the coverage gain is due to the number of antennas that is 97%.
Next, we test the Stadium with Pitch scenario with different parameters. For this purpose, we adopted a higher number of subcarriers Nc combined with more spatial streams and larger transmit powers. 5G NR numerology is 3, with Δf=120KHz, the bandwidth is 400MHz and 112 symbols are transmitted every millisecond. The maximum number of subcarriers is 3300. In this scenario the main KPI is throughput associated to the number of connections per unit of area. In all simulations in the Stadium the transmitter side BS has Ntx= 64, whereas the receiver side user has the same number of antennas namely, Nrx= 16. Different numbers of transmitted symbols per subcarrier were considered in the simulations, Ns=2 and Ns=3. The number of OFDM subcarriers (Nc) used in the evaluations is variable, namely, Nc = 12, 60, 120, 132, 216 or 240. All these numbers are multiples of 12, which corresponds to the number of subcarriers of a PRB in 5G NR.
Users are placed uniformly around TRPs localized close to the four corners of Stadium turf. 46% of users are on turf while the remaining 54% are sit on Stadium stands (see
Figure 4 and
Table 2). The configuration with 100%BS+0%RIS has 4 BS and 12 sectors and matches a typical cellular system. In this case users are uniformly distributed within a radius of 57 m (see
Table 2). On the other hand, when simulating cases with BSs and RIS operating simultaneously, users that are connected to BSs will be uniformly distributed within a circle with a radius of 34 m. RIS panels are uniformly distributed inside the scenario and users served by RISs are distributed uniformly within a radius of 17 m.
In
Figure 10, it is shown the throughput performance versus number of users in the Stadium scenario. The throughput curves of
Figure 10 were obtained with a maximum transmitted power of 250 mW. We keep the maximum number of subcarriers of 3300 to avoid intra cell (sector) interference. When we consider 11 PRB per user, i.e.,
Nc=132 the maximum number of users per sector is 3300/132=25. In the scenario there are 12 sectors, 4 BS each one with 3 sectors. The total number of served users at the same time in the Stadium is 25x12=300. Two different RIS sizes appear since each RIS panel is divided into sub-panels when serving more than 1 user. We have
NRIS=192 and 576 elements. When users are served with a direct link only from BS (100%BS+0%RIS), we have a total average throughput of approximately 48 Gbps for 300 users (represented as a red line). In the Stadium there are 22x3=66 RIS panels and 4x3=12 TRPs. Considering the black and blue lines, when users at cell borders start to be served by RIS panels, 66 users from 300 (i.e., 66/300=0.22) must be subtracted resulting in 300-66=234 users (i.e., 234/300=0.78) are served by 12 TRPs, 78%BS+22%RIS. For the pink and green lines, when 66x3=198 users are served by RIS 198/300=0.66, the remaining 300-198=102 users (i.e., 102/300=0.34) are served by TRPs, 34%BS+66%RIS. We observe that the throughput gain of curves with
Ns=3 compared to
Ns=2, is 1.5. There is no difference of throughput in users served by RIS (0%BS) compared to users served by BS (0%RIS). If more users are served by RIS panels (66%RIS) there is a throughput decrease compared to curves with few users (22%RIS) that is only noticeable for
Ns=3 (curves black and pink). For
Ns=2 there is no decrease of throughput when more users are served by the same RIS panel (curves blue and green). Note that in all curves there is no decrease of throughput for increasing
Nu because C-RAN operation combines transmitted signal power to RIS and BS, decreasing substantially inter cell interference. We are also assuming perfect CSI in link level results and there is no intra cell interference as explained above.
In
Figure 11, it is also shown the throughput performance versus
Nu in the Stadium scenario, with two different maximum transmitted powers from BSs. Low Power BSs (LPBS) transmit 250 mW and High Power BSs (HPBS) transmit 3.16W. We keep the maximum number of subcarriers of 3300 to avoid intra sector interference. Assigning 10 PRB per user, i.e.,
Nc=120, and for 20 PRB per user,
Nc=240. With 10 PRB to each user, 3300/120=27.5 is the number of users per sector. The total number of served users at the same time in the Stadium is 27.5x12=330. With 20 PRB to each user, there are 3300/240=13.75 users per sector. The total number of served users at the same time in the Stadium is 13.75x12=165. When users are exclusively served by the TRPs of BS with direct links 100%BS+0%RIS, we have a total throughput of approximately 53 Gbps for 330 users (pink line), and a total throughput of 48 Gbps for 165 users (red line). Considering the case of 10 PRBs (blue lines), when 66 users at cell borders are served by RIS (i.e., 66/330=0.20), then the remaining 330-66=264 users are served by TRPs, (i.e., 264/330=0.80), 80%BS+20%RIS,10RB. When 66x2=132 users are served by RIS (i.e., 132/330=0.40), the remaining 60% users are served by TRP 60%BS+40%RIS,10RB. Assigning 20 PRB per user (black lines), when 66 users at cell borders are served by RIS (i.e., 66/165=0.40), then 165-66=99 users are served by TRPs, 99/165=0.60, 60%BS+40%RIS,20RB. When 66x2=132 users are served by RIS (i.e., 132/165=0.80), then the remaining 165-132=33 users, (i.e., 33/132=0.20) are served by TRP 20%BS+80%RIS,20RB. We observed only a very slight decrease of the throughput results by increasing the number of users from 1 to 2 served by every RIS panel. There is almost no difference in achieved throughput between HPBS and LPBS. There is a slightly higher throughput for LPBS that is more noticeable for 20 RB. Throughput gains are between 97% for 10RB and 110% for 10RB.
Figure 12 illustrates the average coverage versus transmitted power for both
Nc=120 (10 PRB) and
Nc=240 (20 PRB), corresponding to the throughput performance presented in
Figure 11 considering only HPBS, with a maximum transmitted power of 3.16 W (35dBm). The comparison between
Figure 12 and
Figure 11 indicates that there is a direct correspondence between the throughput performance of
Figure 11 and the associated coverage of
Figure 12. However, in terms of coverage the differences are more noticeable than with throughput. It is noticeable that for the maximum transmitted power, the highest coverage of 99.4% is achieved by the curve 80%BS+20%RIS,10RB followed by the curve 60%BS+40%RIS,10RB. Next follow the curves 60%BS+40%RIS,20RB and 20%BS+80%RIS,20RB. It is obvious that when the number of PRBs per user doubles the coverage decreases. It is also noticeable that doubling the number of users served by RIS panels also decreases the coverage. The lowest coverage values correspond to no RIS panels activated, i.e., 100%BS+0%RIS, 10RB and finally 100%BS+0%RIS,20RB with a coverage of 45.4%. The coverage gain compared to the highest coverage value with 20RB is 111% and 98% for 10RB.
In
Figure 13, we consider again the Stadium scenario, presenting the aggregate throughput for each sector, considering the maximum transmit power of 24 dBm, versus the number of PRBs per user. The maximum number of subcarriers is set as 3300 per sector independently of the number of PRBs per user. We evaluate the sector throughput only for some specific number of PRBs, namely, 1PRB, 5PRB, 10PRB, 15PRB and 20PRB. Assigning 1 PRB per user corresponds to
Nc=12, 5 PRB per user corresponds to
Nc=60, 10 PRB per user corresponds to
Nc=120, 15 PRB per user corresponds to
Nc=180 and finally 20 PRB per user requires
Nc=240 subcarriers. The number of users per sector when each user has 1 PRB is 3300/12=275. In total, there are 12 sectors in the Stadium, thus a total of 3300 active users. For the case where every active user has 20 PRB then there are 3300/240=13.75 users per sector and the total of active users in the Stadium is 165. We observe in
Figure 13 that for the curve BS (100%BS+0%RIS) there is a maximum throughput for 1 PRB, which is the same for the curve BS+RIS. For 5 PRB and 10 PRB per user there is a clear linear decrease of throughput of the red curve (BS) compared to black curve BS+RIS. The decrease of throughput continues with 15 PRB and 20 PRB but in a slower way. The curve BS+RIS has almost the same aggregate throughput per sector, however, it is noticeable a very slight decrease with increasing number of PRBs. When users are transmitting a higher number of PRBs the size of packets transmitted is higher and more likely to suffer deep fades. In this case, the link diversity RIS+BS provided by the C-RAN is able to cope with the expected decrease of throughput.
The third evaluated scenario is the Indoor Factory, which is the only one with a higher frequency band at sub-THz. The system operates at 100 GHz which, due to its shorter wavelength, allows us to work with more elements at the RIS, i.e.,
NRIS=576 up to 2304. The maximum transmit power of BS is 30 mW (14.8dBm). Users are placed uniformly around 2 sites with TRPs equipped with UPA antennas per sector, localized 75 meters close to the walls of the Indoor Factory (see
Figure 5 and
Table 3). We have considered high density clutter inside the Factory, with low BS antenna height and high BS antenna height. For this scenario we keep 5G NR numerology 3, with Δf=120KHz, a bandwidth of 400MHz and 112 symbols transmitted in every millisecond, such as in the Stadium scenario. Thus, it is possible to compare directly the two scenarios in spite of having different frequency bands and areas. In all simulations in the Factory environment, each transmitting TRP has a UPA with
Ntx= 256 antennas, and the user receiver side, has
Nrx= 16. Different numbers of transmitted symbols per subcarrier were considered in the simulations and
Ns=3. When there are no active RIS panels in the Factory there are only 2 BS, each with 3 sectors, making a total of 6 sectors. In this case, users are uniformly distributed within a radius of 70 m (see
Table 3). When simulating performance curves with BSs and RIS operating simultaneously, users that are connected to BSs will be uniformly distributed within a radius of 50 m. 16 RIS panels are uniformly distributed inside the scenario and users served by them are distributed within a radius of 25 m.
In
Figure 14, it is shown the throughput performance versus
Nu in the Factory scenario, with low BS antenna height, for two different resource blocks, and RIS panels. The scenario parameters are
Ntx = 256,
Nrx = 16,
Ns = 3. An assignment of 11 PRB per user, gives
Nc=132, and 25 PRB per user, gives
Nc=300. With 11 PRB per user, 3300/120=25 is the number of users per sector. With 25 PRB per user, 3300/3000=11 is the number of users per sector. As there are 6 sectors in the Factory, the total number of users is 150 or 66 for 11 PRB and 25 PRB, respectively. When users are exclusively served by the TRPs of BS with direct links 100%BS+0%RIS, we have a total throughput of approximately 41 Gbps for 150 users (pink line), and the total throughput of 37 Gbps for 66 users (red line). In the Factory there are 16x3=48 RIS panels and 2x3=6 TRPs (one per sector). When users at cell borders start to be served by RIS panels, 48 users must be subtracted resulting in 150-48=102 users served by 6 TRPs. The 48 users are served by RIS 48/150=0.32 (blue lines), and the remaining 68% users are served by TRPs 68%BS+32%RIS,11RB. In the case of 25 PRB (black lines) the total number of users is 66, we have 48 users served by RIS (i.e., 48/66=0.73) and the remaining 66-48=18 users served by TRPs (i.e.,18/66=0.27), 27%BS+73%RIS,25RB. We observe that by increasing the number of RIS panels from 576 up to 2304 only increases slightly the throughput. This increment is more noticeable for 25RB. For the case of maximum number of users there is almost no difference in the achieved throughput between lines 27%BS+73%,25RB and 68%BS+32%RIS,11RB. When compared to lines 100%BS+0%RIS, the throughput gains are between 29% for 11 PRB and 53% for 25 PRB.
Figure 15 presents the average coverage versus transmitted power for both
Nc=132 (11 PRB) and
Nc=300 (25 PRB), with maximum transmitted power of 30 mW (14.8dBm), corresponding to the throughput performance illustrated in
Figure 14. As expected, there is a direct correspondence between the throughput performance of
Figure 14 and the associated coverage of
Figure 15. In terms of performance of coverage, the differences between 11 PRB, or 25 PRB per user are higher than with throughput. We observe that for 11 PRBs with the maximum transmitted power, the highest coverage of 99% is achieved for both 576 and 2304 RIS panels. However, the coverage for 25 PRBs per user is higher when the RIS panels have 2304 elements compared to 576 elements, not only for 30 mW but also for all transmitted power interval. The lowest coverage values are 67% and 77% corresponding both to no activated RIS panels, i.e., 100%BS+0%RIS,25RB and 100%BS+0%RIS,11RB, respectively. The corresponding coverage gains are between 29% for 11 PRB and 48% for 25 PRB. This is in accordance with throughput values of
Figure 13. The comparison of the average coverage of Factory scenario (
Figure 15) is slightly higher than the coverage of Stadium scenario (
Figure 12).
In
Figure 16, we consider the Factory scenario, illustrating the aggregate throughput for each sector, versus the number of PRBs per user, considering the maximum transmit power of 14.8 dBm. We evaluate the sector throughput for some specific number of PRB, namely, 1PRB, 5PRB, 11PRB, 18PRB and 25PRB. With 1 PRB per user results
Nc=12, with 5 PRB per user results
Nc=60, with 11 PRB per user results
Nc=132, with 18 PRB per user results
Nc=216 and finally 25 PRB per user requires
Nc=300 subcarriers. The number of users per sector when each user has 1 PRB is 3300/12=275. In total, there are 6 sectors in the Stadium, thus there is a maximum of 1650 active users. For the case where every active user has 25 PRB then there are 3300/300=11 users per sector and the total number of active users in the Factory is 66. We observe in
Figure 15 that for the red curve LBS (100%BS+0%RIS) where the BS UPA antenna is at 2 m height, there is a maximum throughput for 1 PRB, which is the same for both curves BS+RIS (black) and HBS (pink). For 5 PRB, 11 PRB, 18 PRB per user there is a clear linear decrease of throughput of the red curve compared to the black curve BS+RIS that has almost the same constant throughput. The decrease in throughput continues with 25 PRB but in a slower way. The pink curve HBS with the BS UPA antenna at 8 m height has also a decreasing aggregate throughput per sector, for an increasing number of PRBs. However, it is a very slight decrease with the increasing number of PRB. When users are served by BS at 8 m height there is almost no interference from the high clutter density because the height of clutter is 6 m. The reason why transmitting higher number of PRBs per user decreases the throughput has to do with packets of large size that are more likely to suffer deep fades. The link diversity of RIS+BS is provided by the C-RAN operation that is capable to avoid the decrease of throughput. Therefore, if possible the antennas of BS should be placed as high as possible to be above the height of high clutter density characteristic of Indoor Factory scenarios.
The comparison of
Figure 16 with
Figure 13 indicates that the aggregate throughput of each sector of the Factory is half of the throughput of the Stadium. This is explained by the number of BSs being half in the Factory compared to the Stadium. Therefore, if we need to increase the aggregate throughput of the Factory scenario, we have two options, such as 1) We increase the number of BS or 2) We use carrier aggregation. Indeed, 5G NR allows for aggregation of 16 carriers each one with 400 MHz. At 100 GHz it is possible to have available a total bandwidth of 6.4 GHz. The first solution involves densification of the network, and it is not possible to be energy efficient with this solution. To be energy efficient and achieve high sector throughput we need to aggregate carriers and place RIS surfaces distributed around the BSs of the scenario.