(Translated by https://www.hiragana.jp/)
A Graph-based Strategic Sensor Deployment Approach for π‘˜-coverage in WSN

A Graph-based Strategic Sensor
Deployment Approach for kπ‘˜kitalic_k-coverage in WSN

Lakshmikanta Sau, Priyadarshi Mukherjee, and Sasthi C. Ghosh Advanced Computing & Microelectronics Unit
Indian Statistical Institute, Kolkata 700108, India
Emails: lakshmikanta_r@isical.ac.in, priyadarshi@ieee.org, sasthi@isical.ac.in
Abstract

This paper studies a graph-based sensor deployment approach in wireless sensor networks (WSNs). Specifically, in today’s world, where sensors are everywhere, detecting various attributes like temperature and movement, their deteriorating lifetime is indeed a very concerning issue. In many scenarios, these sensors are placed in extremely remote areas, where maintenance becomes challenging. As a result, it is not very wise to depend on a single sensor to obtain data from a particular terrain or place. Hence, multiple sensors are deployed in these places, such that no problem arises if one or few of them fail. In this work, this problem of intelligent placement of sensors is modelled from the graph theoretic point of view. We propose a new sensor deployment approach here, which results in lesser sensor density per unit area and less number of sensors as compared to the existing benchmark schemes. Finally, the numerical results also support our claims and provide insights regarding the selection of parameters that enhance the system performance.

Index Terms:
Wireless sensor network, planer graph, polygon tiling, regular hexagon, kπ‘˜kitalic_k-coverage.

I Introduction

With the rapid evolution of applications like Internet of Things (IoT), wireless sensors are everywhere. Wireless sensor network (WSN) is a reality today; the global IoT connections is forecast to reach around 38.938.938.938.9 billion in 2029202920292029 [1]. As a result, the lifetime of these sensors is a very critical issue and hence, they need to be regularly monitored and recharged. In this context, there has been research regarding efficient environment aware wireless transmission techniques, which extend the sensor lifetime by avoiding unwanted transmissions [2],[3]. However, in many applications, they are placed in very remote and challenging areas, where their maintenance is extremely challenging. Thus, it is not wise to rely solely on a single sensor placed in such areas for data. In graph theoretical terms, this is termed as 1111-coverage, i.e., every point in the region is monitored or covered by at least one sensor.

This leads to a more evolved concept of graph theory, which is kπ‘˜kitalic_k-coverage, i.e., all the points in a concerned β€˜region’ is covered or monitored by at least kπ‘˜kitalic_k sensors. Note that, the value of kπ‘˜kitalic_k depends on the application at hand; a higher value of kπ‘˜kitalic_k implies that a higher degree of reliability is required, and vice-versa [4]. The simplest possible strategy to attain a certain coverage is the one, where we deploy all the available sensors in the region of concern, such that the defined set of targets are totally covered. However, as stated earlier, these sensors have a limited lifetime and hence, this method is not desirable from the application point of view [5]. Also, this leads to unwanted wastage of resources. The authors in [6] proposed a bound on the sensor density in order to cover a certain region of interest. The work in [7] proposed a sensor deployment strategy, where they used regular hexagons of a certain side length to cover the entire region and accordingly, place the sensors inside.

In this context, we propose an intelligent graph-based sensor deployment strategy. Specifically, Section II introduces the essential terminology required throughout the work and also the respective system model. Section III discusses the proposed strategy by modeling the sensor deployment problem as a kπ‘˜kitalic_k-coverage problem, where it is assumed that the sensors are placed in the environment as a regular polygon tiling. To do so, we assess and examine various planar convex regular polygons, which possess the ability to tile the entire Euclidean plane. As hexagons prove to be the most efficient of the regular polygons, we assume hexagonal tiling. Based on hexagonal tiling of the region of interest, we strategically place the sensors inside a regular hexagon, and by employing this method for every hexagon, we cover the entire Euclidean plane. Moreover, we characterize the performance of the proposed strategy in terms of sensor density and number of sensors required to cover a certain region of interest. We show that our proposed strategy outperforms the state of the art deployment strategy [7] in terms of both sensor density and total number of sensors required. Finally, Section IV presents the numerical results and Section V concludes our work.

II System Model and Preliminaries

Here we first introduce the essential terminology, which is used throughout the work and then, the network model.

II-A Definitions

Definition 1.

(Regular polygon) : A closed geometric figure formed by edges, not necessarily equal, is called a polygon. Moreover, a polygon with equal edges and identical internal angles is said to be a regular polygon.

Definition 2.

(Regular polygon tiling) : The act of tiles covering up an Euclidean plane is called tiling, where each point on the Euclidean plane is covered by at least one tile. If the Euclidean plane is covered by regular polygons, then it is called a regular polygon tiling.

Definition 3.

(Cover set) : If a set of regular polygons covers an Euclidean plane, then this particular set of polygons is said to be a cover set of this Euclidean plane. Here, an Euclidean plane is said to be covered if each point on the plane is covered by at least one regular polygon.

Definition 4.

(kπ‘˜kitalic_k-coverage set) : A cover set is defined as a kπ‘˜kitalic_k-coverage set if each point of an Euclidean plane is covered by at least kπ‘˜kitalic_k number of regular polygons.

Definition 5.

(Sensing range) : In a WSN, each sensor monitors the physical conditions of the environment and/or the wireless network up to a certain distance, i.e., it can sense its desired attribute within a circle of fixed radius rπ‘Ÿritalic_r. This circle (with radius rπ‘Ÿritalic_r) is defined as the sensing range of the sensor.

Definition 6.

(Communication range) : A sensor monitors the physical conditions within its sensing range and communicates the desired attribute to a destination. In this context, the communication range of a sensor s𝑠sitalic_s is defined as the maximum distance from s𝑠sitalic_s, within which it can communicate directly with another sensor.

Definition 7.

(Adjacent polygons) : Two regular hexagons, of equal sides, are said to be adjacent to each other if they share a common edge.

Definition 8.

(Overlapping polygons) : Two regular hexagons, of equal sides, are said to be overlapping, if they contain at least one common point, which lies inside both the hexagons.

Refer to caption
Figure 1: Two different tilling of Euclidean plane.

II-B Network Model

In a typical WSN, sensors are placed throughout the environment to monitor certain physical conditions. In many cases, they are placed in extremely remote areas to collect useful data and as a result, their maintenance becomes extremely challenging in that scenario. Therefore, it is not wise to depend on a single sensor to obtain the useful data from a particular point. As a result, we want each point of a particular place to be covered by multiple sensors, such that there is no problem in case one or some of them fail. We model this as a kπ‘˜kitalic_k-coverage problem, i.e., each point of a particular area in the Euclidean plane is covered by at least kπ‘˜kitalic_k regular polygons. Here, we assume that each sensor of the WSN have similar characteristics and follow identical lifetime models. We further assume, that all the sensors are aware of their own position with respect to the entire network and also, they possess full knowledge of all the other sensors by means of a common central system. Finally, we characterize the sensing range of sensor sisubscript𝑠𝑖s_{i}italic_s start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT as S⁒Ri𝑆subscript𝑅𝑖SR_{i}italic_S italic_R start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT, which is a disk of radius rπ‘Ÿritalic_r and the corresponding communication range is modeled as C⁒Ri𝐢subscript𝑅𝑖CR_{i}italic_C italic_R start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT (a disk of radius R𝑅Ritalic_R).

III Proposed Strategy

In this section, we discuss the proposed sensor deployment strategy. As stated earlier, we model this scenario as a kπ‘˜kitalic_k-coverage problem, where we assume that the sensors are placed in the environment as a regular polygon tiling. Note that, this is not an ideal solution if the sensors are randomly deployed and thus, we propose to place the sensors strategically, such that the resulting cover set offers kπ‘˜kitalic_k-coverage with fewer sensors. Moreover, as shown in Fig.1, now we examine and assess a collection of planar convex regular polygons, which possess the ability to tile the Euclidean plane.

Now, examining these planar convex tiles is crucial to understand the kπ‘˜kitalic_k-coverage issue in planar WSNs, where Fig. 2 demonstrates a few convex regular polygons. Specifically, Fig. 2(a) represents an equilateral triangle with side rπ‘Ÿritalic_r, i.e., the area of this triangle is 34⁒r234superscriptπ‘Ÿ2\frac{\sqrt{3}}{4}r^{2}divide start_ARG square-root start_ARG 3 end_ARG end_ARG start_ARG 4 end_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT and Fig. 2(b) represents a regular hexagon Hr/2subscriptπ»π‘Ÿ2H_{r/2}italic_H start_POSTSUBSCRIPT italic_r / 2 end_POSTSUBSCRIPT of side r/2π‘Ÿ2r/2italic_r / 2. Since Hr/2subscriptπ»π‘Ÿ2H_{r/2}italic_H start_POSTSUBSCRIPT italic_r / 2 end_POSTSUBSCRIPT consists of six equilateral triangles with side r/2π‘Ÿ2r/2italic_r / 2, its area is obtained as

A⁒(Hr/2)=6Γ—34⁒(r/2)2=3⁒38⁒r2.𝐴subscriptπ»π‘Ÿ2634superscriptπ‘Ÿ22338superscriptπ‘Ÿ2A(H_{r/2})=6\times\frac{\sqrt{3}}{4}(r/2)^{2}=\frac{3\sqrt{3}}{8}r^{2}.italic_A ( italic_H start_POSTSUBSCRIPT italic_r / 2 end_POSTSUBSCRIPT ) = 6 Γ— divide start_ARG square-root start_ARG 3 end_ARG end_ARG start_ARG 4 end_ARG ( italic_r / 2 ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = divide start_ARG 3 square-root start_ARG 3 end_ARG end_ARG start_ARG 8 end_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT . (1)

Fig. 2(c) represents a regular hexagon Hrsubscriptπ»π‘ŸH_{r}italic_H start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT with side rπ‘Ÿritalic_r. Later, we demonstrate that Hrsubscriptπ»π‘ŸH_{r}italic_H start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT can contain a maximum of three non-overlapping Hr/2subscriptπ»π‘Ÿ2H_{r/2}italic_H start_POSTSUBSCRIPT italic_r / 2 end_POSTSUBSCRIPT. Finally, Fig. 2(d) illustrates a Hrsubscriptπ»π‘ŸH_{r}italic_H start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT, which consists of six equilateral triangles of side rπ‘Ÿritalic_r, i.e., regular polygon tiling of six equilateral triangles. Therefore, the area of Hrsubscriptπ»π‘ŸH_{r}italic_H start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT is

A⁒(Hr)=6Γ—34⁒(r)2=3⁒32⁒r2𝐴subscriptπ»π‘Ÿ634superscriptπ‘Ÿ2332superscriptπ‘Ÿ2A(H_{r})=6\times\frac{\sqrt{3}}{4}(r)^{2}=\frac{3\sqrt{3}}{2}r^{2}italic_A ( italic_H start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT ) = 6 Γ— divide start_ARG square-root start_ARG 3 end_ARG end_ARG start_ARG 4 end_ARG ( italic_r ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = divide start_ARG 3 square-root start_ARG 3 end_ARG end_ARG start_ARG 2 end_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT (2)

Though, from (1) and (2), we observe that A⁒(Hr)=4Γ—A⁒(Hr/2)𝐴subscriptπ»π‘Ÿ4𝐴subscriptπ»π‘Ÿ2A(H_{r})=4\times A(H_{r/2})italic_A ( italic_H start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT ) = 4 Γ— italic_A ( italic_H start_POSTSUBSCRIPT italic_r / 2 end_POSTSUBSCRIPT ), in the following lemma, we prove that more than three Hr/2subscriptπ»π‘Ÿ2H_{r/2}italic_H start_POSTSUBSCRIPT italic_r / 2 end_POSTSUBSCRIPT can not completely lie within a Hrsubscriptπ»π‘ŸH_{r}italic_H start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT.

Lemma 1.

A hexagon Hrsubscriptπ»π‘ŸH_{r}italic_H start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT of side rπ‘Ÿritalic_r can contain a maximum of three non-overlapping hexagons Hr/2subscriptπ»π‘Ÿ2H_{r/2}italic_H start_POSTSUBSCRIPT italic_r / 2 end_POSTSUBSCRIPT, each of side r/2π‘Ÿ2r/2italic_r / 2.

Proof.

The length of the largest diagonal of Hrsubscriptπ»π‘ŸH_{r}italic_H start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT and Hr/2subscriptπ»π‘Ÿ2H_{r/2}italic_H start_POSTSUBSCRIPT italic_r / 2 end_POSTSUBSCRIPT is 2⁒r2π‘Ÿ2r2 italic_r and rπ‘Ÿritalic_r, respectively, i.e., a single Hr/2subscriptπ»π‘Ÿ2H_{r/2}italic_H start_POSTSUBSCRIPT italic_r / 2 end_POSTSUBSCRIPT can completely lie within Hrsubscriptπ»π‘ŸH_{r}italic_H start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT. Hereafter, let Hr/21subscriptsuperscript𝐻1π‘Ÿ2H^{1}_{r/2}italic_H start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_r / 2 end_POSTSUBSCRIPT and Hr/22subscriptsuperscript𝐻2π‘Ÿ2H^{2}_{r/2}italic_H start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_r / 2 end_POSTSUBSCRIPT be two non-overlapping hexagons of side r/2π‘Ÿ2r/2italic_r / 2 each, and X=Hr/21βˆͺHr/22𝑋subscriptsuperscript𝐻1π‘Ÿ2subscriptsuperscript𝐻2π‘Ÿ2X=H^{1}_{r/2}\cup H^{2}_{r/2}italic_X = italic_H start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_r / 2 end_POSTSUBSCRIPT βˆͺ italic_H start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_r / 2 end_POSTSUBSCRIPT. Hence, the furthest distance between two points of X𝑋Xitalic_X is at least

l⁒(X)=2⁒(3⁒(r/2))2+(r/4)2=132⁒r<2⁒r.𝑙𝑋2superscript3π‘Ÿ22superscriptπ‘Ÿ42132π‘Ÿ2π‘Ÿl(X)=2\sqrt{(\sqrt{3}(r/2))^{2}+(r/4)^{2}}=\frac{\sqrt{13}}{2}r<2r.italic_l ( italic_X ) = 2 square-root start_ARG ( square-root start_ARG 3 end_ARG ( italic_r / 2 ) ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( italic_r / 4 ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG = divide start_ARG square-root start_ARG 13 end_ARG end_ARG start_ARG 2 end_ARG italic_r < 2 italic_r . (3)

Therefore, two Hr/2subscriptπ»π‘Ÿ2H_{r/2}italic_H start_POSTSUBSCRIPT italic_r / 2 end_POSTSUBSCRIPT can lie completely within Hrsubscriptπ»π‘ŸH_{r}italic_H start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT. Now, if we define X=⋃i=13Hr/2i𝑋superscriptsubscript𝑖13subscriptsuperscriptπ»π‘–π‘Ÿ2X=\bigcup\limits_{i=1}^{3}H^{i}_{r/2}italic_X = ⋃ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_H start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_r / 2 end_POSTSUBSCRIPT where each Hr/2isubscriptsuperscriptπ»π‘–π‘Ÿ2H^{i}_{r/2}italic_H start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_r / 2 end_POSTSUBSCRIPT βˆ€for-all\forallβˆ€ i𝑖iitalic_i is a hexagon of side r/2π‘Ÿ2r/2italic_r / 2, and all three are mutually non-overlapping to each other, the furthest distance, in this case, is at least

l⁒(X)=2⁒(3⁒(r/2))2+(r/4)2=132⁒r<2⁒r.𝑙𝑋2superscript3π‘Ÿ22superscriptπ‘Ÿ42132π‘Ÿ2π‘Ÿl(X)=2\sqrt{(\sqrt{3}(r/2))^{2}+(r/4)^{2}}=\frac{\sqrt{13}}{2}r<2r.italic_l ( italic_X ) = 2 square-root start_ARG ( square-root start_ARG 3 end_ARG ( italic_r / 2 ) ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( italic_r / 4 ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG = divide start_ARG square-root start_ARG 13 end_ARG end_ARG start_ARG 2 end_ARG italic_r < 2 italic_r . (4)

As shown in Fig. 2, this configuration also lies within Hrsubscriptπ»π‘ŸH_{r}italic_H start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT. Note that, we cannot place four Hr/2subscriptπ»π‘Ÿ2H_{r/2}italic_H start_POSTSUBSCRIPT italic_r / 2 end_POSTSUBSCRIPT inside Hrsubscriptπ»π‘ŸH_{r}italic_H start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT, where all of them are mutually non-overlapping to each other. The reason for this is, that in this case, l⁒(X)=r+r+12⁒r=52⁒r>2⁒rπ‘™π‘‹π‘Ÿπ‘Ÿ12π‘Ÿ52π‘Ÿ2π‘Ÿl(X)=r+r+\frac{1}{2}r=\frac{5}{2}r>2ritalic_l ( italic_X ) = italic_r + italic_r + divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_r = divide start_ARG 5 end_ARG start_ARG 2 end_ARG italic_r > 2 italic_r, which follows from Fig. 2. Since we have l⁒(X)>2⁒r𝑙𝑋2π‘Ÿl(X)>2ritalic_l ( italic_X ) > 2 italic_r, we claim that four Hr/2subscriptπ»π‘Ÿ2H_{r/2}italic_H start_POSTSUBSCRIPT italic_r / 2 end_POSTSUBSCRIPT cannot completely lie within Hrsubscriptπ»π‘ŸH_{r}italic_H start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT. ∎

Refer to caption
Figure 2: Different shape tiles
Lemma 2.

If sensor sisubscript𝑠𝑖s_{i}italic_s start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT with sensing range S⁒Ri𝑆subscript𝑅𝑖SR_{i}italic_S italic_R start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT is placed on a vertex of an equilateral triangle of side rπ‘Ÿritalic_r, the entire triangle lies within S⁒Ri𝑆subscript𝑅𝑖SR_{i}italic_S italic_R start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT.

Proof.

Let Δ⁒A⁒B⁒CΔ𝐴𝐡𝐢\Delta ABCroman_Ξ” italic_A italic_B italic_C be an equilateral triangle of side rπ‘Ÿritalic_r and a sensor sisubscript𝑠𝑖s_{i}italic_s start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT with sensing range S⁒Ri𝑆subscript𝑅𝑖SR_{i}italic_S italic_R start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT is placed on A𝐴Aitalic_A. If we draw a circle of radius rπ‘Ÿritalic_r at A𝐴Aitalic_A, the vertices B𝐡Bitalic_B and C𝐢Citalic_C lie within its circumference. Similarly, A,C𝐴𝐢A,Citalic_A , italic_C and A,B𝐴𝐡A,Bitalic_A , italic_B lie on the circles’ circumference, whose center is at B𝐡Bitalic_B and C𝐢Citalic_C, respectively. Since, the circles center at A,B𝐴𝐡A,Bitalic_A , italic_B, and C𝐢Citalic_C are convex, their intersecting region formed by A,B𝐴𝐡A,Bitalic_A , italic_B, and C𝐢Citalic_C is also a convex zone [8]. Moreover, Δ⁒A⁒B⁒CΔ𝐴𝐡𝐢\Delta ABCroman_Ξ” italic_A italic_B italic_C also lies on the convex zone formed by A,B,𝐴𝐡A,B,italic_A , italic_B , and C𝐢Citalic_C. Therefore, Δ⁒A⁒B⁒CΔ𝐴𝐡𝐢\Delta ABCroman_Ξ” italic_A italic_B italic_C lies in S⁒Ri𝑆subscript𝑅𝑖SR_{i}italic_S italic_R start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT. ∎

Note that, for sisubscript𝑠𝑖s_{i}italic_s start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT, both S⁒Ri𝑆subscript𝑅𝑖SR_{i}italic_S italic_R start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT and C⁒Ri𝐢subscript𝑅𝑖CR_{i}italic_C italic_R start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT are circular disks of a certain radius. But we can not cover an Euclidean plane with circular tiles. It is well-studied that a hexagonal tile can cover the greatest area inside a circular disc [7]. As a result, we use the hexagonal tiles to cover the Euclidean plane.

III-A Sensor Deployment Strategy

As both S⁒Ri𝑆subscript𝑅𝑖SR_{i}italic_S italic_R start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT and C⁒Ri𝐢subscript𝑅𝑖CR_{i}italic_C italic_R start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT can be approximated as hexagonal tiles, the sensors can cover the Euclidean plane in a way similar to these tiles. Moreover, here we propose a sensor deployment strategy that enables us to obtain kπ‘˜kitalic_k-coverage, but with fewer sensors. Notably, we strategically place the sensor inside a regular hexagon, and by employing this method for every hexagon, we can cover the entire Euclidean plane. In this context, we classify our strategy into four separate cases, as described below.

Refer to caption
Figure 3: Sensor Deployment Strategy

A-I: For k=1π‘˜1k=1italic_k = 1 (1111-Coverage)

In this case, as shown in Fig.3(a), we deploy the sensor at the center O𝑂Oitalic_O of the hexagon. Note that, as O𝑂Oitalic_O is the common vertex of all the six equilateral triangles contained in the hexagon, this deployment strategy is adequate for providing 1111-coverage to these triangles, which also follows from Lemma 2. Note that, kπ‘˜kitalic_k-coverage can easily be obtained by placing kπ‘˜kitalic_k co-located sensors at O𝑂Oitalic_O.

Before discussing the higher coverage sensor deployment strategy for a hexagon, we demonstrate in the ensuing lemma that, to obtain 2222-coverage, where co-located sensors are prohibited, at least three additional sensors are required.

Lemma 3.

If placing a sensor in the center O𝑂Oitalic_O is prohibited, at least three sensors are required to cover a hexagon.

Proof.

A hexagon consists of six equilateral triangles and two consecutive triangles within the hexagon share a common edge. Moreover, from Lemma 2, it follows that a point on or within a triangle can cover the entire triangle. It is evident that, the entire triangle cannot be covered by a point lying outside of it. Therefore, if we deploy a sensor on the common edge of two triangles it can cover both of them. Moreover, as a hexagon contains six triangles, we can assert that at least three sensors are necessary for covering the hexagon. ∎

A-II: For k=2π‘˜2k=2italic_k = 2 (2222-Coverage)

As shown in Fig.3(b), here also, there are six equilateral triangles within the hexagon. Consequently, to achieve 2222-coverage of all these triangles, we need to deploy the sensors in a way such that each triangle lies inside the coverage area of at least two sensors. From the figure, we observe that O𝑂Oitalic_O is the center, and A,B,C,D,E,F𝐴𝐡𝐢𝐷𝐸𝐹A,B,C,D,E,Fitalic_A , italic_B , italic_C , italic_D , italic_E , italic_F are the vertices of the hexagon. As discussed previously, deploying a sensor at O𝑂Oitalic_O provides 1111-coverage for all the triangles. Therefore, if the sensors are deployed on the alternate vertices, i.e., at B,D,𝐡𝐷B,D,italic_B , italic_D , and F𝐹Fitalic_F; or at A,C,𝐴𝐢A,C,italic_A , italic_C , and E𝐸Eitalic_E, this, along with the sensor placed at O𝑂Oitalic_O, provides 2222-coverage to all the equilateral triangles.

A-III: For k=3π‘˜3k=3italic_k = 3 (3333-Coverage)

Unlike the previous case, here we strategically deploy the sensors at all the vertices of the hexagon, i.e., A,B,C,D,E,𝐴𝐡𝐢𝐷𝐸A,B,C,D,E,italic_A , italic_B , italic_C , italic_D , italic_E , and F𝐹Fitalic_F, and also at the center O𝑂Oitalic_O. As a result of which, Fig.3(c) demonstrates that all the triangles lie inside the coverage area of at least 3333 sensors, i.e., a 3333-coverage scenario.

Note that, up to the case of k=3π‘˜3k=3italic_k = 3-coverage, we deployed the sensor on the vertices and at the center of the hexagon. However, we can not place more than one sensor at the same place. Therefore, for the case of k>3π‘˜3k>3italic_k > 3-coverage scenario, we strategically place the sensors inside the regular hexagon, as discussed next.

A-IV: For k>3π‘˜3k>3italic_k > 3 (k>3π‘˜3k>3italic_k > 3-Coverage)

In Fig.3(d), we observe that A,B,C,D,E,𝐴𝐡𝐢𝐷𝐸A,B,C,D,E,italic_A , italic_B , italic_C , italic_D , italic_E , and F𝐹Fitalic_F are the vertices of the hexagon. Also, O⁒A𝑂𝐴OAitalic_O italic_A and O⁒B𝑂𝐡OBitalic_O italic_B are the sides of Δ⁒O⁒A⁒BΔ𝑂𝐴𝐡\Delta OABroman_Ξ” italic_O italic_A italic_B, O⁒C𝑂𝐢OCitalic_O italic_C and O⁒D𝑂𝐷ODitalic_O italic_D are the sides of Δ⁒O⁒C⁒DΔ𝑂𝐢𝐷\Delta OCDroman_Ξ” italic_O italic_C italic_D, and O⁒E𝑂𝐸OEitalic_O italic_E, O⁒F𝑂𝐹OFitalic_O italic_F are the sides of Δ⁒O⁒E⁒FΔ𝑂𝐸𝐹\Delta OEFroman_Ξ” italic_O italic_E italic_F, respectively. Accordingly, we rename O⁒A,O⁒B,O⁒C,O⁒D,O⁒E𝑂𝐴𝑂𝐡𝑂𝐢𝑂𝐷𝑂𝐸OA,OB,OC,OD,OEitalic_O italic_A , italic_O italic_B , italic_O italic_C , italic_O italic_D , italic_O italic_E, and O⁒F𝑂𝐹OFitalic_O italic_F as l1,β‹―,l6subscript𝑙1β‹―subscript𝑙6l_{1},\cdots,l_{6}italic_l start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , β‹― , italic_l start_POSTSUBSCRIPT 6 end_POSTSUBSCRIPT. Therefore, to get kπ‘˜kitalic_k-coverage (k>3)π‘˜3(k>3)( italic_k > 3 ) from (kβˆ’1)π‘˜1(k-1)( italic_k - 1 )-coverage, we deploy 3333 additional sensors strategically on li,i=1,β‹―,6formulae-sequencesubscript𝑙𝑖𝑖1β‹―6l_{i},i=1,\cdots,6italic_l start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , italic_i = 1 , β‹― , 6. Specifically, for an even (odd) kπ‘˜kitalic_k, we deploy the sensors on those lisubscript𝑙𝑖l_{i}italic_l start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT, whose index is odd (even). For example, for k=4π‘˜4k=4italic_k = 4, we deploy three sensors on l1,l3subscript𝑙1subscript𝑙3l_{1},l_{3}italic_l start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_l start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT and l5subscript𝑙5l_{5}italic_l start_POSTSUBSCRIPT 5 end_POSTSUBSCRIPT, respectively, i.e., on O⁒A,O⁒C𝑂𝐴𝑂𝐢OA,OCitalic_O italic_A , italic_O italic_C and O⁒D𝑂𝐷ODitalic_O italic_D which is shown in Fig.3 (d). Similarly, for k=5π‘˜5k=5italic_k = 5, we deploy three sensors on l2,l4subscript𝑙2subscript𝑙4l_{2},l_{4}italic_l start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , italic_l start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT and l6subscript𝑙6l_{6}italic_l start_POSTSUBSCRIPT 6 end_POSTSUBSCRIPT, respectively, i.e., on O⁒B,O⁒D𝑂𝐡𝑂𝐷OB,ODitalic_O italic_B , italic_O italic_D and O⁒F𝑂𝐹OFitalic_O italic_F. Note that, we can deploy the sensor at any point on the line but doing so at the same location is not allowed.

III-B Required Number of Sensors for kπ‘˜kitalic_k-coverage

In Fig.3(a), we observed that 1111-coverage can be achieved with a single sensor placed at the center of the hexagon. Note that, there are six triangles inside the hexagon and any two consecutive triangles share a common side. Consequently, it is beneficial for both triangles to achieve kπ‘˜kitalic_k-coverage if we place sensors on their common sides. Therefore, four sensors are sufficient for 2222-coverage and eventually, for 3333-coverage, 7777 sensors are adequate. As a result, we can state, that if n1subscript𝑛1n_{1}italic_n start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT sensors are required for kπ‘˜kitalic_k-coverage for k>3π‘˜3k>3italic_k > 3, then n1+3subscript𝑛13n_{1}+3italic_n start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + 3 sensors are adequate for (k+1)π‘˜1(k+1)( italic_k + 1 )-coverage. Therefore, we can express this as an arithmetic progression, whose first term is 1111 and the common difference is 3333. Hence, for the general kπ‘˜kitalic_k-coverage, the required number of sensors are

n1=1+3⁒(kβˆ’1)=3⁒kβˆ’2.subscript𝑛113π‘˜13π‘˜2n_{1}=1+3(k-1)=3k-2.italic_n start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = 1 + 3 ( italic_k - 1 ) = 3 italic_k - 2 . (5)

III-C Sensor Density per Unit Area

Concerning our proposed deployment strategy, here we calculate the respective sensor density per unit area. From (2), we have the area of the hexagon of side rπ‘Ÿritalic_r to be 3⁒32⁒r2332superscriptπ‘Ÿ2\frac{3\sqrt{3}}{2}r^{2}divide start_ARG 3 square-root start_ARG 3 end_ARG end_ARG start_ARG 2 end_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT. Furthermore, from (5), we claim that, 3⁒kβˆ’23π‘˜23k-23 italic_k - 2 sensors are sufficient to support kπ‘˜kitalic_k-coverage within the hexagon. Hence, if sdsubscript𝑠𝑑s_{d}italic_s start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT is the sensor density per unit area, we obtain

sd=3⁒kβˆ’23⁒32⁒r2=2⁒(3⁒kβˆ’2)3⁒3⁒r2.subscript𝑠𝑑3π‘˜2332superscriptπ‘Ÿ223π‘˜233superscriptπ‘Ÿ2s_{d}=\dfrac{3k-2}{\frac{3\sqrt{3}}{2}r^{2}}=\dfrac{2(3k-2)}{3\sqrt{3}r^{2}}.italic_s start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT = divide start_ARG 3 italic_k - 2 end_ARG start_ARG divide start_ARG 3 square-root start_ARG 3 end_ARG end_ARG start_ARG 2 end_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG = divide start_ARG 2 ( 3 italic_k - 2 ) end_ARG start_ARG 3 square-root start_ARG 3 end_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG . (6)

For example, when k=1π‘˜1k=1italic_k = 1, we get sd=23⁒3⁒r2subscript𝑠𝑑233superscriptπ‘Ÿ2s_{d}=\frac{2}{3\sqrt{3}r^{2}}italic_s start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT = divide start_ARG 2 end_ARG start_ARG 3 square-root start_ARG 3 end_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG.

On the other hand, the authors in [7] have discussed another sensor deployment strategy for an Euclidean plane. They used regular hexagons of side r/2π‘Ÿ2r/2italic_r / 2 to cover an Euclidean plane. As a result, for kπ‘˜kitalic_k-coverage, they randomly placed kπ‘˜kitalic_k sensors, each having sensing range rπ‘Ÿritalic_r, inside the hexagon of side r/2π‘Ÿ2r/2italic_r / 2. In their proposed strategy, the respective sensor density per unit area is

sdex=k3⁒32⁒(r/2)2=8⁒k3⁒3⁒r2.superscriptsubscript𝑠𝑑exπ‘˜332superscriptπ‘Ÿ228π‘˜33superscriptπ‘Ÿ2s_{d}^{\rm ex}=\frac{k}{\frac{3\sqrt{3}}{2}(r/2)^{2}}=\frac{8k}{3\sqrt{3}r^{2}}.italic_s start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_ex end_POSTSUPERSCRIPT = divide start_ARG italic_k end_ARG start_ARG divide start_ARG 3 square-root start_ARG 3 end_ARG end_ARG start_ARG 2 end_ARG ( italic_r / 2 ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG = divide start_ARG 8 italic_k end_ARG start_ARG 3 square-root start_ARG 3 end_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG . (7)

Accordingly, by using (6) and (7), we quantify the gain of our proposed strategy as

sdexβˆ’sd=8⁒k3⁒3⁒r2βˆ’2⁒(3⁒kβˆ’2)3⁒3⁒r2=2⁒k+43⁒3⁒r2.superscriptsubscript𝑠𝑑exsubscript𝑠𝑑8π‘˜33superscriptπ‘Ÿ223π‘˜233superscriptπ‘Ÿ22π‘˜433superscriptπ‘Ÿ2s_{d}^{\rm ex}-s_{d}=\frac{8k}{3\sqrt{3}r^{2}}-\frac{2(3k-2)}{3\sqrt{3}r^{2}}=% \frac{2k+4}{3\sqrt{3}r^{2}}.italic_s start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_ex end_POSTSUPERSCRIPT - italic_s start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT = divide start_ARG 8 italic_k end_ARG start_ARG 3 square-root start_ARG 3 end_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG - divide start_ARG 2 ( 3 italic_k - 2 ) end_ARG start_ARG 3 square-root start_ARG 3 end_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG = divide start_ARG 2 italic_k + 4 end_ARG start_ARG 3 square-root start_ARG 3 end_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG . (8)

Thus, we observe that for any positive kπ‘˜kitalic_k, sdexβˆ’sd>0superscriptsubscript𝑠𝑑exsubscript𝑠𝑑0s_{d}^{\rm ex}-s_{d}>0italic_s start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_ex end_POSTSUPERSCRIPT - italic_s start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT > 0, i.e., our strategy outperforms the existing benchmark. Another interesting insight, which can be obtained is

limkβ†’βˆžsdsdex=34,subscriptβ†’π‘˜subscript𝑠𝑑superscriptsubscript𝑠𝑑ex34\lim_{k\rightarrow\infty}\frac{s_{d}}{s_{d}^{\rm ex}}=\frac{3}{4},roman_lim start_POSTSUBSCRIPT italic_k β†’ ∞ end_POSTSUBSCRIPT divide start_ARG italic_s start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT end_ARG start_ARG italic_s start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_ex end_POSTSUPERSCRIPT end_ARG = divide start_ARG 3 end_ARG start_ARG 4 end_ARG , (9)

which demonstrates the advantage of our strategy as kβ†’βˆžβ†’π‘˜k\rightarrow\inftyitalic_k β†’ ∞.

III-D Required Number of Sensor for a Plane

In Section III-A, we discussed our sensor deployment strategy for a single hexagon with side rπ‘Ÿritalic_r. As shown in Fig.1, the entire plane can be tiled using the regular hexagon of side rπ‘Ÿritalic_r. So, we looked into the required number of sensors for kπ‘˜kitalic_k-coverage of the entire plane.

To achieve this, we look into a model with one central hexagon and six more hexagons around it. Each of these second β€˜layer’ of hexagons shares an edge with the six edges of the central one. By continuing with this procedure, we obtain the third, fourth, and consecutive layers of hexagons, with respect to the central one. Let l𝑙litalic_l denote the number of layers that cover the plane. We name this β€˜the solar model’, as is shown in Fig. 1, where O𝑂Oitalic_O denotes the center of the central hexagon. Accordingly, we obtain a closed-form analytical expression for the required number of sensors to support kπ‘˜kitalic_k-coverage upto l𝑙litalic_l layers. In this context, we classify our strategy into four separate cases, as described below.

D-I: For k=1π‘˜1k=1italic_k = 1 (1111-Coverage)

From previous discussions, we observe that in 1111-coverage a single sensor lies at the centre of the hexagon. Thus, the required number of sensors, in this case, is equal to the total number of hexagons in the solar model. Therefore, from Fig.1, we observe that a single hexagon lies in the first layer, 6666 hexagons lie in the second layer, and in general, 6⁒(lβˆ’1)6𝑙16(l-1)6 ( italic_l - 1 ) hexagons lie in the l𝑙litalic_l-th layer, where lβ‰₯2𝑙2l\geq 2italic_l β‰₯ 2. Specifically, the number of hexagons lie in the l𝑙litalic_l-th layer for lβ‰₯2𝑙2l\geq 2italic_l β‰₯ 2, follows an arithmetic progression whose first term as well as common difference is 6666. Hence the total number of hexagons that lie within the solar model is

1+6+12+β‹―+6⁒(lβˆ’1)=1+3⁒l⁒(lβˆ’1).1612β‹―6𝑙113𝑙𝑙1\displaystyle 1+6+12+\cdots+6(l-1)=1+3l(l-1).1 + 6 + 12 + β‹― + 6 ( italic_l - 1 ) = 1 + 3 italic_l ( italic_l - 1 ) . (10)

Therefore, in this case, n⁒(l,k=1)π‘›π‘™π‘˜1n(l,k=1)italic_n ( italic_l , italic_k = 1 ), the required number of sensors to get 1111-coverage of l𝑙litalic_l layers is given by

n(l,k=1)=1+3l(lβˆ’1)).n(l,k=1)=1+3l(l-1)).italic_n ( italic_l , italic_k = 1 ) = 1 + 3 italic_l ( italic_l - 1 ) ) . (11)

D-II: For k=2π‘˜2k=2italic_k = 2 (2222-Coverage)

To achieve 2222-coverage, we consider that n⁒(l,k=1)π‘›π‘™π‘˜1n(l,k=1)italic_n ( italic_l , italic_k = 1 ) sensors are already deployed to support 1111-coverage. Section III-A shows that to support 2222-coverage, we need three more sensors in the central hexagon, placed at the alternate vertices, apart from the one at the center. Similarly, for the second layer of hexagons, we need 9999 sensors and for the third layer, 15151515 sensors are required. Moving forward, we need 3+6⁒(lβˆ’1)36𝑙13+6(l-1)3 + 6 ( italic_l - 1 ) sensors for the l𝑙litalic_l-th layer. Therefore, we obtain the total number of sensors required as

n⁒(l,k=2)π‘›π‘™π‘˜2\displaystyle n(l,k=2)italic_n ( italic_l , italic_k = 2 ) =n(l,k=1)+(3+9+β‹―+(3+6(lβˆ’1))\displaystyle=n(l,k=1)+(3+9+\cdots+(3+6(l-1))= italic_n ( italic_l , italic_k = 1 ) + ( 3 + 9 + β‹― + ( 3 + 6 ( italic_l - 1 ) )
=n⁒(l,k=1)+l2⁒(2Γ—3+(lβˆ’1)⁒6)absentπ‘›π‘™π‘˜1𝑙223𝑙16\displaystyle=n(l,k=1)+\frac{l}{2}\left(2\times 3+(l-1)6\right)= italic_n ( italic_l , italic_k = 1 ) + divide start_ARG italic_l end_ARG start_ARG 2 end_ARG ( 2 Γ— 3 + ( italic_l - 1 ) 6 )
=(a)1+3l(lβˆ’1))+3l2=6l2βˆ’3l+1,\displaystyle\overset{(a)}{=}1+3l(l-1))+3l^{2}=6l^{2}-3l+1,start_OVERACCENT ( italic_a ) end_OVERACCENT start_ARG = end_ARG 1 + 3 italic_l ( italic_l - 1 ) ) + 3 italic_l start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = 6 italic_l start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 3 italic_l + 1 , (12)

where (a)π‘Ž(a)( italic_a ) follows from (11).

D-III: For k=3π‘˜3k=3italic_k = 3 (3333-Coverage)

Similar to the previous discussion, we observed that having 3333-coverage implies placing sensors on all the vertices and at the center of the hexagon. Alternatively, for 3333-coverage scenario up to l𝑙litalic_l layers, as discussed earlier, we consider that n⁒(l,k=2)π‘›π‘™π‘˜2n(l,k=2)italic_n ( italic_l , italic_k = 2 ) sensors are required to support 2222-coverage. Note that if we place the sensors on the remaining sensor-free vertices, the l𝑙litalic_l layers achieve the 3333-coverage. Here, the number of remaining sensor-free vertices is

l2⁒(2Γ—3+(lβˆ’1)⁒6)=3⁒l2.𝑙223𝑙163superscript𝑙2\frac{l}{2}\left(2\times 3+(l-1)6\right)=3l^{2}.divide start_ARG italic_l end_ARG start_ARG 2 end_ARG ( 2 Γ— 3 + ( italic_l - 1 ) 6 ) = 3 italic_l start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT . (13)

This can be proved in a way similar to the process as discussed for the case of 2222-coverage. Thus, the required number of sensors, in this case, is

n⁒(l,k=3)=n⁒(l,k=2)+3⁒l2=9⁒l2βˆ’3⁒l+1.π‘›π‘™π‘˜3π‘›π‘™π‘˜23superscript𝑙29superscript𝑙23𝑙1n(l,k=3)=n(l,k=2)+3l^{2}=9l^{2}-3l+1.italic_n ( italic_l , italic_k = 3 ) = italic_n ( italic_l , italic_k = 2 ) + 3 italic_l start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = 9 italic_l start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 3 italic_l + 1 . (14)

D-IV: For k>3π‘˜3k>3italic_k > 3 (k>3π‘˜3k>3italic_k > 3-Coverage)

Up to 3333-coverage, sensors are strategically deployed on the the vertices, and at the center of the hexagons. But, after 3333-coverage, they are deployed inside the hexagons, and to increase one more coverage, three more sensors are needed for every single hexagon. Hence, for each k>3π‘˜3k>3italic_k > 3, to get kπ‘˜kitalic_k-coverage, 3Γ—(1+3⁒l⁒(lβˆ’1))313𝑙𝑙13\times(1+3l(l-1))3 Γ— ( 1 + 3 italic_l ( italic_l - 1 ) ) more sensors are needed as compared to n⁒(l,kβˆ’1)π‘›π‘™π‘˜1n(l,k-1)italic_n ( italic_l , italic_k - 1 ). Therefore, the total number of sensors required is

n⁒(l,k)π‘›π‘™π‘˜\displaystyle n(l,k)italic_n ( italic_l , italic_k ) =n⁒(l,k=3)+(kβˆ’3)⁒(3Γ—(1+3⁒l⁒(lβˆ’1)))absentπ‘›π‘™π‘˜3π‘˜3313𝑙𝑙1\displaystyle=n(l,k=3)+(k-3)(3\times(1+3l(l-1)))= italic_n ( italic_l , italic_k = 3 ) + ( italic_k - 3 ) ( 3 Γ— ( 1 + 3 italic_l ( italic_l - 1 ) ) )
=9⁒(kβˆ’2)⁒l2βˆ’3⁒(3⁒kβˆ’8)⁒l+3⁒kβˆ’8absent9π‘˜2superscript𝑙233π‘˜8𝑙3π‘˜8\displaystyle=9(k-2)l^{2}-3(3k-8)l+3k-8= 9 ( italic_k - 2 ) italic_l start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 3 ( 3 italic_k - 8 ) italic_l + 3 italic_k - 8 (15)

By combining (11), (III), (14), and (III), we obtain an unified expression for n⁒(l,k)π‘›π‘™π‘˜n(l,k)italic_n ( italic_l , italic_k ) as

n⁒(l,k)={1+3⁒l⁒(lβˆ’1)k=1,6⁒l2βˆ’3⁒l+1k=2,9⁒(kβˆ’2)⁒l2βˆ’3⁒(3⁒kβˆ’8)⁒l+3⁒kβˆ’8kβ‰₯3.π‘›π‘™π‘˜cases13𝑙𝑙1π‘˜16superscript𝑙23𝑙1π‘˜29π‘˜2superscript𝑙233π‘˜8𝑙3π‘˜8π‘˜3\!\!\!\!n(l,k)=\begin{cases}1+3l(l-1)&k=1,\\ 6l^{2}-3l+1&k=2,\\ 9(k-2)l^{2}-3(3k-8)l+3k-8&k\geq 3.\end{cases}italic_n ( italic_l , italic_k ) = { start_ROW start_CELL 1 + 3 italic_l ( italic_l - 1 ) end_CELL start_CELL italic_k = 1 , end_CELL end_ROW start_ROW start_CELL 6 italic_l start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 3 italic_l + 1 end_CELL start_CELL italic_k = 2 , end_CELL end_ROW start_ROW start_CELL 9 ( italic_k - 2 ) italic_l start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 3 ( 3 italic_k - 8 ) italic_l + 3 italic_k - 8 end_CELL start_CELL italic_k β‰₯ 3 . end_CELL end_ROW (16)

On the other hand, by using Lemma 1 and the existing sensor deployment strategy as discussed in [7], the associated total number of sensors required to get kπ‘˜kitalic_k-coverage for l𝑙litalic_l layers is

nex⁒(l,k)={1+3⁒l⁒(lβˆ’1)k=1,k⁒(15⁒l2βˆ’27⁒l+18)kβ‰₯2.superscript𝑛exπ‘™π‘˜cases13𝑙𝑙1π‘˜1π‘˜15superscript𝑙227𝑙18π‘˜2\!\!\!\!n^{\rm ex}(l,k)=\begin{cases}1+3l(l-1)&k=1,\\ k(15l^{2}-27l+18)&k\geq 2.\end{cases}italic_n start_POSTSUPERSCRIPT roman_ex end_POSTSUPERSCRIPT ( italic_l , italic_k ) = { start_ROW start_CELL 1 + 3 italic_l ( italic_l - 1 ) end_CELL start_CELL italic_k = 1 , end_CELL end_ROW start_ROW start_CELL italic_k ( 15 italic_l start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 27 italic_l + 18 ) end_CELL start_CELL italic_k β‰₯ 2 . end_CELL end_ROW (17)

The expression for kβ‰₯2π‘˜2k\geq 2italic_k β‰₯ 2 is obtained based on the fact that 15⁒l2βˆ’27⁒l+1815superscript𝑙227𝑙1815l^{2}-27l+1815 italic_l start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 27 italic_l + 18 is the total number of Hr/2subscriptπ»π‘Ÿ2H_{r/2}italic_H start_POSTSUBSCRIPT italic_r / 2 end_POSTSUBSCRIPT completely lies inside the solar model of l𝑙litalic_l layer where each Hr/2subscriptπ»π‘Ÿ2H_{r/2}italic_H start_POSTSUBSCRIPT italic_r / 2 end_POSTSUBSCRIPT contains kπ‘˜kitalic_k sensors. It is interesting to observe, that irrespective of kπ‘˜kitalic_k and l𝑙litalic_l, we always have nex⁒(l,k)βˆ’n⁒(l,k)β‰₯0superscript𝑛exπ‘™π‘˜π‘›π‘™π‘˜0n^{\rm ex}(l,k)-n(l,k)\geq 0italic_n start_POSTSUPERSCRIPT roman_ex end_POSTSUPERSCRIPT ( italic_l , italic_k ) - italic_n ( italic_l , italic_k ) β‰₯ 0, which demonstrates the benefit of our proposed strategy. Another interesting insight, which can be obtained is

limkβ†’βˆžlimlβ†’βˆžn⁒(l,k)nex⁒(l,k)=35,subscriptβ†’π‘˜subscriptβ†’π‘™π‘›π‘™π‘˜superscript𝑛exπ‘™π‘˜35\lim_{k\rightarrow\infty}\lim_{l\rightarrow\infty}\frac{n(l,k)}{n^{\rm ex}(l,k% )}=\frac{3}{5},roman_lim start_POSTSUBSCRIPT italic_k β†’ ∞ end_POSTSUBSCRIPT roman_lim start_POSTSUBSCRIPT italic_l β†’ ∞ end_POSTSUBSCRIPT divide start_ARG italic_n ( italic_l , italic_k ) end_ARG start_ARG italic_n start_POSTSUPERSCRIPT roman_ex end_POSTSUPERSCRIPT ( italic_l , italic_k ) end_ARG = divide start_ARG 3 end_ARG start_ARG 5 end_ARG , (18)

which demonstrates the advantage of our strategy as k,lβ†’βˆžβ†’π‘˜π‘™k,l\rightarrow\inftyitalic_k , italic_l β†’ ∞.

IV Numerical Results

In this section, we evaluate the proposed sensor deployment strategy and also, compare its performance with that of the existing benchmark scheme in [7]. As stated earlier, the work in [7] follows a computation geometry-based approach (CGA), where the authors suggest to deply the sensors in a certain manner in order to attain kπ‘˜kitalic_k-coverage.

Refer to caption
Figure 4: Impact of sensing radius on sensor density.

Fig. 4 depicts the impact of the sensing radius on the deployed sensor density. Here we demonstrate the effect of both the proposed strategy and CGA specifically to attain kπ‘˜kitalic_k-coverage for k=2,7π‘˜27k=2,7italic_k = 2 , 7. We observe that while the sensor density decreases monotonically with the sensing radius irrespective of the approach being considered, the proposed strategy results in significantly lower sensor density, when compared against CGA. The reason for this is attributed to the intelligent placing of sensors on the vertices and edges of the hexagons.

Fig. 5 illustrates the variation of sensor density with the kπ‘˜kitalic_k-coverage, for both the cases of r=10π‘Ÿ10r=10italic_r = 10 m and r=20π‘Ÿ20r=20italic_r = 20 m, respectively. While we observe that in both the cases, the sensor density increases with the value of kπ‘˜kitalic_k in kπ‘˜kitalic_k-coverage, the reason is intuitive. A higher value of kπ‘˜kitalic_k implies that, at any point of time, any particular area is being covered by kπ‘˜kitalic_k sensors. However, at the same time, we observe that with identical system parameters, our proposed approach results to a lower sensor density as compared to CGA.

Refer to caption
Figure 5: Impact of kπ‘˜kitalic_k-coverage on sensor density.

For a given number of layers l𝑙litalic_l and kπ‘˜kitalic_k-coverage, Section III-D obtains a closed-form expression for the required number of sensors n⁒(l,k)π‘›π‘™π‘˜n(l,k)italic_n ( italic_l , italic_k ). Accordingly, Fig. 6 investigates the impact of l𝑙litalic_l on n⁒(l,k)π‘›π‘™π‘˜n(l,k)italic_n ( italic_l , italic_k ) for k=3π‘˜3k=3italic_k = 3 and k=10π‘˜10k=10italic_k = 10. Similarly, Fig. 7 illustrates the impact of kπ‘˜kitalic_k on n⁒(l,k)π‘›π‘™π‘˜n(l,k)italic_n ( italic_l , italic_k ) for l=3𝑙3l=3italic_l = 3 and l=5𝑙5l=5italic_l = 5, respectively. In both these figures, we observe an increasing trend irrespective of the scheme employed, which is intuitive. Moreover, they also demonstrate the enhanced performance of the proposed sensor deployment scheme against CGA. As explained earlier, here also, this performance gain is attributed to the intelligent placement of sensors to obtain the desired coverage with a certain number of layers.

Refer to caption
Figure 6: Impact of layers on required number of sensors.

Finally, Fig. 8 demonstrates the joint effect of coverage and number of layers on the number of sensors required. Specifically, by varying both l𝑙litalic_l and kπ‘˜kitalic_k, we calculate both n⁒(l,k)π‘›π‘™π‘˜n(l,k)italic_n ( italic_l , italic_k ) and nex⁒(l,k)superscript𝑛exπ‘™π‘˜n^{\rm ex}(l,k)italic_n start_POSTSUPERSCRIPT roman_ex end_POSTSUPERSCRIPT ( italic_l , italic_k ) according to our proposed strategy and CGA, respectively. Henceforth, we plot nex⁒(l,k)βˆ’n⁒(l,k)superscript𝑛exπ‘™π‘˜π‘›π‘™π‘˜n^{\rm ex}(l,k)-n(l,k)italic_n start_POSTSUPERSCRIPT roman_ex end_POSTSUPERSCRIPT ( italic_l , italic_k ) - italic_n ( italic_l , italic_k ) by jointly varying both the parameters. We observe from the figure that this quantity is always positive in nature, i.e., our deployment strategy requires lesser sensors compared to CGA with identical scenarios. Moreover, the variation with l𝑙litalic_l for a given kπ‘˜kitalic_k is relatively sharp compared to the variation with kπ‘˜kitalic_k for a given l𝑙litalic_l. This is because both n⁒(l,k)π‘›π‘™π‘˜n(l,k)italic_n ( italic_l , italic_k ) and nex⁒(l,k)superscript𝑛exπ‘™π‘˜n^{\rm ex}(l,k)italic_n start_POSTSUPERSCRIPT roman_ex end_POSTSUPERSCRIPT ( italic_l , italic_k ) are linear and quadratic with respect to kπ‘˜kitalic_k and l𝑙litalic_l, respectively.

Refer to caption
Figure 7: Impact of kπ‘˜kitalic_k-coverage on required number of sensors.
Refer to caption
Figure 8: Joint effect of kπ‘˜kitalic_k-coverage and number of layers on required number of sensors.

V Conclusion

In this work, we investigated a graph theoretic sensor deployment approach in WSNs. Specifically, we modeled the sensor deployment problem as a kπ‘˜kitalic_k-coverage problem to propose a respective solution. Accordingly, we obtained analytical closed-form expressions of two performance metrics, namely β€˜sensor density per unit area’ and β€˜required number of sensors for a plane’ in terms of the sensing radius rπ‘Ÿritalic_r, coverage kπ‘˜kitalic_k, and level l𝑙litalic_l, respectively. Also, we evaluated these quantities for the existing benchmark schemes to quantify the advantages of our proposed strategy. We demonstrated that the required number of sensors is linear and quadratic with respect to kπ‘˜kitalic_k and l𝑙litalic_l, respectively. Finally, the numerical results validate the advantages of our proposed strategy. Based on the framework presented in this paper, an immediate extension of this work is to investigate the scenario, where tiling does not always necessarily happen by regular hexagons.

References

  • [1] Ericsson Mobility Report, November 2023, [Online]. Available: https://www.ericsson.com/en/reports-and-papers/mobility-report/reports.
  • [2] P. Mukherjee and S. De, β€œDynamic feedback based adaptive modulation for energy-efficient communication,” IEEE Commun. Lett., vol. 23, no. 5, pp. 946–949, May 2019.
  • [3] I. Muhammad, H. Alves, N. H. Mahmood, O. L. A. LΓ³pez, and M. Latva-aho, β€œMission effective capacityβ€”a novel dependability metric: A study case of multiconnectivity-enabled URLLC for IIoT,” IEEE Trans. Ind. Inf., vol. 18, no. 6, pp. 4180–4188, June 2022.
  • [4] Y. Wang, S. Wu, Z. Chen, X. Gao, and G. Chen, β€œCoverage problem with uncertain properties in wireless sensor networks: A survey,” Comput. Netw., vol. 123, pp. 200–232, Aug. 2017.
  • [5] Manju, S. Singh, S. Kumar, A. Nayyar, F. Al-Turjman, and L. Mostarda, β€œProficient QoS-based target coverage problem in wireless sensor networks,” IEEE Access, vol. 8, pp. 74 315–74 325, 2020.
  • [6] H. M. Ammari, β€œInvestigating the energy sink-hole problem in connected kπ‘˜kitalic_k-covered wireless sensor networks,” IEEE Trans. Comput., vol. 63, no. 11, pp. 2729–2742, Nov. 2014.
  • [7] β€”β€”, β€œA computational geometry-based approach for planar kπ‘˜kitalic_k-coverage in wireless sensor networks,” ACM Trans. Sens. Netw., vol. 19, no. 2, pp. 1–42, Feb. 2023.
  • [8] R. T. Rockafellar, Convex Analysis.   Princeton University Press, 1970.