Abstract

Energy hole is an inherent problem caused by heavier traffic loads of sensor nodes nearer the sink because of more frequent data transmission, which is strongly dependent on the topology induced by the sensor deployment. In this paper, we propose an autonomous sensor redeployment algorithm to balance energy consumption and mitigate energy hole for unattended mobile sensor networks. First, with the target area divided into several equal width coronas, we present a mathematical problem modeling sensor node layout as well as transmission pattern to maximize network coverage and reduce communication cost. And then, by calculating the optimal node density for each corona to avoid energy hole, a fully distributed movement algorithm is proposed, which can achieve an optimal distribution quickly only by pushing or pulling its one-hop neighbors. The simulation results demonstrate that our algorithm achieves a much smaller average moving distance and a much longer network lifetime than existing algorithms and can eliminate the energy hole problem effectively.

1. Introduction

A typical wireless sensor network (WSN) is composed of hundreds of sensor nodes reporting their data to the information collector, referred to as the sink node. Sensor nodes are usually of low cost and low power, having limited sensing, computing, and communication capabilities. In recent years, with the rapid progress in advanced VLSI and radio frequency (RF) technologies, WSNs have attracted lots of interest due to their potential use in various applications such as military surveillance, target tracking, emergency navigation, and large scale systems [1, 2].

Sensor deployment is an important issue in designing a WSN since it affects the communication cost, detection capability, coverage, and connectivity [3]. In general, the resource-limited sensor nodes are dropped from airborne vehicles for remote surveying of unattended environment without a preconfigured infrastructure. Those sensor nodes are typically left unattended and remain static after initial deployment. Sometimes, establishing such a fixed sensor network over a hostile or dangerous area to provide complete coverage could be a daunting task. Thus, there exists an urgent need for exploiting sensor mobility in WSNs to improve network performance.

As sensor nodes are usually battery driven, they can survive for only a limited lifetime with nonrenewable batteries. Taking it one step further, the limited constraints of the sensor nodes restrict the use of high complexity algorithms and protocols. How to balance energy consumption is one of the fundamental issues arising in WSN. To address this issue, much work has been done during recent years. Among them, taking advantage of sensor mobility to enhance network lifetime has attracted extensive attention [411]. Typically, most of these traditional approaches aimed at having uniform deployment to achieve full coverage using a minimum number of sensor nodes. However, these approaches do not consider the issue of unbalanced energy depletion with distance to a predetermined sink. As the source nodes sent their sensing messages to the sink via a multihop relay tree in WSN, this many-to-one communication pattern could lead to traffic imbalance among sensor nodes. It has been observed that the closer a sensor is to the sink, the faster its battery exhausts, as time evolves, which will cause an energy hole in the region near the sink [12]. If this happens, no more data could be transmitted to the sink and the network is inactive soon leaving much energy unused, regardless of how many sensor nodes are deployed. Experimental results in [13] show that when the network lifetime ends, nearly 90% of the total initial energy is unused with uniform distribution. Therefore, managing sensor mobility to achieve both optimal energy consumption and full coverage is important to prolong the network lifetime for WSN.

In this paper, we try to solve the energy hole problem by proposing an autonomous coverage-driven sensor redeployment scheme. We first develop the energy hole problem with nonuniform node distribution in WSN theoretically. By importing the energy-aware transmission mechanism and the accessibility condition of energy-balanced depletion in our pervious approach [14], we further propose a fully distributed density control scheme in different coronas to balance energy depletion for the entire network. The rest of the paper is organized as follows. Section 2 reviews the related literature. Section 3 describes the network model, assumptions, and energy-aware transmission mechanism. A novel sensor redeployment strategy is proposed in Section 4. Section 5 presents simulation results for our algorithms, and Section 6 concludes our paper.

More recently, there has been growing interest in optimizing the sensor movement to maintain full coverage and prolong the network lifetime for mobile sensor networks. In [6], the authors proposed a potential field-based deployment algorithm, in which all the nodes explore from a compact region and fill the maximum working area in a way similar to the particles in the microworld [15]. In [7], the authors assumed that there were virtually attractive and repulsive forces among sensors. Using these virtual forces, mobile sensor nodes spread throughout the target area with a uniform distribution such that the network coverage rate is maximized. In [8], the authors proposed a Voronoi diagram-based distribution model, in which each sensor iteratively calculates its Voronoi polygon to detect the coverage holes and moves to a better position to improve the coverage rate. In [9], the authors proposed three independent algorithms (VEC, VOR, and MiniMax), by pushing or pulling nodes to cover the gaps based on virtual forces. These three algorithms have comparable performance in a bounded area, whereas only VEC algorithm can be used in both unbounded and bounded areas. In [10], the authors investigate how to move sensors to some locations while still preserving the degree of coverage under partially controlled placement. In [11], the sensing field was split into grids, and the sensors moved from high-density grids to low-density ones such that the densities remain constant. Overall, most of these traditional algorithms intended to redistribute sensor nodes uniformly and thus to maximize coverage rate, minimize coverage overlap or gap, and reduce the network cost for WSN. However, as the uniform distribution has unbalanced communication traffic, these approaches will cause the network lifetime to end prematurely with a great amount of energy unused.

In [12], the authors investigated the problem of uneven energy consumption in many-to-one WSNs for the first time. Further, they proposed several approaches to mitigate this problem and inferred that simply increasing the number of nodes under a uniform distribution could not prolong the network lifetime [16]. In [17], the authors focused on the nonuniform energy distribution among sensors. It was concluded that the higher the workload was, the higher the initial level of energy would be set. Although such a strategy seems to be promising, its application may be so difficult that it is inconvenient in producing and deploying sensor nodes. Sink mobility was also introduced to tackle the energy hole problem [1820]. The authors in [18] proposed an autonomous sink movement strategy in which mobile sinks move toward half-quadrant zones with abundant energy. Similarly, the authors in [19] drew a conclusion that utilizing mobile sink could prolong network lifetime 3.48 times at most compared to the case with static sinks (when the mobile sink moves around a hexagonal network perimeter and stops at the six corners). They also proposed a sink movement algorithm to keep network connectivity so as to prolong network lifetime. The authors in [20] tried to incorporate both static and mobile sinks to improve the network performance of WSN, in which the static sinks are located at the center of the target area and the mobile sink moves around the network perimeter. Each time when the mobile sink stops at a new location, it only needs to broadcast the location updates messages to a subset of sensor nodes. Despite the fact that mobile sinks would bring some advantages to WSN, some new issues appeared such as the movement control policy, the energy-aware routing protocols, and the efficient path planning. Unfortunately, these issues are very complicated and not readily solvable [21, 22].

Energy-aware sensor redistribution was also proposed to mitigate or avoid the energy hole problem in WSNs. In [23], the authors proposed an efficient node placement, topology control, and a scheduling protocol of the MAC layer to prolong the network lifetime for a grid-based WSN. In [24], the authors developed a mixed-integer linear programming model for determining the locations of sensors and sinks, schedules, and sensor-to-sink data flow routes. Further, they proposed a heuristic algorithm to maximize the network lifetime. In [25], the authors tried to solve the energy provisioning and relay node placement problems simultaneously in a two-tiered WSN. At first, the low-cost sensors sense the surrounding area and forward the sensing data to their cluster head, and then the cluster head forwards the data to the sink. Further, they proposed a heuristic approach to solve the mixed-integer nonlinear programming problem. In [26], the authors investigated sensor self-deployment problem, by constructing focused coverage around a point of interest. The authors in [27] explored variable node distribution density to mitigate the effects of the uneven energy depletion. The authors in [28] proposed a nonuniform distribution algorithm to solve the energy hole problem in a corona-based WSN. With their theoretical analysis, when all the sensors have a constant data acquisition rate, the balanced energy depletion among the whole network is impossible. However, their traffic pattern may not be true for highly dense WSN and the uneven energy depletion still exists between the outermost corona and the inner coronas. In fact, we can prove that balanced energy depletion is achievable with the extra help of energy-aware transmission mechanism in this paper. In [29], the authors investigated the sink-hole problem in duty-cycled connected -covered WSNs, where each point is covered by at least sensor. In [30], the ring-based relays are also used to eliminate energy hole of a WSN; however, how to form such a relay is not discussed in their approach. In our pervious approach [14], we have proposed the centralized algorithm to form an energy-balanced distribution for WSN; however, the exchange of global location information during dynamic sensor movement would put a heavy traffic burden on the network. In [31], the authors proposed a transmission range adjustment approach to tackle the unbalanced energy depletion. However, searching for the optimal transmission ranges among all the coronas is a NP-complete problem.

Overall, though most of the algorithms discussed above intended to maximize coverage rate, minimize deployment density, and eliminate the energy hole, they did not answer a fundamental question in sensor redeployment: what type of node layout and communication pattern that could provide the maximum coverage with the smallest overlap and gap and guarantee that all the working sensors die simultaneously with nearly zero residual energy? We will deal with this issue in the next section.

3. Network Model and Assumptions

In this section, we will present our network model and basic assumptions. Assume that a set of homogeneous sensors with the same initial energy is deployed in a circular area with radius to monitor some physical phenomena and an unlimited amount of energy is set for the sink node. We refer to the set of deployed sensors as , each of which has an ID, a fixed transmission range , a fixed sensing range , and is aware of its location. The only sink node is placed at the center of the circle. We divide the area into adjacent coronas with the equivalent width of , and the th corona is denoted by . Obviously, the corona is composed of nodes whose distances to the sink node are between and .

In this paper, periodic data gathering monitoring is considered, where the network is working in rounds. Each round is further divided into two phases: the node redistribution phase and stability monitoring phase. We will provide a detailed description of the first stage issues in the next section. During the second phase, each working node should send their sensing messages to the sink node per unit time via multihop communication. For the theoretical analysis, we use a simplified power consumption model and an ideal MAC layer with no collisions and retransmissions. The initial energy of each sensor is set as , and the sink has no energy limitation. We further assume that each sensor consumes units of energy when sending one bit, while it depletes units of energy when receiving one bit, where .

When using traditional transmission mechanism, the redundant sensing messages for the same area will be retransmitted by more than one sensor, consuming a considerable amount of energy [26]. In order to save energy, an energy aware data transmission mechanism proposed in our pervious approach [14] is imported. In this mechanism, we first need to build Voronoi graph for each sensor node to determine its own sensing range. After the establishment of the Voronoi graph, the sensing message for each pixel is sent only once. Based on the network model, we can conclude that sensors belonging to corona will forward both the data generated by themselves and the data produced by coronas . Obviously, the sensors located in the outermost corona do not need to forward any data. Assume that the sensors in each corona are distributed uniformly and there is no data aggregation at any forwarding nodes. Define the number of sensors deployed in corona as and the number of pixels in as . When the energy-aware data transmission mechanism is applied, the number of messages for corona to receive and forward is and . As sensing messages are transmitted per working round, the average energy consumption per working round of sensors in corona is where is the node density of corona .

As the sensors in corona only need to send their own sensing messages, the energy depletion of sensors in corona can be calculated as follows:

Definition 1. The optimal energy consumption state means that all the sensors in the network deplete their energy in the constant ratio; namely, all the sensor nodes have the same lifetime which is identical with the corresponding network lifetime. In particular, if the optimal energy depletion is achieved, there is no energy wasted and the network lifetime can be given by

Theorem 2. Optimal energy depletion is possible, if all the working nodes take the energy-aware transmission mechanism, and the node density in corona satisfies [14]

Proof. Without loss of generality, suppose that (4) is true: so (1) can be described as follows: As , we have
Since is a permanent establishment, we can draw a conclusion that . This completes the proof of Theorem 2.

Theorem 2 shows that, in a circular monitored area, based on the energy-aware transmission mechanism, if the sensors in each corona obey a uniform distribution and the node density meets a certain condition, the optimal energy depletion of the whole network can be achieved. Besides, the node density of corona only depends on and the corona number .

Theorem 3. The lifetime comparison of energy–balanced redeployment with traditional uniform approaches is .

Proof. Suppose that these two schemes run with the same initial conditions. The node density of each corona obeys (4) in our proposed scheme, while, in the uniform distribution, the density is equal to . As the innermost corona needs to forward all the sensing messages in the whole network, it consumes the most energy. Therefore, the maximum network lifetime under uniform distribution is determined by the survival time . The network lifetime can be calculated as where is the average energy depletion of per unit time under uniform distribution. Further, we can get the average energy depletion in under optimal energy consumption conditions as Thus the lifetime enhancement is
This completes the proof of Theorem 3.

Therefore, compared with traditional uniform distribution strategy, the network lifetime can raise as much as times by using optimal energy depletion distribution.

4. Optimization of Sensor Redeployment

In this section, we first introduce equivalent sensing radius. Then we develop the uniform sensor distribution for corona-shaped area. Further, we propose a novel autonomous sensor redeployment approach to balance energy depletion.

4.1. Redistribution under Equivalent Sensing Radius

Definition 4. Equivalent sensing radius is defined as the sensing radius when the given distribution density is the lowest one to maintain the full coverage for the target area with size .
According to [9], the equivalent sensing radius and the distribution density satisfy the following:

Theorem 5. Optimal energy depletion can be achieved, where satisfies the following:

Proof. According to the definition of equivalent sensing radius, we can combine it with the energy balance condition. Thus we have After simple transformation, we can get This concludes the proof of Theorem 5.

Since the equivalent sensing radius is only determined by corona number , the redeployment of sensors in corona is similar to the traditional uniform distribution problem. Define the equivalent sensing radius of corona as and the desired number of deployed sensor for as . Note that the boundary effects cannot be ignored simply. In order to cover corona uniformly, the optimal sensor distribution should satisfy the following conditions.(1)If , the optimal deployment can be achieved if all the sensors are uniformly lying on the middle line of corona , and the angle between any two adjacent nodes is .(2)If , as the single ring distribution will cause coverage gaps, the optimal sensor redeployment is equivalent to a multiring distribution. Define the number of rings for as , calculated as . When multiring uniform distribution is achieved for corona , the number of sensors distributed in ring is , and the angle between any two neighbours in is , where is the area of ring . Assuming that all the rings have the same width, , such that can be calculated as

Figure 1 shows the optimal sensor redeployment with different equivalent sensing radii . When , the sensor nodes in obey a single-ring uniform distribution, and the angle between any two neighbours is , as shown in Figure 1(a). When , the sensor nodes in corona obey a two-ring uniform distribution; the adjacent angles of ring and are and , respectively.

Therefore, the optimal sensor redeployment to balance energy depletion can be transformed into a uniform distribution problem with given deployment density and equivalent sensing radius. The novel sensor redeployment algorithm mainly contains two parts: (1) sensor redeployment control among the coronas to regulate the number of nodes for each corona and (2) sensor redeployment control inside coronas to guarantee that each corona achieves the given node density.

4.2. Distributed Sensor Redeployment among Coronas

As sensors are randomly deployed in the target area, the deployment uncertainty may cause the number of deployed nodes to be more or fewer than the corona really needs. Movement control of sensor nodes will satisfy the desired deployment density for each corona and its rings. To avoid consuming much energy during the moving process, the nodes are only allowed to move to the adjacent coronas. After the needed number of sensors is achieved in each corona, a certain number of sensors should move between different rings to achieve uniform distribution. Define the number of sensors deployed in corona as and the desired number of sensors in as . Sensor redeployment control among coronas can be conducted in Algorithm 1.

For each corona ,
   If
 Select nodes close to from ;
 Move these nodes straight to and update ;
 For each ring ,
   Move nodes straight to the middle line of
  Elseif
    Select nodes close to from ;
    Move these nodes straight to and update ;
    For each ring ,
    Move nodes straight to the middle line of .

4.3. Distributed Sensor Redeployment inside Coronas

Sensor redeployment inside coronas mainly focuses on how to redistribute the nodes locating on the median line of these rings in each corona to a perfect layout. Here, it is executed simply with its neighbour nodes. The straightforward idea is to adjust the angle between any two neighbours to . However, as each node may have numerous neighbours, the nearest neighbour is taken in the current working round. Taking Figure 2 as an example, assume that node is the nearest neighbour of node , to form an optimal layout and node should move to or , both of which are the best locations of node . In our algorithm, node moves straight to the nearer location , thus saving energy. If some other sensors existed in this best location, to minimize the total movement among all the sensors, should move with a small deviation along the same orientation to node . The whole process of sensor redeployment inside coronas is shown in Algorithm 2.

For each round
 For each node in ring of corona
Find its nearest neighbor ;
Calculate the angle between the line and ;
Calculate node m’s best location and
  Let ;
  
If distance ( ) < distance ( )
  Move node k to location ( );
  While there exists node with coordinate ( )
  
    Move node to location ( );
Else move node to location ( );
  While there exists node with coordinate ( )
   
    Move node to location ( );

Theorem 6. The sensor redeployment algorithm ends in finite rounds and can achieve uniform layout in the ring.

Proof. As long as the angle between any two nearest neighbours is not equal to , the total covered area with uniform distribution will keep on growing. It suffices to prove that, once the algorithm has reached the stable state when the angle between any two neighbours is equal to , the sensors locating on the ring will be fixed.
Let denote the set of sensors distributed in the ring, whereas is azimuth between the lines and . After reordering these sensors according to , we rename as . Define as the angle of two adjacent nodes in . To ensure that does not exceed , it needs to satisfy the following:
Thus the total covered angle of can be calculated as
If all the sensor nodes are distributed in an optimal layout, we can get the upper bound for as .
Denote the angle between and as . We now consider the changes of in one working round. If , the following three cases may occur.(1)If increases to with no sensors coinciding, the new achieved distribution can be expressed as , as shown in Figure 3(a). Since , we have .(2) If increases to with some sensors coinciding, the coincided sensors would move with a small deviation in our algorithm. Assume that the new achieved distribution is , as shown in Figure 3(b). Therefore, can be calculated as while Thus we have .(3)If increases to with some sensors crossing, the new achieved distribution can be denoted as , as shown in Figure 3(c). Since this new generated distribution is similar to case (2), we can draw the same conclusion of , when becomes the next operating node.
Clearly, in all the cases, is an increasing function on rounds . After a finite number of rounds, would achieve its upper bound. As a result, all the sensors will be redeployed uniformly on the middle line of the ring with the given density. This completes the proof of Theorem 6.

5. Performance Evaluation

In this section, we will present the simulation results of our algorithm for both random and Gaussian deployment models. Three metrics, including coverage rate, the average moving distance, and the network survival lifetime, are imported to evaluate the performance of the algorithm.

First, a random deployment is considered in the circular area of radius 100 with 627 potential sensors. We assume that all the sensors are homogeneous, that is, having the same initial energy reserve of 10000 J, the same transmission radius of 25 m, and the same sensing radius of 9 m, respectively. The values of and are given as follows: /103 J/bit and /103 J/bit. The total number of working rounds for sensor node is determined by , where is the total number of pixels covered by is the total number of messages forwarded by , and is the length of each sensing pixel,  bits.

The sensing data forwarding strategy is similar to [27]. As the nodes obey an approximate uniform distribution in the corona after sensor redeployment, any node in corona can directly communicate with almost nodes in . Among these candidate nodes, the node with most residual energy will be selected as the forwarding one.

Figure 4 depicts the sensor distribution gained in different rounds using random deployment model. The target area is divided into four coronas marked from to , and the equivalent sensing radius from to can be calculated as 1.86, 3.4, 5.11, and 9, respectively, as shown in Figure 4(a). From the simulation results, we can make the observation that our proposed algorithm can converge to the global optimal distribution effectively. What is more, the layout shown in Figure 4(d) is the optimal uniform sensor distribution with the given density, which achieves the coverage rate of 97.6%. And the number of sensors redeployed from corona to is 220, 196, 145, and 66, which approaches to the accessibility condition of optimal energy depletion.

The second assumption examined is the Gaussian distribution, and the simulation environment is the same as the first experiment. In this distribution, each sensor's coordinate obeys a two-dimensional normal Gaussian function: where represents the standard deviation of Gaussian distribution. We set equal to the communication radius to ensure that all the sensors will be deployed in the target area.

Figure 5(a) shows the sensor distribution obtained in different rounds using Gaussian deployment model. From the simulation results, we can draw a conclusion that the total coverage enlarges round by round, as observed from Figure 5(b) to Figure 5(d). At the 28th round, the nodes are uniformly distributed with no sensors moving anymore, which indicates that the algorithm can converge to the optimal solution quickly.

Further, we compare the performance of our algorithm with VEC in terms of coverage rate and average moving. Figure 6 shows the average moving distance and coverage rate for VEC [10] and our approach in different numbers of rounds. From the simulation results in Figure 6, we can conclude that (1) the average moving distance under Gaussian deployment model is shorter than under random model, which is mainly due to having more sensors deployed in the central region under Gaussian deployment model; (2) to achieve the desired coverage rate, many more sensors are needed in VEC than in our proposed algorithm, which shows our algorithm has a shorter average moving distance than VEC under the same condition; (3) our algorithm can get a higher coverage rate than VEC, which means that our algorithm is free of the boundary influence and can achieve the optimal sensor distribution effectively.

We further compare the performance of our algorithm with traditional uniform and nonuniform sensor deployment in terms of network lifetime and energy efficiency. In addition, the performance of our energy-aware transmission mechanism () and traditional transmission mechanism () combined these algorithms is also tested.

Since the node densities for all the coronas are the same in uniform deployment, the number of sensors distributed from to is 274, 196, 118, and 39. While fewer nodes are used to cover the outmost corona in the nonuniform deployment, the ratio is set to 2.12. Therefore, the number of sensors uniformly distributed from corona to is 66, 73, 156, and 332, respectively. For our algorithm, to ensure the optimal energy consumption, the number of sensors distributed from corona to is 66, 145, 196, and 220.

Figure 7 shows the comparisons of average energy depletion per working round. From the simulation results in Figure 7, we can draw the conclusions that (1) the average energy consumed in a round in is much smaller than in . Because the energy-aware transmission can avoid retransmitting the same sensing data, the total messages transferred in each round are much smaller; (2) since the uniform deployment ignores the traffic imbalance, the innermost corona has the largest traffic load and most energy consumption, as shown in Figure 7(a); (3) the nonuniform deployment combined with scheme can achieve a suboptimal distribution; (4) the energy consumption of the whole working set is almost equal in our approach. Although the nodes in inner coronas do not behave as source but router, their sensing pixels are much smaller than those of the outer coronas.

Figure 8 compares the energy unused ratio of each node when the network terminates. Here, the energy unused ratio refers to the ratio of the residual energy to the initial energy at the end of the network lifetime. From Figure 8, clearly, the uniform scheme leaves an abundant amount of nodal energy in both transmission mechanisms as expected, especially for the sensor nodes in the outermost corona. Although the nonuniform deployment combined with scheme achieves much better energy efficiency, it leaves a considerable amount of unused energy in the outermost corona (almost 50% or so of the initial nodal energy). For the nonuniform deployment combined with , as corona has the most transmission load and corona has the largest node density, no doubt, has the largest energy unused ratio. For our algorithm, as the energy consumption of the whole working set is almost equal, all the sensors exhaust their energy simultaneously when the network terminates, as shown in Figure 7. In fact, the energy unused ratio for most sensors is lower than 1%.

Table 1 compares these algorithms in network lifetime. We can see that the algorithm combined with energy-aware transmission mechanism has a longer network lifetime as expected. Because of the imbalance consumption of energy near the sink, the uniform distribution has the shortest network lifetime. The nonuniform distribution combined with achieved a shorter network lifetime than with , because its average energy consumption is larger than in scheme. Our algorithm can achieve the optimal energy consumption among the whole working nodes, thus having the longest network lifetime.

6. Conclusion and Future Work

In this paper, we have investigated the problem of sensor redeployment to achieve optimal energy depletion and minimize sensor movement. We have given a theoretical analysis on energy consumption using nonuniform distribution strategy. Formally, we have proved that, the optimal energy consumption can be achieved through calculating the node densities for different regions of the target area. As a contribution, we have proposed an autonomous coverage-driven sensor redeployment algorithm to produce an optimal solution, which can maximize the network lifetime and minimize total movement of sensors. In addition, extensive simulation results have been presented to demonstrate the effectiveness of our proposed techniques.It also should be noted that we only consider the two dimensional (2D) case in this paper. As part of our future work, we will design new algorithms for the sensor redeployment problem in 3D space.

Conflict of Interests

The authors declare that they have no conflict of interests regarding the publication of this paper.

Acknowledgments

This work is supported by the National Natural Science Foundation of China under Grant no. 61173153 and no. 60903159, the National Science Foundation for Distinguished Young Scholars of China under Grant no. 61225012 and no. 71325002, the Fundamental Research Funds for the Central Universities under Grant nos. N130504007, N110318001, N110204003, N100218001, and N120104001, the China Postdoctoral Science Foundation funded project under Grant no. 20110491508 and no. 2012T50248, the Specialized Research Fund of the Doctoral Program of Higher Education for the Priority Development Areas under Grant no. 20120042130003, and the Specialized Research Fund for the Doctoral Program of Higher Education under Grant no. 20110042110024.