CFD analysis of a turbulent jet behavior induced by a steam jet discharge through a single hole in a subcooled water tank

A computational fluid dynamics (CFD) analysis for a turbulent jet flow induced by a steam jet discharged into a subcooled water pool was performed for 10 s of transients to investigate whether the currently available CFD codes can be suitably used as a tool to investigate the applicability of the existing semi-analytical correlations to a condensing jet-induced turbulent jet and to analyze the thermal-hydraulic behavior, such as global circulation and local hot spot, in a condensation pool for advanced light water reactors. As for the numerical experiment, a series of sensitivity calculations was conducted systematically to elucidate the major factors which can cause different analysis results by varying the mesh distributions, numerical models for a convection term and an eddy viscosity term. The effect of a difference in the velocity and the temperature distribution in a region between the sparger and the pool wall has not been observed in the afore-mentioned sensitivity calculations. The comparison of the CFD results with the test data shows that the CFD analysis does not accurately simulate the local phenomenon of a turbulent jet existing downstream of a steam jet. It was found that the value of the turbulent intensity at the inlet of the turbulent jet region is the most important factor because it can determine the boundary of a turbulent jet through a momentum diffusion process in a radial direction. The comparison of the CFD results with the test data shows that the CFD analysis can accurately simulate the local phenomenon of a turbulent jet existing downstream of a steam jet only when the CFD analysis reflects the physics of a turbulent jet.     

CFD simulation of steam jet-induced thermal mixing in subcooled water pool

This study examined the IRWST thermal mixing phenomena induced by a steam jet in a subcooled water pool. Due to the limitation of the current CFD code to simulate condensation, the steam condensation region model was developed to evaluate the thermal mixing phenomena. Within this region, all the steam was condensed into water, and the steam mass and energy inputs were treated as the source. This calculation was treated using single-phase CFD methods. The benchmark calculation for a thermal mixing experiment in the water tank was performed to develop an optimized 3D evaluation methodology of the thermal-hydraulic behavior in APR1400 IRWST. Steam discharge through the sparger and condensation phenomenon was modeled with the choking flow and thermal mixing model in the quenching tank using CFX11.Three types of thermal mixing experiments, local phenomena test, thermal mixing tests in cylindrical water pool and annulus water pool, were designed to provide data representative of the behavior of the prototype for CFD simulations of the thermal-hydraulic behavior in IRWST. A comparison of the calculated and experimentally measured temperature profiles showed some disagreement particularly around the sparger. The main reason for this disagreement was caused by the difference in the test and simulating conditions at the tank wall. However, moving away from the sparger, the trends of the temperature rise became similar to that in the experiment. Despite these problems, this model is the best way of evaluating the thermal mixing phenomena caused by a steam jet in a subcooled water pool.     

Experimental study on a steam-driven turbulent jet in subcooled water

This study was designed to investigate the behaviour of a steam-driven turbulent jet in subcooled water. The total pressures along the axis of the jet were measured under various steam mass fluxes, water temperatures and nozzle sizes. The results indicated that the lifting-pressure coefficient decreased with inlet steam pressure and water temperature, and an optimal lifting-pressure coefficient could be obtained at pressure ratio 0.228. Then an empirical correlation based on the influencing factors was given to fit the position of the maximum pressure. Moreover, the velocity profiles in the turbulent jet region were obtained, and the results showed the self-similarity property of the velocity profiles. An empirical correlation was suggested to correlate the jet half width obtained from the velocity profiles.     

CFD analysis for thermal mixing in a subcooled water tank under a high steam mass flux discharge condition

A Computational Fluid Dynamics (CFD) analysis for a thermal mixing test was performed for 30 s to develop the methodology for a numerical analysis of the thermal mixing between steam and subcooled water and to apply it to Advanced Power Reactor 1400 MWe (APR1400). In the CFD analysis, the steam condensation phenomenon by a direct contact was simulated by the so-called condensation region model. Thermal mixing phenomenon in the subcooled water tank was treated as an incompressible flow, a free surface flow between the air and the water, and a turbulent flow, which are implemented in the CFX4.4. The comparison of the CFD results with the test data showed a good agreement as a whole, but a small local temperature difference was found at some locations. A sensitivity analysis was performed to find the reason of the temperature difference. The commercial CFD code of CFX4.4 together with the condensation region model can simulate the thermal mixing behavior reasonably well when a sufficient number of mesh distributions and a proper numerical method are selected.     

Condensation regime diagram for supersonic/sonic steam jet in subcooled water

This study was designed to determine the behaviour of the flow patterns of supersonic and sonic steam jet condensation in subcooled water. The effects of steam mass flux, water temperature and pressure ratio on the flow pattern were discussed. The results indicated that the flow pattern was not only affected by steam mass flux and water temperature, but also affected significantly by pressure ratio. The expansion and contraction ratios of the steam plume were influenced significantly by pressure ratio, and a correlation for predicting the expansion and contraction ratio was given based on the theoretical expansion ratio of sonic steam jet. At relatively low steam mass flux and high water temperature, the transition of flow pattern from stable jet to unstable jet was observed, and a criterion of stable-unstable jet transition was given. Moreover, the three-dimensional condensation regime maps, including stable steam plume shapes and unstable jet, for sonic and supersonic steam jet under different conditions were discussed. Four different shapes of the steam plume occurred in sonic steam jet regime map, and six different shapes of the steam plume occurred in supersonic steam jet regime map. The regime maps were also validated against the previous experiments.     

Film boiling from a vertical jet of fuel entering a water pool

The rate of removal of coolant is examined by means of quasi-steady film boiling from a jet of fuel falling into a water pool during a core meltdown accident in a light water reactor. A two-layer turbulence model in the vapor film is used. Entrainment of droplets due to the high vapor velocity is taken into account, which, for a jet of diameter 0.5 m, increases the removal rate by factor of 5 at 0.7 m behind the leading edge. However, the water removal rate is only of the order of kilograms/second, and the total amount of coolant removed is still small compared to the coolant inventory in the lower plenum.     

Analytical tool development for coarse break-up of a molten jet in a deep water pool

A computer code JASMINE-pre was developed for the prediction of premixing conditions of fuel–coolant interactions and debris bed formation behavior relevant to severe accidents of light water reactors. In JASMINE-pre code, a melt model which consists of three components of sub-models for melt jet, melt particles and melt pool, is coupled with a two-phase flow model derived from ACE-3D code developed at JAERI. The melt jet and melt pool models are one-dimensional representations of a molten core stream falling into a water pool and a continuous melt body agglomerated on the bottom, respectively. The melt particles generated by the melt jet break-up are modeled based on a Lagrangian grouped particle concept. Additionally, a simplified model pmjet was developed which considers only steady state break-up of the melt jet, cooling and settlement of particles in a stationary water pool. The FARO corium quenching experiments with a saturation temperature water pool and a subcooled water pool were simulated with JASMINE-pre and pmjet. JASMINE-pre reproduced the pressurization and fragmentation behavior observed in the experiments with a reasonable accuracy. Also, the influences of model parameters on the pressurization and fragmentation were examined. The calculation results showed a quasi-steady state phase of melt jet break-up during which the amount of molten mass contained in the premixture was kept almost constant, and the steady state molten premixed masses evaluated by JASMINE-pre and pmjet agreed well.     

Effect of turbulent natural convection on sodium pool combustion in the steam generator building of a fast breeder reactor

A computational model is proposed to simulate sodium pool combustion considering the effect of turbulent natural convection in a vented enclosure of the steam generator building (SGB) of a fast breeder reactor. The model is validated by comparing the simulated results with the experimental results available in literature for sodium pool combustion in a CSTF vessel. After validation, the effects of vents and the location of the pool on the burning rate of sodium and the associated heat transfer to the walls are studied in an enclosure comparable in size to one floor of the steam generator building. In the presence of ventilation, the burning rate of sodium increases, but the total heat transferred to the walls of the enclosure is reduced. It is also found that the burning rate of sodium pool and the heat transfer to the walls of the enclosures vary significantly with the location of sodium pool.     

Simulating the mixing of helium and air by mixing salt water and fresh water in a scaled enclosure

We are attempting to establish scaling laws to simulate the mixing of helium (a simulant for hydrogen) with air in a large-scale enclosure by mixing salt water and fresh water, in a small-scale enclosure. This will allow us to assess the mixing of gases in a nuclear reactor containment using relatively small-scale liquid-mixing experiments. The scope of our current work does not cover the integrated effects of different mechanisms of gas mixing expected to prevail during postulated LOCA (loss-of-coolant accident) scenarios. The work is limited to mixing caused by jet inertia and buoyancy forces. Within this scope, we have identified the dominant scaling laws, and tested them by conducting gas-mixing experiments in a large-scale enclosure and liquid-mixing experiments in a small-scale enclosure. The experimental results demonstrate the validity of the scaling laws.     

Large-Eddy Simulation study of turbulent mixing in a T-junction

A potential cause of thermal fatigue failures in energy cooling systems is identified with cyclic stresses imposed on a piping system. These are generated due to temperature changes in regions where cold and hot flows are intensively mixed together. A typical situation for such mixing appears in turbulent flow through a T-junction, which is investigated here using Large-Eddy Simulations (LES). In general, LES is well capable in capturing the mixing phenomena and accompanied turbulent flow fluctuations in a T-junction. An assessment of the accuracy of LES predictions is made for the applied Vreman subgrid-scale model through a direct comparison with the available experimental results. In particular, an estimation of the minimal mesh-resolution requirements for LES is examined on the basis of the complementary RANS simulations. This estimation is based on the characteristics turbulent scales (e.g., Taylor micro-scale) that can be computed from LES or RANS simulations.     

Acoustics impact on steam dryer induced by main steam line break in a boiling water reactor

In boiling water reactor (BWR) design, significant acoustic pressure loads impact the steam dryer hood as a result of the main steam line break outside containment (MSLB) event. When a main steam line breaks, it is assumed that the pipe instantaneously breaks completely open to the ambient environment (double-ended guillotine break). Due to the huge pressure difference between the inside of the reactor pressure vessel (RPV) and surrounding ambient environment, a shock wave will form at the break point and burst into the surrounding environment. At the same time, an expansion wave will travel upstream through the main steam line to the RPV, which results in a pressure reduction on the outside of the steam dryer hood. This expansion wave will create a substantial pressure difference between the two sides of the steam dryer hood with a resultant high stress on the hood. This differential pressure load is the acoustic load used in the structure design evaluations for this event. A key design basis requirement for the steam dryer is to maintain structural integrity during transient, and accident conditions. Demonstration that the steam dryers meet this design basis requires a calculation of the magnitude of the acoustic load on the steam dryer during a MSLB. In this study, computational fluid dynamics (CFD) is used as an alternate calculation method to investigate the phenomenon of MSLB. Transient simulations with fine time steps were carried out. The results show that CFD is a useful tool to provide additional information on the acoustic load as compared to the traditional methods. From the CFD results, the minimum pressure value and its distribution area at different flow times was identified. Through the modeling, an understanding of the detailed transient flow field, particularly the acoustic pressure field near the dryer hood during the MSLB was achieved.     

Two-phase turbulent mixing and buoyancy drift in rod bundles

This paper describes the development of generalized relationships for single- and two-phase intersubchannel turbulent mixing in vertical and horizontal flows, and lateral buoyancy drift in horizontal flows.The relationships for turbulent mixing, together with a recommended one for void drift, have been implemented in a subchannel thermalhydraulics code, and assessed using a range of data on enthalpy migration in vertical steam–water flows under BWR and PWR diabatic conditions. The intent of this assessment was to optimize these relationships to give the best agreement with the enthalpy migration data for vertical flows. The optimized turbulent mixing relationships were then used as a basis to benchmark a proposed buoyancy drift model to give the best predictions of void and enthalpy migration data in horizontal flows typical of PHWR CANDU¹ reactor operation under normal and off-normal conditions.Overall, the optimized turbulent mixing and buoyancy drift relationships have been found to predict the available data quite well, and generally better and more consistently than currently used models. This is expected to result in more accurate calculations of subchannel distributions of phasic flows, and hence, in improved predictions of critical heat flux (CHF).     

Generation characteristics of a metal ion plasma jet in vacuum discharge

To improve the performance of a metal ion plasma jet in vacuum discharge,an anode-insulated cone-cylinder electrode with insulating sleeve is proposed in this paper.Discharge characteristics and generation characteristics of plasma of the electrode are investigated,effects of diameter of insulating sleeve,variety of cathode material and length of the insulating sleeve on characteristics of metal ion plasma jet are discussed.Results indicate that a directional and steady plasma jet is formed by using the novel electrode with insulating sleeve under high vacuum conditions.Moreover,the properties of metal ion plasma jet are improved by using the aluminum cathode and thin and long insulating sleeve.The study provides strong support for research of vacuum metal ion plasma thruster and ion implantation technology.     

Pressure distribution due to a steam bubble collapse in a BWR pressure suppression pool

At the end of a steam blowdown into a pressure suppression pool severe steam bubble collapses occur. Because the condenser pipes are arranged in a regular pattern. There are two limit loading conditions: the simultaneous collapse with the same intensity at each condenser pipe which is mostly discussed and a strong collapse only at one pipe or a small group of pipes. Since the influence of the wall elasticity on the pressure input due to the steam bubble collapse is unessentially small, the complete bubble-fluid-structure problem may be solved in two steps: First the bubble-fluid interaction in the rigid containment and later on the fluid-structure interaction in the elastic containment due to the pressure input by the bubble-fluid system. The bubble-fluid problem is also separated in the analysis for the pressure time history and that for the pressure distribution. If compressibility effects are neglected, the steam bubble surface is to idealize as a sphere during the whole collapse time for analyzing the pressure time history. The pressure distribution is fixed then by the potential due to a stationary source in the bubble center. The structure system is solved as usual in structural dynamics by the eigenmode response analysis in which the fluid-structure system is substituted by some uncoupled mass spring oscillators.The aim of this paper is to clarify the behaviour of the bubble-fluid-structure system by simple models and by separating the problem into the important global behaviour and some additional localized effects. This seems to be necessary since the pool geometry especially for steel containments and also the analysis for the structural system, which consists of different shells is too complicated. The influence of different idealisations on the pressure distribution in the rigid containment and on the eigenmodes of the fluid-structure system is presented. Formulas are derived for estimating the eigenfrequencies, the eigenmodes and the pressure distribution as well as the pressure reduction by the wall elasticity. Furthermore the similarity rules for transfering the results from small experimental facilities to the complex ring containment are explained. Analytical and experimental results are compared and how the maximum pressure distribution in steel and concrete containments differ from each other. Finally the following engineering concept for a stepwise analysis is confirmed: The resulting wall loads due to the fluid-pressure and the inertia wall forces may be determined by a simple fluid-structure model and then the stresses due to these loads by a refined shell analysis.     

Degradation of sulfamethoxazole in water by dielectric barrier discharge plasma jet: influencing parameters,degradation pathway,toxicity evaluation

Sulfamethoxazole (SMX) is an antibiotic and widely present in aquatic environments, so it presents a serious threat to human health and sustainable development. A dielectric barrier discharge (DBD) plasma jet was utilized to degrade aqueous SMX, and the effects of various operating parameters (working gas, discharge power, etc) on SMX degradation performance were studied. The experimental results showed that the DBD plasma jet can obtain a relatively high degradation efficiency for SMX when the discharge power is high with an oxygen atmosphere, the initial concentration of SMX is low, and the aqueous solution is under acidic conditions. The reactive species produced in the liquid phase were detected, and OH radicals and O3 were found to play a significant role in the degradation of SMX. Moreover, the process of SMX degradation could be better fitted by the quasi-first-order reaction kinetic equation. The analysis of the SMX degradation process indicated that SMX was gradually decomposed and 4-amino benzene sulfonic acid, benzene sulfonamide, 4-nitro SMX, and phenylsulfinyl acid were detected, and thus three possible degradation pathways were finally proposed. The mineralization degree of SMX reached 90.04% after plasma treatment for 20 min, and the toxicity of the solution fluctuated with the discharge time but eventually decreased.     

Analysis of heat transfer and flow characteristics in turbulent impinging jet

Jet impingement technique is characterized by a high heat removal capability. As such it has been proposed as a cooling method for the helium cooled divertor, a high-heat flux component of the future fusion reactor called DEMO. Since power plant efficiency depends on the divertor’s heat removal capability it has to meet certain demands, i.e. high-heat transfer and low pressure drop.In this paper local heat transfer and flow characteristics of an axis-symmetric impinging jet are analyzed numerically using the RANS approach and eddy viscosity type SST turbulence model. Turbulence models and heat transfer predictions are validated on the free jet impingement experiment (Baughn and Shimizu, 1989). Since the numerical results are affected by the turbulence model the influence of the turbulent production is investigated in particular.The validated numerical model is further applied to analyze the effect of the nozzle inlet shape on the heat removal capability and pressure drop in the confined impinging jet. Two different nozzle inlet parameters are tested; chamfer angle θ and chamfer depth L_ch. The numerical results are compared with the experimental data (Brignoni and Garimella, 2000).     

Low-quality and subcooled film boiling of water at elevated pressures

Stable film boiling heat transfer data have been obtained in an 8.9 mm ID tube at pressures from 2 to 9 MPa. These data were obtained at low-quality and subcooled conditions, over a mass flux range of 0.11 to 2.75 Mg m⁻² s⁻¹. Excessive film boiling surface temperatures were avoided by using the hot patch technique. Contrary to the high-quality data, the low-quality data showed a decrease in heat transfer coefficient with an increase in quality. The film boiling data were compared with existing film boiling correlations. None of these were found to be satisfactory.     

Spectroscopic measurements in Fe-plasma produced by pulsed laser ablation

A Nd:Yag laser is employed to irradiate thick Fe targets placed in vacuum. The obtained non-equilibrium Fe-plasma is investigated with various analytical techniques. An electrostatic ion energy analyzer (IEA) and a ring ion collector (ICR) are employed, in time-of-flight configuration, to monitor in situ the ejected ions from the plasma along the normal direction to the target surface and to determine the plasma core temperature and the ion energy distributions.The visible plasma emission, detected with an optical spectrometer, permitted to evaluate the electronic temperature and density and the fluence threshold of the visible light emission. The spectroscopy measurements of ions and photons and the fast CCD plasma images are employed to evaluate the temperature gradients in the laser-generated plasma plume.     

The investigation of OH radicals produced in a DC glow discharge by laser-induced fluorescence spectrometry

In this paper the OH radicals produced by a needle-plate negative DC discharge in water vapor,N2 + H2O mixture gas and He + H2O mixture gas are investigated by a laser-induced fluorescence (LIF) system.With a ballast resistor in the circuit,the discharge current is limited and the discharges remain in glow.The OH rotation temperature is obtained from fluorescence rotational branch fitting,and is about 350 K in pure water vapor.The effects of the discharge current and gas pressure on the production and quenching processes of OH radicals are investigated.The results show that in water vapor and He + H2O mixture gas the fluorescence intensity of OH stays nearly constant with increasing discharge current,and in N2 + H2O mixture gas the fluorescence intensity of OH increases with increasing discharge current.In water vapor and N2 + H2O mixture gas the fluorescence intensity of OH decreases with increasing gas pressure in the studied pressure range,and in He + H2O mixture gas the fluorescence intensity of OH shows a maximum value within the studied gas pressure range.The physicochemical reactions between electrons,radicals,ground and metastable molecules are discussed.The results in this work contribute to the optimization of plasma reactivity and the establishment of a molecule reaction dynamics model.