Subcooled Water Research Articles

Numerical study on effects of nozzle hole spacing distance on steam submerged jet condensation from triangular holes spray nozzle

Steam-water direct contact condensation (DCC) is a phenomenon which is frequently found in the industrial sector due to its rapid and efficient heat and mass exchange rates. In this work, three dimensional numerical study of steam injection from triangular hole nozzles (with different holes pitch circle diameter, PCD = 6, 8,10 and 12 mm) into subcooled water tank has been carried out using Eulerian multiphase flow model along with realizable k-epsilon turbulence model. The thermal phase change model as User Defined Function (UDF) was employed in Fluent software for capturing submerged steam jet condensation. The simulation results were validated using experimental results as well as previously published numerical study and a fairly good agreement was observed. The effects of nozzle hole spacing distance on distributions of various thermal hydraulics parameters including velocity, pressure, temperature, heat and transfer characteristics have been studied. The axial distribution of the velocity, temperature, heat transfer coefficient and mass transfer computed at the centerline of the nozzle was found weaker than observed at the centerline of the nozzle holes. However, these distributions were not weak at nozzle centerline for nozzle with PCD 6 mm. The heat transfer coefficient has been found in the ranges of 1.75–2.3 MW/m2 K. The results reveal that the heat transfer coefficient for nozzle PCD 8 mm is higher than rate obtained for other nozzles. The rate of mass transfer at the holes centerlines has been found increasing with the holes PCD. The axial distribution of the mass transfer curve shifts downstream for nozzle with smaller PCD. Additionally, the mutual shielding effect of the jets weakens as holes PCD increases.

Passive mitigation of condensation-induced water hammer by converging-diverging structures for offshore nuclear power plants

The equal-height-difference passive heat removal system (EHDPHRS) utilizes the ocean environment as an infinite heat sink that can better suit ocean nuclear power plants (ONPPs). However, a structure with an equal-height-difference can easily cause reverse flow of subcooled water, leading to serious condensation-induced water hammer (CIWH). In this paper, a fully visual experimental facility with same geometric characteristics as the EHDPHRS was established with the aim of studying the CIWH phenomenon. First, two kinds of CIWH modes were observed by the results of the visual experiment, namely the intermittent CIWH and the continuous CIWH. The typical feature of the intermittent CIWH is that the single-phase natural circulation occurs intermittently. Then, the converging-diverging structure (CDS) was designed to avoid the CIWH events in the EHDPHRS. The effectiveness of CIWH prevention was evaluated employing the velocity field, which was obtained by the PIV experimental system. It was found that the CDS can effectively inhibit the further reverse flow of subcooled water and thus eliminate the triggering mechanism of CIWH events. Finally, a dimensionless number was used to balance CIWH prevention and natural circulation capacity. It was shown that when the dimensionless number is in the range of 0.67–0.80, the CDS can eliminate the CIWH events while ensuring the natural circulation capacity.

Numerical study on the interfacial heat and mass transfer characteristics of bubble direct contact condensation with the effect of non-condensable gas

Bubble condensation is one of the fundamental issues in the direct contact condensation, which is also affected significantly by the existence of non-condensable gas. In this study, a further investigation is carried out on air-steam bubble condensation with the Proportion Integration Differentiation (PID) control model which adjusts the mass transfer coefficient of Lee model in each time step by comparing the parameters obtained in the simulation with the ones calculated based on the empirical correlations originated from the experiments. Furthermore, species transport equations are employed to obtain the local mass fractions of steam and air. The proposed simulation schema is verified by the existing experimental data.Results show that the air-steam bubble condensation can be divided into initial and stable stages. The mass transfer coefficient used in Lee model can be regarded as linear growth over time at each stage, but with different slopes. The non-condensable gas can weaken the heat and mass transfer of condensation, but stabilize the bubble deformation. The present study also reveals the influences of subcooling and initial bubble diameter on the air-steam bubble condensation, which is expected to be instructive for better understanding the effect of non-condensable gas on the interfacial characteristics of bubble condensation in subcooled water.

Mechanisms and inhibition method investigation on condensation induced water hammer in an equal-height-difference passive heat removal system

The equal-height-difference passive heat removal system utilizes the ocean as an ultimate heat sink for residual heat removal, demonstrating obvious advantages for ocean nuclear power plant applications. However, the structure with a zero-height difference can trigger the reverse flow of subcooled water. This could cause serious condensation induced water hammer events. In this study, a fully visual experimental setup was established to investigate the occurrence features. First, two kinds of formation mechanisms were defined by the visual experimental results. In this regard, the subcooling and the isolated steam slug are the necessary conditions for the formation and development. In addition, the occurrence positions were found to be concentrated in the test pipe section near the water tank. Then, a butterfly check valve was designed to avoid the condensation induced water hammer events based on the formation mechanisms. The experimental results showed that the check valve usage can effectively inhibit the reverse flow of the subcooled water and thus achieve its purpose of eliminating condensation induced water hammer events. Moreover, the check valve transforms the two-phase natural circulation flow from the condensation oscillation mode to the stable boiling two-phase flow mode. As a result, it enhances the two-phase natural circulation capacity.

Numerical simulation of air discharged in subcooled water pool

Turbulent jet discharges in subcooled water pools are essential for safety systems in nuclear power plants, specifically in the pressure suppression pool of boiling water reactors and In-containment Refueling Water Storage Tank of advanced pressurized water reactors. The gas and liquid flow in these systems is investigated using multiphase flow analysis. This field has been extensively examined using a combination of experiments, theoretical models, and Computational Fluid Dynamics (CFD) simulations. ANSYS CFX offers two approaches to model multiphase flow behavior. The non-homogeneous Eulerian-Eulerian Model has been used in this work; it computes global information and is more convenient to study interpenetrated fluids. This study utilized the Large Eddy Simulation Model as the turbulence model, as it is better suited for non-stationary and buoyant flows. The CFD results of this study were validated with experimental data and theoretical results previously obtained. The figures of merit dimensionless penetration length and the dimensionless buoyancy length show good agreement with the experimental measurements. Correlations for these variables were obtained as a function of dimensionless numbers to give generality using only initial boundary conditions. CFD numerical model developed in this research has the capability to simulate the behavior of non-condensable gases discharged in water.

An experimental study on quenching of metallic spheres in seawater

Motivated by quenching of melt droplets and debris particles in seawater during a hypothetical severe accident of light water reactors (LWRs), the quenching process of stainless-steel spheres in a seawater pool was investigated in the present study. The polished spheres were pre-heated up to 1000℃ in an induction furnace with inserted atmosphere, and then immersed into the subcooled water pool in a chamber made of transparent quartz. A thermocouple was embedded in the center of the sphere to measure the history of the sphere’s temperature, while a high-speed camera was employed to record the quenching process and vapor film dynamics. Quantitative data, e.g. film thickness and oscillation, of the vapor film evolution during the quenching process were obtained through an image processing program developed on the MATLAB platform.The experimental results indicated that the quenching rate was higher in seawater than in deionized water, and the vapor film collapsed at a temperature higher than the Leidenfrost temperature in deionized water. The trend appeared more significant with increasing subcooling of water. The comparison of the quenched spheres suggested the surface of the sphere in seawater achieved higher degree of discoloration and roughening than that in deionized water, probably due to the additives of salt which change water properties. The image processing and analysis revealed that the vapor film had different thickness profile along the upper and lower hemispheres, and the averaged film thickness was smaller in seawater than in de-ionized water during the stage of stable film boiling. The vapor film was thinning and oscillating with time, and its fluctuations appeared different frequencies and amplitudes at the upper and lower locations, which may explain the mechanism of the earlier collapse of vapor film in the quenching process of a high-temperature sphere in seawater.

CFD Simulation on Pressure Profile for Direct Contact Condensation of Steam Jet in a Narrow Pipe

In the published experimental results, it has been observed that when high-speed steam spurt into the subcooled waterflow, the total pressure along the axial direction at trailing edge of the steam plume shows a pressure-lift. To reveal the mechanism behind this phenomenon, this study utilizes a particle model to investigate the pressure profile of steam jet condensation in subcooled water flow in a narrow pipe. A numerical model based on the Eulerian–Eulerian multiphase model has been developed to accurately simulate the characteristics of pressure profile along the axial direction. The model’s validity is established by comparing the steam plume shape and temperature profiles with the experimental data. By analyzing the total pressure profile of the axis and the contours of gas volume fraction, it is found that there exists a pressure-lift phenomenon at trailing edge of the steam plume. The dynamic pressure of the water also shows a pressure-lift at this position, so it can be inferred that the dynamic pressure of the water is the main factor of the total pressure-lift.

Numerical study on the evolution characteristics of phase interface of subsonic and sonic steam jets in transverse subcooled flow

The cooling process of high-temperature steam ejecting into subcooled liquid widely exists in industrial production, with high mass and heat transfer efficiency. This paper establishes an unsteady numerical simulation method for jet condensation based on the two-fluid model and the thermal phase change model, focusing on the phase interface evolution and pressure fluctuation characteristics of steam jets in the crossflow of subcooled water under different mass fluxes. With the increase of mass flux, the ability of the jet plume to resist water flow is enhanced; thus, the heat and mass transfer at the two-phase interface becomes more stable; the length of the plume increases, and the range of oscillation of the plume length in the axial direction and the radial direction decreases. The pressure vortex that generated by the two-phase region and the subcooled water is more stable with reduced total pressure fluctuations. Under subsonic conditions, bubbles are shed on the leeward side, and the ranges of axial oscillation, radial oscillation and total pressure oscillation are 1.57de–3.41de, 0.21de–0.37de, and 0.17de–0.32de, respectively. The steam plume becomes progressively longer in the axial direction as the steam velocity increases from subsonic to sonic and supersonic speeds, and the range of interfacial oscillations and total pressure fluctuations along the radial direction are further reduced. Under the same working condition, with the increase of axial distance, the momentum of the steam jet decreases, the ability to resist the impact of water flow decreases, and the windward side is greatly affected; the axial and radial amplitude fluctuations of the pressure increases, with the amplitude on the windward side being significantly larger than that on the leeward side, and the maximum difference between the fluctuations on both sides increases continuously. The momentum change of the two-phase interface is intense, which is the main factor of pressure fluctuation. The correlation coefficients of pressure oscillation with radial and axial interface oscillations are 0.9110 and 0.7477, respectively.

Thermal stability of 6-filament MgB2 wire with resistive CuNi sheath cooled by liquid He and water ice

In this paper, thermal stability of two 6-filament MgB2 wires with resistive CuNi sheath has been examined by DC transport measurements at the temperature of liquid He and variable external magnetic fields, and at variable temperatures of water ice under self-field conditions. Current-voltage characteristics up to the quench were performed and the operation limits of both wires under different cooling are analyzed and compared. The obtained results show that only 6 % of well conducting copper core allows to improve the thermal stability considerably, and it also limits the wire heating after the quench. MgB2 wires showed more stable operation after the quench under subcooled water ice in comparison to traditional Helium bath cooling, which offer a cheap and safe cooling mode for future devices.

Three-dimensional numerical simulation of the interaction between high-pressure sub-cooled water jet impacting high-temperature liquid lead-bismuth metal

The rupture accident of the steam generator heat transfer tube of the lead-bismuth cooled fast reactor causes a high-pressure sub-cooled water jet impacting high-temperature liquid lead-bismuth, and the interaction mechanism between water and liquid lead-bismuth needs to be studied clearly. In this paper, based on the volume of fluid method, the LES turbulence model and the Lee phase transition model, a three-dimensional numerical model is developed to study the multiphase flow and heat transfer mechanism of water and lead-bismuth interaction. The effects of sub-cooled water temperature, inlet pressure, and nozzle diameter on the distribution of phase state, pressure, temperature, and generated steam volume are studied in detail. The results indicate that the water jet can be divided into the jet water zone, the water-steam mixed zone, and the multiphase flow zone, respectively. Among them, the appearance of water-steam mixed zone depends on the internal boiling of the jet, and this zone tends to appear under the condition of low inlet pressure, small nozzle diameter and high sub-cooled water temperature. For multiphase flow zone, the temperature increases with sub-cooled water temperature, decreases with nozzle diameter and is almost not influenced by inlet pressure; The pressure increases with sub-cooled water temperature, inlet pressure and nozzle diameter. When nozzle diameter is large enough, however, it no longer has any effect on the temperature and pressure distribution within the jet. The volume of steam generated by water boiling is positively correlated with all three factors. Overall, the steam migrates upward through the outer side of the jet, absorbing heat and expanding during the migration process, and forming a high-temperature area at the top of the jet.

Study on thermodynamics and flow characteristics of jet condensate in subcooled flowing liquid

Steam jet condensation in flow water under different conditions was investigated numerically. A cavitation model considering thermodynamic effects was proposed to study the flow and interfacial heat transfer characteristics. The distribution of thermodynamic parameters and the influence of different incident angles on the jet state were studied. The decrease in steam pressure and water subcooling would have a negative impact on the heat transfer coefficient of the interface and the length of plume. Furthermore, with increased water flow, the steam plume was tapered and its length was gradually shortened. Additionally, the turbulent kinetic energy is the maximum at the outlet of the nozzle and decreases with increasing flow distance. The distribution characteristic of the heat transfer coefficient is similar to and turbulence kinetic energy, which increases with the water mass flux, and reaches the maximum at both sides of the axis. Gradually increases in the middle of the axis. Therefore, the improvement of jet condensation should be achieved by adjusting the incidence angle and the water flow velocity.

Study on the Dynamic Generation of Subcooled Water Using a Compact Heat Exchanger

The dynamic generation of ice slurry from subcooled water is one of the most promising ways to make ice; this process is utilized widely in ice storage air-conditioning systems. However, the random occurrence of ice blockage during the generation of subcooled water using conventional heat exchangers prevents the increase in subcooling, thereby reducing the efficiency of the release of the subcooled water and converting it into ice slurry. A more efficient approach to reducing the fluid passage time is to employ a compact heat exchanger with a highly efficient heat transfer performance, a heat transfer length of only 21.5 mm, and a hydraulic diameter of 0.32 mm. A compact heat exchanger was used to build a dynamic generation setup for subcooled water, and 40 wt% of non-freezing liquid and tap water was used as the working fluid for heat exchange to generate subcooled water. The results show that the compact heat exchanger can achieve a greater subcooling degree (3.8 K) and longer duration (108 min). This study further explored the potential for dynamic ice making from deep, subcooled water and improved the overall structure of the compact heat exchanger used. The experimental setup is recommended based on the analysis of the results.

Experimental study on transient flow characteristics in an equal-height-difference passive heat removal system for ocean nuclear power plants

The equal-height-difference passive heat removal system (EHDPHRS) utilizes the ocean as the heat sink for infinite heat dissipation and has broad application prospects in ocean nuclear power plants. However, the equal height difference between the heat source and heat sink is susceptible to severe flow instability. In this study, a fully visualized experimental facility was constructed to analyze the transient flow characteristics in the EHDPHRS. The experimental results revealed three typical two-phase natural circulation flow modes in the EHDPHRS. The formation mechanism of each flow mode was analyzed using thermal–hydraulic parameters. The results revealed that the flow mode is determined by the steam–water direct contact condensation and reverse flow of subcooled water, and both factors promote each other. Moreover, condensation-induced water hammer (CIWH) was found to be the most severe flow instability in the EHDPHRS. CIWH events can lead to violent pressure peaks and reduce the natural circulation capacity. In addition, when the occurrence position of CIWH events is concentrated at the pipe outlet, the influence of CIWH events on the two-phase natural circulation flow capacity gradually decreases. These research results aid in the optimal design of the EHDPHRS and prevention and control of CIWH events.

Gradient Heatmetry in Study of Heat Transfer during Pool Boiling of Subcooled Water and Liquid with the Addition of Al2O3 Microparticles

Gradient Heatmetry in Study of Heat Transfer during Pool Boiling of Subcooled Water and Liquid with the Addition of Al2O3 Microparticles

Experimental investigation on condensation regimes and transition boundary during bubble condensation in narrow rectangular channel

Direct steam contact condensation in subcooled water is widely used in various industries because of high efficiency in heat and mass transfer. A visual experiment is carried out to investigate the condensation process of vapor bubble in narrow rectangular channel. The steam injection mass flux, subcooled water temperature and water mass flow are 0.43–40 kg/m2·s, 50–93 °C, and 0–200 kg/m2·s, respectively. In order to obtain more accurate parameters of vapor bubble in narrow rectangular channel, a new bubble volume model was developed and compared with existed models. A three-dimensional condensation regime map is developed with the coordinates of steam injection flux, water mass flow, and water temperature. Five condensation patterns such as chug, oscillatory bubble, growth bubble, transition region and oscillatory jet are determined based on the dynamic behavior of steam-water interface. The bubble equivalent diameter and aspect ratio during evolution process are calculated, and these two parameters in each pattern exhibit a significant periodicity, except for the oscillatory Jet. The true area of heat transfer interface between bubble and subcooled water is much larger than the actual measured value under high steam injection mass flux. A transition criterion from oscillatory bubble regime to oscillatory jet regime is proposed by the analysis of the bubble critical diameter and the condensation heat transfer of steam-water interface.

Dropwise and liquid-jet laminar flow of subcooled water falling over horizontal tube banks

Falling-film heat exchangers are widely used to promote heat and mass transport within industrial environments. In contemporary condensers, for example, liquid films continuously form over a tube row and impact multiple tubes beneath, resulting in a rich variety of two-phase flow phenomena, including the classical Rayleigh-Plateau instability. It is widely known that altering the contact angle of the tubes influences the heat exchangers’ heat and mass transfer characteristics by changing the falling films’ flow behaviour (i.e. dropwise or jet). Therefore, resolving interfacial flow characteristics, such as the Rayleigh-Plateau instability, is crucial for falling film studies as this instability may change the state of the falling film from being jet flow to dropwise/jet flow or even just dropwise. To explore this aspect, three-dimensional direct numerical simulations are used in the present work to replicate the jet instability behaviour. By employing certain mesh refinement strategies demonstrated from various case studies, the interfacial flow behaviour and jet instability have been captured, showing good agreement with experimental and theoretical data previously published. This is determined from the increased jet bifurcation and a maximum error of 2.2% of the film thickness in comparison to Nusselt's result within the bulk region of the tube. In addition, the average and local heat transfer coefficients across the first tube were compared with numerical solutions to give an average error of 3.2 and 14 %. Jet instability is found to be directly related to the amplitude of the axial pressure variations which is augmented by the growth rate of the capillary wave that travels upstream. Specifically, for contact angles of 30 and 60°, the flow is in a dropwise and dropwise/jet state due to the emergence of the Rayleigh-Plateau instability. Notably, the additional satellite droplets that materialize from the unstable jet lead to an increase in the local heat transfer coefficient of the film compared to without the Rayleigh instability improvement.

Investigation on heat transfer and flow of a sonic steam jet condensation in crossflow of water in pipes

Direct contact condensation is widely used in nuclear reactor thermal-hydraulic systems, underwater propulsion systems, liquid propellant rocket engines and other industrial fields due to its highly efficient heat and mass transfer. Understanding the distribution of the thermal-hydraulic parameters of steam jet condensation systems has become an urgent need for optimal design and safe operation of industrial processes. In this paper, the flow characteristics and heat transfer for the condensation of sonic steam jet in crossflow in a subcooled water pipe are numerically simulated based on experimental results. A two-resistance model under the framework of two-fluid model is employed to consider the mass and heat transfer across the steam-water interface. The steam plume transforms from contraction-cone into expansion-cone when steam mass flux rises from 500 to 900 kg/m2s, and the plume length gradually increases approximately linearly from 2.43de to 4.17de. When the steam mass flux is 900 kg/m2s, the expanding channel of the steam plume allows the jet to accelerate from sonic to supersonic with Ma = 1.43. At the transition of expansion to contraction, the steam dynamic pressure reaches a maximum and the steam total pressure decreases to a minimum. Because the water dynamic pressure is significantly higher than the steam dynamic pressure, the total pressure is maximum at the end of the two-phase region. Mass transfer is much higher in two-phase region than in other regions. The heat transfer coefficient calculated from the interface area decreases from 1.98 to 1.56 MW/m2K, and that calculated from the two-phase region volume increases from 3.37 × 103 to 4.05 × 103 MW/m2K.

Modeling of horizontal steam injection in a water pool

A comprehensive understanding of the direct contact condensation (DCC) phenomena of steam discharged into a subcooled water pool is essential for the nuclear engineering field. DCC is the most rapid and efficient way to condense high-pressure and high-temperature steam with simple engineering systems. However, the pressure oscillation due to the rapid condensation and small length and time scales of turbulent two-phase flow make the investigation of DCC phenomena challenging either by experiments or with numerical simulations.This paper presents the computational fluid dynamics (CFD) simulations of the small-scale separate effect test facility (SEF-POOL) experiment of Lappeenranta-Lahti University of Technology (LUT University). In the SEF-POOL test, the steam was injected using the orifice of a diameter of 16 mm. All simulations were performed by employing the compressible two-phase solver of OpenFOAM which was based on the Eulerian–Eulerian two-fluid approach. The interfacial heat transfer between steam and water was modeled by using the Nusselt number formulation of Coste (2004). The Rayleigh–Taylor interfacial (RTI) area model of Pellegrini et al. (2015) was applied for the interfacial area modeling. In this work, the formation and collapse of the steam bubbles are studied using the extended pattern recognition based image analysis algorithm. The pattern recognition (PR) algorithm is based on the video material recorded during the SEF-INF2 experiment of the SEF-POOL test facility.The presented study show that if the two-phase interface is rough and smear into the subgrid scale, the interfacial area density calculation solely with the volume fraction gradient is inadequate for predicting rapid condensation rates. Results indicate that rapid and high interfacial accelerations in DCC trigger interfacial instabilities which increased the interface roughness and total interfacial area consequently. Thus, interfacial area modeling with the RTI increased the total interfacial heat transfer rate in the simulations. The achieved results exhibit that the selected modeling approach captured the rapid bubbling condensation oscillation mode of the SEF-POOL test. Presented results demonstrate that image analysis and pattern recognition have great potential in thermal-hydraulic research that can be beneficial for direct contact condensation and phase-interface dynamics model development.

Experimental study on condensation and acoustic characteristics of steam underwater injection

Introduction: Bubbling deaeration is widely applied in marine steam power system, and direct contact condensation (DCC) of steam injected into water under low sub-cooled water is the thermodynamic foundation of this technology. The deaeration process involves complex thermal and hydraulic behaviors, such as violent turbulence, interfacial mass transfer, and multiphase flow etc., which induce significant pressure fluctuation and radiated noise. In this paper, the two-phase flow regimes and acoustic characteristics of DCC of steam injected into water were studied for sub-cooled temperature ranging from 3℃ to 17℃.Methods: Using high-speed camera to capture the dynamic plume shape, four typical plumes are found: external chug with encapsulating bubble, oscillating bubble jet, conical and cylindrical jet and cap jet; The underwater acoustic signal is measured synchronously, and the acoustic spectrum is obtained by Fast Fourier transform (FFT).Results and Discussion: Combined with the plume image, the component sources of the acoustic signal are analyzed. It is found that the first dominant frequency and the second dominant frequency of the acoustic are caused by the evolution of the jet main steam zone and the downstream large steam column respectively; By analyzing the influence of injection parameters on acoustic signal, it is found that steam flow and sub-cooled temperatures have a great influence on acoustic root mean square value and spectrum characteristics. With the decrease of sub-cooled temperature, the main frequency of injection noise signal gradually moves to low frequency. Low frequency noise dominates the radiated noise when the flow field approaches the saturation temperature.

Experimental study on the air-steam bubble condensation characteristics in subcooled water

Experimental study on the air-steam bubble condensation characteristics in subcooled water