Experiment and analysis on flow rate and temperature of liquid PbBi in PREKY

Experiment and analysis on flow rate and temperature of liquid PbBi in PREKY 1 2 first-author 2 2 corresp Corresponding author, zhiqiang.zhu@fds.org.cn Corresponding author, zhiqiang.zhu@fds.org.cn zhiqiang.zhu@fds.org.cn 2 2 University of Science and Technology of China, Hefei, Anhui, 230027, China University of Science and Technology of China, Hefei, Anhui, 230027, China Key Laboratory of Neutronics and Radiation Safety, Institute of Nuclear Energy Safety Technology, Chinese Academy of Sciences, Hefei, Anhui, 230031, China Key Laboratory of Neutronics and Radiation Safety, Institute of Nuclear Energy Safety Technology, Chinese Academy of Sciences, Hefei, Anhui, 230031, China Flow rate and temperature are important parameters for design and operation of liquid lead-bismuth (PbBi) experimental loop. The PREKY facility was designed to study test technique of flow rate and local temperature of PbBi loop. In this work, flow rate monitoring of molten PbBi was performed, and temperatures of the PbBi and pipe surface were measured. The results show that the flow rate of the Venturi-nozzle flow meter (VNFM) had an uncertainty of ±5% in the range of 0.6~2.0 m/s, and the maximum temperature difference between the PbBi and pipe surface was about 8°C. PREKY can be used to measure high temperature PbBi. Lead-bismuth PREKY loop Flow rate Temperature measurement 1 Introduction 1 For advanced reactors, the use of molten metal or alloy is an attractive choice [ 1 1 , 2 2 , 3 3 ] . Heavy liquid metal (HLM) such as lead or lead-bismuth (PbBi) is proposed as a promising candidate coolants for next generation reactors due to their favorable properties of good neutron properties, heat transfer performance, chemical inert and security features etc . [ 4 4 , 5 5 ] It can also be used as a spallation target for accelerator-driven systems (ADS), which has been proposed for the transmutation of nuclear waste [ 6 6 ] . The KYLIN-II, a multi-function liquid lead-bismuth eutectic (LBE) experimental platform with liquid PbBi loops, has been built at Institute of Nuclear Energy Safety Technology (INEST) [ 7 7 ] . The planned experimental activities will be focused on compatibility of structure materials, thermal hydraulic characteristics of HLM, fuel assembly testing, and safety study etc . [ 8 8 , 9 9 ] . In addition, measurement technology for liquid PbBi can be studied on the platform. A high temperature liquid PbBi experimental facility, PREKY, was established to study key techniques and construction technology of liquid metal loop, such as flow meters with high precision for input reference of thermal hydraulic analysis and thermocouples with long term lifetime for local temperature etc . Besides the heavy density, corrosiveness and opaqueness, direct measurement of flowing PbBi would face temperatures of 200°C– 450°C or even higher [ 10 10 - 12 12 10-12 ] . So, the measurement components should have the following features. (1) Operator safety shall be of top priority in the experiment; (2) Experimental data shall be repeatable and their accuracy and repeatability can be verified; and (3) Each sensor can be independently calibrated and easily replaced. In this paper, main characteristics of the PREKY loop are presented and the experimental results are discussed. The results can evaluate feasibility of the instrumentations. 2 The PREKY facility 1 PREKY was designed to study and evaluate feasibility of measurement instruments, which include thermal- couples, flow-meters, level meters and differential pressure transducers. They were commonly used in standard industrial applications. The scheme of PREKY is shown in Fig.1 Fig.1 . It consists of the storage vessel and calibration tank, electromagnetic pump, heaters and heat exchanger, electrical system, data acquisition and control unit. Table 1 Table 1 lists the main design parameters, which cover completely the thermal measurement requirements in KYLIN-II loop. 1 1 Parameters 1 1 Values 1 1 Working temperature (°C) 1 1 200~500 1 1 Working pressure (MPa) 1 1 0.5~1.0 1 1 Velocity (m·s −1 ) 1 1 0.5~2.0 1 1 Maximum flow resistance (MPa) 1 1 0.25 The PREKY loop under operation conditions of <1 MPa and <450°C includes the following four tests: (1) flow meter performance test relative to level meters; (2) local temperature comparative experiments for both the pipe and flowing molten PbBi; (3) electromagnetic pump (EMP) performance for the coupling curve of power/efficiency and flow rate; and (4) pressure drop test of valves at different flow rates. Tests (1) and (2) will be described in detail in this paper. 3 Results and Discussion 1 3.1. Flow rate test 2 Owing to the opaque nature of liquid PbBi, especially under high temperature and chemical reaction conditions, traditional measuring technology for water or other flowing medium cannot be applied to the PREKY loop. Therefore, the flow meter should be calibrated in advance. Based on previous experience and test results, the Venturi-nozzle flow meter (VNFM) was chosen for flow rate test of the PREKY loop. The scheme of VNFM is shown in Fig.2 Fig.2 . In principle, VNFM can be used in a wide range of flow rates based on the measuring principle of differential pressures. The flow rate ( Q ) of VNFM is calculated by Eq.(1) [ 13 13 ] : <math display="block" id="M1"><mrow><mi>Q</mi><mo>=</mo><mi>μ</mi><mo>⋅</mo><mi>K</mi><mo>⋅</mo><msqrt><mrow><mi>Δ</mi><mi>h</mi></mrow></msqrt><mtext> ,</mtext></mrow></math> where, μ is flow coefficient of the VNFM, Δ h is the pressure loss between pressure ports located at the entrance 1-1 and throat 2-2 (see Fig.2 Fig.2 ); and K is a constant, determined by Eq.(2) [ 13 13 ] : <math display="block" id="M2"><mrow><mi>K</mi><mo>=</mo><mfrac><mrow><mi>π</mi><msubsup><mi>d</mi><mn>2</mn><mn>2</mn></msubsup></mrow><mn>4</mn></mfrac><mo>⋅</mo><mfrac><mrow><msqrt><mrow><mn>2</mn><mi>g</mi></mrow></msqrt></mrow><mrow><msqrt><mrow><mn>1</mn><mo>−</mo><msup><mrow><mrow><mo>(</mo><mrow><mfrac><mrow><msub><mi>d</mi><mn>2</mn></msub></mrow><mrow><msub><mi>d</mi><mn>1</mn></msub></mrow></mfrac></mrow><mo>)</mo></mrow></mrow><mn>4</mn></msup></mrow></msqrt></mrow></mfrac><mtext> </mtext><mo>,</mo></mrow></math> where, d 1 and d 2 are inner diameters of the pipe at the entrance and throat (see Fig.2 Fig.2 ), respectively; and g is the acceleration of gravity. The flow coefficient μ is determined by Eq.(3): <math display="block" id="M3"><mrow><mi>μ</mi><mo>=</mo><msqrt><mrow><mfrac><mrow><mn>1</mn><mo>−</mo><msup><mrow><mrow><mo>(</mo><mrow><msub><mi>d</mi><mn>2</mn></msub><mo>/</mo><msub><mi>d</mi><mn>1</mn></msub></mrow><mo>)</mo></mrow></mrow><mn>4</mn></msup></mrow><mrow><msub><mi>a</mi><mn>2</mn></msub><mo>+</mo><mi>ζ</mi><mo>−</mo><msub><mi>a</mi><mn>1</mn></msub><mo>⋅</mo><msup><mrow><mrow><mo>(</mo><mrow><msub><mi>d</mi><mn>2</mn></msub><mo>/</mo><msub><mi>d</mi><mn>1</mn></msub></mrow><mo>)</mo></mrow></mrow><mn>4</mn></msup><mo>+</mo><mi>λ</mi><mo>⋅</mo><mi>l</mi><mo>/</mo><msub><mi>d</mi><mn>2</mn></msub></mrow></mfrac></mrow></msqrt><mtext> ,</mtext></mrow></math> where (see Fig.2 Fig.2 ), l is the length between the 1-1 and 2-2; λ is the drag coefficient; a 1 and a 2 are the correction coefficients of kinetic energy at the 1-1 and 2-2, respectively; ζ is the local head loss coefficient from 1-1 to 2-2. In this experiment, the μ is evaluated at 0.9316. Calibration experiment of the VNFM was implemented to obtain the relationship between flow rate Q and pressure loss △ h . The measurement data were compared with the calculation results using the above equations. At beginning of the experiment, the output signal under differential pressures was unstable and changed due to the less wetting of internal surface of the flow channel [ 14 14 , 15 15 ] . This, however, can be sufficiently improved by flowing liquid PbBi at large flow rate and high temperature for tens of hours before the calibration experiment [ 16 16 ] . The flow rate calibration was carried out at different flow rates. Gravity was the driving force of liquid PbBi flow, and different flow rates were obtained by adjusting the valve opening at the calibration tank exit. The calibration procedures are as follows: (1) The system is fully preheated, and evacuated for at least 3 times; (2) Liquid PbBi is driven to the safety level in the calibration tank; and (3) Open the valve quickly to a certain degree, and record the experiment data. In the calibration, the calibration vessel and storage vessel were kept at the same gas pressures. To improve the value precision and extend the range of flow rate, the experiment was repeated for seven times. The maximum flow rate was up to 2 m/s, which was about the maximum flow rate of the KYLIN-II loop, for corrosion and thermal hydraulic experiments. Typically, the experiment was performed at about 350°C and flow rate of 2.0 m/s. Fig.3 Fig.3 shows that basically, the VNFM data agree with the level meter measurement results. The flow rate of VNFM increases stably each time, while the level meter curve fluctuates a little comparing with the VNFM output, and slope of the level meter curve is a little bigger than that of VNFM. The relative errors between VNFM and level meter measurement are shown in Fig.3(b) Fig.3(b) . At flow rates of below about 0.6 m/s, the maximum measurement error is about 20% or even more, so the VNFM is not suitable for measurements at low flow rates. It can also be seen that the flow rate calibration has an uncertainty of ±5% in the range of 0.6‒2.0 m/s. This is because that: (1) the driving force for liquid PbBi to flow through the calibration vessel is the gravity, hence the non-steady flowing state between two level electrodes; and (2) changes in diameter and height of the calibration vessel under different temperatures were not considered in the experiment. 3.2. Local temperature test 2 Temperature controlling and monitoring are of importance for safe, reliable and stable operation of molten PbBi facilities. Accurate measurement of temperature distribution is an important emergency control for safe operation. For studying heat transfer characteristics of PbBi in the channel, temperature tests were implemented by two quick response thermocouples of K-type ( Table 2 Table 2 and Fig. 4 Fig. 4 ) installed inside and outside the flow channel located at the outlet of air-cooled heat exchanger. The inside thermal-couple was directly contacted with flowing PbBi itself under high temperature. 1 1 Parameters 1 1 Values 1 1 Dimension (mm) 1 1 Φ 3×500 1 1 Maximum measuring range (°C) 1 1 ~1000 1 1 Accuracy (°C) 1 1 ±1.5(below 400°C) 1 1 Maximum working pressure (MPa) 1 1 ~1 1 1 Response time (s) 1 1 10 1 1 Material of protection tube 1 1 SS 316L The heat exchanger was used to cool down the PbBi flowing out of the test section. Its inlet and outlet were opened manually. A few minutes were needed to reach stable state of flowing PbBi. The thermal-couples were contacted with the pipe surface, and the heating and controlling system would need a thermal response time of over 60 s. In other words, the heat exchanger would start to stop working 60 seconds later when the temperature of outer surface of flow channel reached the preset value, while the temperature of liquid PbBi had exceeded the preset value at this moment. Therefore there were temperature difference at the same position and time between flowing PbBi and the outer surface of flow channel. Thermal response time of the thermal-couple, which depends on outer diameter of the measuring terminal, can be defined by Eq.(4) [ 13 13 ] : <math display="block" id="M4"><mrow><mi>τ</mi><mo>=</mo><mi>ρ</mi><mi>V</mi><mi>c</mi><mo>/</mo><mi>h</mi><mi>S</mi><mtext> ,</mtext></mrow></math> where, h is heat transfer coefficient; and ρ , V , c and S are density, volume, specific heat and surface area of the thermal-couple, respectively. From Eq. (1), supposed L is the terminal diameter of the thermal-couple, the volume V ∝ L 3 , the surface area S ∝ L 2 , similarly the thermal response time τ ∝( V / S ) ∝ L . Therefore, the thermal response time τ is proportional to L . So, the terminal diameter is 3mm, thus response time is about 10 s. The results of temperature experiments were shown in Fig.5 Fig.5 . It can be seen that (1) the working period of the cooler was about 10 min, the temperature fluctuations of flowing PbBi and the pipe surface were about 35°C and 18°C, respectively; (2) the maximum temperature difference between the PbBi and pipe surface was roughly 8°C； and (3) the temperature peak positions of pipe surface was about 1 min later than that of the PbBi, and the maximum and minimum temperatures of PbBi, and of the pipe surface, had a time difference of about 4.5 min. These can be explained as follows: (1) The two thermocouples were located at the outlet of the heat exchanger. When the heat exchanger stopped working, output of the two thermocouples increased for a short period of time due to continuous heat release from the electromagnetic pump. After a few seconds, temperature went down quickly by heat exchange in the working state. At this time, a working period for the heat exchanger had been completed. (2) When temperature of liquid PbBi began to drop, due to heat transfer of liquid PbBi and thermal response time of thermocouples, temperature of pipe surface started to drop about 60 s later. The maximum temperature difference (about 8°C) between PbBi and pipe surface did not appear in the inflexion of the operating state for the heat exchanger. Following the above results, thermocouples could be in direct contact with flowing PbBi to obtain the accurate and real-time temperature of liquid PbBi. 4 Summary 1 The VNFM experiment results indicate that the measuring technique can meet the experiments requirements. VNFM can be used to measure the flow rate of 0.6‒2.0 m/s. The uncertainty of ±5% in the flow rate measurement was probably due to the driving force of liquid PbBi and temperature difference. The results of temperature monitoring experiment show that the maximum temperature difference between flowing PbBi and pipe surface was about 8°C. The thermal response time of the thermal-couple was about 10 s. These can be considered as an inevitable difference. These results were intended to provide guidance for the planned experiments in the KYLIN-II facility. References Wu Y C , Jiang J Q , Wang M H , et al .. A Fusion-Driven Subcritical System Concept Based on Viable Technologies . Nucl Fus , 2011 , 51 : 103036 . DOI: 10.1088/0029-5515/51/10/103036 http://doi.org/10.1088/0029-5515/51/10/103036 A Fusion-Driven Subcritical System Concept Based on Viable Technologies LU Y , HE J , ZHU Z Q , et al . Preliminary Calibration Test and Analysis of Electromagnetic Flow-Meter in Liquid Lead-Bismuth . Nucl Tech , 2014 , 37 : 080603 . DOI: 10.11889/j.0253-3219.2014.hjs.37.080603 http://doi.org/10.11889/j.0253-3219.2014.hjs.37.080603 Preliminary Calibration Test and Analysis of Electromagnetic Flow-Meter in Liquid Lead-Bismuth HE J , ZHU Z Q , HUANG Q Y , et al .. Test and Analysis on the Flow Resistance of Y-Type Globe Valve in the Liquid Lead-Bismuth Alloy . Nucl Tech , 2015 , 38 : 040604 . DOI: 10.11889/j.0253-3219.2015.hjs.38.040604 http://doi.org/10.11889/j.0253-3219.2015.hjs.38.040604 Test and Analysis on the Flow Resistance of Y-Type Globe Valve in the Liquid Lead-Bismuth Alloy Wu Y C , FDS Team . Conceptual Design Activities of FDS Series Fusion Power Plants in China . Fusion Eng Des , 2006 , 81 : 2713 - 2718 . DOI: 10.1016/j.fusengdes.2006.07.068 http://doi.org/10.1016/j.fusengdes.2006.07.068 Conceptual Design Activities of FDS Series Fusion Power Plants in China Wang Y Q , Huang Q Y , Wu B , et al .. The Performance of Pt/air Oxygen Sensors in Stagnant Pb-Bi Eutectic at High Temperatures . Nucl Sci Tech , 2014 , 25 : 060602 . DOI: 10.13538/j.1001-8042/nst.25.060602 http://doi.org/10.13538/j.1001-8042/nst.25.060602 The Performance of Pt/air Oxygen Sensors in Stagnant Pb-Bi Eutectic at High Temperatures Wang J , Pei X , Cao R F , et al .. A Multiphase Direct Aperture Optimization for Inverse Planning in Radiotherapy . Nucl Sci Tech , 2015 , 26 : 010502 . DOI: 10.13538/j.1001-8042/nst.26.010502 http://doi.org/10.13538/j.1001-8042/nst.26.010502 A Multiphase Direct Aperture Optimization for Inverse Planning in Radiotherapy Chen C , Zou J , Zheng J et al .. Development and Testing of the Code for Automatic Generating of Multi-Temperature Continuous-Energy Neutron Cross Section Libraries . Nucl Sci Tech , 2014 , 25 : 050602 . DOI: 10.13538/j.1001-8042/nst.25.050602 http://doi.org/10.13538/j.1001-8042/nst.25.050602 Development and Testing of the Code for Automatic Generating of Multi-Temperature Continuous-Energy Neutron Cross Section Libraries Wu Y C , FDS Team . Conceptual Design of the China Fusion Power FDS-II . Fusion Eng Des , 2008 , 83 : 1683 - 1689 . DOI: 10.1016/j.fusengdes.2008.06.048 http://doi.org/10.1016/j.fusengdes.2008.06.048 Conceptual Design of the China Fusion Power FDS-II Yao C M , Huang Q Y , Zhu Z Q , et al .. Data Acquisition and Control System for Lead-Bismuth Loop KYLIN-II-M . Nucl Sci Tech , 2015 , 26 : 010401 . DOI: 10.13538/j.1001-8042/nst.26.010401 http://doi.org/10.13538/j.1001-8042/nst.26.010401 Data Acquisition and Control System for Lead-Bismuth Loop KYLIN-II-M Claudia F , Scaddozzo G , Tarantino M , et al .. ENEA experience in LBE technology . Nucl Mater , 2006 , 356 : 264 - 272 . DOI: 10.1016/j.jnucmat.2006.05.020 http://doi.org/10.1016/j.jnucmat.2006.05.020 ENEA experience in LBE technology Gianluca B , Foletti C , Nicola F , et al .. Experimental study on gas-injection enhanced circulation performed with the CIRCE facility . Nuclear Eng Des , 2007 , 237 : 768 - 777 . DOI: 10.1016/j.nucengdes.2006.09.005 http://doi.org/10.1016/j.nucengdes.2006.09.005 Experimental study on gas-injection enhanced circulation performed with the CIRCE facility Zhao K , Cheng M Y , Long P C et al .. A Hybrid Voxel Sampling Method for Constructing Rad-HUMAN Phantom . Nucl Sci Tech , 2014 , 25 : 020503 . DOI: 10.13538/j.1001-8042/nst.25.020503 http://doi.org/10.13538/j.1001-8042/nst.25.020503 A Hybrid Voxel Sampling Method for Constructing Rad-HUMAN Phantom John E , Joseph B . Fliud Mechanics with Engineering Appications . BeiJing (CHINA) : Tsinghua University Press , 2003 , 122 - 156 . Fliud Mechanics with Engineering Appications Jae H , Batta A , Casamassima V , et al .. Benchmarking of thermal hydraulic loop models for Lead-Alloy Cooled Advanced Nuclear Energy System (LACANES), phase-I: Isothermal steady state forced convection . Nucl Mater , 2011 , 415 : 404 - 414 . DOI: 10.1016/j.jnucmat.2011.04.043 http://doi.org/10.1016/j.jnucmat.2011.04.043 Benchmarking of thermal hydraulic loop models for Lead-Alloy Cooled Advanced Nuclear Energy System (LACANES), phase-I: Isothermal steady state forced convection Thomas S , Robert S . Flow measurement techniques in heavy liquid metals . Nuclear Eng Des , 2010 , 240 : 2077 - 2087 . DOI: 10.1016/j.nucengdes.2009.11.017 http://doi.org/10.1016/j.nucengdes.2009.11.017 Flow measurement techniques in heavy liquid metals Ambrosini W , Azzati M , Benamati G , et al .. Testing and qualification of CIRCE instrumentation based on bubble tubes . Nucl Mater , 2004 , 335 : 293 - 298 . DOI: 10.1016/j.jnucmat.2004.07.030 http://doi.org/10.1016/j.jnucmat.2004.07.030 Testing and qualification of CIRCE instrumentation based on bubble tubes This work was supported by ITER 973 project (No. 2014GB112002 and 2014GB116000), Strategic Priority Research Program of Chinese Academy of Sciences (No. XDA03040200) and Foundation of President of Hefei Institute of Physical Science (No. YZJJ201328). The authors would like to thank other members of the FDS Team for their helps in this research. 10.1007/s41365-016-0036-3 2014-11-25 2015-05-04 2015-06-30 2016-04-14 2016-04-20 2021-04-11 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2016 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-016-0036-3 http://dx.doi.org/10.1007/s41365-016-0036-3 . 2016 Nuclear Science and Techniques 02 27