Next Article in Journal
A Lane Detection Method Based on a Ridge Detector and Regional G-RANSAC
Previous Article in Journal
An Adaptive Real-Time Detection Algorithm for Dim and Small Photoelectric GSO Debris
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Message Bundling with Delay and Synchronization Constraints in Wireless Sensor Networks

1
Department of Electrical and Electronic Engineering, Xi’an Jiaotong-Liverpool University (XJTLU), Suzhou 215123, China
2
Department of Architecture, Xi’an Jiaotong-Liverpool University (XJTLU), Suzhou 215123, China
*
Author to whom correspondence should be addressed.
Submission received: 1 August 2019 / Revised: 8 September 2019 / Accepted: 16 September 2019 / Published: 18 September 2019
(This article belongs to the Section Sensor Networks)

Abstract

:
Energy efficiency and end-to-end delay are two of the major requirements for the monitoring and detection applications based on resource-constrained wireless sensor networks (WSNs). As new advanced technologies for accurate monitoring and detection—such as device-free wireless sensing schemes for human activity and gesture recognition—have been developed, time synchronization accuracy becomes an important requirement for those WSN applications too. Message bundling is considered one of the effective methods to reduce the energy consumption for message transmissions in WSNs, but bundling more messages increases the transmission interval of bundled messages and thereby their end-to-end delays; the end-to-end delays need to be maintained within a certain value for time-sensitive applications like factory monitoring and disaster prevention, while the message transmission interval affects time synchronization accuracy when the bundling includes synchronization messages as well. Taking as an example a novel WSN time synchronization scheme recently proposed for energy efficiency, we investigate an optimal approach for message bundling to reduce the number of message transmissions while maintaining the user-defined requirements on end-to-end delay and time synchronization accuracy. Formulating the optimal message bundling problem as integer linear programming, we compute a set of optimal bundling numbers for the sensor nodes to constrain their link-level delays, thereby achieving and maintaining the required end-to-end delay and synchronization accuracy. Extensive experimental results based on a real WSN testbed using TelosB sensor nodes demonstrate that the proposed optimal bundling could reduce the number of message transmissions about 70% while simultaneously maintaining the required end-to-end delay and time synchronization accuracy.

1. Introduction

In a typical wireless sensor network (WSN) based on lots of battery-powered sensor nodes, minimizing the energy consumptions at the sensor nodes is critical to the service time of the WSN. Considering that the radio activities consume the most energy at sensor nodes, reducing the number of message transmissions is key to saving the energy [1,2], and the data bundling (The terms of “data bundling” and “message bundling” are used interchangeably in this paper.). is an efficient way to achieve it [3,4,5]. Depending on the operation locations, the existing data bundling schemes (e.g., [4,5,6,7,8,9,10,11,12]) can be classified into in-node, in-network, and hybrid bundling. However, one major disadvantage of employing data bundling schemes is the increase of end-to-end (E2E) delay, which draws much attention in WSN research [13,14,15,16,17]. By E2E delay, we mean the difference between the measurement time at an originating sensor node and the reception time of the resulting measurement data by the head node via a message. The data bundling has been studied so far in the context of energy efficiency and E2E delay [4,5,11,12], but synchronization accuracy is hardly considered as one of the requirements because conventional energy-efficient time synchronization schemes are mostly designed independent of data bundling.
Note that, for most long-term monitoring and detection applications in WSNs, energy-efficient and high-precision time synchronization is crucial to ordering and analyzing the data measured over a long period of time. For instance, in human activity and gesture recognition based on the novel device-free wireless sensing [18,19], the sensed data must be ordered accurately and processed timely, which requires high-precision time synchronization for the former while demands certain latency—i.e., E2E delay—for the latter. Among many research works on WSN time synchronization aiming at improving time synchronization accuracy while lowering energy consumption [1,20,21,22,23,24,25,26], the energy-efficient time synchronization scheme we proposed in [1], which we call it EE-ASCFR—i.e., short for Energy-Efficient time synchronization scheme based on Asynchronous Source Clock Frequency Recovery—throughout the paper—leverages data bundling to further reduce the energy consumption in terms of the number of synchronization message transmissions; thanks to the reverse two-way message exchange scheme adopted in EE-ASCFR, the time synchronization is done at the head node and the corresponding synchronization data are bundled together with measurements into a bundled message at a sensor node. Consequently, the time synchronization accuracy is tied to the interval of the bundled message transmissions from the sensor node to the head node, i.e., the synchronization interval (SI). This means that the number of measurement data in a bundled message is linked to the interval of the bundled message transmissions and thereby not only affects the aforementioned E2E delay, but also the time synchronization accuracy. In such a case, we need to optimize the number of bundled messages under the constraint of E2E delay and time synchronization accuracy in achieving higher energy efficiency.
In this paper, we formulate the aforementioned bundling optimization problem as integer linear programming (ILP), where the number of bundled messages for each sensor node is optimized while jointly satisfying the user-defined performance requirements on E2E delay and time synchronization accuracy. In this way, we can further reduce the energy consumption through data bundling while meeting the requirements.
The rest of this paper is organized as follows: Section 2 provides the preliminaries of this study where we introduce the time synchronization scheme based on the data bundling procedure and discuss the related conflicts of performance metrics. Section 3 presents our proposed approach for optimal message bundling with delay and synchronization constraints and its formulation as ILP. Section 4 exhibits our system design at both head (The head node and the monitoring station—e.g., PC or server connected to it—are jointly called head throughout this paper.) and sensor node. Section 5 demonstrates the performance of the proposed approach through experimental results on a real WSN testbed. Section 6 concludes our work in this paper and discusses future works.

2. Preliminaries

Most WSNs have the fundamental characteristics of limited resources, multi-hop communication, large scale and dynamic environments [27]. To achieve satisfactory performance in typical WSN applications including environmental monitoring, event detection, and industrial applications, several specific requirements—e.g., E2E delay and synchronization accuracy—have to be met. With consideration of the most basic requirements of energy efficiency, jointly maintaining the three metrics is crucial; however, it has not been addressed in the existing research.

2.1. Existing Methods toward Improving Energy Efficiency

In the literature, diverse methods across different layers have been proposed to improve the energy efficiency: transmission power control and multi-channel communication in the physical layer [28,29,30,31,32,33]; maximum (re-)transmission times and retry delay tuning in the media access control (MAC) layer [17,34,35]; packet payload size and sampling rate adaptation in the application layer [36,37,38], as well as cross-layer methods such as [39].
The data bundling method employed in this paper belongs to the category of packet payload size adaptation since tuning the bundling number directly affects the packet payload size. Note that data bundling is a mature and effective method that has been leveraged in both the conventional optical burst switching [40] and the novel mobile wireless network [41], the latter of which adopts it to minimize the memory and communication usage, thereby the energy efficiency is improved. Note that one notable difference of the data bundling method in this paper compared to those in other areas is that it is investigated as an optimal bundling problem under the constraint of time synchronization accuracy in addition to the E2E delay one.

2.2. Energy-Efficient Time Synchronization Schemes Using Data Bundling

To fulfill the desired requirements of energy efficiency and synchronization accuracy, many schemes have been proposed such as [1,20,21,22,23,24,25,26]. Among those, EE-ASCFR proposed in [1] particularly suits the E2E delay calculation since computing the E2E delay at the head is completely in conformity with the preferential asymmetric scenario of EE-ASCFR. Two approaches of reverse two-way message exchange and data bundling are leveraged to reduce the message transmissions as demonstrated in Figure 1; the latter data bundling is the major contributor for conserving energy. We will discuss this in detail later in this section. Nevertheless, the time synchronization is accomplished through the former reverse two-way message exchange, in which the first-order affine clock model is used to model the hardware clock T i of a sensor node i with respect to the reference clock t of the head node. Since the hardware clocks of the sensor nodes have different clock frequency ratios and clock offsets with respect to the reference clock due to their coarse crystal oscillators, the following equation is employed in EE-ASCFR to represent the hardware clock of the sensor node i: for i 0 , 1 , , N 1 ,
T i ( t ) = 1 + ϵ i t + θ i ,
where 1 + ϵ i R + and θ i R represent the clock frequency ratio and the clock offset, respectively [1]. The corresponding logical clock T i can be described based on Equation (1) as follows: For t k < t t k + 1 ( k = 0 , 1 , ),
T i T i ( t ) = T i T i ( t k ) + T i ( t ) T i ( t k ) 1 + ϵ ^ i , k θ ^ i , k ,
where the reference time for the kth synchronization is denoted by t k , the estimated clock skew and offset from the kth synchronization are described as ϵ ^ i , k and θ ^ i , k . Note that, as mentioned in [1], θ ^ i , k is set to 0 in Equation (2) since its compensation is done at the head node; in addition, the asynchronous SCFR scheme introduced in [42] is employed for the sensor node to synchronize its frequency of the logical clock to that of the reference clock. Furthermore, asymmetric high-precision time synchronization (AHTS) [43] is proposed as an improvement of EE-ASCFR, which is later leveraged in our proposed optimal bundling, to address the practical application and multi-hop scenario.
However, because the synchronization accuracy (SA) is affected by SI as illustrated in Table 1 (i.e., simulation results of EE-ASCFR and practical evaluation results of AHTS), we need to maintain SI to a reasonable threshold for better synchronization accuracy. Based on those evaluation results in Table 1, the relationship ( R ) between SI and synchronization accuracy is represented as follows:
S I = R ( S A ) ,
in which the requirement of SA could be translated to the requirement of SI.
The energy consumption caused by message transmissions could be drastically reduced by the aforementioned data bundling procedure as demonstrated in Figure 1: downstream, the synchronization “Request” message for timestamp T 1 is embedded inside a regular beacon message; upstream, the synchronization “Response” message for timestamps T 1 , T 2 , and T 3 is bundled together with measurement data as shown in Figure 2. Note that bundling more measurement data could lead to better energy efficiency, which, however, increases SI and decreases SA as discovered in the evaluation of [1] and [43]. Nevertheless, as previously mentioned, EE-ASCFR only covers the single-hop scenario, the more asymmetric approach (i.e., AHTS), which is based on EE-ASCFR and extends it to multi-hop resource-constrained sensor networks, is employed in the proposed approach to cover multi-hop scenarios. The data bundling procedure is inherited and extended to the multi-hop case as illustrated in Figure 3. Of particular note is that the bundling number is not discussed either in EE-ASCFR or AHTS since it directly affects the E2E delay and synchronization accuracy, which makes its value selection very complicated.

2.3. Conflicts of Performance Metrics: Energy Efficiency, E2E Delay and Synchronization Accuracy

Simultaneously meeting the requirements for the three performance metrics of energy efficiency, E2E delay, and synchronization accuracy is not possible due to their relationship discussed in Section 2.2. Time synchronization provides the possibility of accurate calculation and maintenance of E2E delay, however, employing the time synchronization scheme requires certain computing and power resources, which produces the conflict between energy efficiency and E2E delay calculation. Although the data bundling procedure introduced in EE-ASCFR could drastically reduce the energy consumption for message transmissions and bundling more data could lead to higher energy efficiency, the data bundling procedure could directly result in high E2E delay as revealed in Figure 3, which is the conflict between energy efficiency and E2E delay maintenance. Specific to the delay calculation illustrated in Figure 3, the E2E delay—i.e., T m r T m s —of measurement m, would be relatively small (e.g., multiples of forwarding delay which typically in milliseconds) in the regular direct forwarding method. However, due to the bundling procedure in the intermediate sensor nodes, the E2E delay could be as large as multiples of measurement interval. Nonetheless, the bundling procedure is shown to be quite efficient in conserving energy by reducing the number of message transmissions as shown in Figure 3, which is critical to low-power sensor nodes. As for synchronization accuracy, SI should be shorter for achieving higher accuracy [1], but shorter SI results in more message transmissions, which again leads to higher energy consumption. This draws forth the conflict between energy efficiency and synchronization accuracy.
To jointly meet the three performance requirements, comprehensive optimizations should be taken so that overall satisfactory performances among those three requirements could be achieved. In this paper, a new approach for the bundling number optimization for each sensor node is proposed to fulfill the requirement of maintaining E2E delay and synchronization accuracy while minimizing the message transmissions.

3. ILP Model for the Optimal Bundling Problem

In a WSN with N sensor nodes, the energy consumption e i t for the message transmissions at sensor node i is modeled as follows: For i 0 , 1 , , N 1 ,
e i t = α i e i m + β i e i s + γ i e i f ,
where e i m and e i s denote the energy consumption for the transmission of a measurement and a synchronization message generated by sensor node i, respectively, and e i f is the energy consumption for forwarding either a synchronization or a measurement message from offspring sensor nodes at sensor node i. The coefficients α i , β i and γ i are the number of transmissions for corresponding messages.
If we apply data bundling, Equation (4) can be modified as follows:
e i t = δ i e i b , all data bundling , δ i e i b + γ i e i f , self data bundling ,
where δ i e i b is the total energy consumption caused by transmitting the bundled messages. Note that, unlike the all data bundling option, which bundles all data into one bundled message, the self data bundling option bundles only the data generated from the sensor node itself, which is the reason the term γ i e i f still remains. From the network-level perspective, the total energy consumption for message transmissions in the network could be described as follows:
E t = i = 0 N 1 e i t .

3.1. Maximization of Bundling Number for Energy Efficiency

With the model for energy consumption caused by message transmissions, we can increase the network energy efficiency by minimizing the total energy consumption for message transmissions (i.e., E t ).
As shown in Equation (5), a large number of bundled message transmissions (i.e., δ i ) could result in more energy consumption, but it could be decreased by bundling more messages in one transmission. Let Γ i be the number of bundled messages at sensor node i. Then, a larger Γ i could lead to a smaller δ i , which could reduce total energy consumption for message transmissions. Consequently, we can minimize E t by maximizing Γ i .
Even though we can increase the amount of bundling for better energy efficiency, we cannot indefinitely because the E2E delay of the measurement data are sensitive to the number of bundling. Therefore, the maximization of bundling number for energy efficiency can be formulated as follows:
maximize Γ subject to χ m i n Γ i χ m a x , i [ 0 , , N 1 ] ,
where
Γ = i = 0 N 1 Γ i .
Note that Γ is the total bundling number in the network and that χ m i n and χ m a x are the lower and upper bounds of the measurement bundling number which are application-specific parameters and could be specified by the user.

3.2. Constraining E2E Delay

We define the E2E delay of sensor node i as the difference between the time of a certain measurement m at the sensor node and the time of the reception of the resulting message by the head:
D e 2 e i l = 0 L 1 D l i = T m i , r T m i , s ,
where D l i is the link delay at link l of L links from sensor node i to the head, T m i , r and T m i , s are the receiving time at the head and the measuring time at the sensor node respectively, the latter of which is a time with respect to the reference clock at the head translated by a time synchronization scheme. Specifically, the D e 2 e i is a path-level delay which consists of several link-level delays in the network. Considering the bundling procedure, the link delay at link l for sensor node i could be described as follows:
D l i = D p r o p i , l + D s e r v i , l + D b u n d i , l ,
where D p r o p i , l denotes the propagation delay, which is typically at the nanosecond level in WSN, D b u n d i , l is the delay caused by the bundling procedure which is in multiples of the measurement interval (e.g., 5 × 1 s for the measurement interval of 1 s and the bundling number of 5). Based on the service time model [39] for TinyOS, which is running on the TelosB sensor node, we can model D s e r v i , l as follows:
D s e r v i , l = D S P I i , l + D s u c c i , l + ( N t r y i , l 1 ) · D r e t r y i , l , N t r y i , l N m a x i , l , D S P I i , l + D f a i l i , l + ( N m a x i , l 1 ) · D r e t r y i , l , N t r y i , l > N m a x i , l ,
where
D s u c c i , l = D M A C i , l + D f r a m e i , l + D A C K i , l , D f a i l i , l = D M A C i , l + D f r a m e i , l + D w a i t A C K i , l , D r e t r y i , l = T r e t r y i , l + D f r a m e i , l + D w a i t A C K i , l .
Note that the delay parameters—i.e., one-time serial-peripheral interface (SPI) bus loading delay D s p i i , l , medium access control (MAC) layer delay D M A C i , l , frame transmission delay D f r a m e i , l , acknowledgment (ACK) transmission delay D A C K i , l and ACK waiting delay D w a i t A C K i , l —in the above equations are platform-dependent values, and their values are typically in the order of milliseconds. In addition, N t r y i , l and N m a x i , l are the current and the maximum allowed number of transmissions for a successful delivery, and T r e t r y i , l is the user-defined backup time of retransmission. For simplicity and energy efficiency, the packet retransmission is not taken into account in our proposed scheme since there are usually not many packet retransmissions in the network with lower traffic. Thus, Equation (10) could be simplified as follows:
D s e r v i , l = D S P I i , l + D M A C i , l + D f r a m e i , l + D A C K i , l ,
where the value of D s e r v i , l is around 10 based on the reference values in [39].
The link delay in Equation (9) could be further simplified when the measurement interval for an application ( I m e a s i , l ) is much larger than the service delay ( D s e r v i , l ): because D b u n d i , l > I m e a s i , l , I m e a s i , l D s e r v i , l also implies D b u n d i , l D s e r v i , l . Hence,
D l i n k i , l D b u n d i , l , if I m e a s i , l D s e r v i , l .
Note that the service delay D s e r v i , l should not be ignored in case the application requires frequent measurements, i.e., I m e a s i , l is comparable to D s e r v i , l .
In typical hierarchical multi-hop WSNs, sensor nodes located in different layers handle different amounts of traffic: for instance, the gateway node in the upper layer has to handle the message traffic from its offspring sensor nodes as well as itself; the higher layer it is located in, the more message traffic it has to handle. This means that even two sensor nodes with the same bundling number (i.e., Γ i ) could have different bundling delays due to the variance in their message traffic. By introducing a message traffic coefficient ( 1 1 + λ i ) for each sensor node that periodically measures data with the same measurement interval (i.e., I m e a s i ), the bundling delay at sensor node i ( D b u n d i ) could be represented as follows:
D b u n d i = Γ i 1 + λ i · I m e a s i ,
where λ i denotes the number of offspring sensor nodes. Then, the D e 2 e i for the applications with normal measurement interval could be modeled as follows:
D e 2 e i = l = 0 L 1 D b u n d i .
With Equation (14), we can also constrain the E2E delay in the optimal bundling problem in Equation (7) with the user-defined E2E delay requirement ( D e 2 e m a x ): for i [ 0 , 1 , , N 1 ] ,
D e 2 e i D e 2 e m a x .

3.3. Constraining Synchronization Accuracy

Since the proposed optimal message bundling approach is based on the reverse asymmetric time synchronization scheme, the synchronization accuracy depends on the report interval of the synchronization messages, which, in turn, are carried by the bundled messages. In such a case, the user-required synchronization accuracy ( SA m i n ) could be maintained through constraining the E2E delay of the bundled message when sensor nodes generate measurement data periodically. Based on the empirical sets—i.e., R —of the relationships between synchronization accuracy and SI previously provided in Section 2, the user-required synchronization accuracy could be translated to the delay requirement as follows:
D e 2 e S A = R ( SA m i n ) .
Combining Equations (14) and (16), the synchronization accuracy could be achieved through constraining the E2E delay as follows: For i [ 0 , 1 , , N 1 ] ,
D e 2 e i D e 2 e S A .

3.4. ILP Model

Combining the objective function Equation (7) and the two constraint sets Equations (15) and (17), we can formulate the optimal bundling problem as the following ILP:
maximize Γ = i = 0 N 1 Γ i , subject to χ m i n Γ i χ m a x , i [ 0 , , N 1 ] , D e 2 e i min D e 2 e m a x , R ( SA m i n ) , i [ 0 , , N 1 ] ,
where the E2E delays of all sensor nodes are jointly constrained by the user-defined E2E delay D e 2 e m a x and the synchronization accuracy SA m i n requirements, and the bundling number of each sensor node is constrained by the user-defined lower χ m i n and upper χ m a x bounds, respectively. Applying the set of optimal bundling numbers computed from this ILP model to the sensor nodes, the required E2E delay and synchronization accuracy could be jointly maintained while the bundled message transmissions are still minimized.

4. System Design

Figure 4 shows a system architecture for the proposed optimal bundling based on the ILP model formulated in Section 3.4, where the two subsystems—i.e., the Performance Maintainer at the head and the Parameter Adapter at each sensor node—are built to achieve the optimization target. This separation of components further releases the computational burden of the sensor node: the head is in charge of monitoring the runtime topology and computing the optimal bundling number for each sensor node dynamically; the sensor node, therefore, is relieved from the complex computation and leverages the optimal result from the head straightforwardly.

4.1. Performance Maintainer at Head

Since the time synchronization operation is centralized at the head in EE-ASCFR and AHTS, we first build the reference time synchronization system at the head as a component of the proposed system. The time synchronization is achieved through the MAC-Layer Time Recorder and the Time Synchronization Maintainer, the latter of which translates timestamps between the head and a sensor node. With the time synchronization component, a measurement timestamp recorded at the sensor node—i.e., T m s in Figure 3—is translated to a timestamp based on the hardware clock of the head. Then, the E2E delay is obtained as a difference between the translated measurement timestamp and the receiving timestamp of T m r through the Delay Calculator. Afterwards, the runtime E2E delay is monitored through the Performance Monitor.
Based on the topology data carried in each bundled message, the routing paths for all sensor nodes are recovered in the Runtime Path Maintainer, and one set of paths is generated and delivered to the Constraint Generator. The user-defined requirements of E2E delay and synchronization accuracy are captured through the user interface of Requirement Input Interface. By combining the performance requirements and the current path information, the Constraint Generator generates a set of constraints and passes it to the Optimal Parameter Generator, where the optimal bundling number for each sensor node is obtained as a solution of the ILP model. Finally, the optimal bundling numbers from the Optimal Parameter Generator are delivered to sensor nodes by the Parameter Disseminator.
Of particular note is that, whenever the topology is updated or the users change their requirements dynamically, the whole process above will be operated to provide up-to-date optimal bundling numbers. In addition, since the proposed system is in a centralized manner, the major trade-off is the communication overhead for disseminating the optimal results to the sensor nodes. This overhead, by the way, would not be a major issue in practice unless topology updates and user requirements changes are too frequent.

4.2. Parameter Adapter at Sensor Nodes

To reduce the computational complexity required by the proposed optimal bundling on the sensor nodes, three lightweight components (including the MAC-Layer Time Recorder from the time synchronization scheme) are implemented at sensor nodes. The Parameter Adapter receives the optimal bundling number and delivers it to the Data Bundler which bundles that number of measurement data temporarily stored in the Queue plus timestamps into one message. Then, the bundled message will be again timestamped for T 3 by the MAC-Layer Time Recorder as shown in Figure 1 and Figure 4.

5. Experimental Results

The proposed approach is implemented on a real three-hop WSN testbed consisting of five TelosB sensor nodes as illustrated in Figure 5. To verify the time synchronization accuracy of the time synchronization approach embedded in the proposed optimal bundling, we establish the time synchronization accuracy measurement based on the setup illustrated in Figure 6. Similar to the approach of Reference-Broadcast Synchronization (RBS) [44], one reference node is employed to broadcast the reference message to the head and sensor node under evaluation to record the receiving timestamps. With neglecting the propagation delay, those receiving timestamps should be at the same time point, thanks to the leveraging of the MAC-layer timestamping [43]. Afterwards, the recorded timestamps in the sensor node under evaluation will be sent to the head for translation. Since the head and sensor node under evaluation are synchronized through the time synchronization approach, so the synchronization error could be computed in the head through comparing the differences of the head’s recorded timestamp and the translated one.
During the experiments, each sensor node generates one measurement per second, and the E2E delays of the latest measurements in bundled messages from all sensor nodes are collected and stored in a time sequence in order of their arrivals at the head, which are shown in Figure 7, Figure 8 and Figure 9. The red horizontal dotted lines indicate the E2E delay requirements for corresponding time periods, and the requirement of synchronization accuracy is set to 5 for all experiments.

5.1. Delay Performance under Static Requirement Setting

We first evaluate the E2E delay performance of the optimal bundling under static requirement setting—e.g., E2E delay requirement of 8 s and 5 s —with different maximum bundling numbers—i.e., χ m a x in Equation (18)—of 15 and 10. We run the experiment for 3600 s to demonstrate the long-term maintenance capability.
The performance of the proposed approach under the static requirement setting is explicitly exhibited through diverse experiments with different parameters and requirements. As shown in Figure 7a,b, the E2E delays can exceed 14 s before the optimal bundling is applied. Once the optimal bundling is applied with the requirement of 8 s and 5 s ; however, we can see that the E2E delay is controlled and kept under the requirements for most of the time. Note that there are few data points crossed the requirement line; because the gateway node serves its own measurement data first, the data from its offspring nodes, sometimes, could be buffered in the queue and sent by the next available message, which would increase the E2E delay of the corresponding message. Next, we change the maximum bundling number from 15 to 10 to evaluate the proposed optimal bundling approach since the maximum bundling number is often limited by the measurement data length and the maximum payload size of the underlying protocol. Figure 8a,b illustrate that the E2E delay requirement can be well fulfilled with applying the proposed optimal bundling approach.

5.2. Delay Performance under Dynamic Requirement Setting

We also evaluate the E2E delay performance of the optimal bundling under dynamic requirement setting with the maximum bundling number of 10 and 15 to demonstrate its run-time maintenance capability.
During the evaluation with maximum bundling number of 10, the E2E delay requirement is dynamically changed from 8 s to 2 s step-by-step as shown in Figure 9a. In addition, for the evaluation with maximum bundling number of 15 shown in Figure 9b, we extend the range of E2E delay requirement to 14 s to 1 s to further illustrate the performance under the dynamic scenario. The results demonstrate that the proposed optimal bundling nicely handles multiple dynamic requirements of E2E delay throughout the experiments. As discussed in Section 5.1, however, some data points slightly cross the requirement line, which is due to the neglect of the service time (i.e., D s e r v i , l ) in Equation (9). When the optimal bundling numbers are too strict, which just fulfills the requirement, the influence of the neglect of the small service time would be notable. Note that, in the extreme case, when the E2E delay is strict to 1 s as shown in Figure 9b, which reaches the measurement interval of our experiment as previously described in Section 5, the proposed optimal bundling could still fulfill the requirement but with assigning no larger optimal bundling number than 1 to all sensor nodes as shown in Table 2. However, when the E2E delay requirement is much stricter, i.e., smaller than the measurement interval, this may lead to an unsatisfactory situation since there is no room for the proposed optimal bundling to play.

5.3. Synchronization Accuracy

With the novel time synchronization schemes of EE-ASCFR and AHTS, the synchronization accuracy could be fulfilled as far as the E2E delay of the bundled message satisfying the synchronization interval: the synchronization interval—i.e., E2E delay—exhibited in the evaluation results as exhibited in Figure 7, Figure 8 and Figure 9, could overfulfill the SI requirement (e.g., 10 s SI could lead to 2 . 3385 μs synchronization accuracy as illustrated in Table 1), intuitively the synchronization accuracy could be strictly followed. Furthermore, we take node 2 in the static requirement setting experiment—i.e., experiment shown in Figure 8a with maximum bundling number of 10 and E2E delay requirement of 8 s —as an example, to practically evaluate the time synchronization accuracy. As illustrated in Figure 10, most of the absolute synchronization errors (i.e., the dotted line) are under 5 μ s , and the mean of absolute synchronization errors is 2 . 0849 μ s , which could overfulfill the 5 synchronization accuracy requirement.

5.4. Energy Efficiency

To evaluate the energy efficiency, we indirectly estimate it by comparing the number of message receptions and transmissions with and without optimal bundling. Taking the path of 0 1 2 4 as an example, the number of message receptions and transmissions of the sensor nodes 1 , 2 , 4 in the experiment—i.e., experiment exhibited in Figure 7a—of static E2E delay requirement setting is illustrated in Figure 11.
The number of message receptions and transmissions of each sensor node is counted every 10 s over the period of 3600 s . With the proposed optimal bundling, all the sensor nodes could maintain their message receptions and transmissions around the number of 10. Without the optimal bundling, on the other hand, the numbers of their message receptions and transmissions are relatively larger, whose average of 35 is more than triple that of the optimal bundling. In other words, around 70% of message transmissions can be reduced using the proposed optimal bundling in this specific case. Of particular note is that, with the increase of the network size and the maximum bundling number, the performance of the proposed optimal bundling would be even better.

5.5. Discussion

As illustrated in Table 2, the optimal bundling numbers that our proposed approach assign to the leaf nodes (i.e., nodes 3 and 4 in our experiment as shown in Figure 5) are relatively small regardless of the variation of the E2E delay requirements. This is because the proposed optimal bundling treats both the gateway node and leaf node equally since they are most likely the same resource-constrained sensor nodes in the real deployment and the gateway nodes undertake more traffic than the leaf ones. However, when the gateway node is resource-abundant rather than resource-constrained, consequently the proposed approach should incline the optimization to the leaf nodes.
In addition, as briefly mentioned in Section 5.4, the maximum energy efficiency performance of the proposed optimal bundling is not exhibited due to the joint constraints of the network size and various parameters such as maximum bundling number. In the follow-up work, a further study will employ a relative larger testbed and proper parameters to evaluate the extreme performance of the proposed optimal bundling.

6. Conclusions

We have proposed an approach to optimize the number of bundled messages at sensor nodes in a WSN under the constraints of time synchronization accuracy and E2E delay. In order to solve the optimal bundling problem, we formulate it as an ILP model and employ the novel asymmetric time synchronization scheme. To the best of the authors’ knowledge, this is the first work to optimize message bundling for energy efficiency under the joint constraints of synchronization accuracy and E2E delay in the context of WSNs. The practical evaluation results on a real testbed composed of TelosB sensor nodes demonstrate the long-term and runtime maintenance capability of the proposed approach in terms of E2E delay. The energy efficiency is also illustrated through a case study, showing about 70% reduction on the overall message transmissions.
Note that bundling a larger number of messages results in a longer payload, which may lead to possible link degradation such as the packet reception ratio degradation. In this regard, link quality requirements could be introduced to the optimization model as additional constraints to take into account more impacts on the overall performance by the bundling procedure. In addition, the proposed approach should be tested with a larger testbed to investigate its scalability and network-wide energy efficiency.

Author Contributions

Methodology, X.H. and K.S.K.; development, X.H.; software, X.H.; validation, X.H. and K.S.K.; writing—original draft preparation, X.H.; writing—review and editing, K.S.K., S.L. and M.K.K.; funding acquisition, K.S.K., S.L. and M.K.K.

Funding

This work was supported in part by the Xi’an Jiaotong-Liverpool University Research Development Fund (under Grant RDF-16-02-39), Key Program Special Fund (under Grant KSF-E-25), and Centre for Smart Grid and Information Convergence (CeSGIC).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kim, K.S.; Lee, S.; Lim, E.G. Energy-Efficient Time Synchronization Based on Asynchronous Source Clock Frequency Recovery and Reverse Two-Way Message Exchanges in Wireless Sensor Networks. IEEE Trans. Commun. 2017, 65, 347–359. [Google Scholar] [CrossRef]
  2. Cui, J.; Valois, F. Data Aggregation in Wireless Sensor Networks: Compressing or Forecasting? In Proceedings of the 2014 IEEE Wireless Communications and Networking Conference (WCNC), Istanbul, Turkey, 6–9 April 2014; pp. 2892–2897. [Google Scholar] [CrossRef]
  3. Rajagopalan, R.; Varshney, P.K. Data-aggregation techniques in sensor networks: A survey. IEEE Commun. Surv. Tutor. 2006, 8, 48–63. [Google Scholar] [CrossRef]
  4. Li, X.; Liu, W.; Xie, M.; Liu, A.; Zhao, M.; Xiong, N.N.; Zhao, M.; Dai, W. Differentiated Data Aggregation Routing Scheme for Energy Conserving and Delay Sensitive Wireless Sensor Networks. Sensors 2018, 18. [Google Scholar] [CrossRef] [PubMed]
  5. Li, X.; Liu, A.; Xie, M.; Xiong, N.N.; Zeng, Z.; Cai, Z. Adaptive Aggregation Routing to Reduce Delay for Multi-Layer Wireless Sensor Networks. Sensors 2018, 18. [Google Scholar] [CrossRef] [PubMed]
  6. Villas, L.A.; Boukerche, A.; Ramos, H.S.; de Oliveira, H.A.B.F.; de Araujo, R.B.; Loureiro, A.A.F. DRINA: A Lightweight and Reliable Routing Approach for In-Network Aggregation in Wireless Sensor Networks. IEEE Trans. Comput. 2013, 62, 676–689. [Google Scholar] [CrossRef]
  7. Bahi, J.M.; Makhoul, A.; Medlej, M. An Optimized In-Network Aggregation Scheme for Data Collection in Periodic Sensor Networks. In Ad-Hoc, Mobile, and Wireless Networks, Proceedings of the 11th International Conference, ADHOC-NOW 2012, Belgrade, Serbia, 9–11 July 2012; Spring: Berlin/Heidelberg, Germany, 2012; pp. 153–166. [Google Scholar]
  8. Intanagonwiwat, C.; Estrin, D.; Govindan, R.; Heidemann, J. Impact of Network Density on Data Aggregation in Wireless Sensor Networks. In Proceedings of the IEEE 22nd International Conference on Distributed Computing Systems, Vienna, Austria, 2–5 July 2002; pp. 457–458. [Google Scholar] [CrossRef]
  9. Madden, S.; Franklin, M.J.; Hellerstein, J.M.; Hong, W. The Design of an Acquisitional Query Processor for Sensor Networks. In Proceedings of the 2003 ACM SIGMOD International Conference on Management of Data, San Diego, CA, USA, 9–12 June 2003; pp. 491–502. [Google Scholar]
  10. Atoui, I.; Ahmad, A.; Medlej, M.; Makhoul, A.; Tawbe, S.; Hijazi, A. Tree-Based Data Aggregation approach in Wireless Sensor Network Using Fitting Functions. In Proceedings of the 2016 Sixth International Conference on Digital Information Processing and Communications (ICDIPC), Beirut, Lebanon, 21–23 April 2016; pp. 146–150. [Google Scholar] [CrossRef]
  11. Xu, X.; Li, X.Y.; Mao, X.; Tang, S.; Wang, S. A Delay-Efficient Algorithm for Data Aggregation in Multihop Wireless Sensor Networks. IEEE Trans. Parallel Distrib. Syst. 2011, 22, 163–175. [Google Scholar] [CrossRef]
  12. YOSHINO, H.; OTA, K.; HIRAGURI, T. Queueing Delay Analysis and Optimization of Statistical Data Aggregation and Transmission Systems. IEICE Trans. Commun. 2018, E101.B, 2186–2195. [Google Scholar] [CrossRef]
  13. Deng, X.; Yang, Y. Online Adaptive Compression in Delay Sensitive Wireless Sensor Networks. IEEE Trans. Comput. 2012, 61, 1429–1442. [Google Scholar] [CrossRef]
  14. Al-Anbagi, I.; Erol-Kantarci, M.; Mouftah, H.T. Delay Critical Smart Grid Applications and Adaptive QoS Provisioning. IEEE Access 2015, 3, 1367–1378. [Google Scholar] [CrossRef]
  15. Doudou, M.; Djenouri, D.; Badache, N. Survey on Latency Issues of Asynchronous MAC Protocols in Delay-Sensitive Wireless Sensor Networks. IEEE Commun. Surv. Tutor. 2013, 15, 528–550. [Google Scholar] [CrossRef]
  16. Wang, Y.; Vuran, M.C.; Goddard, S. Cross-Layer Analysis of the End-to-End Delay Distribution in Wireless Sensor Networks. IEEE/ACM Trans. Netw. 2012, 20, 305–318. [Google Scholar] [CrossRef]
  17. Wang, J.; Dong, W.; Cao, Z.; Liu, Y. On the Delay Performance in a Large-Scale Wireless Sensor Network: Measurement, Analysis, and Implications. IEEE/ACM Trans. Netw. 2015, 23, 186–197. [Google Scholar] [CrossRef]
  18. Liu, Y.; Zhang, D.; Guo, X.; Gao, M.; Ming, Z.; Yang, L.; Ni, L.M. RSS-Based Ranging by Leveraging Frequency Diversity to Distinguish the Multiple Radio Paths. IEEE Trans. Mob. Comput. 2017, 16, 1121–1135. [Google Scholar] [CrossRef]
  19. Zhang, D.; Liu, Y.; Guo, X.; Ni, L.M. RASS: A Real-Time, Accurate, and Scalable System for Tracking Transceiver-Free Objects. IEEE Trans. Parallel Distrib. Syst. 2013, 24, 996–1008. [Google Scholar] [CrossRef]
  20. Kim, B.; Hong, S.; Hur, K.; Eom, D. Energy-efficient and rapid time synchronization for wireless sensor networks. IEEE Trans. Consum. Electron. 2010, 56, 2258–2266. [Google Scholar] [CrossRef]
  21. Gong, F.; Sichitiu, M.L. CESP: A Low-Power High-Accuracy Time Synchronization Protocol. IEEE Trans. Veh. Technol. 2016, 65, 2387–2396. [Google Scholar] [CrossRef]
  22. Nadas, J.P.B.; Souza, R.D.; Pellenz, M.E.; Brante, G.; Braga, S.M. Energy Efficient Beacon Based Synchronization for Alarm Driven Wireless Sensor Networks. IEEE Signal Process. Lett. 2016, 23, 336–340. [Google Scholar] [CrossRef]
  23. Al-Shaikhi, A.; Masoud, A. Efficient, Single Hop Time Synchronization Protocol for Randomly Connected WSNs. IEEE Wirel. Commun. Lett. 2017, 6, 170–173. [Google Scholar] [CrossRef] [Green Version]
  24. Tavares Bruscato, L.; Heimfarth, T.; Pignaton de Freitas, E. Enhancing Time Synchronization Support in Wireless Sensor Networks. Sensors 2017, 17. [Google Scholar] [CrossRef]
  25. Wang, Y.; Nunez, F.; Doyle, F.J. Energy-Efficient Pulse-Coupled Synchronization Strategy Design for Wireless Sensor Networks Through Reduced Idle Listening. IEEE Trans. Signal Process. 2012, 60, 5293–5306. [Google Scholar] [CrossRef] [Green Version]
  26. Akhlaq, M.; Sheltami, T.R. RTSP: An Accurate and Energy-Efficient Protocol for Clock Synchronization in WSNs. IEEE Trans. Instrum. Meas. 2013, 62, 578–589. [Google Scholar] [CrossRef]
  27. Vieira, M.A.M.; Coelho, C.N.; da Silva, D.C.; da Mata, J.M. Survey on Wireless Sensor Network Devices. In Proceedings of the 2003 IEEE Conference on Emerging Technologies and Factory Automation, Lisbon, Portugal, 16–19 September 2003; pp. 537–544. [Google Scholar]
  28. Lin, S.; Miao, F.; Zhang, J.; Zhou, G.; Gu, L.; He, T.; Stankovic, J.A.; Son, S.; Pappas, G.J. ATPC: Adaptive Transmission Power Control for Wireless Sensor Networks. ACM Trans. Sen. Netw. 2016, 12, 6:1–6:31. [Google Scholar] [CrossRef]
  29. Chincoli, M.; Liotta, A. Self-Learning Power Control in Wireless Sensor Networks. Sensors 2018, 18. [Google Scholar] [CrossRef] [PubMed]
  30. Pal, A.; Nasipuri, A. Joint Power Control and Routing for Rechargeable Wireless Sensor Networks. IEEE Access 2019, 123992–124007. [Google Scholar] [CrossRef]
  31. Chipara, O.; He, Z.; Xing, G.; Chen, Q.; Wang, X.; Lu, C.; Stankovic, J.; Abdelzaher, T. Real-Time Power-Aware Routing in Sensor Networks. In Proceedings of the 2006 14th IEEE International Workshop on Quality of Service, New Haven, CT, USA, 19–21 June 2006; pp. 83–92. [Google Scholar] [CrossRef]
  32. Wu, Y.; Stankovic, J.A.; He, T.; Lin, S. Realistic and Efficient Multi-Channel Communications in Wireless Sensor Networks. In Proceedings of the IEEE INFOCOM 2008—The 27th Conference on Computer Communications, Phoenix, AZ, USA, 13–18 April 2008; pp. 1193–1201. [Google Scholar] [CrossRef]
  33. Pal, A.; Nasipuri, A. Distributed Routing and Channel Selection for Multi-Channel Wireless Sensor Networks. J. Sens. Actuator Netw. 2017, 6. [Google Scholar] [CrossRef]
  34. Abdeddaim, N.; Theoleyre, F.; Heusse, M.; Duda, A. Adaptive IEEE 802.15.4 MAC for Throughput and Energy Optimization. In Proceedings of the 2013 IEEE International Conference on Distributed Computing in Sensor Systems, Cambridge, MA, USA, 20–23 May 2013; pp. 223–230. [Google Scholar] [CrossRef]
  35. Zlmmerling, M.; Ferrari, F.; Mottola, L.; Voigt, T.; Thiele, L. pTUNES: Runtime Parameter Adaptation for Low-power MAC Protocols. In Proceedings of the 2012 ACM/IEEE 11th International Conference on Information Processing in Sensor Networks (IPSN), Beijing, China, 16–20 April 2012; pp. 173–184. [Google Scholar] [CrossRef]
  36. Noda, C.; Prabh, S.; Alves, M.; Voigt, T. On Packet Size and Error Correction Optimisations in Low-power Wireless Networks. In Proceedings of the 2013 IEEE International Conference on Sensing, Communications and Networking (SECON), New Orleans, LA, USA, 24–27 June 2013; pp. 212–220. [Google Scholar] [CrossRef]
  37. Dong, W.; Chen, C.; Liu, X.; He, Y.; Liu, Y.; Bu, J.; Xu, X. Dynamic Packet Length Control in Wireless Sensor Networks. IEEE Trans. Wirel. Commun. 2014, 13, 1172–1181. [Google Scholar] [CrossRef]
  38. Pal, A.; Kant, K. On the Feasibility of Distributed Sampling Rate Adaptation in Heterogeneous and Collaborative Wireless Sensor Networks. In Proceedings of the 2016 25th International Conference on Computer Communication and Networks (ICCCN), Waikoloa, HI, USA, 1–4 August 2016; pp. 1–9. [Google Scholar] [CrossRef]
  39. Fu, S.; Zhang, Y.; Jiang, Y.; Hu, C.; Shih, C.; Marrón, P.J. Experimental Study for Multi-layer Parameter Configuration of WSN Links. In Proceedings of the IEEE 35th International Conference on Distributed Computing Systems, Columbus, OH, USA, 29 June–2 July 2015; pp. 369–378. [Google Scholar] [CrossRef]
  40. Battestilli, T.; Perros, H. An introduction to optical burst switching. IEEE Commun. Mag. 2003, 41, S10–S15. [Google Scholar] [CrossRef]
  41. Pal, A.; Kant, K. E-Darwin2: A Smartphone Based Disaster Recovery Network Using WiFi Tethering. In Proceedings of the 2018 15th IEEE Annual Consumer Communications Networking Conference (CCNC), Las Vegas, NV, USA, 12–15 January 2018; pp. 1–5. [Google Scholar] [CrossRef]
  42. Kim, K.S. Asynchronous Source Clock Frequency Recovery through Aperiodic Packet Streams. IEEE Commun. Lett. 2013, 17, 1455–1458. [Google Scholar] [CrossRef] [Green Version]
  43. Huan, X.; Kim, K.S. On the Practical Implementation of Propagation Delay and Clock Skew Compensated High-Precision Time Synchronization Schemes with Resource-Constrained Sensor Nodes in Multi-Hop Wireless Sensor Networks. arXiv 2019, arXiv:1905.00554. [Google Scholar]
  44. Elson, J.; Girod, L.; Estrin, D. Fine-Grained Network Time Synchronization Using Reference Broadcasts. In Proceedings of the 5th Symposium on Operating Systems Design and Implementation, New York, NY, USA, 9–11 December 2002; pp. 147–163. [Google Scholar]
Figure 1. The reverse two-way message exchange and data bundling illustrated in EE-ASCFR [1].
Figure 1. The reverse two-way message exchange and data bundling illustrated in EE-ASCFR [1].
Sensors 19 04027 g001
Figure 2. Payload contents of the data bundling procedure introduced in EE-ASCFR [1].
Figure 2. Payload contents of the data bundling procedure introduced in EE-ASCFR [1].
Sensors 19 04027 g002
Figure 3. Comparison of the message transmissions of the time synchronization schemes with and without data bundling procedure.
Figure 3. Comparison of the message transmissions of the time synchronization schemes with and without data bundling procedure.
Sensors 19 04027 g003
Figure 4. System architecture of the proposed optimal bundling.
Figure 4. System architecture of the proposed optimal bundling.
Sensors 19 04027 g004
Figure 5. Experiment setup of a real three-hop wireless sensor network (WSN) testbed consisting of five TelosB sensor nodes.
Figure 5. Experiment setup of a real three-hop wireless sensor network (WSN) testbed consisting of five TelosB sensor nodes.
Sensors 19 04027 g005
Figure 6. Evaluation of time synchronization accuracy using one additional reference node.
Figure 6. Evaluation of time synchronization accuracy using one additional reference node.
Sensors 19 04027 g006
Figure 7. End-to-end (E2E) delay performance of the optimal bundling under static E2E delay requirement setting of (a) 8 s and (b) 5 s with the maximum bundling number of 15.
Figure 7. End-to-end (E2E) delay performance of the optimal bundling under static E2E delay requirement setting of (a) 8 s and (b) 5 s with the maximum bundling number of 15.
Sensors 19 04027 g007
Figure 8. End-to-end (E2E) delay performance of the optimal bundling under static E2E delay requirement setting of (a) 8 s and (b) 5 s with the maximum bundling number of 10.
Figure 8. End-to-end (E2E) delay performance of the optimal bundling under static E2E delay requirement setting of (a) 8 s and (b) 5 s with the maximum bundling number of 10.
Sensors 19 04027 g008aSensors 19 04027 g008b
Figure 9. End-to-end (E2E) delay performance of the optimal bundling under dynamic E2E delay requirement setting and the maximum bundling number of (a) 10 and (b) 15.
Figure 9. End-to-end (E2E) delay performance of the optimal bundling under dynamic E2E delay requirement setting and the maximum bundling number of (a) 10 and (b) 15.
Sensors 19 04027 g009
Figure 10. Absolute time synchronization error of node 2 during the experiment of static requirement setting with maximum bundling number of 10 and end-to-end (E2E) delay requirement of 8 s illustrated in Figure 8a.
Figure 10. Absolute time synchronization error of node 2 during the experiment of static requirement setting with maximum bundling number of 10 and end-to-end (E2E) delay requirement of 8 s illustrated in Figure 8a.
Sensors 19 04027 g010
Figure 11. Number of message receptions and transmissions of sensor nodes with and without optimal message bundling for 3600   s with sampling interval of 10   s .
Figure 11. Number of message receptions and transmissions of sensor nodes with and without optimal message bundling for 3600   s with sampling interval of 10   s .
Sensors 19 04027 g011
Table 1. MAE and MSE of measurement time estimation of EE-ASCFR and AHTS with different synchronization intervals (SIs) provided in [1] and [43].
Table 1. MAE and MSE of measurement time estimation of EE-ASCFR and AHTS with different synchronization intervals (SIs) provided in [1] and [43].
Synchronization SchemeMAE 1MSE 2
EE-ASCFRSI = 100 s 8.8811 × 10 25 5.8990× 10 19
SI = 1 s 9.1748 × 10 25 5.4210 × 10 19
SI = 10 m s 1.0887 × 10 24 4.7684 × 10 19
AHTSSI = 100 s 8.4225 × 10 6 1.2524 × 10 10
SI = 10 s 2.3385 × 10 6 9.1694 × 10 12
SI = 1 s 1.8166 × 10 6 5.2094 × 10 12
1 MAE is the mean absolute error of measurement time estimation. 2 MSE is the mean square error of measurement time estimation.
Table 2. Optimal bundling number under different end-to-end (E2E) delay requirements during the dynamic experiment as shown in Figure 9b.
Table 2. Optimal bundling number under different end-to-end (E2E) delay requirements during the dynamic experiment as shown in Figure 9b.
Max Bundling NumberE2E Delay Requirement (s)Optimal Bundling Number
Node 1Node 2Node 3Node 4
1514151481
1510151061
15815641
15615221
15414111
1511111

Share and Cite

MDPI and ACS Style

Huan, X.; Kim, K.S.; Lee, S.; Kim, M.K. Optimal Message Bundling with Delay and Synchronization Constraints in Wireless Sensor Networks. Sensors 2019, 19, 4027. https://0-doi-org.brum.beds.ac.uk/10.3390/s19184027

AMA Style

Huan X, Kim KS, Lee S, Kim MK. Optimal Message Bundling with Delay and Synchronization Constraints in Wireless Sensor Networks. Sensors. 2019; 19(18):4027. https://0-doi-org.brum.beds.ac.uk/10.3390/s19184027

Chicago/Turabian Style

Huan, Xintao, Kyeong Soo Kim, Sanghyuk Lee, and Moon Keun Kim. 2019. "Optimal Message Bundling with Delay and Synchronization Constraints in Wireless Sensor Networks" Sensors 19, no. 18: 4027. https://0-doi-org.brum.beds.ac.uk/10.3390/s19184027

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