Next Article in Journal
Biological Scaffolds for Abdominal Wall Repair: Future in Clinical Application?
Next Article in Special Issue
How to Choose the Superconducting Material Law for the Modelling of 2G-HTS Coils
Previous Article in Journal
A Nondestructive Indirect Approach to Long-Term Wood Moisture Monitoring Based on Electrical Methods
Previous Article in Special Issue
Microstructure and Flux Pinning of Reacted-and-Pressed, Polycrystalline Ba0.6K0.4Fe2As2 Powders
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Study on Quenching Characteristics and Resistance Equivalent Estimation Method of Second-Generation High Temperature Superconducting Tape under Different Overcurrent

1
State Key Laboratory of Advanced Electromagnetic Engineering and Technology, School of Electrical and Electronic Engineering, Huazhong University of Science and Technology, Wuhan 430074, China
2
School of Electrical Engineering, Beijing Jiaotong University, Beijing 100044, China
3
Electric Power Research Institute of Guangdong Power Grid Corporation, Guangzhou 510080, China
*
Author to whom correspondence should be addressed.
Materials 2019, 12(15), 2374; https://doi.org/10.3390/ma12152374
Submission received: 6 June 2019 / Revised: 30 June 2019 / Accepted: 24 July 2019 / Published: 25 July 2019
(This article belongs to the Special Issue Engineering Properties of Superconducting Materials)

Abstract

:
In this paper, through AC and DC overcurrent tests on second generation high temperature superconducting tape (2G HTS tape), we respectively summarize the typical types of quenching resistance and corresponding quenching degree, in which there are three types under AC overcurrent and two types under DC overcurrent. According to experimental results, a rule was found that, when 2G HTS tape quenches to normal state, the relationship between quenching resistance and joule heat generated from 2G HTS tape presents a fixed trend line, and the influence of liquid nitrogen can be ignored. Then, the characteristics and rules of quenching resistance found in experiments are well explained and confirmed by a detailed 3D finite element model of 2G HTS tape including electromagnetic field and thermal field. Finally, based on above works, our group proposes a new equivalent method to estimate the quenching resistance, where the results of AC and DC overcurrent experiments can be equivalent to each other within a certain range. Compared with FEM, the method has the following advantages: (i) The method is simple and easy to implement. (ii) This method combines precision and computational efficiency. (iii) With superconducting tape quenching to normal state, this method presents a good consistency with experimental results.

1. Introduction

With the discovery of the cuprate-based high temperature superconductors, first generation high temperature superconducting (1G HTS) tapes represented by BSCCO (Bismuth Strontium Calcium Copper Oxide) Ag-sheathed conductors and second-generation high temperature superconducting (2G HTS) tapes represented by YBCO (Yttrium Barium Copper Oxide) coated conductors have appeared successively [1,2]. Compared with 1G HTS tape, through the improvement of manufacturing process, YBCO HTS tapes have the advantages of higher current density, lower alternating current (AC) loss and lower theoretical cost [3,4,5,6,7], which provide an opportunity for the high power application of superconducting devices applied in the field of electric power.
When current exceeds critical current, superconducting tape will automatically quench and switch to resistance state, which can limit short circuit fault current in power system. Based on the quenching characteristics of YBCO material, resistance type superconducting fault current limiter (R-SFCL) has advantages of simple structure, automatic trigger and fast response. Therefore, R-SFCL can be widely applied to limit short circuit fault current in multiple application scenarios, such as ship power grid [8], railway direct current (DC) traction systems [9], microgrid system [10,11], DC power grid [12,13,14] and so on.
However, HTS tape is very sensitive to magnetic field, temperature field and stress field. The complex working environment formed by multiple physical fields has a great influence on the homogeneity and stability of superconducting tape. Therefore, the applications of R-SFCL in power system requires the manufacture of superconducting tape with excellent performance and obvious quenching characteristics. It is necessary to focus on studying the quenching characteristic of YBCO material to realize the engineering application of R-SFCL. At present, many scholars have conducted extensive and in-depth research on the different characteristics of YBCO superconducting tape including quenching behavior under AC or DC overcurrent [15,16,17,18,19,20,21,22], quench propagation [23,24,25,26], quenching recovery [27,28,29,30], thermal and mechanical properties [31,32], minimum quench energy [33,34], maximum operating condition [35], and the effect of different materials on quenching behavior [36,37,38], which provide a solid foundation for the use of superconducting tape in design and fabrication of R-SFCL. However, there are few studies on common characteristics of quenching resistance of 2G HTS tape under AC and DC overcurrent, which is necessary and basic for R-SFCL design and research.
In addition, it is very difficult to carry out quenching test on a large-scale DC R-SFCL prototype in practical engineering applications before installed in large capacity voltage source converter (VSC) or modular multilevel converter (MMC) high voltage direct current (HVDC) projects. The transmission power of the HVDC projects usually reach hundreds of megawatts. It is very risky to carry out short circuit fault for testing R-SFCL in the existing DC system projects. Meanwhile, the economic and technical cost of manufacturing large-capacity DC test platform for R-SFCL is very high. On the contrary, the technology of large-scale AC experiment platform is mature and its cost is lower. At present, with the support of national key research and development plan, a large-scale DC R-SFCL is being developed in China, which is planned to be installed in Nanao’s ±160 kV MMC three-terminal HVDC system. If the common characteristics and internal relations of 2G HTS tape under different type overcurrent is found, the AC overcurrent experimental results of the R-SFCL can be used to evaluate quenching resistance of R-SFCL under DC overcurrent.
In this study, a kind of YBCO tape produced by Shanghai Superconductor Company for R-SFCL was tested under AC and DC overcurrent. By comparing the experimental results of AC and DC overcurrent test, the common characteristics and laws of quenching resistance variation were found under different overcurrent in 77 K liquid nitrogen immersion environment. The same conclusions were obtained through the simulation of YBCO three-dimensional (3D) finite element model. Based on the findings, a new estimation method is proposed, which can accurately reflect and estimate the quenching resistance within a certain range, whose feasibility and availability was verified by comparing with experimental results. Furthermore, it provides theoretical basis and experimental support for AC equivalent test method before a large-scale DC R-SFCL connected to DC power grid.

2. 2G HTS Tape

The YBCO superconducting tapes used in the experiment are produced by Shanghai Superconductor Company (Shanghai, China) [39], which are mainly adopted in R-SFCL. In general, the typical structure of YBCO tape consists of copper stabilizer, silver stabilizer, YBCO layer and Hastelloy substrate. In addition, a buffer layer exists between YBCO layer and Hastelloy substrate.
The 2G HTS tape is 12 mm wide with a two-core configuration, including two typical structures of YBCO tape with copper stabilizer on both side, which are stacked back-to-back and encapsulated with stainless steel layer, as shown in Figure 1a,b (buffer layer is omitted). This structure can increase self-field critical current, which is suitable for power systems with high current. The corresponding basic parameters of the 2G HTS are listed in Table 1. As shown in Figure 1c, the middle 10 cm area of tape samples are selected for voltage measurement in AC and DC overcurrent tests.

3. AC and DC Overcurrent Experiment

In a strict sense, a current whose magnitude and direction change periodically with time is defined as AC current, while a current whose magnitude and direction remain constant is defined as a DC current. However, in the actual system short circuit fault, the overcurrent will present irregular, even drastic changes. Therefore, in this paper, we define AC and DC overcurrent in a broad sense. The overcurrent in constant direction is DC overcurrent, and the overcurrent in direction varying with time is AC overcurrent. The two kinds of overcurrent analyzed in the experiments are shown in Figure 2.

3.1. AC and DC Overcurrent Experiment Platforms

The AC and DC overcurrent test platforms were, respectively, built for studying quenching characteristics of YBCO tape, as shown in Figure 3 and Figure 4, and their basic corresponding parameters are listed in Table 2. The parameters setting of AC and DC overcurrent experiments are listed in Table 3. The YBCO tape to be tested was immersed in liquid nitrogen and cooled to a superconducting state at 77 K.
The AC overcurrent test platform was equipped with a 20 kVA strong current generator, which can turn the 220 V power supply into a 6 V high-power current source with a maximum output of 3 kA. Output capacity was adjusted by strong current generator. The on–off state of thyristor decides duration of AC overcurrent, which was controlled by the signal generator and controller.
The DC overcurrent test platform was the series circuit of resistance R, inductance L and capacitance C (RLC). The RLC series circuit was used to approximately equivalent and simulate the capacitor discharging process of VSC and MMC systems in the DC side short-circuit fault [40,41]. When S1 is in the closed state and S2 is in the disconnected state, charging circuit begins to charge the capacitor until the voltage UC reaches a set value. When S1 is disconnected and S2 is closed, R0, L0 and C0 form a closed loop and generate a DC overcurrent for the test of YBCO tape. Peak value, rising slope and duration of DC overcurrent waveform vary with adjusting UC, L0 and R0.

3.2. The results and Analysis of AC Overcurrent Experiment

Eight groups of AC overcurrent tests on YBCO tape were carried out by adjusting the amplitude and duration of overcurrent. The amplitude of overcurrent gradually increased from 860 A to 1920 A. The range of overcurrent duration was 50–80 ms. The eight-group experimental results are shown in Figure 5. The AC experimental analysis is summarized in Table 4.
As shown in Figure 5a, continuous half-wave type quenching resistance suggests that the quenching state of superconducting tape cannot be maintained, which is unstable and discontinuous. With AC overcurrent decreasing, the superconducting tape returns back to superconducting state and quenching resistance tends to 0 Ω. Therefore, this quenching state is inadequate. At this time, the state of superconducting tape is defined as partial resistive state.
As shown in Figure 5b,c, there are two types resistance waveform under this AC overcurrent level, continuous half-wave type and incremental curve type. In initial stage, superconducting tape recovers rapidly after quenching. Then, superconducting tape begins to quench continuously and transits from fluctuation to continuous increment. It indicates that the quenching degree of superconducting tape gradually deepens and further transits to normal state. This quenching stage of superconducting tape is defined as transition stage, where the state of superconducting tape starts to transit from partial resistive state to normal state.
As shown in Figure 5d–h, incremental curve type quenching resistance reflects that superconducting tape responds to quench continuously without fast recovering under this AC overcurrent level. The quenching state in superconducting tape is adequate and stable, which is considered as complete quenching. It indicates that the superconducting tape is in normal state.

3.3. The Results and Analysis of DC Overcurrent Experiment

The eight-group DC experimental results are presented in Figure 6. The amplitude of overcurrent gradually increases from 910 A to 2020 A. The range of overcurrent duration is 20–23 ms. The DC experimental analysis is summarized in Table 5.
As shown in Figure 6a, half-wave type quenching resistance suggests that the quenching state of superconducting tape cannot always be maintained under this DC overcurrent level. When DC overcurrent decreases, superconducting tape returns to superconducting state rapidly and quenching resistance tends to 0 Ω. Therefore, this quenching state is inadequate, which is considered as incomplete quenching. Similar to the result of AC experiment in Figure 5a, the superconducting tape is in partial resistive state.
As shown in Figure 6b–h, the quenching degree of superconducting tape is strengthened with the peak value IPmax of DC overcurrent increasing. Because DC overcurrent has no zero crossing before current declines to 0, quenching resistance is incremental curve type. They reflects that superconducting tape quenches continuously without fast recovering under this DC overcurrent level. During this stage, the superconducting tape completely quenches to normal state and presents normal resistance.

3.4. Comparative Analysis of Quenching Resistance under AC and DC Overcurrent

According to the above results of AC and DC overcurrent experiment shown in Figure 5 and Figure 6, the corresponding R-t curves of superconducting tape are obtained, respectively, shown in Figure 7a,b. The R-t curves reflect resistance varying with time under different overcurrent. With the increase of overcurrent amplitude, the quenching degree of superconducting tape is enhanced. However, the rising slope and steady-state value of R-t curves vary with waveform type, amplitude and duration of impact current. Therefore, R-t relationship cannot be used to describe the variation of quenching resistance in superconducting tape.
Temperature plays a key role in quenching resistance and quenching degree in superconducting tape. In a certain cryogenic medium, joule heat accumulation is the fundamental cause of temperature rise. Without considering the variation of cryogenic medium, the main factor affecting the temperature of superconducting tape should be the joule heat generated during the quenching process, which is the result of comprehensive effect from quenching resistance, current and quenching duration. The generated joule heat Q of the superconducting tape varying with time can be calculated by Equation (1), whose unit is joule (J).
Q ( t 2 ) = t 1 t 2 I 2 R ( t )   d t
where t1 is the initial time of quenching, R(t) is quenching resistance corresponding to time, and the range of time t is [t1, t2]. According to sixteen-group data of experimental results, R-Q curves of per meter superconducting tape are obtained, where dependent variable R is the quenching resistance of per meter superconducting tape and independent variable Q is the generated joule heat of per meter superconducting tape, as shown in Figure 8.
Figure 8a is the comparison diagram of R-Q curves with different overcurrent, and Figure 8b is the partial enlarged drawing of Figure 8a. It is quite clear that twelve R-Q curves obtained are basically overlapping and their trend line are consistent under AC and DC overcurrent except the results of DC No. 1, AC No. 1, AC No. 2 and AC No. 3. According to the classification in Section 3.2 and Section 3.3, DC No. 1 and AC No. 1 belong to partial resistive stage of superconducting tape, while AC No. 2 and No. 3 belong to transition stage of superconducting tape.
Therefore, the experimental results present a common trait. When YBCO superconducting tape completely quenches in a cryogenic medium, though superconducting tape is subjected to different overcurrent, the quenching resistance R has an intrinsic fixed relationship with the generated joule heat Q in superconducting tape. In other words, Q is the deciding factor of quenching resistance, while overcurrent is the only inducing factor of quenching resistance. The preliminary finding suggests that, during the complete quenching stage of superconducting tape, quenching resistance can be described by a fixed R-Q curve, which can be obtained from experiments in comprehensive consideration of current, resistance and duration. Therefore, R should satisfy an incremental relationship with Q when superconducting tape tends to complete quenching, described by Equation (2), which also can be reflected by an R-Q curve.
RQ(I, R, t)
Based on experimental phenomena of superconducting tapes, the quenching resistance characteristics and laws of superconducting tapes are summarized. However, experimental data alone are insufficient to provide adequate physical explanations for the phenomena and laws. Therefore, as presented in Section 4, commercial finite element software COMSOL Multiphysics was adopted to establish 3D simulation model of 2G HTS tape for exploring the criteria of quenching resistance classification, validating the R-Q curve and providing corresponding physical explanations.

4. Simulation Study on 2G HTS Tape Quenching under AC and DC Overcurrent

4.1. 3D Finite Element Model of 2G HTS Tape

In software COMSOL Multiphysics, according to the basic parameters of the superconducting tape listed in Table 1, the 3D multilayer structure of the 2G HTS tape was established as superconducting tape domain, whose surface is covered with liquid nitrogen domain. A combination of custom partial differential equation module (PDE module) and heat transfer module were adopted to perform transient coupling calculation of electric field, magnetic field and thermal field. PDE module was used to carry out electromagnetic calculation. Heat transfer module was used to carry out thermal calculation. The two modules were coupled by joule heat and temperature. The basic structure of the simulation model is shown in Figure 9.

4.1.1. 3D Electromagnetic Model

The 3D electromagnetic calculation was based on Maxwell equations and H equation [42] with good convergence. The Faraday’s law and the ampere-law equation, respectively, are described as Equations (3) and (4).
× E = B t = μ 0 μ r H t
× H = J + D t
where E is electric field intensity, B and H , respectively, are magnetic induction intensity and magnetic field intensity, and J and D are conduction current density and displacement current density. Since the conduction current density is much higher than displacement current density in the superconducting tape, the differential term of displacement current with respect to time in Equation (4) can be ignored, thus × H = J .
Because B   = μ 0 μ r H and E = r J , the partial differential Equation (5) can be obtained according to Equations (3) and (4), where μ0 is permeability of vacuum, 4π × 10−7 H/m, μr is relative permeability, and r is electric resistivity.
× ( r × H ) = μ 0 μ r H t
where × H = J x e x + J y e y + J z e z . Jx, Jy and Jz are the components of J in the x, y, and z directions, which can be expressed by Equations (6).
J x = H z y H y z ,   J y = H x z H z x ,   J z = H y x H x y
Combining with Equations (5) and (6), the governing equation of electromagnetic calculation are obtained as follows.
[ 0 E z E y E z 0 E x E y E x 0 ] [ e x e y e z ] + t [ e x e y e z ] [ μ 0 μ r 0 0 0 μ 0 μ r 0 0 0 μ 0 μ r ] [ H x H y H z ] = 0
where = [ x ,   y ,   z ] . e x , e y , and e z are, respectively, the unit direction vector of x, y, and z directions. Ex, Ey, and Ez are the components of E in the x, y, and z directions. Hx, Hy, and Hz are the components of H in the x, y, and z directions, which are set as dependent variables in PDE module.
Overcurrent was used as model excitation in simulation. As a composite conductor, each layer of the superconducting tape was assumed to be parallel structure. Therefore, the current constraint condition of the model should meet Equation (8), which is load on the cross section of superconducting tape by pointwise constraint.
I i n p u t = I S C + I C u + I A g + I S u b + I S t
where Iinput is total current excitation, ISC is the current in YBCO layer, ICu is the current in copper layer, IAg is the current in silver layer, ISub is the current in Hastelloy substrate, and ISt is the current in stainless steel layer.
In the quenching model of superconducting tape, the calculation of YBCO resistivity is crucial; a piecewise formula for YBCO resistivity was adopted, which is expressed as Equations (9). Critical temperature Tc is used as criterion for state transition. When TTc, the superconducting tape is considered as incomplete quenching, whose YBCO resistivity adopts the classical method of parallel equivalent electrical circuit [43]. When T > Tc, the superconducting tape is considered as complete quenching, whose YBCO resistivity adopts constant resistivity rnorm.
r Y B C O = { ( r 1 + r 2 + r 0 ) r n o r m r 0 + r 1 + r 2 + r n o r m    T T c r n o r m    T > T c
where r0 is residual resistivity [44], r1 is expressed as Equation (10), r2 is expressed as Equation (11). According to the authors of [42], n1, n2 and E0 are fitting parameters of two power-law relations, and k can be calculated based on r1(3Jc) = r2(3Jc).
r 1 = { 0 | J | < J c E 0 | J | ( | J | J c k ) n 1 | J | J c
r 2 = { 0 | J | < k J c E 0 | J | ( | J | J c k ) n 2 | J | k J c
According to the authors of [45], when the temperature T of YBCO layer is below critical temperature Tc, the critical current density Jc of YBCO material is affected by temperature T. When T is above Tc, Jc is considered as 0 A/m2. The relationship can be expressed by the improved typical Equation (12) in [43].
J c ( T ) = { J c 0 ( T c T T c T r e f ) α ( T r e f < T < T c ) 0 ( T c T )
where Tref is initial ambient temperature or cryogenic medium temperature of the superconducting tape. Jc0 is critical current density corresponding to Tref. In this simulation, superconducting tapes were cooled by liquid nitrogen immersion. Therefore, Tref was set to 77 K. The relevant simulation parameters of YBCO resistivity are listed in the Table 6.

4.1.2. Heat transfer Model

In heat transfer calculation, because superconducting tape was immersed in liquid nitrogen, which is without obvious flowing, a heat balance Equations (13) in the module of Heat Transfer in Solids in COMSOL Multiphysics was used to calculate the thermal characteristics of the entire superconducting tape region including YBCO layer, Cu layer, Ag layer, stainless steel and Hastelloy substrate.
{ q s = ρ c T t + q q s = E x J x + E y J y + E z J z q = k T
where qs is volume power density, with unit W/m3; ρ is mass density, with unit kg/m3; c is specific heat capacity, with unit J/(kg·K); q is conduction heat flux, with unit W/m2; and k is heat transfer coefficient, with unit W/(m2·K).
In the simulation, the initial temperature of thermal solution domain was set to 77 K, and the boundary condition of the superconducting tape surface was set to Heat Flux. To simulate the heat exchange process between superconducting tape surface and liquid nitrogen, a heat transfer coefficient curve in Heat Flux was set on the interface between equivalent superconducting tape and liquid nitrogen [43,46,47], which is a common method to simulate heat exchange process between liquid nitrogen and superconducting tape surface at different temperatures. As shown in Figure 10, this is a heat transfer coefficient curve of liquid nitrogen including free convection, nuclear boiling, transition boiling and film boiling [48].

4.2. Verification of the 3D Finite Element Model

As shown in Figure 11a–d, four-group quenching data of superconducting tape under different AC overcurrent were selected for comparison between experiment and simulation, covering three types quenching resistance waveforms of superconducting tapes under AC overcurrent. Figure 11e–h compares the results of experiment and simulation, covering two types quenching resistance waveforms of superconducting tape under DC overcurrent.
Through waveform comparison between experiment and simulation, it can be found that the simulation phenomena are basically consistent with the classification of experimental results, where the model can correctly describe the different types of quenching resistance under different overcurrent. Therefore, the 3D finite element model can correctly reflect the quenching state of the superconducting tape, and contributes to research the internal evolution process of superconducting tape under AC and DC overcurrent.

4.3. Thermal Characteristics and Current Distribution Characteristics of 2G HTS Tape Quenching

4.3.1. Resistance, Thermal and Current Distribution Characteristics under AC Overcurrent

To analyze the transient characteristics of superconducting tape under AC overcurrent and obtain corresponding common laws, it was necessary to select and analyze the three typical quenching resistance waveforms of superconducting tape found in Section 3.2 under AC overcurrent, including half-wave type, half-wave and half-incremental curve type and incremental curve type. Therefore, the simulation results with AC 860 A, AC 1200 A and AC 1780 A were selected.
Figure 12 is the simulation results of the superconducting tape quenching under AC 860 A, which represent the characteristics of continuous half-wave type quenching resistance in partial resistive state. According to Figure 12a, when 0 ≤ t < t1, the superconducting tape is in superconducting state with 0 Ω, and current only flows through YBCO layer. When t1t, current begins to exceed the critical current IC which results in YBCO layer quenching. However, the ratio of real-time resistivity to normal resistivity of YBCO is very low, no more than 0.07. It indicates that the superconducting tape is in partial resistive state. As shown in Figure 12b, during t1 to 60 ms, the temperature of YBCO layer rises continuously in fluctuation, but the maximum value of temperature is only 88.7 K, which deso not exceed the critical temperature TC of YBCO layer. Therefore, during the whole quenching process, the quenching resistivity of YBCO layer is in the first stage in Equation (9), which causes the current to be diverted to other layers. The resistivity of YBCO layer fluctuates with the rise and fall of AC overcurrent, which results in the superconducting tape producing the corresponding continuous half-wave type quenching resistance. During initial quenching, the current distribution among material layers is: IYBCO > ICu > IAg > ISt > ISub. With the temperature of YBCO layer rising, the quenching resistivity of YBCO layer increases, and the ordering of current in proportion is: ICu > IYBCO > IAg > ISt > ISub.
Figure 13 is the simulation results of the superconducting tape quenching under AC 1200 A, which represents the characteristics of continuous half-wave and half-incremental curve type quenching resistance in transition state. According to Figure 13a, when 0 ≤ t < t1, the superconducting tape is in superconducting state with 0 Ω, current only flow through YBCO layer. When t1t < t2, the ratio of real-time resistivity to normal resistivity of YBCO is very low, not reaching 1. Therefore, the quenching characteristics and current distribution characteristics of this condition are the same as those in partial resistive state. According to Figure 13b, when t2t, the temperature of YBCO layer exceeds the critical temperature TC. Its maximum temperature is 92.3 K, which results in the quenching resistivity of YBCO transformed into constant resistivity in Equation (9). At this time, the ratio of real-time resistivity to normal resistivity of YBCO is 1, which indicates that YBCO layer transitions from partial resistive state to normal state. Therefore, current mainly flows through other layers except YBCO layer. During this period, the characteristics of current distribution is: ICu > IAg > ISt > ISub > IYBCO. The superconducting tape further presents a composite conventional metal resistivity. However, with the decline of AC overcurrent, the superconducting tape gradually reaches thermal equilibrium. The temperature fluctuation near 90 K causes the quenching degree of YBCO layer to be in a critical state, fluctuating between partial resistive state and normal state.
Figure 14 is the simulation results of the superconducting tape quenching under AC 1780 A, which represent the characteristics of continuous incremental curve type quenching resistance in normal state. In this condition, the duration of partial resistance state is very short. Current rapidly exceeds the critical current of YBCO layer. Temperature correspondingly exceeds critical temperature 90 K within first current half wave, thereby resulting in YBCO layer transformed into normal state within about 5 ms, as shown in Figure 14. With temperature continuously rising, the superconducting tape presents corresponding incremental curve type quenching resistance. The maximum temperature of YBCO layer is 174 K. During most of quenching process, the current distribution in superconducting tape is mainly: ICu > IAg > ISt > ISub > IYBCO.

4.3.2. Resistance, Thermal and Current Distribution Characteristics under DC Overcurrent.

The two typical quenching resistance waveforms of superconducting are found under DC overcurrent in Section 3.3 including half-wave type and incremental curve type. Therefore, the simulation results with DC 910 A and DC 2020 A were selected to analyze the transient characteristics of superconducting tape under DC overcurrent to obtain the corresponding common laws.
Figure 15 is the simulation results of the superconducting tape quenching under DC 910 A, which represent the characteristics of half-wave type quenching resistance in DC partial resistive state, whose quenching characteristics are same as those of AC partial resistive state. When 0 ≤ t < t1, the superconducting tape is in superconducting state with 0 Ω. Therefore, DC overcurrent only flows through YBCO layer. When t1t, DC overcurrent exceeds the critical current of YBCO layer, but the temperature is still below the critical temperature 90 K. During the whole quenching process, the maximum temperature of YBCO layer is 82.9 K. Therefore, the quenching resistivity of YBCO layers is far below normal resistivity, as shown in Figure 16a, resulting in the superconducting tape recovering quickly with DC overcurrent attenuation, thereby forming half-wave type quenching resistance. The current distribution characteristic of superconducting tape is: IYBCO > ICu > IAg > ISt > ISub in Figure 15b).
Figure 16 is the simulation results of the superconducting tape quenching under DC 2020 A, which represent the characteristics of incremental curve type quenching resistance in DC normal state. The quenching characteristics are the same as those of AC normal state. Because the temperature of YBCO layer exceeds the critical temperature and increases continuously, the resistivity of YBCO layer quickly transits from partial resistive state to normal state, as shown in Figure 16a. Therefore, current only flows through other layers except YBCO layer. The equivalent quenching resistance of the superconducting tape is composite resistance of metal materials depending on the temperature, which presents incremental curve type with temperature rising. Although DC overcurrent starts to decay after the first crest, the temperature is maintained at 192 K. Therefore, the superconducting tape remains at normal resistivity state, whose recovery is very slow. In this condition, current distribution characteristics is: ICu > IAg > ISt > ISub > IYBCO (Figure 16b).

4.4. Consistency Analysis of R-Q Curve

According to the simulation results shown in Figure 11, Equation (2) is adopted to plot the R-Q curves shown in Figure 17, which present the same phenomenon as the experimental results. When superconducting tape is in normal state in a cryogenic medium, the variation tendency of R-Q curves under different AC and DC overcurrent is consistent, while the partial resistive state and transition phenomena do not appear.
To explain the phenomena and analyze the consistency in R-Q curves in normal state, the heat transfer process of each simulation result is plotted in Figure 18, including heat power Pheat generated by superconducting tape, cooling power Pcool of liquid nitrogen on superconducting tape surface and temperature T of superconducting tape. In the consideration of the heat exchange process between superconducting tape and liquid nitrogen, accumulated joule heat Q1 can be calculated by Equation (14).
Q 1 ( t 2 ) = t 1 t 2 ( P h e a t P c o o l ) d t
Q 1 = c m Δ T
where c is specific heat capacity, m is mass, and ΔT is equal to T-Tref. According to Equation (15), accumulated joule heat Q1 directly causes the temperature T to rise, which can further aggravate quenching.
When superconducting tape is in partial resistive state or transition, the temperature of superconducting tape is below 100 K. At this time, liquid nitrogen is in free convection, nuclear boiling or transition boiling. During this time, the cooling power of liquid nitrogen is relatively high, which greatly weakens the quench of superconducting tape. Therefore, the effect of Pcool cannot be ignored, as shown in Figure 18a,b,e. In addition, compared with the heating power of superconducting tape, the cooling power of liquid nitrogen shows obvious hysteresis phenomenon. Therefore, heating power increases and decreases rapidly with the change of overcurrent, resulting in temperature rising first and then dropping. During partial resistive state and transition, R-Q curves without considering Pcool cannot correctly reflect the changing rule of quenching resistance, which further validates and explains that the R-Q curves are not consistent at these quenching stages in the experimental results in Figure 8b.
When superconducting tape is in normal state, the surface temperature of the superconducting tape is very high, causing the surrounding liquid nitrogen to be in film boiling. Because the cooling power of liquid nitrogen is much less than the heating power of superconducting tape (Pheat >> Pcool), the thermal equilibrium is quickly broken by a huge thermal shock. Hence, superconducting tape is approximately in an adiabatic environment, whose heat exchange with liquid nitrogen can be ignored, as shown in Figure 18c,d,f–h. At this time, the superconducting tape is in thermal runaway; quenching resistance increases rapidly with temperature rising. Pcool can be omitted, thus the amount of joule heat Q generated by superconducting tape can be approximately equivalent to accumulated joule heat Q1.
In the complete quenching process of superconducting tape, due to T > Tc, the quenching resistivity of YBCO layer is normal resistivity rnorm. The resistivity of other materials have fixed curves varying with temperature T, as shown in the Appendix. The equivalent resistivity of the superconducting tape can be calculated by Equation (16).
r e q = 1 / ( i = 1 n f i r i )
R e q = r e q l S
where ri is the resistivity of each material layer. fi is the volume percentage of each material in the superconducting tape. l is the length of superconducting tape and S is the cross-section area. According to Equations (16) and (17), the equivalent resistance of 1 m superconducting tape varying with temperature is a fixed curve, as shown in Figure 19a.
The specific heat capacity c of each material has a fixed curve varying with temperature T, as shown in the Appendix. The equivalent specific heat capacity ceq of the superconducting tape can be calculated by Equations (18), as shown in Figure 19b.
{ c v = 1 / ( i = 1 n f i / ( γ i c i ) ) c p = i = 1 n f i γ i c i c e q = ( c v + c p ) / 2
where cv is the equivalent specific heat capacity in the vertical direction of the superconducting tape surface, cp is the equivalent specific heat capacity in the parallel direction of the superconducting tape surface, ci is the specific heat capacity of each material, and γi is the mass density ratio of each material layer to the whole superconducting tape.
With the superconducting tape quenching to normal state, the temperature of superconducting tape is mainly affected by the amount of joule heat generated with time. Therefore, during this condition, the R-Q curves of superconducting tape under different AC and DC overcurrent are consistent, which can be used to reflect the change law of quenching resistance. According to Equations (14)–(18), the theoretical R-Q curve in normal state can be calculated by numerical calculation, which is compared with experiment results and FEM results in Figure 20.
The comparison in Figure 20 obviously indicates that, when T > TC, the R-Q curves of superconducting tape are basically consistent with the theoretical result. In normal state, critical temperature TC is the turning point of quenching resistance variation. The accumulated joule heat required is about 92.14 J. Furthermore, it also confirms the correctness of the theoretical analysis of R-Q curves consistency mentioned above. When the superconducting tape enters normal state, it is approximately in an adiabatic environment, and quenching resistance will vary with the accumulated joule heat along a fixed R-Q curve, which is consistent under AC and DC overcurrent.

4.5. The Summary of Quenching Characteristics under AC and DC Overcurrent

Based on the above simulation and experimental analysis, the superconducting tape is divided into three stages: superconducting state, partial resistive state and normal state. Their classified method is shown in Figure 21 and corresponding characteristics are listed in Table 7. Critical current IC(T) is the criterion of whether or not to quench. Critical temperature TC is the criterion of quenching degree. The consistency of R-Q curves is satisfied in normal state.

5. The Equivalent Estimation of Quenching Resistance under AC and DC Overcurrent

5.1. The Basic Principle of the New Method

Combined with experiment and FEM simulation, a general rule is further verified. When superconducting tape quenches to normal state, immersed in liquid nitrogen environment, the superconducting tape has a correspondence relationship between quenching resistance and the accumulated joule heat. The R-Q curve is consistent under different AC and DC overcurrent. Therefore, a new calculation method of quenching resistance is proposed to estimate and predict variation of quenching resistance when superconducting tape is in normal state. The joule heat adopted in the new method can effectively evaluate the quenching resistance of superconducting tape and ignore the difference of current waveforms.
The method is defined as “R-Q curve method”, whose calculation process is shown in Figure 22. In the calculation process, test current is used as input current. IC is quenching criterion. R0 is set as initial quenching resistance for initial energy accumulation. Δt is the time step of calculation. When superconducting tape starts to quench, the accumulated joule heat Q in each time step is calculated by discrete calculation, whose corresponding quenching resistance R can be obtained by interpolation calculation based on R-Q curve. In addition, under overcurrent shock, the effect of liquid nitrogen on superconducting tape can be approximately ignored. Therefore, the temperature variation of superconducting tape can be calculated by Equation (19).
T = Q c m + T r e f

5.2. The Validation of R-Q Curve Method

According to Section 3 and Section 4, under AC and DC overcurrent, the R-Q curves of superconducting tape obtained have the same trend. Therefore, the experimental data with the most extensive coverage are adopted to be converted into the R-Q curve of per meter superconducting tape, as shown in Figure 23, whose data are from AC 1920 A overcurrent test. The simulation parameters setting are listed in Table 8.
According to Figure 24, with comparison between simulations and experiments, the effectiveness of R-Q curve method is well verified, which can be used to simulate the variation process of quenching resistance of small-scale superconducting tape in normal state under AC and DC overcurrent. In addition, the R-Q curve in the simulation comes from the AC experimental data, but the model of R-Q curve method still has good simulation accuracy under DC overcurrent. It also indicates that the quenching resistance of AC and DC overcurrent test can be equivalent by R-Q curve in normal state. Meanwhile R-Q curve method can calculate temperature change in the quenching process, whose results are basically consistent with the results of FEM but computing speed is faster, as shown in Figure 25.
Obviously, in terms of estimating quenching resistance of superconducting tape, R-Q curve method has a good calculation accuracy in normal state. However, compared with FEM, because of R-Q curve method using interpolation calculation while FEM adopting the coupling calculation of multiple physical fields, the model structure of R-Q curve method is simpler and faster, as shown in Table 9. However, it is undeniable that the application range of R-Q curve method is limited.
Having good simulation precision and calculation speed, the R-Q method can be applied to evaluate and predict the quenching resistance of large-scale non-inductive superconducting coil in R-SFCL based on the experimental data of small-scale superconducting tape. The prediction results of the model is compared with the experimental results of the coil to verify this estimation method.
To calculate the quenching resistance of superconducting tape with any length, the model of R-Q curve method is improved in PSCAD, as shown in Figure 26, whose R-Q curve data comes from the experimental results of 10 cm superconducting tape converted to the R-Q curve of 1 m superconducting tape applied in the model. The coefficient K is the length of superconducting tape, which can be set according to the length of experimental sample including the length of superconducting tape in non-inductive superconducting coil.
The non-inductive superconducting coil is shown in Figure 27, which adopts single 12 mm width superconducting tape winding. Its diameter is 1200 mm and its length of superconducting tape is 136 m. Because the superconducting coil adopts non-inductive design, the effect of magnetic field on superconducting tape can be neglected. According to Figure 28, the peak value of DC overcurrent is 3.7 kA, which results in superconducting coil rapidly entering in complete quenching stage. The inductance of superconducting non-inductive coil can be ignored, but its resistance needs to be considered, where resistance and accumulated joule heat are proportional to the length of superconducting tape. Without considering the influence of coil structure, the quenching resistance of 136 m superconducting tape calculated by R-Q curve method is basically consistent with the experimental results of actual superconducting coil. It indicates that the non-inductive coil structure can effectively promote the uniform quench of the superconducting coil.
In addition, Figure 28 clearly shows that R-Q curve method has good simulation accuracy and fast solving speed compared with the experimental results in the time range of tens of milliseconds. The R-Q method can rapidly calculate the quenching resistance and average temperature of large non-inductive superconducting coils during complete quenching, only taking a few seconds. Therefore, the method can be improved to predict the dynamic quenching resistance and temperature of superconducting coil for R-SFCL design.

6. Conclusions

In this paper, through the combination of AC and DC overcurrent experiment and simulation study, the quenching characteristics and common laws of YBCO superconducting tape applied in R-SFCL are summarized as follows:
(1) Under AC overcurrent, the quenching resistance of the HTS tape can be divided into three types: half-wave type, half-wave and half-incremental curve type and incremental curve type. Under the DC overcurrent, there are two quenching resistance types of the HTS tape: half-wave type and incremental curve type.
(2) The quench resistance type is closely related to the current and the temperature of the HTS tape. When the current exceeds the critical current and the temperature is lower than the critical temperature of the HTS tape, the HTS tape is in partial resistance state. When the temperature exceeds the critical temperature, the HTS tape enters normal state.
(3) In the normal state, the relationship between quenching resistance R and joule heat Q of HTS tape and coil is independent of current waveform. In essence, there is a fixed correspondence between R and Q. Therefore, R-Q curves are consistent under AC and DC overcurrent.
Based on these, the R-Q curve method can be used to estimate and predict the quench resistance and temperature change of HTS tapes in completely quenching stage. Compared with FEM, this method can ensure simulation accuracy, greatly shorten calculation time and improve simulation efficiency. It can be used for simulation modeling, prototype design and quench resistance estimation of resistive type SFCL in complex power systems. Furthermore, this method can realize the equivalent calculation of quenching resistance under AC and DC overcurrent, which can provide theoretical support for large-scale DC current limiters to be tested by AC equivalent experiment.

Author Contributions

Conceptualization, S.L. and L.R.; Data curation, J.S.; Formal analysis, S.L. and S.Y.; Funding acquisition, Y.X. and M.S.; Investigation, X.T., Z.L., G.C. and Z.C.; Methodology, S.L. and L.R.; Project administration, T.M. and Y.X.; Resources, T.M. and L.X.; Software, S.L.; Supervision, Y.T.; Validation, X.T., Z.L., G.C., Z.C. and J.S.; Writing—original draft, S.L.; and Writing—review and editing, L.R.

Funding

This research was funded by the Project of the National Natural Science Foundation of China under grant number 51707074 and National Key Research and Development Program 2017YFB0902304.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

BSCCOBismuth Strontium Calcium Copper Oxide
YBCOYttrium Barium Copper Oxide
1GFirst Generation
2GSecond Generation
3DThree Dimensional
HTSHigh Temperature Superconducting
R-SFCLResistance type Superconducting Fault Current Limiter
DCDirect Current
ACAlternating Current
VSCVoltage Source Converter
MMCModular Multilevel Converter
HVDCHigh Voltage Direct Current
FEMFinite element method
PDEPartial Differential Equation
RLCResistance (R), Inductance (L) and Capacitance (C)

Symbols

UVoltage (V)QJoule heat (J)
ICurrent (A)PPower (W)
RResistance (Ω)EElectric field intensity (V/m)
LInductance (H)JCurrent density (A/m2)
CCapacitance (F)rResistivity (Ω·m)
tTime (s)qsVolume power density (W/m)
mMass (kg)kHeat transfer coefficient (W/(m2·K))
ρMass density (kg/m3)μPermeability (H/m)
cSpecific heat capacity (J/(kg·K))HMagnetic field intensity (A/m)
qConduction heat flux (W/m2)BMagnetic induction intensity (T)
TTemperature (K)

Appendix A

The following are the material properties required in the simulation, which mainly come from [46,49,50]. After collection and collation, the properties of materials are shown in Figure A1 and Table A1.
Figure A1. Material properties: (a) specific heat capacity of each layer of material; (b) thermal conductivity of each layer of material; and (c) resistivity of each layer of material.
Figure A1. Material properties: (a) specific heat capacity of each layer of material; (b) thermal conductivity of each layer of material; and (c) resistivity of each layer of material.
Materials 12 02374 g0a1
Table A1. Mass density of each material in superconducting tape.
Table A1. Mass density of each material in superconducting tape.
Mass Density (kg/m3)
Cu8940YBCO5900
Ag10500Stainless steel7930
Substrate8910

References

  1. Malozemoff, A.P. Second-Generation High-Temperature Superconductor Wires for the Electric Power Grid. Annu. Rev. Mater. Res. 2012, 42, 373–397. [Google Scholar] [CrossRef]
  2. Rogalla, H.; Kes, P.H. (Eds.) 100 Years of Superconductivity; CRC Press: London, UK, 2011. [Google Scholar]
  3. Ronald, M.; Alexis, P.; David, C. Superconducting Materials for Large Scale Applications. IEEE Trans. Appl. Supercond. 2004, 92, 1639–1654. [Google Scholar]
  4. Malozemoff, A.P.; Verebelyi, D.T. HTS Wire: Status and prospects. Phys. C. Superconduct. 2003, 386, 424–430. [Google Scholar] [CrossRef]
  5. Chen, Y.; Bian, W.; Huang, W. High critical current density of YBa2Cu3O7-x superconducting films prepared through a DUV-assisted solution deposition process. Sci. Rep. 2016, 6, 38257. [Google Scholar] [CrossRef]
  6. Schoop, U.; Rupich, M.W.; Thieme, C. Second generation HTS wire based on RABiTS substrates and MOD YBCO. IEEE Trans. Appl. Supercond. 2005, 15, 2611–2616. [Google Scholar] [CrossRef]
  7. Rupich, M.W.; Schoop, U.; Verebelyi, D.T. The Development of Second Generation HTS Wire at American Superconductor. IEEE Trans. Appl. Supercond. 2007, 17, 3379–3382. [Google Scholar] [CrossRef]
  8. Pascal, P.T.; Badel, A.; Auran, G. Superconducting Fault Current Limiter for Ship Grid Simulation and Demonstration. IEEE Trans. Appl. Supercond. 2017, 27, 5601705. [Google Scholar] [CrossRef]
  9. Yang, K.; Yang, Y.; Junaid, M. Direct-Current Vacuum Circuit Breaker with Superconducting Fault-Current Limiter. IEEE Trans. Appl. Supercond. 2018, 28, 560010. [Google Scholar] [CrossRef]
  10. Chen, L.; Zhang, X.; Qin, Y. Application and Design of a Resistive-Type Superconducting Fault Current Limiter for Efficient Protection of a DC Microgrid. IEEE Trans. Appl. Supercond. 2019, 29, 5600607. [Google Scholar] [CrossRef]
  11. He, H.; Chen, L.; Yin, T. Application of a SFCL for Fault Ride-Through Capability Enhancement of DG in a Microgrid System and Relay Protection Coordination. IEEE Trans. Appl. Supercond. 2016, 26, 5603608. [Google Scholar] [CrossRef]
  12. Li, B.; He, J. Studies on the Application of R-SFCL in the VSC-Based DC Distribution System. IEEE Trans. Appl. Supercond. 2016, 26, 5601005. [Google Scholar] [CrossRef]
  13. Chen, L.; Chen, H.; Shu, Z. Comparison of Inductive and Resistive SFCL to Robustness Improvement of a VSC-HVDC System with Wind Plants against DC Fault. IEEE Trans. Appl. Supercond. 2016, 26, 5603508. [Google Scholar] [CrossRef]
  14. Xiao, L.; Dai, S.; Lin, L. HTS Power Technology for Future DC Power Grid. IEEE Trans. Appl. Supercond. 2013, 23, 5401506. [Google Scholar] [CrossRef]
  15. Zhou, Y.; Song, Q.; Guo, F. Quench developing process of HTS tapes under sinusoidal over-currents. IEEE Trans. Appl. Supercond. 2005, 15, 1651–1654. [Google Scholar] [CrossRef]
  16. Baldan, C.A.; Lamas, J.S.; Shigue, C.Y. Test of a Modular Fault Current Limiter for 220 V Line Using YBCO Coated Conductor Tapes with Shunt Protection. IEEE Trans. Appl. Supercond. 2011, 21, 1242–1245. [Google Scholar] [CrossRef]
  17. Liu, X.; Wen, J.; Zeng, W. Quenching Characteristics of Different Types of Superconducting Fault Current Limiting Modules. IEEE Trans. Appl. Supercond. 2015, 25, 5601205. [Google Scholar] [CrossRef]
  18. Zhang, Z.; Yang, J.; Qiu, Q. Research on resistance characteristics of YBCO tape under short-time DC large current impact. Cryogenics 2017, 84, 53–59. [Google Scholar] [CrossRef]
  19. Xiang, B.; Zhang, L.; Tan, Y. DC current withstanding characteristics of superconductor. In Proceedings of the 3rd International Conference on Electric Power Equipment—Switching Technology, Busan, Korea, 25–28 October 2015; pp. 217–220. [Google Scholar]
  20. Jiang, Z.; Wang, Y.; Dai, S. Influence of Insulation on Quench and Recovery of YBCO Tape under DC Impact. IEEE Trans. Appl. Supercond. 2018, 29, 7700305. [Google Scholar] [CrossRef]
  21. Xiang, B.; Junaid, M.; Gao, L. Influencing Factors on Quench and Recovery of YBCO Tapes for DC Superconducting Fault Current Limiter. IEEE Trans. Appl. Supercond. 2018, 29, 5600806. [Google Scholar] [CrossRef]
  22. Xiang, B.; Junaid, M.; Gao, L. Effects of Short Circuit Currents on Quench and Recovery Properties of YBCO Tapes for DC SFCL. IEEE Trans. Appl. Supercond. 2018, 29, 5600706. [Google Scholar] [CrossRef]
  23. Badel, A.; Antognazza, L.; Decroux, M. Hybrid Model of Quench Propagation in Coated Conductors Applied to Fault Current Limiter Design. IEEE Trans. Appl. Supercond. 2013, 23, 5603705. [Google Scholar] [CrossRef]
  24. Tong, Y.; Guan, M.; Wang, X. Theoretical estimation of quench occurrence and propagation based on generalized thermoelasticity for LTS/HTS tapes triggered by a spot heater. Supercond. Sci. Technol. 2017, 30, 045002. [Google Scholar] [CrossRef]
  25. Colangelo, D.; Dutoit, B. Analysis of the influence of the normal zone propagation velocity on the design of resistive fault current limiters. Supercond. Sci. Technol. 2014, 27, 124005. [Google Scholar] [CrossRef]
  26. Núñez-Chico, A.B.; Martínez, E.; Angurel, L.A.; Navarro, R. Enhanced quench propagation in 2G-HTS coils co-wound with stainless steel or anodised aluminium tapes. Supercond. Sci. Technol. 2016, 29, 085012. [Google Scholar] [CrossRef] [Green Version]
  27. Hong, Z.; Sheng, J.; Yao, L.; Gu, J.; Jin, Z. The Structure, Performance and Recovery Time of a 10 kV Resistive Type Superconducting Fault Current Limiter. IEEE Trans. Appl. Supercond. 2012, 23, 5601304. [Google Scholar] [CrossRef]
  28. Ge, H.; Yang, K.; Junaid, M. A quenching recovery time test method for resistive type superconducting fault current limiters used in DC circuit. In Proceedings of the 4th International Conference on Electric Power Equipment—Switching Technology, Xi’an, China, 22–25 October 2017; pp. 393–396. [Google Scholar]
  29. Yang, D.G.; Song, J.B.; Choi, Y.H. Quench and Recovery Characteristics of the Zr-Doped (Gd,Y) BCO Coated Conductor Pancake Coils Insulated With Copper and Kapton Tapes. IEEE Trans. Appl. Supercond. 2011, 21, 2415–2419. [Google Scholar] [CrossRef]
  30. Lim, S.H.; Lim, S.T. Current Limiting and Recovery Characteristics of a Trigger-Type SFCL Using Double Quench. IEEE Trans. Appl. Supercond. 2018, 28, 5601305. [Google Scholar] [CrossRef]
  31. Schwarz, M.; Schacherer, C.; Weiss, K.P. Thermodynamic behaviour of a coated conductor for currents above Ic. Supercond. Sci. Technol. 2008, 21, 054008. [Google Scholar] [CrossRef]
  32. Bagrets, N.; Otten, S.; Weiss, K.P. Thermal and mechanical properties of advanced impregnation materials for HTS cables and coils. IOP Conf. Ser. Mater. Sci. Eng. 2015, 102, 012021. [Google Scholar] [CrossRef]
  33. Bae, J.H.; Eom, B.Y.; Sim, K.D. Minimum Quench Energy Characteristic of YBCO Coated Conductor with Different Stabilizer Thickness. IEEE Trans. Appl. Supercond. 2013, 23, 4600404. [Google Scholar]
  34. Falorio, I.; Young, E.A.; Yang, Y. Quench Characteristic and Minimum Quench Energy of 2G YBCO Tapes. IEEE Trans. Appl. Supercond. 2015, 25, 6605505. [Google Scholar] [CrossRef]
  35. Du, H.I.; Kim, Y.J.; Lee, D.H. Study on Maximum Operating Condition of Resistive Type SFCL Using YBCO Coated Conductor. IEEE Trans. Appl. Supercond. 2010, 20, 1238–1241. [Google Scholar]
  36. Levin, G.A.; Jones, W.A.; Novak, K.A. The effects of superconductor-stabilizer interfacial resistance on quenching of a pancake coil made out of coated conductor. Supercond. Sci. Technol. 2011, 24, 035015. [Google Scholar] [CrossRef]
  37. Kwon, N.Y.; Kim, H.S.; Kim, K.L. The Effects of a Stabilizer Thickness of the YBCO Coated Conductor (CC) on the Quench/Recovery Characteristics. IEEE Trans. Appl. Supercond. 2010, 20, 1246–1249. [Google Scholar] [CrossRef]
  38. Du, H. Evaluation on Resistance Tendency and Recovery Characteristics of 2G Wire with Insulation Layer. IEEE Trans. Appl. Supercond. 2013, 23, 6602004. [Google Scholar]
  39. The Shanghai Superconductor Website 2019. Available online: http://www.shsctec.com/en/index (accessed on 20 March 2019).
  40. Yang, J.; Fletcher, J.; O’Reilly, J. Short-Circuit and Ground Fault Analyses and Location in VSC-Based DC Network Cables. IEEE Trans. Ind. Electron. 2012, 59, 3827–3837. [Google Scholar] [CrossRef]
  41. Li, C.; Zhao, C.; Xu, J. A Pole-to-Pole Short-Circuit Fault Current Calculation Method for DC Grids. IEEE Trans. Power Syst. 2017, 32, 4943–4953. [Google Scholar] [CrossRef]
  42. Hong, Z.; Campbell, A.M.; Coombs, T.A. Numerical solution of critical state in superconductivity by finite element software. Supercond. Sci. Technol. 2006, 19, 1246–1252. [Google Scholar] [CrossRef]
  43. Duron, J.; Grilli, F.; Antognazza, L. Finite-element modelling of YBCO fault current limiter with temperature dependent parameters. Supercond. Sci. Technol. 2007, 20, 338–344. [Google Scholar] [CrossRef]
  44. Stavrev, S.; Grilli, F.; Dutoit, B. Comparison of numerical methods for modeling of superconductors. IEEE Trans. Magn. 2002, 38, 849–852. [Google Scholar] [CrossRef]
  45. Curras, S.R.; Vina, J.; Ruibal, M. Normal-state resistivity versus critical current in YBa2Cu3O7−δ thin films at high current densities. Phy. C. Superconduct. 2002, 372, 1095–1098. [Google Scholar] [CrossRef]
  46. Roy, F.; Dutoit, B.; Grilli, F. Magneto-Thermal Modeling of Second-Generation HTS for Resistive Fault Current Limiter Design Purposes. IEEE Trans. Appl. Supercond. 2008, 18, 29–35. [Google Scholar] [CrossRef]
  47. Liang, F.; Yuan, W.; Baldan, C. Modeling and Experiment of the Current Limiting Performance of a Resistive Superconducting Fault Current Limiter in the Experimental System. J. Supercond. Novel Magn. 2015, 28, 2669–2681. [Google Scholar] [CrossRef] [Green Version]
  48. Frost, W.; Harper, W.L. Heat Transfer at Low Temperatures; Plenum: New York, NY, USA, 1975; Chapter 4. [Google Scholar]
  49. De Sousa, W.; Näckel, O.; Noe, M. Transient Simulations of an Air-Coil SFCL. IEEE Trans. Appl. Supercond. 2014, 24, 5601807. [Google Scholar] [CrossRef]
  50. Kalsi, S.S. Applications of High Temperature Superconductors to Electric Power Equipment; John Wiley & Sons: Hoboken, NJ, USA, 2011; pp. 137–138. [Google Scholar]
Figure 1. A 2G HTS tape produced by Shanghai Superconductor Company: (a) 3D structural diagram; (b) sectional view; and (c) the sample prepared for experiment.
Figure 1. A 2G HTS tape produced by Shanghai Superconductor Company: (a) 3D structural diagram; (b) sectional view; and (c) the sample prepared for experiment.
Materials 12 02374 g001
Figure 2. Two type basic waveforms in the experiment: (a) AC overcurrent; and (b) DC overcurrent.
Figure 2. Two type basic waveforms in the experiment: (a) AC overcurrent; and (b) DC overcurrent.
Materials 12 02374 g002
Figure 3. AC overcurrent test platform: (a) circuit diagram of test platform; and (b) physical map of test platform.
Figure 3. AC overcurrent test platform: (a) circuit diagram of test platform; and (b) physical map of test platform.
Materials 12 02374 g003
Figure 4. DC overcurrent test platform: (a) circuit diagram of test platform; and (b) physical map of test platform.
Figure 4. DC overcurrent test platform: (a) circuit diagram of test platform; and (b) physical map of test platform.
Materials 12 02374 g004
Figure 5. The corresponding experimental results of AC overcurrent test with different amplitude of AC overcurrent.
Figure 5. The corresponding experimental results of AC overcurrent test with different amplitude of AC overcurrent.
Materials 12 02374 g005
Figure 6. The corresponding experimental results DC overcurrent test with different amplitude of DC overcurrent.
Figure 6. The corresponding experimental results DC overcurrent test with different amplitude of DC overcurrent.
Materials 12 02374 g006
Figure 7. The R-t curves of superconducting tape under different overcurrent: (a) the R-t curves under AC overcurrent; and (b) the R-t curves under DC overcurrent.
Figure 7. The R-t curves of superconducting tape under different overcurrent: (a) the R-t curves under AC overcurrent; and (b) the R-t curves under DC overcurrent.
Materials 12 02374 g007
Figure 8. The R-Q curves of superconducting tape under AC and DC overcurrent: (a) the comparison diagram of R-Q; and (b) partial enlarged drawing of the R-Q curves.
Figure 8. The R-Q curves of superconducting tape under AC and DC overcurrent: (a) the comparison diagram of R-Q; and (b) partial enlarged drawing of the R-Q curves.
Materials 12 02374 g008
Figure 9. The basic structure of the simulation model: (a) geometric structure of superconducting tape; and (b) the coupling structure of electromagnetic thermal.
Figure 9. The basic structure of the simulation model: (a) geometric structure of superconducting tape; and (b) the coupling structure of electromagnetic thermal.
Materials 12 02374 g009
Figure 10. Heat transfer coefficient curve of liquid nitrogen.
Figure 10. Heat transfer coefficient curve of liquid nitrogen.
Materials 12 02374 g010
Figure 11. Comparison of experimental and simulation results.
Figure 11. Comparison of experimental and simulation results.
Materials 12 02374 g011
Figure 12. The simulation analysis with AC 860 A overcurrent: (a) the quenching resistance characteristics; and (b) the thermal characteristics and current distribution characteristics.
Figure 12. The simulation analysis with AC 860 A overcurrent: (a) the quenching resistance characteristics; and (b) the thermal characteristics and current distribution characteristics.
Materials 12 02374 g012
Figure 13. The simulation analysis with AC 1200 A overcurrent: (a) the quenching resistance characteristics; and (b) the thermal characteristics and current distribution characteristics.
Figure 13. The simulation analysis with AC 1200 A overcurrent: (a) the quenching resistance characteristics; and (b) the thermal characteristics and current distribution characteristics.
Materials 12 02374 g013
Figure 14. The simulation analysis with AC 1780 A overcurrent: (a) the quenching resistance characteristics; and (b) the thermal characteristics and current distribution characteristics.
Figure 14. The simulation analysis with AC 1780 A overcurrent: (a) the quenching resistance characteristics; and (b) the thermal characteristics and current distribution characteristics.
Materials 12 02374 g014
Figure 15. The simulation analysis with DC 910 A overcurrent: (a) the quenching resistance characteristics; and (b) the thermal characteristics and current distribution characteristics.
Figure 15. The simulation analysis with DC 910 A overcurrent: (a) the quenching resistance characteristics; and (b) the thermal characteristics and current distribution characteristics.
Materials 12 02374 g015
Figure 16. The simulation analysis with DC 2020 A overcurrent: (a) the quenching resistance characteristics; and (b) the thermal characteristics and current distribution characteristics.
Figure 16. The simulation analysis with DC 2020 A overcurrent: (a) the quenching resistance characteristics; and (b) the thermal characteristics and current distribution characteristics.
Materials 12 02374 g016
Figure 17. The simulation results of R-Q curves.
Figure 17. The simulation results of R-Q curves.
Materials 12 02374 g017
Figure 18. The simulation analysis of heat transfer process of 1 m superconducting tape under AC and DC overcurrent.
Figure 18. The simulation analysis of heat transfer process of 1 m superconducting tape under AC and DC overcurrent.
Materials 12 02374 g018
Figure 19. Equivalent material properties: (a) the equivalent quenching resistance curve; and (b) the equivalent specific heat capacity.
Figure 19. Equivalent material properties: (a) the equivalent quenching resistance curve; and (b) the equivalent specific heat capacity.
Materials 12 02374 g019
Figure 20. The comparison diagram of R-Q curves from different results.
Figure 20. The comparison diagram of R-Q curves from different results.
Materials 12 02374 g020
Figure 21. Different stage distribution of superconducting tape immersed in liquid nitrogen.
Figure 21. Different stage distribution of superconducting tape immersed in liquid nitrogen.
Materials 12 02374 g021
Figure 22. The calculation process of R-Q curve method.
Figure 22. The calculation process of R-Q curve method.
Materials 12 02374 g022
Figure 23. R-Q curve required in simulation calculation.
Figure 23. R-Q curve required in simulation calculation.
Materials 12 02374 g023
Figure 24. The quenching resistance comparison of FEM method, R-Q curve method and experiment.
Figure 24. The quenching resistance comparison of FEM method, R-Q curve method and experiment.
Materials 12 02374 g024
Figure 25. The temperature comparison of FEM and R-Q curve method.
Figure 25. The temperature comparison of FEM and R-Q curve method.
Materials 12 02374 g025
Figure 26. The model of R-Q curve method improved in PSCAD.
Figure 26. The model of R-Q curve method improved in PSCAD.
Materials 12 02374 g026
Figure 27. The schematic diagram of non-inductive superconducting coil: (a) physical drawing; and (b) structure drawing.
Figure 27. The schematic diagram of non-inductive superconducting coil: (a) physical drawing; and (b) structure drawing.
Materials 12 02374 g027
Figure 28. Comparison between the result of R-Q curve method and experimental result of 136 m non-inductive superconducting coil: (a) quenching resistance; and (b) temperature estimation.
Figure 28. Comparison between the result of R-Q curve method and experimental result of 136 m non-inductive superconducting coil: (a) quenching resistance; and (b) temperature estimation.
Materials 12 02374 g028
Table 1. Material parameters of 2G HTS tape used in experiment.
Table 1. Material parameters of 2G HTS tape used in experiment.
ParametersValue
Width of tape12 mm
Thickness of stainless steel layer80 μm
Thickness of Cu layer3 μm
Thickness of Ag layer1.5 μm
Thickness of YBCO layer1 μm
Thickness of Hastelloy substrate layer50 μm
InsulationNo
Unit resistance at room temperature58.42 mΩ/m
Self-field critical current Ic at 77 K500 A
Effective measurement length of Sample 10 cm
Table 2. Basic parameters of AC overcurrent test platform.
Table 2. Basic parameters of AC overcurrent test platform.
AC Test PlatformDC Test Platform
ParametersValueParametersValue
Input voltage220 VOutput voltage UC1000–2200 V
Frequency50 HzCapacitance C08 mF
Output voltage4–6 VInductance L05 mH
Output current0–3 kAResistance R00.2 Ω/0.5 Ω
Duration50–80 msDuration20–23 ms
Table 3. The parameters setting of AC and DC overcurrent tests.
Table 3. The parameters setting of AC and DC overcurrent tests.
NumberIPmax (A)IPmax/ICDuration (ms)
AC/DCACDCACDCACDC
No. 18609101.721.825022
No. 2106012602.122.526023
No. 3120014002.42.86020
No. 4142015802.843.167020
No. 5158016203.163.248022
No. 6162018603.243.726022
No. 7178019803.563.966022
No. 8192020203.844.047522
Table 4. The summary of AC overcurrent experiments.
Table 4. The summary of AC overcurrent experiments.
Resistance TypeWaveformResultsCharacteristicsState
Half-wave type Materials 12 02374 i001AC No. 1Recovery after quenchingPartial resistive
Half-wave and half-incremental curve type Materials 12 02374 i002AC No. 2
AC No. 3
Initial stage: Recovery after quenching;
Later stage: continuously quenching
Transition:
Partial resistive state transitions to normal state
Incremental curve type Materials 12 02374 i003AC No. 4
AC No. 5
AC No. 6
AC No. 7
AC No. 8
Continuously quenching to normal stateNormal
Table 5. The summary of DC overcurrent experiments.
Table 5. The summary of DC overcurrent experiments.
Resistance TypeWaveformResultsCharacteristicsState
Half-wave type Materials 12 02374 i004DC No. 1Recovery after quenchingPartial resistive
Incremental curve type Materials 12 02374 i005DC No. 2
DC No. 3
DC No. 4
DC No. 5
DC No. 6
DC No. 7
DC No. 8
Continuously quenching to normal stateNormal
Table 6. Simulation parameters of YBCO resistivity.
Table 6. Simulation parameters of YBCO resistivity.
ParametersValueParametersValue
r01 × 10−14 Ω·cmn12.8
rnorm2.5 μΩ·cmn222
E00.5 V/cmk1.92
Tc90 Kα1.5
Jc01.9 × 1010 A/m2Tref77 K
Table 7. The common quenching characteristics of the superconducting tape.
Table 7. The common quenching characteristics of the superconducting tape.
StateBoundaryCurrent Distribution CharacteristicsYBCO ResistivityR-Q Consistency
Superconducting stateI < IC (T)Only flow through superconducting layer0/
QuenchingPartial resistive state:I > IC (T)
T < TC
Initial stage:
IYBCO > ICu > IAg > ISt > ISub
Late stage:
ICu > IYBCO > IAg > ISt > ISub
( r 1 + r 2 + r 0 ) r n o r m r 0 + r 1 + r 2 + r n o r m No
Normal stateTC < TICu > IAg > ISt > ISub > IYBCO r n o r m Yes
Table 8. Simulation parameters setting of R-Q curve method.
Table 8. Simulation parameters setting of R-Q curve method.
ParametersValue
Critical current, IC500 A
Initial quenching resistance, R00.001 Ω
Time step, Δt5 × 10−5 s
Length of superconducting tape1 m
OvercurrentAC 1780 A, AC 1920 A,
DC 1860 A, DC 2020 A,
Duration60 ms, 80 ms, 25 ms, 25 ms
Table 9. The computation time of R-Q curve method and FEM.
Table 9. The computation time of R-Q curve method and FEM.
OvercurrentSoftwareMethodComputation Time
AC 1780A, 60 msPSCADR-Q<3 s
ComsolFEM3.78 h
AC 1920A, 75 msPSCADR-Q< 3 s
ComsolFEM3.85 h
DC 1860A, 22 msPSCADR-Q<3 s
ComsolFEM1.57 h
DC 2020A, 22 msPSCADR-Q<3 s
ComsolFEM1.58 h

Share and Cite

MDPI and ACS Style

Liang, S.; Ren, L.; Ma, T.; Xu, Y.; Tang, Y.; Tan, X.; Li, Z.; Chen, G.; Yan, S.; Cao, Z.; et al. Study on Quenching Characteristics and Resistance Equivalent Estimation Method of Second-Generation High Temperature Superconducting Tape under Different Overcurrent. Materials 2019, 12, 2374. https://doi.org/10.3390/ma12152374

AMA Style

Liang S, Ren L, Ma T, Xu Y, Tang Y, Tan X, Li Z, Chen G, Yan S, Cao Z, et al. Study on Quenching Characteristics and Resistance Equivalent Estimation Method of Second-Generation High Temperature Superconducting Tape under Different Overcurrent. Materials. 2019; 12(15):2374. https://doi.org/10.3390/ma12152374

Chicago/Turabian Style

Liang, Siyuan, Li Ren, Tao Ma, Ying Xu, Yuejin Tang, Xiangyu Tan, Zheng Li, Guilun Chen, Sinian Yan, Zhiwei Cao, and et al. 2019. "Study on Quenching Characteristics and Resistance Equivalent Estimation Method of Second-Generation High Temperature Superconducting Tape under Different Overcurrent" Materials 12, no. 15: 2374. https://doi.org/10.3390/ma12152374

APA Style

Liang, S., Ren, L., Ma, T., Xu, Y., Tang, Y., Tan, X., Li, Z., Chen, G., Yan, S., Cao, Z., Shi, J., Xiao, L., & Song, M. (2019). Study on Quenching Characteristics and Resistance Equivalent Estimation Method of Second-Generation High Temperature Superconducting Tape under Different Overcurrent. Materials, 12(15), 2374. https://doi.org/10.3390/ma12152374

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

Article Metrics

Back to TopTop