Mobile QR Code QR CODE : Journal of the Korean Institute of Illuminating and Electrical Installation Engineers

Journal of the Korean Institute of Illuminating and Electrical Installation Engineers

ISO Journal TitleJ Korean Inst. IIIum. Electr. Install. Eng.

  1. (CEO, CLTech. Co., Ltd., Korea.)



Angular frequency, Damping coefficient, Effective impedance, Gas discharge, Laplace transform, Voltage oscillation wave

1. Introduction

When one device is electrically connected with other devices, if the impedance of the device is known, it is possible to calculate a physical quantity.

Solids with stable physical properties exhibit little change in impedance over time, whereas gases change their impedance due to free molecular motion. The impedance of gas discharge has both resistance and capacitance considerations, so it is very difficult to directly measure these physical quantities. Impedance was measured using the relationship between current and voltage(1,3) or current and power.(3) This paper proposes a method to calculate the effective impedance by applying Laplace transforms of the only the voltage wave form of gas discharge. The impedance is calculated from the attenuation coefficient and the the local oscillation frequency of the voltage wave form. This method can be applied when it is possible to measure the voltage waveform that appears in the case of under-damping.

2. Impedance calculation

2.1 Gas discharge circuit model

The gas discharge lamp circuit model used in this paper is shown in Fig. 1, and the test voltage waveform is shown in Fig. 2.

2.2 Model circuit analysis using Laplace transform

In Fig. 3, the voltage can be assumed since $i_{L}(0)=0$, $i_{R}(0)=0$ and $i_{C}(0)=0$. Therefore, it has an the initial condition of $v_{2}(0)=0$, $v_{2}'(0)=0$.

Equation(1) is established by using the Laplace transform of Fig. 3 considering the initial valuesof $v_{2}$.

(1)
$V_{2}(s)=\dfrac{R}{s^{2}RLC+s L+R}V(s)$

Equation(1) is transformed into equation(2).

Fig. 1. Gas discharge lamp circuit
../../Resources/kiiee/JIEIE.2021.35.7.001/fig1.png

Fig. 2. Voltage source wave form
../../Resources/kiiee/JIEIE.2021.35.7.001/fig2.png

Fig. 3. Voltage and current of model circuit
../../Resources/kiiee/JIEIE.2021.35.7.001/fig3.png

(2)
$\left(s^{2}+s\dfrac{1}{RC}+\dfrac{1}{LC}\right)V_{2}(s)=\dfrac{1}{LC}V(s)$

The voltage source waveform in Fig. 2 is expressed as equation(3).(4)

(3)
$v(t)=V_{0}[u(t)-2u\left(t-\dfrac{T}{2}\right)+2u(t-T)$ $-2u\left(t-\dfrac{3T}{2}\right)+\cdots]$

The Laplace transform of one period of the voltage source waveform is obtained by equation(4).

(4)
$V_{s1}(s)=\dfrac{V_{0}}{s}\left(1-e^{-\dfrac{T}{2}s}\right)^{2}$

Using equation(4), the Laplace transform of the voltage source waveform is represented by equation(5).

(5)
$V(s)=\dfrac{1}{1-e^{-Ts}}F_{1}(s)=\dfrac{V_{0}}{s}\dfrac{1-e^{-\dfrac{T}{2}s}}{1+e^{-\dfrac{T}{2}s}}$

By substituting equation(5) into equation(2), it is summarized as equation.(6)

(6)
$V_{2}(s)=\dfrac{1}{LC}\dfrac{V_{0}}{s\left(s^{2}+s\dfrac{1}{RC}+\dfrac{1}{LC}\right)}\bullet\dfrac{1-2e^{-\dfrac{T}{2}s}+e^{-Ts}}{1-e^{-Ts}}$

(7)
$s^{2}+s\dfrac{1}{RC}+\dfrac{1}{LC}=0$

In order for $v_{2}$ to oscillate, equation(7) must have an imaginary solution, so the condition of equation(8) must be met.(5)

(8)
$D=\left(\dfrac{1}{RC}\right)^{2}-\dfrac{4}{LC}< 0$

The solution of equation(7) is represented by equation.(9)

(9)
$s=-\dfrac{1}{2RC}+-\dfrac{\sqrt{D}}{2}$

In the case of $s=\alpha +- j\beta$ ($\alpha$ and $\beta$ are real numbers, $\beta\ne 0$), $\alpha$ is calculated by equation(10), and $\beta$ is calculated by equation(11).

(10)
$\alpha =-\dfrac{1}{2RC}$

(11)
$\beta =\sqrt{\dfrac{1}{LC}-\left(\dfrac{1}{2RC}\right)^{2}}=\sqrt{\dfrac{1}{LC}-\alpha^{2}}$

Using $\alpha$ and $\beta$, equation(11) becomes equation(12).

(12)
$V_{2}(s)=\dfrac{V_{0}}{LC}\left\{\dfrac{A}{s}+\dfrac{Bs+C}{(s-\alpha)^{2}+\beta^{2}}\right\}$ $=\dfrac{1-2e^{-\dfrac{T}{2}s}+e^{-Ts}}{1-e^{-Ts}}$

Equation(13) is used to analyze the damped oscillating waveform.

(13)
$H_{1}(s)=\dfrac{A}{s}+\dfrac{Bs+C}{(s-\alpha)^{2}+\beta^{2}}$

The inverse Laplace transform of equation(13) becomes equation.(14)

(14)
$h_{1}(t)=L^{-1}\left[H_{1}(s)\right]$ $=A+\dfrac{1}{\beta}e^{-\alpha t}\left(T_{a}\cos\beta t+S_{a}\sin\beta t\right)$

Where, $S_{a}=\dfrac{\alpha}{\alpha^{2}+\beta^{2}}$, $T_{a}=-\dfrac{\beta}{\alpha^{2}+\beta^{2}}$, $A=s H_{1}(s)|_{s=0}=LC$.

Equation(14) is transformed into equation(15).

(15)
$h_{1}(t)=L^{-1}\left[H_{1}(s)\right]$ $=LC+\dfrac{1}{\beta}e^{-\alpha t}\left(-\dfrac{\beta}{\alpha^{2}+\beta^{2}}\cos\beta t+\dfrac{\alpha}{\alpha^{2}+\beta^{2}}\sin\beta t\right)$

Substituting equation(15) into equation(12) to obtain equation.(16)

(16)
$V_{2}(s)=\dfrac{V_{0}}{LC}H_{1}(s)\left(1-2e^{-\dfrac{T}{2}s}+e^{-Ts}\right)\dfrac{1}{1-e^{-Ts}}$

In equation(16), the Laplace transform of the first half-period waveform is represented by equation.(17)

(17)
$V_{2half}(s)=\dfrac{V_{0}}{LC}H_{1}(s)$

Equation(18) can be obtained by inverse Laplace transform of equation.(17)

(18)
$v_{2half}(t)=L^{-1}\left[\dfrac{V_{0}}{LC}H_{1}(s)\right]$ $=V_{0}+\dfrac{V_{0}}{\beta LC}e^{-\alpha t}\left(-\dfrac{\beta}{\alpha^{2}+\beta^{2}}\cos\beta t+\dfrac{\alpha}{\alpha^{2}+\beta^{2}}\sin\beta t\right)$

In equation(18), $\alpha$ is the damping coefficient and $\beta$ is the angular frequency of the oscillation wave.

2.3 Calculationofimpedanceof discharge lamp

The impedance of the discharge lamp is calculated in the following order:

① voltage waveform measurement

② area calculation

③ calculate DC average value

④ calculate the AC value (waveform value minus DC average value)

⑤ absolute value calculation of AC value

⑥ trend line calculation

⑦ damping rate calculation

⑧ average period calculation of oscillation wave

⑨ calculation of effective resistance and effective capacitance

In this paper, a voltage v is 60kHz, 310V square wave, and inductance L is 20μH is applied for the circuit presented in Fig. 1.

Fig. 4 is a measurement of the current and voltage wave forms of the discharge lamp in the Fig. 1 circuit.

Fig. 4. Discharge lamp current and voltage wave forms
../../Resources/kiiee/JIEIE.2021.35.7.001/fig4.png

In Fig. 4, the time and voltage of the peaks and valleys are measured on the half cycle of the voltage waveform. The area of ​​the graph is obtained by using the trapezoidal formula, and the area is divided by time to obtain the average value of voltage.

Fig. 5. Measured crest, valley and average value of voltage wave form
../../Resources/kiiee/JIEIE.2021.35.7.001/fig5.png

Fig. 6. Measured value minus average value of voltage wave form
../../Resources/kiiee/JIEIE.2021.35.7.001/fig6.png

Fig. 7. Trend line of selected data
../../Resources/kiiee/JIEIE.2021.35.7.001/fig7.png

Fig. 8. LTspice simulation circuit
../../Resources/kiiee/JIEIE.2021.35.7.001/fig8.png

Fig. 9. Simulated voltage and current wave forms
../../Resources/kiiee/JIEIE.2021.35.7.001/fig9.png

Fig. 5 shows the peaks and valleys of the measured voltage waveform, and the average voltage value, 310V, of the measured voltage waveform. The calculated average value of the voltage should be the same as the voltage of the applied square wave as shown in equation.(18)

Fig. 6 shows the waveform with the average value (DC component) removed from the measured voltage waveform.

The waveform in Fig. 6 shows an exponential multiply sine wave function form ($v=A e^{\alpha t}\sin\omega t$). In order to obtain the attenuation coefficient α, the absolute value of the waveform in Fig. 6 is taken, and the trend line is calculated by using the first three values. Fig. 7 shows this process.(6)

It is $\alpha =-5.365\times 10^{5}$ calculated from the trend line.

In Fig. 6 the average period is calculated by using the time from the first peak to the second valley. The calculated average period is 1.37μs. From the above relation, it is calculated as $\beta =4.59\times 10^{6}$ by equation.(19)

(19)
$\beta =\omega =\dfrac{2\pi}{T_{avg}}$

From equation(10) and equation(11), the circuit capacitance C is calculated as equation(20), and the circuit resistance R is calculated as equation(21).

(20)
$C=\dfrac{1}{L\left(\alpha^{2}+\beta^{2}\right)}$

(21)
$R=\dfrac{1}{2C\alpha}$

In the circuit used in this paper, the circuit inductance L is 20μH. Therefore the circuit capacitance C is calculated as 2.35nF and the circuit resistance R is calculated as 397Ω.

In order to verify the calculated resistance and capacitance values ​​of the discharge lamp, Fig. 1(b) circuit was simulated by using LTspice.

Fig. 9 shows the simulation result of Figure 8 circuit. It can be seen that the simulation waveform has a similar trend to the actual waveform in Fig. 4.

Fig. 10. Measured and simulated voltage
../../Resources/kiiee/JIEIE.2021.35.7.001/fig10.png

When the measured voltage waveform and the simulated voltage waveform are compared quantitatively, a difference between the two can be seen.

3. Conclusion

In this paper, a gas discharge is modeled as a parallel circuit with resistance and capacitance. In order to obtain the effective resistance and effective capacitance of the gas discharge, the circuit equation is theoretically analyzed based on Laplace transforms. The derived results were applied to an actual discharge lamp system to obtain the effective resistance and effective capacitance.

In order to verify the calculated resistance and capacitance values ​​of the discharge lamp, a simulation was performed by using LTspice, and the results were compared with the measured waveform. When voltage and current wave forms were compared, the macroscopic behavior was similar, but there was a difference in quantitative terms. It is because this paper is based on the assumption that effective resistance and capacitance are constant during discharge, whereas they change over time during actual discharge.

It is necessary to research the dynamic impedance of gas discharge lamps in the future.

References

1 
Bakker L P., Gerrit M., Kroesen W., Frederik J. de H., 1999, RF Discharge Impedance Measurements Using a New Method to Determine the Stray Impedances, IEEE TRANSACTIONS ON PLASMA SCIENCE, Vol. 27, No. 3, pp. 759-765DOI
2 
Overzet L. J., 2010, RF Impedance Measurements of DC Atmospheric Micro-discharges, The European Physical Journal D, Vol. 60, pp. 449-454DOI
3 
Spiliopoulos N., Mataras D., Rapakoulias D. E., 1996, Power Dissipation and Impedance Measurements in Radiofrequency Discharges, Journal of Vacuum Science & Technology A 14, pp. 2757-2765DOI
4 
Lee D. C., Lee N. H., Go U. S., Lee D. I., 1999, Circuit Analysis, Dongilchulpansa, pp. 534-538Google Search
5 
Lee D C., Lee N. H., Go U. S., Lee D. I., 1999, Circuit Analysis, Dongilchulpansa, pp. 558-559Google Search
6 
Chee C. K., 2002, Latest Electrical and Electronic Measurement and Basic Experiments, Moonwoondang, pp. 18-20Google Search

Biography

Chin-Woo Yi
../../Resources/kiiee/JIEIE.2021.35.7.001/au1.png

He received a Ph.D degree from Seoul National University, Korea, in 1990.

He is now a professor at Hoseo University, Korea.

Hyun-Bae Choi
../../Resources/kiiee/JIEIE.2021.35.7.001/au2.png

He is a Ph.D candidate at Hoseo University, Korea.

He is now the CEO of CLTech. Co., Ltd., Korea.