Next Article in Journal
Proximity Environmental Feature Based Tree Health Assessment Scheme Using Internet of Things and Machine Learning Algorithm
Next Article in Special Issue
Bidirectional Piezoelectric Energy Harvester
Previous Article in Journal
Non-Invasive Ambient Intelligence in Real Life: Dealing with Noisy Patterns to Help Older People
Previous Article in Special Issue
Joint Resource Optimization in Simultaneous Wireless Information and Power Transfer (SWIPT) Enabled Multi-Relay Internet of Things (IoT) System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analytical Model of a Wireless Sensor Network (WSN) Node Operation with a Modified Threshold-Type Energy Saving Mechanism

by
Wojciech M. Kempa
Faculty of Applied Mathematics, Silesian University of Technology, 23 Kaszubska Str., 44-100 Gliwice, Poland
Sensors 2019, 19(14), 3114; https://doi.org/10.3390/s19143114
Submission received: 29 April 2019 / Revised: 25 June 2019 / Accepted: 11 July 2019 / Published: 14 July 2019
(This article belongs to the Special Issue Energy Harvesting and Energy-Neutral IoT Devices and Systems)

Abstract

:
In this article, a model of the operation of a wireless sensor network (WSN) node with an energy saving mechanism based on a threshold-controlled multiple vacation policy is considered. When the queue of packets directed to the node becomes empty, a multiple vacation period is started during which the receiving/transmitting of packets is blocked. In such a period, successive vacations of a fixed constant duration are taken until a predetermined number of N packets accumulated in the queue is detected. Then, at the completion epoch of this vacation, the processing restarts normally. The analytic approach is based on the conception of an embedded Markov chain; integral equations and renewal theory are applied to study the queue-size transient behaviour. The representations for the Laplace transforms of the queue-size distribution at an arbitrary fixed time t and on the idle and processing periods are obtained. The compact-form formulae for the distributions of the idle and processing period duration are derived. Numerical examples are attached as well.

1. Introduction

The issue of energy saving in wireless sensor networks (WSNs) is crucial from the point of view of network security, including ensuring the reliability of its operation, but also in securing an adequate level of data transmission quality. One of the important tools supporting these processes are the algorithms used for self-organization of sensor networks. The purpose of this type of algorithm, after the initial organization of the network, is primarily to minimize the risk of data loss during transmission, which is associated with the organization of network node structures, and secondly to maximize the battery life of individual sensors (see e.g., [1,2,3]). Another challenge for ensuring the appropriate level of the quality of service (QoS) is the occurrence of temporal death in some sensor nodes (see e.g., [4,5]). There are two different main approaches to solving the problem of energy saving in WSNs. The first one is based on the selection of the shortest route, so uses various-type shortest path algorithms. However, as it turns out, minimizing the costs of energy consumption can lead to an increase in the ratio of lost packets. Hence, another approach chooses not necessarily the shortest route, but one that gives a relatively small ratio of lost packets. As a consequence, the connection may be established using rather peripherally located network nodes, however, they are those for which the data loss ratio is very small (see e.g., [6,7,8,9,10,11] for proposals and discussion on this topic).
In the analysis of the functioning of packet network nodes and in the process of network optimization, queueing theory is often used. The subjects of its research are stochastic characteristics of models in which the individual elements (jobs, packets, items etc.) of the input stream are being processed, but due to insufficient service speed or limitations in the processing process, there occur phenomena like the accumulation of jobs waiting to start the service and sometimes (in the case of finite accumulating buffer capacity) the loss of incoming jobs.
In [12] (compare also [13,14]) the queueing model of the M / M / 1 -type with a control mechanism based on the N-policy is proposed as a tool for reducing power consumption of nodes in wireless sensor networks. A G e o / G / 1 -type discrete-time model with disasters is considered in [15]. The approach based on Petri nets is proposed in [16] in energy consumption evaluation. A queueing model with priorities is investigated in [17]. Transient (time-dependent, non-stationary) results for queueing systems with multiple (repeated) vacation policy, used as a model of energy saving mechanisms, can be found in [18,19,20,21] (compare also [22,23,24]). In particular, in [18,21] the departure process is investigated in cases of infinite and finite buffer capacities, respectively, while in [19] the queueing delay is investigated (also for the single vacation policy). Similarly, transient analysis of power saving models based on the N-policy can be found in [25] (queueing delay) and [26] (departure process).
In the article we study a model of the operation of a WSN node with an energy saving mechanism based on a threshold-controlled multiple vacation policy. When the queue of incoming packets directed to the node empties, a multiple vacation period is initialized, during which the receiving/transmitting of packets is blocked. In such a period, successive vacations of a fixed constant duration are taken until a predetermined number of N packets accumulated in the queue is detected. Then, at the completion epoch of this vacation, the processing restarts normally.
We propose the analytic approach based on the conception of an embedded Markov chain, integral equations and renewal theory to investigate the queue-size transient distribution. The explicit compact-form representation for the Laplace transform of the queue-size distribution at arbitrary fixed time t and on the idle and processing periods are obtained. The compact-form formulae for the distributions of the idle and processing period duration are derived. A numerical illustration example is attached as well.
The remaining part of the paper is organized as follows. In the next Section 2 we state a detailed mathematical description of the considered queueing system. In Section 3 we analyze the queue-size distribution on the idle period of the system and, moreover, in Section 4 and Section 5 we give explicit results for the probability mass function of the number of packets present in the system at the idle period completion epoch, and the distribution of the idle period duration Laplace transform (LT), respectively. The queue-size distribution on the processing (busy) period is studied in Section 6. The representation for the LT of the distribution of the busy period duration can be found in Section 7. Finally, Section 8 contains the most general result in that the queue-size at an arbitrary epoch is considered.

2. Queueing Model Description

In the article we consider an M / G / 1 -type queueing model with finite buffer capacity and a control (energy saving) mechanism being a mix of the classical multiple vacation policy and the threshold-type policy (N-policy). Successive data packets arrive according to a Poisson process with rate λ and are being processed according to a FIFO (First-In-First-Out) service discipline with generally-distributed processing times with a cumulative distribution function (CDF) F ( · ) and with its Laplace-Stieltjes transform (LST) f ( · ) . The capacity of the accumulating buffer is finite and equals K 1 , so the maximum system size is K .
Every time the queue empties an energy saving mode is initialized, during which the processing is completely blocked (idle period). Namely, successive server vacations of constant length T are started until the number of accumulated (and waiting for processing) packets equals 1 N K or more, where N is fixed. In the occurrence of such a situation, at the moment of this “successful” queue-size measurement the idle period completes and a new processing (busy) period begins immediately.

3. Queue-Size on Processing Suspension Period

Let us investigate, firstly, the queue-size behaviour on the first processing suspension period (idle period) I 1 that starts at time t = 0 and consists of a number of independent constant-length vacations being initialized as far as the buffer queue length at the end of one of them will reach at least the predetermined level 1 N K . Evidently, we have
P { X ( t ) = m , t I 1 } = i = 1 P { X ( t ) = m , t I 1 , i } ,
where I 1 , i stands for the ith successive vacation during the first idle period I 1 .
Let us denote by 1 A ( · ) the characteristic function of the set A, i.e.,
1 A ( x ) = d e f 1 , if x A , 0 , if x A ,
and let δ i , j stand for the Kronecker delta function.
Note that the following representation holds:
P { X ( t ) = m , t I 1 , i } = 1 [ ( i 1 ) T , i T ) ( t ) { 1 { 0 , , N 1 } ( m ) k = 0 m [ λ ( i 1 ) T ] k k ! e λ ( i 1 ) T · [ λ ( t ( i 1 ) T ) ] m k ( m k ) ! e λ [ t ( i 1 ) T ] + 1 { N , , K 1 } ( m ) k = 0 N 1 [ λ ( i 1 ) T ] k k ! e λ ( i 1 ) T [ λ ( t ( i 1 ) T ) ] m k ( m k ) ! e λ [ t ( i 1 ) T ] + δ m , K k = 0 N 1 [ λ ( i 1 ) T ] k k ! e λ ( i 1 ) T j = K k [ λ ( t ( i 1 ) T ) ] j j ! e λ [ t ( i 1 ) T ] } .
Let us briefly comment on (3). Indeed, if t I 1 , i then t [ ( i 1 ) T , i T ) . The first summand on the right side of (3) refers to the situation in which the threshold level N is not reached before time t . In the second summand the number of packets accumulated in the buffer queue reaches at least N at time t but the buffer is still not saturated (clearly, up to time ( i 1 ) T the number of arrivals must equal at most N 1 in order to initialize the ith vacation). Finally, the last summand on the right side of (3) describes the situation, in that the buffer is saturated at time t .
Before we will state the representation for the LT of P { X ( t ) = m , t I 1 , i } , let us observe that for s > 0 we have
( i 1 ) T i T e ( s + λ ) t [ λ ( t ( i 1 ) T ) ] j j ! d t = 0 T e ( s + λ ) [ u + ( i 1 ) T ] ( λ u ) j j ! d u = e ( s + λ ) ( i 1 ) T λ j ( s + λ ) j + 1 E s + λ ( j , T ) ,
where E s + λ ( j , · ) denotes the CDF of the j-Erlang distribution with parameter s + λ .
Referring to (2)–(4), we get now for s > 0
p I ( s , m ) = d e f 0 e s t P { X ( t ) = m , t I 1 } d t = i = 1 e ( s + λ ) ( i 1 ) T { 1 { 0 , , N 1 } ( m ) · k = 0 m λ m k ( s + λ ) m k + 1 · [ λ ( i 1 ) T ] k k ! E s + λ ( m k , T ) + k = 0 N 1 [ λ ( i 1 ) T ] k k ! · 1 { N , , K 1 } ( m ) λ m k ( s + λ ) m k + 1 E s + λ ( m k , T ) + δ m , K j = K k λ j ( s + λ ) j + 1 E s + λ ( j , T ) } .

4. Queue Measurement at the Idle Period Completion Epoch

In this section we are interested in the probability distribution of the queue state X I r = X I at the completion epoch of the arbitrary processing suspension period (idle period) I r .
Denoting
a n = d e f P { X I = n } , n { N , , K } ,
observe that the following equality is true:
a n = i = 1 P { X ( i T ) = n , idle period contains exactly i single vacations } = i = 1 ( 1 { N , , K 1 } ( n ) k = 0 N 1 [ λ ( i 1 ) T ] k k ! e λ ( i 1 ) T · ( λ T ) n k ( n k ) ! e λ T + δ n , K k = 0 N 1 [ λ ( i 1 ) T ] k k ! e λ ( i 1 ) T j = K k ( λ T ) j j ! e λ T ) .
Indeed, the first summand on the right side of (7) refers to the situation where the buffer is not fulfilled before the endpoint of the idle period, while for the second one the processing starts after the idle period with exactly K packets present.

5. Idle Period Duration

In this section, without loss of generality, we identify the rth idle period I r with its duration. Obviously, we have
P { I r i T } = k = 0 N 1 [ λ ( i 1 ) T ] k k ! e λ ( i 1 ) T , r = 1 , 2 , ,
due to the fact that the ith vacation during the idle period is being initialized if and only if the number of accumulated packets at time ( i 1 ) T does not exceed N 1 .
From (8) we obtain immediately
P { I r = i T } = P { I r i T } P { I r ( i + 1 ) T } = k = 0 N 1 [ λ ( i 1 ) T ] k k ! e λ ( i 1 ) T ( λ i T ) k k ! e λ i T .
Introducing the LST f I ( · ) of the CDF of the arbitrary idle period duration, we get
f I ( s ) = d e f i = 0 e s i T P { I r = ( i + 1 ) T } P { I r = i T } = i = 0 e s i T · k = 0 N 1 2 ( λ i T ) k k ! e λ i T [ λ ( i + 1 ) T ] k k ! e λ ( i + 1 ) T [ λ ( i 1 ) T ] k k ! e λ ( i 1 ) T .

6. Queue-Size on Processing Period

We investigate now the case of the queue-size distribution on a single processing period. Obviously, the transient behavior of this stochastic characteristic depends on the system state n at the completion epoch of the preceding idle period, where n { N , , K } . Without loss of generality, let us assume that the processing period starts at time t = 0 and, temporarily, that at the starting moment the number of packets accumulated in the buffer can take on any value from { 1 , , K } . To receive a representation for the LT of the conditional queue-size distribution on a processing period, we will use the method of potential proposed by Korolyuk in [27]. Let us denote by X ( t ) the number of packets present in the system exactly at time t .
Introduce the following notations:
P n ( P ) ( t , m ) = d e f P { X ( t ) = m , t is in a processing period | X ( 0 ) = n } ,
where m , n { 1 , , K } and t > 0 , and
p n ( P ) ( s , m ) = d e f 0 e s t P n ( P ) ( t , m ) d t , s > 0 .
It is well known that, due to the memoryless property of the exponential distribution of interarrival times, the sequence X ( t i ) , i = 1 , 2 , , where t i stands for the ith processing completion epoch after the starting moment, is an embedded Markov chain for the process { X ( t ) , t 0 } . Thus, utilizing the formula of total probability with respect to the first departure moment after t = 0 , we obtain for n = 1
P 1 ( P ) ( t , m ) = i = 1 K 2 0 t ( λ x ) i i ! e λ x P i ( P ) ( t x , m ) d F ( x ) + i = K 1 0 t ( λ x ) i i ! e λ x P K 1 ( P ) ( t x , m ) d F ( x ) + [ 1 F ( t ) ] e λ t 1 { 1 , , K 1 } ( m ) ( λ t ) m 1 ( m 1 ) ! + δ m , K i = K 1 ( λ t ) i i ! .
Similarly, for 2 n K , we have
P n ( P ) ( t , m ) = i = 0 K n 1 0 t ( λ x ) i i ! e λ x P n + i 1 ( P ) ( t x , m ) d F ( x ) + i = K n 0 t ( λ x ) i i ! e λ x P K 1 ( P ) ( t x , m ) d F ( x ) + [ 1 F ( t ) ] e λ t 1 { n , , K 1 } ( m ) ( λ t ) m n ( m n ) ! + δ m , K i = K n ( λ t ) i i ! .
Observe that the first summand on the right side of (13) and (14) describes the case in that the accumulating buffer does not become fulfilled before the first departure occurring precisely at time x < t , while the second one refers to the situation in which the system is saturated before x . The case of the first processing finishing after t is presented in the last summand on the right side of (13) and (14). Let us note that the first sum in (14) is taken from 0 while in (13) from 1. Indeed, if the system contains exactly one packet at the starting moment then no departure occurs before t , otherwise the processing (busy) period will end before t .
Let us define for s > 0 the following functional sequences:
γ n ( s ) = d e f 0 e ( λ + s ) t ( λ t ) n n ! d F ( t )
and
η n ( s , m ) = d e f 0 e ( λ + s ) t [ 1 F ( t ) ] 1 { n , , K 1 } ( m ) ( λ t ) m n ( m n ) ! + δ m , K i = K n ( λ t ) i i ! d t .
Implementing (15) and (16) into (13) and (14) gives
p 1 ( P ) ( s , m ) = i = 1 K 2 γ i ( s ) p i ( P ) ( s , m ) + p K 1 ( P ) ( s , m ) i = K 1 γ i ( s ) + η 1 ( s , m )
and
p n ( P ) ( s , m ) = i = 0 K n 1 γ i ( s ) p n + i 1 ( P ) ( s , m ) + p K 1 ( P ) ( s , m ) i = K n γ i ( s ) + η n ( s , m ) ,
where 2 n K .
Let us rewrite equations of the system (17) and (18) in a another form. Indeed, substituting p ˜ n ( P ) ( s , m ) = d e f p K n ( P ) ( s , m ) , where 0 n K 1 , we transform them as follows:
p ˜ K 1 ( P ) ( s , m ) = i = 1 K 2 γ i ( s ) p ˜ K i ( P ) ( s , m ) + p ˜ 1 ( P ) ( s , m ) i = K 1 γ i ( s ) + η 1 ( s , m ) ,
k = 1 n γ k + 1 ( s ) p ˜ n k ( P ) ( s , m ) p ˜ n ( P ) ( s , m ) = ω n ( s , m ) ,
where 0 n K 2 , and
ω n ( s , m ) = d e f γ n + 1 ( s ) p ˜ 0 ( P ) ( s , m ) p ˜ 1 ( P ) ( s , m ) i = n + 1 γ i ( s ) η K n ( s , m ) .
Each solution of the linear infinite-size system with known coefficients a k ’s and b k ’s
k = 1 n a k + 1 x n k x n = b n , n 0 ,
can be written in the following form (see [27,28]):
x n = c r n + 1 + i = 0 n r n i b i ,
where c is certain constant and the sequence ( r n ) is defined by the following recursive relationship:
r n = 0 , for n = 0 , 1 / a 0 , for n = 1 , ( 1 / a 0 ) ( r n 1 i = 0 n 1 a i + 1 r n 1 i ) , for n 2 .
At this stage let us make some essential observations. Firstly, the systems (20) and (22) are very similar, however in the former, functional sequences occur. Secondly, the number of equations in (20) is finite, while in (22) it is infinite. As it turns out, the representation (23) can be used to find the solution of the system (19) and (20) and, moreover, the fact that the number of equations in (20) is finite allows for finding an explicit representation for c .
So, we can write (compare (23) and (24))
p ˜ n ( P ) ( s , m ) = c ( s , m ) r n + 1 ( s ) + i = 0 n r n i ( s ) ω i ( s , m ) ,
where c ( s , m ) does not depend on the index n and the functional sequence ( r n ( s ) ) is defined via ( γ n ( s ) ) in the following way:
r n ( s ) = 0 , for n = 0 , 1 / γ 0 ( s ) , for n = 1 , ( 1 / γ 0 ( s ) ) [ r n 1 i = 0 n 1 γ i + 1 ( s ) r n 1 i ( s ) ] , for n 2 .
Note that to represent p ˜ n ( P ) ( s , m ) explicitly we need to find the formulae for p ˜ 0 ( P ) ( s , m ) and p ˜ 1 ( P ) ( s , m ) occurring in the definition (21) of ω n ( s , m ) .
Let us start with substituting n = 0 into (25). We obtain
c ( s , m ) = γ 0 ( s ) p ˜ 0 ( P ) ( s , m ) .
Writing (20) for n = 0 leads to
p ˜ 1 ( P ) ( s , m ) = [ γ 0 ( s ) ] 1 [ ω 0 ( s , m ) + p ˜ 0 ( P ) ( s , m ) 1 γ 1 ( s ) ]
and hence, referring to (21), we get
p ˜ 1 ( P ) ( s , m ) = p ˜ 0 ( P ) ( s , m ) η K ( s , m ) f ( s ) .
As a consequence, implementing representations (27) and (29) in the formula (25), we have
p ˜ n ( P ) ( s , m ) = g n ( s ) p ˜ 0 ( P ) ( s , m ) + h n ( s , m ) ,
where n 0 and
g n ( s ) = d e f γ 0 ( s ) r n + 1 ( s ) + i = 0 n r n i ( s ) { γ i + 1 ( s ) [ f ( s ) ] 1 j = i + 1 γ j ( s ) }
and
h n ( s , m ) = d e f i = 0 n r n i ( s ) { [ f ( s ) ] 1 η K ( s , m ) j = i + 1 γ j ( s ) η K i ( s , m ) } .
Applying now the representation (30) in the Equation (19), we can easily eliminate p ˜ 0 ( P ) ( s , m ) in the form
p ˜ 0 ( P ) ( s , m ) = U ( s , m ) V ( s ) ,
where
U ( s , m ) = d e f i = 1 K 2 γ i ( s ) h K i ( s , m ) [ f ( s ) ] 1 η K ( s , m ) i = K 1 γ i ( s ) + η 1 ( s , m ) h K 1 ( s , m ) ,
V ( s ) = d e f h K 1 ( s ) i = 1 K 2 γ i ( s ) g K i ( s ) [ f ( s ) ] 1 i = K 1 γ i ( s ) 1 .
Finally, returning to p n ( P ) ( s , m ) (instead of p ˜ n ( P ) ( s , m ) ) and referring to (21), (25), (27), (29) and (33), we can formulate the following result:
p n ( P ) ( s , m ) = 0 e s t P n ( P ) ( t , m ) d t = γ 0 ( s ) r K n + 1 ( s ) + i = 0 K n r K n i ( s ) γ i + 1 ( s ) ( f ( s ) ) 1 j = i + 1 γ j ( s ) · U ( s , m ) V ( s ) + i = 0 K n r K n i ( s ) ( f ( s ) ) 1 η K ( s , m ) j = i + 1 γ j ( s ) η K i ( s , m ) ,
where s > 0 and the formulae for γ k ( s ) , η k ( s , m ) , r k ( s ) , U ( s , m ) and V ( s ) are given in (15), (16), (26), (34) and (35), respectively.

7. Processing Period Duration

In this section we establish the formula for the LST of the first busy (processing) period duration, conditioned by the number of packets accumulated in the buffer at the starting epoch. The appropriate result we will write in the most general case, namely assuming that the processing period may start with n { 1 , , K } packets present (not necessarily with at least N packets present). Suppose temporarily that the first busy period starts at t = 0 .
Introduce the following notations:
F B P , n ( x ) = d e f P { first busy period duration < x | busy period starts with n packets present }
and
f B P , n ( s ) = d e f 0 e s x d F B P , n ( x ) , s > 0 .
It is easy to note that for f B P , n ( s ) , n { 1 , , K } , the following system of equations is true (compare (17) and (18)):
f B P , 1 ( s ) = i = 1 K 2 γ i ( s ) f B P , i ( s ) + f B P , K 1 ( s ) i = K 1 γ i ( s ) + f ( λ + s ) ,
f B P , n ( s ) = i = 0 K n 1 γ i ( s ) f B P , n + i 1 ( s ) + f B P , K 1 ( s ) i = K n γ i ( s ) ,
where 2 n K .
Let us observe that if we define
η ˜ n ( s , m ) = d e f δ n , 1 f ( λ + s ) ,
then we can obtain the explicit-form representation for f B P , n ( s ) utilizing the Formula (36). Indeed, we obtain
f B P , n ( s , m ) = γ 0 ( s ) r K n + 1 ( s ) + i = 0 K n r K n i ( s ) γ i + 1 ( s ) ( f ( s ) ) 1 j = i + 1 γ j ( s ) · U ˜ ( s , m ) V ( s ) + i = 0 K n r K n i ( s ) ( f ( s ) ) 1 η ˜ K ( s , m ) j = i + 1 γ j ( s ) η ˜ K i ( s , m ) ,
where s > 0 and now
U ˜ ( s , m ) = d e f i = 1 K 2 γ i ( s ) h ˜ K i ( s , m ) [ f ( s ) ] 1 η ˜ K ( s , m ) i = K 1 γ i ( s ) + η ˜ 1 ( s , m ) h ˜ K 1 ( s , m )
and
h ˜ n ( s , m ) = d e f i = 0 n r n i ( s ) [ f ( s ) ] 1 η ˜ K ( s , m ) j = i + 1 γ j ( s ) η ˜ K i ( s , m ) .

8. Queue-Size Distribution at an Arbitrary Time

Consider now the transient queue-size distribution at arbitrary time t . Assume that the system starts the operation at time t = 0 with the first idle period. The evolution of the considered queueing system can be observed on successive idle periods I 1 , I 2 , interspersed with busy periods B P 1 , B P 2 , . The formula of total probability gives
P { X ( t ) = m } = i = 1 P { X ( t ) = m , t I i } + P { X ( t ) = m , t B P i } .
Due to the memoryless property of the exponential distribution of interarrival times, the following representation is true:
P { X ( t ) = m , t I i } = l 1 , , l i 1 { N , , K } j = 1 i 1 a l j 0 t P { X ( t y ) = m , t y I 1 } · d F I ( i 1 ) i * * F B P , a l 1 * * F B P , a l i 1 ( y ) ,
where F I ( · ) stands for the CDF of the idle period duration (all idle periods are independent and identically distributed) and * denotes the Stieltjes convolution. The sum on the right side of (46) is taken over all possible values of l 1 , , l i 1 { N , , K } .
Quite similarly we obtain
P { X ( t ) = m , t B P i } = l 1 , , l i { N , , K } j = 1 i 1 a l j 0 t P { X ( t y ) = m , t y B P 1 | B P 1 starts with l i packets present } · d F I i * * F B P , a l 1 * * F B P , a l i 1 ( y ) .
Taking LTs of both sides in (46) and (47) and referring to (45), we obtain the following representation for the LT of the queue-size distribution in the considered queueing model with a threshold-controlled multiple vacation policy:
0 e s t P { X ( t ) = m } d t = i = 1 { p I ( s , m ) [ f I ( s ) ] i 1 l 1 , , l i 1 { N , , K } j = 1 i 1 a l j f B P , a l j ( s ) + l 1 , , l i { N , , K } a l i p a l i ( P ) ( s , m ) j = 1 i 1 a l j f B P , a l j ( s ) } ,
where s > 0 .

9. Numerical Study

In this section we present some numerical examples illustrating theoretical results obtained in previous sections. Let us investigate the operation of the system on the first idle period. Assume that the intensity of Poisson arrivals equals 1.1 items per time unit (mean interarrival times are equal to 0.909) and that the processing times are distributed according to a mixture of two exponential distributions, namely the probability density function f S ( · ) of the service time has the following form:
f S ( t ) = 0.25 · 2 exp ( 2 t ) + 0.75 · exp ( t ) , t > 0 .
Hence the mean service time of a single packet equals 0.875 and, in consequence, the offered load for the system is ϱ = 0.963 . Let us fix K = 8 , N = 3 and T = 0.8 . In Figure 1 the queue-size distributions P { X ( t ) = m , t I 1 } for m = 0 (solid line), m = 1 (dashed line) and 2 (dotted line) are visualized. As one can observe, for m > 0 the graphs have characteristic shapes in which the character of the function changes every 0.8 time unit. Obviously, it is connected with the fixed duration of a single server vacation equal to 0.8.
In Figure 2 and Figure 3 the appropriate results are presented for two other arrival intensities. In Figure 2 the case of a relatively low arrival intensity λ = 0.7 (mean interarrival time equal to 1.429 and the offered load ϱ = 0.612 ) is shown. Similarly, the case of the overloaded system in that λ = 1.4 (mean interarrival time equal to 0.714 and the offered load ϱ = 1.225 ) is visualized in Figure 3.
We investigate now the behaviour of the probability distribution of the idle period duration. Evidently, this time is a certain multiple of the single vacation duration T . In Figure 4, Figure 5 and Figure 6, respectively, the cases of λ = 1.1 , 0.7 and 1.4 are presented. As it can be observed, the greater arrival intensity, the greater probability that the idle period will finish after one server vacation only (at time T = 0.8 ). Moreover, let us note that (except for the case of one single vacation case) the maximum of the distribution moves to the right with the decrease of λ (or with the decrease of the offered load).
Finally, in Figure 7, Figure 8 and Figure 9 the probability distributions of the number of packets present in the system at the completion epoch of the idle period are shown for the same three values of the arrival intensity, so for the offered traffic load ϱ = 0.963 , 0.612 and 1.225 .

10. Conclusions

In the paper a queueing system with a threshold-type multiple vacation policy is considered as a model of the operation of a WSN node with an energy saving mechanism. The idea of the algorithm is that the node radio transmitter/receiver goes to sleep mode every time the queue of packets directed to the node becomes empty. Being in this mode it takes successive independent vacations of deterministic duration until the predetermined number of packets present is detected. The time-dependent queue-size distribution is studied, applying the analytic approach based on the idea of the embedded Markov chain, integral equations and renewal theory. The compact-form formulae for the Laplace transforms of the queue-size distribution at an arbitrary fixed time t, and on the idle and processing periods are obtained. Moreover, the compact-form formulae for the distributions of the idle and processing period duration are derived. Numerical examples illustrating the theoretical results are obtained.

Funding

This research was funded by Silesian University of Technology, Gliwice, Poland.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Dressler, F. Self-Organization in Sensor and Actor Networks; John Wiley & Sons Ltd.: Chichester, UK, 2007. [Google Scholar]
  2. Pathan, A.S.K. Security of Self-Organizing Networks: MANET, WSN, WMN, VANET; CRC Press: Boca Raton, FL, USA, 2011. [Google Scholar]
  3. Ruiz, L.B.; Braga, T.R.M.; Silva, F.A.; Assunção, H.P.; Nogueira, J.M.S.; Loureiro, A.A. On the design of a self-managed wireless sensor network. IEEE Commun. Mag. 2005, 43, 95–102. [Google Scholar] [CrossRef]
  4. Tan, L.; Tang, S. Energy harvesting wireless sensor node with temporal death: Novel models and analyses. IEEE/ACM Trans. Netw. 2017, 25, 896–909. [Google Scholar] [CrossRef]
  5. Tang, S.; Tan, L. Reward rate maximization and optimal transmission policy of EH device with temporal death in EH-WSNs. IEEE Trans. Wirel. Commun. 2017, 16, 1157–1167. [Google Scholar] [CrossRef]
  6. Hosen, A.S.M.S.; Cho, G.H. An energy centric cluster-based routing protocol for wireless sensor networks. Sensors 2018, 18, 1520. [Google Scholar] [CrossRef] [PubMed]
  7. Huang, Y.; Tian, Y.; Cheng, W. Optimization of energy saving for wireless sensor networks. Environ. Eng. Manag. J. 2014, 13, 1057–1070. [Google Scholar] [CrossRef]
  8. Luo, J.; Hu, J.; Wu, D.; Li, R. Opportunistic routing algorithm for relay node selection in wireless sensor networks. IEEE Trans. Ind. Inform. 2015, 11, 112–121. [Google Scholar] [CrossRef]
  9. Mostafaei, H. Energy-efficient algorithm for reliable routing of wireless sensor networks. IEEE Trans. Ind. Electron. 2019, 66, 5567–5575. [Google Scholar] [CrossRef]
  10. Srbinovska, M.; Cundeva-Blajer, M. Optimization methods for energy consumption estimation in wireless sensor networks. J. Sustain. Dev. Energy Water Environ. Syst. 2019, 7, 261–274. [Google Scholar] [CrossRef]
  11. Zhou, B.; Ahn, D.; Lee, J.; Sun, C.; Ahmed, S.; Kim, Y. A passive tracking system based on geometric constraints in adaptive wireless sensor networks. Sensors 2018, 18, 3276. [Google Scholar] [CrossRef] [PubMed]
  12. Jiang, F.C.; Huang, D.C.; Yang, C.T.; Leu, F.Y. Lifetime elongation for wireless sensor network using queue-based approaches. J. Supercomput. 2012, 59, 1312–1335. [Google Scholar] [CrossRef]
  13. Ghosh, S.; Unnikrishnan, S. Reduced power consumption in wireless sensor networks using queue based approach. In Proceedings of the 5th IEEE International Conference on Advances in Computing, Communication and Control (ICAC3), Mumbai, India, 1–2 December 2017. [Google Scholar] [CrossRef]
  14. Lenin, R.B.; Ramaswamy, S. Performance analysis of wireless sensor networks using queuing networks. Ann. Oper. Res. 2015, 233, 237–261. [Google Scholar] [CrossRef]
  15. Lee, D.H.; Yang, W.S. The N-policy of a discrete time Geo/G/1 queue with disasters and its application to wireless sensor networks. Appl. Math. Model. 2013, 37, 9722–9731. [Google Scholar] [CrossRef]
  16. Li, J.; Zhou, H.Y.; Zuo, D.C.; Hou, K.M.; Xie, H.P.; Zhou, P. Energy consumption evaluation for wireless sensor network nodes based on queuing Petri net. Int. J. Distrib. Sens. Netw. 2014, 10, 262848. [Google Scholar] [CrossRef]
  17. Ke, J.-F.; Chen, W.-J.; Huang, D.-C. A life extend approach based on priority queue N strategy for wireless sensor network. In Proceedings of the 11th EAI International Conference on Heterogeneous Networking for Quality, Reliability, Security and Robustness (QSHINE), Taipei, Taiwan, 19–20 August 2015; pp. 272–279. [Google Scholar]
  18. Kempa, W.M. Analysis of departure process in batch arrival queue with multiple vacations and exhaustive service. Commun. Stat. Theor. Methods 2011, 40, 2856–2865. [Google Scholar] [CrossRef]
  19. Kempa, W.M. Transient workload distribution in the M/G/1 finite-buffer queue with single and multiple vacations. Ann. Oper. Res. 2016, 239, 381–400. [Google Scholar] [CrossRef]
  20. Kempa, W.M. On buffer overflow duration in a finite-capacity queueing system with multiple vacation policy. In Proceedings of the 43rd International Conference on Applications of Mathematics in Engineering and Economics (AMEE), Sozopol, Bulgaria, 7 December 2017. [Google Scholar]
  21. Kempa, W.; Marjasz, R. Departure counting process in a wireless network node with sleep mode modelled via repeated vacations. In Proceedings of the 23rd International Conference on Information and Software Technologies (ICIST), Druskininkai, Lithuania, 12–14 October 2017; pp. 395–407. [Google Scholar]
  22. Kempa, W.M. Buffer overflow duration in a model of WSN node with power saving mechanism based on SV policy. In Proceedings of the 23rd International Conference on Information and Software Technologies (ICIST), Druskininkai, Lithuania, 12–14 October 2017; pp. 385–394. [Google Scholar]
  23. Kempa, W.M. Transient queue-size distribution in a finite-capacity queueing system with server breakdowns and Bernoulli feedback. In Proceedings of the 43rd International Conference on Applications of Mathematics in Engineering and Economics (AMEE), Sozopol, Bulgaria, 7 December 2017. [Google Scholar]
  24. Kempa, W.M. On buffer overflow duration in WSN with a vacation-type power saving mechanism. In Proceedings of the 24th International Conference on Systems, Signals and Image Processing (IWSSIP), Poznań, Poland, 22–24 May 2017. [Google Scholar]
  25. Kempa, W.M. On queueing delay in WSN with energy saving mechanism based on queued wake up. In Proceedings of the 21st International Conference on Systems Signals and Image Processing, Dubrovnik, Croatia, 12–15 May 2014; pp. 187–190. [Google Scholar]
  26. Kempa, W.M. Time-dependent analysis of transmission process in a wireless sensor network with energy saving mechanism based on threshold waking up. In Proceedings of the IEEE 16th International Workshop on Signal Processing Advances in Wireless Communications, Stockholm, Sweden, 28 June–1 July 2015; pp. 26–30. [Google Scholar]
  27. Korolyuk, V.S. Boundary-Value Problems for Compound Poisson Processes; Naukowa Dumka: Kiev, Ukraine, 1975. [Google Scholar]
  28. Kempa, W.M. A comprehensive study on the queue-size distribution in a finite-buffer system with a general independent input flow. Perform. Eval. 2017, 108, 1–15. [Google Scholar] [CrossRef]
Figure 1. Transient probabilities P { X ( t ) = m , t I 1 } for λ = 1.1 and for m = 0 , 1 , 2 .
Figure 1. Transient probabilities P { X ( t ) = m , t I 1 } for λ = 1.1 and for m = 0 , 1 , 2 .
Sensors 19 03114 g001
Figure 2. Transient probabilities P { X ( t ) = m , t I 1 } for λ = 0.7 and for m = 0 , 1 , 2 .
Figure 2. Transient probabilities P { X ( t ) = m , t I 1 } for λ = 0.7 and for m = 0 , 1 , 2 .
Sensors 19 03114 g002
Figure 3. Transient probabilities P { X ( t ) = m , t I 1 } for λ = 1.4 and for m = 0 , 1 , 2 .
Figure 3. Transient probabilities P { X ( t ) = m , t I 1 } for λ = 1.4 and for m = 0 , 1 , 2 .
Sensors 19 03114 g003
Figure 4. Distribution of idle period duration for λ = 1.1 .
Figure 4. Distribution of idle period duration for λ = 1.1 .
Sensors 19 03114 g004
Figure 5. Distribution of idle period duration for λ = 0.7 .
Figure 5. Distribution of idle period duration for λ = 0.7 .
Sensors 19 03114 g005
Figure 6. Distribution of idle period duration for λ = 1.4 .
Figure 6. Distribution of idle period duration for λ = 1.4 .
Sensors 19 03114 g006
Figure 7. Number of packets of at the completion epoch of idle period for λ = 1.1 .
Figure 7. Number of packets of at the completion epoch of idle period for λ = 1.1 .
Sensors 19 03114 g007
Figure 8. Number of packets of at the completion epoch of idle period for λ = 0.7 .
Figure 8. Number of packets of at the completion epoch of idle period for λ = 0.7 .
Sensors 19 03114 g008
Figure 9. Number of packets of at the completion epoch of idle period for λ = 1.4 .
Figure 9. Number of packets of at the completion epoch of idle period for λ = 1.4 .
Sensors 19 03114 g009

Share and Cite

MDPI and ACS Style

Kempa, W.M. Analytical Model of a Wireless Sensor Network (WSN) Node Operation with a Modified Threshold-Type Energy Saving Mechanism. Sensors 2019, 19, 3114. https://doi.org/10.3390/s19143114

AMA Style

Kempa WM. Analytical Model of a Wireless Sensor Network (WSN) Node Operation with a Modified Threshold-Type Energy Saving Mechanism. Sensors. 2019; 19(14):3114. https://doi.org/10.3390/s19143114

Chicago/Turabian Style

Kempa, Wojciech M. 2019. "Analytical Model of a Wireless Sensor Network (WSN) Node Operation with a Modified Threshold-Type Energy Saving Mechanism" Sensors 19, no. 14: 3114. https://doi.org/10.3390/s19143114

APA Style

Kempa, W. M. (2019). Analytical Model of a Wireless Sensor Network (WSN) Node Operation with a Modified Threshold-Type Energy Saving Mechanism. Sensors, 19(14), 3114. https://doi.org/10.3390/s19143114

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop