Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Lizard Body Temperature Acquisition and Lizard Recognition Using Artificial Intelligence
Next Article in Special Issue
Two-Field Excitation for Contactless Inductive Flow Tomography
Previous Article in Journal
An Overall Automated Architecture Based on the Tapping Test Measurement Protocol: Hand Dexterity Assessment through an Innovative Objective Method
Previous Article in Special Issue
Simplified Beam Hardening Correction for Ultrafast X-ray CT Imaging of Binary Granular Mixtures
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Coupled Double-Layer Electrical Impedance Tomography-Based Sensing Skin for Pressure and Leak Detection

Department of Technical Physics, University of Eastern Finland, 70210 Kuopio, Finland
*
Author to whom correspondence should be addressed.
Sensors 2024, 24(13), 4134; https://doi.org/10.3390/s24134134
Submission received: 12 May 2024 / Revised: 20 June 2024 / Accepted: 24 June 2024 / Published: 26 June 2024
(This article belongs to the Special Issue Tomographic and Multi-Dimensional Sensors)

Abstract

:
There is an extensive need for surface sensors for applications such as tactile sensing for robotics, damage and strain detection for structural health monitoring and leak detection for buried structures. One type of surface sensor is electrical impedance tomography (EIT)-based sensing skins, which use electrically conductive coatings applied on the object’s surface to monitor physical or chemical phenomena on the surface. In this article, we propose a sensing skin with two electrically coupled layers separated by an insulator. Based on electrical measurements, the spatial distribution of the electrical coupling between the layers is estimated. This coupling is sensitive to both the pressure distribution on the surface and water entering between the layers through a leak. We present simulations and experimental studies to evaluate the feasibility of the proposed method for pressure sensing and leak detection. The results support the feasibility of the proposed method for both of these applications.

1. Introduction

The need for surface sensors is prevalent in many fields. The quickly developing field of robotics requires novel smart surface sensors for pressure and tactile sensing, e.g., for human–robot interactions and safety concerns [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16]. Similar sensing is required in some biomedical applications [2]. Sensing touch and pressure have been approached by several methods, such as grid-based capacitive sensors [14], functional materials [15], optical wave guides [16] and arrayless supercapacitive sensors [17]. In structural health monitoring (SHM), the structural integrity and lifetime expectation is analyzed based on continuous measurements. While these measurements are not necessarily performed on surface sensors, there are many promising candidates for that, such as fiber optic [18] or piezoresistive sensors [19] and strain sensors [20]. SHM is also used for leak detection in buried structures and monitoring geomembrane integrity in landfills or sewage ponds [21,22,23].
One of the latest developments in smart surface sensing is the so-called sensing skin based on electrical impedance tomography (EIT). It is an imaging method where the conductivity distribution of a target is estimated based on current and voltage measurements on the surface of the target. EIT-based sensing skins have already been studied in detecting the cracking and fatigue of concrete [24] and steel [25], composite material damage and strain [26,27,28,29,30,31,32,33,34], pH [35], surface pressure [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16] and temperature [36].
Although sensing skins show promise in several of the fields that they are studied in, they have certain challenges. One of them is that there can be several different physical and chemical factors contributing to the conductivity of the sensing skin, and, therefore, it may be difficult to distinguish the effect of interest. Hence, there is still need for innovative smart surface sensors for, e.g., detecting surface pressure in robotics and detecting leaks in buried structures.
In publication [37], a functional double-layer sensing skin was proposed for SHM: the two layers of the sensing skin were used separately, yet simultaneously, to distinguish different physical/chemical effects on the surface—particularly, one layer was sensitive to cracking while the other one was sensitive to both cracking and strong electrolytes, such as chlorides dissolved in water in contact with the surface. Another type of multilayer sensing skin proposed for tactile sensing is based on having a second more conductive layer that is weakly separated from the original measured surface [7,38,39]. Applied pressure introduces a contact between the two surfaces, thus effectively increasing the local conductivity of the measurement layer on the region of pressure.
In this paper, we propose a novel EIT-based sensing skin technique that is based on the use of two parallel electrically conductive layers. Similarly to the sensing skins presented in [38,39], these layers are coupled electrically, causing current to flow between the layers. In addition to the pressure on the surface, the electrical coupling is also sensitive to water entering in between the two layers. Unlike in the previous studies, we use two similar conductive layers and have electrodes attached on both of the conductive surfaces. We also measure currents flowing between electrodes of different surfaces. We formulate the mathematical model for the proposed technique and introduce its computational implementation. We study the feasibility of the proposed technique for pressure sensing and leak detection using numerical simulations and experiments.

2. Single-Layer EIT-Based Sensing Skin

EIT-based sensing skins consist of a conductive surface with electrodes attached to it, through which electrical currents can be injected into the surface. Traditionally, using a so-called current injection EIT device, the currents of the electrodes are set to different injection patterns and the resulting potentials are measured. In contrast, so-called potential excitation EIT devices set the potentials of the electrodes to different excitation patterns and measure the resulting currents. In many cases, both the theory and applications work similarly independent of which type of device is used. In this study, we use a potential excitation device for measurements, and hence we also formulate the theory accordingly. We use the conventional, single-layer sensing skin model in the reconstruction of the background conductivities of the two layers before they are installed in the coupled double-layer sensor system.

2.1. Forward Problem

The EIT measurements (whether potentials or currents are measured) are usually modeled by the complete electrode model (CEM), which consists of the following equations [40,41,42]:
· ( σ ( x ) u ( x ) ) = 0 , x Ω
u ( x ) + z l σ ( x ) u ( x ) n ^ = U l , x Ω e l , l = 1 , , N el
Ω e l σ ( x ) u ( x ) n ^ d S = I l , l = 1 , , N el
u ( x ) n ^ = 0 , x Ω l = 1 N el Ω e l ,
where σ and u are the conductivity and potential inside the imaging domain Ω , z l is the contact impedance of the l:th electrode, n ^ is the outward unit normal vector on the boundary Ω , U l and I l are the potential and current on the l:th electrode, Ω e l is the boundary region corresponding to the l:th electrode, x is the spatial variable and N el is the number of electrodes in the system. Either the potentials U l or currents I l can be solved from the CEM as long as the other ones are known. In this article, as we use a potential injection device, correspondingly, we solve the currents from the CEM. In general, the potential u inside the domain has to be solved simultaneously with the electrode potentials or currents.
Finding an analytical solution to the CEM is generally infeasible, and hence a numerical approximation is used instead. In this article, we use a finite element method (FEM) with piecewise linear basis functions to approximate the forward model (for details, see [43]). In FEM, the imaging domain is approximated by a mesh consisting of nodes and elements. Spatially distributed quantities inside the domain are then approximated as a linear combination of basis functions ϕ i . For example, the conductivity on the domain is then approximated as
σ ( x ) = i = 1 N σ i ϕ i ( x ) ,
where σ i R are coefficients of the basis functions and N is the number of basis functions. Using this approximation, we identify the conductivity by the vector of coefficients σ i in (5):
σ ˜ = σ 1 σ 2 σ N .
From the FEM, we obtain the computational model I ˜ σ ( σ ˜ ) mapping the conductivities to a vector containing electrode currents corresponding to all the used excitations.

2.2. Inverse Problem

To solve the coefficients of conductivity σ ˜ , a variety of methods exist. In this article, we use the Bayesian framework, in which unknown or uncertain parameter values are treated as random variables. The Bayesian framework also allows the use of prior knowledge of the measured system and the use of reference measurements. Further information can be found, e.g., in [44].
In this article, we use the maximum a posteriori (MAP) point as the reconstruction. The MAP point can be found by solving the following optimization problem [44]:
arg min σ ˜ F ( σ ˜ ) , F ( σ ˜ ) = | | L I ˜ σ σ ˜ I m | | 2 + R ( σ ˜ ) ,
where I m contains all the measured currents from different excitations stacked into a single vector, R ( σ ˜ ) is the prior potential function (more details in Section 3.2) and L is the Cholesky factor of the data precision matrix.
The optimization problem (7) can be solved by many optimization algorithms. In this study, we use the Gauss–Newton (GN) method [45,46] with a line search to find the step length. For each iteration σ ˜ i , the next GN iteration σ ˜ i + 1 with step length τ is computed by
σ ˜ i + 1 ( τ ) = σ ˜ i τ J T ( σ ˜ i ) Γ 1 J ( σ ˜ i ) + H R ( σ ˜ i ) 1 J T ( σ ˜ i ) Γ 1 I ˜ ( σ ˜ i ) I m + R ( σ ˜ i ) ,
where J is the Jacobian of I ˜ ( σ ˜ ) with regard to σ ˜ and H R and R are the Hessian and the gradient of R. The step length τ is chosen optimally with respect to the objective function
τ = arg min τ F σ ˜ i + 1 ( τ ) .

3. A Coupled Double-Layer Sensing Skin

The double-layer sensing skin presented in this paper consists of two conductive surfaces separated by an insulating layer, as depicted in Figure 1. Although the middle layer is insulating, some current will pass through it as displacement current since alternating potential difference is applied to the electrodes. Alternatively, in case of leak detection, water can form a conductive pathway between the layers. We use single-electrode excitations, i.e., one electrode at a time is excited to a higher potential, while the others remain grounded. Currents passing through the grounded electrodes are then measured.
In this paper, we model the double-layer sensing skin as two two-dimensional conductive surfaces modeled by the CEM and a coupling between the surfaces. More specifically, the coupling is the conductance per area, and thus, together with the local potential difference between the layers, it determines the current density flowing between the layers. This coupling strengthens if water leaks into the space between the layers (a conductive pathway) or if applied pressure pushes the layers closer to each other (increased local capacitance). For the purely capacitive case, the coupling can be computed by assuming that the layers form a parallel plate capacitor with possibly spatially varying distance d ( x ) between the layers. For a small area A having d ( x ) approximately constant, the capacitance between the layers for that area is
C = ϵ A d ( x ) ,
where ϵ is the electrical permittivity of the insulator between the layers. With signal angular frequency ω , this leads to conductance
G = j ω C = j ω ϵ A d ( x ) ,
where j is the imaginary unit. From this, we obtain
c ( x ) = G A = j ω ϵ d ( x ) .
In this article, we neglect the phase of the coupling, i.e., we use only the absolute value of the coupling computed by (12). Due to this approximation, the reconstructed couplings are to be interpreted only qualitatively.
The spatial coordinate system used for the two conductive surfaces consists of the two coordinates of the point along the surface and an additional index distinguishing the two layers. Therefore, the domain for potential u and conductivity σ is Ω R 2 × { 1 , 2 } . On the other hand, the coupling between the layers is shared by the two surfaces, and thus has the domain Ω c R 2 . We approximate the areas under the attached electrodes to remain in the same potential as the electrode, and therefore do not consider these areas as a part of the domain Ω , as the local potential u is not solved from the CEM there. However, the area under an electrode on the surface that the electrode is not attached on is part of domain Ω .
We notate the area where neither of the surfaces has electrodes as Ω 0 , i.e., Ω 0 = ( x , i ) Ω | ( x , j ) Ω j = 1 , 2 . In addition, the area under electrode l on the surface that the electrode is not attached to is notated as Ω e l . Similarly as before, the boundary corresponding to electrode l is denoted by Ω e l . The subdomains are illustrated in Figure 2. It is worth noting that Ω = Ω 0 l = 1 N el Ω e l and Ω c = x R 2 | j { 1 , 2 } such that ( x , j ) Ω .
Adding the coupling into the CEM requires modifying Equation (1), which is the conservation of currents in the continuum. In the coupled double-layer system, the current on each layer does not have to be conserved, but instead, any current sink on one layer has to be matched by a current source on the other layer, and vice versa. Thus, the total current is conserved, although, on the 2D surfaces, it is not.
Adding these sinks and sources to Equation (1) leads to the set of equations
· ( σ ( x , 1 ) u ( x , 1 ) ) = c ( x ) u ( x , 2 ) u ( x , 1 ) , ( x , 1 ) Ω 0 ,
· ( σ ( x , 2 ) u ( x , 2 ) ) = c ( x ) u ( x , 1 ) u ( x , 2 ) , ( x , 2 ) Ω 0 ,
where u ( x , i ) and σ ( x , i ) are the potential and conductivity, respectively, at point ( x , i ) and c ( x ) is the coupling at point x. The coupling c ( x ) , together with the potential difference in the two layers, determines the current density flowing between the surfaces. In Equations (13) and (14), we consider only the area without electrodes on either surface. Wherever one of the surfaces has an electrode, the potential function u is not defined, and instead we use the potential of the electrode
· ( σ ( x , i ) u ( x , i ) ) = c ( x ) U l u ( x , i ) , ( x , i ) Ω e l , l = 1 , , N el ,
where Ω e l is the area under electrode l on the surface that the electrode is not attached to.
Equations (13)–(15) can be combined into a single equation by first defining a function for the potential on the other layer
u 2 ( x , i ) = u ( x , 2 ) , ( x , i ) Ω 0 , i = 1 u ( x , 1 ) , ( x , i ) Ω 0 , i = 2 U l , ( x , i ) Ω e l .
Hence, we can write the CEM extended for the coupled double-layer sensing skin:
· ( σ ( x , i ) u ( x , i ) ) = c ( x ) u 2 ( x , i ) u ( x , i ) , ( x , i ) Ω ,
u ( x , i ) + z l σ ( x , i ) u ( x , i ) n ^ = U l ,
( x , i ) Ω e l , l = 1 , , N el
Ω e l σ ( x , i ) u ( x , i ) n ^ d S + Ω e l c ( x ) U l u ( x , i ) d Ω = I l , l = 1 , , N el
u ( x , i ) n ^ = 0 , ( x , i ) Ω l = 1 N el Ω e l .
Compared to the original CEM, Equation (19) is also modified to take into account the current sources and sinks on Ω e l that contribute to the measured currents as well.

3.1. Forward Model Solution

In this paper, we express the coupling c ( x ) via a separate mesh covering the domain Ω c . The coupling has values c i at nodes x i of this mesh. The function c ( x ) is then approximated in the FEM equations by
c ( x ) = i = 1 N c i ϕ i ( x , j ) , j = 1 , 2
where the coefficients c i are interpolated linearly from values c i based on the two meshes. Further, we notate that
c ˜ = c 1 c 2 c N ,
and
z ˜ = z 1 z 2 z N el .
The finite element formulation for the conventional CEM is well known, but in order to take into account the modifications made to model the double-layer case, we show the full derivation in Appendix A. We derive the FEM for the potential injection case, i.e., we solve the currents from the CEM, similarly as in [43]. We emphasize that the additions presented here work with few changes also using the more conventional case of solving the potentials from the CEM. We denote the computational model of the measurements obtained from FEM by I ˜ ( σ ˜ , c ˜ , z ˜ ) , which is a mapping from the parameter space to the vectors containing currents corresponding to all used potential excitations stacked into a single vector.
The basic principles of finding the coefficients of coupling c ˜ based on the measurements are the same as in Section 2.2. However, as the conductivity distributions of the surfaces and the contact impedances of the electrodes are not known, they are also estimated. Hence, we define a generalized estimate variable Θ , i.e., any one of c ˜ , σ ˜ , z ˜ or a vector obtained by stacking any combination of these. The conductivity estimate σ ˜ in the equations presented in Section 2.2 is then replaced by Θ . The reconstruction process and the constituents of Θ in each step are described in more detail in Section 4.3. For details on the simultaneous estimation of contact impedances, the reader is referred to [47,48].

3.2. Prior Functions

To mitigate the ill-posedness of the inverse problem, prior knowledge of the estimated distributions is used. In this study, we assume the conductivity distributions to have spatially smooth variations throughout the surfaces caused by inhomogeneities of the paint layers. On the other hand, the coupling is assumed to have a sparsely steep gradient caused by either compression of the insulating region under pressure or by water applied in between the two layers. Furthermore, both conductivity and the coupling are known to always be positive.
To promote smooth variations in the conductivity distribution, a Gaussian smoothness prior is used [49]. It has the functional form
R sm ( σ ) = σ σ e T Γ sm 1 σ σ e ,
where σ e is the expected value of σ and Γ sm determines the covariance between conductivities of all pairs of nodes in the FE mesh of each layer of the sensor. Similarly to [49], we use a Gaussian distance covariance, i.e.,
Γ sm ( i , j ) = a sm exp | | x i x j | | 2 2 b sm 2 + c sm δ i j ,
where x i , x j R are the positions of nodes i and j, a sm is the strength of the prior, b sm is a constant related to the length scale of the spatial smoothness, c sm is a small stabilizing constant and δ i j is the Kronecker delta. Throughout this study, we use c sm = 10 4 × a sm . To facilitate choosing the constant b sm , we define correlation length d cor as the length where the correlation of two nodes given by (25) is 0.01. Then, for the constant b sm , we have
b sm = d cor 2 ln 100 .
The total variation prior is used for favouring sparse spatial changes in the coupling. The total variation prior functional is [50]
R TV ( c ) = k = 1 n A k a TV D x c k 2 + D y c k 2 + β ,
where A k is the area of element k, n is the number of elements in the mesh, a TV is a strength parameter and β is a smoothing parameter used to make the functional differentiable. D x and D y are discrete differential operators corresponding to differentiation in x- and y-directions. The strength parameter a TV is chosen based on the procedure outlined in [50]. In short,
a TV = ln 1 p g max ,
where p is the probability of having the maximal gradient norm value of g max on the domain.
The positivity of estimates is enforced by a parabolic barrier function
R pos ( Θ ) = a pos { k | Θ k < Θ m , k } Θ k Θ m , k 2 ,
where a pos is the strength parameter for the constraint, Θ is the generalized estimate variable and Θ m is a vector containing minimum estimate values that are not penalized. Here, we allow for different minimum values for different elements of Θ because the scales of σ ˜ , c ˜ and z ˜ , any of which may be in Θ , may be different.

4. Simulations and Experimental Studies

4.1. Simulation Setup

We use numerical simulations to study the overall feasibility of the proposed method. The simulations were computed using the finite element (FE) approximation of the coupled CEM detailed in Section 3.1. The imaging domain was chosen to correspond to the experimental setup, i.e., two square surfaces with side lengths of 25 cm and 16 electrodes on each layer. Measurement setups were simulated corresponding to both pressure and leak detection. For all targets, three different measurements were simulated:
  • Reference measurements of the layers separately, i.e., with zero coupling between the layers.
  • Reference measurements of the assembled sensing skin, i.e., with homogeneous coupling.
  • Measurements with the effect (i.e., applied pressure or water) in place.
Gaussian noise with a standard deviation of 0.1% of each measurement value was added to all measurements.
In the simulations, the insulating layer had a thickness of 5 mm at rest and relative permittivity ϵ r = 4 . The couplings for simulations were computed using Equation (12), except for the water ingress.
For pressure sensing, two targets (targets I and II) were simulated, both modeling pressing the distance between the layers from 5 mm down to 2 mm on a circular region having a diameter of 5 cm. Target I was simulated with the sensing skin having homogeneous conductivity of 2 × 10 3 Ω 1 , while target II was simulated with the sensing skin having the conductivity distribution drawn from the smoothness prior distribution Equation (24) with parameter values d cor = 10 cm, a sm = 1.6 × 10 7 S 2 and σ e a vector containing values 2 × 10 3 S in each element. For leak detection, water ingress between the layers was simulated by an area of increased coupling (target III). Equation (12) was used for the background coupling, and the peak value on the wet area was chosen to have ten-fold coupling compared to the background. The initial conductivity distribution for this simulation was drawn from the same distribution as for target II but featured an additional increase in conductivity on the region of water ingress. For target I, contact impedances were set to a negligible value of 10 6 Ω . For targets II and III, the contact impedances were drawn from normal distribution N ( 1 , 1 ) Ω , with all values less than 10 6 Ω replaced by the value 10 6 Ω .
A mesh of 47,639 nodes and 92,490 elements was used for the FEM simulations. The element widths varied between 0.5 mm (near the electrodes) and 16 mm (on the central regions). The conductivity was represented in a mesh of 4122 nodes and 7616 elements (element widths ranging from 3 mm to 19 mm) in total for the two surfaces, and the coupling was represented in a mesh of 790 nodes and 1478 elements (element widths approximately 10 mm). The conductivity and coupling values were linearly interpolated to the nodal points of the FE mesh.

4.2. Experimental Setup

For measurements, a commercial potential excitation EIT device manufactured by Rocsole Ltd. (based in Kuopio, Finland) was used. The injection frequency was chosen as 320 kHz and single-electrode excitations were used.
The two conductive surfaces were made of rubber sheets painted with conductive paint. The paint was a mixture of rubber coating (Maston RUBBERcomp) and graphite powder (CRETACOLOR) added for conductivity in a mix ratio of 1:10 by weight. Electrodes were manufactured of copper tape, which was attached to the rubber sheet with the pre-applied adhesive before applying the paint. The pieces of copper tape were longer than the modeled electrode so that the rest of the piece of tape was turned on the back side of the rubber sheet. Short pieces of wire, with SMB connectors on the other end, were soldered on the copper tapes on the back side. The conductive paint was applied by brush so that the layer of paint covers the entire front surface of each rubber sheet, including the parts of the attached electrodes on the front side. Two layers of conductive paint were applied to obtain a more even coating. The first layer was left to dry overnight before applying the second layer of paint. The insulating layer between the conductive surfaces was made of polyester padding when sensing pressure, and polymer foam when sensing water. The sensing skin was assembled by placing the insulating layer between the rubber sheets, both of the rubber sheets having their front sides, i.e., the conductive surfaces, facing the insulating layer. The layers were not attached together but were only placed resting on top of each other in order to facilitate reusability with different insulating materials.
We note that the choice of the materials used in this experiment is not necessarily optimal for the final applications, but they were chosen to facilitate the laboratory experiments. For example, instead of using rubber sheets for the conductive layers for leak detection, a water-transparent material, such as conductive fabric (e.g., the one used in [9]), may be desirable. Otherwise, a rupture in one of the conductive layers is required for the water to enter the insulating layer and be detected.
A similar set of three measurements was taken as described in Section 4.1 for the simulations, except measurement set 1 (layers measured separately), which was taken only once at the beginning and used as the reference data for all the targets. The second set of measurements (assembled sensing skin without pressure or water) was taken separately for each target before applying the effect.
Weights of 11.6 kg and 3.4 kg were placed separately on the sensing skin for the measurements. Target IV was a 11.6 kg weight standing on the sensing skin on a circular surface having a diameter of 4.5 cm. Target V was a 20 cm tall circular cylinder (diameter 10 cm) weighing 3.4 kg resting on its side (see Figure 3). For target VI, 10 ml of water was injected in between the layers for the leak detection measurements.

4.3. Reconstruction

The reference measurements can be used in several ways to mitigate some of the modeling errors present in the system. The reconstruction process used in this study is detailed in this section. During the whole process, the noise of the forward model is assumed to consist of two additive Gaussian components. The first component, a baseline noise, has the same standard deviation s 1 for each measured value, with s 1 being 0.01% of the difference in maximum and minimum measured current. The second component has a standard deviation of 3% of each measured value. We note that the actual measured noise is well below this level, but these values try to approximately include the possible modeling errors as well. For further mitigating the effects of modeling errors, the Bayesian approximation error approach [51,52] could be applied, but it is left outside the scope of this paper.
To use the reference data, solving the inverse problem is divided into three steps:
  • From measurements made of the two layers separately, their initial conductivity distributions σ ˜ 0 and the contact impedances z ˜ 0 of electrodes are estimated as in single-layer EIT described in Section 2.
  • From measurements, where the layers were combined, but no pressure or water had yet been applied, a homogeneous estimate c 0 of the coupling is computed, using the conductivity distribution σ ˜ 0 and contact impedances z ˜ 0 estimated in step 1. Using the estimate c 0 , a correction to the measurements is computed as specified in reference [24].
  • From measurements after applying pressure or water, the coupling distribution is estimated together with re-estimation of the conductivity distribution.
The conductivity distribution is re-estimated at the final step in addition to the coupling distribution since applying pressure or water is also presumed to change the conductivity distribution. The estimations in steps 1 and 3 are performed using the GN algorithm outlined in Section 2.2. In step 2, the estimation is simple, as the homogeneous estimate is characterized by a single parameter. The estimation of contact impedances is performed as presented in [47,48].
In step 1, the generalized estimate variable Θ = σ ˜ T z ˜ T T . The prior potential function for this step is
R ( Θ ) = R sm ( σ ˜ ) + R pos ( Θ ) ,
with parameters d cor = 5 cm, a sm = 4 × 10 10 S 2 , σ e a vector containing values 2 × 10 3 S in each element, Θ m , k = 10 5 S and a pos = 10 5 S 2 for conductivity values and Θ m , k = 10 5 Ω and a pos = 10 3 Ω 2 for contact impedance values.
In step 3, the generalized estimate variable Θ = σ ˜ T z ˜ T T and the prior potential function is
R ( Θ ) = R sm ( σ ˜ ) + R TV ( c ˜ ) + R pos ( Θ ) .
Here, the parameters used are a sm = ( 0.01 × 1 / n i = 1 n σ ˜ 0 , i ) 2 , d cor = 5 cm, σ e = σ ˜ 0 , β = 10 4 . For reconstructions corresponding to measured data, the parameter a TV is chosen by Equation (28), with p = 0.975 and g max = c ˜ 0 / d elem , where d elem is the average edge length of mesh elements. For reconstructions corresponding to simulated data, the same procedure is used, but the parameter value is multiplied by an additional constant of 10 4 for pressure-sensing simulations and 10 2 for leak detection simulations. The positivity constraint parameters are Θ m , k = 10 4 S and a pos = 10 5 S 2 for conductivity values, and Θ m , k = 10 2 Sm 2 and a pos = 10 4 S 2 m 4 for coupling values.
The forward problem is solved by FEM as described in Section 3.1 using a mesh with 4122 nodes and 7616 elements (element widths ranging from 3 mm to 19 mm). The coupling is defined in a mesh containing 790 nodes and 1478 elements (element widths approximately 10 mm).

5. Results and Discussion

5.1. Simulations

The reconstructions computed from the pressure detection simulations (targets I and II) are shown in Figure 4. Both reconstructed locations and shapes of the inclusions are consistent with the targets within an accuracy expected from an EIT-based method. For target I, the value of coupling inside the inclusion (ca. 4.5 × 10 3 Sm 2 ) does not agree with the true value (ca. 5.6 × 10 3 Sm 2 ). The discrepancy of reconstructed and true values is likely caused by the total variation prior, as it penalizes stronger variations more. The homogeneous conductivity distribution of target I is reconstructed well, with the average conductivity of the reconstruction close to the true value of 2 × 10 3 S and variations in the order of 1% in the conductivity distribution of the reconstruction. For target II, the reconstructed coupling value inside the inclusion (ca. 5.2 × 10 3 Sm 2 ) agrees with the true value (ca. 5.6 × 10 3 Sm 2 ) better than for target I, but the shape of the inclusion is more distorted. The distortion of the shape may be caused by the proximity of the boundary and electrodes, increasing the modeling error caused by using a different FE mesh for the reconstruction than that used for the data simulation. The reconstructed conductivity distribution of target II shows all the features of the target within the expected accuracy.
These reconstructions show that, in the absence of large modeling errors, the size and location of single areas of pressure can be reconstructed fairly accurately even with inhomogeneous background conductivity of the sensing skin. The robustness of quantifying the value of coupling inside an inclusion has to be improved for applications requiring the absolute value, e.g., for quantifying the applied pressure on a surface. For many applications, however, the location of the applied pressure is enough, and, for those applications, the proposed method is viable.
In Figure 5, the reconstructions corresponding to water ingress simulation (target III) are shown. The area of water ingress is depicted well in the coupling reconstruction, with the numerical value (ca. 2.1 × 10 2 Sm 2 ) being quite close to correct (ca. 2.4 × 10 2 Sm 2 ). In addition, the conductivities are reconstructed well within the accuracy expected, with also the increased conductivity in the area of water ingress visible.
This evidence supports the feasibility of leak detection by the proposed double-layer sensing skin. The minor discrepancies between the reconstructed coupling and the true distribution are irrelevant for the application of leak detection.

5.2. Experimental Results

The experimental results of pressure detection (targets IV and V) are shown in Figure 6. In this figure, only the reconstructions of the coupling distribution are represented. Both reconstructions show an increased coupling on or close to the correct area. Both reconstructions feature some artifacts of decreased coupling as well. As the applied pressure should not cause a decrease in coupling, these artifacts are easily distinguished from the locations of actual applied pressure.
The reconstruction of target IV in Figure 6a shows the elongated shape of the pressure distribution well. However, the area of increased coupling is reconstructed as continuing through the sensing skin, from one boundary to another, whereas the edges of the object resting on the surface are 2.5 cm from both boundaries. This error is likely caused by the surface deforming further than the exact location of applied pressure.
From the reconstruction of target V in Figure 6b, the size and shape of the inclusion cannot be inferred. The spreading of the increased coupling to a larger area is likely caused by the surface deforming on a larger area than the applied pressure. This deformation may also include wrinkling of the surface, which may explain some artifacts even having decreased electrical coupling. Another possible cause of significant error is the alteration in contact impedances caused by the stretching and wrinkling of the conductive sheets under the electrodes.
Although the proposed method seems prone to having some artifacts and inaccuracies in the location and shape of the observed distribution, the results support the feasibility of detecting the location of applied pressure using the proposed method. In this study, we used weights in the order of several kilograms, but we emphasize that the weight range that the sensing skin is sensitive to can be tuned by choosing the materials of conductive and insulating layers.
In Figure 7, the reconstructed coupling distribution corresponding to the water ingress (target VI) is shown. The reconstruction features increased coupling coinciding well with the approximate location of the water ingress. However, in addition, there is another region of increased coupling with no clear explanation. This may be caused, for example, by variations in contact impedances or by slight bending of the edges of the conductive layers. Despite the artifact, the location of water ingress is clearly visible, and thus the results support the viability of using the proposed method for leak monitoring.

6. Conclusions

In this article, we studied the feasibility of using EIT-based electrically coupled double-layer sensing skin for pressure sensing and leak detection. The aim of the method is to reconstruct the spatial distribution of the electrical coupling between the two layers. This distribution contains information of possible external pressure pushing the layers closer to each other, or a conductive path of water formed between the layers.
The feasibility of the proposed method for pressure sensing and leak detection was tested on both simulated and measured data. The reconstructions corresponding to simulated data localized both the applied pressure and the ingress of water, with slight distortions in the shape or the reconstructed electrical coupling value. The reconstructions from measured data had some artifacts, but the approximate locations of the applied pressure or water leak could still be inferred from the reconstructions. While the reconstruction from a circular pressure distribution failed to show the shape of the applied pressure, the reconstruction of a line-like target showed the elongated shape of the applied pressure.
The results of both the simulations and the experimental studies support the feasibility of the proposed method for pressure sensing and leak detection. Hence, the proposed method may be used in applications such as tactile sensing in robotics and leak monitoring in buried structures.

Author Contributions

Conceptualization, P.K. and A.S.; methodology, P.K.; software, P.K.; validation, P.K. and A.S.; formal analysis, P.K.; investigation, P.K.; resources, A.S.; data curation, P.K.; writing—original draft preparation, P.K.; writing—review and editing, P.K. and A.S.; visualization, P.K.; supervision, A.S.; project administration, P.K. and A.S.; funding acquisition, P.K. and A.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Academy of Finland, Centre of Excellence of Inverse Modelling and Imaging, 2018-2025 (project 353085) and by Auramo-foundation.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original data presented in the study are openly available in Github repository PetriKuusela/DoubleLayerSensingSkin at https://github.com/PetriKuusela/DoubleLayerSensingSkin, accessed on 13 May 2024.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
EITElectrical impedance tomography
SHMStructural health monitoring
CEMComplete electrode model
FEMFinite element method
MAPMaximum a posteriori
GNGauss–Newton
FEFinite element

Appendix A. FEM Formulation for Coupled Double-Layer Sensing Skin

In this appendix, we present the derivation of the finite element method formulation for the coupled double-layer sensing skin. The derivation follows along the lines of previous similar derivations, such as [41,43]. For the derivation, in addition to requiring (17)–(20), i.e.,
· ( σ ( x , i ) u ( x , i ) ) = c ( x ) u 2 ( x , i ) u ( x , i ) , ( x , i ) Ω ,
u ( x , i ) + z l σ ( x , i ) u ( x , i ) n ^ = U l , ( x , i ) Ω e l , l = 1 , , N el
Ω e l σ ( x , i ) u ( x , i ) n ^ d S + Ω e l c ( x ) U l u ( x , i ) d Ω = I l , l = 1 , , N el
u ( x , i ) n ^ = 0 , ( x , i ) Ω l = 1 N el Ω e l ,
we require conservation of currents on the electrodes
l = 1 N el I l = 0 .
To derive the variational form, we multiply both sides of (A1) by a test function v : Ω R and integrate over the domain Ω to obtain
Ω v ( x , i ) · ( σ ( x , i ) u ( x , i ) ) d Ω = Ω v ( x , i ) c ( x ) u 2 ( x , j ) u ( x , i ) d Ω .
Next, we use Green’s formula to obtain
Ω v ( x , i ) σ ( x , i ) u ( x , i ) n ^ d S Ω σ ( x , i ) u ( x , i ) · v ( x , i ) d Ω = Ω v ( x , i ) c ( x ) ( u 2 ( x , i ) u ( x , i ) ) d Ω .
Next, the boundary integral is separated into electrode and non-electrode boundaries and boundary conditions (A2) and (A4) are used. As a result, we obtain
l = 1 N el 1 z l Ω e l v ( x , i ) ( u ( x , i ) U l ) d S Ω σ ( x , i ) u ( x , i ) · v ( x , i ) d Ω = Ω v ( x , i ) c ( x ) ( u 2 ( x , i ) u ( x , i ) ) d Ω .
Now, we split the integral containing u 2 ( x , i ) into electrode and non-electrode parts using (16) to obtain
l = 1 N el 1 z l Ω e l v ( x , i ) ( u ( x , i ) U l ) d S Ω σ ( x , i ) u ( x , i ) · v ( x , i ) d Ω = Ω 0 v ( x , i ) c ( x ) ( u 2 ( x , i ) u ( x , i ) ) d Ω + l = 1 N el Ω e l v ( x , i ) c ( x ) ( U l u ( x , i ) ) d Ω .
Next, we introduce a test vector ι R N el . By multiplying (A2) by ι l and integrating over electrode l, we obtain
Ω e l u ( x , i ) + z l σ ( x , i ) u ( x , i ) n ^ ι l d S = Ω e l U l ι l d S .
To reformulate the second term on the left-hand side, we use (A3), i.e.,
Ω e l σ ( x , i ) u ( x , i ) n ^ d S = I l Ω e l c ( x ) U l u ( x , i ) d Ω
to obtain
Ω e l u ( x , i ) ι l d S Ω e l U l ι l d S z l I l ι l z l ι l Ω e l c ( x ) U l u ( x , i ) d Ω = 0 .
By dividing this by z l and then summing over all electrodes, we obtain
l = 1 N el 1 z l Ω e l u ( x , i ) ι l d S l = 1 N el 1 z l Ω e l U l ι l d S l = 1 N el I l ι l l = 1 N el ι l Ω e l c ( x ) U l u ( x , i ) d Ω = 0 .
Now, by summing (A9) and (A13), we obtain the variational form
l = 1 N el 1 z l Ω e l u ( x , i ) ( v ( x , i ) ι l ) d S Ω σ ( x , i ) u ( x , i ) · v ( x , i ) d Ω + l = 1 N el I l ι l + l = 1 N el Ω e l c ( x ) u ( x , i ) ( v ( x , i ) ι l ) d Ω + Ω 0 c ( x ) ( u ( x , i ) u 2 ( x , i ) ) v ( x , i ) d Ω = l = 1 N el 1 z l Ω e l U l ( v ( x , i ) ι l ) d S + l = 1 N el Ω e l c ( x ) U l ( v ( x , i ) ι l ) d Ω
Next, we use the finite element approximation
u ˜ ( x , i ) = j = 1 N α j ϕ j ( x , i ) , ( x , i ) Ω ,
where u ˜ is the approximated potential, N is the number of basis functions used and α j are coefficients for basis functions ϕ j . In this article, we assume that none of the basis functions ϕ j have support on both layers. In addition, we want to have basis functions satisfying (A5) for the electrode currents, so we define
I ˜ = i = 1 N el 1 β i n i ,
where I ˜ is a vector containing the electrode currents and β i are coefficients for the basis functions n i R N el , which are of the form
n 1 = [ 1 , 1 , 0 , , 0 ] T , n 2 = [ 1 , 0 , 1 , , 0 ] T , n N el 1 = [ 1 , 0 , 0 , , 1 ] T .
Now, by using the Galerkin scheme, i.e., using the basis functions ϕ i and n i as the test functions v and ι l , we obtain the FEM system of equations in matrix form:
A θ = b ,
where
θ = α β ,
where α and β are vectors containing the coefficients α i and β i . Matrix A is of the form
A = B 0 C D ,
where
B ( j , k ) = Ω σ ( x , i ) ϕ j ( x , i ) · ϕ k ( x , i ) + c ( x ) ϕ j ( x , i ) ϕ k ( x , i ) d Ω + l = 1 N el 1 z l Ω e l ϕ j ( x , i ) ϕ k ( x , i ) d S , j , k = 1 , , N nodes j and k on the same layer
B ( j , k ) = Ω c ( x ) ϕ j ( x , i ) ϕ k ( x , i ) d Ω , j , k = 1 , , N nodes j and k on different layers .
C ( j , k ) = l = 1 N el 1 z l Ω e l ϕ k ( x , i ) ( n j ) l d S l = 1 N el ( n j ) l Ω e l c ( x ) ϕ k ( x , i ) d Ω = 1 z j + 1 Ω e j + 1 ϕ k ( x , i ) d S 1 z 1 Ω e 1 ϕ k ( x , i ) d S + Ω e j + 1 c ( x ) ϕ k ( x , i ) d Ω Ω e 1 c ( x ) ϕ k ( x , i ) d Ω k = 1 , , N , j = 1 , , N el 1
D ( j , k ) = l = 1 N el ( n j ) l ( n k ) l = 1 , j k 2 , j = k j , k = 1 , , N el 1 .
Vector b in (A18) has the elements
b ( j ) = l = 1 N el U l z l Ω e l ϕ j ( x , i ) d S + Ω e l c ( x ) U l ϕ j ( x , i ) d Ω , 1 j N l = 1 N el U l z l ( n j ) l | Ω e l | + Ω e l c ( x ) U l ( n j ) l d Ω , N + 1 j N + N el 1
= l = 1 N el U l z l Ω e l ϕ j ( x , i ) d S + Ω e l c ( x ) U l ϕ j ( x , i ) d Ω , 1 j N U j + 1 z j + 1 | Ω e j + 1 | U 1 z 1 | Ω e 1 | + Ω e j + 1 c ( x ) U j + 1 d Ω Ω e 1 c ( x ) U 1 d Ω , N + 1 j N + N el 1 ,
where | Ω e l | denotes the area of electrode l.

Computing the Jacobian

The derivatives of the currents I ˜ with respect to the coupling coefficients c j (defined in (21)) are obtained similarly as the derivatives with respect to conductivity are obtained conventionally. We start by differentiating (A18) with respect to c j to obtain
c j A θ = b c j .
Hence,
A c j θ + A θ c j = b c j
θ c j = A 1 b c j A c j θ .
Now, we may obtain the derivatives of the currents by differentiating (A16)
I ˜ c j = k = 1 N el 1 β k c j n k .
Here, the derivative on the right-hand side is computed by (A28), since
β k = θ N + k .
Next, we compute the derivative of the FEM matrix
A c j = B c j 0 C c j 0 ,
where
B c j ( k , l ) = Ω ϕ j ( x , i ) ϕ k ( x , i ) ϕ l ( x , i ) d Ω , k , l = 1 , , N k , l on the same layer B c j ( k , l ) = Ω ϕ j ( x , i ) ϕ k ( x , i ) ϕ l ( x , i ) d Ω , k , l = 1 , , N k , l on different layers . C c j ( k , l ) = Ω e k + 1 ϕ j ( x , i ) ϕ l ( x , i ) d Ω Ω e 1 ϕ j ( x , i ) ϕ l ( x , i ) d Ω , k = 1 , , N , j = 1 , , N el 1 .
Finally, we compute the elements of the derivative of the vector b
b c j ( k ) = l = 1 N el Ω e l ϕ j ( x , i ) U l ϕ k ( x , i ) d Ω , 1 k N Ω e k + 1 ϕ j ( x , i ) U k + 1 d Ω Ω e 1 ϕ j ( x , i ) U 1 d Ω , N + 1 j N + N el 1 ,
and thus we have all parts required for calculating the Jacobian with respect to coupling c ˜ .

References

  1. Soleimani, M.; Friedrich, M. E-Skin Using Fringing Field Electrical Impedance Tomography with an Ionic Liquid Domain. Sensors 2022, 22, 5040. [Google Scholar] [CrossRef] [PubMed]
  2. Tallman, T.; Gungor, S.; Wang, K.; Bakis, C. Tactile imaging and distributed strain sensing in highly flexible carbon nanofiber/polyurethane nanocomposites. Carbon 2015, 95, 485–493. [Google Scholar] [CrossRef]
  3. Silvera-Tawil, D.; Rye, D.; Soleimani, M.; Velonaki, M. Electrical Impedance Tomography for Artificial Sensitive Robotic Skin: A Review. IEEE Sens. J. 2015, 15, 2001–2016. [Google Scholar] [CrossRef]
  4. Yang, Y.; Zhou, W.; Chen, X.; Ye, J.; Wu, H. A flexible touching sensor with the variation of electrical impedance distribution. Measurement 2021, 183, 109778. [Google Scholar] [CrossRef]
  5. Zhang, H.; Kalra, A.; Lowe, A.; Yu, Y.; Anand, G. A Hydrogel-Based Electronic Skin for Touch Detection Using Electrical Impedance Tomography. Sensors 2023, 23, 1571. [Google Scholar] [CrossRef] [PubMed]
  6. Yao, A.; Soleimani, M. A pressure mapping imaging device based on electrical impedance tomography of conductive fabrics. Sens. Rev. 2012, 32, 310–317. [Google Scholar] [CrossRef]
  7. Liu, K.; Wu, Y.; Wang, S.; Wang, H.; Chen, H.; Chen, B.; Yao, J. Artificial Sensitive Skin for Robotics Based on Electrical Impedance Tomography. Adv. Intell. Syst. 2020, 2, 1900161. [Google Scholar] [CrossRef]
  8. Liu, J.; Liu, N.; Wang, P.; Wang, M.; Guo, S. Array-less touch position identification based on a flexible capacitive tactile sensor for human-robot interactions. In Proceedings of the 2019 IEEE 4th International Conference on Advanced Robotics and Mechatronics (ICARM), Toyonaka, Japan, 3–5 July 2019; pp. 458–462. [Google Scholar] [CrossRef]
  9. Duan, X.; Taurand, S.; Soleimani, M. Artificial skin through super-sensing method and electrical impedance data from conductive fabric with aid of deep learning. Sci. Rep. 2019, 9, 8831. [Google Scholar] [CrossRef]
  10. Lee, H.; Park, H.; Serhat, G.; Sun, H.; Kuchenbecker, K.J. Calibrating a Soft ERT-Based Tactile Sensor with a Multiphysics Model and Sim-to-real Transfer Learning. In Proceedings of the 2020 IEEE International Conference on Robotics and Automation (ICRA), Paris, France, 31 May–31 August 2020; pp. 1632–1638. [Google Scholar] [CrossRef]
  11. Chen, Y.; Yu, M.; Bruck, H.A.; Smela, E. Compliant multi-layer tactile sensing for enhanced identification of human touch. Smart Mater. Struct. 2018, 27, 125009. [Google Scholar] [CrossRef]
  12. Yao, A.; Yang, L.; Seo, K.; Soleimani, M. EIT-Based Fabric Pressure Sensing. In Computational and Mathematical Methods in Medicine; Wiley: Hoboken, NJ, USA, 2013; p. 405325. [Google Scholar] [CrossRef]
  13. Jamshidi, M.; Park, C.B.; Azhari, F. An EIT-based piezoresistive sensing skin with a lattice structure. Mater. Des. 2023, 233, 112227. [Google Scholar] [CrossRef]
  14. Weichart, J.; Roman, C.; Hierold, C. Tactile Sensing with Scalable Capacitive Sensor Arrays on Flexible Substrates. J. Microelectromechanical Syst. 2021, 30, 915–929. [Google Scholar] [CrossRef]
  15. Yang, T.; Xie, D.; Li, Z.; Zhu, H. Recent advances in wearable tactile sensors: Materials, sensing mechanisms, and device performance. Mater. Sci. Eng. R Rep. 2017, 115, 1–37. [Google Scholar] [CrossRef]
  16. Kim, J.T.; Choi, H.; Shin, E.; Park, S.; Kim, I.G. Graphene-based optical waveguide tactile sensor for dynamic response. Sci. Rep. 2018, 8, 16118. [Google Scholar] [CrossRef] [PubMed]
  17. Liu, J.; Liu, N.; Hashimoto, K.; Meng, C.; Guo, S. Touch position identification based on a flexible array-less supercapacitive tactile sensor. AIP Adv. 2019, 9, 015026. [Google Scholar] [CrossRef]
  18. Guo, H.; Xiao, G.; Mrad, N.; Yao, J. Fiber Optic Sensors for Structural Health Monitoring of Air Platforms. Sensors 2011, 11, 3687–3705. [Google Scholar] [CrossRef] [PubMed]
  19. Nauman, S. Piezoresistive Sensing Approaches for Structural Health Monitoring of Polymer Composites—A Review. Eng 2021, 2, 197–226. [Google Scholar] [CrossRef]
  20. Yao, Y.; Glisic, B. Sensing sheets: Optimal arrangement of dense array of sensors for an improved probability of damage detection. Struct. Health Monit. 2015, 14, 513–531. [Google Scholar] [CrossRef]
  21. Wan, Y.; Dong, Z.; Cai, Y.; Xue, Q.; Liu, K.; Liu, L.; Guo, D. Geomembrane leaks detection and leakage correlation factor analysis of composite liner systems for fifty-five (55) solid waste landfills in China. Environ. Technol. Innov. 2023, 32, 103308. [Google Scholar] [CrossRef]
  22. Lugli, F.; Mahler, C.F. Analytical study of the performance of a geomembrane leak detection system. Waste Manag. Res. 2016, 34, 482–486. [Google Scholar] [CrossRef] [PubMed]
  23. Frangos, W. Electrical detection of leaks in lined waste disposal ponds. Geophysics 1997, 62, 1737–1744. [Google Scholar] [CrossRef]
  24. Hallaji, M.; Pour-Ghaz, M. A new sensing skin for qualitative damage detection in concrete elements: Rapid difference imaging with electrical resistance tomography. NDT E Int. 2014, 68, 13–21. [Google Scholar] [CrossRef]
  25. Yao, Y.; Glisic, B. Detection of Steel Fatigue Cracks with Strain Sensing Sheets Based on Large Area Electronics. Sensors 2015, 15, 8088–8108. [Google Scholar] [CrossRef] [PubMed]
  26. Thomas, A.; Kim, J.; Tallman, T.; Bakis, C. Damage detection in self-sensing composite tubes via electrical impedance tomography. Compos. Part B 2019, 177, 107276. [Google Scholar] [CrossRef]
  27. Loh, K.J.; Hou, T.C.; Lynch, J.P.; Kotov, N.A. Carbon Nanotube Sensing Skins for Spatial Strain and Impact Damage Identification. J. Nondestruct. Eval. 2009, 28, 9–25. [Google Scholar] [CrossRef]
  28. Gupta, R.; Mitchell, D.; Blanche, J.; Harper, S.; Tang, W.; Pancholi, K.; Baines, L.; Bucknall, D.G.; Flynn, D. A Review of Sensing Technologies for Non-Destructive Evaluation of Structural Composite Materials. J. Compos. Sci. 2021, 5, 319. [Google Scholar] [CrossRef]
  29. Chung, D. A review to elucidate the multi-faceted science of the electrical-resistance-based strain/temperature/damage self-sensing in continuous carbon fiber polymer-matrix structural composites. J. Mater. Sci. 2023, 58, 483–526. [Google Scholar] [CrossRef]
  30. Tallman, T.N.; Wang, K.W. Damage and strain identification in multifunctional materials via electrical impedance tomography with constrained sine wave solutions. Struct. Health Monit. 2016, 15, 235–244. [Google Scholar] [CrossRef]
  31. Tallman, T.N.; Gungor, S.; Wang, K.W.; Bakis, C.E. Damage detection via electrical impedance tomography in glass fiber/epoxy laminates with carbon black filler. Struct. Health Monit. 2015, 14, 100–109. [Google Scholar] [CrossRef]
  32. Tallman, T.N.; Gungor, S.; Wang, K.W.; Bakis, C.E. Damage detection and conductivity evolution in carbon nanofiber epoxy via electrical impedance tomography. Smart Mater. Struct. 2014, 23, 045034. [Google Scholar] [CrossRef]
  33. Sannamani, M.; Gao, J.; Chen, W.W.; Tallman, T.N. Damage detection in non-planar carbon fiber-reinforced polymer laminates via electrical impedance tomography with surface-mounted electrodes and directional sensitivity matrices. Compos. Sci. Technol. 2022, 224, 109429. [Google Scholar] [CrossRef]
  34. Augustin, T.; Karsten, J.; Fiedler, B. Detection and localization of impact damages in carbon nanotube–modified epoxy adhesive films with printed circuits. Struct. Health Monit. 2018, 17, 1166–1177. [Google Scholar] [CrossRef]
  35. Hou, T.C.; Loh, K.; Lynch, J. Spatial conductivity mapping of carbon nanotube composite thin films by electrical impedance tomography for sensing applications. Nanotechnology 2007, 18, 315501. [Google Scholar] [CrossRef]
  36. Rashetnia, R.; Hallaji, M.; Smyl, D.; Seppänen, A.; Pour-Ghaz, M. Detection and localization of changes in two-dimensional temperature distributions by electrical resistance tomography. Smart Mater. Struct. 2017, 26, 115021. [Google Scholar] [CrossRef]
  37. Seppänen, A.; Hallaji, M.; Pour-Ghaz, M. A functionally layered sensing skin for the detection of corrosive elements and cracking. Struct. Health Monit. 2017, 16, 215–224. [Google Scholar] [CrossRef]
  38. Alirezaei, H.; Nagakubo, A.; Kuniyoshi, Y. A tactile distribution sensor which enables stable measurement under high and dynamic stretch. In Proceedings of the IEEE Symposium on 3D User Interfaces, 3DUI 2009, Lafayette, LA, USA, 14–15 March 2009; IEEE Computer Society: Los Alamitos, CA, USA, 2009; pp. 87–93. [Google Scholar] [CrossRef]
  39. Tawil, D.S.; Rye, D.C.; Velonaki, M. Touch modality interpretation for an EIT-based sensitive skin. In Proceedings of the IEEE International Conference on Robotics and Automation, ICRA 2011, Shanghai, China, 9–13 May 2011; pp. 3770–3776. [Google Scholar] [CrossRef]
  40. Cheng, K.S.; Isaacson, D.; Newell, J.C.; Gisser, D.G. Electrode Models for Electric Current Computed Tomography. IEEE Trans. Biomed. Eng. 1989, 36, 918–924. [Google Scholar] [CrossRef]
  41. Somersalo, E.; Cheney, M.; Isaacson, D. Existence and Uniqueness for Electrode Models for Electric Current Computed Tomography. SIAM J. Appl. Math. 1992, 52, 1023–1040. [Google Scholar] [CrossRef]
  42. Vauhkonen, P.J.; Vauhkonen, M.; Savolainen, T.; Kaipio, J.P. Three-Dimensional Electrical Impedance Tomography Based on the Complete Electrode Model. IEEE Trans. Biomed. Imaging 1999, 46, 1150–1160. [Google Scholar] [CrossRef]
  43. Voss, A. Imaging Moisture Flows in Cement-Based Materials Using Electrical Capacitance Tomography. Ph.D. Thesis, University of Eastern Finland, Kuopio, Finland, 2020. [Google Scholar]
  44. Kaipio, J.; Somersalo, E. Statistical and Computational Inverse Problems; Springer Science+Business Media Inc.: New York, NY, USA, 2004. [Google Scholar]
  45. Nocedal, J.; Wright, S.J. Numerical Optimization; Springer: New York, NY, USA, 2006. [Google Scholar]
  46. Vauhkonen, P. Image Reconstruction in Three-Dimensional Electrical Impedance Tomography. Ph.D. Thesis, University of Kuopio, Kuopio, Finland, 2004. [Google Scholar]
  47. Vilhunen, T.; Kaipio, J.P.; Vauhkonen, P.J.; Savolainen, T.; Vauhkonen, M. Simultaneous reconstruction of electrode contact impedances and internal electrical properties: I. Theory. Meas. Sci. Technol. 2002, 13, 1848. [Google Scholar] [CrossRef]
  48. Heikkinen, L.M.; Vilhunen, T.; West, R.M.; Vauhkonen, M. Simultaneous reconstruction of electrode contact impedances and internal electrical properties: II. Laboratory experiments. Meas. Sci. Technol. 2002, 13, 1855. [Google Scholar] [CrossRef]
  49. Lipponen, A.; Seppänen, A.; Kaipio, J. Electrical impedance tomography imaging with reduced-order model based on proper orthogonal decomposition. J. Electron. Imaging 2013, 22, 023008. [Google Scholar] [CrossRef]
  50. González, G.; Huttunen, J.M.J.; Kolehmainen, V.; Seppänen, A.; Vauhkonen, M. Experimental Evaluation of 3D Electrical Impedance Tomography with Total Variation Prior. Inverse Probl. Sci. Eng. 2016, 24, 1411–1431. [Google Scholar] [CrossRef]
  51. Nissinen, A.; Heikkinen, L.M.; Kaipio, J.P. The Bayesian approximation error approach for electrical impedance tomography—Experimental results. Meas. Sci. Technol. 2007, 19, 015501. [Google Scholar] [CrossRef]
  52. Nissinen, A.; Kolehmainen, V.P.; Kaipio, J.P. Compensation of Modelling Errors Due to Unknown Domain Boundary in Electrical Impedance Tomography. IEEE Trans. Med. Imaging 2011, 30, 231–242. [Google Scholar] [CrossRef] [PubMed]
Figure 1. (a) A schematic diagram showing the layers of the sensing skin and locations of electrodes. The electrodes are drawn in copper brown. (b) The two conductive layers of the sensing skin without the insulating middle layer. The black side of the layer is the conductive paint and the yellow side shows the rubber base of the layer. (c) The assembled sensing skin with the insulating layer in the middle.
Figure 1. (a) A schematic diagram showing the layers of the sensing skin and locations of electrodes. The electrodes are drawn in copper brown. (b) The two conductive layers of the sensing skin without the insulating middle layer. The black side of the layer is the conductive paint and the yellow side shows the rubber base of the layer. (c) The assembled sensing skin with the insulating layer in the middle.
Sensors 24 04134 g001
Figure 2. A schematic diagram of the imaging domain of the double-layer sensing skin. Ω 0 is the area common for both layers. Ω e l is the area under electrode e l on the surface that the electrode is not attached to and Ω e l is the domain of the boundary corresponding to electrode e l . Electrodes are drawn in copper color on the layer that they are attached to. Labels for electrodes 1 and 2 are shown as examples, with blue highlights used to differentiate where the boundary and where the area of the electrode are meant.
Figure 2. A schematic diagram of the imaging domain of the double-layer sensing skin. Ω 0 is the area common for both layers. Ω e l is the area under electrode e l on the surface that the electrode is not attached to and Ω e l is the domain of the boundary corresponding to electrode e l . Electrodes are drawn in copper color on the layer that they are attached to. Labels for electrodes 1 and 2 are shown as examples, with blue highlights used to differentiate where the boundary and where the area of the electrode are meant.
Sensors 24 04134 g002
Figure 3. (a) Target IV: 11.6 kg weight and (b) target V: 3.4 kg weight resting on the sensing skin during measurements.
Figure 3. (a) Target IV: 11.6 kg weight and (b) target V: 3.4 kg weight resting on the sensing skin during measurements.
Sensors 24 04134 g003
Figure 4. (a) The coupling distribution and (b) the conductivity distributions of target I, and (c) The coupling distribution and (d) the conductivity distributions of target II. First row: the true values. Second row: reconstructions. The conductivity distributions of the two layers are shown side by side. Couplings are in units of Sm 2 and conductivities in units of S.
Figure 4. (a) The coupling distribution and (b) the conductivity distributions of target I, and (c) The coupling distribution and (d) the conductivity distributions of target II. First row: the true values. Second row: reconstructions. The conductivity distributions of the two layers are shown side by side. Couplings are in units of Sm 2 and conductivities in units of S.
Sensors 24 04134 g004
Figure 5. (a) The coupling distribution and (b) the conductivity distributions of target III. First row: the true values. Second row: reconstructions. The conductivity distributions of the two layers are shown side by side. Couplings are in units of Sm 2 and conductivities in units of S.
Figure 5. (a) The coupling distribution and (b) the conductivity distributions of target III. First row: the true values. Second row: reconstructions. The conductivity distributions of the two layers are shown side by side. Couplings are in units of Sm 2 and conductivities in units of S.
Sensors 24 04134 g005
Figure 6. Reconstructed coupling distributions computed from measured data. (a) The black line corresponds to the central line of contact of the 3.4 kg cylinder resting on its side. (b) The black circle corresponds to the bottom surface of the 11.6 kg weight. The units are in Sm 2 .
Figure 6. Reconstructed coupling distributions computed from measured data. (a) The black line corresponds to the central line of contact of the 3.4 kg cylinder resting on its side. (b) The black circle corresponds to the bottom surface of the 11.6 kg weight. The units are in Sm 2 .
Sensors 24 04134 g006
Figure 7. Reconstructed coupling distribution computed from measured data corresponding to water injected between the layers. The black circle marks the approximate location of water injection. The units are in Sm 2 .
Figure 7. Reconstructed coupling distribution computed from measured data corresponding to water injected between the layers. The black circle marks the approximate location of water injection. The units are in Sm 2 .
Sensors 24 04134 g007
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kuusela, P.; Seppänen, A. A Coupled Double-Layer Electrical Impedance Tomography-Based Sensing Skin for Pressure and Leak Detection. Sensors 2024, 24, 4134. https://doi.org/10.3390/s24134134

AMA Style

Kuusela P, Seppänen A. A Coupled Double-Layer Electrical Impedance Tomography-Based Sensing Skin for Pressure and Leak Detection. Sensors. 2024; 24(13):4134. https://doi.org/10.3390/s24134134

Chicago/Turabian Style

Kuusela, Petri, and Aku Seppänen. 2024. "A Coupled Double-Layer Electrical Impedance Tomography-Based Sensing Skin for Pressure and Leak Detection" Sensors 24, no. 13: 4134. https://doi.org/10.3390/s24134134

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop