High-Precision Inertial Sensor Charge Ground Measurement Method Based on Phase-Sensitive Demodulation

Liu, Yang; Yu, Tao; Wang, Yuhua; Zhao, Zihan; Wang, Zhi

doi:10.3390/s24031009

Open AccessArticle

High-Precision Inertial Sensor Charge Ground Measurement Method Based on Phase-Sensitive Demodulation

by

Yang Liu

^1,2,

Tao Yu

^1,*

,

Yuhua Wang

^1,2,

Zihan Zhao

^1,2 and

Zhi Wang

^1,3,*

¹

Changchun Institute of Optics, Fine Mechanics and Physics, Chinese Academy of Sciences, Changchun 130033, China

²

University of Chinese Academy of Sciences, Beijing 100049, China

³

School of Fundamental Physics and Mathematical Sciences, Hangzhou Institute for Advanced Study, University of Chinese Academy of Sciences, Hangzhou 310024, China

^*

Authors to whom correspondence should be addressed.

Sensors 2024, 24(3), 1009; https://doi.org/10.3390/s24031009

Submission received: 27 November 2023 / Revised: 22 January 2024 / Accepted: 1 February 2024 / Published: 4 February 2024

(This article belongs to the Section Electronic Sensors)

Download

Browse Figures

Review Reports Versions Notes

Abstract

:

Inertial sensors are the key payloads in space gravitational wave detection missions, and they need to ensure that the test mass (TM), which serves as the inertial reference, freely floats in the spacecraft without contact, so that the TM is not disturbed by the satellite platform and the cosmic environment. Space gravitational wave detection missions require that the residual acceleration of the TM should be less than

3 \times 10^{- 15} {ms}^{- 2} {Hz}^{- 1 / 2}

. However, the TM with charges will interact with surrounding conductors and magnetic fields, introducing acceleration noise such as electrostatic force and Lorentz force. Therefore, it is necessary to carry out charge management on the TM, in which the high-precision measurement of charge is crucial. Space gravitational wave detection missions require a residual charge measurement accuracy of

3 \times 10^{- 13} C

for the TM. In this paper, we design a high-precision inertial sensor charge measurement method based on phase-sensitive demodulation (PSD). By establishing a torsion pendulum rotation model based on the force modulation method, the characteristics of the TM torsion angle signal are analyzed. The PSD is used to extract the amplitude of the specific frequency signal component containing the charge information, and then to calculate the value of the accumulated charges. The method is compared with the Butterworth band-pass filtering method, and the simulation results show that the method has a higher measurement accuracy, shorter settling time, and stronger anti-interference ability, meeting the TM residual charge measurement accuracy index requirement.

Keywords:

space gravitational wave detection; inertial sensor; charge measurement; force modulation method; phase-sensitive demodulation

1. Introduction

Inertial sensors are the key payloads in space gravitational wave detection missions, and their cores are a sensitive structure consisting of the TM and the electrode housing (EH) [1]. The TM serves as the spacecraft’s inertial reference, and the inertial sensors need to ensure that the TM is freely suspended within the spacecraft without contact [2]. Although the spacecraft is able to block most of the external environmental disturbances, some cosmic rays and energetic particles can penetrate the spacecraft and directly interact with the TM, resulting in the accumulation of charges on the TM.

The charges accumulated on the TM interact with the surrounding conductors to generate electrostatic forces. In addition, the accumulated charges also couple with the interplanetary magnetic fields and spacecraft internal magnetic fields to generate the Lorentz force. These all result in additional acceleration noises [3]. Therefore, charge management is necessary for the TM. With charge management, the residual charges on the TM are required to be less than

2 \times 10^{7} e

(where e represents an electron,

1 e = 1.6 \times 10^{- 19} C

and

2 \times 10^{7} e = 3.2 \times 10^{- 12} C

) for LISA (Laser Interferometer Space Antenna) [4], and less than

10^{6} e

(

1.6 \times 10^{- 13} C

) for Tianqin [5,6], so the residual charge measurement accuracy of

3 \times 10^{- 13} C

is required for the TM.

The charge accumulation on the TM is mainly caused by galactic cosmic rays (GCRs) and solar high-energy particles (SEPs). According to the existing studies on the rate of charge accumulation on the TM, the rate of charge accumulation on the TM due to GCRs is in the range of

+ 10 e s^{- 1}

to

+ 100 e s^{- 1}

[7,8,9,10,11,12,13,14,15], assuming

+ 100 e s^{- 1}

as the maximum rate for the calculation, the TM can be charged to a level of

6 \times 10^{7} e

(

9.6 \times 10^{- 12} C

) in about one week. On the other hand, a large SEP event that can reach a maximum charging rate of about +70,000

e s^{- 1}

may only occur once in a solar cycle, with a total charge deposition of

+ 2 \times 10^{9} e

(

3.2 \times 10^{- 10} C

) during a five-day SEP event [13]. Therefore, for short-term charge accumulation in the space environment, the total amount of charge can be focused on less than or equal to

10^{- 10} C

.

The existing research on charge measurement and estimation methods mainly include the force modulation charge measurement method, which has been verified to be feasible by the GP-B (Gravity Probe B), LISA Pathfinder, etc. [16,17], as well as the disturbance-free charge estimation method proposed by the Tianqin team [18]. The force modulation method is a commonly used measurement method in space missions, which is based on the principle of applying signals to the electrodes around the TM along the sensitive axis to stimulate the motion of the TM, and estimating the charge through its position and attitude changes. However, charge measurement using the force modulation method carries the risk of interrupting gravitational wave detection. To address this problem, the Tianqin team proposed a disturbance-free charge estimation method. By analyzing the motion of the TM, the accumulated charge information embedded in its motion signal can be solved. And the method needs to extract specific frequency terms containing charge information from multiple AC components, but due to the similarity of the frequencies of the AC components, it is difficult to accurately extract specific frequency components through a simple filter, which will affect the accuracy of the results. The effectiveness of this method has only been verified through simulations.

To fulfill the ground testing requirements of inertial sensors, this paper integrates the characteristics of the ground testing device, a torsion pendulum, and establishes the torsion pendulum rotation model using the force modulation method. On this basis, a thorough analysis is conducted on the characteristics of the torsion angle signal. This paper has designed a charge measurement method based on PSD, which can resolve the charges accumulated on the TM with a higher accuracy. PSD uses reference signals of the same frequency as the excitation signal to demodulate the amplitude and phase information of the measured signal, which has the advantage of effectively suppressing the interference or noise that is not of the same frequency as the excitation frequency [19]. The amplitude of the specific frequency component containing the charge information can be extracted more accurately by PSD, resulting in a higher accuracy and meeting the ground measurement needs of the charge management system on inertial sensors.

2. Rotational Model of the Torsion Pendulum Based on the Force Modulation Method

In space, the TM inside the inertial sensor is in free suspension, while on the ground due to the influence of Earth’s gravity, such a free state is difficult to realize, so the ground test of the charge measurement system is a very difficult task. A torsion pendulum can simulate the free suspension state of the TM in space due to its special structure, which has a free rotational degree of freedom around the z-axis. Moreover, the torsion pendulum is sensitive to the tiny torques in the horizontal direction, and it can measure very small forces with high accuracy [20]. Therefore, the torsion pendulum can be used to test and evaluate the key technologies and performance of inertial sensors on the ground.

2.1. TM Torque Driver Model

The TM is suspended by a suspension fiber to form the most basic torsion pendulum system. The TM consists of a hollow gold-coated Al cube, which is surrounded by eighteen electrodes, including twelve sensing/actuation electrodes (grey) and six injection electrodes (yellow) distributed as shown in Figure 1, with two sensing/actuation electrodes on each face, one injection electrode on each y face and two injection electrodes on each z face [21]. During operation, the position and attitude of the TM inside the EH are sensed and controlled by the six-degrees-of-freedom differential capacitor pairs consisting of twelve sensing/actuation electrodes [22,23].

The TM of mass m and the electrodes on the EH form parallel capacitors that control the TM by applying the voltages to the sensing electrodes. The electrostatic force between the TM and the electrodes can be described as follows:

F = \frac{1}{2} \frac{\partial C_{i}}{\partial q} {(u_{i} - u_{T M})}^{2}

(1)

where

C_{i}

is the capacitance formed by the TM and the ith electrode, q is a general coordinate,

u_{i}

is the voltage on the ith electrode, and

u_{T M}

is the surface potential of the TM, which can be expressed as follows:

u_{T M} = \frac{Q_{T M}}{C_{T}} + \frac{1}{C_{T}} \sum_{i} C_{i} u_{i}

(2)

where

Q_{T M}

is the charges accumulated on the TM and

C_{T}

is the total capacitance inside the EH.

Let the voltages applied on the four electrodes on the sensitive x axis be

u_{i}

,

i = 1, \dots, 4

, and the electrostatic force between TM and electrode 1 can be expressed as follows:

F_{1} = \frac{ε S}{2 d_{0}^{2}} {(u_{1} - u_{T M})}^{2}

(3)

where

ε

is the vacuum dielectric constant,

ε = 8.85 \times 10^{- 12} F / m

, S is the effective area of the capacitor consisting of the TM and the electrodes, and

d_{0}

is the distance between the TM and the electrodes when the TM is at the center of the EH.

The forces between the TM and the three remaining electrodes can also be obtained by the same method. By applying specific polarity voltages to the four electrodes surrounding the TM along the x-axis direction, the TM can rotate in

φ

degree of freedom, as shown in Figure 2.

For the TM rotation, the voltages applied to each electrode in the x axis are as follows:

\{\begin{matrix} u_{1 φ} = {+ u}_{T S, 1 φ} sin (2 π f_{φ} t) \\ u_{2 φ} = {- u}_{T S, 1 φ} sin (2 π f_{φ} t) \\ u_{3 φ} = {+ u}_{T S, 1 φ} sin (2 π f_{φ} t) \\ u_{4 φ} = {- u}_{T S, 1 φ} sin (2 π f_{φ} t) \end{matrix}

(4)

Assuming that the TM is at the center of the EH, there are the same capacitances formed by the TM and the four electrodes in the x-axis direction, so,

u_{T M} = \frac{Q_{T M}}{C_{T}} + \frac{1}{C_{T}} \sum_{i = 1}^{4} C_{i} u_{i} = \frac{Q_{T M}}{C_{T}}

(5)

The TM is subjected to a combined electrostatic torque of

\begin{matrix} M_{φ} = \frac{b}{2} [(F_{1} - F_{3}) - (F_{2} - F_{4})] = \frac{2 ε b S}{d_{0}^{2}} u_{T M} u_{T S, 1 φ} sin (2 π f_{φ} t) \end{matrix}

(6)

where b is the distance between the centers of the two electrodes on the same side,

A = \frac{2 ε b S}{d_{0}^{2}} u_{T M} u_{T S, 1 φ}

.

2.2. Dynamic Model and Frequency Response Analysis of the Torsion Pendulum

Neglecting the noise interference, the torsion pendulum is subjected to three effects:

The torque $τ (t)$ of interaction between the charged TM and the surrounding conductors;
The damping torque of vibration;
The elastic restoring torque.

Usually, the torsion pendulum system works in a high vacuum environment, the air damping is much smaller than the structural damping of the system, and the dynamics model of the torsion pendulum can be reduced to a second-order mass-spring-damper model, the dynamics equation of which is

I_{r} \ddot{φ} (t) + γ \dot{φ} (t) + Γ φ (t) = τ (t)

(7)

where

I_{r}

is the rotational inertia of the torsion pendulum,

γ

is the damping coefficient,

γ = Γ / (ω_{n} Q)

, Q is the quality factor of the torsion pendulum,

Γ

is the restoring stiffness of the suspension fiber,

ω_{n} = \sqrt{(Γ / I_{r}))}

is the intrinsic frequency of the torsion pendulum,

φ

is the torsion angle, and

τ

is the torque.

The Laplace transform of Equation (7) gives

I_{r} s^{2} φ (s) + γ s φ (s) + Γ φ (s) = τ (s)

(8)

where

φ (s)

and

τ (s)

are the Laplace transforms of the torsion angle and torque, respectively, and

s = j ω

is the Laplace variable.

Then the transfer function

H_{τ, φ} (s)

between the torque and torsion angle is

H_{τ, φ} (s) = \frac{φ (s)}{τ (s)} = \frac{1}{I_{r} s^{2} + γ s + Γ}

(9)

The standard form of the second-order system is

G (s) = \frac{ω_{n}^{2}}{s^{2} + 2 ξ ω_{n} s + ω_{n}^{2}}

, where

ξ = γ / (2 \sqrt{(I_{r} Γ)})

is the damping ratio. Equation (9) can be further transformed into the following form:

H_{τ, φ} (s) = \frac{1}{Γ} \cdot \frac{ω_{n}^{2}}{s^{2} + 2 ξ ω_{n} s + ω_{n}^{2}}

(10)

Let

s = j ω

in Equation (10), then,

H_{τ, φ} (j ω) = \frac{1}{Γ} \cdot \frac{1}{{- (\frac{ω}{ω_{n}})}^{2} + 2 ξ \frac{ω}{ω_{n}} j + 1}

(11)

Then, let

Ω = ω / ω_{n}

, which yields

\begin{matrix} H_{τ, φ} (j ω) & = \frac{1}{Γ} \cdot (\frac{{1 - Ω}^{2}}{{(1 - Ω^{2})}^{2} + 4 ξ^{2} Ω^{2}} - \frac{2 ξ Ω}{{(1 - Ω^{2})}^{2} + 4 ξ^{2} Ω^{2}} j) \end{matrix}

(12)

So, the amplitude response of the torsion pendulum is

|H_{τ, φ} (j ω)| = \frac{1}{Γ} \cdot \sqrt{\frac{1}{{(1 - Ω^{2})}^{2} + 4 ξ^{2} Ω^{2}}}

(13)

where its amplitude response is

∠ H_{τ, φ} (j ω)

.

2.3. Rotational Model of the Torsion Pendulum Based on the Force Modulation Method

The sinusoidal voltages of certain amplitude and frequency are applied to the four sensing electrodes around the TM along the x axis, and due to the existence of the potential difference between the electrodes and the TM, electrostatic interactions will occur, the TM produces a torsional motion, and the torsion angle of the TM can be measured. On the basis of the above mathematical derivation, a model is built in Matlab/Simulink, as shown in Figure 3.

The inputs of the model are the initial set charge value

Q_{1}

, as well as the amplitude

u_{T S, 1 φ}

and frequency

f_{φ}

of the excitation applied to the electrodes, and the output is the torsion angle of the TM, and the parameters in the model are shown in Table 1. Based on the simulation of the model and the analysis of the angle signal both in the time domain and frequency domain, it can be found that the angle signal mainly consists of two frequency components, which are the frequency

f_{φ}

of the applied excitation and the intrinsic frequency

f_{n}

of the torsion pendulum.

According to Equations (6) and (13), if applying the excitation signal with the amplitude of

u_{T S, 1 φ}

and the frequency of

f_{φ}

to each electrode on the x axis, we can obtain the output angle signal ideally as

φ (t) = A_{1} sin (2 π f_{φ} t + φ_{1}) + A_{2} sin (2 π f_{n} t + φ_{2})

(14)

where the signal component with frequency

f_{φ}

has amplitude

A_{1}

, and the signal component with frequency

f_{n}

has amplitude

A_{2}

, with

A_{1} = A \cdot |H_{τ, φ} (j ω)| = \frac{2 ε b S}{d_{0}^{2}} u_{T M} u_{T S, 1 φ} |H_{τ, φ} (j ω)|

.

Also, since

Q_{T M} = u_{T M} C_{T}

, the information about the accumulated charges on the TM is embedded in the amplitude of the signal component with frequency

f_{φ}

in the output angle signal. Therefore, it is necessary to consider extracting the amplitude of the signal component with frequency

f_{φ}

accurately by analyzing and processing the output angle signal.

The realization of the torsion pendulum is to maximize its sensitivity in the millihertz frequency range of interest for space gravitational wave detection, so that its resonance (the maximum angle response of the pendulum to a given torque) is in the millihertz frequency range [24]. In order to obtain a higher sensitivity of the torsion pendulum, the frequency of the applied excitation also needs to be chosen in the millihertz range. Taking the sinusoidal excitation with amplitude

u_{T S, 1 φ} = 3 V

and frequency

f_{φ} = 0.003 Hz

applied to the model as an example, the resulting angle waveform is shown in Figure 4.

Through performing a frequency domain analysis, as shown in Figure 5, the superposition of two sine waves with frequencies of

f_{n} = 0.00177 Hz

and

f_{φ} = 0.003 Hz

can be observed.

Due to the small difference between the frequencies, it is difficult to accurately extract the amplitude of the signal component the frequency of which is

f_{φ}

through a simple filter. Therefore, the PSD method is selected.

3. Charge Calculation Method Based on PSD

3.1. Design of Charge Calculation Method Based on PSD

PSD is based on the principle of signal cross-correlation, which determines the correlation between two signals by calculating their cross-correlation function. When the frequencies of the two signals are the same, the cross-correlation function is the periodic components of the double frequency. When the frequencies of the two signals are different, the cross-correlation function is zero, indicating that the two signals are uncorrelated. The principle diagram is shown in Figure 6.

In Figure 6,

V_{i n}

is the input signal,

r_{s i n}

and

r_{c o s}

are the sine and cosine reference signals,

d_{r}

,

d_{q}

are the in-phase and quadrature components obtained by demodulation, and A and

φ

are the amplitude and phase of the signal components at the desired specific frequency, respectively. The analog multiplier demodulation uses analog multipliers and low-pass filters to perform cross-correlation operations on the signals, while the digital multiplier demodulation method uses multiplication accumulators [25,26].

PSD is performed with the ideal torsion angle signal as the input signal, and the sine and cosine reference signals are introduced as

r_{1} (t) = sin (2 π f_{φ} t - ϕ)

,

r_{2} (t) = cos (2 π f_{φ} t - ϕ)

, so the results of the cross-correlation operations on the torsion angle signal and the sine and cosine reference signals are as follows:

\begin{matrix} d_{r} (t) & = φ (t) \cdot r_{1} (t) \\ = \frac{A_{1}}{2} cos ϕ - \frac{A_{1}}{2} cos ϕ cos (4 π f_{φ} t) - \frac{A_{1}}{2} sin ϕ sin (4 π f_{φ} t) + \\ A_{2} cos ϕ sin (2 π f_{φ} t) sin (2 π f_{n} t) - A_{2} sin ϕ sin (2 π f_{n} t) cos (2 π f_{φ} t) \end{matrix}

(15)

and

\begin{matrix} d_{q} (t) & = φ (t) \cdot r_{2} (t) \\ = \frac{A_{1}}{2} sin ϕ + \frac{A_{1}}{2} cos ϕ sin (4 π f_{φ} t) - \frac{A_{1}}{2} sin ϕ cos (4 π f_{φ} t) + \\ A_{2} cos ϕ sin (2 π f_{n} t) cos (2 π f_{φ} t) + A_{2} sin ϕ sin (2 π f_{n} t) sin (2 π f_{φ} t) \end{matrix}

(16)

From Equations (15) and (16), it can be found that the cross-correlation operation results contain a DC component and AC components, while the DC component amplitude contains the charge information. And the DC component amplitude can be obtained by filtering out the AC components using a low-pass filter, which can then be used to calculate the value of the charges.

After filtering out the AC components using a low-pass filter, we can obtain

\{\begin{matrix} d_{r} = \frac{A_{1}}{2} cos ϕ \\ d_{q} = \frac{A_{1}}{2} sin ϕ \end{matrix}

(17)

where

d_{r}

is called the in-phase component and

d_{q}

is called the quadrature component, from which the amplitude and phase of the signal component at the desired specific frequency can be calculated.

\{\begin{matrix} \hat{A_{1}} = 2 \sqrt{d_{r}^{2} + d_{q}^{2}} \\ \hat{ϕ} = a r c t a n \frac{d_{q}}{d_{r}} \end{matrix}

(18)

The charge

Q_{T M}

can also be solved from

A_{1} = \frac{2 ε b S}{d_{0}^{2}} u_{T M} u_{T S, 1 φ} |H_{τ, φ} (j ω)|

and

Q_{T M} = u_{T M} C_{T}

:

Q = \frac{d_{0}^{2} C_{T} \sqrt{d_{r}^{2} + d_{q}^{2}}}{ε b S u_{T S, 1 φ} |H_{τ, φ} (j ω)|}

(19)

If the phase of the output angle signal is not zero, similarly, it can be deduced that the results of the cross-correlation contain DC components and AC components, and the amplitudes of the DC components can be obtained through a low-pass filter. However, the DC amplitudes also contain the sine and cosine values of the phase of the angle signal, so the phase of the output angle signal needs to be detected first.

3.2. Model of Charge Measurement in the Torsion Pendulum Based on PSD

The following simulation model is established in MATLAB/Simulink, as shown in Figure 7. The inputs of the model are the initial set charge value

Q_{1}

, the amplitude

u_{T S, 1 φ}

, and frequency

f_{φ}

of the excitation applied to the electrodes, and the output is the estimated charge value

Q_{e}

. The TM is twisted by applying excitation to the four electrodes around the TM along the sensitive axis, and the resulting torque is converted into an angle signal by the transfer function of the torsion pendulum, which is then subjected to PSD to calculate the value of the charges.

The initial charge setting value

Q_{1}

is taken as the agreed true value and the simulation value

Q_{e}

is taken as the measurement value, so the absolute error is

Δ Q = Q_{e} - Q_{1}

and the relative error is

δ = \frac{Q_{e} - Q_{1}}{Q_{1}}

.

3.3. Performance Verification

In order to verify the effectiveness and anti-interference performance of the designed charge measurement method based on PSD, we compared the method with three common band-pass filtering methods, namely the Butterworth BPF, the Chebyshev II BPF, and the Bessel BPF in the ideal case and with noise introduced, respectively. We chose two representative charge values as the initial set charge values,

1 \times 10^{- 12} C

and

1 \times 10^{- 13} C

. The simulation results of the different methods in the ideal case and with noise introduced are summarized in Table 2, respectively.

3.3.1. In the Ideal Case

Figure 8 shows the simulation results of these four methods in the ideal case, where the black curve represents the simulation result obtained using the designed charge measurement method based on PSD, the red curve represents Butterworth BPF, the yellow curve represents ChebyshevII BPF and the purple curve represents Bessel BPF, respectively. The results show that in the ideal case, regardless of the initial charge setting value

Q_{1} = 1 \times 10^{- 12} C

or

Q_{1} = 1 \times 10^{- 13} C

, the estimation error of PSD is smaller than that of the other three methods. The time required for the charge simulation results to reach the initial charge setting value and remain within ±2% error is defined as the settling time (

T s

), the

T s

required by PSD is also shorter than that required by Butterworth BPF or ChebyshevII BPF. Therefore, from a comprehensive perspective, PSD has good effects in both Ts and accuracy.

By giving multiple sets of initial set charge values

Q_{1} \in [1, 1 \times 10^{- 25}] C

, multiple sets of simulation results are obtained, and the absolute and relative errors under different initial set charge values are plotted as shown in Figure 9a,b, respectively. The results show that both the absolute and relative errors of the PSD results are significantly better than those of the other three methods results across the full range.

The root mean square error (

R M S E

) is the square root of the ratio of the sum of the squares of the deviations of the measured value from the true value to the number of observations m. It is used to measure the deviation between the measured value and the true value. In addition, the mean absolute error (

M A E

) is the average absolute value of the error, and the mean absolute percentage error (

M A P E

) is the average relative error. The smaller the

R M S E

,

M A E

and

M A P E

, the higher the measurement accuracy. Their formulas are as follows:

R M S E = \sqrt{\frac{\sum_{i = 1}^{m} {(x (i) - y (i))}^{2}}{m}}

(20)

M A E = \frac{1}{m} \sum_{i = 1}^{n} | x (i) - y (i) |

(21)

M A P E = \frac{100 %}{m} \sum_{i = 1}^{n} | \frac{x (i) - y (i)}{y (i)} |

(22)

where

x (i)

is the ith measurement value,

y (i)

is the ith truth value, and m is the number of observations.

For the short-term charge accumulation in the space environment, the total amount of charges that is mainly concerned in the charge measurement should be less than or equal to

1 \times 10^{- 10} C

. In the ideal case, both the MAE and MAPE of the method based on PSD are less than those of the other three methods, as shown in Table 2. Most importantly, the

R M S E

obtained by the PSD method is

1.256 \times 10^{- 13} C

, which is smaller than the TM residual charge measurement accuracy index of

3 \times 10^{- 13} C

. This shows that the designed charge measurement method based on PSD meets the TM charge measurement accuracy index requirement.

3.3.2. With Charge Noise and Angle Sensor Noise Introduced

The charge accumulation on the TM is caused by two main types of particle radiation, GCRs and SEPs, which are driven by certain types of solar events. Although the charging rate caused by GCR events is relatively stable, it also varies with the solar activity cycle. The GCR flux reaches its maximum during the solar minimum and its minimum during the solar maximum. On the other hand, SEP events are short lived and originate from solar eruptions. SEPs vary with solar activity, lasting from a few hours to a few days, but can temporarily increase the cosmic ray flux by several orders of magnitude. SEP events are expected to occur only a few times per year during the solar maximum, and less than once per year during the solar minimum. In addition, due to the complexity of the cosmic environment, there are still unknown sources of charge accumulation in space other than GCRs and SEPs, which can cause variations in the charge accumulation on the TM. Considering these factors that cause charge fluctuations, charge noise is introduced into the model of charge measurement in the torsion pendulum based on PSD. Assuming that the initial set charge value carries a zero-mean Gaussian white noise with variance

σ = 10^{- 13} C

, the noise curve is shown in Figure 10.

In addition, in the TM torsion angle measurement step, noise may be introduced due to the limitations of measurement methods and measurement devices, so the angle sensor noise is introduced in the the model, and the angle signal output from the torsion pendulum rotation model based on the force modulation method is combined with a zero-mean, Gaussian white noise with a variance of

σ = 10^{- 5} rad

as the real angle signal. The noise curve is shown in Figure 11.

With the introduction of charge noise and angle sensor noise, the designed charge measurement method based on PSD is compared with the normal filtering methods using Butterworth BPF, ChebyshevII BPF, and Bessel BPF, similarly. Figure 12 shows the simulation results of these four methods with noise introduced. The black curve represents the simulation result obtained using PSD and the red curve represents Butterworth BPF, the yellow curve represents ChebyshevII BPF, and the purple curve represents Bessel BPF, respectively. The simulation results are shown in Table 2. The results show that when noise is introduced, regardless of the initial charge setting value

Q_{1} = 1 \times 10^{- 12} C

or

Q_{1} = 1 \times 10^{- 13} C

, the estimation error of the PSD is significantly smaller than that of the other three methods, with stronger anti-interference capability and shorter

T s

required.

By giving multiple sets of initial set charge values

Q_{1} \in [1, 1 \times 10^{- 15}] C

, multiple sets of simulation results are obtained, and the absolute and relative errors under different initial set charge values are plotted as shown in Figure 13a,b, respectively. The results show that the results of PSD are significantly better than those of the other three methods in terms of both the absolute and relative errors.

Under the consideration of the charge noise and angle sensor noise, for the charge measurement, the main concern of the charge magnitude, i.e. less than or equal to

1 \times 10^{- 10} C

, both the

M A E

and

M A P E

of the method based on PSD are less than those of the other three methods, as shown in Table 2. Most importantly, the

R M S E

obtained by PSD is

2.051 \times 10^{- 13} C

, which is less than the TM residual charge measurement accuracy index of

3 \times 10^{- 13} C

. This shows that the designed charge measurement method based on PSD meets the TM charge measurement accuracy index requirement.

4. Conclusions

The accurate measurement of the charges accumulated on the TM is of great significance. Firstly, it is conducive to the analysis of the various noises caused by the accumulated charges on the TM, which provides a reference for the construction of high-precision inertial sensors. Secondly, it provides a basis for the subsequent charge control. By accurately measuring the charges accumulated on the TM, we can determine the control mode that the TM charge management system should adopt. This paper studies the ground measurement method of inertial sensor charges and completes the establishment and simulation analysis of the torsion pendulum rotation model based on the force modulation method. By analyzing the characteristics of the TM torsion angle signal, a high-precision inertial sensor charge measurement method based on PSD is designed. The results show that when in the ideal case and considering the charge noise and angle sensor noise, the method meets the TM charge measurement accuracy index requirement, with a high accuracy and strong anti-interference ability.

Author Contributions

Conceptualization, Y.L.; methodology, Y.L. and T.Y.; software, Y.L.; validation, Y.L., T.Y. and Y.W.; formal analysis, Y.L. and Y.W.; investigation, Y.L., Y.W. and Z.Z.; resources, Z.W. and T.Y.; writing—original draft preparation, Y.L.; writing—review and editing, T.Y., Y.W., Z.Z. and Z.W.; supervision, T.Y. and Z.W.; project administration, T.Y. and Z.W.; funding acquisition, Z.W. and T.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the National Key R&D Program of China (2020YFC2200604), the External Cooperation Program of Chinese Academy of Sciences (181722KYSB20190040), and the National Key R&D Program of China (2020YFC2200600).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data and the source code are publicly available at https://github.com/YangLiucc/High-precision-Charge-Measurement-Method-Based-on-PSD.git (accessed on 24 November 2023).

Conflicts of Interest

The authors declare no conflicts of interest.

References

Luo, J.; Chen, L.S.; Duan, H.Z.; Gong, Y.G.; Hu, S.; Ji, J.; Liu, Q.; Mei, J.; Milyukov, V.; Sazhin, M.; et al. TianQin: A space-borne gravitational wave detector. Class. Quantum Gravity 2016, 33, 035010. [Google Scholar] [CrossRef]
Li, H.; Liu, Y.; Wang, C.; Bai, Y.Z.; Liu, L.; Wu, S.; Qu, S.; Tan, D.; Yin, H.; Li, Z.; et al. Preliminary Design Thoughts and Progress of Tianqin Inertial Sensors. J. Sun Yat-Sen Univ. (Nat. Sci. Ed.) 2021, 60, 186–193. [Google Scholar] [CrossRef]
Shaul, D.; Araújo, H.; Rochester, G.; Sumner, T.; Wass, P. Evaluation of disturbances due to test mass charging for LISA. Class. Quantum Gravity 2005, 22, S297. [Google Scholar] [CrossRef]
Sumner, T.; Shaul, D.; Schulte, M.; Waschke, S.; Hollington, D.; Araújo, H. LISA and LISA Pathfinder charging. Class. Quantum Gravity 2009, 26, 094006. [Google Scholar] [CrossRef]
Jun, L.; Linghao, A.; Ai Yanli, E.A. Introduction to the Tianqin Project. J. Sun Yat-Sen Univ. (Nat. Sci. Ed.) 2021, 60, 19. [Google Scholar]
Fangchao, Y. Investigation of Charge Measurement and Control for Space Inertial Sensors. Master’s Thesis, Huazhong University of Science and Technology, Wuhan, China, 2023. [Google Scholar]
Araújo, H.; Wass, P.; Shaul, D.; Rochester, G.; Sumner, T. Detailed calculation of test-mass charging in the LISA mission. Astropart. Phys. 2005, 22, 451–469. [Google Scholar] [CrossRef]
Jafry, Y.; Sumner, T.; Buchman, S. Electrostatic charging of space-borne test bodies used in precision experiments. Class. Quantum Gravity 1996, 13, A97. [Google Scholar] [CrossRef]
Araujo, H.M.; Howard, A.; Davidge, D.; Sumner, T.J. Charging of isolated proof masses in satellite experiments such as LISA. In Proceedings of the Gravitational-Wave Detection, Waikoloa, HI, USA, 22–28 August 2002; Volume 4856, pp. 55–66. [Google Scholar]
Araújo, H.; Howard, A.; Shaul, D.; Sumner, T. Electrostatic charging of cubic test masses in the LISA mission. Class. Quantum Gravity 2003, 20, S311. [Google Scholar] [CrossRef]
Shaul, D.; Sumner, T.; Araújo, H.; Rochester, G.; Wass, P.; Lee, C.G. Unwanted, coherent signals in the LISA bandwidth due to test mass charging. Class. Quantum Gravity 2004, 21, S647. [Google Scholar] [CrossRef]
Simpson, J. Elemental and isotopic composition of the galactic cosmic rays. Annu. Rev. Nucl. Part. Sci. 1983, 33, 323–382. [Google Scholar] [CrossRef]
Wass, P.; Araújo, H.; Shaul, D.; Sumner, T. Test-mass charging simulations for the LISA Pathfinder mission. Class. Quantum Gravity 2005, 22, S311. [Google Scholar] [CrossRef]
Vocca, H.; Grimani, C.; Amico, P.; Bosi, L.; Marchesoni, F.; Punturo, M.; Travasso, F.; Barone, M.; Stanga, R.; Vetrano, F.; et al. Simulation of the charging process of the LISA test masses due to solar flares. Class. Quantum Gravity 2004, 21, S665. [Google Scholar] [CrossRef]
Han, R.; Cai, M.; Yang, T.; Xu, L.; Xia, Q.; Han, J. Mechanism of test mass charging by cosmic ray energetic particles. Acta Phys. Sin. 2021, 70, 366–374. [Google Scholar] [CrossRef]
Buchman, S.; Quinn, T.; Keiser, G.; Gill, D.; Sumner, T. Charge measurement and control for the Gravity Probe B gyroscopes. Rev. Sci. Instrum. 1995, 66, 120–129. [Google Scholar] [CrossRef]
Carbone, L.; Cavalleri, A.; Dolesi, R.; Hoyle, C.; Hueller, M.; Vitale, S.; Weber, W.J. Characterization of disturbance sources for LISA: Torsion pendulum results. Class. Quantum Gravity 2005, 22, S509. [Google Scholar] [CrossRef]
Li, H.; Hong, W.; Xiao, C.; Li, Z.; Bai, Y.; Zhou, Z. A Continuous Charge Estimation for Gravitational Wave Detections. In Proceedings of the 2021 IEEE 16th Conference on Industrial Electronics and Applications (ICIEA), Chengdu, China, 1–4 August 2021; pp. 1160–1163. [Google Scholar]
Fusheng, Y.; Xiuzhen, D.; Xuetao, S.; Feng, F.; Ruigang, L.; Mingxin, Q.; Mengxing, T.; Jiang, L. Research on Bioelectrical Impedance Analog Demodulation Technology. Beijing Biomed. Eng. 2004, 23, 4. [Google Scholar]
Yin, H. Research on conductor surface potential measurement based on electrostatic controlled torsion pendulum scheme. Ph.D. Thesis, Huazhong University of Science and Technology, Wuhan, China, 2015. [Google Scholar]
Yu, T.; Wang, Y.; Liu, Y.; Wang, Z. High-Precision Inertial Sensor Charge Management Based on Ultraviolet Discharge: A Comprehensive Review. Sensors 2023, 23, 7794. [Google Scholar] [CrossRef] [PubMed]
Armano, M.; Audley, H.; Baird, J.; Born, M.; Bortoluzzi, D.; Cardines, N.; Castelli, E.; Cavalleri, A.; Cesarini, A.; Cruise, A.; et al. Analysis of the accuracy of actuation electronics in the laser interferometer space antenna pathfinder. Rev. Sci. Instrum. 2020, 91, 045003. [Google Scholar] [CrossRef] [PubMed]
Mance, D. Development of Electronic System for Sensing and Actuation of Test Mass of the Inertial Sensor LISA. Master’s Thesis, University of Split, Split, Croatia, 2012. [Google Scholar]
Tombolato, D. A Laboratory Study of Force Disturbances for the LISA Pathfinder Free Fall Demonstration Mission. Ph.D. Thesis, University of Trento, Trento, Italy, 2008. [Google Scholar]
Ke, G. Research on a Rapid Acquisition System for Electrical Impedance Imaging. Master’s Thesis, Chongqing University, Chongqing, China, 2009. [Google Scholar]
Chai, M. Research and Application of Capacitance Coupled Resistance Tomography System Based on Digital Phase Sensitive Demodulation Technology. Master’s Thesis, Zhejiang University, Hangzhou, China, 2014. [Google Scholar]

Figure 1. (a) Sensitive structure (including TM and EH), (b) 18 electrodes around the TM.

Figure 2. Actuation of TM in

φ

degree of freedom.

Figure 2. Actuation of TM in

φ

degree of freedom.

Figure 3. Rotational model of the torsion pendulum based on force modulation method.

Figure 4. Time domain waveform of angle signal.

Figure 5. Amplitude spectrum and phase spectrum of angle signal.

Figure 6. Principle diagram of PSD.

Figure 7. Charge measurement model of torsion pendulum based on PSD.

Figure 8. Simulation results of PSD, Butterworth BPF, Chebyshev II BPF and Bessel BPF methods for the ideal case (a) when

Q_{1} = 1 \times 10^{- 12} C

and (b) when

Q_{1} = 1 \times 10^{- 13} C

.

Figure 8. Simulation results of PSD, Butterworth BPF, Chebyshev II BPF and Bessel BPF methods for the ideal case (a) when

Q_{1} = 1 \times 10^{- 12} C

and (b) when

Q_{1} = 1 \times 10^{- 13} C

.

Figure 9. (a) Simulation results and corresponding absolute errors for multiple sets of simulations with PSD, Butterworth BPF, Chebyshev II BPF and Bessel BPF, in the ideal case, respectively. (b) Relative errors of PSD, Butterworth BPF, Chebyshev II BPF and Bessel BPF, for different initial set charge values in the ideal case, respectively.

Figure 10. Charge noise.

Figure 11. Angle sensor noise.

Figure 12. Simulation results of PSD, Butterworth BPF, Chebyshev II BPF and Bessel BPF methods with noise introduced (a) when

Q_{1} = 1 \times 10^{- 12} C

and (b) when

Q_{1} = 1 \times 10^{- 13} C

.

Figure 12. Simulation results of PSD, Butterworth BPF, Chebyshev II BPF and Bessel BPF methods with noise introduced (a) when

Q_{1} = 1 \times 10^{- 12} C

and (b) when

Q_{1} = 1 \times 10^{- 13} C

.

Figure 13. (a) Simulation results and corresponding absolute errors for multiple sets of simulations with PSD, Butterworth BPF, Chebyshev II BPF, and Bessel BPF with noise introduced, respectively. (b) Relative errors of PSD, Butterworth BPF, Chebyshev II BPF, and Bessel BPF for different initial set charge values with noise introduced, respectively.

Table 1. Parameters in the torsion pendulum rotation model based on the force modulation method.

Parameter	Meaning	Value
m	the mass of TM	$1.96 kg$
$C_{T}$	the total capacitance inside the EH	$34.2 \times 10^{- 12} F$
$ε$	the vacuum dielectric constant	$8.85 \times 10^{- 12} F / m$
$d_{0}$	the distance between the TM and the electrodes when the TM is at the center of the EH	$4 \times 10^{- 3} m$
b	the distance between the centers of the two electrodes on the same side	$10.75 \times 10^{- 3} m$
$a_{3}$	factor	$7.7597 \times 10^{- 13}$
$I_{r}$	the rotational inertia of the torsion pendulum	$4.31 \times 10^{- 5} {kgm}^{2}$
$γ$	the damping coefficient	$1.6557 \times 10^{- 10}$
Q	the quality factor of the torsion pendulum	2900
$Γ$	the restoring stiffness of the suspension fiber	$5.3491 \times 10^{- 9} Nm / rad$
$T_{0}$	free period of the torsion pendulum	$564 s$
$ω_{n}$	the intrinsic frequency of the torsion pendulum	$11.12 \times 10^{- 3} Hz$
$u_{T S, 1 φ}$	the absolute value of amplitude of the excitation applied to the electrodes	$3 V$
$f_{φ}$	the frequency of the excitation applied to the electrodes	$3 \times 10^{- 3} Hz$
$ξ$	the damping ratio	$1.7241 \times 10^{- 4}$
$Ω$	frequency normalisation	$1.6949$

Table 2. Comparison of simulation results between PSD, Butterworth BPF, Chebyshev II BPF and Bessel BPF methods in the ideal case and with noise introduced.

		PSD	Butterworth BPF	ChebyshevII BPF	Bessel BPF
in the ideal case	$R M S E$	$1.256 \times 10^{- 13} C$	$4.522 \times 10^{- 13} C$	$6.533 \times 10^{- 13} C$	$2.261 \times 10^{- 13} C$
	$M A E$	$3.472 \times 10^{- 14} C$	$1.250 \times 10^{- 13} C$	$1.806 \times 10^{- 13} C$	$6.25 \times 10^{- 14} C$
	$M A P E$	$0.00794$	$0.02015$	$0.02819$	$0.01107$
with noise introduced	$R M S E$	$2.051 \times 10^{- 13} C$	$9.026 \times 10^{- 13} C$	$1.108 \times 10^{- 12} C$	$3.693 \times 10^{- 13} C$
	$M A E$	$9.310 \times 10^{- 14} C$	$4.081 \times 10^{- 13} C$	$5.004 \times 10^{- 13} C$	$1.674 \times 10^{- 13} C$
	$M A P E$	$0.05575$	$0.10725$	$0.09556$	$0.11406$

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2024 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 (https://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Liu, Y.; Yu, T.; Wang, Y.; Zhao, Z.; Wang, Z. High-Precision Inertial Sensor Charge Ground Measurement Method Based on Phase-Sensitive Demodulation. Sensors 2024, 24, 1009. https://doi.org/10.3390/s24031009

AMA Style

Liu Y, Yu T, Wang Y, Zhao Z, Wang Z. High-Precision Inertial Sensor Charge Ground Measurement Method Based on Phase-Sensitive Demodulation. Sensors. 2024; 24(3):1009. https://doi.org/10.3390/s24031009

Chicago/Turabian Style

Liu, Yang, Tao Yu, Yuhua Wang, Zihan Zhao, and Zhi Wang. 2024. "High-Precision Inertial Sensor Charge Ground Measurement Method Based on Phase-Sensitive Demodulation" Sensors 24, no. 3: 1009. https://doi.org/10.3390/s24031009

APA Style

Liu, Y., Yu, T., Wang, Y., Zhao, Z., & Wang, Z. (2024). High-Precision Inertial Sensor Charge Ground Measurement Method Based on Phase-Sensitive Demodulation. Sensors, 24(3), 1009. https://doi.org/10.3390/s24031009

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

Article Menu

High-Precision Inertial Sensor Charge Ground Measurement Method Based on Phase-Sensitive Demodulation

Abstract

1. Introduction

2. Rotational Model of the Torsion Pendulum Based on the Force Modulation Method

2.1. TM Torque Driver Model

2.2. Dynamic Model and Frequency Response Analysis of the Torsion Pendulum

2.3. Rotational Model of the Torsion Pendulum Based on the Force Modulation Method

3. Charge Calculation Method Based on PSD

3.1. Design of Charge Calculation Method Based on PSD

3.2. Model of Charge Measurement in the Torsion Pendulum Based on PSD

3.3. Performance Verification

3.3.1. In the Ideal Case

3.3.2. With Charge Noise and Angle Sensor Noise Introduced

4. Conclusions

Author Contributions

Funding

Institutional Review Board Statement

Informed Consent Statement

Data Availability Statement

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI