Open AccessArticle

Analysis of the Dynamic Sensitivity of Hemisphere-Shaped Electrostatic Sensors’ Circular Array for Charged Particle Monitoring

Xin Tang

Zhong-Sheng Chen

^*,

Yue Li

Zheng Hu

and

Yong-Min Yang

Science and Technology on Integrated Logistics Support Laboratory, National University of Defense Technology, Changsha 410073, China

Author to whom correspondence should be addressed.

Sensors 2016, 16(9), 1403; https://doi.org/10.3390/s16091403

Submission received: 1 July 2016 / Revised: 26 August 2016 / Accepted: 27 August 2016 / Published: 31 August 2016

(This article belongs to the Section Physical Sensors)

Download

Browse Figures

Figure 1
Schematic of a grounded pipeline-installed HSESCA. "> Figure 2
Schematic of a proportional two-stage signal conditioner channel. "> Figure 3
Schematic of a grounded pipeline-installed sensor unit of the HSESCA. "> Figure 4
Schematic of bi and βi in the new Descartes coordinate system. "> Figure 5
(a–h) Theoretical array signals of the HSESCA from No. 1–No. 8 sensor units. "> Figure 6
Signal flow diagram of the i-th sensor unit of the HSESCA. "> Figure 7
Signal flow diagram of the HSESCA. "> Figure 8
Flow diagram of the component extraction-based array signal processing algorithm. "> Figure 9
The original charged particle q and its imaginary point charges. "> Figure 10
Grid points in the observation cross-section. "> Figure 11
(a–h) Surface diagrams of dynamic sensitivity of the sensor units from No. 1–No. 8. "> Figure 12
Contour diagram of dynamic sensitivity of the No. 7 sensor unit. "> Figure 13
(a) Surface and (b) contour diagrams of the static sensitivity of the HSESCA. "> Figure 14
(a) Surface and (b) contour diagrams of the dynamic sensitivity of the HSESCA. "> Figure 15
Structure of the HSESCA experiment apparatus. "> Figure 16
Test points set in the experiment. "> Figure 17
Array signals of the HSESCA. (a) Released at P1 with −32 pC charge; (b) released at P2 with −20 pC charge; (c) released at P3 with −46 pC charge; (d) released at P4 with −25 pC charge; (e) released at P5 with −29 pC charge. "> Figure 18
Relationships between the absolute charge quantities and the absolute peak voltages of the No. 1 sensor unit. "> Figure 19
Experimental and theoretical values of static sensitivity (a) in Line 1 and (b) in Line 2. "> Figure 20
Experimental and theoretical values of dynamic sensitivity (a) in Line 1 and (b) in Line 2. ">

Versions Notes

Abstract

Electrostatic sensor arrays (ESAs) are promising in industrial applications related to charged particle monitoring. Sensitivity is a fundamental and commonly-used sensing characteristic of an ESA. However, the usually used spatial sensitivity, which is called static sensitivity here, is not proper for moving particles or capable of reflecting array signal processing algorithms integrated in an ESA. Besides, reports on ESAs for intermittent particles are scarce yet, especially lacking suitable array signal processing algorithms. To solve the problems, the dynamic sensitivity of ESA is proposed, and a hemisphere-shaped electrostatic sensors’ circular array (HSESCA) along with its application in intermittent particle monitoring are taken as an example. In detail, a sensing model of the HSESCA is built. On this basis, its array signals are analyzed; the dynamic sensitivity is thereupon defined by analyzing the processing of the array signals. Besides, a component extraction-based array signal processing algorithm for intermittent particles is proposed, and the corresponding dynamic sensitivity is analyzed quantitatively. Moreover, simulated and experimental results are discussed, which validate the accuracy of the models and the effectiveness of the relevant approaches. The proposed dynamic sensitivity of ESA, as well as the array signal processing algorithm are expected to provide references in modeling, designing and using ESAs.

Keywords:

electrostatic sensor array; dynamic sensitivity; hemisphere-shaped electrostatic sensor; circular array; array signal processing algorithm

1. Introduction

Electrostatic monitoring systems (EMSs) are featured with the advantages of robustness and low cost [1,2], making them very promising in parameter monitoring of gas-solid flows [3,4,5], exhaust debris monitoring-based PHM of gas turbines [1,6,7], etc. Electrostatic sensors are key components of EMSs. Among them, hemisphere-shaped electrostatic sensors were proposed by us as a trade-off between intrusive and non-intrusive electrostatic sensors. They combine the advantages of relatively high sensitivity, easy installation, flexible layout, accuracy in theoretical modeling and cause very little disturbance to flows [8,9].

The sensing principle of an electrostatic sensor is electrostatic induction. It determines that even if a particle carries a constant charge, the corresponding induced charge on an electrostatic sensor probe will reduce sharply when the distance between the particle and the probe gets larger. As a result, a near particle with quite weak charge can generate similar signals as a far particle with great charge. That is to say, the charge quantity of a particle can hardly be accurately monitored by a single electrostatic sensor if the position of the particle is not provided. In view of the sensitivity, this means that most electrostatic sensors have the drawback of inhomogeneous and quite localized sensitivity; thus, only particles near the probe can be effectively detected [8,10]. In a word, as for the applications that need the accurate particle charge (e.g., particle size monitoring or gas-solid flow density monitoring), an EMS that contains a single electrostatic sensor is usually of poor monitoring accuracy. Moreover, a total failure of the whole system will occur if only the sensor breaks down. Therefore, in order to obtain better monitoring accuracy and reliability, EMSs composed of multiple electrostatic sensors have become the trend [5,11,12]. They are usually called electrostatic sensor arrays (ESAs), and the electrostatic sensors are called the sensor unit. Two aspects are of high concern in relevant research: the first is the hardware design and related sensing characteristics analysis; the second is the application and related array signal processing algorithms.

Circular ESA is a kind of representative hardware structure, which is made up by placing some identical electrostatic sensors uniformly around a circular pipeline in one cross-section [12,13,14,15,16,17,18]. Xu et al. used arc-sharped electrostatic sensors to construct circular ESAs and used the finite element method to study the sensing characteristics of the sensor units. Some geometric and material influence factors on the sensing characteristics were also analyzed [18]. Thuku et al. studied circular ESAs composed of rod-shaped electrostatic sensors [5,19]. A sensing model of the sensor units was built, and a comparison between four and 16 sensor systems was made [17]. In addition, differently-sized circularplate-shaped and rectangular plate-shaped electrostatic sensors were modeled and analyzed [20]; they are also suitable for circular ESAs. However, existing works on the sensing characteristics were generally focused on individual electrostatic sensors, thus lacking descriptions and models for the sensing characteristics of a whole ESA. In addition, as a fundamental and the most commonly-used sensing characteristic, sensitivity should reflect the monitoring accuracy of an EMS directly. As for an ESA, this means that a proper sensitivity definition should build a direct relationship between the monitored particles and the monitoring parameters of the whole ESA. However, the usually used spatial sensitivity only builds a basic relationship between static charged particles and corresponding induced charge on a sensor probe [10,18,21]; thus, it is called static sensitivity here. In practice, particles to be monitored are always moving. Furthermore, as a three-dimensionally-defined parameter, static sensitivity is difficult to accurately model when real boundary conditions are considered [9]. More importantly, follow-up processing of the induced charge is not considered in the definition, such as charge-voltage transformation and array signal processing algorithms. In a word, in order to build a direct relationship between the monitored particles and the monitoring parameters of a whole ESA, thus reflecting its monitoring accuracy directly, a more practical sensitivity definition needs to be defined in a systemic perspective.

In the aspect of application, a circular ESA is usually used to detect the particle distribution over the cross-section of a pipeline and further infer gas-solid flow parameters, such as particle density and flow regime. In these cases, the particles to be monitored are continuous in the time domain. That is to say, there are always numerous particles passing a cross-section simultaneously, thus forming an obvious profile of the solid-phase, which can be imaged by tomography-based methods [5,11,14]. ESAs with arc-shaped sensor units usually adopt array signal processing algorithms based on electrical capacitance tomography [11,12,16,22,23]. They are suitable for dense gas-solid flows. There are also reports on volume tomography using multi-circular ESAs with arc-shaped sensor units [24,25,26]. ESAs with rod-shaped sensor units usually adopt algorithms based on electrostatic tomography [5,14,18]. They are suitable for both dense and dilute gas-solid flows. However, because of the limited number of sensor units, the tomography-based algorithms are faced with tough ill-conditioned problems, which lead to poor image resolution [11,16]. Improvements can be made by increasing the unit number [17] and adopting some iterative regularization methods, filtering methods, error minimization protocols, etc., which will increase the computation time [11,12,16]. As a result, the tomography-based algorithms only indicate approximate areas where particles concentrate, but they cannot accurately confirm where discrete particles are located. Moreover, the algorithms are more suitable for relatively stable flows, because the monitoring accuracy can be improved by making a compromise on the computation time. In contrast, particles to be monitored could be quite intermittent in some applications. That is to say, there are only one or just a few particles passing a cross-section simultaneously. One representative application is exhaust debris monitoring-based PHM of gas turbines [7,17,18,19,20,21,22,23,24,25,26,27,28,29]. In this case, the occurrences of fault-related particles are accidental and instantaneous. In other words, first, the particles to be monitored are quite intermittent. Second, a fast array signal processing algorithm is needed to avoid missing the fast-passing particles. Besides, the size and charge quantity of the particles, rather than the distribution, are concerned. A measurement technique was studied by Addabbo et al. to estimate the trajectory position, charge and velocity of a moving charged particle. It was based on a small array whose electrostatic sensors have improved low-frequency response [30,31]. However, relevant studies are scarce yet. As for a circular ESA, array signal processing algorithms remain to be designed for the intermittent particle monitoring.

To solve the problems above, the dynamic sensitivity of ESA is proposed to describe the sensitivity of ESAs in a systemic perspective. An ESA called the hemisphere-shaped electrostatic sensors’ circular array (HSESCA) along with its application in intermittent particle monitoring are taken as an example. In detail, a sensing model of the HSESCA, as well as that of the sensor units is built in Section 2. On this basis, array signals of the HSESCA are analyzed, and the dynamic sensitivity of ESA is thereupon defined in Section 3. A component extraction-based array signal processing algorithm for discrete particle monitoring is designed in Section 4. Numerical simulations and experiments are done for validations in Section 5. The main results are concluded in Section 6.

2. Sensing Model of an HSESCA

This section builds a quantitative relationship between a charged particle and the induced charge of an HSESCA, which provides a theoretical foundation for this paper.

2.1. Basic Structure of an HSESCA

An HSESCA with eight sensor units is installed on a grounded pipeline, as shown in Figure 1. It is mainly composed of eight identical hemispherical probes, a multi-channel signal conditioner and an array signal analyzer. Compared to a “flat head” probe leaning inward toward the pipeline, it is much easier to build a theoretical sensing model for a hemispherical probe, because it only contains a simple hemispheric surface. In addition, a hemispherical probe should be more durable than a “flat head” probe under severe working conditions, as a hemispheric surface is free of damageable edges.

The hemispherical probes are uniformly installed in a cross-section around the pipeline. It is called the observation cross-section. Once the charged particles pass, the induced charge will be generated on each probe and then converted into voltage array signals immediately by the multi-channel signal conditioner for post-processing. Every channel of the conditioner is designed as an identical proportional two-stage amplifier as in Figure 2 [9]. The induced charge is proportionally transformed into voltage in the first stage and then negatively amplified in the second stage. Therefore, the total function of each channel is represented by a negative amplification coefficient K whose unit is V/C. A hemispherical probe and its signal conditioner channel compose an independent sensor unit of the HSESCA. The array signal analyzer is composed of a signal acquisition card and a computer. Some specially-designed array signal processing algorithms are usually integrated in the computer to calculate the monitoring parameters (denoted as M) of the whole HSESCA.

As the sensitivity characteristics of an electrostatic sensor are quite localized [8,9,10] and the radius of the probes is much smaller than the distances between them, interactions between the probes can be ignored. Thus, a sensor unit is firstly modeled to lay a foundation for the whole HSESCA.

2.2. Theoretical Modeling of a Hemisphere-Shaped Electrostatic Sensor Unit

As shown in Figure 3, one sensor unit is kept. A Descartes coordinate system is built with the center of the probe’s bottom face set as the origin. The axial direction of the pipeline and that of the probe are set as the x-axis and the z-axis, respectively. Then, the y-axis is determined by the right-hand rule automatically. Accordingly, a point in the pipeline is denoted as P(x, y, z).

Because electrostatic equilibrium states are reached instantaneously (10⁻¹⁹ s), the interaction between the probe and charged particles in the pipeline can be described by a pure electrostatic field [18,21]. It is uniquely determined by the Poisson equation and Dirichlet boundary conditions as:

\nabla^{2} φ_{t} (P) = - \frac{ρ_{t} (P)}{ε}, s . t . {\begin{cases} φ_{t} (P) = 0, P \in Γ_{2} \\ φ_{t} (P) = const (t), P \in Γ_{1} \\ \lim_{P \to \infty} φ_{t} (P) = 0 \end{cases}

(1)

where Γ₁ is the hemispherical surface of the probe and Γ₂ is the inner wall of the pipeline. ρ_t(P) is the volume charge density in the pipeline at time t, and φ_t(P) is the corresponding potential distribution. Besides, const(t) is a constant function of t, describing the equipotential surface of the probe. ε is the permittivity of free space, and

\nabla^{2}

is the Laplace operator.

On account of the pipeline being much bigger than the probe, the electrostatic field on the inner wall is simplified as that on one infinite plane [2,6,17]. In this case, when a charged particle is simplified as a point charge q [8,10], φ_t(P) is firstly solved using the Green function and the method of image charges. Then, the sensing model at P(x, y, z) is obtained according to the relationship between the induced charge and the potential distribution on a conductor surface [8,9]. That is:

\begin{array}{l} Q (x, y, z) = - \frac{q a}{b} - \frac{q z}{2 π} \int_{0}^{2 π} d θ \int_{0}^{a} {(r^{2} + b^{2} - 2 x y \cos θ - 2 y r \sin θ)}^{- \frac{3}{2}} r d r \\ + \frac{q a^{3} z}{2 π b^{3}} \int_{0}^{2 π} d θ \int_{0}^{a} {[r^{2} + \frac{a^{4} - 2 a^{2} (x r \cos θ + y r \sin θ)}{b^{2}}]}^{- \frac{3}{2}} r d r \end{array}

(2)

where Q is the total induced charge on the probe, a is the radius of the probe and

b = \sqrt{x^{2} + y^{2} + z^{2}}

is the distance between the point charge and the origin.

2.3. Theoretical Modeling of an HSESCA

In order to build a sensing model for all of the sensor units in the same coordinate system, firstly, according to the symmetry of the probe, Equation (2) is transformed into:

\begin{array}{l} Q (x, y, z) = - \frac{q a}{b} - \frac{q b \cos β}{2 π} \int_{0}^{2 π} d θ \int_{0}^{a} {(r^{2} + b^{2} - 2 b r \sin β \sin θ)}^{- \frac{3}{2}} r d r \\ + \frac{q a^{3} \cos β}{2 π b^{2}} \int_{0}^{2 π} d θ \int_{0}^{a} {(r^{2} + \frac{a^{4}}{b^{2}} - \frac{2 a^{2} r \sin β \sin θ}{b})}^{- \frac{3}{2}} r d r \end{array}

(3)

where

β = \arccos (z / b)

is the included angle between the line from the point charge to the origin and the axis of the probe. Then, a new Descartes coordinate system is built by translating the original one (see Figure 2) to the center of the observation cross-section. Meanwhile, the sensor units are numbered anticlockwise, with coordinates of the i-th probe denoted as (0, y_i, z_i) and those of the point charge denoted as P(x, y, z) in the new coordinate system; as shown in Figure 4.

Thereupon, the sensing model of the i-th sensor unit is derived from Equation (3) as:

\begin{array}{l} Q_{i} (x, y, z) = - \frac{q a}{b_{i}} - \frac{q b_{i} \cos β_{i}}{2 π} \int_{0}^{2 π} d θ \int_{0}^{a} {(r^{2} + {b_{i}}^{2} - 2 b_{i} r \sin β_{i} \sin θ)}^{- \frac{3}{2}} r d r \\ + \frac{q a^{3} \cos β_{i}}{2 π {b_{i}}^{2}} \int_{0}^{2 π} d θ \int_{0}^{a} {(r^{2} + \frac{a^{4}}{{b_{i}}^{2}} - \frac{2 {a^{2}}_{i} r \sin β_{i} \sin θ}{b_{i}})}^{- \frac{3}{2}} r d r \end{array}

(4)

where b_i is the distance from the charged particle to the i-th probe and β_i is the included angle between the point charge and the axis of the i-th probe. They are determined by geometrical relationships:

b_{i} = \sqrt{x^{2} + {(y - y_{i})}^{2} + {(z - z_{i})}^{2}}

(5)

β_{i} = \arccos [({y_{i}}^{2} + {z_{i}}^{2} - y y_{i} - z z_{i}) / l b_{i}]

(6)

where

l = \sqrt{{y_{i}}^{2} + {z_{i}}^{2}}

is the distance between the origin O and the bottom face of a probe.

Further, the sensing model of the HSESCA is obtained by making a superposition of sensing models of all of the sensor units. That is:

\begin{matrix} Q_{A} (x, y, z) = \sum_{i} Q_{i} (x, y, z) = - \sum_{i} \frac{q a}{b_{i}} - \sum_{i} \frac{{q b}_{i} \cos β_{i}}{2 π} \int_{0}^{2 π} d θ \int_{0}^{a} {(r^{2} + {b_{i}}^{2} - 2 b_{i} r \sin β_{i} \sin θ)}^{- \frac{3}{2}} r d r \\ + \sum_{i} \frac{q a^{3} \cos β_{i}}{2 π {b_{i}}^{2}} \int_{0}^{2 π} d θ \int_{0}^{a} {(r^{2} + \frac{a^{4}}{{b_{i}}^{2}} - \frac{2 {a^{2}}_{i} r \sin β_{i} \sin θ}{b_{i}})}^{- \frac{3}{2}} r d r \end{matrix}

(7)

In addition, an FEM-based calibration method was proposed as a supplement to consider the actual boundary conditions [9]. It is used to improve the accuracy of the sensing models when a charged particle is in the observation cross-section. The main steps are as follows: First, an FEM model of a grounded pipeline-installed probe is built. The model provides simulated induced charge with high fidelity. Then, the relative error between the simulated and the theoretical induced charge is analyzed. On this basis, a calibration item calculated by using fitting methods is added to eliminate the relative error. When the inner radius of the pipeline is 200 mm and the button face of a probe is 2 mm apart away from the inner wall, Equation (4) is calibrated as:

Q_{i, x = 0}^{'} (y, z) = {\begin{cases} Q_{i, x = 0} (y, z) (p_{1} {b_{i}}^{p_{2}} + p_{3}) & , if b_{i} \geq D \\ Q_{i, x = 0} (y, z) & , if b_{i} < D \end{cases}

(8)

where D = −0.0188 β_i + 0.0444 is the calibration boundary that is determined according to the variation trend of the relative error [10]. In the area of b_i < D, the relative error is negligible; thus, the sensing model is only calibrated where b_i ≥ D. The calibration parameters are obtained using fitting methods; they are:

{\begin{cases} p_{1} = - 5.7311 \times 10^{- 8} e^{8.7639 β_{i}^{2}} - 6.103 e^{0.5724 β_{i}^{2}} \\ p_{2} = - 43.9996 e^{0.1261 β_{i}^{2}} + 45.9113 e^{0.113 β_{i}^{2}} \\ p_{3} = 1.0479 e^{7.0643 \times 10^{- 2} β_{i}^{2}} + 1.2218 \times 10^{- 4} e^{4.4059 β_{i}^{2}} \end{cases}

(9)

It is obvious that Equation (7) is calibrated as

Q_{A, x = 0}^{'} (y, z) = \sum_{i} Q_{i, x = 0}^{'} (y, z)

3. Definition of the Dynamic Sensitivity of ESA

3.1. Definition of the Static Sensitivity of ESA

In order to make a comparison with the dynamic sensitivity of ESA, the static sensitivity of ESA is defined here. The static sensitivity of an electrostatic sensor is defined as the absolute induced charge on the probe generated by a unit charge [7,8,21]. Analogously, the static sensitivity of ESA is defined as the absolute induced charge on all of the probes of an ESA generated by a unit charge. That is:

S_{A} (x, y, z) = - \frac{Q_{A} (x, y, z)}{q}

(10)

The minus indicates that the induced charge and the point charge q are of opposite signs. As for an HSESCA, its static sensitivity model is derived from Equations (9) and (10):

\begin{array}{l} S_{A} (x, y, z) = \sum_{i} \frac{a}{b_{i}} + \sum_{i} \frac{b_{i} \cos β_{i}}{2 π} \int_{0}^{2 π} d θ \int_{0}^{a} {(r^{2} + {b_{i}}^{2} - 2 b_{i} r \sin β_{i} \sin θ)}^{- \frac{3}{2}} r d r \\ - \sum_{i} \frac{a^{3} \cos β_{i}}{2 π {b_{i}}^{2}} \int_{0}^{2 π} d θ \int_{0}^{a} {(r^{2} + \frac{a^{4}}{{b_{i}}^{2}} - \frac{2 {a^{2}}_{i} r \sin β_{i} \sin θ}{b_{i}})}^{- \frac{3}{2}} r d r \end{array}

(11)

When actual boundary conditions are considered, the model in the observation cross-section can also be calibrated in a similar way as Equation (8), that is

S_{A, x = 0}^{'} (y, z) = - Q_{A, x = 0}^{'} (y, z) / q

3.2. Theoretical Array Signals of an HSESCA

Array signal processing algorithms are included in the definition of the dynamic sensitivity of ESA. In order to find effective array signal processing algorithms for an HSESCA, the characteristics of its array signal shave to be analyzed.

It is assumed that a charged particle moves along the pipeline’s axial direction at a constant velocity [10,32]. In this case, the x-coordinates of the particle are denoted with velocity v and time t, while the y-coordinates and z-coordinates are regarded as constants. In this way, induced charge on the i-th probe is expressed by Q_i(x₀ + vt, y, z), where x₀ is the initial x-coordinate. After that, the induced charge is transformed and amplified into voltage according to the amplification coefficient K of the signal conditioner (see Figure 4), thus obtaining the out-put signal of the i-th sensor unit. That is:

u_{i} (t) = K Q_{i} (x_{0} + v t, y, z)

(12)

Then, by combining Equations (4) and (12), one can get:

\begin{array}{l} u_{i} (t) = - \frac{q K a}{b_{i}} - \frac{{q K b}_{i} \cos β_{i}}{2 π} \int_{0}^{2 π} d θ \int_{0}^{a} {(r^{2} + {b_{i}}^{2} - 2 b_{i} r \sin β_{i} \sin θ)}^{- \frac{3}{2}} r d r \\ + \frac{q K a^{3} \cos β_{i}}{2 π {b_{i}}^{2}} \int_{0}^{2 π} d θ \int_{0}^{a} {(r^{2} + \frac{a^{4}}{{b_{i}}^{2}} - \frac{2 {a^{2}}_{i} r \sin β_{i} \sin θ}{b_{i}})}^{- \frac{3}{2}} r d r \end{array}

(13)

where

b_{i} = \sqrt{{(x_{0} + v t)}^{2} + {(y - y_{i})}^{2} + {(z - z_{i})}^{2}}

In order to provide a visualized view of the array signals, here, we set K = −1, q = 1C, and the time when the particle reaches the observation cross-section is set as the zero time (i.e., x₀ = 0). The radius of the hemispherical probes is set to 12 mm, and the motion path is set to y = 100 mm and z = 70 mm (see Figure 3). Besides, the velocities of the charged particle are set to 1 m/s, 2 m/s and 3 m/s, respectively. Then, the array signals of the HSESCA are calculated by using Equation (13), as shown in Figure 5.

This shows that all of the signal peaks u_max,i appear when the charged particle reaches the observation cross-section; when q is set as a constant, the i-th peak is just determined by the relative position between the motion path and the i-th sensor unit, but has nothing to do with the particle’s velocity. Moreover, it is obvious that the peaks will change if the position of the charged particle changes. In a word, different combinations of peaks of the array signals indicate different position information of the charged particle in the observation cross-section.

In order to build a quantitative and direct relationship between a charged particle and the out-put signal peaks of a single electrostatic sensor, the dynamic sensitivity of the electrostatic sensor was proposed by us. It is defined as the absolute value of the out-put signal peak of an electrostatic sensor when a unit point charge passes [9]. As for the i-th sensor unit, the relationship is built by its dynamic sensitivity model as:

S_{D, i} (y, z) = \frac{u_{\max, i}}{q} = \frac{K Q_{i, x = 0}^{'} (y, z)}{q}

(14)

where

Q_{i, x = 0}^{'} (y, z)

is shown as Equation (8).

It is worth noting that in many applications, the monitoring parameters are reflected by the charge quantity of the monitored particles directly or indirectly, such as the particle size, concentration, flow rate, etc. [4,33]. Therefore, it is significant to ensure that the charge quantity of any particle can be accurately monitored no matter where it crosses the observation cross-section. Accordingly, a reasonable way of processing the array signals is as follows: first, obtaining the position information of the monitored particles from the array signals; then, weighting the signal components induced by particles in the less sensitive area according to the position information. In this way, the charge quantity of any particle can be effectively reflected.

3.3. Definition and Interpretation of the Dynamic Sensitivity of ESA

The aim of proposing the dynamic sensitivity of ESA is to build a direct relationship between moving charged particles and the monitoring parameters of an ESA, thus reflecting its monitoring accuracy directly. “Dynamic” here has two implications: the first is making a distinction from static sensitivity, that is dynamic sensitivity is proposed for moving charged particles; the second is that the specific expression of the dynamic sensitivity of ESA is changeable according to different array signal processing algorithms. More detailed interpretations are provided by taking intermittent particle monitoring, for example, as follows.

In some applications, the particles to be monitored are intermittent. One representative application is exhaust debris monitoring-based PHM of gas turbines. When faults (e.g., fatigue or carbon deposition) occur on gas path components of gas turbines, some fault-related particles will be produced and charged. Thus, the condition information of the gas path components can be obtained by monitoring the particles [1,6,7,28]. In general, the fault-related particles are intermittent. Moreover, it is considered that larger particles carry greater charge, which indicates more serious faults [1,7]. Therefore, the charge quantity of any fault-related particle should be accurately monitored.

As for a single electrostatic sensor, for example the i-th sensor unit of the HSESCA, its dynamic sensitivity model builds a relationship between moving charged particles and its out-put voltage signals, as Equation (14). This is expressed by Figure 6.

This means that a charged particle contains two kinds of information, charge quantity q and position (y, z). They codetermine the signal peak u_max,i according to the unit’s dynamic sensitivity model S_D,i(y, z), that is u_max,i = qS_D,i(y, z). However, the relationship is irreversible in practical applications. This is because the position information is difficult to obtain by using a single electrostatic sensor; thus, the value of S_D,i(y, z) is not determined. As a result, the charge quantity q cannot be calculated backward for the lack of position information. This is the basic reason for the low monitoring accuracy of a single electrostatic sensor.

It has been mentioned that the array signals of an HSESCA contain the position information of a monitored particle (see Section 3.2), which can be used to improve the monitoring accuracy of the HSESCA. This is expressed by Figure 7.

This means that a charged particle contains two kinds of information, charge quantity q and position (y, z). They codetermine the signal peak u_max,i of each sensor unit according to its dynamic sensitivity model S_D,i(y, z). Then, the peaks are processed by a certain array signal processing algorithm to calculate the monitoring parameter M. In this process, the estimated position information (y’, z’) of the particle is usually obtained and used. Therefore, valves of S_D,i(y, z) are determined by the position information, making it possible to calculate the charge quantity q backward. In fact, the monitoring parameters of many applications are just the charge quantity or some derivative values of it. As a result, by using an HSESCA and corresponding array signal processing algorithms, the monitoring accuracy can be improved.

In order to provide a unified and convenient expression for different array signal processing algorithms, dynamic operator F_D is proposed to denote the array signal processing algorithms. It builds the relationship between the array signals and monitoring parameters of an ESA (see Figure 7).

Further, the dynamic sensitivity of ESA is defined in the observation cross-section as the absolute value of a monitoring parameter when a unit charge passes. As for intermittent particle monitoring, the dynamic sensitivity model of an HSESCA is expressed as:

S_{DA} (y, z) = | \frac{M}{q} | = | \frac{F_{D} [U_{\max} (y, z)]}{q} | = | \frac{F_{D} [q S_{D} (y, z)]}{q} | = | F_{D} [S_{D} (y, z)] |

(15)

where U_max(y, z) is the vector of the signal peaks and S_D(y, z) is the vector of dynamic sensitivity values of the sensor units. The q is reducible because F_D has to be designed as a linear operator of q.

It is seen from Equation (15) that the essence of the dynamic sensitivity of ESA is a recombination of the dynamic sensitivity of the sensor units according to a dynamic operator, which builds a direct relationship between moving charged particles and the system monitor values of an ESA. In other words, the dynamic sensitivity of ESA reflects the monitoring accuracy of an ESA directly by taking array signal processing algorithms into consideration. The concept of dynamic sensitivity has been used by Zhou et al. [34], but the physical meaning is quite different.

4. A Component Extraction-Based Array Signal Processing Algorithm for Intermittent Particles

A proper array signal processing algorithm is significant to ensuring the monitoring accuracy of an ESA. As for intermittent particle monitoring (e.g., exhaust debris monitoring-based PHM of gas turbines), it has been mentioned that the accurate charge quantity of any monitored particle is desired; thus, the position information of the particle should be used. In addition, the demand for real-time monitoring should be met to avoid missing momentary faults; thus, the corresponding array signal processing algorithms have to be fast enough.

To meet the conditions above, a component extraction-based array signal processing algorithm is designed. It is based on a simple idea that a monitored particle is located near the sensor units that produce larger signal peaks. According to this idea, a set of common components is extracted from the peaks of the array signals. The components are considered to be generated by a set of imaginary point charges. They have fixed positions, but their imaginary charge quantities vary with the position of the monitored particle. In this way, the position information of the monitored particle can be used according to different charge quantities of the imaginary point charges. It is worth noting that the accurate position of a particle is not needed in the algorithm, but it is made use of by introducing a component extraction-based method. As a result, the algorithm only contains some simple steps, which are convenient for computer processing. More details are provided as follows.

As shown in Figure 8, three main steps and one pre-step are contained in the component extraction-based array signal processing algorithm. The condition in Figure 4 is taken as an example below, that is to say, the HSESCA contains eight sensor units, and a charged particle moves along a motion path of y = 100 mm and z = 70 mm.

The first step is re-sorting. When a charged particle passes the observation cross-section, the corresponding signal peaks of the eight sensor units are recorded as a vector U_max. Then, U_max is re-sorted into U_r incrementally according to the absolute values of the signal peaks. That is to say, the i-th absolute smallest peak is denoted as u_ri; thus, the new vector is U_r{u_r₁, u_r₂, u_r₃, u_r₄, u_r₅, u_r₆, u_r₇, u_r₈}. For example, the No. 2 probe is the nearest one to the point charge, while the No. 6 probe is the farthest one to the point charge (see Figure 3); thus, according to the fact that nearer probes produce larger signal peaks, one will have u_r₁ = u_max,6, u_r₈ = u_max,2.

The second step is component extraction. It is clear that a monitored particle is located near the sensor units that produce larger signal peaks. Thus, if a signal component is only contained in one signal peak, its value implies how close the particle is to the corresponding sensor probe. Analogously, if a signal component is contained in every signal peak, its value implies how close the particle is to the central point. Accordingly, a set of common signal components is extracted from U_r by using Equation (16). In detail, the absolute smallest signal peak u_r₁ is firstly denoted as c₁; it is a common signal component contained in every signal peak. Then, the second absolute smallest peak u_r₂ is denoted as the sum of c₁ and c₂. It is obvious that c₂ is a common signal component contained in every signal peak, except u_r₁, and so on; a vector of eight common signal components is extracted from U_r, that is C{c₁,c₂, c₃, c₄, c₅, c₆, c₇, c₈}.

{\begin{cases} u_{r 1} = c_{1} \\ u_{r 2} = c_{1} + c_{2} \\ u_{r 3} = c_{1} + c_{2} + c_{3} \\ u_{r 4} = c_{1} + c_{2} + c_{3} + c_{4} \\ u_{r 5} = c_{1} + c_{2} + c_{3} + c_{4} + c_{5} \\ u_{r 6} = c_{1} + c_{2} + c_{3} + c_{4} + c_{5} + c_{6} \\ u_{r 7} = c_{1} + c_{2} + c_{3} + c_{4} + c_{5} + c_{6} + c_{7} \\ u_{r 8} = c_{1} + c_{2} + c_{3} + c_{4} + c_{5} + c_{6} + c_{7} + c_{8} \end{cases}

(16)

The third step is weighting and summation. First of all, the extracted signal components are considered to be generated by a set of imaginary point charges. In detail, based on the symmetry of the circular array and the idea that a charged particle is located near the sensor units that produce larger signal peaks, an imaginary point charge q₁ located in the central point is firstly assumed to generate c₁. Then, because c₂ is contained in every signal peak, except that of the No. 6 sensor unit, an imaginary point charge q₂ is located in the connection line between the origin, and the No. 2 probe is assumed to generate c₂. Analogously, as c₃ is produced by every sensor unit except the No. 5 and No. 6 sensor units, q₃ located in the bisector of the No. 1 and the No. 2 probes is assumed to generate c₃, and so on; eight imaginary point charges are decomposed from the charged particle q, as shown in Figure 9. Among them, the point charge q₈ should be quite close to the No. 2 probe, because its relevant signal component is only produced by that probe. The positions of the imaginary point charges are considered to be fixed due to the symmetry of the HSESCA. However, their imaginary charge quantities, which are determined by the values of the extracted components, vary with the position of the monitored particle. In other words, the charge quantities of the imaginary point charges imply the position information of the monitored particle.

Next, the charge quantities of the imaginary point charges are calculated, which can be regarded as a process of weighting the extracted components according to the position information. According to the dynamic sensitivity model of the sensor units as Equation (14), the i-th imaginary point charge is expressed as q_i = c_i/s_i, where s_i is the dynamic sensitivity of the No. 2 sensor unit at the position of q_i. By recording:

w_{i} = 1 {/ s}_{i}

(17)

one will have q_i = w_ic_i, where w_i is just the weight coefficient of c_i.

Finally, the monitoring parameter M is obtained by summing up the imaginary point charges:

M = \sum_{i} q_{i} = W \cdot C

(18)

where W is the vector of the weight coefficients and C is the vector of the extracted signal components. When q is in some symmetric lines of the HSESCA, some components could be zero, but it does not affect the computing process.

W is determined in the pre-step (see Figure 8). Firstly, according to the estimated positions of the imaginary point charges, an initial value of each weight coefficient is set using Equation (17). Then, the initial values are optimized by numerical simulations to make M and q as close as possible in the whole observation cross-section. In this study, W was finally determined as {238, 232, 208, 122, 119, 98, 31, 0.9} when K = −1. As for a certain HSESCA, the pre-step needs to be done only once, then the weight coefficients are regarded as a priori values in calculating M.

According to the definition of the dynamic sensitivity of ESA, the dynamic sensitivity in accordance with the component extraction-based array signal processing algorithm is modeled as:

S_{DA} (y, z) = | \frac{M}{q} | = \frac{W \cdot C}{q}

(19)

This reflects the monitoring accuracy of the HSESCA at different positions in the observation cross-section after adopting the algorithm.

It is seen from the steps that the accurate position of the monitored particle was not obtained. However, the position information was made use of by using the component extraction-based method. In addition, even accurate positions of the imaginary point charges were not necessary. As a result, the proposed array signal processing algorithm only contains some simple steps, making it promising in applications that need fast processing of the array signals.

5. Numerical Simulation and Experiment

5.1. Simulated Results and Discussion

5.1.1. Simulated Dynamic Sensitivity of the Sensor Units

As mentioned above, an HSESCA with eight sensor units is installed on a grounded pipeline. The pipeline has an inner radius of 200 mm; the hemispherical probes have a radius of 12 mm with their button faces 2 mm far from the inner wall of the pipeline. For the sake of simplicity, the amplification coefficient of the signal conditioner is set to K = −1. Under these conditions, the dynamic sensitivity of the sensor units is calculated by using Equation (14) at the grid points shown in Figure 10. It is seen that the grid points are representative, as they cover almost the whole observation cross-section.

The simulated results are shown by the surface diagrams in Figure 11. It is seen that the sensor units have the same dynamic sensitivity because they are identical in geometry, except that the diagrams rotate 45° in turn, according to the positions of the units. In addition, as for any sensor unit, its dynamic sensitivity reaches the maximum value near the probe surface, and it decreases sharply as the distance to the probe increases. To make a quantitative view, the contour diagram of the No. 7 sensor unit is provided as Figure 12. It shows that the dynamic sensitivity is greater than 0.9 near the probe surface. However, it decreases to below 0.1 just about 30 mm away and decreases to below 0.001 about 150 mm away. In a word, the dynamic sensitivity of the sensor units is inhomogeneous and quite localized, which means that a hemispherical sensor unit can only effectively detect charged particles near its probe. Therefore, it is of practical value to overcome this drawback by using ESAs, thus improving the monitoring accuracy.

5.1.2. Simulated Static Sensitivity of the HSESCA

The static sensitivity of the HSESCA in the observation cross-section is calculated by using Equation (11) and the calibration method. The results are shown in Figure 13.

It is seen from Figure 13 that, in the observation cross-section, the static sensitivity of the HSESCA distributes symmetrically in accordance with the distribution of the probes. In addition, just like the dynamic sensitivity of a certain sensor unit, the static sensitivity of the HSESCA reaches the maximum value greater than 0.9 near the probes, and it decreases sharply as distance to the probes increases. In most of the central area, the values of static sensitivity are below 0.05. That is to say, the static sensitivity of the HSESCA is inhomogeneous and quite localized, which means that an HSESCA will only effectively detect charged particles near the probes if just using the summation of the induced charges to reflect the charge quantity of a monitored particle. Definitely, the static sensitivity in the center region can be improved by increasing the size and number of probes. However, this will cause problems in installation and increase the disturbance to the flow. Therefore, it is advisable to ensure the monitoring accuracy in the center region by using appropriate array signal processing algorithms.

5.1.3. Simulated Dynamic Sensitivity of the HSESCA

For the sake of simplicity, here, we set the amplification coefficient K = −1. Then, by placing a point charge q = 1 C in turn at each grid point shown in Figure 10, the corresponding signal peaks of all of the sensor units are calculated using Equation (14). Next, they are processed by the component extraction-based array signal processing algorithm (see Section 4) to get the monitoring parameter value corresponding to each grid point. In this way, the values of dynamic sensitivity at all of the grid points are obtained. The results are shown in Figure 14.

It is seen from Figure 14 that the dynamic sensitivity of the HSESCA distributes symmetrically in accordance with the distribution of the probes. In addition, it is obvious that the dynamic sensitivity has relatively high values, not only near the probes, but also in most of the area of the observation cross-section. By making a comparison between Figure 13 and Figure 14, it is shown that the dynamic sensitivity and the static sensitivity have close values near the probes and the inner wall of the pipeline. However, the mean value of the static sensitivity is only about 0.04 in most of the central area, while that of the dynamic sensitivity reaches about 0.95. The sensitivity represents charge measurement accuracy, that is to say the closer to one the sensitivity is, the higher the accuracy is. Therefore, in the case of intermittent particle monitoring, the monitoring ability gets enhanced by 23 times in most of the central area because of using the component extraction-based array signal processing algorithm, which implies significant improvements in sensitivity homogeneity and monitoring accuracy of the HSESCA. Furthermore, it shows that there still distribute relatively low values near the inner wall of the pipeline, except where the probes are located. This is unavoidably caused by the influence of the pipeline on the electrostatic field. However, using a strategy to optimize the number of sensor units combined with optimizations on the proposed algorithm is likely to alleviate this drawback, which deserves further studies.

5.2. Experimental Results and Discussion

5.2.1. Experiment Apparatus

Experimental results were obtained by using an eight-channel HSESCA experiment apparatus. As shown in Figure 15, the experiment apparatus is mainly composed of an electrostatic generator, a particle release device, a pipeline, an HSESCA, a Faraday cylinder and a computer. The electrostatic generator produces ionized air in the particle release device by high voltage, which makes released particles charged. An HSESCA that has eight sensor units is installed in the midsection of the grounded pipeline. It produces array signals when a charged particle falls through the pipeline. Then, the signals are processed by using specially-designed software integrated in the computer. The Faraday cylinder connected with a charge meter is placed below the pipeline. A charged particle will finally fall into the Faraday cylinder; thus, the charge quantity of the particle is read from the charge meter with a measurement resolution of 1 pC. Steel balls with a diameter of 2 mm are used as the experimental particles. To create a consistent condition with the numerical simulation, the pipeline has an inner radius of 200 mm; the hemispherical probes have a radius of 12 mm with their button faces 2 mm far from the inner wall of the pipeline.

In the experiment, the particles were released at different test points, as shown in Figure 16. The points were set along two typical lines in the observation cross-section. Line 1 is the axis of the No. 1 sensor unit, and Line 2 is the bisector of the No. 1 and the No. 8 sensor units. Each line contains 12 test points, which are uniformly spaced by 15 mm. In order to reduce the position error of each release and the measurement error of the charge quantity, 10 tests were done at each point. Then, the mean valves of the measured signal peaks and those of the charge quantities were further used to calculate the static sensitivity and the dynamic sensitivity of the HSESCA.

5.2.2. Array Signals of the HSESCA

Five sets of typical array signals with a sample frequency of 4096 Hz are shown in Figure 17. They were acquired when particles were released at P₁, P₂, P₃, P₄ and P₅, respectively (see Figure 16). Among them, P₁ is the origin; P₂ and P₄ are 60 mm away from the origin; P₃ and P₅ are 120 mm away from the origin. Besides, the charge quantities of the particles were measured to be −32 pC, −20 pC, −46 pC, −25 pC and −29 pC, respectively.

Figure 17 shows that the measured signals are similar in shape as those obtained by numerical simulations (see Figure 5), except the sign symbol. The measured ones are minus because of the minus charge quantities of the particles. Besides, a signal variation trend versus the positions in Line 1 can be observed from Figure 17a–c. In detail, the signals of all of the sensor units are almost coincident at P₁, which is explained by the same distance from P₁ to every probe. However, with the release position moving from P₁ to P₃, the signal peak of the No. 1 sensor unit gets larger and larger. Conversely, other signal peaks decrease with respect to that of the No. 1 sensor unit; among them, that of the No. 5 sensor unit decreases the most. In addition, two symmetric sensor units about Line 1 have the same decrease scale in their signal peaks. A similar trend can also be observed when the release position moves from P₁ to P₅ along Line 2. This is because a charged particle induces greater charge on closer probes, and vice versa. In a word, the measured array signals have validated the theoretical ones in form. Moreover, it is validated that different combinations of peaks of the array signals imply different position information of the particles. This agrees with the theoretical analysis and is the precondition of the proposed array signal processing algorithm.

Ten tests were done at the release positions of P₁, P₂, P₃, P₄ and P₅, respectively. The absolute charge quantities of the particles and the absolute signal peaks of the No. 1 sensor unit were measured. Their relationships are shown in Figure 18.

This shows an obvious proportional relationship between the absolute charge quantities and the absolute peak voltages at each position. The different slopes indicate different values of dynamic sensitivity with respect to different positions, and greater slopes represent greater values of dynamic sensitivity. The errors are acceptable and can be explained by the position error of each release and the measurement error of the charge quantity. As a result, the proportional characteristic of the signal conditioner is validated.

5.2.3. Static Sensitivity Test of the HSESCA

To make a comparison between the theoretical and experimental results of the static sensitivity of the HSESCA, the theoretical values in Line 1 and Line 2 were firstly calculated by using Equation (11) and the calibration method when K = −1. Then, according to the definition of static sensitivity as Equation (10), the experimental value at each test point was obtained by dividing the summation of all of the signal peaks by the charge quantity of the corresponding particle, followed by a scaling to suppose K = −1. The results are shown in Figure 19, where the curves represent the theoretical values and the points represent the experimental values.

It can be seen that the experimental values show fine consistencies with the theoretical ones. Further analysis shows that the mean absolute value of relative error from all of the test points in Line 1 is 2.2544%, and that in Line 2 is 1.3516%. The errors are acceptable and can be explained by the errors in controlling the particles’ falling positions and measuring their charge quantities. Therefore, the sensing model of the HSESCA, as well as that of the sensor units are validated.

5.2.4. Dynamic Sensitivity Test of the HSESCA

Theoretical and experimental results of the dynamic sensitivity of the HSESCA were also compared. When the amplification coefficient was set to K = −1, the theoretical values in Line 1 and Line 2 were firstly calculated (see Section 5.1.3). After that, the measured signal peaks were processed by the component extraction-based array signal processing algorithm to obtain the experimental values. In detail, the peaks were decomposed, weighted and summated to obtain the values of the monitoring parameter. Then, the corresponding values of dynamic sensitivity are calculated by using Equation (19), followed by a scaling to suppose K = −1. The results are shown in Figure 20, where the curves represent the theoretical values and the points represent the experimental values.

It is observed that the experimental values match well with the theoretical ones. Further analysis shows that the mean absolute value of relative error from all of the test points in Line 1 is 3.3233%, and that in Line 2 is 0.5936%. Just as the results in the static sensitivity test, the errors here are acceptable and can be explained by the errors in controlling the particles’ falling positions and measuring their charge quantities. In addition, by making a comparison between Figure 19 and Figure 20, it is seen that the values of dynamic sensitivity are larger than 0.8 in most of the section of Line 1 and Line 2, but those of the static sensitivity are almost below 0.1. This demonstrates that the monitoring accuracy of the HSESCA for intermittent particles has been improved significantly by using the component extraction-based signal processing algorithm, which agrees with the theoretical analysis above. Therefore, the proposed algorithm is proven to be effective in actual measurements. More details can be found in Section 5.1.2 and Section 5.1.3.

6. Conclusions

By taking the processing of the array signals into consideration, the dynamic sensitivity of ESA has been defined to describe the sensitivity characteristics of ESAs in a systemic perspective. It builds a direct relationship between the monitored particles and the monitoring parameters of a whole ESA, thus reflecting the monitoring accuracy directly. An HSESCA along with its application in intermittent particle monitoring has been taken as an example. Its dynamic sensitivity in accordance with a proposed component extraction-based array signal processing algorithm has been analyzed quantitatively. Relevant numerical simulations have been made, and experimental validations have been carried out on an eight-channel HSESCA experiment apparatus. Detailed results have been provided, compared and discussed, which are summarized as follows:

The experimental results match well with the theoretical ones, which has validated the accuracy of the theoretical models and effectiveness of the corresponding methods.
Compared with static sensitivity, the dynamic sensitivity of the HSESCA has much greater values in most of the observation cross-section. This demonstrates that the component extraction-based array signal processing algorithm is effective at overcoming the defect of inhomogeneous and localized static sensitivity, thus making a significant improvement on the monitoring accuracy for intermittent particles.
There still exist relatively less sensitive zones near the inner wall of the pipeline after adopting the proposed algorithm. This is caused by the influence of the pipeline on the electrostatic field. A strategy to optimize the number of sensor units combined with optimizations on the proposed algorithm is likely to alleviate the drawback, which deserves further studies.
Ill-conditioning and noise may impact the results of the proposed array signal processing algorithm. This is of great significance, to be discussed in further studies.

Acknowledgments

The authors are grateful to the National Basic Research Program of China (Grant No. 2015CB057400) for supporting this work.

Author Contributions

This work has been done in collaboration between all authors. Xin Tang did the theoretical analyses and prepared the manuscript. Zhong-Sheng Chen designed the experiment and analyzed the experimental results. Yue Li and Zheng Hu did the simulation work and compared the results with those of the experiment. Yong-Min Yang revised the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

EMS	Electrostatic monitoring system
ESA	Electrostatic sensor array
HSESCA	Hemisphere-shaped electrostatic sensors’ circular array
PHM	Prognostics and health management

References

Sun, J.Z.; Zuo, H.F.; Liu, P.P.; Wen, Z.H. Experimental study on engine gas-path component fault monitoring using exhaust gas electrostatic signal. Meas. Sci. Technol. 2013, 24, 125107. [Google Scholar] [CrossRef]
Rahmat, M.F.; Thuku, I.T.; Tajdari, T.; Jusoff, K.; Ghazali, M.R. Sensing and filtering characteristics of electrostatic sensors for pneumatically conveyed particles. Int. J. Phys. Sci. 2011, 6, 5091–5103. [Google Scholar]
Coombes, J.R.; Yan, Y. Experimental investigations into the flow characteristics of pneumatically conveyed biomass particles using an electrostatic sensor array. Fuel 2015, 151, 11–20. [Google Scholar] [CrossRef]
Hu, H.L.; Dong, J.; Zhang, J.; Cheng, Y.J.; Xu, T.M. Identification of gas/solid two-phase flow regimes using electrostatic sensors and neural-network techniques. Flow Meas. Instrum. 2011, 22, 482–487. [Google Scholar] [CrossRef]
Rahmat, M.F.; Isa, M.D.; Rahim, R.A.; Hussin, T.A. Electrodynamics sensor for the image reconstruction process in an electrical charge tomography system. Sensors 2009, 9, 10291–10308. [Google Scholar] [CrossRef] [PubMed]
Addabbo, T.; Fort, A.; Garbin, R.; Mugnaini, M.; Rocchi, S.; Vignoli, V. Theoretical characterization of a gas path debris detection monitoring system based on electrostatic sensors and charge amplifiers. Measurement 2015, 64, 138–146. [Google Scholar] [CrossRef]
Wen, Z.H.; Zuo, H.F.; Pecht, M.G. Electrostatic monitoring of gas path debris for aero-engines. IEEE Trans. Reliab. 2011, 60, 33–40. [Google Scholar] [CrossRef]
Lin, J.; Chen, Z.S.; Hu, Z.; Yang, Y.M.; Tang, X. Analytical and numerical investigations into hemisphere-shaped electrostatic sensors. Sensors 2014, 14, 14021–14037. [Google Scholar] [CrossRef] [PubMed]
Tang, X.; Chen, Z.; Li, Y.; Hu, Z.; Yang, Y. Theoretical analysis and fem simulations on dynamic sensitivity of hemisphere-shaped electrostatic sensors. Adv. Mech. Eng. 2016, 8. [Google Scholar] [CrossRef]
Krabicka, J.; Yan, Y. Finite-element modeling of electrostatic sensors for the flow measurement of particles in pneumatic pipelines. IEEE Trans. Instrum. Meas. 2009, 58, 2730–2736. [Google Scholar] [CrossRef]
Zhang, W.; Wang, C.; Yang, W.; Wang, C.-H. Application of electrical capacitance tomography in particulate process measurement—A review. Adv. Powder Technol. 2014, 25, 174–188. [Google Scholar] [CrossRef]
Zhou, B.; Zhang, J.; Xu, C.; Wang, S. Image reconstruction in electrostatic tomography using a priori knowledge from ect. Nucl. Eng. Des. 2011, 241, 1952–1958. [Google Scholar] [CrossRef]
Gao, H.; Xu, C.; Fu, F.; Wang, S. Effect of particle distribution on electrostatic tomography system. In Proceedings of the 7th International Symposium on Measurement Techniques for Multiphase Flows, Xi’an, China, 26–30 October 2012; pp. 269–276.
Iliya, T.T.; Mohd, F.R. Development of a new method of sensitivity matrix for image reconstruction in electric charge tomography system using finite element method. J. Teknol. 2014, 70, 103–112. [Google Scholar]
Li, Y.; Yang, W. Image reconstruction by nonlinear landweber iteration for complicated distributions. Meas. Sci. Technol. 2008, 19, 094014. [Google Scholar] [CrossRef]
Martínez Olmos, A.; Botella, G.; Castillo, E.; Morales, D.P.; Banqueri, J.; García, A. A reconstruction method for electrical capacitance tomography based on image fusion techniques. Digit. Signal Process. 2012, 22, 885–893. [Google Scholar] [CrossRef]
Tizhe Thuku, I.; Rahmat, M.F.A. Finite-element method modeling in 4 and 16 sensors electric-charge tomography systems for particles moving in pipeline. Flow Meas. Instrum. 2014, 38, 9–20. [Google Scholar] [CrossRef]
Xu, C.; Li, J.; Gao, H.; Wang, S. Investigations into sensing characteristics of electrostatic sensor arrays through computational modelling and practical experimentation. J. Electrost. 2012, 70, 60–71. [Google Scholar] [CrossRef]
Thuku, I.T.; Rahmat, M.F.; Wahab, N.A.; Tajdari, T.; Yusuf, A.A. 2-D finite-element modeling of electrostatic sensor for tomography system. Sens. Rev. 2013, 33, 104–113. [Google Scholar] [CrossRef]
Ghazali, M.R.; Ibrahim, W.I.; Rahmat, M.F. Electrostatic sensor in circular and rectangular shape electrode for solid dry pneumatic conveyor system. In Proceedings of the 2011 International Conference on Circuits, System and Simulation, Bangkok, Thailand, 28–29 May 2011; IACSIT Press: Singapore, 2011; Volume 7, pp. 66–72. [Google Scholar]
Xu, C.; Wang, S.; Tang, G.; Yang, D.; Zhou, B. Sensing characteristics of electrostatic inductive sensor for flow parameters measurement of pneumatically conveyed particles. J. Electrost. 2007, 65, 582–592. [Google Scholar] [CrossRef]
Yang, W.Q.; Conway, W.F. Measurement of sensitivity distributions of capacitance tomography sensors. Rev. Sci. Instrum. 1997, 69, 233–236. [Google Scholar] [CrossRef]
Guo, Z.; Shao, F.; Lv, D. Sensitivity matrix construction for electrical capacitance tomography based on the difference model. Flow Meas. Instrum. 2009, 20, 95–102. [Google Scholar] [CrossRef]
Marashdeh, Q.; Teixeira, F.L. Sensitivity matrix calculation for fast 3-D electrical capacitance tomography (ECT) of flow systems. IEEE Trans. Magnet. 2010, 40, 1204–1207. [Google Scholar] [CrossRef]
Wang, A.; Marashdeh, Q.M.; Teixeira, F.L.; Fan, L.S. Electrical capacitance volume tomography: A comparison between 12- and 24-channels sensor systems. Prog. Electromagn. Res. M 2015, 41, 73–84. [Google Scholar] [CrossRef]
Zeeshan, Z.; Teixeira, F.L.; Marashdeh, Q. Sensitivity map computation in adaptive electrical capacitance volume tomography with multielectrode excitations. Electron. Lett. 2015, 51, 334–336. [Google Scholar] [CrossRef]
Fisher, C.E. Gas path debris monitoring—A 21st century PHM tool. In Proceedings of the 2000 IEEE Aerospace Conference Proceedings, Big Sky, MT, USA, 25 March 2000; Volume 6, pp. 441–448.
Powrie, H.; Novis, A. Gas path debris monitoring for F-35 joint strike fighter propulsion system PHM. In Proceedings of the 2006 IEEE Aerospace Conference, Big Sky, MT, USA, 4–11 March 2006; pp. 4012–4019.
Wen, Z.; Ma, X.; Zuo, H. Characteristics analysis and experiment verification of electrostatic sensor for aero-engine exhaust gas monitoring. Measurement 2014, 47, 633–644. [Google Scholar] [CrossRef]
Addabbo, T.; Fort, A.; Mugnaini, M.; Panzardi, E.; Vignoli, V. Multisensor electrostatic detection of moving charged particles. In Proceedings of the Sensors Applications Symposium (SAS 2016), Catania, Italy, 20–22 April 2016; pp. 445–450.
Addabbo, T.; Fort, A.; Mugnaini, M.; Panzardi, E.; Vignoli, V. A smart measurement system with improved low-frequency response to detect moving charged debris. IEEE Trans. Instrum. Meas. 2016, 65, 1874–1883. [Google Scholar] [CrossRef]
Xu, C.; Zhang, J.; Yang, D.; Zhou, B.; Tang, G.; Wang, S. Dense-phase pneumatically conveyed coal particle velocity measurement using electrostatic probes and cross correlation algorithm. J. Phys. 2009, 147, 012004. [Google Scholar] [CrossRef]
Zhang, J.Q.; Yan, Y. On-line continuous measurement of particle size using electrostatic sensors. Powder Technol. 2003, 135–136, 164–168. [Google Scholar] [CrossRef]
Zhou, B.; Yang, D.; Tang, G.; Xu, C.; Wang, S. Analysis on dynamic sensitivity of ring-shaped electrostatic sensor. Proc. CSEE 2008, 28, 44–49. [Google Scholar]

Figure 1. Schematic of a grounded pipeline-installed HSESCA.

Figure 2. Schematic of a proportional two-stage signal conditioner channel.

Figure 3. Schematic of a grounded pipeline-installed sensor unit of the HSESCA.

Figure 4. Schematic of b_i and β_i in the new Descartes coordinate system.

Figure 5. (a–h) Theoretical array signals of the HSESCA from No. 1–No. 8 sensor units.

Figure 6. Signal flow diagram of the i-th sensor unit of the HSESCA.

Figure 7. Signal flow diagram of the HSESCA.

Figure 8. Flow diagram of the component extraction-based array signal processing algorithm.

Figure 9. The original charged particle q and its imaginary point charges.

Figure 10. Grid points in the observation cross-section.

Figure 11. (a–h) Surface diagrams of dynamic sensitivity of the sensor units from No. 1–No. 8.

Figure 12. Contour diagram of dynamic sensitivity of the No. 7 sensor unit.

Figure 13. (a) Surface and (b) contour diagrams of the static sensitivity of the HSESCA.

Figure 14. (a) Surface and (b) contour diagrams of the dynamic sensitivity of the HSESCA.

Figure 15. Structure of the HSESCA experiment apparatus.

Figure 16. Test points set in the experiment.

Figure 17. Array signals of the HSESCA. (a) Released at P₁ with −32 pC charge; (b) released at P₂ with −20 pC charge; (c) released at P₃ with −46 pC charge; (d) released at P₄ with −25 pC charge; (e) released at P₅ with −29 pC charge.

Figure 18. Relationships between the absolute charge quantities and the absolute peak voltages of the No. 1 sensor unit.

Figure 19. Experimental and theoretical values of static sensitivity (a) in Line 1 and (b) in Line 2.

Figure 20. Experimental and theoretical values of dynamic sensitivity (a) in Line 1 and (b) in Line 2.

© 2016 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC-BY) license (http://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Tang, X.; Chen, Z.-S.; Li, Y.; Hu, Z.; Yang, Y.-M. Analysis of the Dynamic Sensitivity of Hemisphere-Shaped Electrostatic Sensors’ Circular Array for Charged Particle Monitoring. Sensors 2016, 16, 1403. https://doi.org/10.3390/s16091403

AMA Style

Tang X, Chen Z-S, Li Y, Hu Z, Yang Y-M. Analysis of the Dynamic Sensitivity of Hemisphere-Shaped Electrostatic Sensors’ Circular Array for Charged Particle Monitoring. Sensors. 2016; 16(9):1403. https://doi.org/10.3390/s16091403

Chicago/Turabian Style

Tang, Xin, Zhong-Sheng Chen, Yue Li, Zheng Hu, and Yong-Min Yang. 2016. "Analysis of the Dynamic Sensitivity of Hemisphere-Shaped Electrostatic Sensors’ Circular Array for Charged Particle Monitoring" Sensors 16, no. 9: 1403. https://doi.org/10.3390/s16091403

APA Style

Tang, X., Chen, Z. -S., Li, Y., Hu, Z., & Yang, Y. -M. (2016). Analysis of the Dynamic Sensitivity of Hemisphere-Shaped Electrostatic Sensors’ Circular Array for Charged Particle Monitoring. Sensors, 16(9), 1403. https://doi.org/10.3390/s16091403

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

Article Menu

Analysis of the Dynamic Sensitivity of Hemisphere-Shaped Electrostatic Sensors’ Circular Array for Charged Particle Monitoring

Abstract

1. Introduction

2. Sensing Model of an HSESCA

2.1. Basic Structure of an HSESCA

2.2. Theoretical Modeling of a Hemisphere-Shaped Electrostatic Sensor Unit

2.3. Theoretical Modeling of an HSESCA

3. Definition of the Dynamic Sensitivity of ESA

3.1. Definition of the Static Sensitivity of ESA

3.2. Theoretical Array Signals of an HSESCA

3.3. Definition and Interpretation of the Dynamic Sensitivity of ESA

4. A Component Extraction-Based Array Signal Processing Algorithm for Intermittent Particles

5. Numerical Simulation and Experiment

5.1. Simulated Results and Discussion

5.1.1. Simulated Dynamic Sensitivity of the Sensor Units

5.1.2. Simulated Static Sensitivity of the HSESCA

5.1.3. Simulated Dynamic Sensitivity of the HSESCA

5.2. Experimental Results and Discussion

5.2.1. Experiment Apparatus

5.2.2. Array Signals of the HSESCA

5.2.3. Static Sensitivity Test of the HSESCA

5.2.4. Dynamic Sensitivity Test of the HSESCA

6. Conclusions

Acknowledgments

Author Contributions

Conflicts of Interest

Abbreviations

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI