Next Article in Journal
A Photovoltaic Power Curtailment Method for Operation on Both Sides of the Power-Voltage Curve
Previous Article in Journal
Value Co-Creation Behavior in Green Supply Chains: An Empirical Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thermo-Hydraulic Analysis of a Tri-Axial High-Temperature Superconducting Power Cable with Respect to Installation Site Geography

1
Department of Mechanical Engineering, Changwon National University, Changwon 641-773, Korea
2
Korea Electric Power Corporation Research Institute, Daejeon 34056, Korea
*
Author to whom correspondence should be addressed.
Energies 2020, 13(15), 3898; https://doi.org/10.3390/en13153898
Submission received: 8 June 2020 / Revised: 27 July 2020 / Accepted: 30 July 2020 / Published: 30 July 2020

Abstract

:
Various high-temperature superconducting (HTS) power cables are being developed or are ready for commercial operation to help energy suppliers meet the growing power demand in urban areas. Recently, triaxial HTS power cables have been developed by Korea Electric Power Corporation (KEPCO) and LS Cable & System. Further, KEPCO has been planning to install a 2 km long 23 kV/60 MVA triaxial HTS power cable to connect the Munsan and Seonyu substations and increase the stability of the power grid. The HTS power cables should be cooled down to a cryogenic temperature near 77 K. A thermo-hydraulic analysis of the cooling system considering the geographical characteristics of the installation site is essential for long-distance sections. This paper describes the thermo-hydraulic analysis of the triaxial HTS power cable to determine the proper mass flow rates of subcooled liquid nitrogen that meet the operating temperature and pressure of the cable for four configurations of cooling systems: (1) a single cooling system with an external return path, (2) a dual cooling system with an external return path, (3) a single cooling system with an internal return path, and (4) a dual cooling system with internal return path. Since the flow characteristics in a corrugated cable cryostat differ significantly from those in a typical annular tube, a computational fluid dynamics (CFD) analysis was conducted to estimate the pressure drop along the cable cryostat. With the CFD analysis and given heat loads in the cable, the temperature and the pressure variations along the cable were calculated and their pros and cons were compared for each configuration of the cooling system. This thermo-hydraulic analysis will be referenced in the actual installation of the HTS power cable between the Munsan and Seonyu substations.

Graphical Abstract

1. Introduction

High-temperature superconducting (HTS) power cables with low electric loss and large power capacities are being developed or are ready for commercial operation to meet the increasing energy demand in downtown areas [1,2,3,4,5]. Korea Electric power Corporation (KEPCO) and LS Cable & System commercialized a 23 kV, 50 MVA superconducting cable in the Singal-Heungdeok section that is approximately 1 km long. The Jeju smart-grid demonstration center completed the testing of a 1 km, 154 kV, 600 MVA three-phase superconducting cable [6].
Recently, they developed a triaxial HTS cable with enhanced current characteristics, compared to a 3-core-in-1 cryostat type [7], and are now preparing its commercialization connecting the Munsan and Seounyu substations to increase the stability of the power grid [5,8].
The HTS power cable is cooled by the circulation of subcooled liquid nitrogen (LN2) at approximately 77 K in order to remove various heat losses such as AC loss, heat penetration into the cable, and the heat load at terminations. Increasing the mass flow rate of LN2 can reduce the temperature difference between the cable inlet and outlet, reducing the maximum temperature of the cable, but doing so increases the pressure across the cable. That is, the required mass flow rate of LN2 is determined to satisfy the operating temperature while the operating pressure is within the allowable limit.
There have been studies on computational fluid dynamic (CFD) analysis to predict the pressure drop due to the flow of liquid nitrogen in the cross-sectional shape of various types of superconducting cables [9,10,11]. To install a long HTS power cable over several kilometers, it is very important to determine the operating conditions for the “dedicated” HTS power cable in consideration of the heat transfer characteristics of the cable along with the flow characteristics associated with the actual installation environment.
In order to determine appropriate mass flow rate of LN2, it is necessary to investigate the thermo-hydraulic characteristics of the cable cryostat, and the details of the installation environments, such as vertical elevation, the locations and the number of joint boxes and terminals, should also be considered. The Munsan and Seonyu substations are 2.064 km apart from each other; along the installation path there are two vertical sections up to 27 m, four terminals, and four joint boxes.
This paper describes the thermo-hydraulic analysis of the 23 kV/60 MVA tri-axial HTS power cable to determine proper mass flow rates of subcooled liquid nitrogen that meet the operating temperature and pressure of the cable with four configurations of a cooling system: (1) a single cooling system with an external return path, (2) a dual cooling system with an external return path, (3) a single cooling system with an internal return path, and (4) a dual cooling system with internal return path. Since the flow characteristics in a corrugated cable cryostat differ significantly from that of a typical annular tube, a computational fluid dynamics (CFD) analysis was performed to estimate the pressure drop along the cable cryostat [12,13,14]. With the CFD analysis and given heat loads in the cable, the temperature and the pressure variations along the cable were calculated and their pros and cons were compared for each configuration of a cooling system. In addition, an expanded 1-D thermo-hydraulic network analysis was conducted to consider a radial heat transfer in a cable core that has an internal return path.

2. Structure of Tri-Axial HTS Power Cable

Table 1 lists the major specifications of triaxial HTS power cable. Brass laminated coated conductors are used as HTS wires to transport the current and PPLP (Polypropylene Laminated paper) is used for electric insulation in cryogenic environment.
The cable cryostat for the triaxial HTS power cable has a spiral corrugation geometry that makes it easy to bend the cable during installation; the cable core sags down to the cryostat due to gravity, not concentric shape. As shown in Figure 1, two types of the cable structure for flow path were considered: (a) the external return path type, and (b) the internal return path type. For an external return path, LN2 flows outside the cable core and then returns to the cooling system through a separate external return tube. In this case, the flow area of the external return tube is large enough that the pressure drop in the return tube is very small; however, there is additional heat loss in the return tube as in Figure 1a. The cable core consists of three HTS layers for each phase to conduct current, a central copper former to support concentric shape, an insulation layer between each phase of the cable for dielectric strength, and copper shield layers to transport fault current. For an internal return path, LN2 flows inside of the cable core and returns back through the outside of the cable core without requiring a separate external return tube [10,15,16,17]. This prevents additional heat loss of the external return tube and cools the wider surface area of the core. However, the pressure drop tends to increase due to the reduced flow area as in Figure 1b. In lieu of the copper former, the inner flow tube can serve as a concentric support [18].

3. Geographical Characteristics of Installation Path

The Munsan and Seonyu substations are 2.064 km apart from each other and there is also a platform between them. The platform functions as a switching station for terminal distribution lines instead of a convectional substation. The total distance from Munsan to the platform is 983 m, with two terminals, two joint boxes, and a 25.5 m vertical section along the way. The distance from Seonyu to the platform is 1081 m long, with two terminals, two joint boxes, and a 27 m vertical section along the way. Figure 2 illustrates the layout of the triaxial HTS power cable between the Munsan and Seonyu substations.
Δ P = K b m ˙ 2 2 ρ A 2
Figure 3a shows the reference geometry for a cable snake and bending, which was applied for the 154 kV, 600 MVA, 3-phase HTS power cable at the Jeju smart-grid demonstration center. Since a single-phase HTS power cable at the demonstration center has a similar shape to the triaxial HTS power cable, its installation characteristics such as bending and a cable snake were used as reference data. In the demonstration center, we observed a cable snake depth of 15 cm in the middle section of the HTS power cable for every 3.5 m. The bending radius of the cable was 4.5 to 7 m.
By modeling the cable snake as a sin wave as shown in Figure 3b, it was replaced with two 20.16-deg bending section with radius of 5 m. Along the cable installation path of 2.064 km, there are 1180 bends with a 20.16 bend angle. Figure 4 illustrates the minor loss coefficient for a bent round pipe. In this graph, the minor loss coefficient due to the bend becomes smaller as the ratio of the bending radius to the tube diameter increases; if this ratio is greater than 6, the minor loss coefficient only depends on the bending angle. In the actual HTS power cable installation environment, the ratio of bending radius to the tube diameter is much larger than 6. Since the HTS power cable will be installed along the existing underground installation tunnel, we considered eight 90 degree bending sections and one each of 26, 66 and 14 degree bending sections along the installation path referring to the layout of the underground cable installation site provided by KEPCO. Table 2 lists the calculation results of the minor loss coefficient on the cable snake and at the bend as obtained from Figure 4. The additional pressure drop due to the minor loss coefficient is calculated by Equation (1). The total additional pressure drop due to the minor losses from bending of the cable is 1439 Pa at 70 K, 7 bar and 0.5 kg/s, and it seems to be negligible considering the estimated total pressure drop of the cable.

4. Single and Dual Cooling System with External Return Paths

After circulating through the cable and the return tube, the temperature of LN2 increases and the pressure decreases due to heat loads and flow friction. Then, the LN2 is re-cooled using a cryogenic cooling system and pumped up to its original pressure using a cryogenic circulation pump [10]. Since the total length of cable is about 2 km, two types of cryogenic cooling system were proposed considering the operating conditions and facility installation.
A cooling system can only be installed at the substations with sufficient installation space, and the options for the installation are described in Figure 5. Figure 5a,b show the configurations of cooling systems with an external return paths for a single cooling system at Munsan and dual cooling system, respectively. For the single cooling system, the initial installation cost is low; however, the estimated pressure drop is large, and each substation cannot be operated separately. If a fault occurred at one substation, both substations would have to be stopped for maintenance. In contrast, the dual cooling system incurs a high initial installation cost as it consists of two cooling systems; however, the pressure drop is relatively small, and it is easier to maintain because each substation can be operated separately.

5. Temperature Rise and Pressure Drop Analysis Model

5.1. Allowable Operating Conditions

Considering the optimal electrical performance of HTS cables, the maximum temperature of the cooling system was set to 75.5 K and the minimum temperature to 67 K to prevent the freezing of LN2. The operating pressure of the cooling system was 13.5 bar at maximum due to the allowable pressure of the cable cryostat and 5.5 bar at minimum to suppress bubble generation. Therefore, the maximum temperature difference and the maximum pressure differential were 8.5 K and 8 bar, respectively, between the inlet and the outlet. Besides the pressure drop by flow friction, the pressure variations along the vertical section should be considered to ensure the maximum pressure is within the allowable pressure. Table 3 is the limiting operating conditions of the cooling system for the HTS power cable.

5.2. Calculation of LN2 Temperature Variation

Prior to detailed flow and thermal analysis using CFD, the temperature variation and the pressure drop can be estimated with simple equations assuming the 1-D simplification of the cable. Since the LN2 in the cryostat flows in the steady state, the temperature difference ( Δ T ) between the inlet and outlet is related to the energy balance equation of
Δ T = Q ˙ m ˙ C p ,
where m ˙ and C p are the mass flow rate and the specific heat of LN2.
The total heat loss ( Q ˙ t o t a l ) of the cooling system is calculated by
Q ˙ t o t a l = Q ˙ c a b l e + Q ˙ t e r m i n a t i o n + Q ˙ j o i n t b o x + Q ˙ r e t u r n
where Q ˙ c a b l e is the estimated sum of heat losses at the cable including AC loss, dielectric loss and heat penetration into the cable cryostat, referring to the value provided by LS Cable, the cable manufacturer and the values of previous study based on FEM analysis [7]. Q ˙ t e r m i n a t i o n , Q ˙ j o i n t b o x and Q ˙ r e t u r n are the heat loss at terminations, joint boxes and heat penetration into the return tube, respectively. Referring to the previous analysis and experimental value, Table 4 lists the estimated heat losses of the triaxial HTS power cable [7].

5.3. Pressure Drop Model for the Cable Cryostat

With a given pressure loss coefficient, the pressure drop ( Δ P ) in pipes is calculated by
Δ P = f L D h m ˙ 2 2 ρ A 2 ,
where f , D h ,   ρ ,   A are the pressure loss coefficient, the hydraulic diameter of HTS cable and return path, the density of LN2 and the flow area, respectively. Since the cable cryostats have a corrugation shape, the averaged diameter ( D a o ) was used to calculate D h as in Figure 1.
The moody chart generally shows friction coefficient ( f ) from 0.008 to 0.1. However, the f is expected to be greater than 0.1 due to the shape of spiral corrugation. Then, the Colebrook–White equation of Equation (5) can be used to estimate the f with the given relative roughness.
1 f = 2 log ( ε / D h 3.7 + 2.51 R e d f ) ( R e d > 4000 ) ,
where ε ,   ε / D h and R e d are the surface roughness, the relative roughness and the Reynolds number. Table 5 is representative value of flow parameters for each case.
The wave height in a spiral corrugation can be simply assumed as a continuous surface roughness [20]. However, the flow stream is greatly affected by the large spiral pitch of the corrugation and the eccentric sagging of the cable core; further, it is necessary to define an equivalent relative roughness through CFD analysis to use Equations (4) and (5) for the pressure drop.
According to Equations (2) and (4), as the mass flow rate increases, the temperature difference decreases, and the pressure drop increases between the inlet and outlet of the HTS power cable.

5.4. CFD Analysis of Spiral Corrugation Cryostat

The mass flow rate was selected to be in the range 0.28–0.95 kg/s, which is the expected operating range of the cooling system to meet the allowable temperature and pressure limits. The detail process for the determination of appropriate mass flow rates is described in the next section. In the CFD analysis, the roughness of the cable cryostat and the cable core was assumed as smooth wall and 1 mm rough wall, respectively. The mass flow rate was selected to be in the range 0.28–0.95 kg/s, which is the expected operating range of the cooling system. The commonly used κ ε turbulence model was applied using ANSYS CFX ver. 19.2. The length of the model was set to 2 m to ignore the entrance region. Table 6 lists the boundary conditions for the CFD analysis.
Figure 6 shows the streamline of LN2 at a mass flow rate of 0.57 kg/s; it was possible to observe major straight flow and rotating flow along the spiral corrugation in each case.
Figure 7 is the pressure drop per unit length according to the mass flow rate in each flow path of the HTS power cable. The return flow of the internal return path had the greatest pressure drop per unit length. With the obtained pressure drop results, it is possible to determine the equivalent relative roughness for each flow path as in Table 7 by comparing the pressure drop with Equations (4) and (5). Then, the pressure drop can be simply calculated using Equations (4) and (5) instead of time consuming CFD analysis for the long HTS power cable.

6. Pressure and Temperature Analysis Results for External Return Path

The pressure drop for whole HTS power cable was calculated using Equation (4) with the equivalent relative roughness in Table 5. The temperature rise was calculated for the total heat loss in Figure 5 using Equation (2). Figure 8 shows the allowable mass flow rate of single cooling system and dual cooling systems with external return paths. For the single cooling system as in Figure 8a, the minimum mass flow rate was 0.57 kg/s to meet the temperature limit and the maximum mass flow rate was 0.66 kg/s to meet the pressure limit. In case of the dual cooling system on the Munsan side as in Figure 8b, the minimum mass flow rate was 0.28 kg/s to meet the temperature limit and the maximum mass flow rate was 0.98 kg/s to meet the pressure limit. In case of the dual cooling system on the Seonyu side as in Figure 8c, the minimum mass flow rate was 0.3 kg/s to meet the temperature limit and maximum mass flow rate was 0.93 kg/s to meet the pressure limit. For the dual cooling system, the length of the path is shorter than for the single cooling system, resulting in a wide mass flow rate range. However, the pressure of LN2 increases depending on the laying of the actual cable in the underground vertical section [21]. According to Figure 2, there is a 27 m vertical section near the end of the terminal in the Seonyu substation and the pressure rise at this vertical section should be included in the calculated pressure drop in Figure 8.
Figure 9 shows the pressure and temperature of each cooling system at the minimum mass flow rate along with the cable path. The pressure of the cable cryostat reaches its maximum in the vertical section. In particular, on the Seonyu side, the maximum pressure is approximately 2 bar greater than the outlet pressure, as the vertical section is closer to the inlet. The length of the path in the dual cooling system is shorter than that in the single cooling system; hence, the mass flow rate in the former is enough to satisfy the temperature limit.

7. Temperature and Pressure Analysis Results for Internal Return Path

7.1. Expanded One-Dimensional Thermo-Hydraulic Network Analysis Model

For the internal return path, conduction heat transfer ( Q ˙ c o n d u c t ) occurs in the core due to the temperature difference between the inflow and the return flow, as well as the electrical loss ( Q ˙ e l e c t r i c a l ) of the each conducting layer. Heat penetration ( Q ˙ t h e r m a l ) is introduced into the LN2 through the cryostat. Convective heat transfer occurs continuously between the flow and the core. To analyze the internal return path cooling system, we conducted an expanded 1-D network analysis using a commercial software, SINDA/FLUINT [22,23]. The model consists of one node per 50 m cable, and the total number of nodes is approximately 400. Diffusion node means a solid having mass and heat capacity. Junction node means a fluid with temperature and pressure. The connection line between the diffusion nodes is conduction heat transfer and the connection line between the junction node and the diffusion node is convection heat transfer. For thermal conduction inside the cable core, equivalent thermal conductivity and heat capacity for the components in the cable were modeled. The convection by LN2 flow inside and outside of the cable core was considered in the analysis model.
The analysis was performed by referring to the pressure rise in the vertical section and the pressure drop per unit length in the CFD analysis results, and the model was simulated. Figure 10 illustrates the expanded 1-D network model and the heat balance of the internal return path.
Table 8 is thermal properties used in the expanded 1-D network model.

7.2. Analysis Result for Internal Return Path

Figure 11 is the result of the analysis at the maximum temperature condition for single and dual cooling systems. For the single cooling system as in Figure 11a,b, the minimum mass flow rate was 0.8 kg/s to meet the temperature limit. However, due to the large pressure drop that occurs due the narrow flow area and the long path in the cable cryostat, it was not suitable for use. For the dual cooling system on the Munsan side as in Figure 11c, the minimum mass flow rate was 0.35 kg/s to meet the temperature limit. On the Seonyu side as in Figure 11d, the minimum mass flow rate was 0.4 kg/s; however, it could not be applied because the pressure limiting condition was not satisfied. To implement the dual cooling system on the Seonyu side, it was necessary to increase the diameter of the cryostat and move the vertical section away from the inlet. The temperature in each phase of the triaxial HTS power cable was electrically stable as the cable was operated between the temperatures of the inflow and internal return flow.

8. Summary of Analysis Results Using Phase Diagram

Temperature and pressure variations in Figure 9 and Figure 11 are the result of analysis with minimum flow rates that meet the temperature differences. In addition, the possible maximum flow rates can be obtained by considering the pressure limit. The above two cases can be clearly demonstrated in a temperature–pressure phase diagram as in Figure 12 and Figure 13, including pumping and cooling processes. Figure 12a,b are the phase diagrams for the single cooling system which is installed at Munsan and Seonyu with the external return path. The single cooling system at Munsan can operate from 0.57 to 0.66 kg/s. For the single cooling system at Seonyu, the pressure rise, due to the vertical section, occurred near the inlet; thus, the possible operating mass flow rate was reduced to 0.57 to 0.58 kg/s. With the dual cooling system, the allowable mass flow rate was 0.28 to 0.98 kg/s for section between the Munsan station to the platform and 0.3 to 0.8 kg/s for the section from the Seonyu station to the platform as in Figure 12c,d.
For an internal return path with a single cooling system, it was impossible to meet the pressure limit due to the reduced flow area. With the dual cooling system, it was possible to obtain the allowable mass flow rate, 0.35 to 0.41 kg/s from Munsan to the platform, while the pressure limit also exceeded the limit from Seonyu to the platform.

9. Summary and Conclusions

In this study, we conducted a thermo-hydraulic analysis of the single and dual cooling systems applicable to the external and internal return paths with the detailed geographical characteristics of the real installation environment. Through CFD analysis for the cable cryostat with spiral corrugations, it was possible to obtain the equivalent relative roughness to be used in 1-D hydraulic analysis for pressure drop. For a cable cooling system with an internal return path, we performed an expanded 1-D network analysis to consider the radial heat transfer within the cable core due to the temperature difference between the inflow and the return flow.
Table 9 summarizes the resultant allowable mass flow rates for each cooling mode. For the external return path of the single cooling system and the dual cooling system, we provided an adequate mass flow rate range. However, the internal return path cannot be operated in the single cooling system and dual cooling system for the Seonyu to platform section in spite of the low heat loss feature of the internal return path. To operate the cooling system with an internal return path, the flow area should be increased, and the vertical section should be moved away from the inlet, or a cable cryostat with a high allowable pressure should be used. We believe that the thermo-hydraulic analysis model and results will be used as a design and operation guide for a cryogenic cooling system for the future installation of the triaxial HTS power cable between the Munsan and Seonyu substations.
In this paper, we did not consider the fault situation. Prior to installing the triaxial HTS power cable, additional components in previously demonstrated HTS power cable systems, such as a decompressing cooling chamber and a parallel circulation pump, will be considered against the fault condition.

Author Contributions

Y.C. prepared the investigation and performed CFD and the expanded 1-D thermo-hydraulic analysis; C.L. developed CFD analysis model; D.K. developed the expanded 1-D thermo-hydraulic network model; S.K. provided a supervisory role and managed the research; D.W., J.Y. and H.Y., funding acquisition. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by Korea Electric Power Corporation.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chikumoto, N.; Watanabe, H.; Yury, V.I.; Takano, H.; Yamaguchi, S.; Koshizuka, H.; Hayashi, K.; Sawamura, T. Construction and the circulation test of the 500-m and 1000-m DC superconducting power cables in Ishikari. IEEE Trans. Appl. Supercond. 2016, 26, 1–4. [Google Scholar] [CrossRef]
  2. James, F.M.; Schmidt, F.; Bratt, S.; Tom, E.W.; Jie, Y. Installation and testing results of long island. IEEE Trans. Appl. Supercond. 2009, 19, 1692–1697. [Google Scholar]
  3. Yang, H.S.; Kim, D.L.; Sohn, S.H.; Lim, J.H.; Choi, H.O.; Choi, Y.S.; Lee, B.S.; Jung, W.M.; Ryoo, H.S.; Hwang, S.D. Long term performance test of KEPCO HTS power cable. IEEE Trans. Appl. Supercond. 2009, 19, 1782–1784. [Google Scholar] [CrossRef]
  4. Demko, J.A.; Sauers, I.; James, D.R.; Gouge, M.J.; Lindsay, D.; Roden, M.; Tolbert, J.; Willén, D.; Træholt, C.; Nielsen, C.T. Three-phase HTS cable for the AEP bixby project. IEEE Trans. Appl. Supercond. 2007, 17, 2047–2050. [Google Scholar] [CrossRef]
  5. Lee, S.; Kang, S.Y.; Park, M.; Won, D.; Yoo, J.; Yang, H.S. Performance analysis of real-scale 23 kV/60 MVA class tri-axial HTS power cable for real-grid application in Korea. Energies 2020, 13, 2053. [Google Scholar] [CrossRef]
  6. Jeong, K.W.; Moon, B.S.; Park, S.K. Status and Future Direction of HTS Power Application in KEPCO. Available online: http://www.csee.org.cn/data/zt_cigreAorc/doc/PL-18.pdf (accessed on 30 June 2020).
  7. Lee, C.; Kim, D.; Kim, S.; Won, D.Y.; Yang, H.S. Thermo-hydraulic analysis on long three-phase HTS power cable of several kilometers. IEEE Trans. Appl. Supercond. 2019, 29, 5402805. [Google Scholar] [CrossRef]
  8. Takeshi, F.; Jun, N.; Mitsugi, O.; Takao, Y.; Osami, S.T. Analysis of AC loss and current distribution characteristics of multi-layer three-phase HTS cable for 3-phase AC power transmission. IEEE Trans. Appl. Supercond. 2006, 16, 135–138. [Google Scholar]
  9. Sasaki, A.; Ivanov, Y.; Yamaguchi, S. LN2 circulation in cryopipes of superconducting power transmission line. Cryogenics 2011, 51, 471–476. [Google Scholar] [CrossRef]
  10. Yeom, H.-K.; Koh, D.-Y.; Kim, S.-H.; Seo, T.-B. Design of cryogenic systems for 154 kV HTS power cable. IEEE Trans. Appl. Supercond. 2010, 20, 1296–1299. [Google Scholar] [CrossRef]
  11. He, J.; Tang, Y.; Wei, B.; Li, J.; Ren, L.; Shi, J.; Wu, K.; Li, X.; Xu, Y.; Wang, S. Thermal analysis of HTS power cable using 3-D FEM model. IEEE Trans. Appl. Supercond. 2013, 23, 5402404. [Google Scholar]
  12. Thadela, S.; Rao, V.V.; Agarwal, R.; Dondapatib, R.S. Computational investigation on thermohydraulic characteristics of high temperature superconducting (HTS) power cables. Phys. Supercond. Appl. 2019, 559, 25–31. [Google Scholar] [CrossRef]
  13. Zuo, Z.Q.; Jiang, W.B.; Yu, Z.G.; Huang, Y.H. Liquid nitrogen flow in helically corrugated pipes with insertion of high-temperature superconducting power transmission cables. Int. J. Heat Mass Transf. 2019, 140, 88–99. [Google Scholar] [CrossRef]
  14. Dong, Y.; Huixiong, L.; Tingkuan, C. Pressure drop, heat transfer and performance of single-phase turbulent flow in spirally corrugated tubes. Exp. Therm. Fluid Sci. 2001, 24, 131–138. [Google Scholar] [CrossRef]
  15. Fuchino, S.; Furuse, M.; Higuchi, N. Longitudinal temperature distribution in superconducting power cables with counter-flow cooling. IEEE Trans. Appl. Supercond. 2002, 12, 1339–1342. [Google Scholar] [CrossRef]
  16. Zajaczkowski, B.; Giesbers, A.J.M.; Holtrust, M.; Haenen, E.; den Heijer, R. Feasibility of inline cooling in long distance HTS power line. Cryogenics 2011, 51, 180–186. [Google Scholar] [CrossRef]
  17. Miyagi, D.; Sakakibara, R.; Shinozaki, Y.; Tsuda, M.; Hamajima, T. Suitable cable structure of HTS three-phase cable cooled by counter flow cooling method for long-distance power transmission. IEEE Trans. Appl. Supercond. 2018, 28, 5400705. [Google Scholar] [CrossRef]
  18. Lee, S.; Sung, H.; Park, M.; Won, D.; Yoo, J.; Yang, H.S. Analysis of the temperature characteristics of three-phase superconducting power cable according to a liquid nitrogen circulation method for real-grid application in Korea. Energies 2019, 12, 1740. [Google Scholar] [CrossRef] [Green Version]
  19. Stultz, S.C.; Kitto, B. Steam its Generation and Use, 41st ed.; The Babcock & Wilcox Company: Barberton, OH, USA, 2005; pp. 3–14. [Google Scholar]
  20. Marušic-Paloka, E.; Pažanin, I. Effects of boundary roughness and inertia on the fluid flow through a corrugated pipe and the formula for the Darcy—Weisbach friction coefficient. Int. J. Eng. Sci. 2020, 152, 103293. [Google Scholar] [CrossRef]
  21. Chang, H.; Ryu, K.N.; Yang, H.S. Cryogenic design of liquid-nitrogen circulation system for long-length HTS cables with altitude variation. Cryogenics 2017, 83, 50–56. [Google Scholar] [CrossRef]
  22. Heat Transfer and Fluid Analysis Software, SINDA/FLUINT. Available online: crtech.com/products/sindafluint (accessed on 27 June 2020).
  23. Sakowski, B.; Daniel, M.H.; Jason, W.H.; Kassemi, M. Validation of heat transfer correlations in line chill-down tests of cryogenic fluid in SINDA/FLUINT. In Proceedings of the 2018 International Energy Conversion Engineering Conference, Cincinnati, OH, USA, 9–11 July 2018. [Google Scholar]
  24. Choi, Y.S.; Kim, D.L. Thermal property measurement of insulating material used in HTS power device. Cryogenics 2012, 52, 465–470. [Google Scholar] [CrossRef]
Figure 1. Structure of the triaxial HTS power cable with: (a) an external return path; (b) an internal return path.
Figure 1. Structure of the triaxial HTS power cable with: (a) an external return path; (b) an internal return path.
Energies 13 03898 g001
Figure 2. Layout of triaxial HTS power cable between Munsan and Seonyu substations.
Figure 2. Layout of triaxial HTS power cable between Munsan and Seonyu substations.
Energies 13 03898 g002
Figure 3. Installation environment in the Jeju smart-grid demonstration center: (a) Cable snake and cable bending radius; (b) cable snake assumed as bending.
Figure 3. Installation environment in the Jeju smart-grid demonstration center: (a) Cable snake and cable bending radius; (b) cable snake assumed as bending.
Energies 13 03898 g003
Figure 4. Minor loss coefficient for a bended round pipe (Babcock & Wilcox Co., Barberton, OH, USA, 1978) [19].
Figure 4. Minor loss coefficient for a bended round pipe (Babcock & Wilcox Co., Barberton, OH, USA, 1978) [19].
Energies 13 03898 g004
Figure 5. Configurations of cooling system with an external return path: (a) single cooling system at Munsan; (b) dual cooling system.
Figure 5. Configurations of cooling system with an external return path: (a) single cooling system at Munsan; (b) dual cooling system.
Energies 13 03898 g005
Figure 6. Streamlines in the computational fluid dynamic (CFD) analysis result for the: (a) cable with an external return path; (b) cable with internal return path at 0.57 kg/s.
Figure 6. Streamlines in the computational fluid dynamic (CFD) analysis result for the: (a) cable with an external return path; (b) cable with internal return path at 0.57 kg/s.
Energies 13 03898 g006
Figure 7. Pressure drop per unit length according to mass flow rate.
Figure 7. Pressure drop per unit length according to mass flow rate.
Energies 13 03898 g007
Figure 8. Allowable mas flow rates that meets the pressure limit and temperature limit for various cooling modes: (a) single cooling system; (b) dual cooling system, Munsan side; (c) dual cooling system, Seonyu side.
Figure 8. Allowable mas flow rates that meets the pressure limit and temperature limit for various cooling modes: (a) single cooling system; (b) dual cooling system, Munsan side; (c) dual cooling system, Seonyu side.
Energies 13 03898 g008
Figure 9. Temperature and pressure variations along the HTS power cable with an external return path: (a) single cooling system at Munsan; (b) single cooling system at Seonyu; (c) dual cooling system, Munsan side; (d) dual cooling system, Seonyu side. The mass flow rates are minimum values that meet the allowable temperature difference.
Figure 9. Temperature and pressure variations along the HTS power cable with an external return path: (a) single cooling system at Munsan; (b) single cooling system at Seonyu; (c) dual cooling system, Munsan side; (d) dual cooling system, Seonyu side. The mass flow rates are minimum values that meet the allowable temperature difference.
Energies 13 03898 g009aEnergies 13 03898 g009b
Figure 10. Expanded 1-D network model and heat balance diagram for the internal return path.
Figure 10. Expanded 1-D network model and heat balance diagram for the internal return path.
Energies 13 03898 g010
Figure 11. Temperature and pressure result at the position of internal return path: (a) single cooling system at Munsan; (b) single cooling system at Seonyu; (c) dual cooling system on the Munsan side; (d) dual cooling system on the Seonyu side.
Figure 11. Temperature and pressure result at the position of internal return path: (a) single cooling system at Munsan; (b) single cooling system at Seonyu; (c) dual cooling system on the Munsan side; (d) dual cooling system on the Seonyu side.
Energies 13 03898 g011
Figure 12. Possible operating range between pressure and temperature limit with external return path in phase diagram: (a) single cooling system at Munsan; (b) single cooling system at Seonyu; (c) dual cooling system, Munsan side; (d) dual cooling system, Seonyu side.
Figure 12. Possible operating range between pressure and temperature limit with external return path in phase diagram: (a) single cooling system at Munsan; (b) single cooling system at Seonyu; (c) dual cooling system, Munsan side; (d) dual cooling system, Seonyu side.
Energies 13 03898 g012
Figure 13. Possible operating range between pressure and temperature limit with internal return path in phase diagram, Munsan to platform.
Figure 13. Possible operating range between pressure and temperature limit with internal return path in phase diagram, Munsan to platform.
Energies 13 03898 g013
Table 1. Specification of the triaxial high-temperature superconducting (HTS) power cable.
Table 1. Specification of the triaxial high-temperature superconducting (HTS) power cable.
TypeElectric Power CapacityLine-To-Line VoltageCurrentHTS WireInsulation
Tri-axial HTS power cable60 MVA23 kV1.5 kArmsBrass laminated coated conductorPPLP
Table 2. Total minor loss coefficient of the HTS power cable between Munsan and Seonyu. substation.
Table 2. Total minor loss coefficient of the HTS power cable between Munsan and Seonyu. substation.
Angle of Bend (θ)Number of BendsMinor Loss ( K b )
90°80.131 × 8
66°10.114 × 1
26°10.067 × 1
14°10.043 × 1
20.16°11800.056 × 1180
Total minor loss coefficient67.4
Pressure drop (Pa) @70 K, 7 bar and 0.5 kg/s1439
Table 3. Allowable operating range of cooling system for triaxial HTS power cable.
Table 3. Allowable operating range of cooling system for triaxial HTS power cable.
ParameterUnitValueComment
Max. pressure limitbar13.5Allowable pressure of cable cryostat
Min. Pressure limitbar5.5Suppress the bubble
Max. temperatureK75.5Maximum temperature of HTS power cable
Min. temperatureK67.0Preventing freezing
Table 4. Heat loss of the triaxial HTS power cable with external return tube.
Table 4. Heat loss of the triaxial HTS power cable with external return tube.
Munsan to PlatformSeonyu to Platform
ParameterLength, setHeat loss per unitTotal Heat lossLength, setHeat loss per unitTotal Heat loss
HTS power cable ( Q ˙ c a b l e )983 m3.5 W/m3441 W1081 m3.5 W/m3784 W
External return tube
( Q ˙ r e t u r n )
983 m1.0 W/m983 W1081 m1.0 W/m1081 W
Terminal
( Q ˙ t e r m i n a t i o n )
2 set550 W/EA1100 W2 set550 W/EA1100 W
Joint box
( Q ˙ j o i n t b o x )
2 set80 W/EA160 W2 set80 W/EA160 W
Sum5.68 kWSum6.13 kW
Table 5. Representative value of flow parameters for each case.
Table 5. Representative value of flow parameters for each case.
Return TypeFlow PathDh (mm)Red
Mass Flow Rate 0.1 kg/sMass Flow Rate 1 kg/s
ExternalInflow (annular)20.2439643,961
Return flow93.1783678,369
InternalInflow33.619,943199,435
Return flow (annular)15.7428042,801
Table 6. CFD analysis boundary condition.
Table 6. CFD analysis boundary condition.
ParameterCondition
Working fluidSubcooled liquid nitrogen
Mass flow rate0.28–0.95 kg/s
LN2 temperature70 K
Average outlet pressure5.5 bar, g (static gauge pressure)
Model length2 m
Wall roughnessCore–1 mm rough wall,
Cable cryostat–smooth wall
Wall conditionNo-slip condition
Turbulence model κ ε
Table 7. Equivalent relative roughness of the HTS power cable CFD result.
Table 7. Equivalent relative roughness of the HTS power cable CFD result.
Return TypeFlow PathRelative Roughness of CFD
(0.28–0.95 kg/s)
Equivalent Relative Roughness ( ε / D h ) for Equation (5)
ExternalInflow0.238–0.2690.27
Return flow0.114–0.2010.21
InternalInflow0.252–0.2660.27
Return flow0.303–0.3200.32
Table 8. Thermal properties of expanded 1-D network model.
Table 8. Thermal properties of expanded 1-D network model.
Diffusion Node NameMaterialThermal Conductivity
(W/m·K) @ 70 K
Comments
Conducting layer-60.22Averaged value of HTS wire
StabilizerCopper512.96
InsulationPPLP0.05Based on the measurement [24]
Table 9. Allowable mass flow rates for various cooling modes.
Table 9. Allowable mass flow rates for various cooling modes.
Return PathCooling ModeSectionPossible Mass Flow Rate (kg/s)
ExternalSingle cooling at MunsanMunsan to Seonyu0.57–0.66
Single cooling at SeonyuSeonyu to Munsan0.57–0.58
Dual coolingMunsan to platform0.28–0.98
Seonyu to platform0.30–0.80
InternalSingle cooling at MunsanMusan to Seonyuexceed pressure limit
Single cooling at SeonyuSeonyu to Munsanexceed pressure limit
Dual coolingMunsan to platform0.35–0.41
Seonyu to platformexceed pressure limit

Share and Cite

MDPI and ACS Style

Choi, Y.; Kim, D.; Lee, C.; Won, D.; Yoo, J.; Yang, H.; Kim, S. Thermo-Hydraulic Analysis of a Tri-Axial High-Temperature Superconducting Power Cable with Respect to Installation Site Geography. Energies 2020, 13, 3898. https://doi.org/10.3390/en13153898

AMA Style

Choi Y, Kim D, Lee C, Won D, Yoo J, Yang H, Kim S. Thermo-Hydraulic Analysis of a Tri-Axial High-Temperature Superconducting Power Cable with Respect to Installation Site Geography. Energies. 2020; 13(15):3898. https://doi.org/10.3390/en13153898

Chicago/Turabian Style

Choi, Youngjun, Dongmin Kim, Changhyung Lee, Duyeon Won, Jaeun Yoo, Hyungsuk Yang, and Seokho Kim. 2020. "Thermo-Hydraulic Analysis of a Tri-Axial High-Temperature Superconducting Power Cable with Respect to Installation Site Geography" Energies 13, no. 15: 3898. https://doi.org/10.3390/en13153898

APA Style

Choi, Y., Kim, D., Lee, C., Won, D., Yoo, J., Yang, H., & Kim, S. (2020). Thermo-Hydraulic Analysis of a Tri-Axial High-Temperature Superconducting Power Cable with Respect to Installation Site Geography. Energies, 13(15), 3898. https://doi.org/10.3390/en13153898

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