| Citation: |
Wenting WENG, Donghui XIA, Yizhe TIAN, Xixuan CHEN, Zhijiang WANG, Yuan PAN. Investigation on the multi-hole directional coupler for power measurement of the J-TEXT ECRH system[J]. Plasma Science and Technology, 2022, 24(6): 065602. DOI: 10.1088/2058-6272/ac5d7b
|
Power measurement is necessary for an electron cyclotron resonance heating (ECRH) system. The directional coupler method has been put forward to monitor high-power microwave from gyrotrons in real time. A multi-hole directional coupler has been designed and manufactured for the 105 GHz/500 kW ECRH system on the J-TEXT tokamak. During the design process, we established the relationships between hole parameters and coupling characteristics based on the multi-hole coupling method and small-hole coupling theory. High-power tests have been carried out. The results indicated the reasonability of the theoretical design and practicality of the fabricated directional coupler. Sources of test errors have been discussed in detail, and the influences of spurious modes on the directional couplers have been emphatically analyzed.
In order to transfer power energy in large capacity and long-distance, several HVDC lines have been completed. Corona discharge occurs during the operation of the HVDC lines, which leads to a problem of ion flow field [1, 2]. The surrounding environment of the HVDC line is complex, and the temperature and humidity will affect the ion flow field of HVDC lines [3]. Especially in high humidity, suspension droplets will appear in the space and capture the space charge generated by corona discharge. Then the charged suspension droplets will be formed. Charged suspension droplets will cause distortion of the surrounding electric field and complicate the ion flow field. The parameters such as radius and number density of the suspension droplets in space cannot be measured accurately, which makes the charge of suspension droplets difficult to calculate in the ion flow field. Furthermore, the effect of temperature and humidity on the ion flow field is not clear.
Some scholars have studied the charging characteristics of suspension particles in the ion flow field. Dastoori et al used the Faraday cup to measure the charge of dust particles on the bottom of PCB (printed circuit board) [4]. The charge of particles in a fluidized bed was measured by Song D with a Faraday tube, and the effect of charged particles on the fluidized bed was analyzed [5]. Xu et al measured the concentration and size distribution of particles in flue gas using ELPI (electrical low-pressure impactor), and the average charge of a single-particle was calculated [6].
In the aspect of charge of small particles, Hewitt used the coaxial cylinder mobility analyzer to measure the charge of particles with a diameter of 0.07–0.66 μm. The measured charge was consistent with the charging model considering field charging and diffusion charging, and the correctness of the charging model was verified [7]. Luo et al measured the average charge of a single particle under different voltages in DC electric field by using ELPI. The results show that the charging modes of particles with diameters smaller and larger than 0.2 μm are diffusion charging and field charging, respectively [8]. The above literature mainly studied the charging characteristics and models of solid particles but did not involve the charging characteristics of suspension droplets in space.
As for the treatment of suspension droplets in the ion flow field, Ma et al measured the ion flow field of HVDC wire under different temperatures and humidity, and the variation of corona voltage with temperature and humidity is considered to be the main reason for the change of ion flow field, so the influence of suspension droplets was not considered [9]. Li et al considered that the effect of charged suspension droplets on ion flow field should be considered under any humidity condition, and charged suspension droplets were randomly distributed in space [10]. Up to now, the charging characteristics of suspension droplets in the ion flow field under different temperatures and humidity are not clear. It is necessary to study the characterization method of the charging characteristics of suspension droplets in the ion flow field to provide the theoretical basis for calculating the ion flow field of HVDC lines under different temperatures and humidity.
In this paper, firstly, the effective charging factor to characterize the charging characteristics of suspension droplets is introduced. Secondly, the charging factor of suspension droplets under different temperatures and humidity is calculated based on the measurement results of the ion flow field, and the analytic expression of the charging factor considering the influence of temperature and humidity is obtained. Finally, the influence of the charged suspension droplets in the ion flow field is analyzed. The analytic expression can be used to calculate the ion flow field under different temperatures and humidity, which provides a theoretical basis for the construction of HVDC lines.
The charging characteristics of suspension droplets are determined by radius, electrical permittivity, and number density. It is difficult to quantitatively analyze the charging characteristics of suspension droplets because the above parameters are difficult to be measured accurately. In order to calculate the space charge density of suspension droplets, the charging factor is defined. Because the field charging model can be used to calculate the charge of the suspension droplets [10], the expression of effective charging factor is:
| θ=ρwE | (1) |
where E is the electric field intensity, and ρw is the space charge density of charged suspension droplets.
When the charging factor is determined, the calculation model of the ion flow field considering the influence of suspension droplets is as follows:
| ∇2φ=-(ρe+θE)ε0∇·(ρek→E)=0 | (2) |
where φ is the electric potential,
The boundary condition of equation (2) is that the electric potential on the lines and grounding point are operating voltage and 0, respectively.
According to equation (2), when the charging factor is determined, the ion flow field under the current temperature and humidity can be calculated to realize the prediction of the ion flow field.
The analysis is based on the following assumptions in this work:
(1) The ion mobility and corona inception voltage of the conductor are fixed at a certain temperature, humidity, and air pressure.
(2) It can be considered that the suspension droplets are in a static state in the ion flow field [11].
(3) The electric field on the surfaces of the high voltage conductors is the corona inception electric field after the inception of corona discharge.
(4) The suspension droplets are uniformly distributed in space, and their shape is spherical.
It can be seen from equation (2) that the calculated value of electric field Ec is a function of charging factor θ when the ion mobility, corona inception voltage, and electrode structure are determined.
| Ec=f(θ) | (3) |
The error g(Ec) is defined:
| g(Ec)=|Em-Ec|Ec | (4) |
where Em is the measured value of an electric field.
When the error g(Ec) is less than a certain limit, the charging factor can be considered as the charging factor in the current environment. There is no analytical solution for equation (2) because of the introduction of charged suspension droplets. Therefore, the upstream finite element method (FEM) is used to solve equation (2).
According to the above analysis, the calculation process of the charging factor of suspension droplets is as follows:
Step 1. Give the initial values of the charging factor.
Step 2. Give the initial values of charge density.
Step 3. Calculate the space charge density of the suspension droplets and the total space charge density.
Step 4. Calculate the electric field using FEM to solve the Poisson equation.
Step 5. Calculate the space charge density by using the upstream element method to solve the current continuity equation.
Step 6. If the error limiting conditions of space charge density and the electric field in equations (5) and (6) are satisfied, g(Ec) is calculated. Otherwise, the surface charge density value is updated and steps 2–5 are repeated until the error limit condition is satisfied. The error limiting conditions are as follows:
| σρ=|ρn-ρn-1|ρn≤0.01% | (5) |
| σE=|Emax-E0|E0≤0.01% | (6) |
where σρ and σE are the space charge density error and the electric field error on the conductor surface, respectively, ρn and ρn-1 are the surface charge density of the conductor of the nth and (n-1)th iteration, respectively, Emax and E0 are the maximum electric field and corona inception electric field on the conductor surface, respectively.
In order to consider the influence of the suspension droplets in the calculation, the error limit of g(Ec) should be less than the proportion of the electric field component generated by the charged suspension droplets at the cage surface. Therefore, the error limit of g(Ec) in this paper is 0.01%, which can not only ensure the accuracy of the calculation results, but also ensure the calculation speed and convergence of the calculation.
The calculation flow shows that in order to solve the charging factor, the ion flow field needs to be measured, the corona voltage and ion mobility under different temperatures and humidity are determined.
The artificial climate chamber is used to adjust the temperature and humidity. High-precision temperature and humidity sensors are used to measure the temperature and humidity in the climate chamber. The parameters of the artificial climate chamber are shown in table 1.
| Parameter name | Parameter value |
| Range of temperature (℃) | 15–35 |
| Range of humidity (%) | 20–90 |
| Fluctuation of temperature (℃) | ±0.5 |
| Fluctuation of humidity (%) | ±3% |
DownLoad:
CSV
The coaxial cylindrical electrode structure is used to generate the corona discharge, as shown in figure 1. The central electrode is a bare wire with a length of 2.8 m and a radius of 1.1 mm. The grounded corona cage is composed of three parts that are electrically insulated from each other. The measuring section is used to measure the ion flow field of HVDC wire, and the shielding section is used to weaken the end-effect. The length of the measuring and shielding sections is 2 m and 0.3 m, respectively, and their diameter is 0.8 m.
The experimental platform was placed in the artificial climate chamber. Under the conditions of temperature T=30 ℃, 25 ℃, and 20 ℃, the relative humidity RH was adjusted in the range of 30%–90%. It is difficult to clearly reflect the dual impact of temperature and humidity on the measurement and calculation result using the absolute humidity related to both temperature and humidity, so the relative humidity is used in this work. Then, the positive voltage of 0–65 kV was applied to the central electrode by the DC voltage source. The field mill was used to measure the electric field strength, and the ion current plate was used to measure the ion current density.
When the applied voltage on the conductor is 50 kV and the corona discharge occurs, the measurement of the ion flow field is shown in figure 2 under different temperatures and humidity. The measurement shows that the total electric field and ion current density decrease with the relative humidity and increase with the temperature.
Many studies have shown that temperature and humidity have a great influence on corona inception voltage. The influence of temperature on corona inception voltage was consistent, and the corona voltage of the conductor decreases with the increase of temperature [12]. The influence of humidity on the corona discharge characteristics of wire has not been consistent [13–16]. It is necessary to analyze the corona discharge characteristics of wires under different humidity according to the wire structure.
The corresponding voltage at the inflection point of the relation curve between ion current density and conductor voltage is the corona inception voltage of the conductor [17]. The curve of corona inception voltage with temperature and humidity is shown in figure 3. The corona inception voltage increases with the increase of relative humidity and decreases with the increase of temperature.
Some scholars have studied ion mobility under different temperatures and humidity, the calculation models of ion mobility were established. In this work, the calculation model of ion mobility proposed by Zhang Bo et al was used [18]. The calculations of positive ion mobility under different temperatures and humidity are shown in figure 4. From figure 4, the positive ion mobility decreases with the increase of relative humidity and temperature.
In order to ensure the convergence and calculation speed of the calculation results of the charging factor, it is necessary to select a reasonable initial value of the charging factor. The initial value of the charging factor is:
| θ0=9λmε0ερ1r(ε+2ε0) | (7) |
where λ is the correction factor, ρ1 is the mass density of water, m is water content in the air, and the value of r is 0.2 μm.
In equation (7), the value range of λ is [0, 1) (When the relative humidity is 0%, λ is 0). The value of λ represents the ratio of the actual charging factor in the space to the possible maximum charging factor. The value of λ obtained through a large number of calculations is shown in figure 5.
From figure 5, it can be seen that the value of λ is small when the relative humidity is less than 60%, and increases with the relative humidity when the relative humidity is more than 60%. In this work, the average is taken as the value of λ. The value of λ is:
| λ={2.89×10-9RH<60%5.72×10-7RH≥60% | (8) |
When suspension droplets exist in space, charged suspension droplets will strengthen the ground electric field [10], so the iterative formula of charging factor used in this work is:
| θn+1=θn(1+μEm-EcEm+Ec) | (9) |
where θ n+1 and θ n are the charging factors of the (n+1)th and nth iteration, respectively; μ is the iteration coefficient in this paper, and its value is 1.
When T=25 ℃ and RH=70%, the charging factor of suspension droplets is calculated under different iteration coefficients. The calculation time and charging factor are shown in table 2.
| The iteration coefficient | Calculation time (s) | Charging factor (10-15 C/(V·m2)) | |
| 0.05 | 119.46 | 3.24 | |
| 0.1 | 78.03 | 3.24 | |
| 0.5 | 40.96 | 3.24 | |
| 1 | 34.54 | 3.24 | |
| 1.5 | Non convergence | ||
DownLoad:
CSV
As shown in table 2, when the iteration coefficient is less than 1, the calculation time is significantly increased, but the calculation accuracy is not significantly improved. When the iteration coefficient is more than 1, the calculation may not converge. Using the value of μ given in this paper can ensure the calculation speed and accuracy.
A linear function is used to fit the calculated results of the charging factor, and the analytic expression of the fitting curve is shown in figure 6.
The calculated results show that the charging factor is small and increases little with the relative humidity when the relative humidity is less than 60%. When the temperature is 30 ℃, the charging factor at RH=50% is 1.1 times higher than that at RH=40%. When the relative humidity is more than 60%, the charging factor is large and increases obviously with the increase of relative humidity. When the temperature is 30 ℃, the charging factor at RH=70% is 7.8 times higher than that at RH=60%.
The charging factor increases linearly with the increase of temperature. When the relative humidity is 50%, the charging factor at T=30 ℃ is 1.32 times higher than that at T=25 ℃, and the charging factor at T=25 ℃ is 1.33 times higher than that at T=20 ℃. The increase of saturated water content in space is the main reason for the increase of charging factor with temperature.
The intersection points of analytic expression are calculated. When T=30 ℃, 25℃ and 20 ℃, the abscissa of intersection points is about 59%. The results show that when the relative humidity is more than 60%, there will be charged suspension droplets in the ion flow field.
The slope and intercept of analytic expression change with temperature. A linear function is used to fit the slope and intercept, and the fitting curve is shown in figure 7. The curve-fitting results are:
| A={(1.243T-8.638)×10-5RH<60%0.01979T-0.1651RH≥60% | (10) |
| B={(0.911T-6.32)×10-3RH<60%1.218T-10.95RH≥60% | (11) |
where A and B are the slope and intercept of analytic expression, respectively.
In summary, the charging factor of suspension droplets can be obtained as follows:
| θ=A·RH+B | (12) |
In order to verify the validity of the analytic expression, the ion flow field was measured at a temperature is 15 ℃. The charging factor is calculated based on the measurement of the ion flow field and compared with the calculation results of the analytic expression. The comparison results are shown in figure 8.
The comparison results show that the maximum error is 6.4%, which indicates that the analytic expression can be used to calculate the charging factor and space charge density of suspension droplets and model the ion flow field considering the influence of temperature and humidity.
The charging factor is introduced into equation (2) to solve the space charge density under different temperatures and humidity. When the applied voltage on the conductor is 50 kV, the percentages of suspension droplets charge density (ρw) and electric field component (Ew) generated by charged suspension droplets at the cage surface are calculated, and the calculation results are shown in figures 9 and 10.
As seen in figures 9 and 10, ρw and Ew increase with the temperature and humidity. When the relative humidity is less than 60%, ρw and Ew are small and increase little with the relative humidity and temperature. When the relative humidity is more than 60%, ρw and Ew are large and increase obviously with the increase of relative humidity and temperature, which indicates that the electric field distortion caused by charged suspension droplets becomes more serious and charged suspension droplets have an effect on the characteristics of the ion flow field.
When the applied voltage on the conductor is 50 kV, the comparison results of the electric field at the cage surface with and without suspension droplets under different temperatures and humidity are shown in figure 11.
The results show that there is little difference between the electric field with and without suspension droplets at the same temperature and relative humidity of less than 60%. When T=25 ℃ and RH=40%, the difference is 0.5%. The results show that the ion flow field is less affected by suspension droplets, and the influence of suspension droplets can be ignored in engineering calculation.
In an environment of the same temperature and relative humidity of more than 60%, the electric field without suspension droplets is less than that with suspension droplets. The results show that the suspension droplets can enhance the total electric field on the ground. With the increase of relative humidity and temperature, the difference between the electric field with and without suspension droplets is larger. When the temperature is 25 ℃, the differences are 3.5% and 7.6% at RH=70% and RH=90%, respectively. When the relative humidity is 90%, the differences are 5.7% and 11.1% at T=20 ℃ and T=30 ℃, respectively. The above analysis shows that the suspension droplets greatly influence the ion flow field of the HVDC conductor under high temperatures and high humidity.
In this paper, the charging factor used to characterize the charge characteristics of suspension droplets is introduced. After, the calculation results of the charging factor are fitted by a linear function, and the analytic expression of the charging factor is obtained. Then, the influence of charged suspension droplets in the ion flow field is analyzed. The results are as follows:
(1) The charging factor is small and increases little with the relative humidity when the relative humidity is less than 60%. When the relative humidity is more than 60%, the charging factor of suspension droplets is larger and increases rapidly with the relative humidity, which indicates that there will be obvious charged suspension droplets in the ion flow field. The charging factor of suspension droplets increases linearly with the increase of temperature at the same relative humidity.
(2) When the relative humidity is less than 60%, the difference between the electric field at the cage surface with and without suspension droplets is small, so the influence of suspension droplets can be ignored in engineering calculation. When the relative humidity is more than 60%, the difference between the electric field at the cage surface with and without suspension droplets is 11.1% at RH=90% and T=30 ℃, so the influence of suspension droplets cannot be ignored in the calculation of ion flow field.
This work was supported by the National Key Research and Development Program of China (Nos. 2017YFE0300200 and 2017YFE0300204), and in part by National Natural Science Foundation of China (No. 51821005).
| [1] |
Thumm M K and Kasparek W 2002 IEEE Trans. Plasma Sci. 30 755 doi: 10.1109/TPS.2002.801653
|
| [2] |
Lau C et al 2019 IEEE Trans. Microw. Theory Tech. 67 38 doi: 10.1109/TMTT.2018.2879808
|
| [3] |
Zapevalov V E et al 2001 Fusion Eng. Des. 53 377 doi: 10.1016/S0920-3796(00)00513-5
|
| [4] |
Thumm M K A et al 2019 Nucl. Fusion 59 073001 doi: 10.1088/1741-4326/ab2005
|
| [5] |
Clarricoats P J B and Olver A D 1973 Electron. Lett. 9 376 doi: 10.1049/el:19730278
|
| [6] |
Kowalski E J et al 2010 IEEE Trans. Microw. Theory Tech. 58 2772 doi: 10.1109/TMTT.2010.2078972
|
| [7] |
Bruschi A et al 2003 Nucl. Fusion 43 1513 doi: 10.1088/0029-5515/43/11/024
|
| [8] |
Cengher M et al 2009 Fusion Sci. Technol. 55 213 doi: 10.13182/FST09-A4073
|
| [9] |
Michel G, Laqua H and Kasparek W 2001 Fusion Eng. Des. 56–57 655 doi: 10.1016/S0920-3796(01)00383-0
|
| [10] |
Ives R L et al 1999 IEEE Trans. Plasma Sci. 27 531 doi: 10.1109/27.772282
|
| [11] |
Notake T et al 2005 Rev. Sci. Instrum. 76 023504 doi: 10.1063/1.1848657
|
| [12] |
Kasparek W et al 2003 Int. J. Infrared Millimeter Waves 24 451 doi: 10.1023/A:1022468412389
|
| [13] |
Edgcombe C J 1993 Gyrotron Oscillators - Their Principles and Practice (London: Taylor & Francis))
|
| [14] |
Bethe H A 1944 Phys. Rev. 66 163 doi: 10.1103/PhysRev.66.163
|
| [15] |
Cohn S B 1952 Proc. IRE 40 696 doi: 10.1109/JRPROC.1952.274063
|
| [16] |
Mcdonald N A 1972 IEEE Trans. Microw. Theory Tech. 20 689 doi: 10.1109/TMTT.1972.1127844
|
| [17] |
Sporleder F and Unger H G 1979 Waveguide Tapers Transitions and Couplers (London: Peter Peregrinus)
|
| [18] |
Wang W X 1994 Int. J. Electron. 76 131 doi: 10.1080/00207219408925914
|
| [19] |
Wang W X et al 2008 Int. J. Electron. 95 99 doi: 10.1080/00207210701814119
|
| [20] |
Xia D H et al 2014 Plasma Sci. Technol. 16 267 doi: 10.1088/1009-0630/16/3/17
|
| [21] |
Ramo S, Whinnery J R and van Duzer T 1994 Fields and Waves in Communication Electronics (New York: Wiley)
|
| [22] |
Zach A et al 2019 EPJ Web Conf. 203 04018 doi: 10.1051/epjconf/201920304018
|
| [23] |
Liu C H et al 2018 IEEE Trans. Antennas Propag. 66 6043 doi: 10.1109/TAP.2018.2862467
|
| [24] |
Tian Y Z et al 2019 Fusion Eng. Des. 146 269 doi: 10.1016/j.fusengdes.2018.12.043
|
| [25] |
Doane J L et al 2016 Fusion Eng. Des. 102 99 doi: 10.1016/j.fusengdes.2015.11.017
|
| [26] |
Wagner D et al 2017 J Infrared Milli. Terahertz Waves 38 191 doi: 10.1007/s10762-016-0329-y
|
| [27] |
Ruiz J et al 2012 J Infrared Milli. Terahertz Waves 33 491 doi: 10.1007/s10762-012-9883-0
|
| [1] | Jiali CHEN (陈佳丽), Peiyu JI (季佩宇), Chenggang JIN (金成刚), Lanjian ZHUGE (诸葛兰剑), Xuemei WU (吴雪梅). The properties of N-doped diamond-like carbon films prepared by helicon wave plasma chemical vapor deposition[J]. Plasma Science and Technology, 2019, 21(2): 25502-025502. DOI: 10.1088/2058-6272/aaee90 |
| [2] | Gao ZHAO (赵高), Wanying ZHU (朱婉莹), Huihui WANG (王慧慧), Qiang CHEN (陈强), Chang TAN (谭畅), Jiting OUYANG (欧阳吉庭). Study of axial double layer in helicon plasma by optical emission spectroscopy and simple probe[J]. Plasma Science and Technology, 2018, 20(7): 75402-075402. DOI: 10.1088/2058-6272/aab4f1 |
| [3] | Ming SUN (孙明), Zhan TAO (陶瞻), Zhipeng ZHU (朱志鹏), Dong WANG (王东), Wenjun PAN (潘文军). Spectroscopic diagnosis of plasma in atmospheric pressure negative pulsed gas-liquid discharge with nozzle-cylinder electrode[J]. Plasma Science and Technology, 2018, 20(5): 54005-054005. DOI: 10.1088/2058-6272/aab601 |
| [4] | LAN Hui (兰慧), WANG Xinbing (王新兵), ZUO Duluo (左都罗). Time-Resolved Optical Emission Spectroscopy Diagnosis of CO2 Laser-Produced SnO2 Plasma[J]. Plasma Science and Technology, 2016, 18(9): 902-906. DOI: 10.1088/1009-0630/18/9/05 |
| [5] | WU Ding (吴鼎), LIU Ping (刘平), SUN Liying (孙立影), HAI Ran (海然), DING Hongbin (丁洪斌). Influence of a Static Magnetic Field on Laser Induced Tungsten Plasma in Air[J]. Plasma Science and Technology, 2016, 18(4): 364-369. DOI: 10.1088/1009-0630/18/4/06 |
| [6] | LI Cong (李聪), ZHANG Jialiang (张家良), YAO Zhi (姚志), WU Xingwei (吴兴伟), et al.. Diagnosis of Electron, Vibrational and Rotational Temperatures in an Ar/N 2 Shock Plasma Jet Produced by a Low Pressure DC Cascade Arc Discharge[J]. Plasma Science and Technology, 2013, 15(9): 875-880. DOI: 10.1088/1009-0630/15/9/08 |
| [7] | LU Wenqi (陆文琪), JIANG Xiangzhan (蒋相站), LIU Yongxin (刘永新), YANG Shuo (杨烁), et al. Improved Double-Probe Technique for Spatially Resolved Diagnosis of Dual-Frequency Capacitive Plasmas[J]. Plasma Science and Technology, 2013, 15(6): 511-515. DOI: 10.1088/1009-0630/15/6/05 |
| [8] | Miyuki YAJIMA, Masato YAMAGIWA, Shin KAJITA, Noriyasu OHNO, Masayuki TOKITANI, Arimichi TAKAYAM, Seiki SAITO, Atsushi M. ITO, Hiroaki NAKAMURA, Naoaki YOSHIDA. Comparison of Damages on Tungsten Surface Exposed to Noble Gas Plasmas[J]. Plasma Science and Technology, 2013, 15(3): 282-286. DOI: 10.1088/1009-0630/15/3/18 |
| [9] | CHEN Zhaoquan (陈兆权), LIU Minghai (刘明海), HU Yelin (胡业林), ZHENG Xiaoliang (郑晓亮), LI Ping (李平), XIA Guangqing (夏广庆). Character Diagnosis for Surface-Wave Plasmas Excited by Surface Plasmon Polaritons[J]. Plasma Science and Technology, 2012, 14(8): 754-758. DOI: 10.1088/1009-0630/14/8/13 |
| [10] | PANG Jianhua (庞见华), LU Wenqi (陆文琪), XIN Yu (辛煜), WANG Hanghang (王行行), HE Jia (贺佳), XU Jun (徐军). Plasma Diagnosis for Microwave ECR Plasma Enhanced Sputtering Deposition of DLC Films[J]. Plasma Science and Technology, 2012, 14(2): 172-176. DOI: 10.1088/1009-0630/14/2/17 |
Figures(14)
Supported by: Beijing Renhe Information Technology Co., Ltd. 
| Parameter name | Parameter value |
| Range of temperature (℃) | 15–35 |
| Range of humidity (%) | 20–90 |
| Fluctuation of temperature (℃) | ±0.5 |
| Fluctuation of humidity (%) | ±3% |
DownLoad:
CSV
| The iteration coefficient | Calculation time (s) | Charging factor (10-15 C/(V·m2)) | |
| 0.05 | 119.46 | 3.24 | |
| 0.1 | 78.03 | 3.24 | |
| 0.5 | 40.96 | 3.24 | |
| 1 | 34.54 | 3.24 | |
| 1.5 | Non convergence | ||
DownLoad:
CSV
