Development of A Method of Calculating The Melting Characteristics of OPGW Strands Due To DC Arc Simulating Lightning Strike
Development of A Method of Calculating The Melting Characteristics of OPGW Strands Due To DC Arc Simulating Lightning Strike
3, JULY 2013
Abstract—Some strands of composite fiber-optic ground wire lightning impulse current generator cannot supply the current
(OPGW) are sometimes melted and broken by high-energy light- whose duration is so long and the electric charge is so large. In
ning strikes. DC arc tests simulating lightning strikes have been order to clarify the aforementioned characteristics of OPGW
performed to obtain the melting and breaking characteristics
of OPGW strands. The tests have to be performed under many strands, consequently, dc arc tests have been performed and
conditions concerning the arc (e.g., current, duration, polarity, gap test conditions that adequately simulate high-energy lightning
length) and the OPGW (e.g., size, type, and number of strands) to strikes have also been investigated [1], [2], [6], [7]. However,
clarify the melting and breaking characteristics of the strands. In lightning and OPGW conditions vary considerably. The light-
this paper, the calculations regarding the melting characteristics ning conditions are wide ranging as described before. On the
of strands are performed considering the transferred heat and
its area from the arc to the strands under the aforementioned other hand, there are many combinations of the types and num-
conditions. The melting characteristics of strands are calculated bers of strands of OPGW because they are selected according
with an arc current of 1–100 kA considering the measured current to the conditions of the transmission lines (e.g., transmission
of actual lightning. The calculation results of the strand melting voltages and single-line-to-ground fault currents). Figs. 1(a)
duration depending on arc current show good agreement with the and (b) shows the structure of OPGW with cross sections of
experimental values obtained in dc arc tests.
80 and 140 mm , respectively. These are one of the OPGWs
Index Terms—Arc discharges, ground wires, lightning, optical- generally used. The types of strands in Fig. 1(a) and (b) are
fiber cables, power transmission lines.
23AC and 40AC, respectively. “AC” means “Aluminum-Clad
steel”, while “23” and “40” denote that the electric conduc-
I. INTRODUCTION tivity of the AC strand is “23%” and “40%” of the electric
conductivity of copper strand with equivalent cross section to
the AC strand, respectively. If dc arc tests are performed under
forming a new arc root where the arc column is brought into
contact with with the OPGW, insulating tape is wound onto both
sides of the arc-producing (exposed) part of length , about
20–40 mm.
Fig. 1. Examples of the OPGW structure. (a) 80 mm (23AC 6/(3.86)) [2]. III. CALCULATION METHODS AND CONDITIONS
(b) 140 mm (40AC 12/3.0, 40AC 8/(2.86)) [1].
Fig. 3 shows an outline of the flowchart used to calculate
the melting characteristics of OPGW due to the DC arc. By
using the arc conditions (e.g., current, duration, polarity, gap
length) and OPGW conditions (e.g., the size, type and number of
strands), the heat quantity and its area transferred to the OPGW
from the arc are calculated, whereupon the melting quantity of
OPGW due to the arc is calculated. In the calculation of this
paper, the arc current varies from 1–100 kA considering the
measured current of actual lightning, and the objects for the
calculation are the two types of OPGW shown in Fig. 1. This
chapter describes the calculation methods and conditions.
(1)
Fig. 2. Arrangement of dc arc testing.
where is the heat due to electrons emitted from or absorbed
into the OPGW, is the convective heat transfer due to
the arc jet spouted from the counter electrode, is the heat
An arc root on the OPGW tends to be too mobile to melt the
due to the radiation from the arc, and is the joule heating in
strands, due to the electromagnetic force (emf) induced by the
the electrode. As and are relatively smaller than
current flowing through the OPGW. Accordingly, conductors
[8], this small degree of heat is not considered in this paper.
are connected to both ends of the OPGW to divert the arc current
Therefore, is estimated as the following equation:
to the earth side symmetrically so that the arc root cannot easily
move. The current is supplied by a full-wave rectifying three- (2)
phase ac using silicon diodes. The ac is generated using a short-
circuit generator (15 kV, 2500 MVA) and a transformer. The dc 1) Heat Due to Electrons : changes depending on
output can be changed between 120 kV/10 kA to 16 kV/60 kA the electrode polarity are as follows [9]
by changing the ac voltage and combination of diodes. In the dc • In the case of a cathode
arc tests, the supply dc voltage is set at about 10–17 kV so that
the supply voltage is not significantly affected by an increase
in arc voltage. Furthermore, to prevent the arc from moving by (3)
1316 IEEE TRANSACTIONS ON POWER DELIVERY, VOL. 28, NO. 3, JULY 2013
TABLE I
VOLTAGES CONCERNING THE HEAT DUE TO ELECTRONS (SURROUNDING GAS: AIR)
• In the case of an anode (933 K for aluminum and 1809 K for iron [17]) during testing.
for in the room temperature is larger by about
10%–20% than that for the melting point. In the calculation,
(4) is assumed to be 300 K in order to determine the max-
where and are the cathode-fall and the anode-fall volt- imum value of .
ages, respectively; is the work function of the electrode; and The heat-transfer coefficient is estimated using the fol-
is the voltage estimated by dividing the enthalpy due to the lowing equation [16]:
electron drift by the arc current . is given by the following
(8)
equation, in which Bolzmann’s constant and the value of elec-
tronic charge are considered [9]
where is the thermal conductivity of the flowing gas, the air
(5) arc jet is contaminated iron vapor in the calculation, is the
characteristic dimension, the outer diameter of the OPGW in the
where is the arc temperature near the electrode. calculation. Nu is the Nusselt number, which varies depending
These voltages , and vary depending on the elec- on the heat-transfer conditions and is estimated using equations
trode material and the surrounding gas. In this paper, these volt- obtained by experiments. In the calculation, the heat transfer
ages are the ones for the electrode material aluminum and iron due to the flow collision is considered, because the arc jet flow
and the surrounding gas air, because the type of OPGW strand spouting from the counter electrode collides with the OPGW. In
used, as shown in Fig. 1, is AC steel, and dc arc tests are per- this case, Nu varies depending on gap length and the diam-
formed in the air as described in Section II. Although the cen- eter of the counter electrode, and is expressed as the following
tral part material used in the actual strand is steel, the material equation [16] using the Prandtl number (Pr) and the Reynolds
is iron in the calculation because the main component of steel is number (Re). In the case of
iron and its thermal properties are easily obtained from the ref-
erences. Table I shows the voltages used for calculations. is (9)
the voltage estimated by dividing of (3) and (4) by .
is estimated using (5), where the dependence of the arc temper- where C is as follows.
ature on arc current and the electrode material [13], [14] is In the case of
considered. The transferred heat area for is estimated as
an area touched by the arc using the following equation: (10)
(6)
in the case of
where means a function where the smallest value in the
(11)
bracket is picked, is the current density on the OPGW,
is the exposed area of the OPGW mentioned in Section II. ,
in the case of
for the aluminum and iron is 60.0 and 33.8 A/mm , respectively,
and independent of its polarity, which is estimated using the arc
root diameter in the electrode measured in [15].
2) Convective Heat Transfer : As shown in Fig. 2 re-
garding dc arc testing, an arc is ignited between the OPGW and
(12)
the iron counter electrode set above it. A high-current arc can
produce a high-temperature arc jet containing iron vapor from
where Pr and Re are obtained using the following equations:
the counter electrode at high speed from the arc root on the elec-
trode. The convective heat transfer to the OPGW due to the arc (13)
jet flow is calculated using the following [16]:
(7) (14)
where is the heat transfer coefficient, and is the tem- where , and are the specific heat at constant pressure,
perature of the OPGW surface. is estimated assuming the the viscosity, and density and velocity of the flowing gas, which
same value as for . is considered to be room tem- is the air arc jet contaminated with iron vapor in the calculation.
perature (300 K) before dc arc testing, and near melting point These values and are referred to in [13], [14], and [18]–[20].
IWATA et al.: DEVELOPMENT OF A METHOD OF CALCULATING THE MELTING CHARACTERISTICS OF OPGW STRANDS 1317
layer strands of 80 and 140 mm OPGW are investigated con- [8] M. Shibuya, T. Inaba, and Y. Goda, “Molten losses on cut end of rod
sidering the melting and breaking of actual OPGW strands by electrodes in AC high current arc discharge,” Inst. Elect. Eng. Jpn.
Trans., vol. 113-B, pp. 327–334, 1993.
lightning strikes, and the arc current varies from 1–100 kA [9] Y. Kito and T. Mizuno, “An experimental approach to melting of alu-
considering the measured current of actual lightning. In the case minium alloy string at an arc spot formed by surging current,” Inst.
of the actual strand, in the central part, there is steel whose main Elect. Eng. Jpn. Trans., vol. 108-B, pp. 577–584, 1988.
[10] Y. Yokomizu, T. Matsumura, R. Henmi, and Y. Kito, “Total voltage
component is iron, which is covered with aluminum. In the cal- drops in electrode fall regions of , argon and air arcs in current
culation, therefore, the cross-sectional structure of the strand is range from 10 to 20000A,” J. Phys. D, Appl. Phys., vol. 29, pp.
simplified to one which has an aluminum layer on an iron layer. 1260–1267, 1996.
[11] Electrical Discharges Handbook. Chiyoda, Tokyo, Japan: IEEJ,
In the calculation, after a part, which is estimated by mul- 1974, p. 154.
tiplying by the thickness of the aluminum or iron layer, [12] H. B. Michaelson, “The work function of the elements and its period-
is melted, an arc moves to elsewhere and similarly causes the icity,” J. Appl. Phys., vol. 48, no. 11, pp. 4729–4733, 1977.
[13] M. Iwata, S. Tanaka, K. Ikeda, and Y. Goda, “Temperature of arc jet
part to melt. The melting of the upper strands, which are easily of 50 kA class AC Arc in long gap (Part 2)—Temperature and energy
touched by the arc (e.g., the upper three strands in the case of density of arc jet spouted from iron electrode,” Yokosuka, Kanagawa,
80 mm ) is investigated. When parts of the entire cross section Japan, CRIEPI Rep. W00026, 2001.
[14] M. Iwata, K. Ikeda, S. Tanaka, and Y. Goda, “Temperature of arc jet of
region are melted, even if its length along the OPGW axis-di- 50 kA Class AC arc in long gap (Part 3)—Axial distribution of temper-
rection is slight, the melting quantity is considered 100%. ature and energy density,” Yokosuka, Kanagawa, Japan, CRIEPI Rep.
Consequently, the melting time for the cathode is shorter than W01010, 2002.
[15] Y. Goda, M. Iwata, K. Ikeda, and S. Tanaka, “Behaviour of arc jet of 50
that for the anode, because the total heat transferred to kA class high current arc in a long gap,” Yokosuka, Kanagawa, Japan,
the OPGW for the cathode exceeds that for the anode. When CRIEPI Rep. W00034, 2001.
the arc current exceeds 20–30 kA, the gradient of the cal- [16] JSME Data Book: Heat Transfer, 4th ed. Shinjuku, Tokyo, Japan:
JSME, 1986, p. 42, 65.
culation result approaches that of the property where the flowing [17] Tables of Physical Constants, S. Iida, K. Ohno, S. Kamimae, H. Ku-
electric charge is constant, meaning the melting time is in- magaya, and S. Sawada, Eds. Shinjuku, Tokyo, Japan: Asakura Pub.
versely proportional to . Conversely, in the region below the Co. Ltd., 1978, p. 20, 108, 141, and 143.
[18] S. Tanaka, K. Ikeda, M. Iwata, and Y. Goda, “Flow velocity character-
current, the variation of with changing diminishes, which istics of high current arc jet in atmospheric air (Part 2)—Flow velocity
means diminishes more than in the case where is inversely in 50 kA class arc jet spouted from iron and aluminium electrode,”
proportional to . This calculation result shows good agree- Yokosuka, Kanagawa, Japan, CRIEPI Rep. W01021, 2002.
[19] T. Amakawa, “Transport properties of high temperature air with metal
ment with the tendency of the experimental results. These re- vapor,” Yokosuka, Kanagawa, Japan, CRIEPI Rep. W86009, 1986.
sults suggest that the developed calculation method for melting [20] M. Boulos, P. Fauchais, and E. Pfender, Thermal Plasmas—Funda-
characteristics of OPGW strands due to dc arc is a good tool for mentals and Applications. New York: Plenum, 1994, vol. 1, pp.
413–417.
performing the dc arc tests efficiently and clarifying the melting [21] Chronological Scientific Tables, Tokyo Astronomical Observatory
and breaking characteristics of OPGW strands analytically. Fur- ed. Chuo, Tokyo, Japan: Maruzen Co. Ltd., 1988, p. 472.
thermore, the obtained results in this paper also suggest that
Mikimasa Iwata (M’12) was born in Aichi Prefec-
the developed calculation method will be able to estimate the ture, Japan, on May 14, 1964. He received the B.E.,
melting and breaking characteristics of conventional and newly M.E., and Ph.D. degrees in electrical engineering
developed lightning-resistant OPGW strands due to high-en- from Nagoya University, Aichi Prefecture, Japan, in
1988, 1990, and 2000, respectively.
ergy lightning strikes, assuming that the dc arc tests adequately In 1990, he joined the Central Research Institute of
simulate the actual lightning strike. Electric Power Industry (CRIEPI), Yokosuka, Japan.
From 2006 to 2007, he was a Visiting Research As-
sociate of RWTH Aachen University, Aachen, Ger-
many. Currently, he is a Senior Research Scientist of
REFERENCES the Applied High Energy Physics Sector and the High
[1] Y. Goda, S. Yokoyama, S. Watanabe, T. Kawano, and S. Kanda, Power Testing Laboratory, CRIEPI. His research interests are mainly electrode
“Melting and breaking characteristics of OPGW strands by lightning,” erosion and stability of AC arc plasma, high-current arc phenomena on power
IEEE Trans. Power Del., vol. 19, no. 4, pp. 1734–1739, Oct. 2004. equipment, and synthesis of ceramics nanopowder by means of arc plasma.
[2] Y. Goda, M. Shimizu, and A. Matsumoto, “DC arc tests of OPGW Dr. Iwata is a member of the Institute of Electrical Engineers of Japan (IEEJ).
simulating high energy lightning strike,” in Proc. 28th Int. Conf. Light. In 2006, he was awarded the Technical Development Award by the IEEJ.
Protect., 2006, pp. 959–964.
[3] J. Huang, X. Chen, and J. Xu, “Investigation and analysis of an OPGW
lightning strike incident,” in Proc. 58th Int. Wire Cable Symp., 2009,
pp. 389–392. Toshiya Ohtaka (M’03) was born in Kanagawa
[4] Working Group on Lightning Protection for Transmission Systems, Prefecture, Japan, on March 9, 1976. He received
“Lightning observation on Japan sea coast in winter,” Komae, Tokyo, the B.E., M.E., and Ph.D. degrees in electrical
Japan, CRIEPI Rep. T10, 1989. engineering from Waseda University, Tokyo, Japan,
[5] T. Shindo, S. Yokoyama, T. Sunaga, A. Asakawa, H. Motoyama, A. in 1998, 2000, and 2003, respectively.
Wada, and H. Goshima, “Lightning characteristics in winter season to From 2001 to 2003, he was a Research Associate
a high stack—Ten-year observation results from 1989 to 1998 at Fukui at Waseda University, and from 2003 to 2004, he was
Observation Site” Komae, Tokyo, Japan, CRIEPI Rep. T58, 1999. a Visiting Lecturer. In 2004, he joined the Central Re-
[6] M. Böhme and K. Möller, “Arc motion during lightning test of power search Institute of Electric Power Industry (CRIEPI),
ground wires with optical fibers (OPGW),” in Proc. 24th Int. Conf. Yokosuka, Japan, and is currently a Research Scien-
Light. Protect., 1998, pp. 571–575. tist of the Applied High Energy Physics Sector and
[7] S.-H. Xie and R.-S. Yang, “Study of OPGW performance against light- the High Power Testing Laboratory, CRIEPI. His research interests are mainly
ning strikes in UHV transmission lines,” in Proc. 58th Int. Wire Cable high-current arc phenomena and power system analysis.
Symp., 2009, pp. 393–398. Dr. Ohtaka is a member of the Institute of Electrical Engineers of Japan.
IWATA et al.: DEVELOPMENT OF A METHOD OF CALCULATING THE MELTING CHARACTERISTICS OF OPGW STRANDS 1321
Yasukuni Kuzuma was born in Kanagawa Prefec- Yutaka Goda (SM’94) was born in Ehime Prefec-
ture, Japan, on March 29, 1945. He received the B.E. ture, Japan, on April 1, 1956. He received the B.E.
degree in electrical engineering from Kanto Gakuin and Ph.D. degrees in electrical engineering from
University, Kanagawa Prefecture, Japan, in 1969. Nagoya University, Aichi Prefecture, Japan, in 1979
In 1963, he joined the High Voltage Power and 1998, respectively.
Laboratory (HVPL), Yokosuka, Japan. Since 1977, In 1979, he joined the Central Research Institute of
he has been with the Central Research Institute of Electric Power Industry (CRIEPI), Yokosuka, Japan.
Electric Power Industry (CRIEPI), Yokosuka, Japan. From 1991 to 1992, he was a Visiting Research As-
Currently, he is on the staff of the High Power sociate at the University of Connecticut, Storrs. Cur-
Testing Laboratory, CRIEPI. His research interests rently, he is a Senior Research Scientist and Director
are mainly arc phenomena on power equipment. of the High Power Testing Laboratory, CRIEPI. His
Mr. Kuzuma is a member of the Institute of Electrical Engineers of Japan. research interests are mainly arc phenomena on power equipment.
Dr. Goda is a member of the Institute of Electrical Engineers of Japan
(IEEJ) and Working Group 12 of IEC-TC42. He was awarded the Technical
Development Award and the Best Article Award by the IEEJ in 1997 and 1998,
respectively.