Corresponding author: Akram H. Abed ( akraaam82@yahoo.com ) Academic editor: Boris Balakin
© 2019 Akram H. Abed, Sergey E. Shcheklein, Valery M. Pakhaluev.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Abed AH, Shcheklein SE, Pakhaluev VM (2019) Heat transfer intensification in emergency cooling heat exchanger and dry cooling towers on nuclear power plant using airwater mist flow. Nuclear Energy and Technology 5(4): 281287. https://doi.org/10.3897/nucet.5.47972

Advanced nuclear power plants are equipped with passive emergency heat removal systems (PEHRS) for removing the decay heat from reactor equipment in accidents accompanied by primary circuit leakage to the final heat absorber (ambient air). Herein, the intensity of heat dissipation to air from the outer surface of the heat exchanger achieved by buoyancy induced natural convection is extremely low, which need to a large heat exchanger surface area, apply different types of heat transfer intensification including (grooves, ribs and extended surfaces, positioning at higher altitudes, etc.). The intensity of heat removal is also strongly dependent on the ambient air temperature (disposable temperature head).
Construction of nuclear power plants in countries with high ambient temperatures (Iran, Bangladesh, Egypt, Saudi Arabia, and others) which are characterized by a high level of ambient temperature imposes additional requirements on the increase of the heat exchange surfaces.
The experimental investigation results of heat transfer intensification by a low energy ultrasonic which supply a fine liquid droplet (size ~3 µm) in the cooling air are presented in the present paper. In such case, the heat transfer between the surface and cooling flow involves the following three physical effects: convection, conductive heat transfer, and evaporation of water droplets. The last two effects weakly depend on the ambient air temperature and provide an active heat removal in any situation.
The investigation was performed using a highprecision calorimeter with a controlled rate of heat supply (between 7800 and 12831 W/m^{2}) imitating heated surface within the range of Reynolds numbers from 2500 to 55000 and liquid (water) flow rates from 23.39 to 111.68 kg·m^{2}·h^{1}.
The studies demonstrated that the presence of finely dispersed water results in a significant increase in heat transfer compared with the case of using purely aircooling. With a fixed heat flux, the energy efficiency increases with increasing water concentration, reaching the values over 600 W·m^{2}·C^{1} at 111.68 kg·m^{2}·h^{1}, which is 2.8 times higher than for air cooling. With further development of research in order to clarify the optimal areas of intensification, it is possible to use this technology to intensify heat transfer to the air in dry cooling towers of nuclear power plants and thermal power plants used in hot and extreme continental climates.
Nuclear power plants, aerosol cooling, airwater mist, water concentration, heat exchange intensification.
Advanced nuclear power plants are equipped with passive systems for emergency decay heat removal from reactor equipment (PEHRSs) to the final heat absorber (ambient air) in case of development of accidents accompanied with primary cooling circuit leakage and for transferring heat to the final heat absorber (ambient air) (
According to Dalton’s law, the rate of evaporation from the surface of a body W in kg/(m^{2.}s) is directly proportional to the moisture deficit calculated using the temperature of the evaporating surface and inversely proportional to the atmospheric pressure (
W = A (X_{T, sph} – X_{T,in})/P, (1)
where X_{T}, _{sph} is the elasticity of saturation taken using the temperature of the evaporating surface, GPa; X_{T},_{in} is the elasticity of vapor in the ambient air, GPa; Р is the atmospheric pressure, GPa; А is the proportionality factor dependent on the air flow velocity.
It is clear from Dalton’s law that the higher is the difference X_{T},_{sph} – X_{T},_{in} the higher is the evaporation rate. If the evaporating surface is warmer than the air, then X_{T},_{sph} exceeds X_{T},_{in}. In such case, evaporation continues even when air is saturated with water vapor.
General layout of the test stand for investigation of processes of heat dissipation from the surface by airwater mist is presented in Fig.
General layout of the experimental stand: 1 – plexiglass tube; 2 – support; 3 – copper sphere; 4 – blending chamber; 5, 15 – differential pressure gauge; 6 – fan; 7 – tank; 8 – ultrasonic mist generator; 9 – multimeter; 10, 17 – voltage controller; 11 – analogue output module; 12 – data collection module; 13 – computer; 14 – Pitot tube; 16 – air pump.
The test stand consists of the measuring section with the object under measurement (spherical element with embedded heater), air supply, mist generation, phase blending and data collection subsystems (temperature, water and air flowrates, pressure, etc.).
Measurement section is the transparent channel with 50 mm diameter, 2 mm wall thickness and total length of 940 mm. The measured object – a sphere with 34 mm diameter made of copper – is mounted inside the channel and is supported from the rear side with a fabricreinforced laminate rod with 4 mm diameter.
Air supply subsystem incorporates electric vacuum air blowers (3000 rpm, 1000 W) with air flowrate controller, compact Pitot head equipped with digital micromanometer for measuring average airflow velocity.
Spherical element is heated using alternate current electrical heater installed in the stainless steel casing with 8mm diameter and 31mm length. Electrical heater power is controllable for achieving the required heat flux on the surface of the sphere using RNO2501000 autotransformer and is measured using APPA 109N digital multimeter with ± 0.06% precision.
Special hightemperature heatconducting paste is used for reducing thermal resistance between the metal of the sphere and the electrical heater.
Schematic layout of the spherical element with heater and thermocouples is shown in Fig.
All thermocouples are connected to the data collection system consisting of the analogue signal input module of OWEN MV1108A type, MSD200 data registration device and personal computer.
Airwater mist is generated as the result of injection of microsize water droplets in the airflow inside the blending chamber. Ultrasonic mist generator with piezoelectric transducers and generator operated at 1.7MHz frequency (
Mass flow rate of the sprayed water was determined according to the change of water mass in the tank of the mist generation subsystem per unit time. Mixture of finely dispersed water with controlled amount of air is formed in the blending chamber. Application of the process of air humidification under discussion here allows obtaining airwater mist with low mass humidity varying in the course of experiments within the range from 0.1 to 6%.
The precision of measurements of calorimeter temperature amounted to 4%. Uncertainties of measurements of air and water flowrates, heater power and geometrical characteristics ranged from two to 3%. Multiple repetitions of measurements at the preset regimes produced statistical dispersion of results within 5%.
The following dimensionless criteria were used in the processing of experimental results:
Uncertainties of calculation of Nusselt and Reynolds numbers amount to less than 10%, which allows reliable detection of the effect of heat transfer intensification.
Equation for the sum of convection Q_{conv}, conductive Q_{cond.} components of heat transfer and evaporation of water droplets Q_{ev} on the heated surface (
Q = Q_{conv.} + Q_{cond.} + Q_{ev.}, (2)
Q _{conv.} = aF_{sph}(T_{av.sph} – T_{in}), (3)
Q _{cond.} = M_{d}C_{p}_{,w}(T_{av.sph} – T_{in}), (4)
Q _{ev.} = aF_{wet.} (X_{T}_{,sph} – X_{T}_{,in})r/C_{p}_{,w}, (5)
where a is the heat transfer coefficient; F_{sph} is the surface of the sphere; T_{av.sph} is the average temperature of the surface of the sphere; T_{in} is the inlet temperature of the mist; M_{d} is the water mass flowrate; C_{p}_{,w} is the water heat capacity; F_{wet} is the area of the wetted surface of the sphere; X_{T}_{,sph} is the concentration of moisture for the temperature of the sphere; X_{T}_{,in} is the concentration of moisture at the temperature at the mist inlet in the channel.
It is assumed that part of water precipitates on the surface in the process of the mist flowing around the heated sphere. Part of water droplets evaporates, while some of them remain scattered in the main airflow. Experiments were conducted with pure dry air and with airwater mist within the range of water flow densities (j = 23.39, 46.79, 88.31, and 111.68 kg·m^{–2}·hr^{–1}). Range of Reynolds numbers is 2500 – 55000.
Temperature of the surface is an important factor exercising significant influence on the process of heat transfer to airwater mist. When surface temperature is very high, finely dispersed water can completely evaporate before it attains the complete heated surface. On the other hand, wetting of the heated surface and formation of liquid film may take place at low surface temperature. Several values of heat flux were examined for investigating the effects of surface temperature on the heat transfer coefficient. Average temperature of the surface is shown in Fig.
Surface temperature decreases with increasing Reynolds number in the case of singlephase air flow. Because of evaporation of water, surface temperature decreases in the airwater mist with increasing water flow density at constant heat flux. Surface temperature decreases by 4, 17, 37 and 47% as compared to air because of the increase of water flow density. It was found using the methods of direct observation and photography and video filming that the heated surface is covered with thin water film at high density of liquid flow and the values of Re > 17500. At the same time heated surface remains dry at Re < 17500 for all values of heat flux.
Variation of average heat transfer coefficient in the form of dependence Nu = f (Re) is shown in Fig.
It can be noted that average Nusselt number significantly increases with increase of Reynolds number for any water flow density. Value of Nu increases in singlephase flow from 66 to 240 within the range of Re from 2500 to 55000 at constant heat flux (q = 10155 W/m^{2}). Within the investigated range of Reynolds numbers Nusselt numbers for airwater mist are by 1, 19.7, 90.2 and 134% higher than in the air flow for water flow density j = 23.39, 46.79, 88.31 and 111.68 kg·m^{–2}·hr^{–1}, respectively.
High Reynolds numbers correspond to high flow velocity and, consequently, in airwater mist more and more water precipitates from the mist on the surface of the sphere and absorbs significant energy in the process of evaporation.
Finely dispersed water can create filmy layer and enhance the process of heat transfer by increasing contact surface with the main flow. Reduction of Nusselt number due to evaporation of water from the mist without direct wetting of the heated surface occurs with increasing surface temperature (Re < 17500).
Relative change of heat transfer to the mist depending on the aircooling is presented in Fig.
Within the range of small water flow density (j varying from 23.39 to 46.79 kg·m^{–2}·hr^{–1}) with q = const. the factor characterizing intensification of heat transfer undergoes contradictory variation with increasing Re due to the relative reduction of water content in the airwater mist.
Within the whole range of variation of values of moisture concentration there exists the area of optimal intensification of heat transfer determined by the correlation between the mechanisms of evaporation and precipitation of moisture on the surface of the sphere. Mass concentration of moisture in the mist within this optimal zone falls within the range from 0.3 to 0.6%. Dependence of optimal values of the heat transfer intensification factor on the air flowrates (Re number) and on mass moisture concentration x is represented for d/D = const by the following equation:
K ^{max}= 1+ Re^{1.18}×x^{2.2}. (6)
The above equation is valid within the range of Reynolds numbers from 32000 to 60000 and water flow densities from zero to 111.68 kg·m^{–2}·hr^{–1}.
Stabilization and smooth reduction of the degree of intensification is observed for j > 80 kg·m^{–2}·hr^{–1} and high flow velocities of airwater mist after passing through the optimal value. For j < 20 kg·m^{–2}·hr^{–1} intensification of heat transfer is practically nonexistent within the whole range of investigated velocities of the mist flow and heat flux values.
Test facility and methodology for investigation of heat transfer characteristics were developed for the case when the heat transfer surfaces of the model of PEHRS heat exchange element are cooled down by airwater mist. Investigation of heat transfer was performed within wide range of heat flux values, Reynolds numbers and concentrations of water droplets for the case of cooling of object with airwater mist. The following conclusions can be drawn based on the experimental results obtained.
Use of airwater mist is an efficient method for cooling heat exchange surfaces. Temperature of the surface decreases by 4, 17, 37 and 47% as compared with cooling by air flow for water flow densities j = 23.39, 46.79, 88.31 and 111.68 kg·m^{–2}·hr^{–1}, respectively.
Average Nusselt number significantly increases with increased Reynolds number for all values of heat flux and water flow densities. When finely dispersed water is used, the average Nusselt number is found to be 11, 19.7, 90.2 and 134% higher than for the case of cooling with air flow for water flow densities j = 23.39, 46.79, 88.31, 111.68 kg·m^{–2}·hr^{–1}, respectively. Average Nusselt number increases with increased heat flux for fixed values of water flow density and Reynolds numbers.
The value of the factor characterizing heat transfer intensification in the case of cooling with airwater mist flow is significantly higher (up to 300%) than for the case of air cooling and has area of optimal moisture concentration for each value of water flow velocity.
The present study was implemented under the financial support in pursuance with Resolution No. 211of the Government of the Russian Federation under Contract No. 02.A03.21.0006.