Abstract

This paper describes modeling of the available bit rate (ABR) source traffic in asynchronous transfer mode (ATM) network using BLPos/GTEXP traffic generator, which employs Poisson distribution for modeling the burst length (BLPos) and exponential distribution for modeling the gap time (GTEXP). This traffic generator inherits the advantages of both Poisson and exponential distribution functions to achieve enhanced link performance. Analytical and simulation results for BLPos/GTEXP traffic generator have been presented and compared.

1. Introduction

The Poisson process is an extremely useful process for modeling purposes in many practical applications, such as, for example, to model arrival processes for queueing models or demand processes for inventory systems. It is empirically found that in many circumstances the arising stochastic processes can be well approximated by a Poisson process. From the standpoint of the network users, the messages within a session are typically triggered by particular events. But from the standpoint of the network, these message initiations are somewhat arbitrary and unpredictable. Therefore, the sequence of times at which messages arrive for a given session is a random process. Further, the arrival process of one session has to be independent of others. The poisson process has been proved to be the most appropriate for such purpose and therefore has been considered in the present work.

In this paper source traffic modeling and simulation have been carried out for ABR service in ATM networks using Poisson distribution for modeling the burst length (BLPos) and exponential distribution for modeling the gap time (GTExp). The corresponding traffic generator is designated as BLPos/GTExp.

2. BLPos/GTEXP Traffic Generator

BLPos/GTEXP traffic generator generates cells sent at a fixed rate (ACR) during BL and no cells are sent during GT. BL is assumed to be Poisson distributed (BLPos) whereas GT is exponentially distributed (GTExp). BLPos/GTEXP traffic generator was created and parameterized as explained below.

When events occur according to an exponential distribution, they are said to occur completely at random. Thus the arrival time of the next event is not affected by the time elapsed since the previous event [1โ€“5]. An exponential random variable with a mean value of ๐œ†Exp (whereโ€‰โ€‰๐œ†Exp>0) is given by the following probability density function (PDF): ๐‘“(๐‘ฅ)Exp=โŽงโŽชโŽจโŽชโŽฉ๐œ†Expโ‹…eโˆ’๐œ†Expโ‹…๐‘ฅโŽซโŽชโŽฌโŽชโŽญ๐‘ฅโ‰ฅ00๐‘ฅ<0,(1) and the cumulative density function (probability that time between events is x) can be obtained by integrating (1): ๐น(๐‘ฅ)Exp=โŽงโŽชโŽจโŽชโŽฉ1โˆ’eโˆ’๐œ†Expโ‹…๐‘ฅโŽซโŽชโŽฌโŽชโŽญ๐‘ฅโ‰ฅ00๐‘ฅ<0.(2) The mean and variance of the exponential distribution are ๐ธ(๐‘ฅ)Exp=๎€œโˆž0๐‘ฅโ‹…๐œ†Expโ‹…๐‘’โˆ’๐œ†Expโ‹…๐‘ฅ1๐‘‘๐‘ฅ=๐œ†Exp,Var(๐‘ฅ)Exp=๎€œโˆž0๐‘ฅ2โ‹…๐œ†Expโ‹…๐‘’โˆ’๐œ†Expโ‹…๐‘ฅ๎‚ต1๐‘‘๐‘ฅโˆ’๐œ†Exp๎‚ถ2=1๐œ†2Exp.(3) A Poisson random variable ๐‘ฅ with parameter ๐œ‡Posโ€‰โ€‰(whereโ€‰โ€‰๐œ‡Pos>0) has PDF given by [1, 2] ๐‘“(๐‘ฅ)Pos=๐œ‡๐‘ฅPos๐‘’๐‘ฅ!โˆ’๐œ‡Pos,๐‘ฅ=0,1,2,โ€ฆ.(4) The mean and variance are identical for the Poisson distribution [1]: ๐ธ(๐‘‹)Pos=Var(๐‘‹)Pos=๐œ‡Pos,(5) where a Poisson distribution can be developed by involving the definition of a Poisson Process for a fixed time (๐‘ก), so that ๐œ‡Pos=๐œ†Expโ‹…๐‘ก. For Poisson distributed events, the time between successive events is exponentially distributed with mean 1/๐œ‡Pos. For an event that occurs at a time interval with an exponential distribution, the rate of occurrence of the event is Poisson distributed with mean ๐œ†Exp. Required distribution modeling involves a transformation function for converting a random variable of uniform distribution into the required distribution. Considering the fundamental transformation law of probabilities for two probability density functions ๐‘“(๐‘ฅ) and ๐‘(๐‘ข)||||=|||||||๐‘“(๐‘ฅ)๐‘‘๐‘ฅ๐‘(๐‘ข)๐‘‘๐‘ขor๐‘“(๐‘ฅ)=๐‘(๐‘ข)๐‘‘๐‘ข|||,๐‘‘๐‘ฅ(6) where ๐‘(๐‘ข) is the PDF of random variable ๐‘ข and ๐‘“(๐‘ฅ) is another PDF of random variable ๐‘ฅ. Since ๐‘ข is a random variable of a uniform distribution in the range 0 to 1, ๐‘(๐‘ข) is a constant (=1) and hence ๐‘“(๐‘ฅ)=๐‘‘๐‘ข๎€œ๐‘‘๐‘ฅandtherefore๐‘ข=๐น(๐‘ฅ)=๐‘ฅ0๐‘“(๐‘ง)๐‘‘๐‘ง.(7) Equation (7) can be used to find source random variable ๐‘ฅ=๐บ(๐‘ข) through inverse transformation of ๐‘ข=๐น(๐‘ฅ). For the required distribution, the inverse can easily be found from (7) with ๐‘“(๐‘ง) corresponding to the required distribution. ๐‘ข is uniformly distributed in the range 0โ‰ค๐‘ขโ‰ค1. It can be generated by using the function ๐‘Ÿ๐‘Ž๐‘›๐‘‘() provided by the standard Linux library or using Mersenne Twister (MT) [6].

For modeling the GTExp, (7) is used with ๐‘ข=๐น(๐‘ฅ)=1โˆ’๐‘’โˆ’๐œ†Expโ‹…๐‘ฅor1โˆ’๐‘ข=๐‘’โˆ’๐œ†Expโ‹…๐‘ฅ.(8) Therefore, the required transformation is 1๐‘ฅ=โˆ’๐œ†Expโ‹…ln(1โˆ’๐‘ข).(9) By changing the variables ๐‘ฅ and (1โˆ’๐‘ข) of (9) with ๐‘‹Exp and ๐‘ˆ respectively, we get ๐‘‹Exp1=โˆ’๐œ†Expโ‹…ln(๐‘ˆ).(10)

2.1. Estimation of the Load (๐ฟ๐‘–) for the Traffic Generator

The load variation of the traffic can be realized by synthesizing predefined load such that the resulting load โˆ‘๐ฟ=๐‘๐‘–=1๐ฟ๐‘–, where ๐ฟ๐‘– is the traffic load due to ๐‘–th source. Therefore, the aggregate traffic from ๐‘ sources will generate the load ๐ฟ on a link with rate ๐‘… Mbps giving average throughput of ๐‘…โ‹…๐ฟ Mbps. The load ๐ฟ๐‘– generated by an individual source can be expressed as ๐ฟ๐‘–=BLPosโ‹…๐พBLPosโ‹…๎€ท๐พ+๐‘ƒ๐‘Ÿ๎€ธ+GLExp,(11) where BLPos, GLExp, K, and ๐‘ƒ๐‘Ÿ are the mean BLPos, mean gap length, cell size, and minimum intercell gap length (Preamble), respectively in bytes, and then the load ๐ฟ๐‘– can be found from (11).

2.2. Estimation of the Minimum Gap Time (๐‘€GTExp)

The minimum GTExp(๐‘€GTExp) is a secondary parameter dependent on load. Given a desired load, ๐‘€GTExp is calculated by the source automatically using exponential distribution. Using (11) the GTExp can be expressed as GLExp=BLPosโ‹…๎‚ธ๐พโ‹…1โˆ’๐ฟ๐‘–๐ฟ๐‘–โˆ’๐‘ƒ๐‘Ÿ๎‚น.(12) The BLPos can be written as BLPos=๐ธ(๐‘ฅ)BLPos=๐‘€BLPos,(13) where ๐‘€BLPos is the minimum BLPos.

GLExp can be written as GLExp=๐ธ(๐‘ฅ)GLExp=๐‘€GLExp.(14) Substituting the values of BLPos and GLExp from (13) and (14), respectively, in (12) ๐‘€GLExp can be written as ๐‘€GLExp=๐‘€BLPosโ‹…๎‚ธ๐พ๐ฟ๐‘–โˆ’๎€บ๐พ+๐‘ƒ๐‘Ÿ๎€ป๎‚น.(15) Considering the link rate and using the following relation: ByteTime=ByteSize(bits)๐‘LinkRate(bits/sec)orByteTime=๐‘…,(16) the ๐‘€GTExp now can be computed as ๐‘€GTExp=๐‘๐‘…โ‹…๐‘€BLPosโ‹…๎‚ธ๐พ๐ฟ๐‘–โˆ’๎€บ๐พ+๐‘ƒ๐‘Ÿ๎€ป๎‚น.(17) Now the value of ๐‘ƒ๐‘Ÿ=1/ACR is readily available, depending upon the selected value(s) of ACR that can be separately taken as variable, and thus (17) can be rewritten as ๐‘€GTExp=๐พโ‹…๐‘๐‘…๐‘€BLPos๎‚ธ1๐ฟ๐‘–๎‚นโˆ’1.(18) Therefore, (18) can be used for computing the value of ๐‘€GTExp that would result in link load closer to ๐ฟ๐‘– using the selected values of ๐ฟ๐‘– of the ๐‘–th source, ๐พ, and Poisson distribution parameters. Considering the same parameter values for burst and gap lengths and ๐‘€BLPos=1, (18) can be simplified as ๐‘€GTExp=๐พโ‹…๐‘๐‘…๐‘€BLPos๎‚ธ1๐ฟ๐‘–๎‚นโˆ’1.(19)

2.3. Generating the GTEXP

Bearing in mind that ๐œ†Exp=ACR, since in ABR service a source sends its data with a rate equal to ACR, defining ๐œ†GTExp as the exponential mean arrival rate for GTExp, and considering (14) in terms of time, we get (20) ๐ธ(๐‘ฅ)GTExp=1๐œ†GTExp=๐‘€GTExp.(20) Now the BLPos/GTEXP traffic generator can compute the GTExp using the relation GTExp=โˆ’๐‘€GTExpโ‹…ln(๐‘ˆ).(21)

2.4. BLPos Generating

BLPos can be modeled as follows: for ABR service in ATM networks the value of ACR should be within the range MCRโ‰คACRโ‰คPCR. The source first starts transmitting a random-sized burst of cell at ICR. It then waits for a random amount of time, which follows exponential distribution. The source will go on repeating this cycle BLPos/GTExp until completing transmission. In every cycle the number of cells is determined by ACR, which is changed due to the feedback used in the switch. The source repeats the calculation of Poisson distribution, (4) with ๐œ‡Pos=ACR, ๐‘˜ times, where ๐‘˜ is changing from 0 to peak cell rate (PCR). The results will be saved in an array with size ๐‘˜. The number of cells inside the burst can be found by generating a random number u [6], applying the cumulative distribution function (cdf) for all array values, so that the cumulative values are smaller than or equal to a random number u. The following relation gives the BLPos that is used by the BLPos/GTExp traffic generator: BLPos=๐‘ฅ=๐‘ข๎“๐‘ฅ=0๎‚ต๐‘’โˆ’๐œ‡Posโ‹…๐œ‡โˆ’xPos๎‚ถ๐‘ฅ!,๐œ‡Pos=ACR,(22) where the number of cells inside the burst should be at least one (๐‘€BLPos=1).

3. Analytical Results

Consider the values of ๐ฟ๐‘–, ๐‘…, ๐พ, ๐‘€BLPos, Poisson mean arrival rate (๐œ‡BLPos), and ๐œ†GTExp for BLPos/GTExp traffic generator as given in Table 1. The value of ๐‘€GTExpwas determined using (19) (Table 1).

The analytical result of BLPos/GTExp traffic generator for 1000 count values of ๐‘ˆ, generated by the uniform distribution, are shown in Figure 1, and the corresponding computed values of mean, variance, maximum, and minimum values of BL and GT, and are given in Table 2.

The variations in BLPos, GTExp as functions of ๐œ‡BLPos and ๐œ†GTExp for BLPos/GTExp traffic generator for 100 count values of U are shown in Figures 2 and 3 respectively. The increment steps for ๐œ†BLPos(1โ€“110) cells/sec and ๐œ‡GTExp(1โ€“110) cells/sec are 10 for each. Referring to Table 2 it is seen that the minimum values of BLPos/GTExp are greater than their corresponding values of ๐‘€BLPos/๐‘€BLExp. Referring to Figure 2 it can be concluded that the Poisson mean arrival parameter ๐œ‡BLEXP should not be a very large value, because BLPos will, consequently, be very large as well, and the source will spend most of its time sending only the burst cells with a smaller number of gap intervals for BLPos/GTExp traffic generator resulting in less-bursty traffic. Referring to Figure 3 it can be concluded that the exponential mean arrival parameter ๐œ†GTExo should be selected between 2 and 30โ€‰cells/sec for simulation of real bursty traffic because it offers higher peak values of GTExp. This is further supported by the observation that for ๐œ†GTExo in the range 30 to 100โ€‰cells/sec, the peak values of GTExp have the least variation indicating smoothest traffic.

4. Simulation Results

The ATM network simulation was carried out under Linux network programming. The parameters specified in Table 1 were used for this simulation. Six sources ๐‘†๐‘–(๐‘–=1,2,โ€ฆ6) sending their data at the rate ACR๐‘–(๐‘–=1,2,โ€ฆ6) between minimum cell rate (MCR) and peak cell rate (PCR) were considered. The performance of the relative rate marking (RRM) switch was evaluated for traffic generator with respect to the allowed cell rate (ACR), switch input rate (SWIR)/switch output rate (SWOR), memory access time (MAT), queue length (๐‘„), and cell transfer delay (CTD). The initial value of ACR for sources ๐‘†๐‘– was taken as PCR/2 whereas the final ACR value was kept between 200 to 700โ€‰cells/sec in incremental steps of 100 for ๐‘–=1,2,โ€ฆ6 and taking buffer size=1000 cells, higher queue threshold (๐‘„๐ป)=200 cells, lower queue threshold (๐‘„๐ฟ) = 100 cells, and assuming that each source has to send a total of 1000 cells. The variations in ACR, SWIR/SWOR, MAT, ๐‘„, and CTD as a function of time are given, respectively, in Figures 4โ€“8.

The simulation results for mean, variance, maximum and minimum values of ACR, SWIR, SWOR, MAT, ๐‘„ and CTD are given in Tables 3, 4, and 5 respectively.

It can be seen from Figure 4 that ACR changes due to the feed back received from the switch. When the switch is heavily loaded (congestion) there is a decrease in the ACR and when it is lightly loaded (no impeding congestion) the ACR is either kept constant or it is increased. The relative changes between the SWIR and SWOR are shown in Figure 5 are responsible for the variations in the MAT (Figure 6) and ๐‘„ (Figure 7). When the SWIR becomes greater than SWOR the buffer starts filling and reaches a specified threshold level. The switch then signals the source to start reducing its data rate. Consequently, source ACR reduces, and its effect appears at the queue, which causes a reduction in the rate of increase in the queue size. For a SWIR smaller than a SWOR, the buffer starts becoming empty and when ๐‘„ reaches its minimum value, the source is signaled to start increasing its data rate. There is a time lag between the switch experiencing a traffic load variation, effect of switch feedback control, and the occurrence of the new load due to the feedback.

Referring to Figure 8, and Table 5 it is noticed that the CTD between any source and its corresponding destination changes from minimum value of 0.2โ€‰sec to a maximum value of 1.04โ€‰sec.

Examining the simulation results for ACR, SWIR, SWOR, MAT,โ€‰๐‘„ and CTD, which indicate the performance of RRM switch under BLPos/GTExp traffic generator, it is seen that the switch offers the best performance.

5. Conclusion

In this paper a mathematical modeling of the BLPos/GTExp source traffic generator has been carried out. The analytical results showed that the minimum values of BLPos/GTExp are greater than their corresponding values of MBLPos/MBLExp and the Poisson mean arrival parameter ๐œ‡BLExp should not be a very large value, because BLPos will, consequently, be very large as well, and the source will spend most of its time sending only the burst cells with a smaller number of gap intervals for BLPos/GTExp traffic generator resulting in less-bursty traffic. It can be concluded also that the exponential mean arrival parameter ๐œ†GTExo should be selected between 2 and 30โ€‰cells/sec for simulation of real bursty traffic because it offers higher peak values of GTExp. This is further supported by the observation that for ๐œ†GTExo in the range 30 to 100โ€‰cells/sec, the peak values of GTExp have the least variation indicating smoothest traffic. The analytical results have been verified through simulation for six sources.