- 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 383168, 17 pages
Wireless Sensor Network Modeling and Deployment Challenges in Oil and Gas Refinery Plants
1National Research Council (CNR), IEIIT Institute, 20133 Milano, Italy
2Saipem S.p.A. (ENI Group), San Donato, Italy
3DEIB, Politecnico di Milano, 20133 Milano, Italy
Received 1 November 2012; Revised 23 January 2013; Accepted 5 February 2013
Academic Editor: Marc St-Hilaire
Copyright © 2013 Stefano Savazzi 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.
Wireless sensor networks for critical industrial applications are becoming a remarkable technological paradigm. Large-scale adoption of the wireless connectivity in the field of industrial monitoring and process control is mandatorily paired with the development of tools for the prediction of the wireless link quality to mimic network planning procedures similar to conventional wired systems. In industrial sites, the radio signals are prone to blockage due to dense metallic structures. The layout of scattering objects from the existing infrastructure influences the received signal strength observed over the link and thus the quality of service (QoS). This paper surveys the most promising wireless technologies for industrial monitoring and control and proposes a novel channel model specifically tailored to predict the quality of the radio signals in environments affected by highly dense metallic building blockage. The propagation model is based on the diffraction theory, and it makes use of the 3D model of the plant to classify the links based on the number and density of the obstructions surrounding each individual radio device. Accurate link classification opens the way to the optimization of the network deployment to guarantee full end-to-end connectivity with minimal on-site redesign. The link-quality prediction method based on the classification of propagation conditions is validated by experimental measurements in two oil refinery sites using industry standard ISA SP100.11a compliant devices operating at 2.4 GHz.
The increasing demand of oil and gas supplies frequently requires the design of very large production and processing plants over remote locations with harsh environmental conditions and challenging logistics. The adoption of cabling to fully interconnect machines for process monitoring/control lacks flexibility when in large plants, and it is becoming unfeasible due to the increasing fluctuations of wiring costs to high values. The opportunity to replace cabling by deploying a network of wireless sensors is now becoming of strategic interest for several industrial applications ranging from oil and gas refining, smart factories, transport processes , and more recently oil and gas exploration .
The status of current technology allows the deployment of low-power, cost-effective network nodes in a battery-powered configuration that substitute the traditional wired devices in a very cost-effective way . The installation of wireless devices may give significant cost savings for a variety of typical plants . Current wireless networks for industrial control and monitoring are based on the IEEE 802.15.4 standard  and are mostly considered for monitoring tasks and supervised/regulatory control. The typical locations of wireless devices used for remote control and monitoring of industrial oil and gas refinery sites are characterized by harsh environments where radio signals are prone to blockage and multipath fading due to metallic structures (structural pipe racks, metallic towers and buildings, etc.) that obstruct the direct path .
With the widespread use of the wireless technology in industrial environments, the development of virtual (computer aided) network planning software tools is now becoming crucial for accurate system deployment. Inaccuracies during the radio planning design phase will turn into issues during the commissioning phase. As an example, when adding new wired nodes such as gateways and/or access points to improve the coverage, it might be required to reopen excavations along the cable route which is totally unacceptable during the commissioning (or even before the commissioning) phase of the plant. Accurate network planning limits the need to oversize the design of the overall system, which is obviously an extra cost for the contractor. Therefore, it is crucial to develop consistent design guidelines and tools that can guarantee a reasonable accuracy in the prediction of the wireless coverage. Making use of the 3D model of the deploying area (if available) during the design phase is also of utmost importance to achieve this result. An example of a 3D view of two oil refinery sites is illustrated in Figure 1: the wireless end devices (EDs), also referred to as sensors, can be connected by star or mesh mode towards a Gateway device, with the help of intermediate Repeater nodes serving as decode and forward relays. The Gateway device is collecting data and rerouting to a wired network. Network planning is based on the prediction of the pairwise wireless link qualities among all the devices in the distributed network: the link quality is expressed in terms of the strength of the received signal. The prediction can be supported by independent radio measurement campaigns over typical refinery environments and/or by models based on propagation theory and statistical or ray-tracing tools.
Conventional empirical channel models  cannot fully capture the unique propagation characteristics of the industrial environments; in addition, the ray-tracing-based models  turn out to be not practical to process the high number of structures observed in large industrial sites . This motivates the development of accurate site-specific channel models based on a small fraction of measurements taken in the refinery area.
This paper addresses a novel channel model based on the diffraction theory to assess the link quality in radio environments affected by highly dense metallic building blockage. The wireless links are partitioned into mutually exclusive classes: for each class, a separate channel model is proposed to predict the quality of the radio link. The link classification is based on the analysis of the characteristics of the obstructions that impair the wireless propagation. The 3D-CAD model of the refinery site (see Figure 1(b)) is used to identify the structure of the building blockage. Based on link classification, an optimization tool is developed for the prediction of the radio coverage and for wireless connectivity optimization. Although the channel modeling and the classification methodologies proposed in this paper are fairly general and applicable in different scenarios, the model is validated by experimental measurements using industry standard ISA SP100.11a compliant  devices operating at 2.4 GHz based on the IEEE 802.15.4-2011 physical layer. The measurement campaigns have been carried out in two sites located in a large-size oil refinery plant. Different practical deployment cases for coverage testing are discussed in environments characterized by blockage due to a high-density of metallic structures.
1.1. Wireless Industrial Networks: Applications and Technologies
A typical industrial environment shows relevant similarities with dense urban microcellular sites characterized by a harsh environment for short-range (10–50 m) radio-frequency propagation with metallic structures , changing environmental conditions, nonline of sight (NLOS), and possible colocated wireless applications running over unlicensed spectrum . Industrial networks typically require low-jitter sampling period for monitoring, high-integrity data delivery of critical messages, automatic reconfiguration, and usage of redundancy in case of communication failures. The most representative application cases for wireless technology  are commissioning, open-loop maintenance monitoring, closed-loop supervisory, and regulatory remote control. Notice that regulatory control is characterized by stricter reliability and delay requirements compared to supervisory control (some relevant application cases are primary flow and pressure control).
The commercial wireless systems predominantly use the so-called ISM bands at 2.4 GHz. Early experiments for cable replacing in regulatory control applications revealed that the traditional single-hop carrier sense multiple access (CSMA) schemes supported by WiFi (IEEE 802.11) perform poorly when adopted in a factory environment . More recently, wireless extensions to PROFIBUS protocol for critical control have been analyzed by real-time simulations . Today, commercial battery-operated systems are based on the IEEE 802.15.4 standard and enable data to be transmitted at a typical rate of 250 kbit/s, with up to a maximum of 10 dBm output RF power to meet the RF regulations for hazardous environments. The IEEE 802.15.4 physical layer also constitutes the basis for the WirelessHART  and ISA100.11a  industry standard protocols.
2. Wireless Standards for Industrial Monitoring and Control
Low-power wireless architectures and standards widely adopted in industrial automation are reviewed in this section. This introduction is instrumental to the definition of a design tool for coverage prediction and connectivity optimization. Industrial organizations such as HART and the International Society of Automation (ISA) are currently pushing towards the definition of common specifications for wireless industrial monitoring and process automation based on the IEEE 802.15.4 standard. Below we summarize the characteristics of the most relevant network solutions.
WirelessHART has been ratified by the HART Communication Foundation in 2007 as the first open wireless communication standard designed for process control applications and monitoring. Although WirelessHART adopts the IEEE 802.15.4 standard for the physical layer, the MAC layer is slightly modified as it is based on TDMA (while contention access is not allowed ) with guaranteed time slots assigned to the network devices. Frequency hopping spread spectrum access (FHSS) is used as proven technology to provide further improvements in terms of link gain compared to direct sequence spread spectrum (DSSS) option. The adoption of TDMA technology with precisely network-wide time synchronization is the key technology that makes WirelessHART different from other industry standards. Time synchronization is based on the Time Synchronized Mesh Protocol (TSMP). This method allows to synchronize transmitting and receiving node pairs by periodically correcting the relative time offsets misalignments. The offsets corrections are typically transmitted using standard ACK reply messages (with limited extra power consumption). The synchronous TDMA MAC sublayer is built upon the IEEE 802.15.4 physical layer for mesh network communication and defines superframes of 1 sec, fixed timeslot of 10 ms, channel hopping scheme supporting flexible blacklisting options and industry-standard AES-128 block ciphers with related keys.
ISA SP100.11a standard for wireless industrial automation is meant to provide the specifications for reliable and secure wireless operations for monitoring, alerting, open/closed-loop quality control, and predictive monitoring applications . The standard supports the interoperability of multiple radio technologies. The envisioned applications include wireless process control systems (with maximum latencies in the order of 1 sec). The protocol suite, system management, and security specifications are defined for low data-rate wireless connectivity based on IEEE 802.15.4 standard. Network and transport layers are based on UDP with support of IPv6-based solutions (6LoWPAN). Coexistence with other wireless services based on IEEE 802.11x, and IEEE 802.16x standards is also addressed. Although the logical link layer of ISA SP100.11a standard has a similar structure compared to WirelessHART, the standard specifies configurable timeslots with variable durations from 10 ms to 12 ms on a superframe base. Configurable timeslots ease the development of advanced architectures based on duo-casting mechanisms, optimized coexistence, and flexibility. In ISA SP100, a transaction may consist of multiple timeslots; longer transactions can be used to extend the waiting time for multiple consecutive ACKs as required in multicast transmission. The ISA standard supports both slow and fast channel hopping schemes, thus allowing devices with imprecise timing settings to perform resynchronization and neighbor discovery.
The wireless architecture supported by the standard ISA is adopted here as reference for deployment testing. As depicted in Figure 1(a), the network infrastructure consists of the following components.(i)The end devices (ED) are the input/output field instruments with the minimum set of functions that are necessary to join the network. The EDs take the role of reduced-function devices and typically do not provide any mechanism for relaying messages of other devices.(ii)The Repeaters are field EDs specifically configured to serve as relay nodes for other EDs by forming a two-hop (or multihop) mesh network. In typical industrial settings where the real-time responsiveness of the monitoring network is a crucial issue, the number of hops is limited to 2, therefore, the Repeater devices act as ED range extenders.(iii)The Gateways act as access points (or sinks) and collect the measurements acquired by the field devices. In practical settings, the Gateways are connected by cables (or by broadband wireless technology) to a common network manager node and thus also act as a translator between the ISA standard and other wired protocols (Foundation Fieldbus, HART, etc.).
3. Channel Modeling
In this section we introduce the channel model as instrumental to the proposed link classification approach. The wireless links without a clear line-of-sight (LOS) path undergo more severe received signal power attenuations than those where the line-of-sight (LOS) path is fully unobstructed. This additional attenuation is almost uncorrelated from the distance between the transmitter and the receiver . The main scatterers/objects that are responsible for the received signal power attenuation are mostly confined within the first and second Fresnel zones as these can be considered to contribute to the main propagating energy in the wavefield . For a wireless link where the direct path between the transmitter and the receiver has length , the th Fresnel zone is the region inside an ellipsoid with circular cross-section. The radius of the th Fresnel zone at distance is with the signal wavelength.
We assume that any pair of wireless devices connected with an arbitrary link are deployed at fixed locations and distance . The nodes are equipped with radio devices characterized by single omnidirectional antenna transceivers. As for typical scenarios, the Gateway antenna is mounted on an elevated point while flat terrain is assumed.
The propagation model describes the correlation between the size of the (mostly metallic) obstructions located within the Fresnel volumes and the total received signal strength (RSS) experienced along the propagation path. The RSS is thus the metric (in decibel scale) used to assess the quality of the radio link It combines (i) a static component characterized by a distance-dependent LOS term and an excess attenuation that accounts for the building blockage; (ii) a zero-mean random term accounting for the fluctuations of the received power with typical standard deviation around dB in static environments .
In what follows, it is derived a model for the static RSS component . The model is instrumental to the prediction of the average radio link quality for connectivity optimization (see Section 5). The distance-dependent static component describes the channel gain observed over the flat terrain and without obstructions. The term denotes the additional signal attenuation as a function of the size and the density of the metallic objects located within the Fresnel volume (i.e., blocking the LOS path). The model used to characterize the additional attenuation component is derived in Section 3.1. As observed in , the reflection of the radio signals from the flat terrain does not influence the attenuation parameter but only the term and the path-loss exponent . The model is validated based on measurements over the refinery sites (see Section 6).
The distance-dependent loss factor can be modeled as a function of the path-loss exponent (see ): where is the channel gain function of the transmit power and measured at a reference distance ( m typical), while is the Fresnel distance being a function of the antenna heights from the ground and for the pair of devices , respectively. Path loss exponent is typically set to in short-range environments  where ground reflections can be neglected, for . Larger path loss exponents are caused by reflections from the ground and can be experimented in long-range cases for .
The probability of successful communication depends on the random fluctuations of the RSS as in (2). Successful communication is modeled by outage probability such that . The threshold is typically set to dBm such that . Any link experiencing is assumed as unreliable, and thus, it should not be accounted for during network planning.
3.1. Diffraction Model for Prediction of Building Blockage
It is assumed that the additional attenuation in (2) is due to propagating wavefronts diffracting around the building blockage consisting of metallic obstacles with different dimensions. Obstacles are acting as perfectly absorbing interfaces.
The diffraction model for the building blockage term is based on the Fresnel-Kirchhoff method . The attenuation in (2) is obtained as a function of the received electric field : The ratio describes the obstruction loss in excess of the free space field . Large-size metallic objects obstructing the wireless link absorb a large amount of the signal intensity and limit the received field to a small fraction (being ) of the one that would be observed under free-space propagation (without obstructions). A simplified description of the propagation environment (with obstacles blocking the LOS path) is considered in Figure 2. To simplify the reasoning, we assume that the obstacles surrounding the transmitter and the receiver antennas lie in the far-field region. In addition, the shape of the obstacles obstructing the Fresnel zones is squared or rectangular. Shapes that are more typical in refinery sites (tubes, structural pipe racks, etc.) have been approximated by matching a number of rectangles in the plane to get the same shape of the obstructed areas; this is also illustrated in . For the th object, the clearance zone in the plane denotes the region corresponding to the Fresnel volume cross-section that is free from any obstacle. The shaded region in the same plane indicates instead the complementary portion of the surface occupied by the obstacle.
The Fresnel-Kirchhoff approach is used to model the field loss caused by a single th obstacle located at distance . The Huygens principle is used to predict the actual field strength diffracted by one obstacle modeled as a knife edge. The 2D model takes into account both the lateral and the vertical profiles of the obstruction by integrating the exponential phase term of the spherical wavefields over the two dimensions . The electric field measured at the receiver may be interpreted as generated by a virtual array of Huygens sources located in the plane of the single obstacle at distance from the receiver. Considering an object located at distance from the transmitter and occupying an area , the field loss can be approximated for , as  where defined in (1) refers to the radius of the 1st Fresnel volume circular section corresponding to the location of the obstruction. The approximation reasonably fits with the considered environment (see Section 6) as far as the obstacle is confined within the Fresnel volume.
To gain further insight into the interplay between the obstruction size and the corresponding field loss, in what follows we focus on the example of a single obstacle obstructing the LOS path with rectangular cross-section described by lateral and vertical half-dimensions . The loss term in (6) simplifies for the case of large obstacle , as (see the appendix) with
Figure 3 compares the diffraction loss for a single object obstructing the LOS path with varying square cross-sections () measured with respect to the Fresnel radius . Model (6) and approximation (7) are in solid and dashed lines, respectively. The field loss caused by an object fully obstructing the 1st and the 2nd Fresnel circular section such that lies below .
The general model (6) for a single obstacle can be extended to multiple obstacles by following the Deygout approach . For multiple obstacles, the lateral and vertical dimensions of the shaded region for the th obstacle are calculated with respect to the size of the largest obstacle (obstacle in the example of Figure 2). For each dimension, the Deygout method requires to find the th object (edge) with the largest value of parameters () compared to the Fresnel size, such that and and ignoring all the other edges. Based on the selected set of obstacles, a new reference plane is created for each dimension and used to compute the contributions of all the intermediate edges with modified size, , (see Figure 2). The overall obstruction loss for obstacles with meaningful obstructing size at distance is obtained by multiplying each contribution along the LOS path so that  where each term is in (6) or approximated as in (7). In spite of the simplicity of this method, in Section 6, it is proved to be accurate enough for wireless link quality prediction.
4. Wireless Link Classification
The proposed approach for the evaluation of the pairwise link channel qualities is validated by a database of radio measurements taken in different refinery sites to cover the most representative scenarios. Based on the experimental measurements, mutually exclusive link categories have been defined to account for the different sizes and the positions of the most typical obstructions inside the (1st and 2nd) Fresnel volumes surrounding the considered links. The analysis of the building blockage property is based on the inspection of the full 2D/3D model of the plant. Each link type is characterized by a specific configuration of the Fresnel zone clearance that corresponds to a reference value for the obstruction loss according to the model outlined in Section 3. For each link type , the loss is computed as in (5), and it is used to predict the average link quality in (2).
Based on the experimental activity, five different link-types are considered (see Figure 4).
Type I. LOS () link type is characterized by the absence of obstacles (with dimensions larger than the signal wavelength ) within the first and second Fresnel volume, while obstacles might instead occupy the remaining Fresnel volumes. The nominal (such that from (2) ) maximum range to guarantee a reliable connection is found as m (for RSS above dBm). In the worst-case scenario where obstacles completely obstruct the th Fresnel volumes with , the observed received electric field intensity from (6) is the fraction of the one that would be measured in the free-space case (thus corresponding to an attenuation of dB ).
Type II. Near-LOS () link type is observed in environments where the obstacles are located in the first Fresnel outer region at distance from the direct path. The shaded subregion in Figure 4 can be considered as a “forbidden” region: if this region is kept clear, then the total path attenuation will be practically the same as for the unobstructed case (Type I). This clearance zone is thus used here as a criterion to decide whether an object is to be treated as a relevant obstruction. The radio propagation for this link category is characterized by an additional signal energy loss compared to Type I. Based on the radio measurement campaigns and the diffraction model in (6), the Type II links typically retain the of the electric field observed in the free-space case ( dB). The theoretical maximum range reduces to m.
Type III. Obstructed-LOS () link type is observed in environments where the obstacles are located inside the forbidden region, although the direct path connecting the transmitter and the receiver is still unobstructed. The links belonging to this category retain approximately of the electric field that would be measured in the free-space case. The theoretical maximum range is m.
Type IV. NLOS () link type is characterized by large objects obstructing the direct path between transmitter and receiver; therefore, ( dB): the size of those objects is such that a clearance zone is still visible, , suggesting that there might be the possibility of reliable communication. Being the forbidden region and the LOS path both obstructed, the reference value for the field loss is chosen as ( dB). The theoretical maximum range further reduces to m.
Type IV-S. Severe-NLOS () link type refers to a severe NLOS environment where the first and the second Fresnel regions are completely obstructed by one or more obstacles with significant size (and dimensions scaling as ), so that the observed received electric field falls below compared to the one that would be measured in the free-space case ( dB). The theoretical maximum range is m. This model type resembles a propagation environment where the line-of-sight path is blocked by large-size concrete buildings .
5. Radio Planning Optimization
The wireless network deployment problem refers to the determination of the positions of the wireless nodes such that some limiting values of coverage, connectivity, and energy efficiency can be achieved . Wireless device deployment strategies for coverage and connectivity enhancement play a crucial role in providing better quality of service (QoS) to the network. The coverage and the connectivity problems are two fundamental issues that have been widely studied in the literature . In coverage problems, the objective is to deploy wireless sensor devices in strategic ways such that an optimal area coverage is achieved given the requirements of the underlying application . The coverage problem therefore deals with placing a minimum number of nodes so that every measurement point in the sensing field is optimally covered according to application-specific constraints. In industrial monitoring and control applications, the position of the measurement points (sensors or actuators) is constrained by the application; therefore, the coverage optimization is typically carried out based on the structure of the process unit. The focus of this section is thus on connectivity optimization, as this is the most crucial problem for cable replacing in the industrial networking context.
In what follows, we first report on the current state of the research on optimized node placement in wireless sensor networks (Section 5.1). Next, we discuss relevant practical issues and rules that are specifically tailored for network and connectivity optimization in industrial networks (Section 5.2). Finally, we propose an optimization framework tailored for commercially available ISA SP100 two-hop networks that allows the optimal selection of the devices that need to be configured as Repeaters (Section 5.3). The goal is to optimize the number and the position of the infrastructure devices (e.g., the Repeater nodes and/or the Gateways) to guarantee a reliable connection between the measurement points and the control unit with some degree of redundancy . Optimally deployed wireless infrastructure devices guarantee adequate QoS (i.e., outage probability), long network lifetime, and thus reduced costs for network maintenance. The proposed deployment problem is based on the prediction of the RSS for all the pairwise wireless links according to channel modeling and classification outlined in Sections 3 and 4.
5.1. Node Deployment Strategies in Wireless Sensor Networks: A Survey
Extensive work has been reported in the literature relating to wireless sensor and relay node deployment. Deployment of nodes has been considered for targeting connectivity, coverage, node lifetime, and/or QoS. The deployment strategies can be classified into static and dynamic  depending on whether the optimization is performed during network setup or during network operation (for node repositioning, see ). In static environments where data is periodically collected over preset routes, the problem of optimal node placement for connectivity maximization has been proven to be NP-hard for most of the formulations . Several heuristics and rules have been therefore proposed to find suboptimal solutions based on graph theory. Several approaches to the problem of placing nodes are addressed in  to achieve -connectivity at the network setup time so that independent paths are identified for every pair of devices. The majority of published work on sensor network deployment limits its focus on simplified and analytically tractable 1D and 2D environments where connectivity can be considered as a primary/secondary objective or as a constraint in the deployment problem . For example, in  an outdoor random deployment which targets the connectivity as a primary objective in 2D space is investigated. In , a constrained multivariable nonlinear programming problem is analyzed to determine the locations of the sensor nodes to maximize the network lifetime, given a fixed number of sensor nodes with certain coverage and connectivity requirements. A deployment strategy for sensor networks is introduced in  to balance the network lifetime and connectivity goals for single- and two-hop networks.
Focusing on large-scale sensor network applications, controlled placement of nodes is often focused on a subset of network devices (e.g., Repeaters or relays) with the goal of designing the network topology to achieve the desired application requirements . The problem of relay placement in two-hop networks is analyzed in : the objective is to place the fewest number of relay nodes so that each sensor node can communicate with at least one relay node, and the network of relay nodes is connected. The goal is to guarantee a reliable communication between each pairs of sensor nodes while the same reasoning can be extended for sensors communicating with a common Gateway node. Recent literature considers the problem of connectivity in massively dense sensor networks . The problem of deploying relay nodes in heterogeneous sensor network scenarios is considered in  where sensor and relay nodes possess different transmission ranges (e.g., through the use of different hardware, antennas, or high-power radio modules). The work  considers a scenario where sensor devices are equipped with directional antennas: the goal is to find an optimal subset of locations to minimize the total network cost while satisfying the requirements of coverage and connectivity.
The network connectivity problem is mostly considered for 2D planning with the assumption of simple binary communication disk model without looking at site-specific environmental constraints (see also [34–37]). Those approaches are very prone to failure in practical large-scale industrial applications. Some attempts in the literature have been made towards the analysis of deployment and connectivity problems in 3D environments, although the topic is still considered an open issue . The problem of modeling and connectivity optimization in random 3D networks has been recently addressed in [38, 39] where the deployment problem considers the maximization of network connectivity satisfying lifetime constraints.
5.2. Connectivity Optimization in Wireless Industrial Networks
The connectivity optimization for industrial networks can be in general applied to two-hop large-scale networks consisting of Gateways, relays, and sensors, operating in time (and safety) critical applications [36, 37]. Three general practical rules  should be followed during system design and configuration. These are summarized below.
Gateway Deployment Planning. The wireless network is first divided from a single process unit into subsections (subnetworks). Within each subsection, the position of the measurement points, and thus the degree of coverage, is designed to satisfy application-dependent requirements. The devices (or end devices, EDs) are deployed to collect data from the nearby measurement points (depending on the monitored process, EDs might consists of a single or multiple measurement points). Each process unit subsection is served by one Gateway (acting as access point for the corresponding devices). The Gateway should be able to allocate resources for two-way communication in real time with the EDs. For small-size projects (as those analyzed in Section 6), a single Gateway is sufficient if the total number of measurement points is less than the capacity of the Gateway point. Instead, if the project is large with several hundreds of wireless devices and process units, a single network manager should manage multiple Gateways. The required number of Gateways can be defined as a function of the number of measurement points. The following simple calculation can be used in practice to approximate the number of Gateways needed: where is the spare (or residual) fraction of the available capacity to be reserved for emergency signalling with capacity measured in terms of number of measurement points served. is the number of measurement devices assuming that each ED is serving as a single measurement point. A typical design rule prescribes that . The Gateway capacity depends on the wired/wireless protocol used for data transfer towards the network manager.
Connectivity Optimization. The use of site-specific radio propagation models (empirical or ray-tracing based) enables the optimization of the connectivity for virtual network planning. A propagation model can be therefore exploited as instrumental to the prediction of the RSS, from which the quality of the radio link (and of the end-to-end connectivity) can be inferred with some degree of accuracy. Prediction errors are typically caused by modeling mismatches (e.g., link classification errors) or unpredictable RSS fluctuations (see Section 3) due to interference over the 2.4 GHz band or fading induced by objects or people moving in the area. The solution to the connectivity and the Gateway deployment problems is generally well understood in the literature (see, e.g., ). In the industrial context, three practical rules are defined to ensure a sufficiently high link reliability. The rules are summarized as follows (see also ).(i)Rule 1. Every network with more than 5 devices should have a minimum of 25% of devices within the effective range of the Gateway to ensure mesh connection (typically over a maximum number of hops). In any case, every network should have a minimum of 5 infrastructure devices within the effective range of the Gateway. Example: a network consisting of 100 EDs requires 25 EDs at minimum within the effective range of the Gateway (directly connected).(ii)Rule 2. Gateway RF antenna should be mounted at least 2 m from the ground level and should not be surrounded by obstacles. Obstacles should lie at distance from the antenna.(iii)Rule 3. Every device should have a minimum of 3 neighbors in the effective range. This ensures that when implemented, there will be at least one reliable routing path to the Gateway alternative to direct connection (to guarantee connectivity).
On-Site Stress Testing. Stress testing of the deployment design is recommended during an on-site survey to verify potential weaknesses highlighted during the virtual network configuration. Stress testing is performed by altering the position of the EDs from the nominal position and thus by measuring the fluctuations of the RSS field.
Although the context may vary slightly depending on the structure of the environment, almost all these basic steps could be applied regardless of the specific commercial system and standard (i.e., WirelessHART or ISA SP100, see Section 2). The first and the second steps are known to be the most critical for high density applications .
5.3. Optimal Repeater Configuration for Two-Hop ISA Industrial Networks
The wireless network for industrial environment under consideration conforms with the standard ISA SP100.11a and is characterized by one Gateway collecting data from wireless end devices (EDs). A subset of EDs might serve as Repeater nodes acting as range extenders. The Gateway node is an electrically powered device, serving as access point for the EDs. It manages both wireless and wired interfaces. The Repeaters are configured as EDs with superior functionalities: these allow to connect to the Gateway and simultaneously serve as decode and forward relays for extending the range of the neighboring EDs. Repeater nodes are more expensive than standard EDs since they must be preconfigured to multiplex different sensor data and could be more powerful in terms of processing and transmission capabilities.
The connectivity optimization problem is therefore focused on the Repeater configuration. The candidate sites for the deployment of the EDs and of the Gateway node are assumed to be assigned: each candidate site might host either a Repeater node multiplexing sensor data or a standard ED without relaying functionalities. Optimal placement of the Gateway is not addressed in this paper, although we only assume that Gateway locations satisfy connectivity Rule 2 (see Section 5.2). Optimal deployment for the Gateway might be carried out as illustrated in , even if, in practical industrial scenarios, the exact position is subject to stringent environmental constraints.
The optimization approach is based on the selection of the smallest subset of devices that need to be configured as Repeaters to guarantee network connectivity. The optimization jointly minimizes the number of EDs connected to the corresponding Gateway over two hops  and guarantees a minimum quality of service for all links, so that the static RSS component is kept for all the configured links above the system threshold (see Section 3) herein adopted as the minimum tolerable link quality. The static RSS component is predicted based on the 3D model of the plant, as done in Section 3.
Let the wireless network be represented by a set of nodes located at known positions within a specific area of the plant. A sequence of messages is continuously transmitted by the EDs towards a common Gateway node labeled as “0” possibly with the help of one intermediate node serving as a Repeater. To comply with the real-time responsiveness constraints typically required by industrial closed-loop control applications, the maximum number of hops to reach the Gateway node is herein limited to 2: the same constraint is also adopted in recent ISA compliant network implementations. Any wireless node is said to be connected with reasonable quality to the Gateway node “0” if and only if where the indicator for an arbitrary link is defined as while sensitivity threshold accounts for the random fluctuations of the static RSS component . Deployment optimization consists of three phases.
Selection of Candidate Repeaters. First, it is defined the subset of nodes without direct connection to the Gateway The same nodes are preconfigured as EDs, and they should not provide relaying functionalities. The remaining subset of devices observing a reliable connection with the Gateway can be assigned either as Repeaters or EDs. The optimal configuration of devices in subset is carried out in the following steps.
Feasibility Region for the Connectivity Problem. Assuming all nodes be initially configured as Repeaters, a solution to the connectivity problem satisfying Rule 3 for the EDs without a reliable direct connection with the Gateway exists if such that all the EDs can exploit an alternative two-hop link through one Repeater serving as a range extender. In case the condition is not satisfied, the number of candidate points is not sufficient for a feasible solution to the coverage problem; additional candidate sites must be therefore identified. New candidate sites must be assigned during the precommissioning of the plant; the deployment should thus account for application and site-specific environmental constraints.
Repeater Configuration. Among the potential subsets of devices configured as Repeaters, with , the optimal subset is defined as the one satisfying the feasibility region (14) and with the smallest cardinality. By letting be the cardinality of the th subset , the algorithm identifies the optimal th subset of the Repeater devices such that The devices are thus configured as Repeaters while the other devices take the role of EDs. Notice that .
The iterative algorithm described as follows is used to find a solution to problem (15). Let the ordering of the Repeater subsets be such that ; the algorithm starts by picking the largest feasible set of Repeaters, so that and iteratively identifies new feasible subsets with smaller cardinality () by randomly removing nodes from . The optimal subset solution to (15) is such that any smaller subset of Repeaters with is not feasible as or, equivalently, for any Repeater subset with smaller cardinality the sum for some ED without reliable direct connection.
6. Experimental Activity
The experimental validation of network connectivity is based on the link classification and channel modeling described in Sections 3 and 4. The optimization tool used for the optimal selection of the Repeater nodes is described in Section 5. The connectivity design consists of three steps. At first, the candidate positions for the wireless devices are chosen to highlight practical cases of meaningful interest for the deployment of an industrial network. The Gateway node is mounted above ground (according to Rule 2 in Section 5.2) and collects the data received from all the EDs. Second, the pairwise link RSSs are predicted based on channel modeling and classification as outlined in Sections 3 and 4. Finally, the optimal sub-set of the wireless devices that should act as Repeaters is computed based on the connectivity optimization tool illustrated in Section 5.3.
In the proposed experimental set-up, we deployed absolute and gauge pressure transmitters communicating with a Gateway node by star or two-hop mesh topology. Compared to mesh topology, deploying a star topology network should be preferred in practice as it provides better performance in terms of per-link real-time responsiveness that is required for monitoring and control of critical plant parameters. The radio transceivers conform with the ISA SP100.11a protocol  with radio transmit power set to dBm. The experiments have been carried out in two sites within the same oil refinery: the first site is a 100 m × 200 m area around a flare unit; the second one is a 60 m × 30 m area surrounding a furnace structure. All the environments under consideration are characterized by metallic objects and concrete buildings with high-reflectivity surfaces. Before the test, we used a signal analyzer to characterize the interferers in the area. Since no significant activity was detected, the IEEE 802.15.4 channels selected for the experiments have center frequencies 2.405 GHz and 2.480 GHz, corresponding to the ISA SP100.11a channel numbers 1 and 15, respectively.
The static RSS component in (2) characterizing the radio propagation over each link is predicted by following three steps.(i)Step 1. The number and size of the objects blocking the direct path between the transmitter and receiver pair (or the corresponding Fresnel volumes) are identified by analyzing the 3D model of the plant.(ii)Step 2. The link is classified by exploring the 3D maps of the corresponding sites. Based on the link types identified in Section 4, the size of the obstructions is compared with the Fresnel volumes to identify the corresponding clearance zones for the obstacles with relevant size compared to the wavelength . The link category is then selected by comparing the resulting clearance zones with the ones characterizing each link type.(iii)Step 3. The static RSS component is predicted according to the chosen link type. The signal attenuation in (2) for the chosen link category is computed based on the predicted field loss as in (5). The distance-dependent loss factor is defined according to the position of the transmitter and receiver devices. In all of the considered short range cases for , the path loss exponent is . The propagation over long ranges such that suffers from a larger path loss due to ground reflections: for a typical case of m (Gateway at 6 m from the ground), the available measurements indicate a path loss exponent of .
The measurements analyzed in the following sections highlight the accuracy of the proposed channel characterization and modeling approach. Tightness of the proposed model is verified by comparing the predicted RSSs with the corresponding measurements obtained during the on-line testing. The model accuracy is found as reasonably high in all the considered settings (with errors below 4 dB for all the considered cases).
6.1. Site Test No. 1: Flare Unit
In this test, the Gateway is mounted in 4 different locations corresponding to different deployment cases as illustrated in the floor plan maps of Figure 5. For deployment case no. 1 the height from the ground of the Gateway is 1.5 m (m); for the remaining cases; the height is above 6 m (m). For all cases, the 3 EDs labelled as B, C, and D are acting as input/output field devices and are moved in different positions labeled by lowercase letters (, , and ). The corresponding RSS measurements are reported by circle markers for all the deployment cases and analyzed in Figure 6 for devices at ground level and in Figure 7 for devices at 1 m above ground. The markers have different colors to identify the link category while the link classification is based on the inspection of 2D and 3D-CAD maps. The predicted static RSS component is represented by solid lines as a function of the distance and for each link category (see Section 4). The same color code used for the measurements is adopted to highlight the prediction and link classification accuracy.
The effectiveness of the proposed channel characterization and modeling approach can be appreciated in several settings as highlighted in Figure 8. To focus on a relevant example, in the deployment case 3, the ED transmitters located at positions C3a (ground level) and C3b (1 m height from the ground) are hidden behind a big cylindrical vessel that completely obstruct the 1st Fresnel region. The wireless links connected to the Gateway retain the and the of the received field that would be measured in free space, respectively. Therefore, they can be reasonably classified as Type IV-S. As confirmed by measurements, the predicted RSS is below the critical dBm reliability threshold (distance m) suggesting the need for a Repeater device acting as relay. The same transmitter is now moved at position C3c to circumvent the large obstruction and create more favorable propagation conditions. In this case, by analyzing the corresponding 3D map, the 1st Fresnel region is slightly unobstructed: the link retains a larger fraction () of the received electric field and thus can be reasonably classified as Type IV. As confirmed by the chosen model, the connection with the Gateway is now reliable as RSS dBm: this suggests to deploy a Repeater device at position C3c multiplexing the source data received from the devices obstructed by the cylindrical vessel at position C3a and C3b. As confirmed by analysis of the 3D model, the links connecting the Repeater with devices located at the other side of the vessel can be classified as Type II, being the forbidden region free from obstacles.
Figure 8 highlights other relevant deployment cases: the links connecting the Gateway with the EDs at position B4a and C4a are classified as Type IV () and Type III (), respectively. For both links the forbidden region is found as partially obstructed: in addition, at position B4a, the LOS path is blocked by concrete and metallic structures located around the corresponding ED location. For position D4a instead, the forbidden region is found as unobstructed; the corresponding link can be thus classified as Type II ().
6.2. Site Test No. 2: Furnace Structure
In this test, the Gateway is mounted on the stairway in the south-east of the furnace at 10 m above the ground level. In this scenario, devices C and D are moved over four different floors of the furnace structure according to Figure 9. Device B instead is located at ground level, moved in positions in front of the furnace structure. The distance between each device and the Gateway ranges between 14 m and 57 m and is lower than the Fresnel distance m in all cases. Measurements and predicted model for each link category are reported in Figure 10 using the same color code adopted for the flare unit scenario. By exploring the 2D and the 3D maps of the site, the links corresponding to positions B5 ( and ), D5e, and C5a can be reasonably classified as Type III ( or dB), being the forbidden region (see Figure 4) partially obstructed. Measured attenuations are dB and confirm this choice. As highlighted in Figure 11, the NLOS links (Type IV) correspond to positions B5c, D5a, and C5e with observed attenuation ranging from dB. For positions D5d (4rd floor) and C5d (3rd floor), the metallic structure produces a waveguide effect on propagation such that reliable communication occurs even across the whole furnace structure. The wireless signals propagate all around the furnace environment without obstacles and take advantage of the constructive interference. The positions D5b (2rd floor) and C5b (1rd floor) are instead surrounded by the furnace building that fully obstructs the 1st Fresnel volume and absorbs approximately the 84% and the 88% of the free-space field intensity, with and , respectively (Type IV-S). As confirmed by measurements, the predicted RSS for Type IV-S link (with distance m) is below the critical dBm reliability threshold. The installation of one Repeater device located in the example at position C5a is therefore the optimal choice to relay the data acquired by the measurement points located at positions D5b and C5b.
6.3. Long-Range Testing
Although the focus of this paper is mostly on short-range networking modeling and optimization of network deployment in industrial environments, a long-range test have been also carried out as depicted in Figure 12 (deployment case 5) with the Gateway located in the same position of case 4 while the device C at ground level has been moved in two sites. The first one is an open area classified as near LOS environment (Type II) on the right side of the flare unit at distance 109 m from the Gateway, the second site was located at distance 132 m from the Gateway in the southern part of the flare unit where the LOS path is obstructed by a building. The path loss caused by the ground reflections (flat terrain) can be reasonably modeled as in (2) with exponent . For the first test, the wireless link is characterized by . The measured RSS of dBm confirms the predicted range for the corresponding Type II link category (see Section 3). In the second test, the link only retains the of the free-space electric field (the attenuation caused by the building is 7 dB), and it is classified as unreliable with RSS of dBm.
7. Concluding Remarks
Network deployment in industrial settings with dense metallic structures can be based on a simple but effective channel model that makes use of the diffraction theory for 3D environments. The model proves to characterize the wireless propagation in industrial environments with an accuracy that is reasonably high to predict the average quality of the wireless links in different sites. The wireless links are partitioned into mutually exclusive attenuation classes (link types) based on the 3D structure of the building blockage. Each class is characterized by a different amount of obstruction loss; therefore, a separate channel model is proposed to predict the QoS for each link type. The diffraction model is then adopted for virtual planning of two-hop ISA networks: the problem of optimal Repeater configuration of the Repeater devices is addressed to guarantee reliable connectivity between the end devices and the Gateway. The proposed classification approach has been validated by extensive experimental measurements in critical areas within an oil refinery plant characterized by highly dense metallic structure. Industry ISA SP100.11a standard devices operating at 2.4 GHz are adopted. Experimental results from the surveys confirm the effectiveness of the proposed method as it provides a practical tool for virtual network planning with reasonable accuracy that meets the expectations in several industrial settings.
For the symmetric rectangular obstacle case, the loss term in (6) simplifies as where we used the substitutions and . Using an asymptotic expansion  for the integrals in the form valid for large enough with defined in (8), the loss term can be written now as in (7).
This work has been performed in the framework of the European Research Project DIWINE (Dense Cooperative Wireless Cloud Network) under FP7 ICT Objective 1.1—The Network of the Future.
- D. Zuehlke, “Smart factory—towards a factory-of-things,” Annual Reviews in Control, vol. 34, no. 1, pp. 129–138, 2010.
- S. Savazzi, U. Spagnolini, L. Goratti, D. Molteni, M. Latva-aho, and M. Nicoli, “Ultra-wide band sensor networks in oil and gas explorations,” IEEE Communications Magazine. In press.
- A. Willig, “Recent and emerging topics in wireless industrial communications: a selection,” IEEE Transactions on Industrial Informatics, vol. 4, no. 2, pp. 102–124, 2008.
- S. Savazzi, S. Guardiano, and U. Spagnolini, “Wireless critical process control in oil and gas refinery plants,” in Proceedings of the IEEE International Conference on Industrial Technology (ICIT '12), pp. 1003–1008, Athens, Greece, March 2012.
- Standard IEEE 802.15.4-2006, “Part 15.4: Wireless Medium Access Control (MAC) and Physical layer (PHY) specifications for low-rate Wireless Personal Area Networks,” 2006.
- L. Tang, K.-C. Wang, Y. Huang, and F. Gu, “Channel characterization and link quality assessment of IEEE 802.15.4-compliant radio for factory environments,” IEEE Transactions on Industrial Informatics, vol. 3, no. 2, pp. 99–110, 2007.
- V. Erceg, L. J. Greenstein, S. Y. Tjandra et al., “An empirically based path loss model for wireless channels in suburban environments,” IEEE Journal on Selected Areas in Communications, vol. 17, no. 7, pp. 1205–1211, 1999.
- S. J. Fortune, D. M. Gay, B. W. Kernighan, O. Landron, R. A. Valenzuela, and M. H. Wright, “WISE design of indoor wireless systems: practical computation and optimization,” IEEE Computational Science and Engineering, vol. 2, no. 1, pp. 58–68, 1995.
- J. Lei, L. Greenstein, and R. Yates, “Link gain matrix estimation in distributed wireless networks,” in Proceedings of IEEE Global Telecommunications Conference (GLOBECOM '08), pp. 848–852, New Orleans, La, USA, December 2008.
- Standard ISA100.11a-2009, “Wireless systems for industrial automation: process control and related applications,” ISA, July 2009.
- M. Baldi, R. Giacomelli, and G. Marchetto, “Time-driven access and forwarding for industrial wireless multihop networks,” IEEE Transactions on Industrial Informatics, vol. 5, no. 2, pp. 99–112, 2009.
- D. Christin, P. S. Mogre, and M. Hollick, “Survey on wireless sensor network technologies for industrial automation: the security and quality of service perspectives,” Future Internet, no. 2, pp. 96–125, 1999.
- G. Cena, I. C. Bertolotti, A. Valenzano, and C. Zunino, “Evaluation of response times in industrial WLANs,” IEEE Transactions on Industrial Informatics, vol. 3, no. 3, pp. 191–201, 2007.
- P. B. Sousa and L. L. Ferreira, “Hybrid wired/wireless profibus architectures: performance study based on simulation models,” EURASIP Journal on Wireless Communications and Networking, vol. 2010, Article ID 845792, 25 pages, 2010.
- J. Song, S. Han, A. K. Mok et al., “WirelessHART: applying wireless technology in real-time industrial process control,” in Proceedings of the 14th IEEE Real-Time and Embedded Technology and Applications Symposium (RTAS '08), pp. 377–386, April 2008.
- D. J. Y. Lee and W. C. Y. Lee, “Propagation prediction in and through buildings,” IEEE Transactions on Vehicular Technology, vol. 49, no. 5, pp. 1529–1533, 2000.
- S. R. Saunders and A. A. Zavala, Antennas and Propagation for Wireless Communications Systems, Wiley, New York, NY, USA, 2nd edition, 2007.
- W. C. Y. Lee and D. J. Y. Lee, “Microcell prediction in dense urban area,” IEEE Transactions on Vehicular Technology, vol. 47, no. 1, pp. 246–253, 1998.
- H. Mokhtari and P. Lazaridis, “Comparative study of lateral profile knife-edge diffraction and ray tracing technique using GTD in urban environment,” IEEE Transactions on Vehicular Technology, vol. 48, no. 1, pp. 255–261, 1999.
- A. Nyuli and B. Szekeres, “An improved method for calculating the diffraction loss of natural and man made obstacles,” in Proceedings of Personal, Indoor and Mobile Radio Communications, pp. 426–430, Boston, Mass, USA, October 1992.
- C. L. Giovanelli, “An analysis of simplified solutions for multiple knife-edge diffraction,” IEEE Transactions on Antennas and Propagation, vol. 32, no. 3, pp. 297–301, 1984.
- G. E. Athanasiadou, “Incorporating the Fresnel zone theory in ray tracing for propagation modelling of fixed wireless access channels,” in Proceedings of the 18th Annual IEEE International Symposium on Personal, Indoor and Mobile Radio Communications (PIMRC '07), pp. 3–7, September 2007.
- C. F. Huang and Y. C. Tseng, “The coverage problem in a wireless sensor network,” Mobile Networks and Applications, vol. 10, no. 4, pp. 519–528, 2005.
- A. Ghosh and S. K. Das, “Coverage and connectivity issues in wireless sensor networks: a survey,” Pervasive and Mobile Computing, vol. 4, no. 3, pp. 303–334, 2008.
- R. Chandra, L. Qiu, K. Jain, and M. Mahdian, “Optimizing the placement of integration points in multi-hop wireless networks,” in Proceedings of the 12th IEEE International Conference on Network Protocols (ICNP '04), pp. 271–282, Berlin, Germany, 2004.
- M. Younis and K. Akkaya, “Strategies and techniques for node placement in wireless sensor networks: a survey,” Ad Hoc Networks, vol. 6, no. 4, pp. 621–655, 2008.
- X. Cheng, D.-Z. Du, L. Wang, and B. Xu, “Relay sensor placement in wireless sensor networks,” Wireless Networks, vol. 14, no. 3, pp. 347–355, 2008.
- M. Ishizuka and M. Aida, “Performance study of node placement in sensor networks,” in Proceedings of the 24th International Conference on Distributed Computing Systems Workshops (ICDCSW '04), vol. 7, pp. 598–603, Tokyo, Japan, March 2004.
- P. Cheng, C. N. Chuah, and X. Liu, “Energy-aware node placement in wireless sensor networks,” in Proceedings of the 47th IEEE Global Telecommunications Conference (GLOBECOM '04), pp. 3210–3214, Dallas, Tex, USA, December 2004.
- K. Xu, H. Hassanein, and G. Takahara, “Relay node deployment strategies in heterogeneous wireless sensor networks: multiple-hop communication case,” in Proceedings of the 2nd Annual IEEE Communications Society Conference on Sensor and AdHoc Communications and Networks (SECON '05), pp. 575–585, Santa Clara, Calif, USA, September 2005.
- M. Haghpanahi, M. Kalantari, and M. Shayman, “Topology control in large-scale wireless sensor networks between information source and sink,” Ad Hoc Networks, 2012.
- J. Tang, B. Hao, and A. Sen, “Relay node placement in large scale wireless sensor networks,” Computer Communications, vol. 29, no. 4, pp. 490–501, 2006.
- A. Konstantinidis and K. Yang, “Multi-objective energy-efficient dense deployment in wireless sensor networks using a hybrid problem-specific MOEA/D,” Applied Soft Computing, vol. 12, no. 7, pp. 1847–1864, 2012.
- X. Han, X. Cao, E. L. Lloyd, and C. C. Shen, “Fault-tolerant relay node placement in heterogeneous wireless sensor networks,” IEEE Transactions on Mobile Computing, vol. 9, no. 5, pp. 643–656, 2010.
- Y. E. Osais, M. St-Hilaire, and F. R. Yu, “Directional sensor placement with optimal sensing range, field of view and orientation,” Mobile Networks and Applications, vol. 15, no. 2, pp. 216–225, 2010.
- J. Lee, T. Kwon, and J. Song, “Group connectivity model for industrial wireless sensor networks,” IEEE Transactions on Industrial Electronics, vol. 57, no. 5, pp. 1835–1844, 2010.
- T. M. Deyab, U. Baroudi, and S. Z. Selim, “Optimal placement of heterogeneous wireless sensor and relay nodes,” in Proceedings of the 7th International Wireless Communications and Mobile Computing Conference (IWCMC '11), pp. 65–70, July 2011.
- F. M. Al-Turjman, H. S. Hassanein, and M. A. Ibnkahla, “Efficient deployment of wireless sensor networks targeting environment monitoring applications,” Computer Communications, vol. 36, no. 2, pp. 135–148, 2013.
- X. Bai, C. Zhang, D. Xuan, and W. Jia, “Full-coverage and k-connectivity (k =14, 6) three dimensional networks,” in Proceedings of the 28th Conference on Computer Communications (INFOCOM '09), pp. 388–396, April 2009.
- WirelessHART, “IEC 62591, System Engineering Guide,” Revision 2, October 2010.
- B. Aoun, R. Boutaba, Y. Iraqi, and G. Kenward, “Gateway placement optimization in wireless mesh networks with QoS constraints,” IEEE Journal on Selected Areas in Communications, vol. 24, no. 11, pp. 2127–2136, 2006.
- M. Abramovitz and I. A. Stegun, Handbook of Mathematical Functions, Applied Mathematics Series 55, Dover, New York, NY, USA, 1972.