Next Article in Journal
Virtual Antenna Array and Fractional Fourier Transform-Based TOA Estimation for Wireless Positioning
Previous Article in Journal
Embedded Processing and Compression of 3D Sensor Data for Large Scale Industrial Environments
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimating the Volume of Unknown Inclusions in an Electrically Conducting Body with Voltage Measurements

Polytechnical Department of Engineering and Architecture, University of Udine, Via delle Scienze 206, 33100 Udine, Italy
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(3), 637; https://doi.org/10.3390/s19030637
Submission received: 14 December 2018 / Revised: 21 January 2019 / Accepted: 29 January 2019 / Published: 2 February 2019
(This article belongs to the Section Physical Sensors)

Abstract

:
We propose a novel technique to estimate the total volume of unknown insulating inclusions in an electrically conducting body from voltage measurements. Unlike conventional Electrical Impedance Tomography (EIT) systems that usually exhibit low spatial resolution and accuracy, the proposed device is composed of a pair of driving electrodes which, supplied with a known sinusoidal voltage, create a current density field inside a region of interest. The electrodes are designed to generate a current density field in the region of interest that is uniform, to a good approximation, when the inclusions are not present. A set of electrodes with a polygonal geometry is used for four-wires resistance measurements. The proposed technique has been tested designing a low cost prototype, where all electrodes are on the bottom of the conducting body, showing good performances. Such a device may be used to monitor the volume of biological cells inside cell culture dishes or the volume of blood clots in micro-channels in lab-on-a-chip biosensors.

1. Introduction

There is an increasing interest in the measurement of the volume of inclusions growing inside an investigation area, with lots of applications in engineering [1,2,3]; in particular, in this paper we focus on its application in biomedical engineering and lab-on-a-chip biosensors. An appealing method to implement the volume estimation is to use the electrical impedance or voltage measurements which proved to be faster and more cost effective than techniques based on confocal microscopy.
Unlike conventional EIT systems that exhibit low spatial resolution and accuracy [4,5,6,7,8,9,10,11,12,13], we are concerned here only in the estimation of the total volume of unknown insulating inclusions in an electrically conducting body from voltage measurements. In other words, we are not interested in the position and shape of the inclusions and we assume that the electrical conductivity of the inclusions is negligible with respect to the conductivity of the conducting body.
The considered problem is typical in biomedical engineering when, for example, an automated technique to monitor the cell growth in a dish culture is needed. Another application in the medical field is the real time estimation of the blood clots’ volume in flowing blood for studying the thrombotic process in humans [14,15,16]. Either the cell membranes or the blood clot can be considered as insulating inclusions in the hundreds of kilohertz range of frequencies. For applications inside lab-on-a-chip biosensors, like the thrombotic profiling, the physical feasibility of the device inside openable micro-channels is of particular interest [17].
In mathematical literature, there are few papers containing the derivation of theoretical bounds for the inclusion volume from one measurement taken from a pair of electrodes [18,19,20]. While these bounds are very important from the theoretical point of view, in practice they are not very useful mainly because these bounds depend on unknown parameters that render the bounds not computable. In [21], a method based on simulation to obtain the previously introduced theoretical bounds is proposed. The results in [21] show that the obtained bounds are not very tight and therefore do not allow an accurate evaluation of the volume of the inclusion.
The solution adopted in this paper is strongly inspired by the technique proposed in [22,23,24,25]. The appealing feature of this technique is that no forward or inverse numerical solver is needed to compute the volume of the inclusion, allowing its real time estimation. The method determines the volume of the inclusion by applying a predefined pattern of injection currents on a number of electrodes placed in the lateral surface of a cylindrical container. Then, the voltages on these electrodes are sampled with two wire measurements. Finally, the size of the volume is simply estimated by summing for all electrodes the product of the injected current and the voltage readout.
We extend the technique in [25] in various ways. First of all, to estimate the inclusion volume in a three-dimensional container that represents the conductive body, we place all the electrodes only on the bottom wall of the cylindrical container. This planar configuration of electrodes enables the determination of the volume of inclusions flowing in microchannels of lab-on-a-chip devices with low cost prototyping; it would be more expensive to place the electrodes on the lateral walls of the microchannel. Moreover, we propose a technique to avoid the injection of a predefined pattern of currents in the electrodes that requires one independent current source for each electrode and a precise calibration of the value of each current. We avoid this complication by introducing a pair of driving electrodes in which a single current is injected [26]. Voltages are sampled on another set of electrodes disposed on a circle in the bottom of the container in such a way that voltages are measured with a four-wires technique. This has also the advantage of avoiding the so-called double layer at the electrode-fluid interface and therefore does not require the use of the complete model for electrodes [21]. The estimation of the inclusions volume is hence obtained using a formula derived from a simple power balance in the region of interest.
The physical realization of the prototype shows good performances in terms of linearity of the estimated volume with respect to the volume of the inclusions. Its real time operation and the limited cost render it appealing for the estimation of the total volume of inclusions for various practical applications.
The paper is structured as follows. Section 2 describes the sensor, in particular the geometric shape of the driving and the sensing electrodes, and the signal conditioning circuitry. Moreover, the theory behind its operation is explained by using the electrical power balance in the region of interest. In Section 3, for the sake of comparison only, a commercial simulation software (COMSOL MultiPhysics) is used to reproduce the behavior of the sensor. In Section 4, the performances of the sensor are compared to the ones provided by the numerical simulation. Finally, in Section 5, the conclusions are drawn.

2. Material and Methods

2.1. Sensor Description

The developed sensor and conditioning circuitry have been integrated in a single Printed Circuit Board (PCB) with gold plated electrodes, as shown in Figure 1.
The sensing element (shown in Figure 1b) is formed by: a couple of outer, parallel, wide driving electrodes which are responsible for supplying the sensor providing a current density from the top to the bottom of the sensor; an inner, polygonal set of reading electrodes which perform the readout of the potentials in the investigation area. In the present solution, the polygonal shape is a decagon having a width and height of 10 mm. The driving electrodes are 18 mm wide and 600 µm high, spaced 16 mm apart. The reading electrodes are 1.6 mm wide and 600 µm high posed as a decagon centered with respect to driving electrodes. The sensor forms the bottom wall of a circular reservoir having 25 mm diameter.
The schematic of the conditioning circuitry is shown in Figure 2a. A sinusoidal voltage V d having amplitude 100 mV and frequency 100 kHz is generated by a Wien bridge oscillator with amplitude control; the characterization of the oscillator evidenced a total harmonic distortion (THD) lower than 2% and a stability in the order of 0.3%. The amplitude of V d has been chosen to avoid red-ox reactions at the driving electrodes (which would yield to electrodes dissolution), while the frequency of V d has been chosen in order to minimize the double layer effect due to metal-electrolyte interface. Regarding the amplitude in fact, we chose | V d | = 100 mV, one order of magnitude lower with respect to the half-cell standard potential of gold (which is 1.4 V) to assure that no red-ox reaction takes place, thus avoiding the damage of the current electrodes. The frequency of V d has been chosen in order to minimize the double layer effect; when a metal transfers charge to an electrolyte, there is the formation of a narrow metal-electrolyte interface, named double layer, which behaves as a constant phase element (CPE). Its impedance is Z C P E ( ω ) = ( j ω ) ϕ C P E 1 and behaves partially as a resistor (if ϕ = 0 ) and partially as a capacitor (if ϕ = 1 ); this leads to a phase shift from voltage to current whose value is ϕ π / 2 . In previous works, we demonstrated that if the frequency is ≥100 kHz in the present application the double layer is negligible and current is in-phase with voltages. From the experimental characterization, the phase lag of current resulted lower than 1°; for this reason, in Section 2.3 we safely assume the measured impedance to be a resistance in order to obtain the volume information.
The lower driving electrode is connected to a current-to-voltage converter; its potential is fixed to virtual ground and the current I is converted into a voltage by means of a resistor. In the case of the present board, the resistor has been set to 330 Ω in order to obtain that the expected gain of the converter multiplied by the operating frequency results one decade lower than the gain-bandwidth product of the operational amplifier. The characterization of the current to voltage converter provided an uncertainty of conversion in the order of 1%. Each voltage electrode is connected, through a passive first order high-pass filter, to a buffer; the passive filter has input impedance 1.8 MΩ and cut-off frequency at 4 kHz. The aim of the filter is to remove the undesired DC component due to the electrolyte and the power line noise at 50 Hz. The buffer decouples the sensor voltages from the input impedance of the Data Acquisition (DAQ) board connected to the output of the circuit. In this way we minimize the load error during acquisition.

2.2. Data Acquisition

A 16 bits DAQ board is connected to the circuit (National Instruments NI-USB 6353); the DAQ has been set to acquire, at rate 1 MSa/s, a number of 3000 samples (corresponding to 300 periods of the driving voltage) of each differential voltage u i = V i V 0 , where V i is the potential of each electrode surface S i and V 0 is the potential of the lower electrode, chosen as reference, Figure 2a. The current signal I flowing through the outer electrodes is acquired as a single ended trace referred to ground. A dedicated control panel, shown in Figure 3, developed in LabView, acquires the data of the driving current and of u i . During a measurement cycle, the DAQ reads the nine differential voltages ( u i ) and the current signal (I). The measurement time for each channel is 3 ms, thus the acquisition time is 30 ms. After acquisition, the control panel calculates the AC-coupled RMS value of each quantity and saves the results into a text file. The volume of the defects is then computed with a Matlab script which calculates the quantity Δ R ( D ) as described in the next subsection.

2.3. Modelling and Volume Reconstruction Formula

Since (from measurements) voltages u i and current I are in phase sinusoidal signals, a stationary current conduction model is assumed in the domain of interest Ω , Figure 2b,c. At first, a reference configuration is considered in Ω without any insulating domain D, Figure 2b. We denote with I the current in this case, yielding a uniform current density vector J in Ω of amplitude
J = I S
where S denotes the effective cross-section area located between the current electrodes; its precise geometrical extent is not crucial, since S plays the role of an instrumental constant. For example, dashed line in Figure 2a represents the trace of a possible surface S. The electric power P i n entering Ω from its boundary Ω , can be expressed as
P i n = 1 9 u k i k ,
where the current i k on Ω is
i k = J · n k S e = J S e cos θ k ,
where S e is the effective area of the electrodes and θ k is the angle between J and the normal n k to S e , with k = 1 , , 9 . Of course P i n underestimates the actual electrical power flowing through Ω , due to the discrete distribution of the electrodes on Ω . By substituting (3) in (2), we obtain
P i n = J S e 1 9 u k cos θ k
Now, we apply the following approximate power balance in Ω
P i n = R ( Ω , ρ , J ) I 2 ,
where the functional
R ( Ω , ρ , J ) = 1 I 2 Ω ρ J 2 d Ω
is the equivalent resistance of Ω , ρ being the resistivity of the medium filling Ω . By combining (1), (4) and (5), we obtain
R ( Ω , ρ , J ) = S e S 1 9 R k cos θ k ,
where
R k = u k I .
Expression (7) yields an estimate of R ( Ω , ρ , J ) from voltages on Ω and the current I; the quantities R k are deduced directly from resistance measurement at the specified frequency. Moreover, the ratio S / S e plays the role of instrumental constant; the extent of S and S e can be estimated by means of numerical simulations.
Next, we introduce in Ω a number of insulating subregions whose union will be denoted with D, Figure 2c; thus, the new conducting domain of resistivity ρ becomes Ω D , and the corresponding resistance functional becomes
R ( Ω D , ρ , J ) = 1 I 2 Ω D ρ J 2 d Ω ,
where I is the current flowing through the outer electrodes due to the presence of the insulating domain D.
Applying a power balance to Ω D similar to (7), we obtain an estimate of the equivalent resistance of Ω D from voltages u k on Ω
R ( Ω D , ρ , J ) = S e S 1 9 R k cos θ k ,
where
R k = u k I .
Again, the quantities R k are deduced directly from resistance measurement at the specified frequency.
Finally, from (7) and (10), we compute the quantity
Δ R ( D ) = R ( Ω D , ρ , J ) R ( Ω , ρ , J ) = S e S 1 9 ( R k R k ) cos θ k ,
which is proportional to volume of D since the resistivity of the filling conductive medium ρ is an invariant.

3. FEM Simulations

The 3D geometry has been simulated using Comsol Multiphysics with a stationary conduction definition, since all the signals are in-phase. On the circumference that delimits the region we set the boundary condition of insulation, i.e., what happens in experimental setup having a plastic reservoir. We set the upper electrode at potential 100 mV and the lower electrode to ground; we evaluate V i on the voltage electrodes and I as the surface integral of current densities on the boundaries of the upper electrode. The subdomain of water has been set 3 mm high (as in experimental conditions) and its conductivity has been iteratively trimmed to match the current measured in the experiments with no defects in the investigation area. Inside the decagon, we have posed 18 spheres having 1.5 mm diameter leaning on the electrodes plane; their conductivities can assume the saline water value (in the order of 4 S/m) or the glass value (in the order of 10 −7 S/m) depending on the simulation setup.
The meshed geometry consists of 20,577 s order tetrahedra. A simulation script performed 10 series of simulations with random permutation of the spheres conductivities; each series is composed by 18 simulations (for a total amount of 180 simulations), where all the spheres are, progressively one by one, switched from water to glass. In this way, we simulate (with the same meshed geometry) 10 different combinations of possible linear growths of the spheres volume in the investigation area.
Figure 4a shows the simulated current density when all the spheres have the same conductivity of water, while Figure 4b shows the current density distribution when all the spheres have conductivity of glass; it is evident in Figure 4b that the current density avoids the glass insulating spheres, while no discontinuity is manifested in Figure 4a.

4. Experimental Results

In this section, we show the experimental results obtained from measurements and we compare them to FEM simulations. The reservoir of the sensor has been filled with 1.5 × 10−6 m3 (corresponding to 3 mm height) saline solution and we performed 10 measurement series to evaluate the linearity of the overall system. Each measurement series started from the acquisition of the voltages and current with no defects in the reservoir; then, progressively one by one, 18 defects (glass spheres with 1.5 mm diameter) have been added into the solution and, for each additional defect, we measured the respective voltages and current. Figure 5 shows the differential potentials at each electrode u i measured without defects in the decagon; markers represent the acquired experimental data, while the line represents the simulation results.
Figure 6 shows the trend of Δ R ( D ) in (12) for experimental results ( ), simulations (□), and its least squares linear regression (—). As it can be seen, depending on the sequence of insertion of the spheres in the investigation area, there is a small deviation from linearity both in experimental and in simulations.
For this reason, we are interested in evaluating the non-linearity of the sensor and comparing it to the non-linearity manifested in simulations. Evaluating the difference between data and least squares linear regression, we obtain the punctual linearity error of each j th volume acquisition ( v j k ) obtained from the k th series with respect to the least squares linear interpolation. The linearity error thus is evaluated as the difference between linear regression and the mean value v j ¯ of each point among the 10 series of measurements; the uncertainty on the linearity error u ( L i n ) j is thus represented by the standard deviation of the sample means:
u ( L i n ) j = k = 1 10 ( v j k v j ¯ ) 2 10 · 9
Figure 7 shows the plots of linearity error for simulations (top graph) and experimental data (bottom graph); error bars represent the standard deviation of the sample means u ( L i n ) . The results are expressed in terms of number of microspheres, and it is possible to see that both simulations and experimental results agree, showing 0.2 spheres non-linearity. This value corresponds, knowing that each sphere is 1.77 × 10−9 m3, to a non-linearity of 3.53 × 10−10 m3 (i.e., 1% of the measured volume). The agreement of non-linearity between simulated and experimental data is due to the approximation in (5); the accuracy of this approximation depends on the number of electrodes in the polygonal investigation area. With 10 voltage electrodes the non-linearity is in the order of 1%, this could be improved with a higher number of electrodes; obviously, increasing the number of electrodes means introducing more electronic components, a DAQ with more channels and longer acquisition time, thus invalidating the low-cost profile of the proposed solution.
The main difference that emerges from the two plots in Figure 7 is the uncertainty on non-linearity; the error bars in the bottom graph are, predictably, wider with respect to simulations, providing 0.4 sphere uncertainty (corresponding to 7 × 10−10 m3). This higher uncertainty is mainly due to the experimental setup, in particular to the oscillator stability and to the linearity of the current to voltage converter.

5. Conclusions

A technique to estimate the total volume of unknown insulating inclusions in a conducting body has been presented and experimentally validated. The technique turns out to be very cheap both in terms of hardware and software, given that no forward numerical solver solution is needed and expensive reconstruction algorithms are therefore avoided. Only a pair of driving electrodes is used that enable a simplification of the electronic hardware with respect to a similar technique proposed in literature that, in contrast, requires one independent current source for each electrode. The experimental results showed good agreement with simulations and the non-linearity of the proposed device results in the order of 1%.

Author Contributions

Conceptualization, A.A., F.T. and R.S.; methodology, A.A., F.T. and R.S.; software, A.A.; validation, A.A.; formal analysis F.T. and R.S.; investigation, A.A., F.T. and R.S.; resources, F.T. and R.S.; data curation, A.A.; writing—original draft preparation, A.A., F.T. and R.S.; writing—review and editing, A.A., R.S. and F.T.; visualization, A.A.; supervision, F.T.

Funding

This research received no external funding

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rymarczyk, T.; Klosowski, G.; Kozlowski, E. A Non-Destructive System Based on Electrical Tomography and Machine Learning to Analyze the Moisture of Buildings. Sensors 2018, 18, 2285. [Google Scholar] [CrossRef] [PubMed]
  2. Rymarczyk, T.; Tchorzewski, P.; Sikora, J. Implementation of Electrical Impedance Tomography for Analysis of Building Moisture Conditions. Compel Int. J. Comput. Math. Electr. Electron. Eng. 2018, 37, 1837–1861. [Google Scholar] [CrossRef]
  3. Rymarczyk, T.; Tchorzewski, P. Analysis of Historical Wall Dampness Using Electrical Tomography Measuring System. Int. J. Appl. Electromagn. Mech. 2018. [Google Scholar] [CrossRef]
  4. Yorkey, T.J.; Webster, J.G.; Tompkins, W.J. Comparing Reconstruction Algorithms for Electrical Impedance Tomography. IEEE Trans. Biomed. Eng. 1987, 34, 843–852. [Google Scholar] [CrossRef] [PubMed]
  5. Polydorides, N.; Lionheart, W.R.B. A Matlab toolkit for three-dimensional electrical impedance tomography: A contribution to the Electrical Impedance and Diffuse Optical Reconstruction Software project. Meas. Sci. Technol. 2002, 13, 1871. [Google Scholar] [CrossRef]
  6. Henderson, R.P.; Webster, J.G. An Impedance Camera for Spatially Specific Measurements of the Thorax. IEEE Trans. Biomed. Eng. 1978, 25, 250–254. [Google Scholar] [CrossRef] [PubMed]
  7. Barber, D.C.; Brown, B.H. Applied Potential Tomography. J. Phys. E Sci. Instrum. 1984, 17, 723–733. [Google Scholar] [CrossRef]
  8. Barber, D.C.; Brown, B.H.; Freeston, I.L. Imaging Spatial distributions of resistivity using Applied Potential Tomography. Electron. Lett. 1983, 19, 93–95. [Google Scholar] [CrossRef]
  9. Cheney, M.; Isaacson, D.; Newell, J. Electrical Impedance Tomography. SIAM Rev. 1999, 41, 85–101. [Google Scholar] [CrossRef]
  10. Zhang, T.; Zhou, L.; Ammari, H.; Seo, J.K. Electrical Impedance Spectroscopy-Based Defect Sensing Technique in Estimating Cracks. Sensors 2015, 15, 10909–10922. [Google Scholar] [CrossRef]
  11. Avery, J.; Dowrick, T.; Faulkner, M.; Goren, N.; Holder, D. A Versatile and Reproducible Multi-Frequency Electrical Impedance Tomography System. Sensors 2017, 17, 280. [Google Scholar] [CrossRef] [PubMed]
  12. Russo, S.; Nefti-Meziani, S.; Carbonaro, N.; Tognetti, A. A Quantitative Evaluation of Drive Pattern Selection for Optimizing EIT-Based Stretchable Sensors. Sensors 2017, 17, 1999. [Google Scholar] [CrossRef] [PubMed]
  13. Dupré, A.; Mylvaganam, S. A Simultaneous and Continuous Excitation Method for High-Speed Electrical Impedance Tomography with Reduced Transients and Noise Sensitivity. Sensors 2018, 18, 1013. [Google Scholar] [CrossRef] [PubMed]
  14. Affanni, A.; Chiorboli, G.; Codecasa, L.; Cozzi, M.; Marco, L.D.; Mazzucato, M.; Morandi, C.; Specogna, R.; Tartagni, M.; Trevisan, F. A novel inversion technique for imaging thrombus volume in microchannels fusing optical and impedance data. IEEE Trans. Magn. 2014, 50, 1021–1024. [Google Scholar] [CrossRef]
  15. Affanni, A.; Specogna, R.; Trevisan, F. Ex vivo time evolution of thrombus growth through optical and electrical impedance data fusion. J. Phys. Conf. Ser. 2013, 459, 012016. [Google Scholar] [CrossRef]
  16. Affanni, A.; Specogna, R.; Trevisan, F. Combined electro-optical imaging for the time evolution of white thrombus growth in artificial capillaries. IEEE Trans. Instrum. Meas. 2013, 62, 2954–2959. [Google Scholar] [CrossRef]
  17. De Zanet, D.; Battiston, M.; Lombardi, E.; Specogna, R.; Trevisan, F.; Marco, L.D.; Affanni, A.; Mazzucato, M. Impedance biosensor for real-time monitoring and prediction of thrombotic individual profile in flowing blood. PLoS ONE 2017, 12, e018494. [Google Scholar] [CrossRef]
  18. Rosset, E.; Alessandrini, G. The Inverse Conductivity Problem with One Measurement: Bounds on the Size of the Unknown Object. SIAM J. Appl. Math. 1998, 58, 1060–1071. [Google Scholar] [CrossRef]
  19. Alessandrini, G.; Rosset, E.; Seo, J. Optimal size estimates for the inverse conductivity problem with one measurement. Proc. Am. Math. Soc. 2000, 128, 53–64. [Google Scholar] [CrossRef]
  20. Alessandrini, G.; Rosset, E. Volume bounds of inclusions from physical EIT measurements. Inverse Probl. 2004, 20, 575. [Google Scholar] [CrossRef]
  21. Alessandrini, G.; Bilotta, A.; Morassi, A.; Rosset, E.; Turco, E. Computing Volume Bounds of Inclusions by EIT Measurements. J. Sci. Comput. 2007, 33, 293. [Google Scholar] [CrossRef]
  22. Kwon, O.; Seo, J.K. Total size estimation and identification of multiple anomalies in the inverse conductivity problem. Inverse Probl. 2001, 17, 59. [Google Scholar] [CrossRef]
  23. Kwon, O.; Seo, J.K.; Yoon, J.R. A real time algorithm for the location search of discontinuous conductivities with one measurement. Commun. Pure Appl. Math. 2002, 55, 1–29. [Google Scholar] [CrossRef]
  24. Kwon, O.; Seo, J.K.; Woo, E.J.; Yoon, J.R. Electrical impedance imaging for searching anomalies. Commun. Korean Math. Soc. 2001, 16, 459–485. [Google Scholar]
  25. Kwon, O.; Yoon, J.R.; Seo, J.K.; Woo, E.J.; Cho, Y.G. Estimation of anomaly location and size using electrical impedance tomography. IEEE Trans. Biomed. Eng. 2003, 50, 89–96. [Google Scholar] [CrossRef] [PubMed]
  26. Trevisan, F.; Affanni, A.; Specogna, R.; De Marco, L.; Mazzucato, M.; Battiston, M. Apparatus for Analyzing the Process of Formation of Aggregates in a Biological Fluid and Corresponding Method of Analysis. U.S. Patent 20160047827A1, 18 February 2016. [Google Scholar]
Figure 1. (a) overview of the printed circuit board (PCB) with circuitry and sensor (bottom left), and (b) zoom of the sensing element with some defects.
Figure 1. (a) overview of the printed circuit board (PCB) with circuitry and sensor (bottom left), and (b) zoom of the sensing element with some defects.
Sensors 19 00637 g001
Figure 2. (a) Schematic of the board, where potentials V i of the sensing electrodes with i = 0 , , 9 are shown; moreover, the pair of outer current electrodes have applied voltage V d and current I. (b) Current conduction in domain Ω without insulating domain. (c) Current conduction in Ω D , where D is the insulating domain.
Figure 2. (a) Schematic of the board, where potentials V i of the sensing electrodes with i = 0 , , 9 are shown; moreover, the pair of outer current electrodes have applied voltage V d and current I. (b) Current conduction in domain Ω without insulating domain. (c) Current conduction in Ω D , where D is the insulating domain.
Sensors 19 00637 g002
Figure 3. Control panel, developed in LabView, to control the sensor board.
Figure 3. Control panel, developed in LabView, to control the sensor board.
Sensors 19 00637 g003
Figure 4. Current density vector simulated with Comsol Multiphysics; arrows represent the current density direction and the surface color represents the magnitude. (a) When all the spheres are supposed to be water, (b) when all the spheres are supposed to be glass.
Figure 4. Current density vector simulated with Comsol Multiphysics; arrows represent the current density direction and the surface color represents the magnitude. (a) When all the spheres are supposed to be water, (b) when all the spheres are supposed to be glass.
Sensors 19 00637 g004
Figure 5. Differential voltage at i th electrode u i = V i V 0 ; line represents the numerical simulation, markers represent experimental results on 10 different measurements without defects in the sensor.
Figure 5. Differential voltage at i th electrode u i = V i V 0 ; line represents the numerical simulation, markers represent experimental results on 10 different measurements without defects in the sensor.
Sensors 19 00637 g005
Figure 6. Evaluation of the quantity Δ R ( D ) over 10 series of experiments: ( ) Simulation results, (□) Experimental results, (—) least squares fitting.
Figure 6. Evaluation of the quantity Δ R ( D ) over 10 series of experiments: ( ) Simulation results, (□) Experimental results, (—) least squares fitting.
Sensors 19 00637 g006
Figure 7. Evaluation of the non-linearity on the measurement expressed in terms of number of spheres (one sphere is 1.77 × 10−9 m3). Top graph represents the simulations, bottom graph represents experimental results. Error bars are the uncertainty expressed as standard deviation of the sample means.
Figure 7. Evaluation of the non-linearity on the measurement expressed in terms of number of spheres (one sphere is 1.77 × 10−9 m3). Top graph represents the simulations, bottom graph represents experimental results. Error bars are the uncertainty expressed as standard deviation of the sample means.
Sensors 19 00637 g007

Share and Cite

MDPI and ACS Style

Affanni, A.; Specogna, R.; Trevisan, F. Estimating the Volume of Unknown Inclusions in an Electrically Conducting Body with Voltage Measurements. Sensors 2019, 19, 637. https://doi.org/10.3390/s19030637

AMA Style

Affanni A, Specogna R, Trevisan F. Estimating the Volume of Unknown Inclusions in an Electrically Conducting Body with Voltage Measurements. Sensors. 2019; 19(3):637. https://doi.org/10.3390/s19030637

Chicago/Turabian Style

Affanni, Antonio, Ruben Specogna, and Francesco Trevisan. 2019. "Estimating the Volume of Unknown Inclusions in an Electrically Conducting Body with Voltage Measurements" Sensors 19, no. 3: 637. https://doi.org/10.3390/s19030637

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