Investigation of Induced Charge Mechanism on a Rod Electrode

: Rod electrodes based on an electrostatic induction mechanism are widely used in various industrial applications, but the analytic solution of an induced charge mechanism on a metal rod electrode has not yet been systematically established. In this paper, the theoretical model of the induced charge on a rod electrode is obtained through the method of images. Then, the properties of the rod electrode under the action of the point charge are studied, including the induced charge density distribution on the rod electrode, the amount of the induced charge with di ﬀ erent diameters and lengths of the electrode, and the e ﬀ ective space region induced by the electrode. On this basis, a theoretical model of the induced current on a rod electrode is established, which is used to study the induced current properties by a moving point charge. It is found that both the magnitude and bandwidth of the induced current increase with the increased point charge velocity. Finally, three experimental studies are conducted, and the experimental results show good consistency with the analysis of the theoretical model, verifying the correctness, and accuracy of the model. In addition, the induced charge mechanism studied in this paper can act as an e ﬀ ective basis for the rod electrode sensor design in terms of the optimal radius and length.


Introduction
Dust pollution is a common issue in industrial and mining enterprises, such as coal mining, iron mining, etc. [1,2]. Dust moves together with the ventilation air and settles on walls and equipment. Dust in the ventilation air can have negative effects on the working conditions, posing a risk to the health of workers. Many measures have been taken to improve the air quality [3,4]. Real-time and accurate detection of the dust concentration is a basic guarantee to ensure an effective dust removal [5,6].
Based on the electrostatic induction phenomenon, electrodes are widely used to measure the electric charges carried on solid particles, mass flow rate, concentration, volume loading, mean flow velocity, and other electrical and mechanical parameters in two-phase and multiphase gas-solid flows such as those in pneumatic conveyances, in the air, etc. [7]. When charged particles move close to or away from a metal electrode, a charge with the opposite polarity to the dust particles is induced on the surface of the electrode. Hence, the physical parameters of the charged particles can be acquired by measuring the characteristics of the induced charge on the electrode. In practical applications, metal electrodes are fabricated in different shapes to acquire different measurement parameters in various environments. Typical electrodes mainly include the ring, curved, square, and rod electrodes and arrays composed of electrodes of different shapes.
The ring electrode has received wide attention by researchers because of its noninvasive characteristics. It is an electrode form with relatively mature theory and applications. Weinheimer [8] derived a charge numerical solution on the surface of ring electrodes induced by the point charge, which was applied to the measurement of meteorological precipitation charge. Yan, Gajewski, and Woodhead used a correlation method to measure velocity after studying the sensing mechanism, spatial sensitivity, and spatial filtering effect of noninvasive ring induction electrodes [9][10][11]. In recent years, ring charge-sensing sensors have been widely utilized in the measurement of dilute phase/dense phase of gas-solid two-phase flow parameters [12][13][14]. At the same time, some modern signal processing algorithms have further improved the measurement accuracy [15][16][17][18][19]. For example, Wang et al. [20] improved the measurement accuracy by applying the wavelet transform to the multiphase flow parameters. Considering that the signal measured by the ring electrode is an average feature over the entire cross-section, Zhang, Yang, and Dong's research found that the arc-shaped sensing electrode had an advantage in particle velocity and concentration distribution measurement in the monitoring of particle motion in gas-solid fluidized beds [21][22][23]. Qian also combined arc-shaped electrodes with digital images for the measurement of biomass-coal particles in fuel-injection pipelines [24]. Liu and Yao obtained the characteristics of square electrodes through theoretical and simulation studies, and utilized them on a square pneumatic conveying pipeline [25,26]. Compared with the above two models, Zhang used the method of images to obtain a simple analytical solution of the square electrode [27]. With the increasing complexity of the measurement environment, combinations of multistatic sensors and even electrostatic sensor arrays are applied to acquire different parameters of various pipelines [28,29].
Rod sensors based on the electrostatic induction mechanism have been widely used in industrial fields [7,30], mainly due to their simpler installation compared to other types of sensors, its working principle, and its applications in the real world are shown as Figure 1. Since they can accurately reflect pollutant emissions such as the steel mill flue gas, they have been widely applied for the detection of particulate matter concentration and dust in dust collector bags. Although the noninvasiveness is a great advantage of the ring electrodes, it usually takes the form of a spool piece installed in line with the pipe, which leads to an expensive and challenging installation. Moreover, signals collected by the ring electrodes are an overall result of the induced signals, making it difficult to detect local flow regimens. In comparison to the ring electrodes, the installation of rod electrodes is easier, since they only need a suitable drilling hole at any position of the pipe, making it possible to detect local flow regimes [30]. Therefore, rod electrodes are widely used in industrial areas and some research work has been performed. Shao compared the advantages and disadvantages of the ring electrode and rod electrode electrostatic sensors in measuring the pulverized coal speed [31,32], and proved that both types of electrodes achieved the same measurement accuracy in the coal dust measurement.
The ring electrode has received wide attention by researchers because of its noninvasive characteristics. It is an electrode form with relatively mature theory and applications. Weinheimer [8] derived a charge numerical solution on the surface of ring electrodes induced by the point charge, which was applied to the measurement of meteorological precipitation charge. Yan, Gajewski, and Woodhead used a correlation method to measure velocity after studying the sensing mechanism, spatial sensitivity, and spatial filtering effect of noninvasive ring induction electrodes [9][10][11]. In recent years, ring charge-sensing sensors have been widely utilized in the measurement of dilute phase/dense phase of gas-solid two-phase flow parameters [12][13][14]. At the same time, some modern signal processing algorithms have further improved the measurement accuracy [15][16][17][18][19]. For example, Wang et al. [20] improved the measurement accuracy by applying the wavelet transform to the multiphase flow parameters. Considering that the signal measured by the ring electrode is an average feature over the entire cross-section, Zhang, Yang, and Dong's research found that the arc-shaped sensing electrode had an advantage in particle velocity and concentration distribution measurement in the monitoring of particle motion in gas-solid fluidized beds [21][22][23]. Qian also combined arcshaped electrodes with digital images for the measurement of biomass-coal particles in fuel-injection pipelines [24]. Liu and Yao obtained the characteristics of square electrodes through theoretical and simulation studies, and utilized them on a square pneumatic conveying pipeline [25,26]. Compared with the above two models, Zhang used the method of images to obtain a simple analytical solution of the square electrode [27]. With the increasing complexity of the measurement environment, combinations of multistatic sensors and even electrostatic sensor arrays are applied to acquire different parameters of various pipelines [28,29].
Rod sensors based on the electrostatic induction mechanism have been widely used in industrial fields [7,30], mainly due to their simpler installation compared to other types of sensors, its working principle, and its applications in the real world are shown as Figure 1. Since they can accurately reflect pollutant emissions such as the steel mill flue gas, they have been widely applied for the detection of particulate matter concentration and dust in dust collector bags. Although the noninvasiveness is a great advantage of the ring electrodes, it usually takes the form of a spool piece installed in line with the pipe, which leads to an expensive and challenging installation. Moreover, signals collected by the ring electrodes are an overall result of the induced signals, making it difficult to detect local flow regimens. In comparison to the ring electrodes, the installation of rod electrodes is easier, since they only need a suitable drilling hole at any position of the pipe, making it possible to detect local flow regimes [30]. Therefore, rod electrodes are widely used in industrial areas and some research work has been performed. Shao compared the advantages and disadvantages of the ring electrode and rod electrode electrostatic sensors in measuring the pulverized coal speed [31,32], and proved that both types of electrodes achieved the same measurement accuracy in the coal dust measurement. Electrostatic induction is a basic physical phenomenon in which the opposite induced charge is induced on the conductor if the point charge is close to the conductor. It is a qualitative conclusion. However, the analytic solution of the amount of induced charge when the point charge is close to the conductor has always been a difficult problem. It is a physical boundary value problem, which is difficult to represent with an analytical solution. The common method is to solve the Poisson equation through numerical calculations to get the amount of induced charge on the electrode. Krabicka studied the characteristics of electrostatic charges on rod electrodes by the finite element analysis (FEA) method and obtained an approximate solution [30]. However, this method requires remodeling of rod sensing electrodes of different lengths and diameters, which increases the modeling time. The simulation of large sensing electrodes takes too much time and the simulation accuracy is heavily dependent on the simulation software such as COMSOL. Therefore, this method is limited in practical application. Chen [33] established a mathematical model of rod electrodes by a theoretical derivation, but the fact that the induction conductor is a metal conductor was neglected, and the metal conductor was modeled as an insulator.
Analytic solutions of rod electrodes can be used to optimize the sensor design or interpret how particles at various locations influence the signal; for example, the bandwidth and amplitude. In order to establish a simple and easy-to-use mathematical model of rod electrodes, the method of images and the symmetry of rod electrodes are employed in this paper. The mathematical formula of the amount of charge induced by the point charge on the sensing rod electrode is obtained, and the equation is then used to study the physical properties of the sensing electrode under the action of the point charge. Based on this model, the distribution characteristics of the induced charge on the surface of the electrode and the influence of the electrode length on the induced charge density are studied, the induced charge of the induced electrode is simulated when the point charge moves in different directions, the amount of the induced current on the sensing electrode is studied, the spectral characteristics and the influence of the general measurement model of the induced current are analyzed, the experimental model is established, and the validity and accuracy of the model are verified by experiments.

Induced Charge Model on Rod Electrode by Point Charge
The method of images is an indirect method to solve the electrostatic field problem by applying the uniqueness theorem. The electrostatic method can be used to treat the actual partitioned uniform medium as uniform, and replace the actual complex charge distribution on the boundary with a simple charge distribution on the virtual closed setting boundary of the research field for calculation. According to the uniqueness theorem, this result is correct as long as the electric field generated by the imaginary charge together with the actual charge within the boundary satisfies a given boundary condition [34].
In this paper, the relationship between the rod electrode and point charge is established by the method of images, which mainly models the relationship of the induced charge by the point charge with two basic parameters of the rod electrode, length, and radius. The steps of the Electrostatic induction is a basic physical phenomenon in which the opposite induced charge is induced on the conductor if the point charge is close to the conductor. It is a qualitative conclusion. However, the analytic solution of the amount of induced charge when the point charge is close to the conductor has always been a difficult problem. It is a physical boundary value problem, which is difficult to represent with an analytical solution. The common method is to solve the Poisson equation through numerical calculations to get the amount of induced charge on the electrode. Krabicka studied the characteristics of electrostatic charges on rod electrodes by the finite element analysis (FEA) method and obtained an approximate solution [30]. However, this method requires remodeling of rod sensing electrodes of different lengths and diameters, which increases the modeling time. The simulation of large sensing electrodes takes too much time and the simulation accuracy is heavily dependent on the simulation software such as COMSOL. Therefore, this method is limited in practical application. Chen [33] established a mathematical model of rod electrodes by a theoretical derivation, but the fact that the induction conductor is a metal conductor was neglected, and the metal conductor was modeled as an insulator.
Analytic solutions of rod electrodes can be used to optimize the sensor design or interpret how particles at various locations influence the signal; for example, the bandwidth and amplitude. In order to establish a simple and easy-to-use mathematical model of rod electrodes, the method of images and the symmetry of rod electrodes are employed in this paper. The mathematical formula of the amount of charge induced by the point charge on the sensing rod electrode is obtained, and the equation is then used to study the physical properties of the sensing electrode under the action of the point charge. Based on this model, the distribution characteristics of the induced charge on the surface of the electrode and the influence of the electrode length on the induced charge density are studied, the induced charge of the induced electrode is simulated when the point charge moves in different directions, the amount of the induced current on the sensing electrode is studied, the spectral characteristics and the influence of the general measurement model of the induced current are analyzed, the experimental model is established, and the validity and accuracy of the model are verified by experiments.

Induced Charge Model on Rod Electrode by Point Charge
The method of images is an indirect method to solve the electrostatic field problem by applying the uniqueness theorem. The electrostatic method can be used to treat the actual partitioned uniform medium as uniform, and replace the actual complex charge distribution on the boundary with a simple charge distribution on the virtual closed setting boundary of the research field for calculation. According to the uniqueness theorem, this result is correct as long as the electric field generated by the imaginary charge together with the actual charge within the boundary satisfies a given boundary condition [34].
In this paper, the relationship between the rod electrode and point charge q is established by the method of images, which mainly models the relationship of the induced charge by the point charge q with two basic parameters of the rod electrode, length, and radius. The steps of the proposed method to model the induced charge on a rod electrode is shown in Figure 2. The following subsections explains the steps in detail step by step. proposed method to model the induced charge on a rod electrode is shown in Figure 2. The following subsections explains the steps in detail step by step.
Step 1 Select a random point on the electrode surface Step 2 Apply the Method of images Step 3 Verify the suppose satisfies the uniqueness theorem Step 4 Calculate the combined field strength Step 5 Calculate the total amount of charge

Step 1 Select A Random Point On the Electrode Surface
To model the induced current on a rod electrode by a point charge, a cylindrical coordinate system is established, as shown Figure 3a, with the origin point being the center of the electrode. Point C is the intersection of Oq and the electrode's surface, point is one point on the surface of the electrode, and plane OAC is perpendicular to line L. Suppose is the angle between line and line , the coordinate of point can be written as ( , , 0). Point charge is in plane and is located on the extension line of , and the distance to is ℎ and the distance to point is ℎ 1 . Point is a random point on the electrode surface and is the vertical distance from B to plane , hence the coordinates of point are ( , , ).

Step 2 Apply the Method of Images
In addition to the electric field caused by point charge q in the electrolyte, the effect of the induced charge on the rod electrode should also be considered. However, the electric field distribution on the sensing electrode is unknown, which is the problem that this paper needs to solve. As shown in Figure 3b, it is assumed that an infinitely large plane is tangent to the metal rod electrode and the intersection between the rod electrode and plane forms a straight line . The electric field intensity and the electric field line distribution generated on infinite plane can be obtained by the method of images, as shown in Figure 3c. Note that these electric field lines are perpendicular to the surface of the rod electrode. It is assumed that electric field E on line L of infinite plane S formed by point charge q is electric field E on the line of the rod electrode. Then, the image method is assumed to have a point charge at q' with an equivalent charge (-q) to point charge q.

Step 3 Verify the Suppose Satisfies the Uniqueness Theorem
After removing the electrode and conductive plate S, and the following conditions are still satisfied: (1) In addition to the position where the point charge is located, 2 = 0 is satisfied everywhere, and is the potential.
(2) Taking infinity as the reference point, the potential at the interface between the medium and conductor L is zero.
(3) The direction of the electric field line is unchanged, which is the direction of the vertical electrode surface pointing to the axis.
These conditions satisfy the uniqueness theorem.

Step 4 Calculate the Combined Field Strength
By applying the method of images, we can obtain the electric field intensity of point . A cross-section view of plane is illustrated in Figure 3d; the projection of point is point and the projection of plane is line ′ . The connecting line between point charge and image charge ′ intersects with ′ at point , the distance between point and point is , and the distance between point and point is ℎ ′ ; is the distance between point charge and B.

Step 1 Select A Random Point On the Electrode Surface
To model the induced current on a rod electrode by a point charge, a cylindrical coordinate system is established, as shown Figure 3a, with the origin point O being the center of the electrode. Point C is the intersection of Oq and the electrode's surface, point A is one point on the surface of the electrode, and plane OAC is perpendicular to line L. Suppose θ is the angle between line OA and line OC, the coordinate of point A can be written as (r, θ, 0). Point charge q is in plane OAC and is located on the extension line of OC, and the distance to C is h and the distance to point A is h 1 . Point B is a random point on the electrode surface and x is the vertical distance from B to plane OAC, hence the coordinates of point B are (r, θ, x).

Step 2 Apply the Method of Images
In addition to the electric field caused by point charge q in the electrolyte, the effect of the induced charge on the rod electrode should also be considered. However, the electric field distribution on the sensing electrode is unknown, which is the problem that this paper needs to solve. As shown in Figure 3b, it is assumed that an infinitely large plane S is tangent to the metal rod electrode and the intersection between the rod electrode and plane S forms a straight line L. The electric field intensity E and the electric field line distribution generated on infinite plane S can be obtained by the method of images, as shown in Figure 3c. Note that these electric field lines are perpendicular to the surface of the rod electrode. It is assumed that electric field E on line L of infinite plane S formed by point charge q is electric field E on the line of the rod electrode. Then, the image method is assumed to have a point charge at q with an equivalent charge (−q) to point charge q.

Step 3 Verify the Suppose Satisfies the Uniqueness Theorem
After removing the electrode and conductive plate S, and the following conditions are still satisfied: (1) In addition to the position where the point charge is located, ∇ 2 ϕ = 0 is satisfied everywhere, and ϕ is the potential. (2) Taking infinity as the reference point, the potential at the interface between the medium and conductor L is zero. (3) The direction of the electric field line is unchanged, which is the direction of the vertical electrode surface pointing to the axis.
These conditions satisfy the uniqueness theorem.

Step 4 Calculate the Combined Field Strength
By applying the method of images, we can obtain the electric field intensity E of point B. A cross-section view of plane OAC is illustrated in Figure 3d; the projection of point B is point A and the projection of plane S is line S . The connecting line between point charge q and image charge q intersects with S at point D, the distance between point A and point D is l, and the distance between point q and point D is h ; d is the distance between point charge q and B. The mathematical expressions of ℎ ′ , , and ℎ 1 and the range of values of can be obtained, as follows: = (ℎ 1 2 + 2 ) 1/2 . The mathematical expressions of h , l, and h 1 and the range of values of θ can be obtained, as follows: To determine the range of θ, suppose point charge q carries a positive charge and is located at a position with minimum distance h to the surface of the electrode, as shown in Figure 3f. F is at the intersection of the tangential line through point charge q, and the angle θ 2 is equal to acos(r/(r + h)).
For any point M at a position that satisfies the condition θ 1 < acos(r/(r + h)), its electric field strength can be analyzed as follows. According to the method of images, the imaginary charge is q , carrying a negative charge. E Mq is the electric field strength generated by point charge q, and E Mq is the electric field strength generated by point charge q . Their synthetic electric field of strength is E M , which points to the center point O. The orientation of E M conforms with the electric field on the metal being perpendicular to the metal surface. As is known, if the charge of q is positive, the induced charge on the metal surface should be negative. This further confirms the correctness of the orientation of E M .
For any point G at a position that satisfies the condition θ 3 > acos(r/(r + h)), its electric field strength can be analyzed similarly. According to the method of images, the imaginary charge is q , also carrying a negative charge. E Gq is the electric field strength generated by point charge q, and E Gq is the electric field strength generated by point charge q . The synthetic field of strength is E G , which points to a direction off the center O and is perpendicular to the metal surface. This indicates that the charge at point G is positive, which does not comply with the induced charge being negative. Hence, such a point should not have an induced charge generated from point charge q.
According to the above analysis, θ should not be larger than acos(r/(r + h)). Using the same method, we can get that the minimum value of θ is not smaller than −acos(r/(r + h)).
In summary, the range of θ can be determined as the following when the method of images is used: It can be seen from Figure 3e that field strength E qB is the summation result of point charge q and its image q . E qB can be expressed as and E can be expressed as E = E qB cosα.
Electric field strength E of point B based on the method of images can be obtained as where ε 0 is the vacuum permittivity.

Step 5 Calculate the Total Amount of Charge
The charge density ρ of point B is obtained as The dS marked in orange on the surface of the electrode contains point B, as in Figure 3c; the area of dS is expressed as When dθ and dx are infinite to zero, the charge density of dS will be the density of point B.
The total amount of the induced charge on the metal electrode is obtained: where x 1 , x 2 are the x-coordinates of the ends of the metal electrode. By changing the order of the integration, we get The result is obtained as Let the charge distribution function along with x 1 , x 2 and θ be F(x 1 , x 2 , θ): When point charge q moves infinitely away from the electrode, the amount of induced charge is lim h→∞ Q = 0. This means the amount of the induced charge generated by the point charge at infinity on the rod-shaped metal electrode is zero. From the analysis of Section 3.2.1, it can been seen that when point charge q is infinitely close to the electrode, the amount of the induced charge is lim h→0 Q = −q; that is, the amount of the induced charge generated on the rod-shaped metal electrode is −q when the point charge is infinitely close to the metal electrode.

Characteristics of Induced Charge on a Rod Electrode
In this section, the basic characteristics of the induced charge on a rod electrode are analyzed by utilizing the model obtained in Section 2.

Induced Charge Distribution under the Effect of Point Charge
The induced charge distribution with θ is studied. Four cases of h are studied, as shown in Figure 4a, where h is selected as 0.05 m, 0.10 m, 0.15 m, and 0.20 m. The diameter of the electrode is r = 0.005 m, the length of the rod electrode is 0.50 m, and the two ends of the electrode are x 1 = −0.25 m and x 2 = 0.25 m. Charge q carries a charge quantity of −1 C. Figure 4b shows the charge distribution with θ. It can be observed that the charge decreases as distance h increases and the charge reaches a maximum value when θ = 0.
where 1 , 2 are the x-coordinates of the ends of the metal electrode. By changing the order of the integration, we get The result is obtained as Let the charge distribution function along with 1 , 2 and be ( 1 , 2 , ): When point charge moves infinitely away from the electrode, the amount of induced charge is ℎ→∞ = 0. This means the amount of the induced charge generated by the point charge at infinity on the rod-shaped metal electrode is zero. From the analysis of section 3.2.1, it can been seen that when point charge is infinitely close to the electrode, the amount of the induced charge is ℎ→0 = − ; that is, the amount of the induced charge generated on the rod-shaped metal electrode is − when the point charge is infinitely close to the metal electrode.

Characteristics of Induced Charge on a Rod Electrode
In this section, the basic characteristics of the induced charge on a rod electrode are analyzed by utilizing the model obtained in Section 2.

Induced Charge Distribution under the Effect of Point Charge
The induced charge distribution with is studied. Four cases of ℎ are studied, as shown in Figure 4a, where ℎ is selected as 0.05 m, 0.10 m, 0.15 m, and 0.20 m. The diameter of the electrode is = 0.005 m, the length of the rod electrode is 0.50 m, and the two ends of the electrode are 1 = −0.25 m and 2 = 0.25 m. Charge q carries a charge quantity of −1 C. Figure 4b shows the charge distribution with . It can be observed that the charge decreases as distance ℎ increases and the charge reaches a maximum value when = 0. The variation of the induced charge along ( , ) on the electrode is shown in Figure 5. The induced charge reaches its maximum value at the point ( , 0, 0). If the point charge is closer to the electrode, the charge distribution is more concentrated. When the sensing range is smaller, the induced charge is greater. This indicates that the closer the point charge is to the rod electrode, the more charge is induced in the small sensing region, and the charge distribution tends to gradually decrease when the distance from the point ( , 0, 0) increases.

Effect of Distance between Electrode and Point Charge
As shown in Figure 6a, the point charge q carries -1 C charge and moves at a constant velocity of 5 / along the direction parallel to the pipeline and perpendicular to the electrode's axis. The length of the electrode is 0.1 and the radius of electrode is 0.005 m. Suppose the distance between the electrode and point charge is ℎ. We performed simulations on five distances, as presented in Figure 6b  The variation of the induced charge along (θ, x) on the electrode is shown in Figure 5. The induced charge reaches its maximum value at the point (r, 0, 0). If the point charge is closer to the electrode, the charge distribution is more concentrated. When the sensing range x is smaller, the induced charge is greater. This indicates that the closer the point charge is to the rod electrode, the more charge is induced in the small sensing region, and the charge distribution tends to gradually decrease when the distance from the point (r, 0, 0) increases. The variation of the induced charge along ( , ) on the electrode is shown in Figure 5. The induced charge reaches its maximum value at the point ( , 0, 0). If the point charge is closer to the electrode, the charge distribution is more concentrated. When the sensing range is smaller, the induced charge is greater. This indicates that the closer the point charge is to the rod electrode, the more charge is induced in the small sensing region, and the charge distribution tends to gradually decrease when the distance from the point ( , 0, 0) increases.

Effect of Distance between Electrode and Point Charge
As shown in Figure 6a, the point charge q carries -1 C charge and moves at a constant velocity of 5 / along the direction parallel to the pipeline and perpendicular to the electrode's axis. The length of the electrode is 0.1 and the radius of electrode is 0.005 m. Suppose the distance between the electrode and point charge is ℎ. We performed simulations on five distances, as presented in Figure 6b

Effect of Distance between Electrode and Point Charge
As shown in Figure 6a, the point charge q carries −1 C charge and moves at a constant velocity of 5 m/s along the direction parallel to the pipeline and perpendicular to the electrode's axis. The length of the electrode is 0.1 m and the radius of electrode r is 0.005 m.
Suppose the distance between the electrode and point charge is h. We performed simulations on five distances, as presented in Figure 6b It can be observed that the amount of charge decreases when the distance between the electrode and point charge increases. This applies when the distance between the electrode and point charge is within the range of −0.5 m to 0.5 m, and the amount of the induced charge can be neglected when it falls outside of this range Therefore, if there are charged particles uniformly distributed around the electrode, it can be considered that the induced charge on the electrode is generated by charges within a cylindrical section with a radius of 0.5 m.
induced on the electrode by these five simulations are shown in Figure 6c. Here, we define the distance from the point charge before it passes through points A, B, C, D, E as negative and the distance after as positive. It can be observed that the amount of charge decreases when the distance between the electrode and point charge increases. This applies when the distance between the electrode and point charge is within the range of −0.5 m to 0.5 m, and the amount of the induced charge can be neglected when it falls outside of this range Therefore, if there are charged particles uniformly distributed around the electrode, it can be considered that the induced charge on the electrode is generated by charges within a cylindrical section with a radius of 0.5 m. As shown in Figure 7a, the point charge moves from the surface of the electrode to point A and then to point B. The distance between point A and the surface of the electrode is 0.5 m and the distance between point B and the electrode is 0.5 m, the amount of charge is −1 C, the radius of rod electrode is 0.005 m, and the length of the electrode is 0.5 m. As shown in Figure 7b, the closer the point charge is to the electrode, the larger the induced charge. With increased distance from the point charge, the induced charge decreases rapidly. When the distance between the point charge and the electrode reaches a certain value, the change of the induced charge is very small. The trend in Figure 7b can also be observed in Figure 6c. It is also consistent with the conclusion that when the charged particles are evenly distributed around the electrode, the charge on the electrode is induced by a cylindrical region around the electrode with a radius of 0.5 m. As shown in Figure 7a, the point charge q moves from the surface of the electrode to point A and then to point B. The distance between point A and the surface of the electrode is 0.5 m and the distance between point B and the electrode is 0.5 m, the amount of charge Q is −1 C, the radius of rod electrode r is 0.005 m, and the length of the electrode is 0.5 m. As shown in Figure 7b, the closer the point charge is to the electrode, the larger the induced charge. With increased distance from the point charge, the induced charge decreases rapidly. When the distance between the point charge and the electrode reaches a certain value, the change of the induced charge is very small. The trend in Figure 7b can also be observed in Figure 6c. It is also consistent with the conclusion that when the charged particles are evenly distributed around the electrode, the charge on the electrode is induced by a cylindrical region around the electrode with a radius of 0.5 m. From the above analysis, it is concluded that the charged particles capable of inducing a charge on the rod electrode are mostly distributed in the shadow range, as shown in Figure 10. Note that a circular conveying pipe is considered in Figure 10 and the charged particles are uniformly distributed. From the above analysis, it is concluded that the charged particles capable of inducing a charge on the rod electrode are mostly distributed in the shadow range, as shown in Figure 10. Note that a circular conveying pipe is considered in Figure 10 and the charged particles are uniformly distributed. From the above analysis, it is concluded that the charged particles capable of inducing a charge on the rod electrode are mostly distributed in the shadow range, as shown in Figure 10. Note that a circular conveying pipe is considered in Figure 10 and the charged particles are uniformly distributed.

Effect of Electrode Length
The length of the rod electrode is an important parameter that needs to be determined in the sensor design stage. As shown in Figure 11a Figure 11b. It can be observed that as the length of the electrode increases, the amount of the induced charge increases, but the increase rate decreases with the increase of the electrode. It can be concluded that the longer the electrode, the more charge can be induced. However, when the electrode length increases to a certain extent, its growth does not cause a significant change in the amount of the induced charge.
As shown in Figure 12a, the length of the electrode varies from 0 m to 4 m and the distance from the surface of the electrode is ℎ = 0.1 m. The changed amount of the induced charge on the electrode is shown in Figure 12b. The charge of the point charge is = −1 C, and the radius of the rod electrode is = 0.005 m. As shown in Figure 13b, once the length of the electrode reaches 1 m or more, the increase rate of the induced change becomes very small. In other words, increasing the electrode length contributes very little to the amount of the induced charge when the length reaches 1 m.

Effect of Electrode Length
The length of the rod electrode is an important parameter that needs to be determined in the sensor design stage. As shown in Figure 11a Figure 11b. It can be observed that as the length of the electrode increases, the amount of the induced charge increases, but the increase rate decreases with the increase of the electrode. It can be concluded that the longer the electrode, the more charge can be induced. However, when the electrode length increases to a certain extent, its growth does not cause a significant change in the amount of the induced charge.
As shown in Figure 12a, the length of the electrode varies from 0 m to 4 m and the distance from the surface of the electrode is ℎ = 0.1 m. The changed amount of the induced charge on the electrode is shown in Figure 12b. The charge of the point charge is = −1 C, and the radius of the rod electrode is = 0.005 m. As shown in Figure 13b, once the length of the electrode reaches 1 m or more, the increase rate of the induced change becomes very small. In other words, increasing the electrode length contributes very little to the amount of the induced charge when the length reaches 1 m.

Effect of Electrode Length
The length of the rod electrode is an important parameter that needs to be determined in the sensor design stage. As shown in Figure 11a Figure 11b. It can be observed that as the length of the electrode increases, the amount of the induced charge increases, but the increase rate decreases with the increase of the electrode. It can be concluded that the longer the electrode, the more charge can be induced. However, when the electrode length increases to a certain extent, its growth does not cause a significant change in the amount of the induced charge.
As shown in Figure 12a, the length of the electrode varies from 0 m to 4 m and the distance from the surface of the electrode is h = 0.1 m. The changed amount of the induced charge on the electrode is shown in Figure 12b. The charge of the point charge is Q = −1 C, and the radius of the rod electrode is r = 0.005 m. As shown in Figure 13b, once the length of the electrode reaches 1 m or more, the increase rate of the induced change becomes very small. In other words, increasing the electrode length contributes very little to the amount of the induced charge when the length reaches 1 m.

Effect of Electrode Radius
As shown in Figure 13a, the distance of four points J, K, M, and N from the surface of the electrode is 0.05 m, 0.10 m, 0.15 m, and 0.20 m, respectively. The charge amount of the point charge is = −1 C, and the length of the rod electrode is 0.1 m. The radius varies from 0.001 m to 0.2 m. The amount of the induced charge on the electrode is shown in Figure 13b. As the radius of the electrode increases, the amount of the induced charge increases. This means that electrodes with a larger radius can induce more charge, which is beneficial to the design of the detection circuit in later stages. However, the radius of the electrode generally cannot be made too large to adapt to the actual situation in practice.

Induced Current on a Rod Electrode by Moving Point Charge
As shown in Figure 14, point charge q moves with velocity ( ), assuming that its position at time is ( 0 , 0 , 0 ), then its position at time is ( , , ).
ℎ in Equation (1) can be expressed as

Effect of Electrode Radius
As shown in Figure 13a, the distance of four points J, K, M, and N from the surface of the electrode is 0.05 m, 0.10 m, 0.15 m, and 0.20 m, respectively. The charge amount of the point charge is Q = −1 C, and the length of the rod electrode is 0.1 m. The radius varies from 0.001 m to 0.2 m. The amount of the induced charge on the electrode is shown in Figure 13b. As the radius of the electrode increases, the amount of the induced charge increases. This means that electrodes with a larger radius can induce more charge, which is beneficial to the design of the detection circuit in later stages. However, the radius of the electrode generally cannot be made too large to adapt to the actual situation in practice.

Induced Current on a Rod Electrode by Moving Point Charge
As shown in Figure 14, point charge q moves with velocity v(t), assuming that its position at time zero is (x 0 , y 0 , z 0 ), then its position at time t is (x t , y t , z t ).

Effect of Electrode Radius
As shown in Figure 13a Figure 13b. As the radius of the electrode increases, the amount of the induced charge increases. This means that electrodes with a larger radius can induce more charge, which is beneficial to the design of the detection circuit in later stages. However, the radius of the electrode generally cannot be made too large to adapt to the actual situation in practice.

Induced Current on a Rod Electrode by Moving Point Charge
As shown in Figure 14, point charge q moves with velocity ( ), assuming that its position at time is ( 0 , 0 , 0 ), then its position at time is ( , , ).
ℎ in Equation (1) can be expressed as x t , y t , and z t can be expressed as (15) h in Equation (1) can be expressed as Equation (13) can be expressed as The induced current can be calculated by

Simulation for Induced Current
As shown in Figure 14, suppose that point charge q moves from (0, y 0 , z 0 ) to (0, y 0 , −z 0 ) with a different constant velocity v of 2 m/s, 5 m/s, 10 m/s, and 20 m/s. The charge of the point charge is Q = −1 C, y 0 = 0.1 m, the length of the electrode is 0.1 m, and the radius of the electrode is r = 0.005 m. As shown in Figure 15a, the induced current increases with the increased speed, and the range of the induced current is from −0.5 m to 0.5 m, which is consistent with the previous analysis in Section 3.2.1. Figure 15b shows the variation of the maximum induced current with the velocity. Note that the maximum induced current is proportional to the velocity v z of the point charge.
Equation (13) can be expressed as The induced current can be calculated by

Simulation for Induced Current
As shown in Figure 14, suppose that point charge moves from (0, 0 , 0 ) to (0, 0 , − 0 ) with a different constant velocity of 2 m/s, 5 m/s, 10 m/s, and 20 m/s. The charge of the point charge is = −1 C, 0 = 0.1 m, the length of the electrode is 0.1 m, and the radius of the electrode is r = 0.005 m. As shown in Figure 15a, the induced current increases with the increased speed, and the range of the induced current is from −0.5 m to 0.5 m, which is consistent with the previous analysis in Section 3.2.1. Figure 15b shows the variation of the maximum induced current with the velocity. Note that the maximum induced current is proportional to the velocity of the point charge.

Spectrum of Induced Current
Four velocity cases (2 m/s, 5 m/s, 14.3 m/s, and 20 m/s) are selected to study the spectrum of the induced current. As shown in Figure 14, point charge q carries a charge of −1 C and moves from position 0, y 0 , z 0 to position 0, y 0 , −z 0 at velocity v, where z 0 = 1 m and y 0 = 0.0015 m. The length of the electrode is 0.5 m and the radius is r = 0.005 m. The results of the spectrum analysis are shown in Figure 16. It can be seen that the spectrum of the induced current spreads wider as the speed increases. In other words, as the charge moves faster, more abundant frequency components are introduced. This places a requirement on the design of the induced current measurement circuit. In order to capture the signals induced by faster-moving particles, it is necessary to design an acquisition circuit with a sufficiently wide frequency response, which will be discussed in Section 4.4.

Spectrum of Induced Current
Four velocity cases (2 m/s, 5 m/s, 14.3 m/s, and 20 m/s) are selected to study the spectrum of the induced current. As shown in Figure 14, point charge q carries a charge of −1 C and moves from position (0, y 0 , z 0 ) to position (0, y 0 , −z 0 ) at velocity v, where z 0 = 1 m and y 0 = 0.0015 m. The length of the electrode is 0.5 m and the radius is r = 0.005 m. The results of the spectrum analysis are shown in Figure 16. It can be seen that the spectrum of the induced current spreads wider as the speed increases. In other words, as the charge moves faster, more abundant frequency components are introduced. This places a requirement on the design of the induced current measurement circuit. In order to capture the signals induced by faster-moving particles, it is necessary to design an acquisition circuit with a sufficiently wide frequency response, which will be discussed in Section 4.4.

Analysis of the Variation of the Induced Current Spectrum over the Effective Range
To study the variation of the induced current spectrum in the range shown in Figure 10, we take velocity = 20 m/s as an example. As shown in Figure 14, point charge q moves from (0, y 0 , z 0 ) to (0, 0 , − 0 ) at constant velocity parallel to the z-axis direction, where 0 = 1 m, the charge of the point charge is = −1 C, 0 changes from 0.001 m to 0.5 m, the length of the electrode is 0.5 m, and the radius is = 0.005 m. The vertical axis of Figure 17 represents the maximum frequency component. The result indicates that the higher the frequency of the induced current induced by particles closer to the electrode, the lower the frequency of the current induced by particles far from the electrode. In order to acquire high-frequency signals around the electrodes, the detection circuit requires a wide band-pass.

Analysis of the Variation of the Induced Current Spectrum over the Effective Range
To study the variation of the induced current spectrum in the range shown in Figure 10, we take velocity v = 20 m/s as an example. As shown in Figure 14, point charge q moves from 0, y 0 , z 0 to (0, y 0 , −z 0 ) at constant velocity v parallel to the z-axis direction, where z 0 = 1 m, the charge of the point charge is Q = −1 C, y 0 changes from 0.001 m to 0.5 m, the length of the electrode is 0.5 m, and the radius is r = 0.005 m. The vertical axis of Figure 17 represents the maximum frequency component. The result indicates that the higher the frequency of the induced current induced by particles closer to the electrode, the lower the frequency of the current induced by particles far from the electrode. In order to acquire high-frequency signals around the electrodes, the detection circuit requires a wide band-pass.

Spectrum of Induced Current
Four velocity cases (2 m/s, 5 m/s, 14.3 m/s, and 20 m/s) are selected to study the spectrum of the induced current. As shown in Figure 14, point charge q carries a charge of −1 C and moves from position (0, y 0 , z 0 ) to position (0, y 0 , −z 0 ) at velocity v, where z 0 = 1 m and y 0 = 0.0015 m. The length of the electrode is 0.5 m and the radius is r = 0.005 m. The results of the spectrum analysis are shown in Figure 16. It can be seen that the spectrum of the induced current spreads wider as the speed increases. In other words, as the charge moves faster, more abundant frequency components are introduced. This places a requirement on the design of the induced current measurement circuit. In order to capture the signals induced by faster-moving particles, it is necessary to design an acquisition circuit with a sufficiently wide frequency response, which will be discussed in Section 4.4.

Analysis of the Variation of the Induced Current Spectrum over the Effective Range
To study the variation of the induced current spectrum in the range shown in Figure 10, we take velocity = 20 m/s as an example. As shown in Figure 14, point charge q moves from (0, y 0 , z 0 ) to (0, 0 , − 0 ) at constant velocity parallel to the z-axis direction, where 0 = 1 m, the charge of the point charge is = −1 C, 0 changes from 0.001 m to 0.5 m, the length of the electrode is 0.5 m, and the radius is = 0.005 m. The vertical axis of Figure 17 represents the maximum frequency component. The result indicates that the higher the frequency of the induced current induced by particles closer to the electrode, the lower the frequency of the current induced by particles far from the electrode. In order to acquire high-frequency signals around the electrodes, the detection circuit requires a wide band-pass.

Measurement Circuit Analysis
An electrode typically has a high impedance output, hence a charge amplifier is typically added to match the impedance with the input of the data acquisition circuit. A typical charge amplifier is presented in Figure 18a. After passing through the circuit, the induced current signal is converted to a voltage signal. According to the virtual short and virtual break principle of the operational amplifier, it can be obtained that By performing the Fourier transform on Equations (21)-(23), the response of the circuit can be derived as: The two most important features of the measurement circuit are its amplification and frequency characteristics. The amplification function is used to amplify the weak induced current signal to the extent that it can be acquired by an analogue to a digital converter. To study the spectral influence of the measurement circuit over the sensing circuit, its frequency characteristics are examined. From Equation (24), it can be concluded that the amplification factor of the output signal is determined by the resistor R and the cutoff frequency of the output signal U(t) is determined by the capacitor C. The frequency responses by different R and C values are calculated and shown in Figure 18b. It can be observed that the cutoff frequency decreases with the increased C.

Measurement Circuit Analysis
An electrode typically has a high impedance output, hence a charge amplifier is typically added to match the impedance with the input of the data acquisition circuit. A typical charge amplifier is presented in Figure 18a. After passing through the circuit, the induced current signal is converted to a voltage signal. According to the virtual short and virtual break principle of the operational amplifier, it can be obtained that By performing the Fourier transform on Equations (21)-(23), the response of the circuit can be derived as: The two most important features of the measurement circuit are its amplification and frequency characteristics. The amplification function is used to amplify the weak induced current signal to the extent that it can be acquired by an analogue to a digital converter. To study the spectral influence of the measurement circuit over the sensing circuit, its frequency characteristics are examined. From Equation (24), it can be concluded that the amplification factor of the output signal is determined by the resistor and the cutoff frequency of the output signal ( ) is determined by the capacitor C. The frequency responses by different and values are calculated and shown in Figure 18b. It can be observed that the cutoff frequency decreases with the increased .

Experimental Study and Discussion
From the simulation study in Section 4, we can understand the signals of the induced current by a moving point charge. In this section, an experimental system is set up and three experiments are performed to validate the established model.

Experimental Setup
The schematic of the experimental system is shown in Figure 19a. A rod electrode is placed on a supporting holder parallel to the ground, and a funnel is placed on top of the electrode. During the experiment, a charged ball is released from the funnel toward the ground. By adjusting the distance Height1 between supporting rods 1 and 2, the moving speed of the charged ball can be adjusted. The distance Height2 between supporting rod 2 and the ground can decide the final speed of the charging ball. Note that supporting rod 2 is made of insulating material to hold the rod electrode. In addition, supporting rods 1 and 2 are adjustable to a full 360°, hence the horizontal distance between the small ball and the rod electrode can be adjusted by changing the angle between supporting rods 1 and 2; that is, parameter ℎ in Equation (19).
A picture of the rod electrode employed in this paper is shown in Figure 19b. The rod electrode is made of a 316L stainless steel, with a radius of 0.005 m and a length of 0.5 m. Note that the surface is coated with a Teflon insulation 0.3 mm thick to ensure that the current on the electrode is entirely induced by the electrostatic induction.
The output of the rod electrode was connected to a charge amplifier for amplification, as shown in Figure 18a, with = 100 MΩ and = 10 pf. According to the analysis in Section 4.4, the cutoff frequency was set at 159.2 Hz with such a setting. In the experiment, the maximum moving speed of the ball was less than 6 m/s, so the signal of interest could be appropriately magnified based on the analysis in Figure 16. Then, the signal was acquired by a data acquisition board and sent to the computer for further analysis. From the analysis in Section 4, we understand that the highest frequency of the signal is below 3500 Hz, so the sampling rate of the data acquisition board was set at 10 kHz to satisfy the Nyquist sampling theorem, i.e., the sampling frequency should be more than twice the maximum frequency.
The experiments are conducted in an environment with room temperature (about 20-25 °C) and normal pressure (about 101 kPa).

Experimental Study and Discussion
From the simulation study in Section 4, we can understand the signals of the induced current by a moving point charge. In this section, an experimental system is set up and three experiments are performed to validate the established model.

Experimental Setup
The schematic of the experimental system is shown in Figure 19a. A rod electrode is placed on a supporting holder parallel to the ground, and a funnel is placed on top of the electrode. During the experiment, a charged ball is released from the funnel toward the ground. By adjusting the distance Height1 between supporting rods 1 and 2, the moving speed of the charged ball can be adjusted. The distance Height2 between supporting rod 2 and the ground can decide the final speed of the charging ball. Note that supporting rod 2 is made of insulating material to hold the rod electrode. In addition, supporting rods 1 and 2 are adjustable to a full 360 • , hence the horizontal distance between the small ball and the rod electrode can be adjusted by changing the angle between supporting rods 1 and 2; that is, parameter h in Equation (19).
A picture of the rod electrode employed in this paper is shown in Figure 19b. The rod electrode is made of a 316L stainless steel, with a radius of 0.005 m and a length of 0.5 m. Note that the surface is coated with a Teflon insulation 0.3 mm thick to ensure that the current on the electrode is entirely induced by the electrostatic induction.
The output of the rod electrode was connected to a charge amplifier for amplification, as shown in Figure 18a, with R = 100 MΩ and C = 10 pf. According to the analysis in Section 4.4, the cutoff frequency was set at 159.2 Hz with such a setting. In the experiment, the maximum moving speed of the ball was less than 6 m/s, so the signal of interest could be appropriately magnified based on the analysis in Figure 16. Then, the signal was acquired by a data acquisition board and sent to the computer for further analysis. From the analysis in Section 4, we understand that the highest frequency of the signal is below 3500 Hz, so the sampling rate of the data acquisition board was set at 10 kHz to satisfy the Nyquist sampling theorem, i.e., the sampling frequency should be more than twice the maximum frequency.
The experiments are conducted in an environment with room temperature (about 20-25 • C) and normal pressure (about 101 kPa).

Results and Discussion
Three situations were analyzed in this study, and in all these experiments, the gravity acceleration was assumed to be 9.8 m/s 2 and the ball fell from the funnel at an initial velocity of 0 m/s, hence the velocity of the ball at time can be calculated by ( ) = .  Figure 20a,b, respectively. As can be observed, there exists a DC component in the measured signal in Figure 20a. This DC component is from the data acquisition board, not from the true signal, as the data acquisition board has an added offset of 2.5 V to the input signal, converting the input voltage of -2.5 V to 2.5 V to 0 to 5 V. Hence, this DC component should be removed before further processing. It can be observed that the signal interfered with a severe 50 Hz power line signal and the signal of interest mainly falls within 20 Hz. With a proper low-pass digital filter, the power line frequency can easily be removed. The filtered signal and its spectrum are presented in Figure  20c,d, respectively. In Figure 20c, the voltage increases at the beginning when the ball moves nearer to the rod electrode, indicating that the rate of increase of the induced charge is increasing. Then, the voltage decreases to zero, indicating that the rate of increase of the induced charge is decreasing, and the induced charge reaches the maximum value when the ball is nearest to the electrode, at which point the voltage output is equal to zero.
As the ball moves away from the electrode, the voltage becomes negative, indicating that the induced charge on the electrode becomes smaller. The smaller the voltage, the larger the loss rate. When the loss rate becomes the smallest, then it becomes larger, indicating that the loss rate of the induced charge becomes smaller. At the point where the voltage reaches zero again, the induced charge returns to zero. This is in agreement with the fact that when a charged ball moves to an electrode, the electrode will induce charges, and the induced charges will be reduced when the charged ball moves away from the electrode.

Results and Discussion
Three situations were analyzed in this study, and in all these experiments, the gravity acceleration g was assumed to be 9.8 m/s 2 and the ball fell from the funnel at an initial velocity of 0 m/s, hence the velocity of the ball at time t can be calculated by v(t) = gt.  Figure 20a,b, respectively. As can be observed, there exists a DC component in the measured signal in Figure 20a. This DC component is from the data acquisition board, not from the true signal, as the data acquisition board has an added offset of 2.5 V to the input signal, converting the input voltage of −2.5 V to 2.5 V to 0 to 5 V. Hence, this DC component should be removed before further processing. It can be observed that the signal interfered with a severe 50 Hz power line signal and the signal of interest mainly falls within 20 Hz. With a proper low-pass digital filter, the power line frequency can easily be removed. The filtered signal and its spectrum are presented in Figure 20c,d, respectively. In Figure 20c, the voltage increases at the beginning when the ball moves nearer to the rod electrode, indicating that the rate of increase of the induced charge is increasing. Then, the voltage decreases to zero, indicating that the rate of increase of the induced charge is decreasing, and the induced charge reaches the maximum value when the ball is nearest to the electrode, at which point the voltage output is equal to zero.
As the ball moves away from the electrode, the voltage becomes negative, indicating that the induced charge on the electrode becomes smaller. The smaller the voltage, the larger the loss rate. When the loss rate becomes the smallest, then it becomes larger, indicating that the loss rate of the induced charge becomes smaller. At the point where the voltage reaches zero again, the induced charge returns to zero. This is in agreement with the fact that when a charged ball moves to an electrode, the electrode will induce charges, and the induced charges will be reduced when the charged ball moves away from the electrode. By using Equation (19), the induced current change with the distance between the ball and the funnel is calculated and presented in Figure 21a. The simulated signal output after the charge amplifier is shown in Figure 21b. It can be observed that this signal shows good similarity with the experimental output signal (Figure 20c) in both the shape and amplitude. The induced current increases to a maximum value with the ball moving close to the electrode, and then decreases to zero when the ball is closest to the electrode, corresponding with Figure 15a. The validity of the model is verified.
By resetting the time when the ball leaves the funnel at zero and converting the time to distance in Figure 20c, the theoretical output and measured signal are plotted together in Figure 22a and their spectra are compared in Figure 22b. It shows that these two signals have a high similarity in trend and their values are very close to each other, verifying the correctness of the model we established. Specifically, when the charged ball moved with a low velocity (distance Height 1 is less than 0.8 m or velocity is less than 4 m/s), the theoretical and measured waveforms agreed very well. When the charged ball moved with a high velocity, i.e., velocity greater than 4 m/s, the measured signal was smaller than the theoretical one. A possible reason for this phenomenon could be that the charged ball experiences more resistance from the air with the increased velocity, which reduces the velocity of the ball to a level lower than theoretical calculations and hence less current is induced on the rod electrode. By using Equation (19), the induced current change with the distance between the ball and the funnel is calculated and presented in Figure 21a. The simulated signal output after the charge amplifier is shown in Figure 21b. It can be observed that this signal shows good similarity with the experimental output signal (Figure 20c) in both the shape and amplitude. The induced current increases to a maximum value with the ball moving close to the electrode, and then decreases to zero when the ball is closest to the electrode, corresponding with Figure 15a. The validity of the model is verified.
By resetting the time when the ball leaves the funnel at zero and converting the time to distance in Figure 20c, the theoretical output and measured signal are plotted together in Figure 22a and their spectra are compared in Figure 22b. It shows that these two signals have a high similarity in trend and their values are very close to each other, verifying the correctness of the model we established. Specifically, when the charged ball moved with a low velocity (distance Height 1 is less than 0.8 m or velocity is less than 4 m/s), the theoretical and measured waveforms agreed very well. When the charged ball moved with a high velocity, i.e., velocity greater than 4 m/s, the measured signal was smaller than the theoretical one. A possible reason for this phenomenon could be that the charged ball experiences more resistance from the air with the increased velocity, which reduces the velocity of the ball to a level lower than theoretical calculations and hence less current is induced on the rod electrode.  The theoretical and measured output voltages are shown in Figure 23a. As can be observed, the waveforms show similar trends as those in Experiment 1. The induced current increased accordingly with the increased velocity of the ball, which is reflected by the maximum and minimum amplitude of the induced current. By comparing Figure 23a,c, it can be seen that the induced voltage increased significantly, mainly caused by an increased charge carried by the ball. This can be explained by Equation (19), which shows that the induced current is proportional to the amount of charge if other conditions remain unchanged. Taking the maximum value, minimum value, and measured value at four similar distances in Figure 23a,c, the results are summarized in Table 1. It can be observed that the induced current is approximately proportional to the amount of charge.   The theoretical and measured output voltages are shown in Figure 23a. As can be observed, the waveforms show similar trends as those in Experiment 1. The induced current increased accordingly with the increased velocity of the ball, which is reflected by the maximum and minimum amplitude of the induced current. By comparing Figure 23a,c, it can be seen that the induced voltage increased significantly, mainly caused by an increased charge carried by the ball. This can be explained by Equation (19), which shows that the induced current is proportional to the amount of charge if other conditions remain unchanged. Taking the maximum value, minimum value, and measured value at four similar distances in Figure 23a,c, the results are summarized in Table 1. It can be observed that the induced current is approximately proportional to the amount of charge. The theoretical and measured output voltages are shown in Figure 23a. As can be observed, the waveforms show similar trends as those in Experiment 1. The induced current increased accordingly with the increased velocity of the ball, which is reflected by the maximum and minimum amplitude of the induced current. By comparing Figure 23a,c, it can be seen that the induced voltage increased significantly, mainly caused by an increased charge carried by the ball. This can be explained by Equation (19), which shows that the induced current is proportional to the amount of charge if other conditions remain unchanged. Taking the maximum value, minimum value, and measured value at four similar distances in Figure 23a,c, the results are summarized in Table 1. It can be observed that the induced current is approximately proportional to the amount of charge.   The comparison of the output voltage between the theoretical and experimental waveforms is shown in Figure 24a and the spectrum comparison is shown in Figure 24b. By comparing Figures 24b  and 22a, it can be seen that the high-frequency component decreases when ℎ increases, which is consistent with the analysis in Section 4.3.  The comparison of the output voltage between the theoretical and experimental waveforms is shown in Figure 24a and the spectrum comparison is shown in Figure 24b. By comparing Figures 22a  and 24b, it can be seen that the high-frequency component decreases when h increases, which is consistent with the analysis in Section 4.3. Through three experiments with different scenario settings, it can be observed that the established model matches well with the experimental measurements, verifying the correctness and accuracy of the model. In addition, the following observations are also verified: (1) The induced current decreased significantly with the increase of ℎ, which can be observed in Figures 22b and 24a, and this is in line with the result in Figure 7b.
(2) The high-frequency components of the induced current are reduced with the increased ℎ, which can be seen by comparing Figures 24b and 22a, and this is consistent with Figure 17.
(3) The high-frequency components increase with the increased velocity, which can be observed by comparing Figures 23b and 24b. (4) When the velocity and the ℎ value are fixed, the change in the amount of the charge on the ball does not influence the spectral characteristics of the induced current signal, which can be observed in Figure 23b,d.

Conclusions
In this paper, a theoretical model of the induced charge on a rod electrode by a point charge is established by the method of images and the accuracy of the model is verified by three experimental studies. The following conclusions are drawn based on the study.
1. The amount of the induced charge on the rod electrode are mainly determined by the following factors: The distance between the point charge and the electrode, the radius and the length of the electrode.
2. The general model of the relationship between the induced current and velocity is established and the spectrum of the induced current is studied. The induced current increases with the increase of the point charge's velocity and the maximum value of the induced current is linearly proportional to the point charge's velocity. The faster the velocity of the point charge, the wider the spectrum of the induced current. The further the point charge from the rod electrode, the narrower the spectrum of the induced current on the rod electrode.
3. For the measurement circuit, its amplification ratio and pass-band width are determined by the feedback resistance and feedback capacitance, respectively. With the increase of the feedback resistance, the amplification factor of the circuit increases. With the increase of the feedback capacitance, the pass-band width of the measurement circuit becomes narrower.
The phenomena observed in the experiments, which are well explained by the established model, would be not so easy to explain with just the qualitative analysis. This model provides a basis for more complex studies in the future; for example, the sum result of all the particles around the electrode can been studied with this model. In addition, the influence of the length and radius of the rod electrode on the induced charge are also studied, providing theoretical support for the rod electrode sensor development in the future. Through three experiments with different scenario settings, it can be observed that the established model matches well with the experimental measurements, verifying the correctness and accuracy of the model. In addition, the following observations are also verified: (1) The induced current decreased significantly with the increase of h, which can be observed in Figures 22b and 24a, and this is in line with the result in Figure 7b. (2) The high-frequency components of the induced current are reduced with the increased h, which can be seen by comparing Figures 22a and 24b, and this is consistent with Figure 17. (3) The high-frequency components increase with the increased velocity, which can be observed by comparing Figures 23b and 24b. (4) When the velocity and the h value are fixed, the change in the amount of the charge on the ball does not influence the spectral characteristics of the induced current signal, which can be observed in Figure 23b,d.

Conclusions
In this paper, a theoretical model of the induced charge on a rod electrode by a point charge is established by the method of images and the accuracy of the model is verified by three experimental studies. The following conclusions are drawn based on the study.

1.
The amount of the induced charge on the rod electrode are mainly determined by the following factors: The distance between the point charge and the electrode, the radius and the length of the electrode.

2.
The general model of the relationship between the induced current and velocity is established and the spectrum of the induced current is studied. The induced current increases with the increase of the point charge's velocity and the maximum value of the induced current is linearly proportional to the point charge's velocity. The faster the velocity of the point charge, the wider the spectrum of the induced current. The further the point charge from the rod electrode, the narrower the spectrum of the induced current on the rod electrode.

3.
For the measurement circuit, its amplification ratio and pass-band width are determined by the feedback resistance and feedback capacitance, respectively. With the increase of the feedback resistance, the amplification factor of the circuit increases. With the increase of the feedback capacitance, the pass-band width of the measurement circuit becomes narrower.
The phenomena observed in the experiments, which are well explained by the established model, would be not so easy to explain with just the qualitative analysis. This model provides a basis for more complex studies in the future; for example, the sum result of all the particles around the electrode can