- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents

International Journal of Distributed Sensor Networks

Volume 2013 (2013), Article ID 483581, 12 pages

http://dx.doi.org/10.1155/2013/483581

## A Novel *Physarum*-Inspired Routing Protocol for Wireless Sensor Networks

^{1}State Key Laboratory of Networking and Switching Technology, Beijing University of Posts and Telecommunications, Beijing 100876, China^{2}Information Engineering College, Henan University of Science and Technology, Luoyang 471023, China^{3}Institute of Sensing Technology and Business, Beijing University of Posts and Telecommunications, Wuxi 214135, China^{4}National Engineering Laboratory for Next Generation Internet Interconnection Devices, Beijing Jiaotong University, Beijing 100876, China

Received 31 January 2013; Accepted 7 June 2013

Academic Editor: Shuai Li

Copyright © 2013 Mingchuan Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

There is a tradeoff between routing efficiency and
energy equilibrium for sensor nodes in wireless sensor networks
(WSNs). Inspired by the large and single-celled amoeboid
organism, *slime mold Physarum polycephalum*, this paper
presents a novel *Physarum*-inspired routing protocol (P-iRP) for
WSNs to address the above issue. In P-iRP, a sensor node can choose
the proper next hop by using a proposed *Physarum*-inspired
selecting next hop model (P-iSNH), which comprehensively considers
the distance, energy residue, and location of the next hop.
As a result, the P-iRP can get a rather low algorithm complexity
of , which greatly reduces the processing delay and saves
the energy of sensors. Moreover, by theoretical analysis, the
P-iSNH always has an equilibrium solution for multiple next
hop candidates, which is vital factor to the stability of routing
protocol. Finally, simulation results show that P-iRP can perform
better in many scenarios and achieve the effective tradeoff
between routing efficiency and energy equilibrium compared to
other famous algorithms.

#### 1. Introduction

With the development of communication, electron, and sensor technologies, wireless sensor networks (WSNs) have attracted wide concern of both researchers and application providers. WSNs consist of large numbers of sensor nodes deployed over a certain region. Each sensor node is a low-cost, short range wireless transceiver typically equipped with a low-computation processor and a battery operated power supply. Under many scenarios, the sensor nodes need to operate without battery replacement for several years. Thus, there are two questions need to be considered. One is how to achieve energy balance of these nodes to avoid the emergence of energy holes which commonly take place around the sink, since the data traffic follows a many-to-one communication pattern and nodes nearer the sink have to take heavier traffic load. The other is how to obtain high routing efficiency under multihop transmission circumstance, since WSNs can contain hundreds of such low-cost sensor nodes. Therefore, designing such networks should primarily focus on both routing efficiency and energy equilibrium in terms of trade-off.

Location-aware routing protocols seem to possess high routing efficiency, where GPS, phone, or other techniques are used for positioning nodes [1]. However, there are two extremes in location-aware routing—the greedy strategy and the robust strategy. Greedy strategies may suffer failures to route packets to destination, while robust strategies need very high flooding rates to ensure reliability and rapid delivery of data. Thus, many location-aware routing protocols are mostly to propose methods to overcome the mentioned drawbacks [2–4]. The energy-aware routing attracts more researchers’ attentions than that of location-aware routing for the significance of energy. There are many results relating to energy-aware routing recently [5–8] to save energy or prolong WSNs’ lifetime, where energy harvesting [9] is shown to be a promising technique.

Unsurprisingly, the combination between location-aware routing and energy-aware routing becomes another researchers’ focus to balance energy and efficiency for WSNs’ routing protocol [10, 11]. In addition, some researchers focus on other aspects of WSNs’ routing protocol, for example, the distributed characteristic [12, 13] and the trade-off of other two or more indexes [14, 15].

In recent years, the *slime mold Physarum polycephalum*, a large single-celled amoeboid organism, becomes a new researchers’ pet and has shown to be a good technique for solving the shortest-path problem, since it can adapt its organism to forage for patchily distributed food sources, as shown in Figure 1. In this paper, we draw the inspiration from the *Physarum*, introduce the *Physarum* model into WSNs, and improve it through ignoring its dimension and preserving its logical meanings to make it suit for routing selection based on our prior works [16–20]. Our focus is to choose the proper next hops to transmit data to sink in thinking of both routing efficiency and energy equilibrium, which is partially similar to [11].

The rest of this paper is organized as follows. Section 2 gives a brief description of related work. Section 3 formulates the proposed models. Section 4 details the P-iRP. Section 5 discusses the feasibility of P-iSNH. Section 6 evaluates our P-iRP by simulations. Finally, the conclusion is presented in Section 7.

#### 2. Related Work

In the aspect of efficient routing, GPSR [2] is a famous greedy routing protocol, which makes greedy forwarding decisions using only information about a router’s immediate neighbors in the network topology. Li et al. [21] present a neural network approach to plan the shortest path from the target position to the start position in real time. Kuhn et al*.* [3] utilize face (or perimeter) routing to go around voids in the topology. Padmanabh et al. [4] consider unbiased random walk on a regular deployment of nodes, forming a hexagonal lattice pattern.

In the aspect of the energy-aware routing, Trajcevski et al. [5] construct a data aggregation tree that minimizes the total energy cost of data transmission, which is shown as an NP-complete problem, and propose algorithms for addressing it. A battery aware power allocation model was studied in [6] for a single-hop transmission scheme to balance the network energy consumption based on the nonlinear battery parameters proposed in [7]. Chau et al. [8] consider that a portion of the lost charge can be recovered due to the battery’s recovery effect and present a battery model. The approaches propose some of the routes that would otherwise need to bypass the hole along the boundary and should start to deviate from their original path further from the hole instead.

Moreover, distribution and clustering problems are also important branches of WSNs’ routing. Li et al. [13] consider the problem of nonlinear constraints defined on a graph and give a better solution incorporated with Laplacian eigenmap as heuristic information to solve the problem in distributed scenarios. For maximizing the network lifetime, Rao and Fapojuwo [22] present a battery aware distributed clustering and routing protocol which incorporates the state of the battery’s remaining charge and health parameters in computing the charge utility metric at each cluster formation round. Wang and Syue [23] propose a relay selection protocol based on geographical information, in which multihop transmission is realized by concatenation of single cluster-to-cluster hops, where each cluster-to-cluster scheme forms the simplified cooperative network that consists of a single source destination pair and a set of available relays.

However, the trade-off is not comprehensively considered in those papers, which is very necessary for WSNs’ routing due to the features of WSNs, for example, nodes’ failures, limited bandwidth, and power energy. Bai et al. [10] route the connections in a manner that link failure does not shut down the entire stream but allows a continuing flow for a significant portion of the traffic along multiple paths to address the issues of reliability and energy efficiency. Trajcevski et al. [24] present heuristic approaches to relieve some of the routing load of the boundary nodes of energy holes in location-aware WSNs to balance load and latency. Yu et al. [11] use energy aware neighbor selection to route a packet towards the target region and recursive geographic forwarding or restricted flooding algorithm to disseminate the packet inside the destination region. By allowing the battery to rest for certain duration, without being subjected to heavy loads, Yang and Heinzelman [14] propose sleeping multipath routing, which selects the minimum number of disjoint paths to achieve the trade-off of given reliability requirement and energy efficiency. Sivrikaya et al*.* [15] propose randomized routing based on Markov chains to balance the load and routing performance.

In recent years, the *Physarum* becomes a new focus of bioinspired method. It is also the original source of our inspiration. Nakagaki et al. [25] validate that the *Physarum* is apparently able to solve shortest path problems as shown in Figure 1. They build a maze, cover it with pieces of *Physarum* (the *Physarum* can be cut into pieces that will reunite if brought into vicinity), and then feed the *Physarum* with oatmeal at two locations. A few hours later, the *Physarum* retracts to a path that follows the shortest path connecting the food sources in the maze. Tero et al. [26] use *Physarum* forms networks comparable efficiency, fault tolerance, and cost to those of real-world infrastructure networks—Tokyo rail system. Tero et al. [27] propose a mathematical model for the behavior of *Physarum* and argue extensively that the model is adequate.

#### 3. System Models

Our research is built on three assumptions. The first is that all nodes are aware of their locations, which may be achieved through GPS receivers at network deployment time, employing a distributed location discover algorithm shortly after deployment or adopting other positioning methods [1, 28]. The second is that each node is aware of its energy residue [22]. The third is that the link is bidirectional, that is, if a node hears from a neighbor, then its transmission range can reach to the neighbor.

##### 3.1. Typical WSNs Scenario

We consider the large multi-hop WSNs which consist of static sensors. Each node has a fixed circular transmission range which determines the set of sensors in which each node can communicate with node in one hop. We abstract such WSNs using a graph , where each node represents a sensor, and each edge represents the existence of one-hop wireless link between two sensors.

We suppose that node is the source node and node is the sink, as shown in Figure 2. In most cases, the sink is placed in the middle of WSNs field to ease traffic burdens of nodes in the right of the WSNs field. In this paper, we only think of the nodes in left of the WSNs field for simplicity and clarity. The transmission range of is drawn as a dashed circle whose radius is and center is . We call the angle is the angle of deviation of node , which represents a measurement of node deviating from the sink . The Euclidean distance of any two nodes, and , and the angle can be calculated following from (1) and (2), respectively where and are the coordinates of nodes and , respectively.

If the node needs to transmit data to sink , it will select its next hop in the dashed circle. Since our scenario is location aware, we select the next hop in the right semicircle under normal circumstances. Obviously, the smaller the angle is, the closer the next hop is to the sink for a fixed distance. That is to say, we are apt to choose the node whose is smaller as the next hop. In order to avail discussion, we define the is the set of neighbors of node (in the dashed circle), is the set of left neighbors of node (in the left semicircle), and is the set of right neighbors of node (in the right semicircle). Then, , , and meet the following:

In addition, the energy residue of each node is also important for choosing next hop for balancing the energy of WSNs’ nodes. When a node chooses its next hop, it would consider the energy residue of the candidates and be apt to pick the node which has much higher energy residue as the next hop. Therefore, it is important to acquire the energy residue of neighbors.

We think of the basic theory of wireless transmission combined with Figure 2 and the data packet format shown in Figure 3(a). If node transmits a group of data to node , all of the nodes in would receive the wireless radio and check the packet header. The node matches the field NA and receives the packet. Other nodes mismatch the field NA then ignore the packet and go on sleeping.

In order to acquire the energy residue of neighbors, we add a new field ER to the packet header, which is shown in Figure 3(b). When node transmits a group of data to node , all of the nodes in extract the fields of SA and ER from the packet header and save ER in local memory according to SA. Then, the node matches the field NA and receives the packet. Other nodes mismatch the NA then ignore the packet and go on sleeping.

Since each node needs to listen in real time to every packet and try to match its field NA, only adding an operation of saving ER would not add a considerable effect on energy consumption. Therefore, we neglect the cost of acquiring energy residue of neighbors.

##### 3.2. *Physarum*-Inspired Path-Finding Model

Papers in [25–27] exploit the *slime mold Physarum polycephalum* to develop a *Physarum*-inspired path-finding model (PiPf). Suppose that (1) the initial shape of a *Physarum* organism is represented by a graph, (2) the edges represent plasmodial tubes in which protoplasm flows, and nodes are junctions between tubes, (3) the pressures at nodes and are and , respectively, and the two nodes are connected by a cylinder of length and radius , and (4) the flow is laminar and follows the Hagen-Poiseuille equation. Then, the flux through the tube is calculated as in the following
where is the difference of pressures, is the viscosity of the fluid, and is a measure of the conductivity of the tube. As the length is a constant, the behavior of *Physarum* is described by the conductivities, , of the edges.

Equation (4) represents that the flux through the tube is determined by , , and . The better the conductivity of the tube is and the larger the pressure difference is, the more the flux through the tube is, while the longer the length of the tube is, the less the flux through the tube is.

Suppose that the capacity of each node is zero, the conservation law of each node is calculated from the following:
where is the flux flowing from the source node (or into the sink node). It should be noted that is a constant in *Physarum* model, which means that the total flux is fixed constant throughout the process.

Equation (5) illuminates the flux relationship in each node. In the source node , is the flux flowing from it; in the sink node , is the flux flowing into it; and in intermediate nodes, the sum of flowing from and flowing into is zero.

*Physarum* forages for distributed food sources through the adaptive behavior of the plasmodium. The adaptive behavior is illustrated as follows combined with Figure 4(a), where two food sources are connected by two tubes. Because of and , the flux will be greater than from (4). Note that and are kept constant throughout the adaptation process in contrast to ; therefore, the adaptive behavior can be described by the evolution of ,
where is a decay rate of the tube. Equation (6) implies that the conductivity tends to vanish if there is no flux along the edge, while it is enhanced by the flux. It is natural to assume that is a monotonically increasing continuous function satisfying .

Equation (6) illustrates the variation relationship of the conductivity with time to accommodate the flux distribution of the multipath transmission. In equilibrium ( for all edges), the flow through any edge is steady. In nonequilibrium, the diameter grows or shrinks if is larger or smaller than , respectively.

The PiPf consisting of (4), (5), and (6) describes the evolutionary process of *Physarum* to solve the path-finding behavior of self-organized networks.

##### 3.3. *Physarum*-Inspired Selecting Next Hop Model

In this section, we improve the PiPf and make it fit for routing in WSNs based on dimensionless analysis method. That is to say, we improve the PiPf to achieve a *Physarum*-inspired selecting next hop model (P-iSNH). In order to obtain that, there are two problems that need to be solved. The first is which physical quantities are used to replace the conductivity , the length , and the pressure difference . The second is how to select the proper next hop.

Equation (4) derives from fluid dynamics. is a measure of the conductivity of the tube; is the length of the tube; and is the differential pressure of tube on both ends. However, the , , and cannot be directly used in WSNs where we need to consider the link quality, energy residue, transmission direction, and the distance of one hop.

First, because the is an inherent characteristic of the tube, we should replace the by an inherent physical quantity. Apparently, the link quality is an inherent characteristic relating to wireless link, so we replace the by in our model.

Second, the meaning of is the same as in fluid dynamics. However, in wireless communication, there is a path-loss exponent , which has a great effect on transmission. Therefore, we replace the by in our model.

Third, we discuss the combind with Figure 4(b). On one hand, suppose that there is a potential field from node to sink . The potential difference between and can be expressed by . Because parameter is unimportant in the judge process, we use expressing . When node chooses its next hop, it is apt to pick the node whose potential difference is much greater. On the other hand, because the next hop needs to consume energy to deal with data packets, it is apt to pick the node with much higher energy residue as the next hop. Therefore, we replace the by . Since is the base pressure, we replace by through omitting . Using (4), we have where is the virtual communication fluxes through the wireless link ; is the link quality; is the energy residue of node ; is the Euclidean distance of nodes and ; is path-loss exponent; is the angle of deviation and its range is ; is a proportionality factor which uses to adjust the weight of and .

Then, we discuss how to choose the proper next hop. As related in Section 3.2, since and and are kept constant throughout the adaptation process in contrast to , the PiPf can only achieve the adaptation by the evolution of . In our scenario, node chooses the next hop form candidates as shown in Figure 4(b). Since (1) and and are kept constant according to the assumptions and (2) and are different and time-varying, we can achieve the adaptation by the evolution of . If letting the monotonically increasing continuous function , we have where is a decay rate of and is a constant satisfying . We use (8) to determine the next hop in our P-iSNH; namely, we choose the node whose is maximal as the next hop.

#### 4. P-iSNH Based Routing Strategy and Algorithm

##### 4.1. Data Conserved

In this section, we introduce the data which should be conserved in each node. Because of the same characteristic of each node, we suppose that the of each link is the same and ignore it to simplify discussion.

Each node needs to conserve the following information: , , , , and which are shown in Figure 5, where node is the previous hop of node , represents the stored in node , and represents the stored in node . As our WSNs are location aware, the , , and are easily acquired following from (1) and (2). Note that the nodes in our WSNs are static, and we only need to calculate the , and once at WSNs deployment time. For the difference of and , we normalize them to and , respectively. Therefore, we obtain

##### 4.2. Routing Strategy

If node needs to send data to the sink , it searches for a routing in the following method. We illustrate the routing strategy combined with Figure 6.

*Step 1. *Each is calculated following from (10), where , , and are conserved and stored in node beforehand.

*Step 2. *Each node is saved into a temporary array variable in descending order by .

*Step 3. *The first node in is picked as the next hop of the routing.

*Step 4. *If the next hop of node satisfies , namely, there is an energy hole in the right side of node , the node will not send to . Then, the node will trigger a specific processing routine.

*Step 5. *If , node will choose the node as the next hop. Then, the regular processing routine is going on.

*Step 6. *Otherwise, each is calculated following from (11) and the nodes are saved into the in ascending order by . Then, the first node in is chosen as the next hop of the routing and the regular processing routine is going on
where is the angle of line and line . Equation (11) indicates that it tends to choose a node which sharply deviates from the failing node, for example, , as the next hop to avoid entering the energy hole again.

*Step 7. *The process is repeated, like a rolling wheel, until the sink is found.

##### 4.3. Routing Algorithms

Given the data conserved and routing process in the preceding sections, the P-iRP’s algorithms of initialization, regular processing routine, receiving routine, and specific processing routine are described by Algorithm 1, Algorithm 2, Algorithm 3, and Algorithm 4, respectively.

Based on the WSNs scenario in Section 3.1, suppose that the degree of graph is which can be regarded as a constant, the complexity of Algorithm 1, Algorithm 2, Algorithm 3, and Algorithm 4 are , , , and , respectively. From Figure 2, the number of intermediate nodes from node to sink is approximately in the worst case. Note that Algorithm 1 is run only once at WSNs deployment time. Therefore, the complexities of P-iRP is in running time, which greatly reduces the processing delay and saves the energy of sensors.

#### 5. P-iSNH Analysis

In this section, we analyze the feasibility of P-iSNH by mathematical theoretical analysis. We study the cases in which two nodes connected to the same node compete to be the next hop, as shown in Figure 4(b).

There are four nodes , , , and sink . For simplicity, we hereafter replace , , , , , and by , , , , , and , respectively. In multipath routing, the virtual fluxes along each path are calculated as

Since and are nonnegative, adaptation equation (8) becomes

Setting , , and , we have

After some calculations, we obtain

Namely, there is an equilibrium point given by . We perform the simulation using MATLAB by setting the parameters , , , and following from (14), and the solutions are shown in Figure 7, where two curves intersect in a point which superpose on the equilibrium point .

We present a linear stability analysis at the equilibrium point in before parameters. The Jacobi matrix on the right-hand side of (13) is calculated as where and the Jacobi matrix at equilibrium point is denoted . After some calculations, the following formula is obtained where and are virtual communication fluxes along the first and second wireless links at the equilibrium point . Using the relation , we have

Note that is the decay rate of and . If we let , thus,

Since we set , this means that the equilibrium point is stable. Therefore, the routing of WSNs will reach to equilibrium with our P-iSNH, which is very important to a routing strategy.

#### 6. Simulation Results

We design a simulation platform using C++ to validate P-iRP. In the simulation, 441 sensors are relatively regularly deployed in the field of 200 m × 200 m, and the sink node is deployed in the right of the field, shown in Figure 8. The sensing radius of each sensor is 30 m, the original energy of each node is 100, and the energy of sink node is inexhaustible. We suppose that the energy consumption of one transmission is 1, if the transmission distance is 20 m. Therefore, the energy consumption of one transmission of two nodes and is , where is set to 2.

In order to validate the energy equilibrium, we only choose the nodes in the field of to transmit data to the sink. If a chosen node transmits a group of data to sink, the P-iRP is used to choose next hops until the sink is found, which is called a *round*. This iterative process will halt after rounds until WSNs break down. We run GPSR, GEAR () and P-iRP () 10 times, respectively, to acquire their average value and compare them, where we use replace which is used in GEAR to bring into correspondence with P-iRP. If the distance between the nodes and sink is less than 30 m, we let the nodes directly transmit data to sink to quicken convergence of P-iRP, and the energy consumption is set to 1. From Figure 8, there are 18 sensors around the sink. Therefore, the ideal number of rounds of simulation process is .

##### 6.1. Energy Equilibrium of P-iRP

Figure 9 illustrates the energy distribution of GPSR, GEAR (), and P-iRP () in different rounds. We can infer that the energy distributions of GPSR are very imbalanced, the energy distributions of GEAR and P-iRP are rather balanced, and the energy distributions of P-iRP are more balanced than those of GEAR.

Figure 10 illustrates the lifetime of WSNs. In GPSR, the first dead node emerges in round of 192, and the WSNs break down in round of 889. In GEAR (), the first dead node emerges in round of 910, and the WSNs break down in round of 1223. In P-iRP (), the first dead node emerges in round of 1112, and the WSNs break down in round of 1396. In GEAR (), the first dead node emerges in round of 1428, and the WSNs break down in round of 1592. In P-iRP (), the first dead node emerges in round of 1542, and the WSNs break down in round of 1696. Therefore, the lifetime of GEAR () is 48.8% longer than that of GPSR; the lifetime of P-iRP () is 14.2% longer than that of GEAR (); and the lifetime of P-iRP () is 6.5% longer than that of GEAR (). From Figure 9 and Figure 10, we can differ that (1) whether considering energy residue of next hops or not will impacts on the lifetime of WSNs greatly, and (2) in energy balanced WSNs, the time period is very short from emerging dead nodes to networks breaking down because all nodes reach to exhausted status of energy in the same time period.

Figure 11 illustrates the dead nodes distributions of GEAR () and P-iRP () in the rounds of 1380. The results show that P-iRP () has much less dead nodes than GEAR (). We can also differ that the dead nodes of both algorithms are converged on a specific field but do not spread around the entire range of WSNs, which is useful in deploying such WSNs to prolong the lifetime.

The reasons to gain the results of Figures 9, 10, and 11 are that (1) since GPSR does not take energy into account, it utilizes frequently the “hot” nodes to result in imbalanced energy distributions, (2) since GEAR and P-iRP consider both energy and location of nodes, their energy distributions are rather balanced, and (3) since P-iRP is more elaborate in energy utilization than GEAR, the energy distributions of P-iRP are more balanced than those of GEAR.

##### 6.2. Efficiency of P-iRP

Figure 12 illustrates the number of hops that the different algorithms need in different rounds of transmission. By calculating, the average hops of GPSR, GEAR (), P-iRP (), GEAR (), and P-iRP () are 19.4, 24.3, 21.6, 28.9, and 26.8, respectively.

In case of , the average hops of P-iRP are 11.3% more than those of GPSR, and the hops of GEAR are 23.3% more than those of GPSR. Combined with Figure 10, the increment of average hops of 11.3% will lead to the increment of lifetime of more than 60% from GPSR to P-iRP, while the increment of average hops of 23.2% will only lead to the increment of lifetime of about 48% from GPSR to GEAR. Therefore, the P-iRP is more efficient in balance of routing efficiency and energy equilibrium than GEAR.

In case of , the average hops of P-iRP are 38.1% more than those of GPSR, and the hops of GEAR are 49.0% more than those of GPSR. Combined with Figure 10, the increment of 11.3% of average hops will lead to the increment of about 60% of lifetime from GPSR to P-iRP (), while the increment of 24% of average hops will only cause the increment of 21.5% of lifetime from P-iRP () to P-iRP (). That is to say, the larger is, the smaller the increment of impacts on lifetime of WSNs. Therefore, it is improper to set a larger , so does GEAR.

In addition, Figure 12 implies that (1) GPSR can get much higher routing efficiency than those of GEAR and P-iRP at initial periods of time in WSNs’ lifetime, but its efficiency decreases exponentially with the time going, and (2) the routing efficiencies of GEAR and P-iRP are almost the same at initial periods of time in WSNs’ lifetime, while the difference becomes gradually large with the time going. The reasons to gain those results are that (1) since GPSR does not take energy into account, the routing efficiency will decrease exponentially after dead nodes emerge and (2) P-iRP can outperform GEAR in routing efficiency, since P-iRP considers comprehensively the distance, energy residue, and location of the next hop, other than only considering energy and location in GEAR.

#### 7. Conclusion

The *Physarum* forages for patchily distributed food sources through accommodating its body to form networks with comparable efficiency, fault tolerance, and cost, which is the source of P-iRP’s inspiration. For the proposed scenario, the P-iRP ensures the passage of data packets through one by one static sensor nodes to reach the sink. In each intermediate node, the P-iSNH is used to choose the proper next hop. Once an energy hole emerges, a specific processing will be triggered to bypass the hole. The theoretical analysis and simulation results show that the P-iRP possesses many advantages, for example, rather low algorithm complexity for P-iRP, ever-present equilibrium solution for P-iSNH, and effective trade-off between routing efficiency and energy equilibrium compared to other famous algorithms, which greatly reduces the processing delay and saves the sensors’ energy and also demonstrates that the P-iRP is applicable to the proposed scenario. Furthermore, we consider that the model may also provide a useful help to develop the routing protocol in mobile ad hoc networks, which will be our future focus.

#### Acknowledgments

This work was partially supported by the National High-Tech Research and Development Program of China (863 Program) under Grant no. 2011AA010701, in part by the National Basic Research Program of China (973 Program) under Grant no. 2013CB329102, in part by the National Natural Science Foundation of China (NSFC) under Grant nos. 61003283, 61001122, 61232017, U1204614, 61003035, and 61142002, and in part by the Natural Science Foundation of Jiangsu Province under Grant no. BK2011171.

#### References

- S. Li, Y. Lou, and B. Liu, “Bluetooth aided mobile phone localization: a nonlinear neural circuit approach,”
*Transactions on Embedded Computing Systems*, 2013. - B. Karp and H. T. Kung, “GPSR: greedy perimeter stateless routing for wireless networks,” in
*Proceedings of the 6th Annual International Conference on Mobile Computing and Networking (MOBICOM '00)*, pp. 243–254, New York, NY, USA, August 2000. View at Scopus - F. Kuhn, R. Wattenhofer, and A. Zollinger, “An algorithmic approach to geographic routing in ad hoc and sensor networks,”
*IEEE/ACM Transactions on Networking*, vol. 16, no. 1, pp. 51–62, 2008. View at Publisher · View at Google Scholar · View at Scopus - K. Padmanabh, A. M. Reddy V., S. Sen, and P. Gupta, “Random Walk on Random Graph based Outlier detection in wireless sensor networks,” in
*Proceedings of the 3rd International Conference on Wireless Communication and Sensor Networks (WCSN '07)*, pp. 45–49, Allahabad, India, December 2007. View at Publisher · View at Google Scholar · View at Scopus - G. Trajcevski, F. Zhou, R. Tamassia, and B. Avii, “On the construction of data aggregation tree with minimum energy cost in wireless sensor networks: NP-completeness and approximation algorithms,” in
*Proceedings of the INFOCOM*, pp. 2591–2595, Orlando, Fla, USA, 2012. - M. H. Chaudhary and L. Vandendorpe, “Battery-aware power allocation for lifetime maximization of wireless sensor networks,” in
*Proceedings of the IEEE International Conference on Communications (ICC '10)*, Cape Town, South Africa, May 2010. View at Publisher · View at Google Scholar · View at Scopus - J. Zhang, S. Ci, H. Sharif, and M. Alahmad, “Lifetime optimization for wireless sensor networks using the nonlinear battery current effect,” in
*Proceedings of the IEEE International Conference on Communications (ICC '09)*, Dresden, Germany, June 2009. View at Publisher · View at Google Scholar · View at Scopus - C.-K. Chau, F. Qin, S. Sayed, M. H. Wahab, and Y. Yang, “Harnessing battery recovery effect in wireless sensor networks: experiments and analysis,”
*IEEE Journal on Selected Areas in Communications*, vol. 28, no. 7, pp. 1222–1232, 2010. View at Publisher · View at Google Scholar · View at Scopus - J. Yang and S. Ulukus, “Optimal packet scheduling in an energy harvesting communication system,”
*IEEE Transactions on Communications*, vol. 60, no. 1, pp. 220–230, 2012. View at Publisher · View at Google Scholar · View at Scopus - S. Bai, W. Zhang, G. Xue, J. Tang, and C. Wang, “DEAR: delay-bounded energy-constrained adaptive routing in wireless sensor networks,” in
*Proceedings of the INFOCOM*, pp. 1593–1601, Shanghai, China, 2012. - Y. Yu, R. Govindan, and D. Estrin, “Geographical and Energy Aware Routing: a recursive data dissemination protocol for wireless sensor networks,” Tech. Rep. UCLACSD TR-01-0023, UCLA Computer Science Department, 2001.
- N. Roseveare and B. Natarajan, “Energy-aware distributed tracking in wireless sensor networks,” in
*Proceedings of the IEEE Wireless Communications and Networking Conference (WCNC '11)*, pp. 363–368, Cancun, Mexico, March 2011. View at Publisher · View at Google Scholar · View at Scopus - S. Li, Z. Wang, and Y. Li, “Using Laplacian Eigenmap as heuristic information to solve nonlinear constraints defined on a graph and its application in distributed range-free localization of wireless sensor networks,”
*Neural Processing Letters*, vol. 37, no. 3, pp. 411–424, 2012. - O. Yang and W. Heinzelman, “Sleeping multipath routing: a trade-off between reliability and lifetime in wireless sensor networks,” in
*Proceedings of the 54th Annual IEEE Global Telecommunications Conference (GLOBECOM '11)*, Houston, Tex, USA, December 2011. View at Publisher · View at Google Scholar · View at Scopus - F. Sivrikaya, T. Geithner, C. Truong, M. A. Khan, and S. Albayrak, “Stochastic routing in wireless sensor networks,” in
*Proceedings of the IEEE International Conference on Communications Workshops (ICC '09)*, Berlin, Germany, June 2009. View at Publisher · View at Google Scholar · View at Scopus - M. Zhang, C. Xu, J. Guan, R. Zheng, Q. Wu, and H. Zhang, “P-iRP: physarum-inspired routing protocol for wireless sensor networks,” in
*Proceedings of the VTC*, Las Vegas, Nev, USA, 2013. - C. Xu, T. Liu, J. Guan, H. Zhang, and G.-M. Muntean, “CMT-QA: quality-aware adaptive concurrent multipath data transfer in heterogeneous wireless networks,”
*IEEE Transactions on Mobile Computing*, 2012. - Y. Cao, C. Xu, J. Guan, F. Song, and H. Zhang, “Environment-aware CMT for efficient video delivery in wireless multimedia sensor networks,”
*International Journal of Distributed Sensor Networks*, vol. 2012, Article ID 381726, 12 pages, 2012. View at Publisher · View at Google Scholar - R. Zheng, M. Zhang, Q. Wu, S. Sun, and J. Pu, “Analysis and application of bio-inspired multi-net security model,”
*International Journal of Information Security*, vol. 9, no. 1, pp. 1–17, 2010. View at Publisher · View at Google Scholar · View at Scopus - C. Xu, F. Zhao, J. Guan, H. Zhang, and G.-M. Muntean, “QoE-driven user-centric VoD services in urban multi-homed P2P-based vehicular network,”
*IEEE Transactions on Vehicular Technology*, 2012. - S. Li, M. Q. H. Meng, W. Chen et al., “SP-NN: a novel neural network approach for path planning,” in
*Proceedings of the IEEE International Conference on Robotics and Biomimetics (ROBIO '07)*, pp. 1355–1360, Sanya, China, December 2007. View at Publisher · View at Google Scholar · View at Scopus - J. Rao and A. Fapojuwo, “A battery aware distributed clustering and routing protocol for wireless sensor networks,” in
*Proceedings of the WCNC*, pp. 1538–1543, Shanghai, China, 2012. - C.-L. Wang and S.-J. Syue, “An efficient relay selection protocol for cooperative wireless sensor networks,” in
*Proceedings of the IEEE Wireless Communications and Networking Conference (WCNC '09)*, Budapest, Hungary, April 2009. View at Publisher · View at Google Scholar · View at Scopus - G. Trajcevski, F. Zhou, R. Tamassia, B. Avci, P. Scheuermann, and A. Khokhar, “Bypassing holes in sensor networks: load-balance versus latency,” in
*Proceedings of the 54th Annual IEEE Global Telecommunications Conference (GLOBECOM '11)*, Houston, Tex, USA, December 2011. View at Publisher · View at Google Scholar · View at Scopus - T. Nakagaki, H. Yamada, and Á. Tóth, “Maze-solving by an amoeboid organism,”
*Nature*, vol. 407, no. 6803, p. 470, 2000. View at Publisher · View at Google Scholar · View at Scopus - A. Tero, S. Takagi, T. Saigusa et al., “Rules for biologically inspired adaptive network design,”
*Science*, vol. 327, no. 5964, pp. 439–442, 2010. View at Publisher · View at Google Scholar · View at Scopus - A. Tero, R. Kobayashi, and T. Nakagaki, “A mathematical model for adaptive transport network in path finding by true slime mold,”
*Journal of Theoretical Biology*, vol. 244, no. 4, pp. 553–564, 2007. View at Publisher · View at Google Scholar · View at Scopus - A. Savvides, C.-C. Han, and M. B. Strivastava, “Dynamic fine-grained localization in ad-hoc networks of sensors,” in
*Proceedings of the ACM Mobicom*, pp. 166–179, New York, NY, USA, 2001. View at Scopus